(Table) Argon isotope ratios and argon ages of rocks obtained from the Chernorud zone, western Baikal area

Structural-petrologic and isotopic-geochronologic data on magmatic, metamorphic, and metasomatic rocks from the Chernorud zone were used to reproduce the multistage history of their exhumation to upper crustal levels. The process is subdivided into four discrete stages, which corresponded to metamorphism to the granulite facies (500-490 Ma), metamorphism to the amphibolite facies (470-460 Ma), metamorphism to at least the epidote-amphibolite facies (440-430 Ma), and postmetamorphic events (410-400 Ma). The earliest two stages likely corresponded to the tectonic stacking of the backarc basin in response to the collision of the Siberian continent with the Eravninskaya island arc or the Barguzin microcontinent, a process that ended with the extensive generation of synmetamorphic granites. During the third and fourth stages, the granulites of the Chernorud nappe were successively exposed during intense tectonic motions along large deformation zones (Primorskii fault, collision lineament, and Orso Complex). The comparison of the histories of active thermal events for Early Caledonian folded structures in the Central Asian Foldbelt indicates that active thermal events of equal duration are reconstructed for the following five widely spiced accretion-collision structures: the Chernorud granulite zone in the Ol'khon territory, the Slyudyanka crystalline complex in the southwestern Baikal area, the western Sangilen territory in southeastern Tuva, Derbinskii terrane in the Eastern Sayan, and the Bayankhongor ophiolite zone in central Mongolia. The dates obtained by various isotopic techniques are generally consistent with the four discrete stages identified in the Chernorud nappe, whereas the dates corresponding to the island-arc evolutionary stage were obtained only for the western Sangilen and Bayankhongor ophiolite zone.

(Table 1) Basement penetration, recovery, and lithology of crustal DSDP holes in the North Atlantic Ocean

Oceanic crustal drilling by R. V. Glomar Challenger at 15 sites in the North Atlantic has led to a complex picture of the upper half kilometer of the crust. Elements of the picture include the absence of the source for linear magnetic anomalies, marked episodicity of volcanic activity, ubiquitous low temperature alteration and evidence for large scale tectonic disturbance. Comparison sections in the Pacific and much deeper crustal drilling are needed to attack problems arising from the North Atlantic results.

Structural evolution of the Lesser Himalayan fold-thrust belt of west central Nepal

The Lesser Himalayan fold-thrust belt on the south flank of the Jajarkot klippe in west central Nepal was mapped in detail between the Main Central thrust in the north and the Main Boundary thrust in the south. South of the Jajarkot klippe, the fold-thrust belt involves sandstone, shale and carbonate rocks that are unmetamorphosed in the foreland and increase in metamorphic grade with higher structural position to sub-greenschist facies towards the hinterland. The exposed stratigraphy is correlative with the Proterozoic Ranimata, Sangram, Galyang, Syangia Formations and Lakharpata Group of Western Nepal and overlain by the Paleozoic Tansen and Kali Gandaki Groups. Based on field mapping and cross-section construction, three distinct thrust sheets were identified separated by top-to-the-south thrust faults. From the foreland (south) to the hinterland (north), the first thrust sheet in the immediate hanging wall of the Main Boundary thrust defines an open syncline. The second thrust sheet contains a very broad synformal duplex, which is structurally stacked against the third thrust sheet containing a homoclinal panel of the oldest exposed Proterozoic stratigraphy. Outcrop scale folds throughout the study area are predominantly south vergent, open, and asymmetric reflecting the larger regional scale folding style, which corroborate the top-to-the-south deformation style seen in the faults of the region. Field techniques were complemented with microstructural and quartz crystallographic c-axis preferred orientation analyses using a petrographic microscope and a fabric analyzer, respectively. Microstructural analysis identified abundant strain-induced recrystallization textures and occasional occurrences of top-to-the-south shear-sense indicators primarily in the hinterland rocks in the immediate footwall of the Main Central Thrust. Top-to-the-south shearing is also supported by quartz crystallographic c-axis preferred orientations. Quartz recrystallization textures indicate an increase in deformation temperature towards the Main Central thrust. A line balance estimate indicates that approximately 15 km of crustal shortening was accommodated by folding and faulting in the fold-thrust belt south of the Jajarkot klippe. Additionally, estimations of shortening velocity suggest that the shortening velocity operating in this section of the fold-thrust belt between 23 to 14 Ma was slower than what is currently observed as a result of the ongoing deformation of the Sub-Himalayan fold-thrust belt.

