Predicting forest cover in distinct ecosystems: the potential of multi-source sentinel-1 and -2 data fusion

The fusion of microwave and optical data sets is expected to provide great potential for the derivation of forest cover around the globe. As Sentinel-1 and Sentinel-2 are now both operating in twin mode, they can provide an unprecedented data source to build dense spatial and temporal high-resolution time series across a variety of wavelengths. This study investigates (i) the ability of the individual sensors and (ii) their joint potential to delineate forest cover for study sites in two highly varied landscapes located in Germany (temperate dense mixed forests) and South Africa (open savanna woody vegetation and forest plantations). We used multi-temporal Sentinel-1 and single time steps of Sentinel-2 data in combination to derive accurate forest/non-forest (FNF) information via machine-learning classifiers. The forest classification accuracies were 90.9% and 93.2% for South Africa and Thuringia, respectively, estimated while using autocorrelation corrected spatial cross-validation (CV) for the fused data set. Sentinel-1 only classifications provided the lowest overall accuracy of 87.5%, while Sentinel-2 based classifications led to higher accuracies of 91.9%. Sentinel-2 short-wave infrared (SWIR) channels, biophysical parameters (Leaf Area Index (LAI), and Fraction of Absorbed Photosynthetically Active Radiation (FAPAR)) and the lower spectrum of the Sentinel-1 synthetic aperture radar (SAR) time series were found to be most distinctive in the detection of forest cover. In contrast to homogenous forests sites, Sentinel-1 time series information improved forest cover predictions in open savanna-like environments with heterogeneous regional features. The presented approach proved to be robust and it displayed the benefit of fusing optical and SAR data at high spatial resolution

Retrieval of biophysical parameters from multi-sensoral remote sensing data, assimilated into the crop growth model CERES-Wheat

