M. E. Pritchard

A satellite geodetic survey of large-scale deformation of volcanic centres in the central Andes

Surface deformation in volcanic areas usually indicates movement of magma or hydrothermal fluids at depth. Stratovolcanoes tend to exhibit a complex relationship between deformation and eruptive behaviour. The characteristically long time spans between such eruptions requires a long time series of observations to determine whether deformation without an eruption is common at a given edifice. Such studies, however, are logistically difficult to carry out in most volcanic arcs, as these tend to be remote regions with large numbers of volcanoes (hundreds to even thousands). Here we present a satellite-based interferometric synthetic aperture radar (InSAR) survey of the remote central Andes volcanic arc, a region formed by subduction of the Nazca oceanic plate beneath continental South America. Spanning the years 1992 to 2000, our survey reveals the background level of activity of about 900 volcanoes, 50 of which have been classified as potentially active. We find four centres of broad (tens of kilometres wide), roughly axisymmetric surface deformation. None of these centres are at volcanoes currently classified as potentially active, although two lie within about 10 km of volcanoes with known activity. Source depths inferred from the patterns of deformation lie between 5 and 17 km. In contrast to the four new sources found, we do not observe any deformation associated with recent eruptions of Lascar, Chile.

An aseismic slip pulse in northern Chile and along‐strike variations in seismogenic behavior

We use interferometric synthetic aperture radar, GPS, and seismic observations spanning 5 to 18 years to reveal a detailed kinematic picture of the spatiotemporal evolution of fault slip in a region corresponding to the 30 July 1995 M_w 8.1 subduction zone megathrust earthquake in northern Chile. In a single area, we document a complex mosaic of phenomena including large earthquakes, postseismic afterslip with a spatial distribution that appears to be tied to variations in coastal morphology, and a completely aseismic pulse that may have triggered a M_w 7.1 earthquake on 30 January 1998. In contrast to simple models of fault slip behavior, this spatial heterogeneity indicates that frictional parameters on the fault do not have a systematic transition with depth and also vary rapidly along strike. The low amount of afterslip from the M_w 8.1 earthquake relative to other similar events suggests that postseismic behavior may be modulated by the amount of sediment subducted.

Uturuncu volcano, Bolivia: Volcanic unrest due to mid-crustal magma intrusion

Uturuncu volcano, SW Bolivia, is a dormant stratovolcano (∼85 km3) dominated by dacitic lava domes and flows. 39Ar/40Ar ages show that the volcano was active between 890 ka and 271 ka, with the lavas becoming younger and less extensive at higher elevations. There are current signs of unrest. Between 1992 and 2006 geodetic satellite measurements record an ongoing 70 km deformation field with a central uplift rate of 1 to 2 cm/yr. Deformation indicates volume changes of 400 × 108 m3 over 14 years, an average of ∼1 m3/s (10−2 km3/yr). The deformation is attributed to magma intrusion into the Altiplano-Puna regional crustal magma body. Deformation models indicate a source at depths of 17 to 30 km beneath current local relief. In a reconnaissance survey, persistent seismic activity (mean of 2.6 earthquakes per hour with a maximum of 14 per hour) was recorded at about 4 km depth below the center of the uplift, 4 km SW of the volcanos summit. The seismic events have a normal b value (∼1.04) and activity is attributed to brittle deformation in the elastic crust above the active deep magma intrusion. The porphyritic dacite lavas (64−68% SiO2) have a plagioclase-orthopyroxene-biotite-magnetite-ilmenite assemblage and commonly contain juvenile silicic andesite inclusions, cognate norite nodules and crustal xenoliths. Temperature estimates are in the range 805 to 872°C for the dacites and about 980°C for the silicic andesites. The dacite magmas formed by fractional crystallization of andesite forming norite cumulates and involving partial melting of crust. Compositions and zoning patterns of orthopyroxene and plagioclase phenocrysts indicate that compositional variation in the dacites is caused by magma mixing with the silicic andesite. Reversely zoned orthopyroxene phenocrysts in the andesitic end-member are explained by changing oxidation states during crystallization. Fe3+/Fe2+ ratios from orthopyroxene crystals and Fe3+ in plagioclase provide evidence for a relatively reduced melt that subsequently ascended, degassed and became more oxidized as a consequence of degassing. The geophysical and petrological observations suggest that dacite magma is being intruded into the Altiplano-Puna regional crustal magma body at 17 km or more depth, consistent with deformation models. In the Late Pleistocene dacitic and andesitic magmas ascended from the regional crustal magma body to a shallow magma system at a few kilometers depth where they crystallized and mingled together. The current unrest, together with geophysical anomalies and 270 ka of dormancy, indicate that the magmatic system is in a prolonged period of intrusion. Such circumstances might eventually lead to eruption of large volumes of intruded magma with potential for caldera formation.

