Cenozoic climate change as a possible cause for the rise of the Andes

Causal links between the rise of a large mountain range and climate have often been considered to work in one direction, with significant uplift provoking climate change. Here we propose a mechanism by which Cenozoic climate change could have caused the rise of the Andes. Based on considerations of the force balance in the South American lithosphere, we suggest that the height of, and tectonics in, the Andes are strongly controlled both by shear stresses along the plate interface in the subduction zone and by buoyancy stress contrasts between the trench and highlands, and shear stresses in the subduction zone depend on the amount of subducted sediments. We propose that the dynamics of subduction and mountain-building in this region are controlled by the processes of erosion and sediment deposition, and ultimately climate. In central South America, climate-controlled sediment starvation would then cause high shear stress, focusing the plate boundary stresses that support the high Andes.     

Measuring the onset of locking in the Peru-Chile trench with GPS and acoustic measurements

The subduction zone off the west coast of South America marks the convergence of the oceanic Nazca plate and the continental South America plate. Nazca-South America convergence over the past 23 million years has created the 6-km-deep Peru-Chile trench, 150 km offshore. High pressure between the plates creates a locked zone, leading to deformation of the overriding plate. The surface area of this locked zone is thought to control the magnitude of co-seismic release and is limited by pressure, temperature, sediment type and fluid content. Here we present seafloor deformation data from the submerged South America plate obtained from a combination of Global Positioning System (GPS) receivers and acoustic transponders. We estimate that the measured horizontal surface motion perpendicular to the trench is consistent with a model having no slip along the thrust fault between 2 and 40 km depth. A tsunami in 1996, 200 km north of our site, was interpreted as being the result of an anomalously shallow interplate earthquake. Seismic coupling at shallow depths, such as we observe, may explain why co-seismic events in the Peruvian subduction zone create large tsunamis.     

Evolution and diversity of subduction zones controlled by slab width

Subducting slabs provide the main driving force for plate motion and flow in the Earth's mantle, and geodynamic, seismic and geochemical studies offer insight into slab dynamics and subduction-induced flow. Most previous geodynamic studies treat subduction zones as either infinite in trench-parallel extent (that is, two-dimensional) or finite in width but fixed in space. Subduction zones and their associated slabs are, however, limited in lateral extent (250-7,400 km) and their three-dimensional geometry evolves over time. Here we show that slab width controls two first-order features of plate tectonics-the curvature of subduction zones and their tendency to retreat backwards with time. Using three-dimensional numerical simulations of free subduction, we show that trench migration rate is inversely related to slab width and depends on proximity to a lateral slab edge. These results are consistent with retreat velocities observed globally, with maximum velocities (6-16 cm yr(-1)) only observed close to slab edges (<1,200 km), whereas far from edges (>2,000 km) retreat velocities are always slow (<2.0 cm yr(-1)). Models with narrow slabs (< or =1,500 km) retreat fast and develop a curved geometry, concave towards the mantle wedge side. Models with slabs intermediate in width ( approximately 2,000-3,000 km) are sublinear and retreat more slowly. Models with wide slabs (> or =4,000 km) are nearly stationary in the centre and develop a convex geometry, whereas trench retreat increases towards concave-shaped edges. Additionally, we identify periods (5-10 Myr) of slow trench advance at the centre of wide slabs. Such wide-slab behaviour may explain mountain building in the central Andes, as being a consequence of its tectonic setting, far from slab edges.     

Seismic evidence for catastrophic slab loss beneath Kamchatka

In the northwest Pacific Ocean, a sharp corner in the boundary between the Pacific plate and the North American plate joins a subduction zone running along the southern half of the Kamchatka peninsula with a region of transcurrent motion along the western Aleutian arc. Here we present images of the seismic structure beneath the Aleutian-Kamchatka junction and the surrounding region, indicating that: the subducting Pacific lithosphere terminates at the Aleutian-Kamchatka junction; no relict slab underlies the extinct northern Kamchatka volcanic arc; and the upper mantle beneath northern Kamchatka has unusually slow shear wavespeeds. From the tectonic and volcanic evolution of Kamchatka over the past 10 Myr (refs 3-5) we infer that at least two episodes of catastrophic slab loss have occurred. About 5 to 10 Myr ago, catastrophic slab loss shut down island-arc volcanic activity north of the Aleutian-Kamchatka junction. A later episode of slab loss, since about 2 Myr ago, seems to be related to the activity of the world's most productive island-arc volcano, Klyuchevskoy. Removal of lithospheric mantle is commonly discussed in the context of a continental collision, but our findings imply that episodes of slab detachment and loss are also important agents in the evolution of oceanic convergent margins.     

