Forest classification and impact of BIOMASS resolution on forest area and aboveground biomass estimation

The European Space Agency (ESA) is currently implementing the BIOMASS mission as 7th Earth Explorer satellite. BIOMASS will provide for the first time global forest aboveground biomass estimates based on P-band synthetic aperture radar (SAR) imagery. This paper addresses an often overlooked element of the data processing chain required to ensure reliable and accurate forest biomass estimates: accurate identification of forest areas ahead of the inversion of radar data into forest biomass estimates.The use of the P-band data from BIOMASS itself for the classification into forest and non-forest land cover types is assessed in this paper. For airborne data in tropical, hemi-boreal and boreal forests we demonstrate that classification accuracies from 90 up to 97% can be achieved using radar backscatter and phase information. However, spaceborne data will have a lower resolution and higher noise level compared to airborne data and a higher probability of mixed pixels containing multiple land cover types. Therefore, airborne data was reduced to 50 m, 100 m and 200 m resolution. The analysis revealed that about 50–60% of the area within the resolution level must be covered by forest to classify a pixel with higher probability as forest compared to non-forest. This results in forest omission and commission leading to similar forest area estimation over all resolutions. However, the forest omission resulted in a biased underestimated biomass, which was not equaled by the forest commission. The results underline the necessity of a highly accurate pre-classification of SAR data for an accurate unbiased aboveground biomass estimation.     

Estimating temporal changes in soil carbon stocks at ecoregional scale in Madagascar using remote-sensing

Soil organic carbon (SOC) plays an important role in climate change regulation notably through release of CO2 following land use change such a deforestation, but data on stock change levels are lacking. This study aims to empirically assess SOC stocks change between 1991 and 2011 at the landscape scale using easy-to-access spatially-explicit environmental factors. The study area was located in southeast Madagascar, in a region that exhibits very high rate of deforestation and which is characterized by both humid and dry climates. We estimated SOC stock on 0.1 ha plots for 95 different locations in a 43,000 ha reference area covering both dry and humid conditions and representing different land cover including natural forest, cropland, pasture and fallows. We used the Random Forest algorithm to find out the environmental factors explaining the spatial distribution of SOC. We then predicted SOC stocks for two soil layers at 30 cm and 100 cm over a wider area of 395,000 ha. By changing the soil and vegetation indices derived from remote sensing images we were able to produce SOC maps for 1991 and 2011. Those estimates and their related uncertainties where combined in a post-processing step to map estimates of significant SOC variations and we finally compared the SOC change map with published deforestation maps. Results show that the geologic variables, precipitation, temperature, and soil-vegetation status were strong predictors of SOC distribution at regional scale. We estimated an average net loss of 10.7% and 5.2% for the 30 cm and the 100 cm layers respectively for deforested areas in the humid area. Our results also suggest that these losses occur within the first five years following deforestation. No significant variations were observed for the dry region. This study provides new solutions and knowledge for a better integration of soil threats and opportunities in land management policies.     

Above ground biomass and tree species richness estimation with airborne lidar in tropical Ghana forests

Estimates of forest aboveground biomass are fundamental for carbon monitoring and accounting; delivering information at very high spatial resolution is especially valuable for local management, conservation and selective logging purposes. In tropical areas, hosting large biomass and biodiversity resources which are often threatened by unsustainable anthropogenic pressures, frequent forest resources monitoring is needed. Lidar is a powerful tool to estimate aboveground biomass at fine resolution; however its application in tropical forests has been limited, with high variability in the accuracy of results. Lidar pulses scan the forest vertical profile, and can provide structure information which is also linked to biodiversity. In the last decade the remote sensing of biodiversity has received great attention, but few studies focused on the use of lidar for assessing tree species richness in tropical forests.This research aims at estimating aboveground biomass and tree species richness using discrete return airborne lidar in Ghana forests. We tested an advanced statistical technique, Multivariate Adaptive Regression Splines (MARS), which does not require assumptions on data distribution or on the relationships between variables, being suitable for studying ecological variables.We compared the MARS regression results with those obtained by multilinear regression and found that both algorithms were effective, but MARS provided higher accuracy either for biomass (R² = 0.72) and species richness (R² = 0.64). We also noted strong correlation between biodiversity and biomass field values. Even if the forest areas under analysis are limited in extent and represent peculiar ecosystems, the preliminary indications produced by our study suggest that instrument such as lidar, specifically useful for pinpointing forest structure, can also be exploited as a support for tree species richness assessment.     

