Geophysics and Seismology Commons^™

Paleoseismic Investigation Of Rupture On The Dead Coyote Fault In The Kittitas Valley, Washington, Garet Huddleston

All Master's Theses

Recently discovered fault scarps along the Dead Coyote Fault (DCF) in the northern Kittitas Valley (KV) of central Washington suggests active faulting in the late Quaternary, but constraints on the timing and potential magnitude of earthquakes along the fault zone are limited. The KV lies at the northwestern edge of the Yakima Fold and Thrust Belt (YFTB), a low-strain region where individual structures are capable of producing M~7 earthquakes.

This investigation uses surficial geologic mapping and topographic analysis of the DCF scarps, ground penetrating radar transects, and paleoseismic trenching to determine the rupture history of the DCF. The trench was …

Spatial Variations Of Stochastic Noise Properties In Gps Time Series, Xiaoxing He, Michael Simon Bos, Jean-Philippe Montillet, Rui Fernandes, Timothy I. Melbourne, Weiping Jiang, Wudong Li

Geological Sciences Faculty Scholarship

The noise in position time series of 568 GPS (Global Position System) stations across North America with an observation span of ten years has been investigated using solutions from two processing centers, namely, the Pacific Northwest Geodetic Array (PANGA) and New Mexico Tech (NMT). It is well known that in the frequency domain, the noise exhibits a power-law behavior with a spectral index of around −1. By fitting various noise models to the observations and selecting the most likely one, we demonstrate that the spectral index in some regions flattens to zero at long periods while in other regions it …

Stratigraphic Evidence Of Two Historical Tsunamis On The Semi-Arid Coast Of North-Central Chile, Jessica M. Depaolis, Tina Dura, Breanyn Macinnes, Lisa L. Ely, Marco Cisternas, Matías Carvajal, Hui Tang, Hermann M. Fritz, Cyntia Mizobe, Robert L. Wesson, Gino Figueroa, Nicole Brennan, Benjamin P. Horton, Jessica E. Pilarczyk, D. Reide Corbett, Benjamin C. Gill, Robert Weiss

Geological Sciences Faculty Scholarship

On September 16, 2015, a M_w 8.3 earthquake struck the north-central Chile coast, triggering a tsunami observed along 500 km of coastline, between Huasco (28.5°S) and San Antonio (33.5°S). This tsunami provided a unique opportunity to examine the nature of tsunami deposits in a semi-arid, siliciclastic environment where stratigraphic and sedimentological records of past tsunamis are difficult to distinguish. To improve our ability to identify such evidence, we targeted one of the few low-energy, organic-rich depositional environments in north-central Chile: Pachingo marsh in Tongoy Bay (30.3°S).

We found sedimentary evidence of the 2015 and one previous tsunami as tabular …

Toward Near‐Field Tsunami Forecasting Along The Cascadia Subduction Zone Using Rapid Gnss Source Models, Amy L. Williamson, Diego Melgar, Brendan W. Crowell, Diego Argas, Timothy I. Melbourne, Yong Wei, Kevin Kwong

Geological Sciences Faculty Scholarship

Over the past 15 years and through multiple large and devastating earthquakes, tsunami warning systems have grown considerably in their efficacy in providing timely and accurate forecasts to affected communities. However, one part of tsunami warning that still needs improvement is forecasts catered to local, near‐field communities in the time after an earthquake rupture but before coastal inundation. In this study, we test a rapid, Global Navigation Satellite Systems (GNSS)‐driven earthquake characterization model using a large data set of synthetic megathrust ruptures for its near‐field tsunami forecasting potential. We also provide a framework for tsunami forecasting that focuses on the …

Plate Boundary Trench Retreat And Dextral Shear Drive Intracontinental Fault-Slip Histories: Neogene Dextral Faulting Across The Gabbs Valley And Gillis Ranges, Central Walker Lane, Nevada, Jeffrey Lee, Andrew K. R. Hoxey, Andrew Calvert, Peter Dubyoski

All Faculty Scholarship for the College of the Sciences

The spatial-temporal evolution of intracontinental faults and the forces that drive their style, orientation, and timing are central to understanding tectonic processes. Intracontinental NW-striking dextral faults in the Gabbs Valley–Gillis Ranges (hereafter referred to as the GVGR), Nevada, define a structural domain known as the eastern Central Walker Lane located east of the western margin of the North American plate. To consider how changes in boundary type along the western margin of the North American plate influenced both the initiation and continued dextral fault slip to the present day in the GVGR, we combine our new detailed geologic mapping, structural …

