Mapping and estimating land change between 2001 and 2013 in a heterogeneous landscape in West Africa: Loss of forestlands and capacity building opportunities

In West Africa, accurate classification of land cover and land change remains a big challenge due to the patchy and heterogeneous nature of the landscape. Limited data availability, human resources and technical capacities, further exacerbate the challenge. The result is a region that is among the more understudied areas in the world, which in turn has resulted in a lack of appropriate information required for sustainable natural resources management. The objective of this paper is to explore open source software and easy-to-implement approaches to mapping and estimation of land change that are transferrable to local institutions to increase capacity in the region, and to provide updated information on the regional land surface dynamics. To achieve these objectives, stable land cover and land change between 2001 and 2013 in the Kara River Basin in Togo and Benin were mapped by direct multitemporal classification of Landsat data by parameterization and evaluation of two machine-learning algorithms. Areas of land cover and change were estimated by application of an unbiased estimator to sample data following international guidelines. A prerequisite for all tools and methods was implementation in an open source environment, and adherence to international guidelines for reporting land surface activities. Findings include a recommendation of the Random Forests algorithm as implemented in Orfeo Toolbox, and a stratified estimation protocol − all executed in the QGIS graphical use interface. It was found that despite an estimated reforestation of 10,0727 ± 3480 ha (95% confidence interval), the combined rate of forest and savannah loss amounted to 56,271 ± 9405 ha (representing a 16% loss of the forestlands present in 2001), resulting in a rather sharp net loss of forestlands in the study area. These dynamics had not been estimated prior to this study, and the results will provide useful information for decision making pertaining to natural resources management, land management planning, and the implementation of the United Nations Collaborative Programme on Reducing Emissions from Deforestation and Forest Degradation in Developing Countries (UN-REDD).     

Spatial analysis of human-induced vegetation productivity decline over eastern Africa using a decade (2001–2011) of medium resolution MODIS time-series data

Climate variation and land transformations related to exploitative land uses are among the main drivers of vegetation productivity decline and ongoing land degradation in East Africa. We combined analysis of vegetation trends and cumulative rain use efficiency differences (CRD), calculated from 250-m MODIS NDVI time-series data, to map vegetation productivity loss over eastern Africa between 2001 and 2011. The CRD index values were furthermore used to discern areas of particular severe vegetation productivity loss over the observation period. Monthly 25-km Tropical Rainfall Measuring Mission (TRMM) data metrics were used to mask areas of rainfall declines not related to human-induced land productivity loss. To provide insights on the productivity decline, we linked the MODIS-based vegetation productivity map to land transformation processes using very high resolution (VHR) imagery in Google Earth (GE) and a Landsat-based land-cover change map. In total, 3.8 million ha experienced significant vegetation loss over the monitoring period. An overall agreement of 68% was found between the rainfall-corrected MODIS productivity decline map and all reference pixels discernable from GE and the Landsat map. The CRD index showed a good potential to discern areas with ‘severe’ vegetation productivity losses under high land-use intensities.     

Satellite monitoring of urbanization and environmental impacts—A comparison of Stockholm and Shanghai

This study investigates urbanization and its potential environmental consequences in Shanghai and Stockholm metropolitan areas over two decades. Changes in land use/land cover are estimated from support vector machine classifications of Landsat mosaics with grey-level co-occurrence matrix features. Landscape metrics are used to investigate changes in landscape composition and configuration and to draw preliminary conclusions about environmental impacts. Speed and magnitude of urbanization is calculated by urbanization indices and the resulting impacts on the environment are quantified by ecosystem services. Growth of urban areas and urban green spaces occurred at the expense of cropland in both regions. Alongside a decrease in natural land cover, urban areas increased by approximately 120% in Shanghai, nearly ten times as much as in Stockholm, where the most significant land cover change was a 12% urban expansion that mostly replaced agricultural areas. From the landscape metrics results, it appears that fragmentation in both study regions occurred mainly due to the growth of high density built-up areas in previously more natural/agricultural environments, while the expansion of low density built-up areas was for the most part in conjunction with pre-existing patches. Urban growth resulted in ecosystem service value losses of approximately 445 million US dollars in Shanghai, mostly due to the decrease in natural coastal wetlands while in Stockholm the value of ecosystem services changed very little. Total urban growth in Shanghai was 1768 km² and 100 km² in Stockholm. The developed methodology is considered a straight-forward low-cost globally applicable approach to quantitatively and qualitatively evaluate urban growth patterns that could help to address spatial, economic and ecological questions in urban and regional planning.     

