Geological transect across the Northwestern Himalaya in eastern Ladakh and Lahul (A model for the continental collision of India and Asia)

The detailed geological mapping and structural study of a complete transect across the northwestern Himalaya allow to describe the tectonic evolution of the north Indian continental margin during the Tethys ocean opening and the Himalayan Orogeny. The Late Paleozoic Tethys rifting is associated with several tectonomagmatic events. In Upper Lahul and SE Zanskar, this extensional phase is recorded by Lower Carboniferous synsedimentary transtensional faults, a Lower Permian stratigraphic unconformity, a Lower Permian granitic intrusion and middle Permian basaltic extrusions (Panjal Traps). In eastern Ladakh, a Permian listric normal fault is also related to this phase. The scarcity of synsedimentary faults and the gradual increase of the Permian syn-rift sediment thickness towards the NE suggest a flexural type margin. The collision of India and Asia is characterized by a succession of contrasting orogenic phases. South of the Suture Zone, the initiation of the SW vergent Nyimaling-Tsarap Nappe corresponds to an early phase of continental underthrusting. To the S, in Lahul, an opposite underthrusting within the Indian plate is recorded by the NE vergent Tandi Syncline. This structure is associated with the newly defined Shikar Beh Nappe, now partly eroded, which is responsible for the high grade (amphibolite facies) regional metamorphism of South Lahul. The main thrusting of the Nyimaling-Tsarap Nappe followed the formation of the Shikar Beh Nappe. The Nyimaling-Tsarap Nappe developed by ductile shear of the upper part of the subducted Indian continental margin and is responsible for the progressive regional metamorphism of SE Zanskar, reaching amphibolite facies below the frontal part of the nappe, near Sarchu. In Upper Lahul, the frontal parts of the Nyimaling-Tsarap and Shikar Beh nappes are separated by a zone of low grade metamorphic rocks (pumpellyite-actinolite facies to lower greenschist facies). At high structural level, the Nyimaling-Tsarap Nappe is characterized by imbricate structures, which grade into a large ductile shear zone with depth. The related crustal shortening is about 87 km. The root zone and the frontal part of this nappe have been subsequently affected by two zones of dextral transpression and underthrusting: the Nyimaling Shear Zone and the Sarchu Shear Zone. These shear zones are interpreted as consequences of the counterclockwise rotation of the continental underthrusting direction of India relative to Asia, which occurred some 45 and 36 Ma ago, according to plate tectonic models. Later, a phase of NE vergent `'backfolding'' developed on these two zones of dextral transpression, creating isoclinal folds in SE Zanskar and more open folds in the Nyimaling Dome and in the Indus Molasse sediments. During a late stage of the Himalayan Orogeny, the frontal part of the Nyimaling-Tsarap Nappe underwent an extension of about 15 km. This phase is represented by two types of structures, responsible for the tectonic unroofing of the amphibolite facies rocks of the Sarchu area: the Sarchu high angle Normal Fault, cutting a first set of low angle normal faults, which have been created by reactivation of older thrust planes related to the Nyimaling-Tsarap Nappe.

Radiometric age constraints on mineral growth, metamorphism, and tectonism of the Gummfluh klippe, Brianconnais domain of the Prealpes, Switzerland

A combined Ar-40/Ar-39, K/Ar, Rb/Sr and stable isotope study has been made of white micas from the Gummfluh klippe (Brianconnais domain of the Prealpes), Switzerland. The klippe consists mainly of Mesozoic to early Tertiary carbonate rocks metamorphosed from anchizonal to epizonal conditions. At the base of the klippe is a 10-50 m thick, ductilely deformed marble mylonite containing deformed authigenic quartz segregations. Stable isotope measurements of the coexisting calcite (deltaO-18SMOW=24.5) and quartz (deltaO-18SMOW=28.4) from the mylonite indicate relatively low temperatures (< 300-degreesC) during mylonitization. Analyses of white mica separates of varying size fractions from the mylonitic rocks by K/Ar and Rb/Sr techniques yield ages between 57 and 103 Ma. This variation is correlated with two parameters, the size of the mineral fraction, and the proportion of 2M1 (more phengitic) to 1M (more muscovitic) polytype in the sample. The K/Ar and Rb/Sr ages are generally younger in the smaller size fractions, which also containless 2M1 phengite. High precision Ar-40/Ar-39 age spectra from different size fractions of these micas record three distinct components, a small Hercynian component (ca. 200-300 Ma), a significant Eoalpine component (64-80 Ma) forming Ar-40/Ar-39 age plateaus, and a very minor Tertiary component (ca. 20-40 Ma). Characterization of the samples by SEM indicates the presence of two white mica populations, a coarser grained, deformed, detrital mica that probably corresponds to the 2M1 phengite and a finer grained neoformed 1M mica. Collectively these observations suggest that the Gummfluh samples contain a mixture of detrital phengites of Hercynian age together with neocrystallized muscovites grown during the late Eoalpine metamorphic event followed by minor argon loss during the Tertiary. The main geologic episode recorded in the Ar-40/Ar-39 age spectra of white micas in the mylonite is of Late Cretaceous/Early Tertiary age (64-80 Ma), representing the first reliable Eoalpine ages ever to be reported from the Prealpes. Contrary to tectonic models, the marble mylonite at the base of the Gummfluh klippe appears to be a Cretaceous thrust plane and not the thrust surface formed during transport of the klippe into its present position from the Penninic Alps during the Tertiary. The late Cretaceous thrust developed during marine sedimentation at a depth of 800 m below the seafloor at temperatures of approximately 280-degrees-C, facilitated by warm fluids along the tectonic discontinuity.

Late Cretaceous to Eocene geology of the South Central American forearc area (southern Costa Rica and western Panama) : initiation and evolution of an intra-oceanic convergent margin

Résumé : L'arc volcanique du sud de l'Amérique Centrale se situe sur la marge SW de la Plaque Caraïbe, au-dessus des plaques subduites de Cocos et Nazca. Il s'agit de l'un des arcs intra-océaniques les plus étudiés au monde, qui est généralement considéré comme s'étant développé à la fin du Crétacé le long d'un plateau océanique (le Plateau Caraïbe ou CLIP) et se trouvant actuellement dans un régime de subduction érosive. Au cours des dernières décennies, des efforts particuliers ont été faits pour comprendre les processus liés à la subduction sur la base d'études géophysiques et géochimiques. Au sud du Costa Rica et à l'ouest du Panama, des complexes d'accrétions et structures à la base de l'arc volcanique ont été exposés grâce à la subduction de rides asismiques et de failles transformantes. Des affleurements, situés jusqu'à seulement 15 km de la fosse, offrent une possibilité unique de mieux comprendre quelques uns des processus ayant lieu le long de la zone de subduction. Nous présentons de nouvelles contraintes sur l'origine de ces affleurements en alliant une étude de terrain poussée, de nouvelles données géochimiques, sédimentaires et paléontologiques, ainsi que des observations structurales effectuées en télédétection. Une nouvelle stratigraphie tectonique entre le Campanien et l'Éocène est définie pour la région d'avant-arc située entre la Péninsule d'Osa (Costa Rica) et la Péninsule d'Azuero (Panama). Nos résultats montrent que la partie externe de la marge est composée d'un arrangement complexe de roches ignées et de séquences sédimentaires de recouvrement qui comprennent principalement le socle de l'arc, des roches d'arc primitif, des fragments de monts sous-marins accrétés et des mélanges d'accrétion. Des preuves sont données pour le développement de l'arc volcanique du sud de l'Amérique Centrale sur un plateau océanique. Le début de la subduction le long de la marge SW de la Plaque Caraïbe a eu lieu au Campanien et a généré des roches d'arc primitif caractérisées par des affinités géochimiques particulières, globalement intermédiaires entre des affinités de plateau et d'arc insulaire. L'arc était mature au Maastrichtien et formait un isthme essentiellement continu entre l'Amérique du Nord et l'Amérique du Sud. Ceci a permis la migration de faunes terrestres entre les Amériques et pourrait avoir contribué à la crise fin Crétacé -Tertiaire en réduisant les courants océaniques subéquatoriaux entre le Pacifique et l'Atlantique. Plusieurs unités composées de fragments de monts sous-marins accrétés sont définies. La nature et l'arrangement structural de ces unités définissent de nouvelles contraintes sur les modes d'accrétion des monts sous-marins/îles océaniques et sur l'évolution de la marge depuis la formation de la zone de subduction. Entre la fin du Crétacé et l'Éocène moyen, la marge a enregistré plusieurs épisodes ponctuels d'accrétion de monts sous-marins alternant avec de la subduction érosive. A l'Éocène moyen, un événement tectonique régional pourrait avoir causé un fort couplage entre les plaques supérieure et inférieure, menant à des taux plus important d'accrétion de monts sous-marins. Durant cette période, la situation le long de la marge était très semblable à la situation actuelle et caractérisée par la présence de monts sous-marins subductants et l'absence d'accrétion de sédiments. L'enregistrement géologique montre qu'il n'est pas possible d'attribuer une nature érosive ou accrétionnaire à la marge dans le passé ou -par analogie- aujourd'hui, parce que (1) les processus d'accrétion et érosifs varient fortement spatialement et temporellement et (2) il est impossible d'évaluer la quantité exacte de matériel tectoniquement enlevé à la marge depuis le début de la subduction. Au sud du Costa Rica, certains fragments de monts sous-marins accrétés sont représentatifs d'une interaction entre une ride et un point chaud dans le Pacifique au Crétacé terminal/Paléocène. L'existence de ces fragments de monts sous-marins et la morphologie du fond de l'Océan Pacifique indiquent que la formation de la ride de Cocos-Nazca s'est formée au moins ~40 Ma avant l'âge proposé par les modèles tectoniques actuels. Au Panama, nous avons identifié une île océanique d'âge début Éocène qui a été accrétée à l'Éocène moyen. L'accrétion a eu lieu à très faible profondeur par détachement de l'île dans la fosse, et a mené à une exceptionnelle préservation des structures volcaniques. Des affleurement comprenant aussi bien des parties basses et hautes de l'édifice volcanique on été étudiées, depuis la phase sous-marine bouclier jusqu'à la phase subaérienne post-bouclier. La stratigraphie nous a permis de différencier les laves de la phase sous-marine de celles de la phase subaérienne. La composition des laves indique une diminution progressive de l'intensité de la fusion partielle de la source et une diminution de la température des laves produites durant les derniers stades de l'activité volcanique. Nous interprétons ces changements comme étant liés à l'éloignement progressif de l'île océanique de la zone de fusion ou point chaud. Abstract The southern Central American volcanic front lies on the SW edge of the Caribbean Plate, inboard of the subducting Cocos and Nazca Plates. It is one of the most studied intra-oceanic convergent margins around the world, which is generally interpreted to have developed in the late Cretaceous along an oceanic plateau (the Caribbean Large Igneous Province or CLIP) and to be currently undergoing a regime of subduction erosion. In the last decades a particular effort has been made to understand subduction-related processes on the basis of geophysical and geochemical studies. In southern Costa Rica and western Panama accretionary complexes and structures at the base of the volcanic front have been exposed in response to subduction of aseismic ridges and transforms. Onland exposures are located as close as to 15 km from the trench and provide a unique opportunity to better understand some of the processes occurring along the subduction zone. We provide new constraints on the origins of these exposures by integrating a comprehensive field work, new geochemical, sedimentary and paleontological data, as well as structural observations based on remote imaging. A new Campanian to Eocene tectonostratigraphy is defined for the forearc area located between the Osa Peninsula (Costa Rica) and the Azuero Peninsula (Panama). Our results show that the outer margin is composed of a complicated arrangement of igneous complexes and overlapping sedimentary sequences that essentially comprise an arc basement, primitive island-arc rocks, accreted seamount fragments and accretionary mélanges. Evidences are provided for the development of the southern Central American arc on the top an oceanic plateau. The subduction initiation along the SW edge of the Caribbean Plate occurred in the Campanian and led to formation of primitive island-arc rocks characterized by unusual geochemical affinities broadly intermediate between plateau and arc affinities. The arc was mature in the Maastrichtian and was forming a predominantly continuous landbridge between the North and South Americas. This allowed migration of terrestrial fauna between the Americas and may have contributed to the Cretaceous-Tertiary crisis by limiting trans-equatorial oceanic currents between the Pacific and the Atlantic. Several units composed of accreted seamount fragments are defined. The nature of the units and their structural arrangement provide new constraints on the modes of accretion of seamounts/oceanic islands and on the evolution of the margin since subduction initiation. Between the late Cretaceous and the middle Eocene, the margin recorded several local episodes of seamount accretion alternating with tectonic erosion. In the middle Eocene a regional tectonic event may have triggered strong coupling between the overriding and subducting plates, leading to higher rates of seamount accretion. During this period the situation along the margin was very similar to the present and characterized by subducting seamounts and absence of sediment accretion. The geological record shows that it is not possible to ascribe an overall erosive or accretionary nature to the margin in the past and, by analogy, today, because (1) accretionary and erosive processes exhibit significant lateral and temporal variations and (2) it is impossible to estimate the exact amount of material tectonically eroded from the margin since subduction initiation. In southern Costa Rica, accreted seamount fragments point toward a plume-ridge interaction in the Pacific in the late Cretaceous/Paleocene. This occurrence of accreted seamount fragments and morphology of the Pacific Ocean floor is indicative of the formation of the Cocos-Nazca spreading system at least ~40 Ma prior to the age proposed in current tectonic models. In Panama, we identified a remarkably-well preserved early Eocene oceanic island that accreted in the middle Eocene. The accretion probably occurred at very shallow depth by detachment of the island in the trench and led to an exceptional preservation of the volcanic structures. Exposures of both deep and superficial parts of the volcanic edifice have been studied, from the submarine-shield to subaerial-postshield stages. The stratigraphy allowed us to distinguish lavas produced during the submarine and subaerial stages. The lava compositions likely define a progressive diminution of source melting and a decrease in the temperature of erupted melts in the latest stages of volcanic activity. We interpret these changes to primarily reflect the progressive migration of the oceanic island out of the melting region or hotspot.

Paleomagnetic data and K-Ar ages from Mesozoic units of the Santa Marta massif: A preliminary interpretation for block rotation and translations

We report 6 K-Ar ages and paleomagnetic data from 28 sites collected in Jurassic, Lower Cretaceous and Paleocene rocks of the Santa Marta massif, to test previous hypothesis of rotations and translations of this massif, whose rock assemblage differs from other basement-cored ranges adjacent to the Guyana margin. Three magnetic components were identified in this study. A first component has a direction parallel to the present magnetic field and was uncovered in all units (D 352, I = 25.6, k = 57.35, a95 = 5.3, N = 12). A second component was isolated in Cretaceous limestone and Jurassic volcaniclastic rocks (D = 8.8, I = 8.3, k = 24.71, a95 = 13.7, N = 6), and it was interpreted as of Early Cretaceous age. In Jurassic sites with this component, Early Cretaceous K-Ar ages obtained from this and previous studies are interpreted as reset ages. The third component was uncovered in eight sites of Jurassic volcaniclastic rocks, and its direction indicates negative shallow to moderate inclinations and northeastward declinations. K-Ar ages in these sites are of Early (196.5 +/- 4.9 Ma) to early Late Jurassic age (156.6 +/- 8.9 Ma). Due to local structural complexity and too few Cretaceous outcrops to perform a reliable unconformity test, we only used two sites with (1) K-Ar ages, (2) less structural complexity, and (3) reliable structural data for Jurassic and Cretaceous rocks. The mean direction of the Jurassic component is (D = 20.4, I = -18.2, k = 46.9, a95 = 5.1, n = 18 specimens from two sites). These paleomagnetic data support previous models of northward along-margin translations of Grenvillian-cored massifs. Additionally, clockwise vertical-axis rotation of this massif, with respect to the stable craton, is also documented; the sense of rotation is similar to that proposed for the Perija Range and other ranges of the southern Caribbean margin. More data is needed to confirm the magnitudes of rotations and translations. (C) 2009 Elsevier Ltd. All rights reserved.

Mesoproterozoic rapakivi granites of the Rondonia Tin Province, southwestern border of the Amazonian craton, Brazil - I. Reconnaissance U-Pb geochronology and regional implications

Rapakivi granites and associated mafic and ultramafic rocks in the Rondonia Tin Province, southwestern Amazonian craton, Brazil were emplaced during six discrete episodes of magmatism between ca 1600 and 970 Ma. The seven rapakivi granite suites emplaced at this time were the Serra da Providencia Intrusive Suite (U-Pb ages between 1606 and 1532 Ma); Santo Antonio Intrusive Suite(U-Pb age 1406 Ma), Teotonio Intrusive Suite (U-Pb age 1387 Ma); Alto Candeias Intrusive Suite (U-Pb ages between 1346 and 1338 Ma); Sao Lourenco-Caripunas Intrusive Suite (U-Pb ages between 1314 and 1309 Ma); Santa Clara Intrusive Suite (U-Pb ages between 1082 and 1074 Ma); and Younger Granites of Rondonia (U-Pb ages between 998 and 974 Ma). The Serra da Providencia Intrusive Suite intruded the Paleoproterozoic (1.80 to 1.70 Ga) Rio Negro-Juruena crust whereas the other suites were emplaced into the 1.50 to 1.30 Ga Rondonia-San Ignacio crust. Their intrusion was contemporaneous with orogenic activity in other parts of the southwestern Amazonian craton, except for the oldest, Serra da Providencia Intrusive Suite. Orogenic events coeval with emplacement of the Serra da Providencia Intrusive Suite are not clearly recognized in the region. The Santo Antonio, Teotonio, Alto Candeias and Sao Lourenco-Caripunas Intrusive Suites are interpreted to represent extensional anorogenic magmatism associated with the terminal stages of the Rondonian-San Ignacio orogeny. At least the Sao Lourenco-Caripunas rapakivi granites and coeval intra-continental rift sedimentary rocks may, in contrast, represent the products of extensional tectonics and rifting preceding the Sunsas/Aguapei orogeny (1.25 to 1.0 Ga). The two youngest rapakivi suites, the Santa Clara Intrusive Suite and Younger Granites of Rondonia, seemingly represent inboard magmatism in the Rondonian-San Ignacio Province during a younger episode of reworking in the Rio Negro-Juruena Province during the waning stages of the collisional 1.1 to 1.0 Ga Sunsas/Aguapei orogeny. The six intra-plate rapakivi granite episodes in the southwestern part of the Amazonian craton form three broad periods of anorogenic magmatism that have age-correlative events composed of similar rocks and geologic environments in eastern Laurentia and Baltica, although the exact timing of magmatism appears slightly different. Recognition of lithologic and chronological correlations between various cratons provide important constraints to models explaining the interplay between rapakivi granite magmatism and deep crustal evolution of an early Mesoproterozoic supercontinent. They are, furthermore, important to plate tectonic models for the assembly, dispersal and reassembly of Amazonia, Laurentia and Baltica in the Mesoproterozoic and Neoproterozoic. (C) 1999 Elsevier B.V. B.V. All rights reserved.

Levantamento gravimétrico na região de Franca - SP e delfinópolis - MG, no reconhecimento e delimitação de feições estruturais (Zonas de Sutura)

This paper presents the results of a geophysical study carried out in northeastern São Paulo State and southwestern Minas Gerais State over an area 80 km wide and 97.5 km long in SE Brazil. The Bouguer anomaly map, and geological and structural data allowed to identify three different gravity domains - crustal blocks limited by major discontinuities -related to the structural pattern of the area. These discontinuities were interpreted as geosuture zones underlying the Paraná Basin sediments which have extensions in the Pre-cambrian Basement. The crustal discontinuities named Alterosa and Ribeirão Preto may be seen as A - type collision sutures in a triple junction arramgement. Two prominent linear anomalies are recognized in the Bouguer anomaly map, as well as the limit between the Brasilia and São Paulo crustal blocks or paleoplates. The Alterosa suture zone trends NW-SE while the Ribeirão Preto suture has a NE-SW direction. The Bouguer anomaly map provides subsidies and information on new concepts and theories leading to the refinement of tectonic models.

Estudo das feições estruturais e tectônicas do nordeste do estado de São Paulo e sudoeste do estado de Minas Gerais através da gravimetria

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

The Basement of the Deccan Traps and Its Madagascar Connection: Constraints from Xenoliths

Paleogeographic reconstructions of India and Madagascar before their late Cretaceous rifting juxtapose the Antongil Block of Madagascar against the Deccan Traps of India, indicating that the Western Dharwar Craton extends below the Deccan lavas. Some recent studies have suggested that the South Maharashtra Shear Zone along the northern Konkan coast of India limits the northern extent of the Western Dharwar Craton, implying that the craton does not extend below the Deccan Traps, raising a question mark on paleogeographic reconstructions of India and Madagascar. The continuity of the Western Dharwar Craton north of the South Maharashtra Shear Zone below the Deccan Traps—or its lack thereof—is critical for validating tectonic models correlating Madagascar with India. In this study, zircons in tonalitic basement xenoliths hosted in Deccan Trap dykes were dated in situ, using the U-Pb isotope system. The data furnish U-Pb ages that define three populations at 2527 ± 6, 2456 ± 6, and 2379 ± 9 Ma. The 2527 ± 6 Ma ages correspond to the igneous crystallization of the tonalites, whereas the 2456 ± 6 and 2379 ± 9 Ma ages date metamorphic overprints. The results help to establish for the first time that the basement is a part of the Neoarchean granitoid suite of the Western Dharwar Craton, which extends northward up to at least Talvade in central and Kihim beach in the western Deccan. By implication, the South Maharashtra Shear Zone cannot be the northern limit of the Western Dharwar Craton. The granitoids are correlated with the Neoarchean felsic intrusions (2.57–2.49) of the Masaola suite in the Antongil Block of Madagascar, supporting the existence of a Neoarchean Greater Dharwar Craton comprising the Western Dharwar Craton and the Antongil-Masora Block.

40Ar/39Ar single-crystal dating of detrital muscovite and K-feldspar of ODP Holes 116-717A and 116-718C

Detrital K-feldspars and muscovites from Ocean Drilling Program Leg 116 cores that have depositional ages from 0 to 18 Ma have been dated by the 40Ar/39Ar technique. Four to thirteen individual K-feldspars have been dated from seven stratigraphic levels, each of which have a very large range, up to 1660 Ma. At each level investigated, at least one K-feldspar yielded an age minimum which is, within uncertainty, identical to the age of deposition. One to twelve single muscovite crystals from each of six levels have also been studied. The range of muscovite ages is less than that of the K-feldspars and, with one exception, reveal only a 20-Ma spread in ages. As with the K-feldspars, each level investigated contains muscovites with mineral ages essentially identical to depositional ages. These results indicate that a significant portion of the material in the Bengal Fan is first-cycle detritus derived from the Himalayas. Therefore, the significant proportion of sediment deposited in the distal fan in the early to mid Miocene can be ascribed to a significant pulse of uplift and erosion in the collision zone. Moreover, these data indicate that during the entire Neogene, some portion of the Himalayan orogen was experiencing rapid erosion (<= uplift). The lack of granulite facies rocks in the eastern Himalayas and Tibetan Plateau suggests that very rapid uplift must have been distributed in brief pulses in different places in the mountain belt. We suggest that the great majority of the crystals with young apparent ages have been derived from the southern slope of the Himalayas, predominantly from near the main central thrust zone. These data provide further evidence against tectonic models in which the Himalayas and Tibetan plateaus are uplifted either uniformly during the past 40 m.y. or mostly within the last 2 to 5 m.y.

Vertical crustal motion along the Mediterranean and Black Sea coast derived from ocean altimetry and tide gauge data

Tide gauge (TG) data along the northern Mediterranean and Black Sea coasts are compared to the sea-surface height (SSH) anomaly obtained from ocean altimetry (TOPEX/Poseidon and ERS-1/2) for a period of nine years (1993–2001). The TG measures the SSH relative to the ground whereas the altimetry does so with respect to the geocentric reference frame; therefore their difference would be in principle a vertical ground motion of the TG sites, though there are different error sources for this estimate as is discussed in the paper. In this study we estimate such vertical ground motion, for each TG site, from the slope of the SSH time series of the (non-seasonal) difference between the TG record and the altimetry measurement at a point closest to the TG. Where possible, these estimates are further compared with those derived from nearby continuous Global Positioning System (GPS) data series. These results on vertical ground motion along the Mediterranean and Black Sea coasts provide useful source data for studying, contrasting, and constraining tectonic models of the region. For example, in the eastern coast of the Adriatic Sea and in the western coast of Greece, a general subsidence is observed which may be related to the Adriatic lithosphere subducting beneath the Eurasian plate along the Dinarides fault.

Evolution of stress deficit and changing rates of seismicity in cellular automaton models of earthquake faults

We investigate the internal dynamics of two cellular automaton models with heterogeneous strength fields and differing nearest neighbour laws. One model is a crack-like automaton, transferring ail stress from a rupture zone to the surroundings. The other automaton is a partial stress drop automaton, transferring only a fraction of the stress within a rupture zone to the surroundings. To study evolution of stress, the mean spectral density. f(k(r)) of a stress deficit held is: examined prior to, and immediately following ruptures in both models. Both models display a power-law relationship between f(k(r)) and spatial wavenumber (k(r)) of the form f(k(r)) similar tok(r)(-beta). In the crack model, the evolution of stress deficit is consistent with cyclic approach to, and retreat from a critical state in which large events occur. The approach to criticality is driven by tectonic loading. Short-range stress transfer in the model does not affect the approach to criticality of broad regions in the model. The evolution of stress deficit in the partial stress drop model is consistent with small fluctuations about a mean state of high stress, behaviour indicative of a self-organised critical system. Despite statistics similar to natural earthquakes these simplified models lack a physical basis. physically motivated models of earthquakes also display dynamical complexity similar to that of a critical point system. Studies of dynamical complexity in physical models of earthquakes may lead to advancement towards a physical theory for earthquakes.

Long-range automaton models of earthquakes: Power-law accelerations, correlation evolution, and mode-switching

We introduce a conceptual model for the in-plane physics of an earthquake fault. The model employs cellular automaton techniques to simulate tectonic loading, earthquake rupture, and strain redistribution. The impact of a hypothetical crustal elastodynamic Green's function is approximated by a long-range strain redistribution law with a r(-p) dependance. We investigate the influence of the effective elastodynamic interaction range upon the dynamical behaviour of the model by conducting experiments with different values of the exponent (p). The results indicate that this model has two distinct, stable modes of behaviour. The first mode produces a characteristic earthquake distribution with moderate to large events preceeded by an interval of time in which the rate of energy release accelerates. A correlation function analysis reveals that accelerating sequences are associated with a systematic, global evolution of strain energy correlations within the system. The second stable mode produces Gutenberg-Richter statistics, with near-linear energy release and no significant global correlation evolution. A model with effectively short-range interactions preferentially displays Gutenberg-Richter behaviour. However, models with long-range interactions appear to switch between the characteristic and GR modes. As the range of elastodynamic interactions is increased, characteristic behaviour begins to dominate GR behaviour. These models demonstrate that evolution of strain energy correlations may occur within systems with a fixed elastodynamic interaction range. Supposing that similar mode-switching dynamical behaviour occurs within earthquake faults then intermediate-term forecasting of large earthquakes may be feasible for some earthquakes but not for others, in alignment with certain empirical seismological observations. Further numerical investigation of dynamical models of this type may lead to advances in earthquake forecasting research and theoretical seismology.

On the applicability of joint inversion of gravity and resistivity data to the study of a tectonic sedimentary basin in Northern Portugal

The Chaves basin is a pull-apart tectonic depression implanted on granites, schists, and graywackes, and filled with a sedimentary sequence of variable thickness. It is a rather complex structure, as it includes an intricate network of faults and hydrogeological systems. The topography of the basement of the Chaves basin still remains unclear, as no drill hole has ever intersected the bottom of the sediments, and resistivity surveys suffer from severe equivalence issues resulting from the geological setting. In this work, a joint inversion approach of 1D resistivity and gravity data designed for layered environments is used to combine the consistent spatial distribution of the gravity data with the depth sensitivity of the resistivity data. A comparison between the results from the inversion of each data set individually and the results from the joint inversion show that although the joint inversion has more difficulty adjusting to the observed data, it provides more realistic and geologically meaningful models than the ones calculated by the inversion of each data set individually. This work provides a contribution for a better understanding of the Chaves basin, while using the opportunity to study further both the advantages and difficulties comprising the application of the method of joint inversion of gravity and resistivity data.

Depleted arc volcanism in the Alboran Sea and shoshonitic volcanism in Morocco: geochemical and isotopic constraints on Neogene tectonic processes

Samples of volcanic rocks from Alboran Island, the Alboran Sea floor and from the Gourougou volcanic centre in northern Morocco have been analyzed for major and trace elements and Sr-Nd isotopes to test current theories on the tectonic geodynamic evolution of the Alboran Sea. The Alboran Island samples are low-K tholeiitic basaltic andesites whose depleted contents of HFS elements (similar to0.5xN-MORB), especially Nb (similar to0.2xN-MORB), show marked geochemical parallels with volcanics from immature intra-oceanic arcs and back-arc basins. Several of the submarine samples have similar compositions, one showing low-Ca boninite affinity. Nd-143/Nd-144 ratios fall in the same range as many island-arc and back-arc basin samples, whereas Sr-87/Sr-86 ratios (on leached samples) are somewhat more radiogenic. Our data point to active subduction taking place beneath the Alboran region in Miocene times, and imply the presence of an associated back-arc spreading centre. Our sea floor suite includes a few more evolved dacite and rhyolite samples with (Sr-87/Sr-86)(0) up to 0.717 that probably represent varying degrees of crustal melting. The shoshonite and high-K basaltic andesite lavas from Gourougou have comparable normalized incompatible-element enrichment diagrams and Ce/Y ratios to shoshonitic volcanics from oceanic island arcs, though they have less pronounced Nb deficits. They are much less LIL- and LREE-enriched than continental arc analogues and post-collisional shoshonites from Tibet. The magmas probably originated by melting in subcontinental lithospheric mantle that had experienced negligible subduction input. Sr-Nd isotope compositions point to significant crustal contamination which appears to account for the small Nb anomalies. The unmistakable supra-subduction zone (SSZ) signature shown by our Alboran basalts and basaltic andesite samples refutes geodynamic models that attribute all Neogene volcanism in the Alboran domain to decompression melting of upwelling asthenosphere arising from convective thinning of over-thickened lithosphere. Our data support recent models in which subsidence is caused by westward rollback of an eastward-dipping subduction zone beneath the westemmost Mediterranean. Moreover, severance of the lithosphere at the edges of the rolling-back slab provides opportunities for locally melting lithospheric mantle, providing a possible explanation for the shoshonitic volcanism seen in northern Morocco and more sporadically in SE Spain. (C) 2004 Elsevier B.V. All rights reserved.

Quantifying rates of landscape evolution and tectonic processes by thermochronology and numerical modeling of crustal heat transport using PECUBE

PECUBE is a three-dimensional thermal-kinematic code capable of solving the heat production-diffusion-advection equation under a temporally varying surface boundary condition. It was initially developed to assess the effects of time-varying surface topography (relief) on low-temperature thermochronological datasets. Thermochronometric ages are predicted by tracking the time-temperature histories of rock-particles ending up at the surface and by combining these with various age-prediction models. In the decade since its inception, the PECUBE code has been under continuous development as its use became wider and addressed different tectonic-geomorphic problems. This paper describes several major recent improvements in the code, including its integration with an inverse-modeling package based on the Neighborhood Algorithm, the incorporation of fault-controlled kinematics, several different ways to address topographic and drainage change through time, the ability to predict subsurface (tunnel or borehole) data, prediction of detrital thermochronology data and a method to compare these with observations, and the coupling with landscape-evolution (or surface-process) models. Each new development is described together with one or several applications, so that the reader and potential user can clearly assess and make use of the capabilities of PECUBE. We end with describing some developments that are currently underway or should take place in the foreseeable future. (C) 2012 Elsevier B.V. All rights reserved.