Changing Mantle Sources And The Effects Of Crustal Passage On The Steens Basalt, Se Oregon: Chemical And Isotopic Constraints, Nicole E. Moore, A. L. Grunder, Wendy A. Bohrson, Richard W. Carlson, I. N. Bindeman

Geological Sciences Faculty Scholarship

Continental flood basalts are more prone to compositional modification from passage through thicker and (or) more felsic crust in comparison to their oceanic counterparts. The Steens Basalt in southeast Oregon (~17 Ma) is among the oldest and most mafic members of the Columbia River Basalt Group and provides a record of the early stages of flood basalt volcanism. We evaluate the balance of mantle sources in time during the onset of Columbia River Basalt Group magmatism and assess the effect of crustal passage using stratigraphically controlled Sr, Nd, Pb, Hf, Os, and O isotopic compositions, as well as whole rock …

On The Interaction Of The North Andes Plate With The Caribbean And South American Plates In Northwestern South America From Gps Geodesy And Seismic Data, Omar J. Pérez, Steven G. Wesnousky, Roberto De La Rosa, Julio Márquez, Redescal Uzcátegui, Christian Quintero, Luis Liberal, Héctor Mora-Páez, Walter Szeliga

All Faculty Scholarship for the College of the Sciences

We examine the hypocentral distribution of seismicity and a series of geodetic velocity vectors obtained from Global Positioning System observations between 1994 and 2015 both offshore and mainland northwestern South America (66°W–77°W; 8°N–14°N). Our analysis, that includes a kinematic block modelling, shows that east of the Caribbean–South American–North Andes plates’ triple junction at ∼68°W; 10.7°N, right-lateral easterly oriented shear motion (∼19.6 ± 2.0 mm yr⁻¹) between the Caribbean and South America plates is split along two easterly striking, right-lateral strike-slip subparallel fault zones: the San Sebastián fault that runs off-shore the Venezuelan coast and slips about 17.0±0.5mm …

The Uppermost Mantle Seismic Velocity And Viscosity Structure Of Central West Antarctica, J. P. O'Donnell, K. Selway, A. A. Nyblade, R. A. Brazier, D. A. Wiens, S. Anandakrishnan, R. C. Aster, Audrey D. Huerta, T. Wilson, J. Paul Winberry

All Faculty Scholarship for the College of the Sciences

Accurately monitoring and predicting the evolution of the West Antarctic Ice Sheet via secular changes in the Earth's gravity field requires knowledge of the underlying upper mantle viscosity structure. Published seismic models show the West Antarctic lithosphere to be ∼70–100 km thick and underlain by a low velocity zone extending to at least ∼200 km. Mantle viscosity is dependent on factors including temperature, grain size, the hydrogen content of olivine, the presence of partial melt and applied stress. As seismic wave propagation is particularly sensitive to thermal variations, seismic velocity provides a means of gauging mantle temperature. In 2012, a …

Unusual Geologic Evidence Of Coeval Seismic Shaking And Tsunamis Shows Variability In Earthquake Size And Recurrence In The Area Of The Giant 1960 Chile Earthquake, M. Cisternas, E. Garrett, R. Wesson, T. Dura, Lisa L. Ely

All Faculty Scholarship for the College of the Sciences

An uncommon coastal sedimentary record combines evidence for seismic shaking and coincident tsunami inundation since AD 1000 in the region of the largest earthquake recorded instrumentally: the giant 1960 southern Chile earthquake (Mw 9.5). The record reveals significant variability in the size and recurrence of megathrust earthquakes and ensuing tsunamis along this part of the Nazca-South American plate boundary. A 500-m long coastal outcrop on Isla Chiloé, midway along the 1960 rupture, provides continuous exposure of soil horizons buried locally by debris-flow diamicts and extensively by tsunami sand sheets. The diamicts flattened plants that yield geologically precise ages to correlate …

Plate Boundary Observatory And Related Networks: Gps Data Analysis Methods And Geodetic Products, Thomas A. Herring, Timothy I. Melbourne, Mark H. Murray, Michael A. Floyd, Walter M. Szeliga, Robert W. King, David A. Phillips, Christine M. Puskas, Marcelo Santillan, Lei Wang

