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

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

A New Algorithm to Classify the Homogeneity of ERS-2 Wave Mode SAR Imagette

A new classification parameter is developed using 1535 ERS-2 wave mode synthetic aperture radar (SAR) test imagettes to better differentiate homogeneous and inhomogeneous imagettes. The comparison between the new parameter (Min) and the previous one (Inhomo) (Schulz-Stellenfleth and Lehner, 2004) was done under varied threshold values of Inhomo. It is concluded that the performance of ‘Min’ is much better than ‘Inhomo’ when applying to the 1535 test imagettes. Furthermore, both Min and Inhomo are applied to nearly 1 million imagettes collected for the period from 1 September 1998 to 30 November 2000. The comparisons of the global inhomogeneous distribution between ‘Min’ and ‘Inhomo’ reveal that both the areas and percentage of inhomogeneity calculated by ‘Min’ are larger than that calculated by ‘Inhomo’. By analyzing the low wind speed distribution of HOAPS data, we found that low wind speed over the ocean is one of the key reasons for the inhomogeneity of SAR imagettes.     

Global SMOS Soil Moisture Retrievals from The Land Parameter Retrieval Model

A recent study by Van der Schalie et al. (2015) showed good results for applying the Land Parameter Retrieval Model (LPRM) on SMOS observations over southeast Australia and optimizing and evaluating the retrieved soil moisture (θ in m³ m⁻³) against ground measurements from the OzNet sites. In this study, the LPRM parameterization is globally updated for SMOS against modelled θ from MERRA-Land (MERRA) and ERA-Interim/Land (ERA) over the period of July 2010–December 2010, mainly focusing on two parameters: the single scattering albedo (ω) and the roughness (h). The Pearson's coefficient of correlation (r) increased rapidly when increasing the ω up to 0.12 and reached a steady state from thereon, no significant spatial pattern was found in the estimation of the single scattering albedo, which could be an artifact of the used parameter estimation procedure, and a single value of 0.12 was therefore used globally. The h was defined as a function of θ and varied slightly for the different angle bins, with maximum values of 1.1–1.3 as the angle changes from 42.5° to 57.5°.This resulted in an average r of 0.51 and 0.47, with a bias (m³ m⁻³) of −0.02 and −0.01 and an unbiased root mean square error (ubrmse in m³ m⁻³) of 0.054 and 0.056 against MERRA (ascending and descending). For ERA this resulted in an r of 0.61 and 0.53, with a bias of −0.03 and an ubrmse 0.055 and 0.059. The resulting parameterization was then used to run LPRM on SMOS observations over the period of July 2010–December 2013 and evaluated against SMOS Level 3 (L3) θ and available in situ measurements from the International Soil Moisture Network (ISMN). The comparison with L3 shows that the LPRM θ retrievals are very similar, with for the ascending set very high r of over 0.9 in large parts of the globe, with an overall average of 0.85 and the descending set performing less with an average of 0.74, mainly due to the negative r over the Sahara. The mean bias is 0.03, with an ubrmse of 0.038 and 0.044. In this study there are three major areas where the LPRM retrievals do not perform well: very dry sandy areas, densely forested areas and over high latitudes, which are all known limitations of LPRM. The comparison against in situ measurement from the ISMN give very similar results, with average r for LPRM of 0.65 and 0.61 (0.64 and 0.59 for L3) for the ascending and descending sets, while having a comparable bias and ubrmse over the different networks. This shows that LPRM used on SMOS observations produce θ retrievals with a similar quality as the SMOS L3 product.     

Automatic detection of residential buildings using LIDAR data and multispectral imagery