Modeling Mediterranean forest structure using airborne laser scanning data

The conservation of biological diversity is recognized as a fundamental component of sustainable development, and forests contribute greatly to its preservation. Structural complexity increases the potential biological diversity of a forest by creating multiple niches that can host a wide variety of species. To facilitate greater understanding of the contributions of forest structure to forest biological diversity, we modeled relationships between 14 forest structure variables and airborne laser scanning (ALS) data for two Italian study areas representing two common Mediterranean forests, conifer plantations and coppice oaks subjected to irregular intervals of unplanned and non-standard silvicultural interventions. The objectives were twofold: (i) to compare model prediction accuracies when using two types of ALS metrics, echo-based metrics and canopy height model (CHM)-based metrics, and (ii) to construct inferences in the form of confidence intervals for large area structural complexity parameters.Our results showed that the effects of the two study areas on accuracies were greater than the effects of the two types of ALS metrics. In particular, accuracies were less for the more complex study area in terms of species composition and forest structure. However, accuracies achieved using the echo-based metrics were only slightly greater than when using the CHM-based metrics, thus demonstrating that both options yield reliable and comparable results. Accuracies were greatest for dominant height (Hd) (R² = 0.91; RMSE% = 8.2%) and mean height weighted by basal area (R² = 0.83; RMSE% = 10.5%) when using the echo-based metrics, 99th percentile of the echo height distribution and interquantile distance. For the forested area, the generalized regression (GREG) estimate of mean Hd was similar to the simple random sampling (SRS) estimate, 15.5 m for GREG and 16.2 m SRS. Further, the GREG estimator with standard error of 0.10 m was considerable more precise than the SRS estimator with standard error of 0.69 m.     

Augmenting forest inventory attributes with geometric optical modelling in support of regional susceptibility assessments to bark beetle infestations

Assessment of the susceptibility of forests to mountain pine beetle (Dendroctonus ponderosae Hopkins) infestation is based upon an understanding of the characteristics that predispose the stands to attack. These assessments are typically derived from conventional forest inventory data; however, this information often represents only managed forest areas. It does not cover areas such as forest parks or conservation regions and is often not regularly updated resulting in an inability to assess forest susceptibility. To address these shortcomings, we demonstrate how a geometric optical model (GOM) can be applied to Landsat-5 Thematic Mapper (TM) imagery (30 m spatial resolution) to estimate stand-level susceptibility to mountain pine beetle attack. Spectral mixture analysis was used to determine the proportion of sunlit canopy and background, and shadow of each Landsat pixel enabling per pixel estimates of attributes required for model inversion. Stand structural attributes were then derived from inversion of the geometric optical model and used as basis for susceptibility mapping. Mean stand density estimated by the geometric optical model was 2753 (standard deviation ± 308) stems per hectare and mean horizontal crown radius was 2.09 (standard deviation ± 0.11) metres. When compared to equivalent forest inventory attributes, model predictions of stems per hectare and crown radius were shown to be reasonably estimated using a Kruskal–Wallis ANOVA (p < 0.001). These predictions were then used to create a large area map that provided an assessment of the forest area susceptible to mountain pine beetle damage.     

Detecting pruning of individual stems using Airborne Laser Scanning data captured from an Unmanned Aerial Vehicle

Modern forest management involves implementing optimal pruning regimes. These regimes aim to achieve the highest quality timber in the shortest possible rotation period. Although a valuable addition to forest management activities, tracking the application of these treatments in the field to ensure best practice management is not economically viable. This paper describes the use of Airborne Laser Scanner (ALS) data to track the rate of pruning in a Eucalyptus globulus stand. Data is obtained from an Unmanned Aerial Vehicle (UAV) and we describe automated processing routines that provide a cost-effective alternative to field sampling. We manually prune a 500 m² plot to 2.5 m above the ground at rates of between 160 and 660 stems/ha. Utilising the high density ALS data, we first derived crown base height (CBH) with an RMSE of 0.60 m at each stage of pruning. Variability in the measurement of CBH resulted in both false positive (mean rate of 11%) and false negative detection (3.5%), however, detected rates of pruning of between 96% and 125% of the actual rate of pruning were achieved. The successful automated detection of pruning within this study highlights the suitability of UAV laser scanning as a cost-effective tool for monitoring forest management activities.     

