CN112014323A - Remote sensing inversion method and system for surface evapotranspiration - Google Patents

Remote sensing inversion method and system for surface evapotranspiration Download PDF

Info

Publication number
CN112014323A
CN112014323A CN202010895830.0A CN202010895830A CN112014323A CN 112014323 A CN112014323 A CN 112014323A CN 202010895830 A CN202010895830 A CN 202010895830A CN 112014323 A CN112014323 A CN 112014323A
Authority
CN
China
Prior art keywords
soil
temperature
vegetation
relative
limiting factor
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.)
Granted
Application number
CN202010895830.0A
Other languages
Chinese (zh)
Other versions
CN112014323B (en
Inventor
唐荣林
王晟力
姜亚珍
李召良
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Institute of Geographic Sciences and Natural Resources of CAS
Original Assignee
Institute of Geographic Sciences and Natural Resources of CAS
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Institute of Geographic Sciences and Natural Resources of CAS filed Critical Institute of Geographic Sciences and Natural Resources of CAS
Priority to CN202010895830.0A priority Critical patent/CN112014323B/en
Publication of CN112014323A publication Critical patent/CN112014323A/en
Application granted granted Critical
Publication of CN112014323B publication Critical patent/CN112014323B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N21/00Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
    • G01N21/17Systems in which incident light is modified in accordance with the properties of the material investigated
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01WMETEOROLOGY
    • G01W1/00Meteorology
    • G01W1/02Instruments for indicating weather conditions by measuring two or more variables, e.g. humidity, pressure, temperature, cloud cover or wind speed
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N21/00Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
    • G01N21/17Systems in which incident light is modified in accordance with the properties of the material investigated
    • G01N2021/1793Remote sensing
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N21/00Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
    • G01N21/17Systems in which incident light is modified in accordance with the properties of the material investigated
    • G01N2021/1793Remote sensing
    • G01N2021/1795Atmospheric mapping of gases

Landscapes

  • Life Sciences & Earth Sciences (AREA)
  • Environmental & Geological Engineering (AREA)
  • Immunology (AREA)
  • Pathology (AREA)
  • Analytical Chemistry (AREA)
  • Biochemistry (AREA)
  • General Health & Medical Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Physics & Mathematics (AREA)
  • Chemical & Material Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Health & Medical Sciences (AREA)
  • Atmospheric Sciences (AREA)
  • Biodiversity & Conservation Biology (AREA)
  • Ecology (AREA)
  • Environmental Sciences (AREA)
  • Radiation Pyrometers (AREA)

Abstract

The invention discloses a remote sensing inversion method of surface evapotranspiration, which comprises the following steps: acquiring meteorological data and remote sensing data of the earth surface; respectively calculating the soil temperature and the vegetation temperature relative to the soil moisture of 0 and the soil temperature and the vegetation temperature relative to the soil moisture of 1 according to the meteorological data and the remote sensing data; further calculating the soil dry edge temperature, the soil humidity limiting factor, the vegetation humidity limiting factor and the vegetation temperature limiting factor; and calculating the surface evaporation capacity by using a Priestley-Taylor formula according to the soil humidity limiting factor, the vegetation humidity limiting factor and the temperature limiting factor. According to the invention, the calculation limiting factors of the soil temperature and the vegetation temperature relative to the soil moisture of 0 and the calculation limiting factors of the soil temperature and the vegetation temperature relative to the soil moisture of 1 are introduced, so that the water change of the soil surface layer and the vegetation root area on the daily scale can be reflected, and the accuracy of remote sensing inversion of the earth surface evapotranspiration is improved by adopting the inversion method of daily scale earth surface evapotranspiration.

Description

Remote sensing inversion method and system for surface evapotranspiration
Technical Field
The invention relates to the technical field of surface evapotranspiration calculation, in particular to a remote sensing inversion method and system for surface evapotranspiration.
Background
The surface evapotranspiration is an important component of surface energy balance and water circulation, and is also an important index of vegetation growth condition and crop yield. The method has the advantages that regional or global-scale surface evapotranspiration can be accurately estimated, and the method has important significance for global climate evolution research, water resource planning and management, agricultural water conservation and crop yield estimation. However, the evapotranspiration has strong spatial variability, and the remote sensing technology has the characteristics of high speed, high space-time resolution and suitability for large-area long-term observation, and is considered as the most effective means for acquiring the evapotranspiration distribution of the earth surface area or the global scale with high time efficiency and high precision.
Due to the characteristics of instantaneous transit of the remote sensing satellite, the earth surface evapotranspiration obtained based on the thermal infrared remote sensing inversion model is generally an instantaneous value of the transit time of the satellite, and the change rule of the earth surface evapotranspiration in a daily scale and a longer time scale cannot be reflected. In practical application, the daily scale and longer time scale evapotranspiration are more valuable. The ratio of variables closely related to evapotranspiration under clear sky conditions to the evapotranspiration is always assumed to be basically unchanged in the daytime, and the time scale expansion is carried out on the instantaneous evapotranspiration obtained by the inversion of the satellite transit time to obtain the daily scale evapotranspiration. However, the daily scale expansion method is based on the instantaneous evapotranspiration of remote sensing inversion, and the accuracy and the application condition of different remote sensing inversion methods are different. Therefore, the estimation result of the earth surface evapotranspiration of the current day scale mostly has two uncertainties: (1) uncertainty of the remote sensing inversion model; (2) uncertainty of the daily scale extension method. In order to avoid the error accumulation of uncertainty of the two methods, a daily scale earth surface evapotranspiration remote sensing inversion model is directly constructed, and the method has important research significance for improving the accuracy of the daily scale earth surface evapotranspiration of remote sensing inversion.
Disclosure of Invention
The invention aims to provide a remote sensing inversion method and a remote sensing inversion system for ground evapotranspiration, so as to avoid error accumulation of uncertainty of a remote sensing inversion model and uncertainty of a daily scale expansion method.
In order to achieve the purpose, the invention provides the following scheme:
a remote sensing inversion method of ground surface evapotranspiration comprises the following steps:
acquiring meteorological data and remote sensing data of the earth surface; the meteorological data comprise air temperature, atmospheric pressure, air relative humidity, wind speed, uplink short-wave radiation, downlink short-wave radiation and the like; the remote sensing data comprises earth surface temperature, earth surface emissivity, reflectivity, leaf area index, enhanced vegetation index, photosynthetically active radiation and the like;
respectively calculating the soil temperature and the vegetation temperature relative to the soil moisture of 0 and the soil temperature and the vegetation temperature relative to the soil moisture of 1 according to the meteorological data and the remote sensing data of the earth surface;
calculating and distinguishing the soil dry edge temperature of the sufficient water state of the root zone and the water-free state of the soil surface layer according to the soil temperature of 0 relative to the soil moisture, the vegetation temperature of 1 relative to the soil moisture and remote sensing data;
calculating a soil humidity limiting factor, a vegetation humidity limiting factor and a vegetation temperature limiting factor according to the soil temperature and the vegetation temperature relative to the soil moisture of 0, the soil temperature and the vegetation temperature relative to the soil moisture of 1 and the soil dry edge temperature;
and calculating the surface evaporation capacity by using a Priestley-Taylor (Priesterli-Taylor) formula according to the soil humidity limiting factor, the vegetation humidity limiting factor and the vegetation temperature limiting factor.
Optionally, according to meteorological data and the remote sensing data of earth's surface respectively calculate soil temperature and vegetation temperature that relative soil moisture is 0, soil temperature and vegetation temperature that relative soil moisture is 1, specifically include:
using a formula based on the meteorological data and the remote sensing data
Figure BDA0002658402690000021
Calculating the soil temperature relative to the soil moisture of 0;
using a formula based on the meteorological data and the remote sensing data
Figure BDA0002658402690000022
Calculating the vegetation temperature relative to the soil moisture of 0;
using a formula based on the meteorological data and the remote sensing data
Figure BDA0002658402690000023
Calculating the soil temperature relative to the soil moisture of 1;
using a formula based on the meteorological data and the remote sensing data
Figure BDA0002658402690000031
Calculating the vegetation temperature relative to the soil moisture of 1;
wherein, Ts0And Tv0Respectively representing the soil temperature and vegetation temperature, T, relative to the soil moisture of 0s1And Tv1Respectively representing the soil temperature and the vegetation temperature relative to the soil moisture of 1; ρ is the air density; cpIs the specific heat at constant pressure; e.g. of the typesSaturated water vapor pressure; e.g. of the typeaThe actual water vapor pressure of the atmosphere; t isaThe temperature near the surface is obtained; gvwVegetation canopy conductivity sufficient for water supply; gavAnd gasThe aerodynamic conductivity of the vegetation and the upper layer of the soil are respectively; rn,sAnd Rn,vRespectively the soil component net radiation and the vegetation component net radiation; gsIs the soil heat flux.
Optionally, the soil dry edge temperature for distinguishing the sufficient water state of the root zone and the moisture-free state of the soil surface layer is calculated according to the soil temperature relative to the soil moisture of 0 and the vegetation temperature relative to the soil moisture of 1 and the remote sensing data, and specifically includes:
according to the surface remote sensing data, using a formula
Figure BDA0002658402690000032
Calculating the vegetation ratio Fv
Wherein EVI is Enhanced Vegetation Index (EVI)minAnd EVImaxRespectively the minimum enhancement type vegetation index corresponding to bare soil and the maximum enhancement type vegetation index corresponding to full vegetation coverage;
according to relative soilThe soil temperature with 0 soil moisture and the vegetation temperature and vegetation proportion with 1 relative to the soil moisture by using a formula
Figure BDA0002658402690000033
Calculating and distinguishing soil dry edge temperature T of sufficient water state and no water state on soil surface layer of root zone*
Wherein, Ts0And Tv1The soil temperature is 0 relative to the soil moisture and the vegetation temperature is 1 relative to the soil moisture.
Optionally, the soil temperature and vegetation temperature that is 0 according to relative soil moisture, relative soil moisture and the soil temperature and vegetation temperature that is 1 and soil dry edge temperature calculate soil humidity limiting factor, vegetation humidity limiting factor and vegetation temperature limiting factor, specifically include:
judging whether the surface temperature of the remote sensing image pixel is greater than the dry edge temperature of the soil or not to obtain a judgment result;
if the judgment result shows no, the formula f is usedVMCalculating a vegetation moisture limiting factor f 1VM(ii) a Using formulas
Figure BDA0002658402690000041
Calculating the soil humidity limiting factor fSM
If the judgment result shows yes, the formula f is utilizedSMCalculating the soil humidity limiting factor f as 0SM(ii) a Using formulas
Figure BDA0002658402690000042
Calculating vegetation humidity limiting factor fVM
Wherein, Ts0And Ts1Respectively, a soil temperature of 0 relative to the soil moisture and a soil temperature of 1 relative to the soil moisture, TsoilWhich is indicative of the temperature of the soil components,
Figure BDA0002658402690000043
TRfor remote-sensing the surface temperature of image elements, FvIs the vegetation proportion; t isv0And Tv1Respectively vegetation temperature of 0 relative to the soil moisture and vegetation temperature of 1 relative to the soil moisture, TvegIs the temperature of the components of the vegetation,
Figure BDA0002658402690000044
using formulas
Figure BDA0002658402690000045
Calculating a vegetation temperature limiting factor fT
Wherein, TmaxThe maximum atmospheric temperature in one day; t isoptThe temperature is the optimum temperature for vegetation growth.
Optionally, calculating the surface evaporation capacity by using a Priestley-Taylor formula according to the soil humidity limiting factor, the vegetation humidity limiting factor and the vegetation temperature limiting factor specifically includes:
according to the vegetation temperature limiting factor and the vegetation humidity limiting factor, using a formula
Figure BDA0002658402690000046
Calculating the transpiration amount of the vegetation;
according to the soil humidity limiting factor, using the formula
Figure BDA0002658402690000047
Calculating the evaporation capacity of the soil;
by using
Figure BDA0002658402690000048
Calculating the interception evaporation capacity;
according to the vegetation transpiration amount, the soil evaporation amount and the cut-off evaporation amount, utilizing a formula ET ═ ETv+ETs+ETiCalculating the evaporation capacity of the earth surface;
wherein, ETv、ETsAnd ETiThe vegetation transpiration amount, the soil evaporation amount and the closure evaporation amount are respectively, and ET is the earth surface evaporation amount; f. ofwetIs the proportion of the wet leaves of the vegetation,
Figure BDA0002658402690000051
RH is the relative humidity of the atmosphere; f. ofgThe ratio of vegetation to green leaves;
Figure BDA0002658402690000052
fAPARphotosynthetically active radiation absorbed for green plants; f. ofIPARPhotosynthetically active radiation intercepted for green vegetation; f. ofTIs a vegetation temperature limiting factor; f. ofSMIs a soil humidity limiting factor; f. ofVMIs vegetation humidity limiting factor; alpha is the coefficient of a Priestley-Taylor formula, and alpha is 1.26; e.g. of the typesSaturated water vapor pressure; t isaThe temperature near the surface is obtained; rn,vFor canopy net radiation, Rn,sNet radiation for soil; gsIs the soil heat flux.
A remote sensing inversion system of surface evapotranspiration, the inversion system comprising, for example:
the data acquisition module is used for acquiring meteorological data and remote sensing data of the earth surface; the meteorological data comprise air temperature, atmospheric pressure, air relative humidity, wind speed, uplink short-wave radiation, downlink short-wave radiation and the like; the remote sensing data comprises earth surface temperature, earth surface emissivity, reflectivity, leaf area index, enhanced vegetation index, photosynthetically active radiation and the like;
the soil temperature and vegetation temperature calculation module is used for calculating the soil temperature and vegetation temperature relative to the soil moisture of 0 and the soil temperature and vegetation temperature relative to the soil moisture of 1 respectively according to the meteorological data and the remote sensing data of the earth surface;
the soil dry edge temperature calculation module is used for calculating and distinguishing the soil dry edge temperature of the sufficient water state of the root area and the no water state of the soil surface layer according to the soil temperature relative to the soil water content of 0, the vegetation temperature relative to the soil water content of 1 and remote sensing data;
the limiting factor calculation module is used for calculating a soil humidity limiting factor, a vegetation humidity limiting factor and a vegetation temperature limiting factor according to the soil temperature and the vegetation temperature relative to the soil moisture of 0, the soil temperature and the vegetation temperature relative to the soil moisture of 1 and the soil dry edge temperature;
and the earth surface evaporation amount calculation module is used for calculating the earth surface evaporation amount by utilizing a Priestley-Taylor formula according to the soil humidity limiting factor, the vegetation humidity limiting factor and the temperature limiting factor.
Optionally, the soil temperature and vegetation temperature calculation module specifically includes:
a soil temperature calculation submodule with a relative soil moisture of 0 for utilizing a formula according to the meteorological data and the remote sensing data
Figure BDA0002658402690000061
Calculating the soil temperature relative to the soil moisture of 0;
a vegetation temperature calculation submodule with a relative soil moisture of 0 for utilizing a formula according to the meteorological data and the remote sensing data
Figure BDA0002658402690000062
Calculating the vegetation temperature relative to the soil moisture of 0;
a soil temperature calculation submodule with a relative soil moisture of 1 for utilizing a formula according to the meteorological data and the remote sensing data
Figure BDA0002658402690000063
Calculating the soil temperature relative to the soil moisture of 1;
a vegetation temperature calculation submodule with relative soil moisture of 1 for utilizing a formula according to the meteorological data and the remote sensing data
Figure BDA0002658402690000064
Calculating the vegetation temperature relative to the soil moisture of 1;
wherein, Ts0And Tv0Respectively representing the soil temperature and vegetation temperature, T, relative to the soil moisture of 0s1And Tv1Respectively representing the soil temperature and the vegetation temperature relative to the soil moisture of 1; ρ is the air density; cpIs the specific heat at constant pressure; e.g. of the typesSaturated water vapor pressure; e.g. of the typeaThe actual water vapor pressure of the atmosphere; t isaThe temperature near the surface is obtained; gvwVegetation canopy conductivity sufficient for water supply;gavAnd gasThe aerodynamic conductivity of the vegetation and the upper layer of the soil are respectively; rn,sAnd Rn,vRespectively the soil component net radiation and the vegetation component net radiation; gsIs the soil heat flux.
Optionally, the soil dry edge temperature calculation module specifically includes:
a vegetation coverage calculation submodule for utilizing a formula according to the earth surface remote sensing data
Figure BDA0002658402690000065
Calculating the vegetation ratio Fv
Wherein EVI is Enhanced Vegetation Index (EVI)minAnd EVImaxRespectively the minimum enhancement type vegetation index corresponding to bare soil and the maximum enhancement type vegetation index corresponding to full vegetation coverage;
a soil dry edge temperature operator module for utilizing a formula according to the soil temperature of 0 relative to the soil moisture and the vegetation temperature and vegetation proportion of 1 relative to the soil moisture
Figure BDA0002658402690000071
Calculating and distinguishing soil dry edge temperature T of sufficient water state and no water state on soil surface layer of root zone*
Wherein, Ts0And Tv1The soil temperature is 0 relative to the soil moisture and the vegetation temperature is 1 relative to the soil moisture.
Optionally, the limiting factor calculating module specifically includes:
the judgment submodule is used for judging whether the surface temperature of the remote sensing image pixel is greater than the dry edge temperature of the soil or not to obtain a judgment result;
a first limiting factor calculating submodule for utilizing the formula f if the judgment result shows noVMCalculating vegetation humidity limiting factor f 1VM(ii) a Using formulas
Figure BDA0002658402690000072
Calculating the soil humidity limiting factor fSM
A second limiting factor calculating submodule for using the formula f if the judgment result shows yesSMCalculating a vegetation humidity limiting factor fVM(ii) a Using formulas
Figure BDA0002658402690000073
Calculating vegetation humidity limiting factor fVM
Wherein, Ts0And Ts1The soil temperature is 0 relative to the soil moisture and the soil temperature is 1 relative to the soil moisture; t issoilWhich is indicative of the temperature of the soil components,
Figure BDA0002658402690000074
TRfor remote-sensing the surface temperature of image elements, FvIs the vegetation proportion; t isv0And Tv1Respectively vegetation temperature of 0 relative to the soil moisture and vegetation temperature of 1 relative to the soil moisture, TvegIs the temperature of the components of the vegetation,
Figure BDA0002658402690000075
a third limiting factor calculation submodule for utilizing a formula
Figure BDA0002658402690000076
Calculating a vegetation temperature limiting factor fT
Wherein, TmaxThe maximum atmospheric temperature in one day; t isoptThe temperature is the optimum temperature for vegetation growth.
Optionally, the surface evaporation amount calculation module specifically includes:
a vegetation transpiration meter operator module for utilizing a formula according to the vegetation temperature limiting factor and the vegetation humidity limiting factor
Figure BDA0002658402690000081
Calculating the transpiration amount of the vegetation;
a soil evaporation meter operator module for utilizing a formula according to the soil humidity limiting factor
Figure BDA0002658402690000082
Calculating the evaporation capacity of the soil;
a cut-off evaporation amount calculation operator module for utilizing
Figure BDA0002658402690000083
Calculating the interception evaporation capacity;
the earth surface evapotranspiration meter operator module is used for utilizing a formula ET-ET according to the vegetation evapotranspiration, the soil evaporation and the interception evaporationv+ETs+ETiCalculating the evapotranspiration amount of the earth surface;
wherein, ETv、ETsAnd ETiThe vegetation transpiration amount, the soil evaporation amount and the closure evaporation amount are respectively, and ET is the earth surface evaporation amount; f. ofwetIs the proportion of the wet leaves of the vegetation,
Figure BDA0002658402690000084
RH is the relative humidity of the atmosphere; f. ofgThe ratio of vegetation to green leaves;
Figure BDA0002658402690000085
fAPARphotosynthetically active radiation absorbed for green plants; f. ofIPARPhotosynthetically active radiation intercepted for green vegetation; f. ofTIs a vegetation temperature limiting factor; f. ofSMIs a soil humidity limiting factor; f. ofVMIs vegetation humidity limiting factor; alpha is the coefficient of a Priestley-Taylor formula, and alpha is 1.26; e.g. of the typesSaturated water vapor pressure; t isaThe temperature near the surface is obtained; rn,vFor canopy net radiation, Rn,sNet radiation for soil; gsIs the soil heat flux.
According to the specific embodiment provided by the invention, the invention discloses the following technical effects:
the invention discloses a remote sensing inversion method of surface evapotranspiration, which comprises the following steps: acquiring meteorological data and remote sensing data of the earth surface; respectively calculating the soil temperature and the vegetation temperature relative to the soil moisture of 0 and the soil temperature and the vegetation temperature relative to the soil moisture of 1 according to the meteorological data and the remote sensing data of the earth surface; calculating and distinguishing the soil dry edge temperature of the sufficient water state of the root zone and the water-free state of the soil surface layer according to the soil temperature of 0 relative to the soil moisture, the vegetation temperature of 1 relative to the soil moisture and the remote sensing data of the earth surface; calculating a soil humidity limiting factor, a vegetation humidity limiting factor and a vegetation temperature limiting factor according to the soil temperature and the vegetation temperature relative to the soil moisture of 0, the soil temperature and the vegetation temperature relative to the soil moisture of 1 and the soil dry edge temperature; and calculating the surface evaporation capacity by using a Priestley-Taylor formula according to the soil humidity limiting factor, the vegetation humidity limiting factor and the temperature limiting factor. According to the invention, the calculation limiting factors of the soil temperature and the vegetation temperature relative to the soil moisture of 0 and the calculation limiting factors of the soil temperature and the vegetation temperature relative to the soil moisture of 1 are introduced, so that the water change of the soil surface layer and the vegetation root area on the daily scale can be reflected, and the accuracy of remote sensing inversion of the earth surface evapotranspiration is improved by adopting the inversion method of daily scale earth surface evapotranspiration.
Drawings
In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the drawings needed to be used in the embodiments will be briefly described below, and it is obvious that the drawings in the following description are only some embodiments of the present invention, and it is obvious for those skilled in the art to obtain other drawings without inventive exercise.
Fig. 1 is a flowchart of a remote sensing inversion method of surface evapotranspiration provided by the present invention.
Detailed Description
The invention aims to provide a remote sensing inversion method and a remote sensing inversion system for ground evapotranspiration, so as to avoid error accumulation of uncertainty of a remote sensing inversion model and uncertainty of a daily scale expansion method.
In order to make the aforementioned objects, features and advantages of the present invention comprehensible, embodiments accompanied with figures are described in further detail below.
As shown in fig. 1, the present invention provides a remote sensing inversion method of surface evapotranspiration, which comprises the following steps:
step 101, acquiring meteorological data and remote sensing data of the earth surface; the meteorological data comprise air temperature, atmospheric pressure, air relative humidity, wind speed, uplink short-wave radiation, downlink short-wave radiation and the like; the remote sensing data comprises earth surface temperature, earth surface emissivity, reflectivity, leaf area index, enhanced vegetation index, photosynthetically active radiation and the like; and constructing an input data set of the daily scale earth surface evapotranspiration remote sensing inversion method.
And 102, respectively calculating the soil temperature and the vegetation temperature relative to the soil moisture of 0 and the soil temperature and the vegetation temperature relative to the soil moisture of 1 according to the meteorological data and the remote sensing data of the earth surface.
102, respectively calculating the soil temperature and the vegetation temperature of 0 relative to the soil moisture and the soil temperature and the vegetation temperature of 1 relative to the soil moisture according to the meteorological data and the remote sensing data of the earth surface, and specifically comprising the following steps of: using a formula based on the meteorological data and the remote sensing data
Figure BDA0002658402690000101
Calculating the soil temperature relative to the soil moisture of 0; using a formula based on the meteorological data and the remote sensing data
Figure BDA0002658402690000102
Calculating the vegetation temperature relative to the soil moisture of 0; using a formula based on the meteorological data and the remote sensing data
Figure BDA0002658402690000103
Calculating the soil temperature relative to the soil moisture of 1; using a formula based on the meteorological data and the remote sensing data
Figure BDA0002658402690000104
Calculating the vegetation temperature relative to the soil moisture of 1; wherein, Ts0And Tv0Respectively representing the soil temperature and vegetation temperature, T, relative to the soil moisture of 0s1And Tv1Respectively representing the temperature of the soil and the vegetation relative to the water content of the soil of 1(ii) temperature; ρ is the air density; cpIs the specific heat at constant pressure; e.g. of the typesSaturated water vapor pressure; e.g. of the typeaThe actual water vapor pressure of the atmosphere; t isaThe temperature near the surface is obtained; gvwVegetation canopy conductivity sufficient for water supply; gavAnd gasThe aerodynamic conductivity of the vegetation and the upper layer of the soil are respectively; rn,sAnd Rn,vRespectively the soil component net radiation and the vegetation component net radiation; gsIs the soil heat flux.
103, calculating and distinguishing the soil dry edge temperature of the sufficient water state of the root area and the water-free state of the soil surface layer according to the soil temperature relative to the soil water content of 0, the vegetation temperature relative to the soil water content of 1 and the remote sensing data of the earth surface;
step 103, calculating and distinguishing the soil dry edge temperature of the sufficient water state of the root zone and the water-free state of the soil surface layer according to the soil temperature relative to the soil water content of 0, the vegetation temperature relative to the soil water content of 1 and the remote sensing data of the earth surface, and specifically comprising the following steps: according to the surface remote sensing data, using a formula
Figure BDA0002658402690000105
Calculating the vegetation ratio Fv(ii) a Wherein EVI is Enhanced Vegetation Index (EVI)minAnd EVImaxRespectively the minimum enhancement type vegetation index corresponding to bare soil and the maximum enhancement type vegetation index corresponding to full vegetation coverage; according to the soil temperature with the relative soil moisture of 0 and the vegetation temperature with the relative soil moisture of 1 and the vegetation proportion, utilizing a formula
Figure BDA0002658402690000111
Calculating and distinguishing soil dry edge temperature T of sufficient water state and no water state on soil surface layer of root zone*(ii) a Wherein, Ts0And Tv1The soil temperature is 0 relative to the soil moisture and the vegetation temperature is 1 relative to the soil moisture.
And 104, calculating a soil humidity limiting factor, a vegetation humidity limiting factor and a vegetation temperature limiting factor according to the soil temperature and the vegetation temperature relative to the soil moisture of 0, the soil temperature and the vegetation temperature relative to the soil moisture of 1 and the soil dry edge temperature.
Step 104 the soil humidity limiting factor, vegetation humidity limiting factor and vegetation temperature limiting factor are calculated according to the soil temperature and vegetation temperature relative to the soil moisture of 0, the soil temperature and vegetation temperature relative to the soil moisture of 1 and the soil dry edge temperature, and specifically include: judging whether the surface temperature of the remote sensing image pixel is greater than the dry edge temperature of the soil or not to obtain a judgment result; if the judgment result shows no, the formula f is usedVMCalculating a vegetation moisture limiting factor f 1VM(ii) a Using formulas
Figure BDA0002658402690000112
Calculating the soil humidity limiting factor fSM(ii) a If the judgment result shows yes, the formula f is utilizedSMCalculating the soil humidity limiting factor f as 0SM(ii) a Using formulas
Figure BDA0002658402690000113
Calculating vegetation humidity limiting factor fVM(ii) a Wherein, Ts0And Ts1Respectively, a soil temperature of 0 relative to the soil moisture and a soil temperature of 1 relative to the soil moisture, TsoilWhich is indicative of the temperature of the soil components,
Figure BDA0002658402690000114
TRfor remote-sensing the surface temperature of image elements, FvIs the vegetation proportion; t isv0And Tv1Respectively vegetation temperature of 0 relative to the soil moisture and vegetation temperature of 1 relative to the soil moisture, TvegIs the temperature of the components of the vegetation,
Figure BDA0002658402690000115
using formulas
Figure BDA0002658402690000116
Calculating a vegetation temperature limiting factor fT(ii) a Wherein, TmaxThe maximum atmospheric temperature in one day; t isoptThe temperature is the optimum temperature for vegetation growth.
And 105, calculating the surface evaporation capacity by using a Priestley-Taylor formula according to the soil humidity limiting factor, the vegetation humidity limiting factor and the temperature limiting factor.
105, calculating the surface evaporation capacity by using a Priestley-Taylor formula according to the soil humidity limiting factor, the vegetation humidity limiting factor and the temperature limiting factor, wherein the method specifically comprises the following steps: according to the vegetation temperature limiting factor and the vegetation humidity limiting factor, using a formula
Figure BDA0002658402690000121
Calculating the transpiration amount of the vegetation; according to the soil humidity limiting factor, using the formula
Figure BDA0002658402690000122
Calculating the evaporation capacity of the soil; by using
Figure BDA0002658402690000123
Calculating the interception evaporation capacity; according to the vegetation transpiration amount, the soil evaporation amount and the cut-off evaporation amount, utilizing a formula ET ═ ETv+ETs+ETiCalculating the evaporation capacity of the earth surface; wherein, ETv、ETsAnd ETiThe vegetation transpiration amount, the soil evaporation amount and the closure evaporation amount are respectively, and ET is the earth surface evaporation amount; f. ofwetIs the proportion of the wet leaves of the vegetation,
Figure BDA0002658402690000124
RH is the relative humidity of the atmosphere; f. ofgThe ratio of vegetation to green leaves;
Figure BDA0002658402690000125
fAPARphotosynthetically active radiation absorbed for green plants; f. ofIPARPhotosynthetically active radiation intercepted for green vegetation; f. ofTIs a vegetation temperature limiting factor; f. ofSMIs a soil humidity limiting factor; f. ofVMIs vegetation humidity limiting factor; alpha is the coefficient of a Priestley-Taylor formula, and alpha is 1.26; e.g. of the typesSaturated water vapor pressure; t isaThe temperature near the surface is obtained; rn,vFor canopy net radiation, Rn,sNet radiation for soil; gsFor soil heat ventilationAmount of the compound (A).
A remote sensing inversion system of surface evapotranspiration, the inversion system comprising:
the data acquisition module is used for acquiring meteorological data and remote sensing data of the earth surface; the meteorological data comprise air temperature, atmospheric pressure, air relative humidity, wind speed, uplink short-wave radiation, downlink short-wave radiation and the like; the remote sensing data comprises earth surface temperature, earth surface emissivity, reflectivity, leaf area index, enhanced vegetation index, photosynthetically active radiation and the like;
the soil temperature and vegetation temperature calculation module is used for calculating the soil temperature and vegetation temperature relative to the soil moisture of 0 and the soil temperature and vegetation temperature relative to the soil moisture of 1 respectively according to the meteorological data and the remote sensing data of the earth surface;
soil temperature and vegetation temperature calculation module specifically includes: a soil temperature calculation submodule with a relative soil moisture of 0 for utilizing a formula according to the meteorological data and the remote sensing data
Figure BDA0002658402690000131
Calculating the soil temperature relative to the soil moisture of 0; a vegetation temperature calculation submodule with a relative soil moisture of 0 for utilizing a formula according to the meteorological data and the remote sensing data
Figure BDA0002658402690000132
Calculating the vegetation temperature relative to the soil moisture of 0; a soil temperature calculation submodule with a relative soil moisture of 1 for utilizing a formula according to the meteorological data and the remote sensing data
Figure BDA0002658402690000133
Calculating the soil temperature relative to the soil moisture of 1; a vegetation temperature calculation submodule with relative soil moisture of 1 for utilizing a formula according to the meteorological data and the remote sensing data
Figure BDA0002658402690000134
Calculating the vegetation temperature relative to the soil moisture of 1; wherein, Ts0And Tv0Respectively representing the soil temperature and vegetation temperature, T, relative to the soil moisture of 0s1And Tv1Respectively representing the soil temperature and the vegetation temperature relative to the soil moisture of 1; ρ is the air density; cpIs the specific heat at constant pressure; e.g. of the typesSaturated water vapor pressure; e.g. of the typeaThe actual water vapor pressure of the atmosphere; t isaThe temperature near the surface is obtained; gvwVegetation canopy conductivity sufficient for water supply; gavAnd gasThe aerodynamic conductivity of the vegetation and the upper layer of the soil are respectively; rn,sAnd Rn,vRespectively the soil component net radiation and the vegetation component net radiation; gsIs the soil heat flux.
And the soil dry edge temperature calculation module is used for calculating and distinguishing the soil dry edge temperature of the sufficient water state of the root area and the soil surface layer water-free state according to the soil temperature relative to the soil water content of 0, the vegetation temperature relative to the soil water content of 1 and remote sensing data.
The soil dry edge temperature calculation module specifically comprises: a vegetation coverage calculation submodule for utilizing a formula according to the earth surface remote sensing data
Figure BDA0002658402690000135
Calculating the vegetation ratio Fv(ii) a Wherein EVI is Enhanced Vegetation Index (EVI)minAnd EVImaxRespectively the minimum enhancement type vegetation index corresponding to bare soil and the maximum enhancement type vegetation index corresponding to full vegetation coverage; a soil dry edge temperature operator module for utilizing a formula according to the soil temperature of 0 relative to the soil moisture and the vegetation temperature and vegetation proportion of 1 relative to the soil moisture
Figure BDA0002658402690000136
Calculating and distinguishing soil dry edge temperature T of sufficient water state and no water state on soil surface layer of root zone*(ii) a Wherein, Ts0And Tv1The soil temperature is 0 relative to the soil moisture and the vegetation temperature is 1 relative to the soil moisture.
And the limiting factor calculation module is used for calculating a soil humidity limiting factor, a vegetation humidity limiting factor and a vegetation temperature limiting factor according to the soil temperature and the vegetation temperature relative to the soil moisture of 0 and the soil temperature and the vegetation temperature relative to the soil moisture of 1.
The limiting factor calculating module specifically comprises: the judgment submodule is used for judging whether the surface temperature of the remote sensing image pixel is greater than the dry edge temperature of the soil or not to obtain a judgment result; a first limiting factor calculating submodule for utilizing the formula f if the judgment result shows noVMCalculating vegetation humidity limiting factor f 1VM(ii) a Using formulas
Figure BDA0002658402690000141
Calculating the soil humidity limiting factor fSM(ii) a A second limiting factor calculating submodule for using the formula f if the judgment result shows yesSMCalculating a vegetation humidity limiting factor fVM(ii) a Using formulas
Figure BDA0002658402690000142
Calculating vegetation humidity limiting factor fVM(ii) a Wherein, Ts0And Ts1The soil temperature is 0 relative to the soil moisture and the soil temperature is 1 relative to the soil moisture; t issoilWhich is indicative of the temperature of the soil components,
Figure BDA0002658402690000143
TRfor remote-sensing the surface temperature of image elements, FvIs the vegetation proportion; t isv0And Tv1Respectively vegetation temperature of 0 relative to the soil moisture and vegetation temperature of 1 relative to the soil moisture, TvegIs the temperature of the components of the vegetation,
Figure BDA0002658402690000144
a third limiting factor calculation submodule for utilizing a formula
Figure BDA0002658402690000145
Calculating a vegetation temperature limiting factor fT(ii) a Wherein, TmaxThe maximum atmospheric temperature in one day; t isoptFor vegetationLong optimum temperature.
And the earth surface evaporation amount calculation module is used for calculating the earth surface evaporation amount by utilizing a Priestley-Taylor formula according to the soil humidity limiting factor, the vegetation humidity limiting factor and the temperature limiting factor.
The surface evaporation amount calculation module specifically comprises: a vegetation transpiration meter operator module for utilizing a formula according to the vegetation temperature limiting factor and the vegetation humidity limiting factor
Figure BDA0002658402690000146
Calculating the transpiration amount of the vegetation; a soil evaporation meter operator module for utilizing a formula according to the soil humidity limiting factor
Figure BDA0002658402690000151
Calculating the evaporation capacity of the soil; a cut-off evaporation amount calculation operator module for utilizing
Figure BDA0002658402690000152
Calculating the interception evaporation capacity; the earth surface evapotranspiration meter operator module is used for utilizing a formula ET-ET according to the vegetation evapotranspiration, the soil evaporation and the interception evaporationv+ETs+ETiCalculating the evapotranspiration amount of the earth surface; wherein, ETv、ETsAnd ETiThe vegetation transpiration amount, the soil evaporation amount and the closure evaporation amount are respectively, and ET is the earth surface evaporation amount; f. ofwetIs the proportion of the wet leaves of the vegetation,
Figure BDA0002658402690000153
RH is the relative humidity of the atmosphere; f. ofgThe ratio of vegetation to green leaves;
Figure BDA0002658402690000154
fAPARphotosynthetically active radiation absorbed for green plants; f. ofIPARPhotosynthetically active radiation intercepted for green vegetation; f. ofTIs a vegetation temperature limiting factor; f. ofSMIs a soil humidity limiting factor; f. ofVMIs vegetation humidity limiting factor; alpha is the coefficient of a Priestley-Taylor formula, and alpha is 1.26; e.g. of the typesIs saturated waterSteam pressure; t isaThe temperature near the surface is obtained; rn,vFor canopy net radiation, Rn,sNet radiation for soil; gsIs the soil heat flux.
The invention provides a remote sensing inversion method of daily scale earth surface evapotranspiration, and belongs to the technical field of earth surface evapotranspiration remote sensing inversion. The method comprises the following steps: (A) collecting meteorological data and remote sensing data, and constructing an input data set of a daily scale earth surface evapotranspiration remote sensing inversion method; (B) calculating a soil humidity limiting factor and a vegetation humidity limiting factor by calculating the soil temperature and the vegetation temperature of 0 dry relative to the soil moisture, the soil temperature and the vegetation temperature and the soil dry edge temperature of 1 relative to the soil moisture; (C) calculating vegetation green leaf proportion, vegetation wet leaf proportion and vegetation temperature limiting factor according to the vegetation index, relative humidity, atmospheric temperature and vegetation growth optimum temperature; (D) and calculating vegetation transpiration, soil evaporation and cut-off evaporation according to a Priestley-Taylor formula, a vegetation green leaf proportion, a vegetation wet leaf proportion, a vegetation temperature limiting factor, a vegetation humidity limiting factor and a soil humidity limiting factor. The invention can reflect the water change of the soil surface layer and the vegetation root area on the daily scale by calculating the limiting factors of the soil temperature and the vegetation temperature relative to the soil moisture of 0 and the soil temperature and the vegetation temperature relative to the soil moisture of 1, directly provides the inversion method of daily scale surface evapotranspiration, and has important effect on improving the accuracy of remote sensing inversion of the surface evapotranspiration.
The embodiments in the present description are described in a progressive manner, each embodiment focuses on differences from other embodiments, and the same and similar parts among the embodiments are referred to each other.
The principle and the implementation manner of the present invention are explained by applying specific examples, the above description of the embodiments is only used to help understanding the method of the present invention and the core idea thereof, the described embodiments are only a part of the embodiments of the present invention, not all embodiments, and all other embodiments obtained by a person of ordinary skill in the art based on the embodiments of the present invention without creative efforts belong to the protection scope of the present invention.

Claims (10)

1. A remote sensing inversion method of ground surface evapotranspiration is characterized by comprising the following steps:
acquiring meteorological data and remote sensing data of the earth surface;
respectively calculating the soil temperature and the vegetation temperature relative to the soil moisture of 0 and the soil temperature and the vegetation temperature relative to the soil moisture of 1 according to the meteorological data and the remote sensing data of the earth surface;
calculating and distinguishing the soil dry edge temperature of the sufficient water state of the root zone and the water-free state of the soil surface layer according to the soil temperature of 0 relative to the soil moisture, the vegetation temperature of 1 relative to the soil moisture and remote sensing data;
calculating a soil humidity limiting factor, a vegetation humidity limiting factor and a vegetation temperature limiting factor according to the soil temperature and the vegetation temperature relative to the soil moisture of 0, the soil temperature and the vegetation temperature relative to the soil moisture of 1 and the soil dry edge temperature;
and calculating the surface evaporation capacity by using a Priestley-Taylor formula according to the soil humidity limiting factor, the vegetation humidity limiting factor and the vegetation temperature limiting factor.
2. The remote sensing inversion method for surface evapotranspiration according to claim 1, wherein the method for calculating the soil temperature and vegetation temperature relative to the soil moisture of 0 and the soil temperature and vegetation temperature relative to the soil moisture of 1 respectively according to the meteorological data and the remote sensing data of the surface specifically comprises the following steps:
using a formula based on the meteorological data and the remote sensing data
Figure FDA0002658402680000011
Calculating the soil temperature relative to the soil moisture of 0;
using a formula based on the meteorological data and the remote sensing data
Figure FDA0002658402680000012
Calculating relative soil moistureA vegetation temperature of 0;
using a formula based on the meteorological data and the remote sensing data
Figure FDA0002658402680000013
Calculating the soil temperature relative to the soil moisture of 1;
using a formula based on the meteorological data and the remote sensing data
Figure FDA0002658402680000021
Calculating the vegetation temperature relative to the soil moisture of 1;
wherein, Ts0And Tv0Respectively representing the soil temperature and vegetation temperature, T, relative to the soil moisture of 0s1And Tv1Respectively representing the soil temperature and the vegetation temperature relative to the soil moisture of 1; ρ is the air density; cpIs the specific heat at constant pressure; e.g. of the typesSaturated water vapor pressure; e.g. of the typeaThe actual water vapor pressure of the atmosphere; t isaThe temperature near the surface is obtained; gvwVegetation canopy conductivity sufficient for water supply; gavAnd gasThe aerodynamic conductivity of the vegetation and the upper layer of the soil are respectively; rn,sAnd Rn,vRespectively the soil component net radiation and the vegetation component net radiation; gsIs the soil heat flux.
3. The remote sensing inversion method for ground surface evapotranspiration according to claim 1, wherein the soil dry edge temperature for distinguishing the sufficient state of the water content in the root zone and the no water content state of the soil surface layer is calculated according to the soil temperature with the relative soil water content of 0, the vegetation temperature with the relative soil water content of 1 and remote sensing data, and specifically comprises the following steps:
according to the surface remote sensing data, using a formula
Figure FDA0002658402680000022
Calculating the vegetation ratio Fv
Wherein EVI is Enhanced Vegetation Index (EVI)minAnd EVImaxMinimum enhanced vegetation respectively corresponding to bare soilAn index and a maximum enhanced vegetation index corresponding to total vegetation coverage;
according to the soil temperature with the relative soil moisture of 0 and the vegetation temperature with the relative soil moisture of 1 and the vegetation proportion, utilizing a formula
Figure FDA0002658402680000023
Calculating and distinguishing soil dry edge temperature T of sufficient water state and no water state on soil surface layer of root zone*
Wherein, Ts0And Tv1The soil temperature is 0 relative to the soil moisture and the vegetation temperature is 1 relative to the soil moisture.
4. The remote sensing inversion method of surface evapotranspiration according to claim 1, wherein the calculation of the soil humidity limiting factor, the vegetation humidity limiting factor and the vegetation temperature limiting factor according to the soil temperature and vegetation temperature relative to a soil moisture of 0, the soil temperature and vegetation temperature relative to a soil moisture of 1 and the soil dry edge temperature specifically comprises:
judging whether the surface temperature of the remote sensing image pixel is greater than the dry edge temperature of the soil or not to obtain a judgment result;
if the judgment result shows no, the formula f is usedVMCalculating a vegetation moisture limiting factor f 1VM(ii) a Using formulas
Figure FDA0002658402680000031
Calculating the soil humidity limiting factor fSM
If the judgment result shows yes, the formula f is utilizedSMCalculating the soil humidity limiting factor f as 0SM(ii) a Using formulas
Figure FDA0002658402680000032
Calculating vegetation humidity limiting factor fVM
Wherein, Ts0And Ts1Respectively, a soil temperature of 0 relative to the soil moisture and a soil temperature of 1 relative to the soil moisture, TsoilWhich is indicative of the temperature of the soil components,
Figure FDA0002658402680000033
TRfor remote-sensing the surface temperature of image elements, FvIs the vegetation proportion; t isv0And Tv1Respectively vegetation temperature of 0 relative to the soil moisture and vegetation temperature of 1 relative to the soil moisture, TvegIs the temperature of the components of the vegetation,
Figure FDA0002658402680000034
using formulas
Figure FDA0002658402680000035
Calculating a vegetation temperature limiting factor fT
Wherein, TmaxThe maximum atmospheric temperature in one day; t isoptThe temperature is the optimum temperature for vegetation growth.
5. The remote sensing inversion method of surface evapotranspiration according to claim 1, wherein the calculation of the surface evapotranspiration by using a Priestley-Taylor formula according to the soil humidity limiting factor, the vegetation humidity limiting factor and the vegetation temperature limiting factor specifically comprises:
according to the vegetation temperature limiting factor and the vegetation humidity limiting factor, using a formula
Figure FDA0002658402680000036
Calculating the transpiration amount of the vegetation;
according to the soil humidity limiting factor, using the formula
Figure FDA0002658402680000037
Calculating the evaporation capacity of the soil;
by using
Figure FDA0002658402680000038
Calculating the interception evaporation capacity;
steaming according to vegetationThe evaporation capacity, the soil evaporation capacity and the interception evaporation capacity are calculated according to the formula ET ═ ETv+ETs+ETiCalculating the evaporation capacity of the earth surface;
wherein, ETv、ETsAnd ETiThe vegetation transpiration amount, the soil evaporation amount and the closure evaporation amount are respectively, and ET is the earth surface evaporation amount; f. ofwetIs the proportion of the wet leaves of the vegetation,
Figure FDA0002658402680000041
RH is the relative humidity of the atmosphere; f. ofgThe ratio of vegetation to green leaves;
Figure FDA0002658402680000042
fAPARphotosynthetically active radiation absorbed for green plants; f. ofIPARPhotosynthetically active radiation intercepted for green vegetation; f. ofTIs a vegetation temperature limiting factor; f. ofSMIs a soil humidity limiting factor; f. ofVMIs vegetation humidity limiting factor; alpha is the coefficient of a Priestley-Taylor formula, and alpha is 1.26; e.g. of the typesSaturated water vapor pressure; t isaThe temperature near the surface is obtained; rn,vFor canopy net radiation, Rn,sNet radiation for soil; gsIs the soil heat flux.
6. A remote sensing inversion system of surface evapotranspiration, the inversion system comprising:
the data acquisition module is used for acquiring meteorological data and remote sensing data of the earth surface;
the soil temperature and vegetation temperature calculation module is used for calculating the soil temperature and vegetation temperature relative to the soil moisture of 0 and the soil temperature and vegetation temperature relative to the soil moisture of 1 respectively according to the meteorological data and the remote sensing data of the earth surface;
the soil dry edge temperature calculation module is used for calculating and distinguishing the soil dry edge temperature of the sufficient water state of the root area and the no water state of the soil surface layer according to the soil temperature relative to the soil water content of 0, the vegetation temperature relative to the soil water content of 1 and remote sensing data;
the limiting factor calculation module is used for calculating a soil humidity limiting factor, a vegetation humidity limiting factor and a vegetation temperature limiting factor according to the soil temperature and the vegetation temperature relative to the soil moisture of 0, the soil temperature and the vegetation temperature relative to the soil moisture of 1 and the soil dry edge temperature;
and the earth surface evaporation amount calculation module is used for calculating the earth surface evaporation amount by utilizing a Priestley-Taylor formula according to the soil humidity limiting factor, the vegetation humidity limiting factor and the temperature limiting factor.
7. The remote sensing inversion system of surface evapotranspiration of claim 6, wherein the soil temperature and vegetation temperature calculation module specifically comprises:
a soil temperature calculation submodule with a relative soil moisture of 0 for utilizing a formula according to the meteorological data and the remote sensing data
Figure FDA0002658402680000051
Calculating the soil temperature relative to the soil moisture of 0;
a vegetation temperature calculation submodule with a relative soil moisture of 0 for utilizing a formula according to the meteorological data and the remote sensing data
Figure FDA0002658402680000052
Calculating the vegetation temperature relative to the soil moisture of 0;
a soil temperature calculation submodule with a relative soil moisture of 1 for utilizing a formula according to the meteorological data and the remote sensing data
Figure FDA0002658402680000053
Calculating the soil temperature relative to the soil moisture of 1;
a vegetation temperature calculation submodule with relative soil moisture of 1 for utilizing a formula according to the meteorological data and the remote sensing data
Figure FDA0002658402680000054
Calculating the vegetation temperature relative to the soil moisture of 1;
wherein, Ts0And Tv0Respectively representing the soil temperature and vegetation temperature, T, relative to the soil moisture of 0s1And Tv1Respectively representing the soil temperature and the vegetation temperature relative to the soil moisture of 1; ρ is the air density; cpIs the specific heat at constant pressure; e.g. of the typesSaturated water vapor pressure; e.g. of the typeaThe actual water vapor pressure of the atmosphere; t isaThe temperature near the surface is obtained; gvwVegetation canopy conductivity sufficient for water supply; gavAnd gasThe aerodynamic conductivity of the vegetation and the upper layer of the soil are respectively; rn,sAnd Rn,vRespectively the soil component net radiation and the vegetation component net radiation; gsIs the soil heat flux.
8. The remote sensing inversion system of surface evapotranspiration according to claim 6, wherein the soil dry edge temperature calculation module specifically comprises:
a vegetation coverage calculation submodule for utilizing a formula according to the earth surface remote sensing data
Figure FDA0002658402680000055
Calculating the vegetation ratio Fv
Wherein EVI is Enhanced Vegetation Index (EVI)minAnd EVImaxRespectively the minimum enhancement type vegetation index corresponding to bare soil and the maximum enhancement type vegetation index corresponding to full vegetation coverage;
a soil dry edge temperature operator module for utilizing a formula according to the soil temperature of 0 relative to the soil moisture and the vegetation temperature and vegetation proportion of 1 relative to the soil moisture
Figure FDA0002658402680000061
Calculating and distinguishing soil dry edge temperature T of sufficient water state and no water state on soil surface layer of root zone*
Wherein, Ts0And Tv1The soil temperature is 0 relative to the soil moisture and the vegetation temperature is 1 relative to the soil moisture.
9. The remote sensing inversion system of surface evapotranspiration according to claim 6, wherein the limiting factor calculation module specifically comprises:
the judgment submodule is used for judging whether the surface temperature of the remote sensing image pixel is greater than the dry edge temperature of the soil or not to obtain a judgment result;
a first limiting factor calculating submodule for utilizing the formula f if the judgment result shows noVMCalculating vegetation humidity limiting factor f 1VM(ii) a Using formulas
Figure FDA0002658402680000062
Calculating the soil humidity limiting factor fSM
A second limiting factor calculating submodule for using the formula f if the judgment result shows yesSMCalculating a vegetation humidity limiting factor fVM(ii) a Using formulas
Figure FDA0002658402680000063
Calculating vegetation humidity limiting factor fVM
Wherein, Ts0And Ts1The soil temperature is 0 relative to the soil moisture and the soil temperature is 1 relative to the soil moisture; t issoilWhich is indicative of the temperature of the soil components,
Figure FDA0002658402680000064
TRfor remote-sensing the surface temperature of image elements, FvIs the vegetation proportion; t isv0And Tv1Respectively vegetation temperature of 0 relative to the soil moisture and vegetation temperature of 1 relative to the soil moisture, TvegIs the temperature of the components of the vegetation,
Figure FDA0002658402680000065
a third limiting factor calculation submodule for utilizing a formula
Figure FDA0002658402680000066
Calculating a vegetation temperature limiting factor fT
Wherein, TmaxThe maximum atmospheric temperature in one day; t isoptThe temperature is the optimum temperature for vegetation growth.
10. The remote sensing inversion system of surface evapotranspiration according to claim 6, wherein the surface evapotranspiration calculation module specifically comprises:
a vegetation transpiration meter operator module for utilizing a formula according to the vegetation temperature limiting factor and the vegetation humidity limiting factor
Figure FDA0002658402680000071
Calculating the transpiration amount of the vegetation;
a soil evaporation meter operator module for utilizing a formula according to the soil humidity limiting factor
Figure FDA0002658402680000072
Calculating the evaporation capacity of the soil;
a cut-off evaporation amount calculation operator module for utilizing
Figure FDA0002658402680000073
Calculating the interception evaporation capacity;
the earth surface evapotranspiration meter operator module is used for utilizing a formula ET-ET according to the vegetation evapotranspiration, the soil evaporation and the interception evaporationv+ETs+ETiCalculating the evapotranspiration amount of the earth surface;
wherein, ETv、ETsAnd ETiThe vegetation transpiration amount, the soil evaporation amount and the closure evaporation amount are respectively, and ET is the earth surface evaporation amount; f. ofwetIs the proportion of the wet leaves of the vegetation,
Figure FDA0002658402680000074
RH is the relative humidity of the atmosphere; f. ofgThe ratio of vegetation to green leaves;
Figure FDA0002658402680000075
fAPARphotosynthetically active radiation absorbed for green plants; f. ofIPARPhotosynthetically active radiation intercepted for green vegetation; f. ofTIs a vegetation temperature limiting factor; f. ofSMIs a soil humidity limiting factor; f. ofVMIs vegetation humidity limiting factor; alpha is the coefficient of a Priestley-Taylor formula, and alpha is 1.26; e.g. of the typesSaturated water vapor pressure; t isaThe temperature near the surface is obtained; rn,vFor canopy net radiation, Rn,sNet radiation for soil; gsIs the soil heat flux.
CN202010895830.0A 2020-08-31 2020-08-31 Remote sensing inversion method and system for surface evapotranspiration Active CN112014323B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010895830.0A CN112014323B (en) 2020-08-31 2020-08-31 Remote sensing inversion method and system for surface evapotranspiration

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010895830.0A CN112014323B (en) 2020-08-31 2020-08-31 Remote sensing inversion method and system for surface evapotranspiration

Publications (2)

Publication Number Publication Date
CN112014323A true CN112014323A (en) 2020-12-01
CN112014323B CN112014323B (en) 2021-07-20

Family

ID=73503957

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010895830.0A Active CN112014323B (en) 2020-08-31 2020-08-31 Remote sensing inversion method and system for surface evapotranspiration

Country Status (1)

Country Link
CN (1) CN112014323B (en)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112733445A (en) * 2021-01-05 2021-04-30 四川大学 Large-area-scale soil moisture inversion method based on evapotranspiration vegetation index spatial characteristics
CN113075012A (en) * 2021-03-29 2021-07-06 南京晓庄学院 Moisture evaporation device for remotely sensing and monitoring near-surface soil humidity

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106771073A (en) * 2016-12-28 2017-05-31 中国科学院地理科学与资源研究所 A kind of method that soil and vegetation evapotranspiration are estimated based on end member information model
CN109142674A (en) * 2018-08-02 2019-01-04 中国科学院地理科学与资源研究所 The remote sensing inversion method that root zone soil moisture opposite with surface layer is estimated simultaneously
CN110059362A (en) * 2019-03-22 2019-07-26 兰州大学 Consider the construction method of the regional scale double source evapotranspiration model of vegetation dynamic changes
WO2019145895A1 (en) * 2018-01-24 2019-08-01 The State Of Israel, Ministry Of Agriculture & Rural Development Agricultural Research Organization Method and system for estimating crop coefficient and evapotranspiration of crops based on remote sensing
CN110321784A (en) * 2019-05-08 2019-10-11 中国科学院地理科学与资源研究所 Method, apparatus, electronic equipment and the computer media of soil moisture estimation
CN110334381A (en) * 2019-05-08 2019-10-15 中国科学院地理科学与资源研究所 The method, apparatus and computer media of evapotranspiration estimation

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106771073A (en) * 2016-12-28 2017-05-31 中国科学院地理科学与资源研究所 A kind of method that soil and vegetation evapotranspiration are estimated based on end member information model
WO2019145895A1 (en) * 2018-01-24 2019-08-01 The State Of Israel, Ministry Of Agriculture & Rural Development Agricultural Research Organization Method and system for estimating crop coefficient and evapotranspiration of crops based on remote sensing
CN109142674A (en) * 2018-08-02 2019-01-04 中国科学院地理科学与资源研究所 The remote sensing inversion method that root zone soil moisture opposite with surface layer is estimated simultaneously
CN110059362A (en) * 2019-03-22 2019-07-26 兰州大学 Consider the construction method of the regional scale double source evapotranspiration model of vegetation dynamic changes
CN110321784A (en) * 2019-05-08 2019-10-11 中国科学院地理科学与资源研究所 Method, apparatus, electronic equipment and the computer media of soil moisture estimation
CN110334381A (en) * 2019-05-08 2019-10-15 中国科学院地理科学与资源研究所 The method, apparatus and computer media of evapotranspiration estimation

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
王桐等: ""遥感反演蒸散发的日尺度扩展方法研究进展"", 《遥感学报》 *

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112733445A (en) * 2021-01-05 2021-04-30 四川大学 Large-area-scale soil moisture inversion method based on evapotranspiration vegetation index spatial characteristics
CN113075012A (en) * 2021-03-29 2021-07-06 南京晓庄学院 Moisture evaporation device for remotely sensing and monitoring near-surface soil humidity
CN113075012B (en) * 2021-03-29 2021-11-30 南京晓庄学院 Moisture evaporation device for remotely sensing and monitoring near-surface soil humidity

Also Published As

Publication number Publication date
CN112014323B (en) 2021-07-20

Similar Documents

Publication Publication Date Title
Ca et al. Reductions in air conditioning energy caused by a nearby park
Yang et al. Study on the local climatic effects of large photovoltaic solar farms in desert areas
Molion Amazonia rainfall and its variability
CN112014323B (en) Remote sensing inversion method and system for surface evapotranspiration
Nicholson et al. Micrometeorological conditions and surface mass and energy fluxes on Lewis Glacier, Mt Kenya, in relation to other tropical glaciers
Prueger et al. Crop Water Stress Index of an irrigated vineyard in the Central Valley of California
CN111814317B (en) Surface energy balance component estimation method and system based on remote sensing
Ali Khan Niazi et al. Comparative analysis of photovoltaic configurations for agrivoltaic systems in Europe
Linacre Net radiation to various surfaces
Juvik et al. Forest climatology near the tree line in Hawai'i
Iziomon et al. Empirical models for estimating net radiative flux: A case study for three mid-latitude sites with orographic variability
Stanhill et al. The climate of the cotton crop: Physical characteristics and microclimate relationships
Graser et al. Within-canopy temperature patterns of sorghum at two row spacings
Dickinson Evapotranspiration in global climate models
LeDrew et al. A comparison of the radiation and energy balance during the growing season for an arctic and alpine tundra
Hojdová et al. Microclimate of a peat bog and of the forest in different states of damage in the Šumava National Park
Nullet Water balance of pacific atolls 1
Decker Precision of Estimates of Evapotranspiration in Missouri Climate 1
Lockwood The influence of temperature variations on interception loss and water storage in vegetation canopies
Chaowiwat et al. Prediction of future agriculture water demand in Thailand using multi bias corrected climate models
Boyce An historical comparison between the Met Office Surface Exchange Scheme-Probability Distributed Model (MOSES-PDM) and the Met Office Rainfall and Evaporation Calculation System (MORECS)
Skidmore et al. Potential evaporation as influenced by barrier-induced microclimate
Ren et al. Regional yield prediction for winter wheat based on crop biomass estimation using multi-source data
Wijeratne et al. To Assess the Impact of Urbanization-Associated Land Use Changes on Actual Evapotranspiration and Water Balance in the Kelani River Basin, Sri Lanka
CN115796659B (en) Method for estimating urban water storage capacity based on regional scale evaporation decay

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