ALBERT

Description: The objective of this study was to assess the performance of the simulated start of the photosynthetically active season by a large-scale biosphere model in boreal forests in Finland with remote sensing observations. The start of season for two forest types, evergreen needle- and deciduous broad-leaf, was obtained for the period 2003–2011 from regional JSBACH (Jena Scheme for Biosphere–Atmosphere Hamburg) runs, driven with climate variables from a regional climate model. The satellite-derived start of season was determined from daily Moderate Resolution Imaging Spectrometer (MODIS) time series of Fractional Snow Cover and the Normalized Difference Water Index by applying methods that were targeted to the two forest types. The accuracy of the satellite-derived start of season in deciduous forest was assessed with bud break observations of birch and a root mean square error of seven days was obtained. The evaluation of JSBACH modelled start of season dates with satellite observations revealed high spatial correspondence. The bias was less than five days for both forest types but showed regional differences that need further consideration. The agreement with satellite observations was slightly better for the evergreen than for the deciduous forest. Nonetheless, comparison with gross primary production (GPP) determined from CO2 flux measurements at two eddy covariance sites in evergreen forest revealed that the JSBACH-simulated GPP was higher in early spring and led to too-early simulated start of season dates. Photosynthetic activity recovers differently in evergreen and deciduous forests. While for the deciduous forest calibration of phenology alone could improve the performance of JSBACH, for the evergreen forest, changes such as seasonality of temperature response, would need to be introduced to the photosynthetic capacity to improve the temporal development of gross primary production.

Description: An innovative remote sensing approach that combines land-use change and water quality information is proposed in order to investigate if Artisanal and Small-scale Gold Mining (ASGM) area extension is associated with water siltation in the Tapajós River Basin (Brazil), containing the largest small-scale gold mining district in the world. Taking advantage of a 40-year period of the multi-satellite imagery archive, the objective of this paper is to build a normalized time-series in order to evaluate the influence of temporal mining expansion on the water siltation data (TSS, Total Suspended Solids concentration) derived from previous research. The methodological approach was set to deliver a full characterization of the ASGM expansion from its initial stages in the early 1970s to the present. First, based on IRS/LISSIII images acquired in 2012, the historical Landsat image database (1973–2001) was corrected for radiometric and atmospheric effects using dark vegetation as reference to create a normalized time-series. Next, a complete update of the mining areas distribution in 2012 derived from the TerraClass Project (an official land-use classification for the Brazilian Amazon) was conducted having IRS/LISSIII as the base map with the support of auxiliary data and vector editing. Once the ASGM in 2012 was quantified (261.7 km2) and validated with photos, a reverse classification of ASGM in 2001 (171.7 km2), 1993 (166.3 km2), 1984 (47.5 km2), and 1973 (15.4 km2) with the use of Landsat archives was applied. This procedure relies on the assumption that ASGM changes in the land cover are severe and remain detectable from satellite sensors for decades. The mining expansion area over time was then combined with the (TSS) data retrieved from the same atmospherically corrected satellite imagery based on the literature. In terms of gold mining expansion and water siltation effects, four main periods of ASGM activities were identified in the study area: (i) 1958–1977, first occurrence of mining activities and low water impacts; (ii) 1978–1993, introduction of low-budget mechanization associated with very high gold prices resulting in large mining area expansion and high water siltation levels; (iii) 1994–2003, general recession of ASGM activities and exhaustion of easy-access gold deposits, resulting in decreased TSS; (iv) 2004 to present, intensification of ASGM encouraged by high gold prices, resulting in an increase of TSS.

Description: Current earth observation models do not take into account the influence of water salinity on the evaporation rate, even though the salinity influences the evaporation rate by affecting the density and latent heat of vaporization. In this paper, we adapt the SEBS (Surface Energy Balance System) model for large water bodies and add the effect of water salinity to the evaporation rate. Firstly, SEBS is modified for fresh-water whereby new parameterizations of the water heat flux and sensible heat flux are suggested. This is achieved by adapting the roughness heights for momentum and heat transfer. Secondly, a salinity correction factor is integrated into the adapted model. Eddy covariance measurements over Lake IJsselmeer (The Netherlands) are carried out and used to estimate the roughness heights for momentum (~0.0002 m) and heat transfer (~0.0001 m). Application of these values over the Victoria and Tana lakes (freshwater) in Africa showed that the calculated latent heat fluxes agree well with the measurements. The root mean-square of relative-errors (rRMSE) is about 4.1% for Lake Victoria and 4.7%, for Lake Tana. Verification with ECMWF data showed that the salinity reduced the evaporation at varying levels by up to 27% in the Great Salt Lake and by 1% for open ocean. Our results show the importance of salinity to the evaporation rate and the suitability of the adapted-SEBS model (AquaSEBS) for fresh and saline waters.

Description: Area-based analyses of airborne laser scanning (ALS) data are an established approach to obtain wall-to-wall predictions of forest characteristics for vast areas. The analyses of sparse data in particular are based on the height value distributions, which do not produce optimal information on the horizontal forest structure. We evaluated the complementary potential of features quantifying the textural variation of ALS-based canopy height models (CHMs) for both supervised (linear regression) and unsupervised (k-Means clustering) analyses. Based on a comprehensive literature review, we identified a total of four texture analysis methods that produced rotation-invariant features of different order and scale. The CHMs and the textural features were derived from practical sparse-density, leaf-off ALS data originally acquired for ground elevation modeling. The features were extracted from a circular window of 254 m2 and related with boreal forest characteristics observed from altogether 155 field sample plots. Features based on gray-level histograms, distribution of forest patches, and gray-level co-occurrence matrices were related with plot volume, basal area, and mean diameter with coefficients of determination (R2) of up to 0.63–0.70, whereas features that measured the uniformity of local binary patterns of the CHMs performed poorer. Overall, the textural features compared favorably with benchmark features based on the point data, indicating that the textural features contain additional information useful for the prediction of forest characteristics. Due to the developed processing routines for raster data, the CHM features may potentially be extracted with a lower computational burden, which promotes their use for applications such as pre-stratification or guiding the field plot sampling based solely on ALS data.

Description: Various themed streets have recently been developed by local governments in order to stimulate local economies and to establish the identity of the corresponding places. However, the motivations behind the development of some of these themed street projects has been based on profit, without full considerations of people’s perceptions of their local areas, resulting in marginal effects on the local economies concerned. In response to this issue, this study proposed a themed street clustering method to detect the themed streets of a specific region, focusing on the commercial themed street, which is more prevalent than other types of themed streets using location based service data. This study especially uses “the street segment” as a basic unit for analysis. The Sillim and Gangnam areas of Seoul, South Korea were chosen for the evaluation of the adequacy of the proposed method. By comparing trade areas that were sourced from a market analysis report by a reliable agent with the themed streets detected in this study, the experiment results showed high proficiency of the proposed method.

Description: Delimitating trade areas is a major business concern. Today, mobile communication technologies make it possible to use social media data for this purpose. Few studies however, have focused on methods to extract suitable samples from social media data for trade area delimitation. In our case study, we divided Beijing into regular grid cells and extracted activity centers for each social media user. Ten sample sets were obtained by selecting users based on the retail agglomerations they visited and aggregating user activity centers to each grid cell. We calculated distance and visitation frequency attributes for each user and each grid cell. The distance value of a grid cell is the average distance of user activity centers in this grid cell to a retail agglomeration. The visitation frequency of a grid cell refers to the average count of visits to retail agglomerations by user activity centers for a cell. The calculated attribute values of 10 sets were input into a Huff model and the delimitated trade areas were evaluated. Results show that sets obtained by aggregating user activity centers have a better delimitating effect than sets obtained without aggregation. Differences in the distribution and intensity of trade areas also became apparent.

Description: Road information is fundamental not only in the military field but also common daily living. Automatic road extraction from a remote sensing images can provide references for city planning as well as transportation database and map updating. However, owing to the spectral similarity between roads and impervious structures, the current methods solely using spectral characteristics are often ineffective. By contrast, the detailed information discernible from the high-resolution aerial images enables road extraction with spatial texture features. In this study, a knowledge-based method is established and proposed; this method incorporates the spatial texture feature into urban road extraction. The spatial texture feature is initially extracted by the local Moran’s I, and the derived texture is added to the spectral bands of image for image segmentation. Subsequently, features like brightness, standard deviation, rectangularity, aspect ratio, and area are selected to form the hypothesis and verification model based on road knowledge. Finally, roads are extracted by applying the hypothesis and verification model and are post-processed based on the mathematical morphology. The newly proposed method is evaluated by conducting two experiments. Results show that the completeness, correctness, and quality of the results could reach approximately 94%, 90% and 86% respectively, indicating that the proposed method is effective for urban road extraction.

Description: Effectively identifying an airport from satellite and aerial imagery is a challenging task. Traditional methods mainly focus on the use of multiple features for the detection of runways and some also adapt knowledge of airports, but the results are unsatisfactory and the usage limited. A new method is proposed to recognize airports from high-resolution optical images. This method involves the analysis of the saliency distribution and the use of fuzzy rule-based classification. First, a number of images with and without airports are segmented into multiple scales to obtain a saliency distribution map that best highlights the saliency distinction between airports and other objects. Then, on the basis of the segmentation result and the structural information of airports, we analyze the segmentation result to extract and represent the semantic information of each image via the bag-of-visual-words (BOVW) model. The image correlation degree is combined with the BOVW model and fractal dimension calculation to make a more complete description of the airports and to carry out preliminary classification. Finally, the support vector machine (SVM) is adopted for detailed classification to classify the remaining imagery. The experiment shows that the proposed method achieves a precision of 89.47% and a recall of 90.67% and performs better than other state of the art methods on precision and recall.

Description: Inland saline wetlands are complex systems undergoing continuous changes in moisture and salinity and are especially vulnerable to human pressures. Remote sensing is helpful to identify vegetation change in semi-arid wetlands and to assess wetland degradation. Remote sensing-based monitoring requires identification of the spectral characteristics of soils and vegetation and their correspondence with the vegetation cover and soil conditions. We studied the spectral characteristics of soils and vegetation of saline wetlands in Monegros, NE Spain, through field and satellite images. Radiometric and complementary field measurements in two field surveys in 2007 and 2008 were collected in selected sites deemed as representative of different soil moisture, soil color, type of vegetation, and density. Despite the high local variability, we identified good relationships between field spectral data and Quickbird images. A methodology was established for mapping the fraction of vegetation cover in Monegros and other semi-arid areas. Estimating vegetation cover in arid wetlands is conditioned by the soil background and by the occurrence of dry and senescent vegetation accompanying the green component of perennial salt-tolerant plants. Normalized Difference Vegetation Index (NDVI) was appropriate to map the distribution of the vegetation cover if the green and yellow-green parts of the plants are considered.

Description: Extreme heat weather events have received increasing attention and has become of special importance as they can remarkably affect sectors as diverse as public health, energy consumption, water resources, natural biodiversity and agricultural production. In this regard, summer temperatures have become a parameter of essential interest under a framework of a hypothetical increase in the number of intense-heat conditions. Thus, their forecast is a crucial aspect bearing in mind a mitigation of the effects and impacts that these intense-heat situations could produce. The current work tries to reach a better understanding of these sorts of situations that are really common over the Western Mediterranean coast. An extreme heat episode that took place in the Valencia Region in July 2009 is analysed, based on the simulations performed with the Regional Atmospheric Modeling System (RAMS). This event recorded maximum temperatures exceeding 40 °C amply extended over the region besides reaching minimum temperatures up to 25.92 °C. We examine the role of improved skin and soil temperature (ST) initial conditions in the forecast results by means of different modelling and satellite-derived products. The influence of incorporating the Land Surface Temperature (LST) into RAMS is not found to produce a meaningful impact on the simulation results, independently of the resolution of the dataset used in the initial conditions of the model. In contrast, the introduction of the ST in lower levels, not only the skin temperature, has a more marked decisive effect in the simulation. Additionally, we have evaluated the influence of increasing the number of soil levels to spread deeper underground. This sensitivity experiment has revealed that more soil levels do not produce any meaningful impact on the simulation compared to the original one. In any case, RAMS is able to properly capture the observed patterns in those cases where a Western advection is widely extended over the area of study. This region’s variability in orography and in distances to the sea promotes the development of sea-breeze circulations, thus producing a convergence of two opposite wind flows, a Western synoptic advection and a sea-breeze circulation. As a result, the RAMS skill in those cases where a sea breeze is well developed depends on the proper location of the boundary and convergence lines of these two flows.

Description: Spatially explicit knowledge of aboveground biomass (AGB) in large areas is important for accurate carbon accounting and quantifying the effect of forest disturbance on the terrestrial carbon cycle. We estimated AGB from 1990 to 2011 in northern Guangdong, China, based on a spatially explicit dataset derived from six years of national forest inventory (NFI) plots, Landsat time series imagery (1986–2011) and Advanced Land Observing Satellite (ALOS) Phased Array L-band Synthetic Aperture Radars (PALSAR) 25 m mosaic data (2007–2010). Four types of variables were derived for modeling and assessment. The random forest approach was used to seek the optimal variables for mapping and validation. The root mean square error (RMSE) of plot-level validation was between 6.44 and 39.49 (t/ha), the normalized root-mean-square error (NRMSE) was between 7.49% and 19.01% and mean absolute error (MAE) was between 5.06 and 23.84 t/ha. The highest coefficient of determination R2 of 0.8 and the lowest NRMSE of 7.49% were reported in 2006. A clear increasing trend of mean AGB from the lowest value of 13.58 t/ha to the highest value of 66.25 t/ha was witnessed between 1988 and 2000, while after 2000 there was a fluctuating ascending change, with a peak mean AGB of 67.13 t/ha in 2004. By integrating AGB change with forest disturbance, the trend in disturbance area closely corresponded with the trend in AGB decrease. To determine the driving forces of these changes, the correlation analysis was adopted and exploratory factor analysis (EFA) method was used to find a factor rotation that maximizes this variance and represents the dominant factors of nine climate elements and nine human activities elements affecting the AGB dynamics. Overall, human activities contributed more to short-term AGB dynamics than climate data. Harvesting and human-induced fire in combination with rock desertification and global warming made a strong contribution to AGB changes. This study provides valuable information for the relationships between forest AGB and climate as well as forest disturbance in subtropical zones.

Description: Merging satellite and rain gauge data by combining accurate quantitative rainfall from stations with spatial continuous information from remote sensing observations provides a practical method of estimating rainfall. However, generating high spatiotemporal rainfall fields for catchment-distributed hydrological modeling is a problem when only a sparse rain gauge network and coarse spatial resolution of satellite data are available. The objective of the study is to present a satellite and rain gauge data-merging framework adapting for coarse resolution and data-sparse designs. In the framework, a statistical spatial downscaling method based on the relationships among precipitation, topographical features, and weather conditions was used to downscale the 0.25° daily rainfall field derived from the Tropical Rainfall Measuring Mission (TRMM) Multisatellite Precipitation Analysis (TMPA) precipitation product version 7. The nonparametric merging technique of double kernel smoothing, adapting for data-sparse design, was combined with the global optimization method of shuffled complex evolution, to merge the downscaled TRMM and gauged rainfall with minimum cross-validation error. An indicator field representing the presence and absence of rainfall was generated using the indicator kriging technique and applied to the previously merged result to consider the spatial intermittency of daily rainfall. The framework was applied to estimate daily precipitation at a 1 km resolution in the Qinghai Lake Basin, a data-scarce area in the northeast of the Qinghai-Tibet Plateau. The final estimates not only captured the spatial pattern of daily and annual precipitation with a relatively small estimation error, but also performed very well in stream flow simulation when applied to force the geomorphology-based hydrological model (GBHM). The proposed framework thus appears feasible for rainfall estimation at high spatiotemporal resolution in data-scarce areas.

