Middle Miocene peralkaline ignimbrites in the Hermosillo region (Sonora, Mexico): Geodynamic implications

Scattered ignimbritic mesas crop out in the Hermosillo region (Sonora, Mexico). These rocks that have been dated at 12.5 Ma (Middle Miocene) have the petrography and chemical characteristics of comendites. Such a flare-up of peralkaline acidic volcanism, after a long period of subduction-related arc volcanism, emphasises an important change in the source of volcanism. It corresponds to the latest stage of continental extension prior to the marine invasion and the development of spreading centres in the Gulf of California. To cite this article: J. Vidal Solano et al., C. R. Geoscience 337 (2005).     

The Aguablanca Ni–(Cu) sulfide deposit, SW Spain: geologic and geochemical controls and the relationship with a midcrustal layered mafic complex

The Aguablanca Ni–(Cu) sulfide deposit is hosted by a breccia pipe within a gabbro–diorite pluton. The deposit probably formed due to the disruption of a partially crystallized layered mafic complex at about 12–19 km depth and the subsequent emplacement of melts and breccias at shallow levels (<2 km). The ore-hosting breccias are interpreted as fragments of an ultramafic cumulate, which were transported to the near surface along with a molten sulfide melt. Phlogopite Ar–Ar ages are 341–332 Ma in the breccia pipe, and 338–334 Ma in the layered mafic complex, and are similar to recently reported U–Pb ages of the host Aguablanca Stock and other nearby calc-alkaline metaluminous intrusions (ca. 350–330 Ma). Ore deposition resulted from the combination of two critical factors, the emplacement of a layered mafic complex deep in the continental crust and the development of small dilational structures along transcrustal strike-slip faults that triggered the forceful intrusion of magmas to shallow levels. The emplacement of basaltic magmas in the lower middle crust was accompanied by major interaction with the host rocks, immiscibility of a sulfide melt, and the formation of a magma chamber with ultramafic cumulates and sulfide melt at the bottom and a vertically zoned mafic to intermediate magmas above. Dismembered bodies of mafic/ultramafic rocks thought to be parts of the complex crop out about 50 km southwest of the deposit in a tectonically uplifted block (Cortegana Igneous Complex, Aracena Massif). Reactivation of Variscan structures that merged at the depth of the mafic complex led to sequential extraction of melts, cumulates, and sulfide magma. Lithogeochemistry and Sr and Nd isotope data of the Aguablanca Stock reflect the mixing from two distinct reservoirs, i.e., an evolved siliciclastic middle-upper continental crust and a primitive tholeiitic melt. Crustal contamination in the deep magma chamber was so intense that orthopyroxene replaced olivine as the main mineral phase controlling the early fractional crystallization of the melt. Geochemical evidence includes enrichment in SiO₂ and incompatible elements, and Sr and Nd isotope compositions (⁸⁷Sr/⁸⁶Sr_i 0.708–0.710; ¹⁴³Nd/¹⁴⁴Nd_i 0.512–0.513). However, rocks of the Cortegana Igneous Complex have low initial ⁸⁷Sr/⁸⁶Sr and high initial ¹⁴³Nd/¹⁴⁴Nd values suggesting contamination by lower crustal rocks. Comparison of the geochemical and geological features of igneous rocks in the Aguablanca deposit and the Cortegana Igneous Complex indicates that, although probably part of the same magmatic system, they are rather different and the rocks of the Cortegana Igneous Complex were not the direct source of the Aguablanca deposit. Crust–magma interaction was a complex process, and the generation of orebodies was controlled by local but highly variable factors. The model for the formation of the Aguablanca deposit presented in this study implies that dense sulfide melts can effectively travel long distances through the continental crust and that dilational zones within compressional belts can effectively focus such melt transport into shallow environments.Electronic supplementary material Supplementary material is available in the online version of this article at and is accessible for authorized users.     

Visualization of Volcanic Rock Geochemical Data and Classification with Artificial Neural Networks

