Next Article in Journal
Effects of Land Cover and Atmospheric Input on Nutrient Budget in Subtropical Mountainous Rivers, Northeastern Taiwan
Previous Article in Journal
Urban Groundwater Processes and Anthropogenic Interactions (Porto Region, NW Portugal)
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Quantitative Analysis of the Sub-Cloud Evaporation of Atmospheric Precipitation and Its Controlling Factors Calculated By D-Excess in an Inland River Basin of China

1
Yulong Snow Mountain Glacier and Environment Observation and Research Station/State Key Laboratory of Cryospheric Sciences, Northwest Institute of Eco-Environment and Resources, Chinese Academy of Sciences, Lanzhou 730000, China
2
College of Geography and Environment Science, Northwest Normal University, Lanzhou 730070, China
*
Authors to whom correspondence should be addressed.
Water 2020, 12(10), 2798; https://doi.org/10.3390/w12102798
Submission received: 15 July 2020 / Revised: 30 September 2020 / Accepted: 3 October 2020 / Published: 9 October 2020
(This article belongs to the Section Hydrology)

Abstract

:
Atmospheric precipitation is an important part of the water circle in an inland basin. Based on the analytical results of 149 precipitation samples and corresponding surface meteorological data collected at four sampling sites (Lenglong, Ningchang, Huajian and Xiying) at different elevations in the Xiying river basin on the north slope of Qilian Mountains from May to September 2017, the sub-cloud evaporation in precipitation and its controlling factors are analyzed by the Stewart model. The results show that sub-cloud evaporation led to d-excess value in precipitation decrease and d-excess variation from cloud-base to near surface (Δd) increase with decreasing altitude. The remaining evaporation fraction of raindrop (f) decreases with decreasing altitude. The difference of underlying surface led to a difference change of f and Δd in the Xiying sampling site. For every 1% increase in raindrop evaporation, d-excess value in precipitation decreased by about 0.99‰. In an environment of high relative humidity and low temperature, the slope of the linear relationship between f and Δd is less than 0.99. In contrast, in the environment of low relative humidity and high temperature, the slope is higher than 0.99. In this study, set constant raindrop diameter may affect the calculation accuracy. The Stewart model could have different parameter requirements in different study areas. This research is helpful to understand water cycle and land–atmosphere interactions in Qilian Mountains.

1. Introduction

In a cold-alpine region, runoff is mainly fed by Cryosphere meltwater (glaciers, snow, permafrost) and precipitation [1]. Almost all glaciers in the world are shrinking, the thicknesses of permafrost active layers are increasing, and snow cover extent and duration are decreasing [2]. In the alpine region of northwest China, which is the main distribution area of Cryosphere (including glacier, permafrost and snow cover), the annual mean surface temperature has risen by 1.8 °C from 1960 to 2007 [3]. In recent decades, the area of glaciers, permafrost and snow cover in Qilian Mountains has decreased significantly [4,5,6]; on the contrary, precipitation is increasing. In the context of global warming, the contribution of precipitation to runoff has continued to increase, while the contribution of Cryosphere meltwater has gradually decreased [7]. Based on the above background, further research on precipitation can help us better understanding the process of regional water cycle. Hydrogen- and oxygen-stable isotopes are important tracers and environmental indicators which have been widely used in the study of water. Based on that, d-excess is defined as a second level parameter of hydrogen- and oxygen-stable isotopes (d-excess = cycle and hydrologic processes [8,9,10]. Hydrogen- and oxygen-stable isotopes are affected by many processes, such as evaporation, transpiration, condensation and so on, with complex variationsδD−8δ18O) [11]. d-excess is influenced by evaporation, vapor mixing and moisture sources [12]. The influence of transpiration and condensation on d-excess is negligible [13,14]. Therefore, d-excess can be used to calculate the evaporation fraction of the water cycle [15,16].
Sub-cloud evaporation is the process of raindrops from cloud-base to near surface under the condition of unsaturated water vapor pressure, accompanied by the enrichment of 18O and D and the reduction of d-excess [11,17]. In arid areas, sub-cloud evaporation is especially obvious [18,19]. Sub-cloud evaporation is directly affected by the meteorological conditions in the study area [17]. Therefore, the study on sub-cloud evaporation during precipitation is the premise and foundation for the study of watershed hydrological process [20,21]. Stewart et al. (1975) used a model to simulate stable isotopes variation in falling raindrops under the influence of sub-cloud evaporation in different environments in laboratory [22]. According to the cloud physics method, Zhang et al. (1998) improved the Stewart model and found that the stable isotope in raindrops in unsaturated atmosphere were enriched with the increase of falling distance, which was more obvious in the case of air drying [23]. Meanwhile, sub-cloud evaporation also decreases the slope of the local meteoric water line [24,25]. Wang et al. (2016) found that the linear relationship between raindrop evaporation and d-excess is also affected by raindrop diameter, temperature, relative humidity and precipitation in the Tianshan area [26]. However, the Stewart model needs many parameters (temperature, dew-point temperature, precipitation amount, relative humidity, pressure, raindrop diameter, cloud-base height and so on) and the calculation process is complex, so many researches directly quoted Froehlich’s conclusion: constant linear relationship of approximately 1‰ of d-excess value change per 1% evaporation of raindrop [16,18,27].
Qilian Mountains not only ensures the ecological security of Hexi Corridor, but also plays an important role in maintaining the ecological balance of Qinghai-Tibet Plateau, preventing desert expansion and maintaining the stability of oasis. There are three inland river systems in this region: Shiyang river, Heihe river and Shule river from the east to the west in the north slope of Qilian Mountains. Previous studies show that the contributions of precipitation to runoff were 77% in Shiyang river basin, 65% in Heihe river basin and 51% in Shule river basin [1]. With the climate warming and cryosphere shrinking, the contribution of precipitation to runoff will continue to increase in Qilian Mountains, and the sub-cloud evaporation in the water cycle would also gradually increase. As a largest tributary of the Shiyang River, the precipitation across the Xiying river basin is not only affected by large-scale circulation, but also by orographic lifting, causing a complex precipitation pattern. In the past, the research mainly focused on the variation of stable isotopes in precipitation, with little systematic research on sub-cloud evaporation in precipitation and its influencing factors. In this paper, we use the Stewart model to quantitatively study the correlation between sub-cloud evaporation and d-excess at different elevations in the Xiying river basin, and discuss the influence of meteorological factors, underlying surface and raindrop diameter. The above studies will provide a scientific basis for utilization of the water resources in the Qilian Mountains and Hexi corridor.

