scholarly journals Multidataset Study of Optimal Parameter and Uncertainty Estimation of a Land Surface Model with Bayesian Stochastic Inversion and Multicriteria Method

2004 ◽  
Vol 43 (10) ◽  
pp. 1477-1497 ◽  
Author(s):  
Youlong Xia ◽  
Mrinal K. Sen ◽  
Charles S. Jackson ◽  
Paul L. Stoffa

Abstract This study evaluates the ability of Bayesian stochastic inversion (BSI) and multicriteria (MC) methods to search for the optimal parameter sets of the Chameleon Surface Model (CHASM) using prescribed forcing to simulate observed sensible and latent heat fluxes from seven measurement sites representative of six biomes including temperate coniferous forests, tropical forests, temperate and tropical grasslands, temperate crops, and semiarid grasslands. Calibration results with the BSI and MC show that estimated optimal values are very similar for the important parameters that are specific to the CHASM model. The model simulations based on estimated optimal parameter sets perform much better than the default parameter sets. Cross-validations for two tropical forest sites show that the calibrated parameters for one site can be transferred to another site within the same biome. The uncertainties of optimal parameters are obtained through BSI, which estimates a multidimensional posterior probability density function (PPD). Marginal PPD analyses show that nonoptimal choices of stomatal resistance would contribute most to model simulation errors at all sites, followed by ground and vegetation roughness length at six of seven sites. The impact of initial root-zone soil moisture and nonmosaic approach on estimation of optimal parameters and their uncertainties is discussed.

2014 ◽  
Vol 7 (5) ◽  
pp. 6773-6809
Author(s):  
T. Osborne ◽  
J. Gornall ◽  
J. Hooker ◽  
K. Williams ◽  
A. Wiltshire ◽  
...  

Abstract. Studies of climate change impacts on the terrestrial biosphere have been completed without recognition of the integrated nature of the biosphere. Improved assessment of the impacts of climate change on food and water security requires the development and use of models not only representing each component but also their interactions. To meet this requirement the Joint UK Land Environment Simulator (JULES) land surface model has been modified to include a generic parametrisation of annual crops. The new model, JULES-crop, is described and evaluation at global and site levels for the four globally important crops; wheat, soy bean, maize and rice is presented. JULES-crop demonstrates skill in simulating the inter-annual variations of yield for maize and soy bean at the global level, and for wheat for major spring wheat producing countries. The impact of the new parametrisation, compared to the standard configuration, on the simulation of surface heat fluxes is largely an alteration of the partitioning between latent and sensible heat fluxes during the later part of the growing season. Further evaluation at the site level shows the model captures the seasonality of leaf area index and canopy height better than in standard JULES. However, this does not lead to an improvement in the simulation of sensible and latent heat fluxes. The performance of JULES-crop from both an earth system and crop yield model perspective is encouraging however, more effort is needed to develop the parameterisation of the model for specific applications. Key future model developments identified include the specification of the yield gap to enable better representation of the spatial variability in yield.


Author(s):  
N. J. Steinert ◽  
J. F. González-Rouco ◽  
P. de Vrese ◽  
E. García-Bustamante ◽  
S. Hagemann ◽  
...  

AbstractThe impact of various modifications of the JSBACH Land Surface Model to represent soil temperature and cold-region hydro-thermodynamic processes in climate projections of the 21st century is examined. We explore the sensitivity of JSBACH to changes in the soil thermodynamics, energy balance and storage, and the effect of including freezing and thawing processes. The changes involve 1) the net effect of an improved soil physical representation and 2) the sensitivity of our results to changed soil parameter values and their contribution to the simulation of soil temperatures and soil moisture, both aspects being presented in the frame of an increased bottom boundary depth from 9.83 m to 1418.84 m. The implementation of water phase changes and supercooled water in the ground creates a coupling between the soil thermal and hydrological regimes through latent heat exchange. Momentous effects on subsurface temperature of up to ±3 K, together with soil drying in the high northern latitudes, can be found at regional scales when applying improved hydro-thermodynamic soil physics. The sensitivity of the model to different soil parameter datasets occurs to be low but shows important implications for the root zone soil moisture content. The evolution of permafrost under pre-industrial forcing conditions emerges in simulated trajectories of stable states that differ by 4 – 6 • 106 km2 and shows large differences in the spatial extent of 105 –106 km2 by 2100, depending on the model configuration.