Timber production assessment of a plantation forest: An integrated framework with field-based inventory,multi-source remote sensing data and forest management history

Timber production is the purpose for managing plantation forests, and its spatial and quantitative information is critical for advising management strategies. Previous studies have focused on growing stock volume (GSV), which represents the current potential of timber production, yet few studies have investigated historical process-harvested timber. This resulted in a gap in a synthetical ecosystem service assessment of timber production. In this paper, we established a Management Process–based Timber production (MPT) framework to integrate the current GSV and the harvested timber derived from historical logging regimes, trying to synthetically assess timber production for a historical period. In the MPT framework, age-class and current GSV determine the times of historical thinning and the corresponding harvested timber, by using a “space-for-time” substitution. The total timber production can be estimated by the historical harvested timber in each thinning and the current GSV. To test this MPT framework, an empirical study on a larch plantation (LP) with area of 43,946 ha was conducted in North China for a period from 1962 to 2010. Field-based inventory data was integrated with ALOS PALSAR (Advanced Land-Observing Satellite Phased Array L-band Synthetic Aperture Radar) and Landsat-8 OLI (Operational Land Imager) data for estimating the age-class and current GSV of LP. The random forest model with PALSAR backscatter intensity channels and OLI bands as input predictive variables yielded an accuracy of 67.9% with a Kappa coefficient of 0.59 for age-class classification. The regression model using PALSAR data produced a root mean square error (RMSE) of 36.5 m³ ha⁻¹. The total timber production of LP was estimated to be 7.27 × 10⁶ m³, with 4.87 × 10⁶ m³ in current GSV and 2.40 × 10⁶ m³ in harvested timber through historical thinning. The historical process-harvested timber accounts to 33.0% of the total timber production, which component has been neglected in the assessments for current status of plantation forests. Synthetically considering the RMSE for predictive GSV and misclassification of age-class, the error in timber production were supposed to range from −55.2 to 56.3 m³ ha⁻¹. The MPT framework can be used to assess timber production of other tree species at a larger spatial scale, providing crucial information for a better understanding of forest ecosystem service.     

Mapping the expansion of galamsey gold mines in the cocoa growing area of Ghana using optical remote sensing

Artisanal gold mining (galamsey) and cocoa farming are essential sources of income for local populations in Ghana. Unfortunately the former poses serious threats to the environment and human health, and conflicts with cocoa farming and other livelihoods. Timely and spatially referenced information on the extent of galamsey is needed to understand and limit the negative impacts of mining. To address this, we use multi-date UK-DMC2 satellite images to map the extent and expansion of galamsey from 2011 to 2015. We map the total area of galamsey in 2013 over the cocoa growing area, using k-means clustering on a cloud-free 2013 image with strong spectral contrast between galamsey and the surrounding vegetation. We also process a pair of hazy images from 2011 and 2015 with Multivariate Alteration Detection to map the 2011–2015 galamsey expansion in a subset, labelled the change area. We use a set of visually interpreted random sample points to compute bias-corrected area estimates. We also delineate an indicative impact zone of pollution proportional to the density of galamsey, assuming a maximum radius of 10 km. In the cocoa growing area of Ghana, the estimated total area of galamsey in 2013 is 27,839 ha with an impact zone of 551,496 ha. In the change area, galamsey has more than tripled between 2011 and 2015, resulting in 603 ha of direct encroachment into protected forest reserves. Assuming the same growth rate for the rest of the cocoa growing area, the total area of galamsey in 2015 is estimated at 43,879 ha. Galamsey is developing along most of the river network (Offin, Ankobra, Birim, Anum, Tano), with downstream pollution affecting both land and water.     

Forest canopy height estimation using ICESat/GLAS data and error factor analysis in Hokkaido,Japan

