Search

SUMMARY Seamounts are volcanic constructs that litter the seafloor. They are important for understanding numerous aspects of marine science, such as plate tectonics, the volcanic melt budget, oceanic circulation, tsunami wave diffraction, tidal energy dissipation and mass wasting. Geometrically, seamounts come in many sizes and shapes, and for the purpose of modelling them for morphological, gravimetric or isostatic studies it is convenient to have simple analytical models whose properties are well known. Here, we present a family of seamount models that may be used in such studies, covering both the initial construction phase and later mass-wasting by sectoral collapses. We also derive realistic axisymmetric density variations that are compatible with observed first-order structure from seismic tomography studies.

The Hawaii-Emperor seamount chain stretches westward from the “Big Island” of Hawaii for over 6000 km until the oldest part of the Emperor chain is subducted at the Kuril and Aleutian trenches. Still regarded as the iconic hotspot-generated seamount chain it has been sampled, mapped, and studied to give insights into numerous oceanic phenomena, such as seamount and volcano formation and associated intraplate magma budgets, the past absolute motions of the Pacific plate and the drift of the Hawaiian plume, and the thermal and mechanical properties of oceanic lithosphere. Much early work on determining the flexural rigidity and equivalent elastic plate thickness that supports the large volcano loads that comprise the chain was focussed on the Hawaiian Ridge, with a major multichannel seismic expedition to the Hawaiian Islands in 1982 providing clear and direct evidence of plate flexure, as well as the indirect effect this deformation has on Earth’s gravity field. Numerous studies have since followed. However, the older part of the chain, beyond the ~50 Ma “bend”, has been much less well studied due to its remoteness, but recent expeditions have provided new marine seismic data to allow an estimation of elastic thickness along the Emperor chain and how they compare to the information we have along the Hawaiian Ridge. Here, we present preliminary work on determining the elastic thickness beneath the Emperor Seamounts. Unlike the Hawaiian Ridge, where the age of the lithosphere at the time of loading (i.e., the difference in age between the underlying seafloor and the formation age of a seamount or oceanic island) is remarkably constant, along the Emperor chain there are major variations in the age of loading, compounded by higher uncertainty due to limited seamount age sampling and the chain’s location within the Cretaceous Quiet Zone. Thus, models with variable elastic thickness as a function of location along the Emperor chain are required. In this presentation, we discuss several models that seek to account for the new seismic imaging of the top and base of flexed oceanic crust (i.e. Moho) at Jimmu guyot while at the same time honouring the characteristic gravimetric signature of the Emperor seamount edifices and their flanking moats. The Optimal Regional Separation (ORS) method is used to isolate the flexural loads, while seismic tomography and different velocity/density relations are explored for assigning suitable load and infill densities that vary spatially, and we search for optimal density and elastic parameters which minimize the misfit to both the residual gravity as well as the seismically observed flexure in the vicinity of Jimmu guyot. The first-order result is a clear thinning of the elastic thickness as we move from south to north: the implications of which we examine here for the tectonic evolution of the northwest Pacific Ocean and the long-term (>106 a) mechanical properties of oceanic lithosphere.

SUMMARYWe examine the rheology and thermal structure of the oceanic lithosphere, expressed in situ by plate flexure beneath the Hawaiian Ridge, where volcanoes of variable sizes have loaded seafloor of approximately the same age, and thus where the lithosphere is expected to have had an approximately uniform age-dependent thermal structure at the time of loading. Shipboard and satellite-derived gravity, as well as multibeam bathymetry data are used in models of plate flexure with curvature-dependent flexural rigidity, the strength of which is limited, in the shallow lithosphere, by brittle failure, and in the deeper lithosphere, by low-temperature plasticity (LTP). We compute relative likelihoods and posterior probabilities for four model parameters: average crustal density ρc, friction coefficient for brittle failure ${\mu _f}$, a pre-exponential weakening factor F controlling the strength of LTP and lithospheric geotherm age t. Results show that if the lithosphere temperatures were as is expected for normal (t = ) 90-Myr-old seafloor at the time of volcano loading, the rheology must be significantly weaker than expected. Specifically, weak brittle strengths (μf ≤ 0.3) show relatively high probabilities for three of the six published LTP flow laws examined. Alternatively, moderate-to-large brittle strengths (μf ≥ 0.5) require all LTP flow laws to be substantially weakened with F = 102 to > 108 or, equivalently, activation energy reduced by 10–35 per cent. In contrast, if the lithosphere has been moderately reheated by the Hawaiian hotspot, represented by geotherms for t = 50–70 Myr, then the flow laws of Evans & Goetze, Raterron et al. and Krancj et al. require little or no weakening. Such modest thermal rejuvenation is allowed by heatflow constraints, supported by regional mantle seismic tomography imaging as well as compositions of mantle xenoliths, and reconciles previously noted discrepancies between the LTP strengths of lithosphere beneath Hawaii versus that entering the Pacific subduction zones.

