Search

Unravelling causes for absolute plate velocity change and continental dynamic topography change is challenging because of the interdependence of large-scale geodynamic driving processes. Here, we unravel a clear spatio-temporal relation between latest Cretaceous-Early Cenozoic subduction at the northern edge of the Australian plate, Early Cenozoic Australian plate motion changes and Cenozoic topography evolution of the Australian continent. We present evidence for a ~4000 km wide subduction zone, which culminated in ophiolite obduction and arc-continent collision in the New Guinea-Pocklington Trough region during subduction termination, coinciding with cessation of spreading in the Coral Sea, a ~5 cm/yr decrease in northward Australian plate velocity, and slab detachment. Renewed northward motion caused the Australian plate to override the sinking subduction remnant, which we detect with seismic tomography at 800-1200 km depth in the mantle under central-southeast Australia at a position predicted by our absolute plate reconstructions. With a numerical model of slab sinking and mantle flow we predict a long-wavelength subsidence (negative dynamic topography) migrating southward from ~50 Ma to present, explaining Eocene-Oligocene subsidence of the Queensland Plateau, ~330 m of late Eocene-early Oligocene subsidence in the Gulf of Carpentaria, Oligocene-Miocene subsidence of the Marion Plateau, and providing a first-order fit to the present-day, ~200 m deep, topographic depression of the Lake Eyre Basin and Murray-Darling Basin. We propound that dynamic topography evolution provides an independent means to couple geological processes to a mantle reference frame. This is complementary to, and can be integrated with, other approaches such as hotspot and slab reference frames.

The occurrence of salt diapirs is strongly associated with potential geothermal and hydrocarbon energy sources. Many numerical modelling studies of diapirism have been done in the past, though very few of these in fact use geologically realistic settings and materials. Besides, only analogue and structural studies have been done on full scale diapirism during thin-skinned extension. Numerical modelling of this problem, using a code called ELEFANT, could answer the following questions: ”which parameters affect the growth rate and shape of the diapir, and how?” and ”what is the effect of this diapirism on the temperature field and surface heat flux?” The results show that, in compliance with both analogue modelling and structural geological studies, a diapir formed during thin-skinned extension undergoes three phases: reactional piercement, active piercement and passive piercement. Extension rates directly influence the total time required for the diapir to reach the surface, as well as how long the system remains in a state of reactional diapirism, which both affect the shape of the resulting diapir. The mass diffusivity coefficient affects the growth rate of the diapir during its active stage and the total rising time, which affects the shape of the diapir. The density contrast between the salt and the sediments also influences the growth rate during active and passive piercement. Surrounding a rising diapir the temperature is heightened, although not as much as expected. It can, however, still be valuable for hydrocarbon exploration. The heat flux above the diapir is higher than away from the diapir. Besides, both an increased bottom temperature and an increased salt conductivity give higher surface heat fluxes. This makes it an important feature for geothermal energy.

Realistic appraisal of paleoclimatic information obtained from a particular location requires accurate knowledge of its paleolatitude defined relative to the Earth’s spin-axis. This is crucial to, among others, correctly assess the amount of solar energy received at a location at the moment of sediment deposition. The paleolatitude of an arbitrary location can in principle be reconstructed from tectonic plate reconstructions that (1) restore the relative motions between plates based on (marine) magnetic anomalies, and (2) reconstruct all plates relative to the spin axis using a paleomagnetic reference frame based on a global apparent polar wander path. Whereas many studies do employ high-quality relative plate reconstructions, the necessity of using a paleomagnetic reference frame for climate studies rather than a mantle reference frame appears under-appreciated. In this paper, we briefly summarize the theory of plate tectonic reconstructions and their reference frames tailored towards applications of paleoclimate reconstruction, and show that using a mantle reference frame, which defines plate positions relative to the mantle, instead of a paleomagnetic reference frame may introduce errors in paleolatitude of more than 15° (>1500 km). This is because mantle reference frames cannot constrain, or are specifically corrected for the effects of true polar wander. We used the latest, state-of-the-art plate reconstructions to build a global plate circuit, and developed an online, user-friendly paleolatitude calculator for the last 200 million years by placing this plate circuit in three widely used global apparent polar wander paths. As a novelty, this calculator adds error bars to paleolatitude estimates that can be incorporated in climate modeling. The calculator is available at www.paleolatitude.org. We illustrate the use of the paleolatitude calculator by showing how an apparent wide spread in Eocene sea surface temperatures of southern high latitudes may be in part e

