Export 69 results:
Sort by: [ Author  (Asc)] Title Type Year
A B C D E F G H I J K L M N O P Q R [S] T U V W X Y Z   [Show ALL]
Sandwell, DT, Zhang B.  1989.  Global Mesoscale Variability from the Geosat Exact Repeat Mission - Correlation with Ocean Depth. Journal of Geophysical Research-Oceans. 94:17971-17984.   10.1029/JC094iC12p17971   AbstractWebsite

We have developed a new technique for extracting global mesoscale variability from satellite altimeter profiles having large radial orbit error (∼3 m). Long-wavelength radial orbit error, as well as other long-wavelength errors (e.g., tides, ionospheric-atmospheric delay, and electromagnetic bias), are suppressed by taking the derivative (slope) of each altimeter profile. A low-pass filter is used to suppress the short-wavelength altimeter noise (λ<100 km). Twenty-two repeat slope profiles are then averaged to produce a mean sea surface slope profile having a precision of about 0.1 μrad. Variations in sea surface slope, which are proportional to changes in current velocity, are obtained by differencing individual profiles from the average profile. Slopes due to mesoscale dynamic topography are typically 1 μrad (i.e., a 0.1-m change in topography over a 100-km distance). Root-mean-square (rms) slope variability as low as 0.2 μrad are found in the southeast Pacific, and maximum slope variations up to 6–8 μrad are found in major western boundary currents (e.g., Gulf Stream, Kuroshio, Falkland, and Agulhas) and Antarctic Circum-polar Current (ACC) systems. The global rms variability map shows previously unknown spatial details that are highly correlated with seafloor topography. Over most areas, the rms slope variability is less than 1 μrad. However at mid-latitudes, areas of higher variability occur in deep water (>3 km) adjacent to continental shelves, spreading ridges, and oceanic plateaus. Variability is low in shallower areas (<3 km). Along the ACC, the meso-scale variability appears to be organized by the many shallow areas in its path. We do not see convincing evidence that variability is higher downstream from topographic protrusions. Instead, the areas of highest variability occur in the deep basins (>4km).

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.  1984.  Thermomechanical Evolution of Oceanic Fracture-Zones. Journal of Geophysical Research. 89:1401-1413.   10.1029/JB089iB13p11401   AbstractWebsite

A fracture zone (FZ) model is constructed from existing models of the thermal and mechanical evolution of the oceanic lithosphere. As the lithosphere cools by conduction, thermal and mechanical boundary layers develop and increase in thickness as (age)1/2. Surface expressions of this development, such as topography and deflection of the vertical (i.e., gravity field), are most apparent along major FZ's because of the sharp age contrast. A simple model, including the effects of lateral heat transport but with no elastic layer, predicts that variations in seafloor depth and deflection of the vertical will become increasingly smooth and ultimately disappear as the FZ ages. Observations, however, show that both the FZ topography and the deflection of the vertical remain sharp as the FZ evolves. These two observations, as well as the observed asymmetry in deflection of the vertical profiles across the Udintsev, Romanche, and Mendocino FZ's, are explained by including a continuous elastic layer in the model. The asymmetry in deflection of the vertical is a consequence of elastic thickness variations across the FZ. Modeling also shows that the evolution of the FZ topography is extremely sensitive to the initial thermal structure near the ridge-transform intersection. Model geoid steps and their development with age are used to access techniques for measuring geoid offsets across FZ's. Reasonable step estimation techniques will underestimate the overall step amplitude by up to 50%. This implies that abnormally thin thermal boundary layers, derived from studies of geoid height versus age, are not required by the data.

Sandwell, D, Schubert G.  2010.  A contraction model for the flattening and equatorial ridge of Iapetus. Icarus. 210:817-822.   10.1016/j.icarus.2010.06.025   AbstractWebsite

