Earth Sciences Commons^™

Noise Characteristics Of Operational Real‐Time High‐Rate Gnss Positions In A Large Aperture Network, Diego Melgar, Brendan W. Crowell, Timothy I. Melbourne, Walter Szeliga, Marcelo Santillan, Craig Scivner

Geological Sciences Faculty Scholarship

Large earthquakes are difficult to model in real‐time with traditional inertial seismic measurements. Several algorithms that leverage high‐rate real‐time Global Navigation Satellite Systems (HR‐GNSS) positions have been proposed, and it has been shown that they can supplement the earthquake monitoring effort. However, analyses of the long‐term noise behavior of high‐rate real‐time GNSS positions, which are important to understand how the data can be used operationally by monitoring agencies, have been limited to just a few sites and to short time spans. Here, we show results from an analysis of the noise characteristics of 1 year of positions at 213 GNSS …

Dextral, Normal, And Sinistral Faulting Across The Eastern California Shear Zone-Mina Deflection Transition, California-Nevada, Usa, Kevin Delano, Jeffrey Lee, Rachelle Roper, Andrew Calvert

All Faculty Scholarship for the College of the Sciences

Strike-slip faults commonly include extensional and contractional bends and stepovers, whereas rotational stepovers are less common. The Volcanic Tableland, Black Mountain, and River Spring areas (California and Nevada, USA) (hereafter referred to as the VBR region) straddle the transition from the dominantly NW-striking dextral faults that define the northwestern part

of the eastern California shear zone into a rotational stepover characterized by dominantly NE-striking sinistral faults that define the southwestern Mina deflection. New detailed geologic mapping, structural studies, and 40Ar/39Ar geochronology across the VBR region allow us to calculate Pliocene to Pleistocene fault slip rates and test predictions for the …

Seismic Observations Of Crevasse Growth Following Rain-Induced Glacier Acceleration, Haupapa/Tasman Glacier, New Zealand, Samuel Taylor-Offord, Huw Horgan, John Townend, J. Paul Winberry

All Faculty Scholarship for the College of the Sciences

Changing rates of water input can affect both the flow of glaciers and ice sheets and their propensity to crevasse. Here we examine geodetic and seismic observations during two substantial (10–18-times background velocity) rain-induced glacier accelerations at Haupapa/Tasman Glacier, New Zealand. Changes in rain rate result in glacier acceleration and associated uplift, which propagate down-glacier. This pattern of acceleration results in a change to the strain rate field, which correlates with an order of magnitude increase in the apparent seismicity rate and an overall down-glacier migration in located seismicity. After each acceleration event the apparent seismicity rate decreases to below …

Tidal And Spatial Variability Of Flow Speed And Seismicity Near The Grounding Zone Of Beardmore Glacier, Antarctica, Jade Cooley, J. Paul Winberry, Michelle Koutnik, Howard Conway

All Faculty Scholarship for the College of the Sciences

GPS measurements of tidal modulation of ice flow and seismicity within the grounding zone of Beardmore Glacier show that tidally induced fluctuations of horizontal flow are largest near the grounding line and decrease downstream. Seismic activity is continuous, but peaks occur on falling and rising tides. Beamforming methods reveal that most seismic events originate from two distinct locations, one on the grid-north side of the grounding zone, and one on the grid-south side. The broad pattern of deformation generated as Beardmore Glacier merges with the Ross Ice Shelf results in net extension along the grid-north side of the grounding zone …

Seismic Tremor Reveals Spatial Organization And Temporal Changes Of Subglacial Water System, Margot E. Vore, Timothy C. Bartholomaus, J. Paul Winberry, Jacob I. Walter, Jason M. Amundson

All Faculty Scholarship for the College of the Sciences

Subglacial water flow impacts glacier dynamics and shapes the subglacial environment. However, due to the challenges of observing glacier beds, the spatial organization of subglacial water systems and the time scales of conduit evolution and migration are largely unknown. To address these questions, we analyze 1.5‐ to 10‐Hz seismic tremor that we associate with subglacial water flow, hat is, glaciohydraulic tremor, at Taku Glacier, Alaska, throughout the 2016 melt season. We use frequency‐dependent polarization analysis to estimate glaciohydraulic tremor propagation direction (related to the subglacial conduit location) and a degree day melt model to monitor variations in melt‐water input. We …

