[go: up one dir, main page]
More Web Proxy on the site http://driver.im/
You seem to have javascript disabled. Please note that many of the page functionalities won't work as expected without javascript enabled.
 
 
Sign in to use this feature.

Years

Between: -

Subjects

remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline

Journals

Article Types

Countries / Regions

remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline

Search Results (295)

Search Parameters:
Keywords = PS-InSAR

Order results
Result details
Results per page
Select all
Export citation of selected articles as:
17 pages, 12868 KiB  
Article
PSInSAR-Based Time-Series Coastal Deformation Estimation Using Sentinel-1 Data
by Muhammad Ali, Alessandra Budillon, Zeeshan Afzal, Gilda Schirinzi and Sajid Hussain
Land 2025, 14(3), 536; https://doi.org/10.3390/land14030536 - 4 Mar 2025
Viewed by 119
Abstract
Coastal areas are highly dynamic regions where surface deformation due to natural and anthropogenic activities poses significant challenges. Synthetic Aperture Radar (SAR) interferometry techniques, such as Persistent Scatterer Interferometry (PSInSAR), provide advanced capabilities to monitor surface deformation with high precision. This study applies [...] Read more.
Coastal areas are highly dynamic regions where surface deformation due to natural and anthropogenic activities poses significant challenges. Synthetic Aperture Radar (SAR) interferometry techniques, such as Persistent Scatterer Interferometry (PSInSAR), provide advanced capabilities to monitor surface deformation with high precision. This study applies PSInSAR techniques to estimate surface deformation along coastal zones from 2017 to 2020 using Sentinel-1 data. In the densely populated areas of Pasni, an annual subsidence rate of 130 mm is observed, while the northern, less populated region experiences an uplift of 70 mm per year. Seawater intrusion is an emerging issue causing surface deformation in Pasni’s coastal areas. It infiltrates freshwater aquifers, primarily due to excessive groundwater extraction and rising sea levels. Over time, seawater intrusion destabilizes the underlying soil and rock structures, leading to subsidence or gradual sinking of the ground surface. This form of surface deformation poses significant risks to infrastructure, agriculture, and the local ecosystem. Land deformation varies along the study area’s coastline. The eastern region, which is highly reclaimed, is particularly affected by erosion. The results derived from Sentinel-1 SAR data indicate significant subsidence in major urban districts. This information is crucial for coastal management, hazard assessment, and planning sustainable development in the region. Full article
(This article belongs to the Special Issue Assessing Land Subsidence Using Remote Sensing Data)
Show Figures

Figure 1

Figure 1
<p>Location map of the study area (using ArcMap software 10.8 and base map): (<b>a</b>) geographical location of Pakistan; (<b>b</b>) coastline of Pakistan where the white small rectangle indicates the location of the study area along the coastline, the yellow rectangle indicates the Sentinel-1 IW swath, and the red solid line indicates the Makran Subduction Zone, with white arrows showing the relative plate motion; (<b>c</b>) study area with red circular points indicating the major locations along the coast.</p>
Full article ">Figure 2
<p>Geological map of the study area (using ArcMap software) where the lithology of the study area is classified into different colored polygons (@Survey of Pakistan).</p>
Full article ">Figure 3
<p>The line of sight (LOS) velocity map depicts the descending and ascending swath mood of Sentinel-1 data from January 2017 to December 2020.</p>
Full article ">Figure 4
<p>Deformation velocity overlaid on Google Earth (GE) imagery showing the displacement trend of four different locations: Zone A is a zoomed view of the Shadi Khor River Bed, Zone B shows a VLOS overlaid view of Pasni Grid Station, Zone C is an image of Pasni Residential Area, Pasni Harbor, and Pakistan Coast Guard, and Zone D depicts the residential area on the Pasni Jetty.</p>
Full article ">Figure 5
<p>Time-series deformation results for the years (<b>a</b>) 2017, (<b>b</b>) 2018, (<b>c</b>) 2019, and (<b>d</b>) 2020 along the coast of Pasni derived using PSInSAR analysis of Sentinel-1 data. The maps illustrate surface deformation velocity along the LOS, which are shown in red to blue hues.</p>
Full article ">Figure 6
<p>Change in shorelines from 2010 to 2020 overlaid on the latest image using ArcMap software and a base map.</p>
Full article ">Figure 7
<p>Time-series shoreline analysis of the Pasni coastline from 2010 to 2020 showing shoreline changes across the intervals (<b>a</b>) 2010–2014 (red and black lines), (<b>b</b>) 2014–2017 (black and blue lines), (<b>c</b>) 2017–2020 (blue and green lines), and (<b>d</b>) the entire period from 2010 to 2020 (red and green lines).</p>
Full article ">Figure 8
<p>Quantitative analysis of time-series shoreline changes from the year 2010 to the year 2020 in four different zones (A, B, C, and D) of the study area.</p>
Full article ">Figure 9
<p>The relationship between surface deformation from 2017 to 2020 and the change in shorelines from 2010 to 2020 overlaid on the latest image using ArcMap software and a base map.</p>
Full article ">Figure 10
<p>Quantitative results for land use land cover categories and classification (barel land, mud area, built-up area, vegetation, and water) over a ten-year period from the years 2010 to 2020 in Pasni coastal areas.</p>
Full article ">
26 pages, 15590 KiB  
Article
Technical and Policy Analysis: Time Series of Land Subsidence for the Evaluation of the Jakarta Groundwater-Free Zone
by Joko Widodo, Edy Trihatmoko, Nugraheni Setyaningrum, Yuta Izumi, Rendi Handika, Mohammad Ardha, Rahmat Arief, Shinichi Sobue, Nurlinda Nurlinda, Pulung Arya Pranantya, Jovi Rauhillah Wiranu and Muhammad Rokhis Khomarudin
Urban Sci. 2025, 9(3), 67; https://doi.org/10.3390/urbansci9030067 - 4 Mar 2025
Viewed by 251
Abstract
Jakarta faces a critical challenge of extensive land subsidence, ranking prominently globally. This research employs a combined technical and policy evaluation approach to analyze the issue, incorporating sustainability considerations to assess the efficacy of Governor Regulation of Jakarta Number 93 of 2021, focusing [...] Read more.
Jakarta faces a critical challenge of extensive land subsidence, ranking prominently globally. This research employs a combined technical and policy evaluation approach to analyze the issue, incorporating sustainability considerations to assess the efficacy of Governor Regulation of Jakarta Number 93 of 2021, focusing on how the groundwater-free zone relates to land subsidence in the city. We processed 81 ALOS-2 PALSAR-2 synthetic aperture radar (SAR) data using persistent scatterer interferometric synthetic aperture radar (PS-InSAR) with HH polarization from 2017 to 2022 and ground truthing with 255 global positioning system (GPS) real-time kinematic (RTK) validation points. Our findings reveal a significant misalignment in the designated groundwater-free zone in the central part of Jakarta. At the same time, severe land subsidence primarily affects northern and northwestern Jakarta, with an average land subsidence rate of 5–6 cm/year. We strongly advocate for a thorough evaluation to rectify and redefine the boundaries of groundwater-free zones, improve regulatory frameworks, and effectively address land subsidence mitigation in the study area. The impact of domestic water needs on land subsidence highlights the urgency of action. Based on a combination of land subsidence velocity rates and domestic water demand, we have classified the cities in Jakarta into three levels of recommendations for groundwater-free zones. The cities are ranked in order of priority from highest to lowest: (1) West Jakarta, (2) North Jakarta, (3) South Jakarta, (4) East Jakarta, and (5) Central Jakarta, which holds the lowest priority. Full article
Show Figures

Figure 1

Figure 1
<p>The study area of the research.</p>
Full article ">Figure 2
<p>Data statistics and normal baselines of ALOS-2 PALSAR-2 data used in this research (81 images).</p>
Full article ">Figure 2 Cont.
<p>Data statistics and normal baselines of ALOS-2 PALSAR-2 data used in this research (81 images).</p>
Full article ">Figure 3
<p>Research flow.</p>
Full article ">Figure 4
<p>Interferogram of 81 ALOS 2 images and PS-InSAR result for Jakarta and surrounding area.</p>
Full article ">Figure 5
<p>Nash–Sutcliffe efficiency model results with a 0.8 (outstanding) value between PS-InSAR as a model and GNSS data as an observation (Obs).</p>
Full article ">Figure 6
<p>GNSS in conjunction with PS-InSAR data within a 100 m buffer.</p>
Full article ">Figure 7
<p>Jakarta land subsidence, 2017–2022.</p>
Full article ">Figure 8
<p>The average of the velocity subsidence rate (mm) in each city of Jakarta from 2017 to 2022. Interferogram of 81 ALOS 2 images.</p>
Full article ">Figure 9
<p>Time series land subsidence and DWNs (in liters) at the research location.</p>
Full article ">Figure 10
<p>Comparison of annual domestic water needs (DWNs) in the five administrative cities of Jakarta.</p>
Full article ">Figure 11
<p>Groundwater-free zone based on Governor Regulation Number 93 of 2021 with land subsidence rate/year (mm) in Jakarta with 81 scene baselines of ALOS PALSAR data and DWNs (in liters) during the acquisition period from 12 June 2017–06 June 2020.</p>
Full article ">Figure 12
<p>Zoning recommendation map for groundwater-free zones considering land subsidence and domestic water needs in Jakarta.</p>
Full article ">Figure 13
<p>The extent of recommended groundwater-free areas in Jakarta (Ha) based on land subsidence analysis and city-level water demand.</p>
Full article ">
17 pages, 134434 KiB  
Technical Note
The Influence of Parameter Estimation Bounds on Velocity Estimation in PSInSAR
by Timo Balz and Mostafa Ewais
Remote Sens. 2025, 17(5), 788; https://doi.org/10.3390/rs17050788 - 24 Feb 2025
Viewed by 180
Abstract
Permanent scatterer interferometric synthetic aperture radar (PSInSAR) processing requires parameter selection that can significantly impact results, yet these parameters are often not fully disclosed in scientific publications. To demonstrate how different parameter settings in PSInSAR processing affect results, our study analyzed PSInSAR processing [...] Read more.
Permanent scatterer interferometric synthetic aperture radar (PSInSAR) processing requires parameter selection that can significantly impact results, yet these parameters are often not fully disclosed in scientific publications. To demonstrate how different parameter settings in PSInSAR processing affect results, our study analyzed PSInSAR processing results using varying parameters. Results were evaluated both with and without temporal coherence filtering (threshold ≥ 0.8). Parameter variations produced differences that exceeded previously stated accuracy ranges for PSInSAR methods, while overall deformation trends remained similar across parameter sets. This shows that even seemingly minor parameter variations can lead to significant differences in PSInSAR results, exceeding what would be considered acceptable with respect to previously published accuracies. These findings emphasize the need for complete parameter disclosure in scientific publications and suggest more careful interpretation of small velocity differences in PSInSAR results. Full article
Show Figures

Figure 1

