Global Geodynamics Research Articles

Global mantle perturbations following the onset of modern plate tectonics.

Plate tectonics drives the compositional diversity of Earth's convecting mantle through subduction of lithosphere. In this context, the role of evolving global geodynamics and plate (re)organization on the spatial and temporal distribution of compositional heterogeneities in the convecting mantle is poorly understood. Here, using the geochemical compositions of intracontinental basalts formed over the past billion years, we show that intracontinental basalts with subchondritic initial neodymium-144/neodymium-143 values become common only after 300 million years, broadly coeval with the global appearance of kimberlites with geochemically enriched isotopic signatures. These step changes in the sources of intraplate magmatism stem from a rapid increase in the supply of deeply subducted lithosphere during the protracted formation of Pangea following the widespread onset of "modern" (cold and deep) subduction in the late Neoproterozoic. We argue that the delay (~300 million years) in the appearance of enriched intraplate magmas reflects the time required for the sinking and (re)incorporation of slabs into the sources of mantle-derived magmas.

The geodynamics of plume-influenced mid-ocean ridges: insights from the Foundation Segment of the Pacific-Antarctic Ridge

The intersection of the Foundation Plume with the Pacific-Antarctic Ridge is a key location in global geodynamics where a mantle plume is approached by and interacting with a fast-spreading mid-ocean ridge. Here, we discuss a comprehensive major and trace element and Sr-Nd-Pb isotope dataset of new and existing samples from the young Foundation Seamount Chain (<5 Ma) and adjacent section of the Pacific-Antarctic Ridge. We use the geochemistry of axial, off-axis and intraplate lavas to map the spatial extent of plume dispersal underneath the ridge as well as the internal zonation of the upwelling plume. We show that the unusual length, increased crustal thickness and occurrence of silicic rocks on the axis of the Foundation Segment are the direct result of plume being tapped by the axial melting zone. We demonstrate that the plume is not homogeneous but shows a HIMU-like (high time-integrated 238U/204Pb) OIB (Ocean Island Basalt) component characterized by 206Pb/204Pb of up to 20.5 in its center and a more EM1-like (Enriched Mantle one) OIB component characterized by low U/Pb and 206Pb/204Pb but high Rb/Nb and 87Sr/86Sr towards its edges. Plume entrainment leads to a high magma supply rate that fosters the formation of silicic rocks and triggers the lengthening of the segment over time. However, plume dispersal is not symmetric as the geochemical tracers for the OIB component are extending <100 km northwards but >300 km southwards. We relate this to the current plate tectonic framework in which the obliquity between the migrating ridge and the absolute plate motions induces a sub-axial asthenospheric flow that preferentially channels plume material southwards.

Convergent plate boundary environments for formation of ≥ 3800 Ma mafic–ultramafic assemblages (Isua area, Greenland): Implications for early global geodynamics

In the gneiss terrane on the south side of the Eoarchean Isua supracrustal belt, ultramafic rocks with relict abyssal peridotite mineralogy (Bennett et al., 2002; Friend et al., 2002; Nutman et al., 2007; Rollinson, 2007; van de Löcht et al., 2020), layered gabbros with cumulate ultramafic rocks, basalts and associated siliceous sedimentary rocks were tectonically-imbricated, prior to and during intrusion of ca. 3800 Ma tonalites. Together with ≥ 3800 Ma basalts in the Outer Arc Group of the nearby Isua supracrustal belt, the composition of all these mafic rocks (e.g., Th–Hf–Nb systematics, high Th/Yb, Ba/Nb, Ba/Yb ratios and negative Nb and Ti anomalies) shows affinity with modern suprasubduction rocks whose genesis involved fluid fluxing of the upper mantle. However, the majority of these samples have Ba/Nb and Ba/Yb values less than in modern island arc magmas, but similar to many backarc basin magmas (e.g., Pearce and Stern, 2006). It is unknown whether these ca. 3800 Ma mafic rocks are, (i) arc rocks where the Ba/Nb and Ba/Yb signatures reflect lower surficial Ba in Eoarchean oceanic settings, or (ii) in direct comparison with Phanerozoic suites, these signatures reflect a back-arc setting with interplay between fluid fluxing and decompressional melting. The tectonic intercalation of upper mantle with lower and upper crustal rocks, combined with the fluid-fluxing influences seen in chemistry of all the mafic rocks is best accommodated in a compressional Eoarchean convergent plate boundary setting within a mobile-lid regime. Thus stagnant lid scenarios of crust formation, if operative, must have co-existed or alternated with mobile-lid regimes by 3800 Ma.

