Selecting landscape metrics as indicators of spatial heterogeneity—A comparison among Greek landscapes

This paper investigates the spatial heterogeneity of three landscapes along an altitudinal gradient and different human land use. The main aim was the identification of appropriate landscape indicators using different extents. ASTER image was used to create a land cover map consisting of three landscapes which differed in altitude and land use. A number of landscape metrics quantifying patch complexity, configuration, diversity and connectivity were derived from the thematic map at the landscape level. There were significant differences among the three landscapes regarding these four aspects of landscape heterogeneity. The analysis revealed a specific pattern of land use where lowlands are being increasingly utilized by humans (percentage of agricultural land = 65.84%) characterized by physical connectedness (high values of Patch Cohesion Index) and relatively simple geometries (low values of fractal dimension index). The landscape pattern of uplands was found to be highly diverse based upon the Shannon Diversity index. After selecting the scale (600 ha) where metrics values stabilized, it was shown that metrics were more correlated at the small scale of 60 ha. From the original 24 metrics, 14 individual metrics with high Spearman correlation coefficient and Variance Inflation Factor criterion were eliminated, leaving 10 representative metrics for subsequent analysis. Data reduction analysis showed that Patch Density, Area-Weighted Mean Fractal Dimension Index and Patch Cohesion Index are suitable to describe landscape patterns irrespective of the scale. A systematic screening of these metrics could enhance a deeper understanding of the results obtained by them and contribute to a sustainable landscape management of Mediterranean landscapes.     

Methodology for estimating availability of cloud-free image composites: A case study for southern Canada

Image composites are often used for earth surface phenomena studies at regional or national level. The compromise between residual clouds and the length of compositing period is a necessary corollary to the choice of satellite optical data for monitoring earth surface phenomena dynamics. This paper introduced a methodology for estimating availability of cloud-free image composites for optical sensors with various revisiting intervals, using MODIS MOD06 L2 cloud fraction product in the period of 2000–2008. The methodology starts with downscaling of the cloud fraction product to 1 km × 1 km cloud cover binary images. The binary images are then used for the exploration of spatial and temporal characteristics of cloud dynamics, and subsequently for the simulation of cloud-free composite availability with various revisiting intervals of optical sensors. Using Canada's southern provinces as an application case, the study explored several factors important for the design of environmental monitoring system using optical sensors of earth observation, in particular, cloud dynamics and its inter-annual variability, sensors’ revisiting intervals, and cloud-free threshold for targeting composites. While the cloud images used in the analysis are at 1 km × 1 km resolution, our analysis suggests that the simulated availabilities of cloud-free image composites may also provide reasonable estimates for optical sensors with higher than 1 km × 1 km resolution, though the closer to 1 km × 1 km resolution the optical sensor, the more pertinent the application. Also, the methodology can be parameterised to different temporal period and different spatial region, depending on applications.     

A Spatio-Temporal Enhancement Method for medium resolution LAI (STEM-LAI)

