Evaluation of evapotranspiration models over semi‐arid and semi‐humid areas of China

Evapotranspiration (ET) is a critical component in the hydrological cycle. However, its actual values appear to be difficult to obtain, especially in areas in which precipitation has high inter‐annual variability. Here, we evaluated eight commonly used ET models in semi‐arid and semi‐humid areas of China. The order of overall performance from best to worst is as follows: the revised Priestley–Taylor model (PT‐JPL, 0.71, 1.65 [18.37%], 4.72 [49.19%]) a a Statistics (model abbreviation, coefficient of determination, bias [relative value], standard deviation [relative value]).
, the modified PT‐JPL model (M1‐PT‐JPL, 0.67, ?0.68 [7.56%], 3.87 [40.31%]), the Community Land Model (CLM, 0.68, ?2.52 [28.01%], 5.10 [53.17%]), the modified PT‐JPL model (M2‐PT‐JPL, 0.63, 0.57 [6.27%], 5.04 [52.52%]), the revised Penman–Monteith model (RS‐PM, 0.62, 3.56 [37.40%], 6.11 [63.68%]), an empirical model (Wang, 0.59, ?1.04 [11.57%], 5.61 [58.43%]), the advection‐aridity model (AA, 0.55, 5.56 [61.78%], 7.45 [77.60%]), and the energy balance model (SEBS, 0.35, 5.11 [56.72%], 9.43 [98.18%]). The performance of all of the models is comparably poor in winter and summer, except for the PT‐JPL model, and relatively good in spring and autumn. Because of the vegetation control on ET, the Wang, RS‐PM, PT‐JPL, M1‐PT‐JPL, and M2‐PT‐JPL models perform better for cropland, whereas the AA model, SEBS model and CLM perform better for grassland. The CLM, PT‐JPL, and Wang models perform better in semi‐arid region than in semi‐humid region, whereas the opposite is true for SEBS and RS‐PM. The AA, M1‐PT‐JPL, and M2‐PT‐JPL models perform similarly in semi‐arid and semi‐humid regions. When considering the inter‐annual variability in precipitation, the Wang model has relatively good performance under only some annual precipitation conditions; the performance of the PT‐JPL and AA models is reduced under conditions of high precipitation; the two modified PT‐JPL models inherited the steady performance of the PT‐JPL model and improved the performance under conditions of high annual precipitation by the modification of the soil moisture constraint. RS‐PM is more appropriate for humid conditions. CLM and PT‐JPL models could be effectively applied to all precipitation conditions because of their good performance across a wide annual precipitation range. Copyright © 2016 John Wiley & Sons, Ltd.     

Comparison of three models to estimate evapotranspiration for a temperate mixed forest

Forest evapotranspiration is one of the main components in the regional water budget. A comparison between measured and estimated eddy covariance (EC) data, considering the Katerji–Perrier (KP), Todorovic (TD) and Priestley–Taylor (PT) actual evapotranspiration methods, was carried out. These models, relying on more easily obtainable data, are valuable when long‐term direct measurements are not available. The objective of this paper is to compare the effectivity of these three models. In this paper, experimental data were obtained within the temperate mixed forest of broad‐leaved and coniferous trees of the Changbai Mountains in northeastern China during the growing seasons of 2003 to 2005. The KP method gave the most effective values for half‐hourly and daily evapotranspiration computed by summing up half‐hourly estimates, and the TD method overestimated evapotranspiration by about 30%. The diurnal courses of estimated and measured evapotranspiration showed bell curves, similar to that of net radiation, except for a slight increase at about 14:30 solar time due to a peak value of vapour pressure deficit (VPD). For the case of daily evapotranspiration using daily mean micrometeorological variables, the PT method presented the closest values to the measurements. Accuracy of estimation related to VPD negatively (especially for VPD > 1·5 kPa). The KP parameters, considered to be vegetation dependent, were a = 0·545 and b = 1·31 at the experimental site. A constant PT parameter (α = 1·18) was applied to estimated evapotranspiration. Daily values of α responded to VPD (negatively) more strongly than to soil moisture (positively) in this forest. The experiment showed the inherent limits and advantages of the three methods. The KP method, a semi‐empirical approach, was preferred to estimate half‐hourly evapotranspiration. The TD method was a mechanistic approach to estimate reference evapotranspiration and always overestimated actual evapotranspiration. The PT method, being site dependent and the simplest approach, was effective enough to estimate large time‐scale (at least daily) evapotranspiration. Copyright © 2008 John Wiley & Sons, Ltd.     

