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.     

Hydrological simulation of the Betwa River basin (India) using the SWAT model

A seasonal water budget analysis was carried out to quantify various components of the hydrological cycle using the Soil and Water Assessment Tool (SWAT) model for the Betwa River basin (43?500 km²) in central India. The model results were satisfactory in calibration and validation. The seasonal water budget analysis showed that about 90% of annual rainfall and 97% of annual runoff occurred in the monsoon season. A seasonal linear trend analysis was carried out to detect trends in the water balance components of the basin for the period 1973–2001. In the monsoon season, an increasing trend in rainfall and a decreasing trend in ET were observed; this resulted in an increasing trend in groundwater storage and surface runoff. The winter season followed almost the same pattern. A decreasing trend was observed in summer season rainfall. The study evokes the need for conservation structures in the study area to reduce monsoon runoff and conserve it for basin requirements in water-scarce seasons.

EDITOR Z.W. Kundzewicz

ASSOCIATE EDITOR F. Hattermann     

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

Abstract

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

Editor Z.W. Kundzewicz

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

Use of remote sensing for evapotranspiration monitoring over land surfaces

Abstract

Monitoring evapotranspiration (ET) at large scales is important for assessing climate and anthropogenic effects on natural and agricultural ecosystems. This paper describes techniques used in evaluating ET with remote sensing, which is the only technology that can efficiently and economically provide regional and global coverage. Some of the empirical/statistical techniques have been used operationally with satellite data for computing daily ET at regional scales. The more complex numerical simulation models require detailed input parameters that may limit their application to regions containing a large database of soils and vegetation properties. Current efforts are being directed towards simplifying the parameter requirements of these models. Essentially all energy balance models rely on an estimate of the available energy (net radiation less soil heat flux). Net radiation is not easily determined from space, although progress is being made. Simplified approaches for estimating soil heat flux appear promising for operational applications. In addition, most ET models utilize remote sensing data in the shortwave and thermal wavelengths to measure key boundary conditions. Differences between the radiometric surface temperature and aerodynamic temperature can be significant and progress in incorporating this effect is evident. Atmospheric effects on optical data are significant, and optical sensors cannot see through clouds. This has led some to use microwave observations as a surrogate for optical data to provide estimates of surface moisture and surface temperature; preliminary results are encouraging. The approaches that appear most promising use surface temperature and vegetation indices or a time rate of change in surface temperature coupled to an atmospheric boundary layer model. For many of these models, differences with ET observations can be as low as 20% from hourly to daily time scales, approaching the level of uncertainty in the measurement of ET and contradicting some recent pessimistic conclusions concerning the utility of remotely sensed radiometric surface temperature for determining the surface energy balance.     

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

ABSTRACT

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

Assessment of the MODIS global evapotranspiration algorithm using eddy covariance measurements and hydrological modelling in the Rio Grande basin

Abstract

Remote sensing is considered the most effective tool for estimating evapotranspiration (ET) over large spatial scales. Global terrestrial ET estimates over vegetated land surfaces are now operationally produced at 1-km spatial resolution using data from the Moderate Resolution Imaging Spectroradiometer (MODIS) and the MOD16 algorithm. To evaluate the accuracy of this product, ground-based measurements of energy fluxes obtained from eddy covariance sites installed in tropical biomes and from a hydrological model (MGB-IPH) were used to validate MOD16 products at local and regional scales. We examined the accuracy of the MOD16 algorithm at two sites in the Rio Grande basin, Brazil, one characterized by a sugar-cane plantation (USE), the other covered by natural savannah vegetation (PDG) for the year 2001. Inter-comparison between 8-day average MOD16 ET estimates and flux tower measurements yielded correlations of 0.78 to 0.81, with root mean square errors (RMSE) of 0.78 and 0.46 mm d^-1, at PDG and USE, respectively. At the PDG site, the annual ET estimate derived by the MOD16 algorithm was 19% higher than the measured amount. For the average annual ET at the basin-wide scale (over an area of 145 000 km²), MOD16 estimates were 21% lower than those from the hydrological model MGB-IPH. Misclassification of land use and land cover was identified as the largest contributor to the error from the MOD16 algorithm. These estimates improve significantly when results are integrated into monthly or annual time intervals, suggesting that the algorithm has a potential for spatial and temporal monitoring of the ET process, continuously and systematically, through the use of remote sensing data.

Editor D. Koutsoyiannis; Associate editor T. Wagener

