scholarly journals A vertical representation of soil carbon in the JULES land surface scheme (vn4.3_permafrost) with a focus on permafrost regions

2017 ◽  
Vol 10 (2) ◽  
pp. 959-975 ◽  
Author(s):  
Eleanor J. Burke ◽  
Sarah E. Chadburn ◽  
Altug Ekici

Abstract. An improved representation of the carbon cycle in permafrost regions will enable more realistic projections of the future climate–carbon system. Currently JULES (the Joint UK Land Environment Simulator) – the land surface model of the UK Earth System Model (UKESM) – uses the standard four-pool RothC soil carbon model. This paper describes a new version of JULES (vn4.3_permafrost) in which the soil vertical dimension is added to the soil carbon model, with a set of four pools in every soil layer. The respiration rate in each soil layer depends on the temperature and moisture conditions in that layer. Cryoturbation/bioturbation processes, which transfer soil carbon between layers, are represented by diffusive mixing. The litter inputs and the soil respiration are both parametrized to decrease with increasing depth. The model now includes a tracer so that selected soil carbon can be labelled and tracked through a simulation. Simulations show an improvement in the large-scale horizontal and vertical distribution of soil carbon over the standard version of JULES (vn4.3). Like the standard version of JULES, the vertically discretized model is still unable to simulate enough soil carbon in the tundra regions. This is in part because JULES underestimates the plant productivity over the tundra, but also because not all of the processes relevant for the accumulation of permafrost carbon, such as peat development, are included in the model. In comparison with the standard model, the vertically discretized model shows a delay in the onset of soil respiration in the spring, resulting in an increased net uptake of carbon during this time. In order to provide a more suitable representation of permafrost carbon for quantifying the permafrost carbon feedback within UKESM, the deep soil carbon in the permafrost region (below 1 m) was initialized using the observed soil carbon. There is now a slight drift in the soil carbon ( <  0.018 % decade−1), but the change in simulated soil carbon over the 20th century, when there is little climate change, is comparable to the original vertically discretized model and significantly larger than the drift.

2016 ◽  
Author(s):  
Eleanor J. Burke ◽  
Sarah E. Chadburn ◽  
Altug Ekici