Identifying natural and anthropogenically-induced geohazards from satellite ground motion and geospatial data: Stoke-on-Trent,UK

Determining the location and nature of hazardous ground motion resulting from natural and anthropogenic processes such as landslides, tectonic movement and mining is essential for hazard mitigation and sustainable resource use. Ground motion estimates from satellite ERS-1/2 persistent scatterer interferometry (PSI) were combined with geospatial data to identify areas of observed geohazards in Stoke-on-Trent, UK. This investigation was performed within the framework of the EC FP7-SPACE PanGeo project which aimed to provide free and open access to geohazard information for 52 urban areas across Europe. Geohazards identified within the city of Stoke-on-Trent and neighbouring rural areas are presented here alongside an examination of the PanGeo methodology.A total of 14 areas experiencing ground instability caused by natural and anthropogenic processes have been defined, covering 122.35 km². These are attributed to a range of geohazards, including landslides, ground dissolution, made ground and mining activities. The dominant geohazard (by area) is ground movement caused by post-mining groundwater recharge and mining-related subsidence (93.19% of total geohazard area), followed by landsliding (5.81%). Observed ground motions along the satellite line-of-sight reach maxima of +35.23 mm/yr and −22.57 mm/yr. A combination of uplift, subsidence and downslope movement is displayed.‘Construction sites’ and ‘continuous urban fabric’ (European Urban Atlas land use types) form the land uses most affected (by area) by ground motion and ‘discontinuous very low density urban fabric’ the least. Areas of ‘continuous urban fabric’ also show the highest average velocity towards the satellite (5.08 mm/yr) and the highest PS densities (1262.92 points/km²) along with one of the lowest standard deviations. Rural land uses tend to result in lower PS densities and higher standard deviations, a consequence of fewer suitable reflectors in these regions. PSI is also limited in its ability to identify especially rapid ground motion. As a consequence the supporting geospatial data proved especially useful for the identification of landslides and some areas of ground dissolution. The mapped areas of instability are also compared with modelled potential geohazards (the BGS GeoSure dataset).     

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.     

Selection of LiDAR geometric features with adaptive neighborhood size for urban land cover classification

LiDAR has been an effective technology for acquiring urban land cover data in recent decades. Previous studies indicate that geometric features have a strong impact on land cover classification. Here, we analyzed an urban LiDAR dataset to explore the optimal feature subset from 25 geometric features incorporating 25 scales under 6 definitions for urban land cover classification. We performed a feature selection strategy to remove irrelevant or redundant features based on the correlation coefficient between features and classification accuracy of each features. The neighborhood scales were divided into small (0.5–1.5 m), medium (1.5–6 m) and large (>6 m) scale. Combining features with lower correlation coefficient and better classification performance would improve classification accuracy. The feature depicting homogeneity or heterogeneity of points would be calculated at a small scale, and the features to smooth points at a medium scale and the features of height different at large scale. As to the neighborhood definition, cuboid and cylinder were recommended. This study can guide the selection of optimal geometric features with adaptive neighborhood scale for urban land cover classification.     

Mapping the distribution of the main host for plague in a complex landscape in Kazakhstan: An object-based approach using SPOT-5 XS,Landsat 7 ETM+, SRTM and multiple Random Forests