Satellite remote sensing has been used successfully to map leaf area index (LAI) across landscapes, but advances are still needed to exploit multi-scale data streams for producing LAI at both high spatial and temporal resolution. A multi-scale Spatio-Temporal Enhancement Method for medium resolution LAI (STEM-LAI) has been developed to generate 4-day time-series of Landsat-scale LAI from existing medium resolution LAI products. STEM-LAI has been designed to meet the demands of applications requiring frequent and spatially explicit information, such as effectively resolving rapidly evolving vegetation dynamics at sub-field (30 m) scales. In this study, STEM-LAI is applied to Moderate Resolution Imaging Spectroradiometer (MODIS) based LAI data and utilizes a reference-based regression tree approach for producing MODIS-consistent, but Landsat-based, LAI. The Spatial and Temporal Adaptive Reflectance Fusion Model (STARFM) is used to interpolate the downscaled LAI between Landsat acquisition dates, providing a high spatial and temporal resolution improvement over existing LAI products. STARFM predicts high resolution LAI by blending MODIS and Landsat based information from a common acquisition date, with MODIS data from a prediction date. To demonstrate its capacity to reproduce fine-scale spatial features observed in actual Landsat LAI, the STEM-LAI approach is tested over an agricultural region in Nebraska. The implementation of a 250 m resolution LAI product, derived from MODIS 1 km data and using a scale consistent approach based on the Normalized Difference Vegetation Index (NDVI), is found to significantly improve accuracies of spatial pattern prediction, with the coefficient of efficiency (E) ranging from 0.77–0.94 compared to 0.01–0.85 when using 1 km LAI inputs alone. Comparisons against an 11-year record of in-situ measured LAI over maize and soybean highlight the utility of STEM-LAI in reproducing observed LAI dynamics (both characterized by r² = 0.86) over a range of plant development stages. Overall, STEM-LAI represents an effective downscaling and temporal enhancement mechanism that predicts in-situ measured LAI better than estimates derived through linear interpolation between Landsat acquisitions. This is particularly true when the in-situ measurement date is greater than 10 days from the nearest Landsat acquisition, with prediction errors reduced by up to 50%. With a streamlined and completely automated processing interface, STEM-LAI represents a flexible tool for LAI disaggregation in space and time that is adaptable to different land cover types, landscape heterogeneities, and cloud cover conditions.     

SAR interferometry and optical remote sensing for analysis of co-seismic deformation,source characteristics and mass wasting pattern of Lushan (China,April 2013) earthquake

Co-seismic deformation associated with the Lushan (China) earthquake that occurred along the south-western segment of the Longmenshan Fault Zone (LFZ) on the 20th April 2013 has been estimated by differential interferometric SAR (DInSAR) technique using Radarsat-2 data. The Lushan earthquake resulted in the deformation of the Sichuan basin and the Longmenshan ranges in proximity to the LFZ. The line of sight (LOS) displacement values obtained from DInSAR technique mainly range between −4.0 cm to +3.0 cm. The western Sichuan basin shows oblique westward movement with predominant downward component in areas farther from LFZ and predominant westward component over the downward movement in areas closer to the source fault. Inversion modelling has been used to derive the seismic source characteristics from DInSAR derived deformation values using elastic dislocation source type. The linear inversion model converged at a double-fault source solution consisting of a deeper, steep, NW dipping fault plane-1 of 60 km × 16 km dimension and a shallower, gentle, NW dipping fault plane-2 of 60 km × 15 km dimension, with distributed slip values varying between 0 to 2.26 m. These fault planes (fault planes-1 and -2) coincide with the Dachuan-Shuangshi fault and the buried Range Front Fault, respectively. The inversion model gives a moment magnitude of 6.81 and the geodetic moment of 2.07 × 10¹⁹ Nm, comparable to those given in literature, derived using teleseismic body wave data. Thus DInSAR technique helped to quantify the co-seismic deformation and to retrieve the source characteristics from the estimated deformation values. The study also evaluated the distribution pattern of earthquake induced landslides (EIL) triggered fresh or re-activated during the Lushan earthquake and found that they show spatial association with the seismic source zone and also with various pre-conditioning factors of slope instability.     

Mapping seagrass coverage and spatial patterns with high spatial resolution IKONOS imagery

