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.     

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.     

Modelling high arctic percent vegetation cover using field digital images and high resolution satellite data

In this study, digital images collected at a study site in the Canadian High Arctic were processed and classified to examine the spatial-temporal patterns of percent vegetation cover (PVC). To obtain the PVC of different plant functional groups (i.e., forbs, graminoids/sedges and mosses), field near infrared-green-blue (NGB) digital images were classified using an object-based image analysis (OBIA) approach. The PVC analyses comparing different vegetation types confirmed: (i) the polar semi-desert exhibited the lowest PVC with a large proportion of bare soil/rock cover; (ii) the mesic tundra cover consisted of approximately 60% mosses; and (iii) the wet sedge consisted almost exclusively of graminoids and sedges. As expected, the PVC and green normalized difference vegetation index (GNDVI; (R_NIR − R_Green)/(R_NIR + R_Green)), derived from field NGB digital images, increased during the summer growing season for each vegetation type: i.e., ∼5% (0.01) for polar semi-desert; ∼10% (0.04) for mesic tundra; and ∼12% (0.03) for wet sedge respectively. PVC derived from field images was found to be strongly correlated with WorldView-2 derived normalized difference spectral indices (NDSI; (R_x − R_y)/(R_x + R_y)), where R_x is the reflectance of the red edge (724.1 nm) or near infrared (832.9 nm and 949.3 nm) bands; R_y is the reflectance of the yellow (607.7 nm) or red (658.8 nm) bands with R²’s ranging from 0.74 to 0.81. NDSIs that incorporated the yellow band (607.7 nm) performed slightly better than the NDSIs without, indicating that this band may be more useful for investigating Arctic vegetation that often includes large proportions of senescent vegetation throughout the growing season.     

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.     

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

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

Retrieval of leaf water content spanning the visible to thermal infrared spectra

The objective of this study was to investigate the entire spectra (from visible to the thermal infrared; 0.390–14.0 μm) to retrieve leaf water content in a consistent manner. Narrow-band spectral indices (calculated from all possible two band combinations) and a partial least square regression (PLSR) were used to assess the strength of each spectral region. The coefficient of determination (R²) and root mean square error (RMSE) were used to report the prediction accuracy of spectral indices and PLSR models. In the visible-near infrared and shortwave infrared (VNIR–SWIR), the most accurate spectral index yielded R² of 0.89 and RMSE of 7.60%, whereas in the mid infrared (MIR) the highest R² was 0.93 and RMSE of 5.97%. Leaf water content was poorly predicted using two-band indices developed from the thermal infrared (R² = 0.33). The most accurate PLSR model resulted from MIR reflectance spectra (R² = 0.96, RMSE = 4.74% and RMSE cross validation RMSE_CV = 6.17%) followed by VNIR–SWIR reflectance spectra (R² = 0.91, RMSE = 6.90% and RMSE_CV = 7.32%). Using thermal infrared (TIR) spectra, the PLSR model yielded a moderate retrieval accuracy (R² = 0.67, RMSE = 13.27% and RMSE_CV = 16.39%). This study demonstrated that the mid infrared (MIR) and shortwave infrared (SWIR) domains were the most sensitive spectral region for the retrieval of leaf water content.     

Predicting grain protein content of winter wheat using remote sensing data based on nitrogen status and water stress

Advanced site-specific knowledge of grain protein content of winter wheat from remote sensing data would provide opportunities to manage grain harvest differently, and to maximize output by adjusting input in fields. In this study, remote sensing data were utilized to predict grain protein content. Firstly, the leaf nitrogen content at winter wheat anthesis stage was proved to be significantly correlated with grain protein content (R² = 0.36), and spectral indices significantly correlated to leaf nitrogen content at anthesis stage were potential indicators for grain protein content. The vegetation index, VI_green, derived from the canopy spectral reflectance at green and red bands, was significantly correlated to the leaf nitrogen content at anthesis stage, and also highly significantly correlated to the final grain protein content (R² = 0.46). Secondly, the external conditions, such as irrigation, fertilization and temperature, had important influence on grain quality. Water stress at grain filling stage can increase grain protein content, and leaf water content is closely related to irrigation levels, therefore, the spectral indices correlated to leaf water content can be potential indicators for grain protein content. The spectral reflectance of TM channel 5 derived from canopy spectra or image data at grain filling stage was all significantly correlated to grain protein content (R² = 0.31 and 0.37, respectively). Finally, not only this study proved the feasibility of using remote sensing data to predict grain protein content, but it also provided a tentative prediction of the grain protein content in Beijing area using the reflectance image of TM channel 5.     

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).     

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.     