Description: Kernel-based methods and ensemble learning are two important paradigms for the classification of hyperspectral remote sensing images. However, they were developed in parallel with different principles. In this paper, we aim to combine the advantages of kernel and ensemble methods by proposing a kernel supervised ensemble classification method. In particular, the proposed method, namely RoF-KOPLS, combines the merits of ensemble feature learning (i.e., Rotation Forest (RoF)) and kernel supervised learning (i.e., Kernel Orthonormalized Partial Least Square (KOPLS)). In particular, the feature space is randomly split into K disjoint subspace and KOPLS is applied to each subspace to produce the new features set for the training of decision tree classifier. The final classification result is assigned to the corresponding class by the majority voting rule. Experimental results on two hyperspectral airborne images demonstrated that RoF-KOPLS with radial basis function (RBF) kernel yields the best classification accuracies due to the ability of improving the accuracies of base classifiers and the diversity within the ensemble, especially for the very limited training set. Furthermore, our proposed method is insensitive to the number of subsets.

Description: Leaf Area Index (LAI) is a key variable that bridges remote sensing observations to the quantification of agroecosystem processes. In this study, we assessed the universality of the relationships between crop LAI and remotely sensed Vegetation Indices (VIs). We first compiled a global dataset of 1459 in situ quality-controlled crop LAI measurements and collected Landsat satellite images to derive five different VIs including Simple Ratio (SR), Normalized Difference Vegetation Index (NDVI), two versions of the Enhanced Vegetation Index (EVI and EVI2), and Green Chlorophyll Index (CIGreen). Based on this dataset, we developed global LAI-VI relationships for each crop type and VI using symbolic regression and Theil-Sen (TS) robust estimator. Results suggest that the global LAI-VI relationships are statistically significant, crop-specific, and mostly non-linear. These relationships explain more than half of the total variance in ground LAI observations (R2 > 0.5), and provide LAI estimates with RMSE below 1.2 m2/m2. Among the five VIs, EVI/EVI2 are the most effective, and the crop-specific LAI-EVI and LAI-EVI2 relationships constructed by TS, are robust when tested by three independent validation datasets of varied spatial scales. While the heterogeneity of agricultural landscapes leads to a diverse set of local LAI-VI relationships, the relationships provided here represent global universality on an average basis, allowing the generation of large-scale spatial-explicit LAI maps. This study contributes to the operationalization of large-area crop modeling and, by extension, has relevance to both fundamental and applied agroecosystem research.

Description: With its temporal resolution of 10 days (five days with two satellites, and significantly more at high latitudes), its swath width of 290 km, and its 10 m and 20 m spatial resolution bands from the visible to the shortwave infrared, the European Sentinel-2 satellites have significant potential for glacier remote sensing, in particular mapping of glacier outlines and facies, and velocity measurements. Testing Level 1C commissioning and ramp-up phase data for initial sensor quality experiences, we find a high radiometric performance, but with slight striping effects under certain conditions. Through co-registration of repeat Sentinal-2 data we also find lateral offset patterns and noise on the order of a few metres. Neither of these issues will complicate most typical glaciological applications. Absolute geo-location of the data investigated was on the order of one pixel at the time of writing. The most severe geometric problem stems from vertical errors of the DEM used for ortho-rectifying Sentinel-2 data. These errors propagate into locally varying lateral offsets in the images, up to several pixels with respect to other georeferenced data, or between Sentinel-2 data from different orbits. Finally, we characterize the potential and limitations of tracking glacier flow from repeat Sentinel-2 data using a set of typical glaciers in different environments: Aletsch Glacier, Swiss Alps; Fox Glacier, New Zealand; Jakobshavn Isbree, Greenland; Antarctic Peninsula at the Larsen C ice shelf.

