Relationships between riparian evapotranspiration and groundwater depth along a semiarid irrigated river valley

Evapotranspiration (ET) from riparian vegetation can be difficult to estimate due to relatively abundant water supply, spatial vegetation heterogeneity, and interactions with anthropogenic influences such as shallower groundwater tables, increased salinity, and nonpoint source pollution induced by irrigation. In semiarid south-eastern Colorado, reliable ET estimates are scarce for the riparian corridor that borders the Arkansas River. This work investigates relationships between the riparian ecosystem along the Arkansas River and an underlying alluvial aquifer using ET estimates from remotely sensed data and modelled water table depths. Results from a calibrated, finite-difference groundwater model are used to estimate weekly water table fluctuations in the riparian ecosystem from 1999 to 2009, and estimates of ET are calculated using the Operational Simplified Surface Energy Balance (SSEBop) model with over 200 Landsat scenes covering over 30 km² of riparian ecosystem along a 70-km stretch of the river. Comparison of calculated monthly SSEBop ET to estimated alfalfa reference ET from local micrometeorological station data indicated statistically significant high linear correspondence (R² = .87). Daily calculated SSEBop ET showed statistically significant moderate linear correspondence with data from a local weighing lysimeter (R² = .59). Simulated monthly SSEBop ET values were larger in drier years compared with wetter years, and ET variability was also larger in drier years. Peak ET most commonly occurred during the month of June for all 11 years of analysis. Relationships between ET and water table depth showed that peak monthly ET was highest when groundwater depths were less than about 3 m, and ET values were significantly lower for groundwater depths greater than 3 m. Negative sample Spearman correlation highlighted riparian corridor locations where ET increased as a result of decreased groundwater depths across years with different hydroclimatic conditions. This study shows how a combination of remotely sensed riparian ET estimates and a regional groundwater model can improve our understanding of linkages between riparian consumptive use and near-river groundwater conditions influenced by irrigation return flow and different climatic drivers.     

砌体结构抗震抗剪强度分析

进行了墙体抗剪强度理论和公式的研究。将主拉强度理论与剪摩强理论相结合,形成了一种新和强度理论,即拉摩强度理论。根据该理论,导出了墙体抗剪强度公式。给出了公式系数的确定方法,并确定了砖墙体的公式系数。将本文建议的砖墙体抗剪强度公式与抗震和砌体规范采用的抗剪强度公式进行了比较,比较结果表明：抗震规范公式值较小;砌体规范公式值较在,本文公式值介于两者中间。     

On the performance of surface renewal analysis to estimate sensible heat flux over two growing rice fields under the influence of regional advection

High-frequency temperature data were recorded at one height and they were used in Surface Renewal (SR) analysis to estimate sensible heat flux during the full growing season of two rice fields located north–northeast of Colusa, CA (in the Sacramento Valley). One of the fields was seeded into a flooded paddy and the other was drill seeded before flooding. To minimize fetch requirements, the measurement height was selected to be close to the maximum expected canopy height. The roughness sub-layer depth was estimated to discriminate if the temperature data came from the inertial or roughness sub-layer. The equation to estimate the roughness sub-layer depth was derived by combining simple mixing-length theory, mixing-layer analogy, equations to account for stable atmospheric surface layer conditions, and semi-empirical canopy–architecture relationships. The potential for SR analysis as a method that operates in the full surface boundary layer was tested using data collected over growing vegetation at a site influenced by regional advection of sensible heat flux. The inputs used to estimate the sensible heat fluxes included air temperature sampled at 10 Hz, the mean and variance of the horizontal wind speed, the canopy height, and the plant area index for a given intermediate height of the canopy. Regardless of the stability conditions and measurement height above the canopy, sensible heat flux estimates using SR analysis gave results that were similar to those measured with the eddy covariance method. Under unstable cases, it was shown that the performance was sensitive to estimation of the roughness sub-layer depth. However, an expression was provided to select the crucial scale required for its estimation.     