2. Study Area, Date and Method

2.1. Study Area

The Qilian Mountains lie in the intersection belt of Qinghai-Tibet Plateau, Mongolia-Xinjiang Plateau and Loess Plateau (94° E∼104° E and 36° N∼40° N), which are composed of several parallel mountains and valleys (Figure 1). The total area is 18.3 × 104 km2, with an average elevation of 2800 m [28], with some peaks over 4000 m to the cryosphere area. Influenced jointly by westerly, eastern monsoons, plateau monsoons and local moisture recycling, Qilian Mountains have diversified ecosystems such as forests, grasslands, alpine meadows, etc. [29,30]. The Xiying river basin lies in the northern slope of the eastern Qilian Mountains, which is the largest tributary of Shiyang River. It is mainly supplied by glacier-snow meltwater, permafrost meltwater and precipitation. In the Xiying river basin, summer is short and hot and winter is long and cold, with large temperature differences between day and night and strong evaporation. The annual mean temperature is about 5 °C, precipitation amount is about 400 mm/year and the pan evaporation is 850 mm/year.

2.2. Data Sources

Precipitation samples (n = 149) were collected at four sampling sites (Lenglong, Ningchang, Huajian and Xiying) from May to September 2017 (Figure 1, Table 1). The Lenglong sampling site is located at the frozen soil zone in the upper reaches of the Xiying river basin. Ningchang and Huanjian sampling sites are located at the alpine forest zone in the middle reaches of the Xiying river basin. The Xiying sampling site is located at the sub-alpine meadow zone in the lower reaches of the Xiying river basin. Meteorological data are recorded by the automatic weather station there. The sample collector was placed in naturally vegetated land undisturbed, 1.5 m from the surface. Placing a table tennis ball in the collection funnel sealed the collector bottle against evaporation and debris. After each precipitation event, precipitation samples were collected using 50 mL polyethylene bottles. All samples were stored in a refrigerator at −18 °C.
Stable isotopic ratios of precipitation were analyzed at the Laboratory of Stable Isotopes in the College of Geography and Environmental Science at Northwest Normal University by the DLT-100 Liquid Water Isotope Analyzer instrument (developed by the Los Gatos Research company of the United States). Measurement uncertainties were ±0.2‰ for δ18O and ±0.6‰ for δ2H. Normally, all the samples were stored in the freezer to prevent the fractionation of isotopes and they were taken out to thaw at room temperature before analysis. In the test, every test group included three standard samples (standard No.3: δ2H: −96.4 ± 0.5‰, δ18O: −13.10 ± 0.15‰; standard No. 4: δ2H: −51.0 ± 0.5‰, δ18O: −7.69 ± 0.15‰; standard No. 5: δ2H: −9.5 ± 0.5‰, δ18O: −2.80 ± 0.15‰; provided by the LGR company) and six unknown precipitation samples. Every sample was tested for six injections. Data of the first two needles were discarded because of the isotope memory, and values of the last four needles were calculated as final results. Results tested were relative to the Vienna Standard Mean Ocean Water (VSMOW).
δ s a m p l e = ( R s a m p l e R s t a n d a r d ) R s t a n d a r d
Rsample presents the ratio of 2H/1H (18O/16O) in water samples. Rstandard shows the ratio of 2H/1H (18O/16O) in the VSMOW.