This paper presents an automatic building detection technique using LIDAR data and multispectral imagery. Two masks are obtained from the LIDAR data: a ‘primary building mask’ and a ‘secondary building mask’. The primary building mask indicates the void areas where the laser does not reach below a certain height threshold. The secondary building mask indicates the filled areas, from where the laser reflects, above the same threshold. Line segments are extracted from around the void areas in the primary building mask. Line segments around trees are removed using the normalized difference vegetation index derived from the orthorectified multispectral images. The initial building positions are obtained based on the remaining line segments. The complete buildings are detected from their initial positions using the two masks and multispectral images in the YIQ colour system. It is experimentally shown that the proposed technique can successfully detect urban residential buildings, when assessed in terms of 15 indices including completeness, correctness and quality.     

Assessing optical earth observation systems for mapping and monitoring temporary ponds in arid areas

Remote sensing methods for locating and monitoring temporary ponds over large areas in arid lands were tested on a study site in Northern Senegal. Three main results are presented, validated with field data and intended to highlight different spectral, spatial and temporal characteristics of the methods: (1) Among several water indices tested, two Middle Infrared-based indices (MNDWI—Modified Normalized Difference Water Index and NDWI₁—Normalized Difference Water Index) are found to be most efficient; (2) an objective method is given prescribing the necessary sensor spatial resolution in terms of minimal detected pond area; and (3) the potential of multi-temporal MODIS imagery for tracking the filling phases of small ponds is illustrated. These results should assist in epidemiological studies of vector-borne diseases that develop around these ponds, but also more generally for land and water management and preservation of threatened ecosystems in arid areas.     

Assessing reference evapotranspiration at regional scale based on remote sensing,weather forecast and GIS tools

Reference evapotranspiration (ET_o) is a key component in efficient water management, especially in arid and semi-arid environments. However, accurate ET_o assessment at the regional scale is complicated by the limited number of weather stations and the strict requirements in terms of their location and surrounding physical conditions for the collection of valid weather data. In an attempt to overcome this limitation, new approaches based on the use of remote sensing techniques and weather forecast tools have been proposed.Use of the Land Surface Analysis Satellite Application Facility (LSA SAF) tool and Geographic Information Systems (GIS) have allowed the design and development of innovative approaches for ET_o assessment, which are especially useful for areas lacking available weather data from weather stations. Thus, by identifying the best-performing interpolation approaches (such as the Thin Plate Splines, TPS) and by developing new approaches (such as the use of data from the most similar weather station, TS, or spatially distributed correction factors, CITS), errors as low as 1.1% were achieved for ET_o assessment. Spatial and temporal analyses reveal that the generated errors were smaller during spring and summer as well as in homogenous topographic areas.The proposed approaches not only enabled accurate calculations of seasonal and daily ET_o values, but also contributed to the development of a useful methodology for evaluating the optimum number of weather stations to be integrated into a weather station network and the appropriateness of their locations. In addition to ET_o, other variables included in weather forecast datasets (such as temperature or rainfall) could be evaluated using the same innovative methodology proposed in this study.     

A globally applicable, season-specific model for estimating the weighted mean temperature of the atmosphere

In GPS meteorology, the weighted mean temperature is usually obtained by using a linear function of the surface temperature T _s. However, not every GPS station can measure the surface temperature. The current study explores the characteristics of surface temperature and weighted mean temperature based on the global pressure and temperature model (GPT) and the Bevis T _m–T _s relationship (T _m =?a?+?bT _s). A new global weighted mean temperature (GWMT) model has been built which directly uses three-dimensional coordinates and day of the year to calculate the weighted mean temperature. The data of year 2005–2009 from 135 radiosonde stations provided by the Integrated Global Radiosonde Archive were used to calculate the model coefficients, which have been validated through examples. The result shows that the GWMT model is generally better than the existing liner models in most areas according to the statistic indexes (namely, mean absolute error and root mean square). Then we calculated precipitable water vapor, and the result shows that GWMT model can also yield high precision PWV.     

An updated delineation of stand ages of deciduous rubber plantations during 1987-2018 using Landsat-derived bi-temporal thresholds method in an anti-chronological strategy

