Jacob I. Walter

A tremor and slip event on the Cocos-Caribbean subduction zone as measured by a global positioning system (GPS) and seismic network on the Nicoya Peninsula, Costa Rica

In May 2007 a network of global positioning systems (GPS) and seismic stations on the Nicoya Peninsula, of northern Costa Rica, recorded a slow-slip event accompanied by seismic tremor. The close proximity of the Nicoya Peninsula to the seismogenic part of the Cocos-Caribbean subduction plate boundary makes it a good location to study such events. Several centimeters of southwest motion were recorded by the GPS stations over a period of several days to several weeks, and the seismic stations recorded three distinct episodes of tremor during the same time span. Inversion of the surface displacement data for the depth and pattern of slip on the plate interface shows peak slip at a depth of 25–30 km, downdip of the main seismogenic zone. Estimated temperatures here are ∼250°–300°C, lower than in other subduction zones where events of this nature have been previously identified. There may also be a shallower patch of slip at ∼6 km depth. These results are significant in that they are the first to suggest that slow slip can occur at the updip transition from stick slip to stable sliding, and that a critical temperature threshold is not required for slow slip. Tremor and low-frequency earthquake locations are more difficult to determine. Our results suggest they occur on or near the plate interface at the same depth range as the deep slow slip, but not spatially colocated.

Oceanic mechanical forcing of a marine-terminating Greenland glacier

Abstract Dynamics of marine-terminating major outlet glaciers are of high interest because of their potential for drawing down large areas of the Greenland ice sheet. We quantify short-term changes in ice flow speed and calving at a major West Greenland glacier and examine their relationship to the presence of the sea-ice melange and tidal stage. A field campaign at the terminus of Store Gletscher (70.40˚N, 50.55˚W) spanning the spring and summer of 2008 included four broadband seismometers, three time-lapse cameras, a tide gauge, an automatic weather station and an on-ice continuous GPS station. Sub-daily fluctuations in speed coincide with two modes of oceanic forcing: (1) the removal of the ice melange from the terminus front and (2) tidal fluctuations contributing to speed increases following ice melange removal. Tidal fluctuations in ice flow speed were observed 16km from the terminus and possibly extend further. Seismic records suggest that periods of intensive calving activity coincide with ice-flow acceleration following breakup of the melange in spring. A synchronous increase in speed at the front and clearing of the melange suggests that the melange directly resists ice flow. We estimate a buttressing stress (~30–60 kPa) due to the presence of the ice melange that is greater than expected from the range of observed tides, though an order of magnitude less than the driving stress.

Transient slip events from near-field seismic and geodetic data on a glacier fault, Whillans Ice Plain, West Antarctica

[1] Bidaily, tidally modulated stick‐slip speed‐ups of the Whillans Ice Plain (WIP) provide insight into glacier dynamics and failure at a naturally repeating fault asperity. We installed a network of continuously operating GPS receivers in 2007 and deployed on‐ice broadband seismometers during the austral summer of 2008 on Whillans Ice Stream (WIS), West Antarctica, and recorded 26 glacier speed‐up events. Previous work during the 2004 field season suggested that these speed‐ups initiate as failure of an asperity on or near “ice raft A” that triggers rupture across the entire WIP. Our results for 2008 locate the slip initiation farther to the south of this feature, closer to the grounding line and the southernmost extent of the Ross Ice Shelf. The initiation may be controlled by a discontinuity in basal boundary conditions at the suture between two ice streams. A strong correlation between the amplitude of seismic waves generated at the rupture front and the total slip achieved over the duration of the slip event (∼30 min) suggests slip‐predictable behavior, i.e., the ability to forecast the eventual slip based on the first minute of seismic radiation. Successive slip events propagate with different rupture speeds (100–300 m/s) that strongly correlate (R 2 = 0.73) with the interevent duration. In addition, the amount of slip achieved during each event appears to be correlated with the rupture speed. We use these observations to constrain basal shear stress to be 4 kPa by calculating conditions for basal freezing. Our observations yield information regarding mechanics and dynamics of ice streambeds at the scale of tens to hundreds of kilometers. Subglacial processes are notoriously difficult to constrain on these large scales, which are relevant to the understanding of regional and continental ice motion.

Multiscale postseismic behavior on a megathrust: The 2012 Nicoya earthquake, Costa Rica