Atmospheric stability effect on subgrid-scale physics for large-eddy simulation

Field measurements in the atmospheric boundary layer were carried out to identify the effect of atmospheric stability on subgrid-scale physics for large-eddy simulation. The basic instrumentation setup consisted of 12 three-dimensional sonic anemometers arranged in two parallel horizontal arrays (seven sensors in the lower array and five sensors in the upper array). Data from this setup are used to compute the subgrid-scale (SGS) heat fluxes and SGS dissipation of the temperature variance under stable and unstable stability conditions. The relative contribution of the SGS vertical flux to the total turbulent flux increases when going from unstable to stable conditions. The relative importance of negative SGS dissipation (backscatter) events becomes larger under stable conditions. The model coefficients for two well-known SGS models (eddy-viscosity and non-linear) are computed. Model coefficients are found to depend strongly on stability. Under both stable and unstable conditions, large negative SGS dissipation is associated with the onset of ejection events while large positive SGS dissipation tends to occur during the onset of sweep events. These findings are also supported by conditionally sampled 2D velocity and temperature fields obtained using the 12 anemometers placed in a vertical array.     

A method to correct eddy covariance flux underestimates under an advective environment for arid or semi-arid regions

Water scarcity is one of the main factors limiting agricultural development in arid or semi-arid areas. Accurate Evapotranspiration (ET) observations and estimations are crucial in water cycle studies to estimate water losses from the terrestrial surfaces to the atmosphere to close the regional water budget. The eddy covariance (EC) method is an important technique measure ET and other land surface energy fluxes. However, the underestimation of energy fluxes and the problem of EC energy balance non-closure are far from solved. In this study, a new method is proposed to account for advection in order to correct EC data under advective environments. This advection based method was applied to data from Bushland, TX, which is subject to dry air and strong winds. Observations from two identical EC systems as well as two precision monolithic weighing lysimeters were used in this analysis. Both EC sites showed significant underestimates of evapotranspiration (ET) compared with lysimeter measurements. The daily energy balance closure for NE01 and SE02 sites were 0.78 and 0.74 respectively. The advection correction method provided improved performance in daytime, and it is more suitable for ET estimate than forcing closure under the advective environment. For nighttime, two methods (NCM1 and NCM2) were proposed to correct EC underestimates. Finally, all the corrected ET values were compared with the lysimeter measurements. For NE01 site, the MAD (mean absolute deviation) and the RMSD (root mean square deviation) were 47.72 W/m² and 67.66 W/m², respectively; and the r² (coefficient of determination) was 0.85. For SE02 site, the MAD and RMSD were 30.59 W/m² and 44.43 W/m²; and the r² was 0.93. The statistical measures illustrated that the proposed methods are functional and appropriate under an advective environment. The accurate estimate of actual evapotranspiration will benefit both the strategic planning of optimal water uses and the improved understanding the environmental and hydrological processes.     

Multiscale modeling of the atmospheric environment over a forest canopy

Vegetation constitutes one of the fundamental types of land use on Earth. The presence of trees in urban areas can decrease local winds and exchange sensible and latent heat with the surrounding environments, thus exerting notable influences on the urban microenvironment. A better understanding of the turbulent transfer of momentum and scalars around vegetation canopy could significantly contribute to improvement of the urban environment. This work develops a large-eddy simulation(LES) method that is applicable to model the flow and scalar transport over the forest canopy. We study the atmospheric flow over complex forested areas under typical weather conditions by coupling LES to the mesoscale model. Models of radiation and energy balance have been developed with explicit treatment of the vegetation canopy. By examining the flow over a forest canopy under a range of stability conditions, we found that buoyancy enhances or suppresses turbulent mixing in unstable or stable atmosphere respectively, with decreasing or increasing wind shear, respectively. From the multiscale modeling of the Beijing Olympic Forest Park, the present coupling scheme proves to better resolve the diurnal variations in wind speed,temperature, and relative humidity over complex urban terrains. The coupling scheme is superior to the traditional mesoscale model in terms of wind field simulation. This is mainly because the coupling scheme not only takes the influences of external mesoscale flow into consideration, but also resolves the heterogeneous urban surface at a fine scale by downscaling, thus better reproducing the complex flow and turbulent transport in the urban roughness sublayer.     