Spaceborne light detection and ranging (LiDAR) enables us to obtain information about vertical forest structure directly, and it has often been used to measure forest canopy height or above-ground biomass. However, little attention has been given to comparisons of the accuracy of the different estimation methods of canopy height or to the evaluation of the error factors in canopy height estimation. In this study, we tested three methods of estimating canopy height using the Geoscience Laser Altimeter System (GLAS) onboard NASA’s Ice, Cloud, and land Elevation Satellite (ICESat), and evaluated several factors that affected accuracy. Our study areas were Tomakomai and Kushiro, two forested areas on Hokkaido in Japan. The accuracy of the canopy height estimates was verified by ground-based measurements. We also conducted a multivariate analysis using quantification theory type I (multiple-regression analysis of qualitative data) and identified the observation conditions that had a large influence on estimation accuracy. The method using the digital elevation model was the most accurate, with a root-mean-square error (RMSE) of 3.2 m. However, GLAS data with a low signal-to-noise ratio (⩽10.0) and that taken from September to October 2009 had to be excluded from the analysis because the estimation accuracy of canopy height was remarkably low. After these data were excluded, the multivariate analysis showed that surface slope had the greatest effect on estimation accuracy, and the accuracy dropped the most in steeply sloped areas. We developed a second model with two equations to estimate canopy height depending on the surface slope, which improved estimation accuracy (RMSE = 2.8 m). These results should prove useful and provide practical suggestions for estimating forest canopy height using spaceborne LiDAR.     

Historical forest biomass dynamics modelled with Landsat spectral trajectories

Estimation of forest aboveground biomass (AGB) is informative of the role of forest ecosystems in local and global carbon budgets. There is a need to retrospectively estimate biomass in order to establish a historical baseline and enable reporting of change. In this research, we used temporal spectral trajectories to inform on forest successional development status in support of modelling and mapping of historic AGB for Mediterranean pines in central Spain. AGB generated with ground plot data from the Spanish National Forest Inventory (NFI), representing two collection periods (1990 and 2000), are linked with static and dynamic spectral data as captured by Landsat Thematic Mapper (TM) and Enhanced Thematic Mapper Plus (ETM+) sensors over a 25 year period (1984–2009). The importance of forest structural complexity on the relationship between AGB and spectral vegetation indices is revealed by the analysis of wavelet transforms. Two-dimensional (2D) wavelet transforms support the identification of spectral trajectory patterns of forest stands that in turn, are associated with traits of individual NFI plots, using a flexible algorithm sensitive to capturing time series similarity. Single-date spectral indices, temporal trajectories, and temporal derivatives associated with succession are used as input variables to non-parametric decision trees for modelling, estimation, and mapping of AGB and carbon sinks over the entire study area. Results indicate that patterns of change found in Normalized Difference Vegetation Index (NDVI) values are associated and relate well to classes of forest AGB. The Tasseled Cap Angle (TCA) index was found to be strongly related with forest density, although the related patterns of change had little relation with variability in historic AGB. By scaling biomass models through small (∼2.5 ha) spatial objects defined by spectral homogeneity, the AGB dynamics in the period 1990–2000 are mapped (70% accuracy when validated with plot values of change), revealing an increase of 18% in AGB irregularly distributed over 814 km² of pines. The accumulation of C calculated in AGB was on average 0.65 t ha⁻¹ y⁻¹, equivalent to a fixation of 2.38 t ha⁻¹ y⁻¹ of carbon dioxide.     

Aboveground biomass mapping in French Guiana by combining remote sensing,forest inventories and environmental data

Mapping forest aboveground biomass (AGB) has become an important task, particularly for the reporting of carbon stocks and changes. AGB can be mapped using synthetic aperture radar data (SAR) or passive optical data. However, these data are insensitive to high AGB levels (>150 Mg/ha, and >300 Mg/ha for P-band), which are commonly found in tropical forests. Studies have mapped the rough variations in AGB by combining optical and environmental data at regional and global scales. Nevertheless, these maps cannot represent local variations in AGB in tropical forests. In this paper, we hypothesize that the problem of misrepresenting local variations in AGB and AGB estimation with good precision occurs because of both methodological limits (signal saturation or dilution bias) and a lack of adequate calibration data in this range of AGB values. We test this hypothesis by developing a calibrated regression model to predict variations in high AGB values (mean >300 Mg/ha) in French Guiana by a methodological approach for spatial extrapolation with data from the optical geoscience laser altimeter system (GLAS), forest inventories, radar, optics, and environmental variables for spatial inter- and extrapolation. Given their higher point count, GLAS data allow a wider coverage of AGB values. We find that the metrics from GLAS footprints are correlated with field AGB estimations (R² = 0.54, RMSE = 48.3 Mg/ha) with no bias for high values. First, predictive models, including remote-sensing, environmental variables and spatial correlation functions, allow us to obtain “wall-to-wall” AGB maps over French Guiana with an RMSE for the in situ AGB estimates of ∼50 Mg/ha and R² = 0.66 at a 1-km grid size. We conclude that a calibrated regression model based on GLAS with dependent environmental data can produce good AGB predictions even for high AGB values if the calibration data fit the AGB range. We also demonstrate that small temporal and spatial mismatches between field data and GLAS footprints are not a problem for regional and global calibrated regression models because field data aim to predict large and deep tendencies in AGB variations from environmental gradients and do not aim to represent high but stochastic and temporally limited variations from forest dynamics. Thus, we advocate including a greater variety of data, even if less precise and shifted, to better represent high AGB values in global models and to improve the fitting of these models for high values.     