Seagrass habitats in subtidal coastal waters provide a variety of ecosystem functions and services and there is an increasing need to acquire information on spatial and temporal dynamics of this resource. Here, we explored the capability of IKONOS (IKO) data of high resolution (4 m) for mapping seagrass cover [submerged aquatic vegetation (%SAV) cover] along the mid-western coast of Florida, USA. We also compared seagrass maps produced with IKO data with that obtained using the Landsat TM sensor with lower resolution (30 m). Both IKO and TM data, collected in October 2009, were preprocessed to calculate water depth invariant bands to normalize the effect of varying depth on bottom spectra recorded by the two satellite sensors and further the textural information was extracted from IKO data. Our results demonstrate that the high resolution IKO sensor produced a higher accuracy than the TM sensor in a three-class % SAV cover classification. Of note is that the OA of %SAV cover mapping at our study area created with IKO data was 5–20% higher than that from other studies published. We also examined the spatial distribution of seagrass over a spatial range of 4–240 m using the Ripley’s K function [L(d)] and IKO data that represented four different grain sizes [4 m (one IKO pixel), 8 m (2 × 2 IKO pixels), 12 m (3 × 3 IKO pixels), and 16 m (4 × 4 IKO pixels)] from moderate-dense seagrass cover along a set of six transects. The Ripley’s K metric repeatedly indicated that seagrass cover representing 4 m × 4 m pixels displayed a dispersed (or slightly dispersed) pattern over distances of <4–8 m, and a random or slightly clustered pattern of cover over 9–240 m. The spatial pattern of seagrass cover created with the three additional grain sizes (i.e., 2 × 24 m IKO pixels, 3 × 34 m IKO pixels, and 4 × 4 m IKO pixels) show a dispersed (or slightly dispersed) pattern across 4–32 m and a random or slightly clustered pattern across 33–240 m. Given the first report on using satellite observations to quantify seagrass spatial patterns at a spatial scale from 4 m to 240 m, our novel analyses of moderate-dense SAV cover utilizing Ripley’s K function illustrate how data obtained from the IKO sensor revealed seagrass spatial information that would be undetected by the TM sensor with a 30 m pixel size. Use of the seagrass classification scheme here, along with data from the IKO sensor with enhanced resolution, offers an opportunity to synoptically record seagrass cover dynamics at both small and large spatial scales.     

Measuring and monitoring linear woody features in agricultural landscapes through earth observation data as an indicator of habitat availability

The loss of natural habitats and the loss of biological diversity is a global problem affecting all ecosystems including agricultural landscapes. Indicators of biodiversity can provide standardized measures that make it easier to compare and communicate changes to an ecosystem. In agricultural landscapes the amount and variety of available habitat is directly correlated with biodiversity levels. Linear woody features (LWF), including hedgerows, windbreaks, shelterbelts as well as woody shrubs along fields, roads and watercourses, play a vital role in supporting biodiversity as well as serving a wide variety of other purposes in the ecosystem. Earth observation can be used to quantify and monitor LWF across the landscape. While individual features can be manually mapped, this research focused on the development of methods using line intersect sampling (LIS) for estimating LWF as an indicator of habitat availability in agricultural landscapes. The methods are accurate, efficient, repeatable and provide robust results. Methods were tested over 9.5 Mha of agricultural landscape in the Canadian Mixedwood Plains ecozone. Approximately 97,000 km of LWF were estimated across this landscape with results useable both at a regional reporting scale, as well as mapped across space for use in wildlife habitat modelling or other landscape management research. The LIS approach developed here could be employed at a variety of scales in particular for large regions and could be adapted for use as a national scale indicator of habitat availability in heavily disturbed agricultural landscape.     

Determination of the crop row orientations from Formosat-2 multi-temporal and panchromatic images

This paper presents a technique developed for the retrieval of the orientation of crop rows, over anthropic lands dedicated to agriculture in order to further improve estimate of crop production and soil erosion management. Five crop types are considered: wheat, barley, rapeseed, sunflower, corn and hemp. The study is part of the multi-sensor crop-monitoring experiment, conducted in 2010 throughout the agricultural season (MCM’10) over an area located in southwestern France, near Toulouse. The proposed methodology is based on the use of satellite images acquired by Formosat-2, at high spatial resolution in panchromatic and multispectral modes (with spatial resolution of 2 and 8 m, respectively). Orientations are derived and evaluated for each image and for each plot, using directional spatial filters (45° and 135°) and mathematical morphology algorithms. “Single-date” and “multi-temporal” approaches are considered. The single-date analyses confirm the good performances of the proposed method, but emphasize the limitation of the approach for estimating the crop row orientation over the whole landscape with only one date. The multi-date analyses allow (1) determining the most suitable agricultural period for the detection of the row orientations, and (2) extending the estimation to the entire footprint of the study area. For the winter crops (wheat, barley and rapeseed), best results are obtained with images acquired just after harvest, when surfaces are covered by stubbles or during the period of deep tillage (0.27 > R² > 0.99 and 7.15° > RMSE > 43.02°). For the summer crops (sunflower, corn and hemp), results are strongly crop and date dependents (0 > R² > 0.96, 10.22° > RMSE > 80°), with a well-marked impact of flowering, irrigation equipment and/or maximum crop development. Last, the extent of the method to the whole studied zone allows mapping 90% of the crop row orientations (more than 45,000 ha) with an error inferior to 40°, associated to a confidence index ranging from 1 to 5 for each agricultural plot.     