Interception loss from eucalypt forest: Lysimeter determination of hourly rates for long term evaluation

Analyses of the response by a weighing lysimeter in Kioloa State Forest during and after rainfall provided values of interception loss rate. The derived rates for time scales between 0.1 and 1.0mm h^?1 were generally similar throughout storm events to losses determined from throughfall and stemflow observations. During post-rainfall periods of canopy drying, enhanced rates of lysimeter evaporation were consistent with micrometeorological determinations of the partitioning of available radiant energy, based on atmospheric gradients of humidity and temperature. Interception losses from the eucalypt forest, deduced from the lysimeter response, varied between 10 and 15 per cent of gross rainfall in three consecutive 12 month periods whereas the corresponding rainfall ranged between 590 and 1530 mm yr^?1. Daytime losses accounted for about two-thirds of total interception loss with a similar fraction occurring during rain periods. Storage capacity of the evergreen forest canopy was inferred to be 0.35 mm. Hourly loss rates during rainfall ranged up to 0.8 mm h^?1 but with decreasing mean values and variability with increasing time scale resulting in a monthly mean value computed for the total number of hours of rain of approximately 0.1 mm h^?1. A preliminary analysis of loss rate in terms of storm windspeed and rainfall intensity explained about half of its variation in statistically derived relationships. Improved time resolution of the order of seconds was considered a prerequisite to the physical understanding of turbulent transport from saturated canopies. The small value of interception storage capacity was considered in relation to that for pine forest as a basis for explaining observed differences in interception behaviour between eucalypt forest and coniferous plantations in the same area. Large differences in interception losses between the Kioloa site and evergreen forest in the South Island of New Zealand and also eucalypt forest in Western Australia were attributed to dissimilar meteorological conditions at the various sites.     

How representative are local measurements of the surface shear stress for regional values?

Local values of surface shear stress measured by eddy correlation devices were compared to regional ones generated by Monin–Obukhov similarity (MOS) theory in the atmospheric surface layer over the cooperative atmosphere–surface exchange study (CASES) area in eastern Kansas. For this comparison, minisodar average wind speed profiles were used to implement MOS. The local surface shear stress values were found to be representative of the regional ones only under free convective conditions. For weak to strongly neutral atmospheric conditions, regional values of the surface shear stress were underpredicted by a factor of two by the locally measured ones. Empirical interpolation equations are proposed to describe the relationship in the intermediate range. The surface roughness values of 0.12 and 0.40 m, respectively, were determined for the two locations of the minisodar sounders by regional analysis of their mean wind speed measurements under neutral atmospheric conditions. The roughness values found in this present study compared well to past investigations conducted in this region. Minisodar roughness values were found in this present study to be a good alternative to radiosonde determined ones.     

Doppler radar spectral width broadening due to beamwidth and wind shear

The spectral width observed by Doppler radars can be due to several effects including the atmospheric turbulence within the radar sample volume plus effects associated with the background flow and the radar geometry and configuration. This study re-examines simple models for the effects due to finite beam-width and vertical shear of the horizontal wind. Analytic solutions of 1- and 2-dimensional models are presented. Comparisons of the simple 2-dimensional model with numerical integrations of a 3-dimensional model with a symmetrical Gaussian beam show that the 2-dimensional model is usually adequate. The solution of the 2-dimensional model gives a formula that can be applied easily to large data sets. Analysis of the analytic solutions of the 2-dimensional model for off-vertical beams reveals a term that has not been included in mathematical formulas for spectral broadening in the past. This term arises from the simultaneous effects of the changing geometry due to curvature within a finite beamwidth and the vertical wind shear. The magnitude of this effect can be comparable to that of the well-known effects of beam-broadening and wind shear, and since it can have either algebraic sign, it can significantly reduce (or increase) the expected spectral broadening, although under typical conditions it is smaller than the beam-broadening effect. The predictions of this simple model are found to be consistent with observations from the VHP radar at White Sands Missile Range, NM.     