Stress Inversion Using LiDAR Tunnel Deformation Measurements

Far-field stresses are those present in a volume of rock prior to excavations being created. Estimates of the orientation and magnitude of far-field stresses, often used in mine design, are generally obtained by single-point measurements of stress, or large-scale, regional trends. Point measurements can be a poor representation of far-field stresses as a result of excavation-induced stresses and geological structures. For these reasons, far-field stress estimates can be associated with high levels of uncertainty. The purpose of this thesis is to investigate the practical feasibility, applications, and limitations of calibrating far-field stress estimates through tunnel deformation measurements captured using LiDAR imaging. A method that estimates the orientation and magnitude of excavation-induced principal stress changes through back-analysis of deformation measurements from LiDAR imaged tunnels was developed and tested using synthetic data. If excavation-induced stress change orientations and magnitudes can be accurately estimated, they can be used in the calibration of far-field stress input to numerical models. LiDAR point clouds have been proven to have a number of underground applications, thus it is desired to explore their use in numerical model calibration. The back-analysis method is founded on the superposition of stresses and requires a two-dimensional numerical model of the deforming tunnel. Principal stress changes of known orientation and magnitude are applied to the model to create calibration curves. Estimation can then be performed by minimizing squared differences between the measured tunnel and sets of calibration curve deformations. In addition to the back-analysis estimation method, a procedure consisting of previously existing techniques to measure tunnel deformation using LiDAR imaging was documented. Under ideal conditions, the back-analysis method estimated principal stress change orientations within ±5° and magnitudes within ±2 MPa. Results were comparable for four different tunnel profile shapes. Preliminary testing using plastic deformation, a rough tunnel profile, and profile occlusions suggests that the method can work under more realistic conditions. The results from this thesis set the groundwork for the continued development of a new, inexpensive, and efficient far-field stress estimate calibration method.

THE MICROSTRUCTURES OF ZR-EXCEL ALLOY: EFFECTS OF HEAVY ION IRRADIATION, PLASTIC DEFORMATION, THERMAL TREATMENT AND IN-SITU HEATING

Zr-Excel alloy (Zr-3.5Sn-0.8Nb-0.8Mo) is a dual phase (α + β) alloy in the as-received pressure tube condition. It has been proposed to be the pressure tube candidate material for the Generation-IV CANDU-Supercritical Water Reactor (CANDU-SCWR). In this dissertation, the effects of heavy ion irradiation, deformation and heat treatment on the microstructures of the alloy were investigated to enable us to have a better understanding of the potential in-reactor performance of this alloy. In-situ heavy ion (1 MeV) irradiation was performed to study the nucleation and evolution of dislocation loops in both α- and β-Zr. Small and dense type dislocation loops form under irradiation between 80 and 450 °C. The number density tends to saturate at ~ 0.1 dpa. Compared with the α-Zr, the defect yield is much lower in β-Zr. The stabilities of the metastable phases (β-Zr and ω-Zr) and the thermal-dynamically equilibrium phase, fcc Zr(Mo, Nb)2, under irradiation were also studied at different temperatures. Chemi-STEM elemental mapping was carried out to study the elemental redistribution caused by irradiation. The stability of these phases and the elemental redistribution are strongly dependent on irradiation temperature. In-situ time-of-flight neutron diffraction tensile and compressive tests were carried out at different temperatures to monitor lattice strain evolutions of individual grain families during these tests. The β-Zr is the strengthening phase in this alloy in the as-received plate material. Load is transferred to the β-Zr after yielding of the α-Zr grains. The temperature dependence of static strain aging and the yielding sequence of the individual grain families were discussed. Strong tensile/compressive asymmetry was observed in the {0002} grain family at room temperature. The microstructures of the sample deformed at 400 °C and the samples only subjected to heat treatment at the same temperature were characterized with TEM. Concentration of β phase stabilizers in the β grain and the morphology of β grain have significant effect on the stability of β- and ω-Zr under thermal treatment. Applied stress/strain enhances the decomposition of isothermal ω phase but suppresses α precipitation inside the β grains at high temperature. An α → ω/ZrO phase transformation was observed in the thin foils of Zr-Excel alloy and pure Zr during in-situ heating at 700 °C in TEM.

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.

Time-Dependent Deformation of Non-Composite and Composite Sand-Lightweight Prestressed Concrete Structures, Phase 1, HR-137, 1969