Citation Ruhoff, A.L., Paz, A.R., Aragao, L.E.O.C., Mu, Q., Malhi, Y., Collischonn, W., Rocha, H.R., and Running, S.W., 2013. Assessment of the MODIS global evapotranspiration algorithm using eddy covariance measurements and hydrological modelling in the Rio Grande basin. Hydrological Sciences Journal, 58 (8), 1658–1676.     

Up-scaling of SEBAL derived evapotranspiration maps from Landsat (30 m) to MODIS (250 m) scale

Remotely sensed imagery of the Earth’s surface via satellite sensors provides information to estimate the spatial distribution of evapotranspiration (ET). The spatial resolution of ET predictions depends on the sensor type and varies from the 30–60 m Landsat scale to the 250–1000 m MODIS scale. Therefore, for an accurate characterization of the regional distribution of ET, scaling transfer between images of different resolutions is important. Scaling transfer includes both up-scaling (aggregation) and down-scaling (disaggregation). In this paper, we address the up-scaling problem.The Surface Energy Balance Algorithm for Land (SEBAL) was used to derive ET maps from Landsat 7 Enhanced Thematic Mapper Plus (ETM+) and Moderate Resolution Imaging Spectroradiometer (MODIS) images. Landsat 7 bands have spatial resolutions of 30–60 m, while MODIS bands have resolutions of 250, 500 and 1000 m. Evaluations were conducted for both “output” and “input” up-scaling procedures, with aggregation accomplished by both simple averaging and nearest neighboring resampling techniques. Output up-scaling consisted of first applying SEBAL and then aggregating the output variable (daily ET). Input up-scaling consisted of aggregating 30 m Landsat pixels of the input variable (radiance) to obtain pixels at 60, 120, 250, 500 and 1000 m before SEBAL was applied. The objectives of this study were first to test the consistency of SEBAL algorithm for Landsat and MODIS satellite images and second to investigate the effect of the four different up-scaling processes on the spatial distribution of ET.We conclude that good agreement exists between SEBAL estimated ET maps directly derived from Landsat 7 and MODIS images. Among the four up-scaling methods the output simple averaging method produced aggregated data and aggregated differences with the most statistically and spatially predictable behavior. The input nearest neighbor method was the least predictable but was still acceptable. Overall, the daily ET maps over the Middle Rio Grande Basin aggregated from Landsat images were in good agreement with ET maps directly derived from MODIS images.     

Can satellite land surface temperature data be used similarly to river discharge measurements for distributed hydrological model calibration?

Abstract

A new methodology is proposed for the calibration of distributed hydrological models at the basin scale by constraining an internal model variable using satellite data of land surface temperature (LST). The model algorithm solves the system of energy and mass balances in terms of a representative equilibrium temperature that governs the fluxes of energy and mass over the basin domain. This equilibrium surface temperature, which is a critical model state variable, is compared to operational satellite LST, while calibrating soil hydraulic parameters and vegetation variables differently in each pixel, minimizing the errors. This procedure is compared to the traditional calibration using only discharge measurements. The distributed energy water balance model, Flash-flood Event-based Spatially-distributed rainfall–runoff Transformation – Energy Water Balance model (FEST-EWB), is used to test this approach. This methodology is applied to the Upper Yangtze River basin (China) using MODIS LST retrieved from satellite data in the framework of the NRSCC-ESA DRAGON-2 Programme. The calibration procedure based on LST seems to outperform the calibration based on discharge, with lower relative error and higher Nash-Sutcliffe efficiency index on cumulated volume.

Editor D. Koutsoyiannis; Associate editor C. Perrin     

Intercomparison across scales between remotely-sensed land surface temperature and representative equilibrium temperature from a distributed energy water balance model

Abstract

Quantifying the reliability of distributed hydrological models is an important task in hydrology to understand their ability to estimate energy and water fluxes at the agricultural district scale as well the basin scale for water resources management in drought monitoring and flood forecasting. In this context, the paper presents an intercomparison of simulated representative equilibrium temperature (RET) derived from a distributed energy water balance model and remotely-sensed land surface temperature (LST) at spatial scales from the agricultural field to the river basin. The main objective of the study is to evaluate the use of LST retrieved from operational remote sensing data at different spatial and temporal resolutions for the internal validation of a distributed hydrological model to control its mass balance accuracy as a complementary method to traditional calibration with discharge measurements at control river cross-sections. Modelled and observed LST from different radiometric sensors located on the ground surface, on an aeroplane and a satellite are compared for a maize field in Landriano (Italy), the agricultural district of Barrax (Spain) and the Upper Po River basin (Italy). A good ability of the model in reproducing the observed LST values in terms of mean bias error, root mean square error, relative error and Nash-Sutcliffe index is shown.