Multi‐model ensemble prediction of terrestrial evapotranspiration across north China using Bayesian model averaging

Using high‐quality dataset from 12 flux towers in north China, the performance of four evapotranspiration (ET) models and the multi‐model ensemble approaches including the simple averaging (SA) and Bayesian model average (BMA) were systematically evaluated in this study. The four models were the single‐layer Penman–Monteith (P–M) model, the two‐layer Shuttleworthe–Wallace (S–W) model, the advection–aridity (A–A) model, and a modified Priestley–Taylor (PT‐JPL). Based on the mean value of Taylor skill (S) and the regression slope between measured and simulated ET values across all sites, the order of overall performance of the individual models from the best to the worst were: S–W (0.88, 0.87), PT‐JPL (0.80, 1.17), P–M (0.63, 1.73) and A–A (0.60, 1.68) [statistics stated as (Taylor skill, regression slope)]. Here, all models used the same values of parameters, LAI and fractional vegetation cover as well as the forcing meteorological data. Thus, the differences in model performance were mainly attributed to errors in model structure. To the ensemble approach, the BMA method has the advantage of generating more skillful and reliable predictions than the SA scheme. However, successful implementation of BMA requires accurate estimates of its parameters, and some degradation in performance were observed when the BMA parameters generated from the training period were used for the validation period. Thus, it is necessary to explore the seasonal variations of the BMA parameters according the different growth stages. Finally, the optimal conditional density function of half‐hourly ET approximated well by the double‐exponential distribution. Copyright © 2016 John Wiley & Sons, Ltd.     

Evaluating best evaporation estimate model for water surface evaporation in semi‐arid region,India

Evaluating performances of four commonly used evaporation estimate methods, namely; Bowen ratio energy balance (BREB), mass transfer (MT), Priestley–Taylor (PT) and pan evaporation (PE), based on 4 years experimental data, the most effective and the reliable evaporation estimates model for the semi‐arid region of India has been derived. The various goodness‐of‐fit measures, such as; coefficient of determination (R²), index of agreement (D), root mean square error (RMSE), and relative bias (RB) have been chosen for the performance evaluation. Of these models, the PT model has been found most promising when the Bowen ratio, β is known a priori, and based on its limited data requirement. The responses of the BREB, the PT, and the PE models were found comparable to each other, while the response of the MT model differed to match with the responses of the other three models. The coefficients, β of the BREB, µ of the MT, α of the PT and K_P of the PE model were estimated as 0·07, 2·35, 1·31 and 0·65, respectively. The PT model can successfully be extended for free water surface evaporation estimates in semi‐arid India. A linear regression model depicting relationship between daily air and water temperature has been developed using the observed water temperatures and the corresponding air temperatures. The model helped to generate unrecorded water temperatures for the corresponding ambient air temperatures. Copyright © 2007 John Wiley & Sons, Ltd.     

Influences of meteorological and vegetational factors on the partitioning of the energy of a rice paddy field

Observations made in a paddy field were analysed to show the influences of meteorological and vegetational factors on the crop's energy budget. Energy budget in the paddy field was characterized by the major partitioning to latent heat flux LE and by the negative Bowen ratio B mostly in the afternoon. Canopy resistance r_c, estimated with the Penman–Monteith equation, was related to the influences of solar radiation SR, vapour pressure deficit VPD and plant height. The results demonstrated that r_c could not directly account for B but that critical canopy resistance r_cc, defined as the canopy resistance when B = 0, could be used to standardize r_c, and that r_c − r_cc proved to be a good parameter to account for B. Influences of bulk stomatal response on energy partitioning were assessed as follows: the Bowen ratio dropped below zero, while the bulk stomatal aperture dwindled with the increase of VPD. In addition, stomata of a big leaf acted to promote the partitioning to LE against the rise of SR in the condition of higher VPD. Copyright © 2004 John Wiley & Sons, Ltd.     