Abstract. An improved representation of the carbon cycle in permafrost regions will enable more realistic projections of the future climate-carbon system. Currently JULES (the Joint UK Land Environment Simulator) – the land surface model of the UK Earth System Model (UKESM) – uses the standard 4-pool RothC soil carbon model. This paper describes a new version of JULES (vn4.3_permafrost) in which the soil vertical dimension is added to the soil carbon model, with a set of four pools in every soil layer. The respiration rate in each soil layer depends on the temperature and moisture conditions in that layer. Cryoturbation/bioturbation processes, which transfer soil carbon between layers, are represented by diffusive mixing. The litter inputs and the soil respiration are both parameterised to decrease with increasing depth. The model now includes a tracer so that selected soil carbon can be labelled and tracked through a simulation. Simulations show an improvement in the large-scale horizontal and vertical distribution of soil carbon over the standard version of JULES (vn4.3). Like the standard version of JULES, the vertically discretised model is still unable to simulate enough soil carbon in the tundra regions. This is in part because JULES underestimates the plant productivity over the tundra, but also because not all of the processes relevant for the accumulation of permafrost carbon are included in the model. In comparison with the standard model, the vertically discretised model shows a delay in the onset of soil respiration in the spring, resulting in an increased net uptake of carbon during this time. In order to provide a more suitable representation of permafrost carbon for quantifying the permafrost carbon feedback within UKESM, the deep soil carbon in the permafrost region (below 1 m) was initialised using the observed soil carbon. There is now a slight drift in the soil carbon (


2020 ◽  
Author(s):  
Yao Gao ◽  
Eleanor Burke ◽  
Sarah Chadburn ◽  
Maarit Raivonen ◽  
Timo Vesala ◽  
...  

&lt;p&gt;Atmospheric emissions and concentrations of CH&lt;sub&gt;4&lt;/sub&gt; are continuing to increase, making CH&lt;sub&gt;4&lt;/sub&gt; the second most important human-influenced greenhouse gas in terms of climate forcing, after CO&lt;sub&gt;2&lt;/sub&gt;. Previous studies indicated that wetland CH&lt;sub&gt;4&lt;/sub&gt; emission is not only the single largest but also the most uncertain natural source in the global CH&lt;sub&gt;4&lt;/sub&gt; budget. Furthermore, the strong sensitivity of wetland CH&lt;sub&gt;4&lt;/sub&gt; emissions to environmental conditions has raised concerns on potential positive feedbacks to climate change. Therefore, evaluation of the process-based land surface models of earth system models (ESMs) in simulating CH&lt;sub&gt;4&lt;/sub&gt; emission over wetlands is needed for more precise future predictions. In this work, a set of high-latitude wetland sites with various nutrient conditions are studied. The wetland CH&lt;sub&gt;4&lt;/sub&gt; fluxes are simulated by the land surface model JULES of the UK Earth System model and the Helsinki peatland methane emission model (HIMMELI), which is developed at Finnish Meteorological Institute and Helsinki University. The differences between the modelled and observed CH&lt;sub&gt;4&lt;/sub&gt; fluxes are analyzed, complemented with key environmental variables for interpretation (e.g. soil temperature and moisture, vegetation types, snow depth, NPP, soil carbon). In general, the simulated CH&lt;sub&gt;4&lt;/sub&gt; fluxes by HIMMELI is closer to the observed CH&lt;sub&gt;4&lt;/sub&gt; fluxes in magnitude and seasonality at sites than those by JULES. The inter-annual variability of simulated CH&lt;sub&gt;4&lt;/sub&gt; fluxes by HIMMELI depends on the simulated anoxic soil respiration, which serves as the substrate of the CH&lt;sub&gt;4&lt;/sub&gt; fluxes in HIMMELI. The anoxic soil respiration is calculated based on the simulated soil respiration and water table depth in JULES. More accurate simulation of soil carbon pool and water table depth in JULES will lead to improvement in the simulated anoxic soil respiration.&lt;/p&gt;


Forests ◽  
2021 ◽  
Vol 12 (1) ◽  
pp. 89
Author(s):  
Hong Wei ◽  
Xiuling Man

The change of litter input can affect soil respiration (Rs) by influencing the availability of soil organic carbon and nutrients, regulating soil microenvironments, thus resulting in a profound influence on soil carbon cycle of the forest ecosystem. We conducted an aboveground litterfall manipulation experiment in different-aged Betula platyphylla forests (25-, 40- and 61-year-old) of the permafrost region, located in the northeast of China, during May to October in 2018, with each stand treated with doubling litter (litter addition, DL), litter exclusion (no-litter, NL) and control litter (CK). Our results indicated that Rs decreased under NL treatment compared with CK treatment. The effect size lessened with the increase in the stand age; the greatest reduction was found for young Betula platyphylla forest (24.46% for 25-year-old stand) and tended to stabilize with the growth of forest with the reduction of 15.65% and 15.23% for 40-and 61- year-old stands, respectively. Meanwhile, under DL treatment, Rs increased by 27.38%, 23.83% and 23.58% on 25-, 40- and 61-year-old stands, respectively. Our results also showed that the increase caused by DL treatment was larger than the reduction caused by NL treatment, leading to a priming effect, especially on 40- and 61-year-old stands. The change in litter input was the principal factor affecting the change of Rs under litter manipulation. The soil temperature was also a main factor affecting the contribution rate of litter to Rs of different-aged stands, which had a significant positive exponential correlation with Rs. This suggests that there is a significant relationship between litter and Rs, which consequently influences the soil carbon cycle in Betula platyphylla forests of the permafrost region, Northeast China. Our finding indicated the increased litter enhanced the Rs in Betula platyphylla forest, which may consequently increase the carbon emission in a warming climate in the future. It is of great importance for future forest management in the permafrost region, Northeast China.


2017 ◽  
Vol 10 (5) ◽  
pp. 2031-2055 ◽  
Author(s):  
Thomas Schwitalla ◽  
Hans-Stefan Bauer ◽  
Volker Wulfmeyer ◽  
Kirsten Warrach-Sagi

Abstract. Increasing computational resources and the demands of impact modelers, stake holders, and society envision seasonal and climate simulations with the convection-permitting resolution. So far such a resolution is only achieved with a limited-area model whose results are impacted by zonal and meridional boundaries. Here, we present the setup of a latitude-belt domain that reduces disturbances originating from the western and eastern boundaries and therefore allows for studying the impact of model resolution and physical parameterization. The Weather Research and Forecasting (WRF) model coupled to the NOAH land–surface model was operated during July and August 2013 at two different horizontal resolutions, namely 0.03 (HIRES) and 0.12° (LOWRES). Both simulations were forced by the European Centre for Medium-Range Weather Forecasts (ECMWF) operational analysis data at the northern and southern domain boundaries, and the high-resolution Operational Sea Surface Temperature and Sea Ice Analysis (OSTIA) data at the sea surface.The simulations are compared to the operational ECMWF analysis for the representation of large-scale features. To analyze the simulated precipitation, the operational ECMWF forecast, the CPC MORPHing (CMORPH), and the ENSEMBLES gridded observation precipitation data set (E-OBS) were used as references.Analyzing pressure, geopotential height, wind, and temperature fields as well as precipitation revealed (1) a benefit from the higher resolution concerning the reduction of monthly biases, root mean square error, and an improved Pearson skill score, and (2) deficiencies in the physical parameterizations leading to notable biases in distinct regions like the polar Atlantic for the LOWRES simulation, the North Pacific, and Inner Mongolia for both resolutions.In summary, the application of a latitude belt on a convection-permitting resolution shows promising results that are beneficial for future seasonal forecasting.


2017 ◽  
Vol 18 (7) ◽  
pp. 2029-2042
Author(s):  
Tony E. Wong ◽  
William Kleiber ◽  
David C. Noone

Abstract Land surface models are notorious for containing many parameters that control the exchange of heat and moisture between land and atmosphere. Properly modeling the partitioning of total evapotranspiration (ET) between transpiration and evaporation is critical for accurate hydrological modeling, but depends heavily on the treatment of turbulence within and above canopies. Previous work has constrained estimates of evapotranspiration and its partitioning using statistical approaches that calibrate land surface model parameters by assimilating in situ measurements. These studies, however, are silent on the impacts of the accounting of uncertainty within the statistical calibration framework. The present study calibrates the aerodynamic, leaf boundary layer, and stomatal resistance parameters, which partially control canopy turbulent exchange and thus the evapotranspiration flux partitioning. Using an adaptive Metropolis–Hastings algorithm to construct a Markov chain of draws from the joint posterior distribution of these resistance parameters, an ensemble of model realizations is generated, in which latent and sensible heat fluxes and top soil layer temperature are optimized. A set of five calibration experiments demonstrate that model performance is sensitive to the accounting of various sources of uncertainty in the field observations and model output and that it is critical to account for model structural uncertainty. After calibration, the modeled fluxes and top soil layer temperature are largely free from bias, and this calibration approach successfully informs and characterizes uncertainty in these parameters, which is essential for model improvement and development. The key points of this paper are 1) a Markov chain Monte Carlo calibration approach successfully improves modeled turbulent fluxes; 2) ET partitioning estimates hinge on the representation of uncertainties in the model and data; and 3) despite these inherent uncertainties, constrained posterior estimates of ET partitioning emerge.


