Influences of observation errors in eddy flux data on inverse model parameter estimation

Eddy covariance data are increasingly used to estimate parameters of ecosystem models. For proper maximum likelihood parameter estimates the error structure in the observed data has to be fully characterized. In this study we propose a method to characterize the random error of the eddy covariance flux data, and analyse error distribution, standard deviation, cross- and autocorrelation of CO<sub>2</sub> and H<sub>2</sub>O flux errors at four different European eddy covariance flux sites. Moreover, we examine how the treatment of those errors and additional systematic errors influence statistical estimates of parameters and their associated uncertainties with three models of increasing complexity – a hyperbolic light response curve, a light response curve coupled to water fluxes and the SVAT scheme BETHY. In agreement with previous studies we find that the error standard deviation scales with the flux magnitude. The previously found strongly leptokurtic error distribution is revealed to be largely due to a superposition of almost Gaussian distributions with standard deviations varying by flux magnitude. The crosscorrelations of CO<sub>2</sub> and H<sub>2</sub>O fluxes were in all cases negligible (<i>R</i><sup>2</sup> below 0.2), while the autocorrelation is usually below 0.6 at a lag of 0.5 h and decays rapidly at larger time lags. This implies that in these cases the weighted least squares criterion yields maximum likelihood estimates. To study the influence of the observation errors on model parameter estimates we used synthetic datasets, based on observations of two different sites. We first fitted the respective models to observations and then added the random error estimates described above and the systematic error, respectively, to the model output. This strategy enables us to compare the estimated parameters with true parameters. We illustrate that the correct implementation of the random error standard deviation scaling with flux magnitude significantly reduces the parameter uncertainty and often yields parameter retrievals that are closer to the true value, than by using ordinary least squares. The systematic error leads to systematically biased parameter estimates, but its impact varies by parameter. The parameter uncertainty slightly increases, but the true parameter is not within the uncertainty range of the estimate. This means that the uncertainty is underestimated with current approaches that neglect selective systematic errors in flux data. Hence, we conclude that potential systematic errors in flux data need to be addressed more thoroughly in data assimilation approaches since otherwise uncertainties will be vastly underestimated

Recent global and regional trends in burned area and their compensating environmental controls

The apparent decline in the global incidence of fire between 1996 and 2015, as measured by satellite- observations of burned area, has been related to socioeconomic and land use changes. However, recent decades have also seen changes in climate and vegetation that influence fire and fire-enabled vegetation models do not reproduce the apparent decline. Given that the satellite-derived burned area datasets are still relatively short (<20 years), this raises questions both about the robustness of the apparent decline and what causes it. We use two global satellite-derived burned area datasets and a data-driven fire model to (1) assess the spatio-temporal robustness of the burned area trends and (2) to relate the trends to underlying changes in temperature, precipitation, human population density and vegetation conditions. Although the satellite datasets and simulation all show a decline in global burned area over ~20 years, the trend is not significant and is strongly affected by the start and end year chosen for trend analysis and the year-to-year variability in burned area. The global and regional trends shown by the two satellite datasets are poorly correlated for the common overlapping period (2001–2015) and the fire model simulates changes in global and regional burned area that lie within the uncertainties of the satellite datasets. The model simulations show that recent increases in temperature would lead to increased burned area but this effect is compensated by increasing wetness or increases in population, both of which lead to declining burned area. Increases in vegetation cover and density associated with recent greening trends lead to increased burned area in fuel-limited regions. Our analyses show that global and regional burned area trends result from the interaction of compensating trends in controls of wildfire at regional scales

Tropical climate–vegetation–fire relationships: multivariate evaluation of the land surface model JSBACH