The level set method allows for tracking material surfaces in 2-D and 3-D flow modeling and is well suited for applications of multi-material flow modeling. The level set method utilizes smooth level set functions to define material interfaces, which makes the method stable and free of oscillations that are typically observed in case step-like functions parameterize interfaces. By design the level set function is a signed distance function and gives for each point in the domain the exact distance to the interface as well as on which side it is located. In this paper we present four benchmarks which show the validity, accuracy and simplicity of using the level set method for multi-material flow modeling. The benchmarks are simplified setups of dynamical geophysical processes such as the Rayleigh-Taylor instability, post-glacial rebound, subduction and slab detachment. We also demonstrate the benefit of using the level set method for modeling a free surface with the sticky air approach. Our results show that the level set method allows for accurate material flow modeling and that the combination with the sticky air approach works well in mimicking Earth's free surface. Since the level set method tracks material interfaces instead of materials themselves, it has the advantage that the location of these interfaces is accurately known and that it represents a viable alternative to the more commonly used tracer method.

The tightly curved Banda Arc (eastern Indonesia) plays a pivotal role in the convergent interaction of the Eurasian, Australian, and Pacific plates. Historically there have been two opposing hypotheses regarding its origin: either the arc has been formed by two opposite-dipping subduction zones joining at the apex, or a single slab has been wrought in its current position by surface kinematics and mantle processes. Using detailed tectonic reconstructions and mantle tomography, the latter view has recently gained substantial support by the work of Spakman and Hall (2010). Through numerical modelling of the instantaneous dynamics of the intricate subducting system, the research presented here aims to determine the factors that control the evolution of the single slab, but also to provide insight into the complexities of instantaneous numerical modelling of subduction in general. To aid the construction of very complex model geometries, we have developed a tool that allows us to digitally draw faults, slab sections, and plate outlines based on tomography, surface topography, and Wadati-Benioff zone geometry. The tool communicates with the finite element mesh generator, which then gives a composite mesh with highly desirable properties. Our results give a basic one-slab-model that satisfies both the predictions of the Spakman and Hall model and the stress state from CMT solutions and seismic fast directions. Better agreement can be attained by varying the buoyancy structure, to which the behavior of the model is found to be rather sensitive. Additionally, we propose that a mantle flow component is required in order to maintain the northern slab limb in the case that it has detached from the surface lithosphere. In the absence of flow-induced resistance, our results suggests that further convergent motion can effect a transition towards large-scale slab detachment progressing southward if there is little buoyant support at the d660 phase transition.

