三峡水库水面蒸发年内滞后效应分析

水面蒸发在季节变化上相比净辐射等气象要素的相位延迟反映了水体储热对水面蒸发的影响,量化这一对气象要素波动的滞后响应对理解和估算深水水库(湖泊)蒸发非常重要。三峡水库等河道型深水水库的水位和面积具有显著周期性变动,使得其水面蒸发的响应模式更为复杂,而目前对其认识非常薄弱。本文利用2013年8月—2020年7月三峡水库巴东站水面蒸发场和陆面蒸发场的蒸发和气象观测数据,分析了水面蒸发的季节变化及其年内滞后效应。结果表明:水面蒸发场蒸发量在年内的8和12月呈现双峰变化,与只有8月单峰值的陆面蒸发场蒸发量显著不同。水面蒸发场蒸发量相对净辐射、平均气温和水面温度分别存在4、3和2个月的滞后,而陆面蒸发场蒸发量相对滞后时间均在1个月以内;水面与陆面蒸发场相比,水温、蒸发量和水面与大气饱和差之间的滞后时间分别为1、3和4个月,而平均气温和净辐射之间不存在滞后。本文揭示出三峡水库巴东段水面蒸发年内滞后效应主要受到水库水温引起的水面与大气饱和差在季节上滞后的影响,需通过深入分析水温的时空变化来明确整个三峡水库的水面蒸发特征。     

Integration of Penman approach with complementary principle for evaporation research

Natural evaporation occurs with water transportation from an unsaturated land surface into an unsaturated atmosphere. The subprocesses at the land surface and in the atmosphere are one‐sidedly emphasized in the Penman approach and the complementary principle, in which the ratio of actual evaporation to the Penman potential evaporation is expressed as a function of the wetness state of the land surface and the atmosphere, respectively. The Penman approach and complementary principle can be integrated for completely conceptualizing the evaporation process, by expressing the evaporation ratio as a function of both the land surface and atmospheric wetness. The integrated approach has the potential to increase the accuracy of evaporation estimation while reducing the burdens of parameterization.     

Improvement of the Noah land surface model for warm season processes: evaluation of water and energy flux simulation

The Noah model is a land surface model of the National Centers for Environmental Prediction. It has been widely used in regional coupled weather and climate models (i.e. Weather Research and Forecasting Model, Eta Mesoscale Model) and global coupled weather and climate models (i.e. National Centers for Environmental Prediction Global Forecast System, Climate Forecast System). Therefore, its continued improvement and development are keys to enhancing our weather and climate forecast ability and water and energy flux simulation accuracy. North American Land Data Assimilation System phase 1 (NLDAS‐1) experiments indicated that the Noah model exhibited substantial bias in latent heat flux, total runoff and land skin temperature during the warm season, and such bias can significantly affect coupled weather and climate models. This paper presents a study to improve the Noah model by adding model parameterization processes such as including seasonal factor on leaf area index and root distribution and selecting optimal model parameters. We compared simulated latent heat flux, mean annual runoff and land skin temperature from the Noah control and test versions with measured latent heat flux, land surface skin temperature, mean annual runoff and satellite‐retrieved land surface skin temperature. The results show that the test version significantly reduces biases in latent heat, total runoff and land skin temperature simulation. The test version has been used for the NLDAS phase 2 (NLDAS‐2) to produce 30‐year water flux, energy flux and state variable products to support the US drought monitor of National Integrated Drought Information System. Copyright © 2012 John Wiley & Sons, Ltd.     

Simulating daily, monthly and annual water balances in a land surface model using alternative root water uptake schemes