Geodetic, teleseismic, and strong motion constraints on slip from recent southern Peru subduction zone earthquakes

We use seismic and geodetic data both jointly and separately to constrain coseismic slip from the 12 November 1996 M_w 7.7 and 23 June 2001 M_w 8.5 southern Peru subduction zone earthquakes, as well as two large aftershocks following the 2001 earthquake on 26 June and 7 July 2001. We use all available data in our inversions: GPS, interferometric synthetic aperture radar (InSAR) from the ERS-1, ERS-2, JERS, and RADARSAT-1 satellites, and seismic data from teleseismic and strong motion stations. Our two-dimensional slip models derived from only teleseismic body waves from South American subduction zone earthquakes with M_w > 7.5 do not reliably predict available geodetic data. In particular, we find significant differences in the distribution of slip for the 2001 earthquake from models that use only seismic (teleseismic and two strong motion stations) or geodetic (InSAR and GPS) data. The differences might be related to postseismic deformation or, more likely, the different sensitivities of the teleseismic and geodetic data to coseismic rupture properties. The earthquakes studied here follow the pattern of earthquake directivity along the coast of western South America, north of 5°S, earthquakes rupture to the north; south of about 12°S, directivity is southerly; and in between, earthquakes are bilateral. The predicted deformation at the Arequipa GPS station from the seismic-only slip model for the 7 July 2001 aftershock is not consistent with significant preseismic motion.

Global link between deformation and volcanic eruption quantified by satellite imagery

A key challenge for volcanological science and hazard management is that few of the world’s volcanoes are effectively monitored. Satellite imagery covers volcanoes globally throughout their eruptive cycles, independent of ground-based monitoring, providing a multidecadal archive suitable for probabilistic analysis linking deformation with eruption. Here we show that, of the 198 volcanoes systematically observed for the past 18 years, 54 deformed, of which 25 also erupted. For assessing eruption potential, this high proportion of deforming volcanoes that also erupted (46%), together with the proportion of non-deforming volcanoes that did not erupt (94%), jointly represent indicators with ‘strong’ evidential worth. Using a larger catalogue of 540 volcanoes observed for 3 years, we demonstrate how this eruption–deformation relationship is influenced by tectonic, petrological and volcanic factors. Satellite technology is rapidly evolving and routine monitoring of the deformation status of all volcanoes from space is anticipated, meaning probabilistic approaches will increasingly inform hazard decisions and strategic development.

Evidence of Recent Thrust Faulting on the Moon Revealed by the Lunar Reconnaissance Orbiter Camera

Lunar Lobate Scarps Revealed Lunar lobate scarps are relatively small-scale landforms that are thought to be formed by tectonic thrust faulting. Previously, lunar lobate scarps could only be identified clearly in high-resolution Apollo Panoramic Camera images confined to the lunar equatorial zone. Now, an analysis by Watters et al. (p. 936) of images returned by the Lunar Reconnaissance Orbiter Camera reveals 14 previously unknown lobate scarps and shows that lunar lobate scarps may be globally distributed. Their appearance suggests that lunar scarps are relatively young landforms (less than 1 Ga), possibly formed during a recent episode of global lunar radial contraction. The relatively young age of the faults and their distribution suggest global, late-stage contraction of the Moon. Lunar Reconnaissance Orbiter Camera images reveal previously undetected lobate thrust-fault scarps and associated meter-scale secondary tectonic landforms that include narrow extensional troughs or graben, splay faults, and multiple low-relief terraces. Lobate scarps are among the youngest landforms on the Moon, based on their generally crisp appearance, lack of superposed large-diameter impact craters, and the existence of crosscut small-diameter impact craters. Identification of previously known scarps was limited to high-resolution Apollo Panoramic Camera images confined to the equatorial zone. Fourteen lobate scarps were identified, seven of which are at latitudes greater than ±60°, indicating that the thrust faults are globally distributed. This detection, coupled with the very young apparent age of the faults, suggests global late-stage contraction of the Moon.

Distribution of slip from 11 Mw > 6 earthquakes in the northern Chile subduction zone