The Payún Matrú Volcanic Field (Pleistocene–Holocene) is located in the Andean back-arc of the Southern Volcanic Zone, western Argentina, and is contemporaneous with the Andean volcanic arc at the same latitude. It includes two polygenetic, mostly trachytic volcanoes: Payún Matrú (with a summit caldera 8 km wide) and Payún Liso (a smaller stratovolcano). The volcanic field includes about 200 scoria cones and alkali basaltic and trachybasaltic lava flows, forming two basaltic fields around Payún Matrú. New 40Ar–39Ar ages extend the activity of Payún Matrú up to 700 ka. The major and trace element and Sr–Nd isotopic compositions of the basaltic lavas and Payún Matrú rocks indicate that the trachytes of Payún Matrú are the result of fractional crystallization of basaltic parent magmas without significant upper crustal contamination, and that the basalts have a geochemical similarity to ocean island basalt (La/Nb = 0·8–1·5, La/Ba = 0·05–0·08). The Sr–Nd isotopic compositions of the basaltic to trachytic rocks range between 0·703813 and 0·703841 (87Sr/86Sr) and 0·512743 and 0·512834 (143Nd/144Nd). Mass-balance and Rayleigh fractionation models support the proposed origin of the trachytes, and an assimilation–fractional crystallization model indicates a low degree of upper crustal contamination in the youngest trachytes. Magnesium numbers (45–55) and contents of Ni (<20–90 ppm) and Cr (30–180 ppm) in the lavas in the basaltic fields indicate that these are not primary magmas. The data also suggest that the basaltic lavas originated in the asthenospheric mantle, probably within the spinel stability field and beneath an attenuated continental lithosphere in the back-arc area. The lack of a slab-fluid signature in the Payún Matrú Volcanic Field rocks, along with unpublished and published geophysical results (mantle tomography and electrical conductivity anomalies) suggest that magmas were generated by decompression-induced melting of upwelling mantle, The Payún Matrú Volcanic Field (Pleistocene–Holocene) is located in the Andean back-arc of the Southern Volcanic Zone, western Argentina, and is contemporaneous with the Andean volcanic arc at the same latitude. It includes two polygenetic, mostly trachytic volcanoes: Payún Matrú (with a summit caldera 8 km wide) and Payún Liso (a smaller stratovolcano). The volcanic field includes about 200 scoria cones and alkali basaltic and trachybasaltic lava flows, forming two basaltic fields around Payún Matrú. New 40Ar–39Ar ages extend the activity of Payún Matrú up to 700 ka. The major and trace element and Sr–Nd isotopic compositions of the basaltic lavas and Payún Matrú rocks indicate that the trachytes of Payún Matrú are the result of fractional crystallization of basaltic parent magmas without significant upper crustal contamination, and that the basalts have a geochemical similarity to ocean island basalt (La/Nb = 0·8–1·5, La/Ba = 0·05–0·08). The Sr–Nd isotopic compositions of the basaltic to trachytic rocks range between 0·703813 and 0·703841 (87Sr/86Sr) and 0·512743 and 0·512834 (143Nd/144Nd). Mass-balance and Rayleigh fractionation models support the proposed origin of the trachytes, and an assimilation–fractional crystallization model indicates a low degree of upper crustal contamination in the youngest trachytes. Magnesium numbers (45–55) and contents of Ni (<20–90 ppm) and Cr (30–180 ppm) in the lavas in the basaltic fields indicate that these are not primary magmas. The data also suggest that the basaltic lavas originated in the asthenospheric mantle, probably within the spinel stability field and beneath an attenuated continental lithosphere in the back-arc area. The lack of a slab-fluid signature in the Payún Matrú Volcanic Field rocks, along with unpublished and published geophysical results (mantle tomography and electrical conductivity anomalies) suggest that magmas were generated by decom

Cenozoic convergence between the Indian and Asian plates produced the archetypical continental collision zone comprising the Himalaya mountain belt and the Tibetan Plateau. How and where India–Asia convergence was accommodated after collision at or before 52 Ma remains a long-standing controversy. Since 52 Ma, the two plates have converged up to 3,600 35 km, yet the upper crustal shortening documented from the geological record of Asia and the Himalaya is up to approximately 2,350-km less. Here we show that the discrepancy between the convergence and the shortening can be explained by subduction of highly extended continental and oceanic Indian lithosphere within the Himalaya between approximately 50 and 25 Ma. Paleomagnetic data show that this extended continental and oceanic “Greater India” promontory resulted from 2,675 700 km of North–South extension between 120 and 70 Ma, accommodated between the Tibetan Himalaya and cratonic India. We suggest that the approximately 50 Ma “India”– Asia collision was a collision of a Tibetan-Himalayan microcontinent with Asia, followed by subduction of the largely oceanic Greater India Basin along a subduction zone at the location of the Greater Himalaya. The “hard” India–Asia collision with thicker and contiguous Indian continental lithosphere occurred around 25–20 Ma. This hard collision is coincident with far-field deformation in central Asia and rapid exhumation of Greater Himalaya crystalline rocks, and may be linked to intensification of the Asian monsoon system. This two-stage collision between India and Asia is also reflected in the deep mantle remnants of subduction imaged with seismic tomography.