Evapotranspiration models of different complexity for multiple land cover types

A comparison between half‐hourly and daily measured and computed evapotranspiration (ET) using three models of different complexity, namely, the Priestley–Taylor (P‐T), the reference Penman–Monteith (P‐M) and the Common Land Model (CLM), was conducted using three AmeriFlux sites under different land cover and climate conditions (i.e. arid grassland, temperate forest and subhumid cropland). Using the reference P‐M model with a semiempirical soil moisture function to adjust for water‐limiting conditions yielded ET estimates in reasonable agreement with the observations [root mean square error (RMSE) of 64–87 W m^?2 for half‐hourly and RMSE of 0.5–1.9 mm day^?1 for daily] and similar to the complex Common Land Model (RMSE of 60–94 W m^?2 for half‐hourly and RMSE of 0.4–2.1 mm day^?1 for daily) at the grassland and cropland sites. However, the semiempirical soil moisture function was not applicable particularly for the P‐T model at the forest site, suggesting that adjustments to key model variables may be required when applied to diverse land covers. On the other hand, under certain land cover/environmental conditions, the use of microwave‐derived soil moisture information was found to be a reliable metric of regional moisture conditions to adjust simple ET models for water‐limited cases. Further studies are needed to evaluate the utility of the simplified methods for different landscapes. Copyright © 2011 John Wiley & Sons, Ltd.     

Towards a comprehensive approach to parameter estimation in land surface parameterization schemes

In climate models, the land–atmosphere interactions are described numerically by land surface parameterization (LSP) schemes. The continuing improvement in realism in these schemes comes at the expense of the need to specify a large number of parameters that are either directly measured or estimated. Also, an emerging problem is whether the relationships used in LSPs are universal and globally applicable. One plausible approach to evaluate this is to first minimize uncertainty in model parameters by calibration. In this paper, we conduct a comprehensive analysis of some model diagnostics using a slightly modified version of the Simple Biosphere 3 model for a variety of biomes located mainly in the Amazon. First, the degree of influence of each individual parameter in simulating surface fluxes is identified. Next, we estimate parameters using a multi‐operator genetic algorithm applied in a multi‐objective context and evaluate simulations of energy and carbon fluxes against observations. Compared with the default parameter sets, these parameter estimates improve the partitioning of energy fluxes in forest and cropland sites and provide better simulations of daytime increases in assimilation of net carbon during the dry season at forest sites. Finally, a detailed assessment of the parameter estimation problem was performed by accounting for the decomposition of the mean squared error to the total model uncertainty. Analysis of the total prediction uncertainty reveals that the parameter adjustments significantly improve reproduction of the mean and variability of the flux time series at all sites and generally remove seasonality of the errors but do not improve dynamical properties. Our results demonstrate that error decomposition provides a meaningful and intuitive way to understand differences in model performance. To make further advancements in the knowledge of these models, we encourage the LSP community to adopt similar approaches in the future. Copyright © 2012 John Wiley & Sons, Ltd.     

M5 model tree based modelling of reference evapotranspiration

This paper investigates the potential of M5 model tree based regression approach to model daily reference evapotranspiration using climatic data of Davis station maintained by California irrigation Management Information System (CIMIS). Four inputs including solar radiation, average air temperature, average relative humidity, and average wind speed whereas reference evapotranspiration calculated using a relation provided by the CIMIS was used as output. To compare the performance of M5 model tree in predicting the reference evapotranspiration, FAO–56 Penman–Monteith equation and calibrated Hargreaves–Samani relation was used. A comparison of results suggests that M5 model tree approach works well in comparison to both FAO–56 and calibrated Hargreaves–Samani relations. To judge the generalization capability of M5 model tree approach, model created by using the Davis data set was tested with the datasets of four different sites. Results from this part of the study suggest that M5 model tree could successfully be employed in modeling the reference evapotranspiration. Further, sensitivity analysis with M5 model tree approach suggests the suitability of solar radiation, average air temperature, average relative humidity, and average wind speed as input parameters to model the reference evapotranspiration Copyright © 2009 John Wiley & Sons, Ltd.     

