What’s the impact of improved soil representations in the ECMWF land surface model and how does it affect the extremes?

Author(s):  
Souhail Boussetta ◽  
Gabriele Arduini ◽  
Gianpaolo Balsamo ◽  
Emanuel Dutra ◽  
Anna Agusti-Panareda ◽  
...  

<p>With increasingly higher spatial resolution and a broader applications, the importance of soil representation (e.g. soil depth, vertical discretisation, vegetation rooting) within land surface models is enhanced. Those modelling choices actually affects the way land surfaces store and regulate water, energy and also carbon fluxes. Heat and water vapour fluxes towards the atmosphere and deeper soil, exhibit variations spanning a range of time scales from minutes to months in the coupled land-atmosphere system. This is further modulated by the vertical roots' distribution, and soil moisture stress function, which control evapotranspiration under soil moisture stress conditions. Currently in the ECMWF land Surface Scheme the soil column is represented by a fixed 4 layers configuration with a total of approximately 3m depth.</p><p>In the present study we explore new configurations with increased soil depth (up to 8m) and higher vertical discretisation (up to 10 layers) including a dissociation between the treatment of water and heat fluxes. Associated with the soil vertical resolution, the vertical distribution of roots is also investigated. A new scheme that assumes a uniform root distribution with an associated maximum rooting depth is explored. The impact of these new configurations is assessed through surface offline simulations driven by the ERA5 meteorological forcing against in-situ and global products of energy, water and carbon fluxes with a particular focus on the diurnal cycle and extreme events in recent years.</p>

2020 ◽  
Author(s):  
Anna B. Harper ◽  
Karina E. Williams ◽  
Patrick C. McGuire ◽  
Maria Carolina Duran Rojas ◽  
Debbie Hemming ◽  
...  

Abstract. Drought is predicted to increase in the future due to climate change, bringing with it a myriad of impacts on ecosystems. Plants respond to drier soils by reducing stomatal conductance, in order to conserve water and avoid hydraulic damage. Despite the importance of plant drought responses for the global carbon cycle and local/regional climate feedbacks, land surface models are unable to capture observed plant responses to soil moisture stress. We assessed the impact of soil moisture stress on simulated gross primary productivity (GPP) and latent energy flux (LE) in the Joint UK Land Environment Simulator (JULES) vn4.9 on seasonal and annual timescales, and evaluated ten different representations of stress in the model. For the default configuration, GPP was more realistic in temperate biome sites than in the tropics or high latitudes/cold region sites, while LE was best simulated in temperate and high latitude/cold sites. Errors not due to soil moisture stress, possibly linked to phenology, contributed to model biases for GPP in tropical savannah and deciduous forest sites. We found that three alternative approaches to calculating soil moisture stress produced more realistic results than the default parameterization for most biomes and climates. All of these involved increasing the number of soil layers from 4 to 14, and the soil depth from 3m to 10.8m. In addition, we found improvements when soil matric potential replaced volumetric water content in the stress equation, when the onset of stress was delayed, and when roots extended deeper into the soil. For LE, the biases were highest in the default configuration in temperate mixed forests, with overestimation occurring during most of the year. At these sites, reducing soil moisture stress (with the new parameterizations mentioned above) increased LE and made the simulation worse. Further evaluation into the reason for the high bias in LE at many of the sites would enable improvements in both carbon and energy fluxes with new parameterizations for soil moisture stress.


2021 ◽  
Vol 14 (6) ◽  
pp. 3269-3294
Author(s):  
Anna B. Harper ◽  
Karina E. Williams ◽  
Patrick C. McGuire ◽  
Maria Carolina Duran Rojas ◽  
Debbie Hemming ◽  
...  

