Satellite-based automated burned area detection: A performance assessment of the MODIS MCD45A1 in the Brazilian savanna

Burnings, which cause major changes to the environment, can be effectively monitored via satellite data, regarding both the identification of active fires and the estimation of burned areas. Among the many orbital sensors suitable for mapping burned areas on global and regional scales, the moderate resolution imaging spectroradiometer (MODIS), on board the Terra and Aqua platforms, has been the most widely utilized. In this study, the performance of the MODIS MCD45A1 burned area product was thoroughly evaluated in the Brazilian savanna, the second largest biome in South America and a global biodiversity hotspot, characterized by a conspicuous climatic seasonality and the systematic occurrence of natural and anthropogenic fires. Overall, September MCD45A1 polygons (2000–2012) compared well to the Landsat-based reference mapping (r² = 0.92) and were closely accompanied, on a monthly basis, by MOD14 and MYD14 hotspots (r² = 0.89), although large omissions errors, linked to landscape patterns, structures, and overall conditions depicted in each reference image, were observed. In spite of its spatial and temporal limitations, the MCD45A1 product proved instrumental for mapping and understanding fire behavior and impacts on the Cerrado landscapes.     

Semi-automated mapping of burned areas in semi-arid ecosystems using MODIS time-series imagery

Understanding spatial and temporal patterns of burned areas at regional scales, provides a long-term perspective of fire processes and its effects on ecosystems and vegetation recovery patterns, and it is a key factor to design prevention and post-fire restoration plans and strategies. Remote sensing has become the most widely used tool to detect fire affected areas over large tracts of land (e.g., ecosystem, regional and global levels). Standard satellite burned area and active fire products derived from the 500-m Moderate Resolution Imaging Spectroradiometer (MODIS) and the Satellite Pour l’Observation de la Terre (SPOT) are available to this end. However, prior research caution on the use of these global-scale products for regional and sub-regional applications. Consequently, we propose a novel semi-automated algorithm for identification and mapping of burned areas at regional scale. The semi-arid Monte shrublands, a biome covering 240,000 km² in the western part of Argentina, and exposed to seasonal bushfires was selected as the test area. The algorithm uses a set of the normalized burned ratio index products derived from MODIS time series; using a two-phased cycle, it firstly detects potentially burned pixels while keeping a low commission error (false detection of burned areas), and subsequently labels them as seed patches. Region growing image segmentation algorithms are applied to the seed patches in the second-phase, to define the perimeter of fire affected areas while decreasing omission errors (missing real burned areas). Independently-derived Landsat ETM+ burned-area reference data was used for validation purposes. Additionally, the performance of the adaptive algorithm was assessed against standard global fire products derived from MODIS Aqua and Terra satellites, total burned area (MCD45A1), the active fire algorithm (MOD14); and the L3JRC SPOT VEGETATION 1 km GLOBCARBON products. The correlation between the size of burned areas detected by the global fire products and independently-derived Landsat reference data ranged from R² = 0.01–0.28, while our algorithm performed showed a stronger correlation coefficient (R² = 0.96). Our findings confirm prior research calling for caution when using the global fire products locally or regionally.     

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.     

Detecting forest damage after a low-severity fire using remote sensing at multiple scales

