Seismic Rate Changes Associated with Seasonal, Annual, and Decadal Changes in the Cryosphere

Near the Bering Glacier Global Fiducial site in southern Alaska large cryospheric fluctuations occur in a region of upper crustal faulting and folding associated with collision and accretion of the Yakutat terrane. In this study we report constraints on seasonal, annual and decadal cryospheric changes estimated over the last decade from field, aircraft and satellite measurements, and we evaluate the influence of cryospheric changes on the background seismic rate. Multi-year images from the Bering Glacier global fiducial site are available since mid-2003 to constrain changes in extent of the Bering Glacier and to discern feature changes in the glacial surface. Starting around the same time, satellite gravimetric measurements from the Gravity Recovery and Climate experiment (GRACE) commenced. Large spatial-scale mass change calculated from the GRACE 1deg x 1deg mascon solution of Luthcke et al. [2012] indicate a general trend of annual ice mass loss for southern Alaska but with large, variable seasonal mass fluctuations. Since 2007, the station position of a continuous GPS site near Cape Yakataga (Alaska EarthScope PBO site, AB35) has been available as well. In addition to changes in the geodetic position due to tectonic motion, this GPS station shows large seasonal excursions in the detrended vertical and horizontal position components consistent with snow loading in the fall and winter and melt onset/mass decrease in the spring/summer. To better understand the timing of processes responsible for the onset of cryospheric mass loss documented in the GRACE data, we examined changes in the snow cover extent and the onset of melt in the spring. We calculated the surface displacements of the solid Earth and theoretical earthquake failure criteria associated with these annual and seasonal ice and snow changes using layered elastic half-space. Additionally, we compared the seismic rate (M>1.8) from a reference background time period against other time periods with variable ice or tectonic change characteristics to test the significance of seismic rate changes. Our earlier results suggest statistically significant changes in the background seismic rate associated with large seasonal mass changes. INDE

Geodetic measurement of deformation in California

The very long baseline interferometry (VLBI) measurements made in the western U.S. since 1979 as part of the NASA Crustal Dynamics Project provide discrete samples of the temporal and spatial deformation field. The interpretation of the VLBI-derived rates of deformation requires an examination of geologic information and more densely sampled ground-based geodetic data. In the first two of three related studies embodying this thesis triangulation and trilateration data measured on two regional networks are processed, one in the central Mojave Desert and one in the Coast Ranges east of the San Andreas fault. At the spatial scales spanned by these local geodetic networks, auxiliary geologic and geophysical data have been utilized to examine the relation between measured incremental strain and the accommodation of strain seen in local geological structures, strain release in earthquakes, and principal stress directions inferred from in situ measurements. In the third study, VLBI data from stations distributed across the Pacific - North American plate boundary zone in the western United States are processed. The VLBI data have been used to constrain the integrated rate of deformation across portions of the continental plate boundary in California and to provide a tectonic framework to interpret regional geodetic and geologic studies

Geodetic Measurement of Deformation East of the San Andreas Fault in Central California

Triangulation and trilateration data from two geodetic networks located between the western edge of the Great Valley and the San Andreas fault have been used to calculate shear strain rates in the Diablo Range and to estimate the slip rate along the Calaveras and Paicines faults in Central California. Within the Diablo Range the average shear strain rate was determined for the time period between 1962 and 1982 to be 0.15 + or - 0.08 microrad/yr, with the orientation of the most compressive strain at N 16 deg E + or - 14 deg. The orientation of the principal compressive strain predicted from the azimuth of the major structures in the region is N 25 deg E. It is inferred that the measured strain is due to compression across the folds of this area: the average shear straining corresponds to a relative shortening rate of 4.5 + or - 2.4 mm/yr. From an examination of wellbore breakout orientations and the azimuths of P-axes from earthquake focal mechanisms the inferred orientation of maximum compressive stress was found to be similar to the direction of maximum compressive strain implied by the trend of local fold structures. Results do not support the hypothesis of uniform fault-normal compression within the Coast Ranges. From trilateration measurements made between 1972 and 1987 on lines that are within 10 km of the San Andreas fault, a slip rate of 10 to 12 mm/yr was calculated for the Calaveras-Paicines fault south of Hollister. The slip rate of the Paicines fault decreases to 4 mm/yr near Bitter