Timely and accurately monitoring stand ages of deciduous rubber plantations is of great importance for ecological studies and plantations management. The re-establishment of rubber plantations usually experiences a short period (several years) of land clearance and transplantation of rubber seedlings, along with a noticeable landscape change from well-grown forest to bare land and sparse vegetation in situ. With Landsat times series (LTS) data of four commonly-used vegetation indices (VIs), namely the Normalized Difference Vegetation Index (NDVI), Enhanced Vegetation Index (EVI), Normalized Difference Moisture Index (NDMI), and Normalized Burn Ratio (NBR), and three non-visible spectral bands, i.e. the near-infrared (NIR) and shortwave-infrared (SWIR1/2), in this study, an approach by combining the inter-annual defoliating and foliating features of rubber trees and the intra-annual landscape changes of rubber plantations was presented to detect and map stand ages of deciduous rubber plantations in an anti-chronological manner across Xishuangbanna between 1987 and 2018, one of the most intensive regions of deciduous rubber plantations within the tropics. The approach highlighted the repeated distribution of newly-cleared and replanted plot (NCRP) of rubber seedlings due to rubber management. It applied the bi-temporal VIs thresholds of zero of NBR and NDMI during the defoliation to foliation phases to delineate the stand ages of deciduous rubber plantations at an interval of five years, by combining a Landsat-based rubber map in 2018 and 32-year NCRPs as well as quadri-classified age-groups and seven sub-categories (i.e. ≤5 as infantile rubber plantations (IRP), 6–10 as young rubber plantations (YRP), 11–15 and 16–20 as mature rubber plantations (MRP), 21–25, 26–30, and ≥31 years as old rubber plantations (ORP)). The results showed that the areas of IRP, YRP, MRP, and ORP were 19.1 km², 817.1 km², 1681.7 km², and 573.7 km² in 2018, respectively. Spatially, the YRP are mainly around the outskirts of two county-level administrative centers (Jinghong and Mengla), while ORP primarily distributed along main roads. Nearly 53.9% of ORP, 51.8% of IRP, 47.3% of MRP and 46.3% of YRP were in Jinghong City, and Mengla County had 50.5% of YRP, 48.8% of MRP, 42.4% of IRP and 36.3% of ORP. This study demonstrates that the bi-temporal VIs thresholds method (i.e. NBR_defoliation <0, NDMI_defoliation <0, NBR_foliation <0, and NDMI_foliation <0) have great potential for detecting stand ages of deciduous rubber plantations.     

Impact of land-cover layout on particulate matter 2.5 in urban areas of China

ABSTRACT

Urbanization in China is closely connected with ambient particulate matter 2.5 (PM_2.5). However, the potential for altering PM_2.5 through the urban landscape characteristics is uncertain. In this study, we analyzed the urban PM_2.5 pollution situation for 2014–2016 and investigated the impact of landscape factors on urban PM_2.5 in China at the city level. All the prefecture-level cities were stratified by urban population size into small (<500,000), medium (500,000–1,000,000), and large (>1,000,000), and the other second-level administrative cities were assigned as ‘other’ cities. The multivariate regression model including both urban landscape factors and social-economic variables explained 70.0%, 32.8%, 19.2%, and 12.4% of the arithmetic mean PM_2.5 concentration (AMC-PM_2.5) for the other, small, medium, and large cities, respectively. With regard to the configuration of land cover, agricultural activity is a major contributor of PM_2.5 pollution, for which the explanatory power ranged from 7.6% (for the large cities) to 64% (for the other cities). In addition, grassland aggregation also has a limited but negative effect on urban PM_2.5 pollution, despite the negligible effect on dry deposition. Overall, these findings likely reflect the interaction between urban air quality and urbanization, and will have implications for air quality control strategies.     

Step-wise Land-class Elimination Approach for extracting mixed-type built-up areas of Kolkata megacity

