skip to main content

Title: ORCHIDEE-PEAT (revision 4596), a model for northern peatland CO<sub>2</sub>, water, and energy fluxes on daily to annual scales
Peatlands store substantial amounts of carbon and are vulnerable to climate change. We present a modified version of the Organising Carbon and Hydrology In Dynamic Ecosystems (ORCHIDEE) land surface model for simulating the hydrology, surface energy, and CO2 fluxes of peatlands on daily to annual timescales. The model includes a separate soil tile in each 0.5° grid cell, defined from a global peatland map and identified with peat-specific soil hydraulic properties. Runoff from non-peat vegetation within a grid cell containing a fraction of peat is routed to this peat soil tile, which maintains shallow water tables. The water table position separates oxic from anoxic decomposition. The model was evaluated against eddy-covariance (EC) observations from 30 northern peatland sites, with the maximum rate of carboxylation (Vcmax) being optimized at each site. Regarding short-term day-to-day variations, the model performance was good for gross primary production (GPP) (r2 =  0.76; Nash–Sutcliffe modeling efficiency, MEF  =  0.76) and ecosystem respiration (ER, r2 =  0.78, MEF  =  0.75), with lesser accuracy for latent heat fluxes (LE, r2 =  0.42, MEF  =  0.14) and and net ecosystem CO2 exchange (NEE, r2 =  0.38, MEF  =  0.26). Seasonal variations in GPP, ER, NEE, and energy fluxes on monthly scales showed moderate to high r2 values (0.57–0.86). For spatial across-site gradients of annual mean more » GPP, ER, NEE, and LE, r2 values of 0.93, 0.89, 0.27, and 0.71 were achieved, respectively. Water table (WT) variation was not well predicted (r2 &lt; 0.1), likely due to the uncertain water input to the peat from surrounding areas. However, the poor performance of WT simulation did not greatly affect predictions of ER and NEE. We found a significant relationship between optimized Vcmax and latitude (temperature), which better reflects the spatial gradients of annual NEE than using an average Vcmax value. « less
Authors:
; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; more » ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; « less
Award ID(s):
1637459 1026843 1636476
Publication Date:
NSF-PAR ID:
10051341
Journal Name:
Geoscientific Model Development
Volume:
11
Issue:
2
Page Range or eLocation-ID:
497 to 519
ISSN:
1991-9603
Sponsoring Org:
National Science Foundation
More Like this
  1. Sphagnum-dominated peatlands store more carbon than all of Earth’s forests, playing a large role in the balance of carbon dioxide. However, these carbon sinks face an uncertain future as the changing climate is likely to cause water stress, potentially reducing Sphagnum productivity and transitioning peatlands to carbon sources. A mesocosm experiment was performed on thirty-two peat cores collected from two peatland landforms: elevated mounds (hummocks) and lower, flat areas of the peatland (hollows). Both rainfall treatments and water tables were manipulated, and CO2 fluxes were measured. Other studies have observed peat subsiding and tracking the water table downward when experiencing water stress, thought to be a self-preservation technique termed ‘Mire-breathing’. However, we found that hummocks tended to compress inwards, rather than subsiding towards the lowered water table as significantly as hollows. Lower peat height was linearly associated with reduced gross primary production (GPP) in response to lowered water tables, indicating that peat subsidence did not significantly enhance the resistance of GPP to drought. Conversely, Sphagnum peat compression was found to stabilize GPP, indicating that this mechanism of resilience to drought may transmit across the landscape depending on which Sphagnum landform types are dominant. This study draws direct connections between Sphagnummore »traits and peatland hydrology and carbon cycling.« less
  2. We introduce EcH2O-iso, a new development of the physically-based, fully-distributed ecohydrological model EcH2O where the tracking of water isotopic tracers (2H and 18O) and age has been incorporated. EcH2O-iso is evaluated at a montane, low-energy experimental catchment in eastern Scotland using 16 independent isotope time series from various landscape positions and compartments; encompassing soil water, groundwater, stream water, and plant xylem. We find a good model-observation match in most cases, despite having only calibrated the model using hydrometric data and energy fluxes. These results provide further validation of the physical basis of the model for successfully capturing catchment hydrological functioning, both in terms of the celerity in energy propagation (e.g. runoff generation under prevailing hydraulic gradients) and flow velocities of water molecules (e.g., in consistent tracer concentrations at given locations and times). We also show that the spatially-distributed formulation of EcH2O-iso provides a powerful tool for quantitatively linking water stores and fluxes with spatio-temporal patterns of isotopes ratios and water ages. Finally, our study highlights some model development and benchmarking needs, refined using isotope-based calibration, for hypothesis testing and improved simulations of catchment dynamics that is transferable beyond the catchment landscape studied here.
  3. Abstract. Methane (CH4) emissions from natural landscapes constituteroughly half of global CH4 contributions to the atmosphere, yet largeuncertainties remain in the absolute magnitude and the seasonality ofemission quantities and drivers. Eddy covariance (EC) measurements ofCH4 flux are ideal for constraining ecosystem-scale CH4emissions due to quasi-continuous and high-temporal-resolution CH4flux measurements, coincident carbon dioxide, water, and energy fluxmeasurements, lack of ecosystem disturbance, and increased availability ofdatasets over the last decade. Here, we (1) describe the newly publisheddataset, FLUXNET-CH4 Version 1.0, the first open-source global dataset ofCH4 EC measurements (available athttps://fluxnet.org/data/fluxnet-ch4-community-product/, last access: 7 April 2021). FLUXNET-CH4includes half-hourly and daily gap-filled and non-gap-filled aggregatedCH4 fluxes and meteorological data from 79 sites globally: 42freshwater wetlands, 6 brackish and saline wetlands, 7 formerly drainedecosystems, 7 rice paddy sites, 2 lakes, and 15 uplands. Then, we (2) evaluate FLUXNET-CH4 representativeness for freshwater wetland coverageglobally because the majority of sites in FLUXNET-CH4 Version 1.0 arefreshwater wetlands which are a substantial source of total atmosphericCH4 emissions; and (3) we provide the first global estimates of theseasonal variability and seasonality predictors of freshwater wetlandCH4 fluxes. Our representativeness analysis suggests that thefreshwater wetland sites in the dataset cover global wetland bioclimaticattributes (encompassing energy, moisture, and vegetation-relatedparameters) in arctic, boreal, and temperate regions but only sparselycovermore »humid tropical regions. Seasonality metrics of wetland CH4emissions vary considerably across latitudinal bands. In freshwater wetlands(except those between 20∘ S to 20∘ N) the spring onsetof elevated CH4 emissions starts 3 d earlier, and the CH4emission season lasts 4 d longer, for each degree Celsius increase in meanannual air temperature. On average, the spring onset of increasing CH4emissions lags behind soil warming by 1 month, with very few sites experiencingincreased CH4 emissions prior to the onset of soil warming. Incontrast, roughly half of these sites experience the spring onset of risingCH4 emissions prior to the spring increase in gross primaryproductivity (GPP). The timing of peak summer CH4 emissions does notcorrelate with the timing for either peak summer temperature or peak GPP.Our results provide seasonality parameters for CH4 modeling andhighlight seasonality metrics that cannot be predicted by temperature or GPP(i.e., seasonality of CH4 peak). FLUXNET-CH4 is a powerful new resourcefor diagnosing and understanding the role of terrestrial ecosystems andclimate drivers in the global CH4 cycle, and future additions of sitesin tropical ecosystems and site years of data collection will provide addedvalue to this database. All seasonality parameters are available athttps://doi.org/10.5281/zenodo.4672601 (Delwiche et al., 2021).Additionally, raw FLUXNET-CH4 data used to extract seasonality parameterscan be downloaded from https://fluxnet.org/data/fluxnet-ch4-community-product/ (last access: 7 April 2021), and a completelist of the 79 individual site data DOIs is provided in Table 2 of this paper.« less
  4. Marine sediments, speleothems, paleo-lake elevations, and ice core methane and δ18O of O2 (δ18Oatm) records provide ample evidence for repeated abrupt meridional shifts in tropical rainfall belts throughout the last glacial cycle. To improve understanding of the impact of abrupt events on the global terrestrial biosphere, we present composite records of δ18Oatm and inferred changes in fractionation by the global terrestrial biosphere (ΔεLAND) from discrete gas measurements in the WAIS Divide (WD) and Siple Dome (SD) Antarctic ice cores. On the common WD timescale, it is evident that maxima in ΔεLAND are synchronous with or shortly follow small-amplitude WD CH4 peaks that occur within Heinrich stadials 1, 2, 4, and 5 – periods of low atmospheric CH4 concentrations. These local CH4 maxima have been suggested as markers of abrupt climate responses to Heinrich events. Based on our analysis of the modern seasonal cycle of gross primary productivity (GPP)-weighted δ18O of terrestrial precipitation (the source water for atmospheric O2 production), we propose a simple mechanism by which ΔεLAND tracks the centroid latitude of terrestrial oxygen production. As intense rainfall and oxygen production migrate northward, ΔεLAND should decrease due to the underlying meridional gradient in rainfall δ18O. A southward shift should increase ΔεLAND. Monsoon intensitymore »also influences δ18O of precipitation, and although we cannot determine the relative contributions of the two mechanisms, both act in the same direction. Therefore, we suggest that abrupt increases in ΔεLAND unambiguously imply a southward shift of tropical rainfall. The exact magnitude of this shift, however, remains under-constrained by ΔεLAND.« less
  5. null (Ed.)
    Globally important carbon (C) stores in boreal peatlands are vulnerable to altered hydrology through changes in precipitation and runoff patterns, groundwater inputs, and a changing cryosphere. These changes can affect the extent of boreal wetlands and their ability to sequester and transform C and other nutrients. Variation in precipitation patterns has also been increasing, with greater occurrences of both flooding and drought periods. Recent work has pointed to the increasing role of algal production in regulating C cycling during flooded periods in fen peatlands, but exactly how this affects the C sink-strength of these ecosystems is poorly understood. We evaluated temporal trends in algal biomass, ecosystem C uptake and respiration (using static and floating chamber techniques), and spectroscopic indicators of DOM quality (absorbance and fluorescence) in a boreal rich-fen peatland in which water table position had been experimentally manipulated for 13 years. Superimposed on the water table treatments were natural variations in hydrology, including periods of flooding, which allowed us to examine the legacy effects of flooding on C cycling dynamics. We had a particular focus on understanding the role of algae in regulating C cycling, as the relative contribution of algal production was observed to significantly increase with flooding.more »Ecosystem measures of gross primary production (GPP) increased with algal biomass, with higher algal biomass and GPP measured in the lowered water table treatment two years after persistent flooding. Prior to flooding the lowered treatment was the weakest C sink (as CO 2 ), but this treatment became the strongest sink after flooding. The lower degree of humification (lower humification index, HIX) and yet lower bioavailability (higher spectral slope ratio, Sr) of DOM observed in the raised treatment prior to flooding persisted after two years of flooding. An index of free or bound proteins (tyrosine index, TI) increased with algal biomass across all plots during flooding, and was lowest in the raised treatment. As such, antecedent drainage conditions determined the sink-strength of this rich fen—which was also reflected in DOM characteristics. These findings indicate that monitoring flooding history and its effects on algal production could become important to estimates of C balance in northern wetlands.« less