Modelling hourly evapotranspiration and soil water content at the grass-covered boundary-layer field site Falkenberg,Germany

Abstract

The study analyses a 2-year period of hourly rates of real evapotranspiration (ETr) derived from eddy covariance measurements and soil water contents at depths from 8 to 90 cm, monitored by time domain reflectometry probes at the grass-covered boundary-layer field site Falkenberg of the Lindenberg Meteorological Observatory – Richard-Aßmann-Observatory, operated by the German Meteorological Service (DWD). The ETr rates and soil water contents were compared with the results of a modelling approach consisting of the Penman-Monteith equation and the soil water balance model Hydrus-1D using a noncompensatory and a compensatory root-water uptake model. After optimization of soil hydraulic parameters by inverse modelling, using measured soil water contents as the objective function, simulated and measured model outputs showed good agreement for soil water contents above 90 cm depth and for ETr rates simulated by our modelling approaches using noncompensatory root-water uptake. The application of a compensatory root-water uptake model led to a decrease in the simulation quality for the total investigation period.

Editor Z.W. Kundzewicz

Citation Wegehenkel, M. and Beyrich, F., 2014. Modelling of hourly evapotranspiration and soil water content at the grass-covered boundary-layer field site Falkenberg, Germany. Hydrological Sciences Journal, 59 (2), 376–394.     

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

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

Comparison of root water uptake functions to simulate surface energy fluxes within a deep‐rooted desert shrub ecosystem

Root water uptake (RWU) is a unique process whereby plants obtain water from soil, and it is essential for plant survival. The mechanisms of RWU are well understood, but their parameterization and simulation in current Land Surface Models (LSMs) fall short of the requirements of modern hydrological and climatic modelling research. Though various RWU functions have been proposed for potential use in LSMs, none was proven to be applicable for dryland ecosystems where drought was generally the limiting factor for ecosystem functioning. This study investigates the effect of root distribution on the simulated surface energy fluxes by incorporating the observed vertical root distribution. In addition, three different RWU functions were integrated into the Common Land Model (CLM) in place of the default RWU function. A comparison of the modified model's results with the measured surface energy fluxes measured by eddy covariance techniques in a Central Asian desert shrub ecosystem showed that both RWU function and vertical root distribution were able to significantly impact turbulent fluxes. Parameterizing the root distribution based on in‐situ measurement and replacing the default RWU function with a revised version significantly improved the CLM's performance in simulating the latent and sensible heat fluxes. Sensitivity analysis showed that varying the parameter values of the revised RWU function did not significantly impact the CLM's performance, and therefore, this function is recommended for use in the CLM in Central Asian desert ecosystems and, possibly, other similar dryland ecosystems. Copyright © 2013 John Wiley & Sons, Ltd.     

Tree root systems competing for soil moisture in a 3D soil–plant model

Competition for water among multiple tree rooting systems is investigated using a soil–plant model that accounts for soil moisture dynamics and root water uptake (RWU), whole plant transpiration, and leaf-level photosynthesis. The model is based on a numerical solution to the 3D Richards equation modified to account for a 3D RWU, trunk xylem, and stomatal conductances. The stomatal conductance is determined by combining a conventional biochemical demand formulation for photosynthesis with an optimization hypothesis that selects stomatal aperture so as to maximize carbon gain for a given water loss. Model results compare well with measurements of soil moisture throughout the rooting zone, of total sap flow in the trunk xylem, as well as of leaf water potential collected in a Loblolly pine forest. The model is then used to diagnose plant responses to water stress in the presence of competing rooting systems. Unsurprisingly, the overlap between rooting zones is shown to enhance soil drying. However, the 3D spatial model yielded transpiration-bulk root-zone soil moisture relations that do not deviate appreciably from their proto-typical form commonly assumed in lumped eco-hydrological models. The increased overlap among rooting systems primarily alters the timing at which the point of incipient soil moisture stress is reached by the entire soil–plant system.     

Validation of a soil water balance model using soil water content and pressure head data

