Harm Bartholomeus

The response of Arctic vegetation to the summer climate: relation between shrub cover, NDVI, surface albedo and temperature

Recently observed Arctic greening trends from normalized difference vegetation index (NDVI) data suggest that shrub growth is increasing in response to increasing summer temperature. An increase in shrub cover is expected to decrease summer albedo and thus positively feed back to climate warming. However, it is unknown how albedo and NDVI are affected by shrub cover and inter-annual variations in the summer climate. Here, we examine the relationship between deciduous shrub fractional cover, NDVI and albedo using field data collected at a tundra site in NE Siberia. Field data showed that NDVI increased and albedo decreased with increasing deciduous shrub cover. We then selected four Arctic tundra study areas and compiled annual growing season maximum NDVI and minimum albedo maps from MODIS satellite data (2000–10) and related these satellite products to tundra vegetation types (shrub, graminoid, barren and wetland tundra) and regional summer temperature. We observed that maximum NDVI was greatest in shrub tundra and that inter-annual variation was negatively related to summer minimum albedo but showed no consistent relationship with summer temperature. Shrub tundra showed higher albedo than wetland and barren tundra in all four study areas. These results suggest that a northwards shift of shrub tundra might not lead to a decrease in summer minimum albedo during the snow-free season when replacing wetland tundra. A fully integrative study is however needed to link results from satellite data with in situ observations across the Arctic to test the effect of increasing shrub cover on summer albedo in different tundra vegetation types.

Soil Organic Carbon mapping of partially vegetated agricultural fields with imaging spectroscopy

Soil Organic Carbon (SOC) is one of the key soil properties, but the large spatial variation makes continuous mapping a complex task. Imaging spectroscopy has proven to be an useful technique for mapping of soil properties, but the applicability decreases rapidly when fields are partially covered with vegetation. In this paper we show that with only a few percent fractional maize cover the accuracy of a Partial Least Square Regression (PLSR) based SOC prediction model drops dramatically. However, this problem can be solved with the use of spectral unmixing techniques. First, the fractional maize cover is determined with linear spectral unmixing, taking the illumination and observation angles into account. In a next step the influence of maize is filtered out from the spectral signal by a new procedure termed Residual Spectral Unmixing (RSU). The residual soil spectra resulting from this procedure are used for mapping of SOC using PLSR, which could be done with accuracies comparable to studies performed on bare soil surfaces (Root Mean Standard Error of Calibration = 1.34 g/kg and Root Mean Standard Error of Prediction = 1.65 g/kg). With the presented RSU approach it is possible to filter out the influence of maize from the mixed spectra, and the residual soil spectra contain enough information for mapping of the SOC distribution within agricultural fields. This can improve the applicability of airborne imaging spectroscopy for soil studies in temperate climates, since the use of the RSU approach can extend the flight-window which is often constrained by the presence of vegetation.

A Lightweight Hyperspectral Mapping System and Photogrammetric Processing Chain for Unmanned Aerial Vehicles

During the last years commercial hyperspectral imaging sensors have been miniaturized and their performance has been demonstrated on Unmanned Aerial Vehicles (UAV). However currently the commercial hyperspectral systems still require minimum payload capacity of approximately 3 kg, forcing usage of rather large UAVs. In this article we present a lightweight hyperspectral mapping system (HYMSY) for rotor-based UAVs, the novel processing chain for the system, and its potential for agricultural mapping and monitoring applications. The HYMSY consists of a custom-made pushbroom spectrometer (400–950 nm, 9 nm FWHM, 25 lines/s, 328 px/line), a photogrammetric camera, and a miniature GPS-Inertial Navigation System. The weight of HYMSY in ready-to-fly configuration is only 2.0 kg and it has been constructed mostly from off-the-shelf components. The processing chain uses a photogrammetric algorithm to produce a Digital Surface Model (DSM) and provides high accuracy orientation of the system over the DSM. The pushbroom data is georectified by projecting it onto the DSM with the support of photogrammetric orientations and the GPS-INS data. Since an up-to-date DSM is produced internally, no external data are required and the processing chain is capable to georectify pushbroom data fully automatically. The system has been adopted for several experimental flights related to agricultural and habitat monitoring applications. For a typical flight, an area of 2–10 ha was mapped, producing a RGB orthomosaic at 1–5 cm resolution, a DSM at 5–10 cm resolution, and a hyperspectral datacube at 10–50 cm resolution.

