•Variability of soil N and C/N of Italian forests were studied with BRT in 1404 plots.•BRT models explained approximately 50% of the N variability.•Latitude explained a larger share of this ...variability than single climate variables.•Forest category was the main determinant of soil N, explaining 30% of its variability.
N is known to be the most limiting element for vegetation growth in temperate and boreal forests. The expected increases in global temperature are predicted to accelerate N mineralization, therefore incrementing N availability in the soil and affecting the soil C cycle as well. While there is an abundance of C data collected to fulfill the requirements for national GHG accounting, more limited information is available for soil N accumulation and storage in relation to forest categories and altitudinal gradients. The data collected by the second Italian National Forest Inventory, spanning a wide range of temperature and precipitation values (10° latitudinal range), represented a unique opportunity to calculate N content and C/N ratio of the different soil layers to a depth of 30 cm. Boosted Regression Tree (BRT) models were applied to investigate the main determinants of soil N distribution and C/N ratio. Forest category was shown to be the main explanatory factor of soil N variability in seven out of eight models, both for forest floor and mineral soil layers. Moreover latitude explained a larger share of variability than single climate variables. BRT models explained, on average, the 49% of the data variability, with the remaining fraction likely due to soil-related variables that were unaccounted for. Accurate estimations of N pools and their determinants in a climate change perspective are consequently required to predict the potential impact of their degradation on forest soil N pools.
Accurate multi-point monitoring systems are required to derive atmospheric measurements of greenhouse gas concentrations both for the calculation of surface fluxes with inversion transport models and ...for the estimation of non-turbulent components of the mass balance equation (i.e. advection and storage fluxes) at eddy covariance sites. When a single analyser is used to monitor multiple sampling points, the deployment of buffer volumes (BVs) along sampling lines can reduce the uncertainty due to the discrete temporal sampling of the signal. In order to optimize the use of buffer volumes we explored various set-ups by simulating their effect on time series of high-frequency CO2 concentration collected at three Fluxnet sites. Besides, we proposed a novel scheme to calculate half-hourly weighted arithmetic means from discrete point samples, accounting for the probabilistic fraction of the signal generated in the averaging period. Results show that the use of BVs with the new averaging scheme reduces the mean absolute error (MAE) up to 80 % compared to a set-up without BVs and up to 60 % compared to the case with BVs and a standard, non-weighted averaging scheme. The MAE of CO2 concentration measurements was observed to depend on the variability of the concentration field and on the size of BVs, which therefore have to be carefully dimensioned. The optimal volume size depends on two main features of the instrumental set-up: the number of measurement points and the time needed to sample at one point (i.e. line purging plus sampling time). A linear and consistent relationship was observed at all sites between the sampling frequency, which summarizes the two features mentioned above, and the renewal frequency associated with the volume. Ultimately, this empirical relationship can be applied to estimate the optimal volume size according to the technical specifications of the sampling system.
Eddy-covariance measurements, which were performed in the summer 2003 in an alpine grassland site, were used to determine the footprint of CO
2 turbulent fluxes and to investigate its dependence on ...stability. The analysis is based on the spatial variability of carbon sinks generated by the difference in the time of the cutting of two concentric portions of the footprint, located outside and inside the fence surrounding the eddy tower, at approximately 30
m distance. Due to this time difference, turbulent flux measurements were performed both in homogeneous (BEFORE and AFTER cutting) and heterogeneous (BETWEEN cuts) condition of spatial sink distribution. The analysis is based exclusively on daytime measurements. Half-hourly records ranked in two stability classes were used to calculate light response functions separately for three 10-day periods (BEFORE, BETWEEN and AFTER cutting). The contribution of the area inside the fence to the total flux was determined for different photosynthetic photon flux density (PPFD) values considering the BETWEEN cuts light response curve as a weighted average of the BEFORE and AFTER ones. The weight of the BEFORE cutting light response curve has been analytically determined and corresponds to the flux fraction which originates from inside the fence. The experimental estimates of the relative importance of the area inside the fence produced values ranging from 30%, during stable conditions, and up to 80% during unstable conditions. These values derived from observations were later compared with the predictions of three analytical footprint models. One of the models systematically underestimated experimental observations, while observations concurred with the other two, particularly in moderately unstable conditions.
Metabolic theory and body size constraints on biomass production and decomposition suggest that differences in the intrinsic potential net ecosystem production (NEPPOT) should be small among ...contrasting C₃ grasslands and therefore unable to explain the wide range in the annual apparent net ecosystem production (NEPAPP) reported by previous studies. We estimated NEPPOT for nine C₃ grasslands under contrasting climate and management regimes using multiyear eddy covariance data. NEPPOT converged within a narrow range, suggesting little difference in the net carbon dioxide uptake capacity among C₃ grasslands. Our results indicate a unique feature of C₃ grasslands compared with other terrestrial ecosystems and suggest a state of stability in NEPPOT due to tightly coupled production and respiration processes. Consequently, the annual NEPAPP of C₃ grasslands is primarily a function of seasonal and short‐term environmental and management constraints, and therefore especially susceptible to changes in future climate patterns and associated adaptation of management practices.
Abstract Metabolic theory and body size constraints on biomass production and decomposition suggest that differences in the intrinsic potential net ecosystem production ( NEP POT ) should be small ...among contrasting C 3 grasslands and therefore unable to explain the wide range in the annual apparent net ecosystem production ( NEP APP ) reported by previous studies. We estimated NEP POT for nine C 3 grasslands under contrasting climate and management regimes using multiyear eddy covariance data. NEP POT converged within a narrow range, suggesting little difference in the net carbon dioxide uptake capacity among C 3 grasslands. Our results indicate a unique feature of C 3 grasslands compared with other terrestrial ecosystems and suggest a state of stability in NEP POT due to tightly coupled production and respiration processes. Consequently, the annual NEP APP of C 3 grasslands is primarily a function of seasonal and short‐term environmental and management constraints, and therefore especially susceptible to changes in future climate patterns and associated adaptation of management practices.
Winter CO2 fluxes represent an important component of the annual carbon budget in northern ecosystems. Understanding winter respiration processes and their responses to climate change is also central ...to our ability to assess terrestrial carbon cycle and climate feedbacks in the future. However, the factors influencing the spatial and temporal patterns of winter ecosystem respiration (Reco) of northern ecosystems are poorly understood. For this reason, we analyzed eddy covariance flux data from 57 ecosystem sites ranging from ~35° N to ~70° N. Deciduous forests were characterized by the highest winter Reco rates (0.90 0.39 g C m−2 d−1), when winter is defined as the period during which daily air temperature remains below 0 °C. By contrast, arctic wetlands had the lowest winter Reco rates (0.02 0.02 g C m−2 d−1). Mixed forests, evergreen needle-leaved forests, grasslands, croplands and boreal wetlands were characterized by intermediate winter Reco rates (g C m−2 d−1) of 0.70(±0.33), 0.60(±0.38), 0.62(±0.43), 0.49(±0.22) and 0.27(±0.08), respectively. Our cross site analysis showed that winter air (Tair) and soil (Tsoil) temperature played a dominating role in determining the spatial patterns of winter Reco in both forest and managed ecosystems (grasslands and croplands). Besides temperature, the seasonal amplitude of the leaf area index (LAI), inferred from satellite observation, or growing season gross primary productivity, which we use here as a proxy for the amount of recent carbon available for Reco in the subsequent winter, played a marginal role in winter CO2 emissions from forest ecosystems. We found that winter Reco sensitivity to temperature variation across space (QS) was higher than the one over time (interannual, QT). This can be expected because QS not only accounts for climate gradients across sites but also for (positively correlated) the spatial variability of substrate quantity. Thus, if the models estimate future warming impacts on Reco based on QS rather than QT, this could overestimate the impact of temperature changes.
It is well established that warming leads to longer growing seasons in seasonally cold ecosystems. Whether this goes along with an increase in the net ecosystem carbon dioxide (CO2) uptake is much ...more controversial. We studied the effects of warming on the start of the carbon uptake period (CUP) of three mountain grasslands situated along an elevational gradient in the Alps. To this end, we used a simple empirical model of the net ecosystem CO2 exchange, calibrated, and forced with multiyear empirical data from each site. We show that reductions in the quantity and duration of daylight associated with earlier snowmelts were responsible for diminishing returns, in terms of carbon gain, from longer growing seasons caused by reductions in daytime photosynthetic uptake and increases in nighttime losses of CO2. This effect was less pronounced at high, compared to low, elevations, where the start of the CUP occurred closer to the summer solstice when changes in day length and incident radiation are minimal.
Key Points
Effects of warming on the start of the carbon uptake period were simulated
Warming induced earlier snowmelt reduced the quantity and duration of daylight
Earlier snowmelts reduced daytime sink and increased nighttime source of carbon
It is well established that warming leads to longer growing seasons in seasonally cold ecosystems. Whether this goes along with an increase in the net ecosystem carbon dioxide (CO
) uptake is much ...more controversial. We studied the effects of warming on the start of the carbon uptake period (CUP) of three mountain grasslands situated along an elevational gradient in the Alps. To this end we used a simple empirical model of the net ecosystem CO
exchange, calibrated and forced with multi-year empirical data from each site. We show that reductions in the quantity and duration of daylight associated with earlier snowmelts were responsible for diminishing returns, in terms of carbon gain, from longer growing seasons caused by reductions in daytime photosynthetic uptake and increases in nighttime losses of CO
. This effect was less pronounced at high, compared to low, elevations, where the start of the CUP occurred closer to the summer solstice when changes in day length and incident radiation are minimal.
The net ecosystem carbon dioxide (CO₂) exchange (NEE) of nine European mountain grassland ecosystems was measured during 2002-2004 using the eddy covariance method. Overall, the availability of ...photosynthetically active radiation (PPFD) was the single most important abiotic influence factor for NEE. Its role changed markedly during the course of the season, PPFD being a better predictor for NEE during periods favorable for CO₂ uptake, which was spring and autumn for the sites characterized by summer droughts (southern sites) and (peak) summer for the Alpine and northern study sites. This general pattern was interrupted by grassland management practices, that is, mowing and grazing, when the variability in NEE explained by PPFD decreased in concert with the amount of aboveground biomass (BMag). Temperature was the abiotic influence factor that explained most of the variability in ecosystem respiration at the Alpine and northern study sites, but not at the southern sites characterized by a pronounced summer drought, where soil water availability and the amount of aboveground biomass were more or equally important. The amount of assimilating plant area was the single most important biotic variable determining the maximum ecosystem carbon uptake potential, that is, the NEE at saturating PPFD. Good correspondence, in terms of the magnitude of NEE, was observed with many (semi-) natural grasslands around the world, but not with grasslands sown on fertile soils in lowland locations, which exhibited higher maximum carbon gains at lower respiratory costs. It is concluded that, through triggering rapid changes in the amount and area of the aboveground plant matter, the timing and frequency of land management practices is crucial for the short-term sensitivity of the NEE of the investigated mountain grassland ecosystems to climatic drivers.
Basal ecosystem respiration rate (BR), the ecosystem respiration rate at a given temperature, is a common and important parameter in empirical models for quantifying ecosystem respiration (ER) ...globally. Numerous studies have indicated that BR varies in space. However, many empirical ER models still use a global constant BR largely due to the lack of a functional description for BR. In this study, we redefined BR to be ecosystem respiration rate at the mean annual temperature. To test the validity of this concept, we conducted a synthesis analysis using 276 site‐years of eddy covariance data, from 79 research sites located at latitudes ranging from ∼3°S to ∼70°N. Results showed that mean annual ER rate closely matches ER rate at mean annual temperature. Incorporation of site‐specific BR into global ER model substantially improved simulated ER compared to an invariant BR at all sites. These results confirm that ER at the mean annual temperature can be considered as BR in empirical models. A strong correlation was found between the mean annual ER and mean annual gross primary production (GPP). Consequently, GPP, which is typically more accurately modeled, can be used to estimate BR. A light use efficiency GPP model (i.e., EC‐LUE) was applied to estimate global GPP, BR and ER with input data from MERRA (Modern Era Retrospective‐Analysis for Research and Applications) and MODIS (Moderate resolution Imaging Spectroradiometer). The global ER was 103 Pg C yr −1, with the highest respiration rate over tropical forests and the lowest value in dry and high‐latitude areas.
Key Points
Redefine basal ecosystem respiration rate of respiration model
Simulate global basal ecosystem respiration rate
Improve ecosystem respiration model performance