A Joint Analysis of GPS Displacement and GRACE Geopotential Data for Simultaneous Estimation of Geocenter Motion and Gravitational Field

Gravitational potential data from GRACE are being used to study mass redistribution within and between the atmosphere, hydrosphere, cryosphere, and solid Earth. The GRACE data are made available in a reference frame with its origin at the center of mass of the Earth system (geocenter) while many other geophysical models and data sets refer to a reference frame attached to the Earth's surface. Changes in the offset between these reference frames (geocenter motion) must be accounted for when GRACE data are used to quantify surface mass changes. In this study, we developed a technique for coestimation of geocenter motion and gravitational potential field seamlessly from degree 1 to 90 by simultaneously inverting a set of globallydistributed GPS displacement time series and the temporallyvarying GRACE gravity data. We found that the effect of geocenter motion was evident particularly in the GPS time series of horizontal displacements. Our estimates of geocenter motion are most consistent with the Satellite Laser Ranging (SLR) results within 1 mm in X and Z components and a submillimeter in Y component, when compared to monthly variability averaged over the period of 20032016. The overall magnitude of the degree1 (l = 1) surface mass load is estimated to be ~3 cm in equivalent water height annually migrating southwestward from Europe (DecemberJanuary) to the South Pacific (JuneJuly). Our results also show that dense GPS network data improve water storage recovery in major river basins in the United States and Europe by contributing significantly to the recovery of higherdegree (l ~20) geopotential coefficients

Geodetic measurement of deformation in California

Thesis (Ph. D.)--Massachusetts Institute of Technology, Dept. of Earth, Atmospheric, and Planetary Sciences, 1989.Includes bibliographical references (leaves 201-219).by Jeanne Marie Sauber.Ph.D

Seismicity near Palmdale, California, and its relation to strain changes

We evaluate the relationships between the spatio-temporal patterns and faulting mechanisms of small earthquakes and the recent temporal changes in horizontal strain observed along the ‘big bend’ portion of the San Andreas fault near Palmdale, California. Microearthquake activity along the entire big bend of the San Andreas fault increased in November 1976 concurrent with the initiation of an earthquake swarm at Juniper Hills. This activity then decreased abruptly to the northwest and southeast of Juniper Hills during the beginning of 1979. This drop in seismic activity occurred around the time that crustal dilatation was observed on the U. S. Geological Survey Palmdale trilateration network. Focal mechanisms from the study region are predominantly thrust. There are two time periods when the mechanisms are closer to strike slip than to thrust. The first period (December 1976 to February 1977) corresponds to the beginning of the Juniper Hills swarm. The second period (November 1978 to April 1979) approximately coincides with a change in trend of the strain data from uniaxial N-S compression to dilatation

Sea Level Rise in the Samoan Islands Escalated by Viscoelastic Relaxation After the 2009 Samoa-Tonga Earthquake

The Samoan islands are an archipelago hosting a quarter million people mostly residing inthree major islands, Savai'i and Upolu (Samoa), and Tutuila (American Samoa). The islands haveexperienced sea level rise by 2-3 mm/year during the last half century. The rate, however, has dramaticallyincreased following the Mw 8.1 SamoaTonga earthquake doublet (megathrust + normal faulting) inSeptember 2009. Since the earthquake, we found largescale gravity increase (0.5 Gal/year) around theislands and ongoing subsidence (8-16 mm/year) of the islands from our analysis of Gravity Recovery AndClimate Experiment gravity and GPS displacement data. The postseismic horizontal displacement is faster inSamoa, while the postseismic subsidence rate is considerably larger in American Samoa. The analysis oflocal tide gauge records and satellite altimeter data also identified that the relative sea level rise becomesfaster by 7-9 mm/year in American Samoa than Samoa. A simple viscoelastic model with a Maxwellviscosity of 2310(exp 18) Pa s for the asthenosphere explained postseismic deformation at nearby GPS sites aswell as Gravity Recovery And Climate Experiment gravity change. It is found that the constructiveinterference of viscoelastic relaxation from both megathrust and normal faulting has intensified thepostseismic subsidence at American Samoa, causing ~5 times faster sea level rise than the global average.Our model indicates that this trend is likely to continue for decades and result in sea level rise of 30-40 cm,which is independent of and in addition to anticipated climaterelated sea level rise. It will worsen coastalflooding on the islands leading to regular nuisance flooding