The validation of soil water balance models and the evaluation of the quality of the model predictions at field‐scale require time‐series of in situ measured model outputs. In our study, we have validated such a model using a 6‐year period with time‐series of automatically recorded, daily volumetric soil water contents measured with the time‐domain reflectometry with intelligent microelements (TRIME) method and daily pressure heads measured with tensiometers. The comparisons of simulated with measured soil water contents and pressure heads were analysed using the modelling efficiency index (IA) and the square root of the mean square error (RMSE) in order to evaluate the prediction quality of the model. In our study, IA and RMSE, obtained either from the comparison of simulated with measured soil water contents or the comparison of calculated with observed pressure heads, in some cases lead to different results regarding the evaluation of the simulation quality of the soil water balance model. For example, a good fit between simulated and observed soil water contents does not necessarily result in a comparably good fit between the corresponding calculated and measured pressure heads. Therefore, a combined use of both measurement techniques, which takes into account their respective advantages and disadvantages, gives a more complete overview on the simulation quality of the soil water balance model than the single use of one of those techniques. Copyright © 2004 John Wiley & Sons, Ltd.     

Comparative analysis of relationships between NLDAS‐2 forcings and model outputs

This study reports results from an analysis of the relationship between atmospheric forcing and model‐simulated water and energy fluxes for the North American Land Data Assimilation System Project Phase 2 (NLDAS‐2). The relationships between mean monthly precipitation and total runoff are stronger in the Sacramento (SAC) and variable infiltration capacity (VIC) models, which grew out of the hydrological community, than in the Noah and Mosaic models, which grew out of the soil‐vegetation‐atmosphere transfer (SVAT) community. The reverse is true for the relationship between mean monthly precipitation and evapotranspiration. In addition, surface energy fluxes in VIC are less sensitive to model forcing (except for air temperature) than those in the Noah and Mosaic model. Notwithstanding these general conclusions, the relationships between forcings and model‐simulated water and energy fluxes for all models vary for different seasons, variables, and regions. These findings will ultimately inspire a combination of SVAT‐type model energy components with hydrological model water components to develop a SVAT‐hydrology model to improve both evapotranspiration and total runoff simulations. Copyright © 2011 John Wiley & Sons, Ltd.     

A modified ISBA surface scheme for modeling the hydrology of Athabasca River Basin with GCM-scale data

A soil–vegetation–atmosphere transfer model (SVAT), interactions between the soil–biosphere–atmosphere (ISBA) of Météo France, is modified and applied to the Athabasca River Basin (ARB) to model its water and energy fluxes. Two meteorological datasets are used: the archived forecasts from the Meteorological Survey of Canada’s Global Environmental Multiscale Model (GEM) and the European Centre for Mid-range Weather Forecasts global re-analysis (ERA-40), representing spatial scales typical of a weather forecasting model and a global circulation model (GCM), respectively. The original treatment of soil moisture and rainfall in ISBA (OISBA) is modified to statistically account for sub-grid heterogeneity of soil moisture and rainfall to produce new, highly non-linear formulations for surface and sub-surface runoff (MISBA). These new formulations can be readily applied to most existing SVATs. Stand alone mode simulations using the GEM data demonstrate that MISBA significantly improves streamflow predictions despite requiring two fewer parameters than OISBA. Simulations using the ERA-40 data show that it is possible to reproduce the annual variation in monthly, mean annual, and annual minimum flows at GCM scales without using downscaling techniques. Finally, simulations using a simple downscaling scheme show that the better performance of higher resolution datasets can be primarily attributed to improved representation of local variation of land cover, topography, and climate.     

A soil‐water‐balance approach to quantify groundwater recharge from irrigated cropland in the North China Plain

Rapidly depleting unconfined aquifers are the primary source of water for irrigation on the North China Plain. Yet, despite its critical importance, groundwater recharge to the Plain remains an enigma. We introduce a one‐dimensional soil‐water‐balance model to estimate precipitation‐ and irrigation‐generated areal recharge from commonly available crop and soil characteristics and climate data. To limit input data needs and to simplify calculations, the model assumes that water flows vertically downward under a unit gradient; infiltration and evapotranspiration are separate, sequential processes; evapotranspiration is allocated to evaporation and transpiration as a function of leaf‐area index and is limited by soil‐moisture content; and evaporation and transpiration are distributed through the soil profile as exponential functions of soil and root depth, respectively. For calibration, model‐calculated water contents of 11 soil‐depth intervals from 0 to 200 cm were compared with measured water contents of loam soil at four sites in Luancheng County, Hebei Province, over 3 years (1998–2001). Each 50‐m² site was identically cropped with winter wheat and summer maize, but received a different irrigation treatment. Average root mean‐squared error between measured and model‐calculated water content of the top 180 cm was 4·2 cm, or 9·3% of average total water content. In addition, model‐calculated evapotranspiration compared well with that measured by a large‐scale lysimeter. To test the model, 12 additional sites were simulated successfully. Model results demonstrate that drainage from the soil profile is not a constant fraction of precipitation and irrigation inputs, but rather the fraction increases as the inputs increase. Because this drainage recharges the underlying aquifer, improving irrigation efficiency by reducing seepage will not reverse water‐table declines. Copyright © 2003 John Wiley & Sons, Ltd.     

