Volcano Monitoring - Mapping Deformation using Satellite Imagery
This article was originally posted on my firm's website:?https://www.mapmyops.com/volcano-monitoring-mapping-deformation-using-satellite-imagery?and is updated there when an edit needs to be made. Hence, I'd recommend you to read this article from the hyperlink above.
________________________________________________________________
Having never seen one for real, active or otherwise, the first thing that comes to my mind when I think of a Volcano is this depiction from the popular Chacha Chaudhary comics -
For those unfamiliar, Sabu was an alien from Jupiter who was faithful to one of India's most famous comic superheroes - the diminutive yet witty old man - Chacha Chaudhary. Easily winning several strength based duels with villainous thugs, Sabu often faltered due to lack of wit and, on certain occasions, while encountering even more powerful adversaries than him. When the latter occurred, he would be subject to incessant needling which prompted him to get into a Hulk-esque rage- symbolized by a volcano erupting somewhere on his home planet. In 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.
__________________________________________________________________________
Volcanic Eruptions have a mythical aura, a phenomenon which humans have always feared and revered at the same time. Unsurprisingly so, because even in our 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 our national language Hindi - 'Jwalamukhi' which can be translated as 'Mouth of 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 many before it.
For instance, it was quite fascinating to understand that volcanoes do not explode all of a sudden, generally. Rather, to put it in my own way - volcanoes are like 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).
It is the latter technique (InSAR) which we will use 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.
Image: Mt. Nyiragongo eruption in Congo. Source: Captured from video by Raphael Kaliwavyo Raks-Brun
Mapping on-surface deformation phenomenon is relatively easier. The recent and 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 observe and plan before the eventual devastating eruption on 15th Jan 2022. (video depiction below_
__________________________________________________________________________
The processing chain involved in extracting deformation (phase fringes) and precise displacement values (up-swell or 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 to derive and visualize deformation and displacement map of Cumbre Vieja Eruption (Canary Islands) 2021. Do read on in case you do not wish to see the video.
__________________________________________________________________________
To explain the above steps in the processing chain individually in a simple manner-
a) Read & Read(2)- Two sets of imagery are used- one prior to eruption and second one post eruption to do change detection.
Image: You may observe that there are 27 Imagery Sections (3*9) in this SAR Imagery set over Cumbre Vieja on 16th September 2021 (pre-eruption)
__________________________________________________________________________
b) TOPSAR-Split- SAR Imagery is heavy (Each set is around 1 GB in size). To reduce computational load while performing analysis, 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.
Image: TopSAR Splitted image (contains 2 sections only now) over Fogo, Cape Verde on 10th October 2014 (pre-eruption)
__________________________________________________________________________
c) Apply (Updated) Orbit Files- Orbit Files within SAR Imagery contains satellite velocity and position information. This may not be the updated / corrected version. 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 of our Region of Interest. This step also allows us to combine the two imagery sets into a single stack allowing for pixel to pixel comparison.
Image: Back-Geocoded Stack of imagery - contains both pre (19 May 2021) and post eruption (12 June 2021) imagery within over Mt. Nyiragongo. The former is the master image while the latter has been designated as the slave.
__________________________________________________________________________
e) Coregistration - Enhanced Spectral Diversity - This optimization tool helps improve the precision of our geometrically derived Coregistered Stack thereby minimizing pixel overlay errors.
Image: ESD Coregistered image of the output from d) Backgeocoding
__________________________________________________________________________
领英推荐
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 centimetre level accuracy.
Image: The Interferogram output over Fogo contains three outputs - a) intensity (uncalibrated energy which returns to the sensor) b) phase (distance between the satellite antenna and the ground targets) and c) coherence (similarity between the base and overlay pixels). The fringes in the phase diagram are indicative of deformation.
__________________________________________________________________________
g) TOPSAR Deburst - Involves removing the black lines (data discontinuities) between the imagery sections (bursts) thereby making the entire imagery seamless.
Image: Notice that the black line (burst) has been removed from the bottom 'debursted' image of the phase over Fogo..
__________________________________________________________________________
h) TopoPhase Removal - There are two types of phase captured in the Interferogram - a) due to topography and b) due to displacement. This step allows us to subtract a) from the total so that we are only left with displacement phase data values.
Image: Notice that after Topo Phase removal the phase diagram significantly changes in the image on the bottom. This is the more accurate representation of phase.
__________________________________________________________________________
i) Multilooking - This step reduces the inherent speckle noise in the imagery (think of Analog TVs grey screen) and makes the pixels square allowing for easy interpretability.
Image: The phase fringes (deformation) are easy to visualize after Multilooking has been applied, allowing for cleaner, square pixels over Mt. Nyiragongo.
__________________________________________________________________________
j) Goldstein Phase Filtering - The next step is to apply Goldstein Phase Filtering. Like Multilooking, this tool is also applied to reduce 'noise' in the image and make the data ready for the later steps which involve the calculation of precise displacement (deformation) values.
Image: Notice that this phase depiction is even more vivid than the one prior (multilooked) over Mt. Nyiragongo.
__________________________________________________________________________
k) Snaphu Export - Snaphu (Phase) Unwrapping - Snaphu Import - Phase to Displacement - These four steps can be clubbed into one (Lets call the combined step as Phase to Displacement only) for the purpose of our understanding. While the interferometric phase gives us an indication where land subsidence (shrinking) and up-swell (growth) has happened, phase unwrapping using Snaphu tool allows us to calculate precise displacement values between post and pre sets of imagery.
Image: While the displacement depiction is not easy on the eyes, notice that the min and max displacement values are computed on the chart to your bottom left.
__________________________________________________________________________
l) (Range-Doppler) Terrain Correction - Alongside eliminating Non-Land displacement values, 'Terrain Correction' also helps in making the imagery geometrically correct - the tilt of satellite, among other aspects, have inbuilt distortions - so that the end product looks as close as possible to the real world.
Image: Terrain corrected displacement values over Fogo volcano. The blue area has displacement of about 10 cms. This is the cone of the volcano which erupted.
__________________________________________________________________________
m) Write (& Export) - the output for further analysis and visualization purposes
Phase Visualization on Google Earth
Video above: 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 above: 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 above: 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).
Hope you found the process of mapping deformation from SAR imagery to study volcanic activity to be interesting!
__________________________________________________________________________
Intelloc Mapping Services | Mapmyops is engaged in selling products which capture geodata (Drones & Drone Services), process geodata (Geographic Information System) and enhance geodata (Imagery & Location Analytics). Together, these help organizations to benefit from Geo-Intelligence for purposes such as operations improvement and project management.
Write to us on [email protected].
Regards, Arpit