Plague is a zoonotic infectious disease present in great gerbil populations in Kazakhstan. Infectious disease dynamics are influenced by the spatial distribution of the carriers (hosts) of the disease. The great gerbil, the main host in our study area, lives in burrows, which can be recognized on high resolution satellite imagery. In this study, using earth observation data at various spatial scales, we map the spatial distribution of burrows in a semi-desert landscape.The study area consists of various landscape types. To evaluate whether identification of burrows by classification is possible in these landscape types, the study area was subdivided into eight landscape units, on the basis of Landsat 7 ETM+ derived Tasselled Cap Greenness and Brightness, and SRTM derived standard deviation in elevation.In the field, 904 burrows were mapped. Using two segmented 2.5 m resolution SPOT-5 XS satellite scenes, reference object sets were created. Random Forests were built for both SPOT scenes and used to classify the images. Additionally, a stratified classification was carried out, by building separate Random Forests per landscape unit.Burrows were successfully classified in all landscape units. In the ‘steppe on floodplain’ areas, classification worked best: producer's and user's accuracy in those areas reached 88% and 100%, respectively. In the ‘floodplain’ areas with a more heterogeneous vegetation cover, classification worked least well; there, accuracies were 86 and 58% respectively. Stratified classification improved the results in all landscape units where comparison was possible (four), increasing kappa coefficients by 13, 10, 9 and 1%, respectively.In this study, an innovative stratification method using high- and medium resolution imagery was applied in order to map host distribution on a large spatial scale. The burrow maps we developed will help to detect changes in the distribution of great gerbil populations and, moreover, serve as a unique empirical data set which can be used as input for epidemiological plague models. This is an important step in understanding the dynamics of plague.     

Quantifying determinants of cash crop expansion and their relative effects using logistic regression modeling and variance partitioning

Cash crop expansion has been a major land use change in tropical and subtropical regions worldwide. Quantifying the determinants of cash crop expansion should provide deeper spatial insights into the dynamics and ecological consequences of cash crop expansion. This paper investigated the process of cash crop expansion in Hangzhou region (China) from 1985 to 2009 using remotely sensed data. The corresponding determinants (neighborhood, physical, and proximity) and their relative effects during three periods (1985–1994, 1994–2003, and 2003–2009) were quantified by logistic regression modeling and variance partitioning. Results showed that the total area of cash crops increased from 58,874.1 ha in 1985 to 90,375.1 ha in 2009, with a net growth of 53.5%. Cash crops were more likely to grow in loam soils. Steep areas with higher elevation would experience less likelihood of cash crop expansion. A consistently higher probability of cash crop expansion was found on places with abundant farmland and forest cover in the three periods. Besides, distance to river and lake, distance to county center, and distance to provincial road were decisive determinants for farmers’ choice of cash crop plantation. Different categories of determinants and their combinations exerted different influences on cash crop expansion. The joint effects of neighborhood and proximity determinants were the strongest, and the unique effect of physical determinants decreased with time. Our study contributed to understanding of the proximate drivers of cash crop expansion in subtropical regions.     

Monitoring vegetation dynamics with medium resolution MODIS-EVI time series at sub-regional scale in southern Africa

Currently there is a lack of knowledge on spatio-temporal patterns of land surface dynamics at medium spatial scale in southern Africa, even though this information is essential for better understanding of ecosystem response to climatic variability and human-induced land transformations. In this study, we analysed vegetation dynamics across a large area in southern Africa using the 14-years (2000–2013) of medium spatial resolution (250 m) MODIS-EVI time-series data. Specifically, we investigated temporal changes in the time series of key phenometrics including overall greenness, peak and timing of annual greenness over the monitoring period and study region. In order to specifically capture spatial and per pixel vegetation changes over time, we calculated trends in these phenometrics using a robust trend analysis method. The results showed that interannual vegetation dynamics followed precipitation patterns with clearly differentiated seasonality. The earliest peak greenness during 2000–2013 occurred at the end of January in the year 2000 and the latest peak greenness was observed at the mid of March in 2012. Specifically spatial patterns of long-term vegetation trends allowed mapping areas of (i) decrease or increase in overall greenness, (ii) decrease or increase of peak greenness, and (iii) shifts in timing of occurrence of peak greenness over the 14-year monitoring period. The observed vegetation decline in the study area was mainly attributed to human-induced factors. The obtained information is useful to guide selection of field sites for detailed vegetation studies and land rehabilitation interventions and serve as an input for a range of land surface models.     