A simple parameterization of bulk canopy resistance from climatic variables for estimating hourly evapotranspiration

This paper examines a model for estimating canopy resistance r_c and reference evapotranspiration ET_o on an hourly basis. The experimental data refer to grass at two sites in Spain with semiarid and windy conditions in a typical Mediterranean climate. Measured hourly ET_o values were obtained over grass during a 4 year period between 1997 and 2000 using a weighing lysimeter (Zaragoza, northeastern Spain) and an eddy covariance system (Córdoba, southern Spain). The present model is based on the Penman–Monteith (PM) approach, but incorporates a variable canopy resistance r_c as an empirical function of the square root of a climatic resistance r* that depends on climatic variables. Values for the variable r_c were also computed according to two other approaches: with the r_c variable as a straight‐line function of r* (Katerji and Perrier, 1983, Agronomie 3 (6): 513–521) and as a mechanistic function of weather variables as proposed by Todorovic (1999, Journal of Irrigation and Drainage Engineering, ASCE 125 (5): 235–245). In the proposed model, the results show that r_c/r_a (where r_a is the aerodynamic resistance) presents a dependence on the square root of r*/r_a, as the best approach with empirically derived global parameters. When estimating hourly ET_o values, we compared the performance of the PM equation using those estimated variable r_c values with the PM equation as proposed by the Food and Agriculture Organization, with a constant r_c = 70 s m^?1. The results confirmed the relative robustness of the PM method with constant r_c, but also revealed a tendency to underestimate the measured values when ET_o is high. Under the semiarid conditions of the two experimental sites, slightly better estimates of ET_o were obtained when an estimated variable r_c was used. Although the improvement was limited, the best estimates were provided by the Todorovic and the proposed methods. The proposed approach for r_c as a function of the square root of r* may be considered as an alternative for modelling r_c, since the results suggest that the global coefficients of this locally calibrated relationship might be generalized to other climatic regions. It may also be useful to incorporate the effects of variable canopy resistances into other climatic and hydrological models. Copyright © 2005 John Wiley & Sons, Ltd.     

Peat porewater chloride concentration profiles in the Everglades during wet/dry cycles from January 1996 to June 1998: field measurements and theoretical analysis

Water quality is a key aspect of the Everglades Restoration Project, the largest water reclamation and ecosystem management project proposed in the United States. Movement of nutrients and contaminants to and from Everglades peat porewater could have important consequences for Everglades water quality and ecosystem restoration activities. In a study of Everglades porewater, we observed complex, seasonally variable peat porewater chloride concentration profiles at several locations. Analyses and interpretation of these changing peat porewater chloride concentration profiles identifies processes controlling conservative solute movement at the peat–surface water interface, that is, solutes whose transport is minimally affected by chemical and biological reactions. We examine, with an advection–diffusion model, how alternating wet and dry climatic conditions in the Florida Everglades mediate movement of chloride between peat porewater and marsh surface water. Changing surface water–chloride concentrations alter gradients at the interface between peat and overlying water and hence alter chloride flux across that interface. Surface water chloride concentrations at two frequently monitored sites vary with marsh water depth, and a transfer function was developed to describe daily marsh surface water chloride concentration as a function of marsh water depth. Model results demonstrate that porewater chloride concentrations are driven by changing surface water chloride concentrations, and a sensitivity analysis suggests that inclusion of advective transport in the model improves the agreement between the calculated and the observed chloride concentration profiles. Published in 2007 by John Wiley & Sons, Ltd.     

Modelling evapotranspiration in an alpine grassland ecosystem on Qinghai‐Tibetan plateau

