Evaluating the SEBS‐estimated evaporative fraction from MODIS data for a complex underlying surface

With the complex nature of land surfaces, more attention should be paid to the performance of remotely sensed models to estimate evapotranspiration from moderate and low spatial resolution data. Taking into account the characteristic of a stable evaporative fraction (EF) in the daytime, this paper uses the surface energy balance system (SEBS) to estimate the EF from MODIS data for a subtropical evergreen coniferous plantation in southern China and evaluates the stability of the SEBS model in estimating the EF under complex surface conditions. The results show that the SEBS‐estimated EF is larger than the measured EF partly because of the serious lack of energy‐balance closure. This difference can be largely reduced by the residual energy correction method. More evaporative land cover within the MODIS pixel is a main reason for the overestimated EF. SEBS underestimates sensible heat flux, and the underestimation of surface available energy also contributes to the overestimation of the EF. The EF estimated from MODIS/Terra data is in agreement with that from MODIS/Aqua data with a coefficient of determination (R²) of 0.552, a mean bias error (BIAS) of 0.028, and a root mean square error (RMSE) of 0.079, which is consistent with the result from in situ measurements. In addition, the estimation of surface available energy from remotely sensed data is evaluated on this complex underlying surface. Compared with in situ measurements, the available energy is underestimated by 28 W m^?2 with an RMSE = 50 W m^?2 and an R² = 0.87. Copyright © 2012 John Wiley & Sons, Ltd.     

Modelling snowpack surface temperature in the Canadian Prairies using simplified heat flow models

Three practical schemes for computing the snow surface temperature T_s, i.e. the force–restore method (FRM), the surface conductance method (SCM), and the Kondo and Yamazaki method (KYM), were assessed with respect to T_s retrieved from cloud‐free, NOAA‐AVHRR satellite data for three land‐cover types of the Paddle River basin of central Alberta. In terms of R², the mean T_s, the t‐test and F‐test, the FRM generally simulated more accurate T_s than the SCM and KYM. The bias in simulated T_s is usually within several degrees Celsius of the NOAA‐AVHRR T_s for both the calibration and validation periods, but larger errors are encountered occasionally, especially when T_s is substantially above 0 °C. Results show that the simulated T_s of the FRM is more consistent than that of the SCM, which in turn was more consistent than that of the KYM. This is partly because the FRM considers two aspects of heat conduction into snow, a stationary‐mean diurnal (sinusoidal) temperature variation at the surface coupled to a near steady‐state ground heat flux, whereas the SCM assumes a near steady‐state, simple heat conduction, and other simplifying assumptions, and the KYM does not balance the snowpack heat fluxes by assuming the snowpack having a vertical temperature profile that is linear. Copyright © 2005 John Wiley & Sons, Ltd.     

Treatment of anchor pixels in the METRIC model for improved estimation of sensible and latent heat fluxes

Abstract

Reliable estimation of sensible heat flux (H) is important in energy balance models for quantifying evapotranspiration (ET). This study was conducted to evaluate the value of adding the Priestley-Taylor (PT) equation to the METRIC (Mapping Evapotranspiration at high Resolution with Internalized Calibration) model. METRIC was used to estimate energy fluxes for 10 Landsat images from the 2005, 2006 and 2007 crop growing seasons in south-central Nebraska, USA, where each image owing to recent rainfall exhibited high residual moisture content even at the hot pixel. The METRIC model performed satisfactorily for net radiation (R_n ) and soil heat flux (G) estimation with a root mean square error (RMSE) of 52 and 24 W m^-2, respectively. A RMSE of 122 W m^-2 for H indicated the limitation of the METRIC model in estimating H for high residual moisture content of the hot pixel (Alfalfa reference ET fraction, ET _r F > 0.15). The modified METRIC model (wet METRIC or wMETRIC) incorporating the PT equation was applied to calculate H at the anchor pixels (hot and cold) for high residual moisture content of the hot pixel. The α coefficient of the PT equation was locally calibrated using hourly meteorological data from an automatic weather station and R_n and G data from a Bowen ratio flux tower. The mean α coefficient value was 1.14. The wMETRIC model reduced the RMSE of H from 122 to 64 W m^-2 and that of latent heat flux, LE, from 163 to 106 W m^-2. The RMSE of daily ET decreased from 1.7 to 1.1 mm d^-1 with wMETRIC. The results indicate that treatment of anchor pixels for high residual moisture content with the PT approach gives improved estimation of H, LE and daily ET. It is recommended that the wMETRIC model be used for estimating ET if the hot pixel has high residual moisture (i.e. reference ET fraction > 0.15).