2.3. Method

Stewart et al. assumed that the cloud-base vapor reached an isotopic equilibrium state, and the variation of d-excess (Δd) (d-excess values of precipitation at near surface minus d-excess values of precipitation at cloud-base) can be calculated by the following formula [22,26]:
Δ d   =   ( 1     γ 2 α 2   ) ( f β 2     1 )     8 ( 1     γ 18 α 18 ) ( f β 18     1 )
where f is the evaporation remaining fraction. 2γ, 18γ, 2β and 18β are defined by Stewart [22]. 2α and 18α are isotope equilibrium fractionation factor and can be calculated by the following formula [31]:
α 2   =   e x p ( 24.844   ×   10 3 T L C L 2     76.248 T L C L   +   52.612   ×   10 3 )
α 18   =   e x p ( 1.137   ×   10 3 T L C L 2     0.4156 T L C L     2.0667   ×   10 3 )
where TLCL is air temperature (K) at lifting condensation level (LCL). According to Barnes [32], the calculation method is
T L C L = T d ( 0.001296 T d + 0.1963 ) ( T T d )
where Td and T are the dew-point temperature and surface air temperature (°C).
According to Kinzer and Gunn [33], f can be calculated by the following formula:
f = m e n d m e n d + m e v    
where mend is the mass of the raindrop fall to surface and mev is the mass of raindrop evaporated. They can be calculated by the following formula:
m e v = E t
  m e n d = 4 3 π r e n d 3 ρ
where E is evaporation intensity, t is falling time of raindrops. rend is raindrop radius at landing, ρ is the density of water.
According to Kinzer and Gunn [33], evaporation intensity can be calculated by the following formula:
E = Q 1 ( T , D ) Q 2 ( T , h )
where Q1 is controlled by temperature (T) and raindrop diameter (D), Q2 is controlled by temperature (T) and relative humidity (h), and the specific calculation method can be found by Kinzer and Gunn [33].
Based on the calculation method proposed by Best [34], Wang et al. proposed a modified empirical formula to calculate median diameter of the raindrops in a semi-arid cold region [26]:
D = 0.69 n A I P
where I is precipitation intensity (mm·h−1), the parameter of n, A, and P are defined by Wang et al. [26].
Raindrops quickly reach the state of equal velocity motion in the process of falling, and the falling time is:
t = H c b v e n d
where Hcb is cloud-base height, vend is the end velocity of the raindrop which can be calculated as [31].
In the previous researches, the height of cloud-base is set to a constant [19,35]. However, in many cases, the height of cloud-base is often less than 1500 m, affected by topography and meteorological factors in alpine areas. In this study, the height of cloud-base (Hcb) can be calculated by the following formula [36]:
H c b = 18400 ( 1 + T m e a n 273 ) lg   P 0 P L C L
where Tmean is average temperature (°C) between the temperature at lifting condensation level and surface temperature, P0 is surface pressure (hPa), PLCL is pressure (hPa) at lifting condensation level, and PLCL can be calculated by the following formula [32]:
P L C L = P ( T L C L T ) 3.5
where P is pressure at the sampling site (hPa).

3. Result and Discussion

3.1. Variations of d-excess

Table 1 shows d-excess varied from cloud-base to near surface at each sampling site in the Xiying river basin in summer from May to September. With decreasing altitude, d-excess values at near surface decreased from 19.72‰ (Lenglong) to 10.23‰ (Huajian), and then increased to 13.9‰ (Xiying). Δd increased from 20.29‰ to 35.15‰, and then decrease to 29.45‰. In the upper and middle reaches of Xiying river basin (Lenglong, Ningchang and Huajian), d-excess values at near surface decreased with decreasing altitude, and Δd increased with decreasing altitude. In the lower reaches, the change of d-excess at near surface and Δd show opposite trends compared with those of middle and upper reaches, which may be caused by different underlying surfaces. d-excess values at near surface were much lower than that in the cloud-base, indicating that sub-cloud evaporation plays a significant role in precipitation process in semi-arid cold region.

3.2. Variations of Δd

In the middle and upper reaches (Lenglon, Ningchang and Huajian), the evaporation remaining fraction (f) decreased with decreasing altitude, but suddenly increased in the lower reaches (Xiying) (Figure 2a,b). The f fluctuated between 9.71% and 95.42% in May to September 2017. As shown in Table 2, f in May was relatively lower than that in other months. This is mainly because relative humidity is lower and evaporation intensity is higher than that in other months. Between June and August, relative humidity change was not obvious at each sampling site. Sub-cloud evaporation is mainly affected by air temperature, f decreases with the increase of air temperature and evaporation intensity. In September, the air temperature and evaporation intensity weakened, which led to the increase of f. Δd depends on altitude and time. As a whole, Δd increases gradually with decreasing altitude (Figure 2c), and Δd was higher in May than that in June to September (Figure 2d).

3.3. The Correlation between f and Δd