Presented in this report is an investigation of the use of "sand-lightweight" concrete in prestressed concrete structures. The sand-lightweight concrete consists of 100% sand substitution for fines, along with Idealite coarse and medium lightweight aggregate and Type I Portland Cement.

U-Pb SHRIMP zircon dating of high-grade rocks from the Upper Allochthonous Terrane of Bragança and Morais Massifs (NE Portugal); geodynamic consequences

Bragança and Morais Massifs are part of the mega-klippen ensemble of NW Iberia, comprising a tectonic pile of four allochthonous units stacked above the Central-Iberian Zone autochthon. On top of this pile, the Upper Allochthonous Terrane (UAT) includes different high-grade metamorphic series whose age and geodynamic meaning are controversial. Mafic granulites provided U–Pb zircon ages at 399±7 Ma, dating the Variscan emplacement of UAT. In contrast,U–Pb zircon ages of ky- and hb-eclogites, felsic/intermediate HP/HT-granulites and orthogneisses (ca. 500–480 Ma) are identical to those of gabbros (488 ± 10 Ma) and Grt-pyroxenites (495 ± 8 Ma) belonging to a mafic/ultramafic igneous suite that records upper mantle melting and mafic magma crustal underplating at these times. Gabbros intrude the high-grade units of UAT and did not underwent the HP metamorphic event experienced by eclogites and granulites. These features and the zircon dates resemblance among different lithologies, suggest that extensive age resetting of older events may have been correlative with the igneous suite emplacement/crystallisation. Accordingly, reconciliation of structural, petrological and geochronological evidence implies that the development and early deformation of UAT high-grade rocks should be ascribed to an orogenic cycle prior to ≈500 Ma. Undisputable dating of this cycle is impossible, but the sporadic vestiges of Cadomian ages cannot be disregarded. The ca. 500–480 Ma time-window harmonises well with the Lower Palaeozoic continental rifting that trace the VariscanWilson Cycle onset and the Rheic Ocean opening. Subsequent preservation of the high heat-flowregime, possibly related to the Palaeotethys back-arc basin development (ca. 450–420 Ma), would explain the 461 ± 10 Ma age yielded by some zircon domains in felsic granulites, conceivably reflecting zircon dissolution/ recrystallisation till Ordovician times, long before the Variscan paroxysm (ca. 400–390 Ma). This geodynamic scenario suggests also that UAT should have been part of Armorica before its emplacement on top of Iberia after Palaeotethys closure.

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.

Micromechanical insight into the undrained instability of granular materials: deformation characteristics of geomaterials.

There is no agreement between experimental researchers whether the point where a granular material responds with a large change of stresses, strains or excess pore water pressure given a prescribed small input of some of the same variables defines a straight line or a curve in the stress space. This line, known as the instability line, may also vary in shape and position if the onset of instability is measured from drained or undrained triaxial tests. Failure of granular materials, which might be preceded by the onset of instability, is a subject that the geotechnical engineers have to deal with in the daily practice, and generally speaking it is associated to different phenomena observed not only in laboratory tests but also in the field. Examples of this are the liquefaction of loose sands subjected to undrained loading conditions and the diffuse instability under drained loading conditions. This research presents results of DEM simulations of undrained triaxial tests with the aim of studying the influence of stress history and relative density on the onset of instability in granular materials. Micro-mechanical analysis including the evolution of coordination numbers and fabric tensors is performed aiming to gain further insight on the particle-scale interactions that underlie the occurrence of this instability. In addition to provide a greater understanding, the results presented here may be useful as input for macro-scale constitutive models that enable the prediction of the onset of instability in boundary value problems.

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.

Upper crustal structure of Madeira Island revealed from ambient noise tomography

We present the first image of the Madeira upper crustal structure, using ambient seismic noise tomography. 16 months of ambient noise, recorded in a dense network of 26 seismometers deployed across Madeira, allowed reconstructing Rayleigh wave Green's functions between receivers. Dispersion analysis was performed in the short period band from 1.0 to 4.0 s. Group velocity measurements were regionalized to obtain 2D tomographic images, with a lateral resolution of 2.0 km in central Madeira. Afterwards, the dispersion curves, extracted from each cell of the 2D group velocity maps, were inverted as a function of depth to obtain a 3D shear wave velocity model of the upper crust, from the surface to a depth of 2.0 km. The obtained 3D velocity model reveals features throughout the island that correlates well with surface geology and island evolution.