Determining iron content in Mediterranean soils in partly vegetated areas, using spectral reflectance and imaging spectroscopy

The possibility of quantifying iron content in the topsoil of the slopes of the El Hacho Mountain complex in Southern Spain using imaging spectroscopy is investigated. Laboratory, field and airborne spectrometer (ROSIS) data are acquired, in combination with soil samples, which are analysed for dithionite extractable iron (Fed) content. Analysis of the properties of two iron related absorption features present in laboratory spectra demonstrates good relations, especially between the standard deviation (S.D.) of the values in an absorption feature and the Fed content (R2 = 0.67) as well as the ratio based Redness Index (R2 = 0.51). Such derived relations are less strong for the ROSIS data (R2 for S.D. = 0.26 and R2 for Redness Index = 0.22). The spatial distribution of iron in vegetated areas shows a strong sensitivity of these relations with the presence of vegetation. A combination of both methods shows that the overestimation of the Fed content with the one method is (partly) compensated by the underestimation with the other method.

Scaling dimensions in spectroscopy of soil and vegetation

Abstract The paper revises and clarifies definitions of the term scale and scaling conversions for imaging spectroscopy of soil and vegetation. We demonstrate a new four-dimensional scale concept that includes not only spatial but also the spectral, directional and temporal components. Three scaling remote sensing techniques are reviewed: (1) radiative transfer, (2) spectral (un)mixing, and (3) data fusion. Relevant case studies are given in the context of their up- and/or down-scaling abilities over the soil/vegetation surfaces and a multi-source approach is proposed for their integration. Radiative transfer (RT) models are described to show their capacity for spatial, spectral up-scaling, and directional down-scaling within a heterogeneous environment. Spectral information and spectral derivatives, like vegetation indices (e.g. TCARI/OSAVI), can be scaled and even tested by their means. Radiative transfer of an experimental Norway spruce (Picea abies (L.) Karst.) research plot in the Czech Republic was simulated by the Discrete Anisotropic Radiative Transfer (DART) model to prove relevance of the correct object optical properties scaled up to image data at two different spatial resolutions. Interconnection of the successive modelling levels in vegetation is shown. A future development in measurement and simulation of the leaf directional spectral properties is discussed. We describe linear and/or non-linear spectral mixing techniques and unmixing methods that demonstrate spatial down-scaling. Relevance of proper selection or acquisition of the spectral endmembers using spectral libraries, field measurements, and pure pixels of the hyperspectral image is highlighted. An extensive list of advanced unmixing techniques, a particular example of unmixing a reflective optics system imaging spectrometer (ROSIS) image from Spain, and examples of other mixture applications give insight into the present status of scaling capabilities. Simultaneous spatial and temporal down-scaling by means of a data fusion technique is described. A demonstrative example is given for the moderate resolution imaging spectroradiometer (MODIS) and LANDSAT Thematic Mapper (TM) data from Brazil. Corresponding spectral bands of both sensors were fused via a pyramidal wavelet transform in Fourier space. New spectral and temporal information of the resultant image can be used for thematic classification or qualitative mapping. All three described scaling techniques can be integrated as the relevant methodological steps within a complex multi-source approach. We present this concept of combining numerous optical remote sensing data and methods to generate inputs for ecosystem process models.

Estimating salinity stress in sugarcane fields with spaceborne hyperspectral vegetation indices