Towards decadal soil salinity mapping using Landsat time series data

Salinization is one of the major soil problems around the world. However, decadal variation in soil salinization has not yet been extensively reported. This study exploited thirty years (1985–2015) of Landsat sensor data, including Landsat-4/5 TM (Thematic Mapper), Landsat-7 ETM+ (Enhanced Thematic Mapper Plus) and Landsat-8 OLI (Operational Land Imager), for monitoring soil salinity of the Yellow River Delta, China. The data were initially corrected for atmospheric effects, and then matched the spectral bands of EO-1 (Earth Observing One) ALI (Advanced Land Imager). Subsequently, soil salinity maps were derived with a previously developed PLSR (Partial Least Square Regression) model. On intra-annual scale, the retrievals showed that soil salinity increased in February, stabilized in March, and decreased in April. On inter-annual scale, soil salinity decreased within 1985–2000 (−0.74 g kg⁻¹/10a, p < 0.001), and increased within 2000–2015 (0.79 g kg⁻¹/10a, p < 0.001). Our study presents a new perspective for use of multiple Landsat data in soil salinity retrieval, and further the understanding of soil salinization development over the Yellow River Delta.     

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.     

Uncertainty analysis of remote sensing of colored dissolved organic matter: Evaluations and comparisons for three rivers in North America

The uncertainties involved in remote sensing inversion of CDOM (Colored Dissolved Organic Matter) were analyzed in estuarine and coastal regions of three North American rivers: Mississippi, Hudson, and Neponset. Water optical and biogeochemical properties, including CDOM absorption and above-surface spectra, were collected in very high resolution. CDOM’s concentrations (a_g(440), absorption coefficient at 440 nm) were inverted from EO-1 Hyperion images, using a quasi-analytical algorithm for CDOM (QAA-CDOM). Uncertainties are classified to five levels, in which the underwater measurement uncertainty (level 1), image preprocessing uncertainty (level 4) and inverse model uncertainty (level 5) were evaluated. Results indicate that at level 1, in situ CDOM measurement is significant with 0.1 in the unit of QSU and 0.01 in the unit of a_g(440) (m⁻¹). At level 4, surface wave is a potential uncertainty source for high-resolution images in estuarine and coastal regions. The remote sensing reflectance of wavy water is about 10 times of the truth. At level 5, the overall uncertainty of QAA-CDOM inversion is 0.006 m⁻¹, with accuracy R² = 0.77, k = 1.1 and RMSE_log = 0.33 m⁻¹. The correlations between uncertainties and other water properties indicate that the large uncertainty in some rivers, such as the Neponset and Atchafalaya, might be caused by high-concentration chlorophyll or sediments. The relationships among the three level uncertainties show that the level 1 uncertainty generally does not propagate into level 4 and 5, but the large uncertainty at level 4 usually introduce large uncertainty at level 5.     

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

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

Retrieval of tea polyphenol at leaf level using spectral transformation and multi-variate statistical approach

In the present study, field based hyperspectral data was used to estimate the tea (Camellia sinensis L.) polyphenol at Deha Tea garden of Assam state, India. Leaf reflectance spectra were first filtered for noise and then transformed into normalized and first derivative reflectance for further analysis. Stepwise discriminant analysis was carried out to select sensitive bands for a range of polyphenol concentration by minimizing the effects of other factors such as age of the bushes and management practices. The wavelengths at 358, 369, 484, 845, 916, 1387, 1420, 1435, 1621 and 2294 nm were identified as sensitive to tea polyphenol, among which 2294 nm was found to be the most recurring band. The noise removed selected bands, their transformed derivatives and principal components were regressed with the tea polyphenol using univariate and multi-variate analysis. In univariate analysis the correlation was very poor with RMSE more than 3.0. A significant improvement in R² values were observed when multivariate analyses like stepwise multiple linear regression (SMLR) and partial least square regression (PLSR) was carried out. The PLSR of first derivative reflectance was most accurate (R² = 0.81 and RMSE = 1.39 mg g⁻¹) among all the uni- and multivariate analysis for predicting the polyphenol of fresh tea leaves.     

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.     

Evaluating the relationship between biomass,percent groundcover and remote sensing indices across six winter cover crop fields in Maryland,United States