Land subsidence in major cities of Central Mexico: Interpreting InSAR-derived land subsidence mapping with hydrogeological data

Significant structural damages to urban infrastructures caused by compaction of over-exploited aquifers are an important problem in Central Mexico. While the case of Mexico City has been well-documented, insight into land subsidence problems in other cities of Central Mexico is still limited. Among the cities concerned, we present and discuss the cases of five of them, located within the Trans-Mexican Volcanic Belt (TMVB): Toluca, Celaya, Aguascalientes, Morelia, and Queretaro. Applying the SBAS-InSAR method to C-Band RADARSAT-2 data, five high resolution ground motion time-series were produced to monitor the spatio-temporal variations of displacements and fracturing from 2012 to 2014. The study presents recent changes of land subsidence rates along with concordant geological and water data. It aims to provide suggestions to mitigate future damages to infrastructure and to assist in groundwater resources management.Aguascalientes, Celaya, Morelia and Queretaro (respectively in order of decreasing subsidence rates) are typical cases of fault-limited land subsidence of Central Mexico. It occurs as a result of groundwater over-exploitation in lacustrine and alluvial deposits covering highly variable bedrock topography, typical of horst-graben geological settings. Aguascalientes and Toluca show high rates of land subsidence (up to 10 cm/yr), while Celaya and Morelia show lower rates (from 2 to 5 cm/yr). Comparing these results with previous studies, it is inferred that the spatial patterns of land subsidence have changed in the city of Toluca. This change appears to be mainly controlled by the spatial heterogeneity of compressible sediments since no noticeable change occurred in groundwater extraction and related drawdown rates. While land subsidence of up to 8 cm/yr has been reported in the Queretaro Valley before 2011, rates inferior to 1 cm/yr are measured in 2013–2014. The subsidence has been almost entirely mitigated by major changes in the water management practices of the city, i.e., a 122 km long pipeline bringing surface water from an adjacent state allowed to cease pumping in half of the wells.     

Mapping land cover gradients through analysis of hyper-temporal NDVI imagery

The green cover of the earth exhibits various spatial gradients that represent gradual changes in space of vegetation density and/or in species composition. To date, land cover mapping methods differentiate at best, mapping units with different cover densities and/or species compositions, but typically fail to express such differences as gradients. Present interpretation techniques still make insufficient use of freely available spatial-temporal Earth Observation (EO) data that allow detection of existing land cover gradients. This study explores the use of hyper-temporal NDVI imagery to detect and delineate land cover gradients analyzing the temporal behavior of NDVI values. MODIS-Terra MVC-images (250 m, 16-day) of Crete, Greece, from February 2000 to July 2009 are used. The analysis approach uses an ISODATA unsupervised classification in combination with a Hierarchical Clustering Analysis (HCA). Clustering of class-specific temporal NDVI profiles through HCA resulted in the identification of gradients in landcover vegetation growth patterns. The detected gradients were arranged in a relational diagram, and mapped. Three groups of NDVI-classes were evaluated by correlating their class-specific annual average NDVI values with the field data (tree, shrub, grass, bare soil, stone, litter fraction covers). Multiple regression analysis showed that within each NDVI group, the fraction cover data were linearly related with the NDVI data, while NDVI groups were significantly different with respect to tree cover (adj. R² = 0.96), shrub cover (adj. R² = 0.83), grass cover (adj. R² = 0.71), bare soil (adj. R² = 0.88), stone cover (adj. R² = 0.83) and litter cover (adj. R² = 0.69) fractions. Similarly, the mean Sorenson dissimilarity values were found high and significant at confidence interval of 95% in all pairs of three NDVI groups. The study demonstrates that hyper-temporal NDVI imagery can successfully detect and map land cover gradients. The results may improve land cover assessment and aid in agricultural and ecological studies.     

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.     

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.     

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.     

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.     

An easy and low-cost method for preprocessing and matching small-scale amateur aerial photography for assessing agricultural land use in Burkina Faso