The Nicoya Peninsula in northwest Costa Rica overlies a section of the subduction megathrust along the Middle America Trench. On 5 September 2012, a moment magnitude 7.6 megathrust earthquake occurred beneath a dense network of continuous GPS and seismic stations. Many of the GPS stations recorded the event at high rate, 1 Hz or better. We analyze the temporal and spatial evolution of surface deformation after the earthquake. Our results show that the main rupture was followed by significant afterslip within the first 3 h following the main event. The behavior of the surface displacement can be represented by relaxation processes with three characteristic times: 7, 70, and more than 400 days. We assume that the long relaxation time corresponds to viscoelastic relaxation and the intermediate relaxation time corresponds to afterslip on the main fault. The short relaxation time may represent a combination of rapid afterslip, poroelastic adjustment in the upper crust, or other processes. During the first few months that followed the earthquake, afterslip likely released a significant amount of slip deficit still present following the coseismic rupture, in particular updip of the rupture. Afterslip seems to be bounded updip by regions affected by slow slip events prior to the earthquake, suggesting that the two processes are influenced by different frictional properties.

Subglacial discharge at tidewater glaciers revealed by seismic tremor

Abstract Subglacial discharge influences glacier basal motion and erodes and redeposits sediment. At tidewater glacier termini, discharge drives submarine terminus melting, affects fjord circulation, and is a central component of proglacial marine ecosystems. However, our present inability to track subglacial discharge and its variability significantly hinders our understanding of these processes. Here we report observations of hourly to seasonal variations in 1.5–10 Hz seismic tremor that strongly correlate with subglacial discharge but not with basal motion, weather, or discrete icequakes. Our data demonstrate that vigorous discharge occurs from tidewater glaciers during summer, in spite of fast basal motion that could limit the formation of subglacial conduits, and then abates during winter. Furthermore, tremor observations and a melt model demonstrate that drainage efficiency of tidewater glaciers evolves seasonally. Glaciohydraulic tremor provides a means by which to quantify subglacial discharge variations and offers a promising window into otherwise obscured glacierized environments.

Detailed spatiotemporal evolution of microseismicity and repeating earthquakes following the 2012 Mw 7.6 Nicoya earthquake

We apply a waveform matching technique to obtain a detailed earthquake catalog around the rupture zone of the 5 September 2012 moment magnitude 7.6 Nicoya earthquake, with emphasis on its aftershock sequence. Starting from a preliminary catalog, we relocate ~7900 events using TomoDD to better quantify their spatiotemporal behavior. Relocated aftershocks are mostly clustered in two groups. The first is immediately above the major coseismic slip patch, partially overlapping with shallow afterslip. The second one is 50 km SE to the main shock nucleation point and near the terminus of coseismic rupture, in a zone that exhibited little resolvable afterslip. Using the relocated events as templates, we scan through the continuous recording from 29 June 2012 to 30 December 2012, detecting approximately 17 times more than template events. We find 190 aftershocks in the first half hour following the main shock, mostly along the plate interface. Later events become more scattered in location, showing moderate expansion in both along-trench and downdip directions. From the detected catalog we identify 53 repeating aftershock clusters with mean cross-correlation values larger than 0.9, and indistinguishably intracluster event locations, suggesting slip on the same fault patch. Most repeating clusters occurred within the first major aftershock group. Very few repeating clusters were found in the aftershock grouping along the southern edge of the Peninsula, which is not associated with substantial afterslip. Our observations suggest that loading from nearby afterslip along the plate interface drives spatiotemporal evolution of aftershocks just above the main shock rupture patch, while aftershocks in the SE group are to the SE of the observed updip afterslip and poorly constrained.

The Pawnee earthquake as a result of the interplay among injection, faults and foreshocks

The Pawnee M5.8 earthquake is the largest event in Oklahoma instrument recorded history. It occurred near the edge of active seismic zones, similar to other M5+ earthquakes since 2011. It ruptured a previously unmapped fault and triggered aftershocks along a complex conjugate fault system. With a high-resolution earthquake catalog, we observe propagating foreshocks leading to the mainshock within 0.5 km distance, suggesting existence of precursory aseismic slip. At approximately 100 days before the mainshock, two M ≥ 3.5 earthquakes occurred along a mapped fault that is conjugate to the mainshock fault. At about 40 days before, two earthquakes clusters started, with one M3 earthquake occurred two days before the mainshock. The three M ≥ 3 foreshocks all produced positive Coulomb stress at the mainshock hypocenter. These foreshock activities within the conjugate fault system are near-instantaneously responding to variations in injection rates at 95% confidence. The short time delay between injection and seismicity differs from both the hypothetical expected time scale of diffusion process and the long time delay observed in this region prior to 2016, suggesting a possible role of elastic stress transfer and critical stress state of the fault. Our results suggest that the Pawnee earthquake is a result of interplay among injection, tectonic faults, and foreshocks.