Citation Singh, R. K. & Irmak, A. (2011) Treatment of anchor pixels in the METRIC model for improved estimation of sensible and latent heat fluxes. Hydrol. Sci. J. 56(5), 895–906.     

Further evaluation of the Sim-ReSET model for ET estimation driven by only satellite inputs

Abstract

A simple remote sensing evapotranspiration (ET) model (Sim-ReSET) has been proposed but only tested using field measurements at a site with a semi-arid climate. Its performance for mapping ET using only satellite data remained unknown. In this study, the Sim-ReSET model was further evaluated for ET estimation driven by only MODIS data products. The estimated ET rates were compared with ground-based observational data from a variety of ecosystems and climates across China. The results show that MODIS-based ET estimates are consistent with both the ET measurements from eddy covariance flux towers and those from the Penman-Monteith method combined with micrometeorological data. Evaporation fraction (EF) is indicative of land surface moisture. The derivative EF maps demonstrate that the proposed ET data set obtained from the Sim-ReSET model and MODIS data is capable of capturing the spatio-temporal pattern of land surface moisture for different land covers with different climates.

Editor Z.W. Kundzewicz

Citation Sun, Z.G., Wang, Q.X., Matsushita, B., Fukushima, T., Ouyang, Z., Watanabe, M., and Gebremichael, M., 2013. Further evaluation of the Sim-ReSET model for ET estimation driven by only satellite inputs. Hydrological Sciences Journal, 58 (5), 994–1012.     

Spatial variability of the surface energy balance of Lake Kasumigaura and implications for flux measurements

ABSTRACT

The spatial variability of the lake surface energy balance and its causes are not well-understood. Energy balance maps (90 m resolution) of Lake Kasumigaura (172 km²), Japan, obtained by interpolating station data and bulk equations, allowed an investigation of these issues. Due to lake-scale variations in meteorological variables and small-scale fluctuations of surface temperature, T_s, surface heat fluxes differed horizontally at two distinct scales, while radiative fluxes were more uniform. As the key variable to surface flux T_s was only homogeneous for directions with a longer fetch or under calm wind conditions. Using these findings, the suitability of two flux station locations, one at the centre of the lake and another within a cove, was considered. Although both locations satisfied the fetch requirements, T_s was not always found to be homogeneous in the cove, making this location less suitable for flux measurements, an issue that, to date, has been overlooked.     

A modified trapezoid framework model for partitioning regional evapotranspiration

While evapotranspiration (ET) is normally measured as one hydrologic component, evaporation (E), and transpiration (T) result from different physical-biological processes. Using a two-source model, a trapezoid framework has been widely applied in recent years. The key to applying the trapezoid framework model is the determination of the dry/wet boundaries of the land surface temperature-fractional vegetation coverage trapezoid (LST-f_c). Although algorithms have been developed to characterize the two boundaries, there remains a significant uncertainty near the wet boundary which scatters in a discrete and uneven manner. It is therefore difficult to precisely locate the wet boundary. To address this problem, a Wet Boundary Algorithm (WBA) was developed in this study with the algorithm applied in the region of Huang-Huai-Hai plain of China, using the Pixel Component Arranging and Comparing Algorithm (PCACA) to retrieve ET from MODerate-resolution Imaging Spectroradiometer (MODIS) Data. The eddy covariance (EC) measurements from Yucheng station was used to verify the modified model where the root mean square error (RMSE) of 17.8 W/m², Bias of −7.2 W/m² for latent heat flux (LE) simulation in 28 cloudless test days. The ratio of transpiration to evapotranspiration (T/ET) varied between 0.48 and 0.81 over the Huang-Huai-Hai plain. The spatial and temporal distribution of ET revealed that agriculture practices have a significant influence on the hydrological cycle, where crop growth promotes the magnitude of ET. Likewise, harvesting activities significantly reduce ET. The proposed WBA algorithm significantly reduces the uncertainty of the trapezoid ET model caused by wet edge positioning. The analysis of the impact of agricultural activities on ET provide a better understanding how human activities change the hydrological cycle at regional scales.     

Coupling SEBAL with a new radiation module and MODIS products for better estimation of evapotranspiration