Hydrological simulations at multi-temporal time scales by a widely used land surface model (LSM) are investigated under contrasting vegetation and meteorological conditions. Our investigation focuses particularly on the effects of two different representations of root water uptake and root profile on simulated evapotranspiration (ET) and soil moisture by the Integrated BIosphere Simulator (IBIS). For this purpose, multi-year eddy covariance measurements, collected at four flux-tower sites across North America, were used to gauge IBIS simulations with: (a) its standard version (IBIS2.1), in which static root water uptake (RWU) and root profile schemes are incorporated; and (b) a modified version in which dynamic RWU and root profile schemes replaces the static schemes used in the standard version. Overall, our results suggest that the modified version of the model performs more realistically than the standard version, particularly when high atmospheric demand for evaporation is combined with high atmospheric vapour pressure deficit and low soil water availability. The overall correlation between simulated and measured monthly ET rates at the simulated sites reached 0.87 and 0.91 for the standard and the modified versions, respectively. Our results also show that the incorporation of the dynamic RWU in IBIS yields improved simulations of ET under very dry conditions, when soil moisture falls down to very low levels. This suggests that adequate representations of vegetation responses to drought are needed in LSMs as many state of the art climate models projections of future climate indicate more frequent and/or more intense drought events occurring in some regions of the globe. Our analysis also highlighted the urgent need for adequate methodologies to correct field measurements that exhibit energy imbalances in order to provide rigorous assessments of land surface model simulations of heat and mass exchanges between the land surface and the atmosphere.     

A nonlinear function approach for the normalized complementary relationship evaporation model

A nonlinear function approach for the normalized complementary relationship evaporation model that is different from the methodology maintaining the symmetric complementary relationship with appropriate definitions of potential and wet‐environment evaporation is proposed and verified. This approach employs the definitions used in the advection‐aridity model, wherein the potential is estimated using the Penman equation. Normalized by Penman potential evaporation, the complementary relationship model is expressed as a function describing the relationship between the evaporation ratio (the ratio of the actual to the Penman potential evaporation) and the proportion of the radiation term in Penman potential evaporation. The new nonlinear function proposed in the current study is approximately equivalent to the advection‐aridity and the modified Granger models under conditions that are neither too wet nor too dry, but is more reasonable under arid and wet conditions. The new nonlinear function model performs well in estimating actual evaporation, as verified by the observed data from four sites under different land covers. Copyright © 2011 John Wiley & Sons, Ltd.     

Comparison of 15 evaporation methods applied to a small mountain lake in the northeastern USA

Few detailed evaporation studies exist for small lakes or reservoirs in mountainous settings. A detailed evaporation study was conducted at Mirror Lake, a 0.15 km² lake in New Hampshire, northeastern USA, as part of a long-term investigation of lake hydrology. Evaporation was determined using 14 alternate evaporation methods during six open-water seasons and compared with values from the Bowen-ratio energy-budget (BREB) method, considered the standard. Values from the Priestley–Taylor, deBruin–Keijman, and Penman methods compared most favorably with BREB-determined values. Differences from BREB values averaged 0.19, 0.27, and 0.20 mm d⁻¹, respectively, and results were within 20% of BREB values during more than 90% of the 37 monthly comparison periods. All three methods require measurement of net radiation, air temperature, change in heat stored in the lake, and vapor pressure, making them relatively data intensive. Several of the methods had substantial bias when compared with BREB values and were subsequently modified to eliminate bias. Methods that rely only on measurement of air temperature, or air temperature and solar radiation, were relatively cost-effective options for measuring evaporation at this small New England lake, outperforming some methods that require measurement of a greater number of variables. It is likely that the atmosphere above Mirror Lake was affected by occasional formation of separation eddies on the lee side of nearby high terrain, although those influences do not appear to be significant to measured evaporation from the lake when averaged over monthly periods.     

A complementary relationship evaporation model referring to the Granger model and the advection–aridity model