The extraction of urban built-up areas is an important aspect of urban planning and understanding the complex drivers and biophysical mechanism of urban climate processes. However, built-up area extraction using Landsat data is a challenging task due to spatio-temporal dynamics and spatially intermixed nature of Land Use and Land Cover (LULC) in the cities of the developing countries, particularly in tropics. In the light of advantages and drawbacks of the Normalized Difference Built-up Index (NDBI) and Built-up Area Extraction Method (BAEM), a new and simple method i.e. Step-wise Land-class Elimination Approach (SLEA) is proposed for rapid and accurate mapping of urban built-up areas without depending exclusively on the band specific normalized indices, in order to pursue a more generalized approach. It combines the use of a single band layer, Normalized Difference Vegetation Index (NDVI) image and another binary image obtained through Logit model. Based on the spectral designation of the satellite image in use, a particular band is chosen for identification of water pixels. The Double-window Flexible Pace Search (DFPS) approach is employed for finding the optimum threshold value that segments the selected band image into water and non-water categories. The water pixels are then eliminated from the original image. The vegetation pixels are similarly identified using the NDVI image and eliminated. The residual pixels left after elimination of water and vegetation categories belong either to the built-up areas or to bare land categories. Logit model is used for separation of the built-up areas from bare lands. The effectiveness of this method was tested through the mapping of built-up areas of the Kolkata Metropolitan Area (KMA), India from Thematic Mapper (TM) images of 2000, 2005 and 2010, and Operational Land Imager (OLI) image of 2015. Results of the proposed SLEA were 95.33% accurate on the whole, while those derived by the NDBI and BAEM approaches returned an overall accuracy of 83.67% and 89.33%, respectively. Comparisons of the results obtained using this method with those obtained from NDBI and BAEM approaches demonstrate that the proposed approach is quite reliable. The SLEA generates new patterns of evidence and hypotheses for built-up areas extraction research, providing an integral link with statistical science and encouraging trans-disciplinary collaborations to build robust knowledge and problem solving capacity in urban areas. It also brings landscape architecture, urban and regional planning, landscape and ecological engineering, and other practice-oriented fields to bear together in processes for identifying problems and analyzing, synthesizng, and evaluating desirable alternatives for urban change. This method produced very accurate results in a more efficient manner compared to the earlier built-up area extraction approaches for the landscape and urban planning.     

On safe ground? Analysis of European urban geohazards using satellite radar interferometry

Urban geological hazards involving ground instability can be costly, dangerous, and affect many people, yet there is little information about the extent or distribution of geohazards within Europe’s urban areas. A reason for this is the impracticality of measuring ground instability associated with the many geohazard processes that are often hidden beneath buildings and are imperceptible to conventional geological survey detection techniques. Satellite radar interferometry, or InSAR, offers a remote sensing technique to map mm-scale ground deformation over wide areas given an archive of suitable multi-temporal data. The EC FP7 Space project named PanGeo (2011–2014), used InSAR to map areas of unstable ground in 52 of Europe’s cities, representing ∼15% of the EU population. In partnership with Europe’s national geological surveys, the PanGeo project developed a standardised geohazard-mapping methodology and recorded 1286 instances of 19 types of geohazard covering 18,000 km². Presented here is an analysis of the results of the PanGeo-project output data, which provides insights into the distribution of European urban geohazards, their frequency and probability of occurrence. Merging PanGeo data with Eurostat’s GeoStat data provides a systematic estimate of population exposures. Satellite radar interferometry is shown to be as a valuable tool for the systematic detection and mapping of urban geohazard phenomena.     

Determining iron content in Mediterranean soils in partly vegetated areas,using spectral reflectance and imaging spectroscopy