2011 ◽  
Vol 15 (2) ◽  
pp. 647-666 ◽  
Author(s):  
C. Szczypta ◽  
J.-C. Calvet ◽  
C. Albergel ◽  
G. Balsamo ◽  
S. Boussetta ◽  
...  

Abstract. An evaluation of the global ECMWF atmospheric reanalysis ERA-Interim (with a 0.5° grid) is performed over France, based on the high resolution (8 km) SAFRAN atmospheric reanalysis. The ERA-Interim precipitation, Incoming Solar Radiation (ISR), air temperature, air humidity, and wind speed, are compared with their SAFRAN counterparts. Also, interpolated in situ ISR observations are used in order to consolidate the evaluation of this variable. The daily precipitation estimates produced by ERA-Interim over France correlate very well with SAFRAN. However, the values are underestimated by 27%. A GPCP-corrected version of ERA-Interim is less biased (13%). The ERA-Interim estimates of ISR correlate very well with SAFRAN and with in situ observations on a daily basis. Whereas SAFRAN underestimates the ISR by 6 Wm−2, ERA-Interim overestimates the ISR by 10 Wm−2. In order to assess the impact of the ERA-Interim errors, simulations of the ISBA-A-gs land surface model are performed over the SMOSREX grassland site in southwestern France using ERA-Interim (with and without GPCP rescaling) and SAFRAN. Latent and sensible heat fluxes are simulated, together with carbon dioxide fluxes. The rescaled ERA-Interim performs better than the original ERA-Interim and permits to achieve flux scores similar to those obtained with SAFRAN.


2009 ◽  
Vol 9 (1) ◽  
pp. 2319-2380 ◽  
Author(s):  
A. de Meij ◽  
A. Gzella ◽  
P. Thunis ◽  
C. Cuvelier ◽  
B. Bessagnet ◽  
...  

Abstract. The objective of this study is to evaluate the impact of meteorological input data on calculated gas and aerosol concentrations. We use two different meteorological models (MM5 and WRF) together with the chemistry transport model CHIMERE. We focus on the Po valley area (Italy) for January and June 2005. Firstly we evaluate the meteorological parameters with observations. The analysis shows that the performance of both models is similar, however some small differences are still noticeable. Secondly, we analyze the impact of using MM5 and WRF on calculated PM10 and O3 concentrations. In general CHIMERE/MM5 and CHIMERE/WRF underestimate the PM10 concentrations for January. The difference in PM10 concentrations for January between CHIMERE/MM5 and CHIMERE/WRF is around a factor 1.6 (PM10 higher for CHIMERE/MM5). This difference and the larger underestimation in PM10 concentrations by CHIMERE/WRF are related to the differences in heat fluxes and the resulting PBL heights calculated by WRF. In general the PBL height by WRF meteorology is a factor 2.8 higher at noon in January than calculated by MM5. This study showed that the difference in microphysics scheme has an impact on the profile of cloud liquid water (CLW) calculated by the meteorological driver and therefore on the production of SO4 aerosol. A sensitivity analysis shows that changing the Noah Land Surface Model (LSM) for the 5-layer soil temperature model, the calculated monthly mean PM10 concentrations increase by 30%, due to the change in the heat fluxes and the resulting PBL heights. For June, PM10 calculated concentrations by CHIMERE/MM5 and CHIMERE/WRF are similar and agree with the observations. Calculated O3 values for June are in general overestimated by a factor 1.3 by CHIMERE/MM5 and CHIMRE/WRF. The reason for this is that daytime NO2 concentrations are a higher than the observations and nighttime NO concentrations (titration effect) are underestimated.


2015 ◽  
Vol 8 (4) ◽  
pp. 1139-1155 ◽  
Author(s):  
T. Osborne ◽  
J. Gornall ◽  
J. Hooker ◽  
K. Williams ◽  
A. Wiltshire ◽  
...  