ABSTRACT

Evapotranspiration (ET) is an important ecohydrological process especially in arid and semi-arid regions. In this study, a new radiation module based on MODIS data has been coupled with the Surface Energy Balance Algorithms for Land (SEBAL) to better estimate ET. The accuracies of the coupled model for estimating available energy and sensible heat (H) were improved significantly compared with the outputs from the original SEBAL which was based on empirical equations. The coupled SEBAL modelled instantaneous λET agreed much better with observations in the arid land of Central Asia than the original SEBAL, with a bias of ?2.86 W m^-2, root mean square error (RMSE) of 9.75 W m^-2, and normalized RMSE (NRMSE) of 0.13. The accuracy was blurred when scaling ET to a daily or monthly scale, mainly due to the uncertainties associated with temporal upscaling methods that were applied. Sensitivity analysis, which was conducted using numerical variance-based techniques, indicated that the estimated ET is sensitive to the available energy, suggesting the importance of obtaining accurate estimates of net radiation when applying the coupled SEBAL to estimate ET. This study provides a simple and reliable way to utilize MODIS products and contains sensitivity analysis for helping to correctly interpret the outputs, which are both important for large-scale ET estimation.     

CO2 absorption of sandy soil induced by rainfall pulses in a desert ecosystem

Soil CO₂ flux is strongly influenced by precipitation in many ecosystem types, yet knowledge of the effects of precipitation on soil CO₂ flux in semi‐arid desert ecosystems remains insufficient, particularly for sandy soils. To address this, we investigated the response of sandy soil CO₂ flux to rainfall pulses in a desert ecosystem in northern China during August–September 2011. Significant changes (P < 0.05) were found in diel patterns of soil CO₂ flux induced by small (2.1 mm), moderate (12.4 mm) and large (19.7 mm) precipitation events. Further analysis indicated that rainfall pulses modified the response of soil CO₂ flux to soil temperature, including hysteresis between soil CO₂ flux and soil temperature, with F_s higher when T_s was increasing than when T_s was decreasing, and the linear relationship between them. Moreover, our results showed that rainfall could result in absorption of atmospheric CO₂ by soil, possibly owing to mass flow of CO₂ induced by a gradient of gas pressure between atmosphere and soil. After each precipitation event, soil CO₂ flux recovered exponentially to pre‐rainfall levels with time, with the recovery times exhibiting a positive correlation with precipitation amount. On the basis of the amounts of precipitation that occurred at our site during the measurement period (August–September), the accumulated rain‐induced carbon absorption evaluated for rainy days was 1.068 g C m^?2; this corresponds approximately to 0.5–2.1% of the net primary production of a typical desert ecosystem. Thus, our results suggest that rainfall pulses can strongly influence carbon fluxes in desert ecosystems. Copyright © 2014 John Wiley & Sons, Ltd.     

大型浅水湖泊与大气之间的动量和水热交换系数——以太湖为例

湖泊水面与大气之间垂直方向的动量、水汽和热量通量与风速、湿度和温度梯度之间存在比例关系,因此在湖泊水-气相互作用研究中,这比例系数(交换系数)是关键因子.在以往的研究中,交换系数通常直接采用水面梯度观测法或海洋大气近地层的参数化方案进行计算.本文采用涡度相关系统和小气候系统仪器在太湖平台上直接观测的通量和气象要素,对上述交换系数(最小均方差原则)进行优化,结果为:动量交换系数C_D10N=1.52×10^-3、水汽交换系数C_E10N=0.82×10^-3、热量交换系数CH10N=1.02×10^-3,与其他内陆湖泊涡度相关观测数据的推导结果一致.本文的研究结果表明:与海洋参数化方案相比,在相同的风速条件下,湖面的空气动力学粗糙度比海洋高,这可能是由于受到水深的影响;如果采用海洋参数化方案,会导致湖泊年蒸发量的估算值偏大40%.太湖的动量、水汽和热量交换系数可以视为常数,可以不考虑稳定度和风速的影响.这是因为本文中83%的数据为近中性条件.敏感性分析表明:如果考虑稳定度的影响,LE模拟值的平均误差降低了0.5 W/m²,H的平均误差降低了0.4 W/m²,u*的计算值没有变化;如果考虑风速的影响,u*模拟值的平均误差降低了0.004 m/s,LE的平均误差升高了1.3 W/m²,H的模拟结果几乎不受影响.这一结果能为湖气相互作用研究提供参考.     

