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.     

Next generation of global land cover characterization,mapping, and monitoring

Land cover change is increasingly affecting the biophysics, biogeochemistry, and biogeography of the Earth's surface and the atmosphere, with far-reaching consequences to human well-being. However, our scientific understanding of the distribution and dynamics of land cover and land cover change (LCLCC) is limited. Previous global land cover assessments performed using coarse spatial resolution (300 m–1 km) satellite data did not provide enough thematic detail or change information for global change studies and for resource management. High resolution (∼30 m) land cover characterization and monitoring is needed that permits detection of land change at the scale of most human activity and offers the increased flexibility of environmental model parameterization needed for global change studies. However, there are a number of challenges to overcome before producing such data sets including unavailability of consistent global coverage of satellite data, sheer volume of data, unavailability of timely and accurate training and validation data, difficulties in preparing image mosaics, and high performance computing requirements. Integration of remote sensing and information technology is needed for process automation and high-performance computing needs. Recent developments in these areas have created an opportunity for operational high resolution land cover mapping, and monitoring of the world. Here, we report and discuss these advancements and opportunities in producing the next generations of global land cover characterization, mapping, and monitoring at 30-m spatial resolution primarily in the context of United States, Group on Earth Observations Global 30 m land cover initiative (UGLC).     

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

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

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.     

Using remote sensing to monitor the influence of river discharge on watershed outlets and adjacent coral Reefs: Magdalena River and Rosario Islands,Colombia

Worldwide, coral reef ecosystems are being increasingly threatened by sediments loads from river discharges, which in turn are influenced by changing rainfall patterns due to climate change and by growing human activity in their watersheds. In this case study, we explored the applicability of using remote sensing (RS) technology to estimate and monitor the relationship between water quality at the coral reefs around the Rosario Islands, in the Caribbean Sea, and the rainfall patterns in the Magdalena River watershed. From the Moderate Resolution Imaging Spectroradiometer (MODIS), this study used the water surface reflectance product (MOD09GQ) to estimate water surface reflectance as a proxy for sediment concentration and the land cover product (MCD12Q1 V51) to characterize land cover of the watershed. Rainfall was estimated by using the 3B43 V7 product from the Tropical Rainforest Measuring Mission (TRMM). For the first trimester of each year, we investigated the inter-annual temporal variation in water surface reflectance at the Rosario Islands and at the three main mouths of the Magdalena River watershed. No increasing or decreasing trends of water surface reflectance were detected for any of the sites for the study period 2001–2014 (p > 0.05) but significant correlations were detected among the trends of each site at the watershed mouths (r = 0.57–0.90, p < 0.05) and between them and the inter-annual variation in rainfall on the watershed (r = 0.63–0.67, p < 0.05). Those trimesters with above-normal water surface reflectance at the mouths and above-normal rainfall at the watershed coincided with La Niña conditions while the opposite was the case during El Niño conditions. Although, a preliminary analysis of inter-annual land cover trends found only cropland cover in the watershed to be significantly correlated with water surface reflectance at two of the watershed mouths (r = 0.58 and 0.63, p < 0.05), the validation analysis draw only a 40.7% of accuracy in this land cover classification. This requires further analysis to confirm the impact of the cropland on the water quality at the watershed outlets. Spatial analysis with MOD09GQ imagery detected the overpass of river plumes from Barbacoas Bay over the Rosario Islands waters.     

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.     

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.     

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.     

Species’ habitat use inferred from environmental variables at multiple scales: How much we gain from high-resolution vegetation data?