The Southwest Pacific region is tectonically complex and is home to numerous fossil and active subduction zones. At the Earth's surface, there remains a geological controversy regarding the polarity and continuity of fossil subduction zones in New Zealand and New Caledonia, the origin of obducted ophiolites, the presence of high-pressure metamorphism, the occurrence of widespread Cenozoic magmatism, and the potential disappearance of one or more ocean basins. This controversy can be solved by looking at the lower mantle rather then at the Earth's surface. New P-wave and S-wave mantle tomography models from the Southwest Pacific are presented, which identify a previously unrecognized lower-mantle high-velocity anomaly that cannot be linked to Pacific subduction. The anomaly is located below the Tasman Sea at ~ 1100 km depth, strikes NW-SE and is ~ 2200 km by 600-900 km in lateral extent. By combining relative and absolute plate motions it is demonstrated that when the geological structures at the surface are reinterpreted as a single northeast-dipping 2500-km middle Cenozoic subduction zone (the so-called New Caledonia subduction zone) the lower mantle anomaly can be accounted for, as it is found at the predicted location and depth. Discovery of the lower mantle slab anomaly thereby solves a long-standing geological controversy in the New Zealand-New Caledonia region. Finally, reconstructions and analytical calculations predict a lower mantle slab sinking velocity of ~ 1.5 cm/yr and a lower mantle viscosity of ~ 1022 Pa·s.

The main aim of this study is to create a data set of accurate absolute arrival times for stations in Europe which do not report to the International Seismological Centre (ISC). Waveforms were obtained from data centres and temporary experiments and a semi-automatic picking method was applied to determine absolute arrival times for P and S phases. 85 000 arrival times were picked whose distribution of residuals shows generally low standard deviations on the order of 0.5–0.7 s. Furthermore, mean teleseismic station residuals reflect the properties of the underlying crust and uppermost mantle. Comparison to ISC data for matching event-station-phase combinations also confirms the good quality of the new absolute arrival time picks. Most importantly, this data set complements the ISC data as it fills regional data coverage gaps in Europe.

TOPO-EUROPE addresses the 4-D topographic evolution of the orogens and intra-plate regions of Europe through a multidisciplinary approach linking geology, geophysics, geodesy and geotechnology. TOPO-EUROPE integrates monitoring, imaging, reconstruction and modelling of the interplay between processes controlling continental topography and related natural hazards. Until now, research on neotectonics and related topography development of orogens and intra-plate regions has received little attention. TOPO-EUROPE initiates a number of novel studies on the quantification of rates of vertical motions, related tectonically controlled river evolution and land subsidence in carefully selected natural laboratories in Europe. From orogen through platform to continental margin, these natural laboratories include the Alps/Carpathians–Pannonian Basin System, the West and Central European Platform, the Apennines–Aegean–Anatolian region, the Iberian Peninsula, the Scandinavian Continental Margin, the East-European Platform, and the Caucasus–Levant area. TOPO-EUROPE integrates European research facilities and know-how essential to advance the understanding of the role of topography in Environmental Earth System Dynamics. The principal objective of the network is twofold. Namely, to integrate national research programs into a common European network and, furthermore, to integrate activities among TOPO-EUROPE institutes and participants. Key objectives are to provide an interdisciplinary forum to share knowledge and information in the field of the neotectonic and topographic evolution of Europe, Depto. de Geodinámica, Estratigrafía y Paleontología, Fac. de Ciencias Geológicas, TRUE, pub