Mantle convection modeling with viscoelastic/brittle lithosphere: Numerical methodology and plate tectonic modeling

The earth's tectonic plates are strong, viscoelastic shells which make up the outermost part of a thermally convecting, predominantly viscous layer. Brittle failure of the lithosphere occurs when stresses are high. In order to build a realistic simulation of the planet's evolution, the complete viscoelastic/brittle convection system needs to be considered. A particle-in-cell finite element method is demonstrated which can simulate very large deformation viscoelasticity with a strain-dependent yield stress. This is applied to a plate-deformation problem. Numerical accuracy is demonstrated relative to analytic benchmarks, and the characteristics of the method are discussed.

Relative motions of Africa, Iberia and Europe during Alpine orogeny

A revised kinematic model for the motions of Africa and Iberia relative to Europe since the Middle Jurassic is presented in order to provide boundary conditions for Alpine-Mediterranean reconstructions. These motions were calculated using up-to-date kinematic data predominantly based on magnetic isochrons in the Atlantic Ocean and published by various authors during the last 15 years. It is shown that convergence of Africa with respect to Europe commenced during the Cretaceous Normal Superchron (CNS), between chrons MO and 34 (120-83 Ma). This motion was subjected to fluctuations in convergence rates characterised by two periods of relatively rapid convergence (during Late Cretaceous and Eocene-Oligocene times) that alternated with periods of slower convergence (during the Paleocene and since the Early Miocene). Distinct changes in plate kinematics are recognised in the motion of Iberia with respect to Europe, indicated by: (1) a Late Jurassic-Early Cretaceous left-lateral strike-slip motion; (2) Late Cretaceous convergence; (3) Paleocene quiescence; (4) a short period of right-lateral strike-slip motion; and (5) final Eocene-Oligocene convergence. Based on these results, it is speculated that a collisional episode in the Alpine orogeny at ca. 65 Ma resulted in a dramatic decrease in the relative plate motions and that a slower motion since the Early Miocene promoted extension in the Mediterranean back-arc basins. (C) 2002 Elsevier Science B.V. All rights reserved.

Deconvolving the delta O-18 seawater component from subseasonal coral delta O-18 and Sr/Ca at Rarotonga in the southwestern subtropical Pacific for the period 1726 to 1997

To reconstruct oceanographic variations in the subtropical South Pacific, 271-year long subseasonal time series of Sr/Ca and delta(18)O were generated from a coral growing at Rarotonga (21.5degreesS, 159.5degreesW). In this case, coral Sr/Ca appears to be an excellent proxy for sea surface temperature (SST) and coral delta(18)O is a function of both SST and seawater delta(18)O composition (delta(18)O(sw)). Here, we focus on extracting the delta(18)O(sw) signal from these proxy records. A method is presented assuming that coral Sr/Ca is solely a function of SST and that coral delta(18)O is a function of both SST and delta(18)O(sw). This method separates the effects of delta(18)O(sw) from SST by breaking the instantaneous changes of coral delta(18)O into separate contributions by instantaneous SST and delta(18)O(sw) changes, respectively. The results show that on average delta(18)O(sw) at Rarotonga explains similar to39% of the variance in delta(18)O and that variations in SST explains the remaining similar to61% of delta(18)O variance. Reconstructed delta(18)O(sw) shows systematic increases in summer months (December-February) consistent with the regional pattern of variations in precipitation and evaporation. The delta(18)O(sw) also shows a positive linear correlation with satellite-derived estimated salinity for the period 1980 to 1997 (r = 0.72). This linear correlation between reconstructed delta(18)O(sw) and salinity makes it possible to use the reconstructed delta(18)O(sw) to estimate the past interannual and decadal salinity changes in this region. Comparisons of coral delta(18)O and delta(18)O(sw) at Rarotonga with the Pacific decadal oscillation index suggest that the decadal and interdecadal salinity and SST variability at Rarotonga appears to be related to basin-scale decadal variability in the Pacific. Copyright (C) 2002 Elsevier Science Ltd.

Longevity and stability of cratonic lithosphere: Insights from numerical simulations of coupled mantle convection and continental tectonics