Abstract. Drought is predicted to increase in the future due to climate change, bringing with it myriad impacts on ecosystems. Plants respond to drier soils by reducing stomatal conductance in order to conserve water and avoid hydraulic damage. Despite the importance of plant drought responses for the global carbon cycle and local and regional climate feedbacks, land surface models are unable to capture observed plant responses to soil moisture stress. We assessed the impact of soil moisture stress on simulated gross primary productivity (GPP) and latent energy flux (LE) in the Joint UK Land Environment Simulator (JULES) vn4.9 on seasonal and annual timescales and evaluated 10 different representations of soil moisture stress in the model. For the default configuration, GPP was more realistic in temperate biome sites than in the tropics or high-latitude (cold-region) sites, while LE was best simulated in temperate and high-latitude (cold) sites. Errors that were not due to soil moisture stress, possibly linked to phenology, contributed to model biases for GPP in tropical savanna and deciduous forest sites. We found that three alternative approaches to calculating soil moisture stress produced more realistic results than the default parameterization for most biomes and climates. All of these involved increasing the number of soil layers from 4 to 14 and the soil depth from 3.0 to 10.8 m. In addition, we found improvements when soil matric potential replaced volumetric water content in the stress equation (the “soil14_psi” experiments), when the critical threshold value for inducing soil moisture stress was reduced (“soil14_p0”), and when plants were able to access soil moisture in deeper soil layers (“soil14_dr*2”). For LE, the biases were highest in the default configuration in temperate mixed forests, with overestimation occurring during most of the year. At these sites, reducing soil moisture stress (with the new parameterizations mentioned above) increased LE and increased model biases but improved the simulated seasonal cycle and brought the monthly variance closer to the measured variance of LE. Further evaluation of the reason for the high bias in LE at many of the sites would enable improvements in both carbon and energy fluxes with new parameterizations for soil moisture stress. Increasing the soil depth and plant access to deep soil moisture improved many aspects of the simulations, and we recommend these settings in future work using JULES or as a general way to improve land surface carbon and water fluxes in other models. In addition, using soil matric potential presents the opportunity to include plant functional type-specific parameters to further improve modeled fluxes.


2007 ◽  
Vol 8 (1) ◽  
pp. 68-87 ◽  
Author(s):  
Margaret A. LeMone ◽  
Fei Chen ◽  
Joseph G. Alfieri ◽  
Mukul Tewari ◽  
Bart Geerts ◽  
...  

Abstract Analyses of daytime fair-weather aircraft and surface-flux tower data from the May–June 2002 International H2O Project (IHOP_2002) and the April–May 1997 Cooperative Atmosphere Surface Exchange Study (CASES-97) are used to document the role of vegetation, soil moisture, and terrain in determining the horizontal variability of latent heat LE and sensible heat H along a 46-km flight track in southeast Kansas. Combining the two field experiments clearly reveals the strong influence of vegetation cover, with H maxima over sparse/dormant vegetation, and H minima over green vegetation; and, to a lesser extent, LE maxima over green vegetation, and LE minima over sparse/dormant vegetation. If the small number of cases is producing the correct trend, other effects of vegetation and the impact of soil moisture emerge through examining the slope ΔxyLE/ΔxyH for the best-fit straight line for plots of time-averaged LE as a function of time-averaged H over the area. Based on the surface energy balance, H + LE = Rnet − Gsfc, where Rnet is the net radiation and Gsfc is the flux into the soil; Rnet − Gsfc ∼ constant over the area implies an approximately −1 slope. Right after rainfall, H and LE vary too little horizontally to define a slope. After sufficient drying to produce enough horizontal variation to define a slope, a steep (∼−2) slope emerges. The slope becomes shallower and better defined with time as H and LE horizontal variability increases. Similarly, the slope becomes more negative with moister soils. In addition, the slope can change with time of day due to phase differences in H and LE. These trends are based on land surface model (LSM) runs and observations collected under nearly clear skies; the vegetation is unstressed for the days examined. LSM runs suggest terrain may also play a role, but observational support is weak.


Author(s):  
Nemesio Rodriguez-Fernandez ◽  
Patricia de Rosnay ◽  
Clement Albergel ◽  
Philippe Richaume ◽  
Filipe Aires ◽  
...  

The assimilation of Soil Moisture and Ocean Salinity (SMOS) data into the ECMWF (European Centre for Medium Range Weather Forecasts) H-TESSEL (Hydrology revised - Tiled ECMWF Scheme for Surface Exchanges over Land) model is presented. SMOS soil moisture (SM) estimates have been produced specifically by training a neural network with SMOS brightness temperatures as input and H-TESSEL model SM simulations as reference. This can help the assimilation of SMOS information in several ways: (1) the neural network soil moisture (NNSM) data have a similar climatology to the model, (2) no global bias is present with respect to the model even if regional differences can exist. Experiments performing joint data assimilation (DA) of NNSM, 2 metre air temperature and relative humidity or NNSM-only DA are discussed. The resulting SM was evaluated against a large number of in situ measurements of SM obtaining similar results to those of the model with no assimilation, even if significant differences were found from site to site. In addition, atmospheric forecasts initialized with H-TESSEL runs (without DA) or with the analysed SM were compared to measure of the impact of the satellite information. Although, NNSM DA has an overall neutral impact in the forecast in the Tropics, a significant positive impact was found in other areas and periods, especially in regions with limited in situ information. The joint NNSM, T2m and RH2m DA improves the forecast for all the seasons in the Southern Hemisphere. The impact is mostly due to T2m and RH2m, but SMOS NN DA alone also improves the forecast in July- September. In the Northern Hemisphere, the joint NNSM, T2m and RH2m DA improves the forecast in April-September, while NNSM alone has a significant positive effect in July-September. Furthermore, forecasting skill maps show that SMOS NNSM improves the forecast in North America and in Northern Asia for up to 72 hours lead time.


