Arctic‐boreal landscapes are experiencing profound warming, along with changes in ecosystem moisture status and disturbance from fire. This region is of global importance in terms of carbon feedbacks to climate, yet the sign (sink or source) and magnitude of the Arctic‐boreal carbon budget within recent years remains highly uncertain. Here, we provide new estimates of recent (2003–2015) vegetation gross primary productivity (GPP), ecosystem respiration (
The eddy covariance (EC) method is one of the most widely used approaches to quantify surface‐atmosphere fluxes. However, scaling up from a single EC tower to the landscape remains an open challenge. To address this, we used 63 site years of data to examine simulated annual and growing season sums of carbon fluxes from three paired land‐cover type sites of corn, restored‐prairie, and switchgrass ecosystems. This was also done across the landscape by modeling fluxes using different land‐cover type input data. An artificial neural network (ANN) approach was used to model net ecosystem exchange (NEE), ecosystem respiration (
- Award ID(s):
- 1832042
- NSF-PAR ID:
- 10375566
- Publisher / Repository:
- DOI PREFIX: 10.1029
- Date Published:
- Journal Name:
- Journal of Geophysical Research: Biogeosciences
- Volume:
- 126
- Issue:
- 8
- ISSN:
- 2169-8953
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
more » « less
Abstract R eco), net ecosystem CO2exchange (NEE;R eco − GPP), and terrestrial methane (CH4) emissions for the Arctic‐boreal zone using a satellite data‐driven process‐model for northern ecosystems (TCFM‐Arctic), calibrated and evaluated using measurements from >60 tower eddy covariance (EC) sites. We used TCFM‐Arctic to obtain daily 1‐km2flux estimates and annual carbon budgets for the pan‐Arctic‐boreal region. Across the domain, the model indicated an overall average NEE sink of −850 Tg CO2‐C year−1. Eurasian boreal zones, especially those in Siberia, contributed to a majority of the net sink. In contrast, the tundra biome was relatively carbon neutral (ranging from small sink to source). Regional CH4emissions from tundra and boreal wetlands (not accounting for aquatic CH4) were estimated at 35 Tg CH4‐C year−1. Accounting for additional emissions from open water aquatic bodies and from fire, using available estimates from the literature, reduced the total regional NEE sink by 21% and shifted many far northern tundra landscapes, and some boreal forests, to a net carbon source. This assessment, based on in situ observations and models, improves our understanding of the high‐latitude carbon status and also indicates a continued need for integrated site‐to‐regional assessments to monitor the vulnerability of these ecosystems to climate change. -
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 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 < 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.more » « less
-
more » « less
Abstract Permafrost thaw causes the seasonally thawed active layer to deepen, causing the Arctic to shift toward carbon release as soil organic matter becomes susceptible to decomposition. Ground subsidence initiated by ice loss can cause these soils to collapse abruptly, rapidly shifting soil moisture as microtopography changes and also accelerating carbon and nutrient mobilization. The uncertainty of soil moisture trajectories during thaw makes it difficult to predict the role of abrupt thaw in suppressing or exacerbating carbon losses. In this study, we investigated the role of shifting soil moisture conditions on carbon dioxide fluxes during a 13‐year permafrost warming experiment that exhibited abrupt thaw. Warming deepened the active layer differentially across treatments, leading to variable rates of subsidence and formation of thermokarst depressions. In turn, differential subsidence caused a gradient of moisture conditions, with some plots becoming consistently inundated with water within thermokarst depressions and others exhibiting generally dry, but more variable soil moisture conditions outside of thermokarst depressions. Experimentally induced permafrost thaw initially drove increasing rates of growing season gross primary productivity (GPP), ecosystem respiration (
R eco), and net ecosystem exchange (NEE) (higher carbon uptake), but the formation of thermokarst depressions began to reverse this trend with a high level of spatial heterogeneity. Plots that subsided at the slowest rate stayed relatively dry and supported higher CO2fluxes throughout the 13‐year experiment, while plots that subsided very rapidly into the center of a thermokarst feature became consistently wet and experienced a rapid decline in growing season GPP,R eco, and NEE (lower carbon uptake or carbon release). These findings indicate that Earth system models, which do not simulate subsidence and often predict drier active layer conditions, likely overestimate net growing season carbon uptake in abruptly thawing landscapes. -
This paper describes the formation of, and initial results for, a new FLUXNET coordination network for ecosystem-scale methane (CH 4 ) measurements at 60 sites globally, organized by the Global Carbon Project in partnership with other initiatives and regional flux tower networks. The objectives of the effort are presented along with an overview of the coverage of eddy covariance (EC) CH 4 flux measurements globally, initial results comparing CH 4 fluxes across the sites, and future research directions and needs. Annual estimates of net CH 4 fluxes across sites ranged from −0.2 ± 0.02 g C m –2 yr –1 for an upland forest site to 114.9 ± 13.4 g C m –2 yr –1 for an estuarine freshwater marsh, with fluxes exceeding 40 g C m –2 yr –1 at multiple sites. Average annual soil and air temperatures were found to be the strongest predictor of annual CH 4 flux across wetland sites globally. Water table position was positively correlated with annual CH 4 emissions, although only for wetland sites that were not consistently inundated throughout the year. The ratio of annual CH 4 fluxes to ecosystem respiration increased significantly with mean site temperature. Uncertainties in annual CH 4 estimates due to gap-filling and random errors were on average ±1.6 g C m –2 yr –1 at 95% confidence, with the relative error decreasing exponentially with increasing flux magnitude across sites. Through the analysis and synthesis of a growing EC CH 4 flux database, the controls on ecosystem CH 4 fluxes can be better understood, used to inform and validate Earth system models, and reconcile differences between land surface model- and atmospheric-based estimates of CH 4 emissions.more » « less
-
more » « less
Abstract Surface‐atmosphere fluxes and their drivers vary across space and time. A growing area of interest is in downscaling, localizing, and/or resolving sub‐grid scale energy, water, and carbon fluxes and drivers. Existing downscaling methods require inputs of land surface properties at relatively high spatial (e.g., sub‐kilometer) and temporal (e.g., hourly) resolutions, but many observed land surface drivers are not continuously available at these resolutions. We evaluate an approach to overcome this challenge for land surface temperature (LST), a World Meteorological Organization Essential Climate Variable and a key driver for surface heat fluxes. The Chequamegon Heterogenous Ecosystem Energy‐balance Study Enabled by a High‐density Extensive Array of Detectors (CHEESEHEAD19) field experiment provided a scalable testbed. We downscaled LST from satellites (GOES‐16 and ECOsystem Spaceborne Thermal Radiometer Experiment on Space Station [ECOSTRESS]) with further refinement using airborne hyperspectral imagery. Temporally and spatially downscaled LST compared well to independent observations from a network of 20 micrometeorological towers and piloted aircrafts in addition to Landsat‐based LST retrieval and drone‐based LST observed at one tower site. The downscaled 50‐m hourly LST showed good relationships with tower (
r 2 = 0.79, RMSE = 3.5 K) and airborne (r 2 = 0.75, RMSE = 2.4 K) observations over space and time, with precision lower over wetlands and lakes, and some improvement for capturing spatio‐temporal variation compared to a geostationary satellite. Further downscaling to 10 m using hyperspectral imagery resolved hot and cold spots across the landscape as evidenced by independent drone LST, with significant reduction in RMSE by 1.3 K. These results demonstrate a simple pathway for multi‐sensor retrieval of high space and time resolution LST.