A complementary relationship evaporation model has been proposed and verified based on evaluations of the advection–aridity model and the Granger's complementary relationship model (Granger model) in dimensionless forms. Normalized by Penman potential evaporation, the Granger model and the advection–aridity model have been transformed into similar dimensionless forms. Evaporation ratio (ratio of actual evaporation to Penman potential evaporation) has been expressed as a function of dimensionless variable based on radiation and atmospheric conditions. Similar dimensionless variables for the different functions have been used in the two models. By referring to the dimensionless variable from the advection–aridity model and the function from the Granger model, a new model to estimate actual evaporation was proposed. The performance of the new model has been validated by the observed data from four sites under different land covers. The new model is an enhanced Granger model with better evaporation prediction over the aforementioned different land covers. It also offers more stable optimized parameters in a grassland site than the Granger model. The new model somewhat approximates the advection–aridity model under neither too wet nor too dry conditions, but without its system bias. Copyright © 2011 John Wiley & Sons, Ltd.     

Validity of the Bouchet’s complementary relationship at 102 observatories across China

Bouchet in 1963 hypothesized that for large homogeneous land surface with minimum advection of heat and moisture, there exists a 1:1 complementary relationship of potential and actual evaporation coupled through land-atmosphere feedbacks. The complementary relationship has been widely used to estimate regional actual evaporation and explain the pan evaporation paradox. We examine the standardized potential evaporation (potential evaporation divided by wet environment evaporation) at 102 observatories at different elevations across China. Generally, the relationship is appropriate at the low elevations (<1000 m). With the increase of elevation, vapor transfer power becomes much less than radiation energy budget because of lower vapor pressure deficit and stronger global solar radiation. As a result, at the high elevations (over 1000 m), the excess energy resulted by limited moisture availability is not enough to be converted into drying power of the air. This result suggests that the complementary relationship is asymmetric at the high elevations. Supported by the Presidential Special Award Foundation, the Chinese Academy of Sciences (Grant No. O7R70020SD) and the National Key Technology R & D Program (Grant No. 2006BAC08B0408)     

The simulation of heat and water exchange in the boreal spruce forest by the land-surface model SWAP

All previous versions of a physically based land-surface model SWAP have assumed for simplicity that vegetation is fully covered by snow during the cold season. Such assumption is reasonable only for the regions dominated by short vegetation or for warm climates where snow processes are absent. The major goals of this paper are (i) modification of the latest version of SWAP by incorporation of tall vegetation into the cold-season parameterizations to make the model applicable for simulating heat and water transfer within a boreal forest biome and (ii) validation of the modified version using the data from a forested catchment located in the boreal zone. Modification of SWAP required to parameterize radiative and turbulent exchange between the forest crown and forest floor, partitioning snowfall between interception by the canopy (in doing so, snow interception differs from rain interception) and falling to the ground, formation of snow cover on the forest crown and forest floor including snow accumulation (both in solid and liquid fractions), snow evaporation, and snowmelt. The advanced model was validated using a set of hydrometeorological data measured during 18 years (1966–1983) at the Tayozhniy catchment (covered by boreal spruce forest), Valdai, Russia. Simulations of annual and monthly snow/rain interception, daily runoff at the catchment outlet, snow density, snow depth, snow water equivalent, soil water storage in three layers (0–20, 0–50 and 0–100 cm), and monthly evapotranspiration from the catchment were compared with observations. Analysis of the results of validation shows that the new version of the model SWAP reproduces the heat and water exchange processes occurring in mid-latitude boreal forest quite reasonable.     

Assessing the Sensitivity of a Model of Runoff Formation in the Ussuri River Basin

A physically based model of runoff formation with daily resolution has been developed for the upper part of the Ussuri basin with an area of 24400 km² based on ECOMAG hydrological modeling platform. Two versions of the hydrological model have been studied: (1) a crude version with the spatial schematization of the drainage area and river network based on DEM 1 × 1 km with the use of soil and landscape maps at a scale of 1: 2500000 and (2) a detailed version with DEM 80 × 80 m and soil and landscape maps of the scale of 1: 100000. Each version of the model has been tested for two variants of meteorological inputs: (1) meteorological forcing data (temperature, air humidity, precipitation) at eight weather stations and (2) with the involvement of additional data on precipitation collected at 15 gages in the basin. The model has been calibrated and validated over a 34-year period (1979–2012) with the use of runoff data for the Ussuri R. and its tributaries. The results of numerical experiments for assessing the sensitivity of model hydrological response to the spatial resolution of land surface characteristics and the density of precipitation gaging stations are discussed.     

