Export 168 results:
Sort by: Author Title Type [ Year  (Desc)]
Kilb, D, Keen CS, Newman RL, Kent GM, Sandwell DT, Vernon FL, Johnson CL, Orcutt JA.  2003.  The Visualization Center at Scripps Institution of Oceanography: Education and Outreach. Seismological Research Letters. 74:641-648. AbstractWebsite

The immersive environment of the Visualization Center at Scripps, coupled with the presentation of current seismological research, has great education and outreach potential. Since its March 2002 opening, the Visualization Center at Scripps has had more than 2,500 visitors, and numerous virtual visitors have explored our Web pages, which include streaming QuickTime movies of geophysical data, tutorials on how to use SGI/Mac/Windows registered visualization software, and examples of visualizations developed by SIO researchers and faculty members ( We will continue to expand the use of the Visualization Center at Scripps for K-12 and informal education, and to use the center to supply geophysical data sets, movies, and research results to as large a group of educators as possible. Our goal is to develop more sustained K-12 educational programs and to generate assessments of the center's programs and the educational products created at the Center.

Fialko, Y, Sandwell D, Agnew D, Simons M, Shearer P, Minster B.  2002.  Deformation on nearby faults induced by the 1999 Hector Mine earthquake. Science. 297:1858-1862.   10.1126/science.1074671   AbstractWebsite

Interferometric Synthetic Aperture Radar observations of surface deformation due to the 1999 Hector Mine earthquake reveal motion on several nearby faults of the eastern California shear zone. We document both vertical and horizontal displacements of several millimeters to several centimeters across kilometer-wide zones centered on pre-existing faults. Portions of some faults experienced retrograde (that is, opposite to their long-term geologic slip) motion during or shortly after the earthquake. The observed deformation likely represents elastic response of compliant fault zones to the permanent co-seismic stress changes. The induced fault displacements imply decreases in the effective shear modulus within the kilometer-wide fault zones, indicating that the latter are mechanically distinct from the ambient crustal rocks.

Lyons, SN, Bock Y, Sandwell DT.  2002.  Creep along the imperial fault, southern California, from GPS measurements. Journal of Geophysical Research-Solid Earth. 107   10.1029/2001jb000763   AbstractWebsite

[1] In May of 1999 and 2000, we surveyed with Global Positioning System (GPS) 46 geodetic monuments established by Imperial College, London, in a dense grid (half-mile spacing) along the Imperial Fault, with three additional National Geodetic Survey sites serving as base stations. These stations were previously surveyed in 1991 and 1993. The Imperial College sites were surveyed in rapid-static mode (15-20 min occupations), while the NGS sites continuously received data for 10 h d(-1). Site locations were calculated using the method of instantaneous positioning, and velocities were determined relative to one of the NGS base stations. Combining our results with far-field velocities from the Southern California Earthquake Center (SCEC), we fit the data to a simple elastic dislocation model with 35 mm yr(-1) of right-lateral slip below 10 km and 9 mm yr(-1) of creep from the surface down to 3 km. The velocity field is asymmetrical across the fault and could indicate a dipping fault plane to the northeast or a viscosity contrast across the fault.

Jacobs, A, Sandwell D, Fialko Y, Sichoix L.  2002.  The 1999 (M-w 7. 1) Hector Mine, California, earthquake: Near-field postseismic deformation from ERS interferometry. Bulletin of the Seismological Society of America. 92:1433-1442.   10.1785/0120000908   AbstractWebsite

Interferometric synthetic aperture radar (InSAR) data over the area of the Hector Mine earthquake (M-w 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, right-lateral, 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.

Sandwell, DT, Sichoix L, Smith B.  2002.  The 1999 Hector Mine earthquake, southern California: Vector near-field displacements from ERS InSAR. Bulletin of the Seismological Society of America. 92:1341-1354.   10.1785/0120000901   AbstractWebsite

Two components of fault slip are uniquely determined from two line-of-sight (LOS) radar interferograms by assuming that the fault-normal component of displacement is zero. We use this approach with ascending and descending interferograms from the ERS satellites to estimate surface slip along the Hector Mine earthquake rupture. The LOS displacement is determined by visually counting fringes to within I kin of the outboard ruptures. These LOS estimates and uncertainties are then transformed into strike- and dip-slip estimates and uncertainties; the transformation is singular for a N-S oriented fault and optimal for an E-W oriented fault. In contrast to our previous strike-slip estimates, which were based only on a descending interferogram, we now find good agreement with the geological measurements, except at the ends of the rupture. The ascending interferogram reveals significant west-side-down dip-slip (similar to1.0 in) which reduces the strike-slip estimates by I to 2 in, especially along the northern half of the rupture. A spike in the strike-slip displacement of 6 m is observed in central part of the rupture. This large offset is confirmed by subpixel cross correlation of features in the before and after amplitude images. In addition to strike slip and dip slip, we identify uplift and subsidence along the fault, related to the restraining and releasing bends in the fault trace, respectively. Our main conclusion is that at least two look directions are required for accurate estimates of surface slip even along a pure strike-slip fault. Models and results based only on a single look direction could have major errors. Our new estimates of strike slip and dip slip along the rupture provide a boundary condition for dislocation modeling. A simple model, which has uniform slip to a depth of 12 km, shows good agreement with the observed ascending and descending interferograms.

Mellors, RJ, Sichoix L, Sandwell DT.  2002.  Lack of precursory slip to the 1999 Hector Mine, California, earthquake as constrained by InSAR. Bulletin of the Seismological Society of America. 92:1443-1449.   10.1785/0120010244   AbstractWebsite

We looked for evidence of interseismic strain occurring between the 1992 Landers earthquake and the 1999 Hector Mine earthquake near the Lavic Lake and Bullion faults by using interferometric synthetic aperture radar (InSAR). Interferograms covering the Hector Mine epicentral region were studied for possible slip along the Bullion and Lavic Lake faults by both visual inspection and a matched filter technique intended to emphasize slip located at the nucleation point. Some indications of possible deformation associated with the 5 July 1992 M-L 5.4 Pisgah event was observed, but high decorrelation prevented a conclusive determination. We have seen no evidence for precursory slip in the epicentral region up to 30 days before the Hector Mine event. We estimated that the slip equivalent to a M-w 4.5 event would have been observable in the months before the Hector Mine event, and this places an upper bound on the long-term precursory slip, had it occurred. We have noted that InSAR is well suited for detecting precursory slip in general due to the high spatial resolution and the lack of ground instrumentation required but that the detection level depends on the depth and orientation of the slip.

Baer, G, Schattner U, Wachs D, Sandwell D, Wdowinski S, Frydman S.  2002.  The lowest place on Earth is subsiding - An InSAR (interferometric synthetic aperture radar) perspective. Geological Society of America Bulletin. 114:12-23.   10.1130/0016-7606(2002)114<0012:tlpoei>;2   AbstractWebsite

Since the early 1990s, sinkholes and wide, shallow subsidence features (WSSFs) have become major problems along the Dead Sea shores in Israel and Jordan. Sinkholes are readily observed in the field, but their locations and timing are unpredictable. WSSFs are often difficult to observe in the field. However, once identified, they delineate zones of instability and increasing hazard. In this study we identify, characterize, and measure rates of subsidence along the Dead Sea shores by the interferometric synthetic aperture radar (InSAR) technique. We analyze 16 SAR scenes acquired during the years 1992 to 1999 by the European Remote Sensing ERS-1 and ERS-2 satellites. The interferograms span periods of between 2 and 71 months. WSSFs are observed in the Lisan Peninsula and along the Dead Sea shores, in a variety of appearances, including circular and elongate coastal depressions (a few hundred meters to a few kilometers in length), depressions in ancient alluvial fans, and depressions along salt-diapir margins. Phase differences measured in our interferograms correspond to subsidence rates generally in the range of 0-20 mm/yr within the studied period, with exceptional high rates that exceed 60 mm/yr in two specific regions. During the study period, the level of the Dead Sea and of the associated ground water has dropped by similar to6 m. This water-level drop within an aquifer overlying fine-grained, marly layers, would be expected to have caused aquifer-system consolidation resulting in gradual subsidence. Comparison of our InSAR observations with calculations of the expected consolidation shows that in areas where marl layers are known to compose part of the upper 30 m of the profile, estimated consolidation settlements are of the order of the measured subsidence. Our observations also show that in certain locations, subsidence appears to be structurally controlled by faults, seaward landslides, and salt domes. Gradual subsidence is unlikely to be directly related to the sinkholes, excluding the use of the WSSFs features as predictable precursors to sinkhole formation.

Watson, KM, Bock Y, Sandwell DT.  2002.  Satellite interferometric observations of displacements associated with seasonal groundwater in the Los Angeles basin. Journal of Geophysical Research-Solid Earth. 107   10.1029/2001jb000470   AbstractWebsite

[1] The Newport-Inglewood fault zone (NIFZ) displays interferometric synthetic aperture radar (SAR) phase features along most of its length having amplitudes of up to 60 mm. However, interpretation in terms of right-lateral, shallow slip along the fault fails to match the range of geologic estimates of slip. Recently, Bawden et al. [2001] proposed that these phase features, as well as a broader deformation pattern in the Los Angeles basin, are due to vertical motion related to annual variations in the elevation of the water table. We confirm this hypothesis through the analysis of a longer span of data consisting of 26 SAR images collected by the ERS-1 and ERS-2 spacecraft between June 1992 and June 2000. Moreover, we use continuous GPS measurements from 1995 to the present to establish the amplitude and phase of the vertical deformation. The Los Angeles basin becomes most inflated one quarter of the way through the year, which is consistent with water table measurements as well as with the end of the rainy season when the aquifer should be at a maximum. The spatial pattern of the amplitude of the annual signal derived from continuous GPS measurements is consistent with the shape of the interferometric fringes. GPS sites both near the NIFZ and in a 20 by 40 km zone within the basin also show significant N-S annual variations that may be related to the differential expansion across the fault. Since these horizontal signals have peak-to-trough amplitudes of 6 mm, they mask the smaller tectonic signals and need to be taken into account when interpreting GPS time series of site position. Moreover, since the groundwater signal appears to have a long-term vertical trend which varies in sign depending on location, it will be difficult to distinguish interseismic tectonic slip along the NIFZ and within the affected areas in the basin.

Sandwell, DT, Smith WHF.  2001.  Bathymetric Estimation. Satellite altimetry and earth sciences : a handbook of techniques and applications. ( Fu L, Cazenave A, Eds.).:441-457., San Diego, Calif. ; London: Academic Abstract
Baer, G, Shamir G, Sandwell D, Bock Y.  2001.  Crustal deformation during 6 years spanning the M (sub w) = 7.2 1995 Nuweiba earthquake, analyzed by Interferometric Synthetic Aperture Radar. Israel Journal of Earth-Sciences. 50( Baer G, Wdowinski S, Eds.).:9-22., Jerusalem, Israel (ISR): Laser Pages Publishing, Jerusalem AbstractWebsite

The November 22, 1995, M (sub w) = 7.2 Nuweiba earthquake occurred along one of the left-stepping segments of the Dead Sea Transform in the Gulf of Elat (Aqaba). We examine the surface deformation patterns in the region by Interferometric Synthetic Aperture Radar (InSAR) for the period 1993 to 1999, which includes the end of one seismic cycle and the beginning of the next. Because the main rupture was under water, ERS coverage is limited to distances of approximately 5 km or more away from the rupture. Pre-earthquake interferograms do not show any detectable deformation along the Gulf. Coseismic interferograms show deformation at distances of up to 50 km from the main rupture, with the highest fringe rate (strain) NW of the rupture termination. Coseismic phase gradient maps show triggered slip along faults parallel to the main rupture (sinistral or normal with the Gulf side down) along the western shore of the Gulf, and in a belt of extensional faults along the eastern shore, striking at angles of about 30 degrees to the major rupture. Postseismic deformation is observed only in a time window of up to 6 months following the mainshock. It was concentrated in the region of the high coseismic strain, and seems to be related to the M (sub L) <4.5 aftershocks in the respective time window.

Sandwell, DT.  2001.  Plate tectonics; a Martian view. Plate tectonics; an insider's history of the modern theory of the Earth. ( Oreskes N, Le Grand H, Eds.)., Boulder, CO, United States (USA): Westview Press, Boulder, CO AbstractWebsite
Sandwell, DT, Sichoix L, Agnew D, Bock Y, Minster JB.  2000.  Near real-time radar interferometry of the Mw 7.1 Hector Mine Earthquake. Geophysical Research Letters. 27:3101-3104.   10.1029/1999gl011209   AbstractWebsite

The Hector Mine Earthquake (Mw 7.1, 16 October 1999) ruptured 45 km of previously mapped and unmapped faults in the Mojave Desert. The ERS-2 satellite imaged the Mojave Desert on 15 September and again on 20 October, just 4 days after the earthquake. Using a newly-developed ground station we acquired both passes and were able to form an interferogram within 20 hours of the second overflight. Estimates of slip along the main rupture are 1-2 meters greater than slip derived from geological mapping. The gradient of the interferometric phase reveals an interesting pattern of triggered slip on adjacent faults as well as a 30 mm deep sink hole along Interstate 40.

Gille, ST, Yale MM, Sandwell DT.  2000.  Global correlation of mesoscale ocean variability with seafloor roughness from satellite altimetry. Geophysical Research Letters. 27:1251-1254.   10.1029/1999gl007003   AbstractWebsite

Both seafloor bathymetry and eddy kinetic energy at the ocean surface can be estimated by making use of satellite altimeters. Comparing the two quantities shows that in regions of the ocean deeper than about 4800 m, surface eddy kinetic energy is greater over smooth abyssal plains than over rough bathymetry, while the opposite is true in shallower waters. Thus in the deep ocean, bottom roughness may dissipate eddy kinetic energy. A simple model indicates that the dissipation rate increases as root-mean-squared bottom roughness increases from 0 to 250 m and decreases to negative values (implying eddy generation) for higher roughness.

Maia, M, Ackermand D, Dehghani GA, Gente P, Hekinian R, Naar D, O'Connor J, Perrot K, Morgan JP, Ramillien G, Revillon S, Sabetian A, Sandwell D, Stoffers P.  2000.  The Pacific-Antarctic Ridge-Foundation hotspot interaction: a case study of a ridge approaching a hotspot. Marine Geology. 167:61-84.   10.1016/s0025-3227(00)00023-2   AbstractWebsite

The Foundation hotspot-Pacific-Antarctic Ridge (PAI) system is the best documented case of a fast spreading ridge approaching a hotspot and interacting with it. The morphology, crustal structure inferred from gravity anomalies and the chemical composition of the lavas of the axial area of the PAR show evidence of the influence of the hotspot, that is presently located roughly 35 km west of the spreading ridge axis. Along-axis variation in the Mantle Bouguer anomaly is about 28 mGal, corresponding to a crustal thickening of 1.5 km where the hotspot is nearer to the PAR. Anomalous ridge elevation is 650 m and the along-axis width of the chemical anomaly is 200 km. A comparison of these axial parameters with those derived for other ridge-hotspot systems, suggests that the amount of plume material reaching the ridge axis is smaller for the Foundation-PAR system. This implies a weaker connection between the plume and the ridge. Cumulative effects of a fast spreading rate and of a fast ridge-hotspot relative motion can be responsible for this weak plume-ridge flow. The how from the hotspot may be less efficiently channelled towards the ridge axis when a fast ridge is rapidly moving towards a hotspot. (C) 2000 Elsevier Science B.V. All rights reserved.

Lyons, SN, Sandwell DT, Smith WHF.  2000.  Three-dimensional estimation of elastic thickness under the Louisville Ridge. Journal of Geophysical Research-Solid Earth. 105:13239-13252.   10.1029/2000jb900065   AbstractWebsite

A three-dimensional approach to estimating elastic thickness is presented which uses dense satellite altimetry and sparse ship bathymetry. This technique is applied to the Louisville Ridge system to study the tectonic history of the region. The inversion is performed as both a first-order approximation and a nonlinear relationship between gravity and topography based on Parker's [1973] equation. While the higher-order effect on the gravity anomaly is nearly zero for most of the region, the magnitude is significant over the summits of the ridge. Nevertheless, the inclusion of the nonlinear terms has only a minor influence on the elastic thickness estimate within each region, lowering the value by similar to 1-2 km compared with the linear result. The incorrect assumption of two dimensionality for circular features exhibits a marked effect on the gravitational anomaly, resulting in false sidelobe structure of nearly 20 mGal for large seamounts. Our elastic thickness estimates are compared with the contradictory values obtained in previous studies by Cazenave and Dominh [1984] and Watts et al. [1988]. We find an increasing elastic thickness along the chain from southeast to northwest, with a discontinuity along the Wishbone scarp. The jump in elastic thickness values northwest of the scarp appears to be an indication of an age discontinuity caused by an extinct spreading center north of the ridge.

Sandwell, DT, Sichoix L.  2000.  Topographic phase recovery from stacked ERS interferometry and a low-resolution digital elevation model. Journal of Geophysical Research-Solid Earth. 105:28211-28222.   10.1029/2000jb900340   AbstractWebsite

A hybrid approach to topographic recovery from ERS interferometry is developed and assessed. Tropospheric/ionospheric artifacts, imprecise orbital information, and layover are key issues in recovering topography and surface deformation from repeat-pass interferometry. Previously, we developed a phase gradient approach to stacking interferograms to reduce these errors and also to reduce the short-wavelength phase noise (see Sandwell and Pi-ice [1998] and Appendix A). Here the method is extended to use a low-resolution digital elevation model to constrain long-wavelength phase errors and an iteration scheme to minimize errors in the computation of phase gradient. We demonstrate the topographic phase recovery on 16-m postings using 25 ERS synthetic aperture radar images from an area of southern California containing 2700 m of relief. On the basis of a comparison with 81 GPS monuments, the ERS-derived topography has a typical absolute accuracy of better than 10 m except in areas of layover. The resulting topographic phase enables accurate two-pass, real-time interferometry even in mountainous areas where traditional phase unwrapping schemes fail. As an example, we form a topography-free (127-m perpendicular baseline) interferogram spanning 7.5 years; fringes from two major earthquakes and aseismic slip on the San Andreas Fault are clearly isolated.

Yale, MM, Sandwell DT.  1999.  Stacked global satellite gravity profiles. Geophysics. 64:1748-1755.   10.1190/1.1444680   AbstractWebsite

Gravity field recovery from satellite altimetry provides global marine coverage but lacks the accuracy and resolution needed for many exploration geophysics studies. The repeating ground tracks of the ERS-1/2, Geosat, and Topex/Poseidon altimeters offer the possibility of improving the accuracy and resolution of gravity anomalies along widely spaced (similar to 40-km spacing) tracks. However, complete ocean coverage is usually needed to convert the sea-surface height (br along-track slope) measurements into gravity anomalies. Here we develop and test a method for constructing stacked gravity profiles by using a published global gravity grid (Sandwell and Smith, 1997), V7.2, as a reference model for the slope-to-gravity anomaly conversion. The method is applied to stacks (averages) of Geosat/ERM (up to 62 cycles), ERS-1/2 (up to 43 cycles), and Topex (up to 142 cycles) satellite altimeter profiles. We assess the accuracies of the ERS-1/2 profiles through a comparison with a gravity model of the northern Gulf of Mexico (profiles provided by EDCON Inc.). The 40 ERS profiles evaluated have a mean rms difference of 3.77 mGal and full wavelength resolution (0.5 coherence) of 24 km. Our processing retains wavelengths as short as 10 km so smaller, large-amplitude features can be resolved, especially in shallow ocean areas (<1000 m deep). We provide an example of combining these higher resolution profiles with lower resolution gravity data in the Caspian Sea.

Baer, G, Sandwell D, Williams S, Bock Y, Shamir G.  1999.  Coseismic deformation associated with the November 1995, M-w=7.1 Nuweiba earthquake, Gulf of Elat (Aqaba), detected by synthetic aperture radar interferometry. Journal of Geophysical Research-Solid Earth. 104:25221-25232.   10.1029/1999jb900216   AbstractWebsite

The November 22, 1995, M-w=7.1 Nuweiba earthquake occurred along one of the left-stepping segments of the Dead Sea Transform in the Gulf of flat (Aqaba). Although it was the largest earthquake along this fault in the last few centuries, little is yet known about the geometry of the rupture, the slip distribution along it, and the nature of postseismic deformation following the main shock. In this study we examine the surface deformation pattern during the coseismic phase of the earthquake in an attempt to better elucidate the earthquake rupture process. As the entire rupture zone was beneath the waters of the Gulf, and there is very little Global Positioning System (GPS) data available in the region for the period spanning the earthquake, interferometric synthetic aperture radar (INSAR) provides the only source of information of surface deformation associated with this earthquake. We chose four synthetic aperture radar (SAR) scenes of about 90x90 km each spanning the rupture area, imaged by the ERS-1 and ERS-2 satellites. The coseismic interferograms show contours of equal satellite-to-ground range changes that correspond to surface displacements due to the earthquake rupture. Interferograms that span the earthquake by 1 week show similar fringe patterns' as those that span the earthquake by 6 months, suggesting that postseismic deformation is minor or confined to the first week after the earthquake. A high displacement gradient is seen on the western side of the Gulf, 20-40 km south of flat and Aqaba, where the total satellite-to-ground range changes are at least 15 cm. The displacement gradient is relatively uniform on the eastern side of the Gulf and the range changes are less than 10 cm. To interpret these results, we compare them to synthetic interferograms generated by elastic dislocation models with a variety of fault parameters. Although selecting the best fit fault parameters is nonunique, we are able to generate a group of simplified model interferograms that provide a reasonable fit to the coseismic interferogram and serve to constrain the location of the fault. The present analysis shows that if the rupture reached the Gulf-bottom surface, the mean sinistral slip along the fault is constrained to about 1.4 m. If surface rupture did not occur, the average sinistral slip is constrained to the range of 1.4-3 m for a fault patch buried 0-4 km below the Gulf-bottom Surface, respectively, with a minor normal component.

Price, EJ, Sandwell DT.  1998.  Small-scale deformations associated with the 1992 Landers, California, earthquake mapped by synthetic aperture radar interferometry phase gradients. Journal of Geophysical Research-Solid Earth. 103:27001-27016.   10.1029/98jb01821   AbstractWebsite

The Landers earthquake (M-w 7.3) occurred on June 28, 1992, and ruptured nearly 100 km of previously mapped and unmapped faults in the Mojave Desert. We use synthetic aperture radar interferometry (InSAR) to examine the cumulative surface deformation between April 24 and August 7, 1992, in a 100 x 100 km region surrounding the northern portion of the earthquake rupture. Also, we introduce a technique for manipulating SAR interferograms to extract short-wavelength displacement information. This technique involves computation and subsequent combination of interferometric phase gradient maps. The InSAR results show significant deformation signatures associated with faults, fractures, dry lake beds, and mountainous regions within 75-100 km of the main rupture. Using the phase gradient method, we are able to extract small-scale deformation patterns near the main rupture. Many of the preexisting, mapped faults within 50 km of the main rupture experienced triggered slip; these include the Old Woman, Lenwood, Johnson Valley, West Calico, and Calico Faults. The InSAR results also indicate right-lateral offsets along secondary fractures trending N-NE within the left-lateral zone of shear between the main rupture and the Johnson Valley Fault. Additionally, there are interesting interferogram fringe signatures surrounding Troy Dry Lake and Coyote Dry Lake that are related to deformation of dry lake beds.

Sandwell, DT, Price EJ.  1998.  Phase gradient approach to stacking interferograms. Journal of Geophysical Research-Solid Earth. 103:30183-30204.   10.1029/1998jb900008   AbstractWebsite

The phase gradient approach is used to construct averages and differences of interferograms without phase unwrapping. Our objectives for change detection are to increase fringe clarity and decrease errors due to tropospheric and ionospheric delay by averaging many interferograms. The standard approach requires phase unwrapping, scaling the phase according to the ratio of the perpendicular baseline, and finally forming the average or difference; however, unique phase unwrapping is usually not possible. Since the phase gradient due to topography is proportional to the perpendicular baseline, phase unwrapping is unnecessary prior to averaging or differencing. Phase unwrapping may be needed to interpret the results, but it is delayed until all of the largest topographic signals are removed. We demonstrate the method by averaging and differencing six interferograms having a suite of perpendicular baselines ranging from 18 to 406 m. Cross-spectral analysis of the difference between two Tandem interferograms provides estimates of spatial resolution, which are used to design prestack filters. A wide range of perpendicular baselines provides the best topographic recovery in terms of accuracy and coverage. Outside of mountainous areas the topography has a relative accuracy of better than 2 m. Residual interferograms (single interferogram minus stack) have tilts across the unwrapped phase that are typically 50 mm in both range and azimuth, reflecting both orbit error and atmospheric delay. Smaller-scale waves with amplitudes of 15 mm are interpreted as atmospheric lee waves. A few Global Positioning System (GPS) control points within a Game could increase the precision to similar to 20 mm for a single interferogram; further improvements may be achieved by stacking residual interferograms.

Yale, M, Sandwell D, Herring A.  1998.  What are the limitations of satellite altimetry? The Leading Edge. 17:73-76.: Society of Exploration Geophysicists   10.1190/1.1437832   AbstractWebsite

Radar altimeter measurements of the marine geoid collected during the Seasat altimeter mission gave geophysicists hope of uncovering the gravity field over all the ocean basins. However, because of insufficient track density, it has taken 16 years for the full potential of the satellite altimeter to be realized. The high‐density coverage obtained by ERS-1 during its geodetic mapping phase (April 1994–March 1995) prompted the U.S. Navy to declassify all of the Geosat altimeter data in June 1995. The combination of these two high‐density data sets provided the first global view of all the ocean basins at a wavelength resolution of 20–30 km.

Sandwell, DT, Johnson CL, Bilotti F, Suppe J.  1997.  Driving forces for limited tectonics on Venus. Icarus. 129:232-244.   10.1006/icar.1997.5721   AbstractWebsite

The very high correlation of geoid height and topography on Venus, along with the high geoid topography ratio, can be interpreted as local isostatic compensation and/or dynamic compensation of topography at depths ranging from 50 to 350 km. For local compensation within the lithosphere, the swell-push force is proportional to the first moment of the anomalous density. Since the long-wavelength isostatic geoid height is also proportional to the first moment of the anomalous density, the swell push force is equal to the geoid height scaled by -g(2)/2 pi G. Because of this direct relationship, the style (i.e., thermal, Airy, or Pratt compensation) and depth of compensation do not need to be specified and can in fact vary over the surface. Phillips (1990) showed that this simple relationship between swell-push force and geoid also holds for dynamic uplift by shear traction on the base of the lithosphere caused by thermal convection of an isoviscous, infinite half-space mantle. Thus for all reasonable isostatic models and particular classes of dynamic models, the geoid height uniquely determines the magnitude of the swell-push body force that is applied to the venusian lithosphere. Given this body force and assuming Venus can be approximated by a uniform thickness thin elastic shell over an inviscid sphere, we calculate the present-day global strain field using equations given in Banerdt (1986); areas of positive geoid height are in a state of extension while areas of negative geoid height are in a state of compression. The present-day model strain field is compared to global strain patterns inferred from Magellan-derived maps of wrinkle ridges and rift zones. Wrinkle ridges, which are believed to reflect distributed compressive deformation, are generally confined to regions with geoid of less than 20 m while rift zones are found primarily along geoid highs. Moreover, much of the observed deformation matches the present-day model strain orientations suggesting that most of the rifts on Venus and many of the wrinkle ridges formed in a stress field similar to the present one. In several large regions, the present-day model strain pattern does not match the observations. This suggests that either the geoid has changed significantly since most of the strain occurred or our model assumptions are incorrect (e.g., there could be local plate boundaries where the stress pattern is discontinuous). Since the venusian lithosphere shows evidence for limited strain, the calculation also provides an estimate of the overall strength of the lithosphere in compression and extension which can be compared with rheological models of yield strength versus depth. At the crests of the major swells, where evidence for rifting is abundant, we find that the temperature gradient must be at least 7 K/km. (C) 1997 Academic Press.

Smith, WHF, Sandwell DT.  1997.  Global sea floor topography from satellite altimetry and ship depth soundings. Science. 277:1956-1962.   10.1126/science.277.5334.1956   AbstractWebsite

A digital bathymetric map of the oceans with a horizontal resolution of 1 to 12 kilometers was derived by combining available depth soundings with high-resolution marine gravity information from the Geosat and ERS-1 spacecraft. Previous global bathymetric maps lacked features such as the 1600-kilometer-long Foundation Seamounts chain in the South Pacific. This map shows relations among the distributions of depth, sea floor area, and sea floor age that do not fit the predictions of deterministic models of subsidence due to lithosphere cooling but may be explained by a stochastic model in which randomly distributed reheating events warm the lithosphere and raise the ocean floor.

Sandwell, DT, Smith WHF.  1997.  Marine gravity anomaly from Geosat and ERS 1 satellite altimetry. Journal of Geophysical Research-Solid Earth. 102:10039-10054.   10.1029/96jb03223   AbstractWebsite

Closely spaced satellite altimeter profiles collected during the Geosat Geodetic Mission (similar to 6 km) and the ERS 1 Geodetic Phase (8 km) are easily converted to grids of vertical gravity gradient and gravity anomaly. The long-wavelength radial orbit error is suppressed below the noise level of the altimeter by taking the along-track derivative of each profile. Ascending and descending slope profiles are then interpolated onto separate uniform grids. These four grids are combined to form comparable grids of east and north vertical deflection using an iteration scheme that interpolates data gaps with minimum curvature. The vertical gravity gradient is calculated directly from the derivatives of the vertical deflection grids, while Fourier analysis is required to construct gravity anomalies from the two vertical deflection grids. These techniques are applied to a combination of high-density data from the dense mapping phases of Geosat and ERS 1 along with lower-density but higher-accuracy profiles from their repeat orbit phases. A comparison with shipboard gravity data shows the accuracy of the satellite-derived gravity anomaly is about 4-7 mGal for random skip tracks. The accuracy improves to 3 mGal when the ship track follows a Geosat Exact Repeat Mission track line. These data provide the first view of the ocean floor structures in many remote areas of the Earth. Some applications include inertial navigation, prediction of seafloor depth, planning shipboard surveys, plate tectonics, isostasy of volcanoes and spreading ridges, and petroleum exploration.

Phiilips, RJ, Johnson CL, Mackwell SJ, Morgan P, Sandwell DT, Zuber MT.  1997.  Lithospheric Mechanics and Dynamics of Venus. Venus II--geology, geophysics, atmosphere, and solar wind environment. ( Bougher SW, Hunten DM, Phillips RJ, Eds.)., Tucson, Ariz.: University of Arizona Press Abstract