2018 ◽  
Author(s):  
Marwa Tifafi ◽  
Marta Camino-Serrano ◽  
Christine Hatté ◽  
Hector Morras ◽  
Lucas Moretti ◽  
...  

Abstract. Despite the importance of soil as a large component of the terrestrial ecosystems, the soil compartments are not well represented in the Land Surface Models (LSMs). Indeed, soils in current LSMs are generally represented based on a very simplified schema that can induce a misrepresentation of the deep dynamics of soil carbon. Here, we present a new version of the IPSL-Land Surface Model called ORCHIDEE-SOM, incorporating the 14C dynamic in the soil. ORCHIDEE-SOM, first, simulates soil carbon dynamics for different layers, down to 2 m depth. Second, concentration of dissolved organic carbon (DOC) and its transport are modeled. Finally, soil organic carbon (SOC) decomposition is considered taking into account the priming effect. After implementing the 14C in the soil module of the model, we evaluated model outputs against observations of soil organic carbon and 14C activity (F14C) for different sites with different characteristics. The model managed to reproduce the soil organic carbon stocks and the F14C along the vertical profiles. However, an overestimation of the total carbon stock was noted, but was mostly marked on the surface. Then, thanks to the introduction of 14C, it has been possible to highlight an underestimation of the age of carbon in the soil. Thereafter, two different tests on this new version have been established. The first was to increase carbon residence time of the passive pool and decrease the flux from the slow pool to the passive pool. The second was to establish an equation of diffusion, initially constant throughout the profile, making it vary exponentially as a function of depth. The first modifications did not improve the capacity of the model to reproduce observations whereas the second test showed a decrease of the soil carbon stock overestimation, especially at the surface and an improvement of the estimates of the carbon age. This assumes that we should focus more on vertical variation of soil parameters as a function of depth, mainly for diffusion, in order to upgrade the representation of global carbon cycle in LSMs, thereby helping to improve predictions of the future response of soil organic carbon to global warming.