Assessing the utility WorldView-2 imagery for tree species mapping in South African subtropical humid forest and the conservation implications: Dukuduku forest patch as case study

Indigenous forest biome in South Africa is highly fragmented into patches of various sizes (most patches < 1 km²). The utilization of timber and non-timber resources by poor rural communities living around protected forest patches produce subtle changes in the forest canopy which can be hardly detected on a timely manner using traditional field surveys. The aims of this study were to assess: (i) the utility of very high resolution (VHR) remote sensing imagery (WorldView-2, 0.5–2 m spatial resolution) for mapping tree species and canopy gaps in one of the protected subtropical coastal forests in South Africa (the Dukuduku forest patch (ca.3200 ha) located in the province of KwaZulu-Natal) and (ii) the implications of the map products to forest conservation. Three dominant canopy tree species namely, Albizia adianthifolia, Strychnos spp. and Acacia spp., and canopy gap types including bushes (grass/shrubby), bare soil and burnt patches were accurately mapped (overall accuracy = 89.3 ± 2.1%) using WorldView-2 image and support vector machine classifier. The maps revealed subtle forest disturbances such as bush encroachment and edge effects resulting from forest fragmentation by roads and a power-line. In two stakeholders’ workshops organised to assess the implications of the map products to conservation, participants generally agreed amongst others implications that the VHR maps provide valuable information that could be used for implementing and monitoring the effects of rehabilitation measures. The use of VHR imagery is recommended for timely inventorying and monitoring of the small and fragile patches of subtropical forests in Southern Africa.     

Quantifying post-fire fallen trees using multi-temporal lidar

Massive tree-felling due to root damage is a common fire effect on burnt areas in Scandinavia, but has so far not been analyzed in detail. Here we explore if pre- and post-fire lidar data can be used to estimate the proportion of fallen trees. The study was carried out within a large (14,000 ha) area in central Sweden burnt in August 2014, where we had access to airborne lidar data from both 2011 and 2015. Three data-sets of predictor variables were tested: POST (post-fire lidar metrics), DIF (difference between post- and pre-fire lidar metrics) and combination of those two (POST_DIF). Fractional logistic regression was used to predict the proportion of fallen trees. Training data consisted of 61 plots, where the number of fallen and standing trees was calculated both in the field and with interpretation of drone images. The accuracy of the best model was tested based on 100 randomly selected validation plots with a size of 25 × 25 m.Our results showed that multi-temporal lidar together with field-collected training data can be used for quantifying post-fire tree felling over large areas. Several height-, density- and intensity metrics correlated with the proportion of fallen trees. The best model combined metrics from both datasets (POST_DIF), resulting in a RMSE of 0.11. Results were slightly poorer in the validation plots with RMSE of 0.18 using pixel size of 12.5 m and RMSE of 0.15 using pixel size of 6.25 m. Our model performed least well for stands that had been exposed to high-intensity crown fire. This was likely due to the low amount of echoes from the standing black tree skeletons. Wall-to-wall maps produced with this model can be used for landscape level analysis of fire effects and to explore the relationship between fallen trees and forest structure, soil type, fire intensity or topography.     

A Gibbs sampling disaggregation model for orographic precipitation