Absolute plate motion (APM) models derived from hot spot trails must satisfy trail geometries, ages, and paleolatitudes, which requires modeling explicit plume motions. Models lacking plume motions or derived independently from seamounts must also fit these data, provided the implicit plume motions are geodynamically reasonable. We evaluate eight Pacific APM models; three have explicitly modeled plume motions. Seven derive from seamount age progressions; one is a geodynamic model driven by slab pull and ridge push. Using the long‐lived Hawaii‐Emperor and Louisville chains, we derive implicit motions of Hawaii and Louisville plumes for models lacking explicit estimates and compare them with observed paleolatitudes. Inferred plume motions are plausible given rheological constraints on mantle flow, but rates vary considerably and not all models fit the data well. One potential endmember model predicts no APM direction change at 50 Ma, which best explains trails and paleolatitudes, minimizes predicted rotation of Pacific‐Farallon ridge and assumes no true polar motion, yet its implicit plume drift is inconsistent with global circulation models. Alternatively, a global moving hot spot model yields acceptable fits to geometry and ages, implies a major APM change at 50 Ma, but requires significant true polar wander to explain observed paleolatitudes. The inherent inconsistency between age progressions and paleolatitudes may be reconciled by true polar wander, yet questions remain about the accuracy of age progressions for older sections of the Emperor and Louisville chains, the independent geologic evidence for an APM change at 50 Ma, and the uniqueness and relevance of true polar wander estimates.

The Generic Mapping Tools (GMT) software is ubiquitous in the Earth and ocean sciences. As a cross-platform tool producing high-quality maps and figures, it is used by tens of thousands of scientists around the world. The basic syntax of GMT scripts has evolved very slowly since the 1990s, despite the fact that GMT is generally perceived to have a steep learning curve with many pitfalls for beginners and experienced users alike. Reducing these pitfalls means changing the interface, which would break compatibility with thousands of existing scripts. With the latest GMT version 6, we solve this conundrum by introducing a new "modern mode" to complement the interface used in previous versions, which GMT 6 now calls "classic mode." GMT 6 defaults to classic mode and thus is a recommended upgrade for all GMT 5 users. Nonetheless, new users should take advantage of modern mode to make shorter scripts, quickly access commonly used global data sets, and take full advantage of the new tools to draw subplots, place insets, and create animations. Funding Agency National Science Foundation (NSF) Appeared in article as U.S. National Science Foundation MSU Geological Sciences Endowment info:eu-repo/semantics/publishedVersion

We announce a new and integrated system for planning and executing marine geophysical surveys and for scrutinizing and visualizing incoming shipboard data. The system incorporates free software designed for use by scientists and shipboard operators and pertains to underway geophysics and multibeam sonar surveys. Regarding underway data, a crucial first step in the approach is to reduce and merge incoming center beam depth, gravity, and towed magnetic data with navigation, then reformat to the standard exchange format. We are then able to apply established quality control methods including along-track and cross-track analyses to identify error sources and to incrementally build the candidate archive file as new data are acquired. Regarding multibeam data, these are subjected to both an automated error removal scheme for quick visualization and to subsequent ping editing in detail. The candidate archive file and sonar data are automatically and periodically updated and adapted for display in Google Earth, wherein survey planning is also carried out. Data layers are also updated automatically in Google Earth, allowing scientists to focus on visual inspection and interpretation of incoming data. By visualizing underway and sonar data together with reference gravity, magnetic, and bathymetry grids in Google Earth, data familiarity is enhanced and the likelihood of noticing extreme errors increased. We hope scientists will embrace these techniques so that each data set being submitted to a data repository is vetted by the seagoing science party. U.S. National Science FoundationNational Science Foundation (NSF) [1458964, 1558403] Korea Institute of Ocean Science and Technology [PM59941, PM60321]