Impact of spatial variations of land surface parameters on regional evaporation: a case study with remote sensing data

Most precipitation in watersheds is consumed by evaporation, thus techniques to appraise regional evaporation are important to assess the availability of water resources. Many algorithms to estimate evaporation from remotely sensed spectral data have been developed in the recent past. In addition to differences in the physical parameterization of surface fluxes, these algorithms have different solutions for describing spatial variations of the parameters in the soil–vegetation–atmosphere–transfer (SVAT) continuum. In this study, the necessity to spatially distinguish SVAT parameters for computing surface heat fluxes is analysed for the Naivasha watershed in the Kenyan Rift Valley. Landsat Thematic Mapper (TM) spectral data have been used to first delineate the watershed into 15 hydrological units using surface temperature, normalized difference vegetation index and surface albedo as attributes. Thereafter, semi‐empirical relationships between these TM‐based parameters and other SVAT parameters have been applied to compute the spatial variation of SVAT parameters and the associated evaporation from the different hydrological units. The impact of using watershed‐constant or watershed‐distributed SVAT parameters on the fluxes is analysed. The determination of watershed averaged evaporation with area‐aggregated SVAT parameters is feasible without significant loss of accuracy. Distributed evaporation in heterogeneous watersheds, however, can be investigated only with remote sensing flux algorithms that can account for spatially variable air temperature, surface roughness, surface albedo and the stability correction of the temperature profile due to buoyancy. Erroneous results can be expected if area‐aggregated SVAT parameters are used to calculate local evaporation. As most of the recently developed remote sensing flux algorithms are based on areal constant SVAT parameters, direct applications in watersheds are still limited. Copyright © 2001 John Wiley & Sons, Ltd.     

Root controls on water redistribution and carbon uptake in the soil–plant system under current and future climate

Understanding photosynthesis and plant water management as a coupled process remains an open scientific problem. Current eco-hydrologic models characteristically describe plant photosynthetic and hydraulic processes through ad hoc empirical parameterizations with no explicit accounting for the main pathways over which carbon and water uptake interact. Here, a soil–plant-atmosphere continuum model is proposed that mechanistically couples photosynthesis and transpiration rates, including the main leaf physiological controls exerted by stomata. The proposed approach links the soil-to-leaf hydraulic transport to stomatal regulation, and closes the coupled photosynthesis–transpiration problem by maximizing leaf carbon gain subject to a water loss constraint. The approach is evaluated against field data from a grass site and is shown to reproduce the main features of soil moisture dynamics and hydraulic redistribution. In particular, it is shown that the differential soil drying produced by diurnal root water uptake drives a significant upward redistribution of moisture both through a conventional Darcian flow and through the root system, consistent with observations. In a numerical soil drying experiment, it is demonstrated that more than 50% of diurnal transpiration is supplied by nocturnal upward water redistribution, and some 12% is provided directly through root hydraulic redistribution. For a prescribed leaf area density, the model is then used to diagnose how elevated atmospheric CO₂ concentration and increased air temperature jointly impact soil moisture, transpiration, photosynthesis, and whole-plant water use efficiency, along with compensatory mechanisms such as hydraulic lift using several canonical forms of root-density distribution.     

Information and complexity measures applied to observed and simulated soil moisture time series

Abstract

Time series of soil moisture-related parameters provide important insights into the functioning of soil water systems. Analysis of patterns within such time series has been used in several studies. The objective of this work was to compare patterns in observed and simulated soil moisture contents to understand whether modelling leads to a substantial loss of information or complexity. The time series were observed at four plots in sandy soils within the USDA-ARS OPE3 experimental watershed, for a year; precipitation and evapotranspiration (ET) were measured and estimated, respectively, and used for soil water flow simulation with the HYDRUS-1D software. The information content measures are the metric entropy and the mean information gain, and complexity measures are the fluctuation complexity and the effective measure complexity. These measures were computed based on the binary encoding of soil moisture time series, and used probabilities of patterns, i.e. probabilities of joint or sequential appearance of symbol sequences. The information content of daily soil moisture time series was much smaller than that of rainfall data, and had higher complexity, indicating that soil worked essentially as an information filter. Information content and complexity decreased and increased with depth, respectively, demonstrating the increase in the information filtering action of soil. The information measures of simulated soil moisture content were close to those of the measurements, indicating the successful simulation of patterns in the data. The spatial variability of the information measures for simulated soil moisture content at all depths was less pronounced than the one of measured time series. Compared with precipitation and estimated ET, soil moisture time series had more structure and less randomness in this work. The information measures can provide useful complementary knowledge about model performance and patterns in observation and modelling results.

