Search

For carbon (C) trading or any other verifiable C reports, it would be reasonable to identify and quantify continuous changes in carbon stocks at regional scales without high investments into additional C-specific, time- and labor-intensive inventories. Our study demonstrates the potential of using routine stand taxation data from large scale forestry inventories for verifiable quantification of tree biomass C stocks, C stock change rates, and associated uncertainties. Empirical models, parameters, and equations of uncertainty propagation have been assembled and applied to data from a forest management unit in Central Germany (550 000 ha), using stand taxation inventories collected between 1993 and 2006. The study showed: 1) The use of stand taxation data resulted in a verifiable and sufficiently precise (cv = 7%) quantification of tree biomass carbon stocks and their changes at the level of growth-regions (1700 to 140 000Â ha). 2) The forest of the test region accumulated carbon in tree biomass at a mean annual rate of 1.8 (â0.9 to 4.5) tC/ha/yr over the studied period. 3) The taxation inventory data can reveal spatial patterns of rates of C stock changes, specifically low rates of 0.4 tC/ha/yr in the northwest and high rates of 3.0 tC/ha/yr in the south of the study region.

Studies of the carbon sink of forest ecosystems often stratify the studied stands by the dominating species and thereby abstract from differences in the mixed-species, multi-cohort structure of many forests. This case study infers whether the aggregation of forestry data introduces a bias in the estimates of carbon stocks and their changes at the scale of individual stands and the scale of a forest district. The empirical TreeGrOSS-C model was applied to 1616 plots of a forest district in Central Germany to simulate carbon dynamics in biomass, woody debris, and soil. In a first approach each stand was explicitly simulated with all cohorts. In three other approaches the forest inventory data were aggregated in several ways, including a stratification of the stands to 110 classes according to the dominating species, age class, and site conditions. A small but significant bias was confirmed. At stand scale the initial ecosystem carbon stocks by the aggregated approach differed from that of the detailed approach by 2.3%, but at the district scale only by 0.05%. The differences in age between interspersed and dominant cohorts as well as differences in litter production were important for the differences in initial carbon stocks. The amounts of wood extracted by thinning operations were important for the differences in the projection of the carbon stocks over 100 years. Because of the smallness of bias, this case study collects evidence that the approaches, that represent stands or stratums by a single cohort, are valid at the scale of a forest district or larger.

Countries need to assess changes in the carbon stocks of forest soils as a part of national greenhouse gas (GHG) inventories under the United Nations Framework Convention on Climate Change (UNFCCC) and the Kyoto Protocol (KP). Since measuring these changes is expensive, it is likely that many countries will use alternative methods to prepare these estimates. We reviewed seven well-known soil carbon models from the point of view of preparing country-scale soil C change estimates. We first introduced the models and explained how they incorporated the most important input variables. Second, we evaluated their applicability at regional scale considering commonly available data sources. Third, we compiled references to data that exist for evaluation of model performance in forest soils. A range of process-based soil carbon models differing in input data requirements exist, allowing some flexibility to forest soil C accounting. Simple models may be the only reasonable option to estimate soil C changes if available resources are limited. More complex models may be used as integral parts of sophisticated inventories assimilating several data sources. Currently, measurement data for model evaluation are common for agricultural soils, but less data have been collected in forest soils. Definitions of model and measured soil pools often differ, ancillary model inputs require scaling of data, and soil C measurements are uncertain. These issues complicate the preparation of model estimates and their evaluation with empirical data, at large scale. Assessment of uncertainties that accounts for the effect of model choice is important part of inventories estimating large-scale soil C changes. Joint development of models and large-scale soil measurement campaigns could reduce the inconsistencies between models and empirical data, and eventually also the uncertainties of model predictions.

Estimates of litter production are a prerequisite for modeling soil carbon stocks and its changes at regional to national scale. However, the required data on biomass removal is often available only for the recent past. In this study we used yield tables as a source of probable past forest management to drive a single tree based stand growth model. Next, simulated growth and timber volume was converted to tree compartment carbon stocks and biomass turnover. The study explicitly accounted for differences in site quality between stands. In addition we performed a Monte Carlo uncertainty and sensitivity analysis. We exemplify the approach by calculating long-term means of past litter production for 10 species by using yield tables that have been applied in Central Germany during the last century. We found that litter production resulting from harvest residues was almost as large as the one from biomass turnover. Differences in site quality caused large differences in litter production. At a given site quality, the uncertainty in soil carbon inputs were 14%, 17%, and 25% for beech, spruce, and pine stands, respectively. The sensitivity analysis showed that the most influential parameters were associated with foliage biomass and turnover. We conclude that rates of mean past litter production and their uncertainties can reliably be modeled on the basis of yield tables if the model accounts for 1) full rotation length including thinning and final harvest, 2) differences in site quality, and 3) environmental dependency of foliage biomass and foliage turnover.