An updated global bathymetry and topography grid is presented using a spatial sampling interval of 15 arc sec. The bathymetry is produced using a combination of shipboard soundings and depths predicted using satellite altimetry. New data consists of >33.6 million multibeam and singlebeam measurements collated by several institutions, namely, the National Geospatial‐Intelligence Agency, Japan Agency for Marine‐Earth Science and Technology, Geoscience Australia, Center for Coastal and Ocean Mapping, and Scripps Institution of Oceanography. New altimetry data consists of 48, 14, and 12 months of retracked range measurements from Cryosat‐2, SARAL/AltiKa, and Jason‐2, respectively. With respect to SRTM15_PLUS (Olson et al.,), the inclusion of these new data results in a ∼1.4‐km improvement in the minimum wavelength recovered for sea surface free‐air gravity anomalies, a small increase in the accuracy of altimetrically derived predicted depths, and a 1.24% increase, from 9.60% to 10.84%, in the total area of ocean floor that is constrained by shipboard soundings at 15‐arc sec resolution. Bathymetric grid cells constrained by satellite altimetry have estimated uncertainties of ±150 m in the deep oceans and ±180 m between coastlines and the continental rise. Onshore, topography data are sourced from previously published digital elevation models, predominately SRTM‐CGIAR V4.1 between 60°N and 60°S. ArcticDEM is used above 60°N, while Reference Elevation Model of Antarctica is used below 62°S. Auxiliary grids illustrating shipboard data coverage, marine free‐air gravity anomalies, and vertical gradient gradients are also provided in common data formats.

The Generic Mapping Tools (GMT; www.generic-mapping-tools.org) is a well-known set of software for the geosciences, in particular in the marine and solid earth disciplines. GMT is also a prerequisite for many other well-known software infrastructures, including USGS’s ShakeMap for near-real-time maps of ground motion and shaking intensity following significant earthquakes, MBARI/LDEO’s MB-System for multibeam processing and mapping of the seafloor, and Scripps Institution of Oceanography’s GMTSAR for radar interferometric analysis and imaging of crustal deformation. Today, GMT has tens of thousands of users all over the world and remains essential for many terrestrial and planetary data processing and map-making workflows. GMT began its life over 30 years ago when I was a graduate student, and it has enjoyed continuous US National Science Foundation funding since 1993. Leveraging this funding, GMT has succeeded in establishing itself as a collaborative Open Source community resource from the start. Many scientists globally, particularly European scientists, have been instrumental in designing, maintaining, and improving GMT since the early 2000s. As I and several of our core developers approach the end of our academic careers, the GMT team has been pondering how to preserve these collective investments and position GMT to remain an essential geoinformatics infrastructure well into the future. In response, we have made fundamental changes to how GMT works, enabling access to GMT modules from external interfaces such as MATLAB/Octave, Python, and Julia, simplifying user access to large global datasets, and extending our support for the Google Earth platform. However, the biggest impact delivered by the Fall 2019 release of GMT 6 is likely “modern mode”. Modern mode coexists with classic mode (the only previous mode) so that thousands of GMT 4 and 5 scripts will still run as expected. Furthermore, new users will start with modern mode and experience a much-simplified GMT scripting syntax. A new aspect of GMT made possible by modern mode is a greatly simplified animation production. It is clear to all scientists that animations make it easier to elucidate temporal or spatial changes, yet very few scientists create animations as they are traditionally the domain of experts. The GMT team aspires to make animations a task every scientist can do with ease. In this lecture, I will discuss the latest news on GMT, outline modern mode and the external environment access to our modules, highlight a few GMT animations, and present other aspects of our succession planning for strengthening the GMT community.