An unsupervised neural network technique, Growing Cell Structures (GCS) was used to visualize geochemical differences between four different island arc volcanic rock types: basalts, andesites, dacites and rhyolites. The output of the method shows the cluster structure of the dataset clearly, and the relevant geochemical patterns and relationships between its variables. The data can be separated into four clusters, each associated with a specific volcanic rock type (basalt, andesite, dacite and rhyolite), according to a unique combination of major element concentrations. Following clustering, performance of the trained GCS network as a classifier of volcanic rock type was tested using two test datasets with major element concentration data for 312 and 496 island arc volcanic rock samples of known volcanic type. Preliminary classification results are promising. In the first test dataset 94% of basalts, 76% of andesites, 83% of dacites and 100% of the rhyolites were classified correctly. Successful classification rates in the second dataset were 100%, 80%, 77%, and 98% respectively. The success of the analysis suggests that neural networks analysis constitutes a useful analytical tool for identification of natural clusters and examination of the relationships between numeric variables in large datasets, and that can be used for automatic classification of new data.     

Stratigraphy of the Middle Devonian boundary: Formal definition of the susceptibility magnetostratotype in Germany with comparisons to sections in the Czech Republic, Morocco and Spain

One major difficulty in geology is high-resolution correlation among widely separated sections, especially in the Paleozoic where magnetostratigraphy polarity is not well established because rocks are often remagnetized, where critical biostratigraphic zonation may be poor or lacking, or where structural complexities make correlations very difficult. To address this problem, we have been using magnetostratigraphy susceptibility measurements. Here, we report our work from the Middle Devonian in Europe and North Africa. The Middle Devonian (Emsian–Eifelian) global boundary stratotype section and point (GSSP), located in the Eifel Hills, western Germany, was ratified by the International Subcommission on Stratigraphy in 1985, after careful evaluation of the biostratigraphy for this and many other sections. The boundary interval has been characterized using biostratigraphy, and the beginning of the Eifelian stage has been specifically defined by the first occurrence of the conodont Polygnathus costatus partitus. We have collected the Eifel Hills section for magnetic susceptibility (MS) measurement and here we establish it as the magnetostratotype for the Emsian–Eifelian stage boundary, by formally defining the magnetostratigraphy susceptibility for the section. We then collected, measured and compared the magnetostratotype to four other sections for which conodont biostratigraphy has been studied and where P. costatus partitus is present; two Emsian–Eifelian sections in Morocco and two sections in the Czech Republic (including the Emsian–Eifelian parastratotype). Finally, we have measured the MS for the El Puerto Creek section in the Cantabrian Mountains of Spain and identified the location of the Emsian–Eifelian boundary within the section based on MS comparison to the GSSP in conjunction with excellent biostratigraphic indicators, primarily brachiopods. While the conodont zonation in the El Puerto Creek section is poorly defined, we believe that the correspondence between the MS and biostratigraphy in the section allows the identification of the Emsian–Eifelian boundary. These results indicate that this method can be successfully applied to marine sequences where ambiguities in correlation exist.     

Elastic and inelastic triggering of earthquakes in the North Anatolian Fault zone