Winter cover crops are an essential part of managing nutrient and sediment losses from agricultural lands. Cover crops lessen sedimentation by reducing erosion, and the accumulation of nitrogen in aboveground biomass results in reduced nutrient runoff. Winter cover crops are planted in the fall and are usually terminated in early spring, making them susceptible to senescence, frost burn, and leaf yellowing due to wintertime conditions. This study sought to determine to what extent remote sensing indices are capable of accurately estimating the percent groundcover and biomass of winter cover crops, and to analyze under what critical ranges these relationships are strong and under which conditions they break down. Cover crop growth on six fields planted to barley, rye, ryegrass, triticale or wheat was measured over the 2012–2013 winter growing season. Data collection included spectral reflectance measurements, aboveground biomass, and percent groundcover. Ten vegetation indices were evaluated using surface reflectance data from a 16-band CROPSCAN sensor. Restricting analysis to sampling dates before the onset of prolonged freezing temperatures and leaf yellowing resulted in increased estimation accuracy. There was a strong relationship between the normalized difference vegetation index (NDVI) and percent groundcover (r² = 0.93) suggesting that date restrictions effectively eliminate yellowing vegetation from analysis. The triangular vegetation index (TVI) was most accurate in estimating high ranges of biomass (r² = 0.86), while NDVI did not experience a clustering of values in the low and medium biomass ranges but saturated in the higher range (>1500 kg/ha). The results of this study show that accounting for index saturation, senescence, and frost burn on leaves can greatly increase the accuracy of estimates of percent groundcover and biomass for winter cover crops.     

Retrieval of leaf area index in alpine wetlands using a two-layer canopy reflectance model

In this paper, we focused on the retrieval of the LAI in an alpine wetland located in western part of China in late August and early July 2011. A two-layer canopy reflectance model (ACRM) was used to establish the relationships between the LAI and the reflectance of near-infrared (NIR) and red (RED) wavebands. The reflectance data were derived from Landsat TM L1T product and the Terra and Aqua MODIS 16-day and 8-day composite reflectance products (MOD/MYD09) at 250 m resolution. Due to the lack of the information about some major input parameters for ACRM, which are sensitive to model outputs in the reflectance of NIR and RED wavebands, the inverse problem was ill-posed. To overcome this problem, a method of increasing the sensitivity of the LAI while reducing the influence of other model free parameters based on the study of free parameters’ sensitivity to the ACRM outputs and the region’s features was studied. The area of interest was divided into two parts using the approximately statistic normalized difference vegetation index (NDVI) value around 0.5. One part was sparse vegetation (0.1 < NDVI < 0.5), which is more sensitive to soil background effects and less sensitive to the canopy biophysical and biochemical variables. The other part was dense vegetation (0.5 ≤ NDVI < 1.0), which is less sensitive to soil background effects and more sensitive to plant canopies and leaf parameters. Then, the relationships of ρ_nir–LAI and ρ_red–LAI were established using a look-up table algorithm for the two parts. Furthermore, a regularization technique for fast pixel-wise retrieval was introduced to reduce the elements of LUT sets while maintaining a relatively high accuracy. The results were very promising compared to the field measured LAI values that the correlation (R²) of the measured LAI values and retrieved LAI values reached 0.95, and the root-mean-square deviation (RMSD) was 0.33 for late August, 2011, while the R² reached 0.82 and RMSD was 0.25 for early July 2011.     

Empirical models for estimating the suspended sediment concentration in Amazonian white water rivers using Landsat 5/TM

Suspended sediment yield is a very important environmental indicator within Amazonian fluvial systems, especially for rivers dominated by inorganic particles, referred to as white water rivers. For vast portions of Amazonian rivers, suspended sediment concentration (SSC) is measured infrequently or not at all. However, remote sensing techniques have been used to estimate water quality parameters worldwide, from which data for suspended matter is the most successfully retrieved. This paper presents empirical models for SSC retrieval in Amazonian white water rivers using reflectance data derived from Landsat 5/TM. The models use multiple regression for both the entire dataset (global model, N = 504) and for five segmented datasets (regional models) defined by general geological features of drainage basins. The models use VNIR bands, band ratios, and the SWIR band 5 as input. For the global model, the adjusted R² is 0.76, while the adjusted R² values for regional models vary from 0.77 to 0.89, all significant (p-value < 0.0001). The regional models are subject to the leave-one-out cross validation technique, which presents robust results. The findings show that both the average error of estimation and the standard deviation increase as the SSC range increases. Regional models were more accurate when compared with the global model, suggesting changes in optical proprieties of water sampled at different sampling stations. Results confirm the potential for the estimation of SSC from Landsat/TM historical series data for the 1980s and 1990s, for which the in situ database is scarce. Such estimates supplement the SSC temporal series, providing a more comprehensive SSC temporal series which may show environmental dynamics yet unknown.     

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).