Geology topics

Z. Lu

Publications and source records attributed to Z. Lu.

34 records · Page 2Linked to original sources

USGS2006

On the absence of InSAR-detected volcano deformation spanning the 1995-1996 and 1999 eruptions of Shishaldin Volcano, Alaska

Shishaldin Volcano, a large, frequently active basaltic-andesite volcano located on Unimak Island in the Aleutian Arc of Alaska, had a minor eruption in 1995–1996 and a VEI 3 sub-Plinian basaltic eruption in 1999. We used 21 synthetic aperture radar images acquired by ERS-1, ERS-2, JERS-1, and RADARSAT-1 satellites to construct 12 coherent interferograms that span most of the 1993–2003 time interval. All interferograms lack coherence within ∼5 km of the summit, primarily due to persistent snow and ice cover on the edifice. Remarkably, in the 5–15 km distance range where interferograms are coherent, the InSAR images show no intrusion- or withdrawal-related deformation at Shishaldin during this entire time period. However, several InSAR images do show deformation associated with a shallow M L 5.2 earthquake located ∼14 km west of Shishaldin that occurred 6 weeks before the 1999 eruption. We use a theoretical model to predict deformation magnitudes due to a volumetric expansion source having a volume equivalent to the 1999 erupted volume, and find that deformation magnitudes for sources shallower than 10 km are within the expected detection capabilities for interferograms generated from C-band ERS 1/2 and RADARSAT-1 synthetic aperture radar images. We also find that InSAR images cannot resolve relatively shallow deformation sources (1–2 km below sea level) due to spatial gaps in the InSAR images caused by lost coherence. The lack of any deformation, particularly for the 1999 eruption, leads us to speculate that magma feeding eruptions at the summit moves rapidly (at least 80m/day) from > 10 km depth, and that the intrusion–eruption cycle at Shishaldin does not produce significant permanent deformation at the surface.

Alaska↗

USGS2005

Interferometric synthetic aperture radar study of Okmok volcano, Alaska, 1992-2003: Magma supply dynamics and postemplacement lava flow deformation

Okmok volcano, located in the central Aleutian arc, Alaska, is a dominantly basaltic complex topped with a 10-km-wide caldera that formed circa 2.05 ka. Okmok erupted several times during the 20th century, most recently in 1997; eruptions in 1945, 1958, and 1997 produced lava flows within the caldera. We used 80 interferometric synthetic aperture radar (InSAR) images (interferograms) to study transient deformation of the volcano before, during, and after the 1997 eruption. Point source models suggest that a magma reservoir at a depth of 3.2 km below sea level, located beneath the center of the caldera and about 5 km northeast of the 1997 vent, is responsible for observed volcano-wide deformation. The preeruption uplift rate decreased from about 10 cm yr −1 during 1992–1993 to 2 ∼ 3 cm yr −1 during 1993–1995 and then to about −1 ∼ −2 cm yr −1 during 1995–1996. The posteruption inflation rate generally decreased with time during 1997–2001, but increased significantly during 2001–2003. By the summer of 2003, 30 ∼ 60% of the magma volume lost from the reservoir in the 1997 eruption had been replenished. Interferograms for periods before the 1997 eruption indicate consistent subsidence of the surface of the 1958 lava flows, most likely due to thermal contraction. Interferograms for periods after the eruption suggest at least four distinct deformation processes: (1) volcano-wide inflation due to replenishment of the shallow magma reservoir, (2) subsidence of the 1997 lava flows, most likely due to thermal contraction, (3) deformation of the 1958 lava flows due to loading by the 1997 flows, and (4) continuing subsidence of 1958 lava flows buried beneath 1997 flows. Our results provide insights into the postemplacement behavior of lava flows and have cautionary implications for the interpretation of inflation patterns at active volcanoes.

Alaska↗

USGS2005

DEM, tide and velocity over sulzberger ice shelf, West Antarctica

