José Gómez-Dans

Defining pyromes and global syndromes of fire regimes

Fire is a ubiquitous component of the Earth system that is poorly understood. To date, a global-scale understanding of fire is largely limited to the annual extent of burning as detected by satellites. This is problematic because fire is multidimensional, and focus on a single metric belies its complexity and importance within the Earth system. To address this, we identified five key characteristics of fire regimes—size, frequency, intensity, season, and extent—and combined new and existing global datasets to represent each. We assessed how these global fire regime characteristics are related to patterns of climate, vegetation (biomes), and human activity. Cross-correlations demonstrate that only certain combinations of fire characteristics are possible, reflecting fundamental constraints in the types of fire regimes that can exist. A Bayesian clustering algorithm identified five global syndromes of fire regimes, or pyromes. Four pyromes represent distinctions between crown, litter, and grass-fueled fires, and the relationship of these to biomes and climate are not deterministic. Pyromes were partially discriminated on the basis of available moisture and rainfall seasonality. Human impacts also affected pyromes and are globally apparent as the driver of a fifth and unique pyrome that represents human-engineered modifications to fire characteristics. Differing biomes and climates may be represented within the same pyrome, implying that pathways of change in future fire regimes in response to changes in climate and human activity may be difficult to predict.

An Emulator Toolbox to Approximate Radiative Transfer Models with Statistical Learning

Physically-based radiative transfer models (RTMs) help in understanding the processes occurring on the Earth’s surface and their interactions with vegetation and atmosphere. When it comes to studying vegetation properties, RTMs allows us to study light interception by plant canopies and are used in the retrieval of biophysical variables through model inversion. However, advanced RTMs can take a long computational time, which makes them unfeasible in many real applications. To overcome this problem, it has been proposed to substitute RTMs through so-called emulators. Emulators are statistical models that approximate the functioning of RTMs. Emulators are advantageous in real practice because of the computational efficiency and excellent accuracy and flexibility for extrapolation. We hereby present an “Emulator toolbox” that enables analysing multi-output machine learning regression algorithms (MO-MLRAs) on their ability to approximate an RTM. The toolbox is included in the free-access ARTMO’s MATLAB suite for parameter retrieval and model inversion and currently contains both linear and non-linear MO-MLRAs, namely partial least squares regression (PLSR), kernel ridge regression (KRR) and neural networks (NN). These MO-MLRAs have been evaluated on their precision and speed to approximate the soil vegetation atmosphere transfer model SCOPE (Soil Canopy Observation, Photochemistry and Energy balance). SCOPE generates, amongst others, sun-induced chlorophyll fluorescence as the output signal. KRR and NN were evaluated as capable of reconstructing fluorescence spectra with great precision. Relative errors fell below 0.5% when trained with 500 or more samples using cross-validation and principal component analysis to alleviate the underdetermination problem. Moreover, NN reconstructed fluorescence spectra about 50-times faster and KRR about 800-times faster than SCOPE. The Emulator toolbox is foreseen to open new opportunities in the use of advanced RTMs, in which both consistent physical assumptions and data-driven machine learning algorithms live together.

Indoor C-band polarimetric interferometry observations of a mature wheat canopy

We present results from experiments carried out in the ground-based synthetic aperture radar (GB-SAR) facility at the University of Sheffield to ascertain the role of polarimetric interferometry in crop height retrieval. To this end, a mature wheat canopy, grown in outdoor conditions, was reassembled inside the GB-SAR chamber and imaged at C-band using a two-dimensional scan. This allowed fully polarimetric tomography and interferometry. Interferometry using the VV, HH, and VH polarization states shows that the HH and VH interferograms retrieve a height close to the top of the soil layer for all angles of incidence considered, whereas the height retrieved from the VV interferogram increases with angle of incidence. The use of a Pauli basis gives poor results, due to the different location of the scattering phase centers in the VV and HH channels. The use of arbitrary polarization states shows that the top of the soil can be very accurately estimated using left-circular polarization, whereas, for angles of incidence close to 45/spl deg/, a polarization state similar to VV can be used to retrieve the top of the canopy; hence crop height can be recovered as the difference of these two interferometric heights. Polarimetric coherence optimization techniques are also studied. Unconstrained coherence optimization gives very unstable results, due to the small number of available samples. Constrained optimization results in stable retrieved heights, and the retrieved polarization states agree well with the polarization synthesis results.

A Survey on Gaussian Processes for Earth-Observation Data Analysis: A Comprehensive Investigation

