The REFLEX project: Comparing different algorithms and implementations for the inversion of a terrestrial ecosystem model against eddy covariance data

We describe a model-data fusion (MDF) inter-comparison project (REFLEX), which compared various algorithms for estimating carbon (C) model parameters consistent with both measured carbon fluxes and states and a simple C model. Participants were provided with the model and with both synthetic net ecosystem exchange (NEE) of CO2 and leaf area index (LAI) data, generated from the model with added noise, and observed NEE and LAI data from two eddy covariance sites. Participants endeavoured to estimate model parameters and states consistent with the model for all cases over the two years for which data were provided, and generate predictions for one additional year without observations. Nine participants contributed results using Metropolis algorithms, Kalman filters and a genetic algorithm. For the synthetic data case, parameter estimates compared well with the true values. The results of the analyses indicated that parameters linked directly to gross primary production (GPP) and ecosystem respiration, such as those related to foliage allocation and turnover, or temperature sensitivity of heterotrophic respiration, were best constrained and characterised. Poorly estimated parameters were those related to the allocation to and turnover of fine root/wood pools. Estimates of confidence intervals varied among algorithms, but several algorithms successfully located the true values of annual fluxes from synthetic experiments within relatively narrow 90% confidence intervals, achieving >80% success rate and mean NEE confidence intervals <110 gC m−2 year−1 for the synthetic case. Annual C flux estimates generated by participants generally agreed with gap-filling approaches using half-hourly data. The estimation of ecosystem respiration and GPP through MDF agreed well with outputs from partitioning studies using half-hourly data. Confidence limits on annual NEE increased by an average of 88% in the prediction year compared to the previous year, when data were available. Confidence intervals on annual NEE increased by 30% when observed data were used instead of synthetic data, reflecting and quantifying the addition of model error. Finally, our analyses indicated that incorporating additional constraints, using data on C pools (wood, soil and fine roots) would help to reduce uncertainties for model parameters poorly served by eddy covariance data.

Impact of land cover uncertainties on estimates of biospheric carbon fluxes

Large-scale bottom-up estimates of terrestrial carbon fluxes, whether based on models or inventory, are highly dependent on the assumed land cover. Most current land cover and land cover change maps are based on satellite data and are likely to be so for the foreseeable future. However, these maps show large differences, both at the class level and when transformed into Plant Functional Types (PFTs), and these can lead to large differences in terrestrial CO2 fluxes estimated by Dynamic Vegetation Models. In this study the Sheffield Dynamic Global Vegetation Model is used. We compare PFT maps and the resulting fluxes arising from the use of widely available moderate (1 km) resolution satellite-derived land cover maps (the Global Land Cover 2000 and several MODIS classification schemes), with fluxes calculated using a reference high (25 m) resolution land cover map specific to Great Britain (the Land Cover Map 2000). We demonstrate that uncertainty is introduced into carbon flux calculations by (1) incorrect or uncertain assignment of land cover classes to PFTs; (2) information loss at coarser resolutions; (3) difficulty in discriminating some vegetation types from satellite data. When averaged over Great Britain, modeled CO2 fluxes derived using the different 1 km resolution maps differ from estimates made using the reference map. The ranges of these differences are 254 gC m−2 a−1 in Gross Primary Production (GPP); 133 gC m−2 a−1 in Net Primary Production (NPP); and 43 gC m−2 a−1 in Net Ecosystem Production (NEP). In GPP this accounts for differences of −15.8% to 8.8%. Results for living biomass exhibit a range of 1109 gC m−2. The types of uncertainties due to land cover confusion are likely to be representative of many parts of the world, especially heterogeneous landscapes such as those found in western Europe.

Soil organic carbon dynamics under perennial energy crops

