Contribution à la caractérisation de sites sableux : signature spectro-directionnelle, distribution en taille et minéralogie extraites d'échantillons de sables

International audienceThe characterization of sands detailed in this paper has been performed in order to support the in-flight radiometric performance assessment of space-borne optical sensors over so-called Pseudo-Invariant Calibration Sites (PICS). Although the physical properties of PICS surface are fairly stable in time, the signal measured from space varies with the illumination and the viewing geometries. Thus there is a need to characterize the spectro-directional properties of PICS. This can be done, at a broad scale, thanks to multi-spectral multi-directional space-borne sensors such as the POLDER instrument (with old data). However, interpolating or extrapolating the spectro-directional reflectances measured from space to spectral bands of another sensor is not straightforward. The hyperspectral characterization of sand samples collected within or nearby PICS can contribute to a solution. In this context, a set of 31 sand samples was compiled. The BiConical Reflectance Factor (BCRF) was measured between 0.4 and 2.5 µm, over a quarter hemisphere when the amount of sand in the sample was large enough and for only a single fixed angular configuration for small samples. These optical measurements were complemented by grain size distribution measurements and mineralogical analysis and compiled together with previously published measurements in the so-called PICSAND database, freely available on line.La caractérisation des sables détaillée dans cet article a été faite en soutien à l'estimation en vol des performances radiométriques des capteurs optiques spatiaux à partir des sites appelés PICS pour Pseudo-Invariant Calibration Sites. Bien que les propriétés physiques des PICS soient relativement stables dans le temps, le signal mesuré depuis l'espace varie en fonction des géométries d'illumination et d'observation. De ce fait, il est nécessaire de caractériser les propriétés spectro-directionnelles des PICS. Ceci peut être fait, à une grande échelle, à partir de capteurs spatiaux multi-spectraux et multi-directionnels tels que le capteur POLDER (avec des données anciennes). Cependant, l'interpolation ou l'extrapolation des réflectances spectro-directionnelles obtenues depuis l'espace aux bandes spectrales d'un autre capteur est délicate. La caractérisation hyperspectrale d'échantillons de sable issus de PICS ou de leur voisinage peut participer à une solution. Dans ce contexte, 31 échantillons de sable ont été collectés. Le Facteur de Reflectance BiConique (BCRF) a été mesuré entre 0,4 et 2,5 µm, pour une demi-hémisphère lorsque la quantité de sable était suffisante, et pour une seule géométrie pour les échantillons plus petits. Ces mesures optiques ont été complétées par des mesures de distribution en taille et par une analyse minéralogique, et mises dans une base de données appelée PICSAND avec d'autres mesures publiées dans la littérature. Cette base de donnée est en libre accès en ligne

Energy and CO2 fluxes in temperate deciduous broadleaf forests: using local add-covariance data for a global optimization of ORCHIDEE model

Non peer reviewedPublisher PD

Land surface model parameter optimisation using in situ flux data : Comparison of gradient-based versus random search algorithms (a case study using ORCHIDEE v1.9.5.2)

This work used eddy covariance data acquired by the FLUXNET community and in particular by the following networks: AmeriFlux (U.S. Department of Energy, Biological and Environmental Research, Terrestrial Carbon Program; DE-FG02-04ER63917 and DE-FG02-04ER63911), AfriFlux, AsiaFlux, CarboAfrica, CarboEuropeIP, CarboItaly, CarboMont, ChinaFlux, Fluxnet-Canada (supported by CFCAS, NSERC, BIOCAP, Environment Canada, and NRCan), GreenGrass, KoFlux, LBA, NECC, OzFlux, TCOS-Siberia and USCCC. We acknowledge the financial support to the eddy covariance data harmonisation provided by CarboEuropeIP, FAO-GTOS-TCO, iLEAPS, Max Planck Institute for Biogeochemistry, National Science Foundation, University of Tuscia, Universiteì Laval, Environment Canada and US Department of Energy and the database development and technical support from Berkeley Water Center, Lawrence Berkeley National Laboratory, Microsoft Research eScience, Oak Ridge National Laboratory, University of California – Berkeley and the University of Virginia.Peer reviewedPublisher PD

Constraining a global ecosystem model with multi-site eddy-covariance data

Peer reviewedPublisher PD

Improving Estimates of Gross Primary Productivity by Assimilating Solar-Induced Fluorescence Satellite Retrievals in a Terrestrial Biosphere Model Using a Process-Based SIF Model

