Search NASA⌕ Search

Engineering topics

Sandwell, David T.

Publications and source records attributed to Sandwell, David T..

At least 19 records

A detection of the sea level fingerprint of Greenland Ice Sheet melt

Rapid melting of ice sheets and glaciers drives a unique geometry, or fingerprint, of sea level change. However, the detection of individual fingerprints has been challenging because of sparse observations at high latitudes and the difficulty of disentangling ocean dynamic variability from the signal. Here we predict the fingerprint of Greenland Ice Sheet (GrIS) melt using recent ice mass loss estimates from radar altimetry data and model reconstructions of nearby glaciers and compare this prediction to an independent, altimetry-derived sea surface height trend corrected for ocean dynamic variability in the region adjacent to the ice sheet. A statistically significant correlation between the two fields (P < 0.001) provides an unambiguous observational detection of the near-field sea level fingerprint of recent GrIS melting in our warming world.

54 ENVIRONMENTAL SCIENCES↗

The 1999 Hector Mine Earthquake, Southern California: Vector Near-Field Displacements from ERS InSAR

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 1 km 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-sidedown dip-slip (approximately 1.0 m) which reduces the strike-slip estimates by 1 to 2 m, 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.

Sandwell, David T.↗

Global Bathymetric Prediction For Ocean Modeling and Marine Geophysics

We proposed to construct a complete bathymetric map of the oceans at a 3-10 km resolution by combining all of the available depth soundings collected over the past 30 years with high resolution marine gravity information provided by the Geosat, ERS-1/2, and Topex/Poseidon altimeters. Detailed bathymetry is essential for understanding physical oceanography and marine geophysics. Currents and tides are controlled by the overall shapes of the ocean basins as well as the smaller sharp ocean ridges and seamounts. Because erosion rates are low in the deep oceans, detailed bathymetry reveals the mantle convection patterns, the plate boundaries, the cooling/subsidence of the oceanic lithosphere, the oceanic plateaus, and the distribution of off-ridge volcanoes. We proposed to: (1) Accumulate all available depth soundings collected over the past 30 years; (2) Use the short wavelength (< 160 km) satellite gravity information to interpolate between sparse ship soundings; (3) Improve the resolution of the marine gravity field using enhanced estimates along repeat altimeter profiles together with the dense altimeter measurements; (4) Refine/improve bathymetric predictions using the improved resolution gravity field and also by investigating computer-intensive methods for bathymetric prediction such as inverse theory; and (5) Produce a 'Globe of the Earth' similar to the globe of Venus prepared by the NASA Magellan investigation. This will also include the best available digital land data.

Sandwell, David T.↗

