Geophysical time series often feature missing data or data acquired at irregular times. Procedures are needed to either resample these series at systematic time intervals or to generate reasonable ...estimates at specified times in order to meet specific user requirements or to facilitate subsequent analyses. Interpolation methods have long been used to address this problem, taking into account the fact that available measurements also include errors of measurement or uncertainties. This paper inspects some of the currently used approaches to fill gaps and smooth time series (smoothing splines, Singular Spectrum Analysis and Lomb-Scargle) by comparing their performance in either reconstructing the original record or in minimizing the Mean Absolute Error (MAE), Mean Bias Error (MBE), chi-squared test statistics and autocorrelation of residuals between the underlying model and the available data, using both artificially-generated series or well-known publicly available records. Some methods make no assumption on the type of variability in the data while others hypothesize the presence of at least some dominant frequencies. It will be seen that each method exhibits advantages and drawbacks, and that the choice of an approach largely depends on the properties of the underlying time series and the objective of the research.
Monitoring biosphere vegetation 1998-2009 Gobron, N.; Belward, A.; Pinty, B. ...
Geophysical research letters,
August 2010, Letnik:
37, Številka:
15
Journal Article
Recenzirano
Odprti dostop
Earth Observation from space offers the opportunity to produce time‐series of geophysical products that can be used to assess the state and changes of land surfaces. The Fraction of Absorbed ...Photosynthetically Active Radiation (FAPAR) is used to monitor the state and evolution of terrestrial vegetation, and also constitutes a state variable in advanced Earth system models that contain a detailed enough description of the terrestrial biosphere. This present study reports a 12‐year (1998–2009) time series of FAPAR derived from the combination of two satellite‐based sensors. We find that FAPAR exhibits large‐scale inter‐annual variations and multi‐year trends. The fraction of land grid cells showing positive anomalies, as computed by the deviation from the 12‐year climatology, shows a rapid decrease in the early part of the analysis period (until 2004). Large negative anomalies can be associated with previously reported large‐scale climate events, such as global land drying associated with El Niño Southern Oscillation 2000–2003, or the European drought of 2003 or recent Australian droughts The present analysis demonstrates that FAPAR is an important global variable suitable for large‐scale monitoring of climate impacts on the terrestrial biosphere.
Drought and ecosystem carbon cycling van der Molen, M.K.; Dolman, A.J.; Ciais, P. ...
Agricultural and forest meteorology,
07/2011, Letnik:
151, Številka:
7
Journal Article
Recenzirano
► Intermittent droughts interact differently with the carbon cycle than gradual change. ► Direct effects on carbon exchange depend strongly on vegetation species strategie. ► Drought induced ...mortality is species sensitive and gives rise to competition. ► Disturbed nutrient and carbohydrate reservoirs cause carry-over effects. ► Models need to simulate these longer-term climate-carbon cycle feedbacks.
Drought as an intermittent disturbance of the water cycle interacts with the carbon cycle differently than the ‘gradual’ climate change. During drought plants respond physiologically and structurally to prevent excessive water loss according to species-specific water use strategies. This has consequences for carbon uptake by photosynthesis and release by total ecosystem respiration. After a drought the disturbances in the reservoirs of moisture, organic matter and nutrients in the soil and carbohydrates in plants lead to longer-term effects in plant carbon cycling, and potentially mortality. Direct and carry-over effects, mortality and consequently species competition in response to drought are strongly related to the survival strategies of species. Here we review the state of the art of the understanding of the relation between soil moisture drought and the interactions with the carbon cycle of the terrestrial ecosystems. We argue that plant strategies must be given an adequate role in global vegetation models if the effects of drought on the carbon cycle are to be described in a way that justifies the interacting processes.
Current methods for estimating vegetation parameters are generally sub-optimal in the way they exploit information and do not generally consider uncertainties. We look forward to a future where ...operational data assimilation schemes improve estimates by tracking land surface processes and exploiting multiple types of observations. Data assimilation schemes seek to combine observations and models in a statistically optimal way taking into account uncertainty in both, but have not yet been much exploited in this area.
The EO-LDAS scheme and prototype, developed under ESA funding, is designed to exploit the anticipated wealth of data that will be available under GMES missions, such as the Sentinel family of satellites, to provide improved mapping of land surface biophysical parameters. This paper describes the EO-LDAS implementation, and explores some of its core functionality. EO-LDAS is a weak constraint variational data assimilation system. The prototype provides a mechanism for constraint based on a prior estimate of the state vector, a linear dynamic model, and Earth Observation data (top-of-canopy reflectance here). The observation operator is a non-linear optical radiative transfer model for a vegetation canopy with a soil lower boundary, operating over the range 400 to 2500nm. Adjoint codes for all model and operator components are provided in the prototype by automatic differentiation of the computer codes.
In this paper, EO-LDAS is applied to the problem of daily estimation of six of the parameters controlling the radiative transfer operator over the course of a year (>2000 state vector elements). Zero and first order process model constraints are implemented and explored as the dynamic model. The assimilation estimates all state vector elements simultaneously. This is performed in the context of a typical Sentinel-2 MSI operating scenario, using synthetic MSI observations simulated with the observation operator, with uncertainties typical of those achieved by optical sensors supposed for the data.
The experiments consider a baseline state vector estimation case where dynamic constraints are applied, and assess the impact of dynamic constraints on the a posteriori uncertainties. The results demonstrate that reductions in uncertainty by a factor of up to two might be obtained by applying the sorts of dynamic constraints used here. The hyperparameter (dynamic model uncertainty) required to control the assimilation are estimated by a cross-validation exercise. The result of the assimilation is seen to be robust to missing observations with quite large data gaps.
► Monitoring of terrestrial ecosystems using Earth Observation data is often ill-posed. ► This can be improved by Data Assimilation with biophysical/empirical process models. ► A prototype DA system for land surface has been developed: EO-LDAS. ► With synthetic Sentinel2 MSI sampling we show improvements by a factor of 2 using DA. ► These results are quite robust to data gaps from clouds.
A globally integrated carbon observation and analysis system is needed to improve the fundamental understanding of the global carbon cycle, to improve our ability to project future changes, and to ...verify the effectiveness of policies aiming to reduce greenhouse gas emissions and increase carbon sequestration. Building an integrated carbon observation system requires transformational advances from the existing sparse, exploratory framework towards a dense, robust, and sustained system in all components: anthropogenic emissions, the atmosphere, the ocean, and the terrestrial biosphere. The paper is addressed to scientists, policymakers, and funding agencies who need to have a global picture of the current state of the (diverse) carbon observations. We identify the current state of carbon observations, and the needs and notional requirements for a global integrated carbon observation system that can be built in the next decade. A key conclusion is the substantial expansion of the ground-based observation networks required to reach the high spatial resolution for CO2 and CH4 fluxes, and for carbon stocks for addressing policy-relevant objectives, and attributing flux changes to underlying processes in each region. In order to establish flux and stock diagnostics over areas such as the southern oceans, tropical forests, and the Arctic, in situ observations will have to be complemented with remote-sensing measurements. Remote sensing offers the advantage of dense spatial coverage and frequent revisit. A key challenge is to bring remote-sensing measurements to a level of long-term consistency and accuracy so that they can be efficiently combined in models to reduce uncertainties, in synergy with ground-based data. Bringing tight observational constraints on fossil fuel and land use change emissions will be the biggest challenge for deployment of a policy-relevant integrated carbon observation system. This will require in situ and remotely sensed data at much higher resolution and density than currently achieved for natural fluxes, although over a small land area (cities, industrial sites, power plants), as well as the inclusion of fossil fuel CO2 proxy measurements such as radiocarbon in CO2 and carbon-fuel combustion tracers. Additionally, a policy-relevant carbon monitoring system should also provide mechanisms for reconciling regional top-down (atmosphere-based) and bottom-up (surface-based) flux estimates across the range of spatial and temporal scales relevant to mitigation policies. In addition, uncertainties for each observation data-stream should be assessed. The success of the system will rely on long-term commitments to monitoring, on improved international collaboration to fill gaps in the current observations, on sustained efforts to improve access to the different data streams and make databases interoperable, and on the calibration of each component of the system to agreed-upon international scales.
The terrestrial biosphere is currently a strong sink for anthropogenic CO2 emissions. Through the radiative properties of CO2 , the strength of this sink has a direct influence on the radiative ...budget of the global climate system. The accurate assessment of this sink and its evolution under a changing climate is, hence, paramount for any efficient management strategies of the terrestrial carbon sink to avoid dangerous climate change. Unfortunately, simulations of carbon and water fluxes with terrestrial biosphere models exhibit large uncertainties. A considerable fraction of this uncertainty reflects uncertainty in the parameter values of the process formulations within the models. This paper describes the systematic calibration of the process parameters of a terrestrial biosphere model against two observational data streams: remotely sensed FAPAR (fraction of absorbed photosynthetically active radiation) provided by the MERIS (ESA's Medium Resolution Imaging Spectrometer) sensor and in situ measurements of atmospheric CO2 provided by the GLOBALVIEW flask sampling network. We use the Carbon Cycle Data Assimilation System (CCDAS) to systematically calibrate some 70 parameters of the terrestrial BETHY (Biosphere Energy Transfer Hydrology) model. The simultaneous assimilation of all observations provides parameter estimates and uncertainty ranges that are consistent with the observational information. In a subsequent step these parameter uncertainties are propagated through the model to uncertainty ranges for predicted carbon fluxes. We demonstrate the consistent assimilation at global scale, where the global MERIS FAPAR product and atmospheric CO2 are used simultaneously. The assimilation improves the match to independent observations. We quantify how MERIS data improve the accuracy of the current and future (net and gross) carbon flux estimates (within and beyond the assimilation period). We further demonstrate the use of an interactive mission benefit analysis tool built around CCDAS to support the design of future space missions. We find that, for long-term averages, the benefit of FAPAR data is most pronounced for hydrological quantities, and moderate for quantities related to carbon fluxes from ecosystems. The benefit for hydrological quantities is highest for semi-arid tropical or sub-tropical regions. Length of mission or sensor resolution is of minor importance.
Terrestrial productivity in semi-arid woodlands is strongly susceptible to changes in precipitation, and semi-arid woodlands constitute an important element of the global water and carbon cycles. ...Here, we use the Carbon Cycle Data Assimilation System (CCDAS) to investigate the key parameters controlling ecological and hydrological activities for a semi-arid savanna woodland site in Maun, Botswana. Twenty-four eco-hydrological process parameters of a terrestrial ecosystem model are optimized against two data streams separately and simultaneously: daily averaged latent heat flux (LHF) derived from eddy covariance measurements, and decadal fraction of absorbed photosynthetically active radiation (FAPAR) derived from the Sea-viewing Wide Field-of-view Sensor (SeaWiFS). Assimilation of both data streams LHF and FAPAR for the years 2000 and 2001 leads to improved agreement between measured and simulated quantities not only for LHF and FAPAR, but also for photosynthetic CO2 uptake. The mean uncertainty reduction (relative to the prior) over all parameters is 14.9% for the simultaneous assimilation of LHF and FAPAR, 8.5% for assimilating LHF only, and 6.1% for assimilating FAPAR only. The set of parameters with the highest uncertainty reduction is similar between assimilating only FAPAR or only LHF. The highest uncertainty reduction for all three cases is found for a parameter quantifying maximum plant-available soil moisture. This indicates that not only LHF but also satellite-derived FAPAR data can be used to constrain and indirectly observe hydrological quantities.
Long-term precipitation chemistry have been recorded in the rural area of Banizoumbou (Niger), representative of a semi-arid savanna ecosystem. A total of 305 rainfall samples ~90% of the total ...annual rainfall) were collected from June 1994 to September 2005. From ionic chromatography, pH major inorganic and organic ions were detected. Rainwater chemistry is controlled by soil/dust emissions associated with terrigeneous elements represented by SO42−, Ca2+, Carbonates, K+ and Mg2+. It is found that calcium and carbonates represent ~40% of the total ionic charge. The second highest contribution is nitrogenous, with annual Volume Weighed Mean (VWM) for NO3− and NH4+ concentrations of 11.6 and 18.1 μeq.l−1, respectively. This is the signature of ammonia sources from animals and NOx emissions from savannas soil-particles rain-induced. The mean annual NH3 and NO2 air concentration are of 6 ppbv and 2.6 ppbv, respectively. The annual VWM precipitation concentration of sodium and chloride are both of 8.7 μeq.l−1 which reflects the marine signature of monsoonal and humid air masses. The median pH value is of 6.05. Acidity is neutralized by mineral dust, mainly carbonates, and/or dissolved gases such NH3. High level of organic acidity with 8μeq.l−1 and 5.2 μeq.l−1 of formate and acetate were also found. The analysis of monthly Black Carbon emissions and Fraction of Absorbed Photosynthetically Active Radiation (FAPAR) values show that both biogenic emission from vegetation and biomass burning could explain the rainfall organic acidity content. The interannual variability of the VWM concentrations around the mean (1994–2005) is between ±5% and ±30% and mainly due to variations of sources strength and rainfall spatio-temporal distribution. From 1994 to 2005, the total mean wet deposition flux in the Sahelian region is of 60.1 mmol.m−2.yr−1 ±25%. Finally, Banizoumbou measurements are compared to other long-term measurements of precipitation chemistry in the wet savanna of Lamto (Côte d'Ivoire) and in the forested zone of Zoétélé (Cameroon). The total chemical loading presents a maximum in the dry savanna and a minimum in the forest (from 143.7, 100.2 to 86.6 μeq.l−1), associated with the gradient of terrigeneous sources. The wet deposition fluxes present an opposite trend, with 60.0 mmol.m−2.yr−1 in Banizoumbou, 108.6 mmol.m−2.yr−1 in Lamto and 162.9 mmol.m−2.yr−1 in Zoétélé, controlled by rainfall gradient along the ecosystems transect.
Photosynthesis by terrestrial plants is the main driver of the global carbon cycle, and the presence of actively photosynthesizing vegetation can now be observed from space. However, challenges ...remain when translating remotely sensed data into carbon fluxes. One reason is that the Fraction of Absorbed Photosynthetically Active Radiation (FAPAR), which documents the presence of photosynthetically active vegetation, relates more directly to leaf development and leaf phenology than to photosynthetic rates. Here, we present a new approach for linking FAPAR and vegetation‐to‐atmosphere carbon fluxes through variational data assimilation. The scheme extends the Carbon Cycle Data Assimilation System (CCDAS) by a newly developed, globally applicable and generic leaf phenology model, which includes both temperature and water‐driven leaf development. CCDAS is run for seven sites, six of them included in the FLUXNET network. Optimization is carried out simultaneously for all sites against 20 months of daily FAPAR from the Medium Resolution Imaging Spectrometer on board the European Space Agency's ENVISAT platform. Fourteen parameters related to phenology and 24 related to photosynthesis are optimized simultaneously, and their posterior uncertainties are computed. We find that with one parameter set for all sites, the model is able to reproduce the observed FAPAR spanning boreal, temperate, humid‐tropical, and semiarid climates. Assimilation of FAPAR has led to reduced uncertainty (by >10%) of 10 of the 38 parameters, including one parameter related to photosynthesis, and a moderate reduction in net primary productivity uncertainty. The approach can easily be extended to regional or global studies and to the assimilation of further remotely sensed data.
This contribution illustrates results from a large‐scale application of the Joint Research Centre Two‐stream Inversion Package (JRC‐TIP), using MODIS broadband visible and near‐infrared white sky ...surface albedos as inputs. The discussion focuses on products (based on the mean and one‐sigma values of the probability distribution functions (PDFs)) obtained during the summer and winter. This paper discusses the retrieved model parameters including the effective leaf area index (LAI), the background brightness, and the scattering efficiency of the vegetation elements. The similarity between the derived LAI seasonal maps and earlier distributions of this variable comforts us in the quality of the albedo products as well as in the ability of the JRC‐TIP to interpret the latter meaningfully. The opportunity to generate global maps of new products, such as the background albedo, underscores the advantages of using state of the art algorithmic approaches capable of fully exploiting accurate satellite remote sensing data sets. The detailed analyses of the retrieval uncertainties highlight the central role and contribution of the LAI, the main process parameter to interpret radiation transfer observations over vegetated surfaces. The estimation of the radiation fluxes that are absorbed, transmitted, and scattered by the vegetation layer and its background is achieved on the basis of the retrieved PDFs of the model parameters. Results from this latter step are discussed in a companion paper.
Key Points
Consistent estimates of model process parameters from Global MODIS albedos
Novel inversion techniques based on optimal control theory
Remarkable spatiotemporal consistency of the retrievals with uncertainties