2011 ◽  
Vol 8 (2) ◽  
pp. 2555-2608 ◽  
Author(s):  
E. H. Sutanudjaja ◽  
L. P. H. van Beek ◽  
S. M. de Jong ◽  
F. C. van Geer ◽  
M. F. P. Bierkens

Abstract. Large-scale groundwater models involving aquifers and basins of multiple countries are still rare due to a lack of hydrogeological data which are usually only available in developed countries. In this study, we propose a novel approach to construct large-scale groundwater models by using global datasets that are readily available. As the test-bed, we use the combined Rhine-Meuse basin that contains groundwater head data used to verify the model output. We start by building a distributed land surface model (30 arc-second resolution) to estimate groundwater recharge and river discharge. Subsequently, a MODFLOW transient groundwater model is built and forced by the recharge and surface water levels calculated by the land surface model. Although the method that we used to couple the land surface and MODFLOW groundwater model is considered as an offline-coupling procedure (i.e. the simulations of both models were performed separately), results are promising. The simulated river discharges compare well to the observations. Moreover, based on our sensitivity analysis, in which we run several groundwater model scenarios with various hydrogeological parameter settings, we observe that the model can reproduce the observed groundwater head time series reasonably well. However, we note that there are still some limitations in the current approach, specifically because the current offline-coupling technique simplifies dynamic feedbacks between surface water levels and groundwater heads, and between soil moisture states and groundwater heads. Also the current sensitivity analysis ignores the uncertainty of the land surface model output. Despite these limitations, we argue that the results of the current model show a promise for large-scale groundwater modeling practices, including for data-poor environments and at the global scale.


2020 ◽  
Author(s):  
Elizabeth Cooper ◽  
Eleanor Blyth ◽  
Hollie Cooper ◽  
Rich Ellis ◽  
Ewan Pinnington ◽  
...  