Gaussian processes (GPs) have experienced tremendous success in biogeophysical parameter retrieval in the last few years. GPs constitute a solid Bayesian framework to consistently formulate many function approximation problems. This article reviews the main theoretical GP developments in the field, considering new algorithms that respect signal and noise characteristics, extract knowledge via automatic relevance kernels to yield feature rankings automatically, and allow applicability of associated uncertainty intervals to transport GP models in space and time that can be used to uncover causal relations between variables and can encode physically meaningful prior knowledge via radiative transfer model (RTM) emulation. The important issue of computational efficiency will also be addressed. These developments are illustrated in the field of geosciences and remote sensing at local and global scales through a set of illustrative examples. In particular, important problems for land, ocean, and atmosphere monitoring are considered, from accurately estimating oceanic chlorophyll content and pigments to retrieving vegetation properties from multi- and hyperspectral sensors as well as estimating atmospheric parameters (e.g., temperature, moisture, and ozone) from infrared sounders.

African Savanna-Forest Boundary Dynamics: A 20-Year Study

Recent studies show widespread encroachment of forest into savannas with important consequences for the global carbon cycle and land-atmosphere interactions. However, little research has focused on in situ measurements of the successional sequence of savanna to forest in Africa. Using long-term inventory plots we quantify changes in vegetation structure, above-ground biomass (AGB) and biodiversity of trees ≥10 cm diameter over 20 years for five vegetation types: savanna; colonising forest (F1), monodominant Okoume forest (F2); young Marantaceae forest (F3); and mixed Marantaceae forest (F4) in Lopé National Park, central Gabon, plus novel 3D terrestrial laser scanning (TLS) measurements to assess forest structure differences. Over 20 years no plot changed to a new stage in the putative succession, but F1 forests strongly moved towards the structure, AGB and diversity of F2 forests. Overall, savanna plots showed no detectable change in structure, AGB or diversity using this method, with zero trees ≥10 cm diameter in 1993 and 2013. F1 and F2 forests increased in AGB, mainly as a result of adding recruited stems (F1) and increased Basal Area (F2), whereas F3 and F4 forests did not change substantially in structure, AGB or diversity. Critically, the stability of the F3 stage implies that this stage may be maintained for long periods. Soil carbon was low, and did not show a successional gradient as for AGB and diversity. TLS vertical plant profiles showed distinctive differences amongst the vegetation types, indicating that this technique can improve ecological understanding. We highlight two points: (i) as forest colonises, changes in biodiversity are much slower than changes in forest structure or AGB; and (ii) all forest types store substantial quantities of carbon. Multi-decadal monitoring is likely to be required to assess the speed of transition between vegetation types.

Strong constraint on modelled global carbon uptake using solar-induced chlorophyll fluorescence data

Accurate terrestrial biosphere model (TBM) simulations of gross carbon uptake (gross primary productivity – GPP) are essential for reliable future terrestrial carbon sink projections. However, uncertainties in TBM GPP estimates remain. Newly-available satellite-derived sun-induced chlorophyll fluorescence (SIF) data offer a promising direction for addressing this issue by constraining regional-to-global scale modelled GPP. Here, we use monthly 0.5° GOME-2 SIF data from 2007 to 2011 to optimise GPP parameters of the ORCHIDEE TBM. The optimisation reduces GPP magnitude across all vegetation types except C4 plants. Global mean annual GPP therefore decreases from 194 ± 57 PgCyr−1 to 166 ± 10 PgCyr−1, bringing the model more in line with an up-scaled flux tower estimate of 133 PgCyr−1. Strongest reductions in GPP are seen in boreal forests: the result is a shift in global GPP distribution, with a ~50% increase in the tropical to boreal productivity ratio. The optimisation resulted in a greater reduction in GPP than similar ORCHIDEE parameter optimisation studies using satellite-derived NDVI from MODIS and eddy covariance measurements of net CO2 fluxes from the FLUXNET network. Our study shows that SIF data will be instrumental in constraining TBM GPP estimates, with a consequent improvement in global carbon cycle projections.

Estimation of FAPAR over Croplands Using MISR Data and the Earth Observation Land Data Assimilation System (EO-LDAS)