The interactions between climate, vegetation and fire can strongly influence the future trajectories of vegetation in Earth system models. We evaluate the relationships between tropical climate, vegetation and fire in the global vegetation model JSBACH, using a simple fire scheme and the complex fire model SPITFIRE with the aim to identify potential for model improvement. We use two remote-sensing products (based on MODIS and Landsat) in different resolutions to assess the robustness of the obtained observed relationships. We evaluate the model using a multivariate comparison that allows us to focus on the interactions between climate, vegetation and fire and test the influence of land use change on the modelled patterns. Climate–vegetation–fire relationships are known to differ between continents; we therefore perform the analysis for each continent separately.The observed relationships are similar in the two satellite data sets, but maximum tree cover is reached at higher precipitation values for coarser resolution. This shows that the spatial scale of models and data needs to be consistent for meaningful comparisons. The model captures the broad spatial patterns with regional differences, which are partly due to the climate forcing derived from an Earth system model. Compared to the simple fire scheme, SPITFIRE strongly improves the spatial pattern of burned area and the distribution of burned area along increasing precipitation. The correlation between precipitation and tree cover is higher in the observations than in the largely climate-driven vegetation model, with both fire models. The multivariate comparison identifies excessive tree cover in low-precipitation areas and a too-strong relationship between high fire occurrence and low tree cover for the complex fire model. We therefore suggest that drought effects on tree cover and the impact of burned area on tree cover or the adaptation of trees to fire can be improved.The observed variation in the relationship between precipitation and maximum tree cover between continents is higher than the simulated one. Land use contributes to the intercontinental differences in fire regimes with SPITFIRE and strongly overprints the modelled multimodality of tree cover with SPITFIRE.The multivariate model–data comparison used here has several advantages: it improves the attribution of model–data mismatches to model processes, it reduces the impact of biases in the meteorological forcing on the evaluation and it allows us to evaluate not only a specific target variable but also the interactions.</p

No Consistent Evidence for Advancing or Delaying Trends in Spring Phenology on the Tibetan Plateau

Vegetation phenology is a sensitive indicator of climate change and has significant effects on the exchange of carbon, water, and energy between the terrestrial biosphere and the atmosphere. The Tibetan Plateau, the Earth\u27s “third pole,” is a unique region for studying the long‐term trends in vegetation phenology in response to climate change because of the sensitivity of its alpine ecosystems to climate and its low‐level human disturbance. There has been a debate whether the trends in spring phenology over the Tibetan Plateau have been continuously advancing over the last two to three decades. In this study, we examine the trends in the start of growing season (SOS) for alpine meadow and steppe using the Global Inventory Modeling and Mapping Studies (GIMMS)3g normalized difference vegetation index (NDVI) data set (1982–2014), the GIMMS NDVI data set (1982–2006), the Moderate Resolution Imaging Spectroradiometer (MODIS) NDVI data set (2001–2014), the Satellite Pour l\u27Observation de la Terre Vegetation (SPOT‐VEG) NDVI data set (1999–2013), and the Sea‐viewing Wide Field‐of‐View Sensor (SeaWiFS) NDVI data set (1998–2007). Both logistic and polynomial fitting methods are used to retrieve the SOS dates from the NDVI data sets. Our results show that the trends in spring phenology over the Tibetan Plateau depend on both the NDVI data set used and the method for retrieving the SOS date. There are large discrepancies in the SOS trends among the different NDVI data sets and between the two different retrieval methods. There is no consistent evidence that spring phenology (“green‐up” dates) has been advancing or delaying over the Tibetan Plateau during the last two to three decades. Ground‐based budburst data also indicate no consistent trends in spring phenology. The responses of SOS to environmental factors (air temperature, precipitation, soil temperature, and snow depth) also vary among NDVI data sets and phenology retrieval methods. The increases in winter and spring temperature had offsetting effects on spring phenology

Сутність та класифікація ризиків інвестиційної діяльності

Наводиться визначення поняттю "ризики інвестиційної діяльності" за рахунок поєднання його сутнісних характеристик, виконано узагальнення класифікації цих ризиків, запропоновано введення нової класифікаційної групи – "корпоративні ризики", які пов'язані з можливістю втрати контролю над підприємством інвестором-акціонером

BHPMF – a hierarchical Bayesian approach to gap-filling and trait prediction for macroecology and functional biogeography

