David Frantz

An Operational Radiometric Landsat Preprocessing Framework for Large-Area Time Series Applications

We developed a large-area preprocessing framework for multisensor Landsat data, capable of processing large data volumes. Cloud and cloud shadow detection is performed by a modified Fmask code. Surface reflectance is inferred from Tanrés formulation of the radiative transfer, including adjacency effect correction. A precompiled MODIS water vapor database provides daily or climatological fallback estimates. Aerosol optical depth (AOD) is estimated over dark objects (DOs) that are identified in a combined database and image-based approach, where information on their temporal persistency is utilized. AOD is inferred with consideration of the actual target reflectance and background contamination effect. In case of absent DOs in bright scenes, a fallback approach with a modeled AOD climatology is used instead. Topographic normalization is performed by a modified C-correction. The data are projected into a single coordinate system and are organized in a gridded data structure for simplified pixel-based access. We based the assessment of the produced data set on an exhaustive analysis of overlapping pixels: 98.8% of the redundant overlaps are in the range of the expected ±2.5% overall radiometric algorithm accuracy. AOD is in very good agreement with Aerosol Robotic Network sunphotometer data (R2: 0.72 to 0.79, low intercepts, and slopes near unity). The uncertainty in using the water vapor fallback climatology is approximately ±2.8% for the TM SWIR1 band in the wet season. The topographic correction was considered successful by an investigation of the nonrelationship between the illumination angle and the corrected radiance.

Enhanced biomass prediction by assimilating satellite data into a crop growth model

Complex crop growth models (CGM) require a large number of input parameters, which can cause large errors if they are uncertain. Furthermore, they often lack spatial information. The coupling of a CGM with a radiative transfer model offers the possibility to assimilate remote sensing data while taking into account uncertainties in input parameters. A particle filter was used to assimilate satellite data into a CGM coupled with a leaf-canopy radiative transfer model to update biomass simulations of maize. The synthetic experiment set up to test the reliability of the procedure, highlighted the importance of the acquisition time. The real case study with RapidEye observations confirmed these findings. Data assimilation increased the accuracy of biomass predictions in the majority of the six maize fields where biomass validation data was available, with improvements of up to 15%. The smallest and largest errors in biomass prediction after assimilation were 82?kg/ha and 2116?kg/ha, respectively. Furthermore, data assimilation enabled the production of biomass maps showing detailed spatial variability. Data assimilation using a particle filter for biomass estimation was conducted.Proof of concept with synthetic case studies.Multispectral satellite data (visible and near infrared) was found to be suitable for data assimilation.Assimilation of satellite data allowed biomass prediction on a pixel basis.

Atmospheric Correction Inter-Comparison Exercise

The Atmospheric Correction Inter-comparison eXercise (ACIX) is an international initiative with the aim to analyse the Surface Reflectance (SR) products of various state-of-the-art atmospheric correction (AC) processors. The Aerosol Optical Thickness (AOT) and Water Vapour (WV) are also examined in ACIX as additional outputs of an AC processing. In this paper, the general ACIX framework is discussed; special mention is made of the motivation to initiate this challenge, the inter-comparison protocol and the principal results. ACIX is free and open and every developer was welcome to participate. Eventually, 12 participants applied their approaches to various Landsat-8 and Sentinel-2 image datasets acquired over sites around the world. The current results diverge depending on the sensors, products and sites, indicating their strengths and weaknesses. Indeed, this first implementation of processor inter-comparison was proven to be a good lesson for the developers to learn the advantages and limitations of their approaches. Various algorithm improvements are expected, if not already implemented, and the enhanced performances are yet to be investigated in future ACIX experiments.

Evaluating the trade-off between food and timber resulting from the conversion of Miombo forests to agricultural land in Angola using multi-temporal Landsat data.