Topographic Phase Recovery from Stacked ERS Interferometry and a Low-Resolution Digital Elevation Model

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 arid Price [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 a seismic slip on the San Andreas Fault are clearly isolated.

Sandwell, David T.↗

Lithospheric bending at subduction zones based on depth soundings and satellite gravity

A global study of trench flexure was performed by simultaneously modeling 117 bathymetric profiles (original depth soundings) and satellite-derived gravity profiles. A thin, elastic plate flexure model was fit to each bathymetry/gravity profile by minimization of the L(sub 1) norm. The six model parameters were regional depth, regional gravity, trench axis location, flexural wavelength, flexural amplitude, and lithospheric density. A regional tilt parameter was not required after correcting for age-related trend using a new high-resolution age map. Estimates of the density parameter confirm that most outer rises are uncompensated. We find that flexural wavelength is not an accurate estimate of plate thickness because of the high curvatures observed at a majority of trenches. As in previous studies, we find that the gravity data favor a longer-wavelength flexure than the bathymetry data. A joint topography-gravity modeling scheme and fit criteria are used to limit acceptable parameter values to models for which topography and gravity yield consistent results. Even after the elastic thicknesses are converted to mechanical thicknesses using the yield strength envelope model, residual scatter obscures the systematic increase of mechanical thickness with age; perhaps this reflects the combination of uncertainties inherent in estimating flexural wavelength, such as extreme inelastic bending and accumulated thermoelastic stress. The bending moment needed to support the trench and outer rise topography increases by a factor of 10 as lithospheric age increases from 20 to 150 Ma; this reflects the increase in saturation bending moment that the lithosphere can maintain. Using a stiff, dry-olivine rheology, we find that the lithosphere of the GDH1 thermal model (Stein and Stein, 1992) is too hot and thin to maintain the observed bending moments. Moreover, the regional depth seaward of the oldest trenches (approximately 150 Ma) exceeds the GDH1 model depths by about 400 m.

Levitt, Daniel A.↗

Gravity over coronae and chasmata on Venus

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. 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 of underthrusting at these coronae. At arcuate segments of Hecate and Parga Chasmata ADCs are about 125 to 150 km, while those at Fauta Corona, four coronae in Eastern Eistla Regio, and an arcuate segment of Wester Parga Chasmata are about 75 km. The GTRs at Fauta Corona, the four coronae in eastern Eistla Regio, and the accurate segments of Hecate, Parga, and Western Parga Chasmata are about 12 to 21 m/km. 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 Coronam, and elsewhere along Parga and Hecate Chasmata.

Schubert, Gerald↗

Arctic geodynamics: Arctic science and ERS-1 satellite altimetry

A detailed gravity field map of the mid Arctic Ocean, spreading ridge system was produced on the basis of ERS-1 satellite altimetry data. Areas of special concern, the Barents and Kara Seas, and areas surrounding the islands of Svalbard, Frans Josef Land and Novoya Zemlya are reviewed. ERS-1 altimetry covers unique Arctic and Antarctic latitudes above 72 degrees. Before ERS-1 it was not possible to study these areas with satellite altimetry. Gravity field solutions for the Barents Sea, portions of the Arctic Ocean and the Norwegian sea are shown. The largest gravity anomalies occur along the Greenland fracture zone as well as along transform faults near Svalbard.

Anderson, Allen Joel↗

Arctic geodynamics: Continental shelf and deep ocean geophysics. ERS-1 satellite altimetry: A first look

An overall review of the Arctic Geodynamics project is presented. A composite gravity field model of the region based upon altimetry data from ERS-1, Geosat, and Seasat is made. ERS-1 altimetry covers unique Arctic and Antarctic latitudes above 72 deg. Both areas contain large continental shelf areas, passive margins, as well as recently formed deep ocean areas. Until ERS-1 it was not possible to study these areas with satellite altimetry. Gravity field solutions for the Barents sea, portions of the Arctic ocean, and the Norwegian sea north of Iceland are shown. The gravity anomalies around Svalbard (Spitsbergen) and Bear island are particularly large, indicating large isostatic anomalies which remain from the recent breakup of Greenland from Scandinavian. Recently released gravity data from the Armed Forces Topographic Service of Russia cover a portion of the Barents and Kara seas. A comparison of this data with the ERS-1 produced gravity field is shown.

Anderson, Allen Joel↗

Flexural ridges, trenches, and outer rises around coronae on Venus

Flexural signatures outboard of Venusian coronal rims are examined with the purpose of inferring the thickness of the planet's elastic lithosphere. Topographic profiles of several prominent coronae which display clear trench and outer rise signatures are presented. Via a thin elastic plate flexure model to characterize the shape of the trench and outer rise, Venusian flexures are found to be similar in both amplitude and wavelength to lithospheric flexures seaward of subduction zones on earth. It is shown that circumferential fractures are concentrated in areas where the topography is curved downward, in good agreement with the high tensile stress predicted by the flexure models. Two scenarios for the development of the ridge-trench-outer rise flexural topography and circumferential fractures of coronae are presented. The first scenario involves reheating and thermal subsidence of the lithosphere interior to the corona, while the second involves expansion of the corona interior and roll back of the subducting lithosphere exterior to the corona.

Sandwell, David T.↗

Joints in Venusian lava flows

The bright polygonal patterns seen in some Magellan images of the volcanic plains are interpreted as joint networks. Several areas displaying such patterns are identified, noting their areal extent, the presence of volcanic features, and the regional topography. Intersection geometries of the lineations defining the polygonal patterns are typical of those found in tensile networks. In addition, the polygonal patterns exhibit no preferred orientation, implying that they are the result of horizontally isotropic stress fields. It is demonstrated that sufficient tensile stresses develop in a cooling Venusian flow to allow joints to form and that the thermal strain predicted is not incompatible with that observed. It is suggested that the tensioned networks are the result of a thermally stressed layer much thicker than a single lava flow. A second model is also considered in which thermal stresses develop as a result of regional elevation of the temperature gradient.

Johnson, Catherine L.↗

Evidence for retrograde lithospheric subduction on Venus

Annular moats and outer rises around large Venus coronas such as Artemis, Latona, and Eithinoha are similar in arcuate planform and topography to the trenches and outer rises of terrestrial subduction zones. On earth, trenches and outer rises are modeled as the flexural response of a thin elastic lithosphere to the bending moment of the subducted slab; this lithospheric flexure model also accounts for the trenches and outer rises outboard of the major coronas on Venus. Accordingly, it is proposed that retrograde lithospheric subduction may be occurring on the margins of the large Venus coronas while compensating back-arc extension is occurring in the expanding coronas interiors. Similar processes may be taking place at other deep arcuate trenches or chasmata on Venus such as those in the Dali-Diana chasmata area of aestern Aphrodite Terra.

Sandwell, David T.↗

Antarctic marine gravity field from high-density satellite altimetry

High-density (about 2-km profile spacing) Geosat/GM altimetry profiles were obtained for Antarctic waters (6-deg S to 72 deg S) and converted to vertical gravity gradient, using Laplace's equation to directly calculate gravity gradient from vertical deflection grids and Fourier analysis to construct gravity anomalies from two vertical deflection grids. The resultant gravity grids have resolution and accuracy comparable to shipboard gravity profiles. The obtained gravity maps display many interesting and previously uncharted features, such as a propagating rift wake and a large 'leaky transform' along the Pacific-Antarctic Rise.

Sandwell, David T.↗

An analysis of ridge axis gravity roughness and spreading rate

Fast and slow spreading ridges have radically different morphologic and gravimetric characteristics. In this study, altimeter measurements from the Geosat Exact Repeat Mission are used to investigate spreading rate dependence of the ridge axis gravity field. Gravity roughness provides an estimate of the amplitude of the gravity anomaly and is robust to small errors in the location of the ridge axis. Gravity roughness as a weighted root mean square of the vertical deflection at 438 ridge crossings on the mid-ocean ridge system is computed. Ridge axis gravity anomalies show a decrease in amplitude with increasing spreading rate up to an intermediate rate of about 60-80 mm/yr and almost no change at higher rates; overall the roughness decreases by a factor of 10 between the lowest and highest rates. In addition to the amplitude decrease, the range of roughness values observed at a given spreading rate shows a similar order of magnitude decrease with transition between 60 and 80 mm/yr. The transition of ridge axis gravity is most apparent at three relatively unexplored locations on the Southeast Indian Ridge and the Pacific-Antarctic Rise; on these intermediate rate ridges the transition occurs abruptly across transform faults.

Small, Christopher↗

Is the Venusian lithosphere subducting?

Using data collected by the Magellan spacecraft, we are exploring the hypothesis that the cooler and more rigid outer layer of Venus (i.e., the lithosphere) is sinking (subducting) into the interior of Venus. If this process is occurring, it provides a mechanism for cooling the interior of Venus and also for recycling the lighter crustal rocks back into the interior. In addition, since subduction zones drive the plate tectonic motion on the Earth, evidence for lithospheric subduction on Venus raises the possibility of limited plate tectonic-like activity on Venus.

Sandwell, David T.↗

Variations in lithospheric thickness on Venus

Recent analyses of Magellan data have indicated many regions exhibiting topograhic flexure. On Venus, flexure is associated predominantly with coronae and the chasmata with Aphrodite Terra. Modeling of these flexural signatures allows the elastic and mechanical thickness of the lithosphere to be estimated. In areas where the lithosphere is flexed beyond its elastic limit the saturation moment provides information on the strength of the lithosphere. Modeling of 12 flexural features on Venus has indicated lithospheric thicknesses comparable with terrestrial values. This has important implications for the venusian heat budget. Flexure of a thin elastic plate due simultaneously to a line load on a continuous plate and a bending moment applied to the end of a broken plate is considered. The mean radius and regional topographic gradient are also included in the model. Features with a large radius of curvature were selected so that a two-dimensional approximation could be used. Comparisons with an axisymmetric model were made for some features to check the validity of the two-dimensional assumption. The best-fit elastic thickness was found for each profile crossing a given flexural feature. In addition, the surface stress and bending moment at the first zero crossing of each profile were also calculated. Flexural amplitudes and elastic thicknesses obtained for 12 features vary significantly. Three examples of the model fitting procedures are discussed.

Johnson, C. L.↗

Flexural models of trench/outer rise topography of coronae on Venus with axisymmetric spherical shell elastic plates

Magellan altimetry has revealed that many coronae on Venus have trenches or moats around their peripheries and rises outboard of the trenches. This trench/outer rise topographic signature is generally associated with the tectonic annulus of the corona. Sandwell and Schubert have interpreted the trench/outer rise topography and the associated tectonic annulus around coronae to be the result of elastic bending of the Venus lithosphere (though the tectonic structures are consequences of inelastic deformation of the lithosphere). They used two-dimensional elastic plate flexure theory to fit topographic profiles across a number of large coronae and inferred elastic lithosphere thicknesses between about 15 and 40 km, similar to inferred values of elastic thickness for the Earth's lithosphere at subduction zones around the Pacific Ocean. Here, we report the results of using axisymmetric elastic flexure theory for the deformation of thin spherical shell plates to interpret the trench/outer rise topography of the large coronae modeled by Sandwell and Schubert and of coronae as small as 250 km in diameter. In the case of a corona only a few hundred kilometers in diameter, the model accounts for the small planform radius of the moat and the nonradial orientation of altimetric traces across the corona. By fitting the flexural topography of coronae we determine the elastic thickness and loading necessary to account for the observed flexure. We calculate the associated bending moment and determine whether the corona interior topographic load can provide the required moment. We also calculate surface stresses and compare the stress distribution with the location of annular tectonic features.

Moore, W.↗

Evidence for retrograde lithospheric subduction on Venus

Though there is no plate tectonics per se on Venus, recent Magellan radar images and topographic profiles of the planet suggest the occurrence of the plate tectonic processes of lithospheric subduction and back-arc spreading. The perimeters of several large coronae (e.g., Latona, Artemis, and Eithinoha) resemble Earth subduction zones in both their planform and topographic profile. The planform of arcuate structures in Eastern Aphrodite were compared with subduction zones of the East Indies. The venusian structures have radii of curvature that are similar to those of terrestrial subduction zones. Moreover, the topography of the venusian ridge/trench structures is highly asymmetric with a ridge on the concave side and a trough on the convex side; Earth subduction zones generally display the same asymmetry.

Sandwell, David T.↗

Report of the panel on lithospheric structure and evolution, section 3

The panel concluded that NASA can contribute to developing a refined understanding of the compositional, structural, and thermal differences between continental and oceanic lithosphere through a vigorous program in solid Earth science with the following objectives: determine the most fundamental geophysical property of the planet; determine the global gravity field to an accuracy of a few milliGals at wavelengths of 100 km or less; determine the global lithospheric magnetic field to a few nanoTeslas at a wavelength of 100 km; determine how the lithosphere has evolved to its present state via acquiring geologic remote sensing data over all the continents.

Chase, Clement G.↗