The Crust And Upper Mantle Structure Of Central And West Antarctica From Bayesian Inversion Of Rayleigh Wave And Receiver Functions, Weisen Shen, Audrey D. Huerta, J. Paul Winberry

All Faculty Scholarship for the College of the Sciences

We construct a new seismic model for central and West Antarctica by jointly inverting Rayleigh wave phase and group velocities along with P wave receiver functions. Ambient noise tomography exploiting data from more than 200 seismic stations deployed over the past 18 years is used to construct Rayleigh wave phase and group velocity dispersion maps. Comparison between the ambient noise phase velocity maps with those constructed using teleseismic earthquakes confirms the accuracy of both results. These maps, together with P receiver function waveforms, are used to construct a new 3-D shear velocity (Vs) model for the crust and uppermost mantle …

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 …

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 …

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 …

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 …

The Crustal Thickness Of West Antarctica, Julien A. Chaput, Richard C. Aster, Audrey D. Huerta, X. Sun, Andrew Lloyd, Douglas A. Wiens, Andrew Nyblade, Sridhar Anandakrishnan, J. Paul Winberry, Terry Wilson

Geological Sciences Faculty Scholarship

P-to-S receiver functions (PRFs) from the Polar Earth Observing Network (POLENET) GPS and seismic leg of POLENET spanning West Antarctica and the Transantarctic Mountains deployment of seismographic stations provide new estimates of crustal thickness across West Antarctica, including the West Antarctic Rift System (WARS), Marie Byrd Land (MBL) dome, and the Transantarctic Mountains (TAM) margin. We show that complications arising from ice sheet multiples can be effectively managed and further information concerning low-velocity subglacial sediment thickness may be determined, via top-down utilization of synthetic receiver function models. We combine shallow structure constraints with the response of deeper layers using a …

Seismic Slip Deficit In The Kashmir Himalaya From Gps Observations, Celia Schiffman, Bikram Singh Bali, Walter Szeliga, Roger Bilham

Faculty Scholarship for the Cascadia Hazards Institute

GPS measurements in Kashmir Himalaya reveal rangenormal convergence of 11±1 mm/yr with dextral shear of 5±1 mm/yr. The transition from a fully locked 170 km wide décollement to the unrestrained descending Indian plate occurs at ~25 km depth over an ~23 km wide transition zone. The convergence rate is consistent with the lower bounds of geological estimates for the Main Frontal Thrust, Riasi, and Balapora fault systems, on which no surface slip has been reported in the past millennium. Of the 14 damaging Kashmir earthquakes since 1123, none may have exceeded M_w = 7.6. Therefore, either a seismic moment …

Evolution Of Dike Opening During The March 2011 Kamoamoa Fissure Eruption, Kīlauea Volcano, Hawai'i, Paul Lundgren, Michael Poland, Asta Miklius, Tim Orr, Sang-Ho Yun, Eric Fielding, Zhen Liu, Scott Hensley, Susan Owen, Akiko Tanaka, Walter Szeliga

Faculty Scholarship for the Cascadia Hazards Institute