Arctic and Antarctic ice sheets preserve more than 77% of the global fresh water and could raise global sea level by several meters if completely melted. Ocean tides near and under ice shelves shifts the grounding line position significantly and are one of current limitations to study glacier dynamics and mass balance. The Sulzberger ice shelf is an area of ice mass flux change in West Antarctica and has not yet been well studied. In this study, we use repeat-pass synthetic aperture radar (SAR) interferometry data from the ERS-1 and ERS-2 tandem missions for generation of a high-resolution (60-m) Digital Elevation Model (DEM) including tidal deformation detection and ice stream velocity of the Sulzberger Ice Shelf. Other satellite data such as laser altimeter measurements with fine foot-prints (70-m) from NASA's ICESat are used for validation and analyses. The resulting DEM has an accuracy of-0.57??5.88 m and is demonstrated to be useful for grounding line detection and ice mass balance studies. The deformation observed by InSAR is found to be primarily due to ocean tides and atmospheric pressure. The 2-D ice stream velocities computed agree qualitatively with previous methods on part of the Ice Shelf from passive microwave remote-sensing data (i.e., LANDSAT). ?? 2005 IEEE.

Conference Paper↗

USGS2005

Digital elevation model of King Edward VII Peninsula, West Antarctica, from SAR interferometry and ICESat laser altimetry

We present a digital elevation model (DEM) of King Edward VII Peninsula, Sulzberger Bay, West Antarctica, developed using 12 European Remote Sensing (ERS) synthetic aperture radar (SAR) scenes and 24 Ice, Cloud, and land Elevation Satellite (ICESat) laser altimetry profiles. We employ differential interferograms from the ERS tandem mission SAR scenes acquired in the austral fall of 1996, and four selected ICESat laser altimetry profiles acquired in the austral fall of 2004, as ground control points (GCPs) to construct an improved geocentric 60-m resolution DEM over the grounded ice region. We then extend the DEM to include two ice shelves using ICESat profiles via Kriging. Twenty additional ICESat profiles acquired in 2003-2004 are used to assess the accuracy of the DEM. After accounting for radar penetration depth and predicted surface changes, including effects due to ice mass balance, solid Earth tides, and glacial isostatic adjustment, in part to account for the eight-year data acquisition discrepancy, the resulting difference between the DEM and ICESat profiles is -0.57 ?? 5.88 m. After removing the discrepancy between the DEM and ICESat profiles for a final combined DEM using a bicubic spline, the overall difference is 0.05 ?? 1.35 m. ?? 2005 IEEE.

IEEE Geoscience and Remote Sensing Letters↗

USGS2004

Mapping recent lava flows at Westdahl Volcano, Alaska, using radar and optical satellite imagery

Field mapping of young lava flows at Aleutian volcanoes is logistically difficult, and the utility of optical images from aircraft or satellites for this purpose is greatly reduced by persistent cloud cover. These factors have hampered earlier estimates of the areas and volumes of three young lava flows at Westdahl Volcano, including its most recent (1991–1992) flow. We combined information from synthetic aperture radar (SAR) images with multispectral Landsat-7 data to differentiate the 1991–1992 flow from the 1964 flow and a pre-1964 flow, and to calculate the flow areas (8.4, 9.2, and 7.3 km 2 , respectively). By differencing a digital elevation model (DEM) from the 1970–1980s with a DEM from the Shuttle Radar Topography Mission (SRTM) in February 2000, we estimated the average thickness of the 1991–1992 flow to be 13 m, which reasonably agrees with field observations (5–10 m). Lava-flow maps produced in this way can be used to facilitate field mapping and flow-hazards assessment, and to study magma-supply dynamics and thus to anticipate future eruptive activity. Based on the recurrence interval of recent eruptions and the results of this study, the next eruption at Westdahl may occur before the end of this decade.

Alaska↗

USGS2004

Constraining the slip distribution and fault geometry of the M w 7.9, 3 November 2002, Denali fault earthquake with Interferometric Synthetic Aperture Radar and Global Positioning System data