Hydrological applications in complex topographic areas need high spatial resolution precipitation data. Some daily high-resolution products are now available for recent past data, even in complex terrain. While the spatial resolution of Regional Climate Models (RCMs) and operational meteorological models are becoming increasingly fine, there still exists a mismatch between the spatial resolutions of observed or estimated recent past data and simulated or forecasted precipitation.Statistical disaggregation models can generate precipitation on a high-resolution grid using as input a mesoscale precipitation grid (e.g., RCM or meteorological grid). In this paper, a Gibbs sampling disaggregation model previously developed for flat areas is adapted to account for topography. Only one variable, the topographic anomaly, is added to the original model. The model is applied on a 300 km × 300 km area in the northwestern United States, covering the Olympic Mountains and the Cascade Range. Daily high-resolution precipitation data for the 2002–2005 period are used to estimate the model parameters. Using 750 days taken from the 2006–2008 period, 36, 52-km grid boxes are disaggregated on 4.3-, 8.7-, 13-, 17.3- and 26-km grids; each day being simulated nine times. Thank to the Gibbs sampling algorithm, the original model, which does not account for topography, is able to capture the mesoscale topographic structure of the daily precipitation, while the adapted model accounting for topography is better suited to recreate the local impact of topography on interannual means, interday standard deviation, and maximum values. The model outputs could be used by hydrological modelers who need high-resolution precipitation data in complex topographic area application.     

Major forest changes and land cover transitions based on plant functional types derived from the ESA CCI Land Cover product

Land use and land cover change are of prime concern due to their impacts on CO₂ emissions, climate change and ecological services. New global land cover products at 300 m resolution from the European Space Agency (ESA) Climate Change Initiative Land Cover (CCI LC) project for epochs centered around 2000, 2005 and 2010 were analyzed to investigate forest area change and land cover transitions. Plant functional types (PFTs) fractions were derived from these land cover products according to a conversion table. The gross global forest loss between 2000 and 2010 is 172,171 km², accounting for 0.6% of the global forest area in year 2000. The forest changes are mainly distributed in tropical areas such as Brazil and Indonesia. Forest gains were only observed between 2005 and 2010 with a global area of 9844 km², mostly from crops in Southeast Asia and South America. The predominant PFT transition is deforestation from forest to crop, accounting for four-fifths of the total increase of cropland area between 2000 and 2010. The transitions from forest to bare soil, shrub, and grass also contributed strongly to the total areal change in PFTs. Different PFT transition matrices and composition patterns were found in different regions. The highest fractions of forest to bare soil transitions were found in the United States and Canada, reflecting forest management practices. Most of the degradation from grassland and shrubland to bare soil occurred in boreal regions. The areal percentage of forest loss and land cover transitions generally decreased from 2000–2005 to 2005–2010. Different data sources and uncertainty in the conversion factors (converting from original LC classes to PFTs) contribute to the discrepancy in the values of change in absolute forest area.     

Modeling soil parameters using hyperspectral image reflectance in subtropical coastal wetlands

Developing spectral models of soil properties is an important frontier in remote sensing and soil science. Several studies have focused on modeling soil properties such as total pools of soil organic matter and carbon in bare soils. We extended this effort to model soil parameters in areas densely covered with coastal vegetation. Moreover, we investigated soil properties indicative of soil functions such as nutrient and organic matter turnover and storage. These properties include the partitioning of mineral and organic soil between particulate (>53 μm) and fine size classes, and the partitioning of soil carbon and nitrogen pools between stable and labile fractions. Soil samples were obtained from Avicennia germinans mangrove forest and Juncus roemerianus salt marsh plots on the west coast of central Florida. Spectra corresponding to field plot locations from Hyperion hyperspectral image were extracted and analyzed. The spectral information was regressed against the soil variables to determine the best single bands and optimal band combinations for the simple ratio (SR) and normalized difference index (NDI) indices. The regression analysis yielded levels of correlation for soil variables with R² values ranging from 0.21 to 0.47 for best individual bands, 0.28 to 0.81 for two-band indices, and 0.53 to 0.96 for partial least-squares (PLS) regressions for the Hyperion image data. Spectral models using Hyperion data adequately (RPD > 1.4) predicted particulate organic matter (POM), silt + clay, labile carbon (C), and labile nitrogen (N) (where RPD = ratio of standard deviation to root mean square error of cross-validation [RMSECV]). The SR (0.53 μm, 2.11 μm) model of labile N with R² = 0.81, RMSECV= 0.28, and RPD = 1.94 produced the best results in this study. Our results provide optimism that remote-sensing spectral models can successfully predict soil properties indicative of ecosystem nutrient and organic matter turnover and storage, and do so in areas with dense canopy cover.     

Very high spatial resolution optical and radar imagery in tracking water level fluctuations of a small inland reservoir