This study investigated the possibilities and constraints for an integrated use of a crop growth model (CERES-Wheat) and earth observation techniques. The assimilation of information derived from earth observation sensors into crop growth models enables regional applications and may also help to improve the profound knowledge of the different involved processes and interactions. Both techniques can contribute to improved use of resources, reduced crop production risks, minimised environmental degradation, and increased farm income. Up to now, crop growth modelling and remote sensing techniquices mostly have been used separately for the assessment of agricultural applications. Crop growth models have made valuable contributions to, e.g., yield forecasting or to management decision support systems. Likewise, remote sensing techniques were successfully utilized in classification of agricultural areas or in the quantification of vegetation characteristics at various spatial and temporal scales. Multisensoral remote sensing approaches for the quantification biophysical variables are rarely realized. Normally the fusion of the data sources is based on the use of one sensor for classification purposes and the other one for the extraction of the desired parameters, based on the map classified previously. Pixel-based fusions between multispectral and SAR data is seldom realised for the assessment of quantitative parameters. The integration of crop growth models and remote sensing techniques by assimilating remotely sensed parameters into the models, is also still an issue of research. Especially, the integration of, e.g., multi-sensor biophysical parameter time-series for the improvement of the model performance, might feature a high potential. The starting point of the presented study was the question, if it is possible to derive the values of important crop variables from various remote sensing data? For the retrieval of these quantitative parameters by the use of various multispectral remote sensing sensors, intercalibration issues between the different retrieved vegetation indices had to be taken into account, in order to assure the comparability. Features influencing the vegetation indices are, e.g., the sensor geometry (like viewing- and solar-angle), atmospherical conditions, topography and spatial or radiometric resolution. However, the factors taken into account within this study are the spectral characteristics of the different sensors, like band position, bandwidth and centre wavelengths, which are described by the relative spectral response functions. Due to different RSR functions of the sensor bands, measured spectral differences occur, because the sensors record different components of the reflectance’s spectra from the monitored targets. These are then also introduced into the derived vegetation indices. The chosen cross-calibration method, intercalibrated the assessed Normalized Difference Vegetation Index and the Weighted Difference Vegetation Index between the various sensor pairs by regression, based on simulated multispectral sensors. Differences between the various assessed remote sensing sensors decreased form around 7% to below 1%. The intercalibration also had a positive impact on the later biophysical retrieval performance, producing sounder retrieval results. For the retrieval of the biophysical parameters empirical and semi-empirical models were assessed. The results indicate that the semi-empirical CLAIR model outperforms the empirical approaches. Not only for the Leaf Area Index retrieval, but also in the cases of all other assessed parameters. Concerning the other remote sensing data type used, the SAR data, it was analysed what potential different polarizations and incidence angles have for the extraction of the quantitative parameters. It became obvious that especially high incidence angles, as provided by the satellite Envisat ASAR, produce sounder retrieval results than lower incidence angles, due to a smaller amount of received soil signal. In the context of the assessed polarizations, sound results for the VV polarization could only be achieved for the retrieval of fresh biomass and the plant water content. For the ASAR sensor modelling fresh biomass and LAI using the HV polarization or the dry biomass using the ratio (HH/HV) was appropriate. As roughness aspects also have an influence on the retrieval performance from biophysical parameters using SAR data, the impact of soil surface and vegetation roughness was additionally considered. Best results were achieved, when also considering roughness features, however due to the need of regional modelling it is more appropriate not to consider them. For the calibration and re-tuning of crop growth models information about important phenological events such as heading/flowering is rather important. After this stage reproductive growth begins, whereby the number of kernels per plant is often calculated from plant weight at flowering and kernel weight is calculated from time and temperature available for dry matter distribution. By the use of the SAR VV time-series this important stage could be successfully extracted. Further methods for pixel-based fused biophysical parameter estimations, using SAR and multispectral data were analysed. By this approach the different features, being monitored of the two systems, are combined for sounder parameter retrieval. The assessed method of combining the multi-sensoral information by linear regression did not bring sound results and was outperformed by single sensor use, only taking into account the multispectral information. Only for the parameter fresh biomass, modelling based on the NDIV and the ASAR ratio slightly outperformed the single sensor modelling approaches. The complex combined modelling by the use of the CLAIR and the Water Cloud Model featured no valid results. For the combination, by using the CLAIR model and multiple regression slight improvements, in contrast to the single multispectral sensor use, were achieved. Especially, during late phenological stages, the assessed VV information improved the modelling results, in comparison to only using the CLAIR model. All the findings could finally be successfully applied for regional estimations. Only the roughness features could not be applied, due to the fact, that it is hard to regionally assess this needed model input parameter. Regional parameter on the basis of remote sensing data, is the major advantage of this technique, due to the large spatial overview given. The second main question was, if it is possible to integrate the crop variables gained from multisensoral data into a crop growth model, increasing the final yield estimation accuracy. Thus far, beneficial linkages between both techniques have been often limited to land use classification via remote sensing for choosing the adequate model and quantification of crop growth and development curves using biophysical parameters derived from remote sensing images for model calibration. Only a few studies actually considered the potentials of remote sensing for model re-initialization of growth and development characteristics of a specific crop, as the here studied winter wheat. Overall, the integration of remotely sensed variables into the crop growth model CERES-Wheat led to an improved final yield estimation accuracy in comparison to an automatic input parameter setting. The assessed final yield bias for the automatic input parameter setting summed up to 6.6%. When re-initializing the most sensitive input parameters (sowing date and fertilizer application date) by the use of remotely sensed biophysical variables the biases ranged from 0.56% overestimation to 5.4% understimation, in dependence of the data series used for assimilation. Whereby, it was assessed that the combined dense data series, considering SAR and multispectral information, slightly outperformed the performance of the full multispectral data series. However, when analysing the assimilation of the multispectral data series in further detail, it became clear that the actually information from the phenological stage ripening declines the modelling performance and thus the final yield estimation accuracy. When neglecting the information from this phenological stage the reduced multispectral data series performed as sound as the dense data series containing SAR and multispectral information. Thus, when the appropriate phenological stages are monitored by multispectral data, additional SAR information does not lead to a model improvement. However, when important dates are not monitored by multispectral images, e.g., due to cloud coverage, the additionally considered SAR information was not able to appropriatly fill these important multispectral time gaps. They even had a more negeative influence on the modelling performance. Overall, the best results could be obtained by assimilating a multispectral data series, covering the crop development during the important phenological stages stem elongation and flowering (without ripening stage), into the CERES-Wheat model. Finally, the integration of remote sensing data in the point-based crop growth model allowed it‘s spatial application for prediction of wheat production at a more regional scale. This approach also outperformed another evaluated method of direct multi-sensoral regional yield estimation. This study has demonstrated that biophysical parameters can be retrieved from remote sensing data and led, when assimilated into a crop growth model, to an improved final yield estimation. However, overall the SAR information did not really have a significant positive effect on the multi-sensoral biophysical parameter retrieval and on the later assimilation process. Thus, overall SAR information should only be considered, when multispectral data acquisitions are tremendously hampered by cloud coverage. The assessed assimilation of remote sensing information into a crop growth model had a positive effect on the final yield estimation performance. The analysed method, combining remote sensing and crop growth model techniques, was succsessfully demonstrated and will gain even more importance in the future for, e.g., decision support systems fine-tuning fertilizer regimes and thus contributing to more environmentally sound and sustained agricultural production

