Crustal Scale Research Articles

Tectonic anabranching of Brahmaputra River system: implications on survival of world’s largest inhabited river island Majuli

Brahmaputra is a large, tropical, trans-boundary river with high sediment load and has the unique distinction of possessing both largest and smallest inhabited river islands. Majuli is the largest inhabited river island situated in the Indian state of Assam and bounded by Brahmaputra on the south, Subansiri River on the northwest, and the Kherkatia River (anabranch of Brahmaputra) on the northeast. The island was formed during a flooding event in 1750 with a sudden change in the course of the Brahmaputra from a low energy meandering to a high energy braided river and subsequent capture of one of its tributaries. The island is a part of the Brahmaputra alluvium squeezed between active Himalayan and Indo-Burmese orogenic belts in an intricate tectonic regime. The dominant flow path of Brahmaputra is controlled by the crustal scale reorganization in this part of the world. The tectonic factors that controlled the anabranching process and created the island, is at present acting as the destroyer of the same by the means of unabated erosion. The erosion pattern of the two rivers viz. Brahmaputra and Subansiri were studied and areas with marked selective erosion have been found. From a combined study of satellite imaging, bathymetric survey and mapping, it is found that the lobate shaped erosions in the study areas are restricted to traces of the tectonic elements related to the Himalayan thrust system through which the rivers are trying to anabranch.

Geology of the Golobar and Rumaceo faults, Basque-Cantabrian Pyrenees

ABSTRACT The lateral transition between different styles of deformation is crucial for understanding local-scale structures since it influences orogenic reconstructions at the crustal scale. This is especially important at the westernmost Basque-Cantabrian Pyrenees, where the transition from the thin-skinned Basque-Cantabrian Pyrenees to the thick-skinned Cantabrian Mountains remains a subject of debate. To address this, detailed geological mapping and proper fault characterisation are critical. This research focused on the surface characterisation of the Golobar Fault System, one of the faults located at the transition between the two styles of deformation. The results of this study include not only a detailed geological map of the Golobar Fault System but also its integration with neighbouring seismic sections and well data. This approach enhances the understanding of the deeper portions and lateral thickness variations, thereby providing a more comprehensive understanding of the Golobar area within the regional framework.

Shaping the crustal structure of the SW-Alpine Foreland: Insights from 3D Geological modeling

Reactivation processes play a significative role in the localization of deformation but still remain hard to establish at the lithospheric scale. In this work, we built a 3D structural model, which enables to bridge the gap between the main tectonic structures observed at the surface and the geometry of the major interfaces (the Mohorovičić-discontinuity (hereafter Moho) and top of the basement) inferred from geophysical data acquired in the external Western Alps and their foreland. The geometry of these tectonic structures is interpreted in relation to their geodynamic evolution. The main results of this study highlight: (1) a strong contribution of thick-skinned Pyrenean-Provence and Alpine tectonics, (2) a lithospheric rooting of Variscan shear zones and related faults, and (3) the regional-scale influence of these inherited structures on the post-Paleozoic strain localization in the study area. Our 3D model shows that the pattern of Variscan shear zones that were developed at the end of the Paleozoic involved the whole crust, as emphasized by the Moho offsets. These shear zones were reactivated and localized Meso-Cenozoic deformation. The Variscan deformation pattern controlled the geometry of extensional basins, the propagation of Pyrenean-Provence deformation, and finally the Alpine deformation at crustal scale. Our 3D model shows minor crustal thickening (ca. 40 km) located below the Pelvoux External Crystalline Massif, which probably resulted from both Pyrenean and Alpine tectonic phases. In contrast, the southern part of the Alpine front shows a thinned crust (ca. 18 km) resulting from extensional Meso-Cenozoic phases between the Cévennes margin and the Durance basin.

Geologically constrained geometry inversion and null-space navigation to explore alternative geological scenarios: a case study in the Western Pyrenees

