Having never seen one for real - active, dormant or extinct - the first thing that comes to my mind when I think of a Volcano is this popular depiction from Chacha Chaudhary comics-
For those unfamiliar, Sabu was an alien from Jupiter who was residing on Earth with India's most famous comic book superhero - the diminutive yet witty old man - Chacha Chaudhary. Easily winning several strength-based duels with villainous thugs, Sabu sometimes faltered due to his lack of wit or when he encountered adversaries more powerful than himself. When the latter occurred, he would be subject to incessant needling from the antagonist which would prompt him to get into a Hulk-esque rage- symbolized by a volcano erupting on his home planet. In this beast mode, he was indomitable and pulverized the opposition. This is how comic books go, however, I can't help but draw parallels to the current war-like situation unfolding before us in real life.
Volcanic Eruptions have a mythical aura, a phenomenon which humans have always feared and revered at the same time. Unsurprisingly so, because even in my country India, which harbors only 2 active volcanoes and that too in the remote Andaman islands - we have an iconic name for this object of awe. In India's national language Hindi, a volcano is called 'Jwalamukhi' which can be translated as 'the Mouth that spurts Fire'.
Given my limited knowledge about volcanoes or volcanic activity, I am grateful to RUS Copernicus & EO College whose learning resources have been tremendously useful in developing my understanding of this topic, like it has for many topics before this.
For instance, it was quite fascinating to know that volcanoes do not explode all of a sudden, generally. Rather, to put it in my own way, volcanoes behave like living organisms - breathing in, breathing out, expanding and contracting in size, shape and volume, metabolizing energy within and only when the internal pressure is too hard to contain, does it explode in a violent expression of fury.
It is this 'lively' behavior which captures the attention of researchers and they have devised ways to capture it and interpret it in the branch of science known as volcanology. To clarify, this 'liveliness' (the technical term for this is deformation) is not discernible to the naked eye, rather it can be observed using geodetic monitoring techniques such as Leveling, GPS & Interferometric Synthetic Aperture Radar (InSAR).
We will use the InSAR technique on satellite imagery to study three volcanoes - 1) Pico Do Fogo (Cape Verde) 2014, 2) Cumbre Vieja (Canary Islands) 2021 & 3) Mt. Nyiragongo (Congo), 2021.
Before we progress further, I'd recommend you to read this insightful article published by USGS on Volcano deformation.
Mapping on-surface deformation phenomenon is relatively easier. The recent devastating eruption near Tonga happened from the submarine volcano - Hunga Tonga–Hunga Haʻapai below the surface of the ocean which made it difficult to detect and observe before the devastating eruption occurred on 15th January 2022.
The processing chain involved in extracting deformation (phase fringes) and precise displacement values (upswell and subsidence) from SAR Satellite Imagery is fairly extensive. Depiction below -
Below is a detailed video which will take you on a guided recording of performing each of these steps on a satellite imagery processing software (SNAP) to derive and visualize deformation and displacement map of Cumbre Vieja Eruption (Canary Islands) in 2021.
Do read on in case you do not wish to see the video.
Video 1: Mapping Deformation using Satellite Imagery - Cumbre Vieja volcano in 2021
To explain the above steps individually in a simple manner-
a) Read & Read(2)- Two sets of imagery are used - one prior to the eruption and second one post-eruption to do change detection.
b) TOPSAR-Split- SAR Imagery is heavy (Each set is around 1 GB in size). To reduce computational load while performing analysis on it, invariably one has to restrict the image to the area of interest. TOPSAR-Split helps us to select specific imagery sections and data bands from the entire set, eliminating the rest, thereby making the imagery less taxing on the software.
c) Apply (Updated) Orbit Files- Orbit Files within SAR Imagery contains satellite velocity and position information. By default, this may not be the updated / corrected version in the downloaded imagery. We have to update the orbit files to allow for precise pixel overlays in the next steps involving interferometry.
Please note that the visual output remains the same as step b).
d) Coregistration - Back Geocoding - This is a vital step in the Interferometric processing chain. It involves the factoring in of 'Digital Elevation Model' or DEM into the imagery. Interpreting radar reflectance / backscatter is only useful if we were to know the elevation values of the topographic surface that is our Area of Interest (AoI). This step also allows us to combine the two imagery sets into a single stack allowing for pixel to pixel comparison.
e) Coregistration - Enhanced Spectral Diversity - This optimization tool helps improve the precision of our geometrically-derived Coregistered Stack thereby minimizing pixel overlay errors.
f) Interferogram - A visual depiction which allow us to detect and quantify the ground
displacement that occurred between the pre and post volcanic imagery acquisitions with centimeter-level accuracy.
g) TOPSAR Deburst - Involves removing the black lines (data discontinuities) between the imagery sections (bursts) thereby making the entire imagery seamless.
h) Topographic Phase Removal - There are two types of phase information captured in the Interferogram - a) due to the topography and b) due to displacement. This step allows us to eliminate the topo-phase so that we are only left with displacement phase data values.
i) Multilooking - This step reduces the inherent speckle noise in the imagery (think of an Analog TVs dotty-grey screen error) and makes the pixels 'square' in shape allowing for easy interpretability.
j) Goldstein Phase Filtering - The next step is to apply Goldstein Phase Filtering. Like Multilooking, this tool is also applied to reduce the 'noise' in the image and make the data ready for the later steps which involve the calculation of precise displacement (deformation) values.
k) Snaphu Export - Snaphu (Phase) Unwrapping - Snaphu Import - Phase to Displacement - These four steps can be clubbed into one for the purpose of our understanding - let us call the combined step as Phase to Displacement. While the interferometric phase gives us an indication where land subsidence (shrinkage) and upswell (growth) has happened, phase unwrapping using Snaphu tool allows us to calculate precise displacement values between the two sets of imagery (pre & post).
l) (Range-Doppler) Terrain Correction - Alongside eliminating non-Land displacement values, 'Terrain Correction' also helps in making the imagery geometrically correct as the tilt of satellite among other aspects have inbuilt distortions. The correction helps the end product look proportionally as close as possible as it does in the real world.
m) Write (& Export) - the output for further analysis and visualization purposes
Phase Visualization on Google Earth
Video 2: The Foga phase layer is overlayed on Google Earth basemap. Strong deformations (phase fringes) are seen surrounding the volcano indicating loss of surface around the cone as well upswell due to magma accumulation in certain areas.
Video 3: The Mt. Nyiragongo phase layer is overlayed on Google Earth basemap. Strong deformations (phase fringes) are seen, not immediately around the volcano though. These can be attributed to persistent seismic (earthquakes) activity following the eruption - one such as validated here.
Displacement Visualization on Google Earth
Video 4: Fogo displacement values are displayed towards the top left in the depiction overlayed on Google Earth basemap. We can clearly visualize the areas most affected by land deformation (blue and dark green).
With this, we come to the conclusion of this article. Hope you found the process of mapping deformation from SAR imagery to study volcanic activity to be interesting!
ABOUT US
Intelloc Mapping Services | Mapmyops.com is based in Kolkata, India and engages in providing Mapping solutions that can be integrated with Operations Planning, Design and Audit workflows. These include but are not limited to - Drone Services, Subsurface Mapping Services, Location Analytics & App Development, Supply Chain Services & Remote Sensing Services. The services can be rendered pan-India, some even globally, and will aid an organization to meet its stated objectives especially pertaining to Operational Excellence, Cost Reduction, Sustainability and Growth.
Broadly, our area of expertise can be split into two categories - Geographic Mapping and Operations Mapping. The Infographic below highlights our capabilities.
Our 'Mapping for Operations'-themed workflow demonstrations can be accessed from the firm's Website / YouTube Channel and an overview can be obtained from this flyer. Happy to address queries and respond to documented requirements. Custom Demonstration, Training & Trials are facilitated only on a paid-basis. Looking forward to being of service.
Regards,