It is clear that there is a significant correlation between f and Δd, with a slope of about 0.99 (Figure 3). When f is greater than 95%, the correlation between f and Δd is highly significant with a slope of 0.81, indicating that f increase by 1% would lead to Δd decrease about 0.81‰. It shows that when sub-cloud evaporation is weak and f is high (higher relative humidity, lower temperature and larger raindrop radius), there is a significant correlation between f and Δd. The relationship gradually becomes weaker as the f decreases (relative humidity decreases, temperature rises, and raindrop radius decreases). When f is below 40%, there is a stronger scattering around the regression line. Analysis of the relationship between f and Δd at the middle and upper reaches (Lenglong, Ningchang, Huajian) reflects an increasing trend of slope and decreasing trend of correlation as altitude decreases (Table 3). Previous researches [26,35] suggest that there is a significant correlation between f and Δd with a slope about 1 in a context of high f. In a semi-arid cold region, f may be much lower than 95%, with most small raindrops disappearing as they fall to the near surface [19]. Researchers usually assume that there are a 1‰·%−1 correlation between f and Δd [16,37]. However, the results of this paper suggest that this assumption also needs to consider the influence of climate, altitude and other factors.

3.4. Influencing Factors of Sub-cloud Evaporation and Δd

3.4.1. Influence of Artificial Reservoir

The f gradually decreased with decreasing altitude and Δd gradually increased with decreasing altitude, the same as the conclusion in the middle and upper reaches of this study [36]. However, the f suddenly increased and Δd suddenly decreased at the downstream, reflecting the influence of other factors. Through field investigation, we found that there exists an artificial reservoir (Xiying reservoir) near the Xiying sampling site with a distance of 1.4 km between them. As the largest artificial water body in Xiying river basin, with the total storage capacity of 23.5 million m3 [38], the difference of underlying surface may lead to different contribution of local moisture recycling. Because large water body (lakes, reservoirs, etc.) will increase local moisture recycling [26], the monthly average evaporation of Xiying reservoir in summer half year (May to September) can reach 133 mm ± 18.6 mm [38]. In addition, strong evaporation of the Xiying reservoir will lead to decreased temperature, increased relative humidity and decreased vapor pressure deficit in the surrounding area. These changes will lead to a weakened sub-cloud evaporation in precipitation, increased f, and decreased Δd. The temperature, relative humidity and cloud-base height of precipitation at each sampling site were analyzed to verify the above analysis further. As shown in Figure 4, relative humidity, air temperature and cloud-base height show obvious trends of decreasing relative humidity, increasing temperature and cloud-base height corresponding to the decreasing altitude in the middle and upper reaches, but these trends are different downstream. Cloud-base height in precipitation increases with decreasing altitude as mentioned in previous studies [26], further confirming that the sudden increase of f downstream is due to the difference of underlying surface.

3.4.2. Influence of Meteorological Factors

When the temperature is between 0 °C to 8 °C, there is an obvious linear correlation between f and Δd, and the slope is slightly lower than 1. When the temperature is between 8 °C to 15 °C, the linear correlation between them weakens, and the slope is slightly higher than 1. When the temperature is above 15 °C, the linear correlation between them is further weakened with the slope up to 1.24 (Figure 5a–c). When the relative humidity is below 75%, the linear correlation between them is low, and the slope is as high as 1.64. When the relative humidity is between 75% to 90%, the linear correlation is gradually increasing, and the slope sharply decreased to 0.69. When the relative humidity is above 90%, the slope decreased to 0.61 (there are only few points in Figure 5c, the conclusions might be subject to high error) (Figure 5d–f). When the precipitation is below 5 mm, the linear correlation is relatively high and the slope is lower than 1. When the precipitation is between 5 mm~10 mm, the linear correlation gradually decreased, and the slope is slightly higher than 1. When the precipitation is above 10 mm, the slope increased to 1.13 (Figure 5g,h,i). Therefore, it can be concluded that f is high and Δd is close to 0 in the environment of high relative humidity and low temperature. The correlation is significant with a slope lower than 1. By contrast, f is low and Δd is high in the environment of low relative humidity and high temperature. The influence of precipitation amount on the correlation between f and Δd is not obvious. Thus, under different climate background, it is necessary to consider the influence of meteorological factors on sub-cloud evaporation of atmospheric precipitation.

3.4.3. Influence of Raindrop Diameter