Rapid change in drift of the Australian plate records collision with Ontong Java plateau

The subduction of oceanic plateaux, which contain extraordinarily thick basaltic crust and are the marine counterparts of continental flood-basalt provinces, is an important factor in many current models of plate motion and provides a potential mechanism for triggering plate reorganization. To evaluate such models, it is essential to decipher the history of the collision between the largest and thickest of the world's oceanic plateaux, the Ontong Java plateau, and the Australian plate, but this has been hindered by poor constraints for the arrival of the plateau at the Melanesian trench. Here we present (40)Ar-(39)Ar geochronological data on hotspot volcanoes in eastern Australian that reveal a strong link between collision of the Greenland-sized Ontong Java plateau with the Melanesian arc and motion of the Australian plate. The new ages define a short-lived period of reduced northward plate motion between 26 and 23 Myr ago, coincident with an eastward offset in the contemporaneous tracks of seamount chains in the Tasman Sea east of Australia. These features record a brief westward deflection of the Australian plate as the plateau entered and choked the Melanesian trench 26 Myr ago. From 23 Myr ago, Australia returned to a rapid northerly trajectory at roughly the same time that southwest-directed subduction began along the Trobriand trough. The timing and brevity of this collisional event correlate well with offsets in hotspot seamount tracks on the Pacific plate, including the archetypal Hawaiian chain, and thus provide strong evidence that immense oceanic plateaux, like the Ontong Java, can contribute to initiating rapid change in plate boundaries and motions on a global scale.     

Evidence of lower-mantle slab penetration phases in plate motions

It is well accepted that subduction of the cold lithosphere is a crucial component of the Earth's plate tectonic style of mantle convection. But whether and how subducting plates penetrate into the lower mantle is the subject of continuing debate, which has substantial implications for the chemical and thermal evolution of the mantle. Here we identify lower-mantle slab penetration events by comparing Cenozoic plate motions at the Earth's main subduction zones with motions predicted by fully dynamic models of the upper-mantle phase of subduction, driven solely by downgoing plate density. Whereas subduction of older, intrinsically denser, lithosphere occurs at rates consistent with the model, younger lithosphere (of ages less than about 60 Myr) often subducts up to two times faster, while trench motions are very low. We conclude that the most likely explanation is that older lithosphere, subducting under significant trench retreat, tends to lie down flat above the transition to the high-viscosity lower mantle, whereas younger lithosphere, which is less able to drive trench retreat and deforms more readily, buckles and thickens. Slab thickening enhances buoyancy (volume times density) and thereby Stokes sinking velocity, thus facilitating fast lower-mantle penetration. Such an interpretation is consistent with seismic images of the distribution of subducted material in upper and lower mantle. Thus we identify a direct expression of time-dependent flow between the upper and lower mantle.     

Trench-parallel flow and seismic anisotropy in the Mariana and Andean subduction systems

Shear-wave splitting measurements above the mantle wedge of the Mariana and southern Andean subduction zones show trench-parallel seismically fast directions close to the trench and abrupt rotations to trench-perpendicular anisotropy in the back arc. These patterns of seismic anisotropy may be caused by three-dimensional flow associated with along-strike variations in slab geometry. The Mariana and Andean subduction systems are associated with the largest along-strike variations of slab geometry observed on Earth and are ideal for testing the link between slab geometry and solid-state creep processes in the mantle. Here we show, with fully three-dimensional non-newtonian subduction zone models, that the strong curvature of the Mariana slab and the transition to shallow slab dip in the Southern Andes give rise to strong trench-parallel stretching in the warm-arc and warm-back-arc mantle and to abrupt rotations in stretching directions that are accompanied by strong trench-parallel stretching. These models show that the patterns of shear-wave splitting observed in the Mariana and southern Andean systems may be caused by significant three-dimensional flow induced by along-strike variations in slab geometry.     

Low electrical resistivity associated with plunging of the Nazca flat slab beneath Argentina

Beneath much of the Andes, oceanic lithosphere descends eastward into the mantle at an angle of about 30 degrees (ref. 1). A partially molten region is thought to form in a wedge between this descending slab and the overlying continental lithosphere as volatiles given off by the slab lower the melting temperature of mantle material. This wedge is the ultimate source for magma erupted at the active volcanoes that characterize the Andean margin. But between 28 degrees and 33 degrees S the subducted Nazca plate appears to be anomalously buoyant, as it levels out at about 100 km depth and extends nearly horizontally under the continent. Above this 'flat slab', volcanic activity in the main Andean Cordillera terminated about 9 million years ago as the flattening slab presumably squeezed out the mantle wedge. But it is unknown where slab volatiles go once this happens, and why the flat slab finally rolls over to descend steeply into the mantle 600 km further eastward. Here we present results from a magnetotelluric profile in central Argentina, from which we infer enhanced electrical conductivity along the eastern side of the plunging slab, indicative of the presence of partial melt. This conductivity structure may imply that partial melting occurs to at least 250 km and perhaps to more than 400 km depth, or that melt is supplied from the 410 km discontinuity, consistent with the transition-zone 'water-filter' model of Bercovici and Karato.     