Description: Currently there is a lack of quantitative information regarding the driving factors of vegetation dynamics in post-Soviet Central Asia. Insufficient knowledge also exists concerning vegetation variability across sub-humid to arid climatic gradients as well as vegetation response to different land uses, from natural rangelands to intensively irrigated croplands. In this study, we analyzed the environmental drivers of vegetation dynamics in five Central Asian countries by coupling key vegetation parameter “overall greenness” derived from Moderate Resolution Imaging Spectroradiometer (MODIS) Normalized Difference Vegetation Index (NDVI time series data, with its possible factors across various management and climatic gradients. We developed nine generalized least-squares random effect (GLS-RE) models to analyze the relative impact of environmental factors on vegetation dynamics. The obtained results quantitatively indicated the extensive control of climatic factors on managed and unmanaged vegetation cover across Central Asia. The most diverse vegetation dynamics response to climatic variables was observed for “intensively managed irrigated croplands”. Almost no differences in response to these variables were detected for managed non-irrigated vegetation and unmanaged (natural) vegetation across all countries. Natural vegetation and rainfed non-irrigated crop dynamics were principally associated with temperature and precipitation parameters. Variables related to temperature had the greatest relative effect on irrigated croplands and on vegetation cover within the mountainous zone. Further research should focus on incorporating the socio-economic factors discussed here in a similar analysis.

Description: An open, standardized data management and related service infrastructure is a crucial requirement for a seamless storage and exchange of data and information within research projects, for the dissemination of project results and for their application in decision making processes. However, typical project databases often refer to only one research project and are limited to specific purposes. Once implemented, those systems are often not further maintained and updated, rendering the stored information useless once the system stops operating. The River Basin Information System (RBIS) presented here is designed to fit not only the requirements of one research project, but focuses on generic functions, extensibility and standards compliance typically found in interdisciplinary environmental research. Developed throughout more than 10 years of research cooperation worldwide, RBIS is designed to manage different types of environmental data with and without spatial context together with a rich set of metadata. Beside data management and storage, RBIS provides functions for the visualization, linking, analysis and processing of different types of data to support research, decision making, result dissemination and information discovery for all kinds of users. The focus of this paper is on the description of the technical implementation and the presentation of functions. This will be complemented by an overview of example applications and experiences during RBIS development and operation.

Description: In this paper, we consider retrieving individual wave components in a multi-directional sea wave model. To solve this problem, a currently and commonly used method is three-dimensional discrete Fourier transform (3D DFT) on the radar image sequence. However, the uniform frequency and the uniform wavenumber in a wavenumber frequency domain can not always strictly satisfy the dispersion relation, and the spectral leakage in both temporal and spatial domains exists due to the limited analysis area selected from an image sequence. As a result, the DFT method incurs undesirable error performance in retrieving directional wave components. By deeply investigating the data structure of the multi-directional sea wave model, we obtain a new and decomposable matrix representation for processing the wave components. Then, a novel successive cancellation method is proposed to efficiently and effectively extract individual wave components, whose frequency and wavenumber rigorously satisfy the liner dispersion relation. Thus, it avoids spectral leakage in the spatial domain. The algorithm is evaluated by using linear synthetic wave image sequences. The validity of the proposed novel algorithm is verified by comparing the retrieved parameters of amplitude, phase, and direction of the individual wave components with the simulated parameters as well as those obtained by using the 3D DFT method. In addition, the reconstructed sea field using the retrieved wave components is also compared with the simulated remote sensing images as well as those attained using the inverse 3D DFT method. All the simulation results demonstrate that our proposed algorithm is more effective and has better performance for retrieving individual wave components from the spatio-temporal remote sensing image sequences than the 3D DFT method.

Description: Knowing the current phenological state of an agricultural crop is a powerful tool for precision farming applications. In the past, it has been estimated with remote sensing data by exploiting time series of Normalised Difference Vegetation Index (NDVI), but always at the end of the campaign and only providing results for some key states. In this work, a new dynamical framework is proposed to provide real-time estimates in a continuous range of states, for which NDVI images are combined with a prediction model in an optimal way using a particle filter. The methodology is tested over a set of 8 to 13 rice parcels during 2008–2013, achieving a high determination factor R 2 = 0.93 ( n = 379 ) for the complete phenological range. This method is also used to predict the end of season date, obtaining a high accuracy with an anticipation of around 40–60 days. Among the key advantages of this approach, phenology is estimated each time a new observation is available, hence enabling the potential detection of anomalies in real-time during the cultivation. In addition, the estimation procedure is robust in the case of noisy observations, and it is not limited to a few phenological stages.

Description: As a sensitive indicator of climate change, mountain glacier dynamics are of great concern, but the ice motion pattern of an entire glacier surface cannot be accurately and efficiently generated by the use of only phase-based or intensity-based methods with synthetic aperture radar (SAR) imagery. To derive the ice movement of the whole glacier surface with a high accuracy, an integrated approach combining differential interferometric SAR (D-InSAR), multi-aperture interferometry (MAI), and a pixel-tracking (PT) method is proposed, which could fully exploit the phase and intensity information recorded by the SAR sensor. The Chongce Glacier surface flow field is estimated with the proposed integrated approach. Compared with the traditional SAR-based methods, the proposed approach can determine the ice motion over a widely varying range of ice velocities with a relatively high accuracy. Its capability is proved by the detailed ice displacement pattern with the average accuracy of 0.2 m covering the entire Chongce Glacier surface, which shows a maximum ice movement of 4.9 m over 46 days. Furthermore, it is shown that the ice is in a quiescent state in the downstream part of the glacier. Therefore, the integrated approach presented in this paper could present us with a novel way to comprehensively and accurately understand glacier dynamics by overcoming the incoherence phenomenon, and has great potential for glaciology study.

Description: With the upcoming availability of the next generation of high quality orbiting hyperspectral sensors, a major step toward improved regional soil mapping and monitoring and delivery of quantitative soil maps is expected. This study focuses on the determination of the prediction accuracy of spectral models for the mapping of common soil properties based on upcoming EnMAP (Environmental Mapping and Analysis Program) satellite data using semi-operational soil models. Iron oxide (Fed), clay, and soil organic carbon (SOC) content are predicted in test areas in Spain and Luxembourg based on a semi-automatic Partial-Least-Square (PLS) regression approach using airborne hyperspectral, simulated EnMAP, and soil chemical datasets. A variance contribution analysis, accounting for errors in the dependent variables, is used alongside classical error measurements. Results show that EnMAP allows predicting iron oxide, clay, and SOC with an R2 between 0.53 and 0.67 compared to Hyperspectral Mapper (HyMap)/Airborne Hyperspectral System (AHS) imagery with an R2 between 0.64 and 0.74. Although a slight decrease in soil prediction accuracy is observed at the spaceborne scale compared to the airborne scale, the decrease in accuracy is still reasonable. Furthermore, spatial distribution is coherent between the HyMap/AHS mapping and simulated EnMAP mapping as shown with a spatial structure analysis with a systematically lower semivariance at the EnMAP scale.

Description: This study uses a large-scale mobile phone dataset to estimate potential demand of bicycle trips in a city. By identifying two important anchor points (night-time anchor point and day-time anchor point) from individual cellphone trajectories, this study proposes an anchor-point based trajectory segmentation method to partition cellphone trajectories into trip chain segments. By selecting trip chain segments that can potentially be served by bicycles, two indicators (inflow and outflow) are generated at the cellphone tower level to estimate the potential demand of incoming and outgoing bicycle trips at different places in the city and different times of a day. A maximum coverage location-allocation model is used to suggest locations of bike sharing stations based on the total demand generated at each cellphone tower. Two measures are introduced to further understand characteristics of the suggested bike station locations: (1) accessibility; and (2) dynamic relationships between incoming and outgoing trips. The accessibility measure quantifies how well the stations could serve bicycle users to reach other potential activity destinations. The dynamic relationships reflect the asymmetry of human travel patterns at different times of a day. The study indicates the value of mobile phone data to intelligent spatial decision support in public transportation planning.

Description: Recent drastic reductions in the Arctic sea-ice cover have raised an interest in understanding the role of sea ice in the global system as well as pointed out a need to understand the physical processes that lead to such changes. Satellite remote-sensing data provide important information about remote ice areas, and Synthetic Aperture Radar (SAR) data have the advantages of penetration of the omnipresent cloud cover and of high spatial resolution. A challenge addressed in this paper is how to extract information on sea-ice types and sea-ice processes from SAR data. We introduce, validate and apply geostatistical and statistical approaches to automated classification of sea ice from SAR data, to be used as individual tools for mapping sea-ice properties and provinces or in combination. A key concept of the geostatistical classification method is the analysis of spatial surface structures and their anisotropies, more generally, of spatial surface roughness, at variable, intermediate-sized scales. The geostatistical approach utilizes vario parameters extracted from directional vario functions, the parameters can be mapped or combined into feature vectors for classification. The method is flexible with respect to window sizes and parameter types and detects anisotropies. In two applications to RADARSAT and ERS-2 SAR data from the area near Point Barrow, Alaska, it is demonstrated that vario-parameter maps may be utilized to distinguish regions of different sea-ice characteristics in the Beaufort Sea, the Chukchi Sea and in Elson Lagoon. In a third and a fourth case study the analysis is taken further by utilizing multi-parameter feature vectors as inputs for unsupervised and supervised statistical classification. Field measurements and high-resolution aerial observations serve as basis for validation of the geostatistical-statistical classification methods. A combination of supervised classification and vario-parameter mapping yields best results, correctly identifying several sea-ice provinces in the shore-fast ice and the pack ice. Notably, sea ice does not have to be static to be classifiable with respect to spatial structures. In consequence, the geostatistical-statistical classification may be applied to detect changes in ice dynamics, kinematics or environmental changes, such as increased melt ponding, increased snowfall or changes in the equilibrium line.

Description: Urban flooding is a serious natural hazard to many cities all over the world, which has dramatic impacts on the urban environment and human life. Urban flooding mapping has practical significance for the prevention and management of urban flood disasters. Remote sensing images with high temporal resolutions are widely used for urban flooding mapping, but have a limitation of relatively low spatial resolutions. In this study, a new method based on a generalized regression neural network (GRNN) is proposed to achieve improved accuracy in super-resolution mapping of urban flooding (SMUF) from remote sensing images. The GRNN-SMUF algorithm was proposed and then assessed using Landsat 5 and Landsat 8 images of Brisbane city in Australia and Wuhan city in China. Compared to three traditional methods, GRNN-SMUF mapped urban flooding more accurately according to both visual and quantitative assessments. The results of this study will improve the accuracy of urban flooding mapping using easily-available remote sensing images with medium-low spatial resolutions and will be propitious to the prevention and management of urban flood disasters.

Description: The availability of a dense time series of satellite observations at moderate (30 m) spatial resolution is enabling unprecedented opportunities for understanding ecosystems around the world. A time series of data from Landsat was used to generate a series of three maps at decadal time step to show how surface water has changed from 1991 to 2011 in the high northern latitudes of North America. Previous attempts to characterize the change in surface water in this region have been limited in either spatial or temporal resolution, or both. This series of maps was generated for the NASA Arctic and Boreal Vulnerability Experiment (ABoVE), which began in fall 2015. These maps show a “nominal” extent of surface water by using multiple observations to make a single map for each time step. This increases the confidence that any detected changes are related to climate or ecosystem changes not simply caused by short duration weather events such as flood or drought. The methods and comparison to other contemporary maps of the region are presented here. Initial verification results indicate 96% producer accuracy and 54% user accuracy when compared to 2-m resolution WorldView-2 data. All water bodies that were omitted were one Landsat pixel or smaller, hence below detection limits of the instrument.

Description: On the example of the Epembe carbonatite-hosted Nb-Ta-LREE deposit, we demonstrate the use of hyperspectral reflectance data and geomorphic indicators for improving the accuracy of remote sensing exploration data of structurally-controlled critical raw material deposits. The results further show how exploration can benefit from a combination of expert knowledge and remotely-sensed relief, as well as imaging data. In the first stage, multi-source remote sensing data were used in lithological mapping based on Kohonen Self-Organizing Maps (SOM). We exemplify that morphological indices, such as Topographic Position Index (TPI), and spatial coordinates are crucial parameters to improve the accuracy of carbonate classification as much as 10%. The resulting lithological map shows the spatial distribution of the ridge forming carbonatite dyke, the fenitization zone, syenite plugs and mafic intrusions. In a second step, the internal zones of the carbonatite complex were identified using the Multi-Range Spectral Feature Fitting (MRSFF) algorithm and a specific decision tree. This approach allowed detecting potential enrichment zones characterized by an abundance of fluorapatite and pyroxene, as well as dolomite-carbonatite (beforsite). Cross-validation of the mineral map with field observations and radiometric data confirms the accuracy of the proposed method.

Description: Landslide monitoring is important for geological disaster prevention, where Synthetic Aperture Radar (SAR) images have been widely used. Compared with the Interferometric SAR (InSAR) technique, intensity-based offset tracking methods (e.g., Normalized Cross-Correlation method) can overcome the limitation of InSAR’s maximum detectable displacement. The normalized cross-correlation (NCC) method, based on single-channel SAR images, estimates azimuth and range displacement by using statistical correlation between the matching windows of two SAR images. However, the matching windows—especially for the boundary area of landslide—always contain pixels with different moving characteristics, affecting the precision of displacement estimation. Based on the advantages of polarimetric scattering properties, this paper proposes a fully polarimetric SAR (PolSAR) offset tracking method for improvement of the precision of landslide displacement estimation. The proposed method uses the normalized inner product (NIP) of the two temporal PolSAR Pauli scattering vectors to evaluate their similarity, then retrieve the surface displacement of the Slumgullion landslide located in southwestern Colorado, USA. A pair of L-band fully polarimetric SAR images acquired by the Jet Propulsion Laboratory’s Uninhabited Aerial Vehicle Synthetic Aperture Radar (UAVSAR) system are selected for experiment. The results show that the Slumgullion landslide’s moving velocity during the monitoring time ranges between 1.6–10.9 mm/d, with an average velocity of 6.3 mm/d. Compared with the classical NCC method, results of the proposed method present better performance in the sub-pixel estimation. Furthermore, it performs better when estimating displacement in the area around the landslide boundaries.

Description: The savannas of Southern Africa are an important dryland ecosystem as they cover up to 54% of the landscape and support a rich variety of biodiversity. This paper evaluates landscape change in savanna vegetation along Chobe Riverfront within Chobe National Park Botswana, from 1982 to 2011 to understand what change may be occurring in land cover. Classifying land cover in savanna environments is challenging because the vegetation spectral signatures are similar across distinct vegetation covers. With vegetation species and even structural groups having similar signatures in multispectral imagery difficulties exist in making discrete classifications in such landscapes. To address this issue, a Random Forest classification algorithm was applied to predict land-cover classes. Additionally, time series vegetation indices were used to support the findings of the discrete land cover classification. Results indicate that a landscape level vegetation shift has occurred across the Chobe Riverfront, with results highlighting a shift in land cover towards more woody vegetation. This represents a degradation of vegetation cover within this savanna landscape environment, largely due to an increasing number of elephants and other herbivores utilizing the Riverfront. The forested area along roads at a further distance from the River has also had a loss of percent cover. The continuous analysis during 1982–2011, utilizing monthly AVHRR (Advanced Very High Resolution Radiometer) NDVI (Normalized Difference Vegetation Index) values, also verifies this change in amount of vegetation is a continuous and ongoing process in this region. This study provides land use planners and managers with a more reliable, efficient and relatively inexpensive tool for analyzing land-cover change across these highly sensitive regions, and highlights the usefulness of a Random Forest classification in conjunction with time series analysis for monitoring savanna landscapes.

Description: This study is aimed at a better understanding of how upstream runoff formation affected the cropping intensity (CI: number of harvests) in the Aral Sea Basin (ASB) between 2000 and 2012. MODIS 250 m NDVI time series and knowledge-based pixel masking that included settlement layers and topography features enabled to map the irrigated cropland extent (iCE). Random forest models supported the classification of cropland vegetation phenology (CVP: winter/summer crops, double cropping, etc.). CI and the percentage of fallow cropland (PF) were derived from CVP. Spearman’s rho was selected for assessing the statistical relation of CI and PF to runoff formation in the Amu Darya and Syr Darya catchments per hydrological year. Validation in 12 reference sites using multi-annual Landsat-7 ETM+ images revealed an average overall accuracy of 0.85 for the iCE maps. MODIS maps overestimated that based on Landsat by an average factor of ~1.15 (MODIS iCE/Landsat iCE). Exceptional overestimations occurred in case of inaccurate settlement layers. The CVP and CI maps achieved overall accuracies of 0.91 and 0.96, respectively. The Amu Darya catchment disclosed significant positive (negative) relations between upstream runoff with CI (PF) and a high pressure on the river water resources in 2000–2012. Along the Syr Darya, reduced dependencies could be observed, which is potentially linked to the high number of water constructions in that catchment. Intensified double cropping after drought years occurred in Uzbekistan. However, a 10 km × 10 km grid of Spearman’s rho (CI and PF vs. upstream runoff) emphasized locations at different CI levels that are directly affected by runoff fluctuations in both river systems. The resulting maps may thus be supportive on the way to achieve long-term sustainability of crop production and to simultaneously protect the severely threatened environment in the ASB. The gained knowledge can be further used for investigating climatic impacts of irrigation in the region.

Description: The spatial separation of the transmitter and receiver in Bistatic Synthetic Aperture Radar (BiSAR) makes it a promising and useful supplement to a classical Monostatic SAR system (MonoSAR). This paper proposes a novel integrated time and phase synchronization strategy for a multichannel spaceborne-stationary BiSAR system. Firstly, the time synchronization strategy is proposed, which includes Pulse Repetition Frequency (PRF) generation under noisy conditions, multichannel calibration and the alignment of the recorded data with the orbital data. Furthermore, the phase synchronization strategy, which fully considers the deteriorative factors in the BiSAR configuration, is well studied. The contribution of the phase synchronization strategy includes two aspects: it not only compensates the phase error, but also improves the Signal to Noise Ratio (SNR) of the obtained signals. Specifically, all direct signals on different PRF time can be reconstructed with the shift and phase compensation operation using a reference signal. Besides, since the parameters of the reference signal can be estimated only once using the selected practical direct signal and a priori information, the processing complexity is well reduced. Final imaging results with and without compensation for real data are presented to validate the proposed synchronization strategy.

Description: Mangroves are ecologically important ecosystems and globally protected. The purpose of this study was to evaluate the effectiveness of mangrove conservation efforts in two adjacent protected areas in China that were under the management policies of the Ramsar Convention (Mai Po Marshes Nature Reserve (MPMNR), Hong Kong) and China’s National Nature Reserve System (Futian Mangrove National Nature Reserve (FMNNR), Shenzhen). To achieve this goal, eleven Landsat images were chosen and classified, areal extent and landscape metrics were then calculated. The results showed that: from 1973–2015, the areal extent of mangroves in both reserves increased, but the net change for the MPMNR (281.43 hm2) was much higher than those of the FMNNR (101.97 hm2). In general, the area-weighted centroid of the mangroves in FMNNR moved seaward by approximately 120 m, whereas in the MPMNR, the centroid moved seaward even farther (410 m). Although both reserves saw increased integrality and connectivity of the mangrove patches, the patches in the MPMNR always had higher integrality than those in the FMNNR. We concluded that the mangroves in the MPMNR were more effectively protected than those in the FMNNR. This study may provide assistance to the formulation of generally accepted criteria for remote sensing-based evaluation of conservation effectiveness, and may facilitate the development of appropriate mangrove forest conservation and management strategies in other counties.

Description: Monitoring the dynamics of surface water using remotely sensed data generally requires both high spatial and high temporal resolutions. One effective and popular approach for achieving this is image fusion. This study adopts a widely accepted fusion model, the Enhanced Spatial and Temporal Adaptive Reflectance Fusion Model (ESTARFM), for blending the newly available coarse-resolution Suomi NPP-VIIRS data with Landsat data in order to derive water maps at 30 m resolution. The Pan-sharpening technique was applied to preprocessing NPP-VIIRS data to achieve a higher-resolution before blending. The modified Normalized Difference Water Index (mNDWI) was employed for mapping surface water area. Two fusion alternatives, blend-then-index (BI) or index-then-blend (IB), were comparatively analyzed against a Landsat derived water map. A case study of mapping Poyang Lake in China, where water distribution pattern is complex and the water body changes frequently and drastically, was conducted. It has been revealed that the IB method derives more accurate results with less computation time than the BI method. The BI method generally underestimates water distribution, especially when the water area expands radically. The study has demonstrated the feasibility of blending NPP-VIIRS with Landsat for achieving surface water mapping at both high spatial and high temporal resolutions. It suggests that IB is superior to BI for water mapping in terms of efficiency and accuracy. The finding of this study also has important reference values for other blending works, such as image blending for vegetation cover monitoring.

Description: Leaf area index (LAI) is a key parameter used to describe vegetation structures and is widely used in ecosystem biophysical process and vegetation productivity models. Many algorithms have been developed for the estimation of LAI based on remote sensing images. Our goal was to produce accurate and timely predictions of grassland LAI for the meadow steppes of northern China. Here, we compare the predictive power of regression approaches and hybrid geostatistical methods using Chinese Huanjing (HJ) satellite charge coupled device (CCD) data. The regression methods evaluated include partial least squares regression (PLSR), artificial neural networks (ANNs) and random forests (RFs). The two hybrid geostatistical methods were regression kriging (RK) and random forests residuals kriging (RFRK). The predictions were validated for different grassland types and different growing stages, and their performances were also examined by adding several groups of vegetation indices (VIs). The two hybrid geostatistical models (RK and RFRK) yielded the most accurate predictions (root mean squared error (RMSE) = 0.21 m2/m2 and 0.23 m2/m2 for RK and RFRK, respectively), followed by the RF model (RMSE = 0.27 m2/m2), which was the most accurate among the regression models. These three models also exhibited the best temporal performance across the duration of the growing season. The PLSR and ANN models were less accurate (RMSE = 0.33 m2/m2 and 0.35 m2/m2 for ANN and PLSR, respectively), and the PLSR model performed the worst (exhibiting varied temporal performance and unreliable prediction accuracy that was susceptible to ground conditions). By adding VIs to the predictor variables, the predictions of the PLSR and ANN models were obviously improved (RMSE improved from 0.35 m2/m2 to 0.28 m2/m2 for PLSR and from 0.33 m2/m2 to 0.28 m2/m2 for ANN); the RF and RFRK models did not generate more accurate predictions and the performance of the RK model declined (RMSE decreased from 0.21 m2/m2 to 0.32 m2/m2).

Description: Road networks are very important features in geospatial databases. Even though high-resolution optical satellite images have already been acquired for more than a decade, tools for automated extraction of road networks from these images are still rare. One consequence of this is the need for manual interaction which, in turn, is time and cost intensive. In this paper, a multi-stage approach is proposed which integrates structural, spectral, textural, as well as contextual information of objects to extract road networks from very high resolution satellite images. Highlights of the approach are a novel linearity index employed for the discrimination of elongated road segments from other objects and customized tensor voting which is utilized to fill missing parts of the network. Experiments are carried out with different datasets. Comparison of the achieved results with the results of seven state-of-the-art methods demonstrated the efficiency of the proposed approach.

Description: Recently, sparse representation has yielded successful results in hyperspectral image (HSI) classification. In the sparse representation-based classifiers (SRCs), a more discriminative representation that preserves the spectral-spatial information can be exploited by treating the HSI as a whole entity. Based on this observation, a tensor block-sparsity based representation method is proposed for spectral-spatial classification of HSI in this paper. Unlike traditional vector/matrix-based SRCs, the proposed method consists of tensor block-sparsity based dictionary learning and class-dependent block sparse representation. By naturally regarding the HSI cube as a third-order tensor, small local patches centered at the training samples are extracted from the HSI to maintain the structural information. All the patches are then partitioned into a number of groups, on which a dictionary learning model is constructed with a tensor block-sparsity constraint. A test sample is also expressed as a small local patch and the block sparse representation is then performed in a class-wise manner to take advantage of the class label information. Finally, the category of the test sample is determined by using the minimal residual. Experimental results of two real-world HSIs show that our proposed method greatly improves the classification performance of SRC.

Description: With the advance of the World-Wide Web (WWW) technology, people can easily share content on the Web, including geospatial data and web services. Thus, the “big geospatial data management” issues start attracting attention. Among the big geospatial data issues, this research focuses on discovering distributed geospatial resources. As resources are scattered on the WWW, users cannot find resources of their interests efficiently. While the WWW has Web search engines addressing web resource discovery issues, we envision that the geospatial Web (i.e., GeoWeb) also requires GeoWeb search engines. To realize a GeoWeb search engine, one of the first steps is to proactively discover GeoWeb resources on the WWW. Hence, in this study, we propose the GeoWeb Crawler, an extensible Web crawling framework that can find various types of GeoWeb resources, such as Open Geospatial Consortium (OGC) web services, Keyhole Markup Language (KML) and Environmental Systems Research Institute, Inc (ESRI) Shapefiles. In addition, we apply the distributed computing concept to promote the performance of the GeoWeb Crawler. The result shows that for 10 targeted resources types, the GeoWeb Crawler discovered 7351 geospatial services and 194,003 datasets. As a result, the proposed GeoWeb Crawler framework is proven to be extensible and scalable to provide a comprehensive index of GeoWeb.

Description: Water bodies are a fundamental element of urban ecosystems, and water mapping is critical for urban and landscape planning and management. Remote sensing has increasingly been used for water mapping in rural areas; however, when applied to urban areas, this spatially- explicit approach is a challenging task due to the fact that the water bodies are often of a small size and spectral confusion is common between water and the complex features in the urban environment. Water indexes are the most common method of water extraction at the pixel level. More recently, spectral mixture analysis (SMA) has been widely employed in analyzing the urban environment at the subpixel level. The objective of this study is to develop an automatic subpixel water mapping method (ASWM) which can achieve a high accuracy in urban areas. Specifically, we first apply a water index for the automatic extraction of mixed land-water pixels, and the pure water pixels that are generated in this process are exported as the final result. Secondly, the SMA technique is applied to the mixed land-water pixels for water abundance estimation. As for obtaining the most representative endmembers, we propose an adaptive iterative endmember selection method based on the spatial similarity of adjacent ground surfaces. One classical water index method (the modified normalized difference water index (MNDWI)), a pixel-level target detection method (constrained energy minimization (CEM)), and two widely used SMA methods (fully constrained least squares (FCLS) and multiple endmember spectral mixture analysis (MESMA)) were chosen for the water mapping comparison in the experiments. The results indicate that the proposed ASWM was able to detect water pixels more efficiency than other unsupervised water extraction methods, and the water fractions estimated by the proposed ASWM method correspond closely to the reference fractions with the slopes of 0.97, 1.02, 1.04, and 0.98 and the R-squared values of 0.9454, 0.9486, 0.9665, and 0.9607 in regression analysis corresponding to different test regions. In the quantitative accuracy assessment, the ASWM method shows the best performance in water mapping with the mean kappa coefficient of 0.862, mean producer’s accuracy of 82.8%, and mean user’s accuracy of 91.8% for test regions.

Description: Alpine swamp meadow on the Tibetan Plateau is among the most sensitive areas to climate change. Accurate quantification of the GPP in alpine swamp meadow can benefit our understanding of the global carbon cycle. The 8-day MODerate resolution Imaging Spectroradiometer (MODIS) gross primary production (GPP) products (GPP_MOD) provide a pathway to estimate GPP in this remote ecosystem. However, the accuracy of the GPP_MOD estimation in this representative alpine swamp meadow is still unknown. Here five years GPP_MOD was validated using GPP derived from the eddy covariance flux measurements (GPP_EC) from 2009 to 2013. Our results indicated that the GPP_EC was strongly underestimated by GPP_MOD with a daily mean less than 40% of EC measurements. To reduce this error, the ground meteorological and vegetation leaf area index (LAIG) measurements were used to revise the key inputs, the maximum light use efficiency (εmax) and the fractional photosynthetically active radiation (FPARM) in the MOD17 algorithm. Using two approaches to determine the site-specific εmax value, we suggested that the suitable εmax was about 1.61 g C MJ−1 for this alpine swamp meadow which was considerably larger than the default 0.68 g C MJ−1 for grassland. The FPARM underestimated 22.2% of the actual FPAR (FPARG) simulated from the LAIG during the whole study period. Model comparisons showed that the large inaccuracies of GPP_MOD were mainly caused by the underestimation of the εmax and followed by that of the undervalued FPAR. However, the DAO meteorology data in the MOD17 algorithm did not exert a significant affection in the MODIS GPP underestimations. Therefore, site-specific optimized parameters inputs, especially the εmax and FPARG, are necessary to improve the performance of the MOD17 algorithm in GPP estimation, in which the calibrated MOD17A2 algorithm (GPP_MODR3) could explain 91.6% of GPP_EC variance for the alpine swamp meadow.

Description: REDD+ implementation requires robust, consistent, accurate and transparent national land cover historical data and monitoring systems. Satellite imagery is the only data source with enough periodicity to provide consistent land cover information in a cost-effective way. The main aim of this paper is the creation of an operational framework for monitoring land cover dynamics based on Landsat imagery and open-source software. The methodology integrates the entire land cover and land cover change mapping processes to produce a consistent series of Land Cover maps. The consistency of the time series is achieved through the application of a single trained machine learning algorithm to radiometrically normalized imagery using iteratively re-weighted multivariate alteration detection (IR-MAD) across all dates of the historical period. As a result, seven individual Land Cover maps of Costa Rica were produced from 1985/1986 to 2013/2014. Post-classification land cover change detection was performed to evaluate the land cover dynamics in Costa Rica. The validation of the land cover maps showed an overall accuracy of 87% for the 2013/2014 map, 93% for the 2000/2001 map and 89% for the 1985/1986 map. Land cover changes between forest and non-forest classes were validated for the period between 2001 and 2011, obtaining an overall accuracy of 86%. Forest age-classes were generated through a multi-temporal analysis of the maps. By linking deforestation dynamics with forest age, a more accurate discussion of the carbon emissions along the time series can be presented.

Description: This paper addresses how to manage planar spatial data using MongoDB, a popular NoSQL database characterized as a document-oriented, rich query language and high availability. The core idea is to flatten a hierarchical R-tree structure into a tabular MongoDB collection, during which R-tree nodes are represented as collection documents and R-tree pointers are expressed as document identifiers. By following this strategy, a storage schema to support R-tree-based create, read, update, and delete (CRUD) operations is designed and a module to manage planar spatial data by consuming and maintaining flattened R-tree structure is developed. The R-tree module is then seamlessly integrated into MongoDB, so that users could manipulate planar spatial data with existing command interfaces oriented to geodetic spatial data. The experimental evaluation, using real-world datasets with diverse coverage, types, and sizes, shows that planar spatial data can be effectively managed by MongoDB with our flattened R-tree and, therefore, the application extent of MongoDB will be greatly enlarged. Our work resulted in a MongoDB branch with R-tree support, which has been released on GitHub for open access.

Description: Conservation of forests outside protected areas is essential for maintaining forest connectivity, which largely depends on the effectiveness of local institutions. In this study, we use Landsat data to explore the relationship between vegetation structure and forest management institutions, in order to assess the efficacy of local institutions in management of forests outside protected areas. These forests form part of an important tiger corridor in Eastern Maharashtra, India. We assessed forest condition using 450 randomly placed 10 m radius circular plots in forest patches of villages with and without local institutions, to understand the impact of these institutions on forest vegetation. Tree density and species richness were significantly different between villages with and without local forest institutions, but there was no difference in tree biomass. We also found a significant difference in the relationship between tree density and NDVI between villages with and without local forest institutions. However, the relationship between species richness and NDVI did not differ significantly. The methods proposed by this study evaluate the status of forest management in a forest corridor using remotely sensed data and could be effectively used to identify the extent of vegetation health and management status.

Description: Wetlands are important ecosystems playing an essential role for continental water regulation and the hydrologic cycle. Moreover, they are sensitive to climate changes as well as anthropogenic influences, such as land-use or dams. However, the monitoring of these regions is challenging as they are normally located in remote areas without in situ measurement stations. Radar altimetry provides important measurements for monitoring and analyzing water level variations in wetlands and flooded areas. Using the example of the Pantanal region in South America, this study demonstrates the capability and limitations of ENVISAT radar altimeter for monitoring water levels in inundation areas. By applying an innovative processing method consisting of a rigorous data screening by means of radar echo classification as well as an optimized waveform retracking, water level time series with respect to a global reference and with a temporal resolution of about one month are derived. A comparison between altimetry-derived height variations and six in situ time series reveals accuracies of 30 to 50 cm RMS. The derived water level time series document seasonal height variations of up to 1.5 m amplitude with maximum water levels between January and June. Large scale geographical pattern of water heights are visible within the wetland. However, some regions of the Pantanal show water level variations less than a few decimeter, which is below the accuracies of the method. These areas cannot be reliably monitored by ENVISAT.

Description: A new pansharpening method is proposed, based on convolutional neural networks. We adapt a simple and effective three-layer architecture recently proposed for super-resolution to the pansharpening problem. Moreover, to improve performance without increasing complexity, we augment the input by including several maps of nonlinear radiometric indices typical of remote sensing. Experiments on three representative datasets show the proposed method to provide very promising results, largely competitive with the current state of the art in terms of both full-reference and no-reference metrics, and also at a visual inspection.

Description: In the context of climate change, agricultural managers have the imperative to combine sufficient productivity with durability of the resources. Many studies have shown the interest of recent satellite missions as suitable tools for agricultural surveys. Nevertheless, they are not predictive methods. A system able to detect summer crops as early as possible is important in order to obtain valuable information for a better water management strategy. The detection of summer crops before the beginning of the irrigation period is therefore our objective. The study area is located near Toulouse (southwestern France), and is a region of mixed farming with a wide variety of irrigated and non-irrigated crops. Using the reference data for the years concerned, a set of fixed thresholds are applied to a vegetation index (the Normalized Difference Vegetation Index, NDVI) for each agricultural season of multi-spectral satellite optical imagery acquired at decametric spatial resolutions from 2006 to 2013. The performance (i.e., accuracy) is contrasted according to the agricultural practices, the development states of the different crops and the number of acquisition dates (one to three in the results presented here). The detection of summer crops reaches 64% to 88% with a single date, 80% to 88% with two dates and 90% to 99% with three dates. The robustness of this method is tested for several years (showing an impact of meteorological conditions on the actual choice of images), several sensors and several resolutions.

Description: Zebra crossings provide guidance and warning to pedestrians and drivers, thereby playing an important role in traffic safety management. Most previous studies have focused on detecting zebra stripes but have not provided full information about the areas, which is critical to both driver assistance systems and guide systems for blind individuals. This paper presents a stepwise procedure for recognizing and reconstructing zebra crossings using mobile laser scanning data. First, we propose adaptive thresholding based on road surface partitioning to reduce the impact of intensity unevenness and improve the accuracy of road marking extraction. Then, dispersion degree filtering is used to reduce the noise. Finally, zebra stripes are recognized according to the rectangular feature and fixed size, which is followed by area reconstruction according to arrangement patterns. We test our method on three datasets captured by an Optech Lynx mobile mapping system. The total recognition rate of 90.91% demonstrates the effectiveness of the method.

Description: In recent years, with the reforms to the land use system and the development of urbanization in China, land price evaluation has tended towards marketization. Prices are determined by the government, the land transaction market and the public. It is necessary to propose higher standards to be used in the evaluation process. This paper presents an online land price evaluation approach for convenience in evaluation. In a network environment, taking advantage of the data services provided by various departments, we propose two models to assist in decision-making: (1) a geographic information system (GIS)- and fuzzy set-based location factor quantification model, which adopts dynamic data, rules and quantification measures (based on the road network) to dynamically quantify location factors, thus transforming fuzzy sets into appropriate values; and (2) a neartude-based transaction sample push model, which quantifies the similarity between a given land and other samples, thus providing a basis for decision-making by an appraiser. This approach is applied in Shenzhen to evaluate its ability to simplify the work of appraisers and make their decisions more intuitive and objective in a real case.

Description: The archiving of scientific data is a sophisticated mission in nearly all research projects. In this paper, we introduce a new online archive of atmospheric measurement data from the "High definition clouds and precipitation for advancing climate prediction" (HD(CP)2) research initiative. The project data archive is quality managed, easy to use, and is now open for other atmospheric research data. The archive’s creation was already taken into account during the HD(CP)2 project planning phase and the necessary resources were granted. The funding enabled the HD(CP)2 project to build a sound archive structure, which guarantees that the collected data are accessible for all researchers in the project and beyond.

Description: Radar backscatter from forest canopies is related to forest cover, canopy structure and aboveground biomass (AGB). The S-band frequency (3.1–3.3 GHz) lies between the longer L-band (1–2 GHz) and the shorter C-band (5–6 GHz) and has been insufficiently studied for forest applications due to limited data availability. In anticipation of the British built NovaSAR-S satellite mission, this study evaluates the benefits of polarimetric S-band SAR for forest biophysical properties. To understand the scattering mechanisms in forest canopies at S-band the Michigan Microwave Canopy Scattering (MIMICS-I) radiative transfer model was used. S-band backscatter was found to have high sensitivity to the forest canopy characteristics across all polarisations and incidence angles. This sensitivity originates from ground/trunk interaction as the dominant scattering mechanism related to broadleaved species for co-polarised mode and specific incidence angles. The study was carried out in the temperate mixed forest at Savernake Forest and Wytham Woods in southern England, where airborne S-band SAR imagery and field data are available from the recent AirSAR campaign. Field data from the test sites revealed wide ranges of forest parameters, including average canopy height (6–23 m), diameter at breast-height (7–42 cm), basal area (0.2–56 m2/ha), stem density (20–350 trees/ha) and woody biomass density (31–520 t/ha). S-band backscatter-biomass relationships suggest increasing backscatter sensitivity to forest AGB with least error between 90.63 and 99.39 t/ha and coefficient of determination (r2) between 0.42 and 0.47 for the co-polarised channel at 0.25 ha resolution. The conclusion is that S-band SAR data such as from NovaSAR-S is suitable for monitoring forest aboveground biomass less than 100 t/ha at 25 m resolution in low to medium incidence angle range.

Description: The absolute radiometric calibration of a satellite sensor is the critical factor that ensures the usefulness of the acquired data for quantitative applications on remote sensing. This work presents the results of the first cross-calibration of the sensor on board the Sistema Satelital de Observación de la Tierra (SSOT) Chilean satellite or Air Force Satellite FASat-C. RapidEye-MSI was chosen as the reference sensor, and a simultaneous Nadir Overpass Approach (SNO) was applied. The biases caused by differences in the spectral responses of both instruments were compensated through an adjustment factor derived from EO-1 Hyperion data. Through this method, the variations affecting the radiometric response of New AstroSat Optical Modular Instrument (NAOMI-1), have been corrected based on collections over the Frenchman Flat calibration site. The results of a preliminary evaluation of the pre-flight and updated coefficients have shown a significant improvement in the accuracy of at-sensor radiances and TOA reflectances: an average agreement of 2.63% (RMSE) was achieved for the multispectral bands of both instruments. This research will provide a basis for the continuity of calibration and validation tasks of future Chilean space missions.

Description: Satellite-derived precipitation can be a potential source of forcing data for assessing water availability and managing water supply in mountainous regions of East Asia. This study investigates the hydrological utility of satellite-derived precipitation and uncertainties attributed to error propagation of satellite products in hydrological modeling. To this end, four satellite precipitation products (tropical rainfall measuring mission (TRMM) multi-satellite precipitation analysis (TMPA) version 6 (TMPAv6) and version 7 (TMPAv7), the global satellite mapping of precipitation (GSMaP), and the climate prediction center (CPC) morphing technique (CMORPH)) were integrated into a physically-based hydrologic model for the mountainous region of South Korea. The satellite precipitation products displayed different levels of accuracy when compared to the intra- and inter-annual variations of ground-gauged precipitation. As compared to the GSMaP and CMORPH products, superior performances were seen when the TMPA products were used within streamflow simulations. Significant dry (negative) biases in the GSMaP and CMORPH products led to large underestimates of streamflow during wet-summer seasons. Although the TMPA products displayed a good level of performance for hydrologic modeling, there were some over/underestimates of precipitation by satellites during the winter season that were induced by snow accumulation and snowmelt processes. These differences resulted in streamflow simulation uncertainties during the winter and spring seasons. This study highlights the crucial need to understand hydrological uncertainties from satellite-derived precipitation for improved water resource management and planning in mountainous basins. Furthermore, it is suggested that a reliable snowfall detection algorithm is necessary for the new global precipitation measurement (GPM) mission.

Description: In this paper, an automatic approach for zebra crossing extraction and reconstruction from high-resolution aerial images is proposed. In the extraction procedure, zebra crossings are extracted by the JointBoost classifier based on GLCM (Gray Level Co-occurrence Matrix) features and 2D Gabor Features. In the reconstruction procedure, a geometric parameter model based on spatial repeatability relationships is globally fitted to reconstruct the geometric shape of zebra crossings. Additionally, a group of representative experiments is conducted to test the proposed method under interfered conditions, such as zebra crossings covered by pedestrians, shadows and color fading. Furthermore, the performance of the proposed extraction method is compared with the template matching method. Finally, the results show the validation of our proposed method, both in the extraction and reconstruction of zebra crossings.

Description: Simplification of three-dimensional (3D) buildings is critical to improve the efficiency of visualizing urban environments while ensuring realistic urban scenes. Moreover, it underpins the construction of multi-scale 3D city models (3DCMs) which could be applied to study various urban issues. In this paper, we design a generic yet effective approach for simplifying 3D buildings. Instead of relying on both semantic information and geometric information, our approach is based solely on geometric information as many 3D buildings still do not include semantic information. In addition, it provides an integrated means to treat 3D buildings with either sloped or flat roofs. The two case studies, one exploring simplification of individual 3D buildings at varying levels of complexity while the other, investigating the multi-scale simplification of a cityscape, show the effectiveness of our approach.

Description: Ukraine has experienced immense environmental and institutional changes during the last three decades. We have conducted this study to analyze important land surface dynamics and to assess processes underlying the changes. This research was conducted in two consecutive steps. To analyze monotonic changes we first applied a Mann–Kendall trend analysis of the Normalized Difference Vegetation Index (NDVI3g) time series. Gradual and abrupt changes were studied by fitting a seasonal trend model and detecting the breakpoints. Secondly, essential environmental factors were used to quantify their possible relationships with land surface changes. These factors included soil moisture as well as gridded air temperature and precipitation data. This was done using partial rank correlation analysis based on annually aggregated time-series. Our results demonstrate that positive NDVI trends characterize approximately one-third of Ukraine’s land surface, located in the northern and western areas of the country. Negative trends occurred less frequently, covering less than 2% of the area and are distributed irregularly across the country. Monotonic trends were rarely found; shifting trends were identified with a greater frequency. Trend shifts were seen to occur with an increased frequency following the period of the 2000s. We determined that land surface dynamics and climate variability are functionally interdependent; however, the relative influence of the drivers varies in different locations. Among the factors analyzed, the air temperature variable explains the largest portion of NDVI variability. High air temperature/NDVI correlation coefficients (r = 0.36 − 0.77) are observed over the entire country. The soil moisture content is of significant influence in the eastern portion of Ukraine (r = 0.68); precipitation (r = 0.65) was most influential in the central regions of the country. These results increase our understanding of ecosystem responses to climatic changes and anthropogenic activities.

Description: The Biosphere Reserve of La Mancha Húmeda is a wetland-rich area located in central Spain. This reserve comprises a set of temporary lakes, often saline, where water level fluctuates seasonally. Water inflows come mainly from direct precipitation and runoff of small lake watersheds. Most of these lakes lack surface outlets and behave as endorheic systems, where water withdrawal is mainly due to evaporation, causing salt accumulation in the lake beds. Remote sensing was used to estimate the temporal variation of the flooded area in these lakes and their associated hydrological patterns related to the seasonality of precipitation and evapotranspiration. Landsat 7 ETM+ satellite images for the reference period 2013–2015 were jointly used with ground-truth datasets. Several inverse modeling methods, such as two-band and multispectral indices, single-band threshold, classification methods, artificial neural network, support vector machine and genetic programming, were applied to retrieve information on the variation of the flooded areas. Results were compared to ground-truth data, and the classification errors were evaluated by means of the kappa coefficient. Comparative analyses demonstrated that the genetic programming approach yielded the best results, with a kappa value of 0.98 and a total error of omission-commission of 2%. The dependence of the variations in the water-covered area on precipitation and evaporation was also investigated. The results show the potential of the tested techniques to monitor the hydrological patterns of temporary lakes in semiarid areas, which might be useful for management strategy-linked lake conservation and specifically to accomplish the goals of both the European Water Framework Directive and the Habitats Directive.

Description: Spatially and temporally resolved observations of near-surface air temperatures (Ta, 1.5–2 m above ground) are essential for understanding hydrothermal circulation at the land–atmosphere interface. However, the uneven spatial distribution of meteorological stations may not effectively capture the true nature of the overall climate pattern. Several studies have attempted to retrieve spatially continuous Ta from remotely sensed and continuously monitored Land Surface Temperature (LST). However, the topographical control of the relationship between LST and Ta in regions with complex topographies and highly variable weather station densities is poorly understood. The aim of this study is to improve the accuracy of Ta estimations from the Moderate Resolution Imaging Spectroradiometer (MODIS) LST via parameterization of the physiographic variables according to the terrain relief. The performances of both Terra and Aqua MODIS LST in estimating Ta have been explored in China. The results indicated that the best agreement was found between Terra nighttime LST (LSTmodn) and the observed Ta in China. In flat terrain areas, the LSTmodn product is significantly linearly correlated with Ta (R2 > 0.80), while, in mountainous areas, the LSTmodn-Ta relationship differed significantly from simple linear correlation. By taking the physiographic features into account, including the seasonal vegetation cover (NDVI), the altitudinal gradient (RDLS), and the ambient absolute humidity (AH), the accuracy of the estimation was substantially improved. The study results indicated that the relevant environmental factors must be considered when interpreting the spatiotemporal variation of the surface energy flux over complex topography.

Description: Aeolian desertification is a kind of land degradation that is characterized by aeolian activity, resulting from the responses of land ecosystems to climate change and anthropogenic disturbances. The source areas of the Yangtze and Yellow Rivers are typical regions of China’s Tibetan Plateau affected by aeolian desertification. We assessed the vulnerability of these areas to aeolian desertification by combining remote sensing with geographical information system technologies. We developed an assessment model with eight indicators, whose weights were determined by the analytical hierarchy process. Employing this model, we analyzed the spatial distribution of vulnerability to aeolian desertification and its changes from 2000 to 2010, and discuss the implications. Overall, low-vulnerability land was the most widespread, accounting for 64%, 62%, and 71% of the total study area in 2000, 2005, and 2010, respectively. The degree of vulnerability showed regional differences. In the source areas of the Yangtze River, land with high or very high vulnerability accounted for 17.4% of this sub-region in 2010, versus 2.6% in the source areas of the Yellow River. In the Zoige Basin, almost all of the land had very low to low vulnerability. To understand the change in vulnerability to aeolian desertification, we calculated an integrated vulnerability index (IVI). This analysis indicated that the vulnerability to aeolian desertification increased from 2000 to 2005 (IVI increased from 2.1709 to 2.2463), and decreased from 2005 to 2010 (IVI decreased from 2.2463 to 2.0057). Increasing regional temperatures appear to be primarily responsible for the change in vulnerability to aeolian desertification throughout the region. The effects of other factors (climatic variation and human activities) differed among the various sub-regions. The implementation of the ecological restoration project has achieved a noticeable effect since 2005. Our results provide empirical support for effort to protect the ecology of this ecologically fragile region.

Description: In order discuss the impact of land consumption, it is first necessary to localize and quantify the extent of sealed surfaces. Since 2010, the monitoring of land use structures and developments in Germany has been provided by the Monitor of Settlement and Open Space Development at the Leibniz Institute of Ecological Urban and Regional Development (IÖR; IÖR Monitor), a scientific service operated by the Leibniz Institute of Ecological Urban and Regional Development. The IÖR Monitor includes an indicator for soil sealing for the years 2006, 2009 and 2012. Using this new source of data, it is possible for the first time to conduct quantitative studies at the level of Germany’s municipalities with the aim of documenting the extent of soil sealing as a form of spatial classification, as well as to investigate possible correlations with other influential factors. Here, we describe a comprehensive data inspection of soil sealing and potential influential factors. Structural interrelationships are identified under the application of classical and spatial regression methods.

Description: The Göttingen eResearch Alliance is presented as a case study for establishing institutional support for research data management within the context of the Göttingen Campus, a particular alliance of several research institutes at Göttingen. The cross-cutting, “horizontal” approach of the Göttingen eResearch Alliance, established by two research-oriented infrastructure providers, a research library and a computing and IT competence center, aims to coordinate Campus-led activities to establish sustainable and innovative services to support all phases of the research data life cycle. In this article, the core activities of the first phase aimed at developing a modular approach to provide support for research data management to researchers will be described. It closes with lessons learned and an outlook on future activities.

Description: In recent decades, the area and proportion of planted forests have increased; thus, understanding the responses of planted and natural forests to drought are crucial because it forms the basis for forest risk assessments and management strategies. In this study, we combined the moderate-resolution imaging spectroradiometer (MODIS) enhanced vegetation index (EVI), meteorological aridity indices, and standardized precipitation evapotranspiration indices (SPEI) to identify the drought responses of planted and natural forests. In particular, we used the EVI standard anomaly (ESA) as a physiological drought indicator and analyzed the applicability of SPEIs at time scales of 1–30 months, thereby determining the optimal time scale for the SPEI (SPEIopt), i.e., the SPEI that best represents the drought responses of forests in Yunnan. Next, we employed the optimal SPEI and the ESA as indices to statistically analyze the response characteristics of planted and natural forests under different drought intensities. The results indicated the following: (1) The SPEI in June and a time scale of five months (i.e., SPEIJun,5) comprise the optimal meteorological aridity indicator for forests in Yunnan Province, which had the strongest correlation with the EVI standard anomaly (ESAJun). (2) All forest types were affected by drought in Yunnan, but their responses varied according to the forest type, elevation, and drought intensity. In general, natural forests are more vulnerable and sensitive to drought than planted forests, especially natural coniferous forests at low (0–2000 m) and moderate (2000–4000 m) altitudes, and natural mixed forest at low altitudes (0–2000 m). (3) The remote sensing-based ESA (ESAJun) is sensitive to the intensity of water stress, which makes it a good indicator for drought monitoring. In addition, the forests’ inventory survey revealed that 8.05% of forests were affected by drought; thus, we used this as a guide to estimate an approximate threshold to map forest responses to drought across the region. Below this approximate threshold (i.e., ESAJun 〈 −3.85), severe drought-induced effects on forests may occur. Given that natural forests are more vulnerable and sensitive to drought than the planted forests, natural forests need more careful management, especially in the context of projected increases in extreme drought events in the future.

Description: Traditional in situ observation interpolation techniques that provide rainfall data from rain gauges have limitations because they are discrete point-based data records, which may not be sufficient to assess droughts from a spatiotemporal perspective. Considering this limitation, this study has developed a run-off model—a fully satellite-based method for monitoring drought in Peninsular Malaysia. The formulation of the run-off deficit uses a water balance equation based on satellite-based rainfall and evapotranspiration data extracted respectively from calibrated TRMM multi-satellites precipitation analysis data (TMPA) and moderate resolution imaging spectroradiometer data (MODIS). The run-off deficit was calculated based on per pixel spatial scale and allowed to produce the continuous and regular run-off maps. The run-off model was tested and evaluated in a one drought year (2005) within a span of three years (2003–2005) over the Kelantan (3448 km2) and Hulu Perak (3672 km2) catchments of Peninsular Malaysia. The validation results show that (1) monthly TMPA rainfall and MODIS evapotranspiration data significantly improved after calibration; (2) satellite-based run-off data is not only strongly correlated with actual steam flow, but also with spatiotemporal variation of run-off in drought-affected forest catchments. The most severely drought-affected forest catchments that experienced the run-off deficits were Hulu Perak, Ulu Gading, Gunung Stong and Relai over Kelantan. The real time run-off change analysis shows that drought started in January and reached its peak in July of 2005. It was therefore demonstrated that this fully satellite-based run-off deficit model is as good as a conventional drought-monitoring indicator, and can provide not only drought distribution information, but it also can reflect the drought-induced impacts on stream flow, forest catchment and land-use.

Description: Urbanization often occurs in an unplanned and uneven manner, resulting in profound changes in patterns of land cover and land use. Understanding these changes is fundamental for devising environmentally responsible approaches to economic development in the rapidly urbanizing countries of the emerging world. One indicator of urbanization is built-up land cover that can be detected and quantified at scale using satellite imagery and cloud-based computational platforms. This process requires reliable and comprehensive ground-truth data for supervised classification and for validation of classification products. We present a new dataset for India, consisting of 21,030 polygons from across the country that were manually classified as “built-up” or “not built-up,” which we use for supervised image classification and detection of urban areas. As a large and geographically diverse country that has been undergoing an urban transition, India represents an ideal context to develop and test approaches for the detection of features related to urbanization. We perform the analysis in Google Earth Engine (GEE) using three types of classifiers, based on imagery from Landsat 7 and Landsat 8 as inputs. The methodology produces high-quality maps of built-up areas across space and time. Although the dataset can facilitate supervised image classification in any platform, we highlight its potential use in GEE for temporal large-scale analysis of the urbanization process. Our methodology can easily be applied to other countries and regions.

Description: The implementation of sustainable land policies is in need of monitoring methods that go beyond a mere description of the proportion values of land use classes. The annual statistical surface area report on actual land utilization (German: “Bodenfläche nach Art der tatsächlichen Nutzung”), published by the statistical offices of the German federal states and the federation, provides information on a set of pre-defined land use classes for municipalities, districts and federal states. Due to its surveying method of summing up usage information from cadastral registers, it is not possible to determine previous and subsequent usages of land parcels. Hence, it is hard to precisely indicate to what extent particular land use classes contribute to the settlement area increase. Nevertheless, this information is crucial to the understanding of land use change processes, which is needed for a subsequent identification of driving forces. To overcome this lack of information, a method for the spatial and quantitative determination of previous and subsequent land usages has been developed, implemented and tested. It is based on pre-processed land use data for different time slices, which are derived from authoritative geo-topographical base data. The developed method allows for the identification of land use changes considering small geometric shifts and changes in the underlying data model, which can be adaptively excluded from the balance.

Description: Surface urban heat island (SUHI) is defined as the elevated land surface temperature (LST) in urban area in comparison with non-urban areas, and it can influence the energy consumption, comfort and health of urban residents. In this study, the existence of daytime SUHI, in Cairo and its new towns during the summer, is investigated using three different approaches; (1) utilization of pre-urbanization observations as LST references; (2) utilization of rural observations as LST references (urban–rural difference); and (3) utilization of the SIUHI (Surface Intra Urban Heat Island) approach. A time series of Landsat TM & ETM+ data (46 images) from 1984 to 2015 was employed in this study for daytime LST calculation during summer. Different statistical hypothesis tests were utilized for the evaluation of LST and SUHI in the case studies. The results demonstrated that there is no significant LST difference between the urban areas studied, and their corresponding built-up areas. In addition, daytime LST in new towns during the summer is 2 K warmer than in Cairo. Utilization of a pre-urbanization observations approach, alongside an evaluation of the long-term trend, demonstrated that there is no daytime SUHI during the summer in the study areas, and construction activities in the study areas do not result in cooling or warming effects. Utilization of the rural observations approach showed that LST is lower in Cairo than its surrounding areas. This demonstrates why the selection of suitable rural references in SUHI studies is an important and complicated task, and how this approach may lead to misinterpretation in desert city areas with significant landscape and surface difference with their most surrounding areas (e.g., Cairo). Results showed that, although SIUHI technique can be representative for the changes of variance of LST in urban areas, it is not able to identify the changes of mean LST in urban areas.

Description: Three-dimensional information provided by TanDEM-X interferometric phase and airborne Light Detection and Ranging (LiDAR) Digital Elevation Models (DEMs) were used to detect differences in vegetation heterogeneity through a disturbance gradient in Indonesia. The range of vegetation types developed as a consequence of fires during the 1997–1998 El Niño. Two-point statistic (wavelet variance and co-variance) was used to assess the dominant spatial frequencies associated with either topographic features or canopy structure. DEMs wavelet spectra were found to be sensitive to canopy structure at short scales (up to 8 m) but increasingly influenced by topographic structures at longer scales. Analysis also indicates that, at short scale, canopy texture is driven by the distribution of heights. Thematic class separation using the Jeffries–Matusita distance (JM) was greater when using the full wavelet signature (LiDAR: 1.29 ≤ JM ≤ 1.39; TanDEM-X: 1.18 ≤ JM ≤ 1.39) compared to using each decomposition scale individually (LiDAR: 0.1 ≤ JM ≤ 1.26; TanDEM-X: 0.1 ≤ JM ≤ 1.1). In some cases, separability with TanDEM-X was similar to the higher resolution LiDAR. The study highlights the potential of 3D information from TanDEM-X and LiDAR DEMs to explore vegetation disturbance historywhen analyzed using two-point statistics.

Description: The importance of lakes and reservoirs leads to the high need for monitoring lake water quality both at local and global scales. The aim of the study was to test suitability of Sentinel-2 Multispectral Imager’s (MSI) data for mapping different lake water quality parameters. In situ data of chlorophyll a (Chl a), water color, colored dissolved organic matter (CDOM) and dissolved organic carbon (DOC) from nine small and two large lakes were compared with band ratio algorithms derived from Sentinel-2 Level-1C and atmospherically corrected (Sen2cor) Level-2A images. The height of the 705 nm peak was used for estimating Chl a. The suitability of the commonly used green to red band ratio was tested for estimating the CDOM, DOC and water color. Concurrent reflectance measurements were not available. Therefore, we were not able to validate the performance of Sen2cor atmospheric correction available in the Sentinel-2 Toolbox. The shape and magnitude of water reflectance were consistent with our field measurements from previous years. However, the atmospheric correction reduced the correlation between the band ratio algorithms and water quality parameters indicating the need in better atmospheric correction. We were able to show that there is good correlation between band ratio algorithms calculated from Sentinel-2 MSI data and lake water parameters like Chl a (R2 = 0.83), CDOM (R2 = 0.72) and DOC (R2 = 0.92) concentrations as well as water color (R2 = 0.52). The in situ dataset was limited in number, but covered a reasonably wide range of optical water properties. These preliminary results allow us to assume that Sentinel-2 will be a valuable tool for lake monitoring and research, especially taking into account that the data will be available routinely for many years, the imagery will be frequent, and free of charge.

Description: Evapotranspiration (ET) is an essential part of the hydrological cycle and accurately estimating it plays a crucial role in water resource management. Surface energy balance (SEB) models are widely used to estimate regional ET with remote sensing. The presence of horizontal advection, however, perturbs the surface energy balance system and contributes to the uncertainty of energy influxes. Thus, it is vital to consider horizontal advection when applying SEB models to estimate ET. This study proposes an innovative and simplified approach, the surface energy balance-advection (SEB-A) method, which is based on the energy balance theory and also takes into account the horizontal advection to determine ET by remote sensing. The SEB-A method considers that the actual ET consists of two parts: the local ET that is regulated by the energy balance system and the exotic ET that arises from horizontal advection. To evaluate the SEB-A method, it was applied to the middle region of the Heihe River in China. Instantaneous ET for three days were acquired and assessed with ET measurements from eddy covariance (EC) systems. The results demonstrated that the ET estimates had a high accuracy, with a correlation coefficient (R2) of 0.713, a mean average error (MAE) of 39.3 W/m2 and a root mean square error (RMSE) of 54.6 W/m2 between the estimates and corresponding measurements. Percent error was calculated to more rigorously assess the accuracy of these estimates, and it ranged from 0% to 35%, with over 80% of the locations within a 20% error. To better understand the SEB-A method, the relationship between the ET estimates and land use types was analyzed, and the results indicated that the ET estimates had spatial distributions that correlated with vegetation patterns and could well demonstrate the ET differences caused by different land use types. The sensitivity analysis suggested that the SEB-A method requested accurate estimation of the available energy, R n − G , but was less constrained with the difference between ground and air temperature, T 0 − T a – l o c .

Description: Nesting habitat for the federally endangered loggerhead sea turtle (Caretta caretta) were designated as critical in 2014 for beaches along the Atlantic Coast and Gulf of Mexico. Nesting suitability is routinely determined based on site specific information. Given the expansive geographic location of the designated critical C. caretta nesting habitat and the highly dynamic coastal environment, understanding nesting suitability on a regional scale is essential for monitoring the changing status of the coast as a result of hydrodynamic forces and maintenance efforts. The increasing spatial resolution and temporal frequency of remote sensing data offers the opportunity to study this dynamic environment on a regional scale. Remote sensing data were used as input into the spatially-explicit, multi-criteria decision support model to determine nesting habitat suitability. Results from the study indicate that the morphological parameters used as input into the model are well suited to provide a regional level approach with the results from the optimized model having sensitivity and detection prevalence values greater than 80% and the detection rate being greater than 70%. The approach can be implemented in various geographic locations to better communicate priorities and evaluate management strategies as a result of changes to the dynamic coastal environment.

Description: Combining sample plot and image data has been widely used to map forest carbon density at local, regional, national and global scales. When mapping is conducted using multiple spatial resolution images at different scales, field observations have to be collected at the corresponding resolutions to match image values in pixel sizes. Given a study area, however, to save time and cost, field observations are often collected from sample plots having a fixed size. This will lead to inconsistency of spatial resolutions between sample plots and image pixels and impede the mapping and product quality assessment. In this study, a methodological framework was proposed to conduct mapping and accuracy assessment of forest carbon density at four spatial resolutions by combining remotely sensed data and reference values of sample plots from a systematical, nested and clustering sampling design. This design led to one field observation dataset at a 30 m spatial resolution sample plot level and three other reference datasets by averaging the observations from three, five and seven sample plots within each of 250 m and 500 m sub-blocks and 1000 m blocks, respectively. The datasets matched the pixel values of a Landsat 8 image and three MODIS products. A sequential Gaussian co-simulation (SGCS) and a sequential Gaussian block co-simulation (SGBCS), an upscaling algorithm, were employed to map forest carbon density at the spatial resolutions. This methodology was tested for mapping forest carbon density in Huang-Feng-Qiao forest farm of You County in Eastern Hunan of China. The results showed that: First, all of the means of predicted forest carbon density values at four spatial resolutions fell in the confidence intervals of the reference data at a significance level of 0.05. Second, the systematical, nested and clustering sampling design provided the potential to obtain spatial information of forest carbon density at multiple spatial resolutions. Third, the relative root mean square error (RMSE) of predicted values at the plot level was much greater than those at the sub-block and block levels. Moreover, the accuracies of the up-scaled estimates were much higher than those from previous studies. In addition, at the same spatial resolution, SGCSWA (scaling up the SGCS and Landsat derived 30 m resolution map using a window average (WA)) resulted in smallest relative RMSEs of up-scaled predictions, followed by combinations of Landsat images and SGBCS. The accuracies from both methods were significantly greater than those from the combinations of MODIS images and SGCS. Overall, this study implied that the combinations of Landsat 8 images and SGCSWA or SGBCS with the systematical, nested and clustering sampling design provided the potential to formulate a methodological framework to map forest carbon density and conduct accuracy assessment at multiple spatial resolutions. However, this methodology needs to be further refined and examined in other forest landscapes.

Description: Sky View Factor (SVF, a dimensionless value between 0 and 1 representing obstructed and unobstructed sky, respectively) has an important influence on urban energy balance, and is a key contributor to the Urban Heat Island (UHI) effect experienced by heavily built up regions. Continuous urban SVF maps used in modeling the spatial distribution of UHI can be derived analytically using Lidar data; however, Lidar data are costly to obtain and often lack complete coverage of large cities or metropolitan areas. This study develops and validates a method for estimating continuous urban SVF from globally available Landsat TM data, based on the presence of shadows cast by SVF-reducing urban features. SVF and per-pixel shadow proportion (SP) were first calculated for synthetic grid cities to confirm a logarithmic relationship between the two properties; then Lidar data from four US cities were used to determine an empirical regression relating SP to SVF. Spectral Mixture Analysis was then used to estimate per-pixel SP in a Landsat 5 TM image covering the Greater Vancouver Area, Canada, and the empirical regression was used to calculate SVF from per-pixel SP. The accuracy of the resulting SVF map was validated using independent Lidar-derived SVF data (R2 = 0.78; RMSE = 0.056).

Description: Fire activity, in terms of intensity, frequency, and total area burned, is expected to increase with a changing climate. A challenge for landscape-level assessment of fire effects, often termed burn severity, is that current remote sensing assessments provide very little information regarding tree/vegetation physiological performance and recovery, limiting our understanding of fire effects on ecosystem services such as carbon storage/cycling. In this paper, we evaluated whether spectral indicescommoninvegetationstressandburnseverityassessmentscouldaccuratelyquantifypost-fire physiological performance (indicated by net photosynthesis and crown scorch) of two seedling species, Larix occidentalis and Pinus contorta. Seedlings were subjected to increasing fire radiative energy density (FRED) doses through a series of controlled laboratory surface fires. Mortality, physiology, and spectral reflectance were assessed for a month following the fires, and then again at one year post-fire. The differenced Normalized Difference Vegetation Index (dNDVI) spectral index outperformed other spectral indices used for vegetation stress and burn severity characterization in regard to leaf net photosynthesis quantification, indicating that landscape-level quantification of tree physiology may be possible. Additionally, the survival of the majority of seedlings in the low and moderate FRED doses indicates that fire-induced mortality is more complex than the currently accepted binary scenario, where trees survive with no impacts below a certain temperature and duration threshold, and mortality occurs above the threshold.

Description: Mapping of glacier extents from automated classification of optical satellite images has become a major application of the freely available images from Landsat. A widely applied method is based on segmented ratio images from a red and shortwave infrared band. With the now available data from Sentinel-2 (S2) and Landsat 8 (L8) there is high potential to further extend the existing time series (starting with Landsat 4/5 in 1982) and to considerably improve over previous capabilities, thanks to increased spatial resolution and dynamic range, a wider swath width and more frequent coverage. Here, we test and compare a variety of previously used methods to map glacier extents from S2 and L8, and investigate the mapping of snow facies with S2 using top of atmosphere reflectance. Our results confirm that the band ratio method works well with S2 and L8. The 15 m panchromatic band of L8 can be used instead of the red band, resulting in glacier extents similar to S2 (0.7% larger for 155 glaciers). On the other hand, extents derived from the 30 m bands are 4%–5% larger, indicating a more generous interpretation of mixed pixels. Mapping of snow cover with S2 provided accurate results, but the required topographic correction would benefit from a better orthorectification with a more precise DEM than currently used.

Description: The performance of Day-1 Integrated Multi-satellitE Retrievals for the Global Precipitation Measurement (GPM) mission (IMERG) and its predecessor, the Tropical Rainfall Measuring Mission (TRMM) Multisatellite Precipitation Analysis 3B42 Version 7 (3B42V7), was cross-evaluated using data from the best-available hourly gauge network over the Tibetan Plateau (TP). Analyses of three-hourly rainfall estimates in the warm season of 2014 reveal that IMERG shows appreciably better correlations and lower errors than 3B42V7, though with very similar spatial patterns for all assessment indicators. IMERG also appears to detect light rainfall better than 3B42V7. However, IMERG shows slightly lower POD than 3B42V7 for elevations above 4200 m. Both IMERG and 3B42V7 successfully capture the northward dynamic life cycle of the Indian monsoon reasonably well over the TP. In particular, the relatively light rain from early and end Indian monsoon moisture surge events often fails to be captured by the sparsely-distributed gauges. In spite of limited snowfall field observations, IMERG shows the potential of detecting solid precipitation, which cannot be retrieved from the 3B42V7 products.

Description: Land suitability analysis is employed to evaluate the appropriateness of land for a particular purpose whilst integrating both qualitative and quantitative inputs, which can be continuous in nature. However, in agricultural modelling there is often a disregard of this contiguous aspect. Therefore, some parametric procedures for suitability analysis compartmentalise units into defined membership classes. This imposition of crisp boundaries neglects the continuous formations found throughout nature and overlooks differences and inherent uncertainties found in the modelling. This research will compare two approaches to suitability analysis over three differing methods. The primary approach will use an Analytical Hierarchy Process (AHP), while the other approach will use a Fuzzy AHP over two methods; Fitted Fuzzy AHP and Nested Fuzzy AHP. Secondary to this, each method will be assessed into how it behaves in a climate change scenario to understand and highlight the role of uncertainties in model conceptualisation and structure. Outputs and comparisons between each method, in relation to area, proportion of membership classes and spatial representation, showed that fuzzy modelling techniques detailed a more robust and continuous output. In particular the Nested Fuzzy AHP was concluded to be more pertinent, as it incorporated complex modelling techniques, as well as the initial AHP framework. Through this comparison and assessment of model behaviour, an evaluation of each methods predictive capacity and relevance for decision-making purposes in agricultural applications is gained.

Description: The problem of heading drift error using only low cost Micro-Electro-Mechanical (MEMS) Inertial-Measurement-Unit (IMU) has not been well solved. In this paper, a heading estimation method with real-time compensation based on Kalman filter has been proposed, abbreviated as KHD. For the KHD method, a unified heading error model is established for various predictable errors in magnetic compass for pedestrian navigation, and an effective method for solving the model parameters is proposed in the indoor environment with regular structure. In addition, error model parameters are solved by Kalman filtering algorithm with building geometry information in order to achieve real-time heading compensation. The experimental results show that the KHD method can not only effectively correct the original heading information, but also effectively inhibit the accumulation effect of positioning errors. The performance observed in a field experiment performed on the fourth floor of the School of Environmental Science and Spatial Informatics (SESSI) building on the China University of Mining and Technology (CUMT) campus confirms that apply KHD method to PDR(Pedestrian Dead Reckoning) algorithm can reliably achieve meter-level positioning using a low cost MEMS IMU only.

Description: A comprehensive spectral-biogeochemical database was developed for the Wabash River and the Tippecanoe River in Indiana, United States. This database includes spectral measurements of river water, coincident in situ measurements of water quality parameters (chlorophyll (chl), non-algal particles (NAP), and colored dissolved organic matter (CDOM)), nutrients (total nitrogen (TN), total phosphorus (TP), and dissolved organic carbon (DOC)), water-column inherent optical properties (IOPs), water depths, substrate types, and bottom reflectance spectra collected in summer 2014. With this dataset, the temporal variability of water quality observations was first analyzed and studied. Second, radiative transfer models were inverted to retrieve water quality parameters using a look-up table (LUT) based spectrum matching methodology. Results found that the temporal variability of water quality parameters and nutrients in the Wabash River was closely associated with hydrologic conditions. Meanwhile, there were no significant correlations found between these parameters and streamflow for the Tippecanoe River, due to the two upstream reservoirs, which increase the settling of sediment and uptake of nutrients. The poor relationship between CDOM and DOC indicates that most DOC in the rivers was from human sources such as wastewater. It was also found that the source of water (surface runoff or combined sewer overflow (CSO)), water temperature, and nutrients were important factors controlling instream concentrations of phytoplankton. The LUT retrieved NAP concentrations were in good agreement with field measurements with slope close to 1.0 and the average estimation error was 4.1% of independently obtained lab measurements. The error for chl estimation was larger (37.7%), which is attributed to the fact that the specific absorption spectrum of chl was not well represented in this study. The LUT retrievals for CDOM experienced large variability, probably due to the small data range collected in this study and the insensitivity of Rrs to CDOM change. It is concluded that the success of the LUT method requires accurate spectral measurements and enough a priori information of the environment to construct a representative database for water quality retrieval. Therefore, future work will focus on continuing data collection in other seasons of the year and better characterization of the study area.

Description: The northernmost and most arid mangrove ecosystem of West Africa is found in Mauritania, in the Parc National du Banc d’Arguin (PNBA). The existing global and regional maps of Mauritania’s mangroves have little detail, and available estimates of the mangrove area differ among studies. We assessed the use of automated Remote Sensing classification techniques to calculate the extent and map the distribution of the mangrove patches located at Cap Timiris, PNBA, using QuickBird and GeoEye imagery. It was possible to detect the northernmost contiguous mangrove patches of West Africa with an accuracy of 87% ± 2% using the Maximum Likelihood algorithm. The main source of error was the low spectral difference between mangroves and other types of terrestrial vegetation, which resulted in an erroneous classification between these two types of land cover. The most reliable estimate for the mangrove area obtained in this study was 19.48 ± 5.54 ha in 2011. Moreover, we present a special validation procedure that enables a detailed and reliable validation of the land cover maps.

Description: Greenhouse mapping through remote sensing has received extensive attention over the last decades. In this article, the innovative goal relies on mapping greenhouses through the combined use of very high resolution satellite data (WorldView-2) and Landsat 8 Operational Land Imager (OLI) time series within a context of an object-based image analysis (OBIA) and decision tree classification. Thus, WorldView-2 was mainly used to segment the study area focusing on individual greenhouses. Basic spectral information, spectral and vegetation indices, textural features, seasonal statistics and a spectral metric (Moment Distance Index, MDI) derived from Landsat 8 time series and/or WorldView-2 imagery were computed on previously segmented image objects. In order to test its temporal stability, the same approach was applied for two different years, 2014 and 2015. In both years, MDI was pointed out as the most important feature to detect greenhouses. Moreover, the threshold value of this spectral metric turned to be extremely stable for both Landsat 8 and WorldView-2 imagery. A simple decision tree always using the same threshold values for features from Landsat 8 time series and WorldView-2 was finally proposed. Overall accuracies of 93.0% and 93.3% and kappa coefficients of 0.856 and 0.861 were attained for 2014 and 2015 datasets, respectively.

Description: Tree cover maps are used for many purposes, such as vegetation mapping, habitat connectivity and fragmentation studies. Small remnant patches of native vegetation are recognised as ecologically important, yet they are underestimated in remote sensing products derived from Landsat. High spatial resolution sensors are capable of mapping small patches of trees, but their use in large-area mapping has been limited. In this study, multi-temporal Satellite pour l’Observation de la Terre 5 (SPOT5) High Resolution Geometrical data was pan-sharpened to 5 m resolution and used to map tree cover for the Australian state of New South Wales (NSW), an area of over 800,000 km2. Complete coverages of SPOT5 panchromatic and multispectral data over NSW were acquired during four consecutive summers (2008–2011) for a total of 1256 images. After pre-processing, the imagery was used to model foliage projective cover (FPC), a measure of tree canopy density commonly used in Australia. The multi-temporal imagery, FPC models and 26,579 training pixels were used in a binomial logistic regression model to estimate the probability of each pixel containing trees. The probability images were classified into a binary map of tree cover using local thresholds, and then visually edited to reduce errors. The final tree map was then attributed with the mean FPC value from the multi-temporal imagery. Validation of the binary map based on visually assessed high resolution reference imagery revealed an overall accuracy of 88% (±0.51% standard error), while comparison against airborne lidar derived data also resulted in an overall accuracy of 88%. A preliminary assessment of the FPC map by comparing against 76 field measurements showed a very good agreement (r2 = 0.90) with a root mean square error of 8.57%, although this may not be representative due to the opportunistic sampling design. The map represents a regionally consistent and locally relevant record of tree cover for NSW, and is already widely used for natural resource management in the state.

Description: An important component of the AVHRR PATMOS-x climate date record (CDR)—or any satellite cloud climatology—is the performance of its cloud detection scheme and the subsequent quality of its cloud fraction CDR. PATMOS-x employs the NOAA Enterprise Cloud Mask for this, which is based on a naïve Bayesian approach. The goal of this paper is to generate analysis of the PATMOS-x cloud fraction CDR to facilitate its use in climate studies. Performance of PATMOS-x cloud detection is compared to that of the well-established MYD35 and CALIPSO products from the EOS A-Train. Results show the AVHRR PATMOS-x CDR compares well against CALIPSO with most regions showing proportional correct values of 0.90 without any spatial filtering and 0.95 when a spatial filter is applied. Values are similar for the NASA MODIS MYD35 mask. A direct comparison of PATMOS-x and MYD35 from 2003 to 2014 also shows agreement over most regions in terms of mean cloud amount, inter-annual variability, and linear trends. Regional and seasonal differences are discussed. The analysis demonstrates that PATMOS-x cloud amount uncertainty could effectively screen regions where PATMOS-x differs from MYD35.

Description: Land cover mapping (LCM) in complex surface-mined and agricultural landscapes could contribute greatly to regulating mine exploitation and protecting mine geo-environments. However, there are some special and spectrally similar land covers in these landscapes which increase the difficulty in LCM when employing high spatial resolution images. There is currently no research on these mixed complex landscapes. The present study focused on LCM in such a mixed complex landscape located in Wuhan City, China. A procedure combining ZiYuan-3 (ZY-3) stereo satellite imagery, the feature selection (FS) method, and machine learning algorithms (MLAs) (random forest, RF; support vector machine, SVM; artificial neural network, ANN) was proposed and first examined for both LCM of surface-mined and agricultural landscapes (MSMAL) and classification of surface-mined land (CSML), respectively. The mean and standard deviation filters of spectral bands and topographic features derived from ZY-3 stereo images were newly introduced. Comparisons of three MLAs, including their sensitivities to FS and whether FS resulted in significant influences, were conducted for the first time in the present study. The following conclusions are drawn. Textures were of little use, and the novel features contributed to improve classification accuracy. Regarding the influence of FS: FS substantially reduced feature set (by 68% for MSMAL and 87% for CSML), and often improved classification accuracies (with an average value of 4.48% for MSMAL using three MLAs, and 11.39% for CSML using RF and SVM); FS showed statistically significant improvements except for ANN-based MSMAL; SVM was most sensitive to FS, followed by ANN and RF. Regarding comparisons of MLAs: for MSMAL based on feature subset, RF achieved the greatest overall accuracy of 77.57%, followed by SVM and ANN; for CSML, SVM had the highest accuracies (87.34%), followed by RF and ANN; based on the feature subsets, significant differences were observed for MSMAL and CSML using any pair of MLAs. In general, the proposed approach can contribute to LCM in complex surface-mined and agricultural landscapes.

Description: A sequence of Ms ≥ 5.0 earthquakes occurred in 2003 and 2004 in Bange County, Tibet, China, all with similar depths and focal mechanisms. However, the source parameters, kinematics and relationships between these earthquakes are poorly known because of their moderately-sized magnitude and the sparse distribution of seismic stations in the region. We utilize interferometric synthetic aperture radar (InSAR) data from the European Space Agency’s Envisat satellite to determine the location, fault geometry and slip distribution of three large events of the sequence that occurred on 7 July 2003 (Ms 6.0), 27 March 2004 (Ms 6.2), and 3 July 2004 (Ms 5.1). The modeling results indicate that the 7 July 2003 event was a normal-faulting event with a right-lateral slip component, the 27 March 2004 earthquake was associated with a normal fault striking northeast–southwest and dipping northwest with a moderately oblique right-lateral slip, and the 3 July 2004 event was caused by a normal fault. A calculation of the static stress changes on the fault planes demonstrates that the third earthquake may have been triggered by the previous ones.

Description: Information on ice is important for shipping, weather forecasting, and climate monitoring. Historically, ice cover has been detected and ice concentration has been measured using relatively low-resolution space-based passive microwave data. This study presents an algorithm to detect ice and estimate ice concentration in clear-sky areas over the ocean and inland lakes and rivers using high-resolution data from the Visible Infrared Imaging Radiometer Suite (VIIRS) onboard the Suomi National Polar Orbiting Partnership (S-NPP) and on future Joint Polar Satellite System (JPSS) satellites, providing spatial detail that cannot be obtained with passive microwave data. A threshold method is employed with visible and infrared observations to identify ice, then a tie-point algorithm is used to determine the representative reflectance/temperature of pure ice, estimate the ice concentration, and refine the ice cover mask. The VIIRS ice concentration is validated using observations from Landsat 8. Results show that VIIRS has an overall bias of −0.3% compared to Landsat 8 ice concentration, with a precision (uncertainty) of 9.5%. Biases and precision values for different ice concentration subranges from 0% to 100% can be larger.

Description: Optical and microwave images have been combined for land cover monitoring in different agriculture scenarios, providing useful information on qualitative and quantitative land cover changes. This study aims to assess the complementarity and interoperability of optical (SPOT-5 Take-5) and synthetic aperture radar (SAR) (Sentinel-1A) data for crop parameter (basal crop coefficient (Kcb) values and the length of the crop’s development stages) retrieval and crop type classification, with a focus on crop water requirements, for an irrigation perimeter in Angola. SPOT-5 Take-5 images are used as a proxy of Sentinel-2 data to evaluate the potential of their enhanced temporal resolution for agricultural applications. In situ data are also used to complement the Earth O=observation (EO) data. The Normalized Difference Vegetation Index (NDVI) and dual (VV + VH) polarization backscattering time series are used to compute the Kcb curve for four crop types (maize, soybean, bean and pasture) and to estimate the length of each phenological growth stage. The Kcb values are then used to compute the crop’s evapotranspiration and to subsequently estimate the crop irrigation requirements based on a soil water balance model. A significant R2 correlation between NDVI and backscatter time series was observed for all crops, demonstrating that optical data can be replaced by microwave data in the presence of cloud cover. However, it was not possible to properly identify each stage of the crop cycle due to the lack of EO data for the complete growing season.

Description: Determining the relationship between crustal movement and faulting in thrust belts is essential for understanding the growth of geological structures and addressing the proposed models of a potential earthquake hazard. A Mw 5.9 earthquake occurred on 21 January 2016 in Menyuan, NE Qinghai Tibetan plateau. We combined satellite interferometry from Sentinel-1A Terrain Observation with Progressive Scans (TOPS) images, historical earthquake records, aftershock relocations and geological data to determine fault seismogenic structural geometry and its relationship with the Lenglongling faults. The results indicate that the reverse slip of the 2016 earthquake is distributed on a southwest dipping shovel-shaped fault segment. The main shock rupture was initiated at the deeper part of the fault plane. The focal mechanism of the 2016 earthquake is quite different from that of a previous Ms 6.5 earthquake which occurred in 1986. Both earthquakes occurred at the two ends of a secondary fault. Joint analysis of the 1986 and 2016 earthquakes and aftershocks distribution of the 2016 event reveals an intense connection with the tectonic deformation of the Lenglongling faults. Both earthquakes resulted from the left-lateral strike-slip of the Lenglongling fault zone and showed distinct focal mechanism characteristics. Under the shearing influence, the normal component is formed at the releasing bend of the western end of the secondary fault for the left-order alignment of the fault zone, while the thrust component is formed at the restraining bend of the east end for the right-order alignment of the fault zone. Seismic activity of this region suggests that the left-lateral strike-slip of the Lenglongling fault zone plays a significant role in adjustment of the tectonic deformation in the NE Tibetan plateau.

Description: Nitrogen is an essential nutrient element in crop photosynthesis and yield improvement. Thus, it is urgent and important to accurately estimate the leaf nitrogen contents (LNC) of crops for precision nitrogen management. Based on the correlation between LNC and reflectance spectra, the hyperspectral LiDAR (HSL) system can determine three-dimensional structural parameters and biochemical changes of crops. Thereby, HSL technology has been widely used to monitor the LNC of crops at leaf and canopy levels. In addition, the laser-induced fluorescence (LIF) of chlorophyll, related to the histological structure and physiological conditions of green plants, can also be utilized to detect nutrient stress in crops. In this study, four regression algorithms, support vector machines (SVMs), partial least squares (PLS) and two artificial neural networks (ANNs), back propagation NNs (BP-NNs) and radial basic function NNs (RBF-NNs), were selected to estimate rice LNC in booting and heading stages based on reflectance and LIF spectra. These four regression algorithms were used for 36 input variables, including the reflectance spectral variables on 32 wavelengths and four peaks of the LIF spectra. A feature weight algorithm was proposed to select different band combinations for the LNC retrieval models. The determination coefficient (R2) and the root mean square error (RMSE) of the retrieval models were utilized to compare their abilities of estimating the rice LNC. The experimental results demonstrate that (I) these four regression methods are useful for estimating rice LNC in the order of RBF-NNs > SVMs > BP-NNs > PLS; (II) The LIF data in two forms, including peaks and indices, display potential in rice LNC retrieval, especially when using the PLS regression (PLSR) model for the relationship of rice LNC with spectral variables. The feature weighting algorithm is an effective and necessary method to determine appropriate band combinations for rice LNC estimation.

Description: Soil moisture is an important variable in the coupled hydrologic and climate system. In recent years, microwave-based soil moisture products have been shown to be a viable alternative to in situ measurements. A popular way to measure the performance of soil moisture products is to calculate the temporal correlation coefficient (R) against in situ measurements or other appropriate reference datasets. In this study, an existing linear combination method improving R was modified to allow for a non-static or nonstationary model combination as the basis for improving remotely-sensed surface soil moisture. Previous research had noted that two soil moisture products retrieved using the Japan Aerospace Exploration Agency (JAXA) and Land Parameter Retrieval Model (LPRM) algorithms from the same Advanced Microwave Scanning Radiometer 2 (AMSR2) sensor are spatially complementary in terms of R against a suitable reference over a fixed period. Accordingly, a linear combination was proposed to maximize R using a set of spatially-varying, but temporally-fixed weights. Even though this approach showed promising results, there was room for further improvements, in particular using non-static or dynamic weights that take account of the time-varying nature of the combination algorithm being approximated. The dynamic weighting was achieved by using a moving window. A number of different window sizes was investigated. The optimal weighting factors were determined for the data lying within the moving window and then used to dynamically combine the two parent products. We show improved performance for the dynamically-combined product over the static linear combination. Generally, shorter time windows outperform the static approach, and a 60-day time window is suggested to be the optimum. Results were validated against in situ measurements collected from 124 stations over different continents. The mean R of the dynamically-combined products was found to be 0.57 and 0.62 for the cases using the European Centre for Medium-Range Weather Forecasts Reanalysis-Interim (ERA-Interim) and Modern-Era Retrospective Analysis for Research and Applications Land (MERRA-Land) reanalysis products as the reference, respectively, outperforming the statically-combined products (0.55 and 0.54).

Description: Moderate spatial resolution satellite data from the Landsat-8 OLI and Sentinel-2A MSI sensors together offer 10 m to 30 m multi-spectral reflective wavelength global coverage, providing the opportunity for improved combined sensor mapping and monitoring of the Earth’s surface. However, the standard geolocated Landsat-8 OLI L1T and Sentinel-2A MSI L1C data products are currently found to be misaligned. An approach for automated registration of Landsat-8 OLI L1T and Sentinel-2A MSI L1C data is presented and demonstrated using contemporaneous sensor data. The approach is computationally efficient because it implements feature point detection across four image pyramid levels to identify a sparse set of tie-points. Area-based least squares matching around the feature points with mismatch detection across the image pyramid levels is undertaken to provide reliable tie-points. The approach was assessed by examination of extracted tie-point spatial distributions and tie-point mapping transformations (translation, affine and second order polynomial), dense-matching prediction-error assessment, and by visual registration assessment. Two test sites over Cape Town and Limpopo province in South Africa that contained cloud and shadows were selected. A Landsat-8 L1T image and two Sentinel-2A L1C images sensed 16 and 26 days later were registered (Cape Town) to examine the robustness of the algorithm to surface, atmosphere and cloud changes, in addition to the registration of a Landsat-8 L1T and Sentinel-2A L1C image pair sensed 4 days apart (Limpopo province). The automatically extracted tie-points revealed sensor misregistration greater than one 30 m Landsat-8 pixel dimension for the two Cape Town image pairs, and greater than one 10 m Sentinel-2A pixel dimension for the Limpopo image pair. Transformation fitting assessments showed that the misregistration can be effectively characterized by an affine transformation. Hundreds of automatically located tie-points were extracted and had affine-transformation root-mean-square error fits of approximately 0.3 pixels at 10 m resolution and dense-matching prediction errors of similar magnitude. These results and visual assessment of the affine transformed data indicate that the methodology provides sub-pixel registration performance required for meaningful Landsat-8 OLI and Sentinel-2A MSI data comparison and combined data applications.

Description: A Mw 7.1 earthquake struck the Eastern Anatolia, near the city of Van (Turkey), on 23 October 2011. We investigated the coseismic surface displacements using the InSAR technique, exploiting adjacent ENVISAT tracks and COSMO-SkyMed images. Multi aperture interferometry was also applied, measuring ground displacements in the azimuth direction. We solved for the fault geometry and mechanism, and we inverted the slip distribution employing a numerical forward model that includes the available regional structural data. Results show a horizontally elongated high slip area (7–9 m) at 12–17 km depth, while the upper part of the fault results unruptured, enhancing its seismogenic potential. We also investigated the post-seismic phase acquiring most of the available COSMO-SkyMed, ENVISAT and TERRASAR-X SAR images. The computed afterslip distributions show that the shallow section of the fault underwent considerable aseismic slip during the early days after the mainshock, of tens of centimeters. Our results support the hypothesis of a seismogenic potential reduction within the first 8–10 km of the fault through the energy release during the post-seismic phase. Despite non-optimal data coverage and coherence issues, we demonstrate that useful information about the Van earthquake could still be retrieved from SAR data through detailed analysis.

Description: The application of custom classification techniques and posterior probability modeling (PPM) using Worldview-2 multispectral imagery to archaeological field survey is presented in this paper. Research is focused on the identification of Neolithic felsite stone tool workshops in the North Mavine region of the Shetland Islands in Northern Scotland. Sample data from known workshops surveyed using differential GPS are used alongside known non-sites to train a linear discriminant analysis (LDA) classifier based on a combination of datasets including Worldview-2 bands, band difference ratios (BDR) and topographical derivatives. Principal components analysis is further used to test and reduce dimensionality caused by redundant datasets. Probability models were generated by LDA using principal components and tested with sites identified through geological field survey. Testing shows the prospective ability of this technique and significance between 0.05 and 0.01, and gain statistics between 0.90 and 0.94, higher than those obtained using maximum likelihood and random forest classifiers. Results suggest that this approach is best suited to relatively homogenous site types, and performs better with correlated data sources. Finally, by combining posterior probability models and least-cost analysis, a survey least-cost efficacy model is generated showing the utility of such approaches to archaeological field survey.

Description: Multitemporal biodiversity data on a forest ecosystem can provide useful information about the evolution of biodiversity in a territory. The present study describes the recovery of an archive used to determine the main Schmid’s vegetation belts in Trento Province, Italy. The archive covers 20 years, from the 1970s to the 1990s. During the FORCING project (an Italian acronym for Cingoli Forestali, i.e., forest belts), a comprehensive process of database recovering was executed, and missing data were digitized from historical maps, preserving paper-based maps and documents. All of the maps of 16 forest districts, and the related 8000 detected transects, have been georeferenced to make the whole database spatially explicit and to evaluate the possibility of performing comparative samplings on up-to-date datasets. The floristic raw data (approximately 200,000 specific identifications, including frequency indices) still retain an important and irreplaceable information value. The data can now be browsed via a web-GIS. We provide here a set of examples of the use of this type of data, and we highlight the potential and the limits of the specific dataset and of the historical database, in general.

Description: A Bayesian random effects modeling approach was used to examine the influence of neighborhood characteristics on burglary risks in Jianghan District, Wuhan, China. This random effects model is essentially spatial; a spatially structured random effects term and an unstructured random effects term are added to the traditional non-spatial Poisson regression model. Based on social disorganization and routine activity theories, five covariates extracted from the available data at the neighborhood level were used in the modeling. Three regression models were fitted and compared by the deviance information criterion to identify which model best fit our data. A comparison of the results from the three models indicates that the Bayesian random effects model is superior to the non-spatial models in fitting the data and estimating regression coefficients. Our results also show that neighborhoods with above average bar density and department store density have higher burglary risks. Neighborhood-specific burglary risks and posterior probabilities of neighborhoods having a burglary risk greater than 1.0 were mapped, indicating the neighborhoods that should warrant more attention and be prioritized for crime intervention and reduction. Implications and limitations of the study are discussed in our concluding section.

Description: Decision-makers in the fields of urban and regional planning in Germany face new challenges. High rates of urban sprawl need to be reduced by increased inner-urban development while settlements have to adapt to climate change and contribute to the reduction of greenhouse gas emissions at the same time. In this study, we analyze conflicts in the management of urban areas and develop integrated sustainable land use strategies for Germany. The spatial explicit land use change model Land Use Scanner is used to simulate alternative scenarios of land use change for Germany for 2030. A multi-criteria analysis is set up based on these scenarios and based on a set of indicators. They are used to measure whether the mitigation and adaptation objectives can be achieved and to uncover conflicts between these aims. The results show that the built-up and transport area development can be influenced both in terms of magnitude and spatial distribution to contribute to climate change mitigation and adaptation. Strengthening the inner-urban development is particularly effective in terms of reducing built-up and transport area development. It is possible to reduce built-up and transport area development to approximately 30 ha per day in 2030, which matches the sustainability objective of the German Federal Government for the year 2020. In the case of adaptation to climate change, the inclusion of extreme flood events in the context of spatial planning requirements may contribute to a reduction of the damage potential.

Description: Gebhardt et al. (2014) presented the Monitoring Activity Data for the Mexican REDD+ program (MAD-MEX), an automatic nation-wide land cover monitoring system for the Mexican REDD+ MRV. Though MAD-MEX represents a valuable first effort toward establishing a national reference emissions level for the implementation of REDD+ in Mexico, in this paper, we argue that this land cover system has important limitations that may prevent it from becoming operational for REDD+ MRV. Specifically, we show that (1) the accuracy assessment of MAD-MEX land cover maps is optimistically biased; (2) the ability of MAD-MEX to monitor land cover change, including deforestation and forest degradation; is poor and (3) the use of an entirely automatic classification approach, such as that followed by MAD-MEX, is highly problematic in the case of a large and heterogeneous country like Mexico. We discuss these limitations and call into question the ability of a land cover monitoring system, such as MAD-MEX, both to elaborate a national reference emissions level and to monitor future forest cover change, as part of a REDD+ MRV system. We provide some insights with the aim of improving the development of nation-wide land cover monitoring systems in Mexico and elsewhere.

Description: Mas, J.F. et al. have submitted a paper [1] for publication, which aims to respond to a paper published by Gebhardt et al. [2]. Mas, J.F. et al. had received a consultancy in 2013 to assess the quality of the early prototype products partly described in Gebhardt et al. in 2014. This consultancy, although a formal non-disclosure agreement had not been demanded, was awarded under the mutual understanding that the data handed over to Mas et al. constitute the early development phase of the program. Therefore, Mas et al. had been asked to give an assessment on the quality of the prototypes to obtain a proof of concept for the proposed workflow of MAD-Mex. It was clear that this assessment would suffer from limited availability of high quality training and validation data available in 2013. Mas et al. finally did not execute the consultancy due to the limited vector processing capacities in their lab. In October 2014, we sent the latest products, version 4.2 of the MAD-Mex products, including the more than 200,000 validation points gathered from independent expert interpreters of all Mexican ecosystems. Mas et al. did not respond to this transfer or to our request to collaborate in the quality control and assessment of MAD-Mex.

Description: The combined usage of high-resolution satellite images and thermal infrared (TIR) data helps understanding the thermal effect of urban fabric properties and the mechanism of urban heat island (UHI) formation. In this study, three typical urban functional zones (UFZs) of downtown Shanghai were chosen for quantifying the relationship between fine-scale urban fabric properties and their thermal effect. Nine land surfaces and 146 aggregated land parcels extracted from a QuickBird image (dated 14 April 2014) were used to characterize urban fabric properties. The thermal effect was deduced from land surface temperature (LST), intra-UHI intensity, blackbody flux density (BBFD) and blackbody flux (BBF). The net BBF was retrieved from the Landsat 8 TIR band 10 dated 13 August 2013, and 28 May 2014. The products were resampled to fine resolution using a geospatial sharpening approach and further validated. The results show that: (1) On the UFZ level, there is a significant thermal differential among land surfaces. Water, well-vegetated land, high-rises with light color and high-rises with glass curtain walls exhibited relatively low LST, UHI intensity and BBFD. In contrast, mobile homes with light steel roofs, low buildings with bituminous roofs, asphalt roads and composite material pavements showed inverse trends for LST, UHI intensity, and BBFD; (2) It was found that parcel-based per ha net BBF, which offsets the “size-effect” among parcels, is more reasonable and comparable when quantifying excess surface flux emitted by the parcels; (3) When examining the relationship between parcel-level land surfaces and per ha BBF, a partial least squares (PLS) regression model showed that buildings and asphalt roads are major contributors to parcel-based per ha BBF, followed by other impervious surfaces. In contrast, vegetated land and water contribute with a much lower per ha net BBF to parcel warming.

Description: Wetlands worldwide are becoming increasingly degraded, and this has motivated many attempts to manage and restore wetland ecosystems. Restoration actions require a large resource investment, so it is critical to measure the outcomes of these management actions. We evaluated the restoration of floodplain wetland vegetation across a chronosequence of land uses, using remote sensing analyses. We compared the Landsat-based fractional cover of restoration areas with river red gum and lignum reference communities, which functioned as a fixed target for restoration, over three time periods: (i) before agricultural land use (1987–1997); (ii) during the peak of agricultural development (2004–2007); and (iii) post-restoration of flooding (2010–2015). We also developed LiDAR-derived canopy height models (CHMs) for comparison over the second and third time periods. Inundation was crucial for restoration, with many fields showing little sign of similarity to target vegetation until after inundation, even if agricultural land uses had ceased. Fields cleared or cultivated for only one year had greater restoration success compared to areas cultivated for three or more years. Canopy height increased most in the fields that were cleared and cultivated for a short duration, in contrast to those cultivated for >12 years, which showed few signs of recovery. Restoration was most successful in fields with a short development duration after the intervention, but resulting dense monotypic stands of river cooba require future monitoring and possibly intervention to prevent sustained dominance. Fields with intensive land use histories may need to be managed as alternative, drier flood-dependent vegetation communities, such as black box (Eucalyptus largiflorens) grasslands. Remotely-sensed data provided a powerful measurement technique for tracking restoration success over a large floodplain.

Description: The atmospheric correction of satellite observations is crucial for both land and ocean remote sensing. However, the optimal approach for each area is different due to the large spectra difference in the ground reflectance between land and ocean. A unified atmospheric correction (UAC) approach based on a look-up table (LUT) of in situ measurements is developed to remove this difference. The LUT is used to select one spectrum as the in situ ground reflectance needed to obtain the initial aerosol reflectance, which in turn is used for determining the two closest aerosol models. The aerosol reflectance, obtained from these aerosol models, is then used to deduce the estimated ground reflectance. This UAC model is then used to process the Sea-viewing Wide Field-of-view Sensor (SeaWiFS) data, and its performance is validated with a large number of in situ measurements. The mean bias of the land reflectance for this model is 6.59% with a root mean square error (RMSE) of 19.61%. The mean bias and RMSE of the water-leaving reflectance are 7.59% and 17.10% validated by the in situ measurements using the above-water method, while they are 13.60% and 22.53% using the in-water method. The UAC model provides a useful tool for correcting the satellite-received reflectance without separately having to deal with land and ocean pixels. Further, it can seamlessly expand the satellite ocean color data for terrestrial use and improve quantitative remote sensing over land.

Description: The increasing spread and abundance of an invasive perennial grass, buffelgrass (Pennisetum ciliare), represents a critical threat to the native vegetation communities of the Sonoran desert in southern Arizona, USA, where buffelgrass eradication is a high priority for resource managers. Herbicidal treatment of buffelgrass is most effective when the vegetation is actively growing, but the remoteness of infestations and the erratic timing and length of the species’ growth periods confound effective treatment. The goal of our research is to promote buffelgrass management by using remote sensing data to detect where the invasive plants are located and when they are photosynthetically active. We integrated citizen scientist observations of buffelgrass phenology in the Tucson, Arizona area with PRISM precipitation data, eight-day composites of 250-m Moderate-resolution Imaging Spectroradiometer (MODIS) satellite imagery, and aerially-mapped polygons of buffelgrass presence to understand dynamics and relationships between precipitation and the timing and amount of buffelgrass greenness from 2011 to 2013. Our results show that buffelgrass responds quickly to antecedent rainfall: in pixels containing buffelgrass, higher correlations (R2 > 0.5) typically occur after two cumulative eight-day periods of rain, whereas in pixels dominated by native vegetation, four prior 8-day periods are required to reach that threshold. Using the new suite of phenometrics introduced here—Climate Landscape Response metrics—we accurately predicted the location of 49% to 55% of buffelgrass patches in Saguaro National Park. These metrics and the suggested guidelines for their use can be employed by resource managers to treat buffelgrass during optimal time periods.

Description: Recent satellite observations yield estimates of the distribution of sea ice thickness across the entire Arctic Ocean. While these sensors were only placed in operation within the last few years, information from other sensors may assist us with estimating the distribution of sea ice thickness in the Arctic beginning in the 1980s. A previous study found that the age of sea ice is correlated to sea ice thickness from 2003 to 2006, but an extension of the temporal analysis is needed to better quantify this relationship and its variability from year to year. Estimates of the ice age/thickness relationship may allow the thickness record to be extended back to 1985, the beginning of our ice age dataset. Comparisons of ice age and thickness estimates derived from both ICESat (2004–2008) and IceBridge (2009–2015) reveal that the relationship between age and thickness differs between these two campaigns, due in part to the difference in area of coverage. Nonetheless, sea ice thickness and age exhibit a direct relationship when compared on pan-Arctic or regional spatial scales.

Description: The body of scientific literature on slum mapping employing remote sensing methods has increased since the availability of more very-high-resolution (VHR) sensors. This improves the ability to produce information for pro-poor policy development and to build methods capable of supporting systematic global slum monitoring required for international policy development such as the Sustainable Development Goals. This review provides an overview of slum mapping-related remote sensing publications over the period of 2000–2015 regarding four dimensions: contextual factors, physical slum characteristics, data and requirements, and slum extraction methods. The review has shown the following results. First, our contextual knowledge on the diversity of slums across the globe is limited, and slum dynamics are not well captured. Second, a more systematic exploration of physical slum characteristics is required for the development of robust image-based proxies. Third, although the latest commercial sensor technologies provide image data of less than 0.5 m spatial resolution, thereby improving object recognition in slums, the complex and diverse morphology of slums makes extraction through standard methods difficult. Fourth, successful approaches show diversity in terms of extracted information levels (area or object based), implemented indicator sets (single or large sets) and methods employed (e.g., object-based image analysis (OBIA) or machine learning). In the context of a global slum inventory, texture-based methods show good robustness across cities and imagery. Machine-learning algorithms have the highest reported accuracies and allow working with large indicator sets in a computationally efficient manner, while the upscaling of pixel-level information requires further research. For local slum mapping, OBIA approaches show good capabilities of extracting both area- and object-based information. Ultimately, establishing a more systematic relationship between higher-level image elements and slum characteristics is essential to train algorithms able to analyze variations in slum morphologies to facilitate global slum monitoring.