All Faculty Scholarship for the College of the Sciences

The Geodesy Advancing Geosciences and EarthScope (GAGE) Facility Global Positioning System (GPS) Data Analysis Centers produce position time series, velocities, and other parameters for approximately 2000 continuously operating GPS receivers spanning a quadrant of Earth’s surface encompassing the high Arctic, North America, and Caribbean. The purpose of this review is to document the methodology for generating station positions and their evolution over time and to describe the requisite trade-offs involved with combination of results. GAGE GPS analysis involves formal merging within a Kalman filter of two independent, loosely constrained solutions: one is based on precise point positioning produced with the …

2015 Gsa International Distinguished Lecturer Tours Key Locations In Central And South America, Lisa L. Ely

Geological Sciences Faculty Scholarship

As the International Lecturer from North America, I visited 11 destinations in Ecuador, Peru, Chile, Costa Rica, Mexico, and Puerto Rico. The tour was designed to include locations where the topic of my presentation, “Following in the Footsteps of Darwin:
Combining Geological and Historical Evidence to Assess Earthquakes and Tsunami Hazards,” would have direct relevance.

Gps Constraints On Interplate Locking Within The Makran Subduction Zone, Elyse Frohling, Walter Szeliga

All Faculty Scholarship for the College of the Sciences

The Makran subduction zone is one of the last convergent margins to be investigated using space-based geodesy. While there is a lack of historical and modern instrumentation in the region, a sparse sampling of continuous and campaign measurements over the past decade has allowed us to make the first estimates of convergence rates. We combine GPS measurements from 20 stations located in Iran, Pakistan and Oman along with hypocentral locations from the International Seismological Centre to create a preliminary 3-D estimate of the geometry of the megathrust, along with a preliminary fault-coupling model for the Makran subduction zone. Using a …

Crustal And Upper-Mantle Structure Beneath Ice-Covered Regions In Antarctica From S-Wave Receiver Functions And Implications For Heat Flow, C. Ramirez, Andrew A. Nyblade, S. E. Hansen, Douglas A. Wiens, Sridhar Anandakrishnan, Richard C. Aster, Audrey D. Huerta, Partick Shore, Terry Wilson

Geological Sciences Faculty Scholarship

S-wave receiver functions (SRFs) are used to investigate crustal and upper-mantle structure beneath several ice-covered areas of Antarctica. Moho S-to-P (Sp) arrivals are observed at ~6–8 s in SRF stacks for stations in the Gamburtsev Mountains (GAM) and Vostok Highlands (VHIG), ~5–6 s for stations in the Transantarctic Mountains (TAM) and the Wilkes Basin (WILK), and ~3–4 s for stations in the West Antarctic Rift System (WARS) and the Marie Byrd Land Dome (MBLD). A grid search is used to model the Moho Sp conversion time with Rayleigh wave phase velocities from 18 to 30 s period to estimate crustal …

Upper Mantle Structure Of Central And West Antarctica From Array Analysis Of Rayleigh Wave Phase Velocities, David S. Heeszel, Douglas A. Wiens, Sridhar Anandakrishnan, Richard C. Aster, Ian W.D. Dalziel, Audrey D. Huerta, Andrew A. Nyblade, Terry J. Wilson, J. Paul Winberry

Geological Sciences Faculty Scholarship

The seismic velocity structure of Antarctica is important, both as a constraint on the tectonic history of the continent and for understanding solid Earth interactions with the ice sheet. We use Rayleigh wave array analysis methods applied to teleseismic data from recent temporary broadband seismograph deployments to image the upper mantle structure of central and West Antarctica. Phase velocity maps are determined using a two–plane wave tomography method and are inverted for shear velocity using a Monte Carlo approach to estimate three-dimensional velocity structure. Results illuminate the structural dichotomy between the East Antarctic Craton and West Antarctica, with West Antarctica …

Spatiotemporal Slip Rate Variations Along Surprise Valley Fault In Relation To Pleistocene Pluvial Lakes, Brian N. Marion

All Master's Theses

Using mapped paleoshoreline features with high-resolution topographic data and obtained radiocarbon dates on paleoshoreline tufas, I documented precise fault offsets of dated features over the last 25 ka along the Surprise Valley Fault (SVF). Fault offset measured in three lake sections within Surprise Valley ranged from 3.6 m in the southern section to 14.4 m in the central section. The offset paleoshorelines are dated to the late Pleistocene (<22 >ka) and were formed during the latest impoundment of pluvial Lake Surprise since the last glacial maximum. Slip rates vary along strike, assuming a fault dip of 68° with 0.25 ± …