The Fraction of Absorbed Photosynthetically-Active Radiation (FAPAR) is an important parameter in climate and carbon cycle studies. In this paper, we use the Earth Observation Land Data Assimilation System (EO-LDAS) framework to retrieve FAPAR from observations of directional surface reflectance measurements from the Multi-angle Imaging SpectroRadiometer(MISR) instrument. The procedure works by interpreting the reflectance data via the semi-discrete Radiative Transfer (RT) model, supported by a prior parameter distribution and a dynamic regularisation model and resulting in an inference of land surface parameters, such as effective Leaf Area Index (LAI), leaf chlorophyll concentration and fraction of senescent leaves, with full uncertainty quantification. The method is demonstrated over three agricultural FLUXNET sites, and the EO-LDAS results are compared with eight years of in situ measurements of FAPAR and LAI, resulting in a total of 24 site years. We additionally compare three other widely-used EO FAPAR products, namely the MEdium Resolution Imaging Spectrometer (MERIS) Full Resolution, the MISR High Resolution (HR) Joint Research Centre Two-stream Inversion Package (JRC-TIP) and MODIS MCD15 FAPAR products. The EO-LDAS MISR FAPAR retrievals show a high correlation with the ground measurements ( r 2 > 0.8), as well as the lowest average R M S E (0.14), in line with the MODIS product. As the EO-LDAS solution is effectively interpolated, if only measurements that are coincident with MISR observations are considered, the correlation increases ( r 2 > 0.85); the R M S E is lower by 4–5%; and the bias is 2% and 7%. The EO-LDAS MISR LAI estimates show a strong correlation with ground-based LAI (average r 2 = 0.76), but an underestimate of LAI for optically-thick canopies due to saturation (average R M S E = 2.23). These results suggest that the EO-LDAS approach is successful in retrieving both FAPAR and other land surface parameters. A large part of this success is based on the use of a dynamic regularisation model that counteracts the poor temporal sampling from the MISR instrument.

Satellite monitoring of disturbances in Arctic ecosystems

This study explores the capability of satellite remote sensing to detect relatively rapid changes of vegetation cover in northern Fennoscandian regions in response to disturbance more generally, and insect defoliation damage in particular. The data used is a long term time series of leaf area index (LAI) at 8 Km resolution derived from the Advanced Very High Resolution Radiometer (AVHRR) between 1982 and 2006, developed to be structurally consistent with the Moderate Resolution Imaging Spectrometer (MODIS) record. The study explores the potential of frequentist traditional statistics to detect disturbances at this coarse spatial resolution over the 25 years time series, and outlines the possibilities that Bayesian methods offer to improve the detection by including prior information on the profile of such disturbance events.

Modelling the impact of wildfire on spectral reflectance

This paper presents a method to extract information on the impact of wildfire on vegetation canopies. A simple linear model is proposed with a ‘generic’ spectral model of the impacts of wildfire on vegetation canopies. This allows a term related to the projected proportion of a pixel affected by fire (fcc) to be estimated from measurements of pre- and post-fire spectral reflectance. The properties of fcc are investigated using a hybrid radiative transfer model to simulate the impacts of wildfire in a multi-layered canopy. Spectral sampling from MODIS is assumed (7 bands). The fcc is confirmed to relate to the fractional area of a pixel affected by fire, although in multi-layer canopies it is modulated by a term dependent on the contrast between the pre-fire reflectance of areas affected by fire and those unaffected.

Synergies of Spaceborne Imaging Spectroscopy with Other Remote Sensing Approaches

Imaging spectroscopy (IS), also commonly known as hyperspectral remote sensing, is a powerful remote sensing technique for the monitoring of the Earth’s surface and atmosphere. Pixels in optical hyperspectral images consist of continuous reflectance spectra formed by hundreds of narrow spectral channels, allowing an accurate representation of the surface composition through spectroscopic techniques. However, technical constraints in the definition of imaging spectrometers make spectral coverage and resolution to be usually traded by spatial resolution and swath width, as opposed to optical multispectral (MS) systems typically designed to maximize spatial and/or temporal resolution. This complementarity suggests that a synergistic exploitation of spaceborne IS and MS data would be an optimal way to fulfill those remote sensing applications requiring not only high spatial and temporal resolution data, but also rich spectral information. On the other hand, IS has been shown to yield a strong synergistic potential with non-optical remote sensing methods, such as thermal infrared (TIR) and light detection and ranging (LiDAR). In this contribution we review theoretical and methodological aspects of potential synergies between optical IS and other remote sensing techniques. The focus is put on the evaluation of synergies between spaceborne optical IS and MS systems because of the expected availability of the two types of data in the next years. Short reviews of potential synergies of IS with TIR and LiDAR measurements are also provided.