Dana J. Caccamise

The 2010 Maule, Chile earthquake: Downdip rupture limit revealed by space geodesy

Radar interferometry from the ALOS satellite captured the coseismic ground deformation associated with the 2010 Mw 8.8 Maule, Chile earthquake. The ALOS interferograms reveal a sharp transition in fringe pattern at ~150 km from the trench axis that is diagnostic of the downdip rupture limit of the Maule earthquake. An elastic dislocation model based on ascending and descending ALOS interferograms and 13 near-field 3-component GPS measurements reveals that the coseismic slip decreases more or less linearly from a maximum of 17 m (along-strike average of 6.5 m) at 18 km depth to near zero at 43–48 km depth, quantitatively indicating the downdip limit of the seismogenic zone. The depth at which slip drops to near zero appears to be at the intersection of the subducting plate with the continental Moho. Our model also suggests that the depth where coseismic slip vanishes is nearly uniform along the strike direction for a rupture length of ~600 km. The average coseismic slip vector and the interseismic velocity vector are not parallel, which can be interpreted as a deficit in strike-slip moment release.

Bedrock displacements in Greenland manifest ice mass variations, climate cycles and climate change

The Greenland GPS Network (GNET) uses the Global Positioning System (GPS) to measure the displacement of bedrock exposed near the margins of the Greenland ice sheet. The entire network is uplifting in response to past and present-day changes in ice mass. Crustal displacement is largely accounted for by an annual oscillation superimposed on a sustained trend. The oscillation is driven by earth’s elastic response to seasonal variations in ice mass and air mass (i.e., atmospheric pressure). Observed vertical velocities are higher and often much higher than predicted rates of postglacial rebound (PGR), implying that uplift is usually dominated by the solid earth’s instantaneous elastic response to contemporary losses in ice mass rather than PGR. Superimposed on longer-term trends, an anomalous ‘pulse’ of uplift accumulated at many GNET stations during an approximate six-month period in 2010. This anomalous uplift is spatially correlated with the 2010 melting day anomaly.

Geodetic measurements of vertical crustal velocity in West Antarctica and the implications for ice mass balance

We present preliminary geodetic estimates for vertical bedrock velocity at twelve survey GPS stations in the West Antarctic GPS Network, an additional survey station in the northern Antarctic Peninsula, and eleven continuous GPS stations distributed across the continent. The spatial pattern of these velocities is not consistent with any postglacial rebound (PGR) model known to us. Four leading PGR models appear to be overpredicting uplift rates in the Transantarctic Mountains and West Antarctica and underpredicting them in the peninsula north of 65°. This discrepancy cannot be explained in terms of an elastic response to modern ice loss (except, perhaps, in part of the peninsula). Therefore, our initial geodetic results suggest that most GRACE ice mass rate estimates, which are critically dependent on a PGR correction, are systematically biased and are overpredicting ice loss for the continent as a whole.

Local tsunami warnings: Perspectives from recent large events

We demonstrate a flexible strategy for local tsunami warning that relies on regional geodetic and seismic stations. Through retrospective analysis of four recent tsunamigenic events in Japan and Chile, we show that rapid earthquake source information, provided by methodologies developed for earthquake early warning, can be used to generate timely estimates of maximum expected tsunami amplitude with enough accuracy for tsunami warning. We validate the technique by comparing to detailed models of earthquake source and tsunami propagation as well as field surveys of tsunami inundation. Our approach does not require deployment of new geodetic and seismic instrumentation in many subduction zones and could be implemented rapidly by national monitoring and warning agencies. We illustrate the potential impact of our method with a detailed comparison to the actual timeline of events during the recent 2015 Mw8.3 Illapel, Chile, earthquake and tsunami that prompted the evacuation of 1 million people.

Isolating active orogenic wedge deformation in the southern Subandes of Bolivia

A new GPS-derived surface velocity field for the central Andean backarc permits an assessment of orogenic wedge deformation across the southern Subandes of Bolivia, where recent studies suggest that great earthquakes (>Mw 8) are possible. We find that the backarc is not isolated from the main plate boundary seismic cycle. Rather, signals from subduction zone earthquakes contaminate the velocity field at distances greater than 800 km from the Chile trench. Two new wedge-crossing velocity profiles, corrected for seasonal and earthquake affects, reveal distinct regions that reflect (1) locking of the main plate boundary across the high Andes, (2) the location of and loading rate at the back of orogenic wedge, and (3) an east flank velocity gradient indicative of decollement locking beneath the Subandes. Modeling of the Subandean portions of the profiles indicates along-strike variations in the decollement locked width (WL) and wedge loading rate; the northern wedge decollement has a WL of ~100 km while accumulating slip at a rate of ~14 mm/yr, whereas the southern wedge has a WL of ~61 km and a slip rate of ~7 mm/yr. When compared to Quaternary estimates of geologic shortening and evidence for Holocene internal wedge deformation, the new GPS-derived wedge loading rates may indicate that the southern wedge is experiencing a phase of thickening via reactivation of preexisting internal structures. In contrast, we suspect that the northern wedge is undergoing an accretion or widening phase primarily via slip on relatively young thrust-front faults.

Virtual array beamforming of GPS TEC observations of coseismic ionospheric disturbances near the Geomagnetic South Pole triggered by teleseismic megathrusts

We identified co-seismic ionospheric disturbances (CID) in Antarctica generated by the 2010 Maule and the 2011 Tohoku-Oki earthquakes analyzing TEC data with a modified beamforming technique. Beamforming in Antarctica, however, is not straightforward due to the effects of array deformation and atmospheric neutral wave-ionospheric plasma coupling. We take these effects into account and present a method to invert for the seismically generated acoustic wave using TEC observations. The back azimuths, speeds and waveforms obtained by the beamform are in excellent agreement with the hypothesis that the TEC signals are generated by the passage of Rayleigh waves from the Maule and Tohoku-Oki earthquakes. The Tohoku-Oki earthquake is ~12,500 km from Antarctica, making this the farthest observation of CIDs to date using GPS.

Observed rapid bedrock uplift in Amundsen Sea Embayment promotes ice-sheet stability

A quick rebound for Antarctic crust Earths crust deforms under the load of glaciers and ice sheets. When these masses are removed, the crust rebounds at a time scale determined by the viscosity of the upper mantle. Using GPS, Barletta et al. found that the viscosity of the mantle under the West Antarctic Ice Sheet is much lower than expected. This means that as ice is lost, the crust rebounds much faster than previously expected. Although estimates of total ice loss have to be revised upward, the surprising finding indicates that the ice sheet may stabilize against catastrophic collapse. Science, this issue p. 1335 A new viscosity model requires a much lower viscosity under the Amundsen Sea Embayment, stabilizing the Antarctic Ice Sheet. The marine portion of the West Antarctic Ice Sheet (WAIS) in the Amundsen Sea Embayment (ASE) accounts for one-fourth of the cryospheric contribution to global sea-level rise and is vulnerable to catastrophic collapse. The bedrock response to ice mass loss, glacial isostatic adjustment (GIA), was thought to occur on a time scale of 10,000 years. We used new GPS measurements, which show a rapid (41 millimeters per year) uplift of the ASE, to estimate the viscosity of the mantle underneath. We found a much lower viscosity (4 × 1018 pascal-second) than global average, and this shortens the GIA response time scale to decades up to a century. Our finding requires an upward revision of ice mass loss from gravity data of 10% and increases the potential stability of the WAIS against catastrophic collapse.