Explore the vast range of titles available.

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.

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.

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.

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.

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

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.

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 Photsynthetically Active Radiation (FAPAR), which documents the presence of photsynthetically 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 (MERIS) on-board the European Space Agency s ENVISAT platform. 14 parameters related to phenology and 24 related to photosynthesis are optimized simultaneously and their posterior uncertainties 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 semi-arid climates. Assimilation of FAPAR has led to reduced uncertainty (>10%) of 10 of the 38 parameters, including 1 parameter related to photosynthesis, and a moderate reduction in NPP uncertainty. The approach can easily be extended to regional or global studies and to the assimilation of further remotely sensed data sources.

Croplands cover about 12% of the ice-free terrestrial land surface. Compared with natural ecosystems, croplands have distinct characteristics due to anthropogenic influences. Their global gross primary production (GPP) is not well constrained and estimates vary between 8.2 and 14.2 Pg C yr−1. We quantified global cropland GPP using a light use efficiency (LUE) model, employing satellite observations and survey data of crop types and distribution. A novel step in our analysis was to assign a maximum light use efficiency estimate (ϵ*GPP) to each of the 26 different crop types, instead of taking a uniform value as done in the past. These ϵ*GPP values were calculated based on flux tower CO2 exchange measurements and a literature survey of field studies, and ranged from 1.20 g CMJ−1 to 2.96 g CMJ−1. Global cropland GPP was estimated to be 11.05 Pg C yr−1 in the year 2000. Maize contributed most to this (1.55 Pg C yr−1), and the continent of Asia contributed most with 38.9% of global cropland GPP. In the continental United States, annual cropland GPP (1.28 Pg C yr−1) was close to values reported previously (1.24 Pg C yr−1) constrained by harvest records, but our estimates of ϵ*GPP values were much higher. Our results are sensitive to satellite information and survey data on crop type and extent, but provide a consistent and data-driven approach to generate a look-up table of ϵ*GPP for the 26 crop types for potential use in other vegetation models.

The two-stream model parameters and associated uncertainties retrieved by inversion against MODIS broadband visible and near-infrared white sky surface albedos were discussed in a companion paper (Pinty et al., 2010). The present paper concentrates on the partitioning of the solar radiation fluxes delivered by the Joint Research Centre Two-stream Inversion Package (JRC-TIP). The estimation of the various flux fractions related to the vegetation and the background layers separately capitalizes on the PDFs of the model parameters discussed in the companion paper. The propagation of uncertainties from the observations to the model parameters is achieved via the Hessian of the cost function and yields a covariance matrix of posterior parameter uncertainties. This matrix is propagated to the radiation fluxes via the model's Jacobian matrix of first derivatives. Results exhibit a rather good spatio-temporal consistency given that the prior values on the model parameters are not specified as a function of land cover type and/or vegetation phenological states. A specific investigation based on a scenario imposing stringent conditions of leaf absorbing and scattering properties highlights the impact of such constraints that are, as a matter of fact, currently adopted in vegetation index approaches. Special attention is also given to snow covered and snow contaminated areas since these regions encompass significant reflectance changes that strongly affect land surface processes. A definite asset of the JRC-TIP lies in its capability to control and ultimately relax a number of assumptions that are often implicit in traditional approaches. These features greatly help us understand the discrepancies between the different data sets of land surface properties and fluxes that are currently available. Through a series of selected examples, the inverse procedure implemented in the JRC-TIP is shown to be robust, reliable and compliant with large scale processing requirements. Furthermore, this package ensures the physical consistency between the set of observations, the two-stream model parameters and radiation fluxes. It also documents the retrieval of associated uncertainties.