The possibility of quantifying iron content in the topsoil of the slopes of the El Hacho Mountain complex in Southern Spain using imaging spectroscopy is investigated. Laboratory, field and airborne spectrometer (ROSIS) data are acquired, in combination with soil samples, which are analysed for dithionite extractable iron (Fe_d) content. Analysis of the properties of two iron related absorption features present in laboratory spectra demonstrates good relations, especially between the standard deviation (S.D.) of the values in an absorption feature and the Fe_d content (R² = 0.67) as well as the ratio based Redness Index (R² = 0.51). Such derived relations are less strong for the ROSIS data (R² for S.D. = 0.26 and R² for Redness Index = 0.22). The spatial distribution of iron in vegetated areas shows a strong sensitivity of these relations with the presence of vegetation. A combination of both methods shows that the overestimation of the Fe_d content with the one method is (partly) compensated by the underestimation with the other method.     

Normalised difference spectral indices and urban land cover as indicators of land surface temperature (LST)

Land cover changes associated with urbanisation modify microclimate, leading to urban heat islands, whereby cities are warmer than the surrounding countryside. Understanding the factors causing this phenomenon could help urban areas adapt to climate change and improve living conditions of inhabitants. In this study, land surface temperatures (LST) of Aarhus, a city in the high latitudes, are estimated from the reflectance of a thermal band (TIRS1; Band 10; 10.60–11.19 μm) of Landsat 8 on five dates in the summer months (one in 2015, and four in 2018). Spectral indices, modelled on the normalised difference vegetation index (NDVI), using all combinations of the first seven bands of Landsat 8 are calculated and their relationships with LST, analysed. Land cover characteristics, in terms of the percentages of tree cover, building cover and overall vegetation cover are estimated from airborne LiDAR data, building footprints and 4-band aerial imagery, respectively. The correlations between LST, the spectral indices and land cover are estimated.The difference in mean temperature between the rural and urban parts of Aarhus is up to 3.96 °C, while the difference between the warmer and colder zones (based on the mean and SD of LST) is up to 13.26 °C. The spectral index using the near infrared band (NIR; Band 5; 0.85-0.88 μm) and a short-wave infrared band (SWIR2; Band 7; 2.11–2.29 μm) has the strongest correlations (r: 0.62 to 0.89) with LST for the whole study area. This index is the inverse of normalised burn ratio (NBR), which has been used for mapping burnt areas. Spectral indices using different combinations of the infrared bands have stronger correlations with LST than the more widely used vegetation indices such as NDVI. The percentage of tree cover has a higher negative correlation (Pearson’s r: -0.68 to -0.75) with LST than overall vegetation cover (r: -0.45 to -0.63). Tree cover and building cover (r: 0.53 to 0.71) together explain up to 68 % of the variation in LST. Modification of tree and building cover may therefore have the potential to regulate urban LST.     

The comparison index: A tool for assessing the accuracy of image segmentation

Segmentation algorithms applied to remote sensing data provide valuable information about the size, distribution and context of landscape objects at a range of scales. However, there is a need for well-defined and robust validation tools to assessing the reliability of segmentation results. Such tools are required to assess whether image segments are based on ‘real’ objects, such as field boundaries, or on artefacts of the image segmentation algorithm. These tools can be used to improve the reliability of any land-use/land-cover classifications or landscape analyses that is based on the image segments.The validation algorithm developed in this paper aims to: (a) localize and quantify segmentation inaccuracies; and (b) allow the assessment of segmentation results on the whole. The first aim is achieved using object metrics that enable the quantification of topological and geometric object differences. The second aim is achieved by combining these object metrics into a ‘Comparison Index’, which allows a relative comparison of different segmentation results. The approach demonstrates how the Comparison Index CI can be used to guide trial-and-error techniques, enabling the identification of a segmentation scale H that is close to optimal. Once this scale has been identified a more detailed examination of the CI–H- diagrams can be used to identify precisely what H value and associated parameter settings will yield the most accurate image segmentation results.The procedure is applied to segmented Landsat scenes in an agricultural area in Saxony-Anhalt, Germany. The segmentations were generated using the ‘Fractal Net Evolution Approach’, which is implemented in the eCognition software.     

星载微波SSM/Ⅰ遥感在中国东北华北农田的辐射特征分析