The European renewable energy directive 2009/28/EC (E.C. 2009) provides a legislative framework for reducing GHG emissions by 20%, while achieving a 20% share of energy from renewable sources by 2020. Perennial energy crops could significantly contribute to limit GHG emissions through replacing equivalent fossil fuels and by sequestering a considerable amount of carbon into the soil through the large amounts of belowground biomass produced. The objective of this study is to evaluate the effects of land use change that perennial energy crops have on croplands (switchgrass) and marginal grasslands (miscanthus). For that purpose above and belowground biomass, SOC variation and Net Ecosystem Exchange were evaluated after five years of growth. At aboveground level both crops produced high biomass under cropland conditions as well as under marginal soils. At belowground level they also produced large amounts of biomass, but no significant influences on SOC in the upper layer (0-30 cm) were found. This is probably because of the "priming effect" that caused fast carbon substitution. In switchgrass only it was found a significant SOC increase in deeper layers (30-60 cm), while in the whole soil profile (0-60 cm) SOC increased from 42 to 51 ha-1. However, the short experimental periods (for both switchgrass and miscanthus), in which land use change was evaluated, do not permit to determine the real capacity of perennial energy crops to accumulate SOC. In conclusion the large amounts of belowground biomass enhanced the SOC dynamic through the priming effect resulting in increased SOC in cropland but not in marginal grassland.

Effect of different management intensities on natural and human-induced carbon fluxes in a golf course in northern Italy

Turfgrasses are ubiquitous in urban landscape and their role on carbon (C) cycle is increasing important also due to the considerable footprint related to their management practices. It is crucial to understand the mechanisms driving the C assimilation potential of these terrestrial ecosystems Several approaches have been proposed to assess C dynamics: micro-meteorological methods, small-chamber enclosure system (SC), chrono-sequence approach and various models. Natural and human-induced variables influence turfgrasses C fluxes. Species composition, environmental conditions, site characteristics, former land use and agronomic management are the most important factors considered in literature driving C sequestration potential. At the same time different approaches seem to influence C budget estimates. In order to study the effect of different management intensities on turfgrass, we estimated net ecosystem exchange (NEE) through a SC approach in a hole of a golf course in the province of Verona (Italy) for one year. The SC approach presented several advantages but also limits related to the measurement frequency, timing and duration overtime, and to the methodological errors connected to the measuring system. Daily CO2 fluxes changed according to the intensity of maintenance, likely due to different inputs and disturbances affecting biogeochemical cycles, combined also to the different leaf area index (LAI). The annual cumulative NEE decreased with the increase of the intensity of management. NEE was related to the seasonality of turfgrass, following temperatures and physiological activity. Generally on the growing season CO2 fluxes towards atmosphere exceeded C sequestered. The cumulative NEE showed a system near to a steady state for C dynamics. In the final part greenhouse gases (GHGs) emissions due to fossil fuel consumption for turfgrass upkeep were estimated, pinpointing that turfgrass may result a considerable C source. The C potential of trees and shrubs needs to be considered to obtain a complete budget.

Carbon pools and fluxes measured during a field campaign conducted in 2010 on the Tibetan Plateau at Kema

The Tibetan highlands host the largest alpine grassland ecosystems worldwide, bearing soils that store substantial stocks of carbon (C) that are very sensitive to land use changes. This study focuses on the cycling of photoassimilated C within a Kobresia pygmaea pasture, the dominating ecosystems on the Tibetan highlands. We investigated short-term effects of grazing cessation and the role of the characteristic Kobresia root turf on C fluxes and belowground C turnover. By combining eddy-covariance measurements with 13CO2 pulse labeling we applied a powerful new approach to measure absolute fluxes of assimilates within and between various pools of the plant-soil-atmosphere system. The roots and soil each store roughly 50% of the overall C in the system (76 Mg C/ha), with only a minor contribution from shoots, which is also expressed in the root:shoot ratio of 90. During June and July the pasture acted as a weak C sink with a strong uptake of approximately 2 g C/m**2/ in the first half of July. The root turf was the main compartment for the turnover of photoassimilates, with a subset of highly dynamic roots (mean residence time 20 days), and plays a key role for the C cycling and C storage in this ecosystem. The short-term grazing cessation only affected aboveground biomass but not ecosystem scale C exchange or assimilate allocation into roots and soil.

(Table 1) Carbon pools in different vegetation components in grazed and ungrazed plots near Ny-Ålesund, Spitzbergen