Assessment of multi-temporal, multi-sensor radar and ancillary spatial data for grasslands monitoring in Ireland using machine learning approaches

Accurate inventories of grasslands are important for studies of carbon dynamics, biodiversity conservation and agricultural management. For regions with persistent cloud cover the use of multi-temporal synthetic aperture radar (SAR) data provides an attractive solution for generating up-to-date inventories of grasslands. This is even more appealing considering the data that will be available from upcoming missions such as Sentinel-1 and ALOS-2. In this study, the performance of three machine learning algorithms; Random Forests (RF), Support Vector Machines (SVM) and the relatively underused Extremely Randomised Trees (ERT) is evaluated for discriminating between grassland types over two large heterogeneous areas of Ireland using multi-temporal, multi-sensor radar and ancillary spatial datasets. A detailed accuracy assessment shows the efficacy of the three algorithms to classify different types of grasslands. Overall accuracies ≥ 88.7% (with kappa coefficient of 0.87) were achieved for the single frequency classifications and maximum accuracies of 97.9% (kappa coefficient of 0.98) for the combined frequency classifications. For most datasets, the ERT classifier outperforms SVM and RF

Integration of L-band SAR data into Land Surface Models

Abstract—Land surface process modelling might be limited due to lack of reliable model input data. Key surface variables as land cover information or soil moisture conditions have been proven to be observable by remote sensing systems. The integration of remote sensing data into land surface process models might therefore help to improve their simulations results. Longer wavelength SAR data has a higher sensitivity to soil moisture content than higher frequency systems. Recent (ALOS) and planed (e.g. TerraSAR-L) SAR systems are therefore expected to provide valuable information about soil moisture dynamics. The present study investigates the potential to retrieve land cover information and geophysical parameters from L-band SAR data. The retrieval results are assimilated into a state-of-the-art land surface model to evaluate the merit of L-band SAR data assimilation

Coupled modelling of land surface microwave interactions using ENVISAT ASAR data