Thus far, measurements and estimations of actual evapotranspiration (ET) from high‐altitude grassland ecosystems in remote areas like the Qinghai‐Tibetan plateau are still insufficient. To address these issues, a comparison between the results of the eddy covariance (EC) measurements and the estimates, considering the Katerji and Perrier (KP), the Todorovic (TD) and the Priestley–Taylor (PT) models, was carried out over an alpine grassland (38^o03'1.7'' N, 100^o 27’ 26'' E; 3032 m a.s.l.) during the growing seasons in 2008 and 2009. The results indicated that the KP model after a particularly simple calibration gave the most effective ET values in different time scales, the PT model slightly underestimate ET at night and the TD model significantly overestimated ET at noon. In addition, the canopy resistance calculated by the TD model was completely different from that calculated using the inverted EC‐measured data and the KP model, which may be due to some unrealistic assumptions made by the TD model. The KP parameters were a = 0.17 and b = 1.50 for the alpine grassland and appeared to be interannually stable. However, the PT parameter showed some interannual variations (α = 0.83 and 0.74 for 2008 and 2009, respectively). Therefore, the KP model was preferred to estimate the actual ET at both hourly and daily time scales. The PT model, being the simplest approach and field condition dependent, was recommended when available weather data were rare. On the contrary, the TD model always overestimated the actual ET and should be avoided in case of the alpine grassland ecosystems. Copyright © 2012 John Wiley & Sons, Ltd.     

Gullying and erosion control at archaeological sites in Grand Canyon,Arizona

Gully erosion of cultural sites in Grand Canyon National Park is an urgent management problem that has intensified in recent decades, potentially related to the effects of Glen Canyon Dam. We studied 25 gullies at nine sites in Grand Canyon over the 2002 monsoon–erosion season to better understand the geomorphology of the gully erosion and the effectiveness of erosion‐control structures (ECS) installed by the park under the direction of the Zuni Conservation Program. Field results indicate that Hortonian overland flow leads to concentrated flow in gullies and erosion focused at knickpoints along channels as well as at gully heads. Though groundcover type, soil shear strength and permeability vary systemat‐ically across catchments, gradient and, to a lesser degree, contributing drainage area seem to be the first‐order controls on gully extent, location of new knickpoints, and ECS damage. The installed ECS do reduce erosion relative to reaches without them and initial data suggest woody checkdams are preferable to rock linings, but maintenance is essential because damaged structures can exacerbate erosion. Topographic data from intensive field surveys and detailed photogrammetry provide slope–contributing area data for gully heads that have a trend consistent with previous empirical and theoretical formulations from a variety of landscapes. The same scaling holds below gully heads for knickpoint and ECS topographic data, with threshold coefficients the lowest for gully heads, slightly higher for knickpoints, and notably higher for damaged ECS. These topographic thresholds were used with 10‐cm digital elevation models to create simple predictive models for gully extent and structure damage. The model predictions accounted for the observed gullies but there are also many false‐positives. Purely topographical models are probably inadequate at this scale and application, but models that also parameterize the variable soil properties across sites would be useful for predicting erosion problems and ECS failure. Copyright © 2006 John Wiley & Sons, Ltd.     

Comparison of the Penman‐Monteith method and regional calibration of the Hargreaves equation for actual evapotranspiration using SWAT-simulated results in the Seolma-cheon basin,South Korea

ABSTRACT

The Hargreaves method provides reference evapotranspiration (ET_o) estimates when only air temperature data are available, although it requires previous local calibration for an acceptable performance. This method was evaluated using the data from 71 meteorological stations in the Seolma-cheon basin (8.48 km²), South Korea, comparing daily estimates against those from the Penman‐Monteith (PM) method, which was used as the standard. To estimate reference ET_o more exactly, considering the climatological characteristics in South Korea, parameter regionalization of the Hargreaves equation is carried out. First, the modified Hargreaves equation is presented after an analysis of the relationship between solar radiation and temperature. Second, parameter (K_ET) optimization of the regional calibration of the Hargreaves equation (RCH) is performed using the PM method and the modified equation at 71 meteorological stations. Next, an application was carried out to evaluate the evapotranspiration methods (PM, original Hargreaves and RCH) in the SWAT (Soil and Water Assessment Tool) model by comparing these with the measured actual evapotranspiration (AET) in the basin. The SWAT model was calibrated using 3 years (2007–2009) of daily streamflow at the watershed outlet and 3 years (2007–2009) of daily AET measured at a mixed forest. The model was validated with 3 years (2010‐2012) of streamflow and AET. RCH will contribute to a better understanding of evapotranspiration of an ungauged watershed in areas where meteorological information is scarce.

EDITOR D. Koutsoyiannis ASSOCIATE EDITOR Not assigned     

Comparison of automatic procedures for selecting flood peaks over threshold based on goodness‐of‐fit tests

In comparison with the traditional analysis of annual maximums, the peaks over threshold method provides many advantages when performing flood frequency analysis and trend analysis. However, the choice of the threshold remains an important question without definite answers and common visual diagnostic tools are difficult to reproduce on a large scale. This study investigates the behaviour of some automatic methods for threshold selection based on the generalized Pareto model for flood peak exceedances of the threshold and the Anderson–Darling test for fitting this model. In particular, the choice of a critical significance level to define an interval of acceptable values is addressed. First, automatic methods are investigated using a simulation study to assess fitting and prediction performance in a controlled environment. It is shown that P values approximated by an existing table of critical values can speed up computation without affecting the quality of the outcomes. Second, a case study compares automatically and manually selected thresholds for 285 sites across Canada by flood regime and super regions based on site characteristics. Correspondences are examined in terms of prediction of flood quantiles and trend analysis. Results show that trend detection is sensitive to the threshold selection method when studying the evolution of the number of peaks per year. Finally, a hybrid method is developed to combine automatic methods and is calibrated on the basis of super regions. The outcomes of the hybrid method are shown to more closely reproduce the estimates of the manually selected thresholds while reducing the model uncertainty.     

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.     

Stability of the time‐domain analysis method including a frequency‐dependent soil–foundation system

A number of methods have been proposed that utilize the time‐domain transformations of frequency‐dependent dynamic impedance functions to perform a time‐history analysis. Though these methods have been available in literature for a number of years, the methods exhibit stability issues depending on how the model parameters are calibrated. In this study, a novel method is proposed with which the stability of a numerical integration scheme combined with time‐domain representation of a frequency‐dependent dynamic impedance function can be evaluated. The method is verified with three independent recursive parameter models. The proposed method is expected to be a useful tool in evaluating the potential stability issue of a time‐domain analysis before running a full‐fledged nonlinear time‐domain analysis of a soil–structure system in which the dynamic impedance of a soil–foundation system is represented with a recursive parameter model. Copyright © 2015 John Wiley & Sons, Ltd.     

Seasonal variation in the energy and water exchanges above and below a larch forest in eastern Siberia

The water and energy exchanges in forests form one of the most important hydro‐meteorological systems. There have been far fewer investigations of the water and heat exchange in high latitude forests than of those in warm, humid regions. There have been few observations of this system in Siberia for an entire growing season, including the snowmelt and leaf‐fall seasons. In this study, the characteristics of the energy and water budgets in an eastern Siberian larch forest were investigated from the snowmelt season to the leaf‐fall season. The latent heat flux was strongly affected by the transpiration activity of the larch trees and increased quickly as the larch stand began to foliate. The sensible heat dropped at that time, although the net all‐wave radiation increased. Consequently, the seasonal variation in the Bowen ratio was clearly ‘U’‐shaped, and the minimum value (1·0) occurred in June and July. The Bowen ratio was very high (10–25) in early spring, just before leaf opening. The canopy resistance for a big leaf model far exceeded the aerodynamic resistance and fluctuated over a much wider range. The canopy resistance was strongly restricted by the saturation deficit, and its minimum value was 100 s m^?1 (10 mm s^?1 in conductance). This minimum canopy resistance is higher than values obtained for forests in warm, humid regions, but is similar to those measured in other boreal conifer forests. It has been suggested that the senescence of leaves also affects the canopy resistance, which was higher in the leaf‐fall season than in the foliated season. The mean evapotranspiration rate from 21 April 1998 to 7 September 1998 was 1·16 mm day^?1, and the maximum rate, 2·9 mm day^?1, occurred at the beginning of July. For the growing season from 1 June to 31 August, this rate was 1·5 mm day^?1. The total evapotranspiration from the forest (151 mm) exceeded the amount of precipitation (106 mm) and was equal to 73% of the total water input (211 mm), including the snow water equivalent. The understory evapotranspiration reached 35% of the total evapotranspiration, and the interception evaporation was 15% of the gross precipitation. The understory evapotranspiration was high and the interception evaporation was low because the canopy was sparse and the leaf area index was low. Copyright © 2001 John Wiley & Sons, Ltd.