2021 ◽  
Author(s):  
Jaime Gaona ◽  
Pere Quintana-Seguí ◽  
Maria José Escorihuela

<p>Droughts in the Iberian Peninsula are a natural hazard of great relevance due to their recurrence, severity and impact on multiple environmental and socioeconomic aspects. The Ebro Basin, located in the NE of the Iberian Peninsula, is particularly vulnerable to drought with consequences on agriculture, urban water supply and hydropower. This study, performed within the Project HUMID (CGL2017-85687-R), aims at evaluating the influence of the climatic, land cover and soil characteristics on the interactions between rainfall, evapotranspiration and soil moisture anomalies which define the spatio-temporal drought patterns in the basin.</p><p>The onset, propagation and mitigation of droughts in the Iberian Peninsula is driven by anomalies of rainfall, evapotranspiration and soil moisture, which are related by feedback processes. To test the relative importance of such anomalies, we evaluate the contribution of climatic, land-cover and geologic heterogeneity on the definition of the spatio-temporal patterns of drought. We use the Köppen-Geiger climatic classification to assess how the contrasting climatic types within the basin determine differences on drought behavior. Land-cover types that govern the partition between evaporation and transpiration are also of great interest to discern the influence of vegetation and crop types on the anomalies of evapotranspiration across the distinct regions of the basin (e.g. forested mountains vs. crop-dominated areas). The third physical characteristic whose effect on drought we investigate is the impact of soil properties on soil moisture anomalies.</p><p>The maps and time series used for the spatio-temporal analysis are based on drought indices calculated with high-resolution datasets from remote sensing (MOD16A2ET and SMOS1km) and the land-surface model SURFEX-ISBA. The Standardized Precipitation Index (SPI), the EvapoTranspiration Deficit Index (ETDI) and the Soil Moisture Deficit Index (SMDI) are the three indices chosen to characterize the anomalies of the corresponding rainfall (atmospheric), evapotranspiration (atmosphere-land interface) and soil moisture (land) anomalies (components of the water balance). The comparison of the correlations of the indices (with different time lags) between contrasting regions offers insights about the impact of climate, land-cover and soil properties in the dominance, the timing of the response and memory aspects of the interactions. The high spatial and temporal resolution of remote sensing and land-surface model data allows adopting time and spatial scales suitable to investigate the influence of these physical factors with detail beyond comparison with ground-based datasets.</p><p>The spatial and temporal analysis prove useful to investigate the physical factors of influence on the anomalies between rainfall, evapotranspiration and soil moisture. This approach facilitates the physical interpretation of the anomalies of drought indices aiming to improve the characterization of drought in heterogeneous semi-arid areas like the Ebro River Basin.</p>


2006 ◽  
Vol 7 (1) ◽  
pp. 61-80 ◽  
Author(s):  
B. Decharme ◽  
H. Douville ◽  
A. Boone ◽  
F. Habets ◽  
J. Noilhan

Abstract This study focuses on the influence of an exponential profile of saturated hydraulic conductivity, ksat, with soil depth on the water budget simulated by the Interaction Soil Biosphere Atmosphere (ISBA) land surface model over the French Rhône River basin. With this exponential profile, the saturated hydraulic conductivity at the surface increases by approximately a factor of 10, and its mean value increases in the root zone and decreases in the deeper region of the soil in comparison with the values given by Clapp and Hornberger. This new version of ISBA is compared to the original version in offline simulations using the Rhône-Aggregation high-resolution database. Low-resolution simulations, where all atmospheric data and surface parameters have been aggregated, are also performed to test the impact of the modified ksat profile at the typical scale of a climate model. The simulated discharges are compared to observations from a dense network consisting of 88 gauging stations. Results of the high-resolution experiments show that the exponential profile of ksat globally improves the simulated discharges and that the assumption of an increase in saturated hydraulic conductivity from the soil surface to a depth close to the rooting depth in comparison with values given by Clapp and Hornberger is reasonable. Results of the scaling experiments indicate that this parameterization is also suitable for large-scale hydrological applications. Nevertheless, low-resolution simulations with both model versions overestimate evapotranspiration (especially from the plant transpiration and the wet fraction of the canopy) to the detriment of total runoff, which emphasizes the need for implementing subgrid distribution of precipitation and land surface properties in large-scale hydrological applications.