Figure 1
<p>Study area in Wuhan, China (<b>left</b>) and Naples, Italy (<b>right</b>); the red rectangle marks the extent of the TerraSAR-X stack used in this study, and the blue rectangle is the study area. The background image was provided by ESRI.</p>
Full article ">Figure 2
<p>Velocity estimations over Wuhan for parameters <math display="inline"><semantics> <mrow> <mo>±</mo> <mn>40</mn> </mrow> </semantics></math> <math display="inline"><semantics> <mi mathvariant="normal">m</mi> </semantics></math><math display="inline"><semantics> <mi mathvariant="normal">m</mi> </semantics></math>/<math display="inline"><semantics> <mi>yr</mi> </semantics></math> and <math display="inline"><semantics> <mrow> <mo>±</mo> <mn>50</mn> </mrow> </semantics></math> <math display="inline"><semantics> <mi mathvariant="normal">m</mi> </semantics></math> (<b>left</b>) <math display="inline"><semantics> <mrow> <mo>±</mo> <mn>100</mn> </mrow> </semantics></math> <math display="inline"><semantics> <mi mathvariant="normal">m</mi> </semantics></math><math display="inline"><semantics> <mi mathvariant="normal">m</mi> </semantics></math>/<math display="inline"><semantics> <mi>yr</mi> </semantics></math> and <math display="inline"><semantics> <mrow> <mo>±</mo> <mn>200</mn> </mrow> </semantics></math> <math display="inline"><semantics> <mi mathvariant="normal">m</mi> </semantics></math> (<b>right</b>); the red rectangle on the left marks the subset shown in <a href="#remotesensing-17-00788-f004" class="html-fig">Figure 4</a>; background image was provided by ESRI.</p>
Full article ">Figure 3
<p>Velocity estimations over Napoli for parameters <math display="inline"><semantics> <mrow> <mo>±</mo> <mn>40</mn> </mrow> </semantics></math> <math display="inline"><semantics> <mi mathvariant="normal">m</mi> </semantics></math><math display="inline"><semantics> <mi mathvariant="normal">m</mi> </semantics></math>/<math display="inline"><semantics> <mi>yr</mi> </semantics></math> and <math display="inline"><semantics> <mrow> <mo>±</mo> <mn>50</mn> </mrow> </semantics></math> <math display="inline"><semantics> <mi mathvariant="normal">m</mi> </semantics></math> (<b>left</b>) <math display="inline"><semantics> <mrow> <mo>±</mo> <mn>100</mn> </mrow> </semantics></math> <math display="inline"><semantics> <mi mathvariant="normal">m</mi> </semantics></math><math display="inline"><semantics> <mi mathvariant="normal">m</mi> </semantics></math>/<math display="inline"><semantics> <mi>yr</mi> </semantics></math> and <math display="inline"><semantics> <mrow> <mo>±</mo> <mn>200</mn> </mrow> </semantics></math> <math display="inline"><semantics> <mi mathvariant="normal">m</mi> </semantics></math> (<b>right</b>); the red rectangle on the left marks the subset shown in <a href="#remotesensing-17-00788-f005" class="html-fig">Figure 5</a>; background image was provided by ESRI.</p>
Full article ">Figure 4
<p>Velocity estimation (unfiltered) at a subset of the Wuhan test area with different parameters: from ±50 m to ±200 <math display="inline"><semantics> <mi mathvariant="normal">m</mi> </semantics></math> residual height estimation in columns and from ±40 mm/yr to ±100 <math display="inline"><semantics> <mi mathvariant="normal">m</mi> </semantics></math><math display="inline"><semantics> <mi mathvariant="normal">m</mi> </semantics></math>/<math display="inline"><semantics> <mi>yr</mi> </semantics></math> linear velocity in rows.</p>
Full article ">Figure 5
<p>Velocity estimation (unfiltered) at a subset of the Naples test area with different parameters: from ±50 m to ±200 <math display="inline"><semantics> <mi mathvariant="normal">m</mi> </semantics></math> residual height estimation in columns and from ±40 mm/yr to ±100 <math display="inline"><semantics> <mi mathvariant="normal">m</mi> </semantics></math><math display="inline"><semantics> <mi mathvariant="normal">m</mi> </semantics></math>/<math display="inline"><semantics> <mi>yr</mi> </semantics></math> linear velocity in rows.</p>
Full article ">Figure 6
<p>Scatter plot showing the relation between the temporal coherence and the estimated velocity for parameters −60 mm/yr to 60 mm/yr linear velocity and −100 m to 100 m residual height for Wuhan (<b>left</b>) and Napoli (<b>right</b>).</p>
Full article ">Figure 7
<p>Scatter plots for different parameters showing the temporal coherence in relation to the estimated parameters.</p>
Full article ">Figure 8
<p>Scatter plots for different parameters showing the temporal coherence in relation to the estimated parameters in Napoli.</p>
Full article ">Figure 9
<p>Range of velocity differences in mm/yr, i.e., maximum velocity in the vector of 16 results per PS point minus the minimum velocity in the vector, between results from different parameters for the Wuhan test area.</p>
Full article ">Figure 10
<p>Range of velocity differences in mm/yr, i.e., maximum velocity in the vector of 16 results per PS point minus the minimum velocity in the vector, between results from different parameters for the Napels test area.</p>
Full article ">
21 pages, 1385 KiB  
Article
The New Occurrence of Antiphospholipid Syndrome in Severe COVID-19 Cases with Pneumonia and Vascular Thrombosis Could Explain the Post-COVID Syndrome
by Mirjana Zlatković-Švenda, Melanija Rašić, Milica Ovuka, Slavica Pavlov-Dolijanović, Marija Atanasković Popović, Manca Ogrič, Polona Žigon, Snežna Sodin-Šemrl, Marija Zdravković and Goran Radunović
Biomedicines 2025, 13(2), 516; https://doi.org/10.3390/biomedicines13020516 - 19 Feb 2025
Viewed by 398
Abstract
Introduction: The classification of antiphospholipid syndrome (APS) comprises clinical criteria (vascular thrombosis or obstetric complications throughout life) and laboratory criteria (antiphospholipid antibodies (aPLs) positivity, confirmed at least twice at 12-week interval). Methods: In 100 patients admitted to the hospital with COVID-19 pneumonia, thrombosis [...] Read more.
Introduction: The classification of antiphospholipid syndrome (APS) comprises clinical criteria (vascular thrombosis or obstetric complications throughout life) and laboratory criteria (antiphospholipid antibodies (aPLs) positivity, confirmed at least twice at 12-week interval). Methods: In 100 patients admitted to the hospital with COVID-19 pneumonia, thrombosis and pregnancy complications were recorded during the hospital stay and in personal medical history. They were tested for nine types of aPLs at four time points (admission, deterioration, discharge, and 3-month follow-up): anticardiolipin (aCL), anti-β2-glycoproteinI (anti-β2GPI), and antiphosphatidylserine/prothrombin (aPS/PT) isotypes IgM/IgG/IgA. Results: During hospitalization, aPLs were detected at least once in 51% of patients. All 7% of deceased patients tested negative for aPLs upon admission, and only one patient became aCL IgG positive as his condition worsened. In 83.3% of patients, intrahospital thrombosis was not related to aPLs. One patient with pulmonary artery and cerebral artery thrombosis was given an APS diagnosis (triple aPLs positivity on admission, double on follow-up). Personal anamnesis (PA) for thromboembolism was verified in 10 patients, all of whom tested negative for aPLs at admission; however, transition to aPLs positivity at discharge (as the disease subsided) was seen in 60% of patients: three of six with arterial thrombosis (at follow-up, two did not appear, and one was negativized) and three of four with deep vein thrombosis (one was confirmed at follow-up and diagnosed with APS, one was negativized, and one did not appear). At admission, the majority of the aPLs were of the aCL IgG class (58.8%). Unexpectedly, as the COVID-19 disease decreased, anti-β2GPI IgG antibodies (linked with thromboses) became newly positive at discharge (14.9%), as confirmed at follow-up (20.8%). Conclusion: The incidence of APS in our cohort was 2.0%, whereas in the general population, it ranges from 0.001% to 0.002%. The incidence might have increased even more if the four aPLs-positive patients with intrahospital thrombosis/history of thrombosis had attended follow-up. Recommendation: All patients with severe COVID-19 or post-COVID syndrome should be evaluated for current/previous thrombosis and tested for aPLs at least twice: at admission to the hospital and at discharge, then retested 3 months later in positive cases in order to be given the appropriate therapy. Full article
(This article belongs to the Special Issue Emerging Trends in Pathophysiology and Therapy of COVID-19)
Show Figures

Figure 1

Figure 1
<p>Study flow chart with the inclusion/exclusion criteria.</p>
Full article ">Figure 2
<p>Recommendations for the testing of current or previous COVID-19 patients for antiphospholipid antibodies and for their evaluation with regard to current/past thrombosis with the proposed therapy in positive cases. aPLs, antiphospholipid antibodies; aCL, anticardiolipin antibodies; anti-β2GPI, anti-β2-glycoprotein I antibodies; aPS/PT, anti-phosphatidylserine-prothrombin antibodies; LDA, low-dose aspirin; LMWH, low molecular-weight heparin (Fraxiparine); VKA, vitamin K antagonists (warfarin); DOAC (s), direct oral anticoagulants.</p>
Full article ">
22 pages, 15578 KiB  
Article
Analysis of Ground Subsidence Evolution Characteristics and Attribution Along the Beijing–Xiong’an Intercity Railway with Time-Series InSAR and Explainable Machine-Learning Technique
by Xin Liu, Huili Gong, Chaofan Zhou, Beibei Chen, Yanmin Su, Jiajun Zhu and Wei Lu
Land 2025, 14(2), 364; https://doi.org/10.3390/land14020364 - 10 Feb 2025
Viewed by 343
Abstract
The long-term overextraction of groundwater in the Beijing–Tianjin–Hebei region has led to the formation of the world’s largest groundwater depression cone and the most extensive land subsidence zone, posing a potential threat to the operational safety of high-speed railways in the region. As [...] Read more.
The long-term overextraction of groundwater in the Beijing–Tianjin–Hebei region has led to the formation of the world’s largest groundwater depression cone and the most extensive land subsidence zone, posing a potential threat to the operational safety of high-speed railways in the region. As a critical transportation hub connecting Beijing and the Xiong’an New Area, the Beijing–Xiong’an Intercity Railway traverses geologically complex areas with significant ground subsidence issues. Monitoring and analyzing the causes of land subsidence along the railway are essential for ensuring its safe operation. Using Sentinel-1A radar imagery, this study applies PS-InSAR technology to extract the spatiotemporal evolution characteristics of ground subsidence along the railway from 2016 to 2022. By employing a buffer zone analysis and profile analysis, the subsidence patterns at different stages (pre-construction, construction, and operation) are revealed, identifying the major subsidence cones along the Yongding River, Yongqing, Daying, and Shengfang regions, and their impacts on the railway. Furthermore, the XGBoost model and SHAP method are used to quantify the primary influencing factors of land subsidence. The results show that changes in confined water levels are the most significant factor, contributing 34.5%, with strong interactions observed between the compressible layer thickness and confined water levels. The subsidence gradient analysis indicates that the overall subsidence gradient along the Beijing–Xiong’an Intercity Railway currently meets safety standards. This study provides scientific evidence for risk prevention and the control of land subsidence along the railway and holds significant implications for ensuring the safety of high-speed rail operations. Full article
(This article belongs to the Special Issue Assessing Land Subsidence Using Remote Sensing Data)
Show Figures