An accurate and computationally efficient algorithm for ground peak identification in large footprint waveform LiDAR data

Large footprint waveform LiDAR sensors have been widely used for numerous airborne studies. Ground peak identification in a large footprint waveform is a significant bottleneck in exploring full usage of the waveform datasets. In the current study, an accurate and computationally efficient algorithm was developed for ground peak identification, called Filtering and Clustering Algorithm (FICA). The method was evaluated on Land, Vegetation, and Ice Sensor (LVIS) waveform datasets acquired over Central NY. FICA incorporates a set of multi-scale second derivative filters and a k-means clustering algorithm in order to avoid detecting false ground peaks. FICA was tested in five different land cover types (deciduous trees, coniferous trees, shrub, grass and developed area) and showed more accurate results when compared to existing algorithms. More specifically, compared with Gaussian decomposition, the RMSE ground peak identification by FICA was 2.82 m (5.29 m for GD) in deciduous plots, 3.25 m (4.57 m for GD) in coniferous plots, 2.63 m (2.83 m for GD) in shrub plots, 0.82 m (0.93 m for GD) in grass plots, and 0.70 m (0.51 m for GD) in plots of developed areas. FICA performance was also relatively consistent under various slope and canopy coverage (CC) conditions. In addition, FICA showed better computational efficiency compared to existing methods. FICA’s major computational and accuracy advantage is a result of the adopted multi-scale signal processing procedures that concentrate on local portions of the signal as opposed to the Gaussian decomposition that uses a curve-fitting strategy applied in the entire signal. The FICA algorithm is a good candidate for large-scale implementation on future space-borne waveform LiDAR sensors.     

Agriculture pest and disease risk maps considering MSG satellite data and land surface temperature

Pest risk maps for agricultural use are usually constructed from data obtained from in-situ meteorological weather stations, which are relatively sparsely distributed and are often quite expensive to install and difficult to maintain. This leads to the creation of maps with relatively low spatial resolution, which are very much dependent on interpolation methodologies. Considering that agricultural applications typically require a more detailed scale analysis than has traditionally been available, remote sensing technology can offer better monitoring at increasing spatial and temporal resolutions, thereby, improving pest management results and reducing costs. This article uses ground temperature, or land surface temperature (LST), data distributed by EUMETSAT/LSASAF (with a spatial resolution of 3 × 3 km (nadir resolution) and a revisiting time of 15 min) to generate one of the most commonly used parameters in pest modeling and monitoring: “thermal integral over air temperature (accumulated degree-days)”. The results show a clear association between the accumulated LST values over a threshold and the accumulated values computed from meteorological stations over the same threshold (specific to a particular tomato pest). The results are very promising and enable the production of risk maps for agricultural pests with a degree of spatial and temporal detail that is difficult to achieve using in-situ meteorological stations.     

Monitoring land surface albedo and vegetation dynamics using high spatial and temporal resolution synthetic time series from Landsat and the MODIS BRDF/NBAR/albedo product

