Geology Reports⌕ Search

Geology topics

Zhong Lu

Publications and source records attributed to Zhong Lu.

At least 19 recordsLinked to original sources

Satellite interferometry landslide detection and preliminary tsunamigenic plausibility assessment in Prince William Sound, southcentral Alaska

Regional mapping of actively deforming landslides, including measurements of landslide velocity, is integral for hazard assessments in paraglacial environments. These inventories are also critical for describing the potential impacts that the warming effects of climate change have on slope instability in mountainous and cryospheric terrain. The objective of this study is to identify slow-moving landslides in the Prince William Sound region, southcentral Alaska, United States, which has had rapid deglaciation since the mid-1800s, and assess their tsunamigenic plausibility. We use an automated time series persistent scatterer interferometric synthetic aperture radar processing method with 7 years of Sentinel-1 data (2016–22) to identify 43 slow-moving slopes with average velocities ranging from approximately 0.2 to 21 millimeters per year. Landslide presence is confirmed using aerial imagery and previous landslide inventory records. We assess the tsunamigenic plausibility of the landslides using empirically derived estimates of landslide mobility based on modeled landslide volumes. Of the identified landslides, our preliminary analysis suggests that 11 have tsunamigenic potential if they were to fail rapidly and catastrophically. Although our estimate of tsunamigenic plausibility is preliminary and can be refined with additional observations and analyses, it can be used to prioritize ongoing and future hazard assessment, surveillance, and research efforts.

Alaska↗

Spaceborne InSAR mapping of landslides and subsidence in rapidly deglaciating terrain, Glacier Bay National Park and Preserve and vicinity, Alaska and British Columbia

The Glacier Bay area in southeastern Alaska and British Columbia, encompassing Glacier Bay National Park and Preserve, has experienced rapid glacier retreat since the end of the Little Ice Age in the mid-1800s. The impact that rapid deglaciation has had on the slope stability of valley walls and on the sedimentation of fans and deltas adjacent to fjords and inlets is an ongoing research topic. Using 3-year (2018–2020) Sentinel-1 datasets, and an automated time-series persistent scatterer interferometric synthetic aperture radar (PSInSAR) processing method, we detected landslides or delta subsidence at 27 sites within a vast 180 × 180 km remote region encompassing Glacier Bay proper. Most of the sites that we identified had not been previously identified. We categorized the hazard source areas that we identified into three general types:1) slow-moving landslides on steep rocky slopes not near (> 2 km away from) present-day glacier termini (e.g., Tidal Inlet landslide), 2) slow-moving landslides directly adjacent to (< 2 km away from), and associated with glacier thinning and retreat, and 3) subsidence of glacial outwash fan deltas. In categories 1 and 2, we observed 22 landslides moving at velocities ranging from 0.5 to 4 cm/yr. In category 3, we detected five fan deltas subsiding at velocities ranging from 0.5 to 6 cm/yr. Within our measurement error, these velocities were consistent during the monitoring period. Because acceleration was not observed, the issuance of warnings of imminent rapid failure is not warranted, however, continued remote monitoring is warranted. Our interferometric synthetic aperture radar (InSAR) results could be combined with other data sets including field observations, subaerial and submarine landslide inventories, bedrock fabric mapping from newly available light detection and ranging (lidar) data, and geologic maps to produce an inherent susceptibility map for landslides in bedrock and fan deltas. This map could be used to forecast susceptibility for both earthquake and climatically induced landslides.

Alaska, British Columbia↗

Geologic controls of slow-moving landslides near the U.S. West Coast