Figure 1

Figure 1
<p>Study area and the extent of SAR imagery.</p>
Full article ">Figure 2
<p>Ps-InSAR technology roadmap.</p>
Full article ">Figure 3
<p>XGBoost technology roadmap.</p>
Full article ">Figure 4
<p>Time-series evolution trend of subsidence in the background area along the Beijing–Xiong’an Intercity Railway.</p>
Full article ">Figure 5
<p>Accuracy verification of InSAR monitoring subsidence results with leveling precision.</p>
Full article ">Figure 6
<p>Annual average subsidence rates in the buffer zone along the Beijing–Xiong’an Intercity Railway across different periods.</p>
Full article ">Figure 7
<p>Relationship between annual average subsidence rates along the profile of the Beijing–Xiong’an Intercity Railway and the distribution of subsidence funnels.</p>
Full article ">Figure 8
<p>Annual average subsidence rates of the Beijing–Xiong’an Intercity Railway cross-section during different periods (The dashed lines of different colors represent areas where the sedimentation rate varies greatly).</p>
Full article ">Figure 9
<p>Cumulative subsidence at stations along the Beijing–Xiong’an Intercity Railway.</p>
Full article ">Figure 10
<p>Variations in annual average subsidence rates of stations along the Beijing–Xiong’an Intercity Railway across different periods.</p>
Full article ">Figure 11
<p>Distribution of monitoring points on both sides of stations along the Beijing–Xiong’an Intercity Railway.</p>
Full article ">Figure 12
<p>Subsidence rate differences between the eastern and western sides of Xiong’an Station and Bazhou North Station.</p>
Full article ">Figure 13
<p>Subsidence rate differences between the eastern and western sides of Gu’an East Station, Daxing Airport Station, and Beijing Daxing Station.</p>
Full article ">Figure 14
<p>Accuracy validation based on the Extreme Gradient Boosting (XGBoost) model.</p>
Full article ">Figure 15
<p>Distribution of the importance of subsidence-influencing factors for the Beijing–Xiong’an Intercity Railway.</p>
Full article ">Figure 16
<p>SHAP interpretability analysis of overall subsidence characteristics.</p>
Full article ">Figure 17
<p>SHAP interpretability analysis of factor interactions.</p>
Full article ">Figure 18
<p>Variations in subsidence gradient along the Beijing–Xiong’an Intercity Railway.</p>
Full article ">
25 pages, 5948 KiB  
Article
Spatiotemporal Variability of Groundwater Quality for Irrigation: A Case Study in Mimoso Alluvial Valley, Semiarid Region of Brazil
by Thayná A. B. Almeida, Abelardo A. A. Montenegro, João L. M. P. de Lima, Carolyne W. L. A. Farias, Ailton A. Carvalho and Anderson L. R. de Paiva
Water 2025, 17(3), 410; https://doi.org/10.3390/w17030410 - 1 Feb 2025
Viewed by 789
Abstract
Alluvial aquifers are vital for agricultural communities in semiarid regions, where groundwater quality is often constrained by seasonal and spatial salinity variations. This study employed geostatistical methods to analyze the spatial and temporal variability of electrical conductivity (EC) and the sodium adsorption ratio [...] Read more.
Alluvial aquifers are vital for agricultural communities in semiarid regions, where groundwater quality is often constrained by seasonal and spatial salinity variations. This study employed geostatistical methods to analyze the spatial and temporal variability of electrical conductivity (EC) and the sodium adsorption ratio (SAR) and elaborate an indicative quality map in the Mimoso Alluvial Aquifer, Pernambuco, Brazil. Groundwater samples were collected and analyzed for cations, total hardness (TH), and the percentage of sodium (PS). Moreover, the relation between EC and the SAR was used to determine the groundwater quality for irrigation. Cation concentrations followed the order Ca2+ > Mg2+ > Na+ > K+. EC and the SAR exhibited medium to high variability, with spatial dependence ranging from moderate to strong, and presented a strong cross-spatial dependence. Results showed that sequential Gaussian simulation (SGS) provided a more reliable groundwater classification for agricultural purposes compared to kriging methods, enabling a more rigorous evaluation. Based on the strong geostatistical cross correlation between EC and RAS, a novel water quality index was proposed, properly identifying regions with lower groundwater quality. The resulting spatial indicator maps classified groundwater as suitable (64.7%), restricted use (2.08%) and unsuitable (2.38%) for irrigation. The groundwater quality maps indicated that groundwater was mostly suitable for agriculture, except in silty areas, also corresponding to regions with low hydraulic conductivity at the saturated zone. Soil texture, rainfall, and water extraction significantly influenced spatial and temporal patterns of groundwater quality. Such correlations allow a better understanding of the groundwater quality in alluvial valleys, being highly relevant for water resources management in semiarid areas. Full article
(This article belongs to the Special Issue Advance in Groundwater in Arid Areas)
Show Figures

Figure 1

Figure 1
<p>South America map, with location of the semiarid region in Brazil (<b>A</b>); Pernambuco State, and location of Ipanema River Basin (<b>B</b>); Mimoso Alluvial Valley, Alto Ipanema Basin, Pesqueira municipality of Pernambuco State, Brazil (<b>C</b>); saturated hydraulic conductivity map of Mimoso Alluvial Aquifer (<b>D</b>); soil surface texture map of Mimoso Alluvial Aquifer (<b>E</b>); well (<b>F</b>) and piezometers (<b>G</b>) monitoring; and piezometers installation (<b>H</b>,<b>I</b>).</p>
Full article ">Figure 2
<p>Soil map Location of Mimoso Alluvial Valley and locations of irrigated plots (adapted from [<a href="#B31-water-17-00410" class="html-bibr">31</a>]).</p>
Full article ">Figure 3
<p>Historical monthly averages and 2019 data for meteorological parameters as rainfall (<b>A</b>), water exploitation averages (<b>B</b>), and the maximum, minimum, and mean EC values for 2019 (<b>C</b>). (R = accumulated rainfall; Rainfall = historical average rainfall of 20 years; Et0 = reference evapotranspiration; Et0* = historical average evapotranspiration; Qmean = daily exploitation per month; ECmean = mean electrical conductivity; ECmin = minimum electrical conductivity; (<b>D</b>) ECmax = maximum electrical conductivity).</p>
Full article ">Figure 4
<p>Semivariograms scaled by electrical conductivity variances and adjustments of the exponential, gaussian and spherical models, from January to December (<b>A</b>–<b>L</b>) 2019. C0 represents the nugget effect, C1 the sill, A the model range and R<sup>2</sup> the determination coefficient.</p>
Full article ">Figure 5
<p>Semivariograms scaled by the SAR variances and adjustments of the exponential, gaussian and spherical models, from May to November (<b>A</b>–<b>G</b>) 2019. C0 represents the nugget effect, C1 the sill, A the model range and R<sup>2</sup> the determination coefficient.</p>
Full article ">Figure 6
<p>Crossed semivariograms of the SAR × CE and adjustments of the exponential, gaussian and spherical models, from May to November (<b>A</b>–<b>G</b>) 2019. C0 represents the nugget effect, C1 the sill, A the model range and R<sup>2</sup> the determination coefficient.</p>
Full article ">Figure 7
<p>Electrical conductivity (EC) kriging maps from January to December (<b>A</b>–<b>L</b>) 2019 with agricultural plots and saturated electrical conductivity contours.</p>
Full article ">Figure 8
<p>Sodium adsorption ratio (SAR) kriging maps from May to November (<b>A</b>–<b>G</b>) from 2019, with agricultural plots and saturated electrical conductivity contours.</p>
Full article ">Figure 9
<p>Semivariograms of the groundwater adequacy indicators and adjustments of the exponential, gaussian and spherical models, from May to November (<b>A</b>–<b>G</b>) 2019.</p>
Full article ">Figure 10
<p>Kriging maps of groundwater quality for irrigation purposes from May to November (<b>A</b>–<b>G</b>) 2019, with agricultural plots and saturated electrical conductivity contours.</p>
Full article ">Figure 11
<p>Sequential Gaussian simulation maps of groundwater quality for irrigation purposes from May to November (<b>A</b>–<b>G</b>) 2019 with agricultural plots and saturated electrical conductivity contours.</p>
Full article ">Figure 12
<p>Principal Component Analysis (PCA) showing the scores of the first two principal components (<b>A</b>,<b>B</b>), which represent the directions of maximum variance in the dataset.</p>
Full article ">
22 pages, 12094 KiB  
Article
Identification and Analysis on Surface Deformation in the Urban Area of Nanchang Based on PS-InSAR Method
by Mengping Zhang, Jiayi Pan, Peifeng Ma and Hui Lin
Remote Sens. 2025, 17(1), 157; https://doi.org/10.3390/rs17010157 - 5 Jan 2025
Viewed by 722
Abstract
Interferometric Synthetic Aperture Radar (InSAR) technology has emerged as a vital tool for monitoring surface deformation due to its high accuracy and spatial resolution. With the rapid economic development of Nanchang, extensive infrastructure development and construction activities have significantly altered the urban landscape. [...] Read more.
Interferometric Synthetic Aperture Radar (InSAR) technology has emerged as a vital tool for monitoring surface deformation due to its high accuracy and spatial resolution. With the rapid economic development of Nanchang, extensive infrastructure development and construction activities have significantly altered the urban landscape. Underground excavation and groundwater extraction in the region are potential contributors to surface deformation. This study utilized Sentinel-1 satellite data, acquired between September 2018 and May 2023, and applied the Permanent Scatterer Interferometric Synthetic Aperture Radar (PS-InSAR) technique to monitor surface deformation in Nanchang’s urban area. The findings revealed that surface deformation rates in the study area range from −10 mm/a to 6 mm/a, with the majority of regions remaining relatively stable. Approximately 99.9% of the monitored points exhibited deformation rates within −5 mm/a to 5 mm/a. However, four significant subsidence zones were identified along the Gan River and its downstream regions, with a maximum subsidence rate reaching 9.7 mm/a. Historical satellite imagery comparisons indicated that certain subsidence areas are potentially associated with construction activities. Further analysis integrating subsidence data, monthly precipitation, and groundwater depth revealed a negative correlation between surface deformation in Region A and rainfall, with subsidence trends aligning with groundwater level fluctuations. However, such a correlation was not evident in the other three regions. Additionally, water level data from the Xingzi Station of Poyang Lake showed that only Region A’s subsidence trend closely corresponds with water level variations. We conducted a detailed analysis of the spatial distribution of soil types in Nanchang and found that the soil types in areas of surface deformation are primarily Semi-hydromorphic Soils and Anthropogenic Soils. These soils exhibit high compressibility, making them prone to compaction and significantly influencing surface deformation. This study concludes that localized surface deformation in Nanchang is primarily driven by urban construction activities and the compaction of artificial fill soils, while precipitation also has an impact in certain areas. Full article
(This article belongs to the Section Environmental Remote Sensing)
Show Figures