The plate tectonic revolution and decipherment of magnetic isochrons that followed the pioneering work of Marie Tharp and coworkers in visualizing the seafloor has led to a near-complete understanding of the first-order evolution of global seafloor spreading. However, lagging behind in exploration and understanding are areas of seafloor formed during the Cretaceous Normal Superchron (CNS, 121-83 Ma) when no magnetic reversals were recorded to guide investigators. Thus, for such regions tectonic interpretations are largely driven by mapping and identifying seafloor fabric indicators such as fracture zones, abyssal hills, rift propagators, and extinct spreading centers. Here, we focus on the relict spreading system of the Cretaceous Ellice Basin that was apparently formed by seafloor spreading that split the world’s single largest oceanic plateau Ontong Java Nui, composed of present day Ontong Java, Manihiki and Hikurangi plateaus and other (now subducted) fragments. We examine what was known about this basin from historical single and multibeam bathymetry, what was revealed by the advent of satellite altimetry, and why bathymetric mapping is still required to infer short-length-scale tectonic fabric. High-resolution bathymetric data from the central basin were recently acquired by the University of Hawaii’s vessel R/V Kilo Moana. Evolution of the spreading system is characterized by three main stages of spreading based on directional and morphological analyses of the seafloor fabric indicators identified from bathymetry. Spatially conjugate points symmetric about the spreading central zone were identified at the establishment and cessation of each spreading stage and were assumed to be of the same age to form pseudo-isochrons. Pseudo-isochrons were then utilized in reconstructing the basin through time. The earliest Stage 1 fracture zones trend E-W and consist of multiple closely spaced, parallel fault strands that were indistinguishable in satellite altimetry. A clockwise rotation of the spreading direction led to Stage 2 NW-SE trending fracture zones, which splayed from Stage 1 multistrands. An offset between Stage 2 fracture zones indicates a short-lived late Stage 3 that appears to be the result of a counter-clockwise rotation of the spreading direction shortly before spreading ceased. Seafloor evidence for the initial breakup and rifting between Ontong Java and Manihiki plateaus prior to Stage 1 has yet to be mapped. Basaltic rocks dredged from selected locations along the survey track promise to provide tighter temporal constraints on the evolution of Ellice Basin.