Globally, soils store two to three times as much carbon as currently resides in the atmosphere, and it is critical to understand how soil greenhouse gas (GHG) emissions and uptake will respond to ongoing climate change. In particular, the soil-to-atmosphere CO2 flux, commonly though imprecisely termed soil respiration (RS ), is one of the largest carbon fluxes in the Earth system. An increasing number of high-frequency RS measurements (typically, from an automated system with hourly sampling) have been made over the last two decades; an increasing number of methane measurements are being made with such systems as well. Such high frequency data are an invaluable resource for understanding GHG fluxes, but lack a central database or repository. Here we describe the lightweight, open-source COSORE (COntinuous SOil REspiration) database and software, that focuses on automated, continuous and long-term GHG flux datasets, and is intended to serve as a community resource for earth sciences, climate change syntheses and model evaluation. Contributed datasets are mapped to a single, consistent standard, with metadata on contributors, geographic location, measurement conditions and ancillary data. The design emphasizes the importance of reproducibility, scientific transparency and open access to data. While being oriented towards continuously measured RS , the database design accommodates other soil-atmosphere measurements (e.g. ecosystem respiration, chamber-measured net ecosystem exchange, methane fluxes) as well as experimental treatments (heterotrophic only, etc.). We give brief examples of the types of analyses possible using this new community resource and describe its accompanying R software package.

The global land surface absorbs about a third of anthropogenic emissions each year, due to the difference between two key processes: ecosystem photosynthesis and respiration. Despite the importance of these two processes, it is not possible to measure either at the ecosystem scale during the daytime. Eddy-covariance measurements are widely used as the closest 'quasi-direct' ecosystem-scale observation from which to estimate ecosystem photosynthesis and respiration. Recent research, however, suggests that current estimates may be biased by up to 25%, due to a previously unaccounted for process: the inhibition of leaf respiration in the light. Yet the extent of inhibition remains debated, and implications for estimates of ecosystem-scale respiration and photosynthesis remain unquantified. Here, we quantify an apparent inhibition of daytime ecosystem respiration across the global FLUXNET eddy-covariance network and identify a pervasive influence that varies by season and ecosystem type. We develop partitioning methods that can detect an apparent ecosystem-scale inhibition of daytime respiration and find that diurnal patterns of ecosystem respiration might be markedly different than previously thought. The results call for the re-evaluation of global terrestrial carbon cycle models and also suggest that current global estimates of photosynthesis and respiration may be biased, some on the order of magnitude of anthropogenic fossil fuel emissions.

Describing the coupling of nitrogen (N), phosphorus (P), and carbon (C) cycles of land ecosystems requires understanding microbial element use efficiencies of soil organic matter (SOM) decomposition. These efficiencies are studied by the Soil Enzyme Steady Allocation Model (SESAM) at the decadal scale. The model assumes that the soil microbial communities and their element use efficiencies develop towards an optimum where the growth of the entire community is maximized. Specifically, SESAM approximated this growth optimization by allocating resources to several SOM-degrading enzymes proportional to the revenue of these enzymes, called the Relative approach. However, a rigorous mathematical treatment of this approximation has been lacking so far. Therefore, in this study we derive explicit formulas of enzyme allocation that maximize the total return from enzymatic processing, called the Optimal approach. Further, we derive another heuristic approach that prescribes the change of allocation without the need of deriving a formulation for the optimal allocation, called the Derivative approach. When comparing predictions across these approaches, we found that the Relative approach was a special case of the Optimal approach valid at sufficiently high microbial biomass. However, at low microbial biomass, it overestimated allocation to the enzymes having lower revenues compared to the Optimal approach. The Derivative-based allocation closely tracked the Optimal allocation. These findings increase our confidence in conclusions drawn from SESAM studies. Moreover, the new developments extend the range of conditions at which valid conclusions can be drawn. Further, based on these findings we formulated the constrained enzyme hypothesis. This hypothesis provides a complementary explanation why some substrates in soil are preserved over decades, although they are often decomposed within a few years in incubation experiments. This study shows how optimality considerations lead to simplified models, new insights, and new hypotheses. This is another step in deriving a simple representation of an adaptive microbial community, which is required for coupled stoichiometric C–N–P dynamic models that are aimed to study decadal processes beyond the ecosystem scale. [ABSTRACT FROM AUTHOR]