The carbon (C) sink strength of arctic tundra is under pressure from increasing populations of arctic breeding geese. In this study we examined how CO2 and CH4 fluxes, plant biomass and soil C responded to the removal of vertebrate herbivores in a high arctic wet moss meadow that has been intensively used by barnacle geese (Branta leucopsis) for ca. 20 years. We used 4 and 9 years old grazing exclosures to investigate the potential for recovery of ecosystem function during the growing season (July 2007). The results show greater above- and below-ground vascular plant biomass within the grazing exclosures with graminoid biomass being most responsive to the removal of herbivory whilst moss biomass remained unchanged. The changes in biomass switched the system from net emission to net uptake of CO2 (0.47 and -0.77 µmol/m**2/s in grazed and exclosure plots, respectively) during the growing season and doubled the C storage in live biomass. In contrast, the treatment had no impact on the CH4 fluxes, the total litter C pool or the soil C concentration. The rapid recovery of the above ground biomass and CO2 fluxes demonstrates the plasticity of this high arctic ecosystem in terms of response to changing herbivore pressure.

Geographical variation in carbon dioxide fluxes from soils in agro-ecosystems and its implications for life-cycle assessment

1. Exchange of carbon dioxide (CO2) from soils can contribute significantly to the global warming potential (GWP) of agro-ecosystems. Due to variations in soil type, climatic onditions and land management practices, exchange of CO2 can differ markedly in different geographical locations. The food industry is developing carbon footprints for their products necessitating integration of CO2 exchange from soils with other CO2 emissions along the food chain. It may be advantageous to grow certain crops in different geographical locations to minimize CO2 emissions from the soil, and this may provide potential to offset other emissions in the food chain, such as transport. 2. Values are derived for the C balance of soils growing horticultural crops in the UK, Spain and Uganda. Net ecosystem production (NEP) is firstly calculated from the difference in net primary production (NPP) and heterotrophic soil respiration (Rh). Both NPP and Rh were estimated from intensive direct field measurements. Secondly, net biome production (NBP) is calculated by subtracting the crop biomass from NEP to give an indication of C balance. The importance of soil exchange is discussed in the light of recent discussions on carbon footprints and within the context of food life-cycle assessment (LCA). 3. The amount of crop relative to the biomass and the Rh prevailing in the different countries were the dominant factors influencing the magnitude of NEP and NBP. The majority of the biomass for lettuce Lactuca sativa and vining peas Pisum sativum, was removed from the field as crop; therefore, NEP and NBP were mainly negative. This was amplified for lettuces grown in Uganda (-16·5 and -17 t C ha-1 year-1 compared to UK and Spain -4·8 to 7·4 and -5·1 to 6·3 t C ha-1 year-1 for NEP and NBP, respectively) where the climate elevated Rh. 4. Synthesis and applications. This study demonstrates the importance of soil emissions in the overall life cycle of vegetables. Variability in such emissions suggests that assigning a single value to food carbon footprints may not be adequate, even within a country. Locations with high heterotrophic soil respiration, such as Spain and Uganda (21·9 and 21·6 t C ha-1 year-1, respectively), could mitigate the negative effects of climate on the C costs of crop production by growth of crops with greater returns of residue to the soil. This would minimize net CO2 emissions from these agricultural ecosystems.

Arctic ecosystem responses to changes in water availability and warming: Short and long-term responses

Arctic soils store close to 14% of the global soil carbon. Most of arctic carbon is stored below ground in the permafrost. With climate warming the decomposition of the soil carbon could represent a significant positive feedback to global greenhouse warming. Recent evidence has shown that the temperature of the Arctic is already increasing, and this change is associated mostly with anthropogenic activities. Warmer soils will contribute to permafrost degradation and accelerate organic matter decay and thus increase the flux of carbon dioxide and methane into the atmosphere. Temperature and water availability are also important drivers of ecosystem performance, but effects can be complex and in opposition. Temperature and moisture changes can affect ecosystem respiration (ER) and gross primary productivity (GPP) independently; an increase in the net ecosystem exchange can be a result of either a decrease in ER or an increase in GPP. Therefore, understanding the effects of changes in ecosystem water and temperature on the carbon flux components becomes key to predicting the responses of the Arctic to climate change. The overall goal of this work was to determine the response of arctic systems to simulated climate change scenarios with simultaneous changes in temperature and moisture. A temperature and hydrological manipulation in a naturally-drained lakebed was used to assess the short-term effect of changes in water and temperature on the carbon cycle. Also, as part of International Tundra Experiment Network (ITEX), I determined the long-term effect of warming on the carbon cycle in a natural hydrological gradient established in the mid 90's. I found that the carbon balance is highly sensitive to short-term changes in water table and warming. However, over longer time periods, hydrological and temperature changed soil biophysical properties, nutrient cycles, and other ecosystem structural and functional components that down regulated GPP and ER, especially in wet areas.