Subduction and collision processes in the Central Andes constrained by converted seismic phases

The Central Andes are the Earth's highest mountain belt formed by ocean-continent collision. Most of this uplift is thought to have occurred in the past 20 Myr, owing mainly to thickening of the continental crust, dominated by tectonic shortening. Here we use P-to-S (compressional-to-shear) converted teleseismic waves observed on several temporary networks in the Central Andes to image the deep structure associated with these tectonic processes. We find that the Moho (the Mohorovici? discontinuity--generally thought to separate crust from mantle) ranges from a depth of 75 km under the Altiplano plateau to 50 km beneath the 4-km-high Puna plateau. This relatively thin crust below such a high-elevation region indicates that thinning of the lithospheric mantle may have contributed to the uplift of the Puna plateau. We have also imaged the subducted crust of the Nazca oceanic plate down to 120 km depth, where it becomes invisible to converted teleseismic waves, probably owing to completion of the gabbro-eclogite transformation; this is direct evidence for the presence of kinetically delayed metamorphic reactions in subducting plates. Most of the intermediate-depth seismicity in the subducting plate stops at 120 km depth as well, suggesting a relation with this transformation. We see an intracrustal low-velocity zone, 10-20 km thick, below the entire Altiplano and Puna plateaux, which we interpret as a zone of continuing metamorphism and partial melting that decouples upper-crustal imbrication from lower-crustal thickening.     

Late Miocene decoupling of oceanic warmth and atmospheric carbon dioxide forcing

Deep-time palaeoclimate studies are vitally important for developing a complete understanding of climate responses to changes in the atmospheric carbon dioxide concentration (that is, the atmospheric partial pressure of CO(2), p(co(2))). Although past studies have explored these responses during portions of the Cenozoic era (the most recent 65.5 million years (Myr) of Earth history), comparatively little is known about the climate of the late Miocene (～12-5 Myr ago), an interval with p(co(2)) values of only 200-350?parts per million by volume but nearly ice-free conditions in the Northern Hemisphere and warmer-than-modern temperatures on the continents. Here we present quantitative geochemical sea surface temperature estimates from the Miocene mid-latitude North Pacific Ocean, and show that oceanic warmth persisted throughout the interval of low p(co(2)) ～12-5 Myr ago. We also present new stable isotope measurements from the western equatorial Pacific that, in conjunction with previously published data, reveal a long-term trend of thermocline shoaling in the equatorial Pacific since ～13?Myr ago. We propose that a relatively deep global thermocline, reductions in low-latitude gradients in sea surface temperature, and cloud and water vapour feedbacks may help to explain the warmth of the late Miocene. Additional shoaling of the thermocline after 5?Myr ago probably explains the stronger coupling between p(co(2)), sea surface temperatures and climate that is characteristic of the more recent Pliocene and Pleistocene epochs.     

Laboratory models of the thermal evolution of the mantle during rollback subduction

The subduction of oceanic lithosphere plays a key role in plate tectonics, the thermal evolution of the mantle and recycling processes between Earth's interior and surface. Information on mantle flow, thermal conditions and chemical transport in subduction zones come from the geochemistry of arc volcanoes, seismic images and geodynamic models. The majority of this work considers subduction as a two-dimensional process, assuming limited variability in the direction parallel to the trench. In contrast, observationally based models increasingly appeal to three-dimensional flow associated with trench migration and the sinking of oceanic plates with a translational component of motion (rollback). Here we report results from laboratory experiments that reveal fundamental differences in three-dimensional mantle circulation and temperature structure in response to subduction with and without a rollback component. Without rollback motion, flow in the mantle wedge is sluggish, there is no mass flux around the plate and plate edges heat up faster than plate centres. In contrast, during rollback subduction flow is driven around and beneath the sinking plate, velocities increase within the mantle wedge and are focused towards the centre of the plate, and the surface of the plate heats more along the centreline.     

福建平潭岛火山岩及其板块碰撞构造