Hidden Intrabasin Extension: Evidence For Dike-Fault Interaction From Magnetic, Gravity, And Seismic Reflection Data In Surprise Valley, Northeastern California, Noah D. Athens, Jonathan M.G. Glen, Simon L. Klemperer, Anne E. Egger, Valentina C. Fontiveros

All Faculty Scholarship for the College of the Sciences

The relative contributions of tectonic and magmatic processes to continental rifting are highly variable. Magnetic, gravity, and seismic reflection data from Surprise Valley, California, in the northwest Basin and Range, reveal an intrabasin, fault-controlled, ~10-m-thick dike at a depth of ~150 m, providing an excellent example of the interplay between faulting and dike intrusion. The dike, likely a composite structure representing multiple successive intrusions, is inferred from modeling a positive magnetic anomaly that extends ~35 km and parallels the basin-bounding Surprise Valley normal fault on the west side of the valley. A two-dimensional high-resolution seismic reflection profile acquired across the …

Slip Pulse And Resonance Of The Kathmandu Basin During The 2015 Gorkha Earthquake, Nepal, John Galetzka, Walter Szeliga

All Faculty Scholarship for the College of the Sciences

Detailed geodetic imaging of earthquake ruptures enhances our understanding of earthquake physics and associated ground shaking. The 25 April 2015 moment magnitude 7.8 earthquake in Gorkha, Nepal was the first large continental megathrust rupture to have occurred beneath a high-rate (5-hertz) Global Positioning System (GPS) network. We used GPS and interferometric synthetic aperture radar data to model the earthquake rupture as a slip pulse ~20 kilometers in width, ~6 seconds in duration, and with a peak sliding velocity of 1.1 meters per second, which propagated toward the Kathmandu basin at ~3.3 kilometers per second over ~140 kilometers. The smooth slip …

The Mantle Transition Zone Beneath West Antarctica: Seismic Evidence For Hydration And Thermal Upwellings, E. L. Emry, A. A. Nyblade, J. Juliá, S. Anandakrishnan, R. C. Aster, D. A. Wiens, Audrey D. Huerta, T. J. Wilson

All Faculty Scholarship for the College of the Sciences

Although prior work suggests that a mantle plume is associated with Cenozoic rifting and volcanism in West Antarctica, the existence of a plume remains conjectural. Here we use P wave receiver functions (PRFs) from the Antarctic POLENET array to estimate mantle transition zone thickness, which is sensitive to temperature perturbations, throughout previously unstudied parts of West Antarctica. We obtain over 8000 high‐quality PRFs using an iterative, time domain deconvolution method filtered with a Gaussian width of 0.5 and 1.0, corresponding to frequencies less than ∼0.24 and ∼0.48 Hz, respectively. Single‐station and common conversion point stacks, migrated to depth using the …

Estimation Of Offsets In Gps Time-Series And Application To The Detection Of Earthquake Deformation In The Far-Field, Jean-Philippe Montillet, S. D. P. Williams, A. Koulali, S. C. Mcclusky

All Faculty Scholarship for the College of the Sciences

Extracting geophysical signals from Global Positioning System (GPS) coordinate time-series is a well-established practice that has led to great insights into how the Earth deforms. Often small discontinuities are found in such time-series and are traceable to either broad-scale deformation (i.e. earthquakes) or discontinuities due to equipment changes and/or failures. Estimating these offsets accurately enables the identification of coseismic deformation estimates in the former case, and the removal of unwanted signals in the latter case which then allows tectonic rates to be estimated more accurately. We develop a method to estimate accurately discontinuities in time series of GPS positions at …

Clockwise Rotation Of The Brahmaputra Valley Relative To India: Tectonic Convergence In The Eastern Himalaya, Naga Hills, And Shillong Plateau, ‪Philippe Vernant, R. Bilham, Walter Szeliga, D. Drupka, S. Kalita, A. K. Bhattacharyya, V. K. Gaur

All Faculty Scholarship for the College of the Sciences