SUMMARY Reducing the gap between geophysical inversion and geological interpretation can be achieved by integrating geological modelling into geophysical inversion. For this, we use a generalized, iterative level-set gravity inversion scheme in which geological units are deformed automatically. During the inversion process, a regularization term is defined using automated geological modelling to account for geological data and principles. This provides model-dependent geological constraints and encourages geological realism throughout inversion. To alleviate the dependence on the starting model and consider the possibility of features unseen by direct observations, an automated geophysical data-driven method is proposed to insert new rock units in the model. Uncertainty quantification is achieved through the null-space shuttle algorithm, which is used to generate a series of alternative models that are consistent with geophysical data. This methodology is applied to assess the uncertainties of a pre-existing 3-D crustal-scale geological model of the Western Pyrenean orogeny (France, Spain). The area is characterized by a positive gravity anomaly generally attributed to the presence of a shallow mantle body. The impact of variations in shape and density of key crustal and mantle features is investigated. Different scenarios are explored in 3-D space to produce a range of viable, relatively simple crustal-scale models of the area. This application demonstrates the capability and potential of this approach to evaluate alternative interpretations of geophysical data. The results show the plausibility of scenarios with a shorter subducted Iberian lower crust and a denser Axial Zone than in the pre-existing model.

The 3D stress field of Nordland, northern Norway – insights from numerical modelling

Abstract The Nordland area in northern Norway is the seismically most active area on mainland Fennoscandia. It exhibits patterns of coastal extension, which contrasts with the first-order regional stress pattern that reflects compression aligned with the North Atlantic ridge push. The regional stress field has been considered to emanate from the interaction of ridge push and glacial isostatic adjustment; while the local stress pattern can be additionally influenced by gravitational, topographic stresses, as well as the flexural effects of erosion and sediment deposition. We employ finite element numerical models at a crustal scale to study the 3D stress field, using existing geometric constraints from previous geophysical studies. Internal body forces, induced by variations in density, topography or Moho depth, already yield significant deviatoric stresses. In the models tested, these can strongly influence the near-surface stress regime, in particular for the continental margin setting we are investigating. In addition, redistribution of rock mass, which occurred mainly under Pleistocene glaciation, can modify the stress field significantly on a semi-regional scale. We consider this process to be the main driver for the coastal extension, in particular in areas where erosion has been high.

Crucial role of water-present melting in metagranite: Implications for the instigation of crustal-scale shear zones

Abstract Where, when, and why large-scale shear zones nucleate and propagate into the continental lithosphere are critical issues that challenge the research in tectonics. The East Variscan shear zone is one of the crustal-scale strike-slip faults that shaped the Variscan orogenic crust during late Carboniferous time. Field-based structural analysis and petrological observations demonstrate that suprasolidus high-strain deformation zones and metagranite occurrences are spatially correlated. Among the three dominant lithologies forming this orogenic middle crust (metapelite, metagraywacke, and metagranite), petrological observations and phase equilibrium modeling indicate that the latter is the first lithology that melts during collision-induced heating, in response to H2O-fluid-saturated melting. Our field data and modeling suggest that the water-fluxed melting of metagranite has a primary rheological control on the localization, instigation, and growth of crustal-scale shear zones in the middle crust. Thus, the distribution and geometry of metagranite at the crustal scale could be regarded as critical parameters influencing the rheological inheritance governing the tectonic evolution and localization of bulk strain in the continental lithosphere.

Detailed crustal deformation in the Ailaoshan Orogenic Belt revealed by receiver functions from a dense array

The Ailaoshan Orogenic Belt (AOB), located at the southeastern boundary of the Tibetan Plateau, is an ideal place for investigating the mechanisms of lateral growth of Tibet. Using the data recorded by a dense seismic array across the Ailaoshan belt, we investigate the detailed lateral variations of crustal anisotropy on the basis of Pms phase of receiver functions. Remarkable crustal anisotropy is observed throughout this study region with a mean delay time of 0.33 ± 0.19 s, indicating the anisotropy primarily originates in the middle-lower crust. The fast directions beneath the AOB including the Ailaoshan-Red River shear zone (ARRSZ) and its western low-grade metamorphic unit generally align with the NW-SE strike of ARRSZ. The weak anisotropy in the South China Block (SCB) argues that the block is relatively stable, with limited internal deformation. Meanwhile, the anisotropy beneath the western boundary of the SCB is strong, and the N-S oriented fast direction is influenced by both the crustal stress and Xiaojiang Fault. Combining the high Vp/Vs and significant lateral variations of crustal anisotropy parameters, we suggest that the strike-slip motion along the ARRSZ induces the partial melting and pronounced anisotropy in the middle-lower crust of AOB, without the presence of crustal flow. The differences between crustal and mantle anisotropy indicate crust-mantle decoupling deformation of the AOB, supporting the block extrusion model occurring only in the crustal scale as the primary deformation pattern.