The M w 7.9, Denali fault earthquake ( DFE ) is the largest continental strike-slip earthquake to occur since the development of Interferometric Synthetic Aperture Radar (In SAR ). We use five interferograms, constructed using radar images from the Canadian Radarsat-1 satellite, to map the surface deformation at the western end of the fault rupture. Additional geodetic data are provided by displacements observed at 40 campaign and continuous Global Positioning System ( GPS ) sites. We use the data to determine the geometry of the Susitna Glacier fault, thrusting on which initiated the DFE , and to determine a slip model for the entire event that is consistent with both the In SAR and GPS data. We find there was an average of 7.3 ± 0.4 m slip on the Susitna Glacier fault, between 1 and 9.5 km depth on a 29 km long fault that dips north at 41 ± 0.7° and has a surface projection close to the mapped rupture. On the Denali fault, a simple model with large slip patches finds a maximum of 8.7 ± 0.7 m of slip between the surface and 14.3 ± 0.2 km depth. A more complex distributed slip model finds a peak of 12.5 ± 0.8 m in the upper 4 km, significantly higher than the observed surface slip. We estimate a geodetic moment of 670 ± 10 × 10 18 N m ( M w 7.9), consistent with seismic estimates. Lack of preseismic data resulted in an absence of In SAR coverage for the eastern half of the DFE rupture. A dedicated geodetic In SAR mission could obviate coverage problems in the future.

Bulletin of the Seismological Society of America↗

USGS2004

DEM generation and tidal deformation detection for Sulzberger Ice Shelf, West Antarctica using SAR interferometry

In this study we generated a relative Digital Elevation Model (DEM) over the Sulzberger Ice Shelf, West Antarctica using ERS1/2 synthetic aperture radar (SAR) interferometry data. Four repeat pass differential interferograms are used to find the grounding zone and to classify the study area. An interferometrically derived DEM is compared with laser altimetry profile from ICESat. Standard deviation of the relative height difference is 5.12 m and 1.34 m in total length of the profile and at the center of the profile respectively. The magnitude and the direction of tidal changes estimated from interferogram are compared with those predicted tidal differences from four ocean tide models. Tidal deformation measured in InSAR is -16.7 cm and it agrees well within 3 cm with predicted ones from tide models.

Conference Paper↗

USGS2003

Source model for the Mw 6.7, 23 October 2002, Nenana Mountain earthquake (Alaska) from InSAR

The 23 October 2002 Nenana Mountain Earthquake (M w ∼ 6.7) occurred on the Denali Fault (Alaska), to the west of the M w ∼ 7.9 Denali Earthquake that ruptured the same fault 11 days later. We used 6 interferograms, constructed using radar images from the Canadian Radarsat-1 and European ERS-2 satellites, to determine the coseismic surface deformation and a source model. Data were acquired on ascending and descending satellite passes, with incidence angles between 23 and 45 degrees, and time intervals of 72 days or less. Modeling the event as dislocations in an elastic half space suggests that there was nearly 0.9 m of right-lateral strike-slip motion at depth, on a near-vertical fault, and that the maximum slip in the top 4 km of crust was less than 0.2 m. The Nenana Mountain Earthquake increased the Coulomb stress at the future hypocenter of the 3 November 2002, Denali Earthquake by 30–60 kPa.

Alaska↗

USGS2003

Magma supply dynamics at Westdahl volcano, Alaska, modeled from satellite radar interferometry

A group of satellite radar interferograms that span the time period from 1991 to 2000 shows that Westdahl volcano, Alaska, deflated during its 1991-1992 eruption and is reinflating at a rate that could produce another eruption within the next several years. The rates of inflation and deflation are approximated by exponential decay functions having time constants of about 6 years and a few days, respectively. This behavior is consistent with a deep, constant-pressure magma source connected to a shallow reservoir by a magma-filled conduit. An elastic deformation model indicates that the reservoir is located about 6 km below sea level and beneath Westdahl Peak. We propose that the magma flow rate through the conduit is governed by the pressure gradient between the deep source and the reservoir. The pressure gradient, and hence the flow rate, are greatest immediately after eruptions. Pressurization of the reservoir decreases both the pressure gradient and the flow rate, but eventually the reservoir ruptures and an eruption or intrusion ensues. The eruption rate is controlled partly by the pressure gradient between the reservoir and surface, and therefore it, too, decreases with time. When the supply of eruptible magma is exhausted, the eruption stops, the reservoir begins to repressurize at a high rate, and the cycle repeats. This model might also be appropriate for other frequently active volcanoes with stable magma sources and relatively simple magma storage systems.

Journal of Geophysical Research B: Solid Earth↗

USGS2003

The 1997 eruption of Okmok Volcano, Alaska: A synthesis of remotely sensed imagery