Global collections of synthesized flux tower data such as FLUXNET have accelerated scientific progress beyond the eddy covariance community. However, remaining data issues in FLUXNET data pose challenges for users, particularly for multi-site synthesis and modelling activities. Here, we present complementary consistency flags (C2Fs) for flux tower data, which rely on multiple indications of inconsistency among variables, along with a methodology to detect discontinuities in time series. The C2F relates to carbon and energy fluxes, as well as to core meteorological variables, and consists of the following: (1) flags for daily data values, (2) flags for entire-site variables, and (3) flags at time stamps that mark large discontinuities in the time series. The flagging is primarily based on combining outlier scores from a set of predefined relationships among variables. The methodology to detect break points in the time series is based on a non-parametric test for the difference in distributions of model residuals. Applying C2F to the FLUXNET 2015 dataset reveals the following: (1) among the considered variables, gross primary productivity and ecosystem respiration data were flagged most frequently, in particular during rain pulses under dry and hot conditions. This information is useful for modelling and analysing ecohydrological responses. (2) There are elevated flagging frequencies for radiation variables (shortwave, photosynthetically active, and net). This information can improve the interpretation and modelling of ecosystem fluxes with respect to issues in the driver. (3) The majority of long-term sites show temporal discontinuities in the time series of latent energy, net ecosystem exchange, and radiation variables. This should be useful for carefully assessing the results in terms of interannual variations in and trends of ecosystem fluxes. The C2F methodology is flexible for customizing and allows for varying the desired strictness of consistency. We discuss the limitations of the approach that can present starting points for future developments. [ABSTRACT FROM AUTHOR]

Climate change is intensifying the hydrologic cycle and altering ecosystem function, including water flux to the atmosphere through evapotranspiration (ET). ET is made up of evaporation (E) via non‐stomatal surfaces, and transpiration (T) through plant stomata which are impacted by global changes in different ways. E and T are difficult to measure independently at the ecosystem scale, especially across multiple sites that represent different land use and land management strategies. To address this gap in understanding, we applied flux variance similarity (FVS) to quantify how E and T differ across 13 different ecosystems measured using eddy covariance in a 10 × 10 km area from the CHEESEHEAD19 experiment in northern Wisconsin, USA. The study sites included eight forests with a large deciduous broadleaf component, three evergreen needleleaf forests, and two wetlands. Average T/ET for the study period averaged nearly 52% in forested sites and 45% in wetlands, with larger values after excluding periods following rain events when evaporation from canopy interception may be expected. A dominance analysis revealed that environmental variables explained on average 69% of the variance of half‐hourly T, which decreased from summer to autumn. Deciduous and evergreen forests showed similar E trajectories over time despite differences in vegetation phenology, and vapor pressure deficit explained some 13% of the variance E in wetlands but only 5% or less in forests. Retrieval of E and T within a dense network of flux towers lends confidence that FVS is a promising approach for comparing ecosystem hydrology across multiple sites to improve our process‐based understanding of ecosystem water fluxes. Key Points: Eddy covariance‐measured evapotranspiration was partitioned to transpiration and evaporation using flux variance similarity at 14 ecosystemsNet radiation explained most of the variability in transpiration, which was lowest in wetland ecosystems and vapor pressure deficit explained more of the variability of evaporation in wetlands than forestsDaily evaporation varied little over time and amongst sites from summer to autumn [ABSTRACT FROM AUTHOR]

Global collections of synthesized flux tower data such as FLUXNET have accelerated scientific progress beyond the eddy covariance community. However, remaining data issues in FLUXNET data pose challenges for users, particularly for multi-site synthesis and modeling activities. Here we present complementary consistency flags (C2F) for flux tower data, which rely on multiple indications of inconsistency among variables, along with a methodology to detect discontinuities in time series. The C2F relates to carbon and energy fluxes as well as to core meteorological variables, and consists of: (1) flags for daily data values, (2) flags for entire site variables, (3) flags at time stamps that mark large discontinuities in the time series. The flagging is primarily based on combining outlier scores from a set of predefined relationships among variables. The methodology to detect break points in the time series is based on a non-parametric test for the difference of distributions of model residuals. Applying C2F to the FLUXNET 2015 dataset reveals that: (1) Among the considered variables, gross primary productivity and ecosystem respiration data were flagged most frequently, in particular during rain pulses under dry and hot conditions. This information is useful for modelling and analysing ecohydrological responses. (2) there are elevated flagging frequencies for radiation variables (shortwave, photosynthetically active, and net). This information can improve the interpretation and modelling of ecosystem fluxes with respect to issues in the driver. (3) The majority of long-term sites show temporal discontinuities in the time series of latent energy, net ecosystem exchange, and radiation variables. This should be useful for carefully assessing the results on interannual variations and trends of ecosystem fluxes. The C2F methodology is flexible for customizing, and allows for varying the desired strictness of consistency. We discuss the limitations of the approach that can present starting points for future improvements. [ABSTRACT FROM AUTHOR]