Abstract Over the last few years, solar-induced chlorophyll fluorescence (SIF) observations from space have emerged as a promising resource for evaluating the spatio-temporal distribution of gross primary productivity (GPP) simulated by global terrestrial biosphere models. SIF can be used to improve GPP simulations by optimizing critical model parameters through statistical Bayesian data assimilation techniques. A prerequisite is the availability of a functional link between GPP and SIF in terrestrial biosphere models. Here we present the development of a mechanistic SIF observation operator in the ORCHIDEE (Organizing Carbon and Hydrology In Dynamic Ecosystems) terrestrial biosphere model. It simulates the regulation of photosystem II fluorescence quantum yield at the leaf level thanks to a novel parameterization of non-photochemical quenching as a function of temperature, photosynthetically active radiation, and normalized quantum yield of photochemistry. It emulates the radiative transfer of chlorophyll fluorescence to the top of the canopy using a parametric simplification of the SCOPE (Soil Canopy Observation Photosynthesis Energy) model. We assimilate two years of monthly OCO-2 (Orbiting Carbon Observatory-2) SIF product at 0.5° (2015?2016) to optimize ORCHIDEE photosynthesis and phenological parameters over an ensemble of grid points for all plant functional types. The impact on the simulated GPP is considerable with a large decrease of the global scale budget by 28 GtC/year over the period 1990?2009. The optimized GPP budget (134/136 GtC/year over 1990?2009/2001?2009) remarkably agrees with independent GPP estimates, FLUXSAT (137 GtC/year over 2001?2009) in particular and FLUXCOM (121 GtC/year over 1990?2009). Our results also suggest a biome dependency of the SIF-GPP relationship that needs to be improved for some plant functional types.Peer reviewe

A new stepwise carbon cycle data assimilation system using multiple data streams to constrain the simulated land surface carbon cycle

