CN111947707A - Arid and semi-arid region ground surface water circulation full-component monitoring and identifying method - Google Patents
Arid and semi-arid region ground surface water circulation full-component monitoring and identifying method Download PDFInfo
- Publication number
- CN111947707A CN111947707A CN202010630849.2A CN202010630849A CN111947707A CN 111947707 A CN111947707 A CN 111947707A CN 202010630849 A CN202010630849 A CN 202010630849A CN 111947707 A CN111947707 A CN 111947707A
- Authority
- CN
- China
- Prior art keywords
- data
- precipitation
- soil
- moisture
- water
- 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.)
- Pending
Links
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01D—MEASURING NOT SPECIALLY ADAPTED FOR A SPECIFIC VARIABLE; ARRANGEMENTS FOR MEASURING TWO OR MORE VARIABLES NOT COVERED IN A SINGLE OTHER SUBCLASS; TARIFF METERING APPARATUS; MEASURING OR TESTING NOT OTHERWISE PROVIDED FOR
- G01D21/00—Measuring or testing not otherwise provided for
- G01D21/02—Measuring two or more variables by means not covered by a single other subclass
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Investigating Or Analyzing Materials Using Thermal Means (AREA)
Abstract
The invention discloses a method for monitoring and identifying surface water circulation full components in arid and semiarid regions, which relates to the technical field of physics of an atmospheric boundary layer, and comprises the following steps of: acquiring meteorological data including precipitation, multilayer air humidity, multilayer air temperature, multilayer wind speed and the like; and carrying out pretreatment; acquiring soil data including data such as soil humidity, soil temperature, soil heat flux and the like and preprocessing the soil data; acquiring actual evapotranspiration data; the processed data is input to a non-rainfall moisture identification system which will make the following decisions and calculations. The method overcomes the defects that only precipitation and potential evaporation are utilized for rough estimation at present, and the accuracy is greatly improved compared with the past.
Description
Technical Field
The invention relates to the technical field of atmospheric boundary layer physics, in particular to a method for monitoring and identifying the total component of surface water circulation in arid and semi-arid regions.
Background
Global climate change changes global precipitation and distribution patterns, and increases evaporation, which leads to frequent precipitation extreme events, and increases frequency and intensity of disasters such as flooding, drought and the like. Meanwhile, the contradiction between the instability of water resources and supply and demand is also intensified. Therefore, a better understanding of land water circulation is crucial for human sustainability. However, the estimation of land water balance at present has great defects, the measurement of rainfall and evapotranspiration in arid and semi-arid regions is not accurate enough, and a system capable of automatically identifying each component of land water is unavailable. Thus, the land water circulation is not known.
At present, the contribution of evapotranspiration, precipitation and surface runoff or irrigation is mainly considered for the balance of land water balance. For the rain-fed farmland in arid and semi-arid regions, no surface runoff and irrigation factors exist. It is therefore believed that there is a balance between the amount of precipitation (P) and the amount of transpiration (ET) in a region over a period of time. Namely: p = ET. However, many observations have shown that many times the two are not equal. This is because the contribution of non-reducing water content (NRW) is ignored. After considering the influence of non-precipitation water, the two reach the balance that P + NRW = ET. As shown in fig. 1 and 2.
Disclosure of Invention
The invention aims to provide a reasonably designed method for monitoring and identifying the total amount of surface water circulation in arid and semi-arid regions aiming at the defects and shortcomings of the prior art, overcomes the defect that only precipitation and potential evaporation are utilized for rough estimation at present, and greatly improves the accuracy compared with the prior art.
In order to achieve the purpose, the invention adopts the following technical scheme: it comprises the following steps:
1. automatic observation:
1.1, acquiring meteorological data including precipitation, multilayer air humidity, multilayer air temperature, multilayer wind speed and the like; and carrying out pretreatment;
1.2, acquiring soil data including data such as soil humidity, soil temperature, soil heat flux and the like and preprocessing the soil data;
1.3, acquiring actual Evapotranspiration (ET) data;
2. inputting the data processed in the step 1 into a non-precipitation moisture identification system, and judging and calculating by the data identification system as follows:
2.1, firstly, comparing the current data of the lysimeter with the previous data: if less than 0, it can be concluded that transpiration has occurred during the time period and no non-reducing moisture has been produced; if greater than 0, the value may be non-reducing water;
2.2, judging whether non-precipitation water occurs according to precipitation data; on the basis of the step 2.1, precipitation data are introduced: if precipitation occurs within the observation period of step 1, it is assumed that the added value of the lysimeter data is precipitation, not non-precipitation moisture; if there is no precipitation, then the added value is assumed to be probably non-precipitating moisture;
2.3, distinguishing non-precipitation water and dry settlement according to the dust data; on the basis of step 2.2, the dust data observed contemporaneously are imported: if sand dust exists in the observation period, the reason for increasing the data of the lysimeter is considered to be the sand dust; if no sand dust exists, determining the increment value of the lysimeter as non-precipitation moisture;
2.4, secondly, in order to further determine the composition of the non-precipitating moisture, air relative humidity data are introduced on the basis of step 2.3: if the relative humidity of the air is equal to 100%, it is considered fog; otherwise, the soil is regarded as dew or soil adsorption water;
2.5, finally, distinguishing dew from soil adsorbed water with surface temperature: comparing the surface temperature with the dew point temperature, and if the surface temperature is lower than the dew point temperature, determining the surface temperature is the dew point; otherwise, the soil absorbs water; until that, each component of non-precipitating moisture has been determined;
2.6, acquiring the precipitation and the evapotranspiration by using the observation data and the data identification system;
2.7, precipitation, evapotranspiration and non-precipitation water form complete land moisture balance.
After the method is adopted, the invention has the beneficial effects that: the invention provides a method for monitoring and identifying the total amount of circulating earth surface moisture in arid and semi-arid regions, which overcomes the defect that only precipitation and potential evaporation are used for rough estimation at present and greatly improves the accuracy compared with the prior art.
Description of the drawings:
FIG. 1 is a schematic diagram of the land water balance in general in the background art.
FIG. 2 is a diagram of the land moisture balance in the background art.
FIG. 3 is a schematic illustration of land water balance in the present invention.
FIG. 4 is a schematic illustration of the land water balance of the present invention.
FIG. 5 is a flowchart of the operation of an embodiment.
FIG. 6 is a year distribution graph of the non-reducing water content and the respective components of the example.
FIG. 7 is a graph of non-rainfall water content in each year for the example.
Description of reference numerals:
NRW is non-dewatering water; dew: dew formation; WVA: adsorbing water; fog: and (4) atomizing.
The specific implementation mode is as follows:
the technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the drawings in the embodiments of the present invention, and it is obvious that the described embodiments are only a part of the embodiments of the present invention, and not all of the embodiments. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present invention.
As shown in fig. 3 to 5, the following technical solutions are adopted in the present embodiment: it comprises the following steps:
1. automatic observation:
1.1, acquiring meteorological data including precipitation, multilayer air humidity, multilayer air temperature, multilayer wind speed and the like, and preprocessing;
1.2, acquiring soil data including data such as soil humidity, soil temperature, soil heat flux and the like and preprocessing the soil data;
1.3, acquiring actual Evapotranspiration (ET) data;
2. inputting the data processed in the step 1 into a non-precipitation moisture identification system, and judging and calculating by the data identification system as follows:
2.1, firstly, comparing the current data of the lysimeter with the previous data: if less than 0, it can be concluded that transpiration has occurred during the time period and no non-reducing moisture has been produced; if greater than 0, the value may be non-reducing water;
2.2, judging whether non-precipitation water occurs according to precipitation data; on the basis of the step 2.1, precipitation data are introduced: if precipitation occurs within the observation period of step 1, it is assumed that the added value of the lysimeter data is precipitation, not non-precipitation moisture; if there is no precipitation, then the added value is assumed to be probably non-precipitating moisture; the change in the lysimeter data can be expressed as:in the formula (1)
Wherein Δ w is the change value of the evapo-meter [ kg ]](Δ w =0 or, when there is no NRW<0) I is the sequence number of the instantaneous value of the fixed interval; w is aiAnd wi-1Is the weight value [ kg ] of i-hour and i-1 hour measured by a lysimeter];
2.3, distinguishing non-precipitation water and dry settlement according to the dust data; on the basis of step 2.2, the dust data observed contemporaneously are imported: if sand dust exists in the observation period, the reason for increasing the data of the lysimeter is considered to be the sand dust; if no sand dust exists, determining the increment value of the lysimeter as non-precipitation moisture;
2.4, secondly, in order to further determine the composition of the non-precipitating moisture, air relative humidity data are introduced on the basis of step 2.3: if the relative humidity of the air is equal to 100%, it is considered fog; otherwise, the soil is regarded as dew or soil adsorption water;
2.5, finally, distinguishing dew from soil adsorbed water with surface temperature: comparing the surface temperature with the dew point temperature, and if the surface temperature is lower than the dew point temperature, determining the surface temperature is the dew point; otherwise, the soil absorbs water; until that, each component of non-precipitating moisture has been determined; the dew point temperature calculation formula used the calculated dew point temperature given by Michell instruments corporation. The formula is as follows:
In the formula TdThe dew point temperature is the saturated vapor pressure. The application range of the formula is-45 ℃ to +60 ℃, and T in the formula (2)dThe uncertainty of (a) is +/-0.04 ℃;
2.6, acquiring the precipitation and the evapotranspiration by using the observation data and the data identification system;
2.7, precipitation, evapotranspiration and non-precipitation water form complete land moisture balance.
The working principle of the specific embodiment is as follows: the specific embodiment combines the lysimeter, the air temperature sensor, the air humidity sensor, the sand and dust collector and the soil temperature sensor, and obtains the values of non-rainfall moisture, rainfall and evapotranspiration through the data identification system, thereby realizing the calculation of land moisture balance.
After the method is adopted, the beneficial effects of the specific embodiment are as follows: the specific embodiment provides a method for monitoring and identifying the total amount of circulating earth surface moisture in arid and semi-arid regions, which overcomes the defect that only precipitation and potential evaporation are used for rough estimation at present, and greatly improves the accuracy compared with the prior art.
Example (b):
referring to fig. 6 and 7, the embodiment uses the annual observation data of the fixed west weather station to perform annual calculation by using the identification method of the present invention.
Although the present invention has been described in detail with reference to the foregoing embodiments, it will be apparent to those skilled in the art that various changes in the embodiments and/or modifications of the invention can be made, and equivalents and modifications of some features of the invention can be made without departing from the spirit and scope of the invention.
Claims (1)
1. The method for monitoring and identifying the surface water circulation total component of the arid and semi-arid region is characterized by comprising the following steps of: it comprises the following steps:
(1) and automatic observation:
(1.1) acquiring meteorological data including precipitation, multilayer air humidity, multilayer air temperature, multilayer wind speed and the like; and carrying out pretreatment;
(1.2) acquiring soil data including data such as soil humidity, soil temperature, soil heat flux and the like and preprocessing the soil data;
(1.3) acquiring actual evapotranspiration data;
(2) inputting the data processed in the step (1) into a non-precipitation moisture identification system, and judging and calculating by the data identification system as follows:
(2.1), first, comparing the current data of the lysimeter with the previous data: if less than 0, it can be concluded that transpiration has occurred during the time period and no non-reducing moisture has been produced; if greater than 0, the value may be non-reducing water;
(2.2) secondly, judging whether non-precipitation water occurs or not according to precipitation data; on the basis of the step (2.1), precipitation data are introduced: if precipitation occurs within the observation period of step 1, it is assumed that the added value of the lysimeter data is precipitation, not non-precipitation moisture; if there is no precipitation, then the added value is assumed to be probably non-precipitating moisture;
(2.3) distinguishing non-precipitation water from dry settlement according to the dust data; on the basis of the step (2.2), introducing the dust data observed in the same period into: if sand dust exists in the observation period, the reason for increasing the data of the lysimeter is considered to be the sand dust; if no sand dust exists, determining the increment value of the lysimeter as non-precipitation moisture;
(2.4) secondly, in order to further determine the composition of the non-precipitating moisture, air relative humidity data are introduced on the basis of step (2.3): if the relative humidity of the air is equal to 100%, it is considered fog; otherwise, the soil is regarded as dew or soil adsorption water;
(2.5), finally, distinguishing dew from soil adsorbed water using surface temperature: comparing the surface temperature with the dew point temperature, and if the surface temperature is lower than the dew point temperature, determining the surface temperature is the dew point; otherwise, the soil absorbs water; until that, each component of non-precipitating moisture has been determined;
(2.6) obtaining the precipitation and the evapotranspiration by using the observation data and the data identification system;
(2.7), precipitation, evapotranspiration and non-precipitation water constitute a complete land moisture balance.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010630849.2A CN111947707A (en) | 2020-07-03 | 2020-07-03 | Arid and semi-arid region ground surface water circulation full-component monitoring and identifying method |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010630849.2A CN111947707A (en) | 2020-07-03 | 2020-07-03 | Arid and semi-arid region ground surface water circulation full-component monitoring and identifying method |
Publications (1)
Publication Number | Publication Date |
---|---|
CN111947707A true CN111947707A (en) | 2020-11-17 |
Family
ID=73337393
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010630849.2A Pending CN111947707A (en) | 2020-07-03 | 2020-07-03 | Arid and semi-arid region ground surface water circulation full-component monitoring and identifying method |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111947707A (en) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115266596A (en) * | 2022-04-07 | 2022-11-01 | 中国农业大学 | Evaluation method and system for evapotranspiration in arid region |
Citations (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101556219A (en) * | 2009-04-21 | 2009-10-14 | 中国科学院新疆生态与地理研究所 | Dustfall collection method for separating impurities |
CN102253184A (en) * | 2011-06-29 | 2011-11-23 | 南京信息工程大学 | Remote sensing inversion method for land surface evapotranspiration of arid and semi-arid regions |
CN103645295A (en) * | 2013-12-03 | 2014-03-19 | 中国科学院遥感与数字地球研究所 | Multilayer soil moisture simulation method and multilayer soil moisture simulation system |
CN104266926A (en) * | 2014-09-19 | 2015-01-07 | 西北农林科技大学 | Vaporous soil water acquisition measurement system |
CN105260940A (en) * | 2015-10-22 | 2016-01-20 | 南京信息工程大学 | Crop coefficient correction method based on farmland evapotranspiration observation |
CN105644812A (en) * | 2016-02-25 | 2016-06-08 | 上海宇航***工程研究所 | Environmental simulation method and system capable of achieving low dew point in large normal-pressure space |
CN105760814A (en) * | 2016-01-25 | 2016-07-13 | 中国水利水电科学研究院 | Data mining-based drought monitoring method |
CN106845096A (en) * | 2017-01-12 | 2017-06-13 | 安徽省(水利部淮河水利委员会)水利科学研究院(安徽省水利工程质量检测中心站) | A kind of damage caused by a drought comprehensive estimation method based on multi objective |
CN206740568U (en) * | 2017-03-29 | 2017-12-12 | 水利部牧区水利科学研究所 | A kind of note formula miniature evaporator certainly |
CN109839617A (en) * | 2019-03-13 | 2019-06-04 | 中国气象局气象探测中心 | A kind of non-precipitation echo removing method and system of radar network composite product |
CN110334404A (en) * | 2019-06-10 | 2019-10-15 | 淮阴师范学院 | A kind of rapid dry accurate recognition methods of drought of Watershed Scale |
CN110532672A (en) * | 2019-08-27 | 2019-12-03 | 中国农业科学院农业环境与可持续发展研究所 | Basin pollution of area source priority acccess control area recognition methods based on time-frequency coupled simulation |
CN111173075A (en) * | 2020-01-22 | 2020-05-19 | 长安大学 | Device and method for collecting and utilizing dew in arid region |
CN111191673A (en) * | 2019-11-29 | 2020-05-22 | 广州地理研究所 | Ground surface temperature downscaling method and system |
-
2020
- 2020-07-03 CN CN202010630849.2A patent/CN111947707A/en active Pending
Patent Citations (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101556219A (en) * | 2009-04-21 | 2009-10-14 | 中国科学院新疆生态与地理研究所 | Dustfall collection method for separating impurities |
CN102253184A (en) * | 2011-06-29 | 2011-11-23 | 南京信息工程大学 | Remote sensing inversion method for land surface evapotranspiration of arid and semi-arid regions |
CN103645295A (en) * | 2013-12-03 | 2014-03-19 | 中国科学院遥感与数字地球研究所 | Multilayer soil moisture simulation method and multilayer soil moisture simulation system |
CN104266926A (en) * | 2014-09-19 | 2015-01-07 | 西北农林科技大学 | Vaporous soil water acquisition measurement system |
CN105260940A (en) * | 2015-10-22 | 2016-01-20 | 南京信息工程大学 | Crop coefficient correction method based on farmland evapotranspiration observation |
CN105760814A (en) * | 2016-01-25 | 2016-07-13 | 中国水利水电科学研究院 | Data mining-based drought monitoring method |
CN105644812A (en) * | 2016-02-25 | 2016-06-08 | 上海宇航***工程研究所 | Environmental simulation method and system capable of achieving low dew point in large normal-pressure space |
CN106845096A (en) * | 2017-01-12 | 2017-06-13 | 安徽省(水利部淮河水利委员会)水利科学研究院(安徽省水利工程质量检测中心站) | A kind of damage caused by a drought comprehensive estimation method based on multi objective |
CN206740568U (en) * | 2017-03-29 | 2017-12-12 | 水利部牧区水利科学研究所 | A kind of note formula miniature evaporator certainly |
CN109839617A (en) * | 2019-03-13 | 2019-06-04 | 中国气象局气象探测中心 | A kind of non-precipitation echo removing method and system of radar network composite product |
CN110334404A (en) * | 2019-06-10 | 2019-10-15 | 淮阴师范学院 | A kind of rapid dry accurate recognition methods of drought of Watershed Scale |
CN110532672A (en) * | 2019-08-27 | 2019-12-03 | 中国农业科学院农业环境与可持续发展研究所 | Basin pollution of area source priority acccess control area recognition methods based on time-frequency coupled simulation |
CN111191673A (en) * | 2019-11-29 | 2020-05-22 | 广州地理研究所 | Ground surface temperature downscaling method and system |
CN111173075A (en) * | 2020-01-22 | 2020-05-19 | 长安大学 | Device and method for collecting and utilizing dew in arid region |
Non-Patent Citations (6)
Title |
---|
张强等: "关于干旱和半干旱区陆面水分过程的研究", 《干旱气象》 * |
张强等: "关于干旱和半干旱区陆面水分过程的研究", 《干旱气象》, vol. 25, no. 02, 15 June 2007 (2007-06-15) * |
张强等: "半干旱区土壤水汽吸附的影响因素及变化特征", 《中国科学:地球科学》 * |
张强等: "半干旱区土壤水汽吸附的影响因素及变化特征", 《中国科学:地球科学》, vol. 46, no. 11, 20 November 2016 (2016-11-20), pages 1515 - 1527 * |
张强等: "论干旱区非降水性陆面液态水分分量及其与土壤水分的关系", 《干旱区研究》 * |
张强等: "论干旱区非降水性陆面液态水分分量及其与土壤水分的关系", 《干旱区研究》, vol. 27, no. 03, 15 May 2010 (2010-05-15), pages 392 - 400 * |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115266596A (en) * | 2022-04-07 | 2022-11-01 | 中国农业大学 | Evaluation method and system for evapotranspiration in arid region |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Thomas | Spatial and temporal characteristics of potential evapotranspiration trends over China | |
Zhao et al. | Effects of land cover on variations in stable hydrogen and oxygen isotopes in karst groundwater: A comparative study of three karst catchments in Guizhou Province, Southwest China | |
Kato et al. | Effect of canopy interception on spatial variability and isotopic composition of throughfall in Japanese cypress plantations | |
Xuan | Dust emission factors for environment of Northern China | |
Shimizu et al. | Estimation of annual forest evapotranspiration from a coniferous plantation watershed in Japan (2): Comparison of eddy covariance, water budget and sap-flow plus interception loss | |
Shahidian et al. | Parametric calibration of the Hargreaves–Samani equation for use at new locations | |
Yamanaka et al. | Hydrogen and oxygen isotopes in precipitation in the northern part of the North China Plain: climatology and inter‐storm variability | |
Lee | Relative comparison of the local recalibration of the temperature-based evapotranspiration equation for the Korea Peninsula | |
Hudson et al. | Long-term variability in the water balances of the Plynlimon catchments | |
Wang et al. | Dew measurement and estimation of rain-fed jujube (Zizyphus jujube Mill) in a semi-arid loess hilly region of China | |
Meyer et al. | Simulating shrubs and their energy and carbon dioxide fluxes in Canada's Low Arctic with the Canadian Land Surface Scheme Including Biogeochemical Cycles (CLASSIC) | |
Zhang et al. | Influence of climate variability and human activities on stream flow variation in the past 50 years in Taoer River, Northeast China | |
CN111947707A (en) | Arid and semi-arid region ground surface water circulation full-component monitoring and identifying method | |
Li et al. | Aridity trend and response to vegetation restoration in the loess hilly region of northern Shaanxi Province | |
Spieksma et al. | Modelling evaporation from a drained and rewetted peatland | |
Yue et al. | Surface energy-balance closure in a gully region of the Loess Plateau at SACOL on eastern edge of Tibetan Plateau | |
Yuan et al. | A monthly 0.01° terrestrial evapotranspiration product (1982–2018) for the Tibetan Plateau | |
Xi et al. | How changes of groundwater level affect the desert riparian forest ecosystem in the Ejina Oasis, Northwest China | |
Qin et al. | Climatic and biotic influences on isotopic differences among topsoil waters in typical alpine vegetation types | |
Khosravi et al. | Water vapor pressure trends in south and southwest Iran | |
CN117312746B (en) | Rainfall erosion force calculation method based on satellite rainfall data | |
Sen et al. | Use of GIS and remote sensing in identifying recharge zones in an arid catchment: A case study of Roxo river basin, Portugal | |
Aydin et al. | Simulation of soil hydrological components in Chuncheon over 30 years using E-DiGOR model | |
Rončák et al. | Impact of climate change on runoff and development of forest composition in the coming decades in a selected river basin in Slovakia | |
MacBean et al. | Multi-variable, multi-configuration testing of ORCHIDEE land surface model water flux and storage estimates across semi-arid sites in the southwestern US |
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 | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20201117 |
|
RJ01 | Rejection of invention patent application after publication |