Effects of soil moisture aggregation on surface evaporative fluxes

The effects of small-scale heterogeneity in land surface characteristics on the large-scale fluxes of water and energy in the land-atmosphere system has become a central focus of many climatology research experiments. The acquisition of high resolution land surface data through remote sensing and intensive land-climatology field experiments (like HAPEX, FIFE, and BOREAS) has provided data to investigate the interactions between microscale land-atmosphere interactions and macroscale models. To determine the effect of small scale heterogeneities, the spatially averaged evaporative fraction is analytically derived for spatially variable soil moisture and soil-atmospheric controls on evaporation at low soil moisture. This average evaporative fraction is compared with the evaporative fraction determined using the spatially averaged soil moisture, as if from a lumped, or aggregated, land surface model. Results show that the lumped-model based evaporation will over estimate evaporation during periods of low atmospheric demands (early morning/late afternoon, Winter periods, etc.) and under estimate evaporation during periods of high demand (midday Summer periods.) The accuracy of using ‘effective’ parameters in lumped macroscale models depends on the variability of soil moisture and the sensitivity of the soil-vegetation system to low soil moisture.     

Improvement of snowpack simulations in a regional climate model

To improve simulations of regional‐scale snow processes and related cold‐season hydroclimate, the Community Land Model version 3 (CLM3), developed by the National Center for Atmospheric Research (NCAR), was coupled with the Pennsylvania State University/NCAR fifth‐generation Mesoscale Model (MM5). CLM3 physically describes the mass and heat transfer within the snowpack using five snow layers that include liquid water and solid ice. The coupled MM5–CLM3 model performance was evaluated for the snowmelt season in the Columbia River Basin in the Pacific Northwestern United States using gridded temperature and precipitation observations, along with station observations. The results from MM5–CLM3 show a significant improvement in the SWE simulation, which has been underestimated in the original version of MM5 coupled with the Noah land‐surface model. One important cause for the underestimated SWE in Noah is its unrealistic land‐surface structure configuration where vegetation, snow and the topsoil layer are blended when snow is present. This study demonstrates the importance of the sheltering effects of the forest canopy on snow surface energy budgets, which is included in CLM3. Such effects are further seen in the simulations of surface air temperature and precipitation in regional weather and climate models such as MM5. In addition, the snow‐season surface albedo overestimated by MM5–Noah is now more accurately predicted by MM5–CLM3 using a more realistic albedo algorithm that intensifies the solar radiation absorption on the land surface, reducing the strong near‐surface cold bias in MM5–Noah. The cold bias is further alleviated due to a slower snowmelt rate in MM5–CLM3 during the early snowmelt stage, which is closer to observations than the comparable components of MM5–Noah. In addition, the over‐predicted precipitation in the Pacific Northwest as shown in MM5–Noah is significantly decreased in MM5–CLM3 due to the lower evaporation resulting from the longer snow duration. Copyright © 2011 John Wiley & Sons, Ltd.     

Lateral terrestrial water flow contribution to summer precipitation at continental scale – A comparison between Europe and West Africa with WRF-Hydro-tag ensembles