Raindrop diameters are an important parameter to calculate f and Δd in the Stewart model [24]. At present, the field observation and research data of raindrop diameter in Qilian Mountains are very minimal. Limited by technical means and field observation, raindrop diameter is often considered a constant in running the Stewart model [13,14,34]. In order to analyze the influence of raindrop diameter on f and Δd, we set raindrop diameter as a constant from 0.3 to 3 mm with a step-size of 0.1 mm for each sampling site. In order to calculate the value of f and Δd once again by Stewart model, the raindrop diameter is <1.5 mm, the influence on f and Δd is obvious, and the influence is slight when it is more than 2 mm (Figure 6a,b). According to the calculation result of precipitation intensity, the most frequent raindrop diameters in this study are below 1 mm. This was consistent with previous studies which observed the raindrop diameter on the northern slope of Qilian Mountains in 2006 [39], indicating that the raindrop diameter was usually lower than 1 mm with average raindrop diameter of 0.9 mm. In summary, the raindrop diameter, which can reflect spatial variability of precipitation, was a significant parameter for calculating the f and Δd based on the Stewart model. Therefore, setting raindrop diameter as a constant would affect the calculation accuracy in the study of sub-cloud evaporation in precipitation in the Qilian Mountains. In addition, under the conditions of the raindrop diameter being set as constant, the influence of altitude and underlying surface on f and Δd should not be ignored. Based on the analysis of this paper, we can assume that the raindrop diameter can be set as a constant input to the Stewart model when average raindrop diameter in the study area is higher than 2 mm. However, in arid and semi-arid areas, as well as mountainous areas and areas with different underlying surfaces, it is necessary to consider the influence of raindrop diameter on the output of the Stewart model.

4. Conclusions and Prospect

Using Stewart model, this study has quantitatively simulated the variation of f and Δd in precipitation from cloud-base to near surface at different elevation. The influence of meteorological factors, raindrop diameter and underlying surfaces has also discussed. Influenced by sub-cloud evaporation, the analytical results show that the d-excess value in precipitation gradually decreased from cloud-base to near surface. Δd increased with decreasing altitude, and evaporation remaining fraction decreased with decreasing altitude. The change of d-excess and f is most obvious in May. There is an obvious correlation between f and Δd. When f increases by 1%, Δd decreases to 0.99‰. When f is above 60%, the correlation between f and Δd is more significant, and that every 1% increase of could lead to Δd decrease about 0.81‰. However, the application of this method is still limited in an arid environment. Whether the linear relationship can be maintained at the environment of strong evaporation and low relative humidity needs further verification. The existence of the Xiying reservoir led to increased intensity of local moisture recycling and resulted in weakened sub-cloud evaporation in precipitation, increased f, and decreased Δd in the downstream. However, our study didn’t quantitatively analyze the influence of Xiying reservoir on local moisture recycling. In future studies, we will quantitatively analyze the intensity of local moisture recycling at different altitudes and compare it with that in the Xiying reservoir. In the environment of high relative humidity and low temperature, when f is high and Δd is close to 0, the slope of linear relationship between f and Δd is less than 0.99. By contrast, in the environment of low relative humidity and high temperature, f is low and Δd is high. The influence of precipitation amount on the relationship between them is not obvious. In semi-arid cold regions, the influence of raindrop diameter on f and Δd should be noticed. Setting raindrop diameter as a constant may affect the calculation accuracy in the study of sub-cloud evaporation of precipitation. The Stewart model has different parameter requirements in different study areas.
It is necessary to further improve the parameterization scheme in the simulation of sub-cloud evaporation of precipitation using the isotope method. It involves two points: one is whether the required parameters can be measured by technical means to reduce the uncertainty of estimation. The other is whether the algorithm can be optimized, if the original model is not fully suitable for arid environment, then how to debug it. Many parameters such as the raindrop diameter and the cloud-base height are treated by conceptual methods in the model. The differences of these parameters may affect the output results, so it is necessary to further optimize and improve them in the future. Due to the limitation of observation conditions, there are only four sampling sites in this study area, and the observation time is also short (data from May to September). It is necessary to optimize the sampling sites and collect more samples in the future.

Author Contributions

Data curation, X.M.; investigation, G.Z., W.J.; writing—original draft preparation, X.M.; writing—review and editing, W.J., G.Z., and S.W. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by Strategic priority research program of the Chinese Academy of Sciences (XDA19070503), the National Natural Science Foundation of China (41971036, 41867030), and Jiangxi Provincial Natural Science Foundation (20202BAB213013).

Acknowledgments