Deformation models used to explain the triggering mechanism often assume pure elastic behaviour for the crust and upper mantle. In reality however, the mantle and possibly the lower crust behave viscoelastically, particularly over longer time scales. Consequently, the stress field of an earthquake is in general time-dependent. In addition, if the elastic stress increase were enough to trigger a later earthquake, this triggered event should occur instantaneously and not many years after the triggering event. Hence, it is appropriate to include inelastic behaviour when analysing stress transfer and earthquake interaction.In this work, we analyse a sequence of 10 magnitude M_s > 6.5 events along the North Anatolian Fault between 1939 and 1999 to study the evolution of the regional Coulomb stress field. We investigate the triggering of these events by stress transfer, taking viscoelastic relaxation into account. We evaluate the contribution of elastic stress changes, of post-seismic viscoelastic relaxation in the lower crust and mantle, and of steady tectonic loading to the total Coulomb stress field. We analyse the evolution of stress in the region under study, as well as on the rupture surfaces of the considered events and their epicentres. We study the state of the Coulomb stress field before the 1999 İzmit and Düzce earthquakes, as well as in the Marmara Sea region.In general, the Coulomb stress failure criterion offers a plausible explanation for the location of these events. However, we show that using a purely elastic model disregards an important part of the actual stress increase/decrease. In several cases, post-seismic relaxation effects are important and greater in magnitude than the stress changes due to steady tectonic loading. Consequently, viscoelastic relaxation should be considered in any study dealing with Coulomb stress changes.According to our study, and assuming that an important part of the rupture surface must be stressed for an earthquake to occur, the most likely value for the viscosity of the lower crust or mantle in this region is 5 · 10¹⁷–10¹⁸ Pa · s. Our results cannot rule out the possibility of other time-dependent processes involved in the triggering of the 1999 Düzce event. However, the stress increase due to viscoelastic relaxation brought 22% of the 1999 Düzce rupture area over the threshold value of Δσ_c ≥ 0.01 MPa (0.1 bar), and took the whole surface closer to failure by an average of 0.2 MPa. Finally, we argue that the Marmara Sea region is currently being loaded with positive Coulomb stresses at a much faster rate than would arise exclusively from steady tectonic loading on the North Anatolian Fault.     

Geophysical Study of the Valle Fértil Lineament Between 28°45′S and 31°30′S: Boundary Between the Cuyania and Pampia Terranes

A gravimetric and magnetometric study was carried out in the north-eastern portion of the Cuyania terrane and adjacent Pampia terrane. Gravimetric models permitted to interpret the occurrence of dense materials at the suture zone between the latter terranes. Magnetometric models led to propose the existence of different susceptibilities on either side of the suture. The Curie temperature point depth, representing the lower boundary of the magnetised crust, was found to be located at 25 km, consistent with the lower limit of the brittle crust delineated by seismic data; this unusually thick portion of the crust is thought to release stress producing significant seismicity.

Moho depths determined from seismic studies near western Sierras Pampeanas are significantly greater than those obtained from gravimetric crustal models.

Considering mass and gravity changes originated by the flat-slab Nazca plate along Cuyania and western Pampia terranes, it is possible to reconcile Moho thickness obtained either by seismic or by gravity data. Thus, topography and crustal thickness are controlled not only by erosion and shortening but by upper mantle heterogeneities produced by: (a) the oceanic subducted Nazca plate with “normal slope” also including asthenospheric materials between both continental and oceanic lithospheres; (b) flat-slab subducted Nazca plate (as shown in this work) without significant asthenospheric materials between both lithospheres. These changes influence the relationship between topographic altitudes and crustal thickness in different ways, differing from the simple Airy system relationship and modifying the crustal scale shortening calculation. These changes are significantly enlarged in the study area. Future changes in Nazca Plate slope will produce changes in the isostatic balance.     

Optimum seismic zoning for multiple types of structures

We consider zoning for the design criterion that minimizes the expected present value of the total cost, including the initial cost as well as losses due to damage and failure. The problem consists of the following: given the number of zones, their boundaries and design coefficients must be such that they minimize the expected present value of all structures built in the region. We will refer to solutions in one or more dimensions, depending on the number of the types of structures built in the region to be zoned. Two methods are proposed to solve the problems. The first method is based on the different combinations performed in order to attain optimum zoning. The second method uses an analogy to the evolution of biological systems. The work ends by applying the methods developed to a region of known seismicity. Copyright © 2003 John Wiley & Sons, Ltd.     

Lead sources in Mesozoic and Cenozoic Andean ore deposits,north-central Chile (30–34°S)