Figure 1

Figure 1
<p>(<b>a</b>) Map of China, highlighting Nanchang’s location. (<b>b</b>) Map of Jiangxi Province, indicating where Nanchang is situated. (<b>c</b>) Map of Nanchang, showing its geographic features.</p>
Full article ">Figure 2
<p>PS-InSAR technical workflow.</p>
Full article ">Figure 3
<p>Temporal and spatial baseline diagram, with the central image being the master image and the others being the slave images.</p>
Full article ">Figure 4
<p>Surface deformation rate map of Nanchang City along the satellite line of sight from 2018 to 2023. Area A is located in Zhongxu Village, Nanchang County; Area B is situated along the shoreline of Xiazhuang Lake in Xinjian District; Area C is located Along Jiangzhong Avenue in Xihu District and at Sunshine Lighting Plaza; Area D is near the Shiqi Resettlement Housing in Nanchang County.</p>
Full article ">Figure 5
<p>Surface deformation rate distribution.</p>
Full article ">Figure 6
<p>Left (<b>a</b>) is the distribution map of Area A, and right (<b>b</b>) is the time-series subsidence line chart of three points in Area A.</p>
Full article ">Figure 7
<p>Left (<b>a</b>) is the distribution map of Area B, and right (<b>b</b>) is the time-series subsidence line chart of three points in Area B.</p>
Full article ">Figure 8
<p>(<b>a</b>–<b>c</b>) show the Google Earth images of Region B. (<b>a</b>) represents 2 May 2014, (<b>b</b>) represents 14 February 2017, and (<b>c</b>) represents 16 November 2019. Label 1 indicates the location of Baojie Machinery Company and Aonong Central China Science and Technology Park, Label 2 represents the location of Huihua Industrial Company.</p>
Full article ">Figure 9
<p>Left (<b>a</b>) is the distribution map of Area C, and right (<b>b</b>) is the time-series subsidence line chart of three points in Area C.</p>
Full article ">Figure 10
<p>(<b>a</b>–<b>c</b>) show the Google Earth images of Region C. (<b>a</b>) represents 16 November 2019, (<b>b</b>) represents 15 November 2020, and (<b>c</b>) represents 3 March 2022. The red polygon indicates the location of the Oupengwan project.</p>
Full article ">Figure 11
<p>Left (<b>a</b>) is the distribution map of Area D, and right (<b>b</b>) is the time-series subsidence line chart of three points in Area D.</p>
Full article ">Figure 12
<p>(<b>a</b>–<b>c</b>) show the Google Earth images of Region D. (<b>a</b>) represents 2 May 2014, (<b>b</b>) represents 27 March 2017, and (<b>c</b>) represents 15 November 2020. Labels 1, 2, and 3 mark the areas of subsidence corresponding to the three points in Region D.</p>
Full article ">Figure 13
<p>Spatial distribution map of soil types in the Nanchang area.</p>
Full article ">Figure 14
<p>(<b>a</b>–<b>d</b>) Relationship between subsidence at various points and monthly cumulative precipitation. (<b>a</b>) subsidence-precipitation relationship in Region A, (<b>b</b>) subsidence-precipitation relationship in Region B, (<b>c</b>) subsidence-precipitation relationship in Region C, (<b>d</b>) subsidence-precipitation relationship in Region D.</p>
Full article ">Figure 15
<p>(<b>a</b>–<b>d</b>) Relationship between precipitation in four regions and the average depth to groundwater in the Poyang Lake Plain. (<b>a</b>) precipitation-groundwater depth relationship in Region A, (<b>b</b>) precipitation-groundwater depth relationship in Region B, (<b>c</b>) precipitation-groundwater depth relationship in Region C, (<b>d</b>) precipitation-groundwater depth relationship in Region D.</p>
Full article ">Figure 16
<p>(<b>a</b>–<b>d</b>) Relationship between water level of Xingzi Station and the subsidence at various points. (<b>a</b>) subsidence-water level relationship in Region A, (<b>b</b>) subsidence- water level relationship in Region B, (<b>c</b>) subsidence- water level relationship in Region C, (<b>d</b>) subsidence- water level relationship in Region D.</p>
Full article ">
37 pages, 5371 KiB  
Article
Coupling Advanced Geo-Environmental Indices for the Evaluation of Groundwater Quality: A Case Study in NE Peloponnese, Greece
by Panagiotis Papazotos, Maria Vlachomitrou, Despoina Psarraki, Eleni Vasileiou and Maria Perraki
Environments 2025, 12(1), 14; https://doi.org/10.3390/environments12010014 - 4 Jan 2025
Viewed by 1294
Abstract
Water and its management have played a pivotal role in the evolution of organisms and civilizations, fulfilling essential roles in personal use, industry, irrigation, and drinking from ancient times to the present. This study seeks to evaluate groundwater quality for irrigation and drinking [...] Read more.
Water and its management have played a pivotal role in the evolution of organisms and civilizations, fulfilling essential roles in personal use, industry, irrigation, and drinking from ancient times to the present. This study seeks to evaluate groundwater quality for irrigation and drinking in the Northern Peloponnese region, specifically the wells of Loutraki and Schinos areas and the springs of the Gerania Mountains (Mts.), using geo-environmental indices and ionic ratios. For the first time, geo-environmental indices have been applied to a region where groundwater serves multiple purposes, addressing the challenge of understanding their dynamics to optimize their application in environmental science and groundwater pollution research. To achieve this, 68 groundwater samples from the study area were utilized, and a total of 25 geo-environmental indices were calculated to assess water quality. These indices examined: (i) drinking suitability (NPI, RI, PIG, WQI, and WPI), (ii) irrigation suitability (SAR, KR, %Na, PS, MAR, RSC, SSP, TH, PI, IWQI, and TDS), (iii) potentially toxic element (PTE) loadings (Cd, HEI, and HPI), and (iv) major hydrogeochemical processes, expressed as ionic ratios (Ca/Mg, Ca/SO4, Ca/Na, Cl/NO3, Cl/HCO3, and Si/NO3). Data processing involved descriptive statistics, hydrogeochemical bivariate plots, Spearman correlation coefficients, and multivariate statistical analyses, including factor analysis (FA) and R-mode hierarchical cluster analysis (HCA). Results revealed that all groundwater samples (100%) from the Loutraki area and the Gerania Mts. were of good quality for both drinking and irrigation purposes. In contrast, groundwater from the Schinos area exhibited lower quality, with most samples (93.9%) considered suitable only for irrigation. The deterioration in the coastal aquifer of the Schinos area is attributed to elevated concentrations of Cl, Na+, NO3, As, and Cr resulting from salinization and relatively limited anthropogenic influences. The study highlights that relying on individual geo-environmental indices can yield misleading results due to their dependence on factors such as researcher expertise, methodological choices, and the indices’ inherent limitations. Consequently, this research emphasizes the necessity of combining indices to enhance the reliability, accuracy, and robustness of groundwater quality assessments and hydrogeochemical evaluations. Last but not least, the findings demonstrate that calculating all available geo-environmental indices is unnecessary. Instead, selecting a subset of indices that either reflect the impact of specific elemental concentrations or can be effectively integrated with others is sufficient. This streamlined approach addresses challenges in optimizing geo-environmental index applications and contributes to improved groundwater resource management. Full article
(This article belongs to the Special Issue Research Progress in Groundwater Contamination and Treatment)
Show Figures

Graphical abstract

Graphical abstract
Full article ">Figure 1
<p>A simplified geological map of the study area [<a href="#B34-environments-12-00014" class="html-bibr">34</a>,<a href="#B35-environments-12-00014" class="html-bibr">35</a>] with the groundwater sampling sites. An enlarged image of the Schinos area is given in (A) [<a href="#B30-environments-12-00014" class="html-bibr">30</a>].</p>
Full article ">Figure 2
<p>Classification of the 25 calculated geo-environmental indices of this study.</p>
Full article ">Figure 3
<p>Cross-plots of (<b>a</b>) Na<sup>+</sup> vs. Cl<sup>−</sup>, (<b>b</b>) SO<sub>4</sub><sup>2−</sup> vs. Ca<sup>2+</sup>, (<b>c</b>) Na<sup>+</sup> vs. (Ca<sup>2+</sup> + Mg<sup>2+</sup>), (<b>d</b>) EC vs. Na<sup>+</sup>/Cl<sup>−</sup>, (<b>e</b>) (HCO<sub>3<sup>−</sup></sub> + SO<sub>4</sub><sup>2−</sup>) vs. (Ca<sup>2+</sup> + Mg<sup>2+</sup>), and (<b>f</b>) HCO<sub>3<sup>−</sup></sub> vs. Mg<sup>2+</sup> for the 68 groundwater samples from the Loutraki–Schinos–Gerania Mts. region.</p>
Full article ">Figure 4
<p>Ficklin diagram [<a href="#B95-environments-12-00014" class="html-bibr">95</a>] of 68 groundwater samples showing the sum of PTEs vs. pH.</p>
Full article ">Figure 5
<p>Boxplots of ionic ratios (<b>a</b>) Ca/Mg, (<b>b</b>) Ca/SO<sub>4</sub>, (<b>c</b>) Ca/Na, (<b>d</b>) Cl/NO<sub>3</sub>, (<b>e</b>) Cl/HCO<sub>3</sub>, and (<b>f</b>) Si/NO<sub>3</sub> for the groundwater samples from the Loutraki–Schinos–Gerania Mts. region.</p>
Full article ">Figure 6
<p>The Spearman correlation matrix, along with significance levels (<span class="html-italic">p</span>-values), for the 25 calculated geo-environmental indices in the Loutraki–Schinos–Gerania Mts. region (n = 68 groundwater samples).</p>
Full article ">Figure 7
<p>The scree plot of eigenvalues for the components derived from 68 groundwater samples collected from the Loutraki–Schinos–Gerania Mts. region indicates that five components have eigenvalues &gt; 1, signifying their statistical significance within the FA approach.</p>
Full article ">Figure 8
<p>Dendrogram of R-mode HCA for 25 variables calculated in 68 groundwater samples from the Loutraki–Schinos–Gerania Mts. region. The red and yellow dashed lines represent the linkage distances used to create different distinct clusters.</p>
Full article ">Figure 9
<p>Cross-plot of WQI vs. HPI for the 68 groundwater samples from the Loutraki–Schinos–Gerania Mts. Region.</p>
Full article ">Figure 10
<p>Bivariate diagrams of (<b>a</b>) Cl/HCO<sub>3</sub> (molar ratio) vs. Cl (mmol/L) and Cl/HCO<sub>3</sub> (molar ratio) vs. HCO<sub>3</sub> (mmol/L) and (<b>b</b>) Si/NO<sub>3</sub> (molar ratio) vs. NO<sub>3</sub> (mmol/L) and Si/NO<sub>3</sub> (molar ratio) vs. Si (mmol/L) for the 68 groundwater samples from the Loutraki–Schinos–Gerania Mts. region.</p>
Full article ">Figure 11
<p>Evaluation of major hydrogeochemical processes affecting water chemistry using Cl/HCO<sub>3</sub> vs. Si/NO<sub>3</sub> diagram.</p>
Full article ">Figure 12
<p>Bivariate Cl/HCO<sub>3</sub> vs. Si/NO<sub>3</sub> diagrams of (<b>a</b>) As, and (<b>b</b>) Cr (Q1: first quartile; Q2: second quartile or median; Q3: third quartile).</p>
Full article ">
19 pages, 11710 KiB  
Article
Investigating the Structural Health of High-Rise Buildings and Its Influencing Factors Using Sentinel-1 Synthetic Aperture Radar Imagery: A Case Study of the Guangzhou–Foshan Metropolitan Area
by Di Huang, Zhixin Qi, Suya Lin, Yuze Gu, Wenxuan Song and Qianwen Lv
Buildings 2024, 14(12), 4074; https://doi.org/10.3390/buildings14124074 - 21 Dec 2024
Viewed by 1031
Abstract
Urban growth is increasingly shifting from horizontal expansion to vertical development, resulting in skylines dominated by high-rise buildings. The post-construction operations and maintenance of these buildings are critical, requiring regular structural health monitoring (SHM) to proactively identify and address potential safety concerns. Interferometric [...] Read more.
Urban growth is increasingly shifting from horizontal expansion to vertical development, resulting in skylines dominated by high-rise buildings. The post-construction operations and maintenance of these buildings are critical, requiring regular structural health monitoring (SHM) to proactively identify and address potential safety concerns. Interferometric synthetic aperture radar (InSAR) has proven effective for monitoring building safety, but most studies rely on high-resolution synthetic aperture radar (SAR) images. The high cost and limited coverage of these images restrict their use for large-scale monitoring. Sentinel-1 medium-resolution SAR images, which are freely available and offer broad coverage, make large-scale SHM more feasible. However, studies on the use of Sentinel-1 SAR images for structural health monitoring, especially at large spatial scales, remain limited. To address this gap, in this study, Sentinel-1 SAR images and PS-InSAR technology are proposed for performing a comprehensive structural safety assessment of super high-rise buildings in the Guangzhou–Foshan Metropolitan Area (GFMA) and for analyzing the influencing factors. Our assessment shows that while the overall structural safety of these buildings is satisfactory, certain areas, including Pearl River New Town, central Huadu district in Guangzhou, and southeastern Shunde district in Foshan, exhibit suboptimal safety conditions. We verified these findings using GNSS data and on-site investigations, confirming that Sentinel-1 SAR imagery offers reliable accuracy for monitoring building structural health. Furthermore, we identified factors such as settlement in soft soil layers, the construction of surrounding (underground) infrastructure, and building aging, which could potentially impact building structural safety. The results demonstrate that Sentinel-1 SAR images provide a reliable, rapid, and cost-effective method for the large-scale monitoring of building stability, enhancing our understanding of the underlying mechanisms and informing strategies to prevent potential safety crises, and also ensuring the sustainable development of society. Full article
Show Figures

