The GEWEX LandFlux project: evaluation of model evaporation using tower-based and globally-gridded forcing data.

Determining the spatial distribution and temporal development of evaporation at regional and global scales is required to improve our understanding of the coupled water and energy cycles and to better monitor any changes in observed trends and variability of linked hydrological processes. With recent international efforts guiding the development of long-term and globally distributed flux estimates, continued product assessments are required to inform upon the selection of suitable model structures and also to establish the appropriateness of these multi-model simulations for global application. In support of the objectives of the GEWEX LandFlux project, four commonly used evaporation models are evaluated against data from tower-based eddy-covariance observations, distributed across a range of biomes and climate zones. The selected schemes include the Surface Energy Balance System (SEBS) approach, the Priestley-Taylor Jet Propulsion Laboratory (PT-JPL) model, the Penman-Monteith based Mu model (PMMu) and the Global Land Evaporation: the Amsterdam Methodology (GLEAM). Here we seek to examine the fidelity of global evaporation simulations by examining the multimodel response to varying sources of forcing data. To do this, we perform parallel and collocated model simulations using tower-based data together with a global-scale gridbased forcing product. Through quantifying the multi-model response to high-quality tower data, a better understanding of the subsequent model response to coarse-scale globally gridded data that underlies the LandFlux product can be obtained, while also providing a relative evaluation and assessment of model performance. Using surface flux observations from forty-five globally distributed eddy-covariance stations as independent metrics of performance, the tower-based analysis indicated that PT-JPL provided the highest overally statistical performance (0.72; 61Wm^-2; 0.65), followed closely by GLEAM (0.68; 64Wm^-2; 0.62), with values in parenthesis representing the R², RMSD and Nash-Sutclifie Effciency (NSE) and respectively. PM-Mu (0.51; 78Wm^-2; 0.45) tended to underestimate fluxes, while SEBS (0.72; 101Wm^-2; 0.24) overestimated values relative to observations. A focused analysis across specific biome types and climate zones showed considerable variability in the performance of all models, with no single model consistently able to outperform any other. Results also indicated that the global gridded data tended to reduce the performance for all of the studied models when compared to the tower data, likely a response to scale mismatch and issues related to forcing quality. Rather than relying on any single model simulation, the spatial and temporal variability at both the towerand grid-scale highlighted the potential benefits of developing an ensemble or blended evaporation product for global scale LandFlux applications. Challenges related to the robust assessment of the LandFlux product are also discussed. [ABSTRACT FROM AUTHOR]