Acknowledgements. This work was mainly funded by the EU FP7 CARBONES project (contracts FP7-SPACE-2009-1-242316), with also a small contribution from GEOCARBON project (ENV.2011.4.1.1-1-283080). This work used eddy covariance data acquired by the FLUXNET community and in particular by the following networks: AmeriFlux (U.S. Department of Energy, Biological and Environmental Research, Terrestrial Carbon Program; DE-FG02-04ER63917 and DE-FG02-04ER63911), AfriFlux, AsiaFlux, CarboAfrica, CarboEuropeIP, CarboItaly, CarboMont, ChinaFlux, Fluxnet-Canada (supported by CFCAS, NSERC, BIOCAP, Environment Canada, and NRCan), GreenGrass, KoFlux, LBA, NECC, OzFlux, TCOS-Siberia, USCCC. We acknowledge the financial support to the eddy covariance data harmonization provided by CarboEuropeIP, FAO-GTOS-TCO, iLEAPS, Max Planck Institute for Biogeochemistry, National Science Foundation, University of Tuscia, Université Laval and Environment Canada and US Department of Energy and the database development and technical support from Berkeley Water Center, Lawrence Berkeley National Laboratory, Microsoft Research eScience, Oak Ridge National Laboratory, University of California-Berkeley, University of Virginia. Philippe Ciais acknowledges support from the European Research Council through Synergy grant ERC-2013-SyG-610028 “IMBALANCE-P”. The authors wish to thank M. Jung for providing access to the GPP MTE data, which were downloaded from the GEOCARBON data portal (https://www.bgc-jena.mpg.de/geodb/projects/Data.php). The authors are also grateful to computing support and resources provided at LSCE and to the overall ORCHIDEE project that coordinate the development of the code (http://labex.ipsl.fr/orchidee/index.php/about-the-team).Peer reviewedPublisher PD

Land surface model parameter optimisation using in situ flux data: comparison of gradient-based versus random search algorithms (a case study using ORCHIDEE v1.9.5.2)

Land surface models (LSMs), which form the land component of earth system models, rely on numerous processes for describing carbon, water and energy budgets, often associated with highly uncertain parameters. Data assimilation (DA) is a useful approach for optimising the most critical parameters in order to improve model accuracy and refine future climate predictions. In this study, we compare two different DA methods for optimising the parameters of seven plant functional types (PFTs) of the ORCHIDEE LSM using daily averaged eddy-covariance observations of net ecosystem exchange and latent heat flux at 78 sites across the globe. We perform a technical investigation of two classes of minimisation methods – local gradient-based (the L-BFGS-B algorithm, limited memory Broyden–Fletcher–Goldfarb–Shanno algorithm with bound constraints) and global random search (the genetic algorithm) – by evaluating their relative performance in terms of the model–data fit and the difference in retrieved parameter values. We examine the performance of each method for two cases: when optimising parameters at each site independently (“single-site” approach) and when simultaneously optimising the model at all sites for a given PFT using a common set of parameters (“multi-site” approach). We find that for the single site case the random search algorithm results in lower values of the cost function (i.e. lower model–data root mean square differences) than the gradient-based method; the difference between the two methods is smaller for the multi-site optimisation due to a smoothing of the cost function shape with a greater number of observations. The spread of the cost function, when performing the same tests with 16 random first-guess parameters, is much larger with the gradient-based method, due to the higher likelihood of being trapped in local minima. When using pseudo-observation tests, the genetic algorithm results in a closer approximation of the true posterior parameter value in the L-BFGS-B algorithm. We demonstrate the advantages and challenges of different DA techniques and provide some advice on using it for the LSM parameter optimisation.</p

Vegetation Cover Analysis of Hazardous Waste Sites in Utah and Arizona Using Hyperspectral Remote Sensing

This study investigated the usability of hyperspectral remote sensing for characterizing vegetation at hazardous waste sites. The specific objectives of this study were to: (1) estimate leaf-area-index (LAI) of the vegetation using three different methods (i.e., vegetation indices, red-edge positioning (REP), and machine learning regression trees), and (2) map the vegetation cover using machine learning decision trees based on either the scaled reflectance data or mixture tuned matched filtering (MTMF)-derived metrics and vegetation indices. HyMap airborne data (126 bands at 2.3 x 2.3 m spatial resolution), collected over the U. S. Department of Energy uranium processing sites near Monticello, Utah and Monument Valley, Arizona, were used. Grass and shrub species were mixed on an engineered disposal cell cover at the Monticello site while shrub species were dominant in the phytoremediation plantings at the Monument Valley site. Regression trees resulted in the best calibration performance of LAI estimation (R-2 > 0.80. The use of REPs failed to accurately predict LAI (R-2 < 0.2). The use of the MTMF-derived metrics (matched filter scores and infeasibility) and a range of vegetation indices in decision trees improved the vegetation mapping when compared to the decision tree classification using just the scaled reflectance. Results suggest that hyperspectral imagery are useful for characterizing biophysical characteristics (LAI) and vegetation cover on capped hazardous waste sites. However, it is believed that the vegetation mapping would benefit from the use of higher spatial resolution hyperspectral data due to the small size of many of the vegetation patches (<1 m) found on the sites.open111

Investigating the role of prior and observation error correlations in improving a model forecast of forest carbon balance using Four Dimensional Variational data assimilation

Efforts to implement variational data assimilation routines with functional ecology models and land surface models have been limited, with sequential and Markov chain Monte Carlo data assimilation methods being prevalent. When data assimilation has been used with models of carbon balance, prior or “background” errors (in the initial state and parameter values) and observation errors have largely been treated as independent and uncorrelated. Correlations between background errors have long been known to be a key aspect of data assimilation in numerical weather prediction. More recently, it has been shown that accounting for correlated observation errors in the assimilation algorithm can considerably improve data assimilation results and forecasts. In this paper we implement a Four-Dimensional Variational data assimilation (4D-Var) scheme with a simple model of forest carbon balance, for joint parameter and state estimation and assimilate daily observations of Net Ecosystem CO2 Exchange (NEE) taken at the Alice Holt forest CO2 flux site in Hampshire, UK. We then investigate the effect of specifying correlations between parameter and state variables in background error statistics and the effect of specifying correlations in time between observation errors. The idea of including these correlations in time is new and has not been previously explored in carbon balance model data assimilation. In data assimilation, background and observation error statistics are often described by the background error covariance matrix and the observation error covariance matrix. We outline novel methods for creating correlated versions of these matrices, using a set of previously postulated dynamical constraints to include correlations in the background error statistics and a Gaussian correlation function to include time correlations in the observation error statistics. The methods used in this paper will allow the inclusion of time correlations between many different observation types in the assimilation algorithm, meaning that previously neglected information can be accounted for. In our experiments we assimilate a single year of NEE observations and then run a forecast for the next 14 years. We compare the results using our new correlated background and observation error covariance matrices and those using diagonal covariance matrices. We find that using the new correlated matrices reduces the root mean square error in the 14 year forecast of daily NEE by 44% decreasing from 4.22 gCm−2 day−1 to 2.38 gCm−2 day−

Fusing optical and SAR time series for LAI gap filling with multioutput Gaussian processes

The availability of satellite optical information is often hampered by the natural presence of clouds, which can be problematic for many applications. Persistent clouds over agricultural fields can mask key stages of crop growth, leading to unreliable yield predictions. Synthetic Aperture Radar (SAR) provides all-weather imagery which can potentially overcome this limitation, but given its high and distinct sensitivity to different surface properties, the fusion of SAR and optical data still remains an open challenge. In this work, we propose the use of Multi-Output Gaussian Process (MOGP) regression, a machine learning technique that learns automatically the statistical relationships among multisensor time series, to detect vegetated areas over which the synergy between SAR-optical imageries is profitable. For this purpose, we use the Sentinel-1 Radar Vegetation Index (RVI) and Sentinel-2 Leaf Area Index (LAI) time series over a study area in north west of the Iberian peninsula. Through a physical interpretation of MOGP trained models, we show its ability to provide estimations of LAI even over cloudy periods using the information shared with RVI, which guarantees the solution keeps always tied to real measurements. Results demonstrate the advantage of MOGP especially for long data gaps, where optical-based methods notoriously fail. The leave-one-image-out assessment technique applied to the whole vegetation cover shows MOGP predictions improve standard GP estimations over short-time gaps (R 2 of 74% vs 68%, RMSE of 0.4 vs 0.44 [m 2 m −2 ]) and especially over long-time gaps (R 2 of 33% vs 12%, RMSE of 0.5 vs 1.09 [m 2 m −2 ])