Fluid-assisted intra-plate seismicity at the edge of the Gawler Craton, South Australia

The Australian continent, being void of active plate boundaries, is often perceived as seismically quiescent. However, earthquakes of moderate magnitude (M6+) occur on the continent around once per decade. Such intra-plate activity can pose a significant risk as these earthquakes can occur along fault lines that are either unknown or considered inactive, are often non-periodic, and poorly understood. Within Australia, the spatial distribution of intra-plate seismicity is non-uniform, instead tending to concentrate along certain weak zones of increased activity. One such region is the eastern margin of the Gawler Craton in South Australia, one of the oldest building blocks of the continent. Recently, several new temporary seismic arrays have been deployed in the region, transforming data coverage across southern Australia. In total, 139 new local earthquakes have been recorded, most of which went undetected by the national seismic network. Following relocation, the pattern of earthquakes becomes more localised and appears to coincide with the edge of the Gawler Craton. Further, a spatial association was found between earthquakes and mound springs, which act as the discharge point for groundwater migrating from the Great Artesian Basin. Enhanced fluid pressures (mantle-degassing) within permeable crustal scale fault systems, responding to a regional contractional strain field, appear to be a key driver of seismicity in the area.

Regional structural setting of late-orogenic IOCG mineralization along the northern Nautanen deformation zone, Norrbotten, Sweden

The northern Norrbotten ore province in Sweden is one of the most mineralized areas in Europe. Iron, gold and/or copper deposits occur as iron oxide-apatite-style mineralization (IOA, Kiruna-type) as well as iron oxide-copper–gold (IOCG) style. Regardless of mineralization style, most deposits appear to be spatially controlled by a set of crustal-scale Palaeoproterozoic shear zones which share similar structural characteristics and deformation histories.Reappraisal of regional geological and geophysical data, coupled with structural mapping, suggests crustal-scale shear zones form continuous c. N-S-trending zones extending from the Skellefte district in the south into the northern Norrbotten ore province. One example from Norrbotten is a zone that extends SSW from Karesuando in the north towards Svappavaara. While this structure has traditionally been inferred to continue SW towards Arjeplog (i.e. the Karesuando – Arjeplog Deformation Zone; KADZ), we favour its deflection SSE into the Nautanen-Aitik trend, making it a continuous, IOCG-bearing, crustal-scale deformation zone. Similar shear zone geometries can be observed in analogous zones to the west. Most of these crustal scale structures record at least two time-separated deformation events of regional significance. IOA and IOCG deposits form in different tectonic environments, separated in time and overprinting each other.

Lateral changes of crustal extension and passive margin type along the Brazilian southeastern margin

For the last half century, studies based on an increasing, diverse data set, have focused on passive continental margin evolution, as a result of a sequence of tectonic processes that occur in crustal scale and in the geological time. Where extensive subsurface data exists, two distinct endmembers of continental margin architectures were first described along specific margins of the North Atlantic: 1) wide continent-ocean transitions or hyper-extended; 2) narrow continent-ocean transitions or hypo-extended. The lateral transition between these endmember margins, which can occur laterally in short distances, is still not fully understood. The same two endmember margin types are observed across passive margins around the world and notably along the South Atlantic margin. In an area known for controversial interpretations about the crustal nature and the limits of thinned continental crust, this investigation integrates crustal thickness, seismic interpretation, and facies analyses across ∼150,000 km of seismic data along the Brazilian southeastern margin. This study has implications for those investigating the crustal geometry variations for basin analysis and results impact hydrocarbon assessments for the southeastern margin of Brazil. We resume continental margin analyses, in a critical moment as hydrocarbon exploration advances from the continental slope to ultra-deep waters. Results indicate a marked change from a narrow, hypo-extended, sub-aerial, Iceland-like, plume-related volcanic crust in the Pelotas Basin, to a hyper-extended, Iberia-like, magma-rich crust in the Santos Basin, which is separated by a pronounced oceanic transform-fault zone. Dextral movement along this NW-oriented fault zone accommodated differential extensional rates between the hypo- and the hyper-extended margins. Lateral variations in magmatic content within these margin types are interpreted as result from the interaction with mantellic plumes. Margin architecture can locally be affected by pre-existing fault zones and be later modified by oceanic transform faults. The interactions between crustal extensional rates, crustal rheology, mantle underplating, and volcanic material exhumed through mantle-derived plumes, are the key controls for the evolution of continental margins. Tectonic framework classification proposed in this study presents an alternative, original model for continental passive margins evolution.