The repopulation of abandoned areas in Angola after 27years of civil war led to a fast and extensive expansion of agricultural fields to meet the rising food demand. Yet, the increase in crop production at the expense of natural resources carries an inherent potential for conflicts since the demand for timber and wood extraction are also supposed to rise. We use the concept of ecosystem services to evaluate the trade-off between food and woody biomass. Our study area is located in central Angola, in the highlands of the upper Okavango catchment. We used Landsat data (spatial resolution: 30×30m) with a bi-temporal and multi-seasonal change detection approach for five time steps between 1989 and 2013 to estimate the conversion area from woodland to agriculture. Overall accuracy is 95%, users accuracy varies from 89-95% and producers accuracy ranges between 92-99%. To quantify the trade-off between woody biomass and the amount of food, this information was combined with indicator values and we furthermore assessed biomass regrowth on fallows. Our results reveal a constant rise in agricultural expansion from 1989-2013 with the mean annual deforestation rate increasing from roughly 5300ha up to about 12,000ha. Overall, 5.6% of the forested areas were converted to agriculture, whereas the FAO states a national deforestation rate for Angola of 5% from 1990-2010 (FAO, 2010). In the last time step 961,000t per year of woodland were cleared to potentially produce 1240t per year of maize. Current global agro-economical projections forecast increasing pressure on tropical dry forests from large-scale agriculture schemes (Gasparri et al., 2015; Searchinger and Heimlich, 2015). Our study underlines the importance of considering subsistence-related change processes, which may contribute significantly to negative effects associated with deforestation and degradation of these forest ecosystems.

Using Annual Landsat Time Series for the Detection of Dry Forest Degradation Processes in South-Central Angola

Dry tropical forests undergo massive conversion and degradation processes. This also holds true for the extensive Miombo forests that cover large parts of Southern Africa. While the largest proportional area can be found in Angola, the country still struggles with food shortages, insufficient medical and educational supplies, as well as the ongoing reconstruction of infrastructure after 27 years of civil war. Especially in rural areas, the local population is therefore still heavily dependent on the consumption of natural resources, as well as subsistence agriculture. This leads, on one hand, to large areas of Miombo forests being converted for cultivation purposes, but on the other hand, to degradation processes due to the selective use of forest resources. While forest conversion in south-central rural Angola has already been quantitatively described, information about forest degradation is not yet available. This is due to the history of conflicts and the therewith connected research difficulties, as well as the remote location of this area. We apply an annual time series approach using Landsat data in south-central Angola not only to assess the current degradation status of the Miombo forests, but also to derive past developments reaching back to times of armed conflicts. We use the Disturbance Index based on tasseled cap transformation to exclude external influences like inter-annual variation of rainfall. Based on this time series, linear regression is calculated for forest areas unaffected by conversion, but also for the pre-conversion period of those areas that were used for cultivation purposes during the observation time. Metrics derived from linear regression are used to classify the study area according to their dominant modification processes. We compare our results to MODIS latent integral trends and to further products to derive information on underlying drivers. Around 13% of the Miombo forests are affected by degradation processes, especially along streets, in villages, and close to existing agriculture. However, areas in presumably remote and dense forest areas are also affected to a significant extent. A comparison with MODIS derived fire ignition data shows that they are most likely affected by recurring fires and less by selective timber extraction. We confirm that areas that are used for agriculture are more heavily disturbed by selective use beforehand than those that remain unaffected by conversion. The results can be substantiated by the MODIS latent integral trends and we also show that due to extent and location, the assessment of forest conversion is most likely not sufficient to provide good estimates for the loss of natural resources.

Linking Land Surface Phenology and Vegetation-Plot Databases to Model Terrestrial Plant α-Diversity of the Okavango Basin