Others have explained the excess flattening of Iapetus by a model in which the moon formed at a high spin rate, achieved isostatic equilibrium by very rapid interior heating caused by short-lived radioactive isotopes (SLRI), and subsequently cooled, locking in the excess flattening with respect to an equilibrium shape at its present spin rate. Here we propose an alternate model that does not require an unusually high initial spin rate or the SLRI. The initial formation of Iapetus results in a slightly oblate spheroid with porosity >10%. Radioactive heating by long-lived isotopes warms the interior to about 200 K, at which point it becomes ductile and the interior compacts by 10%, while the 120 km-thick exterior shell remains strong. The shell must deform to match the reduced volume of the ductile interior, and we propose that this deformation occurs along the equator, perhaps focused by a thinner equatorial shell. The final shape of the collapsed sphere matches the observed shape of Iapetus today, described as an oblate ellipse, except along the equator where strain concentration forms a broad ridge. To maintain this non-equilibrium shape, the thickness of the shell must exceed 120 km. Testing the equatorial focusing hypothesis will require a model that includes non-linear processes to account for the finite yield strength of the thick lithosphere. Nevertheless, we show that the stress in the lithosphere generated by the contraction of the interior is about 3 times greater than the stress needed to deform the lithosphere, so some type of lithospheric deformation is expected. (C) 2010 Elsevier Inc. All rights reserved.

Sandwell, DT, Lawver LA, Dalziel IWD, Smith WHF, Wiederspahn M.  1992.  ANTARCTICA Gravity Anomaly and Infrared Satellite Image, USGS MAP 1-2284. : U.S. Geol. Survey Abstract
Sandwell, DT, McAdoo DC.  1988.  Marine Gravity of the Southern-Ocean and Antarctic Margin from Geosat. Journal of Geophysical Research-Solid Earth and Planets. 93:10389-&.   10.1029/JB093iB09p10389   AbstractWebsite

In November of 1986 the U.S. Navy satellite Geosat began collecting unclassified (gravity) altimeter data as part of its exact repeat mission (ERM). For national security reasons the Geosat orbit was arranged so that it closely follows the Seasat satellite altimeter ground track. However, there are two advantages of the Geosat data over the Seasat data. First, because of improvements in altimeter design, Geosat profiles are about 3 times more precise than Seasat profiles. This corresponds to an accuracy of 2–3 μrad (i.e., 2–3 mGal) for wavelengths greater than 20 km. Second, the Geosat altimeter data were collected when the Antarctic ice coverage was minimal (February 1987 to March 1987), while Seasat was only active during an Antarctic winter (June 1978 to September 1978). These new data reveal many previously uncharted seamounts and fracture zones in the extreme southern ocean areas adjacent to Antarctica. Seven large age-offset fracture zones, apparent in the Geosat data, record the early breakup of Gondwana. Finally, the new data reveal the detailed gravity signatures of the passive and active continental margins of Antarctica. These data are an important reconnaissance tool for future studies of these remote ocean areas.

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.

Sandwell, D, Smith-Konter B.  2018.  Maxwell: A semi-analytic 4D code for earthquake cycle modeling of transform fault systems. Computers & Geosciences. 114:84-97.   10.1016/j.cageo.2018.01.009   AbstractWebsite

We have developed a semi-analytic approach (and computational code) for rapidly calculating 3D time-dependent deformation and stress caused by screw dislocations imbedded within an elastic layer overlying a Maxwell viscoelastic half-space. The maxwell model is developed in the Fourier domain to exploit the computational advantages of the convolution theorem, hence substantially reducing the computational burden associated with an arbitrarily complex distribution of force couples necessary for fault modeling. The new aspect of this development is the ability to model lateral variations in shear modulus. Ten benchmark examples are provided for testing and verification of the algorithms and code. One final example simulates interseismic deformation along the San Andreas Fault System where lateral variations in shear modulus are included to simulate lateral variations in lithospheric structure.

Sandwell, DT.  1982.  Thermal isostasy; response of a moving lithosphere to a distributed heat source. Journal of Geophysical Research. 87:1001-1014., Washington, DC, United States (USA): American Geophysical Union, Washington, DC   10.1029/JB087iB02p01001   AbstractWebsite