In the last decades microwave remote sensing has proven its capability to provide valuable information about the land surface. New sensor generations as e.g. ENVISAT ASAR are capable to provide frequent imagery with an high information content. To make use of these multiple imaging capabilities, sophisticated parameter inversion and assimilation strategies have to be applied. A profound understanding of the microwave interactions at the land surface is therefore essential. The objective of the presented work is the analysis and quantitative description of the backscattering processes of vegetated areas by means of microwave backscattering models. The effect of changing imaging geometries is investigated and models for the description of bare soil and vegetation backscattering are developed. Spatially distributed model parameterisation is realized by synergistic coupling of the microwave scattering models with a physically based land surface process model. This enables the simulation of realistic SAR images, based on bioand geophysical parameters. The adequate preprocessing of the datasets is crucial for quantitative image analysis. A stringent preprocessing and sophisticated terrain geocoding and correction procedure is therefore suggested. It corrects the geometric and radiometric distortions of the image products and is taken as the basis for further analysis steps. A problem in recently available microwave backscattering models is the inadequate parameterisation of the surface roughness. It is shown, that the use of classical roughness descriptors, as the rms height and autocorrelation length, will lead to ambiguous model parameterisations. A new two parameter bare soil backscattering model is therefore recommended to overcome this drawback. It is derived from theoretical electromagnetic model simulations. The new bare soil surface scattering model allows for the accurate description of the bare soil backscattering coefficients. A new surface roughness parameter is introduced in this context, capable to describe the surface roughness components, affecting the backscattering coefficient. It is shown, that this parameter can be directly related to the intrinsic fractal properties of the surface. Spatially distributed information about the surface roughness is needed to derive land surface parameters from SAR imagery. An algorithm for the derivation of the new surface roughness parameter is therefore suggested. It is shown, that it can be derived directly from multitemporal SAR imagery. Starting from that point, the bare soil backscattering model is used to assess the vegetation influence on the signal. By comparison of the residuals between measured backscattering coefficients and those predicted by the bare soil backscattering model, the vegetation influence on the signal can be quantified. Significant difference between cereals (wheat and triticale) and maize is observed in this context. It is shown, that the vegetation influence on the signal can be directly derived from alternating polarisation data for cereal fields. It is dependant on plant biophysical variables as vegetation biomass and water content. The backscattering behaviour of a maize stand is significantly different from that of other cereals, due to its completely different density and shape of the plants. A dihedral corner reflection between the soil and the stalk is identified as the major source of backscattering from the vegetation. A semiempirical maize backscattering model is suggested to quantify the influences of the canopy over the vegetation period. Thus, the different scattering contributions of the soil and vegetation components are successfully separated. The combination of the bare soil and vegetation backscattering models allows for the accurate prediction of the backscattering coefficient for a wide range of surface conditions and variable incidence angles. To enable the spatially distributed simulation of the SAR backscattering coefficient, an interface to a process oriented land surface model is established, which provides the necessary input variables for the backscattering model. Using this synergistic, coupled modelling approach, a realistic simulation of SAR images becomes possible based on land surface model output variables. It is shown, that this coupled modelling approach leads to promising and accurate estimates of the backscattering coefficients. The remaining residuals between simulated and measured backscatter values are analysed to identify the sources of uncertainty in the model. A detailed field based analysis of the simulation results revealed that imprecise soil moisture predictions by the land surface model are a major source of uncertainty, which can be related to imprecise soil texture distribution and soil hydrological properties. The sensitivity of the backscattering coefficient to the soil moisture content of the upper soil layer can be used to generate soil moisture maps from SAR imagery. An algorithm for the inversion of soil moisture from the upper soil layer is suggested and validated. It makes use of initial soil moisture values, provided by the land surface process model. Soil moisture values are inverted by means of the coupled land surface backscattering model. The retrieved soil moisture results have an RMSE of 3.5 Vol %, which is comparable to the measurement accuracy of the reference field data. The developed models allow for the accurate prediction of the SAR backscattering coefficient. The various soil and vegetation scattering contributions can be separated. The direct interface to a physically based land surface process model allows for the spatially distributed modelling of the backscattering coefficient and the direct assimilation of remote sensing data into a land surface process model. The developed models allow for the derivation of static and dynamic landsurface parameters, as e.g. surface roughness, soil texture, soil moisture and biomass from remote sensing data and their assimilation in process models. They are therefore reliable tools, which can be used for sophisticated practice oriented problem solutions in manifold manner in the earth and environmental sciences

Modelling Net Primary Productivity and Above-Ground Biomass for Mapping of Spatial Biomass Distribution in Kazakhstan