Citation Pan, F., Pachepsky, Y. A., Guber, A. K., & Hill, R. L. (2011) Information and complexity measures applied to observed and simulated soil moisture time series. Hydrol. Sci. J. 56(6), 1027–1039.     

Preferential percolation quantified by large water content sensors with artifactual macroporous envelopes

For many scientific and practical tasks, it is important to estimate the soil–water percolation fluxes. This paper builds on measurements with large horizontal time‐domain reflectometry water content sensors in a loamy Mollisol. The sensors were installed into pre‐drilled holes and the gaps between them, and the soil was filled with a slurry of local soil with water. This gave rise to envelopes around them that contained artificial macropores. The sensors reacted to intensive rains by a rapid increase of their readings, often above the native soil's porosity, followed by an almost equally rapid decrease. The paper explores the feasibility of quantifying the rapid percolation, based on these anomalous water content peaks, and demonstrates that this is possible in principle, if the processes are simulated by a suitable model. A two‐dimensional dual porosity non‐equilibrium (mobile‐immobile) model was tried. The envelope around the sensor was modelled as an annulus with higher porosity and hydraulic conductivity, which attracts preferential flow and amplifies the percolation signal. With the model at hand, the flux hydrographs can be derived from model simulations and measured precipitation. For contrast, the Durner equilibrium dual porosity model was tried but was found little suitable. However, even the mobile‐immobile model did not perform perfectly. Simulated water contents were similar to the measured ones at some depths but not in the others, and the percolation fluxes were overestimated, compared to cumulative soil–water balance. Efforts to improve model performance were not successful. Hence, the model structure needs to be improved. Copyright © 2015 John Wiley & Sons, Ltd.     

Coupling two radar backscattering models to assess soil roughness and surface water content at farm scale

Abstract

Remote sensing techniques are useful for agro-hydrological monitoring at the farm scale because the availability of spatially and temporally distributed data improves agricultural models for irrigation and crop yield optimization under water scarcity conditions. This research focuses on the surface water content retrieval using active microwave data. Two semi-empirical models were chosen as these showed the best performances in simulating cross and co-polarized backscatter. Thus, these models were coupled to obtain reliable assessments of both soil water content and soil roughness. The use of the coupled model enables one to avoid using roughness measured in situ. Remote sensing images and in situ data were collected between April and July 2006 within the European Space Agency-funded project AgriSAR 2006. The images data set includes L-band in HH, VV and VH polarizations acquired from the airborne E-SAR sensor, operated by the German Aerospace Centre. Results were validated using in situ soil water content and roughness measurements. The results show that reliable assessment of both soil roughness (r ² up to ?0.8) and soil water content (r ² ? 0.9) can be retrieved in fields characterized by low fractional coverage.

Editor D. Koutsoyiannis; Associate editor C. Onof

Citation Capodici, F., Maltese, A., Ciraolo, G., La Loggia, G., and D’Urso, G., 2013. Coupling two radar backscattering models to assess soil roughness and surface water content at the farm scale. Hydrological Sciences Journal, 58 (8), 1677–1689.     

Potentials and limitations of modelling vertical distributions of root water uptake of an Austrian pine forest on a sandy soil

