An anisotropic viscous representation of Mohr-Coulomb failure for use is modelling coupled mantle-continent dynamics

In mantle convection models it has become common to make use of a modified (pressure sensitive, Boussinesq) von Mises yield criterion to limit the maximum stress the lithosphere can support. This approach allows the viscous, cool thermal boundary layer to deform in a relatively plate-like mode even in a fully Eulerian representation. In large-scale models with embedded continental crust where the mobile boundary layer represents the oceanic lithosphere, the von Mises yield criterion for the oceans ensures that the continents experience a realistic broad-scale stress regime. In detailed models of crustal deformation it is, however, more appropriate to choose a Mohr-Coulomb yield criterion based upon the idea that frictional slip occurs on whichever one of many randomly oriented planes happens to be favorably oriented with respect to the stress field. As coupled crust/mantle models become more sophisticated it is important to be able to use whichever failure model is appropriate to a given part of the system. We have therefore developed a way to represent Mohr-Coulomb failure within a code which is suited to mantle convection problems coupled to large-scale crustal deformation. Our approach uses an orthotropic viscous rheology (a different viscosity for pure shear to that for simple shear) to define a prefered plane for slip to occur given the local stress field. The simple-shear viscosity and the deformation can then be iterated to ensure that the yield criterion is always satisfied. We again assume the Boussinesq approximation - neglecting any effect of dilatancy on the stress field. An additional criterion is required to ensure that deformation occurs along the plane aligned with maximum shear strain-rate rather than the perpendicular plane which is formally equivalent in any symmetric formulation. It is also important to allow strain-weakening of the material. The material should remember both the accumulated failure history and the direction of failure. We have included this capacity in a Lagrangian-Integration-point finite element code and will show a number of examples of extension and compression of a crustal block with a Mohr-Coulomb failure criterion, and comparisons between mantle convection models using the von Mises versus the Mohr-Coulomb yield criteria. The formulation itself is general and applies to 2D and 3D problems, although it is somewhat more complicated to identify the slip plane in 3D.

Mesozoic to recent evolution of the Andean forearc of northern Chile (22-24 s)

The Andean forearc of northern Chile comprises four morphotectonic units, which include from east to west: 1) The Cordillera de la Costa: composed of Jurassic granites and andesites, thought to represent a volcanic arc, the Mejillones terrane, an accreted allochthonous terrane, and the Lower Cretaceous Coloso basin, which formed through forearc extension along the suture between the Mejillones terrane and the Jurassic arc. Palaeomagnetic studies of the above units have identified approximately 29+/-11 degrees of clockwise rotation. Rotation is due to extension (caused by subduction roll back and slab pull), at an angle to the direction of absolute motion of the South American Plate. 2) The Central Depression: a large arid basin containing isolated fault-bounded blocks of pre-Mesozoic metamorphosed igneous rocks, Triassic sediments and volcanics, and Jurassic carbonates, deposited in a. back-arc basin setting. The isolated blocks formed through extension along previous thrust faults, these originated through compression of the back-arc basin due to accretion of the Jurassic volcanic arc. 3) The Precordillera.: composed of Permian-Triassic rift-related sediments and volcanics, Jurassic continental sediments synchronous with back-arc basin sedimentation, and Cretaceous and Oligo-Miocene continental sediments deposited in foreland basins. Palaeomagnetism has identified clockwise rotation in rocks ranging in age from Jurassic-Miocene. Rotation in the Precordillera. affected larger structural blocks than in the Cordillera de la Costa. 4) The Salar Depression: a. series of arid continental basins developed on continental crust. These basins nay have originated in the Triassic, when rifting of the South American craton is thought to have taken place. In conclusion, palaeomagnetic and geological evidence is consistent with the view that the north Chilean forearc was largely under an extensional stress regime. However, the presence of extensive compressional structures in Palaeocene and older rocks in the forearc together with the currently active foreland thrust belt of Argentina. indicate that throughout the evolution of the Andean Orogen, a delicate balance between compressional and extensional tectonic regimes has existed.

Annotated record of the detailed examination of Mn deposits from R/V "Le Suroit" GEOMANCHE75/1 stations