Figure 1

Figure 1
<p>(<b>a</b>,<b>b</b>) location of the study area; (<b>c</b>) three-dimensional view of buildings in part of the study area based on 3D-GloBFP.</p>
Full article ">Figure 2
<p>(<b>a</b>) Spatial and temporal baseline of Sentinel-1 images from March 2017 to March 2022; (<b>b</b>) diagram of the SHM parameters used in this study.</p>
Full article ">Figure 3
<p>(<b>a</b>) Distribution and (<b>b</b>) height statistics of super high-rise buildings; (<b>c</b>) the GDP per km<sup>2</sup> in 2020. For a clearer display, the area in (<b>a</b>) is the minimum bounding rectangle of the location of super high-rise buildings, rather than the GFMA.</p>
Full article ">Figure 4
<p>(<b>a</b>) Deformation velocity in the GFMA from 2017 to 2022; (<b>b</b>) comparison of the PS-InSAR and GNSS deformation time series; (<b>c</b>) the kernel density of safety level, S.L.K.D is short for safety level kernel density for clearer display.</p>
Full article ">Figure 5
<p>(<b>a</b>) Location and safety level of building B1 and B2; (<b>b</b>) cumulative deformation time series of PS points on building B2.</p>
Full article ">Figure 6
<p>(<b>a</b>) Location and (<b>b</b>) three-dimensional view of building B3 and B4; (<b>c</b>–<b>e</b>) cracks in the walls of B3 and surrounding buildings; (<b>f</b>,<b>g</b>) gaps in the corner of B4.</p>
Full article ">Figure 7
<p>(<b>a</b>) Location of building B5 and B6; (<b>b</b>) cumulative deformation time series of PS points on building B6, the variation amplitude increased in the light blue area after the blue dashed line; (<b>c</b>–<b>e</b>) cracks in building B5 and B6; (<b>f</b>) Huiyangyuan community photo.</p>
Full article ">
34 pages, 90974 KiB  
Article
Multi-Decadal Land Subsidence Risk Assessment at Major Italian Cities by Integrating PSInSAR with Urban Vulnerability
by Michelle Lenardón Sánchez, Celina Anael Farías and Francesca Cigna
Land 2024, 13(12), 2103; https://doi.org/10.3390/land13122103 - 5 Dec 2024
Viewed by 917
Abstract
This study assesses subsidence-induced risk to urban infrastructure in three major Italian cities—Rome, Bologna, and Florence—by integrating satellite-based persistent scatterer interferometric synthetic aperture radar (PSInSAR) ground displacement data with urban vulnerability metrics into a novel risk assessment workflow, incorporating land use and population [...] Read more.
This study assesses subsidence-induced risk to urban infrastructure in three major Italian cities—Rome, Bologna, and Florence—by integrating satellite-based persistent scatterer interferometric synthetic aperture radar (PSInSAR) ground displacement data with urban vulnerability metrics into a novel risk assessment workflow, incorporating land use and population data from the Copernicus Land Monitoring Service (CLMS)—Urban Atlas. This analysis exploits ERS-1/2, ENVISAT, and COSMO-SkyMed PSInSAR datasets from the Italian Extraordinary Plan of Environmental Remote Sensing, plus Sentinel-1 datasets from CLMS—European Ground Motion Service (EGMS), and spans a 30-year period, thus capturing both historical and recent subsidence trends. Angular distortion is introduced as a critical parameter for assessing potential structural damage due to differential settlement, which helps to quantify subsidence-induced hazards more precisely. The results reveal variable subsidence hazard patterns across the three cities, with specific areas exhibiting significant differential ground deformation that poses risks to key infrastructure. A total of 36.15, 11.44, and 0.43 km2 of land at high to very high risk are identified in Rome, Bologna, and Florence, respectively. By integrating geospatial and vulnerability data at the building-block level, this study offers a more comprehensive understanding of subsidence-induced risk, potentially contributing to improved management and mitigation strategies in urban areas. This study contributes to the limited literature on embedding PSInSAR data into urban risk assessment workflows and provides a replicable framework for future applications in other urban areas. Full article
(This article belongs to the Special Issue Assessing Land Subsidence Using Remote Sensing Data)
Show Figures

Figure 1