Abstract. Studies of climate change impacts on the terrestrial biosphere have been completed without recognition of the integrated nature of the biosphere. Improved assessment of the impacts of climate change on food and water security requires the development and use of models not only representing each component but also their interactions. To meet this requirement the Joint UK Land Environment Simulator (JULES) land surface model has been modified to include a generic parametrisation of annual crops. The new model, JULES-crop, is described and evaluation at global and site levels for the four globally important crops; wheat, soybean, maize and rice. JULES-crop demonstrates skill in simulating the inter-annual variations of yield for maize and soybean at the global and country levels, and for wheat for major spring wheat producing countries. The impact of the new parametrisation, compared to the standard configuration, on the simulation of surface heat fluxes is largely an alteration of the partitioning between latent and sensible heat fluxes during the later part of the growing season. Further evaluation at the site level shows the model captures the seasonality of leaf area index, gross primary production and canopy height better than in the standard JULES. However, this does not lead to an improvement in the simulation of sensible and latent heat fluxes. The performance of JULES-crop from both an Earth system and crop yield model perspective is encouraging. However, more effort is needed to develop the parametrisation of the model for specific applications. Key future model developments identified include the introduction of processes such as irrigation and nitrogen limitation which will enable better representation of the spatial variability in yield.


2010 ◽  
Vol 7 (5) ◽  
pp. 7151-7190
Author(s):  
C. Szczypta ◽  
J.-C. Calvet ◽  
C. Albergel ◽  
G. Balsamo ◽  
S. Boussetta ◽  
...  

Abstract. An evaluation of the global ECMWF atmospheric reanalysis ERA-Interim (with a 0.5° grid) is performed over France, based on the high resolution (8 km) SAFRAN atmospheric reanalysis. The ERA-Interim precipitation, Incoming Solar Radiation (ISR), air temperature, air humidity, and wind speed, are compared with their SAFRAN counterparts. Also, interpolated in situ ISR observations are used in order to consolidate the evaluation of this variable. The daily precipitation estimates produced by ERA-Interim over France correlate very well with SAFRAN. However, the values are underestimated by 26%. A GPCP-corrected version of ERA-Interim is less biased (10–15%). The ERA-Interim estimates of ISR correlate very well with SAFRAN and with in situ observations on a daily basis. Whereas SAFRAN underestimates the ISR by 6–8 W m−2, ERA-Interim overestimates the ISR by 9–10 W m−2. In order to assess the impact of the ERA-Interim errors, simulations of the ISBA-A-gs land surface model are performed over the SMOSREX grassland site in southwestern France using ERA-Interim (with and without GPCP rescaling) and SAFRAN. Latent and sensible heat fluxes are simulated, together with carbon dioxide fluxes. The rescaled ERA-Interim performs better than the original ERA-Interim and permits to achieve flux scores similar to those obtained with SAFRAN.


2009 ◽  
Vol 9 (17) ◽  
pp. 6611-6632 ◽  
Author(s):  
A. de Meij ◽  
A. Gzella ◽  
C. Cuvelier ◽  
P. Thunis ◽  
B. Bessagnet ◽  
...  

Abstract. The objective of this study is to evaluate the impact of meteorological input data on calculated gas and aerosol concentrations. We use two different meteorological models (MM5 and WRF) together with the chemistry transport model CHIMERE. We focus on the Po valley area (Italy) for January and June 2005. Firstly we evaluate the meteorological parameters with observations. The analysis shows that the performance of both models in calculating surface parameters is similar, however differences are still observed. Secondly, we analyze the impact of using MM5 and WRF on calculated PM10 and O3 concentrations. In general CHIMERE/MM5 and CHIMERE/WRF underestimate the PMv concentrations for January. The difference in PM10 concentrations for January between CHIMERE/MM5 and CHIMERE/WRF is around a factor 1.6 (PM10 higher for CHIMERE/MM5). This difference and the larger underestimation in PM10 concentrations by CHIMERE/WRF are related to the differences in heat fluxes and the resulting PBL heights calculated by WRF. In general the PBL height by WRF meteorology is a factor 2.8 higher at noon in January than calculated by MM5. This study showed that the difference in microphysics scheme has an impact on the profile of cloud liquid water (CLW) calculated by the meteorological driver and therefore on the production of SO4 aerosol. A sensitivity analysis shows that changing the Noah Land Surface Model (LSM) in our WRF pre-processing for the 5-layer soil temperature model, calculated monthly mean PMv concentrations increase by 30%, due to the change in the heat fluxes and the resulting PBL heights. For June, PM10 calculated concentrations by CHIMERE/MM5 and CHIMERE/WRF are similar and agree with the observations. Calculated O3 values for June are in general overestimated by a factor 1.3 by CHIMERE/MM5 and CHIMERE/WRF. High temporal correlations are found between modeled and observed O3 concentrations.