Seasonal vegetation phenology can significantly alter surface albedo which in turn affects the global energy balance and the albedo warming/cooling feedbacks that impact climate change. To monitor and quantify the surface dynamics of heterogeneous landscapes, high temporal and spatial resolution synthetic time series of albedo and the enhanced vegetation index (EVI) were generated from the 500 m Moderate Resolution Imaging Spectroradiometer (MODIS) operational Collection V006 daily BRDF/NBAR/albedo products and 30 m Landsat 5 albedo and near-nadir reflectance data through the use of the Spatial and Temporal Adaptive Reflectance Fusion Model (STARFM). The traditional Landsat Albedo (Shuai et al., 2011) makes use of the MODIS BRDF/Albedo products (MCD43) by assigning appropriate BRDFs from coincident MODIS products to each Landsat image to generate a 30 m Landsat albedo product for that acquisition date. The available cloud free Landsat 5 albedos (due to clouds, generated every 16 days at best) were used in conjunction with the daily MODIS albedos to determine the appropriate 30 m albedos for the intervening daily time steps in this study. These enhanced daily 30 m spatial resolution synthetic time series were then used to track albedo and vegetation phenology dynamics over three Ameriflux tower sites (Harvard Forest in 2007, Santa Rita in 2011 and Walker Branch in 2005). These Ameriflux sites were chosen as they are all quite nearby new towers coming on line for the National Ecological Observatory Network (NEON), and thus represent locations which will be served by spatially paired albedo measures in the near future. The availability of data from the NEON towers will greatly expand the sources of tower albedometer data available for evaluation of satellite products. At these three Ameriflux tower sites the synthetic time series of broadband shortwave albedos were evaluated using the tower albedo measurements with a Root Mean Square Error (RMSE) less than 0.013 and a bias within the range of ±0.006. These synthetic time series provide much greater spatial detail than the 500 m gridded MODIS data, especially over more heterogeneous surfaces, which improves the efforts to characterize and monitor the spatial variation across species and communities. The mean of the difference between maximum and minimum synthetic time series of albedo within the MODIS pixels over a subset of satellite data of Harvard Forest (16 km by 14 km) was as high as 0.2 during the snow-covered period and reduced to around 0.1 during the snow-free period. Similarly, we have used STARFM to also couple MODIS Nadir BRDF Adjusted Reflectances (NBAR) values with Landsat 5 reflectances to generate daily synthetic times series of NBAR and thus Enhanced Vegetation Index (NBAR-EVI) at a 30 m resolution. While normally STARFM is used with directional reflectances, the use of the view angle corrected daily MODIS NBAR values will provide more consistent time series. These synthetic times series of EVI are shown to capture seasonal vegetation dynamics with finer spatial and temporal details, especially over heterogeneous land surfaces.     

Determination of the factors governing soil erodibility using hyperspectral visible and near-infrared reflectance spectroscopy

Soil erodibility, which is difficult to estimate and upscaling, was determined in this study using multiple spectral models of soil properties (soil organic matter (SOM), water-stable aggregates (WSA) > 0.25 mm, the geometric mean radius (D_g)). Herein, the soil erodibility indicators were calculated, and soil properties were quantitatively analyzed based on laboratory simulation experiments involving two selected contrasting soils. In addition, continuous wavelet transformation was applied to the reflectance spectra (350–2500 nm) of 65 soil samples from the study area. To build the relationship, the soil properties that control erodibility were identified prior to the spectral analysis. In this study, the SOM, D_g and WSA >0.25 mm were selected to represent the most significant soil properties controlling erodibility and describe the erodibility indicator based on a logarithmic regression model as a function of SOM or WSA > 0.25 mm. Five, six and three wavelet features were observed to calibrate the estimated soil properties model, and the best performance was obtained with a combination feature regression model for SOM (R² = 0.86, p < 0.01), D_g (R² = 0.79, p < 0.01) and WSA >0.25 mm (R² = 0.61, p < 0.01), respectively. One part of the wavelet features captured amplitude variations in the broad shape of the reflectance spectra, and another part captured variations in the shape and depth of the soil dry substances. The wavelet features for the validated dataset used to predict the SOM, WSA >0.25 mm and D_g were not significantly different compared with the calibrated dataset. The synthesized spectral models of soil properties, and the formation of a new equation for soil erodibility transformed from the spectral models of soil properties are presented in this study. These results show that a spectral analytical approach can be applied to complex datasets and provide new insights into emerging dynamic variation with erodibility estimation.     

