Examining Land Cover and Greenness Dynamics in Hangzhou Bay in 1985–2016 Using Landsat Time-Series Data

Examining Vegetation Change and Associated Spatial Patterns in Wuyishan National Park at Different Protection Levels

Examining the characteristics of vegetation change and associated spatial patterns under different protection levels can provide a scientific basis for national park protection and management. Based on the dense time-series Landsat enhanced vegetation index (EVI) data between 1986 and 2020, we utilized the Wild Binary Segmentation (WBS) approach to detect spatial and temporal characteristics of abrupt, gradual, and total changes in Wuyishan National Park. The differences in vegetation change in three protection-level areas (strictly protected [Prots], generally protected [Prot], and non-protected [NP]) were examined, and the contributions to their spatial patterns were evaluated through Geodetector. The results showed the following: (1) The highest percentage of area without abrupt change was in Prots (39.89%), and the lowest percentage was in NP (17.44%). The percentage of abrupt change frequency (larger than three times) increased from 4.40% to 9.10% and 12.49% with the decreases in protection. The significance test showed that the difference in changed frequencies was not significant among these regions, but the interannual variation of abrupt change in Prots was significantly different from other areas. (2) The vegetation coverage of the Wuyishan National Park generally improved. The total EVI change (TEVI) showed that the positive percentage of Prots and Prot was 90.43% and 91.71%, respectively, slightly higher than that of NP (88.44%). However, the mean greenness change of NP was higher than that of Prots and Prot. (3) The park’s EVI spatial pattern in 1986 was the strongest factor determining the EVI spatial pattern in 2020; the explanatory power reduced as the protection level decreased. The explanation power (q value) of abrupt vegetation change was lower and increased as the protection level decreased. The interaction detection showed that EVI1986 and TEVI had the strongest explanatory powers, but the explanatory ability gradually weakened from 0.713 to 0.672 to 0.581 in Prots, Prot, and NP, respectively. This study provided a systematic analysis of vegetation changes and their impacts on spatial patterns.

Quantifying Urban Vegetation Dynamics from a Process Perspective Using Temporally Dense Landsat Imagery

Urban vegetation can be highly dynamic due to the complexity of different anthropogenic drivers. Quantifying such dynamics is crucially important as it is a prerequisite to understanding its social and ecological consequences. Previous studies have mostly focused on the urban vegetation dynamics through monotonic trends analysis in certain intervals, but not considered the process which provides important insights for urban vegetation management. Here, we developed an approach that integrates trends with dynamic analysis to measure the vegetation dynamics from the process perspective based on the time-series Landsat imagery and applied it in Shenzhen, a coastal megacity in southern China, as an example. Our results indicated that Shenzhen was turning green from 2000–2020, even though a large-scale urban expansion occurred during this period. Approximately half of the city (49.5%) showed consistent trends in greening, most of which were located in the areas within the ecological protection baseline. We also found that 35.3% of the Shenzhen city experienced at least a one-time change in urban greenness that was mostly caused by changes in land cover types (e.g., vegetation to developed land). Interestingly, 61.5% of these lands showed trends in greening in the recent change period and most of them were distributed in build-up areas. Our approach that integrates trends analysis and dynamic process reveals information that cannot be discovered by monotonic trends analysis alone, and such information can provide insights for urban vegetation planning and management.

Mangrove Forest Cover and Phenology with Landsat Dense Time Series in Central Queensland, Australia

Wetlands are one of the most biologically productive ecosystems. Wetland ecosystem services, ranging from provision of food security to climate change mitigation, are enormous, far outweighing those of dryland ecosystems per hectare. However, land use change and water regulation infrastructure have reduced connectivity in many river systems and with floodplain and estuarine wetlands. Mangrove forests are critical communities for carbon uptake and storage, pollution control and detoxification, and regulation of natural hazards. Although the clearing of mangroves in Australia is strictly regulated, Great Barrier Reef catchments have suffered landscape modifications and hydrological alterations that can kill mangroves. We used remote sensing datasets to investigate land cover change and both intra- and inter-annual seasonality in mangrove forests in a large estuarine region of Central Queensland, Australia, which encompasses a national park and Ramsar Wetland, and is adjacent to the Great Barrier Reef World Heritage site. We built a time series using spectral, auxiliary, and phenology variables with Landsat surface reflectance products, accessed in Google Earth Engine. Two land cover classes were generated (mangrove versus non-mangrove) in a Random Forest classification. Mangroves decreased by 1480 hectares (−2.31%) from 2009 to 2019. The overall classification accuracies and Kappa coefficient for 2008–2010 and 2018–2020 land cover maps were 95% and 95%, respectively. Using an NDVI-based time series we examined intra- and inter-annual seasonality with linear and harmonic regression models, and second with TIMESAT metrics of mangrove forests in three sections of our study region. Our findings suggest a relationship between mangrove growth phenology along with precipitation anomalies and severe tropical cyclone occurrence over the time series. The detection of responses to extreme events is important to improve understanding of the connections between climate, extreme weather events, and biodiversity in estuarine and mangrove ecosystems.