The 5–9 March 2011 Kamoamoa fissure eruption along the east rift zone of Kīlauea Volcano, Hawai`i, followed months of pronounced inflation at Kīlauea summit. We examine dike opening during and after the eruption using a comprehensive interferometric synthetic aperture radar (InSAR) data set in combination with continuous GPS data. We solve for distributed dike displacements using a whole Kīlauea model with dilating rift zones and possibly a deep décollement. Modeled surface dike opening increased from nearly 1.5 m to over 2.8 m from the first day to the end of the eruption, in agreement with field observations of surface fracturing. …

Paleoseismologic Evidence For Holocene Activity On The Pinto Mountain Fault, Twentynine Palms, California, Ana Maria Cadena

All Master's Theses

Excavations across the Pinto Mountain fault in Twentynine Palms, California exposed faulted strata across a 32-m wide zone. Trench wall exposures revealed clear evidence for five ground-rupturing events during the Holocene, and two additional events in the late Pleistocene. Optically stimulated luminescence ages from alluvial sediments suggest that the most recent event occurred between 1.7-2.9 ka B.P. and the penultimate event between 2.7-4.2 ka B.P.. Prior to the penultimate event, there were five ground-rupturing earthquakes on the eastern Pinto Mountain fault between 3.5-13.6 ka B.P.. The average recurrence interval since 13.3-13.6 ka B.P. is 1510-1680 years, and 1200-1500 years in …

Gps Constraints On Indo-Asian Convergence In The Bhutan Himalaya: Segmentation And Potential For A 8.2-8.8 Mw Earthquake, Dowchu Drukpa, Phuntsho Pelgay, Anjan Bhattacharya, Phillipe Vernant, Walter Szeliga, Roser Bilham

Faculty Scholarship for the Cascadia Hazards Institute

The seismogenic setting of Bhutan is unusual due to its lower-than-average 20th century seismic moment release (Drukpa et al. 2006), its absence of a reliable historical record, and its unusual location near the Shillong plateau where a great earthquake in 1897 resulted in ≈10 m of N/S shortening of the Indian plate to its south (Gahalaut et al. 2011). Despite these indicators that lower than normal convergence velocities should currently prevail, the GPS velocity between Shillong and Lhasa suggests that convergence in Sikkim and Bhutan occurs at velocities exceeding 20 mm/yr. GPS points between the Greater Himalaya and the Shillong …

Interseismic Strain Accumulation Along The Western Boundary Of The Indian Subcontinent, Walter Szeliga, Roger Bilham, Din Mohammad Kakar, Sarosh H. Lodi

All Faculty Scholarship for the College of the Sciences

Despite an overall sinistral slip rate of ≈3 cm/yr, few major earthquakes have occurred in the past 200 years along the Chaman fault system, the western boundary of the India Plate with the Eurasia Plate. GPS and InSAR data reported here indicate sinistral shear velocities of 8–17 mm/yr across the westernmost branches of the fault system, suggesting that a significant fraction of the plate boundary slip is distributed in the fold and fault belt to the east. At its southernmost on‐land segment (≈26°N), near the triple junction between the Arabia, Eurasia, and India Plates, we find the velocity across the …

Interseismic Strain Accumulation Along The Western Boundary Of The Indian Subcontinent, Walter Szeliga, Roger Bilham, Din Mohammad Kakar, Sarosh H. Lodi

Faculty Scholarship for the Cascadia Hazards Institute

Despite an overall sinistral slip rate of approximately 3 cm/yr, few major earthquakes have occurred in the past 200 years along the Chaman fault system, the western boundary of the India Plate with the Eurasia Plate. GPS and InSAR data reported here indicate sinistral shear velocities of 8-17 mm/yr across the westernmost branches of the fault system, suggesting that a significant fraction of the plate boundary slip is distributed in the fold and fault belt to the east. At its southernmost on-land segment, near the triple junction between the Arabia, Eurasia, and India Plates, we find the velocity across the …

A Catalog Of Felt Intensity Data For 570 Earthquakes In India From 1636 To 2009, Stacey Martin, Walter Szeliga

Faculty Scholarship for the Cascadia Hazards Institute

Eight thousand three hundred thirty-nine intensity observations have been evaluated for earthquakes that occurred on the Indian subcontinent and surrounding plate boundaries from the seventeenth century to the present. They characterize 570 earthquakes, more than 90% of which occurred in the past two centuries. The electronic supplement to this article lists these data using European Macroseismic Scale (EMS-98) intensities with their geographic coordinates. We summarize these data graphically in the form of a spatially averaged intensity map for the subcontinent, a map that emphasizes the features of many previously published earthquake hazard maps for the Indian plate, but which more …

Intensity, Magnitude, Location, And Attenuation In India For Felt Earthquakes Since 1762, Walter Szeliga, Susan Hough, Stacey Martin, Roger Bilham

Faculty Scholarship for the Cascadia Hazards Institute

A comprehensive, consistently interpreted new catalog of felt intensities for India (Martin and Szeliga, 2010, this issue) includes intensities for 570 earthquakes; instrumental magnitudes and locations are available for 100 of these events. We use the intensity values for 29 of the instrumentally recorded events to develop new intensity versus attenuation relations for the Indian subcontinent and the Himalayan region. We then use these relations to determine the locations and magnitudes of 234 historical events, using the method of Bakun and Wentworth (1997). For the remaining 336 events, intensity distributions are too sparse to determine magnitude or location. We evaluate …

Future Cascadia Megathrust Rupture Delineated By Episodic Tremor And Slip, James S. Chapman, Timothy I. Melbourne

Faculty Scholarship for the Cascadia Hazards Institute

A suite of 15 episodic tremor and slip events imaged between 1997 and 2008 along the northern Cascadia subduction zone suggests future coseismic rupture will extend to 25 km depth, or approximately 60 km inland of the Pacific coast, rather than stopping offshore at 15 km depth. An ETS-derived coupling profile accurately predicts GPS measured interseismic deformation of the overlying North American plate, as measured by approximately 50 continuous GPS stations across western Washington State. When extrapolated over the 550-year average recurrence interval of Cascadia megathrust events, the coupling model also replicates the pattern and amplitude of coseismic coastal subsidence …

Seismic And Geodetic Constraints On Cascadia Slow Slip, Timothy I. Melbourne, Aaron G. Wech, Kenneth C. Creager

Faculty Scholarship for the Cascadia Hazards Institute

Automatically detected and located tremor epicenters from episodic tremor and slip (ETS) episodes in northern Cascadia provide a high-resolution map of Washington’s slow slip region. Thousands of epicenters from the past four ETS events from 2004 to 2008 provide detailed map-view constraints that correlate with geodetic estimates of the simultaneous slow slip. Each of these ETS events exhibits remarkable similarity in the timing and geographic distribution of tremor density and geodetically inferred slip. Analysis of the latest 15-month inter-ETS period also reveals ageodetic tremor activity similar both in duration and extent to ETS tremor. Epicenters from both ETS and inter- …

Fold And Thrust Partitioning In A Contracting Fold Belt: Insights From The 1931 Mach Earthquake In Baluchistan, Walter Szeliga, Roger Bilham, Daniel Schelling, Din Mohamed Kakar, Sarosh Lodi

Faculty Scholarship for the Cascadia Hazards Institute

Surface deformation associated with the 27 August 1931 earthquake near Mach in Baluchistan is quantified from spirit-leveling data and from detailed structural sections of the region interpreted from seismic reflection data constrained by numerous well logs. Mean slip on the west dipping Dezghat/Bannh fault system amounted to 1.2 m on a 42 km x 72 km thrust plane with slip locally attaining 3.2 m up dip of an inferred locking line at approximately 9 km depth. Slip also occurred at depths below the interseismic locking line. In contrast, negligible slip occurred in the 4 km near the interseismic locking line. …

Gps Constraints On 34 Slow Slip Events Within The Cascadia Subduction Zone, 1997–2005, Walter Szeliga, Timothy I. Melbourne, Marcelo Santillan, Meghan Miller

Faculty Scholarship for the Cascadia Hazards Institute

Refinements to GPS analyses in which we factor geodetic time series to better estimate both reference frames and transient deformation resolve 34 slow slip events located throughout the Cascadia subduction zone from 1997 through 2005. Timing of transient onset is determined with wavelet transformation of geodetic time series. Thirty continuous stations are included in this study, ranging from northern California to southwestern British Columbia. Our improvements in analysis better resolve the largest creep events and also identify many smaller events. At 48.5 degrees N latitude, a 14-month average recurrence interval has been observed over eight events since 1997. Farther north …

Extent And Duration Of The 2003 Cascadia Slow Earthquake, Walter Szeliga, Timothy I. Melbourne, M. Meghan Miller, Victor Marcelo Santillan

Faculty Scholarship for the Cascadia Hazards Institute

Inversion of continuous GPS measurements from the Pacific Northwest show the 2003 Cascadia slow earthquake to be among the largest of ten transients recognized here. Twelve stations bracketing slow slip indicate transient slip propagated bi-directionally from initiation in the southern Puget basin, reaching 300 km along-strike over a period of seven weeks. This event produced, for the first time, resolvable vertical subsidence, and horizontal displacement reaching six mm in southern Washington State. Inverted for non-negative thrust slip, a maximum of 3.8 cm of slip is inferred, centered at 28 km depth near the sharp arch in the subducting Juan de …

Southern Cascadia Episodic Slow Earthquakes, Walter Szeliga, Timothy I. Melbourne, M. Meghan Miller, Marcelo Santillan

Faculty Scholarship for the Cascadia Hazards Institute

Continuous GPS and seismic data from northern California show that slow earthquakes periodically rupture the Gorda-North America plate interface within southern Cascadia. On average, these creep events have occurred every 10.9±1.2 months since at least 1998. Appearing as week-long GPS extensional transients that reverse secular forearc contraction, the data show a recurrence interval 22% shorter than slow events recognized to the north. Seismic tremor here accompanies the GPS reversals, correlated across as many as 5 northern California seismometers. Tremor occurs sporadically throughout the year, but increases in duration and intensity by a factor of about 10 simultaneous with the GPS …

Precursory Transient Slip During The 2001 MW = 8.4 Peru Earthquake Sequence From Continuous Gps, Timothy I. Melbourne, Frank H. Webb

Faculty Scholarship for the Cascadia Hazards Institute

Two-hour position estimates from a continuous GPS station located at Arequipa, Peru, document precursory deformation beginning 18 hours prior to an M_w = 7.6 aftershock of the June 23rd 2001 M_w = 8.4 earthquake. This preseismic signal appears on the north and east components as a slow displacement with an amplitude twice that of the subsequent coseismic. Analysis of three years of 18-hour rate measurement shows this signal to be unprecedented and beyond four standard deviations from the mean rate. The best fitting centroid is directionally consistent with slow slip along the plate interface and suggests the preseismic …

Mantle Control Of Plate Boundary Deformation, Timothy I. Melbourne, Don Helmberger

Faculty Scholarship for the Cascadia Hazards Institute

The seismic wavefield propagating along the recently instrumented Pacific-North American plate boundary (California) displays remarkable variation, with regional shear waves arriving at coastal stations up to 20 seconds earlier than equidistant stations in eastern California. Broadband modeling of this data reveals that coastal paths sample fast upper mantle typical of Miocene-aged ocean plate (> 50 Km thickness). Inland paths sample slower uppermost mantle, with the seismic lithosphere, or lid, measuring less than 5 Km thick, characteristic of the Basin and Range extensional province. The boundary in the uppermost mantle between these provinces is sharp, expressing the juxtaposition of the stronger …

Surficial Slip Distribution On The Central Emerson Fault During The June 28, 1992, Landers Earthquake, California, Sally F. Mcgill, Charles M. Rubin

All Faculty Scholarship for the College of the Sciences

We present the results of our mapping of a 5.6‐km length of the central Emerson fault that ruptured during the 1992 Landers earthquake in the southwestern Mojave Desert, California. The right‐lateral slip along this portion of the rupture varied from about 150 to 530 cm along the main rupture zone. In some locations a total of up to 110 cm of additional right‐lateral slip occurred on secondary faults up to 1.7 km away from the main rupture zone. Other secondary faults carried up to several tens of centimeters of left‐lateral or thrust displacement. The maximum net vertical displacement was 175 …

Long Dormancy, Low Slip Rate, And Similar Slip‐Per‐Event For The Emerson Fault, Eastern California Shear Zone, Charles M. Rubin, Kerry Sieh

All Faculty Scholarship for the College of the Sciences

Excavations in a playa along the 1992 rupture of the Emerson fault reveal evidence of two paleoseismic events, with only one large prehistoric rupture in the past 15 millennia. Accelerator mass spectrometer radiocarbon ages of charcoal from playa sediments and from fault‐scarp colluvium directly beneath the playa beds indicate that the last large prehistoric slip event occurred about 9000 ka. Trench‐wall exposures revealed clear evidence of at least one pre‐9 ka rupture at the playa site. The event horizon of this earthquake is between two pedogenic carbonate layers that have radiocarbon ages of 14.8 ka and 24.1 ka, implying that …