Figure 1
<p>Location and land use/cover of the cities of (<b>a</b>) Rome; (<b>b</b>) Bologna, and (<b>c</b>) Florence in Italy, according to the Copernicus Urban Atlas (UA) 2018 dataset [<a href="#B24-land-13-02103" class="html-bibr">24</a>]. Land use/cover types associated with UA codes are provided in <a href="#land-13-02103-t001" class="html-table">Table 1</a>.</p>
Full article ">Figure 2
<p>Comparison of the vertical displacement velocity derived using the Sentinel-1 datasets, for (<b>a</b>) Rome, (<b>b</b>) Bologna, and (<b>c</b>) Florence. The linear regression between ascending and descending mode geometries is represented with a dashed yellow line. Notation: ‘Asc.’, ascending; ‘Desc.’, descending; ‘Sum’, combined ascending and descending.</p>
Full article ">Figure 3
<p>Vertical displacement velocity (<span class="html-italic">V<sub>U</sub></span>) maps of the 2018–2022 period for the city of Rome (<b>a</b>), Bologna (<b>b</b>), and Florence (<b>c</b>). PSInSAR data are overlapped onto Google Earth imagery.</p>
Full article ">Figure 3 Cont.
<p>Vertical displacement velocity (<span class="html-italic">V<sub>U</sub></span>) maps of the 2018–2022 period for the city of Rome (<b>a</b>), Bologna (<b>b</b>), and Florence (<b>c</b>). PSInSAR data are overlapped onto Google Earth imagery.</p>
Full article ">Figure 4
<p>Vertical displacement (<span class="html-italic">d<sub>U</sub></span>) and yearly vertical displacement velocities (<span class="html-italic">V<sub>U</sub></span>) time series corresponding to the sample locations selected for (<b>a</b>) Rome, (<b>b</b>) Bologna, and (<b>c</b>) Florence.</p>
Full article ">Figure 5
<p>Rome hazard maps, representing (<b>a</b>) the 30-year period angular distortion, <span class="html-italic">β</span><sub><span class="html-italic">T</span></sub>, and (<b>b</b>) the 2018–2022 period angular distortion, <span class="html-italic">β</span><sub>2018–2022</sub>. The 2018 Urban Atlas layer is superimposed and zooms (I) and (II) are overlapped onto the OpenStreetMap on grey scale.</p>
Full article ">Figure 6
<p>Bologna hazard maps, representing (<b>a</b>) the 30-year period angular distortion, <span class="html-italic">β</span><sub><span class="html-italic">T</span></sub>, and (<b>b</b>) the 2018–2022 period angular distortion, <span class="html-italic">β</span><sub>2018–2022</sub>. The 2018 Urban Atlas layer is superimposed and zooms (I) and (II) are overlapped onto the OpenStreetMap on grey scale.</p>
Full article ">Figure 7
<p>Florence hazard maps, representing (<b>a</b>) the 30-year period angular distortion, <span class="html-italic">β</span><sub><span class="html-italic">T</span></sub>, and (<b>b</b>) the 2018–2022 period angular distortion, <span class="html-italic">β</span><sub>2018–2022</sub>. The 2018 Urban Atlas layer is superimposed and zooms (I) and (II) are overlapped onto the OpenStreetMap on grey scale.</p>
Full article ">Figure 8
<p>Exposure–vulnerability of urban infrastructure in (<b>a</b>) Rome, (<b>b</b>) Bologna, and (<b>c</b>) Florence, based on scoring of UA 2018 classes.</p>
Full article ">Figure 9
<p>Subsidence-induced risk in Rome, assessed by exploiting the risk matrix combining angular distortion and exposure–vulnerability scores: risk mapping referred to (<b>a</b>) the 30-year period 1992–2022, and (<b>b</b>) the 2018–2022 period. Risk maps are overlapped onto Google Earth imagery.</p>
Full article ">Figure 10
<p>Subsidence-induced risk in Bologna, assessed by exploiting the risk matrix combining angular distortion and exposure–vulnerability scores: risk mapping referred to (<b>a</b>) the 30-year period 1992–2022, and (<b>b</b>) the 2018–2022 period. Risk maps are overlapped onto Google Earth imagery.</p>
Full article ">Figure 11
<p>Subsidence-induced risk in Florence, assessed by exploiting the risk matrix combining angular distortion and exposure–vulnerability scores: risk mapping referred to (<b>a</b>) the 30-year period 1992–2022, and (<b>b</b>) the 2018–2022 period. Risk maps are overlapped onto Google Earth imagery.</p>
Full article ">Figure 12
<p>Overview of the subsidence-induced risk mapping statistics for Rome, Bologna, and Florence in 1992–2022, highlighting the amount of square kilometers of land and population involved in each risk category. Note that each bar represents 100% of the area of each city.</p>
Full article ">Figure 13
<p>Example of correlation between lithology, subsidence, and resulting risk in Rome. (<b>a</b>) Lithology map representing (1) alluvial fan deposits, (2) pozzolana, (3) marginal sandy facies, (4) anthropogenic deposits, and (5) river, obtained from Lazio Region open data catalog; (<b>b</b>) total angular distortion values with Urban Atlas polygons superimposed; and (<b>c</b>) risk map for the 1992–2022 period. Notation: VH, very high; H, high; M, medium; L, low; ND, no data; NA, not applicable.</p>
Full article ">Figure 14
<p>Example of correlation between groundwater extraction, subsidence, and derived risk in Bologna: (<b>a</b>) 2018–2022 angular distortion values with Urban Atlas polygons superimposed. The orange dot represents the position of ARPAE’S groundwater monitoring well (modified from [<a href="#B38-land-13-02103" class="html-bibr">38</a>]). The recorded change in its piezometric level is −11.45 m. (<b>b</b>) Exposure–vulnerability of urban infrastructure. (<b>c</b>) Risk map for the 2018–2022 period.</p>
Full article ">Figure 15
<p>Example of correlation between landslide processes and ground displacement measured in Florence. (<b>a</b>) Shows the landslides in the study area, obtained from the Italian Landslide Inventory (IFFI project), (<b>b</b>) total angular distortion map, and (<b>c</b>) the risk matrix map.</p>
Full article ">Figure A1
<p>Vertical displacement velocities of Rome in: (<b>a</b>) 1992–2000 (ERS-1/2 datasets), (<b>b</b>) 2002–2010 (ENVISAT datasets), and (<b>c</b>) 2013–2014 (COSMO-SkyMed dataset). PSInSAR data are overlapped onto Google Earth imagery.</p>
Full article ">Figure A1 Cont.
<p>Vertical displacement velocities of Rome in: (<b>a</b>) 1992–2000 (ERS-1/2 datasets), (<b>b</b>) 2002–2010 (ENVISAT datasets), and (<b>c</b>) 2013–2014 (COSMO-SkyMed dataset). PSInSAR data are overlapped onto Google Earth imagery.</p>
Full article ">Figure A2
<p>Vertical displacement velocities of Bologna in: (<b>a</b>) 1992–2000 (ERS-1/2 datasets), and (<b>b</b>) 2002–2010 (ENVISAT datasets). PSInSAR data are overlapped onto Google Earth imagery.</p>
Full article ">Figure A2 Cont.
<p>Vertical displacement velocities of Bologna in: (<b>a</b>) 1992–2000 (ERS-1/2 datasets), and (<b>b</b>) 2002–2010 (ENVISAT datasets). PSInSAR data are overlapped onto Google Earth imagery.</p>
Full article ">Figure A3
<p>Vertical displacement velocities of Florence in: (<b>a</b>) 1992–2000 (ERS-1/2 datasets), and (<b>b</b>) 2002–2010 (ENVISAT datasets). PSInSAR data are overlapped onto Google Earth imagery.</p>
Full article ">Figure A3 Cont.
<p>Vertical displacement velocities of Florence in: (<b>a</b>) 1992–2000 (ERS-1/2 datasets), and (<b>b</b>) 2002–2010 (ENVISAT datasets). PSInSAR data are overlapped onto Google Earth imagery.</p>
Full article ">Figure A4
<p>Rome hazard maps, representing (<b>a</b>) the 1992–2000 period angular distortion, <span class="html-italic">β</span><sub>1992–2000</sub>, (<b>b</b>) the 2003–2010 period angular distortion, <span class="html-italic">β</span><sub>2003–2010</sub>, and (<b>c</b>) the 2011–2014 period angular distortion, <span class="html-italic">β</span><sub>2011–2014</sub>. The 2018 Urban Atlas layer is superimposed and zooms (I) and (II) are overlapped onto the OpenStreetMap on grey scale.</p>
Full article ">Figure A4 Cont.
<p>Rome hazard maps, representing (<b>a</b>) the 1992–2000 period angular distortion, <span class="html-italic">β</span><sub>1992–2000</sub>, (<b>b</b>) the 2003–2010 period angular distortion, <span class="html-italic">β</span><sub>2003–2010</sub>, and (<b>c</b>) the 2011–2014 period angular distortion, <span class="html-italic">β</span><sub>2011–2014</sub>. The 2018 Urban Atlas layer is superimposed and zooms (I) and (II) are overlapped onto the OpenStreetMap on grey scale.</p>
Full article ">Figure A5
<p>Bologna hazard maps, representing (<b>a</b>) the 1992–2000 period angular distortion, <span class="html-italic">β</span><sub>1992–2000</sub> and (<b>b</b>) the 2003–2010 period angular distortion, <span class="html-italic">β</span><sub>2003–2010</sub>. The 2018 Urban Atlas layer is superimposed and zooms (I) and (II) are overlapped onto the OpenStreetMap on grey scale.</p>
Full article ">Figure A6
<p>Florence hazard maps, representing (<b>a</b>) the 1992–2000 period angular distortion, <span class="html-italic">β</span><sub>1992–2000</sub> and (<b>b</b>) the 2003–2010 period angular distortion, <span class="html-italic">β</span><sub>2003–2010</sub>. The 2018 Urban Atlas layer is superimposed and zooms (I) and (II) are overlapped onto the OpenStreetMap on grey scale.</p>
Full article ">Figure A6 Cont.
<p>Florence hazard maps, representing (<b>a</b>) the 1992–2000 period angular distortion, <span class="html-italic">β</span><sub>1992–2000</sub> and (<b>b</b>) the 2003–2010 period angular distortion, <span class="html-italic">β</span><sub>2003–2010</sub>. The 2018 Urban Atlas layer is superimposed and zooms (I) and (II) are overlapped onto the OpenStreetMap on grey scale.</p>
Full article ">
20 pages, 10429 KiB  
Article
Dynamic Geo-Visualization of Urban Land Subsidence and Land Cover Data Using PS-InSAR and Google Earth Engine (GEE) for Spatial Planning Assessment
by Joko Widodo, Edy Trihatmoko, Muhammad Rokhis Khomarudin, Mohammad Ardha, Udhi Catur Nugroho, Nugraheni Setyaningrum, Galih Prasetya Dinanta, Rahmat Arief, Andie Setiyoko, Dandy Aditya Novresiandi, Rendi Handika, Muhammad Priyatna, Shinichi Sobue, Dwi Sarah and Wawan Hermawan
Urban Sci. 2024, 8(4), 234; https://doi.org/10.3390/urbansci8040234 - 1 Dec 2024
Viewed by 1798
Abstract
The North Java coastal area, known as the Pantura region, is experiencing significant land subsidence, with certain areas sinking up to 10 cm per year. Pekalongan is among the most affected, with subsidence rates between 10 and 19 cm annually, mainly due to [...] Read more.
The North Java coastal area, known as the Pantura region, is experiencing significant land subsidence, with certain areas sinking up to 10 cm per year. Pekalongan is among the most affected, with subsidence rates between 10 and 19 cm annually, mainly due to groundwater extraction, sediment compaction, and coastal erosion. Other coastal cities, like Semarang and Demak, show rates averaging 4 to 10 cm per year. This rapid subsidence is due to favorable geological conditions and ongoing urban development. This study investigates land subsidence in Pekalongan using the PS-InSAR method and dynamic visualization of time-series land cover data. PS-InSAR was applied to 45 scenes from ALOS-2 PALSAR-2 to monitor subsidence from 2014 to 2022. The results were validated with in situ subsidence benchmarks. Urban development dynamics were analyzed through land cover and land use change (LULC) and population density over the same period, using the GLC_FCS30D dataset in the GEE to detect non-natural LULC. The PS-InSAR results indicated that over 60.9% of investigation points experienced subsidence, up to 100 cm between 2014 and 2022. Ground validation showed an 83% agreement with PS-InSAR results. A statistical analysis of LULC from 2014 to 2022 did not show significant built-up area development, but the extension of salt marshes and water bodies indicated subsidence expansion. The population density reached 6873 people per square km by 2022, causing extensive groundwater use for domestic and industrial purposes, further aggravating the subsidence. Full article
Show Figures

Figure 1

Figure 1
<p>An example of urban change in Pekalongan City is derived from LANDSAT data comparing 1993 (<b>a</b>) and 2019 (<b>b</b>).</p>
Full article ">Figure 2
<p>Research location (land area clipped from the ALOS-2 scene utilized in the study).</p>
Full article ">Figure 3
<p>PS-InSAR records distribution and its velocity (in cm) during 2014–2022.</p>
Full article ">Figure 4
<p>The velocity of the land subsidence within the research period using 17,121 records.</p>
Full article ">Figure 5
<p>A graph of the accuracy test results for ALOS-2 data and field observations shows an <span class="html-italic">NSE</span> value of 0.83.</p>
Full article ">Figure 6
<p>LULC within the research period (2014–2022) in the research location.</p>
Full article ">Figure 7
<p>The LULC comparison in 2014 and 2022.</p>
Full article ">Figure 8
<p>Pixel number changes for LULC in 2014 and 2022.</p>
Full article ">Figure 9
<p>Pekalongan city’s spatial plan.</p>
Full article ">
25 pages, 41258 KiB  
Article
The Deformation Monitoring Capability of Fucheng-1 Time-Series InSAR
by Zhouhang Wu, Wenjun Zhang, Jialun Cai, Hongyao Xiang, Jing Fan and Xiaomeng Wang
Sensors 2024, 24(23), 7604; https://doi.org/10.3390/s24237604 - 28 Nov 2024
Cited by 1 | Viewed by 833
Abstract
The Fucheng-1 (FC-1) satellite has successfully transitioned from its initial operational phase and is now undergoing a detailed performance assessment for time-series deformation monitoring. This study evaluates the surface deformation monitoring capabilities of the newly launched FC-1 satellite using the interferometric synthetic aperture [...] Read more.
The Fucheng-1 (FC-1) satellite has successfully transitioned from its initial operational phase and is now undergoing a detailed performance assessment for time-series deformation monitoring. This study evaluates the surface deformation monitoring capabilities of the newly launched FC-1 satellite using the interferometric synthetic aperture radar (InSAR) technique, particularly in urban applications. By analyzing the observation data from 20 FC-1 scenes and 20 Sentinel-1 scenes, deformation velocity maps of a university in Mianyang city were obtained using persistent scatterer interferometry (PSI) and distributed scatterer interferometry (DSI) techniques. The results show that thanks to the high resolution of 3 × 3 m of the FC-1 satellite, significantly more PS points and DS points were detected than those detected by Sentinel-1, by 13.4 times and 17.9 times, respectively. The distribution of the major deformation areas detected by both satellites in the velocity maps is generally consistent. FC-1 performs better than Sentinel-1 in monitoring densely structured and vegetation-covered areas. Its deformation monitoring capability at the millimeter level was further validated through comparison with leveling measurements, with average errors and root mean square errors of 1.761 mm and 2.172 mm, respectively. Its high-resolution and high-precision interferometry capabilities make it particularly promising in the commercial remote sensing market. Full article
(This article belongs to the Special Issue Recent Advances in Synthetic Aperture Radar (SAR) Remote Sensing)
Show Figures

