Geology Reports⌕ Search

Geology topics

Zhong Lu

Publications and source records attributed to Zhong Lu.

At least 73 records · Page 4Linked to original sources
USGS2007

Interferometric synthetic-aperature radar (InSAR): Chapter 5

Geodesists are, for the most part, a patient and hardworking lot. A day spent hiking to a distant peak, hours spent waiting for clouds to clear a line-of-sight between observation points, weeks spent moving methodically along a level line – such is the normal pulse of the geodetic profession. The fruits of such labors are all the more precious because they are so scarce. A good day spent with an electronic distance meter (EDM) or level typically produces fewer than a dozen data points. A year of tiltmeter output sampled at ten-minute intervals constitutes less than half a megabyte of data. All of the leveling data ever collected at Yellowstone Caldera fit comfortably on a single PC diskette. These quantities are trivial by modern data-storage standards, in spite of the considerable efforts expended to produce them.

Book chapter↗
USGS2006

Satellite radar interferometry measures deformation at Okmok Volcano

The center of the Okmok caldera in Alaska subsided 140 cm as a result of its February– April 1997 eruption, according to satellite data from ERS-1 and ERS-2 synthetic aperture radar (SAR) interferometry. The inferred deflationary source was located 2.7 km beneath the approximate center of the caldera using a point source deflation model. Researchers believe this source is a magma chamber about 5 km from the eruptive source vent. During the 3 years before the eruption, the center of the caldera uplifted by about 23 cm, which researchers believe was a pre-emptive inflation of the magma chamber. Scientists say such measurements demonstrate that radar interferometry is a promising spaceborne technique for monitoring remote volcanoes. Frequent, routine acquisition of images with SAR interferometry could make near realtime monitoring at such volcanoes the rule, aiding in eruption forecasting.

Eos, Transactions, American Geophysical Union↗
USGS2006

Inflation model of Uzon caldera, Kamchatka, constrained by satellite radar interferometry observations

We analyzed RADARSAT-1 synthetic aperture radar (SAR) data to compute interferometric SAR (InSAR) images of surface deformation at Uzon caldera, Kamchatka, Russia. From 2000 to 2003 approximately 0.15 m of inflation occurred at Uzon caldera, extending beneath adjacent Kikhpinych volcano. This contrasts with InSAR data showing no significant deformation during either the 1999 to 2000, or 2003 to 2004, time periods. We performed three sets of numerical source inversions to fit InSAR data from three different swaths spanning 2000 to 2003. The preferred source model is an irregularly shaped, pressurized crack, dipping ∼20° to the NW, 4 km below the surface. The geometry of this solution is similar to the upper boundary of the geologically inferred magma chamber. Extension of the surface deformation and source to adjacent Kikhpinych volcano, without an eruption, suggests that the deformation is more likely of hydrothermal origin, possibly driven by recharge of the magma chamber.

Kamchatka↗
USGS2006

Thickness distribution of a cooling pyroclastic flow deposit on Augustine Volcano, Alaska: Optimization using InSAR, FEMs, and an adaptive mesh algorithm

Interferometric synthetic aperture radar (InSAR) imagery documents the consistent subsidence, during the interval 1992–1999, of a pyroclastic flow deposit (PFD) emplaced during the 1986 eruption of Augustine Volcano, Alaska. We construct finite element models (FEMs) that simulate thermoelastic contraction of the PFD to account for the observed subsidence. Three-dimensional problem domains of the FEMs include a thermoelastic PFD embedded in an elastic substrate. The thickness of the PFD is initially determined from the difference between post- and pre-eruption digital elevation models (DEMs). The initial excess temperature of the PFD at the time of deposition, 640 °C, is estimated from FEM predictions and an InSAR image via standard least-squares inverse methods. Although the FEM predicts the major features of the observed transient deformation, systematic prediction errors (RMSE = 2.2 cm) are most likely associated with errors in the a priori PFD thickness distribution estimated from the DEM differences. We combine an InSAR image, FEMs, and an adaptive mesh algorithm to iteratively optimize the geometry of the PFD with respect to a minimized misfit between the predicted thermoelastic deformation and observed deformation. Prediction errors from an FEM, which includes an optimized PFD geometry and the initial excess PFD temperature estimated from the least-squares analysis, are sub-millimeter (RMSE = 0.3 mm). The average thickness (9.3 m), maximum thickness (126 m), and volume (2.1 × 10 7 m 3 ) of the PFD, estimated using the adaptive mesh algorithm, are about twice as large as the respective estimations for the a priori PFD geometry. Sensitivity analyses suggest unrealistic PFD thickness distributions are required for initial excess PFD temperatures outside of the range 500–800 °C.

Alaska↗
USGS2006