Spatial resolution of environmental data may influence the results of habitat selection models. As high-resolution data are usually expensive, an assessment of their contribution to the reliability of habitat models is of interest for both researchers and managers. We evaluated how vegetation cover datasets of different spatial resolutions influence the inferences and predictive power of multi-scale habitat selection models for the endangered brown bear populations in the Cantabrian Range (NW Spain). We quantified the relative performance of three types of datasets: (i) coarse resolution data from Corine Land Cover (minimum mapping unit of 25 ha), (ii) medium resolution data from the Forest Map of Spain (minimum mapping unit of 2.25 ha and information on forest canopy cover and tree species present in each polygon), and (iii) high-resolution Lidar data (about 0.5 points/m²) providing a much finer information on forest canopy cover and height. Despite all the models performed well (AUC > 0.80), the predictive ability of multi-scale models significantly increased with spatial resolution, particularly when other predictors of habitat suitability (e.g. human pressure) were not used to indirectly filter out areas with a more degraded vegetation cover. The addition of fine grain information on forest structure (LiDAR) led to a better understanding of landscape use and a more accurate spatial representation of habitat suitability, even for a species with large spatial requirements as the brown bear, which will result in the development of more effective measures to assist endangered species conservation.     

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.     

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.     

Comparative assessment of thematic accuracy of GLC maps for specific applications using existing reference data

Inputs to various applications and models, current global land cover (GLC) maps are based on different data sources and methods. Therefore, comparing GLC maps is challenging. Statistical comparison of GLC maps is further complicated by the lack of a reference dataset that is suitable for validating multiple maps. This study utilizes the existing Globcover-2005 reference dataset to compare thematic accuracies of three GLC maps for the year 2005 (Globcover, LC-CCI and MODIS). We translated and reinterpreted the LCCS (land cover classification system) classifier information of the reference dataset into the different map legends. The three maps were evaluated for a variety of applications, i.e., general circulation models, dynamic global vegetation models, agriculture assessments, carbon estimation and biodiversity assessments, using weighted accuracy assessment. Based on the impact of land cover confusions on the overall weighted accuracy of the GLC maps, we identified map improvement priorities. Overall accuracies were 70.8 ± 1.4%, 71.4 ± 1.3%, and 61.3 ± 1.5% for LC-CCI, MODIS, and Globcover, respectively. Weighted accuracy assessments produced increased overall accuracies (80–93%) since not all class confusion errors are important for specific applications. As a common denominator for all applications, the classes mixed trees, shrubs, grasses, and cropland were identified as improvement priorities. The results demonstrate the necessity of accounting for dissimilarities in the importance of map classification errors for different user application. To determine the fitness of use of GLC maps, accuracy of GLC maps should be assessed per application; there is no single-figure accuracy estimate expressing map fitness for all purposes.     

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.     

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

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

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.     

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.     

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.     

Resolution vs. image quality in pre-tsunami imagery used for tsunami impact models in Aceh,Indonesia

Land cover roughness coefficients (LCRs) have been used in multivariate spatial models to test the mitigation potential of coastal vegetation to reduce impacts of the 2004 tsunami in Aceh, Indonesia. Previously, a Landsat 2002 satellite imagery was employed to derive land cover maps, which were then combined with vegetation characteristics, i.e., stand height, stem diameter and planting density to obtain LCRs. The present study tested LCRs extracted from 2003 and 2004 Landsat (30 m) images as well as a combination of 2003 and 2004 higher spatial resolution SPOT (10 m) imagery, while keeping the previous vegetation characteristics. Transects along the coast were used to extract land cover, whenever availability and visibility allowed. These new LCRs applied in previously developed tsunami impact models on wave outreach, casualties and damages confirmed previous findings regarding distance to the shoreline as a main factor reducing tsunami impacts. Nevertheless, the models using the new LCRs did not perform better than the original one. Particularly casualties models using 2002 LCRs performed better (δAIC > 2) than the more recent Landsat and SPOT counterparts. Cloud cover at image acquisition for Landsat and low area coverage for SPOT images decreased statistical predictive power (fewer observations). Due to the large spatial heterogeneity of tsunami characteristics as well as topographic and land-use features, it was more important to cover a larger area. Nevertheless, if more land cover classes would be referenced and high resolution imagery with low cloud cover would be available, the full benefits of higher spatial resolution imagery used to extract more precise land use roughness coefficients could be exploited.     

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.