GPS data reveal that the Brahmaputra Valley has broken from the Indian Plate and rotates clockwise relative to India about a point a few hundred kilometers west of the Shillong Plateau. The GPS velocity vectors define two distinct blocks separated by the Kopili fault upon which 2–3 mm/yr of dextral slip is observed: the Shillong block between longitudes 89 and 93°E rotating clockwise at 1.15°/Myr and the Assam block from 93.5°E to 97°E rotating at ≈1.13°/Myr. These two blocks are more than 120 km wide in a north‐south sense, but they extend locally a similar distance beneath the Himalaya and …

Thermodynamic Model For Energy-Constrained Open-System Evolution Of Crustal Magma Bodies Undergoing Simultaneous Recharge, Assimilation And Crystallization: The Magma Chamber Simulator, Wendy A. Bohrson, Frank J. Spera, Mark S. Ghiorso, Guy Brown, Jeffrey Creamer, Aaron Mayfield

All Faculty Scholarship for the College of the Sciences

The Magma Chamber Simulator quantifies the impact of simultaneous recharge, assimilation and crystallization through mass and enthalpy balance in a multicomponent–multiphase (melt + solids ± fluid) composite system. As a rigorous thermodynamic model, the Magma Chamber Simulator computes phase equilibria and geochemical evolution self-consistently in resident magma, recharge magma and wallrock, all of which are connected by specified thermodynamic boundaries, to model an evolving open-system magma body. In a simulation, magma cools from its liquidus temperature, and crystals ± fluid are incrementally fractionated to a separate cumulate reservoir. Enthalpy from cooling, crystallization, and possible magma recharge heats wallrock from its …

Seismic And Geodetic Evidence For Grounding-Line Control Of Whillans Ice Stream Stick-Slip Events, Martin J. Pratt, J. Paul Winberry, Douglas A. Wiens, Sridhar Anandakrishnan, Richard B. Alley

All Faculty Scholarship for the College of the Sciences

The tidally modulated, stick‐slip events of Whillans Ice Stream in West Antarctica produce seismic energy from three locations near the grounding line. Using ice velocity records obtained by combining time series from colocated broadband seismometers and GPS receivers installed on the ice stream during the 2010–2011 and 2011–2012 austral summers, along with far‐field seismic recordings of elastic waves, we locate regions of high rupture velocity and stress drop. These regions, which are analogous to “asperities” in traditional seismic fault studies, are areas of elevated friction at the base of the ice stream. Slip events consistently initiate at one of two …

Tidal Pacing, Skipped Slips And The Slowdown Of Whillans Ice Stream, Antarctica, J. Paul Winberry, Sridhar Anandakrishnan, Richard B. Alley, Douglas A. Wiens, Martin J. Pratt

All Faculty Scholarship for the College of the Sciences

We summarize new observations of the deceleration and stick–slip motion of Whillans Ice Stream (WIS), Antarctica. We refine the location of the large sticky spots that resist motion between slip events, the locations of which are controlled by the patterns of subglacial water flow. Our examination of the long-term velocity time series for the ice stream reveals that the decadal-scale deceleration is not occurring at a steady rate, but varies at the sub-decadal timescale. This unsteady deceleration modulates the temporal evolution of a broad (~50 km across) surface-elevation bulge forming at the junction between the relatively narrow upstream portion of …

Thermobarometry And Geochronology In The Dulan Region, North Qaidam Ultrahigh-Pressure Terrane: Resolving Spatial Variation Of Ages, Temperatures, And Pressures, Megan Eileen Regel

All Master's Theses

The Dulan ultrahigh-pressure (UHP) terrane exposes 449-422 Ma eclogites in the east and 449-418 Ma high-pressure granulites to the west, but the relationship between the units is poorly understood. Zirconium-in-rutile temperatures show an east-west 3- 5°C/km gradient in four UHP eclogites: 662-716±10°C; garnet-omphacite-phengite thermobarometry of peak conditions yields 659-730±82°C, 25.3-27.1±3.2 kbars. Zircon U-Pb ages of four eclogites are the same within uncertainty, averaging 434±2.9 Ma. Tiin- zircon thermometry and rare earth element patterns indicate these ages represent eclogite facies metamorphism. One high-pressure/high-temperature (HP/HT) granulite yields a minimum rutile temperature of 773±4°C. Previous thermobarometry of highpressure granulites yields ≥800°C and 14-17 …

Proxy-Based Reconstructions Of Earthquakes And Tsunamis At Quidico, South-Central Chile, Isabel J. Hong

