Search for: All records

Environmental observation networks, such as AmeriFlux, are foundational for monitoring ecosystem response to climate change, management practices, and natural disturbances; however, their effectiveness depends on their representativeness for the regions or continents. We proposed an empirical, time series approach to quantify the similarity of ecosystem fluxes across AmeriFlux sites. We extracted the diel and seasonal characteristics (i.e., amplitudes, phases) from carbon dioxide, water vapor, energy, and momentum fluxes, which reflect the effects of climate, plant phenology, and ecophysiology on the observations, and explored the potential aggregations of AmeriFlux sites through hierarchical clustering. While net radiation and temperature showed latitudinal clustering as expected, flux variables revealed a more uneven clustering with many small (number of sites < 5), unique groups and a few large (> 100) to intermediate (15–70) groups, highlighting the significant ecological regulations of ecosystem fluxes. Many identified unique groups were from under-sampled ecoregions and biome types of the International Geosphere-Biosphere Programme (IGBP), with distinct flux dynamics compared to the rest of the network. At the finer spatial scale, local topography, disturbance, management, edaphic, and hydrological regimes further enlarge the difference in flux dynamics within the groups. Nonetheless, our clustering approach is a data-driven method to interpret the AmeriFlux network, informing future cross-site syntheses, upscaling, and model-data benchmarking research. Finally, we highlighted the unique and underrepresented sites in the AmeriFlux network, which were found mainly in Hawaii and Latin America, mountains, and at under- sampled IGBP types (e.g., urban, open water), motivating the incorporation of new/unregistered sites from these groups. 

Abstract Northern peatlands are a globally significant source of methane (CH₄), and emissions are projected to increase due to warming and permafrost loss. Understanding the microbial mechanisms behind patterns in CH₄production in peatlands will be key to predicting annual emissions changes, with stable carbon isotopes (δ¹³C‐CH₄) being a powerful tool for characterizing these drivers. Given that δ¹³C‐CH₄is used in top‐down atmospheric inversion models to partition sources, our ability to model CH₄production pathways and associated δ¹³C‐CH₄values is critical. We sought to characterize the role of environmental conditions, including hydrologic and vegetation patterns associated with permafrost thaw, on δ¹³C‐CH₄values from high‐latitude peatlands. We measured porewater and emitted CH₄stable isotopes, pH, and vegetation composition from five boreal‐Arctic peatlands. Porewater δ¹³C‐CH₄was strongly associated with peatland type, with δ¹³C enriched values obtained from more minerotrophic fens (−61.2 ± 9.1‰) compared to permafrost‐free bogs (−74.1 ± 9.4‰) and raised permafrost bogs (−81.6 ± 11.5‰). Variation in porewater δ¹³C‐CH₄was best explained by sedge cover, CH₄concentration, and the interactive effect of peatland type and pH (r² = 0.50,p < 0.001). Emitted δ¹³C‐CH₄varied greatly but was positively correlated with porewater δ¹³C‐CH₄. We calculated a mixed atmospheric δ¹³C‐CH₄value for northern peatlands of −65.3 ± 7‰ and show that this value is more sensitive to landscape drying than wetting under permafrost thaw scenarios. Our results suggest northern peatland δ¹³C‐CH₄values are likely to shift in the future which has important implications for source partitioning in atmospheric inversion models. 

Abstract Climate change is rapidly altering composition, structure, and functioning of the boreal biome, across North America often broadly categorized into ecoregions. The resulting complex changes in different ecoregions present a challenge for efforts to accurately simulate carbon dioxide (CO₂) and energy exchanges between boreal forests and the atmosphere with terrestrial ecosystem models (TEMs). Eddy covariance measurements provide valuable information for evaluating the performance of TEMs and guiding their development. Here, we compiled a boreal forest model benchmarking dataset for North America by harmonizing eddy covariance and supporting measurements from eight black spruce (Picea mariana)-dominated, mature forest stands. The eight forest stands, located in six boreal ecoregions of North America, differ in stand characteristics, disturbance history, climate, permafrost conditions and soil properties. By compiling various data streams, the benchmarking dataset comprises data to parameterize, force, and evaluate TEMs. Specifically, it includes half-hourly, gap-filled meteorological forcing data, ancillary data essential for model parameterization, and half-hourly, gap-filled or partitioned component flux data on CO₂(net ecosystem production, gross primary production [GPP], and ecosystem respiration [ER]) and energy (latent [LE] and sensible heat [H]) and their daily aggregates screened based on half-hourly gap-filling quality criteria. We present a case study with the Canadian Land Surface Scheme Including Biogeochemical Cycles (CLASSIC) to: (1) demonstrate the utility of our dataset to benchmark TEMs and (2) provide guidance for model development and refinement. Model skill was evaluated using several statistical metrics and further examined through the flux responses to their environmental controls. Our results suggest that CLASSIC tended to overestimate GPP and ER among all stands. Model performance regarding the energy fluxes (i.e., LE and H) varied greatly among the stands and exhibited a moderate correlation with latitude. We identified strong relationships between simulated fluxes and their environmental controls except for H, thus highlighting current strengths and limitations of CLASSIC. 