[1] The physical conditions required to provide for the tectonic stability of cratonic crust and for the relative longevity of deep cratonic lithosphere within a dynamic, convecting mantle are explored through a suite of numerical simulations. The simulations allow chemically distinct continents to reside within the upper thermal boundary layer of a thermally convecting mantle layer. A rheologic formulation, which models both brittle and ductile behavior, is incorporated to allow for plate-like behavior and the associated subduction of oceanic lithosphere. Several mechanisms that may stabilize cratons are considered. The two most often invoked mechanisms, chemical buoyancy and/or high viscosity of cratonic root material, are found to be relatively ineffective if cratons come into contact with subduction zones. High root viscosity can provide for stability and longevity but only within a thick root limit in which the thickness of chemically distinct, high-viscosity cratonic lithosphere exceeds the thickness of old oceanic lithosphere by at least a factor of 2. This end-member implies a very thick mechanical lithosphere for cratons. A high brittle yield stress for cratonic lithosphere as a whole, relative to oceanic lithosphere, is found to be an effective and robust means for providing stability and lithospheric longevity. This mode does not require exceedingly deep strength within cratons. A high yield stress for only the crustal or mantle component of the cratonic lithosphere is found to be less effective as detachment zones can then form at the crust-mantle interface which decreases the longevity potential of cratonic roots. The degree of yield stress variations between cratonic and oceanic lithosphere required for stability and longevity can be decreased if cratons are bordered by continental lithosphere that has a relatively low yield stress, i.e., mobile belts. Simulations that combine all the mechanisms can lead to crustal stability and deep root longevity for model cratons over several mantle overturn times, but the dominant stabilizing factor remains a relatively high brittle yield stress for cratonic lithosphere.

Convective instability of 3-D fluid-saturated geological fault zones heated from below

We conduct a theoretical analysis to investigate the convective instability of 3-D fluid-saturated geological fault zones when they are heated uniformly from below. In particular, we have derived exact analytical solutions for the critical Rayleigh numbers of different convective flow structures. Using these critical Rayleigh numbers, three interesting convective flow structures have been identified in a geological fault zone system. It has been recognized that the critical Rayleigh numbers of the system have a minimum value only for the fault zone of infinite length, in which the corresponding convective flow structure is a 2-D slender-circle flow. However, if the length of the fault zone is finite, the convective flow in the system must be 3-D. Even if the length of the fault zone is infinite, since the minimum critical Rayleigh number for the 2-D slender-circle flow structure is so close to that for the 3-D convective flow structure, the system may have almost the same chance to pick up the 3-D convective flow structures. Also, because the convection modes are so close for the 3-D convective flow structures, the convective flow may evolve into the 3-D finger-like structures, especially for the case of the fault thickness to height ratio approaching zero. This understanding demonstrates the beautiful aspects of the present analytical solution for the convective instability of 3-D geological fault zones, because the present analytical solution is valid for any value of the ratio of the fault height to thickness. Using the present analytical solution, the conditions, under which different convective flow structures may take place, can be easily determined.

Neogene and Quaternary rollback evolution of the Tyrrhenian Sea, the Apennines, and the Sicilian Maghrebides

Reconstruction of the evolution of the Tyrrhenian Sea shows that the major stage of rifting associated with the opening of this basin began at similar to10 Ma. It involved two episodes of back arc extension, which were induced by the rollback of a west dipping subducting slab. The first period of extension (10-6 Ma) was prominent in the northern Tyrrhenian Sea and in the western part of the southern Tyrrhenian Sea. The second period of extension, mainly affected the southern Tyrrhenian Sea, began in the latest Messinian (6-5 Ma) and has been accompanied by subduction rollback at rates of 60-100 km Myr(-1). Slab reconstruction, combined with paleomagnetic and paleogeographic constraints, indicates that in the central Apennines, the latest Messinian (6-5 Ma) arrival of a carbonate platform at the subduction zone impeded subduction and initiated a slab tear and major strike-slip faults. These processes resulted in the formation of a narrow subducting slab beneath the Ionian Sea that has undergone faster subduction rollback and induced extreme rates of back arc extension.

5th International Symposium on Applied Isotope Geochemistry - Preface

Selected Papers from the 5th International Symposium on Applied Isotope Geochemistry, Heron Island, Great Barrier Reef, Australia, 26–30 May 2003

Subduction of the Nazca Ridge and the Inca Plateau: Insights into the formation of ore deposits in Peru

A large number of ore deposits that formed in the Peruvian Andes during the Miocene (15-5 Ma) are related to the subduction of the Nazea plate beneath the South American plate. Here we show that the spatial and temporal distribution of these deposits correspond with the arrival of relatively buoyant topographic anomalies, namely the Nazca Ridge in central Peru and the now-consumed Inca Plateau in northern Peru, at the subduction zone. Plate reconstruction shows a rapid metallogenic response to the arrival of the topographic anomalies at the subduction trench. This is indicated by clusters of ore deposits situated within the proximity of the laterally migrating zones of ridge subduction. It is accordingly suggested that tectonic changes associated with impingement of the aseismic ridge into the subduction zone may trigger the formation of ore deposits in metallogenically fertile suprasubduction environments. (c) 2005 Elsevier B.V All rights reserved.