Impacts of Urban Expansion Forms on Ecosystem Services in Urban Agglomerations: A Case Study of Shanghai-Hangzhou Bay Urban Agglomeration

Exploring impacts of urban expansion on ecosystem services has become a hot topic for regional sustainable development, while analyzing the ecological effects of urban expansion forms under different expansion intensities and city sizes is relatively rare. Therefore, taking a typical urban agglomeration, Shanghai-Hangzhou Bay Urban Agglomeration, as a case study, this study first analyzed the dynamics of urban expansion forms (leapfrogging, edge-expansion, and infilling) and four critical ecosystem services (carbon sequestration, food supply, habitat quality, and soil retention) in three periods from 1990 to 2019. The multiple linear regression model and zonal statistics analysis model were used to quantitatively identify the impacts of urban expansion forms on ecosystem services, taking into account different expansion intensities and city sizes. The results showed that the urban expansion trend in the study area experienced a morphological change from integration to diffusion and then to integration in 1990–2019; edge-expansion was the dominant expansion form. Food supply decreased continuously while other ecosystem services had fluctuating changes, and they all had spatial heterogeneity. The leapfrogging, edge-expansion, and infilling all had negative impacts on ecosystem services, and among them, the edge-expansion intensity had the highest influence degree in the early expansion, and the leapfrogging intensity occupied the dominant position in all influences with the expansion of urban scales. For different city sizes, the impact of edge-expansion in large-scale cities was greater than in small-scale cities in the early expansion, and the impact of leapfrogging in large-scale cities exceeded the edge-expansion in the subsequent expansion. These findings will help further understand the influential mechanisms between urban expansion and ecosystem services and provide a scientific basis for formulating reasonable urban planning.

Decision-Level and Feature-Level Integration of Remote Sensing and Geospatial Big Data for Urban Land Use Mapping

Information about urban land use is important for urban planning and sustainable development. The emergence of geospatial big data (GBD), increased the availability of remotely sensed (RS) data and the development of new methods for data integration to provide new opportunities for mapping types of urban land use. However, the modes of RS and GBD integration are diverse due to the differences in data, study areas, classifiers, etc. In this context, this study aims to summarize the main methods of data integration and evaluate them via a case study of urban land use mapping in Hangzhou, China. We first categorized the RS and GBD integration methods into decision-level integration (DI) and feature-level integration (FI) and analyzed their main differences by reviewing the existing literature. The two methods were then applied for mapping urban land use types in Hangzhou city, based on urban parcels derived from the OpenStreetMap (OSM) road network, 10 m Sentinel-2A images, and points of interest (POI). The corresponding classification results were validated quantitatively and qualitatively using the same testing dataset. Finally, we illustrated the advantages and disadvantages of both approaches via bibliographic evidence and quantitative analysis. The results showed that: (1) The visual comparison indicates a generally better performance of DI-based classification than FI-based classification; (2) DI-based urban land use mapping is easy to implement, while FI-based land use mapping enables the mixture of features; (3) DI-based and FI-based methods can be used together to improve urban land use mapping, as they have different performances when classifying different types of land use. This study provides an improved understanding of urban land use mapping in terms of the RS and GBD integration strategy.

Impacts of the Tropical Cyclone Idai in Mozambique: A Multi-Temporal Landsat Satellite Imagery Analysis