Mesozoic and Cenozoic ore deposits in the Chilean Andes between La Serena (~30°S) and Santiago (~34°S) include polymetallic vein, low- and high-sulfidation epithermal vein, skarn, porphyry copper-molybdenum and porphyry copper-gold. These deposits are associated with volcanic and plutonic complexes emplaced in eastward-migrating longitudinal arcs which formed during subduction along the continental margin of South America since the Middle Jurassic. Stratabound, but epigenetic, volcanic rock- and sedimentary rock-hosted manto deposits contain additional copper resources. Lead isotopic compositions in ore minerals from 29 deposits vary with age and geographic location, and hence with basement and host rocks. Lead in most ore deposits is derived from temporally related igneous rocks, except for the manto deposits whose lead is derived from host volcanic and sedimentary rock sequences. Lead in the ore deposits is dominated by two crustal sources. Low ²⁰⁷Pb/²⁰⁴Pb characterizes one source whereas high ²⁰⁷Pb/²⁰⁴Pb characterizes the second source. Lead isotopic compositions of Jurassic and Miocene ore minerals (²⁰⁶Pb/²⁰⁴Pb>18.50; ²⁰⁷Pb/²⁰⁴Pb>15.61) lie along the average crustal growth curve. By contrast, most Cretaceous deposits have ore minerals with lower ²⁰⁶Pb/²⁰⁴Pb (<18.39) and ²⁰⁷Pb/²⁰⁴Pb (<15.58) than Jurassic ore minerals. The shift in lead isotopic composition to lower lead isotopic values precludes derivation of lead from a source of similar composition to those in the Jurassic or Tertiary deposits. For Cretaceous deposits, polymetallic and low-sulfidation epithermal veins and a skarn have lower ²⁰⁶Pb/²⁰⁴Pb than a porphyry copper-gold system and peripheral gold veins at Andacollo (18.43-18.50). Late Cretaceous veins from the Bellavista deposit have the lowest ²⁰⁶Pb/²⁰⁴Pb (18.33) of all deposits. Ore minerals in Miocene and Pliocene porphyry copper-molybdenum deposits have higher ²⁰⁶Pb/²⁰⁴Pb (18.58-18.67) than Cretaceous deposits, consistent with their age being younger. The Miocene and Pliocene ore minerals also have higher ²⁰⁷Pb/²⁰⁴Pb (15.58-15.66) than Cretaceous ore minerals, thereby requiring an additional input from the high-²⁰⁷Pb/²⁰⁴Pb source into the younger deposits. Miocene auriferous deposits in the north have similar ²⁰⁶Pb/²⁰⁴Pb values as the Miocene and Pliocene porphyry copper-molybdenum deposits in the south, but they are distinguished by higher and variable ²⁰⁷Pb/²⁰⁴Pb (15.61-15.66) and ²⁰⁸Pb/²⁰⁴Pb (38.54-39.01), which are arrayed along steep mixing trends. These ore minerals have the largest input of high-²⁰⁷Pb/²⁰⁴Pb material in the deposits studied. By contrast, lead in the epigenetic manto deposits appears to be derived from the host volcanic or sedimentary rock-dominated sequences, and locally exhibits large-scale isotopic heterogeneity within a deposit. Overall, the lead isotopic compositions of ore minerals mimic the values and variations established in age-equivalent rock sequences. The low-²⁰⁷Pb/²⁰⁴Pb material in the deposits is derived from Cretaceous igneous rocks or their sources as they evolved with time; low ²⁰⁷Pb/²⁰⁴Pb characterizes these rocks. By contrast, high-²⁰⁷Pb/²⁰⁴Pb material is likely derived from Carboniferous to Triassic igneous rocks or their sources, as this lead isotopic characteristic dominates these rocks.     

Large Scale Rockfall Reach Susceptibility Maps in La Cabrera Sierra (Madrid) performed with GIS and Dynamic Analysis at 1:5,000

For the testing of the effect on the tsunami prevention facilities, a simplified methodfor tsunami risk assessment was suggested without wave run-up analysis. This methodis proposed using calculated offshore tsunami waveform and field reconnaissance suchas the seawall height, time necessary for residents' evacuation and tsunami warninginsurance. Then, two normalized values are evaluated; one is the ratio of calculatedmaximum tsunami height to seawall height, the other is the ratio of time betweentsunami over-topping and evacuation completion to total time required for evacuation.These two values are used to qualitatively estimate the safety of residents and the effectof tsunami prevention facilities, eliminating the necessity to compute complicatedtsunami run-up onshore.