The eddy covariance (EC) method is a standard micrometeorological technique for monitoring the exchange rate of the main greenhouse gases across the interface between the atmosphere and ecosystems. ...One of the first EC data processing steps is the temporal alignment of the raw, high frequency measurements collected by the sonic anemometer and gas analyser. While different methods have been proposed and are currently applied, the application of the EC method to trace gases measurements highlighted the difficulty of a correct time lag detection when the fluxes are small in magnitude. Failure to correctly synchronise the time series entails a systematic error on covariance estimates and can introduce large uncertainties and biases in the calculated fluxes. This work aims at overcoming these issues by introducing a new time lag detection procedure based on the assessment of the cross-correlation function (CCF) between variables subject to (i) a pre-whitening based on autoregressive filters and (ii) a resampling technique based on block-bootstrapping. Combining pre-whitening and block-bootstrapping facilitates the assessment of the CCF, enhancing the accuracy of time lag detection between variables with correlation of low order of magnitude (i.e. lower than
-
1
) and allowing for a proper estimate of the associated uncertainty. We expect the proposed procedure to significantly improve the temporal alignment of the EC time-series measured by two physically separate sensors, and to be particularly beneficial in centralised data processing pipelines of research infrastructures (e.g. the Integrated Carbon Observation System, ICOS-RI) where the use of robust and fully data-driven methods, like the one we propose, constitutes an essential prerequisite.
Eddy covariance sites are ideally suited for the study of extreme events on ecosystems as they allow the exchange of trace gases and energy fluxes between ecosystems and the lower atmosphere to be ...directly measured on a continuous basis. However, standardized definitions of hydroclimatic extremes are needed to render studies of extreme events comparable across sites. This requires longer datasets than are available from on-site measurements in order to capture the full range of climatic variability. We present a dataset of drought indices based on precipitation (Standardized Precipitation Index, SPI), atmospheric water balance (Standardized Precipitation Evapotranspiration Index, SPEI), and soil moisture (Standardized Soil Moisture Index, SSMI) for 101 ecosystem sites from the Integrated Carbon Observation System (ICOS) with daily temporal resolution from 1950 to 2021. Additionally, we provide simulated soil moisture and evapotranspiration for each site from the Mesoscale Hydrological Model (mHM). These could be utilised for gap-filling or long-term research, among other applications. We validate our data set with measurements from ICOS and discuss potential research avenues.
•Three different understory leaf area index (LAIu) retrieval methods are compared.•The three methods show different strengths/weaknesses depending on LAI scale.•Variation in LAIu is significantly ...related to vegetation diversity.
Leaf area index (LAI) is a key ecological indicator for describing the structure of canopies and for modelling energy exchange between atmosphere and biosphere. While LAI of the forest overstory can be accurately assessed over large spatial scales via remote sensing, LAI of the forest understory (LAIu) is still largely ignored in ecological studies and ecosystem modelling due to the fact that it is often too complex to be destructively sampled or approximated by other site parameters. Additionally, so far only few attempts have been made to retrieve understory LAI via remote sensing, because dense canopies with high LAI are often hindering retrieval algorithms to produce meaningful estimates for understory LAI. Consequently, the forest understory still constitutes a poorly investigated research realm impeding ecological studies to properly account for its contribution to the energy absorption capacity of forest stands. This study aims to compare three conceptually different indirect retrieval methodologies for LAIu over a diverse panel of forest understory types distributed across Europe. For this we carried out near-to-surface measurements of understory reflectance spectra as well as digital surface photography over the extended network of Integrated Carbon Observation System (ICOS) forest ecosystem sites. LAIu was assessed by exploiting the empirical relationship between vegetation cover and light absorption (Beer-Lambert- Bouguer law) as well as by utilizing proposed relationships with two prominent vegetation indices: normalized difference vegetation index (NDVI) and simple ratio (SR). Retrievals from the three methods were significantly correlated with each other (r = 0.63–0.99, RMSE = 0.53–0.72), but exhibited also significant bias depending on the LAI scale. The NDVI based retrieval approach most likely overestimates LAI at productive sites when LAIu > 2, while the simple ratio algorithm overestimates LAIu at sites with sparse understory vegetation and presence of litter or bare soil. The purely empirical method based on the Beer-Lambert law of light absorption seems to offer a good compromise, since it provides reasonable LAIu values at both low and higher LAI ranges. Surprisingly, LAIu variation among sites seems to be largely decoupled from differences in climate and light permeability of the overstory, but significantly increased with vegetation diversity (expressed as species richness) and hence proposes new applications of LAIu in ecological modelling.
For an assessment of the roles of soil and vegetation in the climate system, a further understanding of the flux components of H2O and CO2 (e.g., transpiration, soil respiration) and their ...interaction with physical conditions and physiological functioning of plants and ecosystems is necessary. To obtain magnitudes of these flux components, we applied source partitioning approaches after Scanlon and Kustas (2010; SK10) and after Thomas et al. (2008; TH08) to high-frequency eddy covariance measurements of 12 study sites covering different ecosystems (croplands, grasslands, and forests) in different climatic regions. Both partitioning methods are based on higher-order statistics of the H2O and CO2 fluctuations, but proceed differently to estimate transpiration, evaporation, net primary production, and soil respiration. We compared and evaluated the partitioning results obtained with SK10 and TH08, including slight modifications of both approaches. Further, we analyzed the interrelations among the performance of the partitioning methods, turbulence characteristics, and site characteristics (such as plant cover type, canopy height, canopy density, and measurement height). We were able to identify characteristics of a data set that are prerequisites for adequate performance of the partitioning methods. SK10 had the tendency to overestimate and TH08 to underestimate soil flux components. For both methods, the partitioning of CO2 fluxes was less robust than for H2O fluxes. Results derived with SK10 showed relatively large dependencies on estimated water use efficiency (WUE) at the leaf level, which is a required input. Measurements of outgoing longwave radiation used for the estimation of foliage temperature (used in WUE) could slightly increase the quality of the partitioning results. A modification of the TH08 approach, by applying a cluster analysis for the conditional sampling of respiration–evaporation events, performed satisfactorily, but did not result in significant advantages compared to the original method versions developed by Thomas et al. (2008). The performance of each partitioning approach was dependent on meteorological conditions, plant development, canopy height, canopy density, and measurement height. Foremost, the performance of SK10 correlated negatively with the ratio between measurement height and canopy height. The performance of TH08 was more dependent on canopy height and leaf area index. In general, all site characteristics that increase dissimilarities between scalars appeared to enhance partitioning performance for SK10 and TH08.
Background: Capturing the response of forest ecosystems to inter-annual climate variability is a great challenge.In this study, we tested the capability of an individual-based forest gap model to ...display carbon fluxes at yearly and daily time scales.The forest model was applied to a spruce forest to simulate the gross primary production(GPP), respiration and net ecosystem exchange(NEE).We analyzed how the variability in climate affected simulated carbon fluxes at the scale of the forest model.Results: Six years were simulated at a daily time scale and compared to the observed eddy covariance(EC) data.In general, the seasonal cycle of the individual carbon fluxes was correctly described by the forest model.However, the estimated GPP differed from the observed data on the days of extreme climatic conditions.Two new parameterizations were developed: one resulting from a numerical calibration, and the other resulting from a filtering method.We suggest new parameter values and even a new function for the temperature limitation of photosynthesis.Conclusions: The forest model reproduced the observed carbon fluxes of a forest ecosystem quite wel.Of the three parameterizations, the calibrated model version performed best.However, the filtering approach showed that calibrated parameter values do not necessarily correctly display the individual functional relations.The concept of simulating forest dynamics at the individual base is a valuable tool for simulating the NEE, GPP and respiration of forest ecosystems.
Extensive field measurements have been performed at three CarboEurope-Integrated Project forest sites with different topography (Renon/Ritten, Italian Alps, Italy; Wetzstein, Thuringia, Germany; ...Norunda, Uppland, Sweden) to evaluate the relevant terms of the carbon balance by measuring CO₂ concentrations CO₂ and the wind field in a 3D multi-tower cube setup. The same experimental setup (geometry and instrumentation) and the same methodology were applied to all the three experiments. It is shown that all sites are affected by advection in different ways and strengths. Everywhere, vertical advection (F VA) occurred only at night. During the day, F VA disappeared because of turbulent mixing, leading to a uniform vertical profile of CO₂. Mean F VA was nearly zero at the hilly site (Wetzstein) and at the flat site (Norunda). However, large, momentary positive or negative contributions occurred at the flat site, whereas vertical non-turbulent fluxes were generally very small at the hilly site. At the slope site (Renon), F VA was always positive at night because of the permanently negative mean vertical wind component resulting from downslope winds. Horizontal advection also occurred mainly at night. It was positive at the slope site and negative at the flat site in the mean diurnal course. The size of the averaged non-turbulent advective fluxes was of the same order of magnitude as the turbulent flux measured by eddy-covariance technique, but the scatter was very high. This implies that it is not advisable to use directly measured quantities of the non-turbulent advective fluxes for the estimation of net ecosystem exchange (NEE) on e.g. an hourly basis. However, situations with and without advection were closely related to local or synoptic meteorological conditions. Thus, it is possible to separate advection affected NEE estimates from fluxes which are representative of the source term. However, the development of a robust correction scheme for advection requires a more detailed site-specific analysis of single events for the identification of the relevant processes. This paper presents mean characteristics of the advective CO₂ fluxes in a first site-to-site comparison and evaluates the main problems for future research.
Differences in the seasonal pattern of assimilatory and respiratory processes are responsible for divergences in seasonal net carbon exchange among ecosystems. Using FLUXNET data (
...http://www.eosdis.ornl.gov/FLUXNET) we have analyzed seasonal patterns of gross primary productivity (
F
GPP), and ecosystem respiration (
F
RE) of boreal and temperate, deciduous and coniferous forests, Mediterranean evergreen systems, a rainforest, temperate grasslands, and C
3 and C
4 crops. Based on generalized seasonal patterns classifications of ecosystems into vegetation functional types can be evaluated for use in global productivity and climate change models. The results of this study contribute to our understanding of respiratory costs of assimilated carbon in various ecosystems.
Seasonal variability of
F
GPP and
F
RE of the investigated sites increased in the order
tropical<
Mediterranean<
temperate
coniferous
<
temperate
deciduous
<
boreal
forests
. Together with the boreal forest sites, the managed grasslands and crops show the largest seasonal variability. In the temperate coniferous forests, seasonal patterns of
F
GPP and
F
RE are in phase, in the temperate deciduous and boreal coniferous forests
F
RE was delayed compared to
F
GPP, resulting in the greatest imbalance between respiratory and assimilatory fluxes early in the growing season.
F
GPP adjusted for the length of the carbon uptake period decreased at the sampling sites across functional types in the order C
4 crops, temperate and boreal deciduous forests
(7.5–8.3
g
C
m
−2
per
day
)>
temperate
conifers, C
3 grassland and crops
(5.7–6.9
g
C
m
−2
per
day
)>
boreal
conifers (4.6
g
C
m
−2
per
day). Annual
F
GPP and net ecosystem productivity (
F
NEP) decreased across climate zones in the order tropical>temperate>boreal. However, the decrease in
F
NEP with latitude was greater than the decrease in
F
GPP, indicating a larger contribution of respiratory (especially heterotrophic) processes in boreal systems.
Both carbon dioxide uptake and albedo of the land surface affect global climate. However, climate change mitigation by increasing carbon uptake can cause a warming trade-off by decreasing albedo, ...with most research focusing on afforestation and its interaction with snow. Here, we present carbon uptake and albedo observations from 176 globally distributed flux stations. We demonstrate a gradual decline in maximum achievable annual albedo as carbon uptake increases, even within subgroups of non-forest and snow-free ecosystems. Based on a paired-site permutation approach, we quantify the likely impact of land use on carbon uptake and albedo. Shifting to the maximum attainable carbon uptake at each site would likely cause moderate net global warming for the first approximately 20 years, followed by a strong cooling effect. A balanced policy co-optimizing carbon uptake and albedo is possible that avoids warming on any timescale, but results in a weaker long-term cooling effect.
The presence or absence of leaves within plant canopies exert a strong influence on the carbon, water and energy balance of ecosystems. Identifying key changes in the timing of leaf elongation and ...senescence during the year can help to understand the sensitivity of different plant functional types to changes in temperature. When recorded over many years these data can provide information on the response of ecosystems to long-term changes in climate. The installation of digital cameras that take images at regular intervals of plant canopies across the Integrated Carbon Observation System ecosystem stations will provide a reliable and important record of variations in canopy state, colour and the timing of key phenological events. Here, we detail the procedure for the implementation of cameras on Integrated Carbon Observation System flux towers and how these images will help us understand the impact of leaf phenology and ecosystem function, distinguish changes in canopy structure from leaf physiology and at larger scales will assist in the validation of (future) remote sensing products. These data will help us improve the representation of phenological responses to climatic variability across Integrated Carbon Observation System stations and the terrestrial biosphere through the improvement of model algorithms and the provision of validation datasets.