提出在福建东南沿海地区晚白垩世除了产出有拉张地球动力环境下火山（石帽山群火山岩）和深成（魁岐碱性花岗岩类）产物外，在长乐－南澳断裂带以东的平潭岛还产出有挤压地球动力环境下火山产物．其化学成分和岩系演化趋势显著不同于闽东南南园组和石帽山群火山岩，为独立的钙碱性角闪石英安岩系．其成因可能与库拉板块沿台湾中央山脉东侧向台湾海峡陆块的俯冲有关．闽东南地区早白垩世挤压环境火山产物（南园组火山岩）和晚自垩世拉张环境火山－深层岩系则可能与燕山早期台湾海峡陆块向闽东南大陆边缘的碰撞俯冲，而在晚白垩世停止俯冲，同时在台湾中央山脉东侧发生的俯冲作用所造成的“弧后扩张”有关．     

Deep roots of the Messinian salinity crisis

The Messinian salinity crisis--the desiccation of the Mediterranean Sea between 5.96 and 5.33 million years (Myr) ago--was one of the most dramatic events on Earth during the Cenozoic era. It resulted from the closure of marine gateways between the Atlantic Ocean and the Mediterranean Sea, the causes of which remain enigmatic. Here we use the age and composition of volcanic rocks to reconstruct the geodynamic evolution of the westernmost Mediterranean from the Middle Miocene epoch to the Pleistocene epoch (about 12.1-0.65 Myr ago). Our data show that a marked shift in the geochemistry of mantle-derived volcanic rocks, reflecting a change from subduction-related to intraplate-type volcanism, occurred between 6.3 and 4.8 Myr ago, largely synchronous with the Messinian salinity crisis. Using a thermomechanical model, we show that westward roll back of subducted Tethys oceanic lithosphere and associated asthenospheric upwelling provides a plausible mechanism for producing the shift in magma chemistry and the necessary uplift (approximately 1 km) along the African and Iberian continental margins to close the Miocene marine gateways, thereby causing the Messinian salinity crisis.     

Subduction erosion along the Middle America convergent margin

'Subduction erosion' has been invoked to explain material missing from some continents along convergent margins. It has been suggested that this form of tectonic erosion removes continental material at the front of the margin or along the underside of the upper (continental) plate. Frontal erosion is interpreted from disrupted topography at the base of a slope and is most evident in the wake of subducting seamounts. In contrast, structures resulting from erosion at the base of a continental plate are seldom recognized in seismic reflection images because such images typically have poor resolution at distances greater than approximately 5 km from the trench axis. Basal erosion from seamounts and ridges has been inferred, but few large subducted bodies--let alone the eroded base of the upper plate--are imaged convincingly. From seismic images we identify here two mechanisms of basal erosion: erosion by seamount tunnelling and removal of large rock lenses of a distending upper plate. Seismic cross-sections from Costa Rica to Nicaragua indicate that erosion may extend along much of the Middle America convergent margin.     

Unusually large earthquakes inferred from tsunami deposits along the Kuril trench

The Pacific plate converges with northeastern Eurasia at a rate of 8-9 m per century along the Kamchatka, Kuril and Japan trenches. Along the southern Kuril trench, which faces the Japanese island of Hokkaido, this fast subduction has recurrently generated earthquakes with magnitudes of up to approximately 8 over the past two centuries. These historical events, on rupture segments 100-200 km long, have been considered characteristic of Hokkaido's plate-boundary earthquakes. But here we use deposits of prehistoric tsunamis to infer the infrequent occurrence of larger earthquakes generated from longer ruptures. Many of these tsunami deposits form sheets of sand that extend kilometres inland from the deposits of historical tsunamis. Stratigraphic series of extensive sand sheets, intercalated with dated volcanic-ash layers, show that such unusually large tsunamis occurred about every 500 years on average over the past 2,000-7,000 years, most recently approximately 350 years ago. Numerical simulations of these tsunamis are best explained by earthquakes that individually rupture multiple segments along the southern Kuril trench. We infer that such multi-segment earthquakes persistently recur among a larger number of single-segment events.     

A satellite geodetic survey of large-scale deformation of volcanic centres in the central Andes

Surface deformation in volcanic areas usually indicates movement of magma or hydrothermal fluids at depth. Stratovolcanoes tend to exhibit a complex relationship between deformation and eruptive behaviour. The characteristically long time spans between such eruptions requires a long time series of observations to determine whether deformation without an eruption is common at a given edifice. Such studies, however, are logistically difficult to carry out in most volcanic arcs, as these tend to be remote regions with large numbers of volcanoes (hundreds to even thousands). Here we present a satellite-based interferometric synthetic aperture radar (InSAR) survey of the remote central Andes volcanic arc, a region formed by subduction of the Nazca oceanic plate beneath continental South America. Spanning the years 1992 to 2000, our survey reveals the background level of activity of about 900 volcanoes, 50 of which have been classified as potentially active. We find four centres of broad (tens of kilometres wide), roughly axisymmetric surface deformation. None of these centres are at volcanoes currently classified as potentially active, although two lie within about 10 km of volcanoes with known activity. Source depths inferred from the patterns of deformation lie between 5 and 17 km. In contrast to the four new sources found, we do not observe any deformation associated with recent eruptions of Lascar, Chile.     