Abstract At the edge of alpine and Arctic ecosystems all over the world, a transition zone exists beyond which it is either infeasible or unfavorable for trees to exist, colloquially identified as the treeline. We explore the possibility of a thermodynamic basis behind this demarcation in vegetation by considering ecosystems as open systems driven by thermodynamic advantage—defined by vegetation’s ability to dissipate heat from the earth’s surface to the air above the canopy. To deduce whether forests would be more thermodynamically advantageous than existing ecosystems beyond treelines, we construct and examine counterfactual scenarios in which trees exist beyond a treeline instead of the existing alpine meadow or Arctic tundra. Meteorological data from the Italian Alps, United States Rocky Mountains, and Western Canadian Taiga-Tundra are used as forcing for model computation of ecosystem work and temperature gradients at sites on both sides of each treeline with and without trees. Model results indicate that the alpine sites do not support trees beyond the treeline, as their presence would result in excessive CO$$_2$$ ${}_2$loss and extended periods of snowpack due to temperature inversions (i.e., positive temperature gradient from the earth surface to the atmosphere). Further, both Arctic and alpine sites exhibit negative work resulting in positive feedback between vegetation heat dissipation and temperature gradient, thereby extending the duration of temperature inversions. These conditions demonstrate thermodynamic infeasibility associated with the counterfactual scenario of trees existing beyond a treeline. Thus, we conclude that, in addition to resource constraints, a treeline is an outcome of an ecosystem’s ability to self-organize towards the most advantageous vegetation structure facilitated by thermodynamic feasibility. 

Abstract The Arctic–Boreal Zone is rapidly warming, impacting its large soil carbon stocks. Here we use a new compilation of terrestrial ecosystem CO₂fluxes, geospatial datasets and random forest models to show that although the Arctic–Boreal Zone was overall an increasing terrestrial CO₂sink from 2001 to 2020 (mean ± standard deviation in net ecosystem exchange, −548 ± 140 Tg C yr⁻¹; trend, −14 Tg C yr⁻¹;P < 0.001), more than 30% of the region was a net CO₂source. Tundra regions may have already started to function on average as CO₂sources, demonstrating a shift in carbon dynamics. When fire emissions are factored in, the increasing Arctic–Boreal Zone sink is no longer statistically significant (budget, −319 ± 140 Tg C yr⁻¹; trend, −9 Tg C yr⁻¹), and the permafrost region becomes CO₂neutral (budget, −24 ± 123 Tg C yr⁻¹; trend, −3 Tg C yr⁻¹), underscoring the importance of fire in this region. 

Abstract Arctic wetlands are known methane (CH₄) emitters but recent studies suggest that the Arctic CH₄sink strength may be underestimated. Here we explore the capacity of well-drained Arctic soils to consume atmospheric CH₄using >40,000 hourly flux observations and spatially distributed flux measurements from 4 sites and 14 surface types. While consumption of atmospheric CH₄occurred at all sites at rates of 0.092 ± 0.011 mgCH₄ m⁻² h⁻¹(mean ± s.e.), CH₄uptake displayed distinct diel and seasonal patterns reflecting ecosystem respiration. Combining in situ flux data with laboratory investigations and a machine learning approach, we find biotic drivers to be highly important. Soil moisture outweighed temperature as an abiotic control and higher CH₄uptake was linked to increased availability of labile carbon. Our findings imply that soil drying and enhanced nutrient supply will promote CH₄uptake by Arctic soils, providing a negative feedback to global climate change. 

Abstract Process‐based land surface models are important tools for estimating global wetland methane (CH₄) emissions and projecting their behavior across space and time. So far there are no performance assessments of model responses to drivers at multiple time scales. In this study, we apply wavelet analysis to identify the dominant time scales contributing to model uncertainty in the frequency domain. We evaluate seven wetland models at 23 eddy covariance tower sites. Our study first characterizes site‐level patterns of freshwater wetland CH₄fluxes (FCH₄) at different time scales. A Monte Carlo approach was developed to incorporate flux observation error to avoid misidentification of the time scales that dominate model error. Our results suggest that (a) significant model‐observation disagreements are mainly at multi‐day time scales (<15 days); (b) most of the models can capture the CH₄variability at monthly and seasonal time scales (>32 days) for the boreal and Arctic tundra wetland sites but have significant bias in variability at seasonal time scales for temperate and tropical/subtropical sites; (c) model errors exhibit increasing power spectrum as time scale increases, indicating that biases at time scales <5 days could contribute to persistent systematic biases on longer time scales; and (d) differences in error pattern are related to model structure (e.g., proxy of CH₄production). Our evaluation suggests the need to accurately replicate FCH₄variability, especially at short time scales, in future wetland CH₄model developments. 

Abstract Photosynthesis of terrestrial ecosystems in the Arctic-Boreal region is a critical part of the global carbon cycle. Solar-induced chlorophyll Fluorescence (SIF), a promising proxy for photosynthesis with physiological insight, has been used to track gross primary production (GPP) at regional scales. Recent studies have constructed empirical relationships between SIF and eddy covariance-derived GPP as a first step to predicting global GPP. However, high latitudes pose two specific challenges: (a) Unique plant species and land cover types in the Arctic–Boreal region are not included in the generalized SIF-GPP relationship from lower latitudes, and (b) the complex terrain and sub-pixel land cover further complicate the interpretation of the SIF-GPP relationship. In this study, we focused on the Arctic-Boreal vulnerability experiment (ABoVE) domain and evaluated the empirical relationships between SIF for high latitudes from the TROPOspheric Monitoring Instrument (TROPOMI) and a state-of-the-art machine learning GPP product (FluxCom). For the first time, we report the regression slope, linear correlation coefficient, and the goodness of the fit of SIF-GPP relationships for Arctic-Boreal land cover types with extensive spatial coverage. We found several potential issues specific to the Arctic-Boreal region that should be considered: (a) unrealistically high FluxCom GPP due to the presence of snow and water at the subpixel scale; (b) changing biomass distribution and SIF-GPP relationship along elevational gradients, and (c) limited perspective and misrepresentation of heterogeneous land cover across spatial resolutions. Taken together, our results will help improve the estimation of GPP using SIF in terrestrial biosphere models and cope with model-data uncertainties in the Arctic-Boreal region.