Quiescent deformation of the Aniakchak Caldera, Alaska mapped by InSAR

The 10-km-wide caldera of the historically active Aniakchak volcano, Alaska, subsides ∼13 mm/yr, based on data from 19 European Remote Sensing Satellite (ERS-1 and ERS-2) interferometric synthetic aperture radar (InSAR) images from 1992 through 2002. The pattern of subsidence does not reflect the distribution of pyroclastic deposits from the last eruption in 1931 and therefore is not related to compaction of fragmental debris. Weighted least-squares inversion of the deformation maps indicates a relatively constant subsidence rate. Modeling the deformation with a Mogi point source locates the source of subsidence at ∼4 km below the central caldera floor, which is consistent with the inferred depth of magma storage before the 1931 eruption. Magmatic CO 2 and He have been measured at a warm soda spring within the caldera, and several sub-boiling fumaroles persist elsewhere in the caldera. These observations suggest that recent subsidence can be explained by the cooling or degassing of a shallow magma body (∼4 km deep), and/or the reduction of the pore-fluid pressure of a cooling hydrothermal system. Ongoing deformation of the volcano detected by InSAR, in combination with magmatic gas output from at least one warm spring, and infrequent low-level bursts of seismicity below the caldera, indicate that the volcanic system is still active and requires close attention for the timely detection of possible hazards.

Alaska↗
USGS2006

Interferometric synthetic aperture radar: Building tomorrow's tools today

A synthetic aperture radar (SAR) system transmits electromagnetic (EM) waves at a wavelength that can range from a few millimeters to tens of centimeters. The radar wave propagates through the atmosphere and interacts with the Earth’s surface. Part of the energy is reflected back to the SAR system and recorded. Using a sophisticated image processing technique, called SAR processing (Curlander and McDonough, 1991), both the intensity and phase of the reflected (or backscattered) signal of each ground resolution element (a few meters to tens of meters) can be calculated in the form of a complex-valued SAR image representing the reflectivity of the ground surface. The amplitude or intensity of the SAR image is determined primarily by terrain slope, surface roughness, and dielectric constants, whereas the phase of the SAR image is determined primarily by the distance between the satellite antenna and the ground targets, slowing of the signal by the atmosphere, and the interaction of EM waves with ground surface. Interferometric SAR (InSAR) imaging, a recently developed remote sensing technique, utilizes the interaction of EM waves, referred to as interference, to measure precise distances. Very simply, InSAR involves the use of two or more SAR images of the same area to extract landscape topography and its deformation patterns.

Alaska Satellite Facility News and Notes↗
USGS2006

Interferometric synthetic aperture radar (InSAR) and its applications to study volcanoes, part 1: Principles of InSAR

Interferometric synthetic aperture radar is an ability to measure the surface deformation of remote sensing technology, in a huge area, its deformation measurement with sub-centimeter accuracy, and spatial resolution in the tens of meters or less. In this paper, the basic theory of InSAR technology is reviewed, its working principle is clarified, and the related problems of surface deformation measurement using InSAR technology are discussed.

Science of Surveying and Mapping↗
USGS2006

Interferometric synthetic aperture radar (InSAR) and its applications to study volcanoes, part 2: InSAR imaging of Alaskan Volcanoes

Interferometric synthetic aperture radar (InSAR) is a remote sensing technique which can measure ground surface deformation with sub-centimeter precision and spatial resolution in tens-of-meters over a large region. This paper summarizes our recent InSAR studies of Alaskan volcanoes, associated with both eruptive and non-eruptive activity. It shows that InSAR can improve our understanding of how the Alaskan volcanoes work and enhance our capability to predict future eruptions and the associated hazards.

Science of Surveying and Mapping↗
USGS2005

InSAR studies of Alaska volcanoes

Interferometric synthetic aperture radar (InSAR) is a remote sensing technique capable of measuring ground surface deformation with sub-centimeter precision and spatial resolution in tens-of­meters over a large region. This paper describes basics of InSAR and highlights our studies of Alaskan volcanoes with InSAR images acquired from European ERS-l and ERS-2, Canadian Radarsat-l, and Japanese JERS-l satellites.

Alaska↗
USGS2004

Toward mapping surface deformation in three dimensions using InSAR

One of the limitations of deformation measurements made with interferometric synthetic aperture radar (InSAR) is that an interferogram only measures one component of the surface deformation — in the satellite's line of sight. We investigate strategies for mapping surface deformation in three dimensions by using multiple interferograms, with different imaging geometries. Geometries for both current and future missions are evaluated, and their abilities to resolve the displacement vector are compared. The north component is always the most difficult to determine using data from near-polar orbiting satellites. However, a satellite with an inclination of about 60°/120° would enable all three components to be well resolved. We attempt to resolve the 3D displacements for the 23 October 2002 Nenana Mountain (Alaska) Earthquake. The north component's error is much larger than the signal, but proxies for eastward and vertical motion can be determined if the north component is assumed negligible. Inversions of hypothetical coseismic interferograms demonstrate that earthquake model parameters can be well recovered from two interferograms, acquired on ascending and descending tracks.