We use interferometric synthetic aperture radar, GPS, and teleseismic data to constrain the relative location of coseismic slip from 11 earthquakes on the subduction interface in northern Chile (23°–25°S) between the years 1993 and 2000. We invert body wave waveforms and geodetic data both jointly and separately for the four largest earthquakes during this time period (1993 M_w 6.8; 1995 M_w 8.1; 1996 M_w 6.7; 1998 M_w 7.1). While the location of slip in the teleseismic-only, geodetic-only, and joint slip inversions is similar for the small earthquakes, there are differences for the 1995 M_w 8.1 event, probably related to nonuniqueness of models that fit the teleseismic data. There is a consistent mislocation of the Harvard centroid moment tensor locations of many of the 6 6 earthquakes, as well as three M_w > 7 events from the 1980s. All of these earthquakes appear to rupture different portions of the fault interface and do not rerupture a limited number of asperities.

Surface cracks record long-term seismic segmentation of the Andean margin

Understanding the long-term patterns of great earthquake rupture along a subduction zone provides a framework for assessing modern seismic hazard. However, evidence that can be used to infer the size and location of past earthquakes is typically erased by erosion after a few thousand years. Meter-scale cracks that cut the surface of coastal areas in northern Chile and southern Peru preserve a record of earthquakes spanning several hundred thousand years owing to the hyperarid climate of the region. These cracks have been observed to form during and/or shortly after strong subduction earthquakes, are preserved for long time periods throughout the Atacama Desert, demonstrate evidence for multiple episodes of reactivation, and show changes in orientation over spatial scales similar to the size of earthquake segments. Our observations and models show that crack orientations are consistent with dynamic and static stress fi elds generated by recent earthquakes. While localized structural and topographic processes infl uence some cracks, the strong preferred orientation over large regions indicates that cracks are primarily formed by plate boundary‐scale stresses, namely repeated earthquakes. We invert the crack-based strain data for slip along the well-known Iquique seismic gap segment of the margin and fi nd consistency with gravity anomaly‐based inferences of long-term earthquake slip patterns, as well as the magnitude and location of the November 2007 Tocopilla earthquake. We suggest that the meter-scale cracks can be used to map characteristic earthquake rupture segments that persist over many seismic cycles, which encourages future study of cracks and other small-scale structures to better constrain the persistence of asperities in other arid, tectonically active regions.

Accounting for Atmospheric Delays in InSAR Data in a Search for Long-Wavelength Deformation in South America

InSAR has been successfully used to observe the deformation of the Earths surface from many processes, but mostly dealing with relatively large signals (>;1 cm) over short wavelengths (<; 100 km). We use interferometric synthetic aperture radar (InSAR) data from two orbital tracks in northern Chile to study the feasibility of imaging the broad interseismic ground deformation signal from the Nazca Plate subduction. In order to measure ~1.5 cm/year of ground motion across ~1000 km of satellite track length due to interseismic loading of the subduction interface, the atmospheric contribution cannot be ignored. We attempt to remove the atmospheric signal using global weather models and by estimating atmospheric parameters directly from the InSAR data. Due to the poor temporal and spatial resolutions of the weather model, this method fails to produce reliable results. The empirical model reduces the phase variance in the interferograms but leaves a residual signal that continues to mask the interseismic signal, which demonstrates the importance of carefully applying corrections to the data as they can significantly affect any interpretation that is based on the corrected observations. Although the methods presented here are not suited for removing all atmospheric path delays, this paper does provide suggestions about improvements that can be made to corrective techniques. Methods that should be further developed are the following: 1) corrections with independent and direct observations of atmospheric properties, e.g., continuous GPS or satellite observations (e.g., the MERIS sensor); 2) improvements using empirical corrections, either in conjunction with a deformation model or constrained by real atmospheric structures; and 3) further work with high resolution and improved weather models.

Magnitude and duration of surface uplift above the Socorro magma body

We use interferometric synthetic aperture radar and geomorphic data to constrain the magnitude and duration of uplift driven by the magma body beneath Socorro, New Mexico, United States. Interferometry spanning 1992–2006 confirms uplift of the Socorro magma body at a rate of ~2.5 mm/yr. However, we find no clear evidence for volcanic uplift after an examination of three rivers (Rio Salado, Rio Puerco, Rio Grande), and two terraces (Llano de Manzano, Llano de Albuquerque) crossing the Socorro magma body. Our geomorphic measurements permit at most 25–50 m of cumulative surface uplift above the Socorro magma body since the middle Pleistocene, but require no long-term uplift. Given previously articulated thermal arguments for the Socorro magma body, we therefore suggest either a recent (within the last few centuries) initiation of uplift at Socorro, or that long-term uplift and subsidence have been essentially equal.