{"id":1389,"date":"2026-10-23T08:18:59","date_gmt":"2026-10-23T08:18:59","guid":{"rendered":"https:\/\/helpwithmygisassignment.com\/blog\/?p=1389"},"modified":"2026-10-23T08:18:59","modified_gmt":"2026-10-23T08:18:59","slug":"insar-surface-deformation-mapping-sentinel-1-subsidence","status":"publish","type":"post","link":"https:\/\/helpwithmygisassignment.com\/blog\/insar-surface-deformation-mapping-sentinel-1-subsidence\/","title":{"rendered":"InSAR Surface Deformation Mapping: Tracking Ground Subsidence with Sentinel-1"},"content":{"rendered":"<div><\/div>\n<div>\n<div id=\"ez-toc-container\" class=\"ez-toc-v2_0_89 counter-hierarchy ez-toc-counter ez-toc-grey ez-toc-container-direction\">\n<div class=\"ez-toc-title-container\">\n<p class=\"ez-toc-title\" style=\"cursor:inherit\">Table of Contents<\/p>\n<span class=\"ez-toc-title-toggle\"><a href=\"#\" class=\"ez-toc-pull-right ez-toc-btn ez-toc-btn-xs ez-toc-btn-default ez-toc-toggle\" aria-label=\"Toggle Table of Content\"><span class=\"ez-toc-js-icon-con\"><span class=\"\"><span class=\"eztoc-hide\" style=\"display:none;\">Toggle<\/span><span class=\"ez-toc-icon-toggle-span\"><svg style=\"fill: #999;color:#999\" xmlns=\"http:\/\/www.w3.org\/2000\/svg\" class=\"list-377408\" width=\"20px\" height=\"20px\" viewBox=\"0 0 24 24\" fill=\"none\"><path d=\"M6 6H4v2h2V6zm14 0H8v2h12V6zM4 11h2v2H4v-2zm16 0H8v2h12v-2zM4 16h2v2H4v-2zm16 0H8v2h12v-2z\" fill=\"currentColor\"><\/path><\/svg><svg style=\"fill: #999;color:#999\" class=\"arrow-unsorted-368013\" xmlns=\"http:\/\/www.w3.org\/2000\/svg\" width=\"10px\" height=\"10px\" viewBox=\"0 0 24 24\" version=\"1.2\" baseProfile=\"tiny\"><path d=\"M18.2 9.3l-6.2-6.3-6.2 6.3c-.2.2-.3.4-.3.7s.1.5.3.7c.2.2.4.3.7.3h11c.3 0 .5-.1.7-.3.2-.2.3-.5.3-.7s-.1-.5-.3-.7zM5.8 14.7l6.2 6.3 6.2-6.3c.2-.2.3-.5.3-.7s-.1-.5-.3-.7c-.2-.2-.4-.3-.7-.3h-11c-.3 0-.5.1-.7.3-.2.2-.3.5-.3.7s.1.5.3.7z\"\/><\/svg><\/span><\/span><\/span><\/a><\/span><\/div>\n<nav><ul class='ez-toc-list ez-toc-list-level-1 ' ><li class='ez-toc-page-1 ez-toc-heading-level-1'><a class=\"ez-toc-link ez-toc-heading-1\" href=\"https:\/\/helpwithmygisassignment.com\/blog\/insar-surface-deformation-mapping-sentinel-1-subsidence\/#InSAR_Surface_Deformation_Mapping_Tracking_Ground_Subsidence_with_Sentinel-1\" >InSAR Surface Deformation Mapping: Tracking Ground Subsidence with Sentinel-1<\/a><ul class='ez-toc-list-level-2' ><li class='ez-toc-heading-level-2'><a class=\"ez-toc-link ez-toc-heading-2\" href=\"https:\/\/helpwithmygisassignment.com\/blog\/insar-surface-deformation-mapping-sentinel-1-subsidence\/#1_Introduction_How_InSAR_detects_ground_subsidence\" >1. Introduction: How InSAR detects ground subsidence<\/a><ul class='ez-toc-list-level-3' ><li class='ez-toc-heading-level-3'><a class=\"ez-toc-link ez-toc-heading-3\" href=\"https:\/\/helpwithmygisassignment.com\/blog\/insar-surface-deformation-mapping-sentinel-1-subsidence\/#11_Understanding_interferometric_phase\" >1.1 Understanding interferometric phase<\/a><\/li><li class='ez-toc-page-1 ez-toc-heading-level-3'><a class=\"ez-toc-link ez-toc-heading-4\" href=\"https:\/\/helpwithmygisassignment.com\/blog\/insar-surface-deformation-mapping-sentinel-1-subsidence\/#12_Converting_phase_into_displacement\" >1.2 Converting phase into displacement<\/a><\/li><li class='ez-toc-page-1 ez-toc-heading-level-3'><a class=\"ez-toc-link ez-toc-heading-5\" href=\"https:\/\/helpwithmygisassignment.com\/blog\/insar-surface-deformation-mapping-sentinel-1-subsidence\/#13_The_processing_workflow_at_a_glance\" >1.3 The processing workflow at a glance<\/a><\/li><\/ul><\/li><li class='ez-toc-page-1 ez-toc-heading-level-2'><a class=\"ez-toc-link ez-toc-heading-6\" href=\"https:\/\/helpwithmygisassignment.com\/blog\/insar-surface-deformation-mapping-sentinel-1-subsidence\/#2_Data_requirements_Sentinel-1_SLC_image_pairs\" >2. Data requirements: Sentinel-1 SLC image pairs<\/a><ul class='ez-toc-list-level-3' ><li class='ez-toc-heading-level-3'><a class=\"ez-toc-link ez-toc-heading-7\" href=\"https:\/\/helpwithmygisassignment.com\/blog\/insar-surface-deformation-mapping-sentinel-1-subsidence\/#21_Download_suitable_Sentinel-1_acquisitions\" >2.1 Download suitable Sentinel-1 acquisitions<\/a><\/li><li class='ez-toc-page-1 ez-toc-heading-level-3'><a class=\"ez-toc-link ez-toc-heading-8\" href=\"https:\/\/helpwithmygisassignment.com\/blog\/insar-surface-deformation-mapping-sentinel-1-subsidence\/#22_Why_SLC_is_essential\" >2.2 Why SLC is essential<\/a><\/li><li class='ez-toc-page-1 ez-toc-heading-level-3'><a class=\"ez-toc-link ez-toc-heading-9\" href=\"https:\/\/helpwithmygisassignment.com\/blog\/insar-surface-deformation-mapping-sentinel-1-subsidence\/#23_Prepare_the_supporting_datasets\" >2.3 Prepare the supporting datasets<\/a><\/li><\/ul><\/li><li class='ez-toc-page-1 ez-toc-heading-level-2'><a class=\"ez-toc-link ez-toc-heading-10\" href=\"https:\/\/helpwithmygisassignment.com\/blog\/insar-surface-deformation-mapping-sentinel-1-subsidence\/#3_Step-by-step_processing_workflow_in_SNAP\" >3. Step-by-step processing workflow in SNAP<\/a><ul class='ez-toc-list-level-3' ><li class='ez-toc-heading-level-3'><a class=\"ez-toc-link ez-toc-heading-11\" href=\"https:\/\/helpwithmygisassignment.com\/blog\/insar-surface-deformation-mapping-sentinel-1-subsidence\/#Step_1_Apply_orbit_corrections_and_select_matching_bursts\" >Step 1: Apply orbit corrections and select matching bursts<\/a><\/li><li class='ez-toc-page-1 ez-toc-heading-level-3'><a class=\"ez-toc-link ez-toc-heading-12\" href=\"https:\/\/helpwithmygisassignment.com\/blog\/insar-surface-deformation-mapping-sentinel-1-subsidence\/#Step_2_Coregister_the_SLC_images\" >Step 2: Coregister the SLC images<\/a><\/li><li class='ez-toc-page-1 ez-toc-heading-level-3'><a class=\"ez-toc-link ez-toc-heading-13\" href=\"https:\/\/helpwithmygisassignment.com\/blog\/insar-surface-deformation-mapping-sentinel-1-subsidence\/#Step_3_Generate_the_interferogram_and_coherence_raster\" >Step 3: Generate the interferogram and coherence raster<\/a><\/li><li class='ez-toc-page-1 ez-toc-heading-level-3'><a class=\"ez-toc-link ez-toc-heading-14\" href=\"https:\/\/helpwithmygisassignment.com\/blog\/insar-surface-deformation-mapping-sentinel-1-subsidence\/#Step_4_Remove_the_topographic_phase\" >Step 4: Remove the topographic phase<\/a><\/li><li class='ez-toc-page-1 ez-toc-heading-level-3'><a class=\"ez-toc-link ez-toc-heading-15\" href=\"https:\/\/helpwithmygisassignment.com\/blog\/insar-surface-deformation-mapping-sentinel-1-subsidence\/#Step_5_Filter_the_phase_and_prepare_for_unwrapping\" >Step 5: Filter the phase and prepare for unwrapping<\/a><\/li><li class='ez-toc-page-1 ez-toc-heading-level-3'><a class=\"ez-toc-link ez-toc-heading-16\" href=\"https:\/\/helpwithmygisassignment.com\/blog\/insar-surface-deformation-mapping-sentinel-1-subsidence\/#Step_6_Unwrap_the_interferometric_phase_with_SNAPHU\" >Step 6: Unwrap the interferometric phase with SNAPHU<\/a><\/li><li class='ez-toc-page-1 ez-toc-heading-level-3'><a class=\"ez-toc-link ez-toc-heading-17\" href=\"https:\/\/helpwithmygisassignment.com\/blog\/insar-surface-deformation-mapping-sentinel-1-subsidence\/#Step_7_Convert_unwrapped_phase_to_displacement\" >Step 7: Convert unwrapped phase to displacement<\/a><\/li><\/ul><\/li><li class='ez-toc-page-1 ez-toc-heading-level-2'><a class=\"ez-toc-link ez-toc-heading-18\" href=\"https:\/\/helpwithmygisassignment.com\/blog\/insar-surface-deformation-mapping-sentinel-1-subsidence\/#4_Track_land_subsidence_in_QGIS\" >4. Track land subsidence in QGIS<\/a><ul class='ez-toc-list-level-3' ><li class='ez-toc-heading-level-3'><a class=\"ez-toc-link ez-toc-heading-19\" href=\"https:\/\/helpwithmygisassignment.com\/blog\/insar-surface-deformation-mapping-sentinel-1-subsidence\/#Step_1_Load_the_rasters\" >Step 1: Load the rasters<\/a><\/li><li class='ez-toc-page-1 ez-toc-heading-level-3'><a class=\"ez-toc-link ez-toc-heading-20\" href=\"https:\/\/helpwithmygisassignment.com\/blog\/insar-surface-deformation-mapping-sentinel-1-subsidence\/#Step_2_Calculate_an_annualized_LOS_displacement_rate\" >Step 2: Calculate an annualized LOS displacement rate<\/a><\/li><li class='ez-toc-page-1 ez-toc-heading-level-3'><a class=\"ez-toc-link ez-toc-heading-21\" href=\"https:\/\/helpwithmygisassignment.com\/blog\/insar-surface-deformation-mapping-sentinel-1-subsidence\/#Step_3_Apply_a_coherence_mask\" >Step 3: Apply a coherence mask<\/a><\/li><li class='ez-toc-page-1 ez-toc-heading-level-3'><a class=\"ez-toc-link ez-toc-heading-22\" href=\"https:\/\/helpwithmygisassignment.com\/blog\/insar-surface-deformation-mapping-sentinel-1-subsidence\/#Step_4_Style_the_deformation_map\" >Step 4: Style the deformation map<\/a><\/li><\/ul><\/li><li class='ez-toc-page-1 ez-toc-heading-level-2'><a class=\"ez-toc-link ez-toc-heading-23\" href=\"https:\/\/helpwithmygisassignment.com\/blog\/insar-surface-deformation-mapping-sentinel-1-subsidence\/#5_Limitations_and_quality_control\" >5. Limitations and quality control<\/a><ul class='ez-toc-list-level-3' ><li class='ez-toc-heading-level-3'><a class=\"ez-toc-link ez-toc-heading-24\" href=\"https:\/\/helpwithmygisassignment.com\/blog\/insar-surface-deformation-mapping-sentinel-1-subsidence\/#Can_LOS_displacement_be_converted_to_vertical_subsidence\" >Can LOS displacement be converted to vertical subsidence?<\/a><\/li><\/ul><\/li><li class='ez-toc-page-1 ez-toc-heading-level-2'><a class=\"ez-toc-link ez-toc-heading-25\" href=\"https:\/\/helpwithmygisassignment.com\/blog\/insar-surface-deformation-mapping-sentinel-1-subsidence\/#6_Alternative_workflow_using_GMTSAR\" >6. Alternative workflow using GMTSAR<\/a><\/li><li class='ez-toc-page-1 ez-toc-heading-level-2'><a class=\"ez-toc-link ez-toc-heading-26\" href=\"https:\/\/helpwithmygisassignment.com\/blog\/insar-surface-deformation-mapping-sentinel-1-subsidence\/#7_Summary\" >7. Summary<\/a><\/li><\/ul><\/li><\/ul><\/nav><\/div>\n<h1><span class=\"ez-toc-section\" id=\"InSAR_Surface_Deformation_Mapping_Tracking_Ground_Subsidence_with_Sentinel-1\"><\/span>InSAR Surface Deformation Mapping: Tracking Ground Subsidence with Sentinel-1<span class=\"ez-toc-section-end\"><\/span><\/h1>\n<h2><span class=\"ez-toc-section\" id=\"1_Introduction_How_InSAR_detects_ground_subsidence\"><\/span>1. Introduction: How InSAR detects ground subsidence<span class=\"ez-toc-section-end\"><\/span><\/h2>\n<p><img loading=\"lazy\" decoding=\"async\" class=\"aligncenter size-large wp-image-1542\" src=\"https:\/\/helpwithmygisassignment.com\/blog\/wp-content\/uploads\/2026\/10\/Gemini_Generated_Image_zhj019zhj019zhj0-1024x572.jpg\" alt=\"InSAR Surface Deformation Mapping Sentinel 1 Subsidence\" width=\"1024\" height=\"572\" srcset=\"https:\/\/helpwithmygisassignment.com\/blog\/wp-content\/uploads\/2026\/10\/Gemini_Generated_Image_zhj019zhj019zhj0-1024x572.jpg 1024w, https:\/\/helpwithmygisassignment.com\/blog\/wp-content\/uploads\/2026\/10\/Gemini_Generated_Image_zhj019zhj019zhj0-300x167.jpg 300w, https:\/\/helpwithmygisassignment.com\/blog\/wp-content\/uploads\/2026\/10\/Gemini_Generated_Image_zhj019zhj019zhj0-768x429.jpg 768w, https:\/\/helpwithmygisassignment.com\/blog\/wp-content\/uploads\/2026\/10\/Gemini_Generated_Image_zhj019zhj019zhj0.jpg 1376w\" sizes=\"auto, (max-width: 1024px) 100vw, 1024px\" \/><\/p>\n<p>When performing InSAR surface deformation mapping Sentinel 1 subsidence, detecting millimeters of ground displacement from satellite radar requires careful processing. Interferometric Synthetic Aperture Radar (InSAR) is a remote sensing technique that measures changes in the distance between a radar satellite and the Earth&#8217;s surface by comparing the phase of radar signals acquired at different times.<\/p>\n<p>Unlike optical remote sensing, Synthetic Aperture Radar (SAR) actively transmits microwave pulses and records the returning signal. This makes Sentinel-1 useful for monitoring ground deformation in areas where clouds, darkness, or frequent optical image gaps would complicate conventional satellite mapping.<\/p>\n<div>\n<div><\/div>\n<\/div>\n<p>InSAR surface deformation mapping with Sentinel-1 is particularly valuable for identifying gradual land subsidence associated with groundwater extraction, underground mining, urban development, and natural geological processes. It can also help monitor volcanic inflation, earthquake-related displacement, and unstable slopes.<\/p>\n<p>The central principle is that a radar signal records both amplitude and phase. Amplitude describes the strength of the returning signal, while phase describes the position of the signal within its wave cycle. When two compatible radar acquisitions cover the same ground, changes in their phase can reveal subtle changes in the radar-to-ground path length.<\/p>\n<p>Under favorable conditions, carefully processed InSAR observations can resolve millimeter-scale changes. That does not mean every Sentinel-1 interferogram has millimeter accuracy: coherence, atmospheric conditions, processing choices, and reference stability all affect the reliability of a measurement.<\/p>\n<p>&nbsp;<\/p>\n<h3><span class=\"ez-toc-section\" id=\"11_Understanding_interferometric_phase\"><\/span>1.1 Understanding interferometric phase<span class=\"ez-toc-section-end\"><\/span><\/h3>\n<p>A Sentinel-1 Single Look Complex (SLC) image stores complex-valued radar samples. A simplified representation of a pixel is:<\/p>\n<div><span data-math-source=\"S=Ae^{i\\phi}\" data-math-display=\"true\">\\[ S=Ae^{i\\phi} \\]<\/span><\/div>\n<p>where:<\/p>\n<ul>\n<li>\n<div>\n<div>\n<p><span data-math-source=\"S\" data-math-display=\"false\">\\(S\\)<\/span> is the complex radar signal.<\/p>\n<\/div>\n<\/div>\n<\/li>\n<li>\n<div>\n<div>\n<p><span data-math-source=\"A\" data-math-display=\"false\">\\(A\\)<\/span> is its amplitude.<\/p>\n<\/div>\n<\/div>\n<\/li>\n<li>\n<div>\n<div>\n<p><span data-math-source=\"\\phi\" data-math-display=\"false\">\\(\\phi\\)<\/span> is its phase, measured in radians.<\/p>\n<\/div>\n<\/div>\n<\/li>\n<li>\n<div>\n<div>\n<p><span data-math-source=\"i=\\sqrt{-1}\" data-math-display=\"false\">\\(i=\\sqrt{-1}\\)<\/span>.<\/p>\n<\/div>\n<\/div>\n<\/li>\n<\/ul>\n<p>Suppose two SLC acquisitions observe the same area: a master image acquired at time <span data-math-source=\"t_1\" data-math-display=\"false\">\\(t_1\\)<\/span> and a secondary image acquired at time <span data-math-source=\"t_2\" data-math-display=\"false\">\\(t_2\\)<\/span>.<\/p>\n<p>Their complex interferogram is commonly formed as:<\/p>\n<div><span data-math-source=\"I=S_1S_2^*\" data-math-display=\"true\">\\[ I=S_1S_2^* \\]<\/span><\/div>\n<p>where <span data-math-source=\"S_2^*\" data-math-display=\"false\">\\(S_2^*\\)<\/span> is the complex conjugate of the secondary image.<\/p>\n<p>The wrapped interferometric phase is:<\/p>\n<div><span data-math-source=\"\\Delta\\phi_{\\mathrm{wrapped}} =\\arg(I)\" data-math-display=\"true\">\\[ \\Delta\\phi_{\\mathrm{wrapped}} =\\arg(I) \\]<\/span><\/div>\n<p>This phase difference is initially restricted to a principal interval such as <span data-math-source=\"[-\\pi,\\pi]\" data-math-display=\"false\">\\([-\\pi,\\pi]\\)<\/span>. It contains contributions from ground displacement, terrain elevation, atmospheric propagation delays, orbit inaccuracies, and noise.<\/p>\n<p>The processing objective is to isolate the deformation-related component from these other contributions.<\/p>\n<h3><span class=\"ez-toc-section\" id=\"12_Converting_phase_into_displacement\"><\/span>1.2 Converting phase into displacement<span class=\"ez-toc-section-end\"><\/span><\/h3>\n<p>Sentinel-1 operates in the C-band, with a radar wavelength of approximately 5.6 cm. Because the radar signal travels from the satellite to the ground and back again, a change in the one-way line-of-sight distance of half a wavelength produces approximately one complete <span data-math-source=\"2\\pi\" data-math-display=\"false\">\\(2\\pi\\)<\/span> phase cycle.<\/p>\n<p>After removing the modeled topographic contribution and unwrapping the phase, the line-of-sight displacement can be expressed as:<\/p>\n<div><span data-math-source=\"d_{\\mathrm{LOS}} =-\\frac{\\lambda}{4\\pi}\\Delta\\phi_{\\mathrm{def}}\" data-math-display=\"true\">\\[ d_{\\mathrm{LOS}} =-\\frac{\\lambda}{4\\pi}\\Delta\\phi_{\\mathrm{def}} \\]<\/span><\/div>\n<p>where:<\/p>\n<ul>\n<li>\n<div>\n<div>\n<p><span data-math-source=\"d_{\\mathrm{LOS}}\" data-math-display=\"false\">\\(d_{\\mathrm{LOS}}\\)<\/span> is the displacement along the satellite&#8217;s radar line of sight, in metres.<\/p>\n<\/div>\n<\/div>\n<\/li>\n<li>\n<div>\n<div>\n<p><span data-math-source=\"\\lambda\" data-math-display=\"false\">\\(\\lambda\\)<\/span> is the radar wavelength, in metres.<\/p>\n<\/div>\n<\/div>\n<\/li>\n<li>\n<div>\n<div>\n<p><span data-math-source=\"\\Delta\\phi_{\\mathrm{def}}\" data-math-display=\"false\">\\(\\Delta\\phi_{\\mathrm{def}}\\)<\/span> is the deformation-related unwrapped phase, in radians.<\/p>\n<\/div>\n<\/div>\n<\/li>\n<\/ul>\n<p>The negative sign assumes a particular master-secondary interferogram convention. Different software or image ordering can reverse the phase sign, so verify the convention before interpreting positive and negative displacement.<\/p>\n<p>For a wavelength of approximately 5.6 cm, one full <span data-math-source=\"2\\pi\" data-math-display=\"false\">\\(2\\pi\\)<\/span> phase cycle represents about 2.8 cm of line-of-sight displacement. Phase unwrapping is essential when accumulated phase changes exceed one cycle.<\/p>\n<p>Important distinction: InSAR measures displacement projected onto the satellite&#8217;s line of sight, not automatically vertical subsidence. Horizontal ground movement can contribute to the same measurement.<\/p>\n<h3><span class=\"ez-toc-section\" id=\"13_The_processing_workflow_at_a_glance\"><\/span>1.3 The processing workflow at a glance<span class=\"ez-toc-section-end\"><\/span><\/h3>\n<div>\n<div>\n<p>1. Acquire SLC data<\/p>\n<p>Select compatible Sentinel-1 acquisitions and a DEM<\/p>\n<\/div>\n<div><\/div>\n<div>\n<p>2. Coregister acquisitions<\/p>\n<p>Orbit correction \u2192 TOPS split \u2192 Back-geocoding \u2192 ESD<\/p>\n<\/div>\n<div><\/div>\n<div>\n<p>3. Generate differential interferogram<\/p>\n<p>Interferogram \u2192 Coherence \u2192 Topographic phase removal<\/p>\n<\/div>\n<div><\/div>\n<div>\n<p>4. Unwrap and calculate displacement<\/p>\n<p>Phase filtering \u2192 SNAPHU \u2192 LOS displacement<\/p>\n<\/div>\n<div><\/div>\n<div>\n<p>5. Map and assess<\/p>\n<p>Geocode \u2192 Export GeoTIFF \u2192 <a href=\"https:\/\/helpwithmygisassignment.com\/blog\/building-real-time-spatial-dashboards-in-qgis-with-geonode-and-postgis\/\">QGIS<\/a> \u2192 Validate subsidence hotspots<\/p>\n<\/div>\n<\/div>\n<p>Figure 1. A differential InSAR workflow for Sentinel-1 Interferometric Wide Swath (IW) TOPS data.<\/p>\n<h2><span class=\"ez-toc-section\" id=\"2_Data_requirements_Sentinel-1_SLC_image_pairs\"><\/span>2. Data requirements: Sentinel-1 SLC image pairs<span class=\"ez-toc-section-end\"><\/span><\/h2>\n<h3><span class=\"ez-toc-section\" id=\"21_Download_suitable_Sentinel-1_acquisitions\"><\/span>2.1 Download suitable Sentinel-1 acquisitions<span class=\"ez-toc-section-end\"><\/span><\/h3>\n<p>Obtain Sentinel-1 Level-1 Single Look Complex products from the Copernicus Data Space Ecosystem data catalogue.<\/p>\n<p>Open the <a href=\"https:\/\/dataspace.copernicus.eu\/?utm_source=chatgpt.com\" target=\"_blank\" rel=\"noopener\">Copernicus Data Space Ecosystem<\/a> and search for products covering the area of interest. The catalogue and its OData API support searching Sentinel-1 products by acquisition time, product type, and other metadata.<\/p>\n<div>\n<div><\/div>\n<\/div>\n<p>For conventional land subsidence mapping, prioritize the following settings:<\/p>\n<div>\n<table>\n<tbody>\n<tr>\n<th>Selection criterion<\/th>\n<th>Recommended choice<\/th>\n<\/tr>\n<tr>\n<td>Product type<\/td>\n<td>Level-1 SLC<\/td>\n<\/tr>\n<tr>\n<td>Acquisition mode<\/td>\n<td>IW (Interferometric Wide Swath), where available<\/td>\n<\/tr>\n<tr>\n<td>Polarization<\/td>\n<td>A common polarization, often VV<\/td>\n<\/tr>\n<tr>\n<td>Orbit direction<\/td>\n<td>Same ascending or descending pass<\/td>\n<\/tr>\n<tr>\n<td>Relative orbit<\/td>\n<td>Same track and acquisition geometry<\/td>\n<\/tr>\n<tr>\n<td>Spatial coverage<\/td>\n<td>Both acquisitions must cover the same target area and compatible bursts<\/td>\n<\/tr>\n<tr>\n<td>Acquisition interval<\/td>\n<td>Prefer a short temporal separation when possible<\/td>\n<\/tr>\n<tr>\n<td>Perpendicular baseline<\/td>\n<td>Prefer a sufficiently small baseline to limit geometric decorrelation and topographic sensitivity<\/td>\n<\/tr>\n<\/tbody>\n<\/table>\n<\/div>\n<p>Short temporal separation helps preserve coherence because the surface has less time to change. A smaller perpendicular baseline generally reduces geometric decorrelation and sensitivity to elevation differences. However, the best pair also depends on the landscape, acquisition geometry, and objectives; there is no single baseline threshold that works for every site.<\/p>\n<p>For long-term subsidence monitoring, plan a multi-date stack rather than relying on one image pair. A pair can estimate displacement between two acquisition dates, but it cannot by itself establish whether the motion followed a steady annual trend.<\/p>\n<h3><span class=\"ez-toc-section\" id=\"22_Why_SLC_is_essential\"><\/span>2.2 Why SLC is essential<span class=\"ez-toc-section-end\"><\/span><\/h3>\n<p>Sentinel-1 products come in different processing levels. For interferometry, SLC preserves complex phase information, whereas Ground Range Detected (GRD) products are detected products that do not retain the complex phase required for conventional interferometric processing.<\/p>\n<p>An SLC product typically includes a <code>.SAFE<\/code> product directory containing measurement data and metadata. Sentinel-1 IW SLC products are the standard input for the TOPS interferometric workflow described below.<\/p>\n<p>Useful official resources include the <a href=\"https:\/\/documentation.dataspace.copernicus.eu\/APIs\/OData.html?utm_source=chatgpt.com\" target=\"_blank\" rel=\"noopener\">Copernicus data access documentation<\/a> and the <a href=\"https:\/\/step.esa.int\/main\/doc\/tutorials\/?utm_source=chatgpt.com\" target=\"_blank\" rel=\"noopener\">ESA SNAP tutorials<\/a>.<\/p>\n<h3><span class=\"ez-toc-section\" id=\"23_Prepare_the_supporting_datasets\"><\/span>2.3 Prepare the supporting datasets<span class=\"ez-toc-section-end\"><\/span><\/h3>\n<p>In addition to the SLC pair, prepare:<\/p>\n<ul>\n<li>\n<div>\n<div>\n<p>A digital elevation model (DEM): Use a suitable source such as Copernicus DEM or SRTM to model and remove topographic phase. DEM accuracy affects the residual topographic signal.<\/p>\n<\/div>\n<\/div>\n<\/li>\n<li>\n<div>\n<div>\n<p>Orbit files: Apply the best suitable orbit information available for each acquisition, typically precise orbit files when available.<\/p>\n<\/div>\n<\/div>\n<\/li>\n<li>\n<div>\n<div>\n<p>A study-area boundary: Optional, but useful for subsetting the images and limiting processing costs.<\/p>\n<\/div>\n<\/div>\n<\/li>\n<li>\n<div>\n<div>\n<p>A reference area or stable reference point: Needed to interpret displacement relative to a consistent location.<\/p>\n<\/div>\n<\/div>\n<\/li>\n<li>\n<div>\n<div>\n<p>Independent validation data: GNSS measurements, leveling observations, or other suitable ground evidence are valuable for checking the inferred displacement.<\/p>\n<\/div>\n<\/div>\n<\/li>\n<\/ul>\n<p>Organize the inputs by acquisition date and preserve the original products. Avoid mixing incompatible tracks or silently changing the selected burst extent between the two acquisitions.<\/p>\n<h2><span class=\"ez-toc-section\" id=\"3_Step-by-step_processing_workflow_in_SNAP\"><\/span>3. Step-by-step processing workflow in SNAP<span class=\"ez-toc-section-end\"><\/span><\/h2>\n<p>The following procedure uses the classic Sentinel-1 IW TOPS differential InSAR workflow in ESA&#8217;s SNAP desktop application. Menu labels may vary by SNAP version; where necessary, use the operator search or Graph Builder to configure the processing chain.<\/p>\n<h3><span class=\"ez-toc-section\" id=\"Step_1_Apply_orbit_corrections_and_select_matching_bursts\"><\/span>Step 1: Apply orbit corrections and select matching bursts<span class=\"ez-toc-section-end\"><\/span><\/h3>\n<p>Purpose: Improve the satellite positioning information and isolate the same ground coverage in both acquisitions.<\/p>\n<ol>\n<li>\n<div>\n<div>\n<p>Open SNAP and load both Sentinel-1 SLC <code>.SAFE<\/code> products using File \u2192 Open Product.<\/p>\n<\/div>\n<\/div>\n<\/li>\n<li>\n<div>\n<div>\n<p>Run Radar \u2192 Apply Orbit File on each acquisition. Select the best appropriate orbit file available.<\/p>\n<\/div>\n<\/div>\n<\/li>\n<li>\n<div>\n<div>\n<p>Run the <code>TOPSAR-Split<\/code> operator on both products.<\/p>\n<\/div>\n<\/div>\n<\/li>\n<li>\n<div>\n<div>\n<p>Select the same subswath, such as IW2, and the same polarization, such as VV.<\/p>\n<\/div>\n<\/div>\n<\/li>\n<li>\n<div>\n<div>\n<p>Select matching bursts that cover the study area. Apply the same selection to the master and secondary products.<\/p>\n<\/div>\n<\/div>\n<\/li>\n<\/ol>\n<p>The master provides the reference geometry. The secondary acquisition will subsequently be resampled into alignment with it.<\/p>\n<p>Do not select burst ranges independently just because they cover similar-looking areas. The same radar coverage must be available in both acquisitions for valid interferometric processing.<\/p>\n<h3><span class=\"ez-toc-section\" id=\"Step_2_Coregister_the_SLC_images\"><\/span>Step 2: Coregister the SLC images<span class=\"ez-toc-section-end\"><\/span><\/h3>\n<p>Coregistration aligns corresponding pixels from the two acquisitions. Small misalignments can contaminate the interferometric phase, so this step is particularly important for Sentinel-1 TOPS data.<\/p>\n<p>In SNAP&#8217;s processing chain, use the following operators:<\/p>\n<p><code>Back-Geocoding<\/code> \u2192 <code>Enhanced-Spectral-Diversity<\/code> (ESD)<\/p>\n<ol>\n<li>\n<div>\n<div>\n<p>Run <code>Back-Geocoding<\/code> with the master and secondary split products.<\/p>\n<\/div>\n<\/div>\n<\/li>\n<li>\n<div>\n<div>\n<p>Select an appropriate DEM for the processing area.<\/p>\n<\/div>\n<\/div>\n<\/li>\n<li>\n<div>\n<div>\n<p>Confirm that the operator completes successfully and generates a coregistered product.<\/p>\n<\/div>\n<\/div>\n<\/li>\n<li>\n<div>\n<div>\n<p>Run <code>Enhanced-Spectral-Diversity<\/code> on the coregistered product to estimate and correct residual azimuth misregistration between TOPS bursts.<\/p>\n<\/div>\n<\/div>\n<\/li>\n<\/ol>\n<p>ESD is an important part of the conventional Sentinel-1 TOPS interferometry workflow because residual azimuth misregistration can introduce phase errors at burst boundaries. ESA&#8217;s Sentinel-1 tutorial covers TOPS coregistration and the subsequent interferometric operators.<\/p>\n<p>Quality check: Inspect the coregistration results before continuing. Failed ESD convergence, incompatible bursts, or poor overlap can compromise the entire deformation product.<\/p>\n<h3><span class=\"ez-toc-section\" id=\"Step_3_Generate_the_interferogram_and_coherence_raster\"><\/span>Step 3: Generate the interferogram and coherence raster<span class=\"ez-toc-section-end\"><\/span><\/h3>\n<p>Purpose: Calculate the wrapped phase difference and determine where the two acquisitions retain a sufficiently consistent radar signal.<\/p>\n<p>Run SNAP&#8217;s <code>Interferogram<\/code> operator on the coregistered product.<\/p>\n<p>Configure the available options to:<\/p>\n<ul>\n<li>\n<div>\n<div>\n<p>Subtract the modeled flat-earth phase.<\/p>\n<\/div>\n<\/div>\n<\/li>\n<li>\n<div>\n<div>\n<p>Generate the coherence band.<\/p>\n<\/div>\n<\/div>\n<\/li>\n<li>\n<div>\n<div>\n<p>Preserve the interferometric phase output.<\/p>\n<\/div>\n<\/div>\n<\/li>\n<\/ul>\n<p>The interferogram operator combines the complex radar observations to calculate the phase difference between the master and secondary acquisitions.<\/p>\n<p>Coherence measures the consistency of the complex signals within a local estimation window. It is commonly represented by values between 0 and 1, where higher values indicate more reliable phase information and lower values indicate greater decorrelation.<\/p>\n<p>Next, run <code>TOPSAR-Deburst<\/code> if your workflow requires a continuous image across the selected TOPS bursts.<\/p>\n<p>Inspect both the wrapped phase and coherence layers. Coherent built-up areas may provide useful measurement points, while dense vegetation, changing agricultural fields, and water often produce less reliable phase measurements.<\/p>\n<div>\n<div><\/div>\n<\/div>\n<p>A wrapped interferogram should not be interpreted as a displacement map yet. Its colors represent cyclic phase values, and those values still contain topographic and other unwanted contributions.<\/p>\n<h3><span class=\"ez-toc-section\" id=\"Step_4_Remove_the_topographic_phase\"><\/span>Step 4: Remove the topographic phase<span class=\"ez-toc-section-end\"><\/span><\/h3>\n<p>Purpose: Isolate phase variations associated with ground deformation by subtracting the phase predicted from terrain elevation and acquisition geometry.<\/p>\n<p>Run SNAP&#8217;s <code>TopoPhaseRemoval<\/code> operator using the selected DEM and interferometric product.<\/p>\n<p>Conceptually, the observed phase can be decomposed as:<\/p>\n<div><span data-math-source=\"\\phi_{\\mathrm{obs}} =\\phi_{\\mathrm{flat}} +\\phi_{\\mathrm{topo}} +\\phi_{\\mathrm{def}} +\\phi_{\\mathrm{atm}} +\\phi_{\\mathrm{orb}} +\\phi_{\\mathrm{noise}}\" data-math-display=\"true\">\\[ \\phi_{\\mathrm{obs}} =\\phi_{\\mathrm{flat}} +\\phi_{\\mathrm{topo}} +\\phi_{\\mathrm{def}} +\\phi_{\\mathrm{atm}} +\\phi_{\\mathrm{orb}} +\\phi_{\\mathrm{noise}} \\]<\/span><\/div>\n<p>where the terms represent flat-earth geometry, topography, deformation, atmospheric effects, orbit-related errors, and noise.<\/p>\n<p>After flat-earth and modeled topographic phase corrections, the residual interferogram should primarily represent deformation, atmospheric effects, residual orbit or elevation errors, and noise.<\/p>\n<p>Important: Topographic phase removal is not optional for this differential deformation workflow. Omitting it can cause terrain-related phase patterns to be mistaken for subsidence or uplift.<\/p>\n<p>A higher-resolution DEM may improve the terrain representation, but DEM quality, alignment, and height reference also matter. Residual elevation errors can still generate spurious phase signals.<\/p>\n<h3><span class=\"ez-toc-section\" id=\"Step_5_Filter_the_phase_and_prepare_for_unwrapping\"><\/span>Step 5: Filter the phase and prepare for unwrapping<span class=\"ez-toc-section-end\"><\/span><\/h3>\n<p>Purpose: Reduce phase noise before reconstructing a continuous phase surface.<\/p>\n<p>Run SNAP&#8217;s <code>GoldsteinPhaseFiltering<\/code> operator on the differential interferogram.<\/p>\n<p>Phase filtering can improve the local continuity of coherent phase patterns, but aggressive filtering may erase small-scale deformation features or smooth genuine displacement gradients.<\/p>\n<p>Keep the coherence layer available so that poorly measured pixels can be excluded or flagged.<\/p>\n<p>For initial investigations, a coherence threshold around 0.3\u20130.4 can be a useful trial setting, not a universal scientific cutoff. Select the actual threshold based on the land cover, coherence distribution, and acceptable measurement uncertainty.<\/p>\n<h3><span class=\"ez-toc-section\" id=\"Step_6_Unwrap_the_interferometric_phase_with_SNAPHU\"><\/span>Step 6: Unwrap the interferometric phase with SNAPHU<span class=\"ez-toc-section-end\"><\/span><\/h3>\n<p>The wrapped phase only records values modulo <span data-math-source=\"2\\pi\" data-math-display=\"false\">\\(2\\pi\\)<\/span>. Phase unwrapping estimates the integer number of cycles that must be added or subtracted to construct a continuous phase field.<\/p>\n<p>SNAPHU is a statistical-cost, network-flow algorithm designed for two-dimensional phase unwrapping. Its statistical model and parameter configuration should reflect whether the interferogram is being processed for topography, deformation, or another application.<\/p>\n<p>Install the SNAPHU plugin in SNAP if it is not already available:<\/p>\n<p>Tools \u2192 Plugins \u2192 Available Plugins \u2192 SNAPHU Unwrapping<\/p>\n<p>Then follow this sequence:<\/p>\n<ol>\n<li>\n<div>\n<div>\n<p>Run Radar \u2192 Interferometric \u2192 Unwrapping \u2192 Snaphu Export on the filtered differential interferogram.<\/p>\n<\/div>\n<\/div>\n<\/li>\n<li>\n<div>\n<div>\n<p>Select the appropriate deformation processing mode where the operator offers a statistical-cost mode, and configure the output folder and other required parameters.<\/p>\n<\/div>\n<\/div>\n<\/li>\n<li>\n<div>\n<div>\n<p>Run Snaphu-unwrapping using the exported input and configuration files.<\/p>\n<\/div>\n<\/div>\n<\/li>\n<li>\n<div>\n<div>\n<p>Preserve the exported files and unwrapping logs.<\/p>\n<\/div>\n<\/div>\n<\/li>\n<li>\n<div>\n<div>\n<p>Run Snaphu Import to bring the unwrapped result back into SNAP.<\/p>\n<\/div>\n<\/div>\n<\/li>\n<\/ol>\n<div>\n<div>\n<div><\/div>\n<\/div>\n<\/div>\n<p>Inspect the unwrapped phase for abrupt jumps, large isolated patches, and discontinuities across low-coherence regions. Unwrapping algorithms can produce plausible-looking surfaces even when the solution is wrong, particularly where phase information is sparse or noisy.<\/p>\n<p>Do not attempt to interpret disconnected, unreliable areas as measured displacement. Mask or flag those areas and document their limitations.<\/p>\n<h3><span class=\"ez-toc-section\" id=\"Step_7_Convert_unwrapped_phase_to_displacement\"><\/span>Step 7: Convert unwrapped phase to displacement<span class=\"ez-toc-section-end\"><\/span><\/h3>\n<p>Once the unwrapped phase has been imported, use SNAP&#8217;s <code>Phase to Displacement<\/code> operator.<\/p>\n<p>The operator converts phase to line-of-sight displacement using the radar wavelength. The output represents a relative displacement measurement between acquisitions, subject to the interferogram convention and reference choice.<\/p>\n<p>For the conventional phase-to-displacement relationship:<\/p>\n<div><span data-math-source=\"d_{\\mathrm{LOS}} =-\\frac{\\lambda}{4\\pi}\\phi_{\\mathrm{unwrapped}}\" data-math-display=\"true\">\\[ d_{\\mathrm{LOS}} =-\\frac{\\lambda}{4\\pi}\\phi_{\\mathrm{unwrapped}} \\]<\/span><\/div>\n<p>Here, the input must be the appropriate unwrapped deformation phase\u2014not the original wrapped phase.<\/p>\n<p>After conversion:<\/p>\n<ol>\n<li>\n<div>\n<div>\n<p>Inspect the displacement raster and its units.<\/p>\n<\/div>\n<\/div>\n<\/li>\n<li>\n<div>\n<div>\n<p>Identify an appropriate stable reference location.<\/p>\n<\/div>\n<\/div>\n<\/li>\n<li>\n<div>\n<div>\n<p>Reference the displacement to that location if needed.<\/p>\n<\/div>\n<\/div>\n<\/li>\n<li>\n<div>\n<div>\n<p>Run Radar \u2192 Geometric \u2192 Terrain Correction \u2192 Range-Doppler Terrain Correction to geocode the displacement and coherence products.<\/p>\n<\/div>\n<\/div>\n<\/li>\n<\/ol>\n<p>Select the desired map CRS and pixel spacing, and export the geocoded results as GeoTIFFs.<\/p>\n<p>The resulting deformation raster is a line-of-sight displacement map, not automatically a vertical elevation-change map. This distinction must remain explicit in the final GIS layer and map legend.<\/p>\n<h2><span class=\"ez-toc-section\" id=\"4_Track_land_subsidence_in_QGIS\"><\/span>4. Track land subsidence in QGIS<span class=\"ez-toc-section-end\"><\/span><\/h2>\n<p>Once the geocoded GeoTIFFs are ready, QGIS can be used to visualize displacement, calculate an indicative rate, and delineate potential subsidence hotspots.<\/p>\n<h3><span class=\"ez-toc-section\" id=\"Step_1_Load_the_rasters\"><\/span>Step 1: Load the rasters<span class=\"ez-toc-section-end\"><\/span><\/h3>\n<p>Open QGIS and navigate to Layer \u2192 Add Layer \u2192 Add Raster Layer.<\/p>\n<p>Load the following:<\/p>\n<ul>\n<li>\n<div>\n<div>\n<p>The geocoded LOS displacement raster.<\/p>\n<\/div>\n<\/div>\n<\/li>\n<li>\n<div>\n<div>\n<p>The geocoded coherence raster.<\/p>\n<\/div>\n<\/div>\n<\/li>\n<li>\n<div>\n<div>\n<p>An optional administrative boundary, groundwater extraction dataset, mine boundary, or other reference layer.<\/p>\n<\/div>\n<\/div>\n<\/li>\n<\/ul>\n<p>Check the raster CRS, pixel dimensions, NoData values, and reported displacement units before performing calculations.<\/p>\n<h3><span class=\"ez-toc-section\" id=\"Step_2_Calculate_an_annualized_LOS_displacement_rate\"><\/span>Step 2: Calculate an annualized LOS displacement rate<span class=\"ez-toc-section-end\"><\/span><\/h3>\n<p>A single interferometric pair measures displacement accumulated between two acquisition dates. To express that displacement as an indicative annualized rate, divide it by the time interval.<\/p>\n<p>If the displacement raster stores metres, calculate:<\/p>\n<div><span data-math-source=\"v_{\\mathrm{LOS}} =\\frac{100d_{\\mathrm{LOS}}}{\\Delta t} \\times365.25\" data-math-display=\"true\">\\[ v_{\\mathrm{LOS}} =\\frac{100d_{\\mathrm{LOS}}}{\\Delta t} \\times365.25 \\]<\/span><\/div>\n<p>where:<\/p>\n<ul>\n<li>\n<div>\n<div>\n<p><span data-math-source=\"v_{\\mathrm{LOS}}\" data-math-display=\"false\">\\(v_{\\mathrm{LOS}}\\)<\/span> is in centimetres per year.<\/p>\n<\/div>\n<\/div>\n<\/li>\n<li>\n<div>\n<div>\n<p><span data-math-source=\"d_{\\mathrm{LOS}}\" data-math-display=\"false\">\\(d_{\\mathrm{LOS}}\\)<\/span> is in metres.<\/p>\n<\/div>\n<\/div>\n<\/li>\n<li>\n<div>\n<div>\n<p><span data-math-source=\"\\Delta t\" data-math-display=\"false\">\\(\\Delta t\\)<\/span> is the number of days between acquisitions.<\/p>\n<\/div>\n<\/div>\n<\/li>\n<\/ul>\n<p>In Raster \u2192 Raster Calculator, an example for a 12-day interval is:<\/p>\n<div>\n<div>\n<pre dir=\"ltr\"><code>\"displacement_m@1\" * 100 * 365.25 \/ 12<\/code><\/pre>\n<\/div>\n<\/div>\n<p>Replace <code>displacement_m@1<\/code> with the actual band name and <code>12<\/code> with the true acquisition interval in days.<\/p>\n<p>This calculation assumes that the motion was approximately linear throughout the period. It does not establish a long-term annual trend from a single pair. For defensible annual velocity estimates, process a multi-date stack using a time-series method, such as Small Baseline Subset (SBAS) or Persistent Scatterer InSAR (PSInSAR), where appropriate. GMTSAR provides documented workflows for batch interferograms and SBAS time-series processing.<\/p>\n<h3><span class=\"ez-toc-section\" id=\"Step_3_Apply_a_coherence_mask\"><\/span>Step 3: Apply a coherence mask<span class=\"ez-toc-section-end\"><\/span><\/h3>\n<p>Poor-coherence pixels should not be treated as reliable displacement measurements.<\/p>\n<p>Use the coherence raster to mask low-quality areas before calculating hotspot statistics. For example, test a threshold of 0.35, inspect the resulting mask, and adjust it to match the study area&#8217;s land cover and measurement requirements.<\/p>\n<p>Do not confuse a low coherence value with a subsidence signal. Coherence describes measurement reliability, not the amount of ground movement.<\/p>\n<h3><span class=\"ez-toc-section\" id=\"Step_4_Style_the_deformation_map\"><\/span>Step 4: Style the deformation map<span class=\"ez-toc-section-end\"><\/span><\/h3>\n<p>In QGIS, open Layer Properties \u2192 Symbology and select a single-band pseudocolor renderer.<\/p>\n<p>Use a diverging color ramp centered around zero, with one side representing motion toward the satellite and the other representing motion away from it. Confirm the processor&#8217;s sign convention before assigning those labels; do not automatically label negative values as subsidence.<\/p>\n<p>Overlay roads, buildings, mine boundaries, geological units, and groundwater extraction locations where suitable data are available.<\/p>\n<p>For hotspot identification, consider both the magnitude and spatial continuity of the signal. A coherent deformation bowl may be more informative than isolated extreme pixels, which can reflect phase-unwrapping errors or atmospheric artifacts.<\/p>\n<p>Finally, export the styled layer or layout as a map for reporting. Keep the displacement units, acquisition dates, reference point, coherence threshold, and line-of-sight interpretation in the map metadata or accompanying notes.<\/p>\n<h2><span class=\"ez-toc-section\" id=\"5_Limitations_and_quality_control\"><\/span>5. Limitations and quality control<span class=\"ez-toc-section-end\"><\/span><\/h2>\n<div>\n<table>\n<tbody>\n<tr>\n<th>Limitation<\/th>\n<th>Effect on the result<\/th>\n<th>Recommended response<\/th>\n<\/tr>\n<tr>\n<td>Temporal decorrelation<\/td>\n<td>Surface changes between acquisitions reduce phase reliability<\/td>\n<td>Prefer shorter time intervals and coherent targets<\/td>\n<\/tr>\n<tr>\n<td>Vegetation<\/td>\n<td>Changing leaves, soil moisture, and canopy structure can reduce coherence<\/td>\n<td>Mask unreliable areas and consider multi-temporal approaches<\/td>\n<\/tr>\n<tr>\n<td>Atmospheric delay<\/td>\n<td>Differences in atmospheric conditions introduce phase patterns unrelated to ground motion<\/td>\n<td>Examine multiple dates and use suitable atmospheric correction methods<\/td>\n<\/tr>\n<tr>\n<td>DEM error<\/td>\n<td>Residual topographic phase can resemble deformation<\/td>\n<td>Use a suitable DEM and inspect residual terrain-related patterns<\/td>\n<\/tr>\n<tr>\n<td>Phase-unwrapping errors<\/td>\n<td>Incorrect cycle counts cause false displacement gradients or offsets<\/td>\n<td>Inspect coherence, unwrapped-phase continuity, and reference areas<\/td>\n<\/tr>\n<tr>\n<td>Orbit errors<\/td>\n<td>Long-wavelength phase ramps can bias displacement<\/td>\n<td>Use suitable orbit products and assess spatial trends<\/td>\n<\/tr>\n<tr>\n<td>LOS geometry<\/td>\n<td>Horizontal and vertical motions contribute differently to the measured displacement<\/td>\n<td>Use ascending and descending observations or independent measurements when possible<\/td>\n<\/tr>\n<\/tbody>\n<\/table>\n<\/div>\n<h3><span class=\"ez-toc-section\" id=\"Can_LOS_displacement_be_converted_to_vertical_subsidence\"><\/span>Can LOS displacement be converted to vertical subsidence?<span class=\"ez-toc-section-end\"><\/span><\/h3>\n<p>Only under additional assumptions.<\/p>\n<p>If movement is assumed to be purely vertical, a simplified approximation is:<\/p>\n<div><span data-math-source=\"d_{\\mathrm{vertical}} \\approx\\frac{d_{\\mathrm{LOS}}}{\\cos\\theta}\" data-math-display=\"true\">\\[ d_{\\mathrm{vertical}} \\approx\\frac{d_{\\mathrm{LOS}}}{\\cos\\theta} \\]<\/span><\/div>\n<p>where <span data-math-source=\"\\theta\" data-math-display=\"false\">\\(\\theta\\)<\/span> is the radar incidence angle measured from the vertical. The sign still depends on the chosen convention.<\/p>\n<p>This approximation is not valid when significant horizontal movement contributes to the LOS measurement. Combining ascending and descending observations can help distinguish vertical from east-west motion, although north-south sensitivity remains limited for many satellite viewing geometries.<\/p>\n<p>For groundwater or mining investigations, an InSAR hotspot is evidence of a surface-motion pattern, not proof of its cause. Compare the mapped pattern with groundwater levels, extraction histories, mine records, local geology, and independent ground observations before attributing the deformation to a specific process.<\/p>\n<h2><span class=\"ez-toc-section\" id=\"6_Alternative_workflow_using_GMTSAR\"><\/span>6. Alternative workflow using GMTSAR<span class=\"ez-toc-section-end\"><\/span><\/h2>\n<p>GMTSAR is an open-source SAR processing system with command-line scripts for individual interferograms and multi-date time series. It is particularly useful when<\/p>\n<p>After installing GMTSAR and its dependencies, create a project directory:<\/p>\n<div>\n<div>\n<pre dir=\"ltr\"><code>mkdir -p insar_project\/{raw,topo}\r\ncd insar_project<\/code><\/pre>\n<\/div>\n<\/div>\n<p>Place the two Sentinel-1 SLC products and their orbit files in <code>raw\/<\/code>, and prepare the DEM in the format expected by GMTSAR in <code>topo\/<\/code>.<\/p>\n<p>A pair-processing command follows this general pattern:<\/p>\n<div>\n<div>\n<pre dir=\"ltr\"><code>p2p_processing.csh \\\r\n  raw\/MASTER.SAFE \\\r\n  raw\/SECONDARY.SAFE \\\r\n  config.tops.txt \\\r\n  &gt; p2p_processing.log 2&gt;&amp;1<\/code><\/pre>\n<\/div>\n<\/div>\n<p>Replace the product names and configuration file with the actual files. The example assumes a GMTSAR installation in which the relevant script and dependencies are available.<\/p>\n<p>Consult the <a href=\"https:\/\/github.com\/gmtsar\/documentation\/blob\/main\/source\/Quick_Start_Guide.rst?utm_source=chatgpt.com\" target=\"_blank\" rel=\"noopener\">GMTSAR Quick Start Guide<\/a> for directory structure, orbit files, DEM preparation, and Sentinel-1 pair processing. The configuration controls options such as phase unwrapping and geocoding, so review it before interpreting the output.<\/p>\n<p>For multi-date monitoring, use the documented Sentinel-1 batch preprocessing and SBAS workflows instead of simply averaging a few independently generated interferograms.<\/p>\n<h2><span class=\"ez-toc-section\" id=\"7_Summary\"><\/span>7. Summary<span class=\"ez-toc-section-end\"><\/span><\/h2>\n<p>InSAR surface deformation mapping with Sentinel-1 provides a repeatable way to identify potential ground subsidence using radar phase differences.<\/p>\n<p>The essential workflow is to obtain compatible SLC acquisitions, coregister them accurately, generate and correct the interferogram, unwrap its phase, convert the result to LOS displacement, and export a geocoded raster for analysis in QGIS.<\/p>\n<p>Reliable results depend on more than completing the processing steps. Coherence assessment, suitable DEM and orbit data, atmospheric-error analysis, careful reference selection, and independent validation are all necessary to distinguish genuine surface motion from processing artifacts.<\/p>\n<p>For initial screening, a carefully processed pair can reveal where deformation may be occurring. For long-term monitoring of groundwater extraction, urban subsidence, or mining hazards, a multi-temporal InSAR time series provides a stronger basis for estimating displacement rates and evaluating how the ground changes through time.<\/p>\n<\/div>\n","protected":false},"excerpt":{"rendered":"<p>InSAR Surface Deformation Mapping: Tracking Ground Subsidence with Sentinel-1 1. Introduction: How InSAR detects ground subsidence When performing InSAR surface deformation mapping Sentinel 1 subsidence, detecting millimeters of ground displacement from satellite radar requires careful processing. Interferometric Synthetic Aperture Radar (InSAR) is a remote sensing technique that measures changes in the distance between a radar &#8230; <a title=\"InSAR Surface Deformation Mapping: Tracking Ground Subsidence with Sentinel-1\" class=\"read-more\" href=\"https:\/\/helpwithmygisassignment.com\/blog\/insar-surface-deformation-mapping-sentinel-1-subsidence\/\" aria-label=\"Read more about InSAR Surface Deformation Mapping: Tracking Ground Subsidence with Sentinel-1\">Read more<\/a><\/p>\n","protected":false},"author":1,"featured_media":1542,"comment_status":"closed","ping_status":"closed","sticky":false,"template":"","format":"standard","meta":{"footnotes":""},"categories":[1,4],"tags":[],"class_list":["post-1389","post","type-post","status-publish","format-standard","has-post-thumbnail","hentry","category-uncategorized","category-gis-articles"],"_links":{"self":[{"href":"https:\/\/helpwithmygisassignment.com\/blog\/wp-json\/wp\/v2\/posts\/1389","targetHints":{"allow":["GET"]}}],"collection":[{"href":"https:\/\/helpwithmygisassignment.com\/blog\/wp-json\/wp\/v2\/posts"}],"about":[{"href":"https:\/\/helpwithmygisassignment.com\/blog\/wp-json\/wp\/v2\/types\/post"}],"author":[{"embeddable":true,"href":"https:\/\/helpwithmygisassignment.com\/blog\/wp-json\/wp\/v2\/users\/1"}],"replies":[{"embeddable":true,"href":"https:\/\/helpwithmygisassignment.com\/blog\/wp-json\/wp\/v2\/comments?post=1389"}],"version-history":[{"count":4,"href":"https:\/\/helpwithmygisassignment.com\/blog\/wp-json\/wp\/v2\/posts\/1389\/revisions"}],"predecessor-version":[{"id":1543,"href":"https:\/\/helpwithmygisassignment.com\/blog\/wp-json\/wp\/v2\/posts\/1389\/revisions\/1543"}],"wp:featuredmedia":[{"embeddable":true,"href":"https:\/\/helpwithmygisassignment.com\/blog\/wp-json\/wp\/v2\/media\/1542"}],"wp:attachment":[{"href":"https:\/\/helpwithmygisassignment.com\/blog\/wp-json\/wp\/v2\/media?parent=1389"}],"wp:term":[{"taxonomy":"category","embeddable":true,"href":"https:\/\/helpwithmygisassignment.com\/blog\/wp-json\/wp\/v2\/categories?post=1389"},{"taxonomy":"post_tag","embeddable":true,"href":"https:\/\/helpwithmygisassignment.com\/blog\/wp-json\/wp\/v2\/tags?post=1389"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}