Export 38 results:
Sort by: Author Title Type [ Year  (Desc)]
Muller, RD, Qin XD, Sandwell DT, Dutkiewicz A, Williams SE, Flament N, Maus S, Seton M.  2016.  The GPlates Portal: Cloud-based interactive 3D visualization of global geophysical and geological data in a web browser. Plos One. 11   10.1371/journal.pone.0150883   AbstractWebsite

The pace of scientific discovery is being transformed by the availability of 'big data' and open access, open source software tools. These innovations open up new avenues for how scientists communicate and share data and ideas with each other and with the general public. Here, we describe our efforts to bring to life our studies of the Earth system, both at present day and through deep geological time. The GPlates Portal ( is a gateway to a series of virtual globes based on the Cesium Javascript library. The portal allows fast interactive visualization of global geophysical and geological data sets, draped over digital terrain models. The globes use WebGL for hardware-accelerated graphics and are cross-platform and cross-browser compatible with complete camera control. The globes include a visualization of a high-resolution global digital elevation model and the vertical gradient of the global gravity field, highlighting small-scale seafloor fabric such as abyssal hills, fracture zones and seamounts in unprecedented detail. The portal also features globes portraying seafloor geology and a global data set of marine magnetic anomaly identifications. The portal is specifically designed to visualize models of the Earth through geological time. These space-time globes include tectonic reconstructions of the Earth's gravity and magnetic fields, and several models of long-wavelength surface dynamic topography through time, including the interactive plotting of vertical motion histories at selected locations. The globes put the on-the-fly visualization of massive data sets at the fingertips of end-users to stimulate teaching and learning and novel avenues of inquiry.

Xu, X, Tong X, Sandwell DT, Milliner CWD, Dolan JF, Hollingsworth J, Leprince S, Ayoub F.  2016.  Refining the shallow slip deficit. Geophysical Journal International. 204:1867-1886.   10.1093/gji/ggv563   Abstract

Geodetic slip inversions for three major (Mw > 7) strike-slip earthquakes (1992 Landers, 1999 Hector Mine and 2010 El Mayor–Cucapah) show a 15–60 per cent reduction in slip near the surface (depth < 2 km) relative to the slip at deeper depths (4–6 km). This significant difference between surface coseismic slip and slip at depth has been termed the shallow slip deficit (SSD). The large magnitude of this deficit has been an enigma since it cannot be explained by shallow creep during the interseismic period or by triggered slip from nearby earthquakes. One potential explanation for the SSD is that the previous geodetic inversions lack data coverage close to surface rupture such that the shallow portions of the slip models are poorly resolved and generally underestimated. In this study, we improve the static coseismic slip inversion for these three earthquakes, especially at shallow depths, by: (1) including data capturing the near-fault deformation from optical imagery and SAR azimuth offsets; (2) refining the interferometric synthetic aperture radar processing with non-boxcar phase filtering, model-dependent range corrections, more complete phase unwrapping by SNAPHU (Statistical Non-linear Approach for Phase Unwrapping) assuming a maximum discontinuity and an on-fault correlation mask; (3) using more detailed, geologically constrained fault geometries and (4) incorporating additional campaign global positioning system (GPS) data. The refined slip models result in much smaller SSDs of 3–19 per cent. We suspect that the remaining minor SSD for these earthquakes likely reflects a combination of our elastic model's inability to fully account for near-surface deformation, which will render our estimates of shallow slip minima, and potentially small amounts of interseismic fault creep or triggered slip, which could ‘make up’ a small percentages of the coseismic SSD during the interseismic period. Our results indicate that it is imperative that slip inversions include accurate measurements of near-fault surface deformation to reliably constrain spatial patterns of slip during major strike-slip earthquakes.

Matthews, KJ, Mullner RD, Sandwell DT.  2016.  Oceanic microplate formation records the onset of India-Eurasia collision. Earth and Planetary Science Letters. 433:204-214.   10.1016/j.epsl.2015.10.040   AbstractWebsite

Mapping of seafloor tectonic fabric in the Indian Ocean, using high-resolution satellite-derived vertical gravity gradient data, reveals an extinct Pacific-style oceanic microplate ('Mammerickx Microplate') west of the Ninetyeast Ridge. It is one of the first Pacific-style microplates to be mapped outside the Pacific basin, suggesting that geophysical conditions during formation probably resembled those that have dominated at eastern Pacific ridges. The microplate formed at the Indian-Antarctic ridge and is bordered by an extinct ridge in the north and pseudofault in the south, whose conjugate is located north of the Kerguelen Plateau. Independent microplate rotation is indicated by asymmetric pseudofaults and rotated abyssal hill fabric, also seen in multibeam data. Magnetic anomaly picks and age estimates calculated from published spreading rates suggest formation during chron 21o (similar to 47.3 Ma). Plate reorganizations can trigger ridge propagation and microplate development, and we propose that Mammerickx Microplate formation is linked with the India-Eurasia collision (initial 'soft' collision). The collision altered the stress regime at the Indian-Antarctic ridge, leading to a change in segmentation and ridge propagation from an establishing transform. Fast Indian-Antarctic spreading that preceded microplate formation, and Kerguelen Plume activity, may have facilitated ridge propagation via the production of thin and weak lithosphere; however both factors had been present for tens of millions of years and are therefore unlikely to have triggered the event. Prior to the collision, the combination of fast spreading and plume activity was responsible for the production of a wide region of undulate seafloor to the north of the extinct ridge and 'W' shaped lineations that record back and forth ridge propagation. Microplate formation provides a precise means of dating the onset of the India-Eurasia collision, and is completely independent of and complementary to timing constraints derived from continental geology or convergence histories. (C) 2015 Elsevier B.V. All rights reserved.

O'Connor, JM, Hoernle K, Muller RD, Morgan JP, Butterworth NP, Hau F, Sandwell DT, Jokat W, Wijbrans JR, Stoffers P.  2015.  Deformation-related volcanism in the Pacific Ocean linked to the Hawaiian-Emperor bend. Nature Geoscience. 8:393-397.   10.1038/ngeo2416   AbstractWebsite

Ocean islands, seamounts and volcanic ridges are thought to form above mantle plumes. Yet, this mechanism cannot explain many volcanic features on the Pacific Ocean floor(1) and some might instead be caused by cracks in the oceanic crust linked to the reorganization of plate motions(1-3). A distinctive bend in the Hawaiian-Emperor volcanic chain has been linked to changes in the direction of motion of the Pacific Plate(4,5), movement of the Hawaiian plume(6-8), or a combination of both(9). However, these links are uncertain because there is no independent record that precisely dates tectonic events that affected the Pacific Plate. Here we analyse the geochemical characteristics of lava samples collected from the Musicians Ridges, lines of volcanic seamounts formed close to the Hawaiian-Emperor bend. We find that the geochemical signature of these lavas is unlike typical ocean island basalts and instead resembles mid-ocean ridge basalts. We infer that the seamounts are unrelated to mantle plume activity and instead formed in an extensional setting, due to deformation of the Pacific Plate. Ar-40/Ar-39 dating reveals that the Musicians Ridges formed during two time windows that bracket the time of formation of the Hawaiian-Emperor bend, 53-52 and 48-47 million years ago. We conclude that the Hawaiian-Emperor bend was formed by plate-mantle reorganization, potentially triggered by a series of subduction events at the Pacific Plate margins.

Malinverni, ES, Sandwell DT, Tassetti AN, Cappelletti L.  2014.  InSAR decorrelation to assess and prevent volcanic risk. European Journal of Remote Sensing. 47:537-556.   10.5721/EuJRS20144730   AbstractWebsite

SAR can be invaluable describing pre-eruption surface deformation and improving the understanding of volcanic processes. This work studies correlation of pairs of SAR images focusing on the influence of surface, climate conditions and acquisition band. Chosen L-band and C-band images (ENVISAT, ERS and ALOS) cover most of the Yellowstone caldera (USA) over a span of 4 years, sampling all the seasons. Interferograms and correlation maps are generated and studied in relation to snow depth and temperature. To isolate temporal decorrelation pairs of images with the shortest baseline are chosen. Results show good performance during winter, bad attitude towards wet snow and good coherence during summer with L-band performing better over vegetation.

Sandwell, DT, Müller DR, Smith WHF, Garcia E, Francis R.  2014.  New global marine gravity model from CryoSat-2 and Jason-1 reveals buried tectonic structure. Science. 346:65-67.   10.1126/science.1258213   AbstractWebsite

Gravity models are powerful tools for mapping tectonic structures, especially in the deep ocean basins where the topography remains unmapped by ships or is buried by thick sediment. We combined new radar altimeter measurements from satellites CryoSat-2 and Jason-1 with existing data to construct a global marine gravity model that is two times more accurate than previous models. We found an extinct spreading ridge in the Gulf of Mexico, a major propagating rift in the South Atlantic Ocean, abyssal hill fabric on slow-spreading ridges, and thousands of previously uncharted seamounts. These discoveries allow us to understand regional tectonic processes and highlight the importance of satellite-derived gravity models as one of the primary tools for the investigation of remote ocean basins.

Marks, KM, Smith WHF, Sandwell DT.  2013.  Significant improvements in marine gravity from ongoing satellite missions. Marine Geophysical Research. 34:137-146.   10.1007/s11001-013-9190-8   AbstractWebsite

Incorporating new altimeter data from CryoSat-2 (30 months), Envisat (18 months), and Jason-1 (7 months) satellites into an updated marine gravity field yields significant reduction in noise and improved resolution. Compared to an older gravity field that did not include the new altimeter data, incoherent power is reduced globally by approximately 2.9 dB at 15 km, 1.6 dB at 20 km, and 1.0 dB at 25 km wavelengths. Coherence analyses between the updated gravity and recent multibeam surveys distributed throughout the world's oceans shows an average increase of similar to 0.023 in mean coherence in the 20-160 km waveband, with the biggest increase (> 0.08) over fast spreading ridges and smallest (< 0.02) over slow spreading ridges and continental shelves. The shortest wavelength at which coherence is above 0.5 decreased globally by similar to 2 km wavelength, with the biggest decrease (> 3.5 km) over fast spreading ridges and smallest (< 1.5 km) over slow spreading ridges and continental shelves. In the Clipperton fracture zone area these improvements result in seamounts that are more accurately located, the detection of smaller seamounts, and the expression of north-south trending abyssal hill fabric. As more altimeter data from the ongoing satellite missions are incorporated into future gravity field updates, finer-scale details of the seafloor will continue to emerge.

Marks, KM, Smith WHF, Sandwell DT.  2010.  Evolution of errors in the altimetric bathymetry model used by Google Earth and GEBCO. Marine Geophysical Research. 31:223-238.   10.1007/s11001-010-9102-0   AbstractWebsite

We analyze errors in the global bathymetry models of Smith and Sandwell that combine satellite altimetry with acoustic soundings and shorelines to estimate depths. Versions of these models have been incorporated into Google Earth and the General Bathymetric Chart of the Oceans (GEBCO). We use Japan Agency for Marine-Earth Science and Technology (JAMSTEC) multibeam surveys not previously incorporated into the models as "ground truth" to compare against model versions 7.2 through 12.1, defining vertical differences as "errors." Overall error statistics improve over time: 50th percentile errors declined from 57 to 55 to 49 m, and 90th percentile errors declined from 257 to 235 to 219 m, in versions 8.2, 11.1 and 12.1. This improvement is partly due to an increasing number of soundings incorporated into successive models, and partly to improvements in the satellite gravity model. Inspection of specific sites reveals that changes in the algorithms used to interpolate across survey gaps with altimetry have affected some errors. Versions 9.1 through 11.1 show a bias in the scaling from gravity in milliGals to topography in meters that affected the 15-160 km wavelength band. Regionally averaged (> 160 km wavelength) depths have accumulated error over successive versions 9 through 11. These problems have been mitigated in version 12.1, which shows no systematic variation of errors with depth. Even so, version 12.1 is in some respects not as good as version 8.2, which employed a different algorithm.

Becker, JJ, Sandwell DT, Smith WHF, Braud J, Binder B, Depner J, Fabre D, Factor J, Ingalls S, Kim SH, Ladner R, Marks K, Nelson S, Pharaoh A, Trimmer R, Von Rosenberg J, Wallace G, Weatherall P.  2009.  Global Bathymetry and Elevation Data at 30 Arc Seconds Resolution: SRTM30_PLUS. Marine Geodesy. 32:355-371.   10.1080/01490410903297766   AbstractWebsite

A new 30-arc second resolution global topography/bathymetry grid (SRTM30_PLUS) has been developed from a wide variety of data sources. Land and ice topography comes from the SRTM30 and ICESat topography, respectively. Ocean bathymetry is based on a new satellite-gravity model where the gravity-to-topography ratio is calibrated using 298 million edited soundings. The main contribution of this study is the compilation and editing of the raw soundings, which come from NOAA, individual scientists, SIO, NGA, JAMSTEC, IFREMER, GEBCO, and NAVOCEANO. The gridded bathymetry is available for ftp download in the same format as the 33 tiles of SRTM30 topography. There are 33 matching tiles of source identification number to convey the provenance of every grid cell. The raw sounding data, converted to a simple common format, are also available for ftp download.

Myer, D, Sandwell D, Brooks B, Foster J, Shimada M.  2008.  Inflation along Kilauea's Southwest Rift Zone in 2006. Journal of Volcanology and Geothermal Research. 177:418-424.   10.1016/j.jvolgeores.2008.06.006   AbstractWebsite

We report on InSAR and GPS results showing the first crustal inflation along the southwest rift zone at Kilauea volcano in over 20 years. Two independent interferograms (May 2-August 2, 2006 and June 22-Nov 7, 2006) from the ALOS PALSAR instrument reveal domal uplift located southwest of the main caldera. The uplift is bounded on the northeast by the caldera and follows the southwest rift zone for about 12 km. It is approximately 8 km wide. We use data derived from permanent GPS stations to calibrate the InSAR displacement data and estimate uplift of 7.7 cm during the first interferogram and 8.9 cm during the second with line-of-sight volumes of 2.8 x 10(6) m(3) and 3.0 X 10(6) m(3) respectively. The earthquake record for the periods before, during, and after inflation shows that a swarm of shallow earthquakes (z<5 km) signaled the beginning of the uplift and that elevated levels of shallow seismicity along the rift zones occurred throughout the uplift period. GPS data indicate that the inflation occurred steadily over nine months between mid-January and mid-October, 2006 making injection of a sill unlikely. We attribute the inflation to recharge of a shallow ductile area under the SWRZ. (c) 2008 Elsevier B.V. All rights reserved.

Sandwell, DT, Myer D, Mellors R, Shimada M, Brooks B, Foster J.  2008.  Accuracy and Resolution of ALOS Interferometry: Vector Deformation Maps of the Father's Day Intrusion at Kilauea. Ieee Transactions on Geoscience and Remote Sensing. 46:3524-3534.   10.1109/tgrs.2008.2000634   AbstractWebsite

We assess the spatial resolution and phase noise of interferograms made from L-band Advanced Land Observing Satellite (ALOS) synthetic-aperture-radar (SAR) data and compare these results with corresponding C-band measurements from European Space Agency Remote Sensing Satellite (ERS). Based on cross-spectral analysis of phase gradients, we find that the spatial resolution of ALOS interferograms is 1.3x better than ERS interferograms. The phase noise of ALOS (i.e., line-of-sight precision in the 100-5000-m wavelength band) is 1.6x worse than ERS (3.3 mm versus 2.1 mm). In both cases, the largest source of error is tropospheric phase delay. Vector deformation maps associated with the June 17, 2007 (Father's day) intrusion along the east rift zone of the Kilauea Volcano were recovered using just four ALOS SAR images from two look directions. Comparisons with deformation vectors from 19 continuous GPS sites show rms line-of-site precision of 14 mm and rms azimuth precision (flight direction) of 71 mm. This azimuth precision is at least 4x better than the corresponding measurements made at C-band. Phase coherence is high even in heavily vegetated areas in agreement with previous results. This improved coherence combined with similar or better accuracy and resolution suggests that L-band ALOS will outperform C-band ERS in the recovery of slow crustal deformation.

Brooks, BA, Foster J, Sandwell D, Wolfe CJ, Okubo P, Poland M, Myer D.  2008.  Magmatically triggered slow slip at Kilauea volcano, Hawaii. Science. 321:1177-1177.   10.1126/science.1159007   AbstractWebsite

We demonstrate that a recent dike intrusion probably triggered a slow fault-slip event (SSE) on Kilauea volcano's mobile south flank. Our analysis combined models of Advanced Land Observing Satellite interferometric dike-intrusion displacement maps with continuous Global Positioning System (GPS) displacement vectors to show that deformation nearly identical to four previous SSEs at Kilauea occurred at far-field sites shortly after the intrusion. We model stress changes because of both secular deformation and the intrusion and find that both would increase the Coulomb failure stress on possible SSE slip surfaces by roughly the same amount. These results, in concert with the observation that none of the previous SSEs at Kilauea was directly preceded by intrusions but rather occurred during times of normal background deformation, suggest that both extrinsic (intrusion-triggering) and intrinsic (secular fault creep) fault processes can lead to SSEs.

Sandwell, D, Rosen P, Moore W, Gurrola E.  2004.  Radar interferometry for measuring tidal strains across cracks on Europa. Journal of Geophysical Research-Planets. 109   10.1029/2004je002276   AbstractWebsite

A major uncertainty in understanding the interaction between the surface of Europa and its ocean below is the present-day activity of fractures. Using well-constrained models for tidal strain and a force balance in a cracked shell, we estimate the shear and normal displacement of cracks that penetrate upward from the base of the shell. If more than half of the plate is fractured, then surface displacements having amplitudes of 3 to 30 cm will be localized in a band 1 to 100 km from the crack. Plate spreading will occur if more than similar to85% of the plate is fractured. The pattern of deformation is sensitive to both the percentage of plate that is cracked and the total thickness of the shell. Repeat-pass radar interferometry could easily detect and map the activity of the cracks during a short experiment from a variety of suitable orbits with repeating ground tracks.

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.

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.

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.

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.

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
Sandwell, DT, Winterer EL, Mammerickx J, Duncan RA, Lynch MA, Levitt DA, Johnson CL.  1995.  Evidence for Diffuse Extension of the Pacific Plate from Pukapuka Ridges and Cross-Grain Gravity Lineations. Journal of Geophysical Research-Solid Earth. 100:15087-15099.   10.1029/95jb00156   AbstractWebsite

Satellite altimeter measurements of marine gravity reveal 100 to 200-km wavelength lineations over a wide area of the Pacific plate oriented roughly in the direction of absolute plate motion. At least three mechanisms have been proposed for their origin: small-scale convective rolls aligned in the direction of absolute plate motion by shear in the asthenosphere; diffuse N-S extension of the lithosphere resulting in lineated zones of extension (boudins); and minihotspots that move slowly with respect to major hotspots and produce intermittent volcanism. Recently, several chains of linear volcanic ridges have been found to be associated with the gravity lineations. Following ridgelike gravity signatures apparent in high-resolution Geosat gravity measurements, we surveyed a series of volcanic ridges that extend northwest from the East Pacific Rise flank for 2600 km onto 40 Ma seafloor. Our survey data, as well as radiometric dates on samples we collected from the ridges, provide tight constraints on their origin: (1) Individual ridge segments and sets of ridges are highly elongate in the direction of present absolute plate motion. (2) The ridges formed along a band 50 to 70-km-wide in the trough of one of the more prominent gravity lineations. (3) Radiometric dates of the largest ridges show no hotspot age progression. Moreover, the directions predicted for minihotspot traces older than 24 Ma do not match observed directions of either the gravity lineations or the ridges. Based on this last observation, we reject the minihotspot model. The occurrence of the ridges in the trough of the gravity lineation is incompatible with the small-scale convection model which would predict increased volcanism above the convective upwelling. We favor the diffuse extension model because it is consistent with the occurrence of ridges in the trough above the more highly extended lithosphere. However, the multibeam data show no evidence for widespread normal faulting of the crust as predicted by the model. Perhaps the fault scarps are buried under more than 30 m of sediments and/or covered by the elongated ridges. Finally, we note that if ridge-push force is much smaller than trench-pull force, then near the ridge axis the direction of maximum tensile stress must be perpendicular to the direction of absolute plate motion.

Schubert, G, Moore WB, Sandwell DT.  1994.  Gravity over Coronae and Chasmata on Venus. Icarus. 112:130-146.   10.1006/icar.1994.1174   AbstractWebsite

The global spherical harmonic model of Venus' gravity field MGNP60FSAAP, with horizontal resolution of about 600 km, shows that most coronae have little or no signature in the gravity field. Nevertheless, some coronae and some segments of chasmata are associated with distinct positive gravity anomalies. No corona has been found to have a negative gravity anomaly. The spatial coincidence of the gravity highs over four closely spaced 300- to 400-km-diameter coronae in Eastern Eistla Regio with the structures themselves is remarkable and argues for a near-surface or lithospheric origin of the gravity signals over such relatively small features. Apparent depths of compensation (ADCs) of the prominent gravity anomalies at Artemis, Latona, and Heng-o Coronae are about 150 to 200 km. The geoid/topography ratios (GTRs) at Artemis, Latona, and Heng-o Coronae lie in the range 32 to 35 m km(-1). The large ADCs and GTRs of Artemis, Latona, and Heng-o Coronae are consistent with topographically related gravity and a thick Venus lithosphere or shallowly compensated topography and deep positive mass anomalies due to subduction or underthrusting at these coronae. At arcuate segments of Hecate and Parga Chasmata ADCs are about 125 to 150 km, while those at Fatua Corona, four coronae in Eastern Eistla Regio, and an arcuate segment of Western Parga Chasma are about 75 km. The GTRs at Fatua Corona, the four coronae in eastern Eistla Regio, and the arcuate segments of Hecate, Parga, and Western Parga Chasmata are about 12 to 21 m km(-1). The ADCs and GTRs of these coronae and arcuate chasmata segments are generally too large to reflect compensation by crustal thickness variations. Instead, they suggest compensation by thermally induced thickness variations in a moderately thick (approximate to 100 km) lithosphere. Alternatively, the gravity signals at these sites could originate from deep positive mass anomalies due to subduction or underthrusting. Weighted linear least squares fits to GTR vs h (long-wavelength topography) data from Heng-o and Fatua Coronae, the four coronae in eastern Eistla Regio, and the arcuate segments of Hecate, Parga, and western Parga Chasmata are consistent with compensation by thermally induced thickness variations of a dense lithosphere above a less dense mantle; the fits imply an average lithosphere thickness of about 180 km and an excess lithospheric density of about 0.5 to 0.7%. Gravity anomalies at the arcuate segments of Dali and Diana Chasmata that form Latona Corona, at Artemis Chasma, and other arcuate segments of Parga and Hecate Chasmata occur on the concave sides of the arcs. By analogy with gravity anomalies of similar horizontal scale (600 km-several thousand kilometers) on the concave sides of terrestrial subduction zone arcs, which are due in large part to subducted lithosphere, it is inferred that the gravity anomalies on Venus are consistent with retrograde subduction at Artemis Chasma, along the northern and southern margins of Latona Corona, and elsewhere along Parga and Hecate Chasmata. (C) 1994 Academic Press, Inc.

Atwater, T, Sclater J, Sandwell D, Severinghaus J, Marlow M.  1993.  Fracture zone traces across the North Pacific Cretaceous Quiet Zone and their tectonic implications. The Mesozoic Pacific : geology, tectonics, and volcanism : a volume in memory of Sy Schlanger. ( Pringle MS, Sager WW, Sliter WV, Stein S, Eds.).:137-154., Washington, DC: American Geophysical Union Abstract
McKenzie, D, Ford PG, Johnson C, Parsons B, Sandwell D, Saunders S, Solomon SC.  1992.  Features on Venus Generated by Plate Boundary Processes. Journal of Geophysical Research-Planets. 97:13533-13544.   10.1029/92JE01350   AbstractWebsite

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.

Marks, KM, Sandwell DT.  1991.  Analysis of Geoid Height Versus Topography for Oceanic Plateaus and Swells Using Nonbiased Linear-Regression. Journal of Geophysical Research-Solid Earth and Planets. 96:8045-8055.   10.1029/91jb00240   AbstractWebsite

We have investigated the relationship between geoid height and topography for 53 oceanic plateaus and swells to determine the mode of compensation. The ratio of geoid height to topography was obtained from the slope of a best line fit by functional analysis (i.e. nonbiased linear regression), a method that minimizes both geoid height and topography residuals. This method is more appropriate than traditional least squares analysis that minimizes only geoid height residuals, because uncertainties are present in both data types. We find that approximately half of the oceanic and continental plateaus analyzed have low ratios that are consistent with Airy-compensated crustal thickening. The remaining plateaus, however, have higher geoid/topography ratios than predicted by the simple Airy model, and the seismically determined Moho depths beneath some of these features are too shallow for crustal thickening alone. A two-layer Airy compensation model, composed of thickened crust underlain by an anomalously low density "mantle root", is used to explain these observations. The Walvis Ridge, and the Agulhas, Crozet, and north Kerguelen plateaus have geoid/topography ratios and Moho depths that are consistent with the two-layer Airy model. The proximity of the Agulhas Plateau to a RRR triple junction during its early development, and the excessive volcanism at active spreading ridges that created the Crozet and north Kerguelen plateaus and the Walvis Ridge, may have produced regions of enhanced depletion and hence the low-density mantle anomalies. If this explanation is correct, then the low-density mantle anomaly persists over time and remains embedded in the lithosphere beneath the oceanic feature.

Marks, KM, Sandwell DT, Vogt PR, Hall SA.  1991.  Mantle Downwelling beneath the Australian-Antarctic Discordance Zone - Evidence from Geoid Height Versus Topography. Earth and Planetary Science Letters. 103:325-338.   10.1016/0012-821x(91)90170-m   AbstractWebsite

The Australian-Antarctic discordance zone (AAD) is an anomalously deep and rough segment of the Southeast Indian Ridge between 120-degrees and 128-degrees-E. A large, negative (deeper than predicted) depth anomaly is centered on the discordance, and a geoid low is evident upon removal of a low-order geoid model and the geoid height-age relation. We investigate two models that may explain these anomalies: a deficiency in ridge-axis magma supply that produces thin oceanic crust (i.e. shallow Airy compensation), and a downwelling and/or cooler mantle beneath the AAD that results in deeper convective-type compensation. To distinguish between these models, we have calculated the ratio of geoid height to topography from the slope of a best line fit by functional analysis (i.e. non-biased linear regression), a method that minimizes both geoid height and topography residuals. Geoid/topography ratios of 2.1 +/- 0.9 m/km for the entire study area (38-degrees-60-degrees-S, 105-degrees-140-degrees-E), 2.3 +/- 1.8 m/km for a subset comprising crust less-than-or-equal-to 25 Ma, and 2.7 +/- 2.0 m/km for a smaller area centered on the AAD were obtained. These ratios are significantly larger than predicted for thin oceanic crust (0.4 m/km), and 2.7 m/km is consistent with downwelling convection beneath young lithosphere. Average compensation depths of 27, 29, and 34 km, respectively, estimated from these ratios suggest a mantle structure that deepens towards the AAD. The deepest compensation (34 km) of the AAD is below the average depth of the base of the young lithosphere (approximately 30 km), and a downwelling of asthenospheric material is implied. The observed geoid height-age slope over the discordance is unusually gradual at -0.133 m/m.y. We calculate that an upper mantle 170-degrees-C cooler and 0.02 g/cm3 denser than normal can explain the shallow slope. Unusually fast shear velocities in the upper 200 km of mantle beneath the discordance, and major-element geochemical trends consistent with small amounts of melting at shallow depths, provide strong evidence for cooler temperatures beneath the AAD.

Mueller, D, Sandwell DT, Tucholke BE, Sclater JG, Shaw PR.  1991.  Depth to basement and geoid expression of the Kane Fracture Zone: A comparison. Marine Geophysical Researches. 13:105-129. AbstractWebsite

Geoid data from Geosat and subsatellite basement depth profiles of the Kane Fracture Zone in the central North Atlantic were used to examine the correlation between the short-wavelength geoid ( lambda = 25-100 km) and the uncompensated basement topography. The processing technique we apply allows the stacking of geoid profiles, although each repeat cycle has an unknown long-wavelength bias. We first formed the derivative of individual profiles, stacked up to 22 repeat cycles, and then integrated the average-slope profile to reconstruct the geoid height. The stacked, filtered geoid profiles have a noise level of about 7 mm in geoid height. The subsatellite basement topography was obtained from a recent compilation of structure contours on basement along the entire length of the Kane Fracture Zone.