Quantifying net ecosystem exchange by multilevel ecophysiological and turbulent transport models

To quantify the interplay between scalar sources and sinks (S_c) and net ecosystem exchange (NEE), “forward” and “inverse” approaches have been proposed. The canonical form of forward approaches is a one-dimensional ecophysiological-radiative transfer scheme coupled to turbulent transport theory. In contrast, inverse approaches strictly rely on turbulent transport theory and mean scalar concentration as their primary input to infer S_c and NEE. While the formulation of both approaches have evolved over the past decade, no systematic comparison between them was undertaken for the same data set, and over a wide range of atmospheric conditions. Our objective is to compare the predicted S_c and NEE from these two approaches with eddy-covariance measurements. The results show that the forward method outperformed all three inverse methods for unstable and neutral conditions on short time scales (30 min) but yielded comparable results at longer time scales. Poor agreement was obtained under stable conditions for all models. Hence, for modeling event-based flux variations, forward models are preferred. Since the forward method requires detailed knowledge of ecophysiological, drag, radiative transfer and other canopy attributes, all of which are difficult to obtain on a routine basis, a symbiotic use of forward and inverse approaches is most advantageous.     

The current-voltage relationship revisited: exact and approximate formulas with almost general validity for hot magnetospheric electrons for bi-Maxwellian and kappa distributions

We derive the current-voltage relationship in the auroral region taking into account magnetospheric electrons for the bi-Maxwellian and kappa source plasma distribution functions. The current-voltage formulas have in principle been well known for a long time, but the kappa energy flux formulas have not appeared in the literature before. We give a unified treatment of the bi-Maxwellian and kappa distributions, correcting some errors in previous work. We give both exact results and two kinds of approximate formulas for the current density and the energy flux. The first approximation is almost generally valid and is practical to compute. The first approximation formulas are therefore suitable for use in simulations. In the second approximation we assume in addition that the thermal energy is small compared to the potential drop. This yields even simpler linear formulas which are suitable for many types of event studies and which have a more transparent physical interpretation than the first approximation formulas. We also show how it is possible to derive the first approximation formulas even for those distributions for which the exact results can not be computed analytically. The kappa field-aligned conductance value turns out always to be smaller than the corresponding Maxwellian conductance. We also verify that the obtained kappa current density and energy flux formulas go to Maxwellian results when k .     

重力场向下延拓Milne法

重力场向下延拓能够突出局部和浅部的异常信息,分离叠加的异常特征.但是向下延拓通常具有过程不稳定、下延深度小、结果不准确等问题.针对向下延拓所存在的不足,本文利用重力场及其垂向一阶导数,基于辛普森(Simpson)求积公式,推导出重力场向下延拓米尔尼(Milne)公式.将本文向下延拓方法应用于模型数据,向下延拓模型结果及误差曲线表明,相对于向下延拓快速傅里叶变换(FFT)法和积分迭代法,向下延拓Milne法的深度更大,相对误差更小;相对模型值,向下延拓Milne法能够获得稳定且准确的结果.对加拿大乃查科(Nechako)盆地地区实测航空重力数据进行本文方法向下延拓验证,处理结果表明,相对于实测异常,本文方法向下延拓结果能够很好还原实测数据,并且在进一步向下延拓中反映原始异常的趋势,增强局部和细小异常信息.     

Hydrologic response of a Hawaiian watershed to future climate change scenarios