Re-dating mid-Holocene betelnut (Areca catechu L.) and other plant use at Dongan, Papua New Guinea

Direct accelerator mass spectrometry (AMS) dating of anaerobically preserved plant remains from the Dongan site in New Guinea, combined with assessment of preservation condition, confirms earlier doubts about the antiquity of betel-nut (Areca catechu L.) found at the site. A possible sago leaf fragment is also identified as a modem contaminant. The mid-Holocene age of other fruit and nut remains is verified using these methods. The utility of AMS dating in combination with detailed archaeobotanical assessment is demonstrated, thus improving chronometric hygiene and with it knowledge of past plant use in Oceania.

Timing and duration of the last interglacial inferred from high resolution U-series chronology of stalagmite growth in Southern Hemisphere

High-precision Th-230-U-238 ages for a stalagmite from Newdegate Cave in southern Tasmania, Australia define a rare record of precipitation between 100 and 155 ka before the present. The fastest stalagmite growth occurred between 129.2 +/- 1.6 and 122.1 +/- 2.0 ka (similar to 61.5 mm/ka), coinciding with a time of prolific coral growth from Western Australia (128-122 ka). This is the first high-resolution continental record in the Southern Hemisphere that can be compared and correlated with the marine record. Such correlation shows that in southern Australia the onset of full interglacial sea level and the initiation of highest precipitation on land were synchronous. The stalagmite growth rate between 129.2 and 142.2 ka (similar to 5.9 mm/ka) was lower than that between 142.2 and 154.5 ka (similar to 18.7 mm/ka), implying drier conditions during the Penultimate Deglaciation, despite rising temperature and sea level. This asymmetrical precipitation pattern is caused by latitudinal movement of subtropical highs and an associated Westerly circulation, in response to a changing Equator-to-Pole temperature gradient. Both marine and continental records in Australia strongly suggest that the insolation maximum between 126 and 128 ka at 65 degreesN was directly responsible for the maintenance of full Last Interglacial conditions, although the triggers that initiated Penultimate Deglaciation (at similar to 142 ka) remain unsolved. (C) 2001 Elsevier Science B.V. All rights reserved.

Early-Mid Holocene climatic variations in Tasmania, Australia: multi-proxy records in a stalagmite from Lynds Cave

Mass spectrometric uranium-series dating and C-O isotopic analysis of a stalagmite from Lynds Cave, northern Tasmania, Australia provide a high-resolution record of regional climate change between 5100 and 9200 yr before present (BP). Combined delta(18)O, delta(13)C, growth rate, initial U-234/U-238 and physical property (color, transparency and porosity) records allow recognition of seven climatic stages: Stage I ( > 9080 yr BP) - a relatively dry period at the beginning of stalagmite growth evidenced by elevated U-234/U-238; Stage II (9080-8600 yr BP) - a period of unstable climate characterized by high-frequency variability in temperature and bio-productivity; Stage 111 (8600-8000 yr BP) - a period of stable and moderate precipitation and stable and high bio-productivity, with a continuously rising temperature; Stage IV (8000-7400 yr BP) - the warmest period with high evaporation and low effective precipitation (rainfall less evaporation); Stage V (7400-7000 yr BP) - the wettest period with highest stalagmite growth and enhanced but unstable bio-productivity; Stage VI (7000-6600 yr BP) - a period with a significantly reduced precipitation and bio-productivity without noticeable change in temperature; Stage VII (6600-5100 yr BP) - a period of lowest temperature and precipitation marking a significant climatic deterioration. Overall, the records suggest that the warmest climate occurred between 8000 and 7400 yr BP, followed by a wettest period between 7400 and 7000 yr BP. These are broadly correlated with the so-called 'Mid Holocene optimum' previously proposed using pollen and lake level records. However, the timing and resolution of the speleothem. record from Lynds Cave are significantly higher than in both the pollen and lake level records. This allows us to correlate the abrupt change in physical property, delta(18)O, delta(13)C, growth rate, and initial U-234/U-238 of the stalagmite at similar to8000 yr BP with a global climatic event at Early-Mid Holocene transition. (C) 2001 Elsevier Science B.V. All rights reserved.