Although partial melt in the asthenosphere is important geodynamically, geophysical constraints on its abundance remain ambiguous. We use a database of seamounts detected using satellite altimetry to constrain the temporal history of erupted asthenospheric melt. We find that intraplate volcanism on young seafloor (

Terrestrial laser scanning (TLS), also known as ground-based Light Detection and Ranging (LiDAR), has been frequently applied to build bare-earth digital elevation models (DEMs) for high-accuracy geomorphology studies. The point clouds acquired from TLS often achieve a spatial resolution at fingerprint (e.g., 3cm3cm) to handprint (e.g., 10cm10cm) level. A downsampling process has to be applied to decimate the massive point clouds and obtain manageable DEMs. It is well known that downsampling can result in aliasing that causes different signal components to become indistinguishable when the signal is reconstructed from the datasets with a lower sampling rate. Conventional DEMs are mainly the results of upsampling of sparse elevation measurements from land surveying, satellite remote sensing, and aerial photography. As a consequence, the effects of aliasing caused by downsampling have not been fully investigated in the open literature of DEMs. This study aims to investigate the spatial aliasing problem of regridding dense TLS data. The TLS data collected from the beach and dune area near Freeport, Texas in the summer of 2015 are used for this study. The core idea of the anti-aliasing procedure is to apply a low-pass spatial filter prior to conducting downsampling. This article describes the successful use of a fourth-order Butterworth low-pass spatial filter employed in the Generic Mapping Tools (GMT) software package as an anti-aliasing filter. The filter can be applied as an isotropic filter with a single cutoff wavelength or as an anisotropic filter with two different cutoff wavelengths in the X and Y directions. The cutoff wavelength for the isotropic filter is recommended to be three times the grid size of the target DEM. An anti-aliasing method for using TLS data has been introduced.A fourth-order Butterworth low-pass spatial filter is applied.The cutoff wavelength is three times the grid size of the target DEM.

The GMT/MATLAB toolbox is a basic interface between MATLAB® (or Octave) and GMT, the Generic Mapping Tools, which allows MATLAB users full access to all GMT modules. Data may be passed between the two programs using intermediate MATLAB structures that organize the metadata needed; these are produced when GMT modules are run. In addition, standard MATLAB matrix data can be used directly as input to GMT modules. The toolbox improves interoperability between two widely used tools in the geosciences and extends the capability of both tools: GMT gains access to the powerful computational capabilities of MATLAB while the latter gains the ability to access specialized gridding algorithms and can produce publication-quality PostScript-based illustrations. The toolbox is available on all platforms and may be downloaded from the GMT website.

Preserving costly marine geophysical trackline data is of paramount importance but variability in priorities, funding, personnel, and technology impact our data archival capacity. We have addressed one crucial facet of this dilemma by devising an open source approach to merge and reduce underway geophysical data and to generate marine geophysical archive files using common command line programs along with the Generic Mapping Tools and its mgd77 supplement. Archive files generated using this approach retain full precision and may be converted automatically to MGD77T, MGD77+, as well as MGD77 formats. We successfully applied the approach to 340 geophysical data sets acquired by R/V Kilo Moana from 2002 to 2018 and in the near term we plan to submit the non-proprietary archive files to the National Centers for Environmental Information's trackline geophysics archive. We encourage international oceanographic communities to explore our methodology as a larger user-base will strengthen the software and the procedures. University of Hawaii at Manoa United States's School of Ocean and Earth Science and Technology U.S. National Science FoundationNational Science Foundation (NSF) [1458964, 1558403] info:eu-repo/semantics/publishedVersion

We present a method for interpolation of sparse two-dimensional vector data. The method is based on the Green's functions of an elastic body subjected to in-plane forces. This approach ensures elastic coupling between the two components of the interpolation. Users may adjust the coupling by varying Poisson's ratio. Smoothing can be achieved by ignoring the smallest eigenvalues in the matrix solution for the strengths of the unknown body forces. We demonstrate the method using irregularly distributed GPS velocities from southern California. Our technique has been implemented in both GMT and MATLAB®.

Relative plate motions provide high-resolution descriptions of motions of plates relative to other plates. Yet geodynamically, motions of plates relative to the mantle are required since such motions can be attributed to forces (e.g., slab pull and ridge push) acting upon the plates. Various reference frames have been proposed, such as the hot spot reference frame, to link plate motions to a mantle framework. Unfortunately, both accuracy and precision of absolute plate motion models lag behind those of relative plate motion models. Consequently, it is paramount to use relative plate motions in improving our understanding of absolute plate motions. A new technique called “ridge-spotting” combines absolute and relative plate motions and examines the viability of proposed absolute plate motion models. We test the method on six published Pacific absolute plate motions models, including fixed and moving hot spot models as well as a geodynamically derived model. Ridge-spotting reconstructs the Pacific-Farallon and Pacific-Antarctica ridge systems over the last 80 Myr. All six absolute plate motion models predict large amounts of northward migration and monotonic clockwise rotation for the Pacific-Farallon ridge. A geodynamic implication of our ridge migration predictions is that the suggestion that the Pacific-Farallon ridge may have been pinned by a large mantle upwelling is not supported. Unexpected or erratic ridge behaviors may be tied to limitations in the models themselves or (for Indo-Atlantic models) discrepancies in the plate circuits used to project models into the Pacific realm. Ridge-spotting is promising and will be extended to include more plates and other ocean basins.

Modern processing of satellite altimetry for use in marine gravimetry involves computing the along-track slopes of observed sea-surface heights, projecting them into east-west and north-south deflection of the vertical grids, and using Laplace's equation to algebraically obtain a grid of the vertical gravity gradient (VGG). The VGG grid is then integrated via overlapping, flat Earth Fourier transforms to yield a free-air anomaly grid. Because of this integration and associated edge effects, the VGG grid retains more short-wavelength information (e.g., fracture zone and seamount signatures) that is of particular importance for plate tectonic investigations. While modeling of gravity anomalies over arbitrary bodies has long been a standard undertaking, similar modeling of VGG anomalies over oceanic features is not commonplace yet. Here we derive analytic solutions for VGG anomalies over simple bodies and arbitrary 2-D and 3-D sources. We demonstrate their usability in determining mass excess and deficiency across the Mendocino fracture zone (a 2-D feature) and find the best bulk density estimate for Jasper seamount (a 3-D feature). The methodologies used herein are implemented in the Generic Mapping Tools, available from gmt.soest.hawaii.edu.

New bathymetric and gravity mapping, refined volume calculations and petrologic analyses show that the Hawaiian volcano Pūhāhonu is the largest and hottest shield volcano on Earth. This ∼12.5-14.1 Ma volcano in the northwest Hawaiian Ridge (NWHR) is twice the size of Mauna Loa volcano (148 ± 29 vs. 74.0 × 10 3 km 3 ), which was assumed to be not only the largest Hawaiian volcano but also the largest known shield volcano. We considered four testable mechanisms to increase magma production, including 1) thinner lithosphere, 2) slower propagation rate, 3) more fertile source, and 4) hotter mantle. The first three of these have been ruled out. The lithosphere was old (∼88 Myrs) when Pūhāhonu was formed, and thus, too thick and cold to allow for greater extents of partial melting. The propagation rate was relatively fast when it erupted (87 km/Myr), so this is another unlikely reason. Source fertility was Kea-like and no more fertile than for other much smaller NWHR volcanoes. A hotter mantle remains the best mechanism to produce the large magma volumes and is consistent with the high forsteritic olivine phenocryst compositions (up to 91.8%) and the calculated high percent of melting (24%). Thus, the gargantuan size of Pūhāhonu reflects its high melting temperature, the highest reported for any Cenozoic basalt. A solitary wave within the Hawaiian plume is the probable cause of Pūhāhonu's higher melting temperature and the resulting increased volume flux given the absence of a more fertile source for Pūhāhonu basalts, as found for many basalts from the Hawaiian Islands.

