We present four companion digital models of the age, age uncertainty, spreading rates, and spreading asymmetries of the world's ocean basins as geographic and Mercator grids with 2 arc min resolution. The grids include data from all the major ocean basins as well as detailed reconstructions of back-arc basins. The age, spreading rate, and asymmetry at each grid node are determined by linear interpolation between adjacent seafloor isochrons in the direction of spreading. Ages for ocean floor between the oldest identified magnetic anomalies and continental crust are interpolated by geological estimates of the ages of passive continental margin segments. The age uncertainties for grid cells coinciding with marine magnetic anomaly identifications, observed or rotated to their conjugate ridge flanks, are based on the difference between gridded age and observed age. The uncertainties are also a function of the distance of a given grid cell to the nearest age observation and the proximity to fracture zones or other age discontinuities. Asymmetries in crustal accretion appear to be frequently related to asthenospheric flow from mantle plumes to spreading ridges, resulting in ridge jumps toward hot spots. We also use the new age grid to compute global residual basement depth grids from the difference between observed oceanic basement depth and predicted depth using three alternative age-depth relationships. The new set of grids helps to investigate prominent negative depth anomalies, which may be alternatively related to subducted slab material descending in the mantle or to asthenospheric flow. A combination of our digital grids and the associated relative and absolute plate motion model with seismic tomography and mantle convection model outputs represents a valuable set of tools to investigate geodynamic problems.

/pdf/age-spreading-rates-and-spreading-asymmetry-of-the-world-s-5bj3otyvfc.pdf

Age, spreading rates, and spreading asymmetry of the world's ocean crust

/pdf/towards-the-standardization-of-sequence-stratigraphy-b61e2b7dm7.pdf

Towards the standardization of sequence stratigraphy

The sudden, widespread glaciation of Antarctica and the associated shift towards colder temperatures at the Eocene/Oligocene boundary (approximately 34 million years ago) (refs 1-4) is one of the most fundamental reorganizations of global climate known in the geologic record. The glaciation of Antarctica has hitherto been thought to result from the tectonic opening of Southern Ocean gateways, which enabled the formation of the Antarctic Circumpolar Current and the subsequent thermal isolation of the Antarctic continent. Here we simulate the glacial inception and early growth of the East Antarctic Ice Sheet using a general circulation model with coupled components for atmosphere, ocean, ice sheet and sediment, and which incorporates palaeogeography, greenhouse gas, changing orbital parameters, and varying ocean heat transport. In our model, declining Cenozoic CO2 first leads to the formation of small, highly dynamic ice caps on high Antarctic plateaux. At a later time, a CO2 threshold is crossed, initiating ice-sheet height/mass-balance feedbacks that cause the ice caps to expand rapidly with large orbital variations, eventually coalescing into a continental-scale East Antarctic Ice Sheet. According to our simulation the opening of Southern Ocean gateways plays a secondary role in this transition, relative to CO2 concentration.

/pdf/rapid-cenozoic-glaciation-of-antarctica-induced-by-declining-4yuos133vq.pdf

Rapid Cenozoic glaciation of Antarctica induced by declining atmospheric CO2.

Clues to the dynamics of the subduction process are found in the many measurable parameters of modern subduction zones. Based on a critical appraisal of the geophysical and geological literature, 26 parameters are estimated for each of 39 modern subduction zones. To isolate causal relationships among these parameters, multivariate analysis is applied to this data set. This analysis yields empirical quantitative relations that predict strain regime and strike-slip faulting in the overriding plate, maximum earthquake magnitude, Benioff zone length, slab dip, arc-trench gap, and maximum trench depth. Excellent correlation is found between length of the Benioff zone and the product of convergence rate and age of the downgoing slab. This relationship is consistent with the conductive heating model of Molnar et al. (1979), if the model is modified in one respect. The rate of heating of the slab is not constant; it is substantially slower during passage of the slab beneath the accretionary prism and overriding plate. The structural style in the overriding plate is determined by its stress state. Though the stress state of overriding plates cannot be quantified, one can classify each individual subduction zone into one of seven semiquantitative strain classes that form a continuum from strongly extensional (class 1, back-arc spreading) to strongly compressional (class 7, active folding and thrusting). This analysis indicates that strain class is probably determined by a linear combination of convergence rate, slab age, and shallow slab dip. Interplate coupling, controlled by convergence rate and slab age, is an important control on strain regime and the primary control on earthquake magnitude. Arc-parallel strike-slip faulting is a common feature of convergent margins, forming a forearc sliver between the strike-slip fault and trench. Optimum conditions for the development of forearc slivers are oblique convergence, a compressional environment, and a continental overriding plate. The primary factor controlling presence of strike-slip faulting is coupling; strongly oblique convergence is not required. The rate of strike-slip faulting is affected by both convergence obliquity and convergence rate. Maximum trench depth is a response to flexure of the underthrusting plate. The amount of flexural deflection at the trench depends on the vertical component of slab pull force, which is very sensitive to slab age and shallow slab dip. Shallow slab dip conforms to the cross-sectional shape of the overriding plate, which is controlled by width of the accretionary prism and duration of subduction. Deep slab dip is affected by the mantle trajectory established at shallow depth but may be modified by mantle flow. Much of the structural complexity of convergent margins is probably attributable to terrane juxtaposition associated with temporal changes in both forearc strike-slip faulting and strain regime. Empirical equations relating subduction parameters can provide both a focus for future theoretical studies and a conceptual and kinematic link between plate tectonics and the geology of subduction zones.