Subseasonal changes observed in subglacial channel pressure, size, and sediment transport

Water that pressurizes the base of glaciers and ice sheets enhances glacier velocities and modulates glacial erosion. Predicting ice flow and erosion therefore requires knowledge of subglacial channel evolution, which remains observationally limited. Here we demonstrate that detailed analysis of seismic ground motion caused by subglacial water flow at Mendenhall Glacier (Alaska) allows for continuous measurement of daily to subseasonal changes in basal water pressure gradient, channel size, and sediment transport. We observe intermittent subglacial water pressure gradient changes during the melt season, at odds with common assumptions of slowly varying, low-pressure channels. These observations indicate that changes in channel size do not keep pace with changes in discharge. This behavior strongly affects glacier dynamics and subglacial channel erosion at Mendenhall Glacier, where episodic periods of high water pressure gradients enhance glacier surface velocity and channel sediment transport by up to 30% and 50%, respectively. We expect the application of this framework to future seismic observations acquired at glaciers worldwide to improve our understanding of subglacial processes.

A Tremor and Slip Event on the Cocos-Caribbean Subduction zone as measured by a GPS and Seismic Network on the Nicoya Peninsula, Costa Rica

In May 2007 a network of global positioning systems (GPS) and seismic stations on the Nicoya Peninsula, of northern Costa Rica, recorded a slow-slip event accompanied by seismic tremor. The close proximity of the Nicoya Peninsula to the seismogenic part of the Cocos-Caribbean subduction plate boundary makes it a good location to study such events. Several centimeters of southwest motion were recorded by the GPS stations over a period of several days to several weeks, and the seismic stations recorded three distinct episodes of tremor during the same time span. Inversion of the surface displacement data for the depth and pattern of slip on the plate interface shows peak slip at a depth of 25–30 km, downdip of the main seismogenic zone. Estimated temperatures here are ∼250°–300°C, lower than in other subduction zones where events of this nature have been previously identified. There may also be a shallower patch of slip at ∼6 km depth. These results are significant in that they are the first to suggest that slow slip can occur at the updip transition from stick slip to stable sliding, and that a critical temperature threshold is not required for slow slip. Tremor and low-frequency earthquake locations are more difficult to determine. Our results suggest they occur on or near the plate interface at the same depth range as the deep slow slip, but not spatially colocated.

A Decade of Induced Slip on the Causative Fault of the 2015 Mw 4.0 Venus Earthquake, Northeast Johnson County, Texas

On 7 May 2015, a MW 4.0 earthquake occurred near Venus, northeast Johnson County, Texas, in an area of the Bend Arch-Fort Worth Basin that reports long-term, high-volume wastewater disposal and that has hosted felt earthquakes since 2009. In the weeks following the MW 4.0 earthquake, we deployed a local seismic network and purchased nearby active-source seismic reflection data to capture additional events, characterize the causative fault, and explore potential links between ongoing industry activity and seismicity. Hypocenter relocations of the resulting local earthquake catalog span ~4-6 km depth and indicate a fault striking ~230°, dipping to the west, consistent with a nodal plane of the MW 4.0 regional moment tensor. Fault plane solutions indicate normal faulting, with B-axes striking parallel to maximum horizontal compressive stress. Seismic reflection data image the reactivated basement fault penetrating the Ordovician disposal layer and Mississippian production layer, but not displacing post-Lower Pennsylvanian units. Template matching at regional seismic stations indicates that low magnitude earthquakes with similar waveforms began in April 2008, with increasing magnitude over time. Pressure data from five saltwater disposal wells within 5 km of the active fault indicate a disposal formation that is 0.9-4.8 MPa above hydrostatic. We suggest that the injection of 28,000,000 m3 of wastewater between 2006 and 2015 at these wells led to an increase in subsurface pore fluid pressure that contributed to inducing this long-lived earthquake sequence. The 2015 MW 4.0 event represents the largest event in the continuing evolution of slip on the causative fault.