Spreading ridges and hot spot swells are identified by their high surface heat flow, shallow seafloor, and high geopotential. To understand these and other thermotectonic features, the oceanic lithosphere is modeled as a thermomechanical boundary layer moving through a three-dimensional, time-independent heat source. The heat source mimics the heat advection associated with a spreading ridge or hot spot without introducing the nonlinearities of these flow processes. The Fourier transforms of three Green's functions (response functions), which relate the three observable fields to their common heat source, are determined analytically. Each of these reponse functions is highly anisotropic because the lithosphere is moving with respect to the source. However, the ratio of the gravity response function to the topography response function (i.e., gravity/topography transfer function) is nearly isotropic and has a maximum lying between the flexural wavelength and 2pi times the thickness of the thermal boundary layer. The response functions are most useful for determining the surface heat flow, seafloor topography, and geopotential for complex lithospheric thermal structures. In practice, these three observables are calculated by multiplying the Fourier transform of the heat source by the appropriate response function and inverse transforming the products. Almost any time-independent thermotectonic feature can be modeled using this technique. Included in this report are examples of spreading ridges and thermal swells, although more complex geometries such as ridges offset by transform faults and RRR-type triple junctions can also be modeled. Because forward modeling is both linear and computationally simple, the inverse of this technique could be used to infer some basic characteristics of the heat source directly from the observed fields.

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.

Sandwell, DT.  1986.  Thermal-Stress and the Spacings of Transform Faults. Journal of Geophysical Research-Solid Earth and Planets. 91:6405-6417.   10.1029/JB091iB06p06405   AbstractWebsite

Bathymetric charts are used with satellite altimeter profiles to locate major ridge-transform intersections along five spreading ridges. The ridges are the Mid-Atlantic Ridge, the East Pacific Rise, the Chile Rise, the Pacific-Antarctic Rise, and the Southeast Indian Ridge. Analysis of these data show spacings between transform faults W increase linearly with spreading rate ν (W/ν = 6.28 m.y.). This linear correlation is explained by a thermoelastic model of a cooling strip of lithosphere spreading at a rate ν. The traction-free boundaries of the thin elastic strip simulate cracks in the lithosphere at transform faults. A two-dimensional thermoelastic solution for the in-plane stress shows the largest stress component is tensional and parallel to the ridge. Stresses are zero at the ridge and increase as (age)½ to a maximum value at an age of W/4ν. All stress components are small for ages greater than W/ν. When the transform spacing is large compared with the spreading rate (W/ν > 100 m.y.) thermal stresses exceed the strength of the lithosphere for ages between 0 and 30 Ma. The observed maximum ratio of transform spacing to spreading rate (W/ν = 10 m.y.) results in low thermal stresses that only exceed the strength of the lithosphere for ages less than 1 Ma. Thus transform faults relieve most of the thermal stress. Model predictions also agree with earthquake studies showing that normal faults in young lithosphere have tensional axes aligned with the ridge. Moreover, oceanic intraplate earthquakes rarely occur in lithosphere older than 30 Ma as predicted by the model. These and other geophysical observations confirm Turcotte's hypothesis that transform faults are thermal contraction cracks.

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.

Sandwell, DT, Poehls KA.  1980.  A Compensation Mechanism for the Central Pacific. Journal of Geophysical Research. 85:3751-3758.   10.1029/JB085iB07p03751   AbstractWebsite