Oceanic fracture zone traces are widely used in studies of seafloor morphology and plate kinematics. Satellite altimetry missions have resulted in high-resolution gravity maps in which all major fracture zones and other tectonic fabric can be identified, and numerous scientists have digitized such lineaments. We have initiated a community effort to maintain low-cost infrastructure that allows seafloor fabric lineaments to be stored, accessed, and updated. A key improvement over past efforts is our processing software (released as a GMT5 supplement) that allows for semiautomatic corrections to previously digitized fracture zone traces given improved gridded data sets. Here we report on our seafloor fabric processing tools, which complement our database of seafloor fabric lineations, magnetic anomaly identifications, and plate kinematic models.

Tectonic plate and lower-mantle motions are often considered independent. Plate tectonic reconstructions reveal long-lived interactions between mantle plumes and mid-ocean ridges that imply feedback between plate boundaries and the deep mantle.

More than 60% of the Earth's land and shallow marine areas are covered by > 2 km of sediments and sedimentary rocks, with the thickest accumulations on rifted continental margins (Figure 1). Free-air marine gravity anomalies derived from Geosat and ERS-1 satellite altimetry (Fairhead et al., 2001; Sandwell and Smith, 2009; Andersen et al., 2009) outline most of these major basins with remarkable precision. Moreover, gravity and bathymetry data derived from altimetry are used to identify current and paleo-submarine canyons, faults, and local recent uplifts. These geomorphic features provide clues to where to look for large deposits of sediments. While current altimeter data delineate large offshore basins and major structures, they do not resolve some of the smaller geomorphic features and basins (Yale et al., 1998; Fairhead et al., 2001). Improved accuracy and resolution is desirable: to facilitate comparisons between continental margins; as an exploration tool and to permit extrapolation of known structures from well-surveyed areas; to follow fracture zones out of the deep-ocean basin into antecedent continental structures, to define and compare segmentation of margins along strike and identify the position of the continent-ocean boundary; and to study mass anomalies (e.g., sediment type and distribution) and isostatic compensation at continental margins. In this article, we assess the accuracy of a new global marine gravity model based on a wealth of new radar altimetry data and demonstrate that these gravity data are superior in quality to the majority of publicly available academic and government ship gravity data.