Relations among subduction parameters

[1] We present a new global model for Recent plate velocities, REVEL, describing the relative velocities of 19 plates and continental blocks. The model is derived from publicly available space geodetic (primarily GPS) data for the period 1993–2000. We include an independent and rigorous estimate for GPS velocity uncertainties to assess plate rigidity and propagate these uncertainties to the velocity estimates. The velocity fields for North America, Eurasia, and Antarctica clearly show the effects of glacial isostatic adjustment, and Australia appears to depart from rigid plate behavior in a manner consistent with the mapped intraplate stress field. Two thirds of tested plate pairs agree with the NUVEL-1A geologic (3 Myr average) velocities within uncertainties. Three plate pairs (Caribbean–North America, Caribbean–South America, and North America–Pacific) exhibit significant differences between the geodetic and geologic model that may reflect systematic errors in NUVEL-1A due to the use of seafloor magnetic rate data that do not reflect the full plate rate because of tectonic complexities. Most other differences probably reflect real velocity changes over the last few million years. Several plate pairs (Arabia–Eurasia, Arabia–Nubia, Eurasia–India) move more slowly than the 3 Myr NUVEL-1A average, perhaps reflecting long-term deceleration associated with continental collision. Several other plate pairs, including Nazca–Pacific, Nazca–South America and Nubia–South America, are experiencing slowing that began ∼25 Ma, the beginning of the current phase of Andean crustal shortening.

/pdf/revel-a-model-for-recent-plate-velocities-from-space-geodesy-18oq1lu9ey.pdf

REVEL: A model for Recent plate velocities from space geodesy

Summary. The active Australian-Pacific plate boundary passes through New Zealand. In the north, the Pacific plate subducts beneath the Australian plate with an accretionary wedge forming the eastern continental (Hikurangi) margin of the North Island. The structure of the region behind the Hikurangi margin changes from the extensional back-arc basin under central North Island to a postulated crustal downwarp under the southern North Island. A 100 km long multichannel seismic reflection profile was recorded across the region of crustal downwarp. The data show discontinuous coherent reflectors dipping westwards at the east end of the profile, and east dipping reflectors at the west end, from depths of 9 to 15 s two way time. Simple hand migration of these events indicate that the east dipping reflectors, interpreted as the base of the Australian plate crust, abut against the west dipping reflectors which are interpreted as marking the top of the subducted Pacific plate. Detailed earthquake hypocentre locations in the area show a dipping zone of high seismicity, the top of which coincides closely with the west dipping events, thus supporting this interpretation.

/pdf/seismic-reflection-measurements-behind-the-hikurangi-4omt451izr.pdf

Seismic reflection measurements behind the Hikurangi convergent margin, southern North Island, New Zealand.

Summary results from CRP-1, Cape Roberts Project, Antarctica

https://homepages.dias.ie/sergei/reprints/Fry_et_al_EPSL2014.pdf

Fred Davey

Papers

Seismic reflection measurements behind the Hikurangi convergent margin, southern North Island, New Zealand.

Summary results from CRP-1, Cape Roberts Project, Antarctica

Depth variable crustal anisotropy, patterns of crustal weakness, and destructive earthquakes in Canterbury, New Zealand

Crustal Structure Studies in New Zealand

The Campbell Plateau and its relationship with the Ross Sea, Antarctica