Comprehensive examination on codon usage bias pattern of the Bovine Ephemeral fever virus

Bovine Ephemeral Fever Virus (BEFV) is a non-contagious virus that commonly infects cattle and water buffalo, reduces milk productivity, decreases the quality of beef, and causes an adverse economic impact on the global livestock industry. However, the evolution of BEFV is unclear, and uncertainty exists regarding its global geodynamics. Consequently, this study aims to comprehend the pattern of viral evolution and gene expression in the BEFV genes G, M, N, and P, including synonymous codons. Additionally, we performed recombination analyses, which exclusively detected recombination signals in the G- and P-genes. Subsequently, a phylogenetic tree was constructed to validate and support these findings. The codon usage bias results showed that the BEFV-selected genes were influenced by both natural and mutation pressure. Furthermore, nucleotide A is more abundant in all the selected genes. The eNC values, ranging from 42.99 to 47.10, revealed the presence of moderate codon usage bias, where gene P exhibited the highest and gene G had the lowest codon usage bias. The neutrality and PR-2 plots, specified codon usage patterns of the genes, are also being shaped by strong selectional pressure. This comprehensive analysis of BEFV genes (G, M, N, and P) sheds light on the molecular evolutionary patterns, co-adaptation, and different genes expression in diverse regions, facilitating the development of preventative programs and insights into viral pathogenesis and vaccine design. Communicated by Ramaswamy H. Sarma

The Superconducting Gravimeter CD-034 at Moxa Observatory: More than 20 Years of Scientific Experience and a Reanimation

In this work, advanced methods and processing schemes for the analysis of data from a Superconducting Gravimeter (SG) will be introduced and their relevance on acquired data proved. The SG CD-034 was installed on Easter of 1999 in the Geodynamic Observatory Moxa of the Friedrich Schiller University Jena, Germany. Initially, the quality of the recorded data was examined, spectra for the detection of the parasitic modes were calculated and the calibration values for the two sensors were determined. Ever since very high-quality gravity data of this SG and most of the other worldwide SGs were made available through the storage archive of the Global Geodynamics Project (GGP later changed to IGETS, International Geodynamics and Earth Tide Service) for global scientific investigations at that time. SG’s such as the one in Moxa (Germany) still deliver significant scientific value for global gravitational field studies as well as for regional/local studies which will be shortly reviewed. Examples are the detection of polar motion, the influence of continental water loading in general and in particular river basin loads, the gravimetric effect of North Sea storm surges and the study of hydro-gravimetric signals, which could be compared with satellite observations and global hydraulic models. The long-term, low-noise operation of complex SG’s requires some effort on maintenance. In order to evaluate the correct operation of the SG, new data processing steps were introduced to assist in the analysis of the data in case of issues with the instrumentation. For example, in 2012/2013 and 2020/2021 severe interference in the gravimeter electronics in Moxa led to a significant loss of data. In both cases, however, the cause could be determined, and the corresponding electronic components renewed. Since July 2021, the SG in Moxa registers again with high data quality comparable or slightly better than before the incident. Initial tests and tidal analyses confirm the validity of the old calibration factors, and the authors now look forward to the re-established long-term recording with excitement and confidence.

Tectonic setting of metamorphism and exhumation of eclogite-facies rocks in the South Beishan orogen, northwestern China

