Premium
Modeling soil evaporation efficiency in a range of soil and atmospheric conditions using a meta‐analysis approach
Author(s) -
Merlin O.,
Stefan V. G.,
Amazirh A.,
Chanzy A.,
Ceschia E.,
ErRaki S.,
Gentine P.,
Tallec T.,
Ezzahar J.,
Bircher S.,
Beringer J.,
Khabba S.
Publication year - 2016
Publication title -
water resources research
Language(s) - English
Resource type - Journals
SCImago Journal Rank - 1.863
H-Index - 217
eISSN - 1944-7973
pISSN - 0043-1397
DOI - 10.1002/2015wr018233
Subject(s) - soil texture , soil science , environmental science , range (aeronautics) , forcing (mathematics) , water content , evaporation , mathematics , soil water , atmospheric sciences , meteorology , geotechnical engineering , geology , materials science , geography , composite material
A meta‐analysis data‐driven approach is developed to represent the soil evaporative efficiency (SEE) defined as the ratio of actual to potential soil evaporation. The new model is tested across a bare soil database composed of more than 30 sites around the world, a clay fraction range of 0.02–0.56, a sand fraction range of 0.05–0.92, and about 30,000 acquisition times. SEE is modeled using a soil resistance ( r ss ) formulation based on surface soil moisture ( θ ) and two resistance parametersr s s , r e fand θ efolding . The data‐driven approach aims to express both parameters as a function of observable data including meteorological forcing, cut‐off soil moisture valueθ 1 / 2at which SEE=0.5, and first derivative of SEE atθ 1 / 2, named Δ θ 1 / 2 − 1. An analytical relationship between ( r s s , r e f ; θ e f o l d i n g ) and ( θ 1 / 2 ; Δ θ 1 / 2 − 1 ) is first built by running a soil energy balance model for two extreme conditions with r ss = 0 andr s s ∼ ∞ using meteorological forcing solely, and by approaching the middle point from the two (wet and dry) reference points. Two different methods are then investigated to estimate the pair ( θ 1 / 2 ; Δ θ 1 / 2 − 1 ) either from the time series of SEE and θ observations for a given site, or using the soil texture information for all sites. The first method is based on an algorithm specifically designed to accomodate for strongly nonlinear SEE ( θ ) relationships and potentially large random deviations of observed SEE from the mean observed SEE ( θ ) . The second method parameterizesθ 1 / 2as a multi‐linear regression of clay and sand percentages, and sets Δ θ 1 / 2 − 1to a constant mean value for all sites. The new model significantly outperformed the evaporation modules of ISBA (Interaction Sol‐Biosphère‐Atmosphère), H‐TESSEL (Hydrology‐Tiled ECMWF Scheme for Surface Exchange over Land), and CLM (Community Land Model). It has potential for integration in various land‐surface schemes, and real calibration capabilities using combined thermal and microwave remote sensing data.