Geos 3 derived geoid heights and sea floor topography were averaged into 256 square areas (203 km on a side) for a region in the central Pacific containing a large portion of the Hawaiian Island chain. The whole region is about 500 m shallower than normal sea floor of the same age. The major portion of the depth anomaly is the Hawaiian swell. Data were analyzed using a two-dimensional fast Fourier transform. A transfer function was computed to determine the part of the observed geoid height that is coherent and in phase with the topography. A number of compensation models were tested against this function. Of these models no single physically reasonable model was found to have an acceptable fit. Accordingly, two models were introduced, one compensating short-wavelength topography at a shallow depth (14 km) and the other compensating the longer wavelengths by a deep mechanism. Acceptable deep compensation models include Airy-Heiskanen type compensation at depths between 40 and 80 km. Using the transition wavelength between the two models (1100 km), an estimate is made of the amplitude and shape of the heat anomaly needed to uplift the Hawaiian swell. The peak of the anomaly has an amplitude of 530 mW m−2 and is located 275 km east of Hawaii.

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
Sandwell, DT.  1984.  Along-track deflection of the vertical from Seasat : GEBCO overlays. , Rockville, Md.: U.S. Dept. of Commerce, National Oceanic and Atmospheric Administration, National Ocean Service, National Geodetic Survey, Charting and Geodetic Services : For sale by the National Geodetic Information Center Abstract
Sandwell, DT.  1996.  Exploration of the remote ocean basins with satellite altimeters. McGraw-Hill 1996 yearbook of science & technology. :178-182., Maidenhead: McGraw-Hill Abstract
Sandwell, DT, Wessel P.  2010.  Seamount Discovery Tool Aids Navigation to Uncharted Seafloor Features. Oceanography. 23:34-36. AbstractWebsite
Sandwell, DT, Ruiz MB.  1992.  Along-Track Gravity-Anomalies from Geostat and Seasat Altimetry - Gebco Overlays. Marine Geophysical Researches. 14:165-205.   10.1007/bf01270629   AbstractWebsite

To provide easy access to the large number of Seastat and Geosat altimeter observations collected over the last decade, we have plotted these satellite altimeter profiles as overlays to the General Bathymetric Chart of the Oceans (GEBCO). Each of the 32 overlays displays along-track gravity anomalies for either ascending (southeast to northwest) or descending (northeast to southwest) altimeter passes. Where Seasat and Geosat profiles coincide, only the more accurate Geosat profiles were plotted. In poorly charted southern ocean areas, satellite altimeter profiles reveal many previously undetected features of the seafloor.

Sandwell, DT, Anderson D, Wessel P.  2005.  Plates, plumes, and paradigms. Plates, plumes, and paradigms. ( Foulger GR, Natland JH, Presnall DC, Anderson DL, Eds.).:1-10., Boulder, Colo.: Geological Society of America Abstract
Sandwell, DT, Renkin ML.  1988.  Compensation of Swells and Plateaus in the North Pacific - No Direct Evidence for Mantle Convection. Journal of Geophysical Research-Solid Earth and Planets. 93:2775-2783.   10.1029/JB093iB04p02775   AbstractWebsite

At intermediate and long wavelengths the ratio of geoid height to topography is sensitive to the depth and mode of compensation. A low geoid/topography ratio (<2 m/km) signifies shallow Airy compensation. A higher ratio (∼6 m/km) signifies thermal isostasy and/or dynamic uplift from a mantle plume. A very high geoid/topography ratio (>8 m/km) in conjunction with a poor correlation between geoid height and topography is evidence of mantle convection. After subtracting a reference geoid from the observed geoid, previous studies have found a regular pattern of geoid highs and lows with a characteristic wavelength of 3000–4000 km. Since these geoid highs and lows were poorly correlated with topography and resulted in very high geoid/topography ratios (10–20 m/km), they were believed to reflect the planform of mantle convection. We show that the regular pattern of geoid highs and lows is an artifact caused by truncating the reference geoid at spherical harmonic degree 10. Since the geoid spectrum is “red,” the residual geoid is dominated by degree 11. When the harmonics of the reference geoid are rolled off gradually, the regular pattern of geoid highs and lows disappears. In the Northeast Pacific, the new residual geoid reflects the lithosphere age offsets across the major fracture zones. In the Northwest Pacific, the residual geoid corresponds to isostatically compensated swells and plateaus. We have calculated the geoid/topography ratio for 10 swells and plateaus and have found a range of compensation depths. The highest geoidAopography ratio of 5.5 m/km occurs on the flanks of the Hawaiian Swell. Intermediate ratios occur in four areas, including the Midway Swell. These intermediate ratios reflect a linear combination of the decaying thermal swell and the increasing volume of Airy-compensated seamounts. Low geoid/topography ratios occur over the remaining five areas (e.g., Emperor Seamounts), reflecting the absence of a thermal swell. Our findings do not support the hypothesis that the planform of mantle convection is evident in the geoid. We see only indirect evidence of thermal plumes reheating the lower lithosphere.

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.