Biomass is an important ecological variable for understanding the responses of vegetation to the currently observed global change. The impact of changes in vegetation biomass on the global ecosystem is also of high relevance. The vegetation in the arid and semi-arid environments of Kazakhstan is expected to be affected particularly strongly by future climate change. Therefore, it is of great interest to observe large-scale vegetation dynamics and biomass distribution in Kazakhstan. At the beginning of this dissertation, previous research activities and remote-sensing-based methods for biomass estimation in semi-arid regions have been comprehensively reviewed for the first time. The review revealed that the biggest challenge is the transferability of methods in time and space. Empirical approaches, which are predominantly applied, proved to be hardly transferable. Remote-sensing-based Net Primary Productivity (NPP) models, on the other hand, allow for regional to continental modelling of NPP time-series and are potentially transferable to new regions. This thesis thus deals with modelling and analysis of NPP time-series for Kazakhstan and presents a methodological concept for derivation of above-ground biomass estimates based on NPP data. For validation of the results, biomass field data were collected in three study areas in Kazakhstan. For the selection of an appropriate model, two remote-sensing-based NPP models were applied to a study area in Central Kazakhstan. The first is the Regional Biomass Model (RBM). The second is the Biosphere Energy Transfer Hydrology Model (BETHY/DLR). Both models were applied to Kazakhstan for the first time in this dissertation. Differences in the modelling approaches, intermediate products, and calculated NPP, as well as their temporal characteristics were analysed and discussed. The model BETHY/DLR was then used to calculate NPP for Kazakhstan for 2003–2011. The results were analysed regarding spatial, intra-annual, and inter-annual variations. In addition, the correlation between NPP and meteorological parameters was analysed. In the last part of this dissertation, a methodological concept for derivation of above-ground biomass estimates of natural vegetation from NPP time-series has been developed. The concept is based on the NPP time-series, information about fractional cover of herbaceous and woody vegetation, and plants’ relative growth rates (RGRs). It has been the first time that these parameters are combined for biomass estimation in semi-arid regions. The developed approach was finally applied to estimate biomass for the three study areas in Kazakhstan and validated with field data. The results of this dissertation provide information about the vegetation dynamics in Kazakhstan for 2003–2011. This is valuable information for a sustainable land management and the identification of regions that are potentially affected by a changing climate. Furthermore, a methodological concept for the estimation of biomass based on NPP time-series is presented. The developed method is potentially transferable. Providing that the required information regarding vegetation distribution and fractional cover is available, the method will allow for repeated and large-area biomass estimation for natural vegetation in Kazakhstan and other semi-arid environments

Application of RADARSAT-2 Polarimetric Data for Land Use and Land Cover Classification and Crop monitoring in Southwestern Ontario

Timely and accurate information of land surfaces is desirable for land change detection and crop condition monitoring. Optical data have been widely used in Land Use and Land Cover (LU/LC) mapping and crop condition monitoring. However, due to unfavorable weather conditions, high quality optical images are not always available. Synthetic Aperture Radar (SAR) sensors, such as RADARSAT-2, are able to transmit microwaves through cloud cover and light rain, and thus offer an alternative data source. This study investigates the potential of multi-temporal polarimetric RADARSAT-2 data for LU/LC classification and crop monitoring in the urban rural fringe areas of London, Ontario. Nine LU/LC classes were identified with a high overall accuracy of 91.0%. Also, high correlations have been found within the corn and soybean fields between some polarimetric parameters and Normalized Difference Vegetation Index (NDVI). The results demonstrate the capability of RADARSAT-2 in LU/LC classification and crop condition monitoring

Evaluation of Sentinel-2 Red-Edge Bands for Empirical Estimation of Green LAI and Chlorophyll Content

ESA’s upcoming satellite Sentinel-2 will provide Earth images of high spatial, spectral and temporal resolution and aims to ensure continuity for Landsat and SPOT observations. In comparison to the latter sensors, Sentinel-2 incorporates three new spectral bands in the red-edge region, which are centered at 705, 740 and 783 nm. This study addresses the importance of these new bands for the retrieval and monitoring of two important biophysical parameters: green leaf area index (LAI) and chlorophyll content (Ch). With data from several ESA field campaigns over agricultural sites (SPARC, AgriSAR, CEFLES2) we have evaluated the efficacy of two empirical methods that specifically make use of the new Sentinel-2 bands. First, it was shown that LAI can be derived from a generic normalized difference index (NDI) using hyperspectral data, with 674 nm with 712 nm as best performing bands. These bands are positioned closely to the Sentinel-2 B4 (665 nm) and the new red-edge B5 (705 nm) band. The method has been applied to simulated Sentinel-2 data. The resulting green LAI map was validated against field data of various crop types, thereby spanning a LAI between 0 and 6, and yielded a RMSE of 0.6. Second, the recently developed “Normalized Area Over reflectance Curve” (NAOC), an index that derives Ch from hyperspectral data, was studied on its compatibility with simulated Sentinel-2 data. This index integrates the reflectance curve between 643 and 795 nm, thereby including the new Sentinel-2 bands in the red-edge region. We found that these new bands significantly improve the accuracy of Ch estimation. Both methods emphasize the importance of red-edge bands for operational estimation of biophysical parameters from Sentinel-2