A detailed energy budget analysis of river supercooling and the importance of accurately quantifying net radiation to predict ice formation

River supercooling and ice formation is a regular occurrence throughout the winter in northern countries. The resulting frazil ice production can obstruct the flow through intakes along the river, causing major problems for hydropower and water treatment facilities, among others. Therefore, river ice modellers attempt to calculate the river energy budget and predict when supercooling will occur in order to anticipate and mitigate the effects of potential intake blockages. Despite this, very few energy budget studies have taken place during freeze-up, and none have specifically analysed individual supercooling events. To improve our understanding of the freeze-up energy budget detailed measurements of air temperature, relative humidity, barometric pressure, wind speed and direction, short- and longwave radiation, and water temperature were made on the Dauphin River in Manitoba. During the river freeze-up period of late October to early November 2019, a total of six supercooling events were recorded. Analysis of the energy budget throughout the supercooling period revealed that the most significant heat source was net shortwave radiation, reaching up to 298 W/m², while the most significant heat loss was net longwave radiation, accounting for losses of up to 135 W/m². Longwave radiation was also the most significant heat flux overall during the individual supercooling events, accounting for up to 84% of the total heat flux irrespective of flux direction, highlighting the importance of properly quantifying this flux during energy budget calculations. Five different sensible (Q_h) and latent (Q_e) heat flux calculations were also compared, using the bulk aerodynamic method as the baseline. It was found that the Priestley and Taylor method most-closely matched the bulk aerodynamic method on a daily timescale with an average offset of 8.5 W/m² for Q_h and 10.1 W/m² for Q_e, while a Dalton-type equation provided by Webb and Zhang was the most similar on a sub-daily timescale with average offsets of 20.0 and 14.7 W/m² for Q_h and Q_e, respectively.     

On the heat flux related to stretching in the NW-Mediterranean continental margins

Summary The surface thermal flux of the continental margins of the northwestern Mediterranean Sea is interpreted on the basis of a 1-D instantaneous pure shear stretching model of the lithosphere in terms of three components: the background heat flowing out from the asthenosphere (38 mW m^–2), the transient contribution depending on the rift age and extension amount (35 mW m^–2 at the most), and the contribution due to the radiogenic elements of the lithosphere. The radiogenic component is estimated at the continental margins of the Ligurian-Provençal basin and Valencia trough, and in the surrounding mainland areas by means of available data of surface heat generation from Variscan Corsica, Maures-Estérel and the Central Massif along with a geophysical-petrological relationship between heat production and seismic velocity. The lithosphere radiogenic heat contribution q_l decreases with the thinning factor according to the exponential law: q_l() = a exp(-b), in which factor b is greater for that part of the lithosphere below the uppermost 10 km. Considering also the heat generated by radioactive isotopes in sediments, the stable Variscan lithosphere produces an average thermal flux of 30 mW m^–2 which decreases by about one half where the lithosphere is thinned by one third. Although the surface heat generation is 2·1 – 3·3 µW m^–3 in the Maures-Estérel massif — excepting small outcrops of dioritic rocks with lower heat production — and 1·8 µW m^–3 for most of Corsica, the radiogenic heating within the lithosphere for such areas is nearly the same and does not explain the higher heat flux of the Corsica margin. This asymmetric thermal pattern with surface heat flux which is 10 – 15 mW m^–2 higher than predictions is probably of upper mantle origin, or can be ascribed to penetrative magmatism.     

A simple method using climatic variables to estimate canopy temperature,sensible and latent heat fluxes in a winter wheat field on the North China Plain

Estimation of evapotranspiration from a crop field is of great importance for detecting crop water status and proper irrigation scheduling. The Penman–Monteith equation is widely viewed as the best method to estimate evapotranspiration but it requires canopy resistance, which is very difficult to determine in practice. This paper presents a simple method simplified from the Penman–Monteith equation for estimating canopy temperature (T_c). The proposed method is a biophysically‐sound extended version of that proposed by Todorovic. The estimated canopy temperature is used to calculate sensible heat flux, and then latent heat flux is calculated as the residual of the surface energy balance. An eddy covariance (EC) system and an infrared thermometer (IRT) were installed in an irrigated winter wheat field on the North China Plain in 2004 and 2005, to measure T_c, and sensible and latent heat fluxes were used to test the modified Todorovic model (MTD). The results indicate that the original Todorovic model (TD) severely underestimates T_c and sensible heat flux, and hence severely overestimates the latent heat flux. However, the MTD model has good capability for estimating T_c, and gives acceptable results for latent heat flux at both half‐hourly and daily scales. The MTD model results also agreed well with the evapotranspiration calculated from the measured T_c. Copyright © 2008 John Wiley & Sons, Ltd.     