AbstractHigh-pressure metamorphic rocks occur as distinct belts along subduction zones and collisional orogens or as isolated blocks within orogens or mélanges and represent continental materials that were subducted to deep depths and subsequently exhumed to the shallow crust. Understanding the burial and exhumation processes and the sizes and shapes of the high-pressure blocks is important for providing insight into global geodynamics and plate tectonic processes. The South Beishan orogen of northwestern China is notable for the exposure of early Paleozoic high-pressure (HP), eclogite-facies metamorphic rocks, yet the tectonism associated with the HP metamorphism and mechanism of exhumation are poorly understood despite being key to understanding the tectonic evolution of the larger Central Asian Orogenic System. To address this issue, we examined the geometries, kinematics, and overprinting relationships of structures and determined the temperatures and timings of deformation and metamorphism of the HP rocks of the South Beishan orogen. Geochronological results show that the South Beishan orogen contains ca. 1.55–1.35 Ga basement metamorphic rocks and ca. 970–866 Ma granitoids generated during a regional tectono-magmatic event. Ca. 500–450 Ma crustal thickening and HP metamorphism may have been related to regional contraction in the South Beishan orogen. Ca. 900–800 Ma protoliths experienced eclogite-facies metamorphism (~1.2–2.1 GPa and ~700–800 °C) in thickened lower crust. These HP rocks were subsequently exhumed after ca. 450 Ma to mid-crustal depths in the footwall of a regional detachment fault during southeast-northwest–oriented crustal extension, possibly as the result of rollback of a subducted oceanic slab. Prior to ca. 438 Ma, north-south–oriented contraction resulted in isoclinal folding of the detachment fault and HP rocks. Following this contractional phase in the middle Mesozoic, the South Beishan orogen experienced thrusting interpreted to be the response to the closure of the Tethyan and Paleo-Asian Ocean domains. This contractional phase was followed by late Mesozoic extension and subsequent surface erosion that controlled exhumation of the HP rocks.

Desilicification Rims of Zircon Xenocrysts Record the Timing of Kimberlite Emplacement

AbstractKimberlites are deep and relatively rare mantle‐derived igneous rocks, which provide key insights into the lithospheric mantle composition beneath continents (i.e., from the entrained mantle materials) at the time of magma emplacement. In addition, the temporal distribution of kimberlites throughout Earth's history may reveal secular changes in global geodynamics. However, alteration and abundant exotic material (i.e., xenoliths and xenocrysts) of different ages, both typical of kimberlites, often hinder reliable dating of these rocks. In this study, we investigated zircon xenocrysts and their replacement by desilicification reaction rims (DSRs), which consist mainly of baddeleyite, in the metamorphosed Kimozero kimberlites from the Karelian Craton, Russia. Similar DSRs have been identified previously on mantle‐derived zircons in numerous kimberlites, but the precise timing and location of formation of these rims remain unclear. U–Pb isotope dating of mantle zircon xenocrysts in the Kimozero kimberlites indicates that these rocks were emplaced at ca. 1.97 Ga. DSRs with similar textures and baddeleyite chemistry occur on both the 1.97 Ga mantle and 2.70–2.38 Ga crustal zircon xenocrysts in these kimberlites. Pb–Pb isotope dating of the DSRs yields the same age (i.e., ca. 1.97 Ga) as the mantle zircons. The data constrain the formation of DSRs to crustal levels and, more importantly, indicate that these rims record the timing of kimberlite emplacement. Therefore, we suggest that such baddeleyite reaction rims are a new and promising material for kimberlite geochronology.

Absolute Paleolatitude of Northern Zealandia From the Middle Eocene to the Early Miocene