Testing estimation of water surface in Italian rice district from MODIS satellite data

Recent changes in rice crop management within Northern Italy rice district led to a reduction of seeding in flooding condition, which may have an impact on reservoir water management and on the animal and plant communities that depend on the flooded paddies. Therefore, monitoring and quantifying the spatial and temporal variability of water presence in paddy fields is becoming important. In this study we present a method to estimate dynamics of presence of standing water (i.e. fraction of flooded area) in rice fields using MODIS data. First, we produced high resolution water presence maps from Landsat by thresholding the Normalised Difference Flood Index (NDFI) made: we made it by comparing five Landsat 8 images with field-obtained information about rice field status and water presence. Using these data we developed an empirical model to estimate the flooding fraction of each MODIS cell. Finally we validated the MODIS-based flooding maps with both Landsat and ground information. Results showed a good predictability of water surface from Landsat (OA = 92%) and a robust usability of MODIS data to predict water fraction (R² = 0.73, EF = 0.57, RMSE = 0.13 at 1 × 1 km resolution). Analysis showed that the predictive ability of the model decreases with the greening up of rice, so we used NDVI to automatically discriminate estimations for inaccurate cells in order to provide the water maps with a reliability flag. Results demonstrate that it is possible to monitor water dynamics in rice paddies using moderate resolution multispectral satellite data. The achievement is a proof of concept for the analysis of MODIS archives to investigate irrigation dynamics in the last 15 years to retrieve information for ecological and hydrological studies.     

Spectral detection of stress-related pigments in salt-lake succulent halophytic shrubs

The spectral detection of vegetation pigment concentrations has a high potential value, but it is still underdeveloped, especially for pigments other than chlorophylls. In this study, the seasonal pigment dynamics of two Tecticornia species (samphires; halophytic shrubs) from north-western Australia were correlated with spectral indices that best document the pigment changes over time. Pigment dynamics were assessed by analysing betacyanin, chlorophyll and carotenoid concentrations at plant level and by measuring reflectance at contrasting seasonal dates. Plant reflectance was used to define a new reflectance index that was most sensitive to the seasonal shifts in Tecticornia pigment concentrations. The two Tecticornia species turned from green to red-pinkish for the period March–August 2012 when betacyanins increased almost nine times in both species. Chlorophyll levels showed the opposite pattern to that of betacyanins, whereas carotenoid levels were relatively stable. Normalised difference indices correlated well with betacyanin (r = 0.805, using bands at 600 and 620 nm) and chlorophyll (r = 0.809, using bands at 737 and 726 nm). Using knowledge of chlorophyll concentrations slightly improved the ability of the spectral index to predict betacyanin concentration (r = 0.822 at bands 606 and 620 nm, in the case of chemically determined chlorophyll, r = 0.809 when using remotely sensed chlorophyll). Our results suggest that this new spectral index can reliably detect changes in betacyanin concentrations in vegetation, with potential applications in ecological studies and environmental impact monitoring.     

Estimating crop net primary production using national inventory data and MODIS-derived parameters

