Abstract. Methane (CH4) emissions from natural landscapes constitute roughly half of global CH4 contributions to the atmosphere, yet large uncertainties remain in the absolute magnitude and the seasonality of emission quantities and drivers. Eddy covariance (EC) measurements of CH4 flux are ideal for constraining ecosystem-scale CH4 emissions due to quasi-continuous and high-temporal-resolution CH4 flux measurements, coincident carbon dioxide, water, and energy flux measurements, lack of ecosystem disturbance, and increased availability of datasets over the last decade. Here, we (1) describe the newly published dataset, FLUXNET-CH4 Version 1.0, the first open-source global dataset of CH4 EC measurements (available at https://fluxnet.org/data/fluxnet-ch4-community-product/, last access: 7 April 2021). FLUXNET-CH4 includes half-hourly and daily gap-filled and non-gap-filled aggregated CH4 fluxes and meteorological data from 79 sites globally: 42 freshwater wetlands, 6 brackish and saline wetlands, 7 formerly drained ecosystems, 7 rice paddy sites, 2 lakes, and 15 uplands. Then, we (2) evaluate FLUXNET-CH4 representativeness for freshwater wetland coverage globally because the majority of sites in FLUXNET-CH4 Version 1.0 are freshwater wetlands which are a substantial source of total atmospheric CH4 emissions; and (3) we provide the first global estimates of the seasonal variability and seasonality predictors of freshwater wetland CH4 fluxes. Our representativeness analysis suggests that the freshwater wetland sites in the dataset cover global wetland bioclimatic attributes (encompassing energy, moisture, and vegetation-related parameters) in arctic, boreal, and temperate regions but only sparsely cover humid tropical regions. Seasonality metrics of wetland CH4 emissions vary considerably across latitudinal bands. In freshwater wetlands (except those between 20∘ S to 20∘ N) the spring onset of elevated CH4 emissions starts 3 d earlier, and the CH4 emission season lasts 4 d longer, for each degree Celsius increase in mean annual air temperature. On average, the spring onset of increasing CH4 emissions lags behind soil warming by 1 month, with very few sites experiencing increased CH4 emissions prior to the onset of soil warming. In contrast, roughly half of these sites experience the spring onset of rising CH4 emissions prior to the spring increase in gross primary productivity (GPP). The timing of peak summer CH4 emissions does not correlate with the timing for either peak summer temperature or peak GPP. Our results provide seasonality parameters for CH4 modeling and highlight seasonality metrics that cannot be predicted by temperature or GPP (i.e., seasonality of CH4 peak). FLUXNET-CH4 is a powerful new resource for diagnosing and understanding the role of terrestrial ecosystems and climate drivers in the global CH4 cycle, and future additions of sites in tropical ecosystems and site years of data collection will provide added value to this database. All seasonality parameters are available at https://doi.org/10.5281/zenodo.4672601 (Delwiche et al., 2021). Additionally, raw FLUXNET-CH4 data used to extract seasonality parameters can be downloaded from https://fluxnet.org/data/fluxnet-ch4-community-product/ (last access: 7 April 2021), and a complete list of the 79 individual site data DOIs is provided in Table 2 of this paper.

Abstract. Methane (CH4) emissions from natural landscapes constitute roughly half of global CH4 contributions to the atmosphere, yet large uncertainties remain in the absolute magnitude and the seasonality of emission quantities and drivers. Eddy covariance (EC) measurements of CH4 flux are ideal for constraining ecosystem-scale CH4 emissions due to quasi-continuous and high-temporal-resolution CH4 flux measurements, coincident carbon dioxide, water, and energy flux measurements, lack of ecosystem disturbance, and increased availability of datasets over the last decade. Here, we (1) describe the newly published dataset, FLUXNET-CH4 Version 1.0, the first open-source global dataset of CH4 EC measurements (available at https://fluxnet.org/data/fluxnet-ch4-community-product/, last access: 7 April 2021). FLUXNET-CH4 includes half-hourly and daily gap-filled and non-gap-filled aggregated CH4 fluxes and meteorological data from 79 sites globally: 42 freshwater wetlands, 6 brackish and saline wetlands, 7 formerly drained ecosystems, 7 rice paddy sites, 2 lakes, and 15 uplands. Then, we (2) evaluate FLUXNET-CH4 representativeness for freshwater wetland coverage globally because the majority of sites in FLUXNET-CH4 Version 1.0 are freshwater wetlands which are a substantial source of total atmospheric CH4 emissions; and (3) we provide the first global estimates of the seasonal variability and seasonality predictors of freshwater wetland CH4 fluxes. Our representativeness analysis suggests that the freshwater wetland sites in the dataset cover global wetland bioclimatic attributes (encompassing energy, moisture, and vegetation-related parameters) in arctic, boreal, and temperate regions but only sparsely cover humid tropical regions. Seasonality metrics of wetland CH4 emissions vary considerably across latitudinal bands. In freshwater wetlands (except those between 20∘ S to 20∘ N) the spring onset of elevated CH4 emissions starts 3 d earlier, and the CH4 emission season lasts 4 d longer, for each degree Celsius increase in mean annual air temperature. On average, the spring onset of increasing CH4 emissions lags behind soil warming by 1 month, with very few sites experiencing increased CH4 emissions prior to the onset of soil warming. In contrast, roughly half of these sites experience the spring onset of rising CH4 emissions prior to the spring increase in gross primary productivity (GPP). The timing of peak summer CH4 emissions does not correlate with the timing for either peak summer temperature or peak GPP. Our results provide seasonality parameters for CH4 modeling and highlight seasonality metrics that cannot be predicted by temperature or GPP (i.e., seasonality of CH4 peak). FLUXNET-CH4 is a powerful new resource for diagnosing and understanding the role of terrestrial ecosystems and climate drivers in the global CH4 cycle, and future additions of sites in tropical ecosystems and site years of data collection will provide added value to this database. All seasonality parameters are available at https://doi.org/10.5281/zenodo.4672601 (Delwiche et al., 2021). Additionally, raw FLUXNET-CH4 data used to extract seasonality parameters can be downloaded from https://fluxnet.org/data/fluxnet-ch4-community-product/ (last access: 7 April 2021), and a complete list of the 79 individual site data DOIs is provided in Table 2 of this paper.

• Phenological controls over aerodynamic resistance ( R ah ) were investigated. • R ah exhibits significant seasonal variability across a wide range of sites. • These shifts in R ah were caused by phenology in some ecosystems. • Accounting for variation in kB −1 is important for improving predictions of H . Surface roughness – a key control on land-atmosphere exchanges of heat and momentum – differs between dormant and growing seasons. However, how surface roughness shifts seasonally at fine time scales (e.g., days) in response to changing canopy conditions is not well understood. This study: (1) explores how aerodynamic resistance changes seasonally; (2) investigates what drives these seasonal shifts, including the role of vegetation phenology; and (3) quantifies the importance of including seasonal changes of aerodynamic resistance in “big leaf” models of sensible heat flux ( H ). We evaluated aerodynamic resistance and surface roughness lengths for momentum ( z 0m ) and heat ( z 0h ) using the kB −1 parameter (ln( z 0m / z 0h )). We used AmeriFlux data to obtain surface-roughness estimates, and PhenoCam greenness data for phenology. This analysis included 23 sites and ∼190 site years from deciduous broadleaf, evergreen needleleaf, woody savanna, cropland, grassland, and shrubland plant-functional types (PFTs). Results indicated clear seasonal patterns in aerodynamic resistance to sensible heat transfer ( R ah ). This seasonality tracked PhenoCam-derived start-of-season green-up transitions in PFTs displaying the most significant seasonal changes in canopy structure, with R ah decreasing near green-up transitions. Conversely, in woody savanna sites and evergreen needleleaf forests, patterns in R ah were not linked to green-up. Our findings highlight that decreases in kB −1 are an important control over R ah , explaining > 50% of seasonal variation in R ah across most sites. Decreases in kB −1 during green-up are likely caused by increasing z 0h in response to higher leaf area index. Accounting for seasonal variation in kB −1 is key for predicting H as well; assuming kB −1 to be constant resulted in significant biases that also exhibited strong seasonal patterns. Overall, we found that aerodynamic resistance can be sensitive to phenology in ecosystems having strong seasonality in leaf area, and this linkage is critical for understanding land-atmosphere interactions at seasonal time scales.

• Phenological controls over aerodynamic resistance ( R ah ) were investigated. • R ah exhibits significant seasonal variability across a wide range of sites. • These shifts in R ah were caused by phenology in some ecosystems. • Accounting for variation in kB −1 is important for improving predictions of H . Surface roughness – a key control on land-atmosphere exchanges of heat and momentum – differs between dormant and growing seasons. However, how surface roughness shifts seasonally at fine time scales (e.g., days) in response to changing canopy conditions is not well understood. This study: (1) explores how aerodynamic resistance changes seasonally; (2) investigates what drives these seasonal shifts, including the role of vegetation phenology; and (3) quantifies the importance of including seasonal changes of aerodynamic resistance in “big leaf” models of sensible heat flux ( H ). We evaluated aerodynamic resistance and surface roughness lengths for momentum ( z 0m ) and heat ( z 0h ) using the kB −1 parameter (ln( z 0m / z 0h )). We used AmeriFlux data to obtain surface-roughness estimates, and PhenoCam greenness data for phenology. This analysis included 23 sites and ∼190 site years from deciduous broadleaf, evergreen needleleaf, woody savanna, cropland, grassland, and shrubland plant-functional types (PFTs). Results indicated clear seasonal patterns in aerodynamic resistance to sensible heat transfer ( R ah ). This seasonality tracked PhenoCam-derived start-of-season green-up transitions in PFTs displaying the most significant seasonal changes in canopy structure, with R ah decreasing near green-up transitions. Conversely, in woody savanna sites and evergreen needleleaf forests, patterns in R ah were not linked to green-up. Our findings highlight that decreases in kB −1 are an important control over R ah , explaining > 50% of seasonal variation in R ah across most sites. Decreases in kB −1 during green-up are likely caused by increasing z 0h in response to higher leaf area index. Accounting for seasonal variation in kB −1 is key for predicting H as well; assuming kB −1 to be constant resulted in significant biases that also exhibited strong seasonal patterns. Overall, we found that aerodynamic resistance can be sensitive to phenology in ecosystems having strong seasonality in leaf area, and this linkage is critical for understanding land-atmosphere interactions at seasonal time scales.

• Large-scale eddy-covariance flux datasets need to be used with footprint-awareness • Using a fixed-extent target area across sites can bias model-data integration • Most sites do not represent the dominant land-cover type at a larger spatial extent • A representativeness index provides general guidance for site selection and data use Large datasets of greenhouse gas and energy surface-atmosphere fluxes measured with the eddy-covariance technique (e.g., FLUXNET2015, AmeriFlux BASE) are widely used to benchmark models and remote-sensing products. This study addresses one of the major challenges facing model-data integration: To what spatial extent do flux measurements taken at individual eddy-covariance sites reflect model- or satellite-based grid cells? We evaluate flux footprints—the temporally dynamic source areas that contribute to measured fluxes—and the representativeness of these footprints for target areas (e.g., within 250–3000 m radii around flux towers) that are often used in flux-data synthesis and modeling studies. We examine the land-cover composition and vegetation characteristics, represented here by the Enhanced Vegetation Index (EVI), in the flux footprints and target areas across 214 AmeriFlux sites, and evaluate potential biases as a consequence of the footprint-to-target-area mismatch. Monthly 80% footprint climatologies vary across sites and through time ranging four orders of magnitude from 10 3 to 10 7 m 2 due to the measurement heights, underlying vegetation- and ground-surface characteristics, wind directions, and turbulent state of the atmosphere. Few eddy-covariance sites are located in a truly homogeneous landscape. Thus, the common model-data integration approaches that use a fixed-extent target area across sites introduce biases on the order of 4%–20% for EVI and 6%–20% for the dominant land cover percentage. These biases are site-specific functions of measurement heights, target area extents, and land-surface characteristics. We advocate that flux datasets need to be used with footprint awareness, especially in research and applications that benchmark against models and data products with explicit spatial information. We propose a simple representativeness index based on our evaluations that can be used as a guide to identify site-periods suitable for specific applications and to provide general guidance for data use.

• Large-scale eddy-covariance flux datasets need to be used with footprint-awareness • Using a fixed-extent target area across sites can bias model-data integration • Most sites do not represent the dominant land-cover type at a larger spatial extent • A representativeness index provides general guidance for site selection and data use Large datasets of greenhouse gas and energy surface-atmosphere fluxes measured with the eddy-covariance technique (e.g., FLUXNET2015, AmeriFlux BASE) are widely used to benchmark models and remote-sensing products. This study addresses one of the major challenges facing model-data integration: To what spatial extent do flux measurements taken at individual eddy-covariance sites reflect model- or satellite-based grid cells? We evaluate flux footprints—the temporally dynamic source areas that contribute to measured fluxes—and the representativeness of these footprints for target areas (e.g., within 250–3000 m radii around flux towers) that are often used in flux-data synthesis and modeling studies. We examine the land-cover composition and vegetation characteristics, represented here by the Enhanced Vegetation Index (EVI), in the flux footprints and target areas across 214 AmeriFlux sites, and evaluate potential biases as a consequence of the footprint-to-target-area mismatch. Monthly 80% footprint climatologies vary across sites and through time ranging four orders of magnitude from 10 3 to 10 7 m 2 due to the measurement heights, underlying vegetation- and ground-surface characteristics, wind directions, and turbulent state of the atmosphere. Few eddy-covariance sites are located in a truly homogeneous landscape. Thus, the common model-data integration approaches that use a fixed-extent target area across sites introduce biases on the order of 4%–20% for EVI and 6%–20% for the dominant land cover percentage. These biases are site-specific functions of measurement heights, target area extents, and land-surface characteristics. We advocate that flux datasets need to be used with footprint awareness, especially in research and applications that benchmark against models and data products with explicit spatial information. We propose a simple representativeness index based on our evaluations that can be used as a guide to identify site-periods suitable for specific applications and to provide general guidance for data use.

Evergreen conifer forests are the most prevalent land cover type in North America. Seasonal changes in the color of evergreen forest canopies have been documented with near-surface remote sensing, but the physiological mechanisms underlying these changes, and the implications for photosynthetic uptake, have not been fully elucidated. Here, we integrate on-the-ground phenological observations, leaf-level physiological measurements, near surface hyperspectral remote sensing and digital camera imagery, tower-based CO2 flux measurements, and a predictive model to simulate seasonal canopy color dynamics. We show that seasonal changes in canopy color occur independently of new leaf production, but track changes in chlorophyll fluorescence, the photochemical reflectance index, and leaf pigmentation. We demonstrate that at winter-dormant sites, seasonal changes in canopy color can be used to predict the onset of canopy-level photosynthesis in spring, and its cessation in autumn. Finally, we parameterize a simple temperature-based model to predict the seasonal cycle of canopy greenness, and we show that the model successfully simulates interannual variation in the timing of changes in canopy color. These results provide mechanistic insight into the factors driving seasonal changes in evergreen canopy color and provide opportunities to monitor and model seasonal variation in photosynthetic activity using color-based vegetation indices.

In our recent study in Global Change Biology (Li et al., ), we examined the relationship between solar-induced chlorophyll fluorescence (SIF) measured from the Orbiting Carbon Observatory-2 (OCO-2) and gross primary productivity (GPP) derived from eddy covariance flux towers across the globe, and we discovered that there is a nearly universal relationship between SIF and GPP across a wide variety of biomes. This finding reveals the tremendous potential of SIF for accurately mapping terrestrial photosynthesis globally.

Solar-induced chlorophyll fluorescence (SIF) has been increasingly used as a proxy for terrestrial gross primary productivity (GPP). Previous work mainly evaluated the relationship between satellite-observed SIF and gridded GPP products both based on coarse spatial resolutions. Finer resolution SIF (1.3 km × 2.25 km) measured from the Orbiting Carbon Observatory-2 (OCO-2) provides the first opportunity to examine the SIF–GPP relationship at the ecosystem scale using flux tower GPP data. However, it remains unclear how strong the relationship is for each biome and whether a robust, universal relationship exists across a variety of biomes. Here we conducted the first global analysis of the relationship between OCO-2 SIF and tower GPP for a total of 64 flux sites across the globe encompassing eight major biomes. OCO-2 SIF showed strong correlations with tower GPP at both midday and daily timescales, with the strongest relationship observed for daily SIF at the 757 nm (R2 = 0.72, p < 0.0001). Strong linear relationships between SIF and GPP were consistently found for all biomes (R2 = 0.57–0.79, p < 0.0001) except evergreen broadleaf forests (R2 = 0.16, p < 0.05) at the daily timescale. A higher slope was found for C4 grasslands and croplands than for C3 ecosystems. The generally consistent slope of the relationship among biomes suggests a nearly universal rather than biome-specific SIF–GPP relationship, and this finding is an important distinction and simplification compared to previous results. SIF was mainly driven by absorbed photosynthetically active radiation and was also influenced by environmental stresses (temperature and water stresses) that determine photosynthetic light use efficiency. OCO-2 SIF generally had a better performance for predicting GPP than satellite-derived vegetation indices and a light use efficiency model. The universal SIF–GPP relationship can potentially lead to more accurate GPP estimates regionally or globally. Our findings revealed the remarkable ability of finer resolution SIF observations from OCO-2 and other new or future missions (e.g., TROPOMI, FLEX) for estimating terrestrial photosynthesis across a wide variety of biomes and identified their potential and limitations for ecosystem functioning and carbon cycle studies.