We report here preliminary results from the detailed sampling of granitic basement at a depth of 4.000 m at the lower edge of the Armorican continental margin near 48°N and 12°W. We consider that we have sampled here the extreme lateral limit of the continental crust where local exposures emerge from under a cover of Cainozoic,Mesozoic and probable Palaeozoic sediments.

Chemical composition of Late Cenozoic volcanic and volcaniclastic deposits from the western Woodlark Basin area, SW Pacific

Tephra fallout layers and volcaniclastic deposits, derived from volcanic sources around and on the Papuan Peninsula, form a substantial part of the Woodlark Basin marine sedimentary succession. Sampling by the Ocean Drilling Program Leg 180 in the western Woodlark Basin provides the opportunity to document the distribution of the volcanically-derived components as well as to evaluate their chronology, chemistry, and isotope compositions in order to gain information on the volcanic sources and original magmatic systems. Glass shards selected from 57 volcanogenic layers within the sampled Pliocene-Pleistocene sedimentary sequence show predominantly rhyolitic compositions, with subordinate basaltic andesites, basaltic trachy-andesites, andesites, trachy-andesites, dacites, and phonolites. It was possible to correlate only a few of the volcanogenic layers between sites using geochemical and age information apparently because of the formation of strongly compartmentalised sedimentary realms on this actively rifting margin. In many cases it was possible to correlate Leg 180 volcanic components with their eruption source areas based on chemical and isotope compositions. Likely sources for a considerable number of the volcanogenic deposits are Moresby and Dawson Strait volcanoes (D'Entrecasteaux Islands region) for high-K calc-alkaline glasses. The Dawson Strait volcanoes appear to represent the source for five peralkaline tephra layers. One basaltic andesitic volcaniclastic layer shows affinities to basaltic andesites from the Woodlark spreading tip and Cheshire Seamount. For other layers, a clear identification of the sources proved impossible, although their isotope and chemical signatures suggest similarities to south-west Pacific subduction volcanism, e.g. New Britain and Tonga- Kermadec island arcs. Volcanic islands in the Trobriand Arc (for example, Woodlark Island Amphlett Islands and/or Egum Atoll) are probable sources for several volcaniclastic layers with ages between 1.5 to 3 Ma. The Lusancay Islands can be excluded as a source for the volcanogenic layers found during Leg 180. Generally, the volcanogenic layers indicate much calc-alkaline rhyolitic volcanism in eastern Papua since 3.8 Ma. Starting at 135 ka, however, peralkaline tephra layers appear. This geochemical change in source characteristics might reflect the onset of a change in geotectonic regime, from crustal subduction to spreading, affecting the D'Entrecasteaux Islands region. Initial 143Nd/144Nd ratios as low as 0.5121 and 0.5127 for two of the tephra layers are interpreted as indicating that D'Entrecasteaux Islands volcanism younger than 2.9 Ma occasionally interacted with the Late Archean basement, possibly reflecting the mobilisation of the deep continental crust during active rift propagation.

(Table 1) Sample locations, Hf and Nd isotopic compositions and concentrations of Atlantic seawater samples with hydrographic details

The first full water column hafnium isotopic compositions of Atlantic seawater have been obtained at seven locations from the Labrador Sea to the Drake Passage. Despite subpicomolar concentrations in seawater, a precision of the Hf isotopic measurements of <0.7 epsilon-Hf units was achieved. An overall epsilon-Hf range between -3.1 in the Labrador Sea and +4.4 in Antarctic bottom water was determined, the distribution of which broadly reflects continental weathering inputs. Within particular water column profiles, significant differences of up to 4 epsilon-Hf units occur. Combined with Nd isotope data of the same samples, it is evident that the Hf isotopic composition of seawater is too radiogenic for a given Nd isotopic composition and that the largest difference between expected and measured Hf isotopic compositions in seawater occurs near the oldest continental crust in the Labrador Sea. This corroborates the previous proposition, which was mainly based on ferromanganese crust data, that the Hf isotopic composition of seawater is controlled by incongruent weathering of continental crust and possibly, to some extent, by hydrothermal contributions. Hafnium concentrations in the ocean do not increase along the deep ocean conveyer indicating an oceanic residence time of only a few hundred years, which is significantly shorter than previously assumed. The Hf isotopic composition of past seawater can therefore serve as a proxy for short distance, basin scale mixing processes and the regime and intensity of nearby continental weathering processes.

The geochemistry, mineralogy, and petrology of basalts at DSDP Leg 59 Holes

Legs 59 and 60 of the International Phase of Oceanic Drilling (IPOD) were designed to study the nature and history of volcanism of the active Mariana arc, its currently spreading inter-arc basin (the Mariana Trough), and the series of inactive basins and intervening ridges that lie to the west. The older basins and ridges were drilled during Leg 59 as the first part of a transect of single-bit holes drilled in each major basin and ridge. The eastern part of the transect - the technically active region - was drilled during Leg 60. The evolution of island-arc volcanos and magma genesis associated with lithospheric subduction remain some of the most complex petrologic problems confronting us. Many types of source material (mantle, oceanic crust, continental crust) and an unusually wide range of possible physical conditions at the time of magma genesis must be identified even before the roles of partial melting and subsequent magma fractionation, mixing, and contamination can be assessed.

Isotope analysis and heat flow measurements of Maria S. Merian cruise MSM21/4 on the western Svalbard margin

Methane hydrate is an ice-like substance that is stable at high-pressure and low temperature in continental margin sediments. Since the discovery of a large number of gas flares at the landward termination of the gas hydrate stability zone off Svalbard, there has been concern that warming bottom waters have started to dissociate large amounts of gas hydrate and that the resulting methane release may possibly accelerate global warming. Here, we can corroborate that hydrates play a role in the observed seepage of gas, but we present evidence that seepage off Svalbard has been ongoing for at least three thousand years and that seasonal fluctuations of 1-2°C in the bottom-water temperature cause periodic gas hydrate formation and dissociation, which focus seepage at the observed sites.

Geochemical composition of minerals obtained from ODP Sites 180-1109, 180-1117 and 180-1118

During the last 8 m.y. the Papuan Peninsula region of Papua New Guinea has been affected by extension which opened the Woodlark Basin. The present-day spreading tip is located at the foot of the Moresby Seamount, a crustal block whose northern flank is an active low-angle normal fault related to this extension. During Ocean Drilling Program Leg 180 (7 June-11 August 1998), 11 sites (1108-1118) were drilled along a north-south-trending transect across the Woodlark Basin just ahead of the spreading tip. Four of these sites (1118, 1109, 1114, and 1117) reached the crystalline basement, which is composed of diabase and gabbro. Sites 1118 and 1109, located on the Woodlark Rise, belong to the hanging wall block, and Sites 1114 and 1117, located on the crest of the Moresby Seamount, belong to the footwall block and the fault zone itself. Most of the basalt, diabase, and gabbro that were recovered show a well-preserved magmatic texture. The diabase, which is the most abundant rock type, has a coarse-grained ophitic texture composed of poikilitic clinopyroxene including radiating, locally skeletal plagioclase laths with interstitial iron oxide grains. Secondary mineralogy consists of chlorite, zeolite, calcite, albite, and quartz. The gabbro shows a medium-grained granular texture. The magmatic mineralogy consists of euhedral laths of plagioclase and anhedral interstitial clinopyroxene. Secondary mineralogy consists of a magnesio to actinolitic hornblende, chlorite, clinozoisite, zeolite, quartz, and calcite. The retrograde metamorphic evolution of both gabbro and diabase occurred under low amphibolite to subgreenschist facies conditions associated mainly with brittle deformation and the development of a local low-temperature shear zone. This shows no evidence for high thermal gradient in the crust during the continental rifting.

Radionuclides of bulk and terrigenous surface sediments in the Southern Ocean

Sr and Nd isotopic compositions of Late Quaternary surface sediment and sediment cores from the south Atlantic and southeast Pacific sectors of the Southern Ocean are used to constrain the provenance and transport mechanisms of their terrigenous component. We report isotopic and mineralogical data for core samples from three localities, the Mid-Atlantic Ridge at 41°S and the northern and southern Scotia Sea. In addition, data for surface sediment samples from the south Atlantic and southeast Pacific sectors of the Southern Ocean are presented. The variations of Sr and Nd isotopic compositions of the bulk sediment samples in all cores were correlated with the magnetic susceptibility of the sediment and with the inferred glacial-interglacial stages. The isotopic data indicate that, during glacial periods, sediment was delivered from continental crust with a shorter residence time than that supplying material during interglacial periods. At the core site near the Mid-Atlantic Ridge, Nd isotopic, combined with mineralogical evidence indicates interglacial period deposition of a relatively high amount of kaolinite and silt with low epsilon-Nd values < -8. The material was probably supplied by North Atlantic Deep Water from low latitudes. For glacial periods, a high contribution of silt and clay with epsilon-Nd > -4.5, probably derived from southern South America, was indicated. The glacial-interglacial shift in sources may be due to either a decreasing influence of North Atlantic Deep Water during glacial times or by a larger contribution of glaciogenic detritus from southern South America. At the core site in the northern Scotia Sea, sediment of interglacial periods is dominated by smectite with epsilon-Nd < - 6 and silt with epsilon-Nd > -4. We suggest that smectite was derived from the Falkland shelf and silt was derived from the Argentinian shelf. During glacial periods, the Argentinian shelf was an important source for silt and chlorite with epsilon-Nd > -4. The contribution from the Falkland shelf seems to have remained similar during glacial and interglacial periods. Hydrographic transport by bottom currents and turbidites could account for the high glacial detrital flux. An evaluation of the significance of an aeolian contribution to deep sea sediment suggests that it plays only a minor role. In the southern Scotia Sea, the Antarctic Peninsula is considered an important source for young material with epsilon-Nd > -4, in particular during glacial periods. During interglacial periods, sediment supply from the Antarctic Peninsula was lower than during glacial times, resulting in a relatively high contribution of old material (epsilon-Nd < -8) from East Antarctica. Deep water currents and icebergs could account for the transport of the old component to the southern Scotia Sea. The accumulation rates of material from the various source regions for glacial times are in agreement with an increase in the strength of the Antarctic Circumpolar Current. The production rate and the circulation pattern of bottom water in the Weddell Sea appear to have remained similar over most of the last 150 kyr.

Elemental composition and petroleum-generating potential of the methylene chloride extracts from DSDP Leg 66 Holes

Organic geochemical studies on samples from Holes 487, 488, and 490 in the southern Mexico Middle America Trench provided an opportunity to characterize the organic fraction of the sedimentary section in an active trench environment and to project the petroleum-producing potential of the extracted lipid fractions. The samples were geologically young and of shallow burial history. Samples from Hole 487, located on the oceanic plate, range in age from late Miocene to middlelate Pleistocene. Samples from Hole 488, representing undifferentiated Quaternary sediment, were collected on the landward side of the lower trench slope. Miocene(?) to Quaternary sediments from Hole 490 were obtained from the upper slope immediately seaward of the inferred location of the continental crust.

Strength of the Iberian intraplate lithosphere: Cenozoic deformations and seismicity

Previous studies about the strength of the lithosphere in the Iberia centre fail to resolve the depth of earthquakes because of the rheological uncertainties. Therefore, new contributions are considered (the crustal structure from a density model) and several parameters (tectonic regime, mantle rheology, strain rate) are checked in this paper to properly examine the role of lithospheric strength in the intraplate seismicity and the Cenozoic evolution. The strength distribution with depth, the integrated strength, the effective elastic thickness and the seismogenic thickness have been calculated by a finite element modelling of the lithosphere across the Central System mountain range and the bordering Duero and Madrid sedimentary basins. Only a dry mantle under strike-slip/extension and a strain rate of 10-15 s-1, or under extension and 10-16 s-1, causes a strong lithosphere. The integrated strength and the elastic thickness are lower in the mountain chain than in the basins. These anisotropies have been maintained since the Cenozoic and determine the mountain uplift and the biharmonic folding of the Iberian lithosphere during the Alpine deformations. The seismogenic thickness bounds the seismic activity in the upper–middle crust, and the decreasing crustal strength from the Duero Basin towards the Madrid Basin is related to a parallel increase in Plio–Quaternary deformations and seismicity. However, elasto–plastic modelling shows that current African–Eurasian convergence is resolved elastically or ductilely, which accounts for the low seismicity recorded in this region.

Age and nature of 560–520 Ma calc-alkaline granitoids of Biarjmand, northeast Iran: insights into Cadomian arc magmatism in northern Gondwana

The Biarjmand granitoids and granitic gneisses in northeast Iran are part of the Torud–Biarjmand metamorphic complex, where previous zircon U–Pb geochronology show ages of ca. 554–530 Ma for orthogneissic rocks. Our new U–Pb zircon ages confirm a Cadomian age and show that the granitic gneiss is ~30 million years older (561.3 ± 4.7 Ma) than intruding granitoids(522.3 ± 4.2 Ma; 537.7 ± 4.7 Ma). Cadomian magmatism in Iran was part of an approximately 100-million-year-long episode of subduction-related arc and back-arc magmatism, which dominated the whole northern Gondwana margin, from Iberia to Turkey and Iran. Major REE and trace element data show that these granitoids have calc-alkaline signatures. Their zircon O (δ18O = 6.2–8.9‰) and Hf (–7.9 to +5.5; one point with εHf ~ –17.4) as well as bulk rock Nd isotopes (εNd(t)= –3 to –6.2) show that these magmas were generated via mixing of juvenile magmas with an older crust and/or melting of middle continental crust. Whole-rock Nd and zircon Hf model ages (1.3–1.6 Ga) suggest that this older continental crust was likely to have been Mesoproterozoic or even older. Our results, including variable zircon εHf(t) values, inheritance of old zircons and lack of evidence for juvenile Cadomian igneous rocks anywhere in Iran, suggest that the geotectonic setting during late Ediacaran and early Cambrian time was a continental magmatic arc rather than back-arc for the evolution of northeast Iran Cadomian igneous rocks.

Estimation of the Crustal Bulk Properties Beneath Mainland Portugal from P-Wave Teleseismic Receiver Functions

In this work, we present results from teleseismic P-wave receiver functions (PRFs) obtained in Portugal, Western Iberia. A dense seismic station deployment conducted between 2010 and 2012, in the scope of the WILAS project and covering the entire country, allowed the most spatially extensive probing on the bulk crustal seismic properties of Portugal up to date. The application of the H-kappa stacking algorithm to the PRFs enabled us to estimate the crustal thickness (H) and the average crustal ratio of the P- and S-waves velocities V (p)/V (s) (kappa) for the region. Observations of Moho conversions indicate that this interface is relatively smooth with the crustal thickness ranging between 24 and 34 km, with an average of 30 km. The highest V (p)/V (s) values are found on the Mesozoic-Cenozoic crust beneath the western and southern coastal domain of Portugal, whereas the lowest values correspond to Palaeozoic crust underlying the remaining part of the subject area. An average V (p)/V (s) is found to be 1.72, ranging 1.63-1.86 across the study area, indicating a predominantly felsic composition. Overall, we systematically observe a decrease of V (p)/V (s) with increasing crustal thickness. Taken as a whole, our results indicate a clear distinction between the geological zones of the Variscan Iberian Massif in Portugal, the overall shape of the anomalies conditioned by the shape of the Ibero-Armorican Arc, and associated Late Paleozoic suture zones, and the Meso-Cenozoic basin associated with Atlantic rifting stages. Thickened crust (30-34 km) across the studied region may be inherited from continental collision during the Paleozoic Variscan orogeny. An anomalous crustal thinning to around 28 km is observed beneath the central part of the Central Iberian Zone and the eastern part of South Portuguese Zone.

Geochemistry of crustal formation and evolution

Terrestrial planets produce crusts as they differentiate. The Earth’s bi-modal crust, with a high-standing granitic continental crust and a low-standing basaltic oceanic crust, is unique in our solar system and links the evolution of the interior and exterior of this planet. Here I present geochemical observations to constrain processes accompanying crustal formation and evolution. My approach includes geochemical analyses, quantitative modeling, and experimental studies. The Archean crustal evolution project represents my perspective on when Earth’s continental crust began forming. In this project, I utilized critical element ratios in sedimentary records to track the evolution of the MgO content in the upper continental crust as a function time. The early Archean subaerial crust had >11 wt. % MgO, whereas by the end of Archean its composition had evolved to about 4 wt. % MgO, suggesting a transition of the upper crust from a basalt-like to a more granite-like bulk composition. Driving this fundamental change of the upper crustal composition is the widespread operation of subduction processes, suggesting the onset of global plate tectonics at ~ 3 Ga (Abstract figure). Three of the chapters in this dissertation leverage the use of Eu anomalies to track the recycling of crustal materials back into the mantle, where Eu anomaly is a sensitive measure of the element’s behavior relative to neighboring lanthanoids (Sm and Gd) during crustal differentiation. My compilation of Sm-Eu-Gd data for the continental crust shows that the average crust has a net negative Eu anomaly. This result requires recycling of Eu-enriched lower continental crust to the mantle. Mass balance calculations require that about three times the mass of the modern continental crust was returned into the mantle over Earth history, possibly via density-driven recycling. High precision measurements of Eu/Eu* in selected primitive glasses of mid-ocean ridge basalt (MORB) from global MORs, combined with numerical modeling, suggests that the recycled lower crustal materials are not found within the MORB source and may have at least partially sank into the lower mantle where they can be sampled by hot spot volcanoes. The Lesser Antilles Li isotope project provides insights into the Li systematics of this young island arc, a representative section of proto-continental crust. Martinique Island lavas, to my knowledge, represent the only clear case in which crustal Li is recycled back into their mantle source, as documented by the isotopically light Li isotopes in Lesser Antilles sediments that feed into the fore arc subduction trench. By corollary, the mantle-like Li signal in global arc lavas is likely the result of broadly similar Li isotopic compositions between the upper mantle and bulk subducting sediments in most arcs. My PhD project on Li diffusion mechanism in zircon is being carried out in extensive collaboration with multiple institutes and employs analytical, experimental and modeling studies. This ongoing project, finds that REE and Y play an important role in controlling Li diffusion in natural zircons, with Li partially coupling to REE and Y to maintain charge balance. Access to state-of-art instrumentation presented critical opportunities to identify the mechanisms that cause elemental fractionation during laser ablation inductively coupled plasma mass spectrometry (LA-ICP-MS) analysis. My work here elucidates the elemental fractionation associated with plasma plume condensation during laser ablation and particle-ion conversion in the ICP.

Deep Crustal Structure of the North-West African Margin from Combined Wide-angle and Reflection Seismic Data (MIRROR Seismic Survey)

The structure of the Moroccan and Nova Scotia conjugate rifted margins is of key importance for understanding the Mesozoic break-up and evolution of the northern central Atlantic Ocean basin. Seven combined multichannel reflection (MCS) and wide-angle seismic (OBS) data profiles were acquired along the Atlantic Moroccan margin between the latitudes of 31.5° and 33° N during the MIRROR seismic survey in 2011, in order to image the transition from continental to oceanic crust, to study the variation in crustal structure and to characterize the crust under the West African Coast Magnetic Anomaly (WACMA). The data were modeled using a forward modeling approach. The final models image crustal thinning from 36 km thickness below the continent to approximately 8 km in the oceanic domain. A 100 km wide zone characterized by rough basement topography and high seismic velocities up to 7.4 km/s in the lower crust is observed westward of the West African Coast Magnetic Anomaly. No basin underlain by continental crust has been imaged in this region, as has been identified north of our study area. Comparison to the conjugate Nova Scotian margin shows a similar continental crustal thickness and layer geometry, and the existence of exhumed and serpentinized upper mantle material on the Canadian side only. The oceanic crustal thickness is lower on the Canadian margin.