The presence of salt in the soil profile negatively affects the growth and development of vegetation. As a result, the spectral reflectance of vegetation canopies varies for different salinity levels. This research was conducted to (1) investigate the capability of satellite-based hyperspectral vegetation indices (VIs) for estimating soil salinity in agricultural fields, (2) evaluate the performance of 21 existing VIs and (3) develop new VIs based on a combination of wavelengths sensitive for multiple stresses and find the best one for estimating soil salinity. For this purpose a Hyperion image of September 2, 2010, and data on soil salinity at 108 locations in sugarcane (Saccharum officina L.) fields were used. Results show that soil salinity could well be estimated by some of these VIs. Indices related to chlorophyll absorption bands or based on a combination of chlorophyll and water absorption bands had the highest correlation with soil salinity. In contrast, indices that are only based on water absorption bands had low to medium correlations, while indices that use only visible bands did not perform well. From the investigated indices the optimized soil-adjusted vegetation index (OSAVI) had the strongest relationship (R2 = 0.69) with soil salinity for the training data, but it did not perform well in the validation phase. The validation procedure showed that the new salinity and water stress indices (SWSI) implemented in this study (SWSI-1, SWSI-2, SWSI-3) and the Vogelmann red edge index yielded the best results for estimating soil salinity for independent fields with root mean square errors of 1.14, 1.15, 1.17 and 1.15 dS/m, respectively. Our results show that soil salinity could be estimated by satellite-based hyperspectral VIs, but validation of obtained models for independent data is essential for selecting the best model.

Using MERIS on Envisat for land cover mapping in the Netherlands

This paper describes the results of a feasibility study to test the usefulness of MERIS for land cover mapping. The Netherlands was used as a test site because of its highly fragmented landscape. Results showed that the geometric and radiometric properties of the studied MERIS images of the Netherlands are suitable for land applications. Calculation of principal components and correlation coefficients revealed that the 15 MERIS bands provided a lot of redundant spectral information. For land applications, information came from the visible part of the spectrum on the one hand and from the near‐infrared part on the other hand. In addition, the red‐edge slope of the reflectance curve (in particular MERIS band 9 at about 708 nm) provided supplementary information. The Dutch land use database LGN5 was used as a reference for classifications in this study after aggregation from 25 m to 300 m and recoding to 7 relevant land cover classes. For land cover classification best results in terms of classification accuracies were obtained for the image of 14 July 2003. For the seven land cover classes selected the overall classification accuracy was 67.2%. A multitemporal classification did not improve the overall classification accuracy.

Measuring and modeling the effect of surface moisture on the spectral reflectance of coastal beach sand

Surface moisture is an important supply limiting factor for aeolian sand transport, which is the primary driver of coastal dune development. As such, it is critical to account for the control of surface moisture on available sand for dune building. Optical remote sensing has the potential to measure surface moisture at a high spatio-temporal resolution. It is based on the principle that wet sand appears darker than dry sand: it is less reflective. The goals of this study are (1) to measure and model reflectance under controlled laboratory conditions as function of wavelength () and surface moisture () over the optical domain of 350–2500 nm, and (2) to explore the implications of our laboratory findings for accurately mapping the distribution of surface moisture under natural conditions. A laboratory spectroscopy experiment was conducted to measure spectral reflectance (1 nm interval) under different surface moisture conditions using beach sand. A non-linear increase of reflectance upon drying was observed over the full range of wavelengths. Two models were developed and tested. The first model is grounded in optics and describes the proportional contribution of scattering and absorption of light by pore water in an unsaturated sand matrix. The second model is grounded in soil physics and links the hydraulic behaviour of pore water in an unsaturated sand matrix to its optical properties. The optical model performed well for volumetric moisture content 24% ( 0.97), but underestimated reflectance for between 24–30% ( 0.92), most notable around the 1940 nm water absorption peak. The soil-physical model performed very well ( 0.99) but is limited to 4% 24%. Results from a field experiment show that a short-wave infrared terrestrial laser scanner ( = 1550 nm) can accurately relate surface moisture to reflectance (standard error 2.6%), demonstrating its potential to derive spatially extensive surface moisture maps of a natural coastal beach.