Sandwell, DT, Schubert G.  1982.  Geoid Height-Age Relation from Seasat Altimeter Profiles across the Mendocino Fracture-Zone. Journal of Geophysical Research. 87:3949-3958.   10.1029/JB087iB05p03949   AbstractWebsite

Twenty-eight SEASAT altimeter profiles crossing the Mendocino Fracture Zone are used together with seafloor ages determined from magnetic lineations to estimate the change in oceanic geoid height with age, between ages of 15 and 135 m.y. An unbiased estimate of the overall geoid offset along each profile is determined from a least-squares fit of the along-track derivative of the geoid to the geoid slope predicted from a simple two-layer gravitational edge effect model. Uncertainties based upon the statistical properties of each profile are also determined. A geoid slope-age relation is constructed by normalizing the geoid offsets and uncertainties by the age offsets. The results are in agreement with geoid slope-age relations determined from symmetrically spreading ridges (Sandwell and Schubert, 1980). However, the fracture zone estimates have smaller uncertainties and show less scatter. A comparison of these results with the geoid slope-age prediction of the boundary layer cooling model shows that the thermal structure begins to deviate from this model at an early age (20–40 m.y.). A plate cooling model with a thickness of 125 km is most compatible with the geoid slope-age estimates, although significant deviations occur; these may indicate that the lithospheric thermal structure is not entirely age dependent.

Sandwell, DT, Smith WHF.  2009.  Global marine gravity from retracked Geosat and ERS-1 altimetry: Ridge segmentation versus spreading rate. Journal of Geophysical Research-Solid Earth. 114   10.1029/2008jb006008   AbstractWebsite

Three approaches are used to reduce the error in the satellite-derived marine gravity anomalies. First, we have retracked the raw waveforms from the ERS-1 and Geosat/GM missions resulting in improvements in range precision of 40% and 27%, respectively. Second, we have used the recently published EGM2008 global gravity model as a reference field to provide a seamless gravity transition from land to ocean. Third, we have used a biharmonic spline interpolation method to construct residual vertical deflection grids. Comparisons between shipboard gravity and the global gravity grid show errors ranging from 2.0 mGal in the Gulf of Mexico to 4.0 mGal in areas with rugged seafloor topography. The largest errors of up to 20 mGal occur on the crests of narrow large seamounts. The global spreading ridges are well resolved and show variations in ridge axis morphology and segmentation with spreading rate. For rates less than about 60 mm/a the typical ridge segment is 50-80 km long while it increases dramatically at higher rates (100-1000 km). This transition spreading rate of 60 mm/a also marks the transition from axial valley to axial high. We speculate that a single mechanism controls both transitions; candidates include both lithospheric and asthenospheric processes.

Sandwell, DT.  1991.  Geophysical Applications of Satellite Altimetry. Reviews of Geophysics. 29:132-137. AbstractWebsite
Sandwell, D, Fialko Y.  2004.  Warping and cracking of the Pacific plate by thermal contraction. Journal of Geophysical Research-Solid Earth. 109   10.1029/2004jb003091   AbstractWebsite

Lineaments in the gravity field and associated chains of volcanic ridges are widespread on the Pacific plate but are not yet explained by plate tectonics. Recent studies have proposed that they are warps and cracks in the plate caused by uneven thermal contraction of the cooling lithosphere. We show that the large thermoelastic stress produced by top-down cooling is optimally released by lithospheric flexure between regularly spaced parallel cracks. Both the crack spacing and approximate gravity amplitude are predicted by elastic plate theory and variational principle. Cracks along the troughs of the gravity lineaments provide conduits for the generation of volcanic ridges in agreement with new observations from satellite-derived gravity. Our model suggests that gravity lineaments are a natural consequence of lithospheric cooling so that convective rolls or mantle plumes are not required.