The impact of potential future climate change scenarios on streamflow and evapotranspiration (ET) in a mountainous Hawaii watershed was studied using the distributed hydrology soil vegetation model (DHSVM). The hydrologic response of the watershed was simulated for 43 years for different levels of atmospheric CO₂ (330, 550, 710 and 970 ppm), temperature (+1.1 and + 6.4 °C) and precipitation (±5%, ±10% and ±20%) on the basis of the Intergovernmental Panel on Climate Change (IPCC) AR4 projections under current, B1, A1B1 and A1F1 emission scenarios. Vegetation leaf conductance and leaf area index were modified to reflect the increase in CO₂ concentration. The relative departure of streamflow and ET from their levels during the reference scenarios was calculated on a monthly and annual basis. Results of this study indicate that the streamflow and ET are less sensitive to changes in temperature compared with changes in precipitation. However, temperature increase coupled with precipitation showed significant effect on ET and streamflow. Changes in leaf conductance and leaf area index with increasing CO₂ concentration under A1F1 scenario had a significant effect on ET and subsequently on streamflow. Evapotranspiration is less sensitive than streamflow for a similar level of change in precipitation. On the basis of a range of climate change scenarios, DHSVM predicted a change in ET by ±10% and streamflow between ?51% and 90%. From the six ensemble mean scenarios for AR4 A1B, simulations suggest reduction in streamflow by 6.7% to 17.2%. These reductions would produce severe impact on water availability in the region. Copyright © 2011 John Wiley & Sons, Ltd.     

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

湍流运动是大气边界层的本质特征,是地表与大气之间能量和物质交换的主要方式.本文利用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²）.     

Micrometeorological measurements and vapour pressure deficit relations under in-field rainwater harvesting

In a cropped field, microclimate and thermal stability conditions depend on the canopy structures and the prevailing weather. The main aim of the study therefore was to characterize the vertical profiles of weather variables within and above a maize (Zea mays L.) canopy and to describe the water vapour pressure deficit (VPD) under different atmospheric and soil surface conditions for both wide and narrow runoff strips with the in-field rainwater harvesting (IRWH) system. Micrometeorological measurements of wind, temperature and relative humidity were performed at eight levels, within canopy (1.8 and 2.1 m), and just above the canopy (2.4, 2.7, 3.0, and 3.3 m) up to reference levels (3.9 and 4.5 m) when the maize reached a maximum height of 2.2 m. Under incomplete canopy cover of the IRWH system, two important factors complicated evapotranspiration estimation, namely the local advection and high temperatures of the bare soil between adjacent plant rows. Diurnal variations of water vapour related to turbulence at each locality and its position in the thermal internal boundary layers. Generally, advection was more pronounced in wide runoff strips than narrow strips. On wide runoff strips the wind was more effective in replacing the air between the rows and maintained a higher driving force for evaporation. The maximum VPD over the narrow strips was observed at reference level during a dry day, at about 2.2 kPa in the afternoon, while wet day VPD reached a maximum of 1.8 kPa. The VPD of the wide runoff strips correlated negatively with wind speed, but showed a fairly positive correlation with some scattered values on wet days after rain. Therefore, profile characteristics within and above plant canopies played a key role in determining the VPD and consequently, could help to explain transpiration rates of crops. Hence, VPD relations enhanced the understanding of the heat energy exchange processes under the heterogeneous nature of maize canopy of the IRWH tillage system.     

Hydro-gravimetry in West-Africa: First results from the Djougou (Benin) superconducting gravimeter