国家遥感应用工程技术研究中心ＮａｔｉｏｎａｌＥｎｇｉｎｅｅｒｉｎｇＲｅｓｅａｒｃｈＣｅｎｔｅｒｆｏｒＧｅｏｍａｔｉｃｓ（ＮＣＧ）国家遥感应用工程技术研究中心于１９９６年１２月２５日由国家科委正式批准组建（国科发计字［１９９６］６０３号文件）。中国科学...     

Imaging and visualizing the spectral signatures from Landsat TM and ‘τ’ value-based surface soil microzonation mapping at and around Agartala (India)

This work aims to assess the soil microzonation of Agartala city and its surrounding areas based on spectral geophysical signatures. Different spectral resolutions of Landsat TM have been used for assessing the Normalized Difference Vegetative Index, spatial thermal emission representation and plant water moisture representation. Normalized Difference Vegetative Index (NDVI) was measured from band 4 (near-infrared (NIR)) and band 3 (photosynthetically active radiation (PAR)). The Digital Number (DN) values of thermal infrared band (TIR) were used for measuring spatial variation of thermal representation in the city area. A very simple model was developed for measuring thermal emission representative index from NDVI and automated classified TIR band. Overlaid NDVI and classified TIR shows the spatial distribution of thermal emission representative values. Classified mid-wave infrared band (MWIR) was used for measuring the surface geotherm units (τ) which are related with different types of soil. On the basis of spatial distribution of τ value which is clearly visible in a thermal emission representative map overlaid by classified MWIR, the soil microzonation map of the study area was prepared. This soil microzonation map shows that Agartala and its surrounding areas are characterized by four types of soil which are related to different geomorphic and geological units. The soil of this area is classified as dry sandy soil and sandy clay soil of the highland areas and wet sandy alluvium and clayey alluvium of the flood plain area.     

A New 3‐D Solar Radiation Model for 3‐D City Models

Estimates of solar radiation distribution in urban areas are often limited by the complexity of urban environments. These limitations arise from spatial structures such as buildings and trees that affect spatial and temporal distributions of solar fluxes over urban surfaces. The traditional solar radiation models implemented in GIS can address this problem only partially. They can be adequately used only for 2‐D surfaces such as terrain and rooftops. However, vertical surfaces, such as facades, require a 3‐D approach. This study presents a new 3‐D solar radiation model for urban areas represented by 3‐D city models. The v.sun module implemented in GRASS GIS is based on the existing solar radiation methodology used in the topographic r.sun model with a new capability to process 3‐D vector data representing complex urban environments. The calculation procedure is based on the combined vector‐voxel approach segmenting the 3‐D vector objects to smaller polygon elements according to a voxel data structure of the volume region. The shadowing effects of surrounding objects are considered using a unique shadowing algorithm. The proposed model has been applied to the sample urban area with results showing strong spatial and temporal variations of solar radiation flows over complex urban surfaces.     

The urban greenness score: A satellite-based metric for multi-decadal characterization of urban land dynamics

Canada’s urban areas have experienced extensive growth over the past quarter century; however, there has been no consistent, spatially explicit approach for quantifying the loss and gain of greenness in cities nationally. Herein, we developed a novel urban greenness score metric using greenness fractions from a multi-decadal time series (1984–2016) of spectrally unmixed annual Landsat satellite image composites to characterize final year (2016) greenness and its overall change for 18 major Canadian urban areas, summarized by census dissemination area (DA). The applied validation procedure confirmed correlation coefficients (ρ) ranging from 0.67 – 0.85 between reference and estimated greenness fractions, indicating that spectral unmixing is an appropriate method for extracting urban greenness from a time series of medium spatial resolution satellite imagery. Most DAs across Canada sustained a moderate (∼20 % – 40 %) or low (≲ 20 %) level of greenness between 1984 and 2016, but overall there was a decreasing trend in greenness. Eastern urban areas maintained the most greenness over time, while urban areas in the Prairies had the greatest increase in greenness. Densely populated urban areas experienced the greatest loss in greenness (16 % of DAs); whereas, urban areas with a moderately-low density experienced the greatest increase (14 % of DAs). In agreement with previous studies, we found that greenness was negatively associated with urban infilling, with lower greenness levels typically found in urban cores, and greenness loss most often found in the urban periphery in conjunction with urban expansion. Methods presented in this analysis take advantage of the open and longstanding Landsat archive, as well as multiple spatial scales, including sub-pixel unmixing techniques, pixel level greenness faction data summarized for management units, and analysis conducted nationally. The developed urban greenness score provides a comprehensive framework to understand current urban greenness and relate it to its recent past, which supports long-term strategic planning, and can be transferred to other regions across spatial and temporal scales.     