It is well accepted that summer precipitation can be altered by soil moisture condition. Coupled land surface – atmospheric models have been routinely used to quantify soil moisture – precipitation feedback processes. However, most of the land surface models (LSMs) assume a vertical soil water transport and neglect lateral terrestrial water flow at the surface and in the subsurface, which potentially reduces the realism of the simulated soil moisture – precipitation feedback. In this study, the contribution of lateral terrestrial water flow to summer precipitation is assessed in two different climatic regions, Europe and West Africa, for the period June–September 2008. A version of the coupled atmospheric-hydrological model WRF-Hydro with an option to tag and trace land surface evaporation in the modelled atmosphere, named WRF-Hydro-tag, is employed. An ensemble of 30 simulations with terrestrial routing and 30 simulations without terrestrial routing is generated with random realizations of turbulent energy with the stochastic kinetic energy backscatter scheme, for both Europe and West Africa. The ensemble size allows to extract random noise from continental-scale averaged modelled precipitation. It is found that lateral terrestrial water flow increases the relative contribution of land surface evaporation to precipitation by 3.6% in Europe and 5.6% in West Africa, which enhances a positive soil moisture – precipitation feedback and generates more uncertainty in modelled precipitation, as diagnosed by a slight increase in normalized ensemble spread. This study demonstrates the small but non-negligible contribution of lateral terrestrial water flow to precipitation at continental scale.     

Analysis of long‐term water balance in the source area of the Yellow River basin

To analyse the long‐term water balance of the Yellow River basin, a new hydrological model was developed and applied to the source area of the basin. The analysis involved 41 years (1960–2000) of daily observation data from 16 meteorological stations. The model is composed of the following three sub‐models: a heat balance model, a runoff formation model and a river‐routing network model. To understand the heat and water balances more precisely, the original model was modified as follows. First, the land surface was classified into five types (bare, grassland, forest, irrigation area and water surface) using a high‐resolution land‐use map. Potential evaporation was then calculated using land‐surface temperatures estimated by the heat balance model. The maximum evapotranspiration of each land surface was calculated from potential evaporation using functions of the leaf area index (LAI). Finally, actual evapotranspiration was estimated by regulating the maximum evapotranspiration using functions of soil moisture content. The river discharge estimated by the model agreed well with the observed data in most years. However, relatively large errors, which may have been caused by the overestimation of surface flow, appeared in some summer periods. The rapid decrease of river discharge in recent years in the source area of the Yellow River basin depended primarily on the decrease in precipitation. Furthermore, the results suggested that the long‐term water balance in the source area of the Yellow River basin is influenced by land‐use changes. Copyright © 2007 John Wiley & Sons, Ltd.     

水面蒸发与散热系数公式研究(一)

根据自1976年以来全国水面蒸发与散热研究协作组在我国各典型地区的原体与室内实验资料和大量水文站历史资料,通过理论分析和统计检验,确定丁影响水面蒸发的诸因子及其非线性相互作用,引入了新的无量纲参数(w_e、Pv、Pe)和公式结构,用实测资料统计确定厂公式中的常系数,得到了用开敞湖面一般水文气象资料计算逐日蒸发和散热系数的公式。经全国各典型气候带内各季节湖泊(水库)和受热污染水体上原体观测和室内专题实验共1860组口平均资料检验,公式的精度高于现有其他公式。全文分两部分,这是第一部分,内容包括:影响水面蒸发的土要无量纲参数;感热输送和大气饱和度对蒸发影响的修正;水面蒸发计算公式的结构及其经验系数。     

Modeling the components of heat and water balance for the land surface of the globe

Physically-based model SWAP, developed by the authors earlier and describing the processes of heat and water exchange between the land surface and the atmosphere was adapted for calculating the components of heat and water balance for the entire land surface of the globe. An information base for the model was prepared as a version of global dataset with one-degree spatial resolution for three-hour hydrometeorological data and land surface parameters. The dynamics of various parameters of heat and water regimes of the soil-vegetation (snow) cover-atmosphere system was calculated by using a new version of the land surface model SWAP with a three-hour time step from July 1, 1982, to December 31, 1995. Calculation results were compared with estimates available from the literature.     

Examination of Morton's CRAE model for estimating daily evaporation from field-sized areas