A New Lithospheric Density and Magnetic Susceptibility Model of Iran, Starting From High‐Resolution Seismic Tomography

AbstractWe propose a new model for the crust and upper mantle in Iran by joint inversion of gravity and magnetic fields, constrained with a seismic tomography model. We then calculate shear modulus from the Vs velocities and densities. The crust and mantle tomography model is first converted to a density cube through empirical and petrological velocity‐density relations. The starting susceptibility is assigned to a two‐layer homogeneous model, above a heat flow‐derived Curie depth. Considering the uncertainties in the density‐velocity relations, and the starting layered susceptibility variation, we refine the model by a constrained inversion of the gravity and magnetic fields with a Bayesian approach, producing the final 3D density and susceptibility model. The area is tectonically active with high seismicity and active faulting which are regulated by the crustal density and rigidity variations. Higher rigidity matches lower seismicity and extended deserts and basins, suggesting the control of their development. The Neo‐Tethys suture, extending ∼1,500 km long, as well as the Paleo‐Tethys suture match crustal scale density variations, defining characteristics of the lighter Arabian plate and denser Eurasian crust. The South Caspian Basin is enigmatic, due to focusing on the seismicity along all its borders, but with relatively low average rigidity, which is contrary to what is observed for Iran, where the reduced rigidity correlates with higher seismicity. The 3D density model will be useful for numerical geodynamic models and obtaining geologic inferences from the crustal‐scale units.

A database of detrital zircon U–Pb geochronology and Hf isotopes from the Songpan–Ganzi and Western Qinling terranes

AbstractSince the advent of detrital geochronology, significant progress has been made in the study of sedimentary sources. The detrital zircon U–Pb age spectrum, especially when combined with zircon Hf isotopes, provides a powerful tool for determining the provenance location, nature, and evolutionary history of sediments. The Western Qinling–Songpan tectonic junction, located in central China, is a giant structure on the crustal scale. A total of 20,640 zircon U–Pb ages of 281 sandstone samples from the Songpan–Ganzi and Western Qinling terranes were obtained. Among these, 1,636 zircons from 33 samples contained Hf isotopic data. Moreover, we compiled data on papers, samples, analytical methods, and results and obtained a total of 76 zircon U–Pb and 60 Hf isotope data heads. The obtained dataset has multiple data interfaces, which can be well associated with research or applications in different fields, for provenance analysis, climate change, and paleogeographic reconstruction studies.

Deep crustal structures with reverse time migration applied to offshore wide-angle seismic data: Equatorial and North-West Brazilian margins

We image the Moho discontinuity and deep crustal layers by applying the Reverse Time Migration (RTM) method to wide-angle seismic (WAS) data that were acquired along two different profiles in NW and SE offshore Brazil, where the spacing between OBS is 12.5 km and 150 m between shots. The application of this method is quite uncommon to ocean bottom seismometers (OBS) data due to the OBS wide spacing deployment and low folds. We analyze the effectiveness of the RTM method when applied to the reflectivity of the WAS data. We imaged the structures by cross-correlating the forward and backward wavefields given by the acoustic seismic equation. This allowed us to use each interface crossed by a ray as both a source and a receiver. The velocity models used to perform RTM were previously obtained by applying a procedure of two-dimensional forward ray-tracing followed by a damped least-squares travel time inversion. The results obtained have an unexpected large contribution from the wavefield traveling as refractions within the Earth and, because of that, we can talk about recovering refractors from the RTM results. We obtained strong and continuous refractors for depths of 7–15 km that correspond to the basement and the Moho discontinuity. As we move landwards, the refractors that correspond to the Moho discontinuity disappear due to the deepening of the surface. Finally, we investigate how the slight change on the velocity model influences the result of the RTM method. The obtained results are promising for a wide range of applications at a crustal scale seismic exploration, with wide-angle seismic data.