AbstractThe absolute position during the Cenozoic of northern Zealandia, a continent that lies more than 90% submerged in the southwest Pacific Ocean, is inferred from global plate motion models, because local paleomagnetic constraints are virtually absent. We present new paleolatitude constraints using paleomagnetic data from International Ocean Discovery Program Site U1507 on northern Zealandia and Site U1511 drilled in the adjacent Tasman Sea Basin. After correcting for inclination shallowing, five paleolatitude estimates provide a trajectory of northern Zealandia past position from the middle Eocene to the early Miocene, spanning geomagnetic polarity chrons C21n to C5Er (∼48–18 Ma). The paleolatitude estimates support previous works on global absolute plate motion where northern Zealandia migrated 6° northward between the early Oligocene and early Miocene, but with lower absolute paleolatitudes, particularly in the Bartonian and Priabonian (C18n–C13r). True polar wander (solid Earth rotation with respect to the spin axis), which only can be resolved using paleomagnetic data, may explain the discrepancy. This new paleomagnetic information anchors past latitudes of Zealandia to Earth's spin axis, with implications not only for global geodynamics, but also for addressing paleoceanographic and paleoclimate problems, which generally require precise paleolatitude placement of proxy data.

On computing viscoelastic Love numbers for general planetary models: the ALMA3 code

SUMMARY The computation of the Love numbers (LNs) for a spherically symmetric self-gravitating viscoelastic Earth is a classical problem in global geodynamics. Here we revisit the problem of the numerical evaluation of loading and tidal LNs in the static limit for an incompressible planetary body, adopting a Laplace inversion scheme based upon the Post-Widder formula as an alternative to the traditional viscoelastic normal modes method. We also consider, within the same framework, complex-valued, frequency-dependent LNs that describe the response to a periodic forcing, which are paramount in the study of the tidal deformation of planets. Furthermore, we numerically obtain the time-derivatives of LNs, suitable for modelling geodetic signals in response to surface loads variations. A number of examples are shown, in which time and frequency-dependent LNs are evaluated for the Earth and planets adopting realistic rheological profiles. The numerical solution scheme is implemented in ALMA3 (the plAnetary Love nuMbers cAlculator, version 3), an upgraded open-source Fortran 90 program that computes the LNs for radially layered planetary bodies with a wide range of rheologies, including transient laws like Andrade or Burgers.

Structure of the Earth's crust and upper mantle and problems of global geodynamics

This work presents findings of deep geological and geophysical studies, which have revealed new and often unexpected features of the structure of the Earth's crust and the upper mantle of the continents and oceans, not easily explainable by current geodynamic concepts. For example, deepwater drilling has shown the old age of the oceanic crust and its inhomogeneous composition. Seismic studies of long-range profiles with nuclear explosions and laboratory studies of mantle xenoliths at high pressures and temperatures have revealed that the advection of energy-intensive and geochemically active deep fluids is central to the formation of the upper mantle structure. A comprehensive analysis of these data provides a geodynamic model for the natural transformation of the upper mantle substance influenced by deep fluids and the formation of continents and oceans as a result of unevenly distributed degassing of the Earth. The global variability of the deep fluid advection intensity led to formation of two hemispheres: the Indo Atlantic Hemisphere where almost all continents are located and the Pacific Hemisphere that consists of oceans. This formation of the hemispheres can be explained by the specific features of the Earth's rotation in the Earth-Moon-Sun system. This work provides a review and analysis of these data and then offers a possible explanation based on an integrated geodynamic model.

Fact-condition statements and super relation extraction for geothermic knowledge graphs construction

Researchers utilize information from the geoscience literature to deduce the regional or global geological evolution. Traditionally this process has relied on the labor of researchers. As the number of papers continues to increase, acquiring domain-specific knowledge becomes a heavy burden. Knowledge Graph (KG) is proposed as a new knowledge representation technology to change this situation. However, the super relation is not considered in the previous KG, which bridges the geological phenomenon (fact) and its precondition (condition). For instance, in the statement (“the late Archean was a crucial transition period in the history of global geodynamics”), the condition statement (“crucial transition for global geodynamics”) works as the complementary fact statement (“the late Archean was a crucial transition period”), which defines the scale of crucial transition accurately in the late Archean. In this study, fact-condition statement extraction is introduced to construct a geological knowledge graph. A rule-based multi-input multi-output model (R-MIMO) is proposed for information extraction. In the R-MIMO, fact-condition statements and their super relation are considered and extracted for the first time. To verify its performances, a GeothCF dataset with 1455 fact tuples and 789 condition tuples is constructed. In experiments, the R-MIMO model achieves the best performance by using BERT as encoder and LSTM-d as decoder, achieving F1 80.24% in tuple extraction and F1 70.03% in tag prediction task. Furthermore, the geothermic KG with super relation is automatically constructed for the first time by trained R-MIMO, which can provide structured data for further geothermic research.