Editor Z.W. Kundzewicz; Associate editor D. Gerten     

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.     

Event-based estimation of water budget components using a network of multi-sensor capacitance probes

Abstract

A new approach was developed for estimating vertical soil water fluxes using soil water content time series data. Instead of a traditional fixed time interval, this approach utilizes the time interval between two sequential minima of the soil water storage time series to identify groundwater recharge events and calculate components of the soil water budget. We calculated water budget components: surface-water excess (Sw), infiltration less evapotranspiration (I – ET) and groundwater recharge (R) from May 2001 to January 2003 at eight locations at the USDA Agricultural Research Center, Beltsville, Maryland, USA. High uncertainty was observed for all budget components. This uncertainty was attributed to spatial and temporal variation in Sw, I – ET and R, and was caused by nonuniform rainfall distributions during recharge events, variability in the profile water content, and spatial variability in soil hydraulic properties. The proposed event-based approach allows estimating water budget components when profile water content monitoring data are available.

Citation Guber, A., Gish, T., Pachepsky, Y., McKee, L., Nicholson, T. & Cady, R. (2011) Event-based estimation of water budget components using a network of multi-sensor capacitance probes. Hydrol. Sci. J. 56(7), 1227–1241.     

Energy fluxes and melt rate of a seasonal snow cover in the Moroccan High Atlas

ABSTRACT

In this study, we used an energy balance model and two simple methods based on readily available data to identify the processes driving the point-scale energy and mass balance of the snowpack. Data were provided from an experimental site located at 3200 m. All models were evaluated by comparing observed and modelled snow water equivalents. Performances are variable from one season to the next and the energy balance model gives better results (mean of root mean square error, RMSE = 25 mm and r² = 0.90) than the two simplified approaches (mean of RMSE = 54 mm and r² = 0.70). There are significant amounts of snow sublimation but they are highly variable from season to season, depending on wind conditions (between 7 and 20% of the total). While the main source of energy for melting is net radiation, the amount of heat brought by sensible heat flux is significant for two of the most windy snow seasons.

Editor Z.W. Kundzewicz Associate editor not assigned     

Trends in land surface evapotranspiration across China with remotely sensed NDVI and climatological data for 1981–2010

Abstract

Using satellite observations of Normalized Difference Vegetation Index (NDVI) from NOAA-AVHRR and Terra-MODIS, together with climatic data in a physical evapotranspiration (ET) model, the spatio-temporal variability of ET is investigated in terrestrial China from 1981 to 2010. The model predictions of actual ET (ET_a) are validated with ET values from in situ eddy covariance flux measurements and from basin water balance calculations. The national averaged crop reference ET (ET_p) and ET_a values are 916 ± 21 and 415 ± 12 mm year^-1, respectively. The annual ET_a pattern is closely associated with vegetation conditions in the eastern part of China, whereas ET_a is low in the sparsely-vegetated areas and deserts in the northwestern region, corresponding to scarce rainfall events and amounts. The trends of ET_p and ET_a are remarkably different over the country, and the complementary relationship between ET_p and ET_a is revealed for the study period. Averaged over the whole country, ET_a showed an increasing trend from the 1980s to the mid-1990s, followed by a decreasing trend, consistent with the precipitation anomaly. Across the main vegetation types, annual ET_a amounts are found to correspond clearly with the bands of precipitation and ET_p.     

Verification of land surface evapotranspiration estimation from remote sensing spatial contextual information

Estimation of evapotranspiration (ET) is of great significance in modeling the water and energy interactions between land and atmosphere. Negative correlation of surface temperature (T_s) versus vegetation index (VI) from remote sensing data provides diagnosis on the spatial pattern of surface soil moisture and ET. This study further examined the applicability of T_s–VI triangle method with a newly developed edges determination technique in estimating regional evaporative fraction (EF) and ET at MODIS pixel scale through comparison with large aperture scintillometer (LAS) and high‐level eddy covariance measurements collected at Changwu agro‐ecological experiment station from late June to late October, 2009. An algorithm with merely land and atmosphere products from MODIS onboard Terra satellite was used to estimate the surface net radiation (R_n) and soil heat flux. In most cases, the estimated instantaneous R_n was in good agreement with surface measurement with slight overestimation by 12 W/m². Validation results from LAS measurement showed that the root mean square error is 0.097 for instantaneous EF, 48 W/m² for instantaneous sensible heat flux, and 30 W/m² for daily latent heat flux. This paper successfully presents a miniature of the overall capability of T_s–VI triangle in estimating regional EF and ET from limited number of data. For a thorough interpretation, further comprehensive investigation needs to be done with more integration of remote sensing data and in‐situ surface measurements. Copyright © 2011 John Wiley & Sons, Ltd.     