Okmok Volcano, in the eastern Aleutian Islands, erupted in February and March of 1997 producing a 6-km-long lava flow and low-level ash plumes. This caldera is one of the most active in the Aleutian Arc, and is now the focus of international multidisciplinary studies. A synthesis of remotely sensed data (AirSAR, derived DEMs, Landsat MSS and ETM+ data, AVHRR, ERS, JERS, Radarsat) has given a sequence of events for the virtually unobserved 1997 eruption. Elevation data from the AirSAR sensor acquired in October 2000 over Okmok were used to create a 5-m resolution DEM mosaic of Okmok Volcano. AVHRR nighttime imagery has been analyzed between February 13 and April 11, 1997. Landsat imagery and SAR data recorded prior to and after the eruption allowed us to accurately determine the extent of the new flow. The flow was first observed on February 13 without precursory thermal anomalies. At this time, the flow was a large single lobe flowing north. According to AVHRR Band 3 and 4 radiance data and ground observations, the first lobe continued growing until mid to late March, while a second, smaller lobe began to form sometime between March 11 and 12. This is based on a jump in the thermal and volumetric flux determined from the imagery, and the physical size of the thermal anomalies. Total radiance values waned after March 26, indicating lava effusion had ended and a cooling crust was growing. The total area (8.9 km 2 ), thickness (up to 50 m) and volume (1.54×10 8 m 3 ) of the new lava flow were determined by combining observations from SAR, Landsat ETM+, and AirSAR DEM data. While the first lobe of the flow ponded in a pre-eruption depression, our data suggest the second lobe was volume-limited. Remote sensing has become an integral part of the Alaska Volcano Observatory’s monitoring and hazard mitigation efforts. Studies like this allow access to remote volcanoes, and provide methods to monitor potentially dangerous ones.

Journal of Volcanology and Geothermal Research↗

USGS2003

Estimating lava volume by precision combination of multiple baseline spaceborne and airborne interferometric synthetic aperture radar: The 1997 eruption of Okmok Volcano, Alaska

Interferometric synthetic aperture radar (InSAR) techniques are used to calculate the volume of extrusion at Okmok volcano, Alaska by constructing precise digital elevation models (DEMs) that represent volcano topography before and after the 1997 eruption. The posteruption DEM is generated using airborne topographic synthetic aperture radar (TOPSAR) data where a three-dimensional affine transformation is used to account for the misalignments between different DEM patches. The preeruption DEM is produced using repeat-pass European Remote Sensing satellite data; multiple interferograms are combined to reduce errors due to atmospheric variations, and deformation rates are estimated independently and removed from the interferograms used for DEM generation. The extrusive flow volume associated with the 1997 eruption of Okmok volcano is 0.154 ?? 0.025 km3. The thickest portion is approximately 50 m, although field measurements of the flow margin's height do not exceed 20 m. The in situ measurements at lava edges are not representative of the total thickness, and precise DEM data are absolutely essential to calculate eruption volume based on lava thickness estimations. This study is an example that demonstrates how InSAR will play a significant role in studying volcanoes in remote areas.

IEEE Transactions on Geoscience and Remote Sensing↗

USGS2002

Preeruptive inflation and surface interferometric coherence characteristics revealed by satellite radar interferometry at Makushin Volcano, Alaska: 1993-2000

Pilot reports in January 1995 and geologic field observations from the summer of 1996 indicate that a relatively small explosive eruption of Makushin, one of the more frequently active volcanoes in the Aleutian arc of Alaska, occurred on 30 January 1995. Several independent radar interferograms that each span the time period from October 1993 to September 1995 show evidence of ???7 cm of uplift centered on the volcano's east flank, which we interpret as preeruptive inflation of a ???7-km-deep magma source (??V = 0.022 km3). Subsequent interferograms for 1995-2000, a period that included no reported eruptive activity, show no evidence of additional ground deformation. Interferometric coherence at C band is found to persist for 3 years or more on lava flow and other rocky surfaces covered with short grass and sparsely distributed tall grass and for at least 1 year on most pyroclastic deposits. On lava flow and rocky surfaces with dense tall grass and on alluvium, coherence lasts for a few months. Snow and ice surfaces lose coherence within a few days. This extended timeframe of coherence over a variety of surface materials makes C band radar interferometry an effective tool for studying volcano deformation in Alaska and other similar high-latitude regions.

Journal of Geophysical Research B: Solid Earth↗

USGS2002

Interferometric synthetic aperture radar studies of Alaska volcanoes

Interferometric synthetic aperture radar (InSAR) imaging is a recently developed geodetic technique capable of measuring ground-surface deformation with centimeter to subcentimeter vertical precision and spatial resolution of tens-of-meter over a relatively large region (/spl sim/10/sup 4/ km/sup 2/). The spatial distribution of surface deformation data, derived from InSAR images, enables the construction of detailed mechanical models to enhance the study of magmatic and tectonic processes associated with volcanoes. This paper summarizes our recent InSAR studies of several Alaska volcanoes, which include Okmok, Akutan, Kiska, Augustine, Westdahl, and Peulik volcanoes.

Alaska↗

USGS2000

Aseismic inflation of Westdahl volcano, Alaska, revealed by satellite radar interferometry

Westdahl volcano, located at the west end of Unimak Island in the central Aleutian volcanic arc, Alaska, is a broad shield that produced moderate-sized eruptions in 1964, 1978-79, and 1991-92. Satellite radar interferometry detected about 17 cm of volcano-wide inflation from September 1993 to October 1998. Multiple independent interferograms reveal that the deformation rate has not been steady; more inflation occurred from 1993 to 1995 than from 1995 to 1998. Numerical modeling indicates that a source located about 9 km beneath the center of the volcano inflated by about 0.05 km3 from 1993 to 1998. On the basis of the timing and volume of recent eruptions at Westdahl and the fact that it has been inflating for more than 5 years, the next eruption can be expected within the next several years.

Geophysical Research Letters↗

USGS2000

Studies of volcanoes of Alaska by satellite radar interferometry

Interferometric synthetic aperture radar (InSAR) has provided a new imaging geodesy technique to measure the deformation of volcanoes at tens-of-meter horizontal resolution with centimeter to subcentimeter vertical precision. The two-dimensional surface deformation data enables the construction of detailed numerical models allowing the study of magmatic and tectonic processes beneath volcanoes. This paper summarizes our recent: InSAR studies over the Alaska-Aleutian volcanoes, which include New Trident, Okmok, Akutan, Augustine, Shishaldin, and Westdahl volcanoes. The first InSAR surface deformation over the Alaska volcanoes was applied to New Trident. Preliminary InSAR study suggested that New Trident volcano experienced several centimeters inflation from 1993 to 1995. Using the InSAR technique, we studied the 1997 eruption of Okmok. We have measured ???1.4 m deflation during the eruption, ???20 cm pre-eruptive inflation during 1992 to 1995, and >10 cm post-eruptive inflation within a year after the eruption, and modeled the deformations using Mogi sources. We imaged the ground surface deformation associated with the 1996 seismic crisis over Akutan volcano. Although seismic swarm did not result in an eruption, we found that the western part of the volcano uplifted ???60 cm while the eastern part of the island subsided. The majority of the complex deformation field at the Akutan volcano was modeled by dike intrusion and Mogi inflation sources. Our InSAR results also indicate that the pyroclastic flows from last the last eruption have been undergoing contraction/subsidence at a rate of about 3 cm per year since 1992. InSAR measured no surface deformation before and during the 1999 eruption of Shishaldin and suggested the eruption may be a type of open system. Finally, we applied satellite radar interferometry to Westdahl volcano which erupted 1991 and has been quiet since. We discovered this volcano had inflated about 15 cm from 1993 to 1998. In summary, satellite radar interferometry can not only be used to study a volcanic eruption, but also to detect aseismic deformation at quiescent volcanoes preceding a seismic swarm; it is a useful technique to study volcanic eruptions as well as to guide scientists to better focus their monitoring efforts.

Alaska↗

USGS1996

Use of SAR data to study active volcanoes in Alaska

Synthetic Aperture Radar (SAR) data of the Westdahl, Veniaminof, and Novarupta volcanoes in the Aleutian Arc of Alaska were analysed to investigate recent surface volcanic processes. These studies support ongoing monitoring and research by the Alaska Volcano Observatory (AVO) in the North Pacific Ocean Region. Landforms and possible crustal deformation before, during, or after eruptions were detected and analysed using data from the European Remote Sensing Satellites (ERS), the Japanese Earth Resources Satellite (JERS) and the US Seasat platforms. Field observations collected by scientists from the AVO were used to verify the results from the analysis of SAR data.

Earth Observation Quarterly↗