Controls on mangrove forest-atmosphere carbon dioxide exchanges in western Everglades National Park

We report on net ecosystem production (NEP) and key environmental controls on net ecosystem exchange (NEE) of carbon dioxide (CO2) between a mangrove forest and the atmosphere in the coastal Florida Everglades. An eddy covariance system deployed above the canopy was used to determine NEE during January 2004 through August 2005. Maximum daytime NEE ranged from −20 to −25 mmol (CO2) m−2 s−1 between March and May. Respiration (Rd) was highly variable (2.81 ± 2.41 mmol (CO2) m−2 s−1), reaching peak values during the summer wet season. During the winter dry season, forest CO2 assimilation increased with the proportion of diffuse solar irradiance in response to greater radiative transfer in the forest canopy. Surface water salinity and tidal activity were also important controls on NEE. Daily light use efficiency was reduced at high (>34 parts per thousand (ppt)) compared to low (ppt) salinity by 46%. Tidal inundation lowered daytime Rd by ∼0.9 mmol (CO2) m−2 s−1 and nighttime Rd by ∼0.5 mmol (CO2) m−2 s−1. The forest was a sink for atmospheric CO2, with an annual NEP of 1170 ± 127 g C m−2 during 2004. This unusually high NEP was attributed to year‐round productivity and low ecosystem respiration which reached a maximum of only 3 g C m−2 d−1. Tidal export of dissolved inorganic carbon derived from belowground respiration likely lowered the estimates of mangrove forest respiration. These results suggest that carbon balance in mangrove coastal systems will change in response to variable salinity and inundation patterns, possibly resulting from secular sea level rise and climate change. Citation: Barr, J. G., V. Engel, J. D. Fuentes,

Effects of Simulated Drought on the Carbon Balance of Everglades Short-Hydroperiod Marsh

Hydrology drives the carbon balance of wetlands by controlling the uptake and release of CO2 and CH4. Longer dry periods in between heavier precipitation events predicted for the Everglades region, may alter the stability of large carbon pools in this wetland's ecosystems. To determine the effects of drought on CO2 fluxes and CH4 emissions, we simulated changes in hydroperiod with three scenarios that differed in the onset rate of drought (gradual, intermediate, and rapid transition into drought) on 18 freshwater wetland monoliths collected from an Everglades short-hydroperiod marsh. Simulated drought, regardless of the onset rate, resulted in higher net CO2 losses net ecosystem exchange (NEE) over the 22-week manipulation. Drought caused extensive vegetation dieback, increased ecosystem respiration (Reco), and reduced carbon uptake gross ecosystem exchange (GEE). Photosynthetic potential measured by reflective indices (photochemical reflectance index, water index, normalized phaeophytinization index, and the normalized difference vegetation index) indicated that water stress limited GEE and inhibited Reco. As a result of drought-induced dieback, NEE did not offset methane production during periods of inundation. The average ratio of net CH4 to NEE over the study period was 0.06, surpassing the 100-year greenhouse warming compensation point for CH4 (0.04). Drought-induced diebacks of sawgrass (C3) led to the establishment of the invasive species torpedograss (C4) when water was resupplied. These changes in the structure and function indicate that freshwater marsh ecosystems can become a net source of CO2 and CH4 to the atmosphere, even following an extended drought. Future changes in precipitation patterns and drought occurrence/duration can change the carbon storage capacity of freshwater marshes from sinks to sources of carbon to the atmosphere. Therefore, climate change will impact the carbon storage capacity of freshwater marshes by influencing water availability and the potential for positive feedbacks on radiative forcing.