2016 ◽  
Vol 16 (13) ◽  
pp. 8375-8387 ◽  
Author(s):  
Liang Chen ◽  
Yanping Li ◽  
Fei Chen ◽  
Alan Barr ◽  
Michael Barlage ◽  
...  

Abstract. A thick top layer of organic matter is a dominant feature in boreal forests and can impact land–atmosphere interactions. In this study, the multi-parameterization version of the Noah land surface model (Noah-MP) was used to investigate the impact of incorporating a forest-floor organic soil layer on the simulated surface energy and water cycle components at the BERMS Old Aspen site (OAS) field station in central Saskatchewan, Canada. Compared to a simulation without an organic soil parameterization (CTL), the Noah-MP simulation with an organic soil (OGN) improved Noah-MP-simulated soil temperature profiles and soil moisture at 40–100 cm, especially the phase and amplitude (Seasonal cycle) of soil temperature below 10 cm. OGN also enhanced the simulation of sensible and latent heat fluxes in spring, especially in wet years, which is mostly related to the timing of spring soil thaw and warming. Simulated top-layer soil moisture is better in OGN than that in CTL. The effects of including an organic soil layer on soil temperature are not uniform throughout the soil depth and are more prominent in summer. For drought years, the OGN simulation substantially modified the partitioning of water between direct soil evaporation and vegetation transpiration. For wet years, the OGN-simulated latent heat fluxes are similar to CTL except for the spring season when OGN produced less evaporation, which was closer to observations. Including organic soil produced more subsurface runoff and resulted in much higher runoff throughout the freezing periods in wet years.


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.


2015 ◽  
Vol 8 (6) ◽  
pp. 1857-1876 ◽  
Author(s):  
J. J. Guerrette ◽  
D. K. Henze

Abstract. Here we present the online meteorology and chemistry adjoint and tangent linear model, WRFPLUS-Chem (Weather Research and Forecasting plus chemistry), which incorporates modules to treat boundary layer mixing, emission, aging, dry deposition, and advection of black carbon aerosol. We also develop land surface and surface layer adjoints to account for coupling between radiation and vertical mixing. Model performance is verified against finite difference derivative approximations. A second-order checkpointing scheme is created to reduce computational costs and enable simulations longer than 6 h. The adjoint is coupled to WRFDA-Chem, in order to conduct a sensitivity study of anthropogenic and biomass burning sources throughout California during the 2008 Arctic Research of the Composition of the Troposphere from Aircraft and Satellites (ARCTAS) field campaign. A cost-function weighting scheme was devised to reduce the impact of statistically insignificant residual errors in future inverse modeling studies. Results of the sensitivity study show that, for this domain and time period, anthropogenic emissions are overpredicted, while wildfire emission error signs vary spatially. We consider the diurnal variation in emission sensitivities to determine at what time sources should be scaled up or down. Also, adjoint sensitivities for two choices of land surface model (LSM) indicate that emission inversion results would be sensitive to forward model configuration. The tools described here are the first step in conducting four-dimensional variational data assimilation in a coupled meteorology–chemistry model, which will potentially provide new constraints on aerosol precursor emissions and their distributions. Such analyses will be invaluable to assessments of particulate matter health and climate impacts.


Sign in / Sign up

Export Citation Format

Share Document