Geospatial scenario based modelling of urban and agricultural intrusions in Ramsar wetland Deepor Beel in Northeast India using a multi-layer perceptron neural network

In recent decades, the world has experienced unprecedented urban growth which endangers the green environment in and around urban areas. In this work, an artificial neural network (ANN) based model is developed to predict future impacts of urban and agricultural expansion on the uplands of Deepor Beel, a Ramsar wetland in the city area of Guwahati, Assam, India, by 2025 and 2035 respectively. Simulations were carried out for three different transition rates as determined from the changes during 2001–2011, namely simple extrapolation, Markov Chain (MC), and system dynamic (SD) modelling, using projected population growth, which were further investigated based on three different zoning policies. The first zoning policy employed no restriction while the second conversion restriction zoning policy restricted urban-agricultural expansion in the Guwahati Municipal Development Authority (GMDA) proposed green belt, extending to a third zoning policy providing wetland restoration in the proposed green belt. The prediction maps were found to be greatly influenced by the transition rates and the allowed transitions from one class to another within each sub-model. The model outputs were compared with GMDA land demand as proposed for 2025 whereby the land demand as produced by MC was found to best match the projected demand. Regarding the conservation of Deepor Beel, the Landscape Development Intensity (LDI) Index revealed that wetland restoration zoning policies may reduce the impact of urban growth on a local scale, but none of the zoning policies was found to minimize the impact on a broader base. The results from this study may assist the planning and reviewing of land use allocation within Guwahati city to secure ecological sustainability of the wetlands.     

Using a stochastic gradient boosting algorithm to analyse the effectiveness of Landsat 8 data for montado land cover mapping: Application in southern Portugal

This study aims to develop and propose a methodological approach for montado ecosystem mapping using Landsat 8 multi-spectral data, vegetation indices, and the Stochastic Gradient Boosting (SGB) algorithm. Two Landsat 8 scenes (images from spring and summer 2014) of the same area in southern Portugal were acquired. Six vegetation indices were calculated for each scene: the Enhanced Vegetation Index (EVI), the Short-Wave Infrared Ratio (SWIR32), the Carotenoid Reflectance Index 1 (CRI1), the Green Chlorophyll Index (CI_green), the Normalised Multi-band Drought Index (NMDI), and the Soil-Adjusted Total Vegetation Index (SATVI). Based on this information, two datasets were prepared: (i) Dataset I only included multi-temporal Landsat 8 spectral bands (LS8), and (ii) Dataset II included the same information as Dataset I plus vegetation indices (LS8 + VIs). The integration of the vegetation indices into the classification scheme resulted in a significant improvement in the accuracy of Dataset II’s classifications when compared to Dataset I (McNemar test: Z-value = 4.50), leading to a difference of 4.90% in overall accuracy and 0.06 in the Kappa value. For the montado ecosystem, adding vegetation indices in the classification process showed a relevant increment in producer and user accuracies of 3.64% and 6.26%, respectively. By using the variable importance function from the SGB algorithm, it was found that the six most prominent variables (from a total of 24 tested variables) were the following: EVI_summer; CRI1_spring; SWIR32_spring; B6_summer; B5_summer; and CI_green_summer.