We would like to thank the colleagues in the Northwest Normal University for their help in writing process. We would like to thank to those who helped us in the experiment. We are grateful to anonymous reviewers and editorial staff for their constructive and helpful suggestions.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Li, Z.; Yuan, R.; Feng, Q.; Zhang, B.; Lv, Y.; Li, Y.; Wei, W.; Chen, W.; Ning, T.; Gui, J.; et al. Climate background, relative rate, and runoff effect of multiphase water transformation in Qilian Mountains, the third pole region. Sci. Total Environ. 2019, 663, 315–328. [Google Scholar] [CrossRef] [PubMed]
  2. IPCC. Summary for policymakers, in climate change 2013: The physical science basis. In Contribution of Working Group I to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change; Thomas, S., Stocker, F., Eds.; Cambridge University Press: Cambridge, UK, 2013; pp. 5–28. [Google Scholar]
  3. Ding, Y.H.; Wang, H.J. New understanding of the scientific issues of climate change in China during the past 100 years. Chin. Sci. Bull. 2016, 61, 1029–1041. (In Chinese) [Google Scholar] [CrossRef] [Green Version]
  4. Zhang, W.J.; Chen, W.M.; Li, B.L.; Tong, C.M.; Zhao, M.; Wang, N. Simulation of the permafrost distribution on Qilian Mountains over past 40 years under the influence of climate change. Geogr. Res. 2014, 33, 1275–1284. (In Chinese) [Google Scholar]
  5. Liang, P.B.; Li, Z.Q.; Zhang, H. Temporal-spatial variation characteristics of snow cover in Qilian Mountains from 2001 to 2017. Arid Zone Res. 2019, 42, 56–66. (In Chinese) [Google Scholar]
  6. Gao, Y.P.; Yao, X.J.; An, L.N.; Li, X.; Gong, P.; Qi, M. College of Geography and Environment Sciences, Northwest Normal University. Change of ice volume in the Qilian Mountains during the period from 2000 to 2010. Arid Zone Res. 2018, 35, 325–333. (In Chinese) [Google Scholar]
  7. Li, Z.X.; Feng, Q.; Li, Z.J.; Yuan, R.F.; Gui, J.; Lv, Y.M. Climate background, fact and hydrological effect of multiphase water transformation in cold regions of the Western China: A review. Earth-Sci. Rev. 2019, 190, 33–57. [Google Scholar]
  8. Dansgaard, W. The abundance of 18O in atmospheric water and water vapor. Tellus B 1953, 5, 461–469. [Google Scholar]
  9. Aggarwal, P.K.; Romatschke, U.; Araguas-Araguas, L.; Belachew, D.; Longstaffe, F.J.; Berg, P.; Schumacher, C.; Funk, A. Proportions of convective and stratiform precipitation revealed in water isotope ratios. Nat. Geosci. 2016, 9, 624–629. [Google Scholar] [CrossRef] [Green Version]
  10. Kern, Z.; Harmon, R.S.; Fórizs, I. Stable isotope signatures of seasonal precipitation on the Pacific coast of central Panama. Isot. Environ. Health Stud. 2015, 52, 1–13. [Google Scholar] [CrossRef] [Green Version]
  11. Dansgaard, W. Stable isotopes in precipitation. Tellus 1964, 16, 436–468. [Google Scholar] [CrossRef]
  12. Bershaw, J. Controls on deuterium excess across Asia. Geoscience 2018, 8, 257. [Google Scholar] [CrossRef] [Green Version]
  13. Salati, E.; Dall’Olio, A.; Matsui, E.; Gat, J.R. Recycling of water in the Amazon Basin: An isotopic study. Water Resour. Res. 1979, 15, 1250–1258. [Google Scholar] [CrossRef]
  14. Gat, J.R.; Bowser, C.J.; Kendall, C. The contribution of evaporation from the Great Lakes to the continental atmosphere: Estimate based on stable isotope data. Geophys. Res. Lett. 1994, 21, 557–560. [Google Scholar] [CrossRef]
  15. Froehlich, K.; Kralik, M.; Papesch, W.; Rank, D.; Scheifinger, H.; Stichler, W. Deuterium excess in precipitation of Alpine regions—Moisture recycling. Isot. Environ. Health Stud. 2008, 44, 61–70. [Google Scholar] [CrossRef] [PubMed]
  16. Peng, T.-R.; Liu, K.-K.; Wang, C.-H.; Chuang, K.-H. A water isotope approach to assessing moisture recycling in the island-based precipitation of Taiwan: A case study in the western Pacific. Water Resour. Res. 2011, 47, 08507. [Google Scholar] [CrossRef]
  17. Zhu, G.-F.; Li, J.; Shi, P.-J.; He, Y.-Q.; Cai, A.; Tong, H.-L.; Liu, Y.-F.; Yang, L. Relationship between sub-cloud secondary evaporation and stable isotope in precipitation in different regions of China. Environ. Earth Sci. 2016, 75, 876. [Google Scholar] [CrossRef]
  18. Pang, Z.; Kong, Y.; Froehlich, K.; Huang, T.; Yuan, L.; Li, Z.; Wang, F. Processes affecting isotopes in precipitation of an arid region. Tellus B 2011, 63, 352–359. [Google Scholar] [CrossRef] [Green Version]
  19. Salamalikis, V.; Argiriou, A.; Dotsika, E. Isotopic modeling of the sub-cloud evaporation effect in precipitation. Sci. Total Environ. 2016, 544, 1059–1072. [Google Scholar] [CrossRef]
  20. Breitenbach, S.F.M.; Adkins, J.F.; Mayer, H.; Marwan, N.; Kumar, K.K.; Haug, G.H. Strong influence of water vapor source dynamics on stable isotopes in precipitation observed in Southern Meghalaya, NE India. Earth Planet. Sci. Lett. 2010, 292, 212–220. [Google Scholar] [CrossRef]
  21. Aggarwal, P.K.; Alduchov, O.A.; Froehlich, K.O.; Araguas-Araguas, L.J.; Sturchio, N.C.; Kurita, N. Stable isotopes in global precipitation: A unified interpretation based on atmospheric moisture residence time. Geophys. Res. Lett. 2012, 39, 11705. [Google Scholar] [CrossRef]
  22. Stewart, M.K. Stable isotope fractionation due to evaporation and isotopic exchange of falling water drops: Applications to atmospheric processes and evaporation of lakes. J. Geophys. Res. 1975, 80, 1133–1146. [Google Scholar] [CrossRef]
  23. Zhang, X.; Xie, Z.; Yao, T. Mathematical modeling of variations on stable isotopic ratios in falling raindrops. Acta Meteorol. Sin. 1998, 12, 213–220. [Google Scholar]
  24. Peng, H.; Mayer, B.; Harris, S.; Krouse, H.R. The influence of below-cloud secondary effects on the stable isotope composition of hydrogen and oxygen in precipitation at Calgary, Alberta, Canada. Tellus B 2007, 59, 698–704. [Google Scholar] [CrossRef]
  25. Wu, H.; Zhang, X.; Xiaoyan, L.; Li, G.; Huang, Y. Seasonal variations of deuterium and oxygen-18 isotopes and their response to moisture source for precipitation events in the subtropical monsoon region. Hydrol. Process. 2014, 29, 90–102. [Google Scholar] [CrossRef]
  26. Wang, S.; Zhang, M.; Che, Y.; Zhu, X.; Liu, X. Influence of below-cloud evaporation on deuterium excess in precipitation of arid central asia and its meteorological controls. J. Hydrometeorol. 2016, 17, 1973–1984. [Google Scholar] [CrossRef]
  27. Chen, F.; Zhang, M.; Wang, S.; Ma, Q.; Zhu, X.; Dong, L. Relationship between sub-cloud secondary evaporation and stable isotopes in precipitation of Lanzhou and surrounding area. Quat. Int. 2015, 68–74. [Google Scholar] [CrossRef]
  28. McVicar, T.R.; Körner, C. On the use of elevation, altitude, and height in the ecological and climatological literature. Oecologia 2012, 171, 335–337. [Google Scholar] [CrossRef]
  29. Ma, X.G.; Jia, W.X.; Zhu, G.F.; Ding, D.; Pan, H.; Xu, X.; Guo, H.; Zhang, Y.; Yuan, R. Stable isotope composition of precipitation at different elevations in the monsoon marginal zone. Quat. Int. 2018, 493, 86–95. [Google Scholar]
  30. Zhu, G.; Guo, H.; Qin, D.; Pan, H.; Zhang, Y.; Jia, W.; Ma, X. Contribution of recycled moisture to precipitation in the monsoon marginal zone: Estimate based on stable isotope data. J. Hydrol. 2019, 569, 423–435. [Google Scholar] [CrossRef]
  31. Criss, R.E. Principles of Stable Isotope Distribution; Oxford University Press: New York, NY, USA, 1999. [Google Scholar]
  32. Barnes, S.L. An empirical shortcut to the calculation of temperature and pressure at the lifted condensation level. J. Appl. Meteorol. 1968, 7, 511. [Google Scholar] [CrossRef] [Green Version]
  33. Kinzer, G.D.; Gunn, R. The evaporation, temperature and thermal relaxation-time of freely falling waterdrops. J. Meteorol. 1951, 8, 71–83. [Google Scholar] [CrossRef] [Green Version]
  34. Best, A.C. Empirical formulae for the terminal velocity of water drops falling through the atmosphere. Q. J. R. Meteorol. Soc. 1950, 76, 302–311. [Google Scholar] [CrossRef]
  35. Kong, Y.; Pang, Z.; Froehlich, K. Quantifying recycled moisture fraction in precipitation of an arid region using deuterium excess. Tellus B Chem. Phys. Meteorol. 2013, 65, 9. [Google Scholar] [CrossRef]
  36. Berberan-Santos, M.N.; Bodunov, E.N.; Pogliani, L. On the barometric formula. Am. J. Phys. 1997, 65, 404–412. [Google Scholar] [CrossRef]
  37. Ma, Q.; Zhang, M.; Wang, S.; Wang, Q.; Liu, W.; Li, F.; Chen, F. An investigation of moisture sources and secondary evaporation in Lanzhou, Northwest China. Environ. Earth Sci. 2013, 71, 3375–3385. [Google Scholar] [CrossRef]
  38. Wang, H. Demonstration of limited water level of Xiying reservoir before flood. Gansu Agric. 2006, 11, 388–389. (In Chinese) [Google Scholar]
  39. Yu, G. Analysis of Raindrop Size-Distribution in Orographic Clouds in Qilian Mountains, China; Peking University: Beijing, China, 2008. (In Chinese) [Google Scholar]
