The California Mediterranean savanna has harsh summer conditions with minimal soil moisture, high temperature, high incoming solar radiation and little or no precipitation. Deciduous blue oaks, Quercus douglasii Hook. and Arn., are winter-deciduous obligate phreatophytes, transpiring mostly groundwater throughout the summer drought. The objective of this work is to fully characterize the seasonal trends of photosynthesis in blue oaks as well as the mechanistic relationships between leaf structure and function. We estimate radiative load of the leaves via the FLiES model and perform in situ measurements of leaf water potential, leaf nitrogen content, an index of chlorophyll content (SPAD readings), photosynthetic and electron transport capacity, and instantaneous rates of CO2 assimilation and electron transport. We measured multiple trees over 3 years providing data from a range of conditions. Our study included one individual that demonstrated strong drought stress as indicated by changes in SPAD readings, leaf nitrogen and all measures of leaf functioning. In the year following severe environmental stress, one individual altered foliation patterns on the crown but did not die. In all other individuals, we found that net carbon assimilation and photosynthetic capacity decreased during the summer drought. SPAD values, electron transport rate (ETR) and quantum yield of photosystem II (PSII) did not show a strong decrease during the summer drought. In most individuals, PSII activity and SPAD readings did not indicate leaf structural or functional damage throughout the season. While net carbon assimilation was tightly coupled to stomatal conductance, the coupling was not as tight with ETR possibly due to contributions from photorespiration or other protective processes. Our work demonstrates that the blue oaks avoid structural damage by maintaining the capacity to convert and dissipate incoming solar radiation during the hot summer drought and are effective at fixing carbon by maximizing rates during the mild spring conditions.

Water availability is the dominant control of global terrestrial primary productivity with concurrent effects on evapotranspiration and ecosystem respiration, especially in water-limited ecosystems. Process-oriented ecosystem models are critical tools for understanding land-atmosphere exchanges and for up-scaling this information to regional and global scales. Thus, it is important to understand how ecosystem models simulate ecosystem fluxes under changing weather conditions. Here, we applied both time-series analysis and meta-analysis techniques to study how five ecosystem process-oriented models-simulated gross primary production (GPP), ecosystem respiration (Reco), and evapotranspiration (ET). Ecosystem fluxes were simulated for 3 years at a daily time step from four evergreen and three deciduous Mediterranean oak woodlands (21 site-year measurements; 105 site-year-simulations). Mediterranean ecosystems are important test-beds for studying the interannual dynamics of soil moisture on ecosystem mass and energy exchange as they experience cool, wet winters with hot, dry summers and are typically subject to drought. Results show data-model disagreements at multiple temporal scales for GPP, Reco, and ET at both plant functional types. Overall there was a systematic underestimation of the temporal variation of Reco at both plant functional types at temporal scales between weeks and months, and an overestimation at the yearly scale. Modeled Reco was systematically overestimated during drought for all sites, but daily GPP was systematically underestimated only for deciduous sites during drought. In contrast, daily estimates of ET showed good data-model agreement even during drought conditions. This meta-analysis brings attention to the importance of drought conditions for modeling purposes in representing forest dynamics in water-limited ecosystems. © 2013 Springer Science+Business Media New York.

