Global Sea Floor Topography from Satellite Altimetry and Ship Depth Soundings
26 Sep 1997-Science (American Association for the Advancement of Science)-Vol. 277, Iss: 5334, pp 1956-1962
TL;DR: In this paper, a digital bathymetric map of the oceans with a horizontal resolution of 1 to 12 kilometers was derived by combining available depth soundings with high-resolution marine gravity information from the Geosat and ERS-1 spacecraft.
Abstract: A digital bathymetric map of the oceans with a horizontal resolution of 1 to 12 kilometers was derived by combining available depth soundings with high-resolution marine gravity information from the Geosat and ERS-1 spacecraft. Previous global bathymetric maps lacked features such as the 1600-kilometer-long Foundation Seamounts chain in the South Pacific. This map shows relations among the distributions of depth, sea floor area, and sea floor age that do not fit the predictions of deterministic models of subsidence due to lithosphere cooling but may be explained by a stochastic model in which randomly distributed reheating events warm the lithosphere and raise the ocean floor.
Citations
More filters
TL;DR: In this paper, a relocatable system for generalized inverse (GI) modeling of barotropic ocean tides is described, where the GI penalty functional is minimized using a representer method, which requires repeated solution of the forward and adjoint linearized shallow water equations.
Abstract: A computationally efficient relocatable system for generalized inverse (GI) modeling of barotropic ocean tides is described. The GI penalty functional is minimized using a representer method, which requires repeated solution of the forward and adjoint linearized shallow water equations (SWEs). To make representer computations efficient, the SWEs are solved in the frequency domain by factoring the coefficient matrix for a finite-difference discretization of the second-order wave equation in elevation. Once this matrix is factored representers can be calculated rapidly. By retaining the first-order SWE system (defined in terms of both elevations and currents) in the definition of the discretized GI penalty functional, complete generality in the choice of dynamical error covariances is retained. This allows rational assumptions about errors in the SWE, with soft momentum balance constraints (e.g., to account for inaccurate parameterization of dissipation), but holds mass conservation constraints. Wh...
3,133 citations
TL;DR: MORVEL as discussed by the authors is a new closure-enforced set of angular velocities for the geologically current motions of 25 tectonic plates that collectively occupy 97 per cent of Earth's surface.
Abstract: SUMMARY
We describe best-fitting angular velocities and MORVEL, a new closure-enforced set of angular velocities for the geologically current motions of 25 tectonic plates that collectively occupy 97 per cent of Earth's surface. Seafloor spreading rates and fault azimuths are used to determine the motions of 19 plates bordered by mid-ocean ridges, including all the major plates. Six smaller plates with little or no connection to the mid-ocean ridges are linked to MORVEL with GPS station velocities and azimuthal data. By design, almost no kinematic information is exchanged between the geologically determined and geodetically constrained subsets of the global circuit—MORVEL thus averages motion over geological intervals for all the major plates. Plate geometry changes relative to NUVEL-1A include the incorporation of Nubia, Lwandle and Somalia plates for the former Africa plate, Capricorn, Australia and Macquarie plates for the former Australia plate, and Sur and South America plates for the former South America plate. MORVEL also includes Amur, Philippine Sea, Sundaland and Yangtze plates, making it more useful than NUVEL-1A for studies of deformation in Asia and the western Pacific. Seafloor spreading rates are estimated over the past 0.78 Myr for intermediate and fast spreading centres and since 3.16 Ma for slow and ultraslow spreading centres. Rates are adjusted downward by 0.6–2.6 mm yr−1 to compensate for the several kilometre width of magnetic reversal zones. Nearly all the NUVEL-1A angular velocities differ significantly from the MORVEL angular velocities. The many new data, revised plate geometries, and correction for outward displacement thus significantly modify our knowledge of geologically current plate motions. MORVEL indicates significantly slower 0.78-Myr-average motion across the Nazca–Antarctic and Nazca–Pacific boundaries than does NUVEL-1A, consistent with a progressive slowdown in the eastward component of Nazca plate motion since 3.16 Ma. It also indicates that motions across the Caribbean–North America and Caribbean–South America plate boundaries are twice as fast as given by NUVEL-1A. Summed, least-squares differences between angular velocities estimated from GPS and those for MORVEL, NUVEL-1 and NUVEL-1A are, respectively, 260 per cent larger for NUVEL-1 and 50 per cent larger for NUVEL-1A than for MORVEL, suggesting that MORVEL more accurately describes historically current plate motions. Significant differences between geological and GPS estimates of Nazca plate motion and Arabia–Eurasia and India–Eurasia motion are reduced but not eliminated when using MORVEL instead of NUVEL-1A, possibly indicating that changes have occurred in those plate motions since 3.16 Ma. The MORVEL and GPS estimates of Pacific–North America plate motion in western North America differ by only 2.6 ± 1.7 mm yr−1, ≈25 per cent smaller than for NUVEL-1A. The remaining difference for this plate pair, assuming there are no unrecognized systematic errors and no measurable change in Pacific–North America motion over the past 1–3 Myr, indicates deformation of one or more plates in the global circuit. Tests for closure of six three-plate circuits indicate that two, Pacific–Cocos–Nazca and Sur–Nubia–Antarctic, fail closure, with respective linear velocities of non-closure of 14 ± 5 and 3 ± 1 mm yr−1 (95 per cent confidence limits) at their triple junctions. We conclude that the rigid plate approximation continues to be tremendously useful, but—absent any unrecognized systematic errors—the plates deform measurably, possibly by thermal contraction and wide plate boundaries with deformation rates near or beneath the level of noise in plate kinematic data.
2,089 citations
TL;DR: Using the Levitus climatology, the authors showed that 2.1 TW (terawatts) is required to maintain the global abyssal density distribution against 30 Sverdrups of deep water formation.
1,958 citations
TL;DR: EGM2008 as mentioned in this paper is a spherical harmonic model of the Earth's gravitational potential, developed by a least squares combination of the ITG-GRACE03S gravitational model and its associated error covariance matrix, with the gravitational information obtained from a global set of area-mean free-air gravity anomalies defined on a 5 arc-minute equiangular grid.
Abstract: [1] EGM2008 is a spherical harmonic model of the Earth's gravitational potential, developed by a least squares combination of the ITG-GRACE03S gravitational model and its associated error covariance matrix, with the gravitational information obtained from a global set of area-mean free-air gravity anomalies defined on a 5 arc-minute equiangular grid This grid was formed by merging terrestrial, altimetry-derived, and airborne gravity data Over areas where only lower resolution gravity data were available, their spectral content was supplemented with gravitational information implied by the topography EGM2008 is complete to degree and order 2159, and contains additional coefficients up to degree 2190 and order 2159 Over areas covered with high quality gravity data, the discrepancies between EGM2008 geoid undulations and independent GPS/Leveling values are on the order of ±5 to ±10 cm EGM2008 vertical deflections over USA and Australia are within ±11 to ±13 arc-seconds of independent astrogeodetic values These results indicate that EGM2008 performs comparably with contemporary detailed regional geoid models EGM2008 performs equally well with other GRACE-based gravitational models in orbit computations Over EGM96, EGM2008 represents improvement by a factor of six in resolution, and by factors of three to six in accuracy, depending on gravitational quantity and geographic area EGM2008 represents a milestone and a new paradigm in global gravity field modeling, by demonstrating for the first time ever, that given accurate and detailed gravimetric data, asingle global model may satisfy the requirements of a very wide range of applications
1,755 citations
British Antarctic Survey1, University of Bristol2, Columbia University3, National Institute of Geophysics and Volcanology4, University of Aberdeen5, University of Texas at Austin6, Centro de Estudios Científicos7, Université libre de Bruxelles8, University of Washington9, Swansea University10, Institute for Geosciences and Natural Resources11, Technical University of Denmark12, National Institute of Polar Research13, California Institute of Technology14, University of Kansas15, Stockholm University16, St. Olaf College17, Norwegian Polar Institute18, Wallops Flight Facility19, University of Canterbury20, University of Oslo21, University of California, Santa Barbara22, University of California, Irvine23, University of York24, Australian Antarctic Division25, Newcastle University26, Goddard Space Flight Center27, Polar Research Institute of China28
TL;DR: Bedmap2 as discussed by the authors is a suite of gridded products describing surface elevation, ice-thickness and the seafloor and subglacial bed elevation of the Antarctic south of 60° S. In particular, the Bedmap2 ice thickness grid is made from 25 million measurements, over two orders of magnitude more than were used in Bedmap1.
Abstract: We present Bedmap2, a new suite of gridded products describing surface elevation, ice-thickness and the seafloor and subglacial bed elevation of the Antarctic south of 60° S. We derived these products using data from a variety of sources, including many substantial surveys completed since the original Bedmap compilation (Bedmap1) in 2001. In particular, the Bedmap2 ice thickness grid is made from 25 million measurements, over two orders of magnitude more than were used in Bedmap1. In most parts of Antarctica the subglacial landscape is visible in much greater detail than was previously available and the improved data-coverage has in many areas revealed the full scale of mountain ranges, valleys, basins and troughs, only fragments of which were previously indicated in local surveys. The derived statistics for Bedmap2 show that the volume of ice contained in the Antarctic ice sheet (27 million km3) and its potential contribution to sea-level rise (58 m) are similar to those of Bedmap1, but the mean thickness of the ice sheet is 4.6% greater, the mean depth of the bed beneath the grounded ice sheet is 72 m lower and the area of ice sheet grounded on bed below sea level is increased by 10%. The Bedmap2 compilation highlights several areas beneath the ice sheet where the bed elevation is substantially lower than the deepest bed indicated by Bedmap1. These products, along with grids of data coverage and uncertainty, provide new opportunities for detailed modelling of the past and future evolution of the Antarctic ice sheets.
1,678 citations
References
More filters
Book•
01 Jan 1999
TL;DR: Cotton and Wilkinson's Advanced Inorganic Chemistry (AIC) as discussed by the authors is one of the most widely used inorganic chemistry books and has been used for more than a quarter century.
Abstract: For more than a quarter century, Cotton and Wilkinson's Advanced Inorganic Chemistry has been the source that students and professional chemists have turned to for the background needed to understand current research literature in inorganic chemistry and aspects of organometallic chemistry. Like its predecessors, this updated Sixth Edition is organized around the periodic table of elements and provides a systematic treatment of the chemistry of all chemical elements and their compounds. It incorporates important recent developments with an emphasis on advances in the interpretation of structure, bonding, and reactivity.From the reviews of the Fifth Edition:* "The first place to go when seeking general information about the chemistry of a particular element, especially when up-to-date, authoritative information is desired." -Journal of the American Chemical Society.* "Every student with a serious interest in inorganic chemistry should have [this book]." -Journal of Chemical Education.* "A mine of information . . . an invaluable guide." -Nature.* "The standard by which all other inorganic chemistry books are judged."-Nouveau Journal de Chimie.* "A masterly overview of the chemistry of the elements."-The Times of London Higher Education Supplement.* "A bonanza of information on important results and developments which could otherwise easily be overlooked in the general deluge of publications." -Angewandte Chemie.
12,231 citations
TL;DR: The Generic Mapping Tools (GMT) is introduced, which is a free, public domain software package that can be used to manipulate columns of tabular data, time series, and gridded data sets and to display these data in a variety of forms ranging from simple x-y plots to maps and color, perspective, and shaded-relief illustrations.
Abstract: When creating camera-ready figures, most scientists are familiar with the sequence of raw data → processing → final illustration and with the spending of large sums of money to finalize papers for submission to scientific journals, prepare proposals, and create overheads and slides for various presentations. This process can be tedious and is often done manually, since available commercial or in-house software usually can do only part of the job.
To expedite this process, we introduce the Generic Mapping Tools (GMT), which is a free, public domain software package that can be used to manipulate columns of tabular data, time series, and gridded data sets and to display these data in a variety of forms ranging from simple x-y plots to maps and color, perspective, and shaded-relief illustrations. GMT uses the PostScript page description language, which can create arbitrarily complex images in gray tones or 24-bit true color by superimposing multiple plot files. Line drawings, bitmapped images, and text can be easily combined in one illustration. PostScript plot files are device-independent, meaning the same file can be printed at 300 dots per inch (dpi) on an ordinary laserwriter or at 2470 dpi on a phototypesetter when ultimate quality is needed. GMT software is written as a set of UNIX tools and is totally self contained and fully documented. The system is offered free of charge to federal agencies and nonprofit educational organizations worldwide and is distributed over the computer network Internet.
4,128 citations
TL;DR: In this paper, a simple cooling model and the plate model were proposed to account for the variation in depth and heat flow with increasing age of the ocean floor. But the results were limited to the North Pacific and North Atlantic basins.
Abstract: Two models, a simple cooling model and the plate model, have been advanced to account for the variation in depth and heat flow with increasing age of the ocean floor. The simple cooling model predicts a linear relation between depth and t½, and heat flow and 1/t½, where t is the age of the ocean floor. We show that the same t½ dependence is implicit in the solutions for the plate model for sufficiently young ocean floor. For larger ages these relations break down, and depth and heat flow decay exponentially to constant values. The two forms of the solution are developed to provide a simple method of inverting the data to give the model parameters. The empirical depth versus age relation for the North Pacific and North Atlantic has been extended out to 160 m.y. B.P. The depth initially increases as t½, but between 60 and 80 m.y. B.P. the variation of depth with age departs from this simple relation. For older ocean floor the depth decays exponentially with age toward a constant asymptotic value. Such characteristics would be produced by a thermal structure close to that of the plate model. Inverting the data gives a plate thickness of 125±10 km, a bottom boundary temperature of 1350°±275°C, and a thermal expansion coefficient of (3.2±1.1) × 10−5°C−1. Between 0 and 70 m.y. B.P. the depth can be represented by the relation d(t) = 2500 + 350t½ m, with t in m.y. B.P., and for regions older than 20 m.y. B.P. by the relation d(t) = 6400 - 3200 exp (−t/62.8) m. The heat flow data were treated in a similar, but less extensive manner. Although the data are compatible with the same model that accounts for the topography, their scatter prevents their use in the same quantitative fashion. Our analysis shows that the heat flow only responds to the bottom boundary at approximately twice the age at which the depth does. Within the scatter of the data, from 0 to 120 m.y. B.P., the heat flow pan be represented by the relation q(t) = 11.3/t½ μcal cm−2s−1. The previously accepted view that the heat flow observations approach a constant asymptotic value in the old ocean basins needs to be tested more stringently. The above results imply that a mechanism is required to supply heat at the base of the plate.
2,667 citations
TL;DR: In this article, 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 is used to construct gravity anomalies from the two vertical deflection grids.
Abstract: 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.
1,695 citations
TL;DR: In this paper, the authors proposed a joint fitting of heat flow and bathymetry with a model with a hotter, thinner oceanic lithosphere than in previous models, including those from older lithosphere previously treated as anomalous.
Abstract: Variations in sea-floor depth and heat flow with age provide the main constraints on the thermal structure and evolution of the oceanic lithosphere. Joint fitting of heat flow and bathymetry yields a model with a hotter, thinner lithosphere than in previous models. The new model provides a significantly better fit to the data, including those from older lithosphere previously treated as anomalous. This will facilitate the analysis of lithospheric processes, including the effects of mid-plate volcanism and swells, regional subsidence, and hydrothermal circulation near spreading centers.
1,379 citations