Identification of Hercynian shoshonitic intrusive rocks in central Hainan Island and its geotectonic implications

Hercynian (Variscan) orogenic belts, including the European-NW African orogen, the Appalachian orogen of North America and the central Asia-Mongolia- Hinggan orogen, etc., are widely distributed in the world. Their extensions are often several thousand ki…     

Seismic tomography of the northwest Pacific and its geodynamic implications

High-resolution tomographic images across Japan Trenh-Changhai Mountains-lDong Ujimqinqi are displayed, showing the morphological feature of the subducted slab in the norhwestem Pacific margin and the eharaeter istics of lithosphere stmctures under the Changhai Mountains and the Da Hinggan Mnuntains. The Pacific plate began to penetrate into the deeper mantle after it subducted to the 660 km discontinuity with an underthmsting angle of 26°, but did not continue to mnve furrther westward. In contrast, there appeared a remarkable thermal upwelling zone to the west of the downward plate. In addition, the evidence frnm the subduction time and time lag between the subduetion and eon sequent magmatism indicates that there is no direct genetic correlatiom between the Mesoznic magmatism in eastern China ami subduction of the Pacific plate. In this work. we also emphasize that what the tomographic images reflect is the pre sent structure in the deep earth interior, which should preserve some Mesozoic lithospheric structure characteristics. In summary, we attribute the Mesozoic intense magmatic evolution in north China to the intraplate asthenosphere upwelling zone.     

The Cenozoic palaeoenvironment of the Arctic Ocean

The history of the Arctic Ocean during the Cenozoic era (0-65 million years ago) is largely unknown from direct evidence. Here we present a Cenozoic palaeoceanographic record constructed from >400 m of sediment core from a recent drilling expedition to the Lomonosov ridge in the Arctic Ocean. Our record shows a palaeoenvironmental transition from a warm 'greenhouse' world, during the late Palaeocene and early Eocene epochs, to a colder 'icehouse' world influenced by sea ice and icebergs from the middle Eocene epoch to the present. For the most recent approximately 14 Myr, we find sedimentation rates of 1-2 cm per thousand years, in stark contrast to the substantially lower rates proposed in earlier studies; this record of the Neogene reveals cooling of the Arctic that was synchronous with the expansion of Greenland ice (approximately 3.2 Myr ago) and East Antarctic ice (approximately 14 Myr ago). We find evidence for the first occurrence of ice-rafted debris in the middle Eocene epoch (approximately 45 Myr ago), some 35 Myr earlier than previously thought; fresh surface waters were present at approximately 49 Myr ago, before the onset of ice-rafted debris. Also, the temperatures of surface waters during the Palaeocene/Eocene thermal maximum (approximately 55 Myr ago) appear to have been substantially warmer than previously estimated. The revised timing of the earliest Arctic cooling events coincides with those from Antarctica, supporting arguments for bipolar symmetry in climate change.     

Orbitally induced oscillations in the East Antarctic ice sheet at the Oligocene/Miocene boundary

Between 34 and 15 million years (Myr) ago, when planetary temperatures were 3-4 degrees C warmer than at present and atmospheric CO2 concentrations were twice as high as today, the Antarctic ice sheets may have been unstable. Oxygen isotope records from deep-sea sediment cores suggest that during this time fluctuations in global temperatures and high-latitude continental ice volumes were influenced by orbital cycles. But it has hitherto not been possible to calibrate the inferred changes in ice volume with direct evidence for oscillations of the Antarctic ice sheets. Here we present sediment data from shallow marine cores in the western Ross Sea that exhibit well dated cyclic variations, and which link the extent of the East Antarctic ice sheet directly to orbital cycles during the Oligocene/Miocene transition (24.1-23.7 Myr ago). Three rapidly deposited glacimarine sequences are constrained to a period of less than 450 kyr by our age model, suggesting that orbital influences at the frequencies of obliquity (40 kyr) and eccentricity (125 kyr) controlled the oscillations of the ice margin at that time. An erosional hiatus covering 250 kyr provides direct evidence for a major episode of global cooling and ice-sheet expansion about 23.7 Myr ago, which had previously been inferred from oxygen isotope data (Mi1 event).