The Central Region of Mozambique (Sofala Province) bordering on the active cyclone area of the southwestern Indian Ocean has been particularly affected by climate hazards. The Cyclone Idai, which hit the region in March 2019 with strong winds causing extensive flooding and a massive loss of life, was the strongest recorded tropical cyclone in the Southern Hemisphere. The aim of this study was to use pre- and post-cyclone Idai Landsat satellite images to analyze temporal changes in Land Use and Land Cover (LULC) across the Sofala Province. Specifically, we aimed—(i) to quantify and map the changes in LULC between 2012 and 2019; (ii) to investigate the correlation between the distance to Idai’s trajectory and the degree of vegetation damage, and (iii) to determine the damage caused by Idai on different LULC. We used Landsat 7 and 8 images (with 30 m resolution) taken during the month of April for the 8-year period. The April Average Normalized Difference Vegetation Index (NDVI) over the aforementioned period (2012–2018, pre-cyclone) was compared with the values of April 2019 (post-cyclone). The results showed a decreasing trend of the productivity (NDVI 0.5 to 0.8) and an abrupt decrease after the cyclone. The most devastated land use classes were dense vegetation (decreased by 59%), followed by wetland vegetation (−57%) and shrub land (−56%). The least damaged areas were barren land (−23%), barren vegetation (−27%), and grassland and dambos (−27%). The Northeastern, Central and Southern regions of Sofala were the most devastated areas. The Pearson Correlation Coefficient between the relative vegetation change activity after Idai (NDVI%) and the distance to Idai’s trajectory was 0.95 (R-square 0.91), suggesting a strong positive linear correlation. Our study also indicated that the LULC type (vegetation physiognomy) might have influenced the degree of LULC damage. This study provides new insights for the management and conservation of natural habitats threatened by climate hazards and human factors and might accelerate ongoing recovery processes in the Sofala Province.

Examining Water Area Changes Accompanying Dam Construction in the Madeira River in the Brazilian Amazon

Two recently constructed run-of-the-river dams (Santo Antônio and Jirau), along the Madeira River in Brazil, have been controversial due to their large unquantified impacts on (1) land use and land cover (LULC) and (2) on the area that would be flooded. Based on annual LULC data from 1985 to 2017, this study integrated intensity analysis and difference components methods to analyze the impacts of the two dams on the annual flooded area in upstream, midstream, and downstream regions of the Madeira River. The dam construction significantly influenced LULC change intensity in the upstream and midstream regions since 2011 and 2010, respectively. An increase of 18.5% of the newly flooded area (462.58 km2) in the post-dam construction period was observed. The water gross gain intensity was active during 2011–2017 and 2011–2014 in upstream and midstream, respectively. The dominant difference components of water change were exchanged in the pre-dam period and became quantity in the post-dam period for both upstream and midstream regions. Forest was the major land category replaced by water; however, the highest gain intensities occurred in other non-vegetated areas in upstream and midstream. This study provided a useful approach for characterizing impacts of dam construction on water area change.

Remote Sensing of Mangroves and Estuarine Communities in Central Queensland, Australia

Great Barrier Reef catchments are under pressure from the effects of climate change, landscape modifications, and hydrology alterations. With the use of remote sensing datasets covering large areas, conventional methods of change detection can expose broad transitions, whereas workflows that excerpt data for time-series trends divulge more subtle transformations of land cover modification. Here, we combine both these approaches to investigate change and trends in a large estuarine region of Central Queensland, Australia, that encompasses a national park and is adjacent to the Great Barrier Reef World Heritage site. Nine information classes were compiled in a maximum likelihood post classification change analysis in 2004–2017. Mangroves decreased (1146 hectares), as was the case with estuarine wetland (1495 hectares), and saltmarsh grass (1546 hectares). The overall classification accuracies and Kappa coefficient for 2004, 2006, 2009, 2013, 2015, and 2017 land cover maps were 85%, 88%, 88%, 89%, 81%, and 92%, respectively. The cumulative area of open forest, estuarine wetland, and saltmarsh grass (1628 hectares) was converted to pasture in a thematic change analysis showing the “from–to” change. We generated linear regression relationships to examine trends in pixel values across the time series. Our findings from a trend analysis showed a decreasing trend (p value range = 0.001–0.099) in the vegetation extent of open forest, fringing mangroves, estuarine wetlands, saltmarsh grass, and grazing areas, but this was inconsistent across the study site. Similar to reports from tropical regions elsewhere, saltmarsh grass is poorly represented in the national park. A severe tropical cyclone preceding the capture of the 2017 Landsat 8 Operational Land Imager (OLI) image was likely the main driver for reduced areas of shoreline and stream vegetation. Our research contributes to the body of knowledge on coastal ecosystem dynamics to enable planning to achieve more effective conservation outcomes.

Sparse Unmixing for Hyperspectral Image with Nonlocal Low-Rank Prior

Hyperspectral unmixing is a key preprocessing technique for hyperspectral image analysis. To further improve the unmixing performance, in this paper, a nonlocal low-rank prior associated with spatial smoothness and spectral collaborative sparsity are integrated together for unmixing the hyperspectral data. The proposed method is based on a fact that hyperspectral images have self-similarity in nonlocal sense and smoothness in local sense. To explore the spatial self-similarity, nonlocal cubic patches are grouped together to compose a low-rank matrix. Then, based on the linear mixed model framework, the nuclear norm is constrained to the abundance matrix of these similar patches to enforce low-rank property. In addition, the local spatial information and spectral characteristic are also taken into account by introducing TV regularization and collaborative sparse terms, respectively. Finally, the results of the experiments on two simulated data sets and two real data sets show that the proposed algorithm produces better performance than other state-of-the-art algorithms.

