CN109598082A - The calculation method of lake evaporation capacity and crucial hydrographic information based on Hydrogen-oxygen Isotope - Google Patents
The calculation method of lake evaporation capacity and crucial hydrographic information based on Hydrogen-oxygen Isotope Download PDFInfo
- Publication number
- CN109598082A CN109598082A CN201811542949.9A CN201811542949A CN109598082A CN 109598082 A CN109598082 A CN 109598082A CN 201811542949 A CN201811542949 A CN 201811542949A CN 109598082 A CN109598082 A CN 109598082A
- Authority
- CN
- China
- Prior art keywords
- lake
- isotope
- item
- formula
- crucial
- 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
Links
- 238000001704 evaporation Methods 0.000 title claims abstract description 65
- 230000008020 evaporation Effects 0.000 title claims abstract description 63
- 238000004364 calculation method Methods 0.000 title claims abstract description 31
- 239000001301 oxygen Substances 0.000 title claims abstract description 27
- 229910052760 oxygen Inorganic materials 0.000 title claims abstract description 27
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims abstract description 48
- 238000001556 precipitation Methods 0.000 claims abstract description 28
- 230000000155 isotopic effect Effects 0.000 claims abstract description 13
- 230000035764 nutrition Effects 0.000 claims abstract description 7
- 235000016709 nutrition Nutrition 0.000 claims abstract description 7
- 238000004088 simulation Methods 0.000 claims abstract description 4
- 230000029142 excretion Effects 0.000 claims description 9
- 238000000926 separation method Methods 0.000 claims description 6
- 238000005070 sampling Methods 0.000 claims description 5
- 238000005194 fractionation Methods 0.000 claims description 3
- 230000010354 integration Effects 0.000 claims description 3
- 238000000034 method Methods 0.000 abstract description 15
- 238000012544 monitoring process Methods 0.000 abstract description 2
- 238000011160 research Methods 0.000 description 5
- 230000000694 effects Effects 0.000 description 3
- 230000015556 catabolic process Effects 0.000 description 2
- 230000008859 change Effects 0.000 description 2
- 238000006731 degradation reaction Methods 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 description 1
- UFHFLCQGNIYNRP-UHFFFAOYSA-N Hydrogen Chemical compound [H][H] UFHFLCQGNIYNRP-UHFFFAOYSA-N 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 239000012141 concentrate Substances 0.000 description 1
- 230000008878 coupling Effects 0.000 description 1
- 238000010168 coupling process Methods 0.000 description 1
- 238000005859 coupling reaction Methods 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 230000002950 deficient Effects 0.000 description 1
- 238000003745 diagnosis Methods 0.000 description 1
- 239000001257 hydrogen Substances 0.000 description 1
- 229910052739 hydrogen Inorganic materials 0.000 description 1
- 230000002452 interceptive effect Effects 0.000 description 1
- 238000011835 investigation Methods 0.000 description 1
- 239000007788 liquid Substances 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000002360 preparation method Methods 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 238000010025 steaming Methods 0.000 description 1
- 239000002352 surface water Substances 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 238000010792 warming Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Other Investigation Or Analysis Of Materials By Electrical Means (AREA)
Abstract
The calculation method of lake evaporation capacity provided by the invention based on Hydrogen-oxygen Isotope and crucial hydrographic information, comprising steps of (1) obtains the mean annual precipitation and many years average evaporation capacity of lake their location, it is exported lake as sub-basin, the corresponding sub-basin in lake is divided using digital elevation, estimates lake and sub-basin area；(2) lake water, the temperature and relative humidity of precipitation Hydrogen-oxygen Isotope abundance and nutrition lake surface are acquired and measured by certain frequency；(3) the water balance equation and isotopic mass conservation equation for proposing lake, establish computation model；(4) simulation calculates the isotope abundance of steam above lake evaporation item and lake；(5) simultaneous isotope abundance information and model equation calculate the evaporation ratio and crucial hydrographic information in lake.The present invention obtains for the hydrology Ecological information of remote Cross Some Region Without Data and provides new method with field monitoring；Its analysis result credibility is higher, is suitable for a plurality of types of lakes.
Description
Technical field
The invention belongs to hydrographic water resource field, concretely relate to construct conservation of mass model based on Hydrogen-oxygen Isotope
To calculate the technical method of lake evaporation capacity and crucial hydrographic information.
Background technique
In recent years, hydrographic information is combined closely with water isotope conservation of mass model in each aquifer System becomes
Research hotspot.In particular, isotopic mass conservation model be theoretically quite suitable for quantitative output lake hydrological information and
Hydraulic connection signal.Using Hydrogen-oxygen Isotope (18O,2H) constructing simple quality model can effectively identify that lake is equilibrated
All multi-parameters of journey, for remote mountain area and headwaters region, hydrology basic data is quite deficient, however the Lake Water in these areas
Literary situation is very crucial, controls the status of water resources environment in region and affects the following hydrological variation situation and therefore seeks
Rationally effective observation method understands that region lake hydrological element seems very urgent.
Think currently, having some researchs, if isotopic mass conservation model can be applied to the research in these lakes
In, the method for proposing a kind of lake isotope sampling of effective certain frequency constructs Water Balance of Lakes model to obtain lake
The method for mooring hydrographic features status, will solve continuously not observing position for a long time, observation method is single, and observed result is not true
It is qualitative big, the problems such as precision is inadequate.In addition, it is even primary representative same to seek a kind of short-term quickly lake isotope monitoring
Position element sampling and the technical method analyzed highly to represent lake hydrological present characteristics are also huge innovation and challenge.
The research of the Isotope Information diagnosis in many specific lakes is quite lacked, current existing lake isotope mould
Type is more dull, and lake input and output are changed, and lake and surrounding water vapour environment reciprocation excessively simplify, and rests on mostly
In terms of qualitative or sxemiquantitative estimation lake evaporation capacity, it is difficult to it is objective comprehensively show lake in hydrologic cycle role and
Quantitative effect.
Summary of the invention
In order to solve the deficiencies in the prior art, lake evaporation capacity provided by the invention based on Hydrogen-oxygen Isotope and
The calculation method of crucial hydrographic information is to utilize the hydrogen in lake based on a kind of lake isotope conservation model being simple and efficient
The water balance equation in O isotopes combination lake carrys out inverting and explains the hydrologic regime in lake, quantitatively calculates lake
Participate in each element effect of hydrologic cycle, the input and output item including lake, evaporation capacity estimation and the water yield in lake etc.
Crucial hydrographic information provides new technical method to study the lake hydrological effect in basin.
To solve the above problems, the present invention specifically uses following technical scheme:
The calculation method of lake evaporation capacity based on Hydrogen-oxygen Isotope and crucial hydrographic information, which is characterized in that including with
Lower step:
Step 1, the mean annual precipitation p and many years average evaporation capacity e for obtaining lake their location, using lake as son
Basin outlet divides the corresponding sub-basin in lake using digital elevation model, estimates area of lake and sub-basin area；
Step 2, by certain frequency acquisition lake water, the precipitation near lake, and lake water is measured, Hydrogen-oxygen Isotope in precipitation
(2H,18O) the temperature T on abundance and nutrition lake surface, relative humidity h；
Step 3, the water balance equation and isotopic mass conservation equation for proposing lake, establish computation model；
Step 4, simulation calculates the isotope abundance of steam and lake evaporation item above lake；
Step 5, the isotope abundance information and computation model of simultaneous step 2, step 3 and step 4, calculates the pass in lake
Key hydrographic information.
The calculation method of lake evaporation capacity above-mentioned based on Hydrogen-oxygen Isotope and crucial hydrographic information, which is characterized in that
Acquire and measure lake water by certain frequency in the step 2, Hydrogen-oxygen Isotope abundance in precipitation, frequency can be it is primary,
It can sample for several times；When sampling for several times, the isotope abundance of lake water uses the maximum of the isotope abundance of lake water sample for several times
Value, the isotope abundance of precipitation use the weighting isotope abundance mean value of precipitation weight.Preferably, in general practical study,
Recommend using the Lake water isotope abundance in 1 year when 9-10 month.
The calculation method of lake evaporation capacity above-mentioned based on Hydrogen-oxygen Isotope and crucial hydrographic information, which is characterized in that
In the step 3,
The water balance equation in lake is as follows:
I=Q+E (1)
In formula, I, Q, E are input item, excretion item and the evaporation item in lake respectively；
I=P+R (2)
In formula, P is the precipitation for dropping to lake surface, and R is the inbound traffics that do not observe；
The isotopic mass conservation equation in lake is as follows:
IδI=Q δQ+EδE (3)
In formula, δI、δQ、δEIt is input item I, the excretion item Q and the corresponding isotope abundance of evaporation item E in lake respectively；
By (1), (2), (3) integration, obtain:
E/I=(δI-δQ)/(δE-δQ) (4)
In formula, E/I is the ratio of lake evaporation item and input item；
Assuming that water body in lake is sufficiently mixed, the δ of the excretion item Q in lakeQUsually use the isotope of Lake Water
Abundance δLInstead of；Normally due to the inbound traffics R not observed is difficult to determine its isotope abundance, it is assumed that δI≈δp, δpIt is the same of precipitation
The plain abundance in position.
The calculation method of lake evaporation capacity above-mentioned based on Hydrogen-oxygen Isotope and crucial hydrographic information, which is characterized in that
In the step 4, the isotope abundance δ of item is evaporated in lakeECalculation method it is as follows:
δE=((δL-ε+)/α+-hδA-εK)/(1-h+10-3εK) (5)
In formula, h is the relative humidity on nutrition lake surface, ε+It is the equilibrium separation value of isotope, α+It is isotope
Equilibrium fractionation factor, ε+=α+- 1, εKIt is isotopic power separation value, δAFor the isotope abundance of steam above lake；
Wherein, α+It is determined by temperature, α+(18O)=exp [- 7.685/10-3+6.7123/(273.15+T)-1666.4/
(273.15+T)2+350410/(273.15+T)3]) (17)
α+(2H)=exp [1158.8 (273.15+T)3/1012)-1620.1×((273.15+T)2/109)+794.84
((273.15+T)/106)-161.04/103+2999200/(273.15+T)3] (18)
Therefore, formula (4) is rewritten are as follows: E/I=(δL-δI)/(m(δ*-δL)) (6)
In formula,
M=(h-10-3×(εK+ε+/α+))/(1-h+10-3εK) (7)
δ*=(h δA+εK+ε+/α+)/(h-10-3×(εK+ε+/α+)) (8)。
The calculation method of lake evaporation capacity above-mentioned based on Hydrogen-oxygen Isotope and crucial hydrographic information, which is characterized in that
In the step 4, the isotope abundance δ of steam above lakeACalculation method it is as follows:
(1) when area of lake is less than 1km2, and many years average evaporation capacity be less than 1000mm when,
δA=(δP-ε+)/(1+10-3ε+) (9)
(2) when area of lake is less than 1km2, and many years average evaporation capacity be greater than or equal to 1000mm when,
δA=(δP-kε+)/(1+10-3·kε+) (10)
Wherein, k=0.5+ (e-1000)/2e；
(3) when area of lake is greater than 1km2, consider that the evaporation in lake itself can have an impact to steam, at this time water above lake
The isotope abundance of vapour is denoted as δ 'A,
δ′A=(1-f) δA+f·δE (11)
Wherein, f=(1-h),
The isotope abundance of such lake evaporation item is denoted as δ 'E, therefore
δ′E=((δL-ε+)/α+-hδ′A-εK)/(1-h+10-3εK) (12)。
The calculation method of lake evaporation capacity above-mentioned based on Hydrogen-oxygen Isotope and crucial hydrographic information, which is characterized in that
The crucial hydrology information in lake includes in step 5: the ratio E/I of item and input item is evaporated in lake, and what evaporation lake was not observed enters
Flow R, the runoff yield WY of lake basins, the runoff coefficient Z in lake, the residence time г in lake；
The lake evaporation item and the ratio E/I of input item are calculated by formula (6)；
The inbound traffics R calculation method that the lake is not observed is as follows:
R=eLA/ (E/I)-pLA (13)
In formula, e and p are the year evaporation capacity average for many years annual precipitation average with many years respectively, and LA is area of lake；
The runoff yield WY calculation method of the lake basins is as follows:
WY=R/WA (14)
In formula, WA is drainage area；
The runoff coefficient Z in the lake calculates as follows:
Z=R/ (pLA) (15)
Residence time г in the lake calculates as follows:
г=(E/IV)/e (16)
In formula, V is the volume in lake.
Beneficial effects of the present invention:
1, the present invention only by one acquisition lake water and precipitation isotope sample and tests its isotope abundance, constructs lake
Water balance and isotopic mass conservation equation, so that it may which the evaporation capacity and a variety of crucial hydrographic informations for obtaining lake do not need
Complexity is carried out to lake continuously to monitor, it is time saving and energy saving, it is that the hydrology Ecological information acquisition of remote Cross Some Region Without Data is supervised with field
Survey provides new method；
2, the analysis result credibility of this method is higher, effective coupling based on water isotope phase-state change and hydrological model
It closes, has physical basis, a variety of crucial lake hydrological information can be obtained, can be adapted for a plurality of types of lakes, such as large-scale/
Small lakes, seasonality/Out of season lake etc. provide the multi-disciplinary interactive development such as hydrometeorology and chemistry new
Thinking and important science refer to.
Detailed description of the invention
Fig. 1 is the techniqueflow chart of the inventive method.
Specific embodiment
The invention will be further described in the following with reference to the drawings and specific embodiments.
As shown in Figure 1, the present embodiment is many thermokarst lake pools the inventive method to be applied to Northeast Tibetan Plateau
Come obtain evaporation capacity and key lake hydrological information, choose 1 thermokarst lake pool and carry out representative calculating.The region be monsoon with it is non-
The band across of monsoon region, Cold and drought Region ecological environment is significant, and the lake origin cause of formation and variation tendency are complicated, the presence on the thermokarst lake pool and
Development accelerates permafrost degradation around it, profound impact surface water resources distribution.Under Dry-warming trends background, expand extensively at present
The thermokarst lake pool is probably in the following significant, rapid drop, or even disappears.The correlative study in China concentrates on Qinghai-Tibet Platean at present
Change and feedback of the permafrost degradation to the thermokarst lake pool, and control the lower thermokarst lake pool in the work of hydrologic cycle process with ever-frozen ground
It is few for Journal of Sex Research.
1. the lake Calculation of Water producing capacity method based on Hydrogen-oxygen Isotope, preparation needed for early period is as follows, i.e. step 1, step 2:
Consult this area Water Year Book data, it is known that the mean annual precipitation p (322.3mm) of lake their location and
Many years average evaporation capacity e (1354.1mm)；
Exported lake as sub-basin, divide the corresponding sub-basin in lake using digital elevation model, estimate lake and
Sub-basin area, area of lake, which simplifies, is less than 1km2, drainage area 6km2；
Measure lake water, precipitation Hydrogen-oxygen Isotope (2H,18O) abundance, the water surface temperature T and relative humidity h of nutrition,
It carries out a field sampling in April, 2014, and utilizes liquid water isotope analysis instrument (model: Picarro L- in laboratory
2130i) measure,18O guarantees precision 0.025 ‰,2H guarantees precision 0.1 ‰, measurement result are as follows: lake water δL:δ18O=-0.88 ‰,
δ2H=-26.87 ‰, precipitation δP:δ18O=-12.20 ‰, δ2H=-86.40 ‰, T=272.9K, h=0.588.
Step 3: proposing the water balance equation and isotopic mass conservation equation in lake, establish computation model；
Wherein, the water balance equation in lake is as follows:
I=Q+E (1)
In formula, I, Q, E are input item, excretion item and the evaporation item in lake respectively；
I=P+R (2)
In formula, P is the precipitation for dropping to lake surface, and R is the inbound traffics that do not observe；
The isotopic mass conservation equation in lake is as follows:
IδI=Q δQ+EδE (3)
In formula, δI、δQ、δEIt is input item I, the excretion item Q and the corresponding isotope abundance of evaporation item E in lake respectively；
By (1), (2), (3) integration, obtain:
E/I=(δI-δQ)/(δE-δQ) (4)
In formula, E/I is the ratio for evaporating item and output item；
Assuming that the δ of item Q is drained in lake in the case where water body in lake is sufficiently mixedQIt is usually rich with the isotope of Lake Water
Spend δLInstead of；Normally due to the inbound traffics R not observed is difficult to determine its isotope abundance, it is assumed that δI≈δp, δpIt is the same position of precipitation
Plain abundance.
Step 4, simulation calculates the isotope abundance of steam above lake evaporation item and lake；
The isotope abundance δ of lake evaporation itemECalculation method it is as follows:
δE=((δL-ε+)/α+-hδA-εK)/(1-h+10-3εK) (5)
In formula, h is the relative humidity on nutrition lake surface, ε+It is the equilibrium separation value of isotope, α+It is isotope
Equilibrium fractionation factor, ε+=α+- 1, εKIt is isotopic power separation value, δAFor the isotope abundance of steam above lake；
Wherein, α+(18O)=exp [- 7.685/10-3+6.7123/(273.15+T)-1666.4/(273.15+T)2+
350410/(273.15+T)3])=1.01184 (17)
α+(2H)=exp [1158.8 (273.15+T)3/1012)-1620.1×((273.15+T)2/109)+794.84
((273.15+T)/106)-161.04/103+2999200/(273.15+T)3]=1.112 (18)
Formula (4) is rewritten are as follows:
E/I=(δL-δI)/(m(δ*-δL)) (6)
In formula,
M=(h-10-3×(εK+ε+/α+))/(1-h+10-3εK) (7)
δ*=(h δA+εK+ε+/α+)/(h-10-3×(εK+ε+/α+)) (8)；
The isotope abundance δ of steam above lakeACalculation method it is as follows:
(1) when area of lake is less than 1km2, and many years average evaporation capacity be less than 1000mm when,
δA=(δP-ε+)/(1+10-3ε+) (9)
(2) when area of lake is less than 1km2, and many years average evaporation capacity be greater than or equal to 1000mm when,
δA=(δP-kε+)/(1+10-3·kε+) (10)
Wherein, k=0.5+ (e-1000)/2e；
(3) when area of lake is greater than 1km2, consider that the evaporation in lake itself can have an impact to steam, at this time water above lake
The isotope abundance of vapour is denoted as δ 'A,
δ′A=(1-f) δA+f·δE (11)
Wherein, f=(1-h),
The isotope abundance of such lake evaporation item is denoted as δ 'E, therefore
δ′E=((δ L- ε+)/α+-hδ′A-εK)/(1-h+10-3εK) (12),
And since area of lake is less than 1km in the present embodiment2, and many years average evaporation capacity is greater than or equal to 1000mm,
Therefore, using formula (10) δA=(δP-kε+)/(1+10-3·kε+) calculate δA, wherein k=0.5+ (e-1000)/2e=0.631.
Step 5, the isotope abundance information and computation model of simultaneous step 2, step 3 and step 4, calculates the steaming in lake
Hair amount and crucial hydrographic information:
Wherein, the crucial hydrology information in lake includes: the ratio E/I of item and input item is evaporated in lake, and evaporation lake is not seen
The inbound traffics R of survey, the runoff yield WY of lake basins, the runoff coefficient Z in lake, the residence time г in lake；
The lake evaporation item and the ratio E/I of input item are calculated by formula (6),；
The inbound traffics R calculation method that the lake is not observed is as follows:
R=eLA/ (E/I)-pLA (13)
Wherein, e and p is the year evaporation capacity average for many years annual precipitation average with many years respectively, and LA is area of lake；
The runoff yield WY calculation method of the lake basins is as follows:
WY=R/WA (14)
In formula, WA is drainage area；
The runoff coefficient Z in the lake calculates as follows:
Z=R/ (pLA) (15)
Residence time г in the lake calculates as follows:
г=(E/IV)/e (16).
It is as shown in Table 1 that for crucial hydrographic information, (item is evaporated in lake and what the ratio E/I of input item, lake were not observed enters
Flow R, the runoff yield WY of lake basins, the runoff coefficient Z in lake, lake residence time г) result:
The crucial hydrographic information calculated result in 1 lake of table
| E/I | R(mm) | WY(mm) | Z | г(d) | |
| 18O calculated result | 1.06 | 1113.87 | 168.86 | 0.346 | 61.8 |
| 2H calculated result | 0.81 | 776.50 | 87.09 | 0.241 | 76.3 |
| Average result | 0.94 | 945.185 | 127.975 | 0.2935 | 69.05 |
As seen from table:
Lake evaporation, close to 1, shows that Water Balance of Lakes is relatively stable, lake receives increment with the ratio fed
More than half be non-precipitation source, the mean residence time of Lake Water reaches 70 days.The Field Geology Investigations in the lake in early days
The runoff yield of report display, the lake basins is 125mm, and lake increment is almost balanced with excretion, with calculated result (E/I
=0.94, WY=127mm) it is very close.
Basic principles and main features and advantage of the invention have been shown and described above.The technical staff of the industry should
Understand, the present invention is not limited to the above embodiments, and the above embodiments and description only describe originals of the invention
Reason, without departing from the spirit and scope of the present invention, various changes and improvements may be made to the invention, these changes and improvements
It all fall within the protetion scope of the claimed invention.The claimed scope of the invention is by appended claims and its equivalent circle
It is fixed.
Claims (8)
1. the calculation method of lake evaporation capacity and crucial hydrographic information based on Hydrogen-oxygen Isotope, which is characterized in that including following
Step:
Step 1, the mean annual precipitation p and many years average evaporation capacity e for obtaining lake their location, using lake as sub-basin
Outlet divides the corresponding sub-basin in lake using digital elevation model, estimates area of lake and sub-basin area；
Step 2, lake water, precipitation are acquired by certain frequency, and measures lake water, Hydrogen-oxygen Isotope abundance and nutrition lake in precipitation
Moor temperature T, the relative humidity h on surface；
Step 3, the water balance equation and isotopic mass conservation equation for proposing lake, establish computation model；
Step 4, simulation calculates the isotope abundance of steam above lake evaporation item and lake；
Step 5, the isotope abundance information and computation model of simultaneous step 2, step 3 and step 4, calculates the evaporation capacity in lake
And crucial hydrographic information.
2. the calculation method of the lake evaporation capacity according to claim 1 based on Hydrogen-oxygen Isotope and crucial hydrographic information,
It is characterized in that, acquiring and measuring lake water by certain frequency in the step 2, Hydrogen-oxygen Isotope abundance in precipitation, frequency is
Once or for several times sample；When sampling for several times, the isotope abundance of lake water uses the isotope abundance of lake water sample for several times most
Big value, the isotope abundance of precipitation use the weighting isotope abundance mean value of precipitation weight.
3. the calculation method of the lake evaporation capacity according to claim 2 based on Hydrogen-oxygen Isotope and crucial hydrographic information,
It is characterized in that, in step 3,
The water balance equation in lake is as follows:
I=Q+E (1)
In formula, I, Q, E are input item, excretion item and the evaporation item in lake respectively；
I=P+R (2)
In formula, P is the precipitation for dropping to lake surface, and R is the inbound traffics that do not observe；
The isotopic mass conservation equation in lake is as follows:
IδI=Q δQ+EδE (3)
In formula, δI、δQ、δEIt is input item I, the excretion item Q and the corresponding isotope abundance of evaporation item E in lake respectively；
By formula (1), formula (2), formula (3) integration, obtain:
E/I=(δI-δQ)/(δE-δQ) (4)
In formula, E/I is the ratio for evaporating item and output item；
In the case where water body in lake is sufficiently mixed, the δ of the excretion item Q in lakeQWith the isotope abundance δ of Lake WaterLInstead of；And
δI≈δp, δpIt is the isotope abundance of precipitation.
4. the calculation method of the lake evaporation capacity according to claim 3 based on Hydrogen-oxygen Isotope and crucial hydrographic information,
It is characterized in that, the isotope abundance δ of item is evaporated in lake in step 4ECalculation method it is as follows:
δE=((δL-ε+)/α+-hδA-εK)/(1-h+10-3εK) (5)
In formula, h is the relative humidity on nutrition lake surface, ε+It is the equilibrium separation value of isotope, α+It is the balance of isotope
Fractionation factor, ε+=α+- 1, εKIt is isotopic power separation value, δAFor the isotope abundance of steam above lake；
Wherein, formula (4) is rewritten as
E/I=(δL-δI)/(m(δ*-δL)) (6)
In formula,
M=(h-10-3×(εK+ε+/α+))/(1-h+10-3εK) (7)
δ*=(h δA+εK+ε+/α+)/(h-10-3×(εK+ε+/α+)) (8)。
5. the calculation method of the lake evaporation capacity according to claim 4 based on Hydrogen-oxygen Isotope and crucial hydrographic information,
It is characterized in that, in step 4, the isotope abundance δ of steam above lakeACalculation method it is as follows:
When area of lake is less than 1km2, and many years average evaporation capacity be less than 1000mm when,
δA=(δP-ε+)/(1+10-3ε+) (9)
When area of lake is less than 1km2, and many years average evaporation capacity be greater than or equal to 1000mm when,
δA=(δP-kε+)/(1+10-3·kε+) (10)
In formula, k=0.5+ (e-1000)/2e；
When area of lake is greater than 1km2, consider that the evaporation in lake itself can have an impact to steam, the isotope of steam is rich above lake
Degree is denoted as δ 'A,
δ′A=(1-f) δA+f·δE (11)
In formula, f=(1-h)；
Then the isotope abundance of such lake evaporation item is denoted as δ 'E, as follows
δ′E=((δL-ε+)/α+-hδ′A-εK)/(1-h+10-3εK) (12)。
6. the calculation method of the lake evaporation capacity according to claim 5 based on Hydrogen-oxygen Isotope and crucial hydrographic information,
It is characterized in that, the crucial hydrology information in lake includes in step 5: the ratio E/I of item and input item is evaporated in lake, evaporates lake
The inbound traffics R not observed, the runoff yield WY of lake basins, the runoff coefficient Z in lake and the residence time г in lake；
The lake evaporation item and the ratio E/I of input item are calculated by formula (6)；
The inbound traffics R calculation method that the lake is not observed is as follows:
R=eLA/ (E/I)-pLA (13)
In formula, e and p are the year evaporation capacity average for many years annual precipitation average with many years respectively, and LA is area of lake；
The runoff yield WY calculation method of the lake basins is as follows:
WY=R/WA (14)
In formula, WA is drainage area；
The runoff coefficient Z in the lake calculates as follows:
Z=R/ (pLA) (15)
Residence time г in the lake calculates as follows:
г=(E/IV)/e (16)
In formula, V is the volume in lake.
7. the calculation method of the lake evaporation capacity according to claim 2 based on Hydrogen-oxygen Isotope and crucial hydrographic information,
It is characterized in that, the Lake water isotope abundance is using the Lake water isotope abundance in 1 year when 9-10 month.
8. the calculation method of the lake evaporation capacity according to claim 4 based on Hydrogen-oxygen Isotope and crucial hydrographic information,
It is characterized in that, the α+It is determined by temperature, α+(18O)=exp [- 7.685/10-3+6.7123/(273.15+T)-1666.4/
(273.15+T)2+350410/(273.15+T)3]) (17)
α+(2H)=exp [1158.8 (273.15+T)3/1012)-1620.1×((273.15+T)2/109)+794.84
((273.15+T)/106)-161.04/103+2999200/(273.15+T)3] (18)。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811542949.9A CN109598082B (en) | 2018-12-17 | 2018-12-17 | Calculation method of lake evaporation capacity and key hydrological information based on hydrogen and oxygen isotopes |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811542949.9A CN109598082B (en) | 2018-12-17 | 2018-12-17 | Calculation method of lake evaporation capacity and key hydrological information based on hydrogen and oxygen isotopes |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN109598082A true CN109598082A (en) | 2019-04-09 |
| CN109598082B CN109598082B (en) | 2020-06-02 |
Family
ID=65963043
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201811542949.9A Expired - Fee Related CN109598082B (en) | 2018-12-17 | 2018-12-17 | Calculation method of lake evaporation capacity and key hydrological information based on hydrogen and oxygen isotopes |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109598082B (en) |
Cited By (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111504277A (en) * | 2020-03-20 | 2020-08-07 | 河海大学 | Lake water supply tracing method utilizing hydrogen and oxygen isotopes |
| CN111539090A (en) * | 2020-03-20 | 2020-08-14 | 河海大学 | Watershed evapotranspiration estimation method combining water body hydrogen and oxygen isotopes and chloride ions |
| CN113256101A (en) * | 2021-05-20 | 2021-08-13 | 东南大学 | Key driving force analysis method for water storage capacity change of urban and rural lakes |
| CN113724802A (en) * | 2021-09-06 | 2021-11-30 | 长安大学 | Calculation method for influence of decoupling evaporation on salt content of lake and reservoir water body |
| CN115144556A (en) * | 2022-07-05 | 2022-10-04 | 重庆交通大学 | Method and device for quantifying rainfall water vapor source based on isotope and leaf area index |
| CN115420811A (en) * | 2022-07-21 | 2022-12-02 | 河海大学 | Method for analyzing fractionation effect of hydrogen and oxygen isotopes in evaporated water vapor |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20130146542A1 (en) * | 2010-08-06 | 2013-06-13 | Brown University | Functionalized chromatographic materials and methods of making and using therefor |
| CN106370802A (en) * | 2016-08-16 | 2017-02-01 | 华南理工大学 | Method for determining water loss of reservoir by using hydrogen and oxygen isotopes |
| CN106885892A (en) * | 2017-02-22 | 2017-06-23 | 北京林业大学 | The method and device split to forest ecosystem evapotranspiration |
-
2018
- 2018-12-17 CN CN201811542949.9A patent/CN109598082B/en not_active Expired - Fee Related
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20130146542A1 (en) * | 2010-08-06 | 2013-06-13 | Brown University | Functionalized chromatographic materials and methods of making and using therefor |
| CN106370802A (en) * | 2016-08-16 | 2017-02-01 | 华南理工大学 | Method for determining water loss of reservoir by using hydrogen and oxygen isotopes |
| CN106885892A (en) * | 2017-02-22 | 2017-06-23 | 北京林业大学 | The method and device split to forest ecosystem evapotranspiration |
Cited By (12)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111504277A (en) * | 2020-03-20 | 2020-08-07 | 河海大学 | Lake water supply tracing method utilizing hydrogen and oxygen isotopes |
| CN111539090A (en) * | 2020-03-20 | 2020-08-14 | 河海大学 | Watershed evapotranspiration estimation method combining water body hydrogen and oxygen isotopes and chloride ions |
| CN111504277B (en) * | 2020-03-20 | 2021-05-04 | 河海大学 | Lake water supply tracing method utilizing hydrogen and oxygen isotopes |
| CN111539090B (en) * | 2020-03-20 | 2023-04-07 | 河海大学 | Watershed evapotranspiration estimation method combining water body hydrogen and oxygen isotopes and chloride ions |
| CN113256101A (en) * | 2021-05-20 | 2021-08-13 | 东南大学 | Key driving force analysis method for water storage capacity change of urban and rural lakes |
| CN113256101B (en) * | 2021-05-20 | 2022-10-28 | 东南大学 | Key driving force analysis method for water storage capacity change of urban and rural lakes |
| CN113724802A (en) * | 2021-09-06 | 2021-11-30 | 长安大学 | Calculation method for influence of decoupling evaporation on salt content of lake and reservoir water body |
| CN113724802B (en) * | 2021-09-06 | 2023-08-11 | 长安大学 | Calculation method for influence of decoupling evaporation on salt content of lake and reservoir water body |
| CN115144556A (en) * | 2022-07-05 | 2022-10-04 | 重庆交通大学 | Method and device for quantifying rainfall water vapor source based on isotope and leaf area index |
| CN115144556B (en) * | 2022-07-05 | 2024-03-26 | 重庆交通大学 | Method and device for quantifying rainfall vapor source based on isotope and leaf area index |
| CN115420811A (en) * | 2022-07-21 | 2022-12-02 | 河海大学 | Method for analyzing fractionation effect of hydrogen and oxygen isotopes in evaporated water vapor |
| CN115420811B (en) * | 2022-07-21 | 2024-05-10 | 河海大学 | Analysis method for oxygen isotope fractionation effect in evaporated water vapor |
Also Published As
| Publication number | Publication date |
|---|---|
| CN109598082B (en) | 2020-06-02 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN109598082A (en) | The calculation method of lake evaporation capacity and crucial hydrographic information based on Hydrogen-oxygen Isotope | |
| Dutkiewicz et al. | Interannual variability of phytoplankton abundances in the North Atlantic | |
| Cui et al. | Application of remote sensing to water environmental processes under a changing climate | |
| Liu et al. | Stable isotopic compositions of precipitation in China | |
| Horita et al. | Isotope effects in the evaporation of water: a status report of the Craig–Gordon model | |
| Rϋhland et al. | Reorganization of algal communities in the Lake of the Woods (Ontario, Canada) in response to turn‐of‐the‐century damming and recent warming | |
| Jasechko et al. | Stable isotope mass balance of the Laurentian Great Lakes | |
| Dong et al. | The role of interocean exchanges on decadal variations of the meridional heat transport in the South Atlantic | |
| Zareian et al. | The effect of climate change on watershed water balance | |
| Shi et al. | Evaluation of the parameter sensitivities of a coupled land surface hydrologic model at a critical zone observatory | |
| Crespin et al. | Re-examination of the temperature-dependent relationship between δ 18 O diatoms and δ 18 O lake water and implications for paleoclimate inferences | |
| Huang et al. | Net community production in the South China Sea Basin estimated from in situ O2 measurements on an Argo profiling float | |
| Seginer et al. | Modeling solar drying rate of wastewater sludge | |
| Xia et al. | Quantifying the effects of climate change and human activities on runoff in the water source area of Beijing, China | |
| Lamane et al. | Soil water erosion assessment in Morocco through modeling and fingerprinting applications: A review | |
| Rice et al. | The influence of watershed characteristics on spatial patterns of trends in annual scale streamflow variability in the continental US | |
| Sun et al. | Quantitative evaluation of the rainfall influence on streamflow in an inland mountainous river basin within Central Asia | |
| Pan et al. | Using data‐driven methods to explore the predictability of surface soil moisture with FLUXNET site data | |
| Rosero-López et al. | Standardization of instantaneous fluoroprobe measurements of benthic algal biomass and composition in streams | |
| Anoh et al. | Modelling water flow in a complex watershed in humid a tropical area using SWAT: a case study of Taabo watershed in Ivory Coast | |
| Zhang et al. | The temporal variations in runoff-generation parameters of the Xinanjiang model due to human activities: A case study in the upper Yangtze River Basin, China | |
| Wang et al. | Comparison of boosted regression trees vs WA-PLS regression on diatom-inferred glacial-interglacial climate reconstruction in Lake Tiancai (southwest China) | |
| Xie et al. | Striking seasonal pattern of primary production in the river-dominated ocean margin of the northern South China Sea (NSCS-RiOMar) revealed by new field and remotely sensed data | |
| Cai et al. | Stable water isotope and surface heat flux simulation using ISOLSM: Evaluation against in-situ measurements | |
| Iraheta et al. | A preliminary isotope‐based evapotranspiration partitioning approach for tropical Costa Rica |
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 | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20200602 |