Root water uptake patterns are often studied with simulation models of the unsaturated soil water flow, as they are difficult to measure directly. Calibration of these models is not straightforward and causes uncertainties in simulated uptake distributions. In this paper we study how uncertainties in the calibration of the SWIF model affect uncertainty intervals in simulated uptake patterns of an Austrian pine stand (Pinus nigra var. nigra) on a sandy soil. After calibrating and validating SWIF with a large data set of more than 125 000 measured soil water contents over a three year period, uncertainty ranges in simulated soil water dynamics and root water uptake distributions were estimated with a Monte Carlo analysis. In general, uncertainties in root uptake patterns were small (typically <2 10⁻⁴ m³ m⁻³ day⁻¹) and were higher for trees with a shallow rooting system (0·8 m) than for trees with a deep rooting system (2·5 m). Uncertainties arose mainly from uncertainties in simulated soil water fluxes and from variations in the reduction of uptake during periods of drought. Uncertainties in soil water contents were far higher (typically 0·01 m³ m⁻³) than uncertainties in uptake, illustrating that uncertainties in uptake parameters and those in the distribution of water uptake hardly affect the modelling of soil water dynamics. Root water uptake models should therefore be validated against measured uptake distributions, which can be determined on sandy soils during dry periods with a high water use when soil fluxes are negligible to uptake. Copyright © 2000 John Wiley & Sons, Ltd.     

Estimating evapotranspiration from small watersheds using a water and energy balance approach

Meteorological and environmental data measured in semiarid watersheds during the summer monsoon and winter periods were used to study the interrelationships among flux, meteorological and soil water variables, and to evaluate the effects of these variables on the daily estimation of actual evapotranspiration (AET). The relationship between AET and potential evapotranspiration (PET) as a function of soil water content, as suggested by Thornthwaite–Mather and by Morton, was studied to determine its applicability to the study area. Furthermore, multiple linear regression (MLR) analysis was employed to evaluate the order of importance of the meteorological and soil water factors involved. The results of MLR analysis showed that the combined effects of available energy, soil water content and wind speed were responsible for more than 70% of the observed variations in AET during the summer monsoon period. The analyses also indicate that the combined effects of available energy, vapour pressure deficit and wind speed were responsible for more than 70% of the observed variations in AET during the winter period. However, the test results of two different approaches, using the relationships between AET and PET as a function of soil water content, indicated some inadequacy. Copyright © 2007 John Wiley & Sons, Ltd.     

Ammonium Transport Simulation in Horizontal Soil Columns from a Natural Inland Alkaline Wetland

Ammonium transport was simulated in horizontal soil columns from an inland alkaline wetland (Fulaowenpao wetland) of Northeast China. The primary objectives of this work are to investigate the changes in ammonium transport rate with increasing distances along horizontal soil column and to determine the effects of water diffusion rate and volumetric water content on ammonium transport rate. Our results showed that water diffusion coefficient was the lowest at the soil layer of 10–20 cm, followed by the 0–10 cm soil layer, and the highest value occured at the soil layer of 20–60 cm. The highest ammonium transport rate also appeared at the soil layer of 20–60 cm, while the lowest value was observed at the soil layer of 10–20 cm. Ammonium transport rates decreased with increasing distances along horizontal soil columns. The ammonium transport rates showed higher values at the distance from 0 to 6 cm and then decreased rapidly from 6 to 18 cm. However, they nearly kept constant and approached to zero after exceeding the distance of 18 cm. Ammonium transport rates increased exponentially with increasing volumetric water contents and water diffusion rates. The alkaline wetland soils prevented ammonium from horizontal diffusion at all soil layers under drying conditions.     

Multi‐scale assimilation of root zone soil water predictions

When hydrology model parameters are determined, a traditional data assimilation method (such as Kalman filter) and a hydrology model can estimate the root zone soil water with uncertain state variables (such as initial soil water content). The simulated result can be quite good. However, when a key soil hydraulic property, such as the saturated hydraulic conductivity, is overestimated or underestimated, the traditional soil water assimilation process will produce a persistent bias in its predictions. In this paper, we present and demonstrate a new multi‐scale assimilation method by combining the direct insertion assimilation method, particle swarm optimisation (PSO) algorithm and Richards equation. We study the possibility of estimating root zone soil water with a multi‐scale assimilation method by using observed in situ data from the Wudaogou experiment station, Huaihe River Basin, China. The results indicate there is a persistent bias between simulated and observed values when the direct insertion assimilation surface soil water content is used to estimate root zone soil water contents. Using a multi‐scale assimilation method (PSO algorithm and direct insertion assimilation) and an assumed bottom boundary condition, the results show some obvious improvement, but the root mean square error is still relatively large. When the bottom boundary condition is similar to the actual situation, the multi‐scale assimilation method can well represent the root zone soil water content. The results indicate that the method is useful in estimating root zone soil water when available soil water data are limited to the surface layer and the initial soil water content even when the soil hydraulic conductivities are uncertain. Copyright © 2011 John Wiley & Sons, Ltd.     