Distribution of magnetic energy in αΩ-dynamos,I: The method

Abstract

In this paper a method for solving the equation for the mean magnetic energy <BB> of a solar type dynamo with an axisymmetric convection zone geometry is developed and the main features of the method are described. This method is referred to as the finite magnetic energy method since it is based on the idea that the real magnetic field B of the dynamo remains finite only if <BB> remains finite. Ensemble averaging is used, which implies that fields of all spatial scales are included, small-scale as well as large-scale fields. The method yields an energy balance for the mean energy density ε ≡ B ²/8π of the dynamo, from which the relative energy production rates by the different dynamo processes can be inferred. An estimate for the r.m.s. field strength at the surface and at the base of the convection zone can be found by comparing the magnetic energy density and the outgoing flux at the surface with the observed values. We neglect resistive effects and present arguments indicating that this is a fair assumption for the solar convection zone. The model considerations and examples presented indicate that (1) the energy loss at the solar surface is almost instantaneous; (2) the convection in the convection zone takes place in the form of giant cells; (3) the r.m.s. field strength at the base of the solar convection zone is no more than a few hundred gauss; (4) the turbulent diffusion coefficient within the bulk of the convection zone is about 10¹⁴cm²s^?1, which is an order of magnitude larger than usually adopted in solar mean field models.     

Evapotranspiration estimation considering anthropogenic heat based on remote sensing in urban area

Urbanization influences hydrologic cycle significantly on local, regional even global scale. With urbanization the water resources demand for dense population sharpened, thus it is a great challenge to ensure water supply for some metropolises such as Beijing. Urban area is traditionally considered as the area with lower evapotranspiration (ET) on account of the impervious surface and the lower wind speed. For most remote sensing models, the ET, defined as latent heat in energy budget, is estimated as the difference between net radiation and sensible heat. The sensible heat is generally higher in urban area due to the high surface temperature caused by heat island, therefore the latent heat (i.e. the ET) in urban area is lower than that in other region. We estimated water consumption from 2003 to 2012 in Beijing based on water balance method and found that the annual mean ET in urban area was about 654 mm. However, using Surface Energy Balance System (SEBS) model, the annual mean ET in urban area was only 348 mm. We attributed this inconsistence to the impact of anthropogenic heat and quantified this impact on the basis of the night-light maps. Therefore, a new model SEBS-Urban, coupling SEBS model and anthropogenic heat was developed to estimate the ET in urban area. The ET in urban area of Beijing estimated by SEBS-Urban showed a good agreement with the ET from water balance method. The findings from this study highlighted that anthropogenic heat should be included in the surface energy budget for a highly urbanized area.     

Effects of land-use/cover change on hydrological processes using a GIS/RS-based integrated hydrological model: case study of the East River,China

Abstract

An integrated model, combining a surface energy balance system, an LAI-based interception model and a distributed monthly water balance model, was developed to predict hydrological impacts of land-use/land-cover change (LUCC) in the East River basin, China, with the aid of GIS/RS. The integrated model is a distributed model that not only accounts for spatial variations in basin terrain, rainfall and soil moisture, but also considers spatial and temporal variation of vegetation cover and evapotranspiration (ET), in particular, thus providing a powerful tool for investigating the hydrological impact of LUCC. The model was constructed using spatial data on topography, soil types and vegetation characteristics together with time series of precipitation from 170 stations in the basin. The model was calibrated and validated based on river discharge data from three stations in the basin for 21 years. The calibration and validation results suggested that the model is suitable for application in the basin. The results show that ET has a positive relationship with LAI (leaf area index), while runoff has a negative relationship with LAI in the same climatic zone that can be described by the surface energy balance and water balance equation. It was found that deforestation would cause an increase in annual runoff and a decrease in annual ET in southern China. Monthly runoff for different land-cover types was found to be inversely related to ET. Also, for most of the scenarios, and particularly for grassland and cropland, the most significant changes occurred in the rainy season, indicating that deforestation would cause a significant increase in monthly runoff in that season in the East River basin. These results are important for water resources management and environmental change monitoring.

Editor Z.W. Kundzewicz     

Assessment of the crop coefficient for saltgrass under native riparian field conditions in the desert southwest