Abstract. Soil moisture predictions from land surface models are important in hydrological, ecological and meteorological applications. In recent years the availability of wide-area soil-moisture measurements has increased, but few studies have combined model-based soil moisture predictions with in-situ observations beyond the point scale. Here we show that we can markedly improve soil moisture estimates from the JULES land surface model using field scale observations and data assimilation techniques. Rather than directly updating soil moisture estimates towards observed values, we optimize constants in the underlying pedotransfer functions, which relate soil texture to JULES soil physics parameters. In this way we generate a single set of newly calibrated pedotransfer functions based on observations from a number of UK sites with different soil textures. We demonstrate that calibrating a pedotransfer function in this way can improve the performance of land surface models, leading to the potential for better flood, drought and climate projections.


2008 ◽  
Vol 5 (5) ◽  
pp. 4161-4207 ◽  
Author(s):  
H. W. Ter Maat ◽  
R. W. A. Hutjes

Abstract. A large scale mismatch exists between our understanding and quantification of ecosystem atmosphere exchange of carbon dioxide at local scale and continental scales. This paper will focus on the carbon exchange on the regional scale to address the following question: What are the main controlling factors determining atmospheric carbon dioxide content at a regional scale? We use the Regional Atmospheric Modelling System (RAMS), coupled with a land surface scheme simulating carbon, heat and momentum fluxes (SWAPS-C), and including also sub models for urban and marine fluxes, which in principle include the main controlling mechanisms and capture the relevant dynamics of the system. To validate the model, observations are used which were taken during an intensive observational campaign in the central Netherlands in summer 2002. These included flux-site observations, vertical profiles at tall towers and spatial fluxes of various variables taken by aircraft. The coupled regional model (RAMS-SWAPS-C) generally does a good job in simulating results close to reality. The validation of the model demonstrates that surface fluxes of heat, water and CO2 are reasonably well simulated. The comparison against aircraft data shows that the regional meteorology is captured by the model. Comparing spatially explicit simulated and observed fluxes we conclude that in general simulated latent heat fluxes are underestimated by the model to the observations which exhibit large standard deviation for all flights. Sensitivity experiments demonstrated the relevance of the urban emissions of carbon dioxide for the carbon balance in this particular region. The same test also show the relation between uncertainties in surface fluxes and those in atmospheric concentrations.


2021 ◽  
Author(s):  
Adam Pasik ◽  
Wolfgang Preimesberger ◽  
Bernhard Bauer-Marschallinger ◽  
Wouter Dorigo

&lt;p&gt;Multiple satellite-based global surface soil moisture (SSM) datasets are presently available, these however, address exclusively the top layer of the soil (0-5cm). Meanwhile, root-zone soil moisture cannot be directly quantified with remote sensing but can be estimated from SSM using a land surface model. Alternatively, soil water index (SWI; calculated from SSM as a function of time needed for infiltration) can be used as a simple approximation of root-zone conditions. SWI is a proxy for deeper layers of the soil profile which control evapotranspiration, and is hence especially important for studying hydrological processes over vegetation-covered areas and meteorological modelling.&lt;/p&gt;&lt;p&gt;Here we introduce the advances in our work on the first operationally capable SWI-based root-zone soil moisture dataset from C3S Soil Moisture v201912 COMBINED product, spanning the period 2002-2020. The uniqueness of this dataset lies in the fact that T-values (temporal lengths ruling the infiltration) characteristic of SWI were translated into particular soil depths making it much more intuitive, user-friendly and easily applicable. Available are volumetric soil moisture values for the top 1 m of the soil profile at 10 cm intervals, where the optimal T-value (T-best) for each soil layer is selected based on a range of correlation metrics with in situ measurements from the International Soil Moisture Network (ISMN) and the relevant soil and climatic parameters.&lt;br&gt;Additionally we present the results of an extensive global validation against in situ measurements (ISMN) as well as the results of investigations into the relationship between a range of soil and climate characteristics and the optimal T-values for particular soil depths.&lt;/p&gt;


Sign in / Sign up

Export Citation Format

Share Document