Aim: Functional traits of organisms are key to understanding and predicting biodiversity and ecological change, which motivates continuous collection of traits and their integration into global databases. Such trait matrices are inherently sparse, severely limiting their usefulness for further analyses. On the other hand, traits are characterized by the phylogenetic trait signal, trait–trait correlations and environmental constraints, all of which provide information that could be used to statistically fill gaps. We propose the application of probabilistic models which, for the first time, utilize all three characteristics to fill gaps in trait databases and predict trait values at larger spatial scales. Innovation: For this purpose we introduce BHPMF, a ierarchical Bayesian extension of probabilistic matrix actorization (PMF). PMF is a machine learning technique which exploits the correlation structure of sparse matrices to impute missing entries. BHPMF additionally utilizes the taxonomic hierarchy for trait prediction and provides uncertainty estimates for each imputation. In combination with multiple regression against environmental information, BHPMF allows for extrapolation frompoint measurements to larger spatial scales.We demonstrate the applicability of BHPMF in ecological contexts, using different plant functional trait datasets, also comparing results to taking the species mean and PMF. Main conclusions: Sensitivity analyses validate the robustness and accuracy of BHPMF: our method captures the correlation structure of the trait matrix as well as the phylogenetic trait signal – also for extremely sparse trait matrices – and provides a robust measure of confidence in prediction accuracy for each missing entry. The combination of BHPMF with environmental constraints provides a promising concept to extrapolate traits beyond sampled regions, accounting for intraspecific trait variability. We conclude that BHPMF and its derivatives have a high potential to support future trait-based research in macroecology and functional biogeography

What eddy-covariance measurements tell us about prior land flux errors in CO2-flux inversion schemes

0.2 after 200 km). Separating out the plant functional types did not increase the spatial correlations, except for the deciduous broad-leaved forests. Using the statistics of the flux measurements as a proxy for the statistics of the prior flux errors was shown not to be a viable approach. A statistical model allowed us to upscale the site-level flux error statistics to the coarser spatial and temporal resolutions used in regional or global models. This approach allowed us to quantify how aggregation reduces error variances, while increasing correlations. As an example, for a typical inversion of grid point (300 km × 300 km) monthly fluxes, we found that the prior flux error follows an approximate e-folding correlation length of 500 km only, with correlations from one month to the next as large as 0.6

Geologic carbon sources may confound ecosystem carbon balance estimates: Evidence from a semiarid steppe in the southeast of Spain

At a semiarid steppe site located in the SE of Spain, relatively large CO2 emissions were measured that could not be attributed to the ecosystem activity alone. Since the study site was located in a tectonically active area, it was hypothesized that a part of the measured CO2 was of geologic origin. This investigation included a survey of soil CO2 efflux, together with carbon isotope analyses of the CO2 in the soil atmosphere, soil CO2 efflux (ie, Keeling plots), groundwater and local thermal springs. These measurements confirmed the ... Articoli in Schola

The importance of mesophyll conductance in regulating forest ecosystem productivity during drought periods

Water availability is the most limiting factor to global plant productivity, yet photosynthetic responses to seasonal drought cycles are poorly understood, with conflicting reports on which limiting process is the most important during drought. We address the problem using a model-data synthesis approach to look at canopy level fluxes, integrating twenty years of half hour data gathered by the FLUXNET network across six Mediterranean sites. The measured canopy level, water and carbon fluxes were used, together with an inverse canopy ecophysiological model, to estimate the bulk canopy conductance, bulk mesophyll conductance, and the canopy scale carbon pools in both the intercellular spaces and at the site of carboxylation in the chloroplasts. Thus the roles of stomatal and mesophyll conductance in the regulation of internal carbon pools and photosynthesis could be separated. A quantitative limitation analysis allowed for the relative seasonal responses of stomatal, mesophyll, and biochemical limitations to be gauged. The concentration of carbon in the chloroplast was shown to be a potentially more reliable estimator of assimilation rates than the intercellular carbon concentration. Both stomatal conductance limitations and mesophyll conductance limitations were observed to regulate the response of photosynthesis to water stress in each of the six species studied. The results suggest that mesophyll conductance could bridge the gap between conflicting reports on plant responses to soil water stress, and that the inclusion of mesophyll conductance in biosphere-atmosphere transfer models may improve their performance, in particular their ability to accurately capture the response of terrestrial vegetation productivity to drought