Fold interference pattern and crustal decoupling in northern Tanzania

High-temperature granulite facies metamorphism and deformation evolved through Neoproterozoic tectonic processes in large parts of southeastern Africa. The Eastern Granulites of Tanzania, one example for these extreme conditions, show progressive decoupling between lower and middle crust together with progressive cooling and rheology stratification. Coupled deformation during earliest stages of deformation arises from extremely small viscosity contrast of all minerals and rocks at granulite facies metamorphic conditions around 850 °C. With cooling of the orogen rheology stratification increases and deformation decouples into horizontal stretching through gravity driven channel flow in the lower crust and thickening and horizontal west-east shortening in the middle/upper crust. As soon as a stiff and cold indenter, represented by external parts of the Archean Tanzania Craton (Western Granulites) enters the system, the indenter acts as lower crustal ramp along which deep crustal units are exhumed and displaced towards the foreland. During this orogenic period the crust was fully stratified and a second phase of decoupling, detachment over crustal scale ramp, occurs. Superposition with late orogenic north-south shortening structures resulted in type 1 fold interference pattern and formation a of dome-and-basin geometry in northeastern Tanzania. Structural modelling of fold superposition helps to explain the complex structural pattern and to predict sites of decollement horizons. By comparison with fold interference pattern exposed in southern Madagascar we define different types of dome-and-basin formation. Those with neglectable rheology contrast are viscosity determined where mixing of upper and lower crustal material is enabled by vertical flow folding. Those with strong viscosity contrast are mechanically determined where domes and basins are formed by changes in regional principal stress directions.

Timescales of ultra-high temperature metamorphism and crustal differentiation: Zircon petrochronology from granulite xenoliths of the Variscan French Massif Central

Lower crustal ultra-high temperature (UHT) metamorphism and the resulting production and transfer of granitic magmas represent key processes of intracrustal differentiation. The timescales of these phenomena are debated, due to the complexity of the granulite zircon U-Pb record and because direct links between granulites and granites are difficult to establish. To address these issues, we present the results of zircon petrochronology (coupled U-Pb dating, trace element and Lu-Hf isotopic analyses) for lower crustal felsic and mafic granulite xenoliths from the Variscan eastern French Massif Central and compare them with data from well-characterized mid-/upper crustal migmatites and granites. The felsic and mafic granulites represent pre-Variscan meta-sedimentary and meta-igneous mafic rocks, respectively; which experienced Variscan UHT peak metamorphism at 940–970°C and 8 ± 2 kbar. Zircons from the felsic granulites show U-Pb dates spreading over ∼50 Myr between ∼315 and ∼265 Ma, correlated to Ti-in-zircon temperatures decreasing from 940–970°C to 800°C and REE contents consistent with growth in equilibrium with garnet at (U)HT conditions. This is best explained by continuous crystallization of zircon upon cooling from the thermal peak owing to decreasing Zr solubility of residual melt, despite significant prograde melt loss. Zircons from the mafic granulite only record the last stage of this time-temperature evolution (295–265 Ma; 800–900°C) due to later zircon saturation. Upper crustal granite emplacement started at 340 Ma and culminated at the age of the lower crustal thermal peak of 313 ± 3 Ma (defined by the highest-temperature zircons from felsic granulites), reflecting melt extraction along the ∼27 Myr prograde path of the lower crust. In turn, the crystallization ages of mid-crustal migmatites (315–300 Ma) and the lower crustal granulites (315–265 Ma) are consistent with slow cooling in the presence of melt. These results provide a direct assessment of the timescales of melt production and residence at the crustal scale; and validate the granulite-granite connection in the framework of the melt loss theory in migmatitic systems.