Tracking water level fluctuations in small lakes and reservoirs is important in order to better understand and manage these ecosystems. A geographic object-based image analysis (GEOBIA) method using very high spatial and temporal resolution optical (Pléiades) and radar (COSMO-SkyMed and TerraSAR-X) remote sensing imagery is presented here which (1) tracks water level fluctuations via variations in water surface area and (2) avoids common difficulties found in using single-band radar images for water-land image classification. Results are robust, with over 98% of image surface area correctly classified into land or water, R² = 0.963 and RMSE = 0.42 m for a total water level fluctuation range of 5.94 m. Multispectral optical imagery is found to be more straightforward in producing results than single-band radar imagery, but the latter crucially increase temporal resolution to the point where fluctuations can be satisfactorily tracked in time. Moreover, an analysis suggests that high and medium spatial resolution imagery is sufficient, in at least some cases, in tracking the water level fluctuations of small inland reservoirs. Finally, limitations of the methodology presented here are briefly discussed along with potential solutions to overcome them.     

Spatially detailed retrievals of spring phenology from single-season high-resolution image time series

Vegetation indices derived from satellite image time series have been extensively used to estimate the timing of phenological events like season onset. Medium spatial resolution (≥250 m) satellite sensors with daily revisit capability are typically employed for this purpose. In recent years, phenology is being retrieved at higher resolution (≤30 m) in response to increasing availability of high-resolution satellite data. To overcome the reduced acquisition frequency of such data, previous attempts involved fusion between high- and medium-resolution data, or combinations of multi-year acquisitions in a single phenological reconstruction. The objectives of this study are to demonstrate that phenological parameters can now be retrieved from single-season high-resolution time series, and to compare these retrievals against those derived from multi-year high-resolution and single-season medium-resolution satellite data. The study focuses on the island of Schiermonnikoog, the Netherlands, which comprises a highly-dynamic saltmarsh, dune vegetation, and agricultural land. Combining NDVI series derived from atmospherically-corrected images from RapidEye (5 m-resolution) and the SPOT5 Take5 experiment (10m-resolution) acquired between March and August 2015, phenological parameters were estimated using a function fitting approach. We then compared results with phenology retrieved from four years of 30 m Landsat 8 OLI data, and single-year 100 m Proba-V and 250 m MODIS temporal composites of the same period. Retrieved phenological parameters from combined RapidEye/SPOT5 displayed spatially consistent results and a large spatial variability, providing complementary information to existing vegetation community maps. Retrievals that combined four years of Landsat observations into a single synthetic year were affected by the inclusion of years with warmer spring temperatures, whereas adjustment of the average phenology to 2015 observations was only feasible for a few pixels due to cloud cover around phenological transition dates. The Proba-V and MODIS phenology retrievals scaled poorly relative to their high-resolution equivalents, indicating that medium-resolution phenology retrievals need to be interpreted with care, particularly in landscapes with fine-scale land cover variability.     

Object-based class modelling for multi-scale riparian forest habitat mapping

Object-based class modelling allows for mapping complex, hierarchical habitat systems. The riparian zone, including forests, represents such a complex ecosystem. Forests within riparian zones are biologically high productive and characterized by a rich biodiversity; thus considered of high community interest with an imperative to be protected and regularly monitored. Satellite earth observation (EO) provides tools for capturing the current state of forest habitats such as forest composition including intermixture of non-native tree species. Here we present a semi-automated object based image analysis (OBIA) approach for the mapping of riparian forests by applying class modelling of habitats based on the European Nature Information System (EUNIS) habitat classifications and the European Habitats Directive (HabDir) Annex 1. A very high resolution (VHR) WorldView-2 satellite image provided the required spatial and spectral details for a multi-scale image segmentation and rule-base composition to generate a six-level hierarchical representation of riparian forest habitats. Thereby habitats were hierarchically represented within an image object hierarchy as forest stands, stands of homogenous tree species and single trees represented by sunlit tree crowns. 522 EUNIS level 3 (EUNIS-3) habitat patches with a mean patch size (MPS) of 12,349.64 m² were modelled from 938 forest stand patches (MPS = 6868.20 m²) and 43,742 tree stand patches (MPS = 140.79 m²). The delineation quality of the modelled EUNIS-3 habitats (focal level) was quantitatively assessed to an expert-based visual interpretation showing a mean deviation of 11.71%.