Energy partitioning over a semi‐arid shrubland in northern China

During the last decade, the widely distributed shrublands in northern China have shown significant signs of recovery from desertification, the result of widespread conservation practices. However, to support the current efforts in conservation, more knowledge is needed on surface energy partitioning and its biophysical controls. Using eddy‐covariance measurements made over a semi‐arid shrubland in northwest China in 2012, we examined how surface energy‐balance components vary on diurnal and seasonal scales, and how biophysical factors control bulk surface parameters and energy exchange. Sensible heat flux (H) exceeded latent heat flux (λE) during most of the year, resulting in an annual Bowen ratio (β, i.e. H/λE) of 2.0. λE exceeded H only in mid‐summer when frequent rainfall co‐occurred with the seasonal peak in leaf area index (LAI). Evapotranspiration reached a daily maximum of 3.3 mm day^?1, and summed to 283 mm yr^?1. The evaporative fraction (EF, i.e. λE/R_n), Priestley–Taylor coefficient (α), surface conductance (g_s) and decoupling coefficient (Ω) were all positively correlated with soil water content (SWC) and LAI. The direct enhancement of λE by high vapour pressure deficit (VPD) was buffered by a concurrent suppression of g_s. The g_s played a direct role in controlling EF and α by mediating the effects of LAI, SWC and VPD. Our results highlight the importance of adaptive plant responses to water scarcity in regulating ecosystem energy partitioning, and suggest an important role for revegetation in the reversal of desertification in semi‐arid areas. Copyright © 2015 John Wiley & Sons, Ltd.     

Estimation of spatially distributed surface energy fluxes using remotely‐sensed data for agricultural fields

Land surface energy fluxes are required in many environmental studies, including hydrology, agronomy and meteorology. Surface energy balance models simulate microscale energy exchange processes between the ground surface and the atmospheric layer near ground level. Spatial variability of energy fluxes limits point measurements to be used for larger areas. Remote sensing provides the basis for spatial mapping of energy fluxes. Remote‐sensing‐based surface energy flux‐mapping was conducted using seven Landsat images from 1997 to 2002 at four contiguous crop fields located in Polk County, northwestern Minnesota. Spatially distributed surface energy fluxes were estimated and mapped at 30 m pixel level from Landsat Thematic Mapper and Enhanced Thematic Mapper images and weather information. Net radiation was determined using the surface energy balance algorithm for land (SEBAL) procedure. Applying the two‐source energy balance (TSEB) model, the surface temperature and the latent and sensible heat fluxes were partitioned into vegetation and soil components and estimated at the pixel level. Yield data for wheat and soybean from 1997 to 2002 were mapped and compared with latent heat (evapotranspiration) for four of the fields at pixel level. The spatial distribution and the relation of latent heat flux and Bowen ratio (ratio of sensible heat to latent heat) to crop yield were studied. The root‐mean‐square error and the mean absolute percentage of error between the observed and predicted energy fluxes were between 7 and 22 W m⁻² and 12 and 24% respectively. Results show that latent heat flux and Bowen ratio were correlated (positive and negative) to the yield data. Wheat and soybean yields were predicted using latent heat flux with mean R² = 0·67 and 0·70 respectively, average residual means of −4·2 bushels/acre and 0·11 bushels/acre respectively, and average residual standard deviations of 16·2 bushels/acre and 16·6 bushels/acre respectively (1 bushel/acre ≈ 0·087 m³ ha⁻¹). The flux estimation procedure from the SEBAL‐TSEB model was useful and applicable to agricultural fields. Copyright © 2005 John Wiley & Sons, Ltd.     

Energy balance and evaporation loss of an agricultural reservoir in a semi‐arid climate (south‐eastern Spain)