The lithospheric S-wave velocity structure beneath the NE Tibetan Plateau and its surrounding craton basins

It is essential to investigate the spatial distribution of the lithosphere and asthenosphere in detail, to further obtain the understanding of the effect of plate collision and the process of orogenic movement. From the joint inversion of receiver functions and surface waves, the three-dimensional S-wave velocity structure results down to 200 km depth in the study area were obtained at 1,843 seismic stations. Analysis was performed on the sedimentary thickness, crustal thickness, lower crustal wave velocity, and lithospheric thickness. According to the crustal thickness, we evaluated the distribution of low-velocity zones in the lower crust. The results show that there are low-velocity bodies in the lower crust in the Qinling tectonic belt, but they are not connected, indicating that they may not be able to be used as a channel for material extrusion from the NE Tibetan Plateau at the crustal scale. According to the section results and the depth distribution of the lithosphere-astenosphere boundary, a relatively thick lithosphere exists below the Sichuan Basin and Ordos Basin, and the lithosphere in the east of the study area is relatively thin with a thickness of about 60–80 km, indicating that the lithosphere in the east of the study area has been severely destructed and restructured. The delamination has been observed in the lithosphere under the Songpan-Ganzi Block, showing characteristics of vertical movement of asthenosphere materials. There is a relatively thick low-velocity zone at the top of the mantle lithosphere of the NE plateau; however, it does not exist under the relatively stable Sichuan Basin and the Ordos Block. Compared with the Sichuan Basin and the Ordos Basin at both sides, the Qinling tectonic belt has a low-velocity zone at the depth of 100–160 km, which may be asthenosphere material. In combination with the polarization direction characteristics of the SKS wave, it is clearly observed that asthenospheric material movement exists in an approximate east-west direction beneath the Qinling tectonic belt. Therefore, the asthenosphere beneath the Qinling tectonic belt may serve as an important channel for material extrusion in the NE Tibetan Plateau.

Rapid construction of Rayleigh wave dispersion curve based on deep learning

Introduction:The dispersion curve of the Rayleigh-wave phase velocity (VR) is widely utilized to determine site shear-wave velocity (Vs) structures from a distance of a few metres to hundreds of metres, even on a ten-kilometre crustal scale. However, the traditional theoretical-analytical methods for calculating VRs of a wide frequency range are time-consuming because numerous extensive matrix multiplications, transfer matrix iterations and the root searching of the secular dispersion equation are involved. It is very difficult to model site structures with many layers and apply them to a population-based inversion algorithm for which many populations of multilayers forward modelling and many generations of iterations are essential.Method:In this study, we propose a deep learning method for constructing the VR dispersion curve in a horizontally layered site with great efficiency. A deep neural network (DNN) based on the fully connected dense neural network is designed and trained to directly learn the relationships between Vs structures and dispersion curves. First, the training and validation sets are generated randomly according to a truncated Gaussian distribution, in which the mean and variance of the Vs models are statistically analysed from different regions’ empirical relationships between soil Vs and its depth. To be the supervised dataset, the corresponding VRs are calculated by the generalized reflection-transmission (R/T) coefficient method. Then, the Bayesian optimization (BO) is designed and trained to seek the optimal architecture of the deep neural network, such as the number of neurons and hidden layers and their combinations. Once the network is trained, the dispersion curve of VR can be constructed instantaneously without building and solving the secular equation.Results and Discussion:The results show that the DNN-BO achieves a coefficient of determination (R2) and MAE for the training and validation sets of 0.98 and 8.30 and 0.97 and 8.94, respectively, which suggests that the rapid method has satisfactory generalizability and stability. The DNN-BO method accelerates the dispersion curve calculation by at least 400 times, and there is almost no increase in computation expense with an increase in soil layers.

From floor to roof of a batholith: geology and petrography of the north-eastern Serre Batholith (Calabria, southern Italy)