Ecosystem health monitoring in the Shanghai-Hangzhou Bay Metropolitan Area: A hidden Markov modeling approach

Ecosystem health assessment is an important method for obtaining information on ecosystem conditions, and it plays a vital role in preserving and enhancing ecosystem health status. In addition, it provides useful information and knowledge for urban agglomeration development decision makers. However, ecological phenomena often vary considerably from one observation to the next, which makes it difficult to distinguish different status of the ecosystem health. In this study, hidden Markov model (HMM) was employed to simulate the internal-external correlations of ecosystem status through establishing the relationships between internal ecological health level and combination state of external observation. Based on the statistics and land use data in 2001, 2007 and 2013, the Vigor-Organization-Resilience (VOR) framework was employed to identify the ecosystem health in Shanghai-Hangzhou Bay Metropolitan (SHBM), in which the ecosystem health state was considered as a hidden state that could be estimated according to the conditions of vigor, organization and resilience. In addition, two parameter learning cases including mathematical statistics and extensible sequence method were employed to solve the iterative convergence problem of parameters in short-time series of ecosystem health simulation. Results show that HMM not only provides a comparable descriptive ability to that of the VOR model, but also can monitor ecosystem health at the optimal grid scale in SHBM. The combination of HMM and VOR greatly expands the spatiotemporal characteristics and provides a new research approach for the study of ecosystem health assessment of urban agglomerations.

Normalized Difference Vegetation Index Continuity of the Landsat 4-5 MSS and TM: Investigations Based on Simulation

Landsat 4-5, built at the same time and with the same design, carrying the Multispectral Scanner System (MSS) and the Thematic Mapper (TM) simultaneously, jointly provided observation service for about 30 years (1982–2013). Considering the importance of data continuity for time series analyses, investigations on the continuity of the Landsat 4-5 MSS and TM are required. In this paper, characterization differences between the Landsat 4-5 MSS and TM were initially discussed using the synthesized reflectance records generated from a collection of Hyperion hyperspectral profiles which were well calibrated and widely distributed. The difference in near-infrared region mostly contributed to the difference in normalized difference vegetation index (NDVI) between MSS and TM, while the between-sensor difference in red spectrum was relatively minor. Models for transforming MSS NDVI to TM NDVI were proposed, and validated subsequently through cross-validation tests. Furthermore, effectiveness of the transformation models was investigated using eight synchronous observation pairs of the Landsat 5 MSS and TM. On average, the univariate models through ordinary least squares regression (OLS) regression resulted in a decrease about 10% of the median relative difference (MdRD). Meanwhile, the bivariate models improved the NDVI comparability in most cases, especially when the transformation models through ridge regression were implemented. The univariate model through OLS regression could be the only solution for cases when problems of data quality are encountered (e.g., problem in the MSS near-infrared channel (800–1000 nm)). In conclusion, the findings on NDVI transformation models from MSS to TM are valuable for reference, because of the collection of diverse Hyperion hyperspectral profiles used.

Editorial for the Special Issue “Remote Sensing in Coastal Zone Monitoring and Management—How Can Remote Sensing Challenge the Broad Spectrum of Temporal and Spatial Scales in Coastal Zone Dynamic?”

Coastal zones are sensitive areas responding at various scales (events to long-term trends) where the monitoring and management of physico-chemical, biological, morphological processes, and fluxes are highly challenging [...]

A Model for Estimating the Vegetation Cover in the High-Altitude Wetlands of the Andes (HAWA)

The natural salt meadows of Tilopozo in the hyperarid, Atacama Desert of northern Chile, which are located at approximately 2800 m above sea level, are under pressure from industrial activity, and cultivation and grazing by local communities. In this research, the land surface covered by salt meadow vegetation was estimated from normalized difference vegetation indices (NDVI) derived from Landsat Thematic Mapper (TM), Enhanced Thematic Mapper (ETM+) and Operational Land Imager (OLI) data from 1985 to 2016. The vegetated area of the Tilopozo salt meadows decreased by 34 ha over the 32-year period studied. Multiple regression models of the area covered by vegetation and climate data and groundwater depths were derived on an annual basis, as well as for both the dry and wet seasons and had R2 values of 83.0%, 72.8% and 92.4% respectively between the vegetated areas modeled and those estimated from remotely sensed data. These models are potentially useful tools for studies into the conservation of the Tilopozo salt meadows, as they provide relevant information on the state of vegetation and enable changes in vegetation in response to fluctuations in climate parameters and groundwater depths to be predicted.