Slow-moving landslides, often with nearly imperceptible creeping motion, are an important landscape shaper and a dangerous natural hazard across the globe, yet their spatial distribution and geologic controls are still poorly known owing to a paucity of detailed, large-area observations. Here, we use interferometry of L-band satellite radar images to reveal 617 spatially large (4 &#x00D7; "> × 10 4 – 13 &#x00D7; "> × 10 6 m 2 ) and presently active (2007 – 2019) slow-moving landslides near the populous US West Coast (only 4.6% of these slides were previously known) and provide evidence for their fundamental controls by bedrock lithology and vertical land motion. We found that slow-moving landslides are generally larger and more spatially frequent in homogeneous bedrock with low rock strength, and they are preferentially located on hillslopes with geologically recent uplift. Notably, landslide size and spatial density in the relatively weak metamorphic rocks and mélange (due to pervasive tectonically sheared discontinuities, foliation, and abundant clay minerals) were two times larger than those in sedimentary and igneous rocks, and the hillslopes with landslides were found to be uplifting approximately three times faster than the average for the whole region. These results suggest that slow-moving landslides can be effectively uncovered by satellite radar imagery and their occurrence and character may be anticipated from vertical land uplift and bedrock lithology. Hence, our study provides understanding critical for reducing landslide hazards and quantifying landslide impacts on landscape change.

California, Oregon, Washington↗

Quantifying eruptive and background seismicity, deformation, degassing, and thermal emissions at volcanoes in the United States during 1978–2020

An important aspect of volcanic hazard assessment is determination of the level and character of background activity at a volcano so that deviations from background (called unrest) can be identified. Here, we compile the instrumentally recorded eruptive and noneruptive activity for 161 US volcanoes between 1978 and 2020. We combine monitoring data from four techniques: seismicity, ground deformation, degassing, and thermal emissions. To previous work, we add the first comprehensive survey of US volcanoes using medium-spatial resolution satellite thermal observations, newly available field surveys of degassing, and new compilations of seismic and deformation data. We report previously undocumented thermal activity at 30 volcanoes using data from the spaceborne ASTER sensor during 2000–2020. To facilitate comparison of activity levels for all US volcanoes, we assign a numerical classification of the Activity Intensity Level for each monitoring technique, with the highest ranking corresponding to an eruption. There are 96 US volcanoes (59%) with at least one type of detected activity, but this represents a lower bound: For example, there are 12 volcanoes where degassing has been observed but has not yet been quantified. We identify dozens of volcanoes where volcanic activity is only measured by satellite (45% of all thermal observations), and other volcanoes where only ground-based sensors have detected activity (e.g., all seismic and 62% of measured degassing observations). Our compilation provides a baseline against which future measurements can be compared, demonstrates the need for both ground-based and remote observations, and serves as a guide for prioritizing future monitoring efforts.

Journal of Geophysical Research↗

Landslide monitoring and runout hazard assessment by integrating multi-source remote sensing and numerical models: An application to the Gold Basin landslide complex, northern Washington

The landslide complex at Gold Basin, Washington, has been drawing considerable attention after a catastrophic runout of the nearby landslide at Oso, Washington, in 2014. To evaluate potential threats of the Gold Basin landslide to the campground down the slope, remote sensing and numerical modeling were integrated to monitor recent landslide activity and simulate hypothetical runout scenarios. Bare-earth LiDAR DEM (digital elevation model) differencing, InSAR (Interferometric Synthetic Aperture Radar), and offset tracking of SAR images reveal that localized collapses at the headscarps have been the primary type of landslide activity at Gold Basin from 2005 to 2019, and currently no signs indicative of movement of a large centralized block or a deep-seated main body were detected. The maximum horizontal deformation rate is 5 m/year occurring primarily from headscarp recession of the middle lobe, and the annual landsliding volume of the whole landslide complex averages 1.03 × 10 5 m 3 . From three-dimensional limit equilibrium analysis of generalized terrace structures, the maximum landslide volume is estimated as 2.0 × 10 6 m 3 . Simulations of hypothetical runout scenarios were carried out using the depth-averaged two-phase model D-claw with above-obtained landslide geometry constraints. The simulation results demonstrate that debris flows with volume less than 10 5 m 3 only pose limited threats to the campground, while volumes over 10 6 m 3 could cause severe damages. Consequently, the estimated maximum landslide volume of 2.0 × 10 6 m 3 suggests a potential risk to the campground nearby. Adaption of our methodology could prove useful for evaluating other similar landslides globally for hazards prevention and mitigation.

Washington↗

High rates of inflation during a noneruptive episode of seismic unrest at Semisopochnoi Volcano, Alaska in 2014–2015

Magma intrusion rate is a key parameter in eruption triggering but is poorly quantified in existing geodetic studies. Here we examine two episodes of rapid inflation in this context. Two noneruptive microseismic swarms were recorded at Semisopochnoi Volcano, Alaska in 2014–2015. We use differential SAR techniques and TerraSAR‐X images to document surface deformation from 2011 to 2015, which comprises island‐wide radial inflation totaling ~25 cm (+/−1 cm) line of sight displacement in 2014–2015. Multiple source geometries are tested in an inversion of the deformation data, and InSAR data are best fit by a spheroid trending to the northeast and plunging to the southeast, with a major axis of ~4 km and minor axes of ~1 km, directly under the central caldera of Semisopochnoi. In 2014, a modeled influx of 0.043 km 3 of magma caused line of sight displacement of ~17 cm. This magma was stored at a depth of ~8 km, until 2015 when 0.029 km 3 was added. Along with the definition of inflation source parameters, the recorded seismic events are relocated using differential travel times. These relocated events outline a linear aseismic area within a larger zone of shallow (<10 km) seismicity. This aseismic region aligns with the centroid of the deformation model. Based on these geodetic and seismic models, the plumbing system at Semisopochnoi is interpreted as a spheroidal magma storage zone at a depth of ˜8 km below a linear feature of partial melt. The observed deformation and seismicity appear to result from rapid injection into this main storage region.

Alaska↗

Space-based imaging radar studies of U.S. volcanoes

The arrival of space-based imaging radar as a revolutionary land-surface mapping and monitoring tool little more than a quarter century ago enabled a spate of innovative volcano research worldwide. Soon after launch of European Space Agency’s ERS-1 spacecraft in 1991, the U.S. Geological Survey began SAR and InSAR studies of volcanoes in the Aleutian and Cascades arcs, in Hawai’i, and elsewhere in the western U.S. including the Yellowstone and Long Valley calderas. This paper summarizes results of that effort and presents new findings concerning: (1) prevalence of volcano deformation in the Aleutian and Cascade arcs; (2) surface-change detection and hazard assessment during eruptions at Aleutian and Hawaiian volcanoes; (3) geodetic imaging of magma storage and transport systems in Hawai’i; and (4) deformation sources and processes at the Yellowstone and Long Valley calderas. Surface deformation caused by a variety of processes is common in arc settings and could easily escape detection without systematic InSAR surveillance. Space-based SAR imaging of active lava flows and domes in remote or heavily vegetated settings, including during periods of bad weather and darkness, extends land-based monitoring capabilities and improves hazards assessments. At Kīlauea Volcano, comprehensive SAR and InSAR observations identify multiple magma storage zones beneath the summit area and along the East Rift Zone, and illuminate magma transport pathways. The same approach at Yellowstone tracks the ascent of magmatic volatiles from a mid-crustal intrusion to shallow depth and relates that process to increased hydrothermal activity at the surface. Together with recent and planned launches of highly capable imaging-radar satellites, these findings support an optimistic outlook for near-real time surveillance of volcanoes at global scale in the coming decade.

Frontiers in Earth Science↗

Interferometric synthetic aperture radar study of recent eruptive activity at Shrub mud volcano, Alaska

Shrub mud volcano is one of three large mud volcanoes that comprise the Klawasi Group in the Copper River Basin of southcentral Alaska. Except for minor discharges in the mid-1950s when the group was first described, Shrub was dormant prior to its reactivation in summer 1996. From 1997 to 1999, Shrub vigorously erupted more than 5 x 105 cubic meters of saline mud and carbon dioxide-rich gas at temperatures as high as 54 degrees C. Thereafter, activity waned but continued at least through 2015. We analyzed 192 interferograms derived from 106 synthetic aperture radar (SAR) images acquired by the JERS-1 (L-band), ERS-1/2 (C-band), RADARSAT-1 (C-band), and ALOS PALSAR (L-band) satellites to characterize ground deformation at Shrub before, during, and after its reactivation. Collectively, the interferograms span 1992–2000 and 2006–2011. We fit the observations with two deformation sources: a deflating, steeply dipping, pipe-like body under the summit area and an inflating, shallow-dipping, sill-like body under the southwest flank. Both sources are shallow, with centroids less than 1 km beneath the summit. Prior to reactivation, the flank source inflated ~0.35 x 105 cubic meters per year from July 1992 to May 1996. During eruptive activity, the summit source deflated at higher rates that peaked at ~8.71 x 105 cubic meters per year during May–November 1997 and continued at ~0.95 x 105 cubic meters per year during the 2006–2011 observation window. Cumulative source-volume loss is comparable to the volume of mud erupted. We interpret the summit source as the volcano’s feeder conduit that pressurized prior to the first SAR observation in 1992. Also before 1992, the conduit ruptured to feed a lateral intrusion of mud under the southwest flank, perhaps along a bedding plane in underlying glaciolacustrine deposits. The growing sill caused the southwest flank to inflate while it accommodated the mud supply from depth, which explains why we observed pre-eruptive inflation of the flank but not the summit. The summit began deflating when the conduit ruptured to the surface at the onset of eruptive activity. The flank source did not deflate concurrently because the weight of the thin overburden was insufficient to collapse the sill. There is a suggestion in the modern topography that lateral intrusions under Shrub’s southwest flank are a common feature of activity there.

Alaska↗

Modeling wildfire-induced permafrost deformation in an Alaskan boreal forest using InSAR observations

The discontinuous permafrost zone is one of the world’s most sensitive areas to climate change. Alaskan boreal forest is underlain by discontinuous permafrost, and wildfires are one of the most influential agents negatively impacting the condition of permafrost in the arctic region. Using interferometric synthetic aperture radar (InSAR) of Advanced Land Observation Satellite (ALOS) Phased Array type L-band Synthetic Aperture Radar (PALSAR) images, we mapped extensive permafrost degradation over interior Alaskan boreal forest in Yukon Flats, induced by the 2009 Big Creek wildfire. Our analyses showed that fire-induced permafrost degradation in the second post-fire thawing season contributed up to 20 cm of ground surface subsidence. We generated post-fire deformation time series and introduced a model that exploited the deformation time series to estimate fire-induced permafrost degradation and changes in active layer thickness. The model showed a wildfire-induced increase of up to 80 cm in active layer thickness in the second post-fire year due to pore-ice permafrost thawing. The model also showed up to 15 cm of permafrost degradation due to excess-ice thawing with little or no increase in active layer thickness. The uncertainties of the estimated change in active layer thickness and the thickness of thawed excess ice permafrost are 27.77 and 1.50 cm, respectively. Our results demonstrate that InSAR-derived deformation measurements along with physics models are capable of quantifying fire-induced permafrost degradation in Alaskan boreal forests underlain by discontinuous permafrost. Our results also have illustrated that fire-induced increase of active layer thickness and excess ice thawing contributed to ground surface subsidence.

Alaska↗

Combining InSAR and GPS to determine transient movement and thickness of a seasonally active low-gradient translational landslide

The combined application of continuous Global Positioning System data (high temporal resolution) with spaceborne interferometric synthetic aperture radar data (high spatial resolution) can reveal much more about the complexity of large landslide movement than is possible with geodetic measurements tied to only a few specific measurement sites. This approach is applied to an ~4 km 2 reactivated translational landslide in the Columbia River Gorge (Washington State), which moves mainly during the winter rainy season. Results reveal the complex three-dimensional shape of the landslide mass, how onset of sliding relates to cumulative rainfall, how surface velocity during sliding varies with location on the topographically complex landslide surface, and how the ground surface subsides slightly in weeks prior to downslope sliding.

Geophysical Research Letters↗

Rapid late Pleistocene/Holocene uplift and coastal evolution of the southern Arabian (Persian) Gulf

The coastline along the southern Arabian Gulf between Al Jubail, Kingdom of Saudi Arabia, and Dubai, UAE, appears to have risen at least 125 m in the last 18,000 years. Dating and topographic surveying of paleo-dunes (43–53 ka), paleo-marine terraces (17–30 ka), and paleo-marine shorelines (3.3–5.5 ka) document a rapid, > 1 mm/a subsidence, followed by a 6 mm/a uplift that is decreasing with time. The mechanism causing this movement remains elusive but may be related to the translation of the coastal area through the backbasin to forebulge hinge line movement of the Arabian plate or, alternatively, by movement of the underlying Infracambrian-age Hormuz salt in response to sea-level changes associated with continental glaciation. Independent of the mechanism, rapid and episodic uplift may impact the design of engineering projects such as nuclear power plants, airports, and artificial islands as well as the interpretation of sedimentation and archeology of the area.

Quaternary Research↗

Detecting seasonal landslide movement within the Cascade landslide complex (Washington) using time-series SAR imagery

Detection of slow or limited landslide movement within broad areas of forested terrain has long been problematic, particularly for the Cascade landslide complex (Washington) located along the Columbia River Gorge. Although parts of the landslide complex have been found reactivated in recent years, the timing and magnitude of motion have not been systematically monitored or interpreted. Here we apply novel time-series strategies to study the spatial distribution and temporal behavior of the landslide movement between 2007 and 2011 using InSAR images from two overlapping L-band ALOS PALSAR-1 satellite tracks. Our results show that the reactivated part has moved approximately 700 mm downslope during the 4-year observation period, while other parts of the landslide complex have generally remained stable. However, we also detect about 300 mm of seasonal downslope creep in a terrain block upslope of the Cascade landslide complex—terrain previously thought to be stable. The temporal oscillation of the seasonal movement can be correlated with precipitation, implying that seasonal movement here is hydrology-driven. The seasonal movement also has a frequency similar to GPS-derived regional ground oscillations due to mass loading by stored rainfall and subsequent rebound but with much smaller magnitude, suggesting different hydrological loading effects. From the time-series amplitude information on terrain upslope of the headscarp, we also re-evaluate the incipient motion related to the 2008 Greenleaf Basin rock avalanche, not previously recognized by traditional SAR/InSAR methods. The approach used in this study can be used to identify active landslides in forested terrain, to track the seasonal movement of landslides, and to identify previously unknown landslide hazards.

Washington↗

Dome growth at Mount Cleveland, Aleutian Arc, quantified by time-series TerraSAR-X imagery

Synthetic aperture radar imagery is widely used to study surface deformation induced by volcanic activity; however, it is rarely applied to quantify the evolution of lava domes, which is important for understanding hazards and magmatic system characteristics. We studied dome formation associated with eruptive activity at Mount Cleveland, Aleutian Volcanic Arc, in 2011&ndash;2012 using TerraSAR-X imagery. Interferometry and offset tracking show no consistent deformation and only motion of the crater rim, suggesting that ascending magma may pass through a preexisting conduit system without causing appreciable surface deformation. Amplitude imagery has proven useful for quantifying rates of vertical and areal growth of the lava dome within the crater from formation to removal by explosive activity to rebirth. We expect that this approach can be applied at other volcanoes that host growing lava domes and where hazards are highly dependent on dome geometry and growth rates.

Alaska↗

Post-eruptive inflation of Okmok Volcano, Alaska, from InSAR, 2008–2014

Okmok, a ~10-km wide caldera that occupies most of the northeastern end of Umnak Island, is one of the most active volcanoes in the Aleutian arc. The most recent eruption at Okmok during July-August 2008 was by far its largest and most explosive since at least the early 19th century. We investigate post-eruptive magma supply and storage at the volcano during 2008&ndash;2014 by analyzing all available synthetic aperture radar (SAR) images of Okmok acquired during that time period using the multi-temporal InSAR technique. Data from the C-band Envisat and X-band TerraSAR-X satellites indicate that Okmok started inflating very soon after the end of 2008 eruption at a time-variable rate of 48-130 mm/y, consistent with GPS measurements. The &ldquo;model-assisted&rdquo; phase unwrapping method is applied to improve the phase unwrapping operation for long temporal baseline pairs. The InSAR time-series is used as input for deformation source modeling, which suggests magma accumulating at variable rates in a shallow storage zone at ~3.9 km below sea level beneath the summit caldera, consistent with previous studies. The modeled volume accumulation in the 6 years following the 2008 eruption is ~75% of the 1997 eruption volume and ~25% of the 2008 eruption volume.

Alaska↗

Measurement of slow-moving along-track displacement from an efficient multiple-aperture SAR interferometry (MAI) stacking

Multiple-aperture SAR interferometry (MAI) has demonstrated outstanding measurement accuracy of along-track displacement when compared to pixel-offset-tracking methods; however, measuring slow-moving (cm/year) surface displacement remains a challenge. Stacking of multi-temporal observations is a potential approach to reducing noise and increasing measurement accuracy, but it is difficult to achieve a significant improvement by applying traditional stacking methods to multi-temporal MAI interferograms. This paper proposes an efficient MAI stacking method, where multi-temporal forward- and backward-looking residual interferograms are individually stacked before the MAI interferogram is generated. We tested the performance of this method using ENVISAT data from Kīlauea Volcano, Hawai&lsquo;i, where displacement on the order of several centimeters per year is common. By comparing results from the proposed stacking methods with displacements from GPS data, we documented measurement accuracies of about 1.03 and 1.07 cm/year for the descending and ascending tracks, respectively&mdash;an improvement of about a factor of two when compared with that from the conventional stacking approach. Three-dimensional surface-displacement maps can be constructed by combining stacked InSAR and MAI observations, which will contribute to a better understanding of a variety of geological phenomena.

Hawaii↗

The 2008 phreatomagmatic eruption of Okmok volcano, Aleutian Islands, Alaska: Chronology, deposits, and landform changes

Okmok volcano, Aleutian Islands, Alaska, explosively erupted over a five-week period between July 12 and August 23, 2008. The eruption was predominantly phreatomagmatic, producing fine-grained tephra that covered most of northeastern Umnak Island. The eruption had a maximum Volcanic Explosivity Index (VEI) of 4, with eruption column heights up to 16 km during the opening phase. Several craters and a master tuff cone formed in the caldera as a result of phreatomagmatic explosions and accumulated tephra-fall and surge deposits. Ascending magma continuously interacted with an extensive shallow groundwater table in the caldera, resulting in the phreatomagmatic character of the eruption. Syneruptive explosion and collapse processes enlarged a pre-existing lake, created a second, entirely new lake, and formed new, deep craters. A field of ephemeral collapse pits and collapse escarpments formed where rapid groundwater withdrawal removed material from beneath capping lava flows. This was the first significant phreatomagmatic event in the U.S. since the Ukinrek Maars eruption in 1977.

Alaska↗

InSAR detects possible thaw settlement in the Alaskan Arctic Coastal Plain

Satellite interferometric synthetic aperture radar (InSAR) has proven to be an effective tool for monitoring surface deformation from volcanoes, earthquakes, landslides, and groundwater withdrawal. This paper seeks to expand the list of applications of InSAR data to include monitoring subsidence possibly associated with thaw settlement over the Alaskan Arctic Coastal Plain. To test our hypothesis that InSAR data are sufficiently sensitive to detect subsidence associated with thaw settlement, we acquired all Japanese Earth Resources Satellite-1 (JERS-1) L-band data available for the summers of 1996, 1997, and 1998 over two sites on the Alaska North Slope. The least amount of subsidence for both study sites was detected in the interferograms covering the summer of 1996 (2-3 cm), interferograms from 1997 and 1998 revealed that about 3 cm of subsidence occurred at the northern Cache One Lake site, and about 5 cm of subsidence was detected at the southern Kaparuk River site. These preliminary results illustrate the capacity of the L-band (24 cm) wavelength JERS-1 radar data to penetrate the short Arctic vegetation to monitor subsidence possibly associated with thaw settlement of the active layer and (or) other hydrologic changes over relatively large areas.

Canadian Journal of Remote Sensing↗