Relationships between declining summer sea ice, increasing temperatures and changing vegetation in the Siberian Arctic tundra from MODIS time series (2000–11)

The concern about Arctic greening has grown recently as the phenomenon is thought to have significant influence on global climate via atmospheric carbon emissions. Earlier work on Arctic vegetation highlighted the role of summer sea ice decline in the enhanced warming and greening phenomena observed in the region, but did not contain enough details for spatially characterizing the interactions between sea ice, temperature and vegetation photosynthetic absorption. By using 1 km resolution data from the Moderate Resolution Imaging Spectrometer (MODIS) as a primary data source, this study presents detailed maps of vegetation and temperature trends for the Siberian Arctic region, using the time integrated normalized difference vegetation index (TI-NDVI) and summer warmth index (SWI) calculated for the period 2000‐11 to represent vegetation greenness and temperature respectively. Spatio-temporal relationships between the two indices and summer sea ice conditions were investigated with transects at eight locations using sea ice concentration data from the Special Sensor Microwave/Imager (SSM/I). In addition, the derived vegetation and temperature trends were compared among major Arctic vegetation types and bioclimate subzones. The fine resolution trend map produced confirms the overall greening (C1% yr 1 ) and warming (C0.27% yr 1 ) of the region, reported in previous studies, but also reveals browning areas. The causes of such local decreases in vegetation, while surrounding areas are experiencing the opposite reaction to changing conditions, are still unclear. Overall correlations between sea ice concentration and SWI as well as TI-NDVI decreased in strength with increasing distance from the coast, with a particularly pronounced pattern in the case of SWI. SWI appears to be driving TI-NDVI in many cases, but not systematically, highlighting the presence of limiting factors other than temperature for plant growth in the region. Further unravelling those limiting factors constitutes a priority in future research. This study demonstrates the use of medium resolution remotely sensed data for studying the complexity of spatio-temporal vegetation dynamics in the Arctic.

Prediction of soil properties using imaging spectroscopy: Considering fractional vegetation cover to improve accuracy

Spectroscopic techniques have become attractive to assess soil properties because they are fast, require little labor and may reduce the amount of laboratory waste produced when compared to conventional methods. Imaging spectroscopy (IS) can have further advantages compared to laboratory or field proximal spectroscopic approaches such as providing spatially continuous information with a high density. However, the accuracy of IS derived predictions decreases when the spectral mixture of soil with other targets occurs. This paper evaluates the use of spectral data obtained by an airborne hyperspectral sensor (ProSpecTIR-VS – Aisa dual sensor) for prediction of physical and chemical properties of Brazilian highly weathered soils (i.e., Oxisols). A methodology to assess the soil spectral mixture is adapted and a progressive spectral dataset selection procedure, based on bare soil fractional cover, is proposed and tested. Satisfactory performances are obtained specially for the quantification of clay, sand and CEC using airborne sensor data (R2 of 0.77, 0.79 and 0.54; RPD of 2.14, 2.22 and 1.50, respectively), after spectral data selection is performed; although results obtained for laboratory data are more accurate (R2 of 0.92, 0.85 and 0.75; RPD of 3.52, 2.62 and 2.04, for clay, sand and CEC, respectively). Most importantly, predictions based on airborne-derived spectra for which the bare soil fractional cover is not taken into account show considerable lower accuracy, for example for clay, sand and CEC (RPD of 1.52, 1.64 and 1.16, respectively). Therefore, hyperspectral remotely sensed data can be used to predict topsoil properties of highly weathered soils, although spectral mixture of bare soil with vegetation must be considered in order to achieve an improved prediction accuracy.