Morton's complementary relationship areal evapotranspiration (CRAE) model was originally designed to provide regional estimates of monthly evapotranspiration. Often, however, hydrologists and others require estimates of evapotranspiration for field-sized land units under a specific land use, for shorter intervals of time. This paper examines CRAE with respect to the algorithms used to describe different terms and its applicability to reduced spatial and temporal scales.

Daily estimates by CRAE of atmospheric radiation fluxes during the summer months are compared with monitored values. It is shown that errors in estimation of the extra-terrestrial flux, the transmittancy of clouds to short-wave radiation, the surface albedo and the net long-wave flux result in standard deviations of the difference between ‘modelled’ and ‘measured’ net all-wave radiation for 1-, 5- and 10-day periods of 2.58, 1.8 and 1.50 MJm⁻² day⁻¹ respectively.

The assumption in CRAE that the vapour transfer coefficient is independent of wind speed may lead to appreciable error in computing evapotranspiration. A procedure for incorporating a wind correction factor is described and the improvement in estimating regional evaporation is illustrated.

Comparisons of evapotranspiration estimates by CRAE and measurements obtained from soil moisture and precipitation observations in the semi-arid, cold-climate Prairie region of western Canada demonstrate that the assumptions that the soil heat flux and storage terms are negligible, lead to large overestimation by the model during periods of soil thaw.     

A refinement of the combination equations for evaporation

Most combination equations for evaporation rely on a linear expansion of the saturation vapor-pressure curve around the air temperature. Because the temperature at the surface may differ from this temperature by several degrees, and because the saturation vapor-pressure curve is nonlinear, this approximation leads to a certain degree of error in those evaporation equations. It is possible, however, to introduce higher-order polynomial approximations for the saturation vapor-pressure curve and to derive a family of explicit equations for evaporation, having any desired degree of accuracy. Under the linear approximation, the new family of equations for evaporation reduces, in particular cases, to the combination equations of H. L. Penman (Natural evaporation from open water, bare soil and grass,Proc. R. Soc. London, Ser. A 193, 120–145, 1948) and of subsequent workers. Comparison of the linear and quadratic approximations leads to a simple approximate expression for the error associated with the linear case. Equations based on the conventional linear approximation consistently underestimate evaporation, sometimes by a substantial amount.     

Influence of condensation and latent heat release upon barotropic and baroclinic instabilities of vortices in a rotating shallow water f-plane model

Analysis of the influence of condensation and related latent heat release upon developing barotropic and baroclinic instabilities of large-scale low Rossby-number shielded vortices on the f-plane is performed within the moist-convective rotating shallow water model, in its barotropic (one-layer) and baroclinic (two-layer) versions. Numerical simulations with a high-resolution well-balanced finite-volume code, using a relaxation parameterisation for condensation, are made. Evolution of the instability in four different environments, with humidity (i) behaving as passive scalar, (ii) subject to condensation beyond a saturation threshold, (iii) subject to condensation and evaporation, with three different parameterisations of the latter, are inter-compared. The simulations are initialised with unstable modes determined from the detailed linear stability analysis in the “dry” version of the model. In a configuration corresponding to low-level mid-latitude atmospheric vortices, it is shown that the known scenario of evolution of barotropically unstable vortices, consisting in formation of a pair of dipoles (dipolar breakdown) is substantially modified by condensation and related moist convection, especially in the presence of surface evaporation. No enhancement of the instability due to precipitation was detected in this case. Cyclone-anticyclone asymmetry with respect to sensitivity to the moist effects is evidenced. It is shown that inertia-gravity wave emission during the vortex evolution is enhanced by the moist effects. In the baroclinic configuration corresponding to idealised cut-off lows in the atmosphere, it is shown that the azimuthal structure of the leading unstable mode is sensitive to the details of stratification. Scenarios of evolution are completely different for different azimuthal structures, one leading to dipolar breaking, and another to tripole formation. The effects of moisture considerably enhance the perturbations in the lower layer, especially in the tripole formation scenario.