The increasing number of hydro-gravimetry studies proves the rising interest of the hydrology community toward this monitoring method. The accuracy of superconducting gravimeters (SG) potentially allows the retrieval of small water storage changes (WSC) down to a few millimeters of equivalent water thickness. However, the importance of corrections applied to SG data to achieve such a precision in gravity residuals should be recalled. The Djougou permanent gravity station presented in this paper and located in northern Benin, West-Africa, provides a good opportunity to review these considerations. This station is equipped since July 2010 with the superconducting gravimeter SG-060 aimed at deriving WSC at different time-scales, daily to inter-annual. In this area, WSC are (1) part of the control system for evapotranspiration (ET) process, a key variable of the West-African monsoon cycle and (2) the state variable for resource management, a critical issue in storage-poor hard rock basement contexts such as in northern Benin. The potential for deriving WSC from time-lapse gravity data partly depends on environmental features such as topography and the instrument shelter. Therefore, this issue is addressed first, with the background idea that such sensitivity analysis should be undertaken before setting up any new instrument. In Djougou, local topography is quite flat leading to a theoretical straightforward relationship between gravity changes and WSC, close to the standard Bouguer value. However, the shelter plays a significant masking role, which is the principal limitation to the retrieval of fast hydrological processes such as ET following a rain event. Several issues concerning classical gravity corrections are also addressed in the paper. These include gap-filling procedures during rain-events and drift estimates for short time series. Special attention is provided to atmospheric corrections, and different approaches are tested: a simple scalar admittance, a filtered scalar admittance, a frequency-dependent admittance and direct atmospheric loading calculations. It is shown that the physically based approach of direct loading calculations performs better in both residual minimization and ET retrieval. Moreover, non-local hydrological effects are investigated and account for about 20% of the gravity residuals. Finally, gravity residuals are briefly analyzed at two distinct time scales: rapid (up to a few days) and seasonal. At the rapid time-scale, it is shown that ET retrieval is hardly achievable given shelter size and state-of-the-art atmospheric corrections. Still, mean values retrieved from this study are in accordance with known values of potential ET and lateral flow. Direct comparison of gravity changes with hydrological data (neutron probe monitoring and water table levels) show some discrepancies, particularly for the hydrological year of 2011, for which all hydrological data show a deficit, but SG and FG5 data do not. This preliminary analysis both provides a basis and call for further hydro-gravity modeling, to comprehensively investigate the water-cycle at the Djougou station.     

Lysimetric evaluation of SEBAL using high resolution airborne imagery from BEAREX08

In this study, Surface Energy Balance Algorithm for Land (SEBAL) was evaluated for its ability to derive aerodynamic components and surface energy fluxes from very high resolution airborne remote sensing data acquired during the Bushland Evapotranspiration and Agricultural Remote Sensing Experiment 2008 (BEAREX08) in Texas, USA. Issues related to hot and cold pixel selection and the underlying assumptions of difference between air and surface temperature (dT) being linearly related to the surface temperature were also addressed. Estimated instantaneous evapotranspiration (ET) and other components of the surface energy balance were compared with measured data from four large precision weighing lysimeter fields, two each managed under irrigation and dryland conditions. Instantaneous ET was estimated with overall mean bias error and root mean square error (RMSE) of 0.13 and 0.15 mm h⁻¹ (23.8 and 28.2%) respectively, where relatively large RMSE was contributed by dryland field. Sensitivity analysis of the hot and cold pixel selection indicated that up to 20% of the variability in ET estimates could be attributed to differences in the surface energy balance and roughness properties of the anchor pixels. Adoption of an excess resistance to heat transfer parameter model into SEBAL significantly improved the instantaneous ET estimates.     

Photosynthesis-transpiration coupling model at canopy scale in terrestrial ecosystem

