Engineering topics
Sandwell, David
Publications and source records attributed to Sandwell, David.
Grounding Zone and Tidal Response of the Amery Ice Shelf, East Antarctica
This report summarizes the main findings of the research project. Unfortunately, it turned out that there was not a great deal of SAR data over the Amery Ice Shelf that we were able to work with on the project; nevertheless, we did make considerable progress on this project, with both the existing SAR data and new field measurements that were collected under this grant. In total we had constructed two SAR interferograms (SSIs), and four SSIs. The latter were combined them to construct two differential SAR interferograms (DSIs;). DSIs are useful because the contribution to the SAR phase from horizontal ice motion is eliminated, since the time difference between the first and second pass within both image pairs used to make the DSI is the same for each pair. The SSIs and DSIs have revealed several interesting glaciological features, and have added to our knowledge of the Amery Ice Shelf (AIS).
Coulomb Stress Accumulation along the San Andreas Fault System
Stress accumulation rates along the primary segments of the San Andreas Fault system are computed using a three-dimensional (3-D) elastic half-space model with realistic fault geometry. The model is developed in the Fourier domain by solving for the response of an elastic half-space due to a point vector body force and analytically integrating the force from a locking depth to infinite depth. This approach is then applied to the San Andreas Fault system using published slip rates along 18 major fault strands of the fault zone. GPS-derived horizontal velocity measurements spanning the entire 1700 x 200 km region are then used to solve for apparent locking depth along each primary fault segment. This simple model fits remarkably well (2.43 mm/yr RMS misfit), although some discrepancies occur in the Eastern California Shear Zone. The model also predicts vertical uplift and subsidence rates that are in agreement with independent geologic and geodetic estimates. In addition, shear and normal stresses along the major fault strands are used to compute Coulomb stress accumulation rate. As a result, we find earthquake recurrence intervals along the San Andreas Fault system to be inversely proportional to Coulomb stress accumulation rate, in agreement with typical coseismic stress drops of 1 - 10 MPa. This 3-D deformation model can ultimately be extended to include both time-dependent forcing and viscoelastic response.
Fault Creep along the Southern San Andreas from Interferometric Synthetic Aperture Radar, Permanent Scatterers, and Stacking
Interferometric synthetic aperture radar (InSAR) provides a practical means of mapping creep along major strike-slip faults. The small amplitude of the creep signal (less than 10 mm/yr), combined with its short wavelength, makes it difficult to extract from long time span interferograms, especially in agricultural or heavily vegetated areas. We utilize two approaches to extract the fault creep signal from 37 ERS SAR images along the southem San Andreas Fault. First, amplitude stacking is utilized to identify permanent scatterers, which are then used to weight the interferogram prior to spatial filtering. This weighting improves correlation and also provides a mask for poorly correlated areas. Second, the unwrapped phase is stacked to reduce tropospheric and other short-wavelength noise. This combined processing enables us to recover the near-field (approximately 200 m) slip signal across the fault due to shallow creep. Displacement maps fiom 60 interferograms reveal a diffuse secular strain buildup, punctuated by localized interseismic creep of 4-6 mm/yr line of sight (LOS, 12-18 mm/yr horizontal). With the exception of Durmid Hill, this entire segment of the southern San Andreas experienced right-lateral triggered slip of up to 10 cm during the 3.5-year period spanning the 1992 Landers earthquake. The deformation change following the 1999 Hector Mine earthquake was much smaller (4 cm) and broader than for the Landers event. Profiles across the fault during the interseismic phase show peak-to-trough amplitude ranging from 15 to 25 mm/yr (horizontal component) and the minimum misfit models show a range of creeping/locking depth values that fit the data.
The 1999 (Mw 7.1) Hector Mine, California, Earthquake: Near-Field Postseismic Deformation from ERS Interferometry
Interferometric synthetic aperture radar (InSAR) data over the area of the Hector Mine earthquake (Mw 7.1, 16 October 1999) reveal postseismic deformation of several centimeters over a spatial scale of 0.5 to 50 km. We analyzed seven SAR acquisitions to form interferograms over four time periods after the event. The main deformations seen in the line-of-sight (LOS) displacement maps are a region of subsidence (60 mm LOS increase) on the northern end of the fault, a region of uplift (45 mm LOS decrease) located to the northeast of the primary fault bend, and a linear trough running along the main rupture having a depth of up to 15 mm and a width of about 2 km. We correlate these features with a double left-bending, rightlateral, strike-slip fault that exhibits contraction on the restraining side and extension along the releasing side of the fault bends. The temporal variations in the near-fault postseismic deformation are consistent with a characteristic time scale of 135 + 42 or - 25 days, which is similar to the relaxation times following the 1992 Landers earthquake. High gradients in the LOS displacements occur on the fault trace, consistent with afterslip on the earthquake rupture. We derive an afterslip model by inverting the LOS data from both the ascending and descending orbits. Our model indicates that much of the afterslip occurs at depths of less than 3 to 4 km.
Bounds on Lithospheric Thickness on Venus from Magellan Gravity and Topography Data
The primary objective of the work executed under NAGW-4784 is to provide constraints on the thermal and tectonic evolution of Venus. Establishing thermal and tectonic evolution models requires not only geological, but geophysical constraints, in particular the nature of temporal and spatial variations in crustal and lithospheric thickness. The major topics of study completed under NAGW-4784 (described more fully below) are: (1) detailed analyses of the resolution of Magellan Line-Of-Site (LOS) Doppler data to establish the minimum resolvable wavelength in the gravity data; (2) calculations of the global strain field in the venusian lithosphere and comparisons with global strain patterns from geological mapping; (3) study of the geological history of coronae at E. Eistla Regio; (4) estimation of crustal and lithospheric thickness by modeling of topography at asymmetric and symmetric rift-like chasmata; (5) preliminary investigations of spatial versus temporal variations in lithospheric thickness. Both the PI and Co-I have presented papers based on these topics at national and international meetings (American Geophysical Union Meetings, Lunar and Planetary Science Conferences, Chapman Conference on the Geodynamics of Venus).
Comparison of marine gravity from shipboard and high-density satellite altimetry along the Mid-Atlantic Ridge, 30.5-35.5 deg S
We compare new marine gravity fields derived from satellite altimetry with shipboard measurements over a region of more than 120,000 sq km in the central South Atlantic. Newly declassified satellite data were employed to construct free-air anomaly maps on 0.05 degree grids. An extensive gravity and bathymetry data set from four cruises along the Mid-Atlantic Ridge from 30.5-35.5 deg S provides a benchmark for testing the 2D resolution and accuracy of the satellite measurements where their crosstrack spacing is near their widest. The satellite gravity signal is coherent with bathymetry in this region down to wavelengths of 26 km, compared to 12.5 km for shipboard gravity. Residuals between the shipboard and satellite data sets have a roughly normal distribution. The standard deviation of satellite gravity with respect to shipboard measurements is nearly 7 mGal in a region of 140 mGal total variation, whereas the internal standard deviation at crossovers for GPS-navigated shipboard data is 1.8 mGal. The differences between shipboard and satellite data are too large to use satellite gravity to determine crustal thickness variations within a typical ridge segment.
Features on Venus generated by plate boundary processes
Various observations suggest that there are processes on Venus that produce features similar to those associated with plate boundaries on earth. Synthetic aperture radar images of Venus, taken with a radar whose wavelength is 12.6 cm, are compared with GLORIA images of active plate boundaries, obtained with a sound source whose wavelength is 23 cm. Features similar to transform faults and to abyssal hills on slow and fast spreading ridges can be recognized within the Artemis region of Venus but are not clearly visible elsewhere. The composition of the basalts measured by the Venera 13 and 14 and the Vega 2 spacecraft corresponds to that expected from adiabatic decompression, like that which occurs beneath spreading ridges on earth. Structures that resemble trenches are widespread on Venus and show the same curvature and asymmetry as they do on earth. These observations suggest that the same simple geophysical models that have been so successfully used to understand the tectonics of earth can also be applied to Venus.