SHARED RESEARCH FACILITIES "PETROPHYSICS, GEOMECHANICS AND PALEOMAGNETISM" OF THE SCHMIDT INSTITUTE OF PHYSICS OF THE EARTH RAS

Shared Research Facilities "Petrophysics, Geomechanics and Paleomagnetism", created based on actively developing divisions of the Schmidt Institute of Physics of the Earth RAS, has advanced research equipment and unique methods, which determines the possibility of solving a wide range of applied and fundamental problems in Earth sciences. The main activities of the Shared Research Facilities "Petrophysics, Geomechanics and Paleomagnetism" are in the field of geomechanical modelling and petrophysics, paleomagnetism, rock magnetism and applied analytical studies of the composition and structure of minerals, rocks and materials. The research results carried out at the Center are used to solve actual problems of prospecting and exploration of oil and gas basins, mineral deposits, tectonics, global geodynamics and tectonophysics. The article presents the scientific equipment of the Center and its key characteristics, describes the methodological and analytical capabilities, and the most significant results of research groups.

Discussion of: “A geological map of the Scotia Sea area constrained by bathymetry, geological data, geophysical data and seismic tomography models from the deep mantle”, by A. Beniest and W.P. Schellart (https://doi.org/10.1016/j.earscirevv.2020.103391)

Global geophysical networks provide powerful databases to infer globally coherent signals, and array processing techniques are useful for inferring them. In this study, we comprehensively analyze seven spherical harmonic-based array processing techniques: spherical harmonic stacking (SHS) as well as its gridded form (SHS_GT) and grid-interval weighted forms (SHS_GK1, SHS_GK2); matrix SHS (MSHS); multi-station experiment (MSE); and optimal sequence estimation (OSE). We first use more specific synthetic tests to evaluate the pros and cons of these techniques, and estimate bias in solutions caused by the station distributions. These methods are applied to four global observation networks, the Global Geodynamics Project (GGP) Network, the Global Seismographic Network (GSN), the Global Geomagnetic (GGM) Network and the Global GNSS Network. For the first time, we restored a much cleaner sequence for one singlet of the 0S2 mode based on the GGP network, and restored similar result for one singlet of the 3S1 mode based on the GSN network. We further isolate different Ylm-related tidal signals from the GGM network for the first time. Moreover, based on global GNSS observations, we estimate the Love number h21 = 0.6243(±7e−4) − 0.01(±6e−3)i at the Chandler Wobble (CW) frequency with OSE/MSHS (The accuracy of the estimate is an order of magnitude higher than the previous results), and further obtain the corresponding lower-mantle anelasticity (fr(ω) = −29.5 ± 0.9, fi(ω) = 12.0 ± 7.2). Our findings confirm that OSE and MSHS methods can more accurately obtain the complex amplitude of any Ylm-related signal, which is not possible with other methods (and we do obtain more precise results than previous studies upon using them); besides, we also confirm that OSE and MSHS methods can greatly reduce the interference of other signals to the target signals. Hence, we believe the results obtained from the OSE/MSHS will helpful for obtaining reasonable geophysical explanations.

Array processing techniques for analyzing global geophysical observation networks: Comprehensive analysis and typical applications

Global geophysical networks provide powerful databases to infer globally coherent signals, and array processing techniques are useful for inferring them. In this study, we comprehensively analyze seven spherical harmonic-based array processing techniques: spherical harmonic stacking (SHS) as well as its gridded form (SHS_GT) and grid-interval weighted forms (SHS_GK1, SHS_GK2); matrix SHS (MSHS); multi-station experiment (MSE); and optimal sequence estimation (OSE). We first use more specific synthetic tests to evaluate the pros and cons of these techniques, and estimate bias in solutions caused by the station distributions. These methods are applied to four global observation networks, the Global Geodynamics Project (GGP) Network, the Global Seismographic Network (GSN), the Global Geomagnetic (GGM) Network and the Global GNSS Network. For the first time, we restored a much cleaner sequence for one singlet of the 0S2 mode based on the GGP network, and restored similar result for one singlet of the 3S1 mode based on the GSN network. We further isolate different Ylm-related tidal signals from the GGM network for the first time. Moreover, based on global GNSS observations, we estimate the Love number h21 = 0.6243(±7e−4) − 0.01(±6e−3)i at the Chandler Wobble (CW) frequency with OSE/MSHS (The accuracy of the estimate is an order of magnitude higher than the previous results), and further obtain the corresponding lower-mantle anelasticity (fr(ω) = −29.5 ± 0.9, fi(ω) = 12.0 ± 7.2). Our findings confirm that OSE and MSHS methods can more accurately obtain the complex amplitude of any Ylm-related signal, which is not possible with other methods (and we do obtain more precise results than previous studies upon using them); besides, we also confirm that OSE and MSHS methods can greatly reduce the interference of other signals to the target signals. Hence, we believe the results obtained from the OSE/MSHS will helpful for obtaining reasonable geophysical explanations.

Mantle metasomatic influence on water contents in continental lithosphere: New constraints from garnet pyroxenite xenoliths (France & Cameroon volcanic provinces)

Quantifying water contents in the lithospheric mantle is key to our understanding of global geodynamics, mantle composition, and related physical properties. Most mantle lithologies (peridotite) contain little water (~50 ppm), but petrological heterogeneities such as pyroxenites are more hydrous (~300 ppm) relative to the mantle rocks. Pyroxenites also melt at lower temperatures than peridotites and are thus important to magma genesis. Thus, quantifying pyroxenite water contents provides new information on the distribution of water in the mantle. Here, we present phase-specific FTIR measurements of the water contents in pyroxenite mantle xenoliths from two continental lithospheric domains that experienced intense metasomatism: the French Massif Central (FMC, France) and the Adamawa Volcanic Plateau (AVP, Cameroon). The AVP garnet pyroxenites are more hydrated ([H2O]Clinopyroxene = 386–685 ppm; [H2O]Orthopyroxene = 124–155 ppm; [H2O]Garnet < 0.5 ppm) than FMC ones ([H2O]Clinopyroxene = 112–465 ppm; [H2O]Orthopyroxene = 61–104 ppm; [H2O]Garnet < 0.5 ppm). These water concentrations are homogenous at the grain and correlate with equilibrated major element concentrations, indicating that they are representative of lithospheric water, although the FMC pyroxenites were dehydrated during metasomatism by a carbonatitic fluid (based on the correlation between LaN/SmN and Ti/Eu ratios); the water contents of AVP pyroxenites were likely not affected by metasomatism. FMC pyroxenites show peculiar FTIR spectra that may reflect the preferential dehydration of specific sites in the pyroxene structure. In both regions, metasomatism modified the light rare Earth element contents (e.g., Ce) of the pyroxenites, resulting in highly variable H2O/Ce ratios. Therefore, we conclude that the utility of the H2O/Ce ratio to identify the involvement of pyroxenites in magmas genesis is limited.

Two contrasting Neoarchean metavolcanic rock suites in eastern Hebei and their geodynamic implications for the northern North China Craton

The late Archean was a crucial transition period for global geodynamics. Metavolcanic rocks, as one of the predominant components of Archean cratons, are key to deciphering the geodynamic regimes of early Earth crust-mantle interactions. Neoarchean supracrustal rocks in the central K-rich granitoid belt of the eastern Hebei-western Liaoning Province, northern North China Craton, may be subdivided into a lower metavolcanic sequence and an upper metavolcanic-sedimentary sequence by an angular unconformity.The lower sequence metavolcanic rocks consist primarily of 2558–2536 Ma amphibolite to granulite facies metamorphosed basaltic-andesitic rocks and high- and low-MgO dacitic rocks. The basaltic-andesitic samples display calc-alkaline affinity, high MgO contents, low (Nb/La)N and moderate (Hf/Sm)N values and were most likely derived from partial melting of mantle materials metasomatized by subducted slab- and sediment-derived melts and fluids. The high-MgO dacitic samples are characterized by strongly fractionated rare earth element (REE) patterns and relatively high MgO contents, which resemble typical high-SiO2 adakites. The magmatic precursors originated from subducted oceanic slabs and sediments. The samples with low-MgO dacite compositions have the highest SiO2 and lowest MgO contents, and the magmatic precursors were generated by partial melting of metamorphosed high-K mafic rocks and sedimentary rocks in the lower crustal level.The upper sequence metavolcanic rocks chiefly comprise greenschist to amphibolite facies metamorphosed tholeiites, tholeiitic to transitional basalts, calc-alkaline basalts, andesites, and dacites with magmatic crystallization ages of 2522–2511 Ma. The tholeiite samples are characterized by left-declined REE patterns and high (Nb/La)N and chondrite-like (Hf/Sm)N values resembling typical N-MORBs and were formed by partial melting of a depleted mantle source. The tholeiitic to transitional basalt samples resemble back-arc basin basalts and display weakly right-declined REE patterns, mildly negative Nb and Ti anomalies, and low (Nb/La)N and chondrite-like (Hf/Sm)N values. Their magmatic precursors were derived from partial melting of a slightly fluid- and melt-metasomatized depleted mantle source. The calc-alkaline basalt and andesite samples exhibit fractionated REE patterns and apparent Nb, Ta, and Ti depletions and originated from a strongly fluid- and melt-metasomatized depleted mantle source. The dacite samples show the highest SiO2 and lowest MgO contents in the upper sequence metavolcanic rocks and were derived from partial melting of lower crustal metasedimentary rocks. Systematic petrological, geochemical, and petrogenetic studies indicate that a mature continental arc setting and a back-arc basin setting may account for the formation of the lower and upper sequence metavolcanic rocks, respectively. Three stages of deformation can be identified in the K-rich granitoid belt, including 2525–2505 Ma extensional deformation in a back-arc basin setting (D1), ca. 2500 Ma top-to-the-NEE-SSE thrusting in a collisional orogenic setting (D2) and 2500–1800 Ma nearly N-S-oriented strike-slip shearing (D3), with D2 deformation accompanied by the synchronous occurrence of paired high-pressure and high-temperature medium-pressure granulite facies metamorphism. Integrating the petrogenesis, spatiotemporal relationships, deformational and metamorphic features of the Archean lithological assemblages, we consider that the late Neoarchean K-rich granitoid belt was formed in an active continental margin setting with a southeastward subduction polarity and experienced multi-stage trench retreats and a subsequent collision between the arc and continent.

Time Between 3 and 2 Ga: Transitional Events in the Earth’s History

Abstract —The time span between 3 and 2 Ga in the geologic history encompassed a number of key events on the cooling Earth. The cooling interrupted heat transfer within and across the mantle, which caused changes in Earth’s major spheres and in the mechanisms of their interaction. The great thermal divergence at 2.5 Ga and differentiation into the depleted upper asthenospheric and primitive lower mantle affected the compositions of oceanic basalts. The lower mantle cooling recorded by a systematic decrease in the temperature of komatiite magma generation at the respective depths began at 2.5 Ga and was accompanied by increasing abundance of arc basalts and by changes in the behavior of the Sr, Nd, and O isotope systems. It was the time when the continental lithosphere consisting of subcontinental lithospheric mantle and crust began its rapid growth, while the crust became enriched in felsic material with high contents of lithophile elements. Magmatism of the 3–2 Ga time span acquired more diverse major-element chemistry, with calc-alkaline and alkaline lithologies like carbonatite and kimberlite. The dramatic changes were driven by subduction processes, whereby the crust became recycled in the mantle and the double layer (D”) formed at the core–mantle boundary. The events of the 3–2 Ga interval created prerequisites for redox changes on the surface and release of free oxygen into the atmosphere. In terms of global geodynamics, it was transition from stagnantlid tectonics to plate tectonic regime, which approached the present-day style about 2.0–1.8 Ga.

The Caribbean and Farallon Plates Connected: Constraints From Stratigraphy and Paleomagnetism of the Nicoya Peninsula, Costa Rica

AbstractPlate kinematic reconstructions play an essential role in our understanding of global geodynamics, but become increasingly difficult to constrain back in geological time due to the subduction of oceanic lithosphere. Here, we attempt to kinematically reconstruct the Cretaceous and older plate tectonic history of the Caribbean Plate within the Mesozoic Panthalassa (paleo‐Pacific) Ocean. To this end, we present new paleomagnetic data from Jurassic and Cretaceous oceanic sedimentary and volcanic Large Igneous Province‐related rocks of the Nicoya Peninsula and Murciélago Islands of northwestern Costa Rica. We use these data, in combination with constraints from marine magnetic anomalies to infer the age of the lithospheric basement, seismic tomography to locate deep‐mantle plume generation zones, and general kinematic feasibility, to test different reconstruction scenarios connecting the Caribbean Plate to the Farallon Plate as restored from Pacific spreading records. Our resulting reconstruction implies that the western Caribbean subduction zone initiated around 100 Ma, in an intraoceanic setting, breaking up oceanic lithosphere of at least 70 Myr old.

Deep hydrous mantle reservoir provides evidence for crustal recycling before 3.3 billion years ago.

Water strongly influences the physical properties of the mantle and enhances its ability to melt or convect. Its presence can also be used to trace recycling of surface reservoirs down to the deep mantle1, which makes knowledge of the water content in the Earth's interior and its evolution crucial for understanding global geodynamics. Komatiites (MgO-rich ultramafic magmas) result from a high degree of mantle melting at high pressures2 and thus are excellent probes of the chemical composition and water contents of the deep mantle. An excess of water over elements that show similar geochemical behaviour during mantle melting (for example, cerium) was recently found in melt inclusions in the most magnesium-rich olivine in 2.7-billion-year-old komatiites from Canada3 and Zimbabwe4. These data were taken as evidence for a deep hydrated mantle reservoir, probably the transition zone, in the Neoarchaean era(2.8 to 2.5 billion years ago). Here we confirm the mantle source of this water by measuring deuterium-to-hydrogen ratios in these melt inclusions and present similar data for 3.3-billion-year-old komatiites from the Barberton greenstone belt. From the hydrogen isotope ratios, we show that the mantle sources of these melts contained excess water, which implies that a deep hydrous mantle reservoir has been present in the Earth's interior since at least the Palaeoarchaean era(3.6 to 3.2 billion years ago). The reconstructed initial hydrogen isotope composition of komatiites is more depleted in deuterium than surface reservoirs or typical mantle but resembles that of oceanic crust that was initially altered by seawater andthen dehydrated during subduction. Together with an excess of chlorine and depletion of lead in the mantle sources of komatiites, these results indicate that seawater-alteredlithosphere recycling into the deep mantle, arguably by subduction, started before 3.3 billion years ago.

Problems of global geodynamics

Global geodynamics is determined by thermal convection in the mantle which manifests itself on the surface by movements, relief, heat flow, and volcanism. Thermal convection in the Earth is complicated by the fact that the lithosphere is broken into rigid plates, the crust is broken into six separate floating continents and a number of islands, on the mantle bottom there are two giant piles of heavy material, at high convection intensity the ascending convective flows acquire a plume shape, and phase transformations take place in the mantle. The impacts of many factors on the mantle structure have been thoroughly studied and fairly well understood. It is pertinent to reconcile the new data on phase transformations at depths of 650 to 700 km with the seismic data on the positions of these boundaries. The ultimate problem of global geodynamics has not yet been solved; the three-dimensional structure of the whole-mantle flows, consistent with the observations in geophysics, geochemistry, geology, and numerical modeling, is not known even on a semischematic level.