2009 ◽  
Vol 10 (6) ◽  
pp. 1548-1560 ◽  
Author(s):  
Richard J. Ellis ◽  
Christopher M. Taylor ◽  
Graham P. Weedon ◽  
Nicola Gedney ◽  
Douglas B. Clark ◽  
...  

Abstract A critical function of a land surface scheme, used in climate and weather prediction models, is to partition the energy from insolation into sensible and latent heat fluxes. Many use a soil moisture function to control the surface moisture fluxes through the transpiration. The validity and global distribution of the parameters used to calculate this soil moisture stress function are difficult to assess. This work presents a method to map soil moisture stress globally from an earth observation vegetation index and precipitation data, and it compares the resulting distributions with output from the Joint U.K. Land Environment Simulator (JULES) land surface scheme. A number of model runs with different soil and vegetation parameters are compared. These examine the sensitivity of the seasonality of soil moisture stress, within the model, to the parameterization of soil hydraulic properties and the seasonality of leaf area index in the vegetation. It is found that the seasonality of soil moisture within the model is more sensitive to the soil hydraulic properties than the leaf area index. The partitioning of throughfall into evaporation and runoff, in the model, is the dominant factor in determining the timing of soil moisture stress.


2020 ◽  
Author(s):  
Leqiang Sun ◽  
Stéphane Belair ◽  
Marco Carrera ◽  
Bernard Bilodeau

<p>Canadian Space Agency (CSA) has recently started receiving and processing the images from the recently launched C-band RADARSAT Constellation Mission (RCM). The backscatter and soil moisture retrievals products from the previously launched RADARSAT-2 agree well with both in-situ measurements and surface soil moisture modeled with land surface model Soil, Vegetation, and Snow (SVS). RCM will provide those products at an even better spatial coverage and temporal resolution. In preparation of the potential operational application of RCM products in Canadian Meteorological Center (CMC), this paper presents the scenarios of assimilating either soil moisture retrieval or outright backscatter signal in a 100-meter resolution version of the Canadian Land Data Assimilation System (CaLDAS) on field scale with time interval of three hours. The soil moisture retrieval map was synthesized by extrapolating the regression relationship between in-situ measurements and open loop model output based on soil texture lookup table. Based on this, the backscatter map was then generated with the surface roughness retrieved from RADARSAT-2 images using a modified Integral Equation Model (IEM) model. Bias correction was applied to the Ensemble Kalman filter (EnKF) to mitigate the impact of nonlinear errors introduced by multi-sourced perturbations. Initial results show that the assimilation of backscatter is as effective as assimilating soil moisture retrievals. Compared to open loop, both can improve the analysis of surface moisture, particularly in terms of reducing bias.  </p>


2020 ◽  
Author(s):  
Daniel Abel ◽  
Felix Pollinger ◽  
Katrin Ziegler ◽  
Heiko Paeth

<p>The EFRE-Project BigData@Geo, founded by the European Union, aims to create highly resolved climate projections for the model region of Lower Franconia in Bavaria, Germany. These projections are analyzed and made available to local stakeholders of agriculture, forestry, and viniculture as well as the public. As recent regional climate models are not dealing with the necessary spatiotemporal resolution the model REMO will be developed in the project‘s frame in cooperation with the Climate Service Center Germany (GERICS).</p><p>For these very high resolutions, besides improvements like the non-hydrostatic atmosphere, higher resolved static land surface parameters, and land use land cover changes, etc., realistic modeling of the soil hydrology becomes absolutely necessary. Therefore, REMO is extended by a 5-layer soil scheme which is a first step to overcome restrictions of the recently used soil hydrology scheme due to the included vertical water flow. Furthermore, the current work also aims to implement lateral water flows between grid cells because this is the only way to model the soil hydrology appropriate to the project‘s question.</p><p>The current model version of REMO includes a bucket scheme that treats the soil hydrology as a single layer. The soil depth is equal to the rooting depth and, thus, depends on the overlying vegetation class. Consequently, the whole soil moisture of the soil is available for transpiration. Evaporation only occurs if the soil moisture reaches at least 90 % of the field capacity.</p><p>The 5-layer scheme has 5 layers with increasing thicknesses for deeper layers. The maximum depth of the soil is at approximately 10 m or the depth of the bedrock. Due to the existence of water below the rooting zone and the processes of capillary rise and percolation more water becomes available for transpiration compared to the bucket scheme. Furthermore, evaporation only occurs if the uppermost layer contains soil moisture which is a more realistic process representation as well.</p><p>First results of the comparison of the two schemes and with observation data in the EURO-CORDEX region and a german subregion are presented. We also show some sensitivity studies of the current improvements to the parameterizations of the 5-layer scheme which are necessary for the goal of incorporation of the lateral flow.</p>


Sign in / Sign up

Export Citation Format

Share Document