Hammond Hill Research Catchment: Supporting hydrologic investigations of rooting zone and vegetation water dynamics under climate change

The Hammond Hill Research Catchment (HH) is a small (120 ha), temperate, second order tributary to Six Mile Creek, Cayuga Lake, and the Great Lakes (42.42°, −76.32°). The HH has been monitored since January 2017 for the purpose of understanding how recent infiltration mixes with antecedent soil water on hillslope forest floors and the spatial and temporal patterns of Root Water Uptake (RWU) by temperate northeastern US tree species (eastern hemlock [Tsuga canadensis], American beech [Fagus grandifolia], and sugar maple [Acer saccharum]). These data are informing us about the hydrologic consequences of anticipated tree species composition change and supporting the development of more refined ecohydrological models. The glaciated catchment is underlain by a shallow confining siltstone layer (1–1.5 m depth) and densely covered with an approximately 60 year old regrowth mixed species forest of hemlock, beech, and other deciduous tree species common to the northeastern US. Current datasets from the HH include precipitation snow water equivalent, discharge, and associated isotopic water compositions, δ²H & δ¹⁸O. Measurements of (top 10 cm) soil water content, as well as bulk soil water and hemlock and beech xylem isotopic compositions are made at several locations across a topographic wetness gradient. The near-term role of the HH is to support an understanding of the environmental and ecological drivers of plant RWU competition. All data from the HH are publicly available.     

Simulation of crop growth and energy and carbon dioxide fluxes at different time steps from hourly to daily

Understanding the exchange processes of energy and carbon dioxide (CO₂) in the soil–vegetation–atmosphere system is important for assessing the role of the terrestrial ecosystem in the global water and carbon cycle and in climate change. We present a soil–vegetation–atmosphere integrated model (ChinaAgrosys) for simulating energy, water and CO₂ fluxes, crop growth and development, with ample supply of nutrients and in the absence of pests, diseases and weed damage. Furthermore, we test the hypotheses of whether there is any significant difference between simulations over different time steps. CO₂, water and heat fluxes were estimated by the improving parameterization method of the coupled photosynthesis–stomatal conductance–transpiration model. Soil water evaporation and plant transpiration were calculated using a multilayer water and heat‐transfer model. Field experiments were conducted in the Yucheng Integrated Agricultural Experimental Station on the North China Plain. Daily weather and crop growth variables were observed during 1998–2001, and hourly weather variables and water and heat fluxes were measured using the eddy covariance method during 2002–2003. The results showed that the model could effectively simulate diurnal and seasonal changes of net radiation, sensible and latent heat flux, soil heat flux and CO₂ fluxes. The processes of evapotranspiration, soil temperature and leaf area index agree well with the measured values. Midday depression of canopy photosynthesis could be simulated by assessing the diurnal change in canopy water potential. Moreover, the comparisons of simulated daily evapotranspiration and net ecosystem exchange (NEE) under different time steps indicated that time steps used by a model affect the simulated results. There is no significant difference between simulated evapotranspiration using the model under different time steps. However, simulated NEE produces large differences in the response to different time steps. Therefore, the accurate calculation of average absorbed photosynthetic active radiation is important for the scaling of the model from hourly steps to daily steps in simulating energy and CO₂ flux exchanges between winter wheat and the atmosphere. Copyright © 2007 John Wiley & Sons, Ltd.     

High-resolution simulation of the spatial pattern of water use in continental China

ABSTRACT

High-resolution data on the spatial pattern of water use are a prerequisite for appropriate and sustainable water management. Based on one well-validated hydrological model, the Distributed Time Variant Gains Model (DTVGM), this paper obtains reliable high-resolution spatial patterns of irrigation, industrial and domestic water use in continental China. During the validation periods, ranges of correlation coefficient (R) and Nash-Sutcliffe efficiency (NSE) coefficient are 0.67–0.96 and 0.51–0.84, respectively, between the observed and simulated streamflow of six hydrological stations, indicating model applicability to simulate the distribution of water use. The simulated water use quantities have relative errors (RE) less than 5% compared with the observed. In addition, the changes in streamflow discharge were also correctly simulated by our model, such as the Zhangjiafen station in the Hai River basin with a dramatic decrease in streamflow, and the Makou station in the Pearl River basin with no significant changes. These changes are combined results of basin available water resources and water use. The obtained high-resolution spatial pattern of water use could decrease uncertainty of hydrological simulation and guide water management efficiently.

Editor M.C. Acreman; Associate editor X. Fang