A typical agricultural water reservoir (AWR) of 2400 m² area and 5 m depth, located in a semi‐arid area (southern Spain), was surveyed on a daily basis for 1 year. The annual evaporation flux was 102·7 W m^?2, equivalent to an evaporated water depth of 1310 mm year^?1. The heat storage rate G exhibited a clear annual cycle with a peak gain in April (G ～ 45 W m^?2) and a peak loss in November (G ～ 40 W m^?2), leading to a marked annual hysteretic trend when evaporation (λE) was related to net radiation (R_n). λE was strongly correlated with the available energy A, representing 91% of the annual AWR energy loss. The sensible heat flux H accounted for the remaining 9%, leading to an annual Bowen ratio in the order of 0·10. The equilibrium and advective evaporation terms of the Penman formula represented 76 and 24%, respectively, of the total evaporation, corresponding to a annual value of the Priestley–Taylor (P–T) coefficient (α) of 1·32. The P–T coefficient presented a clear seasonal pattern, with a minimum of 1·23 (July) and a maximum of 1·65 (December), indicating that, during periods of limited available energy, AWR evaporation increased above the potential evaporation as a result of the advection process. Overall, the results stressed that accurate prediction of monthly evaporation by means of the P–T formula requires accounting for both the annual cycle of storage and the advective component. Some alternative approaches to estimating R_n, G and α are proposed and discussed. Copyright © 2009 John Wiley & Sons, Ltd.     

Variability of tree transpiration across three zones in a southeastern U.S. Piedmont watershed

