Engineering topics
Sauber, Jeanne
Publications and source records attributed to Sauber, Jeanne.
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.
Gravitational Changes of the Earth's Free Oscillation From Earthquakes: Theory and Feasibility Study Using GRACE Inter‐Satellite Tracking
GRACE satellites have detected regional‐scale preseismic, coseismic, and postseismic gravity changes associated with great earthquakes during the GRACE era (2002‐2017). Earthquakes also excite global‐scale transient gravity changes associated with free oscillations that may be discerned for a few days. In this study, we examine such global gravity changes due to Earth's free oscillations and quantify how they affect GRACE measurements. We employ the normal mode formalism to synthesize the global gravity changes after the 2004 Sumatra earthquake and simulate the (gravitational) free oscillation signals manifested in the GRACE K‐band ranging (KBR) measurements. Using the Kaula orbit perturbation theory, we show how GRACE inter‐satellite distances are perturbed through a complex coupling of eigenfrequencies of the normal modes with the Earth's rotation rate and the GRACE satellites' orbital frequency. It is found that a few gravest normal modes can generate range‐rate perturbations as large as 0.2 μm/s, which are comparable to actual errors of GRACE KBR ranging and accelerometer instruments. Wavelet time‐frequency analysis of the GRACE KBR residual data in December 2004 reveals the existence of a significant transient signal after the 2004 Sumatra earthquake. This transient signal is characterized by a frequency of ~0.022 mHz that could be potentially associated with the largest excitation due to the “football” mode of the Earth's free oscillation. However, the results are also affected by low‐frequency noise of the GRACE accelerometers. Improved space‐borne gravitational instrumentation may open new opportunities to study the Earth's interior and earthquakes independently from global seismological analysis.
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 2–3×10(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.
Estimating Viscoelastic Deformation Due to Seasonal Loading
Scientists have been making summer-‐time geodetic measurements in south central Alaska for decades to estimate the rate at which a continental-‐ocean terrane is accreting to the North American continent. Southern Alaska has big earthquakes every century and large, rapidly changing glaciers. In the last decade, primarily as part of the EarthScope Plate Boundary Observatory project, continuous GPS measurements have recorded the response of sites such as the near-‐coastal geodetic site, AB35 to competing processes: uplift and movement to the northwest due to tectonic forces and the response of the solid Earth to seasonal and longer-‐term changes in the cryosphere (snow and ice) surrounding the site. Which process causes the largest displacements of the site? Figure 1 (Blewitt, Nevada Geodetic Lab, 2015) shows the Northward, Eastward, and Upward motion of AB35 between 2007 and 2015. The site is moving rapidly to the north and west reflecting the tectonic convergence of site toward interior Alaska but there is small wiggle on the North component reflecting seasonal displacements of the site associated with snow loading and unloading. However, the Up component, shows a large seasonal signal due to snow loading in the winter (down) and ice and snow melting in the warmer months (site goes up). Between 2007 and the present, the site position is slowly moving upward, due to tectonic forcing but probably associated with longer-‐ term ice melting as well. We are using the CIG finite element modeling (FEM) program Pylith to estimate the surface displacements and stresses associated with seasonal loading changes (top figure and Figure 2 far right) for water year 2012, 2011.8 - 2012.8) and the longer-‐term retreat of the surrounding glaciers.
Broadscale Postseismic Gravity Change Following the 2011 Tohoku-Oki Earthquake and Implication for Deformation by Viscoelastic Relaxation and Afterslip
The analysis of GRACE gravity data revealed post-seismic gravity increase by 6 micro-Gal over a 500 km scale within a couple of years after the 2011 Tohoku-Oki earthquake, which is nearly 40-50% of the co-seismic gravity change. It originates mostly from changes in the isotropic component corresponding to the M(sub rr) moment tensor element. The exponential decay with rapid change in a year and gradual change afterward is a characteristic temporal pattern. Both viscoelastic relaxation and afterslip models produce reasonable agreement with the GRACE free-air gravity observation, while their Bouguer gravity patterns and seafloor vertical deformations are distinctly different. The post-seismic gravity variation is best modeled by the bi-viscous relaxation with a transient and steady state viscosity of 10(exp 18) and 10(exp 19) Pa s, respectively, for the asthenosphere. Our calculated higher-resolution viscoelastic relaxation model, underlying the partially ruptured elastic lithosphere, yields the localized post-seismic subsidence above the hypocenter reported from the GPS-acoustic seafloor surveying.
Ice Surface Morphology and Flow on Malaspina Glacier, Alaska: Implications for Regional Tectonics in the Saint Elias Orogen
The Saint Elias Mountains in southern Alaska are located at a structural syntaxis where the coastal thrust and fold belt of the Fairweather plate boundary intersects thrust faults and folds generated by collision of the Yakutat Terrane. The axial trace of this syntaxis extends southeastward out of the Saint Elias Mountains and beneath Malaspina Glacier where it is hidden from view and cannot be mapped using conventional methods. Here we examine the surface morphology and flow patterns of Malaspina Glacier to infer characteristics of the bedrock topography and organization of the syntaxis. Faults and folds beneath the eastern part of the glacier trend northwest and reflect dextral transpression near the terminus of the Fairweather fault system. Those beneath the western part of the glacier trend northeast and accommodate folding and thrust faulting during collision and accretion of the Yakutat Terrane. Mapping the location and geometry of the structural syntaxis provides important constraints on spatial variations in seismicity, fault kinematics, and crustal shortening beneath Malaspina Glacier, as well as the position of the collisional deformation front within the Yakutat Terrane. We also speculate that the geometrical complexity of intersecting faults within the syntaxis formed a barrier to rupture propagation during two regional Mw 8.1earthquakes in September 1899.
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.
Source Parameter Inversion for Recent Great Earthquakes from a Decade-long Observation of Global Gravity Fields
We quantify gravity changes after great earthquakes present within the 10 year long time series of monthly Gravity Recovery and Climate Experiment (GRACE) gravity fields. Using spherical harmonic normal-mode formulation, the respective source parameters of moment tensor and double-couple were estimated. For the 2004 Sumatra-Andaman earthquake, the gravity data indicate a composite moment of 1.2x10(exp 23)Nm with a dip of 10deg, in agreement with the estimate obtained at ultralong seismic periods. For the 2010 Maule earthquake, the GRACE solutions range from 2.0 to 2.7x10(exp 22)Nm for dips of 12deg-24deg and centroid depths within the lower crust. For the 2011 Tohoku-Oki earthquake, the estimated scalar moments range from 4.1 to 6.1x10(exp 22)Nm, with dips of 9deg-19deg and centroid depths within the lower crust. For the 2012 Indian Ocean strike-slip earthquakes, the gravity data delineate a composite moment of 1.9x10(exp 22)Nm regardless of the centroid depth, comparing favorably with the total moment of the main ruptures and aftershocks. The smallest event we successfully analyzed with GRACE was the 2007 Bengkulu earthquake with M(sub 0) approx. 5.0x10(exp 21)Nm. We found that the gravity data constrain the focal mechanism with the centroid only within the upper and lower crustal layers for thrust events. Deeper sources (i.e., in the upper mantle) could not reproduce the gravity observation as the larger rigidity and bulk modulus at mantle depths inhibit the interior from changing its volume, thus reducing the negative gravity component. Focal mechanisms and seismic moments obtained in this study represent the behavior of the sources on temporal and spatial scales exceeding the seismic and geodetic spectrum.
Effects of Bedrock Lithology and Subglacial Till on the Motion of Ruth Glacier, Alaska, Deduced from Five Pulses from 1973-2012
A pulse is a type of unstable glacier flow intermediate between normal flow and surging. Using Landsat MSS, TM, and ETM+ imagery and feature tracking software, a time-series of mostly annual velocity maps from 1973 to 2012 was produced that reveals five pulses of Ruth Glacier, Alaska. Peaks in ice velocity were found in the 1981, 1989, 1997, 2003, and 2010; approximately every 7 years. During these peak years the ice velocity increased 300%, from approximately 40 m/yr to 160 m/yr, and occurred in an area of the glacier underlain by sedimentary bedrock. Based on the spatio-temporal behavior of Ruth Glacier during the pulse cycles, we suggest the pulses are due to enhanced basal motion via deformation of a subglacial till. The cyclical nature of the pulses is theorized to be due to a thin till, with low permeability, that causes incomplete drainage of the till between the pulses, followed by eventual recharge and dilation of the till. These findings suggest care is needed when attempting to correlate changes in regional climate with decadal-scale changes in velocity, because in some instances basal conditions may have a greater influence on ice dynamics than climate.
Plate Margin Deformation and Active Tectonics Along the Northern Edge of the Yakutat Terrane in the Saint Elias Orogen, Alaska and Yukon, Canada
The northwest directed motion of the Pacific plate is accompanied by migration and collision of the Yakutat terrane into the cusp of southern Alaska. The nature and magnitude of accretion and translation on upper crustal faults and folds is poorly constrained, however, due to pervasive glaciation. In this study we used high-resolution topography, geodetic imaging, seismic, and geologic data to advance understanding of the transition from strike-slip motion on the Fairweather fault to plate margin deformation on the Bagley fault, which cuts through the upper plate of the collisional suture above the subduction megathrust. The Fairweather fault terminates by oblique-extensional splay faulting within a structural syntaxis, allowing rapid tectonic upwelling of rocks driven by thrust faulting and crustal contraction. Plate motion is partly transferred from the Fairweather to the Bagley fault, which extends 125 km farther west as a dextral shear zone that is partly reactivated by reverse faulting. The Bagley fault dips steeply through the upper plate to intersect the subduction megathrust at depth, forming a narrow fault-bounded crustal sliver in the obliquely convergent plate margin. Since . 20 Ma the Bagley fault has accommodated more than 50 km of dextral shearing and several kilometers of reverse motion along its southern flank during terrane accretion. The fault is considered capable of generating earthquakes because it is linked to faults that generated large historic earthquakes, suitably oriented for reactivation in the contemporary stress field, and locally marked by seismicity. The fault may generate earthquakes of Mw <= 7.5.
Contribution of Satellite Gravimetry to Understanding Seismic Source Processes of the 2011 Tohoku-Oki Earthquake
The 2011 great Tohoku-Oki earthquake, apart from shaking the ground, perturbed the motions of satellites orbiting some hundreds km away above the ground, such as GRACE, due to coseismic change in the gravity field. Significant changes in inter-satellite distance were observed after the earthquake. These unconventional satellite measurements were inverted to examine the earthquake source processes from a radically different perspective that complements the analyses of seismic and geodetic ground recordings. We found the average slip located up-dip of the hypocenter but within the lower crust, as characterized by a limited range of bulk and shear moduli. The GRACE data constrained a group of earthquake source parameters that yield increasing dip (7-16 degrees plus or minus 2 degrees) and, simultaneously, decreasing moment magnitude (9.17-9.02 plus or minus 0.04) with increasing source depth (15-24 kilometers). The GRACE solution includes the cumulative moment released over a month and demonstrates a unique view of the long-wavelength gravimetric response to all mass redistribution processes associated with the dynamic rupture and short-term postseismic mechanisms to improve our understanding of the physics of megathrusts.
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.
Crustal Deformation and the Seismic Cycle across the Kodiak Islands, Alaska
The Kodiak Islands are located approx.120 to 250 km from the Alaska-Aleutian Trench and are within the southern extent of the 1964 Prince William Sound (M(sub W) = 9.2) earthquake rupture and aftershock zone. Here we report new campaign GPS results (1993-2001) from northeastern Kodiak and reprocessed GPS results (1993-1997) from southwestern Kodiak. The rate and orientation of the horizontal velocities, relative to a fixed North America, range from 29.7 +/- 1.7 mm/yr at N30.3degW +/- 3.3deg, located approx.120 km from the deepest point of the trench, to 8.0 +/- 1.3 mm/yr at N62.4degW +/- 9.3deg, located approx.230 km from the trench. We evaluated alternate models of coseismic and interseismic slip to test the importance of the mechanisms that account for surface deformation rates. Near the Gulf of Alaska coastal region of Kodiak the horizontal velocity can be accounted for primarily by the viscoelastic response to plate motion and a locked main thrust zone (MTZ), down-dip creep, and to a lesser extent, slip in the 1964 earthquake. Further inland the dominant mechanisms that account for post-1964 uplift rates are time-dependent, down-dip creep and a locked MTZ; for the horizontal velocity component southwest translation of western Kodiak may be important as well. Based on the pre-1964 and post-1964 earthquake pattern of interseismic earthquakes, we suggest that between the occurrences of great earthquakes like the 1964 event, more moderate to large earthquakes occur in the southwestern Kodiak region than near northeastern Kodiak .