All Master's Theses

We utilized geomorphic, microfossil, sedimentological, and stratigraphic methods to investigate the paleotsunami history at Quidico, Chile (38.1° S, 73.2° W). A combination of pits, cores, and riverbank exposures provide stratigraphic evidence of 8 sand layers at Quidico, including tsunami deposits from 2010 and 1960. The sands are laterally continuous but display landward thinning, disappearing completely from the stratigraphy 1.2 km inland from the coast. Intervening sediment layers consist of organicrich silts and peat. Grain size analysis shows the sand units are bimodal in distribution whereas the organic-rich silts and peat have a broader, bimodal distribution. Lithological analysis showed similarities among …

Controls On Fault Geometry During Early Stages Of Extension In The Larkspur Hills, Northwest Basin And Range, Diana Jean Strickley

All Master's Theses

Detailed analyses of normal faults in the Larkspur Hills, CA-NV, northwest Basin and Range, offer insight into factors controlling normal fault initiation, growth, and distribution. N-trending faults in the southern portion of the study area share trends of major range-bounding structures and Pliocene linear volcanic vents; in contrast, NNW- and NNE- trending faults dominate further north and into south-central Oregon. Stress analyses and comparison with experimental and field data suggest that preexisting structures control faults in the northern Larkspur Hills, while faults form perpendicular to σ3 in the southern hills. The change in fault orientations is abrupt, occurring across a …

Spatial Relationship Between Gps Slip And Seismic Tremor During Cascadia Slow Slip Events, Hillary Lynn Goodner

All Master's Theses

We model GPS deformation and timing of seismic tremor associated with transient deformation in Cascadia to test the hypothesis that tremor and slip occur synchronously but are spatially offset. For the period 2010–2013, we use seismic tremor data with a duration-moment relationship to predict GPS time series and compare them to observations. We find that observed GPS displacements are best predicted when tremor locations on the plate interface are shifted 15 km up-dip of their published epicenter. To test whether the spatial offset of tremor and slip is due to systematic mislocation of published epicenters, we attempt to identify individual …

Glacier Slip And Seismicity Induced By Surface Melt, Peter L. Moore, J. Paul Winberry, Neal R. Iverson, Knut A. Christianson, Sridhar Anandakrishnan, Miriam Jackson, Mark E. Mathison, Denis Cohen

All Faculty Scholarship for the College of the Sciences

Many of the key processes governing fast glacier flow involve interaction between a glacier and its basal hydrological system, which is hidden from direct observation. Passive seismic monitoring has shown promise as a tool for remotely monitoring basal processes, but lack of glacier-bed access prevents clear understanding of the relationships between subglacial processes and corresponding seismic emissions. Here we describe direct measurements of basal hydrology, sliding, and broadband seismicity made in a unique subglacial facility in Norway during the onset of two summer melt seasons. In the most pronounced of these episodes, rapid delivery of surface meltwater to the bed …

2012 Haida Gwaii Quake: Insight Into Cascadia's Subduction Extent, Walter Szeliga

Geological Sciences Faculty Scholarship

The limits of Cascadia were first defined to contain nearly the entire margin of the Pacific Northwest, from Cape Mendocino through the Alaska Panhandle [Schuchert, 1910; Schuchert and Barrell, 1914]. Since that time, the boundary of Cascadia has shrunk to become essentially synonymous with the region where the Juan de Fuca plate subducts beneath the North American plate. As a consequence, seismic hazard assessments in the Pacific Northwest have conventionally focused on the potential for large megathrust earthquakes along the interface of the Juan de Fuca and North American plates.

Nucleation And Seismic Tremor Associated With The Glacial Earthquakes Of Whillans Ice Stream, Antarctica, J. Paul Winberry, Sridhar Anandakrishnan, Douglas A. Wiens, Richard B. Alley

All Faculty Scholarship for the College of the Sciences

The ability to monitor transient motion along faults is critical to improving our ability to understand many natural phenomena such as landslides and earthquakes. Here, we usedata from a GPS and seismometer network that were deployed to monitor the regularly repeating glacial earthquakes of Whillans Ice Stream, West Antarctica to show that a unique pattern of precursory slip precedes complete rupture along the bed of the ice stream. Additionally, we show that rupture can be independently tracked by increased levels of microseismic activity, including harmonic tremor, that are coincident with the onset of slip at any location, thus providing a …