The Propagation of a Surge Front on Bering Glacier, Alaska, 2001-2011

Bering Glacier, Alaska, USA, has a 20 year surge cycle, with its most recent surge reaching the terminus in 2011. To study this most recent activity a time series of ice velocity maps was produced by applying optical feature-tracking methods to Landsat-7 ETM+ imagery spanning 2001-11. The velocity maps show a yearly increase in ice surface velocity associated with the down-glacier movement of a surge front. In 2008/09 the maximum ice surface velocity was 1.5 plus or minus 0.017 kilometers per a in the mid-ablation zone, which decreased to 1.2 plus or minus 0.015 kilometers per a in 2009/10 in the lower ablation zone, and then increased to nearly 4.4 plus or minus 0.03 kilometers per a in summer 2011 when the surge front reached the glacier terminus. The surge front propagated down-glacier as a kinematic wave at an average rate of 4.4 plus or minus 2.0 kilometers per a between September 2002 and April 2009, then accelerated to 13.9 plus or minus 2.0 kilometers per a as it entered the piedmont lobe between April 2009 and September 2010. Thewave seems to have initiated near the confluence of Bering Glacier and Bagley Ice Valley as early as 2001, and the surge was triggered in 2008 further down-glacier in the mid-ablation zone after the wave passed an ice reservoir area

DESDynI Lidar for Solid Earth Applications

As part of the NASA's DESDynI mission, global elevation profiles from contiguous 25 m footprint Lidar measurements will be made. Here we present results of a performance simulation of a single pass of the multi-beam Lidar instrument over uplifted marine terraces in southern Alaska. The significance of the Lidar simulations is that surface topography would be captured at sufficient resolution for mapping uplifted terraces features but it will be hard to discern I-2m topographic change over features less than tens of meters in width. Since Lidar would penetrate most vegetation, the accurate bald Earth elevation profiles will give new elevation information beyond the standard 30-m OEM

A Joint Analysis of GPS Displacement and GRACE Geopotential Data for Simultaneous Estimation of Geocenter Motion and Gravitational Field

Gravitational potential data from GRACE are being used to study mass redistribution within and between the atmosphere, hydrosphere, cryosphere, and solid Earth. The GRACE data are made available in a reference frame with its origin at the center of mass of the Earth system (geocenter) while many other geophysical models and data sets refer to a reference frame attached to the Earth's surface. Changes in the offset between these reference frames (geocenter motion) must be accounted for when GRACE data are used to quantify surface mass changes. In this study, we developed a technique for co‐estimation of geocenter motion and gravitational potential field seamlessly from degree 1 to 90 by simultaneously inverting a set of globally‐distributed GPS displacement time series and the temporally‐varying GRACE gravity data. We found that the effect of geocenter motion was evident particularly in the GPS time series of horizontal displacements. Our estimates of geocenter motion are most consistent with the Satellite Laser Ranging (SLR) results within 1 mm in X and Z components and a submillimeter in Y component, when compared to monthly variability averaged over the period of 2003–2016. The overall magnitude of the degree‐1 (l = 1) surface mass load is estimated to be ~3 cm in equivalent water height annually migrating south‐westward from Europe (December–January) to the South Pacific (June–July). Our results also show that dense GPS network data improve water storage recovery in major river basins in the United States and Europe by contributing significantly to the recovery of higher‐degree (l ≥ ~20) geopotential coefficients.This work is funded by The University of Newcastle to support NASA's GRACE and GRACE Follow‐On projects as an international science team member to the missions and by Australian Research Council (DP170100224) as well as NASA support to J. Sauber as a science team membe