Remote sensing technologies are an ideal platform to examine the extent and impact of fire on the landscape. In this study we assess that capacity of the RapidEye constellation and Landsat (Thematic Mapper and Operational Land Imager to map fine-scale burn attributes for a small, low severity prescribed fire in a dry Western Canadian forest. Estimates of burn severity from field data were collated into a simple burn index and correlated with a selected suite of common spectral vegetation indices. Burn severity classes were then derived to map fire impacts and estimate consumed woody surface fuels (diameter ≥2.6 cm). All correlations between the simple burn index and vegetation indices produced significant results (p < 0.01), but varied substantially in their overall accuracy. Although the Landsat Soil Adjusted Vegetation Index provided the best regression fit (R² = 0.56), results suggested that RapidEye provided much more spatially detailed estimates of tree damage (Soil Adjusted Vegetation Index, R² = 0.51). Consumption estimates of woody surface fuels ranged from 3.38 ± 1.03 Mg ha⁻¹ to 11.73 ± 1.84 Mg ha⁻¹, across four derived severity classes with uncertainties likely a result of changing foliage moisture between the before and after fire images. While not containing spectral information in the short wave infrared, the spatial variability provided by the RapidEye imagery has potential for mapping and monitoring fine scale forest attributes, as well as the potential to resolve fire damage at the individual tree level.     

Aboveground total and green biomass of dryland shrub derived from terrestrial laser scanning

Sagebrush (Artemisia tridentata), a dominant shrub species in the sagebrush-steppe ecosystem of the western US, is declining from its historical distribution due to feedbacks between climate and land use change, fire, and invasive species. Quantifying aboveground biomass of sagebrush is important for assessing carbon storage and monitoring the presence and distribution of this rapidly changing dryland ecosystem. Models of shrub canopy volume, derived from terrestrial laser scanning (TLS) point clouds, were used to accurately estimate aboveground sagebrush biomass. Ninety-one sagebrush plants were scanned and sampled across three study sites in the Great Basin, USA. Half of the plants were scanned and destructively sampled in the spring (n = 46), while the other half were scanned again in the fall before destructive sampling (n = 45). The latter set of sagebrush plants was scanned during both spring and fall to further test the ability of the TLS to quantify seasonal changes in green biomass. Sagebrush biomass was estimated using both a voxel and a 3-D convex hull approach applied to TLS point cloud data. The 3-D convex hull model estimated total and green biomass more accurately (R² = 0.92 and R² = 0.83, respectively) than the voxel-based method (R² = 0.86 and R² = 0.73, respectively). Seasonal differences in TLS-predicted green biomass were detected at two of the sites (p < 0.001 and p = 0.029), elucidating the amount of ephemeral leaf loss in the face of summer drought. The methods presented herein are directly transferable to other dryland shrubs, and implementation of the convex hull model with similar sagebrush species is straightforward.     

An approach of surface coal fire detection from ASTER and Landsat-8 thermal data: Jharia coal field,India

Radiant temperature images from thermal remote sensing sensors are used to delineate surface coal fires, by deriving a cut-off temperature to separate coal-fire from non-fire pixels. Temperature contrast of coal fire and background elements (rocks and vegetation etc.) controls this cut-off temperature. This contrast varies across the coal field, as it is influenced by variability of associated rock types, proportion of vegetation cover and intensity of coal fires etc. We have delineated coal fires from background, based on separation in data clusters in maximum v/s mean radiant temperature (13th band of ASTER and 10th band of Landsat-8) scatter-plot, derived using randomly distributed homogeneous pixel-blocks (9 × 9 pixels for ASTER and 27 × 27 pixels for Landsat-8), covering the entire coal bearing geological formation. It is seen that, for both the datasets, overall temperature variability of background and fires can be addressed using this regional cut-off. However, the summer time ASTER data could not delineate fire pixels for one specific mine (Bhulanbararee) as opposed to the winter time Landsat-8 data. The contrast of radiant temperature of fire and background terrain elements, specific to this mine, is different from the regional contrast of fire and background, during summer. This is due to the higher solar heating of background rocky outcrops, thus, reducing their temperature contrast with fire. The specific cut-off temperature determined for this mine, to extract this fire, differs from the regional cut-off. This is derived by reducing the pixel-block size of the temperature data. It is seen that, summer-time ASTER image is useful for fire detection but required additional processing to determine a local threshold, along with the regional threshold to capture all the fires. However, the winter Landsat-8 data was better for fire detection with a regional threshold.     

Predictive fire occurrence modelling to improve burned area estimation at a regional scale: A case study in East Caprivi,Namibia

Fires threaten human lives, property and natural resources in Southern African savannas. Due to warming climate, fire occurrence may increase and fires become more intense. It is crucial, therefore, to understand the complexity of spatiotemporal and probabilistic characteristics of fires. This study scrutinizes spatiotemporal characteristics of fires and the role played by abiotic, biotic and anthropogenic factors for fire probability modelling in a semiarid Southern African savanna environment. The MODIS fire products: fire hot spots (MOD14A2 and MYD14A2) and burned area product MODIS (MCD45A1), and GIS derived data were used in analysis. Fire hot spots occurrence was first analysed, and spatial autocorrelation for fires investigated, using Moran's I correlograms. Fire probability models were created using generalized linear models (GLMs). Separate models were produced for abiotic, biotic, anthropogenic and combined factors and an autocovariate variable was tested for model improvement. The hierarchical partitioning method was used to determine independent effects of explanatory variables. The discriminating ability of models was evaluated using area under the curve (AUC) from the receiver operating characteristic (ROC) plot. The results showed that 19.2–24.4% of East Caprivi burned when detected using MODIS hot spots fire data and these fires were strongly spatially autocorrelated. Therefore, the autocovariate variable significantly improved fire probability models when added to them. For autologistic models, i.e. models accounting for spatial autocorrelation, discrimination was good to excellent (AUC 0.858–0.942). For models not counting spatial autocorrelation, prediction success was poor to moderate (AUC 0.542–0.745). The results of this study clearly showed that spatial autocorrelation has to be taken in to account in the fire probability model building process when using remotely sensed and GIS derived data. This study also showed that fire probability models accounting for spatial autocorrelation proved to be superior in regional scale burned area estimation when compared with MODIS burned area product (MCD45A1).     

Influences of forest roads and their edge effects on the spatial pattern of burn severity

Previous research has shown that forest roads are an important feature in many landscapes and have significant effects on wildfire ignition and cessation. However, forest road effects on burn severity have not been studied at the landscape level. Therefore, the overarching goal of our study is to identify the influences of road edge effects on the spatial patterns of burn severity. We analyzed six fires within the Okanogan–Wenatchee National Forest on the eastern slope of the Cascades mountain range of central Washington.We generated two categories for assessing road variables: (1) Primary Road Effect Zone (area within 150 m of the nearest road) and (2) Secondary Road Effect Zone (area from 150 m to 300 m to the nearest road). A regular sampling grid including one out of every 9 cells was created for each fire.These grids were intersected with burn severity data in the form of the Relative Differenced Normalized Burn Ratio (RdNBR), road distance category, stream distance, elevation, slope, terrain shape index, heat load index, canopy cover, and fuel type. We fit spatial regression models with RdNBR as the dependent variable.We found that high burn severity is less likely to occur in the Primary Road Effect Zone for most fires, although one fire exhibited the opposite relationship. Forest road edge effects were hypothesized to be an important determinant of burn severity because fragmentation created by roads alters the roadside fuel profile and environment and because road corridors create barriers to fire spread. Recognizing roadside effects on burn severity patterns highlights the need for further study of the range of effects that roads have on fuels and the fire environment and the potential for incorporating road effects into landscape-level assessments of fire risk.     

Improving the prediction of arsenic contents in agricultural soils by combining the reflectance spectroscopy of soils and rice plants

Visible and near-infrared reflectance spectroscopy provides a beneficial tool for investigating soil heavy metal contamination. This study aimed to investigate mechanisms of soil arsenic prediction using laboratory based soil and leaf spectra, compare the prediction of arsenic content using soil spectra with that using rice plant spectra, and determine whether the combination of both could improve the prediction of soil arsenic content. A total of 100 samples were collected and the reflectance spectra of soils and rice plants were measured using a FieldSpec3 portable spectroradiometer (350–2500 nm). After eliminating spectral outliers, the reflectance spectra were divided into calibration (n = 62) and validation (n = 32) data sets using the Kennard-Stone algorithm. Genetic algorithm (GA) was used to select useful spectral variables for soil arsenic prediction. Thereafter, the GA-selected spectral variables of the soil and leaf spectra were individually and jointly employed to calibrate the partial least squares regression (PLSR) models using the calibration data set. The regression models were validated and compared using independent validation data set. Furthermore, the correlation coefficients of soil arsenic against soil organic matter, leaf arsenic and leaf chlorophyll were calculated, and the important wavelengths for PLSR modeling were extracted. Results showed that arsenic prediction using the leaf spectra (coefficient of determination in validation, R_v² = 0.54; root mean square error in validation, RMSE_v = 12.99 mg kg⁻¹; and residual prediction deviation in validation, RPD_v = 1.35) was slightly better than using the soil spectra (R_v² = 0.42, RMSE_v = 13.35 mg kg⁻¹, and RPD_v = 1.31). However, results also showed that the combinational use of soil and leaf spectra resulted in higher arsenic prediction (R_v² = 0.63, RMSE_v = 11.94 mg kg⁻¹, RPD_v = 1.47) compared with either soil or leaf spectra alone. Soil spectral bands near 480, 600, 670, 810, 1980, 2050 and 2290 nm, leaf spectral bands near 700, 890 and 900 nm in PLSR models were important wavelengths for soil arsenic prediction. Moreover, soil arsenic showed significantly positive correlations with soil organic matter (r = 0.62, p < 0.01) and leaf arsenic (r = 0.77, p < 0.01), and a significantly negative correlation with leaf chlorophyll (r = −0.67, p < 0.01). The results showed that the prediction of arsenic contents using soil and leaf spectra may be based on their relationships with soil organic matter and leaf chlorophyll contents, respectively. Although RPD of 1.47 was below the recommended RPD of >2 for soil analysis, arsenic prediction in agricultural soils can be improved by combining the leaf and soil spectra.     

Dry-season vegetation mass and cover fraction from SWIR1.6 and SWIR2.1 band ratio: Ground-radiometer and MODIS data in the Sahel

The potential of the short-wave infrared (SWIR) bands to detect dry-season vegetation mass and cover fraction is investigated with ground radiometry and MODIS data, confronted to vegetation data collected in rangeland and cropland sites in the Sahel (Senegal, Niger, Mali). The ratio of the 1.6 and 2.1 μm bands (called STI) acquired with a ground radiometer proved well suited for grassland mass estimation up to 2500 kg/ha with a linear relation (r² = 0.89). A curvilinear regression is accurate for masses ranging up to 3500 kg/ha. STI proved also well suited to retrieve vegetation cover fraction in crop fields, fallows and rangelands. Such dry-season monitoring, with either ground or satellite data, has important applications for forage, erosion risk and fire risk assessment in semi-arid areas.     

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.     

An integrated remote sensing and GIS approach for monitoring areas affected by selective logging: A case study in northern Mato Grosso,Brazilian Amazon

Forest cover disturbances due to processes such as logging and forest fires are a widespread issue especially in the tropics, and have heavily affected forest biomass and functioning in the Brazilian Amazon in the past decades. Satellite remote sensing has played a key role for assessing logging activities in this region; however, there are still remaining challenges regarding the quantification and monitoring of these processes affecting forested lands. In this study, we propose a new method for monitoring areas affected by selective logging in one of the hotspots of Mato Grosso state in the Brazilian Amazon, based on a combination of object-based and pixel-based classification approaches applied on remote sensing data. Logging intensity and changes over time are assessed within grid cells of 300 m × 300 m spatial resolution. Our method encompassed three main steps: (1) mapping forest/non-forest areas through an object-based classification approach applied to a temporal series of Landsat images during the period 2000–2015, (2) mapping yearly logging activities from soil fraction images on the same Landsat data series, and (3) integrating information from previous steps within a regular grid-cell of 300 m × 300 m in order to monitor disturbance intensities over this 15-years period. The overall accuracy of the baseline forest/non-forest mask (year 2000) and of the undisturbed vs disturbed forest (for selected years) were 93% and 84% respectively. Our results indicate that annual forest disturbance rates, mainly due to logging activities, were higher than annual deforestation rates during the whole period of study. The deforested areas correspond to circa 25% of the areas affected by forest disturbances. Deforestation rates were highest from 2001 to 2005 and then decreased considerably after 2006. In contrast, the annual forest disturbance rates show high temporal variability with a slow decrease over the 15-year period, resulting in a significant increase of the ratio between disturbed and deforested areas. Although the majority of the areas, which have been affected by selective logging during the period 2000–2014, were not deforested by 2015, more than 70% of the deforested areas in 2015 had been at least once identified as disturbed forest during that period.     

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.     

Non-linear partial least square regression increases the estimation accuracy of grass nitrogen and phosphorus using in situ hyperspectral and environmental data

Grass nitrogen (N) and phosphorus (P) concentrations are direct indicators of rangeland quality and provide imperative information for sound management of wildlife and livestock. It is challenging to estimate grass N and P concentrations using remote sensing in the savanna ecosystems. These areas are diverse and heterogeneous in soil and plant moisture, soil nutrients, grazing pressures, and human activities. The objective of the study is to test the performance of non-linear partial least squares regression (PLSR) for predicting grass N and P concentrations through integrating in situ hyperspectral remote sensing and environmental variables (climatic, edaphic and topographic). Data were collected along a land use gradient in the greater Kruger National Park region. The data consisted of: (i) in situ-measured hyperspectral spectra, (ii) environmental variables and measured grass N and P concentrations. The hyperspectral variables included published starch, N and protein spectral absorption features, red edge position, narrow-band indices such as simple ratio (SR) and normalized difference vegetation index (NDVI). The results of the non-linear PLSR were compared to those of conventional linear PLSR. Using non-linear PLSR, integrating in situ hyperspectral and environmental variables yielded the highest grass N and P estimation accuracy (R² = 0.81, root mean square error (RMSE) = 0.08, and R² = 0.80, RMSE = 0.03, respectively) as compared to using remote sensing variables only, and conventional PLSR. The study demonstrates the importance of an integrated modeling approach for estimating grass quality which is a crucial effort towards effective management and planning of protected and communal savanna ecosystems.     

Characterization of aboveground biomass in an unmanaged boreal forest using Landsat temporal segmentation metrics

Canada is dominated by forested ecosystems which are subject to various inventory and management practices, with more northern boreal forests subject to neither. Our objectives were to measure the capacity of temporal trajectory metrics for estimating selected forest attributes in a northern Canadian boreal forest context using Landsat imagery and investigate the importance of different types of temporal trajectory metrics. Results indicated that Wetness was the best Tasseled Cap (TC) component for aboveground biomass estimation (R² = 50%, RMSE% = 56%), and the combination of simple and complex metrics from all TC components produced the highest R² (62%) and lowest RMSE% (49%). Using a similar combination of variables, other forest attributes were estimated equally reliably with lower RMSE% values. The most important temporal trajectory metrics were simple and described TC component values at each point of change in the temporal trajectory, however the most important variables overall were environmental variables.     

Sensitivity of the MODIS fire detection algorithm (MOD14) in the savanna region of the Northern Territory,Australia

Fire detection using satellites is an important source of information for fire management, ecological studies and emission estimates. However, little is known about the minimum sizes of fires that are being detected. This paper presents an approach using fire radiative power estimated from MODIS satellite data to determine the detection threshold for fire-prone savannas in Northern Australia. The results indicate that fires with an active flaming area 100–300 m² can be detected in the study region. It is also shown that the algorithm is slightly more sensitive at night. As expected the detection threshold shows strong view angle dependence. While this study has been undertaken in the savannas of Northern Australia, the results should be transferable to other savanna regions worldwide and other areas where fires are not obscured by a dense tree canopy.     

A new 500-m resolution map of canopy height for Amazon forest using spaceborne LiDAR and cloud-free MODIS imagery

In the present study, we aimed to map canopy heights in the Brazilian Amazon mainly on the basis of spaceborne LiDAR and cloud-free MODIS imagery with a new method (the Self-Organizing Relationships method) for spatial modeling of the LiDAR footprint. To evaluate the general versatility, we compared the created canopy height map with two different canopy height estimates on the basis of our original field study plots (799 plots located in eight study sites) and a previously developed canopy height map. The compared canopy height estimates were obtained by: (1) a stem diameter at breast height (D) – tree height (H) relationship specific to each site on the basis of our original field study, (2) a previously developed D–H model involving environmental and structural factors as explanatory variables (Feldpausch et al., 2011), and (3) a previously developed canopy height map derived from the spaceborne LiDAR data with different spatial modeling method and explanatory variables (Simard et al., 2011). As a result, our canopy height map successfully detected a spatial distribution pattern in canopy height estimates based on our original field study data (r = 0.845, p = 8.31 × 10⁻³) though our canopy height map showed a poor correlation (r = 0.563, p = 0.146) with the canopy height estimate based on a previously developed model by Feldpausch et al. (2011). We also confirmed that the created canopy height map showed a similar pattern with the previously developed canopy height map by Simard et al. (2011). It was concluded that the use of the spaceborne LiDAR data provides a sufficient accuracy in estimating the canopy height at regional scale.     

A practical algorithm for estimating surface soil moisture using combined optical and thermal infrared data

Surface soil moisture (SSM) is a critical variable for understanding the energy and water exchange between the land and atmosphere. A multi-linear model was recently developed to determine SSM using ellipse variables, namely, the center horizontal coordinate (x₀), center vertical coordinate (y₀), semi-major axis (a) and rotation angle (θ), derived from the elliptical relationship between diurnal cycles of land surface temperature (LST) and net surface shortwave radiation (NSSR). However, the multi-linear model has a major disadvantage. The model coefficients are calculated based on simulated data produced by a land surface model simulation that requires sufficient meteorological measurements. This study aims to determine the model coefficients directly using limited meteorological parameters rather than via the complicated simulation process, decreasing the dependence of the model coefficients on meteorological measurements. With the simulated data, a practical algorithm was developed to estimate SSM based on combined optical and thermal infrared data. The results suggest that the proposed approach can be used to determine the coefficients associated with all ellipse variables based on historical meteorological records, whereas the constant term varies daily and can only be determined using the daily maximum solar radiation in a prediction model. Simulated results from three FLUXNET sites over 30 cloud-free days revealed an average root mean square error (RMSE) of 0.042 m³/m³ when historical meteorological records were used to synchronously determine the model coefficients. In addition, estimated SSM values exhibited generally moderate accuracies (coefficient of determination R² = 0.395, RMSE = 0.061 m³/m³) compared to SSM measurements at the Yucheng Comprehensive Experimental Station.     

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.     

A new calibration of the effective scattering albedo and soil roughness parameters in the SMOS SM retrieval algorithm

This study focuses on the calibration of the effective vegetation scattering albedo (ω) and surface soil roughness parameters (H_R, and N_Rp, p = H,V) in the Soil Moisture (SM) retrieval from L-band passive microwave observations using the L-band Microwave Emission of the Biosphere (L-MEB) model. In the current Soil Moisture and Ocean Salinity (SMOS) Level 2 (L2), v620, and Level 3 (L3), v300, SM retrieval algorithms, low vegetated areas are parameterized by ω = 0 and H_R = 0.1, whereas values of ω = 0.06 − 0.08 and H_R = 0.3 are used for forests. Several parameterizations of the vegetation and soil roughness parameters (ω, H_R and N_Rp, p = H,V) were tested in this study, treating SMOS SM retrievals as homogeneous over each pixel instead of retrieving SM over a representative fraction of the pixel, as implemented in the operational SMOS L2 and L3 algorithms. Globally-constant values of ω = 0.10, H_R = 0.4 and N_Rp = −1 (p = H,V) were found to yield SM retrievals that compared best with in situ SM data measured at many sites worldwide from the International Soil Moisture Network (ISMN). The calibration was repeated for collections of in situ sites classified in different land cover categories based on the International Geosphere-Biosphere Programme (IGBP) scheme. Depending on the IGBP land cover class, values of ω and H_R varied, respectively, in the range 0.08–0.12 and 0.1–0.5. A validation exercise based on in situ measurements confirmed that using either a global or an IGBP-based calibration, there was an improvement in the accuracy of the SM retrievals compared to the SMOS L3 SM product considering all statistical metrics (R = 0.61, bias = −0.019 m³ m⁻³, ubRMSE = 0.062 m³ m⁻³ for the IGBP-based calibration; against R = 0.54, bias = −0.034 m³ m⁻³ and ubRMSE = 0.070 m³ m⁻³ for the SMOS L3 SM product). This result is a key step in the calibration of the roughness and vegetation parameters in the operational SMOS retrieval algorithm. The approach presented here is the core of a new forthcoming SMOS optimized SM product.