Integrated Carbon Budget Models for the Everglades Terrestrial-Coastal-Oceanic Gradient: Current Status and Needs for Inter-Site Comparisons

Recent studies suggest that coastal ecosystems can bury significantly more C than tropical forests, indicating that continued coastal development and exposure to sea level rise and storms will have global biogeochemical consequences. The Florida Coastal Everglades Long Term Ecological Research (FCE LTER) site provides an excellent subtropical system for examining carbon (C) balance because of its exposure to historical changes in freshwater distribution and sea level rise and its history of significant long-term carbon-cycling studies. FCE LTER scientists used net ecosystem C balance and net ecosystem exchange data to estimate C budgets for riverine mangrove, freshwater marsh, and seagrass meadows, providing insights into the magnitude of C accumulation and lateral aquatic C transport. Rates of net C production in the riverine mangrove forest exceeded those reported for many tropical systems, including terrestrial forests, but there are considerable uncertainties around those estimates due to the high potential for gain and loss of C through aquatic fluxes. C production was approximately balanced between gain and loss in Everglades marshes; however, the contribution of periphyton increases uncertainty in these estimates. Moreover, while the approaches used for these initial estimates were informative, a resolved approach for addressing areas of uncertainty is critically needed for coastal wetland ecosystems. Once resolved, these C balance estimates, in conjunction with an understanding of drivers and key ecosystem feedbacks, can inform cross-system studies of ecosystem response to long-term changes in climate, hydrologic management, and other land use along coastlines.

Seawater carbonate chemistry and net calcification during experiments with coral communities, 2009

Acidification of seawater owing to oceanic uptake of atmospheric CO2 originating from human activities such as burning of fossil fuels and land-use changes has raised serious concerns regarding its adverse effects on corals and calcifying communities. Here we demonstrate a net loss of calcium carbonate (CaCO3) material as a result of decreased calcification and increased carbonate dissolution from replicated subtropical coral reef communities (n=3) incubated in continuous-flow mesocosms subject to future seawater conditions. The calcifying community was dominated by the coral Montipora capitata. Daily average community calcification or Net Ecosystem Calcification (NEC=CaCO3 production - dissolution) was positive at 3.3 mmol CaCO3 m-2 h-1 under ambient seawater pCO2 conditions as opposed to negative at -0.04 mmol CaCO3 m-2 h-1 under seawater conditions of double the ambient pCO2. These experimental results provide support for the conclusion that some net calcifying communities could become subject to net dissolution in response to anthropogenic ocean acidification within this century. Nevertheless, individual corals remained healthy, actively calcified (albeit slower than at present rates), and deposited significant amounts of CaCO3 under the prevailing experimental seawater conditions of elevated pCO2.

Physical controls on carbon dioxide transfer velocity and flux in low-gradient river systems and implications for regional carbon budgets

Outgassing of carbon dioxide (CO(2)) from rivers and streams to the atmosphere is a major loss term in the coupled terrestrial-aquatic carbon cycle of major low-gradient river systems (the term ""river system"" encompasses the rivers and streams of all sizes that compose the drainage network in a river basin). However, the magnitude and controls on this important carbon flux are not well quantified. We measured carbon dioxide flux rates (F(CO2)), gas transfer velocity (k), and partial pressures (p(CO2)) in rivers and streams of the Amazon and Mekong river systems in South America and Southeast Asia, respectively. F(CO2) and k values were significantly higher in small rivers and streams (channels <100 m wide) than in large rivers (channels >100 m wide). Small rivers and streams also had substantially higher variability in k values than large rivers. Observed F(CO2) and k values suggest that previous estimates of basinwide CO(2) evasion from tropical rivers and wetlands have been conservative and are likely to be revised upward substantially in the future. Data from the present study combined with data compiled from the literature collectively suggest that the physical control of gas exchange velocities and fluxes in low-gradient river systems makes a transition from the dominance of wind control at the largest spatial scales (in estuaries and river mainstems) toward increasing importance of water current velocity and depth at progressively smaller channel dimensions upstream. These results highlight the importance of incorporating scale-appropriate k values into basinwide models of whole ecosystem carbon balance.