In many parts of Africa, spatially-explicit information on plant α-diversity, i.e., the number of species in a given area, is missing as baseline information for spatial planning. We present an approach on how to combine vegetation-plot databases and remotely-sensed land surface phenology (LSP) metrics to predict plant α-diversity on a regional scale. We gathered data on plant α-diversity, measured as species density, from 999 vegetation plots sized 20 m × 50 m covering all major vegetation units of the Okavango basin in the countries of Angola, Namibia and Botswana. As predictor variables, we used MODIS LSP metrics averaged over 12 years (250-m spatial resolution) and three topographic attributes calculated from the SRTM digital elevation model. Furthermore, we tested whether additional climatic data could improve predictions. We tested three predictor subsets: (1) remote sensing variables; (2) climatic variables; and (3) all variables combined. We used two statistical modeling approaches, random forests and boosted regression trees, to predict vascular plant α-diversity. The resulting maps showed that the Miombo woodlands of the Angolan Central Plateau featured the highest diversity, and the lowest values were predicted for the thornbush savanna in the Okavango Delta area. Models built on the entire dataset exhibited the best performance followed by climate-only models and remote sensing-only models. However, models including climate data showed artifacts. In spite of lower model performance, models based only on LSP metrics produced the most realistic maps. Furthermore, they revealed local differences in plant diversity of the landscape mosaic that were blurred by homogenous belts as predicted by climate-based models. This study pinpoints the high potential of LSP metrics used in conjunction with biodiversity data derived from vegetation-plot databases to produce spatial information on a regional scale that is urgently needed for basic natural resource management applications.

Improving the Spatial Resolution of Land Surface Phenology by Fusing Medium- and Coarse-Resolution Inputs

Satellite-derived land surface phenology (LSP) serves as a valuable input source for many environmental applications such as land cover classifications and global change studies. Commonly, LSP is derived from coarse-resolution (CR) sensors due to their well-suited temporal resolution. However, LSP is increasingly demanded at medium resolution (MR), but inferring LSP directly from MR imagery remains a challenging task (e.g., due to acquisition frequency). As such, we present a methodology that directly predicts MR LSP on the basis of the respective CR LSP and MR reflectance imagery. The approach considers information from the local pixel neighborhood at both resolutions by utilizing several prediction proxies, including spectral distance and multiscale heterogeneity metrics. The prediction performs well with simulated data

Forest Disturbance Mapping Using Dense Synthetic Landsat/MODIS Time-Series and Permutation-Based Disturbance Index Detection

(R^{2} = 0.84)

Fire spread from MODIS burned area data: obtaining fire dynamics information for every single fire

, and the approach substantially reduces noise. The size of the smallest reliably predicted object coincides with the effective CR pixel size (i.e., field-of-view). Nevertheless, even subpixel objects can be reliably predicted provided that pure CR pixels are located within the search radius. The application to real MODIS LSP and Landsat reflectance well preserves the phenological landscape composition, and the spatial refinement is especially striking in heterogeneous agricultural areas, where, for example, the circular shape of center pivot irrigation schemes is successfully restored at MR.

On the derivation of a spatially distributed aerosol climatology for its incorporation in a radiometric Landsat pre-processing framework

Spatio-temporal information on process-based forest loss is essential for a wide range of applications. Despite remote sensing being the only feasible means of monitoring forest change at regional or greater scales, there is no retrospectively available remote sensor that meets the demand of monitoring forests with the required spatial detail and guaranteed high temporal frequency. As an alternative, we employed the Spatial and Temporal Adaptive Reflectance Fusion Model (STARFM) to produce a dense synthetic time series by fusing Landsat and Moderate Resolution Imaging Spectroradiometer (MODIS) nadir Bidirectional Reflectance Distribution Function (BRDF) adjusted reflectance. Forest loss was detected by applying a multi-temporal disturbance detection approach implementing a Disturbance Index-based detection strategy. The detection thresholds were permutated with random numbers for the normal distribution in order to generate a multi-dimensional threshold confidence area. As a result, a more robust parameterization and a spatially more coherent detection could be achieved. (i) The original Landsat time series; (ii) synthetic time series; and a (iii) combined hybrid approach were used to identify the timing and extent of disturbances. The identified clearings in the Landsat detection were verified using an annual woodland clearing dataset from Queensland’s Statewide Landcover and Trees Study. Disturbances caused by stand-replacing events were successfully identified. The increased temporal resolution of the synthetic time series indicated promising additional information on disturbance timing. The results of the hybrid detection unified the benefits of both approaches, i.e., the spatial quality and general accuracy of the Landsat detection and the increased temporal information of synthetic time series. Results indicated that a temporal improvement in the detection of the disturbance date could be achieved relative to the irregularly spaced Landsat data for sufficiently large patches.