Figure 1

Figure 1
<p>(<b>a</b>) Coverage areas of Sentinel-1 (purple) and FC-1 (brown), study area location marked by a five-pointed star, and COPDEM topographic map. (<b>b</b>) Google Maps image of the study area.</p>
Full article ">Figure 2
<p>Flow chart of DSI and PSI.</p>
Full article ">Figure 3
<p>(<b>a</b>) Spatio-temporal baseline map of FC-1 single master image. (<b>b</b>) Spatio-temporal baseline map of Sentinel-1 single master image.</p>
Full article ">Figure 4
<p>(<b>a</b>,<b>b</b>) Vertical deformation velocity maps from FC-1 using the DSI and PSI methods. (<b>c</b>,<b>d</b>) Vertical deformation velocity maps from Sentinel-1 using the DSI and PSI methods. (<b>e</b>) Drone orthophoto of the reference point.</p>
Full article ">Figure 5
<p>(<b>a</b>,<b>b</b>) Histograms of deformation velocity from FC-1 using the DSI and PSI methods. (<b>c</b>,<b>d</b>) Histograms of deformation velocity from Sentinel-1 using the DSI and PSI methods.</p>
Full article ">Figure 6
<p>(<b>a</b>) Schematic diagram of the research area on Google Earth. (<b>b</b>,<b>c</b>) Deformation rate maps of region R1 obtained by FC-1 and Sentinel-1 using the PSI method, with a drone image as the base map. (<b>d</b>–<b>g</b>) Deformation rate maps of regions R2 and R3 obtained by FC-1 and Sentinel-1 using the PSI method, with Google Earth as the base map. (<b>h</b>–<b>k</b>) Deformation rate maps of regions R4 and R5 obtained by FC-1 and Sentinel-1 using the DSI method, with Google Earth or a drone image as the base map.</p>
Full article ">Figure 7
<p>(<b>a</b>) Deformation velocity points obtained by FC-1 using the PSI method overlaid onto a drone image. (<b>b</b>) Deformation velocity points obtained by Sentinel-1 using the PSI method overlaid onto a drone image.</p>
Full article ">Figure 8
<p>(<b>a</b>,<b>d</b>) Deformation velocity maps from FC-1 and Sentinel-1 using the PSI method, with schematic maps of ZZ1 and ZZ2 locations. (<b>b</b>,<b>c</b>) PS deformation points from FC-1 overlaid onto drone oblique images of ZZ1 and ZZ2. (<b>e</b>,<b>f</b>) PS deformation points from Sentinel-1 overlaid onto drone oblique images of ZZ1 and ZZ2.</p>
Full article ">Figure 9
<p>(<b>a</b>,<b>b</b>) Deformation velocity points obtained by FC-1 using DSI and PSI methods overlaid onto Google imagery. (<b>c</b>,<b>d</b>) Deformation velocity points obtained by Sentinel-1 using DSI and PSI methods overlaid onto Google imagery.</p>
Full article ">Figure 10
<p>Diagram of road profile location.</p>
Full article ">Figure 11
<p>(<b>a</b>,<b>b</b>) Deformation velocity profile of FC-1 under the DSI and PSI methods. (<b>c</b>,<b>d</b>) Deformation velocity profile of Sentinel-1 under the DSI and PSI methods.</p>
Full article ">Figure 12
<p>(<b>a</b>,<b>b</b>) Deformation rate profiles of Sentinel-1 and FC-1 under the DSI method. (<b>c</b>) Diagram of position of vegetation section line. (<b>d</b>) UAV 3D model of vegetation area.</p>
Full article ">Figure 13
<p>(<b>a</b>,<b>b</b>) Coherence histograms and average coherence values for the PSI method with FC-1 and Sentinel-1. (<b>c</b>,<b>d</b>) Coherence histograms and average coherence values for the DSI method with FC-1 and Sentinel-1.</p>
Full article ">Figure 14
<p>(<b>a</b>,<b>b</b>) Standard deviation maps of deformation velocity for Sentinel-1 using PSI and DSI methods. (<b>c</b>,<b>d</b>) Standard deviation maps of deformation velocity for FC-1 using PSI and DSI methods.</p>
Full article ">Figure 15
<p>(<b>a</b>) Diagram of locations of four regions A, B, C and D. (<b>b</b>–<b>e</b>) Time-series settlement maps of FC-1 and Sentinel-1 using the DSI and PSI methods in four regions.</p>
Full article ">Figure 16
<p>(<b>a</b>–<b>d</b>) Spearman’s correlation matrix heatmaps of the time-series settlement amounts obtained by FC-1 and Sentinel-1 using the DSI and PSI methods in four regions.</p>
Full article ">Figure 17
<p>(<b>a</b>–<b>d</b>) Pearson’s correlation matrix plots of the time-series subsidence values between FC-1 and Sentinel-1 using the DSI and PSI methods in four regions.</p>
Full article ">Figure 18
<p>Illustrative Google Earth map showing the locations of level points.</p>
Full article ">Figure 19
<p>The subsidence measured by FC-1 using the DSI method compared to the subsidence measured by leveling.</p>
Full article ">Figure 20
<p>The subsidence measured by FC-1 and Sentinel-1 using the DSI method compared to the subsidence measured by leveling.</p>
Full article ">
33 pages, 4423 KiB  
Article
Harnessing Deep Learning for Real-Time Water Quality Assessment: A Sustainable Solution
by Selma Toumi, Sabrina Lekmine, Nabil Touzout, Hamza Moussa, Noureddine Elboughdiri, Reguia Boudraa, Ouided Benslama, Mohammed Kebir, Subhan Danish, Jie Zhang, Abdeltif Amrane and Hichem Tahraoui
Water 2024, 16(23), 3380; https://doi.org/10.3390/w16233380 - 24 Nov 2024
Cited by 3 | Viewed by 1402
Abstract
This study presents an innovative approach utilizing artificial intelligence (AI) for the prediction and classification of water quality parameters based on physico-chemical measurements. The primary objective was to enhance the accuracy, speed, and accessibility of water quality monitoring. Data collected from various water [...] Read more.
This study presents an innovative approach utilizing artificial intelligence (AI) for the prediction and classification of water quality parameters based on physico-chemical measurements. The primary objective was to enhance the accuracy, speed, and accessibility of water quality monitoring. Data collected from various water samples in Algeria were analyzed to determine key parameters such as conductivity, turbidity, pH, and total dissolved solids (TDS). These measurements were integrated into deep neural networks (DNNs) to predict indices such as the sodium adsorption ratio (SAR), magnesium hazard (MH), sodium percentage (SP), Kelley’s ratio (KR), potential salinity (PS), exchangeable sodium percentage (ESP), as well as Water Quality Index (WQI) and Irrigation Water Quality Index (IWQI). The DNNs model, optimized through the selection of various activation functions and hidden layers, demonstrated high precision, with a correlation coefficient (R) of 0.9994 and a low root mean square error (RMSE) of 0.0020. This AI-driven methodology significantly reduces the reliance on traditional laboratory analyses, offering real-time water quality assessments that are adaptable to local conditions and environmentally sustainable. This approach provides a practical solution for water resource managers, particularly in resource-limited regions, to efficiently monitor water quality and make informed decisions for public health and agricultural applications. Full article
Show Figures

Figure 1

Figure 1
<p>Mapping of sampling points in Médéa, Algeria.</p>
Full article ">Figure 2
<p>Detailed diagram illustrating the development and optimization process of Deep Neural Networks (DNNs).</p>
Full article ">Figure 3
<p>The architecture of the optimal DNNs model.</p>
Full article ">Figure 4
<p>Comparison between experimental and predicted values: (<b>a</b>) Training phase, (<b>b</b>) Validation phase, and (<b>c</b>) All phases.</p>
Full article ">Figure 4 Cont.
<p>Comparison between experimental and predicted values: (<b>a</b>) Training phase, (<b>b</b>) Validation phase, and (<b>c</b>) All phases.</p>
Full article ">Figure 5
<p>Comparison between experimental and predicted values for test data.</p>
Full article ">Figure 6
<p>Analysis of residuals using multiple techniques based on estimated values: (<b>a</b>) Comparison of experimental data with predicted values, and (<b>b</b>) Histogram of the frequency distribution of residuals.</p>
Full article ">Figure 6 Cont.
<p>Analysis of residuals using multiple techniques based on estimated values: (<b>a</b>) Comparison of experimental data with predicted values, and (<b>b</b>) Histogram of the frequency distribution of residuals.</p>
Full article ">Figure 7
<p>Application for prediction and classification of water quality.</p>
Full article ">
22 pages, 42906 KiB  
Article
Monitoring the Subsidence in Wan’an Town of Deyang Based on PS-InSAR Technology (Sichuan, China)
by Hongyi Guo, Antonio Miguel Martínez-Graña and José Angel González-Delgado
Sustainability 2024, 16(22), 10010; https://doi.org/10.3390/su162210010 - 16 Nov 2024
Viewed by 947
Abstract
In recent years, land subsidence has become a crucial factor affecting urban safety and sustainable development, especially in Wan’an Town. To accurately monitor and analyze the land subsidence in Wan’an Town, this study uses the PS-InSAR technique combined with an improved DEM for [...] Read more.
In recent years, land subsidence has become a crucial factor affecting urban safety and sustainable development, especially in Wan’an Town. To accurately monitor and analyze the land subsidence in Wan’an Town, this study uses the PS-InSAR technique combined with an improved DEM for detailed research on land subsidence in Wan’an Town. PS-InSAR, or Permanent Scatterer Interferometric SAR, is suitable for high-precision monitoring of surface deformation. The natural neighbor interpolation method optimizes DEM data, improving its spatial resolution and accuracy. In this study, multiple periods of SAR imagery data of Wan’an Town were collected and preprocessed through radiometric calibration, phase unwrapping, and other steps. Using the PS-InSAR technique, the phase information of permanent scatterers (PS points) on the surface was extracted to establish a deformation model and preliminarily analyze the land subsidence in Wan’an Town. Concurrently, the DEM data were optimized using the natural neighbor interpolation method to enhance its accuracy. Finally, the optimized DEM data were combined with the surface deformation information extracted through the PS-InSAR technique for a detailed analysis of the land subsidence in Wan’an Town. The research results indicate that the DEM data optimized by the natural neighbor interpolation method have higher accuracy and spatial resolution, providing a more accurate reflection of the topographical features of Wan’an Town. The research found that the optimized DEM provided a more accurate reflection of Wan’an Town’s topographical features. By combining PS-InSAR data, subsidence information from 2016 to 2024 was calculated. The study area showed varying degrees of subsidence, with rates ranging from 6 mm/year to 10 mm/year. Four characteristic deformation areas were analyzed for causes and influencing factors. The findings contribute to understanding urban land subsidence, guiding urban planning, and providing data support for geological disaster warning and prevention. Full article
Show Figures