In recent decades, the Kou watershed in south-western Burkina Faso has suffered from poor water management. Despite the abundance of water, most water users regularly face water shortages because of the increase in the amount of land under irrigation. To help them achieve a more equitable allocation of irrigated land, local stakeholders need an easily managed low-cost tool for monitoring and mapping these irrigated zones. The aim of this study was to develop a fast and low-cost procedure for mosaicing and geo referencing amateur small-scale aerial photographs for land-use surveys. Sets of tens (2009) and hundreds (2007) of low-altitude aerial photographs, with a resolution of 0.4 m and 0.8 m, respectively, were used to create a detailed land-cover map of typical African small-scale irrigated agriculture. A commercially available stitching tool and GIS allowed geo referenced 'mono-images’ to be constructed; both mosaics were warped on a high-resolution SPOT image with a horizontal root mean square error (RMSE) of about 11 m. The RMSE between the two image datasets was 2 m. This approach is less sensitive to atmospheric conditions that are non-predictable in programming satellite imagery.     

Developing a temporally land cover-based look-up table (TL-LUT) method for estimating land surface temperature based on AMSR-E data over the Chinese landmass

The land surface temperature (LST) is an important parameter when studying the interface between the atmosphere and the Earth's surface. Compared to satellite thermal infrared (TIR) remote sensing, passive microwave (PMW) remote sensing is better able to overcome atmospheric influences and to estimate the LST, especially in cloudy regions. However, methods for estimating PMW LSTs at the country and continental scales are still rare. The necessity of training such methods from a temporally dynamic perspective also needs further investigations. Here, a temporally land cover based look-up table (TL-LUT) method is proposed to estimate the LSTs from AMSR-E data over the Chinese landmass. In this method, the synergies between observations from MODIS (Moderate Resolution Imaging Spectroradiometer) and AMSR-E (Advanced Microwave Scanning Radiometer for EOS), which are onboard the same Aqua satellite, are explored. Validation with the synchronous MODIS LSTs demonstrates that the TL-LUT method has better performances in retrieving LSTs with AMSR-E data than the method that uses a single brightness temperature in 36.5 GHz vertical polarization channel. The accuracy of the TL-LUT method is better than 2.7 K for forest and 3.2 K for cropland. Its accuracy varies according to land cover type, time of day, and season. When compared with the in-situ measured LSTs at four sites without urban warming in the Tibet Plateau, the standard errors of estimation between the estimated AMSR-E LST and in-situ measured LST are from 5.1 K to 6.0 K in the daytime and 3.1 K to 4.5 K in the nighttime. Further comparison with the in-situ measured air temperatures at 24 meteorological stations confirms the good performance of the TL-LUT method. The feasibility of PMW remote sensing in estimating the LST for China can complement the TIR data and can, therefore, aid in the generation of daily LST maps for the entire country. Further study of the penetration of PMW radiation would benefit the LST estimations in barren and other sparsely vegetated environments.     

Object-based land-cover classification for metropolitan Phoenix,Arizona, using aerial photography

Detailed land-cover mapping is essential for a range of research issues addressed by the sustainability and land system sciences and planning. This study uses an object-based approach to create a 1 m land-cover classification map of the expansive Phoenix metropolitan area through the use of high spatial resolution aerial photography from National Agricultural Imagery Program. It employs an expert knowledge decision rule set and incorporates the cadastral GIS vector layer as auxiliary data. The classification rule was established on a hierarchical image object network, and the properties of parcels in the vector layer were used to establish land cover types. Image segmentations were initially utilized to separate the aerial photos into parcel sized objects, and were further used for detailed land type identification within the parcels. Characteristics of image objects from contextual and geometrical aspects were used in the decision rule set to reduce the spectral limitation of the four-band aerial photography. Classification results include 12 land-cover classes and subclasses that may be assessed from the sub-parcel to the landscape scales, facilitating examination of scale dynamics. The proposed object-based classification method provides robust results, uses minimal and readily available ancillary data, and reduces computational time.     

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.     

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.