At the hypothesis of big leaf, an ecosystem photosynthesis-transpiration coupling cycle model was established by the scaled SMPT-SB model from single leaf to canopy, and model parameterization methods were discussed. Through simulating the canopy light distribution, canopy internal conductance to CO2 can be scaled from single leaf to canopy by integrating to canopy using the relationship between single internal conductance and photosynthetic photon flux density. Using the data observed by eddy covariance method from the Changbai Mountains site of ChinaFLUX, the application of the model at the canopy scale was examined. Under no water stress, the simulated net ecosystem photosynthesis rate fitted with the observed data very well, the slope and R2 of the line regression equation of the observed and simulated values were 0.7977 and 0.8892, respectively (n = 752), and average absolute error was 3.78 μmol CO2 m-2s-1; the slope, R2 and average absolute error of transpiration rate were 0.7314, 0.4355 and 1.60mmol H2O m-2 s-1, respectively (n = 752). The relationship between canopy photosynthesis,transpiration and external environmental conditions was discussed by treating the canopy as a whole and neglecting the comprehensive feedback mechanism within canopy, and it was noted that the precipitation course affected the transpiration rate simulation badly. Compared to the models based on eco-physiological processes, the SMPT-SB model was simple and easy to be used. And it can be used as a basic carbon and water coupling model of soil-plant-atmosphere continuum.     

Responses of canopy transpiration and canopy conductance of peach (Prunus persica) trees to alternate partial root zone drip irrigation

We investigated canopy transpiration and canopy conductance of peach trees under three irrigation patterns: fixed 1/2 partial root zone drip irrigation (FPRDI), alternate 1/2 partial root zone drip irrigation (APRDI) and full root zone drip irrigation (FDI). Canopy transpiration was measured using heat pulse sensors, and canopy conductance was calculated using the Jarvis model and the inversion of the Penman–Monteith equation. Results showed that the transpiration rate and canopy conductance in FPRDI and APRDI were smaller than those in FDI. More significantly, the total irrigation amount was greatly reduced, by 34·7% and 39·6%, respectively for APRDI and FPRDI in the PRDI (partial root zone drip irrigation) treatment period. The daily transpiration was linearly related to the reference evapotranspiration in the three treatments, but daily transpiration of FDI is more than that of APRDI and FPRDI under the same evaporation demand, suggesting a restriction of transpiration water loss in the APRDI and FPRDI trees. FDI needed a higher soil water content to carry the same amount of transpiration as the APRDI and FPRDI trees, suggesting the hydraulic conductance of roots of APRDI and FPRDI trees was enhanced, and the roots had a greater water uptake than in FDI when the average soil water content in the root zone was the same. By a comparison between the transpiration rates predicted by the Penman–Monteith equation and the measured canopy transpiration rates for 60 days during the experimental period, an excellent correlation along the 1:1 line was found for all the treatments (R² > 0·80), proving the reliability of the methodology. Copyright © 2005 John Wiley & Sons, Ltd.     

Canopy interception of apple orchards should not be ignored when assessing evapotranspiration partitioning on the Loess Plateau in China

Canopy interception is one of the most important processes in an ecosystem, but it is still neglected when assessing evapotranspiration (ET) partitioning in apple orchards on the Loess Plateau in China. To explore the importance of canopy interception, we monitored two neighbouring apple orchards on the Loess Plateau in China, one 8‐year‐old and the other 18‐years old at the start of the study, from May to September for four consecutive years (2013–2016). We measured parameters of canopy interception (I) including precipitation, throughfall, stemflow, leaf area index, transpiration (T), and soil evaporation (S) to quantify ET. The importance of canopy interception was then assessed by comparing the relationship between water supply (precipitation) and water demand (ET), calculated with and without considering canopy interception (T + S and T + S + I, respectively). Tree age clearly influenced canopy interception, as estimates of annual canopy interception during the study years in the younger and older orchards amounted to 22.2–29.4 mm and 26.8–39.9 mm, respectively. Daily incident rainfall and rainfall intensity in both orchards were significantly positively correlated with daily canopy interception in each year. The relationship between annual precipitation and annual ET (calculated with and without consideration of canopy interception) in the younger orchard differed during 2015 and 2016. Ignoring canopy interception would result in underestimation of annual ET in both apple orchards and hence incorrect evaluation of the relationship between water supply and water demand, particularly for the younger orchard during 2015 and 2016. Thus, for a complete understanding of water consumption in apple orchards in this and similar regions, canopy interception should not be ignored when assessing ET partitioning.