SUMMARY We have discovered evidence of a previously unrecognized, large-scale rotation of the Ontong Java Plateau (OJP) recorded in its basement palaeolatitudes. When palaeolatitude differences computed among Ocean Drilling Program Sites 807 and 1183–1187 are plotted versus their present-day site latitude differences, a systematic 2:1 slope bias is evident. While it is possible to resolve this bias by introducing ad hoc tilt corrections at all six sites, drilling records indicate relatively undisturbed conditions at Sites 1183 and 1185–1187. Of the possible causes of the bias, only whole plateau rotation resolves it while honouring the majority of published palaeolatitudes. This implies that only Sites 807 and 1184 palaeolatitudes, both questioned in the literature, are erroneous. A 9 ◦ northward dip previously reported at Site 1184 appears to stem from inclined deposition rather than post-emplacement deformation. We also estimate an 8 ◦ southward tilt correction at Site 807 to make the data set self-consistent. Based on the six sites analysed, we find that OJP may have experienced ∼40 ◦ of clockwise rotation since its formation at ∼123 Ma. In contrast, available Pacific absolute plate motion (APM) models predict less than 10 ◦ of rotation. If our analysis is correct, it suggests that the plateau moved independently of the Pacific Plate early in its history or that Pacific APM models for the Lower Cretaceous are unreliable. While our corrections to Sites 807 and 1184 combined with ∼40 ◦ rotation resolve the internal inconsistencies, the mean palaeolatitude value of Ontong Java remains largely unchanged and is still anomalous with respect to the Pacific apparent polar wander path at ∼123 Ma.

The Taylor (2006) hypothesis suggesting a common origin for the Ontong Java, Manihiki, and Hikurangi large igneous provinces provides an opportunity for a quantitative reconstruction and reassessment of the Ontong Java–Louisville hotspot connection. Our plate tectonic reconstructions of the three plateaus into Ontong Java Nui, or greater Ontong Java, combined with models for Pacific absolute plate motion (APM), allow an analysis of this connection. A new survey of the central Ellice Basin confirms easterly fracture zones, northerly abyssal hill fabric, as well as an area of sigmoidally-southeast-trending fracture zones associated with a late-stage spreading reorientation. From the fracture zone trends we derive new rotation poles for a two-stage model of Ellice Basin opening between the Ontong Java and Manihiki Plateaus. We use these and a single stage pole for separation of the Manihiki and Hikurangi Plateaus, together with three different Pacific APMs, to reconstruct the Ontong Java Nui super plateau back to 123 Ma and compare its predicted location with paleolatitude data obtained from the Ontong Java and Manihiki plateaus. Discrepancies between our Ontong Java Nui reconstructions and Ontong Java and Manihiki paleolatitudes are largest for the fixed Pacific hotspot APM. Assuming a Louisville hotspot source for Ontong Java Nui, remaining disparity between Ontong Java Nui's paleo-location at 123 Ma and published paleomagnetic latitudes for Ontong Java plateau imply that 8°–19° of Louisville hotspot drift or true polar wander may have occurred since the formation of Ontong Java Nui. However, the older portions of the Pacific APMs could easily be biased by a similar amount, making a firm identification of the dominant source of misfit difficult. Prior studies required a combined 26° of hotspot drift, octupole bias effects, and true polar wander just to link the Ontong Java Plateau to Louisville. Consequently, we suggest the super plateau hypothesis and our new reconstructions have considerably strengthened the case for a Louisville plume origin for Ontong Java Nui.

This is the publisher’s final pdf. The published article is copyrighted by Oxford University Press and can be found at: http://petrology.oxfordjournals.org/.