Describing the coupling of nitrogen (N), phosphorus (P), and carbon (C) cycles of land ecosystems requires understanding microbial element use efficiencies of soil organic matter (SOM) decomposition. These efficiencies are studied by the soil enzyme steady allocation model (SESAM) at decadal scale. The model assumes that the soil microbial communities and their element use efficiencies develop towards an optimum where the growth of the entire community is maximized. Specifically, SESAM approximated this growth optimization by allocating resources to several SOM degrading enzymes pro- portional to the revenue of these enzymes, called the Relative approach. However, a rigorous mathematical treatment of this approximation has been lacking so far. Therefore, in this study we derive explicit formulas of enzyme allocation that maximize total return from enzymatic processing, called the Optimal approach. Further, we derive another heuristic approach that prescribes the change of allocation without the need of deriving a formulation for the optimal allocation, called the Derivative approach. When comparing pre- dictions across these approaches, we found that the Relative approach was a special case of the Optimal approach valid at sufficiently high microbial biomass. However, at low microbial biomass, it overestimated allocation to the enzymes having lower revenues compared to the Optimal approach. The Derivative-based allocation closely tracked the Optimal allocation. The model finding that the Relative approach was a special case of the more rigorous Optimal approach together with observing the same patterns across optimization approaches increases our confidence into conclusions drawn from SESAM studies. Moreover, the new developments extend the range of conditions at which valid conclusions can be drawn. The new model finding that a smaller set of enzyme types was expressed at low microbial biomass led us to formulate the constrained enzyme hypothesis, which provides a complementary explanation why some substrates in soil are preserved over decades although often being decomposed within a few years in incubation experiments. This study shows how optimality considerations lead to simplified models, new insights and new hypotheses. It is another step in deriving a simple representation of an adaptive micro- bial community, which is required for coupled stoichiometric CNP dynamic models that are aimed to study decadal processes beyond ecosystem scale. [ABSTRACT FROM AUTHOR]

Understanding the coupling of nitrogen (N) and carbon (C) cycles of land ecosystems, requires understanding microbial element use efficiencies of soil organic matter (SOM) decomposition. Whereas important controls of those efficiencies by microbial community adaptations have been shown at the scale of a soil pore, a simplified representation of those controls is needed at the ecosystem scale. However, without abstracting from the many details, models are not identifiable, i.e. can not be fitted without ambiguities to observations. There is a need to find, implement, and validate abstract simplified formulations of theses processes. Therefore, we developed the SESAM model as an abstraction of the more detailed soil enzyme allocation model (SEAM) model and tested, whether it can provide the same decadal-term predictions. SEAM explicitly models community adaptation strategies of resource allocation to extracellular enzymes and enzyme limitations on SOM decomposition. It thus provides a scaling from representing several microbial functional groups to a single holistic microbial community. Here we further abstracted the model using quasi-steady-state assumption for extracellular enzyme pools to derive the SESAM model. SESAM reproduced the priming effect, the SOM banking mechanism, and the damping of fluctuations of carbon use efficiency with microbial competition as predicted by SEAM and other more detailed models. This development is an important step towards more parsimonious representation of soil microbial effects in global land surface models. [ABSTRACT FROM AUTHOR]

Phosphorus (P) availability affects the response of terrestrial ecosystems to environmental and climate change (e.g. elevated CO2), yet the magnitude of this effect remains uncertain. This uncertainty arises mainly from a lack of quantitative understanding of the soil biological and geochemical P cycling processes, particularly the P exchange with soil mineral surfaces, which is often described by a Langmuir sorption isotherm. We first conducted a literature review on P sorption experiments and terrestrial biosphere models (TBMs) using Langmuir isotherm. We then developed a new algorithm to describe the inorganic P exchange between soil solution and soil matrix based on the double-surface Langmuir isotherm and extracted empirical equations to calculate the sorption capacity and Langmuir coefficient. We finally tested the conventional and new models of P sorption at five beech forest sites in Germany along a soil P stock gradient using the QUINCY (Quantifying Interactions between terrestrial Nutrient CYcles and the climate system) TBM. We found that the conventional (single-surface) Langmuir isotherm approach in most TBMs largely differed from P sorption experiments regarding the sorption capacities and Langmuir coefficients, and simulated a too low soil P buffering capacity. Conversely, the double-surface Langmuir isotherm approach adequately reproduced the observed patterns of soil inorganic P pools. The better representation of inorganic P cycling using the double Langmuir approach also improved simulated foliar N and P concentrations, and the patterns of gross primary production and vegetation carbon across the soil P gradient. The novel model generally reduces the estimates of P limitation compared to the conventional model, particularly at the low-P site, as the model constraint of slow inorganic P exchange on plant productivity is reduced. [ABSTRACT FROM AUTHOR]

Abstract: This study provides a comprehensive set of functions for predicting biomass for Common beech (Fagus sylvatica L.) in Central Europe for all major tree compartments. The equations are based [...]