Figure 1

Figure 1
<p>Digital elevation model of the study area.</p>
Full article ">Figure 2
<p>Topography of the study area and radar image coverage area.</p>
Full article ">Figure 3
<p>Geology map of the study area.</p>
Full article ">Figure 4
<p>Elevation contrast chart.</p>
Full article ">Figure 5
<p>Workflow of PS processing.</p>
Full article ">Figure 6
<p>Spatial and temporal baseline distribution map.</p>
Full article ">Figure 7
<p>Differential interferogram.</p>
Full article ">Figure 8
<p>Spatial distribution of the average subsidence rate in the study area from 2016 to 2024.</p>
Full article ">Figure 9
<p>Settlement comparison diagram.</p>
Full article ">Figure 10
<p>Total subsidence in the study area from 2016 to 2024.</p>
Full article ">Figure 11
<p>Natural neighbor interpolation.</p>
Full article ">Figure 12
<p>Time-series deformation map of the study area from 2016 to 2024.</p>
Full article ">Figure 13
<p>GPS survey map. (<b>A</b>) Field survey map of target A, (<b>B</b>) field survey map of target B, (<b>C</b>) field survey map of target C, (<b>D</b>) field survey map of target D.</p>
Full article ">Figure 14
<p>(<b>A</b>) Deformation analysis of Deformation Target Area A. (<b>B</b>) Deformation analysis of Deformation Target Area B. (<b>C</b>) Deformation analysis of Deformation Target Area C. (<b>D</b>) Deformation analysis of Deformation Target Area D.</p>
Full article ">Figure 14 Cont.
<p>(<b>A</b>) Deformation analysis of Deformation Target Area A. (<b>B</b>) Deformation analysis of Deformation Target Area B. (<b>C</b>) Deformation analysis of Deformation Target Area C. (<b>D</b>) Deformation analysis of Deformation Target Area D.</p>
Full article ">
28 pages, 6037 KiB  
Article
Statistical and Independent Component Analysis of Sentinel-1 InSAR Time Series to Assess Land Subsidence Trends
by Celina Anael Farías, Michelle Lenardón Sánchez, Roberta Bonì and Francesca Cigna
Remote Sens. 2024, 16(21), 4066; https://doi.org/10.3390/rs16214066 - 31 Oct 2024
Cited by 2 | Viewed by 1882
Abstract
Advanced statistics can enable the detailed characterization of ground deformation time series, which is a fundamental step for thoroughly understanding the phenomena of land subsidence and their main drivers. This study presents a novel methodological approach based on pre-existing open-access statistical tools to [...] Read more.
Advanced statistics can enable the detailed characterization of ground deformation time series, which is a fundamental step for thoroughly understanding the phenomena of land subsidence and their main drivers. This study presents a novel methodological approach based on pre-existing open-access statistical tools to exploit satellite differential interferometric synthetic aperture radar (DInSAR) data to investigate land subsidence processes, using European Ground Motion Service (EGMS) Sentinel-1 DInSAR 2018?2022 datasets. The workflow involves the implementation of Persistent Scatterers (PS) time series classification through the PS-Time tool, deformation signal decomposition via independent component analysis (ICA), and drivers’ investigation through spatio-temporal correlation with geospatial and monitoring data. Subsidence time series at the three demonstration sites of Bologna, Ravenna and Carpi (Po Plain, Italy) were classified into linear and nonlinear (quadratic, discontinuous, uncorrelated) categories, and the mixed deformation signal of each PS was decomposed into independent components, allowing the identification of new spatial clusters with linear, accelerating/decelerating, and seasonal trends. The relationship between the different independent components and DInSAR-derived displacement velocity, acceleration, and seasonality was also analyzed via regression analysis. Correlation with geological and groundwater monitoring data supported the investigation of the relationship between the observed deformation and subsidence drivers, such as aquifer resource exploitation, local geological setting, and gas extraction/reinjection. Full article
(This article belongs to the Special Issue Monitoring Geohazard from Synthetic Aperture Radar Interferometry)
Show Figures

Figure 1

Figure 1
<p>The study areas of (1) Ravenna, (2) Bologna, and (3) Carpi–Correggio–Soliera: (<b>a</b>) geographical location in Italy; (<b>b</b>) extent of the European Ground Motion Service (EGMS) Level-3 (L3) and Level-2b (L2b) dataset footprints used for the statistical analysis, overlapped onto the Copernicus Global Digital Elevation Model [<a href="#B44-remotesensing-16-04066" class="html-bibr">44</a>]; and (<b>c</b>) detail of the mean vertical deformation velocity from EGMS L3 datasets, overlapped onto a Google satellite imagery basemap.</p>
Full article ">Figure 2
<p>(<b>a</b>) Mean LOS deformation velocity; (<b>b</b>) acceleration; (<b>c</b>) annual seasonality amplitude; and (<b>d</b>) PS-Time classification maps for Ravenna, overlapped onto Google satellite imagery. The area selected for the following ICA analysis is highlighted on (<b>a</b>). DCV = discontinuous with constant velocity; DVV = discontinuous with variable velocity.</p>
Full article ">Figure 3
<p>(<b>a</b>) Mean LOS deformation velocity; (<b>b</b>) acceleration; (<b>c</b>) annual seasonality amplitude; and (<b>d</b>) PS-Time classification maps for Bologna, overlapped onto Google satellite imagery, with indication of the administrative boundary of the city of Bologna (black polygon). The rectangles (i.e., 1 in (<b>c</b>), and 2 in (<b>b</b>)) indicate the testing areas utilized in the following ICA analysis. DCV = discontinuous with constant velocity; DVV = discontinuous with variable velocity.</p>
Full article ">Figure 4
<p>(<b>a</b>) Mean LOS deformation velocity; (<b>b</b>) acceleration; (<b>c</b>) annual seasonality amplitude; and (<b>d</b>) PS-Time classification maps in the Carpi–Correggio–Soliera area, overlapped onto Google satellite imagery. The area selected for the following ICA analysis is highlighted on (<b>a</b>). DCV = discontinuous with constant velocity; DVV = discontinuous with variable velocity.</p>
Full article ">Figure 5
<p>Independent components identified in Ravenna (Ra) testing area, overlapped onto Google satellite imagery.</p>
Full article ">Figure 6
<p>Independent components identified in Bologna (Bo), covering (<b>a</b>) Area 1, and (<b>b</b>) Area 2, overlapped onto Google satellite imagery.</p>
Full article ">Figure 7
<p>Independent components identified in Soliera (So), overlapped onto Google satellite imagery.</p>
Full article ">Figure 8
<p>Correlation between mean deformation velocity, acceleration and amplitude of the APC, from EGMS products and PS-Time analysis, and the ICs retained for the area of Ravenna. Linear or bilinear fitting (red lines) and R<sup>2</sup> values are shown in the graphs that show the best correlation.</p>
Full article ">Figure 9
<p>Correlation between mean deformation velocity, acceleration, and amplitude of the APC, from EGMS products and PS-Time analysis, and the ICs retained for Area 1 in Bologna. Linear or quadratic fitting (red lines) and R2 values are shown in the graphs that show the best correlation.</p>
Full article ">Figure 10
<p>Correlation between mean deformation velocity, acceleration, and amplitude of the APC, from EGMS products and PS-Time analysis, and the ICs retained for Area 2 in Bologna. Linear or quadratic fitting (red lines) and R<sup>2</sup> values are shown in the graphs that show the best correlation.</p>
Full article ">Figure 11
<p>Correlation between mean deformation velocity, acceleration, and amplitude of the APC, from EGMS products and PS-Time analysis, and the ICs retained for the area of Soliera. Linear or quadratic fitting (red lines) and R<sup>2</sup> values are shown in the graphs that show the best correlation.</p>
Full article ">Figure 12
<p>(<b>a</b>) Geological map with the location of the principal gas fields operating near the coast of Ravenna and two ARPAE groundwater monitoring wells, overlapped onto Google satellite imagery; (<b>b</b>) Comparison between piezometric level variations in ARPAE’s monitoring wells RA49-00 and RA29-00 and the deformation time series of contiguous points; (<b>c</b>) Deformation velocities observed within each lithological unit, expressed in [mm/year]. Gas exploitation data in (<b>a</b>) is made available by the Italian Ministry of Environment and Energy Security [<a href="#B58-remotesensing-16-04066" class="html-bibr">58</a>], while the location of the monitoring wells and the geological layers were downloaded from the MinERva Portal, managed by the Emilia-Romagna Region service [<a href="#B46-remotesensing-16-04066" class="html-bibr">46</a>].</p>
Full article ">Figure 13
<p>(<b>a</b>) Position of ARPAE’S groundwater monitoring wells and the recorded change in piezometric levels (Δ<span class="html-italic">h<sub>i</sub></span>) for the area of Bologna during the studied time period (2018−2022), overlapped onto Google Satellite imagery; (<b>b</b>) Comparison between piezometric level variations in three of the wells and the deformation time series of contiguous PS–DS points; (<b>c</b>) Geological map of Bologna; (<b>d</b>) Deformation velocities observed within each lithological unit, ex-pressed in [mm/year]. Geological layers used in (<b>c</b>) were downloaded from the MinERva Portal, managed by the Emilia-Romagna Region service [<a href="#B46-remotesensing-16-04066" class="html-bibr">46</a>].</p>
Full article ">Figure 14
<p>(<b>a</b>) Geological map of Carpi–Correggio–Soliera subsidence hotspot, overlapped onto Google satellite imagery; (<b>b</b>) Comparison between piezometric level variations in MO10-01 ARPAE’s monitoring well and a deformation time series of a contiguous PS–DS point; (<b>c</b>) Deformation velocities observed within each lithology, expressed in [mm/year]. Geological layers used were downloaded from the MinERva Portal, managed by the Emilia-Romagna Region service [<a href="#B46-remotesensing-16-04066" class="html-bibr">46</a>].</p>
Full article ">Figure A1
<p>Example of a time series classified as “Bilinear” by PS-Time automatic classification algorithm, in the southern area of Soliera.</p>
Full article ">Figure A2
<p>Time series of one of the PS–DS points scored positively for Bo2–IC2 seasonal component.</p>
Full article ">Figure A3
<p>Acceleration variations vs. buffer distances from Angela Angelina reinjection well.</p>
Full article ">
Back to TopTop