ABSTRACT We present an original geological map of the north-eastern Serre Batholith (southern Italy), together with new field and petrographic data of its main five granitoid units. Our study provides an overall picture of a c. 13-thick floor-to-roof batholith exposure, focusing on the relationships between the magmatic units, which were sequentially emplaced at depth from c. 23 to c. 6 km, in a time interval from c. 297 Ma to c. 292 Ma. Indeed, this composite and zoned batholith, with its crustal scale exposure and large compositional and structural diversity, is a real natural laboratory where to test models of granitoid magma production and batholith construction. The new geological map and related field and petrographic data provide a valuable addition to the existing knowledge of the Serre Batholith and, at the same time, a new starting point for further in-depth multidisciplinary investigations aimed to better understand its architecture and build-up mechanisms.

Crustal Seismic Anisotropy Along the Continental Margin in Western Canada From Receiver Function Analysis

AbstractCoastal British Columbia (BC), Canada, has the highest seismic hazard in the country due to convergent and transpressive deformation at offshore plate boundaries between the Pacific, Juan de Fuca (JdF) and North American (NA) plates. Further landward, the crust of the NA plate is made up of several geologically unique exotic terranes and is unusually thin. Investigating the geophysical features in this area can help us better constrain its tectonic history and the geophysical processes that are currently underway. Here, we conduct an analysis of teleseismic body‐wave scattering data (i.e., receiver functions [RFs]) recorded at stations across western coastal BC including northern Vancouver Island and southeastern Alaska. Using these RFs, we perform a harmonic decomposition with respect to earthquake back‐azimuths to determine the orientation of seismic anisotropy over a series of depth ranges, attributable to either mineral alignment or dipping structures. We find a coherent pattern of margin‐parallel orientations at upper crustal depths that persist onto the mainland at distances 420 km from the margin. Furthermore, dominant receiver function orientations at depth are attributed to dipping faults and interfaces, and fabrics due to lower crustal shearing or inherited from tectonic assembly along the margin. This work supports models for the tectonic assembly of this region that involve a combination of plate subduction and transpressive motion along crustal scale faults that pervade a wide portion of the margin. This work also helps to constrain the current geometries of the subducting Pacific and JdF slabs.

Long-term magmatism in the Ribeira Orogen, Western Gondwana: SHRIMP U-Pb and O isotope data for two Ediacaran-Cambrian magmatic events in the Rio de Janeiro area

In the southeastern Brazilian margin, a Neoproterozoic to the Early Paleozoic long-term magmatic interval is recorded within the Oriental Terrane, part of the Ribeira Orogen, characterizing it as a hot orogenic belt. Metaigneous and igneous rock units of ca. 570–540 Ma and 510–480 Ma are coeval with two main collisional orogenies related to the Western Gondwana assembly: the Ribeira Orogeny (590–565 Ma) and the Búzios Orogeny (540–490 Ma). The Oriental Terrane, interpreted as a continental microplate amalgamated during the convergence between the São Francisco and Congo cratons, comprises Cryogenian to Ediacaran metaigneous rocks and Cambro-Ordovician plutons intrusive in metasedimentary sequences. Zircons from an orthogneiss and a porphyritic leucogneiss sample gave SHRIMP U-Pb crystallization ages of 558 ± 2 Ma and 560 ± 2 Ma, respectively, with high δ18O values from 7 ‰ to 10 ‰. Zircons from a tonalite pluton and tonalitic dikes, gave SHRIMP U-Pb crystallization ages of 483 ± 2 Ma and 492 ± 2 Ma, with lower δ18O values, between 6 ‰ to 8 ‰. These two magmatic pulses, integrated with a compiled database from other plutons in this terrane, corroborate with the hypothesis that they represent crustal scale magmatic provinces here interpreted as pre- and post-collisional in relation with the Búzios Orogeny. Their distribution within the tectonic map of the Oriental terrane reflect a paleogeography linked to their tectonic origin. The 570–540 Ma plutons are aligned ENE -WSW, parallel and close to the suture zone with the Cabo Frio Tectonic Domain, representing a crustal melt related to a NW-subduction. A subsequent magmatic hiatus of 30 ma is interpreted as the collision of the Cabo Frio Tectonic Domain. The post collisional period is marked the intrusion of non– linear 510–480 Ma pluton, widespread in the Oriental terrane.