Figure 1. The location of Xiying river basin and the distribution of sampling sites.
Figure 1. The location of Xiying river basin and the distribution of sampling sites.
Water 12 02798 g001
Figure 2. The variation of f and Δd of precipitation with altitude and time in the Xiying river basin from May to September 2017. (a): f varies with altitude; (b): f varies with time; (c): Δd varies with altitude; (d): Δd varies with time.
Figure 2. The variation of f and Δd of precipitation with altitude and time in the Xiying river basin from May to September 2017. (a): f varies with altitude; (b): f varies with time; (c): Δd varies with altitude; (d): Δd varies with time.
Water 12 02798 g002
Figure 3. Relationship between f and Δd in Xiying river basin from May to September 2017.
Figure 3. Relationship between f and Δd in Xiying river basin from May to September 2017.
Water 12 02798 g003
Figure 4. Changes of average cloud-base height, temperature and relative humidity for each sampling site in Xiying river basin from May to September 2017.
Figure 4. Changes of average cloud-base height, temperature and relative humidity for each sampling site in Xiying river basin from May to September 2017.
Water 12 02798 g004
Figure 5. Comparing the correlation between f and Δd at different environments in the Xiying river basin from May to September 2017. (ac): different temperature conditions; (df): different relative humidity conditions; (gi): different precipitation amount conditions. (T is temperature, RH is relative humidity, P is precipitation amount).
Figure 5. Comparing the correlation between f and Δd at different environments in the Xiying river basin from May to September 2017. (ac): different temperature conditions; (df): different relative humidity conditions; (gi): different precipitation amount conditions. (T is temperature, RH is relative humidity, P is precipitation amount).
Water 12 02798 g005
Figure 6. Relationship between raindrop diameter (D) and Δd or f at each sampling site in Xiying river basin. (a): diameter (D) and Δd. (b): raindrop diameter (D) and f.
Figure 6. Relationship between raindrop diameter (D) and Δd or f at each sampling site in Xiying river basin. (a): diameter (D) and Δd. (b): raindrop diameter (D) and f.
Water 12 02798 g006
Table 1. The variation of d-excess of precipitation at each sampling sites in the Xiying river basin from May to September 2017.
Table 1. The variation of d-excess of precipitation at each sampling sites in the Xiying river basin from May to September 2017.
Sampling SitesLatitudeLongitudeElevation (m)d-excess in Near Surface (‰)d-excess in Cloud-base (‰)Δd
(‰)
Lenglong37.55° N101.85° E360019.7240.01−20.29
Ningchang37.70° N101.89° E272113.1740.14−26.97
Huajian37.82° N102.01° E239010.2345.74−35.15
Xiying37.89° N102.18° E209713.9043.35−29.45
Table 2. Comparison of temperature, relative humidity and evaporation intensity at each sampling sites from May to September 2017.
Table 2. Comparison of temperature, relative humidity and evaporation intensity at each sampling sites from May to September 2017.
Surface Air Temperature (°C)Relative Humidity (%)Evaporation Intensity (g·s−1)
LenglongNingchangHuanjianXiyingLenglongNingchangHuanjianXiyingLenglongNingchangHuanjianXiying
May3.18.413.78.269.763.347.563.80.490.871.600.67
June5.38.312.878.277.659.30.440.500.95
July8.212.617.817.181.972.171.981.70.450.751.010.60
August7.613.017.116.584.879.766.477.40.370.561.010.65
September4.29.514.510.780.979.965.471.40.340.500.890.67
Table 3. Relationship of evaporation remaining fraction and Δd for each sampling site.
Table 3. Relationship of evaporation remaining fraction and Δd for each sampling site.
Sampling SiteLinear RegressionCoefficient of Determination (R2)
LenglongΔd = 0.83 f − 83.320.96
NingchangΔd = 0.98 f − 95.840.95
HuajianΔd = 1.38 f − 117.540.76
XiyingΔd = 0.82 f − 82.070.98

Share and Cite

MDPI and ACS Style

Ma, X.; Jia, W.; Zhu, G.; Wang, S. Quantitative Analysis of the Sub-Cloud Evaporation of Atmospheric Precipitation and Its Controlling Factors Calculated By D-Excess in an Inland River Basin of China. Water 2020, 12, 2798. https://doi.org/10.3390/w12102798

AMA Style

Ma X, Jia W, Zhu G, Wang S. Quantitative Analysis of the Sub-Cloud Evaporation of Atmospheric Precipitation and Its Controlling Factors Calculated By D-Excess in an Inland River Basin of China. Water. 2020; 12(10):2798. https://doi.org/10.3390/w12102798

Chicago/Turabian Style

Ma, Xinggang, Wenxiong Jia, Guofeng Zhu, and Shijin Wang. 2020. "Quantitative Analysis of the Sub-Cloud Evaporation of Atmospheric Precipitation and Its Controlling Factors Calculated By D-Excess in an Inland River Basin of China" Water 12, no. 10: 2798. https://doi.org/10.3390/w12102798

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop