Geology Reports⌕ Search

Geology topics

Zhong Lu

Publications and source records attributed to Zhong Lu.

At least 37 records · Page 2Linked to original sources
USGS2014

InSAR imaging of aleutian volcanoes: Monitoring a volcanic arc from space

Interferometric synthetic aperture radar (InSAR) is a relatively new remote sensing tool that is capable of measuring ground-surface deformation with centimeter-to-subcentimeter precision at a spatial resolution of tens of meters over an area of hundreds to thousands of square kilometers. With its global coverage and all-weather imaging capability, InSAR has become an increasingly important technique for studying volcanoes in remote regions such as the Aleutian Islands. The spatial distribution of surface deformation data derived from InSAR images enables the construction of detailed mechanical models to enhance the study of magmatic processes.

Book↗
USGS2014

Monitoring Everglades freshwater marsh water level using L-band synthetic aperture radar backscatter

The Florida Everglades plays a significant role in controlling floods, improving water quality, supporting ecosystems, and maintaining biodiversity in south Florida. Adaptive restoration and management of the Everglades requires the best information possible regarding wetland hydrology. We developed a new and innovative approach to quantify spatial and temporal variations in wetland water levels within the Everglades, Florida. We observed high correlations between water level measured at in situ gages and L-band SAR backscatter coefficients in the freshwater marsh, though C-band SAR backscatter has no close relationship with water level. Here we illustrate the complementarity of SAR backscatter coefficient differencing and interferometry (InSAR) for improved estimation of high spatial resolution water level variations in the Everglades. This technique has a certain limitation in applying to swamp forests with dense vegetation cover, but we conclude that this new method is promising in future applications to wetland hydrology research.

Florida↗
USGS2013

The utility of atmospheric analyses for the mitigation of artifacts in InSAR

The numerical weather models (NWMs) developed by the meteorological community are able to provide accurate analyses of the current state of the atmosphere in addition to the predictions of the future state. To date, most attempts to apply the NWMs to estimate the refractivity of the atmosphere at the time of satellite synthetic aperture radar (SAR) data acquisitions have relied on predictive models. We test the hypothesis that performing a final assimilative routine, ingesting all available meteorological observations for the times of SAR acquisitions, and generating customized analyses of the atmosphere at those times will better mitigate atmospheric artifacts in differential interferograms. We find that, for our study area around Mount St. Helens (Amboy, Washington, USA), this approach is unable to model the refractive changes and provides no mean benefit for interferogram analysis. The performance is improved slightly by ingesting atmospheric delay estimates derived from the limited local GPS network; however, the addition of water vapor products from the GOES satellites reduces the quality of the corrections. We interpret our results to indicate that, even with this advanced approach, NWMs are not a reliable mitigation technique for regions such as Mount St. Helens with highly variable moisture fields and complex topography and atmospheric dynamics. It is possible, however, that the addition of more spatially dense meteorological data to constrain the analyses might significantly improve the performance of weather modeling of atmospheric artifacts in satellite radar interferograms.

Journal of Geophysical Research B: Solid Earth↗
USGS2013

Pre-eruption deformation caused by dike intrusion beneath Kizimen volcano, Kamchatka, Russia, observed by InSAR

Interferometric synthetic aperture radar (InSAR) images reveal a pre-eruption deformation signal at Kizimen volcano, Kamchatka, Russia, where an ongoing eruption began in mid-November, 2010. The previous eruption of this basaltic andesite-to-dacite stratovolcano occurred in 1927–1928. InSAR images from both ascending and descending orbital passes of Envisat and ALOS PALSAR satellites show as much as 6 cm of line-of-sight shortening from September 2008 to September 2010 in a broad area centered at Kizimen. About 20 cm of opening of a nearly vertical dike provides an adequate fit to the surface deformation pattern. The model dike is approximately 14 km long, 10 km high, centered 13 km beneath Kizimen, and strikes NE–SW. Time-series analysis of multi-temporal interferograms indicates that (1) intrusion started sometime between late 2008 and July 2009, (2) continued at a nearly constant rate, and (3) resulted in a volume expansion of 3.2 × 10 7 m 3 by September 2010, i.e., about two months before the onset of the 2010 eruption. Earthquakes located above the tip of the dike accompanied the intrusion. Eventually, magma pressure in the dike exceeded the confining strength of the host rock, triggering the 2010 eruption. Our results provide insight into the intrusion process that preceded an explosive eruption at a Pacific Rim stratovolcano following nearly a century of quiescence, and therefore have implications for monitoring and hazards assessment at similar volcanoes elsewhere.

Kamchatka, Kizimen Volcano↗
USGS2013

Dynamic deformation of Seguam Island, Alaska, 1992--2008, from multi-interferogram InSAR processing

We generated a time-series of ERS-1/2 and ENVISAT interferometric synthetic aperture radar (InSAR) images to study ground surface deformation at Seguam Island from 1992 to 2008. We used the small baseline subset (SBAS) technique to reduce artifacts associated with baseline uncertainties and atmospheric delay anomalies, and processed images from two adjacent tracks to validate our results. Seguam Island comprises the remnants of two late Quaternary calderas, one in the western caldera of the island and one in the eastern part of the island. The western caldera subsided at a constant rate of ~ 1.6 cm/yr throughout the study period, while the eastern caldera experienced alternating periods of subsidence and uplift: ~ 5 cm/year uplift during January 1993–October 1993 (stage 1), ~ 1.6 cm/year subsidence during October 1993–November 1998 (stage 2), ~ 2.0 cm/year uplift during November 1998–September 2000 (stage 3), ~ 1.4 cm/year subsidence during September 2000–November 2005 (stage 4), and ~ 0.8 cm/year uplift during November 2005– July 2007 (stage 5). Source modeling indicates a deflationary source less than 2 km below sea level (BSL) beneath the western caldera and two sources beneath the eastern caldera: an inflationary source 2.5–6.0 km BSL and a deflationary source less than 2 km BSL. We suggest that uplift of the eastern caldera is driven by episodic intrusions of basaltic magma into a poroelastic reservoir 2.5–6.0 km BSL beneath the caldera. Cooling and degassing of the reservoir between intrusions results in steady subsidence of the overlying surface. Although we found no evidence of magma intrusion beneath the western caldera during the study period, it is the site (Pyre Peak) of all historical eruptions on the island and therefore cooling and degassing of intrusions presumably contributes to subsidence there as well. Another likely subsidence mechanism in the western caldera is thermoelastic contraction of lava flows emplaced near Pyre Peak during several historical eruptions, most recently in 1977 and 1992–93.

Alaska↗
USGS2012

Aftershocks of the 13 May 1993 Shumagin Alaska earthquake

The 13 May 1993 M s 6.9 Shumagin earthquake had an aftershock sequence of 247 earthquakes with magnitudes greater than or equal to 1.5 by 1 June 1993. Of these aftershocks, 79 were located by using S-P travel times at the only two stations within 570 km of the mainshock epicenter. The rupture area inferred from the aftershocks is about 600 km 2 and we estimate for the mainshock a mean fault displacement of 1.0 m and a 28 bar stress drop. The magnitude-frequency plots give a b-value for the aftershock sequence of about 0.4, which is low compared to the background value of approximately 0.8. The decay of the aftershock sequence followed the modified Omori law with a p-value of 0.79, which is also lower than the typical values of about 1.1 observed in Alaska. Both of these facts can be interpreted as indicating relatively high ambient stress in the Shumagin seismic gap and the possibility that the 13 May earthquake was a foreshock to a larger gap-filling event to occur within the next few years.

Geophysical Research Letters↗
USGS2012

Mapping ground surface deformation using temporarily coherent point SAR interferometry: Application to Los Angeles Basin

Multi-temporal interferometric synthetic aperture radar (InSAR) is an effective tool to detect long-term seismotectonic motions by reducing the atmospheric artifacts, thereby providing more precise deformation signal. The commonly used approaches such as persistent scatterer InSAR (PSInSAR) and small baseline subset (SBAS) algorithms need to resolve the phase ambiguities in interferogram stacks either by searching a predefined solution space or by sparse phase unwrapping methods; however the efficiency and the success of phase unwrapping cannot be guaranteed. We present here an alternative approach – temporarily coherent point (TCP) InSAR (TCPInSAR) – to estimate the long term deformation rate without the need of phase unwrapping. The proposed approach has a series of innovations including TCP identification, TCP network and TCP least squares estimator. We apply the proposed method to the Los Angeles Basin in southern California where structurally active faults are believed capable of generating damaging earthquakes. The analysis is based on 55 interferograms from 32 ERS-1/2 images acquired during Oct. 1995 and Dec. 2000. To evaluate the performance of TCPInSAR on a small set of observations, a test with half of interferometric pairs is also performed. The retrieved TCPInSAR measurements have been validated by a comparison with GPS observations from Southern California Integrated GPS Network. Our result presents a similar deformation pattern as shown in past InSAR studies but with a smaller average standard deviation (4.6 mm) compared with GPS observations, indicating that TCPInSAR is a promising alternative for efficiently mapping ground deformation even from a relatively smaller set of interferograms.

Remote Sensing of Environment↗
USGS2012

Large-area landslide detection and monitoring with ALOS/PALSAR imagery data over Northern California and Southern Oregon, USA

Multi-temporal ALOS/PALSAR images are used to automatically investigate landslide activity over an area of ~ 200 km by ~ 350 km in northern California and southern Oregon. Interferometric synthetic aperture radar (InSAR) deformation images, InSAR coherence maps, SAR backscattering intensity images, and a DEM gradient map are combined to detect active landslides by setting individual thresholds. More than 50 active landslides covering a total of about 40 km 2 area are detected. Then the short baseline subsets (SBAS) InSAR method is applied to retrieve time-series deformation patterns of individual detected landslides. Down-slope landslide motions observed from adjacent satellite tracks with slightly different radar look angles are used to verify InSAR results and measurement accuracy. Comparison of the landslide motion with the precipitation record suggests that the landslide deformation correlates with the rainfall rate, with a lag time of around 1–2 months between the precipitation peak and the maximum landslide displacement. The results will provide new insights into landslide mechanisms in the Pacific Northwest, and facilitate development of early warning systems for landslides under abnormal rainfall conditions. Additionally, this method will allow identification of active landslides in broad areas of the Pacific Northwest in an efficient and systematic manner, including remote and heavily vegetated areas difficult to inventory by traditional methods.

California, Oregon↗
USGS2012

Digital elevation model generation from satellite interferometric synthetic aperture radar: Chapter 5

An accurate digital elevation model (DEM) is a critical data set for characterizing the natural landscape, monitoring natural hazards, and georeferencing satellite imagery. The ideal interferometric synthetic aperture radar (InSAR) configuration for DEM production is a single-pass two-antenna system. Repeat-pass single-antenna satellite InSAR imagery, however, also can be used to produce useful DEMs. DEM generation from InSAR is advantageous in remote areas where the photogrammetric approach to DEM generation is hindered by inclement weather conditions. There are many sources of errors in DEM generation from repeat-pass InSAR imagery, for example, inaccurate determination of the InSAR baseline, atmospheric delay anomalies, and possible surface deformation because of tectonic, volcanic, or other sources during the time interval spanned by the images. This chapter presents practical solutions to identify and remove various artifacts in repeat-pass satellite InSAR images to generate a high-quality DEM.

Book chapter↗
USGS2012

Limitations and potential of satellite imagery to monitor environmental response to coastal flooding

Storm-surge flooding and marsh response throughout the coastal wetlands of Louisiana were mapped using several types of remote sensing data collected before and after Hurricanes Gustav and Ike in 2008. These included synthetic aperture radar (SAR) data obtained from the (1) C-band advance SAR (ASAR) aboard the Environmental Satellite, (2) phased-array type L-band SAR (PALSAR) aboard the Advanced Land Observing Satellite, and (3) optical data obtained from Thematic Mapper (TM) sensor aboard the Land Satellite (Landsat). In estuarine marshes, L-band SAR and C-band ASAR provided accurate flood extent information when depths averaged at least 80 cm, but only L-band SAR provided consistent subcanopy detection when depths averaged 50 cm or less. Low performance of inundation mapping based on C-band ASAR was attributed to an apparent inundation detection limit (>30 cm deep) in tall Spartina alterniflora marshes, a possible canopy collapse of shoreline fresh marsh exposed to repeated storm-surge inundations, wind-roughened water surfaces where water levels reached marsh canopy heights, and relatively high backscatter in the near-range portion of the SAR imagery. A TM-based vegetation index of live biomass indicated that the severity of marsh dieback was linked to differences in dominant species. The severest impacts were not necessarily caused by longer inundation but rather could be caused by repeated exposure of the palustrine marsh to elevated salinity floodwaters. Differential impacts occurred in estuarine marshes. The more brackish marshes on average suffered higher impacts than the more saline marshes, particularly the nearshore coastal marshes occupied by S. alterniflora.

Journal of Coastal Research↗
USGS2011

Deformation of the 2002 Denali Fault earthquakes, mapped by Radarsat-1 interferometry

The magnitude 7.9 earthquake that struck central Alaska on 3 November 2002 was the largest strike-slip earthquake in North America for more than 150 years. The earthquake ruptured about 340 km of the Denali Fault system with observed right-lateral offsets of up to 9 m [ Eberhart-Phillips et al. , 2003] (Figure l). The rupture initiated with slip on a previously unknown thrust fault, the 40-km-long Susitna Glacier Fault. The rupture propagated eastward for about 220 km along the right-lateral Denali Fault where right-lateral slip averaged ˜5 m, before stepping southeastward onto the Totschunda Fault for about 70 km, with offsets as large as 2 m. The 3 November earthquake was preceded by a magnitude 6.7 shock on 23 October—the Nenana Mountain Earthquake—which was located about 25 km to the west of the 3 November earthquake.

Eos, Transactions, American Geophysical Union↗
USGS2011

C-band radar observes water level change in swamp forests

C-band radar pulses backscatter from the upper canopy of swamp forests, and consequently interferometric synthetic aperture radar (InSAR) analysis of C-band imagery has not been exploited to study water level changes in swamp forests. This article explores C-band ERS-1 (European Remote Sensing Satellite) and ERS-2 InSAR data over swamp forests composed of moderately dense trees with a medium-low canopy closure in southeastern Louisiana to measure water level changes beneath tree cover. Wetlands cover more than 4% of the Earth's land surface and interact with hydrologic, biogeochemical, and sediment transport processes that are fundamental in understanding ecological and climatic changes [ Alsdorf et al , 2003; Prigent et al ., 2001 ; Melack and Forsberg , 2000; Dunne et al ., 1998]. Measurement of water level changes in wetlands, and consequently of changes in water storage capacity, provides a required input for hydrologic models, and is required to comprehensively assess flood hazards [e.g., Coe , 1998].

Eos, Transactions, American Geophysical Union↗
USGS2011

Volcano deformation and gravity workshop synopsis and outcomes: The 2008 volcano deformation and temporal gravity change workshop

A volcano workshop was held in Washington State, near the U.S. Geological Survey (USGS) Cascades Volcano Observatory. The workshop, hosted by the USGS Volcano Hazards Program (VHP), included more than 40 participants from the United States, the European Union, and Canada. Goals were to promote (1) collaboration among scientists working on active volcanoes and (2) development of new tools for studying volcano deformation. The workshop focused on conventional and emerging techniques, including the Global Positioning System (GPS), borehole strain, interferometric synthetic aperture radar (InSAR), gravity, and electromagnetic imaging, and on the roles of aqueous and magmatic fluids.

Eos, Transactions, American Geophysical Union↗
USGS2011

Monitoring a boreal wildfire using multi-temporal Radarsat-1 intensity and coherence images

Twenty-five C-band Radarsat-1 synthetic aperture radar (SAR) images acquired from the summer of 2002 to the summer of 2005 are used to map a 2003 boreal wildfire (B346) in the Yukon Flats National Wildlife Refuge, Alaska under conditions of near-persistent cloud cover. Our analysis is primarily based on the 15 SAR scenes acquired during arctic growing seasons. The Radarsat-1 intensity data are used to map the onset and progression of the fire, and interferometric coherence images are used to qualify burn severity and monitor post-fire recovery. We base our analysis of the fire on three test sites, two from within the fire and one unburned site. The B346 fire increased backscattered intensity values for the two burn study sites by approximately 5–6 dB and substantially reduced coherence from background levels of approximately 0.8 in unburned background forested areas to approximately 0.2 in the burned area. Using ancillary vegetation information from the National Land Cover Database (NLCD) and information on burn severity from Normalized Burn Ratio (NBR) data, we conclude that burn site 2 was more severely burned than burn site 1 and that C-band interferometric coherence data are useful for mapping landscape changes due to fire. Differences in burn severity and topography are determined to be the likely reasons for the observed differences in post-fire intensity and coherence trends between burn sites.

Geomatics, Natural Hazards and Risk↗
USGS2011

Monitoring of urban subsidence with SAR interferometric point target analysis: A case study in Suzhou, China

Ground subsidence, mainly caused by over exploitation of groundwater and other underground resources, such as oil, gas and coal, occurs in many cities in China. The annual direct loss associated with subsidence across the country is estimated to exceed 100 million US dollar. Interferometric SAR (InSAR) is a powerful tool to map ground deformation at an unprecedented level of spatial detail. It has been widely used to investigate the deformation resulting from earthquakes, volcanoes and subsidence. Repeat-pass InSAR, however, may fail due to impacts of spatial decorrelation, temporal decorrelation and heterogeneous refractivity of atmosphere. In urban areas, a large amount of natural stable radar reflectors exists, such as buildings and engineering structures, at which radar signals can remain coherent during a long time interval. Interferometric point target analysis (IPTA) technique, also known as persistent scatterers (PS) InSAR is based on these reflectors. It overcomes the shortfalls in conventional InSAR. This paper presents a procedure for urban subsidence monitoring with IPTA. Calculation of linear deformation rate and height residual, and the non-linear deformation estimate, respectively, are discussed in detail. Especially, the former is highlighted by a novel and easily implemented 2-dimensional spatial search algorithm. Practically useful solutions that can significantly improve the robustness of IPTA, are recommended. Finally, the proposed procedure is applied to mapping the ground subsidence in Suzhou city, Jiangsu province, China. Thirty-four ERS-1/2 SAR scenes are analyzed, and the deformation information over 38,881 point targets between 1992 and 2000 are generated. The IPTA-derived deformation estimates correspond well with leveling measurements, demonstrating the potential of the proposed subsidence monitoring procedure based on IPTA technique. Two shortcomings of the IPTA-based procedure, e.g., the requirement of large number of SAR images and assumed linear plus non-linear deformation model, are discussed as the topics of further research.

Jaingsu Province↗
USGS2011

Monitoring duration and extent of storm-surge and flooding in Western Coastal Louisiana marshes with Envisat ASAR data

Inundation maps of coastal marshes in western Louisiana were created with multitemporal Envisat Advanced Synthetic Aperture (ASAR) scenes collected before and during the three months after Hurricane Rita landfall in September 2005. Corroborated by inland water-levels, 7 days after landfall, 48% of coastal estuarine and palustrine marshes remained inundated by storm-surge waters. Forty-five days after landfall, storm-surge inundated 20% of those marshes. The end of the storm-surge flooding was marked by an abrupt decrease in water levels following the passage of a storm front and persistent offshore winds. A complementary dramatic decrease in flood extent was confirmed by an ASAR-derived inundation map. In nonimpounded marshes at elevations <;80 cm, storm-surge waters rapidly receded while slower recession was dominantly associated with impounded marshes at elevations >;80 cm during the first month after Rita landfall. After this initial period, drainage from marshes-especially impounded marshes-was hastened by the onset of offshore winds. Following the abrupt drops in inland water levels and flood extent, rainfall events coinciding with increased water levels were recorded as inundation re-expansion. This postsurge flooding decreased until only isolated impounded and palustrine marshes remained inundated. Changing flood extents were correlated to inland water levels and largely occurred within the same marsh regions. Trends related to incremental threshold increases used in the ASAR change-detection analyses seemed related to the preceding hydraulic and hydrologic events, and VV and HH threshold differences supported their relationship to the overall wetland hydraulic condition.

Louisiana↗
USGS2010

Broad accommodation of rift-related extension recorded by dyke intrusion in Saudi Arabia

The extensive harrat lava province of Arabia formed during the past 30 million years in response to Red Sea rifting and mantle upwelling. The area was regarded as seismically quiet, but between April and June 2009 a swarm of more than 30,000 earthquakes struck one of the lava fields in the province, Harrat Lunayyir, northwest Saudi Arabia. Concerned that larger damaging earthquakes might occur, the Saudi Arabian government evacuated 40,000 people from the region. Here we use geologic, geodetic and seismic data to show that the earthquake swarm resulted from magmatic dyke intrusion. We document a surface fault rupture that is 8 km long with 91 cm of offset. Surface deformation is best modelled by the shallow intrusion of a north-west trending dyke that is about 10 km long. Seismic waves generated during the earthquakes exhibit overlapping very low- and high-frequency components. We interpret the low frequencies to represent intrusion of magma and the high frequencies to represent fracturing of the crystalline basement rocks. Rather than extension being accommodated entirely by the central Red Sea rift axis, we suggest that the broad deformation observed in Harrat Lunayyir indicates that rift margins can remain as active sites of extension throughout rifting. Our analyses allowed us to forecast the likelihood of a future eruption or large earthquake in the region and informed the decisions made by the Saudi Arabian government to return the evacuees.

Nature Geoscience↗