Limited urban water supplies in southwestern USA cause water managers and planners to re‐assess water losses and needs from consumptive water use by riparian vegetation. Here, we report on field measurements of evapotranspiration (ET) for inland saltgrass [Distichlis spicata var. stricta (L.) Greene]; a once common riparian plant native to the desert southwest. The objective was to develop a saltgrass crop coefficient, Kc, similar to agricultural crop coefficients commonly used in irrigation water management. The developed Kc, in conjunction with the local climate, can then be used to assess the water savings that may be achieved in riparian zones for saltgrass versus invasive species and for use in irrigation management and scheduling of saltgrass in urban setting. The ET of saltgrass was measured in its native riparian setting located in the flood plain of the Rio Grande, north of Caballo Lake, New Mexico, in 2011 using an eddy covariance technique in the energy budget method. Total ET of 692 mm was measured during the growing season (n = 241 days) and 837 mm during the year. The American Society of Civil Engineers standardized ET for short crop (ETso) was calculated using climate data measured at the study site as 1560 mm during the growing season and 1870 mm during the year. Crop coefficients (ET/ETso) were fitted with a polynomial equation as a function of day of the year to develop saltgrass Kc function. A graphical and simplified method of computing Kc as a function of day of the year and crop season was also developed as an alternative method. Copyright © 2013 John Wiley & Sons, Ltd.     

Evaluating actual evapotranspiration and impacts of groundwater storage change in the North China Plain

As a critical water discharge term in basin‐scale water balance, accurate estimation of evapotranspiration (ET) is therefore important for sustainable water resources management. The understanding of the relationship between ET and groundwater storage change can improve our knowledge on the hydrological cycle in such regions with intensive agricultural land usage. Since the 1960s, the North China Plain (NCP) has experienced groundwater depletion because of overexploitation of groundwater for agriculture and urban development. Using meteorological data from 23 stations, the complementary relationship areal evapotranspiration model is evaluated against estimates of ET derived from regional water balance in the NCP during the period 1993–2008. The discrepancies between calculated ET and that derived by basin water balance indicate seasonal and interannual variations in model parameters. The monthly actual ET variations during the period from 1960 to 2008 are investigated by the calibrated model and then are used to derive groundwater storage change. The estimated actual ET is positively correlated with precipitation, and the general higher ET than precipitation indicates the contributions of groundwater irrigation to the total water supply. The long term decreasing trend in the actual ET can be explained by declining in precipitation, sunshine duration and wind speed. Over the past ~50 years, the calculated average annual water storage change, represented by the difference between actual ET and precipitation, was approximately 36 mm, or 4.8 km³; and the cumulative groundwater storage depletion was approximately 1700 mm, or 220 km³ in the NCP. The significantly groundwater storage depletion conversely affects the seasonal and interannual variations of ET. Irrigation especially during spring cause a marked increase in seasonal ET, whereas the rapid increasing of agricultural coverage over the NCP reduces the annual ET and is the primary control factor of the strong linear relationship between actual and potential ET. Copyright © 2013 John Wiley & Sons, Ltd.     

On Creating Global Gridded Terrestrial Water Budget Estimates from Satellite Remote Sensing

The increasing availability and reliability of satellite remote sensing products [e.g., precipitation (P), evapotranspiration (ET), and the total water storage change (TWSC)] make it feasible to estimate the global terrestrial water budget at fine spatial resolution. In this study, we start from a reference water budget dataset that combines all available data sources, including satellite remote sensing, land surface model (LSM) and reanalysis, and investigate the roles of different non-satellite remote sensing products in closing the terrestrial water budget through a sensitivity analysis by removing/replacing one or more categories of products during the budget estimation. We also study the differences made by various satellite products for the same budget variable. We find that the gradual removal of non-satellite data sources will generally worsen the closure errors in the budget estimates, and remote sensing retrievals of P, ET, and TWSC together with runoff (R) from LSM give the worst closure errors. The gauge-corrected satellite precipitation helps to improve the budget closure (4.2–9 % non-closure errors of annual mean precipitation) against using the non-gauge-corrected precipitation (7.6–10.4 % non-closure errors). At last, a data assimilation technique, the constrained Kalman filter, is applied to enforce the water balance, and it is found that the satellite remote sensing products, though with worst closure, yield comparable budget estimates in the constrained system to the reference data. Overall, this study provides a first comparison between the water budget closure using the satellite remote sensing products and a full combination of remote sensing, LSM, and reanalysis products on a quasi-global basis. This study showcases the capability and potential of the satellite remote sensing in closing the terrestrial water budget at fine spatial resolution if properly constrained.