National estimates of spatially-resolved cropland net primary production (NPP) are needed for diagnostic and prognostic modeling of carbon sources, sinks, and net carbon flux between land and atmosphere. Cropland NPP estimates that correspond with existing cropland cover maps are needed to drive biogeochemical models at the local scale as well as national and continental scales. Existing satellite-based NPP products tend to underestimate NPP on croplands. An Agricultural Inventory-based Light Use Efficiency (AgI-LUE) framework was developed to estimate individual crop biophysical parameters for use in estimating crop-specific NPP over large multi-state regions. The method is documented here and evaluated for corn (Zea mays L.) and soybean (Glycine max L. Merr.) in Iowa and Illinois in 2006 and 2007. The method includes a crop-specific Enhanced Vegetation Index (EVI), shortwave radiation data estimated using the Mountain Climate Simulator (MTCLIM) algorithm, and crop-specific LUE per county. The combined aforementioned variables were used to generate spatially-resolved, crop-specific NPP that corresponds to the Cropland Data Layer (CDL) land cover product. Results from the modeling framework captured the spatial NPP gradient across croplands of Iowa and Illinois, and also represented the difference in NPP between years 2006 and 2007. Average corn and soybean NPP from AgI-LUE was 917 g C m⁻² yr⁻¹ and 409 g C m⁻² yr⁻¹, respectively. This was 2.4 and 1.1 times higher, respectively, for corn and soybean compared to the MOD17A3 NPP product. Site comparisons with flux tower data show AgI-LUE NPP in close agreement with tower-derived NPP, lower than inventory-based NPP, and higher than MOD17A3 NPP. The combination of new inputs and improved datasets enabled the development of spatially explicit and reliable NPP estimates for individual crops over large regional extents.     

Modeling BVOC isoprene emissions based on a GIS and remote sensing database

This paper presents a geographic information systems (GIS) model to relate biogenic volatile organic compounds (BVOCs) isoprene emissions to ecosystem type, as well as environmental drivers such as light intensity, temperature, landscape factor and foliar density. Data and techniques have recently become available which can permit new improved estimates of isoprene emissions over Hong Kong. The techniques are based on Guenther et al., 1993, Guenther et al., 1999 model. The spatially detailed mapping of isoprene emissions over Hong Kong at a resolution of 100 m and a database has been constructed for retrieval of the isoprene maps from February 2007 to January 2008. This approach assigns emission rates directly to ecosystem types not to individual species, since unlike in temperate regions where one or two single species may dominate over large regions, Hong Kong's vegetation is extremely diverse with up to 300 different species in 1 ha. Field measurements of emissions by canister sampling obtained a range of ambient emissions according to different climatic conditions for Hong Kong's main ecosystem types in both urban and rural areas, and these were used for model validation. Results show the model-derived isoprene flux to have high to moderate correlations with field observations (i.e. r² = 0.77, r² = 0.63, r² = 0.37 for all 24 field measurements, subset for summer, and winter data, respectively) which indicate the robustness of the approach when applied to tropical forests at detailed level, as well as the promising role of remote sensing in isoprene mapping. The GIS model and raster database provide a simple and low cost estimation of the BVOC isoprene in Hong Kong at detailed level. City planners and environmental authorities may use the derived models for estimating isoprene transportation, and its interaction with anthropogenic pollutants in urban areas.     

Texture augmented detection of macrophyte species using decision trees

Image classification using multispectral sensors has shown good performance in detecting macrophytes at the species level. However, species level classification often does not utilize the texture information provided by high resolution images. This study investigated whether image texture provides useful vector(s) for the discrimination of monospecific stands of three floating macrophyte species in Quickbird imagery of the South Nation River. Semivariograms indicated that window sizes of 5 × 5 and 13 × 13 pixels were the most appropriate spatial scales for calculation of the grey level co-occurrence matrix and subsequent texture attributes from the multispectral and panchromatic bands. Of the 214 investigated vectors (13 Haralick texture attributes * 15 bands + 9 spectral bands + 10 transformations/indices), feature selection determined which combination of spectral and textural vectors had the greatest class separability based on the Mann–Whitney U-test and Jefferies–Matusita distance. While multispectral red and near infrared (NIR) performed satisfactorily, the addition of panchromatic-dissimilarity slightly improved class separability and the accuracy of a decision tree classifier (Kappa: red/NIR/panchromatic-dissimilarity – 93.2% versus red/NIR – 90.4%). Class separability improved by incorporating a second texture attribute, but resulted in a decrease in classification accuracy. The results suggest that incorporating image texture may be beneficial for separating stands with high spatial heterogeneity. However, the benefits may be limited and must be weighed against the increased complexity of the classifier.     

