Walter H. F. Smith

New, improved version of generic mapping tools released

Version 3.1 of the Generic Mapping Tools (GMT) has been released. More than 6000 scientists worldwide are currently using this free, public domain collection of UNIX tools that contains programs serving a variety of research functions. GMT allows users to manipulate (x,y) and (x,y,z) data, and generate PostScript illustrations, including simple x-y diagrams, contour maps, color images, and artificially illuminated, perspective, and/or shaded-relief plots using a variety of map projections (see Wessel and Smith [1991] and Wessel and Smith [1995], for details.). GMT has been installed under UNIX on most types of workstations and both IBM-compatible and Macintosh personal computers.

New version of the generic mapping tools

An updated, new version (3.0) of the Generic Mapping Tools (GMT) has just been released. GMT is a public domain collection of UNIX tools that contains programs to manipulate (x,y) and (x,y,z) data and to generate PostScript illustrations, including simple x-y diagrams, contour maps, color images, and artificially illuminated, perspective, shaded-relief plots using a variety of map projections [Wessel and Smith, 1991]. GMT has been installed on super computers, workstations and personal computers, all running some flavor of UNIX. We estimate that approximately 5000 scientists worldwide are currently using GMT in their work.

Marine gravity anomaly from Geosat and ERS 1 satellite altimetry

Closely spaced satellite altimeter profiles collected during the Geosat Geodetic Mission (-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 ship 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.

Gridding with continuous curvature splines in tension

A gridding method commonly called minimum curvature is widely used in the earth sciences. The method interpolates the data to be gridded with a surface having continuous second derivatives and minimal total squared curvature. The minimum-curvature surface has an analogy in elastic plate flexure and approximates the shape adopted by a thin plate flexed to pass through the data points. Minimum-curvature surfaces may have large oscillations and extraneous inflection points which make them unsuitable for gridding in many of the applications where they are commonly used. These extraneous inflection points can be eliminated by adding tension to the elastic-plate flexure equation. It is straightforward to generalize minimum-curvature gridding algorithms to include a tension parameter; the same system of equations must be solved in either case and only the relative weights of the coefficients change. Therefore, solutions under tension require no more computational effort than minimum-curvature solutions, and any algorithm which can solve the minimum-curvature equations can solve the more general system. We give common geologic examples where minimum-curvature gridding produces erroneous results but gridding with tension yields a good solution. We also outline how to improve the convergence of an iterative method of solution for the gridding equations.

Generic Mapping Tools: Improved Version Released

Generic Mapping Tools (GMT) is an open-source software package for the analysis and display of geoscience data, helping scientists to analyze, interpolate, filter, manipulate, project, and plot time series and gridded data sets. The GMT toolbox includes about 80 core and 40 supplemental program modules sharing a common set of command options, file structures, and documentation. Its power to process data and produce publication-quality graphic presentations has made it vital to a large scientific community that now includes more than 25,000 individual users. GMTs website (http://gmt.soest.hawaii.edu/) exceeds 20,000 visits per month, and server logs show roughly 2000 monthly downloads.

Global Bathymetry and Elevation Data at 30 Arc Seconds Resolution: SRTM30_PLUS

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.

Global marine gravity from retracked Geosat and ERS‐1 altimetry: Ridge segmentation versus spreading rate

[1] 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.

A global, self-consistent, hierarchical, high-resolution shoreline database

We present a high-resolution shoreline data set amalgamated from two databases in the public domain. The data have undergone extensive processing and are free of internal inconsistencies such as erratic points and crossing segments. The shorelines are constructed entirely from hierarchically arranged closed polygons. The data can be used to simplify data searches and data selections or to study the statistical characteristics of shorelines and landmasses. The data set can be accessed both electronically over Internet and from the National Geophysical Data Center, Boulder, Colorado; it comes with access software and routines to facilitate decimation based on a standard line-reduction algorithm.

New global marine gravity model from CryoSat-2 and Jason-1 reveals buried tectonic structure

High-resolution tectonic solutions Detailed topographic maps are available for only a small fraction of the ocean floor, severely limited by the number of ship crossings. Global maps constructed using satellite-derived gravity data, in contrast, are limited in the size of features they can resolve. Sandwell et al. present a new marine gravity model that greatly improves this resolution (see the Perspective by Hwang and Chang). They identify several previously unknown tectonic features, including extinct spreading ridges in the Gulf of Mexico and numerous uncharted seamounts. Science, this issue p. 65; see also p. 32 A high-resolution marine gravity model shows buried tectonic features and ocean-floor topography. [Also see Perspective by Hwang and Chang] 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.

Bathymetric prediction from dense satellite altimetry and sparse shipboard bathymetry

The southern oceans (south of 30°S) are densely covered with satellite-derived gravity data (track spacing 2–4 km) and sparsely covered with shipboard depth soundings (hundreds of kilometers between tracks in some areas). Flexural isostatic compensation theory suggests that bathymetry and downward continued gravity data may show linear correlation in a band of wave-lengths 15–160 km, if sediment cover is thin and seafloor relief is moderate. At shorter wave-lengths, the gravity field is insensitive to seafloor topography because of upward continuation from the seafloor to the sea surface; at longer wavelengths, isostatic compensation cancels out most of the gravity field due to the seafloor topography. We combine this theory with Wiener optimization theory and empirical evidence for gravity noise-to-signal ratios to design low-pass and band-pass filters to use in predicting bathymetry from gravity. The prediction combines long wavelengths (>160 km) from low-pass-filtered soundings with an intermediate-wavelength solution obtained from multiplying downward continued, band-pass-filtered (15–160 km) gravity data by a scaling factor S. S is empirically determined from the correlation between gravity data and existing soundings in the 15–160 km band by robust regression and varies at long wave-lengths. We find that areas with less than 200 m of sediment cover show correlation between gravity and bathymetry significant at the 99% level, and S may be related to the density of seafloor materials in these areas. The prediction has a horizontal resolution limit of 5–10 km in position and is within 100 m of actual soundings at 50% of grid points and within 240 m at 80% of these. In areas of very rugged topography the prediction underestimates the peak amplitudes of seafloor features. Images of the prediction reveal many tectonic features not seen on any existing bathymetrie charts. Because the prediction relies on the gravity field at wavelengths <160 km, it is insensitive to errors in the navigation of sounding lines but also cannot completely reproduce them. Therefore it may be used to locate tectonic features but should not be used to assess hazards to navigation. The prediction is available from the National Geophysical Data Center in both digital and printed form.