SUMMARY Recent revisions to the satellite-derived vertical gravity gradient (VGG) data reveal more detail of the ocean bottom and have allowed us to develop a non-linear inversion method to detect seamounts in VGG data. We approximate VGG anomalies over seamounts as sums of individual, partially overlapping, elliptical polynomial functions, which allows us to form a non-linear inverse problem by fitting the polynomial model to the observations. Model parameters for a potential seamount include geographical location, peak VGG amplitude, major and minor axes of the elliptical base, and the azimuth of the major axis. The non-linear inversion is very sensitive to the initial values for the location and amplitude; hence, they are constrained by the centre and amplitude of the uppermost contours obtained with a 1-Eotvos contour interval. With these initial conditions from contouring, we execute a step-wise and fully automated inversion and obtain optimal model estimates for potential seamounts; these are statistically evaluated for significance using the Akaike Information Criterion and F tests. A logarithmic barrier technique is applied to ensure positivity of all seamount amplitudes. After automatic and manual inspections of the model parameters we estimate actual heights and basal ellipses of the inspected potential seamounts directly from the predicted bathymetry grid. In this study, we find globally 24 643 potential seamounts (h ≥ 0.1 km) that are located away from continental margins; 8458 potential seamounts are taller than 1 km. Although our global estimate is significantly lower than predictions from previous studies, a first-order reconciliation of the size-frequency statistics obtained from those studies reveals that the previous counts are systematically overestimated. Because of the ambiguity of gravity signals due to small seamounts of h < 1 km and the overlap with abyssal hills, we estimate the global seamount census to lie in the 40 000–55 000 range. The seamount data from this study are accessible from http://www.soest.hawaii.edu/PT/SMTS.

Volcanism observed far from plate boundaries, in the interior of oceanic and continental plates, may result from flow in the underlying mantle. Comparison between a numerical model of mantle flow and the spatial distribution of intraplate volcanism indicates that rapid shear motion in the mantle may drive melting that causes intraplate eruptions. Most of Earth’s volcanism occurs at plate boundaries, in association with subduction or rifting. A few high-volume volcanic fields are observed both at plate boundaries and within plates, fed by plumes upwelling from the deep mantle1. The remaining volcanism is observed away from plate boundaries. It is typically basaltic, effusive and low volume, occurring within continental interiors2,3,4,5,6,7 or creating seamounts on the ocean floor8,9,10,11. This intraplate volcanism has been attributed to various localized processes12 such as cracking of the lithosphere8,13,14, small-scale convection in the mantle beneath the lithosphere15,16,17 or shear-induced melting of low-viscosity pockets of asthenospheric mantle that have become embedded along the base of the lithosphere18. Here we compare the locations of observed intraplate volcanism with global patterns of mantle flow from a numerical model. We find a correlation between recent continental and oceanic intraplate volcanism and areas of the asthenosphere that are experiencing rapid shear due to mantle convection. We detect particularly high correlations in the interior of the continents of western North America, eastern Australia, southern Europe and Antarctica, as well as west of the East Pacific Rise in the Pacific Ocean. We conclude that intraplate volcanism associated with mantle convection is best explained by melting caused by shear flow within the asthenosphere, whereas other localized processes are less important.

We examine how bathymetric mapping coverage varies with distance from the coastline, here a proxy for the effort involved in collecting the data. Distances to the nearest coastline were evaluated on a 1′ × 1′ global grid. We evaluate the density of marine survey track lines, which falls off with increasing distance from the coastline and drops off precipitously for the most remote regions. Bathymetric coverage shows a marked asymmetry between the southern and northern hemispheres, the latter having a factor of 2–4 denser coverage. We find a rapid decrease in data acquisition for previously unexplored regions beginning in 1973–1975. This rate change may reflect a transition from serendipitous exploration to more targeted investigations as the plate tectonics hypothesis became accepted, but it could also reflect the 1970s oil shocks. Coverage of the seafloor varies logarithmically with mapping resolution. At 0.5° resolution, only ∼60% of the seafloor has been mapped; the 50% mark was reached in 1979 and coverage of unexplored seafloor has since been less rapid. For comparison, at 1′ resolution less than 10% of the seafloor has been mapped. Given rising fuel costs we predict the most remote areas will see a decline in future surveys. Better coordination of exploration among agencies and nations could mitigate this concern and improve global coverage, as could future altimetric mapping dedicated to bathymetric prediction.