Quantifying the spatial variability of species-specific tree transpiration across hillslopes is important for estimating watershed-scale evapotranspiration (ET) and predicting spatial drought effects on vegetation. The objectives of this study are to (1) assess sap flux density (J_s) and tree-level transpiration (T_s) across three contrasting zones a (riparian buffer, mid-hillslope and upland-hillslope, (2) determine how species-specific J_s responds to vapour pressure deficit (VPD) and (3) estimate watershed-level transpiration (T_w) using T_s derived from each zone. During 2015 and 2016, we measured J_s in eight tree species in the three topographic zones in a small 12-ha forested watershed in the Piedmont region of central North Carolina. In the dry year of 2015, loblolly pine (Pinus taeda), Virginia pine (Pinus virginiana) and sweetgum (Liquidambar styraciflua) J_s rates were significantly higher in the riparian buffer when compared to the other two zones. In contrast, J_s rates in tulip poplar (Liriodendron tulipifera) and red maple (Acer rubrum) were significantly lower in the buffer than in the mid-hillslope. Daily T_s varied by zone and ranged from 10 to 93 L/day in the dry year and from 9 to 122 L/day in the wet year (2016). J_s responded nonlinearly to VPD in all species and zones. Annual T_w was 447, 377 and 340 mm based on scaled-J_s data for the buffer, mid-hillslope and upland-hillslope, respectively. We conclude that large spatial variability in J_s and scaled T_w was driven by differences in soil moisture at each zone and forest composition. Consequently, spatial heterogeneity of vegetation and soil moisture must be considered when accurately quantifying watershed level ET.     

藏南羊卓雍错湖面大气湍流特征观测分析

湍流运动是大气边界层的本质特征,是地表与大气之间能量和物质交换的主要方式.本文利用2016和2017年4-10月藏南羊卓雍错湖泊涡动观测资料,分析了湖面大气湍流方差和湍流特征量的统计和变化特征.结果表明：（1）不稳定层结下,三维风速分量和超声虚温、水汽密度、CO₂密度的无量纲标准差随稳定度变化符合Monin-Obukhov相似理论的"1/3"或"-1/3"次幂律,垂直风速的拟合效果最好;稳定层结下,除CO₂密度无量纲标准差与稳定度无明显关系外,其他量基本上满足相似性规律;中性条件下,以上物理量的无量纲标准差分别趋近常数：3.57、3.93、0.77、20.91、6.35和11.96.（2）水平方向平均湍流强度（0.60和0.58）大于垂直方向（0.13）,三维方向湍流强度与平均风速的变化呈显著负相关,相关系数分别为-0.39、-0.42和-0.34.（3）湖面湍流动能随风速呈线性增长,增长率达0.45 m/s;近中性层结时湍流动能最大,层结越稳定或不稳定湍流动能均减小.（4）湖泊下午到傍晚动量输送较强,13：00-22：30时间段平均动量通量达0.091 kg/（m·s²）;热量输送以潜热为主,潜热通量日平均值（77.3 W/m²）是感热通量（14.6 W/m²）的5.3倍,感热和潜热通量日变化峰值分别出现在5：30（22.4 W/m²）和16：00（106.6 W/m²）.     

Short‐term stream water temperature observations permit rapid assessment of potential climate change impacts

Assessment of potential climate change impacts on stream water temperature (T_s) across large scales remains challenging for resource managers because energy exchange processes between the atmosphere and the stream environment are complex and uncertain, and few long‐term datasets are available to evaluate changes over time. In this study, we demonstrate how simple monthly linear regression models based on short‐term historical T_s observations and readily available interpolated air temperature (T_a) estimates can be used for rapid assessment of historical and future changes in T_s. Models were developed for 61 sites in the southeastern USA using ≥18 months of observations and were validated at sites with longer periods of record. The T_s models were then used to estimate temporal changes in T_s at each site using both historical estimates and future T_a projections. Results suggested that the linear regression models adequately explained the variability in T_s across sites, and the relationships between T_s and T_a remained consistent over 37 years. We estimated that most sites had increases in historical annual mean T_s between 1961 and 2010 (mean of +0.11 °C decade^?1). All 61 sites were projected to experience increases in T_s from 2011 to 2060 under the three climate projections evaluated (mean of +0.41 °C decade^?1). Several of the sites with the largest historical and future T_s changes were located in ecoregions home to temperature‐sensitive fish species. This methodology can be used by resource managers for rapid assessment of potential climate change impacts on stream water temperature. Copyright © 2014 John Wiley & Sons, Ltd.     

It goes both ways: measurements of simultaneous evapotranspiration and fog droplet deposition at a montane cloud forest

Fluxes of latent heat, sensible heat, and water vapor, including turbulent deposition of fog droplets, were measured for two months in autumn 2005 within a subtropical montane cypress forest in Taiwan. The goal of the study was to determine whether significant evapotranspiration can occur during foggy conditions. Water vapor fluxes, Q_W, as determined with the Bowen Ratio method, were compared to those simultaneously measured with the eddy covariance method. The median Bowen Ratio was 1.06, and the median Q_W flux was 5 · 2 × 10^?5 kg m^?2 s^?1. The vertical gradients of temperature and specific humidity over the forest, ΔT and Δq, peaked around noon during days without fog, and were reduced during foggy conditions. For 66% of the data points, ΔT and Δq were negative, corresponding to positive (upward) fluxes of sensible heat Q_H and latent heat Q_E. A Monte Carlo simulation proved that statistically significant evapotranspiration rates, i.e., upward water vapor fluxes, occurred during fog. At the same time, deposition fluxes of fog droplets occurred. Our results show that even during fog events, significant evapotranspiration may occur. Copyright © 2008 John Wiley & Sons, Ltd.     

Mechanics of interrill erosion with wind‐driven rain

The vector physics of wind‐driven rain (WDR) differs from that of wind‐free rain, and the interrill soil detachment equations in the Water Erosion Prediction Project (WEPP) model were not originally developed to deal with this phenomenon. This article provides an evaluation of the performance of the interrill component of the WEPP model for WDR events. The interrill delivery rates were measured in the wind tunnel facility of the International Center for Eremology (ICE), Ghent University, Belgium with an experimental setup to study different raindrop impact velocity vectors. Synchronized wind and rain simulations with wind velocities of 6, 10 and 14 m s^–1 were applied to a test surface placed on windward and leeward slopes of 7, 15 and 20%. Since both rainfall intensity and raindrop impact velocity varied greatly depending on differences in the horizontal wind velocity under WDRs, the resultant kinetic energy flux (KE_r, in J m^–2 s^–1) was initially used in place of the WEPP model intensity term in order to incorporate the effect of wind on impact velocity and frequency of raindrops. However, our results showed only minor improvement in the model predictions. For all research data, the model Coefficients of Determination (r²) were 0·63 and 0·71, when using the WEPP and the KE_r approaches, respectively. Alternately, integrating the angle of rain incidence into the model by vectorally partitioning normal kinetic energy flux (KE_rn, in J m^–2 s^–1) from the KE_r greatly improved the model's ability to estimate the interrill sediment delivery rates (r² = 0·91). This finding suggested that along with the fall trajectory of wind‐driven raindrops with a given frequency, raindrop velocity and direction at the point of impact onto the soil surface provided sufficient physical information to improve WEPP sediment delivery rate predictions under WDR. Copyright © 2012 John Wiley & Sons, Ltd.