Intercomparison and validation of MODIS and GLASS leaf area index (LAI) products over mountain areas: A case study in southwestern China

The validation study of leaf area index (LAI) products over rugged surfaces not only gives additional insights into data quality of LAI products, but deepens understanding of uncertainties regarding land surface process models depended on LAI data over complex terrain. This study evaluated the performance of MODIS and GLASS LAI products using the intercomparison and direct validation methods over southwestern China. The spatio-temporal consistencies, such as the spatial distributions of LAI products and their statistical relationship as a function of topographic indices, time, and vegetation types, respectively, were investigated through intercomparison between MODIS and GLASS products during the period 2011–2013. The accuracies and change ranges of these two products were evaluated against available LAI reference maps over 10 sampling regions which standed for typical vegetation types and topographic gradients in southwestern China.The results show that GLASS LAI exhibits higher percentage of good quality data (i.e. successful retrievals) and smoother temporal profiles than MODIS LAI. The percentage of successful retrievals for MODIS and GLASS is vulnerable to topographic indices, especially to relief amplitude. Besides, the two products do not capture seasonal dynamics of crop, especially in spring over heterogeneously hilly regions. The yearly mean LAI differences between MODIS and GLASS are within ±0.5 for 64.70% of the total retrieval pixels over southwestern China. The spatial distribution of mean differences and temporal profiles of these two products are inclined to be dominated by vegetation types other than topographic indices. The spatial and temporal consistency of these two products is good over most area of grasses/cereal crops; however, it is poor for evergreen broadleaf forest. MODIS presents more reliable change range of LAI than GLASS through comparison with fine resolution reference maps over most of sampling regions. The accuracies of direct validation are obtained for GLASS LAI (r = 0.35, RMSE = 1.72, mean bias = −0.71) and MODIS LAI (r = 0.49, RMSE = 1.75, mean bias = −0.67). GLASS performs similarly to MODIS, but may be marginally inferior to MODIS based on our direct validation results. The validation experience demonstrates the necessity and importance of topographic consideration for LAI estimation over mountain areas. Considerable attention will be paid to the improvements of surface reflectance, retrieval algorithm and land cover types so as to enhance the quality of LAI products in topographically complex terrain.     

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.     

Comparative analysis of different retrieval methods for mapping grassland leaf area index using airborne imaging spectroscopy

Fine scale maps of vegetation biophysical variables are useful status indicators for monitoring and managing national parks and endangered habitats. Here, we assess in a comparative way four different retrieval methods for estimating leaf area index (LAI) in grassland: two radiative transfer model (RTM) inversion methods (one based on look-up-tables (LUT) and one based on predictive equations) and two statistical modelling methods (one partly, the other entirely based on in situ data). For prediction, spectral data were used that had been acquired over Majella National Park in Italy by the airborne hyperspectral HyMap instrument. To assess the performance of the four investigated models, the normalized root mean squared error (nRMSE) and coefficient of determination (R²) between estimates and in situ LAI measurements are reported (n = 41). Using a jackknife approach, we also quantified the accuracy and robustness of empirical models as a function of the size of the available calibration data set. The results of the study demonstrate that the LUT-based RTM inversion yields higher accuracies for LAI estimation (R² = 0.91, nRMSE = 0.18) as compared to RTM inversions based on predictive equations (R² = 0.79, nRMSE = 0.38). The two statistical methods yield accuracies similar to the LUT method. However, as expected, the accuracy and robustness of the statistical models decrease when the size of the calibration database is reduced to fewer samples. The results of this study are of interest for the remote sensing community developing improved inversion schemes for spaceborne hyperspectral sensors applicable to different vegetation types. The examples provided in this paper may also serve as illustrations for the drawbacks and advantages of physical and empirical models.