Geophysical Research Letters↗
USGS2004

Transient volcano deformation sources imaged with interferometric synthetic aperture radar: Application to Seguam Island, Alaska

Thirty interferometric synthetic aperture radar (InSAR) images, spanning various intervals during 1992–2000, document coeruptive and posteruptive deformation of the 1992–1993 eruption on Seguam Island, Alaska. A procedure that combines standard damped least squares inverse methods and collective surfaces, identifies three dominant amorphous clusters of deformation point sources. Predictions generated from these three point source clusters account for both the spatial and temporal complexity of the deformation patterns of the InSAR data. Regularized time series of source strength attribute a distinctive transient behavior to each of the three source clusters. A model that combines magma influx, thermoelastic relaxation, poroelastic effects, and petrologic data accounts for the transient, interrelated behavior of the source clusters and the observed deformation. Basaltic magma pulses, which flow into a storage chamber residing in the lower crust, drive this deformational system. A portion of a magma pulse is injected into the upper crust and remains in storage during both coeruption and posteruption intervals. This injected magma degasses and the volatile products accumulate in a shallow poroelastic storage chamber. During the eruption, another portion of the magma pulse is transported directly to the surface via a conduit roughly centered beneath Pyre Peak on the west side of the island. A small amount of this magma remains in storage during the eruption, and posteruption thermoelastic contraction ensues. This model, made possible by the excellent spatial and temporal coverage of the InSAR data, reveals a relatively simple system of interrelated predictable processes driven by magma dynamics.

Journal of Geophysical Research B: Solid Earth↗
USGS2003

Studies of Alaskan volcanoes using synthetic aperature radar and Landsat imagery

Approximately 10 percent of the world’s active volcanoes are located in the Alaskan Aleutian arc and produce about 3-4 explosive eruptions per year. Even with this high amount of volcanic activity, the remote locations and harsh environments of the Aleutian volcanoes conspire to keep them among some of the most poorly studied volcanoes in the world. Space-borne remote sensed imagery can play a significant role in improving our understanding of activity at these volcanoes. Synthetic aperture radar (SAR), Landsat imagery, and Digital Elevation Models (DEMs) derived from SRTM and the National Elevation Database (NED) are used to study several Alaskan volcanoes. Interferometric SAR (InSAR) techniques with ERS-1 and ERS-2 SAR imagery are used to measure ground-surface deformation, which enables the construction of detailed mechanical models that enhance the study of magmatic and tectonic processes. The 30-year historical archive of Landsat data is used to study land cover change, visualize the ash plumes of Aleutian volcanic eruptions, and to map the extent of lava flows. Differencing two DEMs that represent volcano topography before and after an eruption makes it possible to calculate the volume of extruded materials. This paper provides a progress report on how InSAR, Landsat imagery and digital elevation data can be used to better understand the volcanic processes at three Aleutian volcanoes.

Alaska↗
USGS2003

Deriving rangeland structural attributes using Landsat ETM+, ERS-1/ERS-2

The purpose of this study is to determine if Synthetic Aperture Radar (SAR) can be used independently, or in conjunction with Landsat Enhanced Thematic Mapper Plus (ETM+) to improve the classification accuracy of structural attributes of rangeland vegetation, particularly percent shrub cover and top shrub canopy height. Such information, if mapped accurately, can be used in models to better characterize fuel conditions and fire regimes, as well as to evaluate fire hazard status, called for by the U.S. National Fire Plan. The input datasets utilized in this investigation included eighteen bands of Landsat ETM+ path 38 / row 32 (three image dates, six bands each), backscattering and interferometic data derived from tandem ERS-1/2 SAR image pairs (C-band), and extensive field point data. The results showed the use of SAR data provided no significant improvement over the ETM+ data for estimating percent cover or shrub canopy height. The lack of improvement in classification accuracy is possibly due to the influence of topography on the radar backscattering signal. Additional results demonstrated improved model accuracies when a 3x3-averaging filter was applied to the eighteen bands of ETM+ imagery.

Conference Paper↗
USGS2002

Magmatic inflation at a dormant stratovolcano: 1996-1998 activity at Mount Peulik volcano, Alaska, revealed by satellite radar interferometry