Reliable and fine resolution estimates of surface net-radiation are required for estimating latent and sensible heat fluxes between the land surface and the atmosphere. However, currently, fine resolution estimates of net-radiation are not available and consequently it is challenging to develop multi-year estimates of evapotranspiration at scales that can capture land surface heterogeneity and are relevant for policy and decision-making. We developed and evaluated a global net-radiation product at 5 km and 8-day resolution by combining mutually consistent atmosphere and land data from the Moderate Resolution Imaging Spectroradiometer (MODIS) on board Terra. Comparison with net-radiation measurements from 154 globally distributed sites (414 site-years) from the FLUXNET and Surface Radiation budget network (SURFRAD) showed that the net-radiation product agreed well with measurements across seasons and climate types in the extratropics (Wilmott's index ranged from 0.74 for boreal to 0.63 for Mediterranean sites). Mean absolute deviation between the MODIS and measured net-radiation ranged from 38.0 ± 1.8 W.m-2 in boreal to 72.0 ± 4.1 W.m-2 in the tropical climates. The mean bias was small and constituted only 11%, 0.7%, 8.4%, 4.2%, 13.3%, and 5.4% of the mean absolute error in daytime net-radiation in boreal, Mediterranean, temperate-continental, temperate, semi-arid, and tropical climate, respectively. To assess the accuracy of the broader spatiotemporal patterns, we upscaled error-quantified MODIS net-radiation and compared it with the net-radiation estimates from the coarse spatial (1° x 1°) but high temporal resolution gridded net-radiation product from the Clouds and Earth's Radiant Energy System (CERES). Our estimates agreed closely with the net-radiation estimates from the CERES. Difference between the two was less than 10W.m-22 in 94% of the total land area. MODIS net-radiation product will be a valuable resource for the science community studying turbulent fluxes and energy budget at the Earth's surface.

The leaf economics spectrum^1,2 and the global spectrum of plant forms and functions³ revealed fundamental axes of variation in plant traits, which represent different ecological strategies that are shaped by the evolutionary development of plant species². Ecosystem functions depend on environmental conditions and the traits of species that comprise the ecological communities⁴. However, the axes of variation of ecosystem functions are largely unknown, which limits our understanding of how ecosystems respond as a whole to anthropogenic drivers, climate and environmental variability^4,5. Here we derive a set of ecosystem functions⁶ from a dataset of surface gas exchange measurements across major terrestrial biomes. We find that most of the variability within ecosystem functions (71.8%) is captured by three key axes. The first axis reflects maximum ecosystem productivity and is mostly explained by vegetation structure. The second axis reflects ecosystem water-use strategies and is jointly explained by variation in vegetation height and climate. The third axis, which represents ecosystem carbon-use efficiency, features a gradient related to aridity, and is explained primarily by variation in vegetation structure. We show that two state-of-the-art land surface models reproduce the first and most important axis of ecosystem functions. However, the models tend to simulate more strongly correlated functions than those observed, which limits their ability to accurately predict the full range of responses to environmental changes in carbon, water and energy cycling in terrestrial ecosystems^7,8.

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 103 to 107 m2 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.

Past efforts to synthesize and quantify the magnitude and change in carbon dioxide (CO2) fluxes in terrestrial ecosystems across the rapidly warming Arctic-boreal zone (ABZ) have provided valuable information but were limited in their geographical and temporal coverage. Furthermore, these efforts have been based on data aggregated over varying time periods, often with only minimal site ancillary data, thus limiting their potential to be used in large-scale carbon budget assessments. To bridge these gaps, we developed a standardized monthly database of Arctic-boreal CO2 fluxes (ABCflux) that aggregates in situ measurements of terrestrial net ecosystem CO2 exchange and its derived partitioned component fluxes: gross primary productivity and ecosystem respiration. The data span from 1989 to 2020 with over 70 supporting variables that describe key site conditions (e.g., vegetation and disturbance type), micrometeorological and environmental measurements (e.g., air and soil temperatures), and flux measurement techniques. Here, we describe these variables, the spatial and temporal distribution of observations, the main strengths and limitations of the database, and the potential research opportunities it enables. In total, ABCflux includes 244 sites and 6309 monthly observations; 136 sites and 2217 monthly observations represent tundra, and 108 sites and 4092 observations represent the boreal biome. The database includes fluxes estimated with chamber (19% of the monthly observations), snow diffusion (3%) and eddy covariance (78%) techniques. The largest number of observations were collected during the climatological summer (June-August; 32%), and fewer observations were available for autumn (September-October; 25%), winter (December-February; 18%), and spring (March-May; 25%). ABCflux can be used in a wide array of empirical, remote sensing and modeling studies to improve understanding of the regional and temporal variability in CO2 fluxes and to better estimate the terrestrial ABZ CO2 budget. ABCflux is openly and freely available online (Virkkala et al., 2021b, 10.3334/ORNLDAAC/1934). Copyright:

Tundra and boreal ecosystems encompass the northern circumpolar permafrost region and are experiencing rapid environmental change with important implications for the global carbon (C) budget. We analysed multi-decadal time series containing 302 annual estimates of carbon dioxide (CO2) flux across 70 permafrost and non-permafrost ecosystems, and 672 estimates of summer CO2 flux across 181 ecosystems. We find an increase in the annual CO2 sink across non-permafrost ecosystems but not permafrost ecosystems, despite similar increases in summer uptake. Thus, recent non-growing-season CO2 losses have substantially impacted the CO2 balance of permafrost ecosystems. Furthermore, analysis of interannual variability reveals warmer summers amplify the C cycle (increase productivity and respiration) at putatively nitrogen-limited sites and at sites less reliant on summer precipitation for water use. Our findings suggest that water and nutrient availability will be important predictors of the C-cycle response of these ecosystems to future warming.

The FLUXNET2015 dataset provides ecosystem-scale data on CO₂, water, and energy exchange between the biosphere and the atmosphere, and other meteorological and biological measurements, from 212 sites around the globe (over 1500 site-years, up to and including year 2014). These sites, independently managed and operated, voluntarily contributed their data to create global datasets. Data were quality controlled and processed using uniform methods, to improve consistency and intercomparability across sites. The dataset is already being used in a number of applications, including ecophysiology studies, remote sensing studies, and development of ecosystem and Earth system models. FLUXNET2015 includes derived-data products, such as gap-filled time series, ecosystem respiration and photosynthetic uptake estimates, estimation of uncertainties, and metadata about the measurements, presented for the first time in this paper. In addition, 206 of these sites are for the first time distributed under a Creative Commons (CC-BY 4.0) license. This paper details this enhanced dataset and the processing methods, now made available as open-source codes, making the dataset more accessible, transparent, and reproducible.

The following authors were omitted from the original version of this Data Descriptor: Markus Reichstein and Nicolas Vuichard. Both contributed to the code development and N. Vuichard contributed to the processing of the ERA-Interim data downscaling. Furthermore, the contribution of the co-author Frank Tiedemann was re-evaluated relative to the colleague Corinna Rebmann, both working at the same sites, and based on this re-evaluation a substitution in the co-author list is implemented (with Rebmann replacing Tiedemann). Finally, two affiliations were listed incorrectly and are corrected here (entries 190 and 193). The author list and affiliations have been amended to address these omissions in both the HTML and PDF versions.