Emma M. Hill

The 2010 Mw 7.8 Mentawai earthquake: Very shallow source of a rare tsunami earthquake determined from tsunami field survey and near-field GPS data

[1] The Mw 7.8 October 2010Mentawai, Indonesia, earthquake was a“tsunami earthquake,” a rare type of earthquake that generates a tsunami much larger than expected based on the seismicmagnitude.Itproducedalocallydevastatingtsunami,withrunupcommonlyinexcess of 6 m. We examine this event using a combination of high-rate GPS data, from instruments located on the nearby islands, and a tsunami field survey. The GPS displacement time series are deficient in high-frequency energy, and show small coseismic displacements ( 16 m. Our modeling results show that the combination of the small GPS displacements and large tsunami can only be explained by high fault slip at very shallow depths, far from the islands and close to the oceanic trench. Inelastic uplift of trench sediments likely contributed to the size of the tsunami. Recent results for the 2011 Mw 9.0 Tohoko-Oki earthquake have also shown shallow fault slip, but the results from our study, which involves a smaller earthquake, provide much stronger constraints on how shallow the rupture can be, with the majority of slip for the Mentawai earthquake occurring at depths of <6 km. This result challenges the conventional wisdom that the shallow tips of subduction megathrusts are aseismic, and therefore raises important questions both about the mechanical properties of the shallow fault zone and the potential seismic and tsunami hazard of this shallow region.

Impact of self-attraction and loading on the annual cycle in sea level

The annual exchange of water between the continents and oceans is observed by GPS, gravimetry, and altimetry. However, the global average amplitude of this annual cycle (observed amplitude of ∼8 mm) is not representative of the effects that would be observed at individual tide gauges or at ocean bottom pressure recorders because of self-attraction and loading effects (SAL). In this paper, we examine the spatial variation of sea level change caused by the three main components that load the Earth and contribute to the water cycle: hydrology (including snow), the atmosphere, and the dynamic ocean. The SAL effects cause annual amplitudes at tide gauges (modeled here with a global average of ∼9 mm) to vary from less than 2 mm to more than 18 mm. We find a variance reduction (global average of 3 to 4%) after removing the modeled time series from a global set of tide gauges. We conclude that SAL effects are significant in places (e.g., the south central Pacific and coastal regions in Southeast Asia and west central Africa) and should be considered when interpreting these data sets and using them to constrain ocean circulation models.

Rupture process of the 2010 Mw 7.8 Mentawai tsunami earthquake from joint inversion of near-field hr-GPS and teleseismic body wave recordings constrained by tsunami observations

The 25 October 2010 Mentawai tsunami earthquake (Mw 7.8) ruptured the shallow portion of the Sunda megathrust seaward of the Mentawai Islands, offshore of Sumatra, Indonesia, generating a strong tsunami that took 509 lives. The rupture zone was updip of those of the 12 September 2007 Mw 8.5 and 7.9 underthrusting earthquakes. High-rate (1 s sampling) GPS instruments of the Sumatra GPS Array network deployed on the Mentawai Islands and Sumatra mainland recorded time-varying and static ground displacements at epicentral distances from 49 to 322 km. Azimuthally distributed tsunami recordings from two deepwater sensors and two tide gauges that have local high-resolution bathymetric information provide additional constraints on the source process. Finite-fault rupture models, obtained by joint inversion of the high-rate (hr)-GPS time series and numerous teleseismic broadband P and S wave seismograms together with iterative forward modeling of the tsunami recordings, indicate rupture propagation ~50 km up dip and ~100 km northwest along strike from the hypocenter, with a rupture velocity of ~1.8 km/s. Subregions with large slip extend from 7 to 10 km depth ~80 km northwest from the hypocenter with a maximum slip of 8 m and from ~5 km depth to beneath thin horizontal sedimentary layers beyond the prism deformation front for ~100 km along strike, with a localized region having >15 m of slip. The seismic moment is 7.2 × 1020 N m. The rupture model indicates that local heterogeneities in the shallow megathrust can accumulate strain that allows some regions near the toe of accretionary prisms to fail in tsunami earthquakes.

Earthquake magnitude calculation without saturation from the scaling of peak ground displacement

GPS instruments are noninertial and directly measure displacements with respect to a global reference frame, while inertial sensors are affected by systematic offsets—primarily tilting—that adversely impact integration to displacement. We study the magnitude scaling properties of peak ground displacement (PGD) from high-rate GPS networks at near-source to regional distances (~10–1000 km), from earthquakes between Mw6 and 9. We conclude that real-time GPS seismic waveforms can be used to rapidly determine magnitude, typically within the first minute of rupture initiation and in many cases before the rupture is complete. While slower than earthquake early warning methods that rely on the first few seconds of P wave arrival, our approach does not suffer from the saturation effects experienced with seismic sensors at large magnitudes. Rapid magnitude estimation is useful for generating rapid earthquake source models, tsunami prediction, and ground motion studies that require accurate information on long-period displacements.

Characterization of site‐specific GPS errors using a short‐baseline network of braced monuments at Yucca Mountain, southern Nevada

[1] We use a short-baseline network of braced monuments to investigate site-specific GPS effects. The network has baseline lengths of 10, 100, and 1000 m. Baseline time series have root mean square (RMS) residuals, about a model for the seasonal cycle, of 0.05–0.24 mm for the horizontal components and 0.20–0.72 mm for the radial. Seasonal cycles occur, with amplitudes of 0.04–0.60 mm, even for the horizontal components and even for the shortest baselines. For many time series these lag seasonal cycles in local temperature measurements by 23–43 days. This could suggest that they are related to bedrock thermal expansion. Both shorter-period signals and seasonal cycles for shorter baselines to REP2, the one short-braced monument in our network, are correlated with temperature, with no lag time. Differences between REP2 and the other stations, which are deep-braced, should reflect processes occurring in the upper few meters of the ground. These correlations may be related to thermal expansion of these upper ground layers, and/or thermal expansion of the monuments themselves. Even over these short distances we see a systematic increase in RMS values with increasing baseline length. This, and the low RMS levels, suggests that site-specific effects are unlikely to be the limiting factor in the use of similar GPS sites for geophysical investigations.

The 2012 Mw 8.6 Wharton Basin sequence: A cascade of great earthquakes generated by near-orthogonal, young, oceanic mantle faults

We improve constraints on the slip distribution and geometry of faults involved in the complex, multisegment, Mw 8.6 April 2012 Wharton Basin earthquake sequence by joint inversion of high-rate GPS data from the Sumatran GPS Array (SuGAr), teleseismic observations, source time functions from broadband surface waves, and far-field static GPS displacements. This sequence occurred under the Indian Ocean, ∼400 km offshore Sumatra. The events are extraordinary for their unprecedented rupture of multiple cross faults, deep slip, large strike-slip magnitude, and potential role in the formation of a discrete plate boundary between the Indian and Australian plates. The SuGAr recorded static displacements of up to ∼22 cm, along with time-varying arrivals from the complex faulting, which indicate that the majority of moment release was on young, WNW trending, right-lateral faults, counter to initial expectations that an old, lithospheric, NNE trending fracture zone played the primary role. The new faults are optimally oriented to accommodate the present-day stress field. Not only was the greatest moment released on the younger faults, but it was these that sustained very deep slip and high stress drop (>20 MPa). The rupture may have extended to depths of up to 60 km, suggesting that the oceanic lithosphere in the northern Wharton Basin may be cold and strong enough to sustain brittle failure at such depths. Alternatively, the rupture may have occurred with an alternative weakening mechanism, such as thermal runaway.

Combination of geodetic observations and models for glacial isostatic adjustment fields in Fennoscandia

We demonstrate a new technique for using geodetic data to update a priori predictions for Glacial Isostatic Adjustment (GIA) in the Fennoscandia region. Global Positioning System (GPS), tide gauge, and Gravity Recovery and Climate Experiment (GRACE) gravity rates are assimilated into our model. The technique allows us to investigate the individual contributions from these data sets to the output GIA model in a self‐consistent manner. Another benefit of the technique is that we are able to estimate uncertainties for the output model. These are reduced with each data set assimilated. Any uncertainties in the GPS reference frame are absorbed by reference frame adjustments that are estimated as part of the assimilation. Our updated model shows a spatial pattern and magnitude of peak uplift that is consistent with previous models, but our location of peak uplift is slightly to the east of many of these. We also simultaneously estimate a spatially averaged rate of local sea level rise. This regional rate (∼1.5 mm/yr) is consistent for all solutions, regardless of which data sets are assimilated or the magnitude of a priori GPS reference frame constraints. However, this is only the case if a uniform regional gravity rate, probably representing errors in, or unmodeled contributions to, the low‐degree harmonic terms from GRACE, is also estimated for the assimilated GRACE data. Our estimated sea level rate is consistent with estimates obtained using a more traditional approach of direct “correction” using collocated GPS and tide gauge sites.

The mechanism of partial rupture of a locked megathrust: The role of fault morphology

Assessment of seismic hazard relies on estimates of how large an area of a tectonic fault could potentially rupture in a single earthquake. Vital information for these forecasts includes which areas of a fault are locked and how the fault is segmented. Much research has focused on exploring downdip limits to fault rupture from chemical and thermal boundaries, and along-strike barriers from subducted structural features, yet we regularly see only partial rupture of fully locked fault patches that could have ruptured as a whole in a larger earthquake. Here we draw insight into this conundrum from the 25 April 2015 M w 7.8 Gorkha (Nepal) earthquake. We invert geodetic data with a structural model of the Main Himalayan thrust in the region of Kathmandu, Nepal, showing that this event was generated by rupture of a decollement bounded on all sides by more steeply dipping ramps. The morphological bounds explain why the event ruptured only a small piece of a large fully locked seismic gap. We then use dynamic earthquake cycle modeling on the same fault geometry to reveal that such events are predicted by the physics. Depending on the earthquake history and the details of rupture dynamics, however, great earthquakes that rupture the entire seismogenic zone are also possible. These insights from Nepal should be applicable to understanding bounds on earthquake size on megathrusts worldwide.

An estimate of increases in storm surge risk to property from sea level rise in the first half of the twenty-first century.

Abstract Sea level is rising as the World Ocean warms and ice caps and glaciers melt. Published estimates based on data from satellite altimeters, beginning in late 1992, suggest that the global mean sea level has been rising on the order of 3 mm yr−1. Local processes, including ocean currents and land motions due to a variety of causes, modulate the global signal spatially and temporally. These local signals can be much larger than the global signal, and especially so on annual or shorter time scales. Even increases on the order of 10 cm in sea level can amplify the already devastating losses that occur when a hurricane-driven storm surge coincides with an astronomical high tide. To quantify the sensitivity of property risk to increasing sea level, changes in expected annual losses to property along the U.S. Gulf and East Coasts are calculated as follows. First, observed trends in sea level rise from tide gauges are extrapolated to the year 2030, and these changes are interpolated to all coastal location...

An ancient shallow slip event on the Mentawai segment of the Sunda megathrust, Sumatra

The outer-arc islands of western Sumatra rise during great megathrust earthquakes, due to large slip on the underlying megathrust. In contrast, the islands subsided up to a few centimeters during the recent tsunamigenic earthquake of October 2010, due to slip far updip, near the trench. Coral microatolls on one of the islands recorded a much larger subsidence, at least 35 cm, during an event in approximately A.D. 1314. We calculate a suite of slip models, slightly deeper and/or larger than the 2010 event, that are consistent with this large amount of subsidence. Sea level records from older coral microatolls suggest that these events occur at least once every millennium, but likely far less frequently than their great downdip neighbors. The revelation that shallow slip events are important contributors to the seismic cycle of the Mentawai segment further complicates our understanding of this subduction megathrust and our assessment of the regions exposure to seismic and tsunami hazards.