A series of ERS radar interferograms that collectively span the time interval from July 1992 to August 2000 reveal that a presumed magma body located 6.6 ??? 0.5 km beneath the southwest flank of the Mount Peulik volcano inflated 0.051 ??? 0.005 km3 between October 1996 and September 1998. Peulik has been active only twice during historical time, in 1814 and 1852, and the volcano was otherwise quiescent during the 1990s. The inflation episode spanned at least several months because separate interferograms show that the associated ground deformation was progressive. The average inflation rate of the magma body was ???0.003 km3/month from October 1996 to September 1997, peaked at 0.005 km3/month from 26 June to 9 October 1997, and dropped to ???0.001 km3/month from October 1997 to September 1998. An intense earthquake swarm, including three ML 4.8 - 5.2 events, began on 8 May 1998 near Becharof Lake, ???30 km northwest of Peulik. More than 400 earthquakes with a cumulative moment of 7.15 ?? 1017 N m were recorded in the area through 19 October 1998. Although the inflation and earthquake swarm occurred at about the same time, the static stress changes that we calculated in the epicentral area due to inflation beneath Peulik appear too small to provide a causal link. The 1996-1998 inflation episode at Peulik confirms that satellite radar interferometry can be used to detect magma accumulation beneath dormant volcanoes at least several months before other signs of unrest are apparent. This application represents a first step toward understanding the eruption cycle at Peulik and other stratovolcanoes with characteristically long repose periods.

Alaska↗
USGS2002

Magmatic activity beneath the quiescent Three Sisters volcanic center, central Oregon Cascade Range, USA

Images from satellite interferometric synthetic aperture radar (InSAR) reveal uplift of a broad ~10 km by 20 km area in the Three Sisters volcanic center of the central Oregon Cascade Range, ~130 km south of Mt. St. Helens. The last eruption in the volcanic center occurred ~1500 years ago. Multiple satellite images from 1992 through 2000 indicate that most if not all of ~100 mm of observed uplift occurred between September 1998 and October 2000. Geochemical (water chemistry) anomalies, first noted during 1990, coincide with the area of uplift and suggest the existence of a crustal magma reservoir prior to the uplift. We interpret the uplift as inflation caused by an ongoing episode of magma intrusion at a depth of ~6.5 km.

Geophysical Research Letters↗
USGS2001

Elastic rebound following the Kocaeli earthquake, Turkey, recorded using synthetic aperture radar interferometry

A basic model incorporating satellite synthetic aperture radar (SAR) interferometry of the fault rupture zone that formed during the Kocaeli earthquake of August 17, 1999, documents the elastic rebound that resulted from the concomitant elastic strain release along the North Anatolian fault. For pure strike-slip faults, the elastic rebound function derived from SAR interferometry is directly invertible from the distribution of elastic strain on the fault at criticality, just before the critical shear stress was exceeded and the fault ruptured. The Kocaeli earthquake, which was accompanied by as much as ∼5 m of surface displacement, distributed strain ∼110 km around the fault prior to faulting, although most of it was concentrated in a narrower and asymmetric 10-km-wide zone on either side of the fault. The use of SAR interferometry to document the distribution of elastic strain at the critical condition for faulting is clearly a valuable tool, both for scientific investigation and for the effective management of earthquake hazard.

Geology↗
USGS2000

Synthetic aperture radar interferometry of Okmok volcano, Alaska: radar observations

ERS-1/ERS-2 synthetic aperture radar interferometry was used to study the 1997 eruption of Okmok volcano in Alaska. First, we derived an accurate digital elevation model (DEM) using a tandem ERS-1/ERS-2 image pair and the preexisting DEM. Second, by studying changes in interferometric coherence we found that the newly erupted lava lost radar coherence for 5-17 months after the eruption. This suggests changes in the surface backscattering characteristics and was probably related to cooling and compaction processes. Third, the atmospheric delay anomalies in the deformation interferograms were quantitatively assessed. Atmospheric delay anomalies in some of the interferograms were significant and consistently smaller than one to two fringes in magnitude. For this reason, repeat observations are important to confidently interpret small geophysical signals related to volcanic activities. Finally, using two-pass differential interferometry, we analyzed the preemptive inflation, coeruptive deflation, and posteruptive inflation and confirmed the observations using independent image pairs. We observed more than 140 cm of subsidence associated with the 1997 eruption. This subsidence occurred between 16 months before the eruption and 5 months after the eruption, was preceded by ∼18 cm of uplift between 1992 and 1995 centered in the same location, and was followed by ∼10 cm of uplift between September 1997 and 1998. The best fitting model suggests the magma reservoir resided at 2.7 km depth beneath the center of the caldera, which was ∼5 km from the eruptive vent. We estimated the volume of the erupted material to be 0.055 km 3 and the average thickness of the erupted lava to be ∼7.4 m. Copyright 2000 by the American Geophysical Union.

Journal of Geophysical Research B: Solid Earth↗