Decratonization by rifting enables orogenic reworking and transcurrent dispersal of old terranes in NE Brazil

Page created by Carmen Martin
 
CONTINUE READING
Decratonization by rifting enables orogenic reworking and transcurrent dispersal of old terranes in NE Brazil
www.nature.com/scientificreports

                OPEN             Decratonization by rifting
                                 enables orogenic reworking
                                 and transcurrent dispersal of old
                                 terranes in NE Brazil
                                 Carlos E. Ganade1*, Roberto F. Weinberg2, Fabricio A. Caxito3, Leonardo B. L. Lopes1,
                                 Lucas R. Tesser4 & Iago S. Costa1
                                 Dispersion and deformation of cratonic fragments within orogens require weakening of the craton
                                 margins in a process of decratonization. The orogenic Borborema Province, in NE Brazil, is one
                                 of several Brasiliano/Pan-African late Neoproterozoic orogens that led to the amalgamation of
                                 Gondwana. A common feature of these orogens is that a period of extension and opening of narrow
                                 oceans preceded inversion and collision. For the case of the Borborema Province, the São Francisco
                                 Craton was pulled away from its other half, the Benino-Nigerian Shield, during an intermittent
                                 extension event between 1.0–0.92 and 0.9–0.82 Ga. This was followed by inversion of an embryonic
                                 and confined oceanic basin at ca. 0.60 Ga and transpressional orogeny from ca. 0.59 Ga onwards.
                                 Here we investigate the boundary region between the north São Francisco Craton and the Borborema
                                 Province and demonstrate how cratonic blocks became physically involved in the orogeny. We
                                 combine these results with a wide compilation of U–Pb and Nd-isotopic model ages to show that
                                 the Borborema Province consists of up to 65% of strongly sheared ancient rocks affiliated with the
                                 São Francisco/Benino-Nigerian Craton, separated by major transcurrent shear zones, with only ≈
                                 15% addition of juvenile material during the Neoproterozoic orogeny. This evolution is repeated
                                 across a number of Brasiliano/Pan-African orogens, with significant local variations, and indicate
                                 that extension weakened cratonic regions in a process of decratonization that prepared them for
                                 involvement in the orogenies, that led to the amalgamation of Gondwana.

                                  Cratons are old, stable continental regions with thick buoyant keels that resist deformation. Their thick and
                                  cold refractory lithospheric keel is taken to be responsible for their integrity, protecting them from tectonic
                                  ­reworking1. However, during continental collisions, colliding cratonic blocks can be partially reworked to gen-
                                   erate continental tracts that are no longer cratons but that are not typical orogens e­ ither2. Thus, in order to
                                   understand the origin of these old continental crustal tracts within orogenies, it is necessary to address the
                                   specific mechanisms for both the weakening of the craton lithosphere (c.f. decratonization) and their subsequent
                                   reworking. First order controls in the reworking of these strong lithospheres are temporal and spatial variations
                                   of their thermal state and pre-existing mechanical anisotropies or major compositional ­boundaries3–8. The North
                                   China Craton is a classic example of decratonization where lithospheric thinning, asthenospheric decompres-
                                   sion and magmatism changed the thermal and chemical state of the cratonic lithosphere resulting in wholesale
                                 ­decratonization9–12.
                                       Here we describe the process of decratonization of the São Francisco Craton (SFC) during the Neoprote-
                                   rozoic in northeast Brazil. This craton is bound to the north by the late Neoproterozoic orogenic Borborema
                                   Province (BP). In the province, pre-orogenic extension started as early as 1.0 Ga and resulted in the develop-
                                   ment of intracontinental extensional b  ­ asins13, passive margin basins flanking the c­ raton14–16, and separation of
                                   the SFC from the African Benino-Nigerian ­Shield17,18. Remarkably, this orogenic province includes a number
                                   of reworked and deformed Archean-Paleoproterozoic terranes among the Neoproterozoic sedimentary basins
                                   and granitic ­intrusions19. These old terranes are relatively small (10,000 to 40,000 k­ m2), bounded by large-scale
                                   continental shear ­zones19,20 (Fig. 1A and B), and their affinity with the SFC have been previously ­proposed21.

                                 1
                                  Geological Survey of Brazil/CEDES, Rio de Janeiro, Brazil. 2School of Earth, Atmosphere and Environment, Monash
                                 University, Clayton, Australia. 3Universidade Federal de Minas Gerais, Belo Horizonte, Brazil. 4Universidade de São
                                 Paulo, São Paulo, Brazil. *email: carlos.ganade@cprm.gov.br

Scientific Reports |   (2021) 11:5719                 | https://doi.org/10.1038/s41598-021-84703-x                                                     1

                                                                                                                                                  Vol.:(0123456789)
Decratonization by rifting enables orogenic reworking and transcurrent dispersal of old terranes in NE Brazil
www.nature.com/scientificreports/

                                            Figure 1.  Age of lithological units and distribution of model ages in the Borborema Province. (A) Interpreted
                                            geological-age map of the Borborema Province. Labels correspond to names of shear zones. Geological base
                                            maps are from the Geological Survey of Brazil and freely available at http://geosg​b.cprm.gov.br. (B) Sm–Nd T ­ DM
                                            age map of the Borborema Province and northern São Francisco Craton with magnetic first derivative image as
                                            background (see Fig. S1 for data point location). Note the Entremontes promontory (labelled E), a part of the
                                            SFC separating two fold-and-thrust belts of the Southern Borborema Orogen. SP Senador Pompeo shear zone,
                                            PT Patos shear zone, PE Pernambuco shear zone, BV Boa Vista shear zone, SQ Santa Quitéria arc, TP Tróia-
                                            Pedra Branca terrane, NBB Northern Borborema block, SC São José do Campestre terrane, AM Alto Moxotó
                                            terrane, E Entremontes block, SBO Southern Borborema Orogen, RP Riacho do Pontal fold-and-thrust belt,
                                            SER Sergipano fold-and-thrust belt, C Cachoeirinha Belt, S Seridó Belt. (C) P-wave tomography model for the
                                            BP at 100 km ­depth35 showing a markedly high velocity under the Alto Moxotó terrane. (D) Spatial distribution
                                            of zircon U–Pb crystallization ages of igneous and metaigneous rocks of the Borborema Province. (E) Spatial
                                            distribution of εNd(t) values calculated for the time of crystallization shown in (D). The Sm–Nd T
                                                                                                                                             ­ DM ages, U–Pb
                                            ages and εNd(t) values were gridded in the ArcGIS software using Inverse-Distance-Weighted Interpolation
                                            (IDW).

                                            Although no quantitative crustal growth curves of the orogenic Borborema Province have been reported, pre-
                                            sent models range from dominantly intracontinental r­ eworking21,22, to arc ­accretion23,24,37 and docking followed
                                            by reworking of unrelated t­ erranes25, or even to questionable autochthonous accretion surrounding old small
                                            reworked Archean ­crust26.
                                               The orogenic Borborema Province (BP) comprises two separate and interacting collisions: one on the west
                                            side, part of the vast West Gondwana O­ rogen27, and one in the south, against the SFC, and termed here the South-
                                            ern Borborema Orogen. These two collisions were partly contemporaneous and interacted involving the entire
                                            BP. In this paper we will use the term BP to refer to the whole orogenic area. The current tectonic configuration

          Scientific Reports |   (2021) 11:5719 |                https://doi.org/10.1038/s41598-021-84703-x                                                  2

Vol:.(1234567890)
Decratonization by rifting enables orogenic reworking and transcurrent dispersal of old terranes in NE Brazil
www.nature.com/scientificreports/

                                  of the province is controlled by these ­collisions19. Their tectonic nature is indicated by the record of UHP meta-
                                  morphism and ­eclogites27,28, marginal ­ophiolites29,30, and arc ­sequences16,24,37, that together suggest the West
                                  Gondwana Orogen resulted from oceanic subduction followed by collision between the Benino-Nigerian Shield
                                  and the West African C ­ raton19,27, while the Southern Borborema Orogen resulted from collision of the BP against
                                  the SFC. The two different orogens interacted in a complex collisional zone, giving rise to an intricate network
                                  of continental shear zones that controlled d   ­ eformation19. In this context, the large Transbrasiliano-Kandi shear
                                  zone, that developed as a result of the collision along the West Gondwana Orogen, acted as a dextral transfer
                                  zone leading to BP collision against the SFC and the development of the Southern Borborema ­Orogen17.
                                      In this paper, we compile large volumes of whole-rock Nd isotope data, zircon U–Pb geochronology data, as
                                  well as geological and geophysical data from the BP and northern part of the SFC and the Benino-Nigerian Shield.
                                  We use the data to first quantify the Neoproterozoic continental growth of the BP and establish the dominance
                                  of recycled cratonic material, and then to investigate how the decratonization allowed former cratonic Archean-
                                  Paleoproterozoic terranes to be entrained in the BP during the Neoproterozoic Orogeny.

                                  Geological setting of the Borborema Province
                                    The distribution of rocks in the triangular wedge-shaped orogenic BP (Fig. 1) was ultimately controlled by a
                                    set of Neoproterozoic strike-slip shear z­ ones19,20. These shear zones bound the north, central and south sub-
                                    provinces and within them several other informal d      ­ omains32. Shearing was accompanied by granitoids with
                                   intrusions concentrated at 0.59–0.56 ­Ga19,22. Archean and Paleoproterozoic rocks are found as basement inliers,
                                    often referred as blocks and/or terranes, all over the BP, and always bounded by Neoproterozoic shear zones
                                    (Fig. 1). Early Neoproterozoic extensional events starting as early as 1.0 Ga are marked by granitoid intrusions,
                                    bimodal volcanism and deposition of immature terrigenous and minor carbonatic sediments during the so-called
                                    Cariris Velhos ­event13. This event was initially defined as ­orogenic33, however no associated Tonian deforma-
                                   tion and metamorphism have been reported, which led to the proposition that this was an extensional ­event21.
                                    In the southern BP, this event culminated in the development of a passive margin sequence along the edge of
                                    the SFC associated with 0.9–0.82 Ga mafic–ultramafic intrusions, proximal margin-type ophiolites and A-type
                                  ­granites14,16. Early Neoproterozoic passive margin and intracontinental sedimentation was overlain by orogenic
                                    sedimentary successions as young as 0.58 ­Ga16 contemporaneous with other basins within the central and north
                                   ­BP34. The lithological boundary of the SFC and BP is defined in many geological maps (Fig. 1A), however seismic
                                   data indicate that the lithospheric cratonic signature at 100 km depth extends further north ending within the
                                   central portion of ­BP35,36 (Fig. 1C).
                                        During the Late Neoproterozoic, from 0.65 to 0.61 Ga, predating continental collisional events, continental
                                    magmatic arcs developed in the w   ­ est24 and s­ outh37 parts of the province, closing intervening oceanic basins.
                                    During the collisions that ensued, the province was squeezed between two impinging continents, one continent
                                    coming in from the west and the other, the São Francisco Craton, coming in from the ­south19. As a consequence,
                                    the BP tectonically “escaped” obliquely with the development of a throughgoing network of transcurrent shear
                                  ­zones19,20 that both shear and bound Archean-Paleoproterozoic ­terranes38–40, such as the Alto Moxotó t­ errane25
                                    (Fig. 1). Just how these old terranes became involved in the orogen remains unclear.

                                  Results
                                   U–Pb ages, Nd isotopic spatial distribution and crustal growth. The U–Pb ages of igneous rocks in
                                   the northern SFC and Benino-Nigerian Shield show a similar pattern to those of the ancient basement terranes
                                   in the BP, recording similar episodes of magma production during the Archean and most intensely during the
                                   Paleoproterozoic (Fig. 2A, Supplementary Data 1). In the Neoproterozoic, U–Pb ages of igneous rocks from the
                                   BP show two intervals of magma production at 1.0–0.92 Ga and 0.67–0.52 Ga, with a few ages < 0.50 Ga (Fig. 2B).
                                   Both intervals are characterized by production of new juvenile crust combined with recycling of old pre-existing
                                   crust. The spatial distribution of the interpolated U–Pb ages and εNd(t) values (Fig. 1D and E), excluding the areas
                                   covered by sedimentary and metasedimentary rocks, was used to estimate the crustal growth during the evolu-
                                   tion of the Borborema Province (see “Methods” section, Supplementary Fig. S1).
                                       The gridded values resulted in a map with 46,426 pixels with a pixel size of 3.2 × 3.2 km (pixel area = 10.24
                                  ­km2). These pixels where clipped to the crystalline basement area of the Borborema Province, excluding the
                                   Neoproterozoic metasedimentary belts, Phanerozoic basins and more recent cover, resulting in a map with
                                   17,235 pixels with a total area of 176,486 ­km2. Pixel values of Nd model ages (Fig. 1B) were binned into 100 Myr
                                   interval, from 0.5 to 3.4 Ga, and so were the U–Pb crystallization ages from individual rock samples. Positive
                                   εNd(t) values were used to distinguish juvenile addition of new crust to the Borborema Province through time.
                                       The results based on the exposed surface area, incorporating both zircon U–Pb and Nd isotopes, indicate that
                                   65–60% of the province was already formed by the end of the Paleoproterozoic, with a rapid growth rate between
                                   2.3 to 2.0 Ga, when 55–50% of the continental crust of the Borborema Province formed (Fig. 2C). During the
                                   Early Neoproterozoic, from 1.0 to 0.7 Ga, new juvenile additions accounts for ~ 20% of the crust. This period
                                   includes the early stages of the Santa Quitéria arc at ca. 0.85 G­ a24 and magmatism related to the Cariris Velhos
                                                               41
                                   event from 1.0 to 0.92 G ­ a . The production of new continental crust during the main Neoproterozoic orogenic
                                   period, between 0.67 and 0.52 Ga, accounts for only ~ 15% of the new continental crust in the province, sug-
                                   gesting the orogeny was characterized by whole-sale lithospheric reworking with minor juvenile magma input.
                                       The compiled dataset of detrital zircons, including 2147 grains from 60 samples (Supplementary Data 2) from
                                   the central and southern Borborema Province, defines three main groups (G1, G2 and G3) based on the young-
                                   est zircon age and source areas given by the older zircon dates in each group (Fig. 3A,B). The first group (G1) is
                                   dominated by samples with youngest zircon ages > 1.4 Ga and source areas with dominantly Paleoproterozoic
                                   and Archean ages. The second group (G2) is characterized by samples with youngest zircon grains between 0.7

Scientific Reports |   (2021) 11:5719 |                https://doi.org/10.1038/s41598-021-84703-x                                                    3

                                                                                                                                                Vol.:(0123456789)
Decratonization by rifting enables orogenic reworking and transcurrent dispersal of old terranes in NE Brazil
www.nature.com/scientificreports/

                                            Figure 2.  Zircon U–Pb, Nd whole-rock data and crustal growth. (A) Similarity of zircon U–Pb crystallization
                                            age distribution for the North São Francisco Craton, Benino-Nigerian Shield and Borborema Province
                                            (see “Methods” section for explanations on data compilation). (B) Zircon U–Pb crystallization ages for
                                            the Neoproterozoic magmatism in the Borborema Province and their εNd(t) values and ­TDM (given by the colour
                                            of the circles). Data show increased involvement of old crust during the Brasiliano Orogen, as oceanic basins
                                            closed, and collision progressed (see “Methods” section for data compilation). (C) Continental growth for the
                                            Borborema Province using data from samples of igneous and metaigneous rocks (solid lines) and exposed
                                            area coverd by these rocks based on pixel values from the maps in Fig. 1D, E (dashed lines). Colour of the lines
                                            represent cumulative U–Pb ages (red), εNd(t) < 0 (black) and Nd model ages (green), binned for each 100 Myr
                                            interval from 3.4 to 0.5 Ga. Growth is estimated from the cumulative proportion of juvenile crustal addition
                                            using the εNd(t) cut-off value of zero for the igneous and metaigneous rock samples and gridded pixels for the
                                            crystalline igneous and metaigneous areas from the maps in Fig. 1D, E (see “Methods” section).

                                            and 1.0 Ga with source areas dominated by rocks ranging from 0.9 to 1.0 Ga. A subgroup of samples in G2 also
                                            shows source areas dominated by Paleoproterozoic and Archean rocks such as those in G1. The last group (G3)
                                            comprises samples with youngest zircon ages between 0.5 and 0.7 Ga, representing the youngest sedimentary
                                            record in the province, associated with syn-orogenic sedimentation. The source areas of G3 vary from rocks with
                                            ages in the range of 0.65 to 0.52 Ga, such as the foreland strata of the Sergipano fold-and-thrust ­belt16, to rocks
                                            with source components from 1.0 to 0.65 Ga and even older (Fig. 3A and B).

                                            Transpression and correlation of terranes. Another aspect revealed by our data compilation and geo-
                                            physical image investigation is the link between old Paleoproterozoic and Archean terranes across the province.
                                            The Southern Borborema Orogen (Figs. 1A,B, 4), is split into the Sergipano and Riacho do Pontal fold-and-thrust
                                            belts by a promontory of the São Francisco Craton, well-defined in the Nd isotopic model age map (Fig. 1B) and
                                            geological ­maps15,16. This promontory forms the Entremontes block and may mark an inherited feature of the
                                            paleocontinental passive margin inverted during onset of the Southern Borborema Orogen. It comprises mainly
                                            Archean and Paleoproterozoic gneisses, migmatites and supracrustal ­rocks42.
                                                The combination of our field data, geological maps, isotopic data and geophysical images show that the
                                            Entremontes block is bound by the Pernambuco shear zone in the north and the Boa Vista shear zone in the
                                            south (Fig. 4A). Movement on these shear zones forced rotation and internal deformation of the block under a
                                            transpressive regime (Fig. 4B) leading to the folds documented in the aeromagnetic images at wavelengths of 5
                                            to 20 km with 2D axial planes subparallel to the shear zones (Fig. 4B). In outcrop, folds are tight and fold a pre-
                                            vious foliation (Sn), generating a steep axial plane (Sn + 1) permeated by leucosomes indicative of syn-tectonic
                                            partial melting (Fig. 4B). Stretching lineation has low rake and is commonly parallel to the fold axes, plunging

          Scientific Reports |   (2021) 11:5719 |                https://doi.org/10.1038/s41598-021-84703-x                                                   4

Vol:.(1234567890)
www.nature.com/scientificreports/

                                  Figure 3.  Provenance patterns of sedimentary rocks from central and south Borborema Province. (A)
                                  Cumulative probability plot showing the distribution of U–Pb ages from detrital zircon grains for samples
                                  from the central Borborema Province (region between the Pernambuco and Patos shear zones in Fig. 1). (B)
                                  Same as A, but for the southern Borborema Province (region between Pernambuco and the northern margin
                                  of the São Francisco Craton in Fig. 1). The distributions show three different patterns denoting different source
                                  areas: G1 cratonic provenance only, G2 Tonian extensional magmatic rocks plus cratonic provenance, and G3
                                  Neoproterozoic orogenic magmatic rocks plus a variety of older sources. N is the number of rock samples; n
                                  is the total number of zircon grains analyzed. The bars on top of the diagrams depict the range of the youngest
                                  zircon grains found in each group.

                                  at low angles to WSW or ENE. Kinematic indicators on Sn + 1 planes indicate consistent top-to-NE or E trans-
                                  port in the Entremontes block (yellow arrows in Fig. 4A). This direction could either be a result of rotation of
                                  early-formed shear zones, or formed contemporaneously with transcurrent movement, defining a constrictional
                                  transpressional belt. The metasedimentary Riacho do Pontal fold-and-thrust belt, west of the Entremontes block,
                                  includes the Afeição domain comprising Tonian metasedimentary and metaigneous rocks. In contrast to the
                                  Entremontes, the direction of tectonic transport for the entire Riacho do Pontal belt is dominantly top-to SW,
                                  thrusting the belt over the ­SFC43 (Fig. 4A).
                                      The geophysical and Nd isotopic characteristics of the ancient Entremontes block and the Tonian-age Afeição
                                  domain can be recognized in the central Borborema Province, north of the Pernambuco shear zone, displaced
                                  dextrally by ≈ 200 km and deformed into sigmoidal terranes characteristic of the BP (Fig. 4, Supplementary
                                  Fig. S1). U–Pb ages and εNd(t), as well as K-Th-U contents, as recorded in the radiometric images (Supplementary
                                  Fig. S1), show that the Afeição d­ omain43 can be linked with the contemporaneous Alto Pajeú ­terrane13, and the
                                  Entremontes block can be linked with the Alto Moxotó terrane, also composed of Paleoproterozoic and Archean
                                  gneisses and ­migmatites25,42. The Alto Moxotó terrane includes high-grade Late Neoproterozoic metasedimentary
                                  rocks of the Surubim Complex that can be correlated to the lower-grade equivalents of the Riacho do Pontal
                                  fold-and-thrust belt and of the Cabrobó Complex overthrusting the Entremontes block. The Surubim Complex
                                  and these lower grade-rocks share similar detrital zircon age ­signatures22,42,43.
                                      In summary, much of the Borborema Province has old Nd model ages and zircon population signatures
                                  similar to the São Francisco Craton/Benino-Nigerian Shield (Figs. 1 and 2A), suggesting direct or indirect deri-
                                  vation from these areas. The stepwise process of breakdown and involvement of cratonic blocks is preserved at
                                  the margin of the craton, where increased deformation of cratonic blocks is recorded across the Pernambuco
                                  shear zone and associated with a dextral displacement of ≈ 200 km.

                                  Discussion: decratonization, terrane dispersion and reworking
                                  In the Borborema Province, crustal extension is marked by mafic–ultramafic intrusions (e.g., the ca. 0.90 Ga
                                  Brejo Seco ­Unit14 and the ca. 0.82 Ga Monte Orebe ­Unit29), continental rift-like basic volcanic rocks (e.g., the
                                  ca. 0.88 Ga Paulistana ­Complex15) and A-type orthogneisses (e.g., the ca. 0.87 Ga Pinhões ­pluton44), follow-
                                  ing intraplate 1.0–0.92 Ga A-type granitoids and bimodal v­ olcanism21,41. Together this magmatism defines the
                                  Cariris Velhos event that records the initial extension of cratonic l­ithosphere21 (Fig. 5A,B). Mafic dikes of the
                                  Bahia-Gangila LIP intruding the north São Francisco Craton at 0.92–0.90 G      ­ a45 further support the interpreta-
                                  tion that this was a widespread extensional event. Depocenters in this dynamic extensional environment were
                                  favorably filled with detritus shed by the topographic highs of the cratonic lithosphere and/or by local supply of
                                  the 1.0–0.87 Ga igneous rocks, giving rise to the detrital zircon pattern observed in groups G1 and G2 (Figs. 3
                                  and 5B). However, sampling of older Paleoproterozoic metasedimentary rocks can also generate the G1 pattern.
                                  The geochemical ­signature14,29 of the Brejo Seco and Monte Orebe units, tectonically emplaced in the distal por-
                                  tion of the SFC passive margin, points to a continental rather than an oceanic origin, typical of continental margin
                                  ophiolites of Ligurian and Western Alpine e­ xamples46. The last pulses of rifting are recorded by the 0.72–0.68 Ga
                                  mafic–ultramafic intrusive rocks of the Canindé Group in the Southern Borborema ­Province16, but the relation

Scientific Reports |   (2021) 11:5719 |                https://doi.org/10.1038/s41598-021-84703-x                                                     5

                                                                                                                                               Vol.:(0123456789)
www.nature.com/scientificreports/

                                            Figure 4.  Transpression and correlation of terranes in the central and southern Borborema Province. (A)
                                            Sm–Nd ­TDM age map of the south Borborema Province and northern edge of the São Francisco Craton with
                                            magnetic first derivative image as background. Inset illustrates dextral displacements of the 1.0 to 0.8 Ga and
                                            the older 2.8–2.0 Ga terranes across the Pernambuco shear zone, resulting from transcurrent tectonics. The
                                            yellow arrows indicate the direction of tectonic transport. BV Boa Vista shear zone, PE Pernambuco shear zone,
                                            E Entremontes block, AF Afeição domain; RP Riacho do Pontal fold-and-thrust belt; AP Alto Pajeú terrane;
                                            AM Alto Moxotó terrane, SFC São Francisco Craton. The stereograms report the Sn + 1 foliation along with
                                            its stretching lineation (red arrows) and shear sense when available. The Sm–Nd T ­ DM ages were gridded in
                                            the ArcGIS software using Inverse-Distance-Weighted Interpolation (IDW). (B) Block diagram illustrating
                                            the transpressive regime in the Entremontes block with progressive displacement and increased deformation
                                            from the craton to the Entremontes block and to the Alto Moxotó terrane across the Boa Vista (BV) and the
                                            Pernambuco (PE) shear zones. Note how the Tonian rocks follow the same increase in deformation from the
                                            Afeição (AF) domain to the Alto Pajeú (AP) terrane. Together, they indicate a dextral displacement across the
                                            Pernambuco shear zone of ≈ 200 km.

          Scientific Reports |   (2021) 11:5719 |               https://doi.org/10.1038/s41598-021-84703-x                                                6

Vol:.(1234567890)
www.nature.com/scientificreports/

                                  Figure 5.  Tectonic evolution and formation of the northern São Francisco passive margin. (A) Initial extension
                                  and partial melting of the sub-continental lithospheric mantle (SCLM) generating mafic magmas associated
                                  with the early Cariris Velhos granitoids and filling of depocenters with sediments shed by cratonic sources
                                  (G1 pattern of detrital zircon ages in Fig. 3). (B) Continuous extension and necking of the cratonic lithosphere
                                  associated with intrusion of late Cariris Velhos granitoids, extrusion of volcanic rocks and sedimentation of
                                  strata derived from erosion of Cariris Velhos magmatic rocks as well as cratonic rocks (G2 pattern of detrital
                                  zircon ages in Fig. 3). (C) Breakup of the cratonic lithosphere and formation of the PEAL continental ribbon
                                  and separation of the Benino-Nigerian Shield (BNS). This stage marks development of the São Francisco Craton
                                  passive margin accompanied by marginal mafic–ultramafic rocks of the Brejo S­ eco14 and Monte ­Orebe29 units
                                  interpreted as exhumed sub-continental lithospheric mantle. At this stage the margin is supplied by mixture
                                  of detrital material from the craton and Cariris Velhos magmatism. (D) Inversion of the passive margin
                                  basin and initial collision of the São Francisco Craton against the PEAL. At this stage sedimentary strata are
                                  dominated by Brasiliano-age orogenic detritus as well as detritus from the older sequences (G3 pattern of
                                  detrital zircon ages in Fig. 3). (E) Yield stress ­envelope54 for continental lithosphere as a function of the time
                                  since the last thermotectonic event (1.0 Ga: dashed blue lines) and (0.25 Ga: dashed yellow lines, corresponding
                                  approximately to the time gap between Tonian extension and start of the orogeny in the BP). The stress
                                  envelopes show that during the Brasiliano orogenic event, the cratonic area would be stronger than the regions
                                  that underwent the more recent Tonian thermotectonic event. This implies that the orogenic deformation would
                                  focus on the weaker decratonized lithosphere.

                                  to the development of the SFC passive margin is still contentious. The absence of suprasubduction-type ophi-
                                                                     ­ yprus46), suggest that extension did not lead to the formation of a large ocean
                                  olites (e.g., Troodos ophiolite in C
                                  separating the rifted blocks. Instead, the break-up of the greater São Francisco Craton/Benino-Nigerian Shield
                                  resulted in local sub-continental lithospheric mantle (SCLM) exhumation followed by formation of an embryonic
                                  narrow ocean, the Sergipano ­Ocean16 (Fig. 5B and C). This is similar to the Alpine-Apennine poorly evolved
                                  oceanic basins, as also described in other Neoproterozoic passive m    ­ argins47. If this is the case, the extensional
                                  event between 1.0 and 0.82 Ga must have been intermittent so as to prevent the opening of a large ocean. This
                                  intermittent extension pulled away crustal ribbons from the cratonic margin, such as the Pernambuco-Alagoas
                                  terrane (PEAL, Fig. 5C), that were subsequently deformed during inversion and final collision (Fig. 5D).
                                      The break-up of the conjoined São Francisco Craton/Benino-Nigerian Shield during extension indicates that
                                  decratonization started already in the Early Neoproterozoic, in a manner akin to the Mesozoic decratonization of
                                  the North China ­Craton9, where extension associated with ­metasomatism10 reactivated and replaced lithospheric

Scientific Reports |   (2021) 11:5719 |                https://doi.org/10.1038/s41598-021-84703-x                                                      7

                                                                                                                                                  Vol.:(0123456789)
www.nature.com/scientificreports/

                                            mantle by a­ sthenosphere48 and facilitated continental b   ­ reakup49. This hypothesis is supported by alkaline volcan-
                                            ism (e.g., kimberlites) in the northern SFC from 1.15 to 0.68 G      ­ a50–52.
                                                  The long-term effect of extension on lithospheric strength depends on the relative thinning of the crust in rela-
                                             tion to the sub-continental lithospheric mantle (SCLM)53,54. Starting from a thick Archean lithosphere, as for the
                                            ­SFC55,56, Tonian thinning and refertilization of the SCLM (indicated by the alkaline volcanism), would have led
                                            to weakening that persisted long-after extension ended. Only after 0.40 to 0.75 Gyr from a given thermotectonic
                                            event, the temperature distribution in the continental lithosphere approaches equilibrium and stops evolving
                                            with ­time53 (Fig. 5E). Thus, the Tonian extensional event would have made sections of the cratonic lithosphere
                                             amenable to subsequent reworking within the orogenic realm.
                                                  This weakened, decratonized lithosphere between the Benino-Nigerian Shield and the São Francisco Craton
                                             (Fig. 6A) was reworked in two steps during the construction of the Borborema Province Orogen. The first step
                                             was a result of oblique continental collision with the West African Craton being thrust under the Benino-Nigerian
                                             ­Shield17,27,57,58. This was preceded by oceanic lithosphere subduction and magmatism in the aforementioned Santa
                                            Quitéria arc, and terminal collision led to the development of the dextral Transbrasiliano-Kandi strike-slip belt
                                            (Fig. 6B and C). The second step was a result of dextral movement on the Transbrasiliano-Kandi strike-slip belt.
                                            This movement brought the Benino-Nigerian Shield closer to the São Francisco Craton deforming the weakened
                                            decratonized area in between these two stiffer lithospheric domains and leading to the transpressional Borborema
                                            Province Orogeny (see cross-sections in Fig. 6D). The closure of the small intervening Sergipano Ocean (Fig. 5C),
                                            prior to collision and transpression, resulted in voluminous magmatism between 0.64 and 0.60 Ga, especially in
                                            the PEAL terrane and Sergipano fold-and-thrust ­belt16,37. This abundant magmatism could have resulted from
                                             increased ­H2O flux released from the subduction of sediments and serpentinized mantle from the extended
                                             margin of the approaching SFC, in a similar manner to that described in the transition from subduction to
                                             continental collision in the West Gondwana ­Orogen59. Erosion of this orogenic system contributed to the 0.65-
                                            0.52 Ga sedimentary input in the Southern Borborema Orogen, as observed in detrital zircon pattern of group
                                            G3 (Figs. 4 and 5B), while orogenic contribution in the central Borborema Province could represent a mixing
                                              between the Southern Borborema and West Gondwana Orogens.
                                                  The major transcurrent shear zones that characterize the Borborema Province splay out of the Transbrasi-
                                              liano-Kandi shear z­ one19 and deformed the weakened cratonic margin. Deformation was a result of these two
                                             quasi-contemporaneous ­collisions19,27 (Fig. 6D). In the north part (north of the Patos shear zone in Fig. 1B),
                                             old Paleoproterozoic and minor volumes of Archean rocks dominate and represent the reworked southern con-
                                             tinuation of the Benino-Nigerian S­ hield17. To the south, the Borborema Province records the interaction of the
                                            decratonized lithosphere and the São Francisco Craton. In this region, the transcurrent shear zones wrenched the
                                            blocks of the northern margin of the São Francisco Craton previously weakened by the Cariris Velhos extension.
                                            Deformation included a number of continental ribbons (e.g., PEAL terrane) pulled away from the craton, and
                                            the intervening Tonian sedimentary basins (Fig. 6C and D).
                                                  The exact geometry of the decratonized terrane at the start of wrenching, and the source region of individual
                                            blocks now embedded in the BP remain undetermined. However, using rock ages, isotopic and geophysical sig-
                                            natures, we have linked the Tonian-age Afeição domain and the ancient Entremontes block south of the Pernam-
                                            buco shear zone, with the Tonian-age Alto Pajeú and the ancient Alto Moxotó terranes, north of the shear zone.
                                            This implies ≈ 200 km dextral wrenching across the Pernambuco shear zone (Fig. 4). The rocks south of the shear
                                            zone record only incipient deformation, whereas those to the north are intensely strained into sigmoidal terranes,
                                            in harmony with the regional transcurrent deformation and the geological-geophysical structure of this region.
                                                  It is also possible to recognize signatures of specific sections of the SFC in the old blocks now within the BP.
                                            For example, the geological features of the Gavião Block in the SFC and of the Kaduna Massif in the Benino-
                                            Nigerian ­Shield62 can be found in blocks within the Central and Northern Borborema ­Province38,40,60. These
                                            features include rare occurrences of Paleoarchean rocks (> 3.4 Ga) embedded in Neoarchean (2.8–2.6 Ga) to
                                            Paleoproterozoic (2.1–2.0 Ga) rocks imprinted by ca. 2.0 Ga high-grade m         ­ etamorphism60,61. We conclude that
                                             the stepwise increase in deformation intensity of craton margin blocks, illustrated by the Entremontes block and
                                             the sigmoidal-shaped Alto Moxotó terrane, illustrates how a number of other Archean-Paleoproterozoic blocks
                                             may have been pulled away from the original craton to form inliers within the orogeny. It is important to note
                                             that dispersal of decratonized blocks from the SFC was more effective along the southern and central zones of
                                             the Borborema Province, where evidence for Cariris Velhos extensional events has been better defined.
                                                  Decratonization related to Neoproterozoic extension and juvenile magmatism could have been widespread
                                             throughout West Gondwana, related to the break-up of Rodinia, although the position of the SFC in this super-
                                             continent is c­ ontentious63. For example, it might account for thinning of the cratonic lithosphere in the Ribeira
                                            and Araçuaí belts of southeastern ­Brazil64, decratonizing the eastern margin of the SFC, preparing it for later
                                            involvement in the transcurrent Brasiliano t­ ectonics65. In the Araçuaí Belt fissural mafic magmatism, associated
                                             with the Bahia-Gangila LIP event, preceded the opening of the Adamastor Ocean starting at ca. 0.9 G         ­ a45,66 with
                                            subsequent ca. 0.87 Ga rift-related, A-type continental ­plutonism . Further south, in the Dom Feliciano-Kaoko
                                                                                                                     31

                                            belt, rift-related siliciclastic and bimodal volcanic rocks preserved in the Neoproterozoic schist belts, from both
                                            the Rio de la Plata/Paranapanema and Congo/Kalahari cratonic margins, suggest continental rifting between 0.9
                                            and 0.78 G    ­ a67. In the African side of these orogens (e.g., Gariep, Kaoko and Damara-Lufilian belts) extensional
                                            tectonics and breakup of surrounding cratons are also constrained to between 0.85 to 0.77 G         ­ a68–70. These pro-
                                             tracted extensional events preceding the Brasiliano/Pan-African orogens of coastal South America and African
                                             equivalents, disrupted and weakened the surrounding cratons, enabling orogenic reworking and transcurrent
                                             dispersal of old terranes such as the São Luiz, Curitiba, and Cabo Frio t­ erranes71,72. The same could account for
                                            Archean/Paleoproterozoic lithosphere in the West Gondwana Orogen in Africa where 1.1–1.0 Ga extension in
                                            the Tuareg Shield allowed reworking during the Pan-African ­Orogeny2.

          Scientific Reports |   (2021) 11:5719 |                 https://doi.org/10.1038/s41598-021-84703-x                                                        8

Vol:.(1234567890)
www.nature.com/scientificreports/

                                  Figure 6.  Neoproterozoic decratonization, transcurrent tectonics and terrane dispersion in NE Brazil. (A) Overall extension
                                  of the São Francisco Craton/Benino-Nigerian Shield (SFC and BNS), following intrusion of intraplate Cariris Velhos granitoids
                                  (1.0–0.92 Ga in orange) and leading to development of Tonian rift basins in restricted oceanic realms at 0.85–0.70 Ga (c.f. Sergipano
                                  Ocean recorded by the ca. 0.82 Ga Monte Orebe o        ­ phiolite28 (blue star in C) and aulacogens such as the 0.9–0.5 Ga Paramirim (P).
                                  Inception of peripherical intraoceanic and transitional arcs in the west along the Goiás-Pharusian Ocean. The Gavião Block (GB)
                                  is the old Archean cratonic core of the conjoined São Francisco Craton (SFC), Congo Craton (CC) and the Benino-Nigerian Shield
                                   (BNS). The Pernambuco-Alagoas terrane (PEAL) represents a continental ribbon of SFC rocks within the rift, that was extensively
                                   intruded by Neoproterozoic granitoids. Note the Entremontes (E) promontory on the paleomargin of the SFC. (B) Subduction stage
                                   marked by intrusion of continental arc related granitoids such as the Santa Quitéria and Major Isidoro/Betania arcs in the western and
                                   southern Borborema Province, respectively. Also shown are the Archean domains of Kaduna (K), Campo Grande (C) and São José
                                   do Campestre (SJ) and the Nigerian schist belts of the Benino-Nigerian Shield. (C) Overall shortening of the interior basins (Tonian
                                   rift + Brasiliano/Pan-African Ediacaran orogenic-related basins) as a response to the oblique continental collision along the West
                                   Gondwana Orogen (WGO) marked by ultra-high-pressure (UHP) metamorphism, and nucleation of the dextral Transbrasiliano-
                                   Kandi shear zone (TK)17,19 as a result of continental collision in the west. This stage is also marked by continental arc magmatism
                                  (0.63-0.60 Ga) in the Southern Borborema ­Orogen16 and shearing along the inactive Santa Quitéria a­ rc17. (D) Zoom on NE Brazil and
                                  NW Africa illustrating final shortening, inversion of rift basins and development of the Borborema transcurrent shear zone system
                                  as a consequence of the double cratonic indentation associated with continued convergence of the West African Craton against the
                                  greater São Francisco-Saharan paleocontinent in the west and a counterclockwise rotation of the São Francisco Craton in the south.
                                  Transpression deformed and dispersed terranes sourced from the SFC within the Borborema Province. Transcurrent shear zones
                                  are associated with crustal heating and magmatism associated with delamination of the lithospheric mantle and asthenospheric
                                  ­upwelling80,81 that facilitated lateral crustal flow. Position of the cross-sections (on the right-hand-side of D) are indicated by pink
                                  dashed lines in A to D. Arrows and their size indicate direction and relative magnitude of the stresses, respectively. Shear zones: TK
                                  Transbrasiliano-Kandi, PE Pernambuco, PT Patos, SP Senador Pompeu, BV Boa Vista. Disrupted Tonian and cratonic domains: AF
                                  Afeição, AP Alto Pajeú, AM Alto Moxotó, T Tróia, G Granjeiro, PI Cristalândia do Piauí. PM in cross section (g) to (h) is partial
                                  melting.

Scientific Reports |   (2021) 11:5719 |                   https://doi.org/10.1038/s41598-021-84703-x                                                                     9

                                                                                                                                                                   Vol.:(0123456789)
www.nature.com/scientificreports/

                                                 Global estimates for the construction of Gondwana between 0.6 and 0.5 Ga indicate only minor mantle
                                            ­addition73, in accordance with our observations in the Borborema Province. The final stage of reworking and
                                             transpression was accompanied by voluminous syn-transcurrent high-K calc-alkaline magmatism ranging from
                                             mafic to felsic rocks from 0.59 to 0.56 G   ­ a74–76, peaking at ca. 0.58 Ga (Fig. 2B). Their strongly negative εNd
                                             values (Fig. 2B) and radiogenic Sr isotopes suggest a metasomatized, enriched lithospheric mantle s­ ource74.
                                             Lead isotopic ratios provide complementary information, suggesting involvement of asthenospheric fluids, pos-
                                             sibly responsible for triggering voluminous lithospheric ­melting76. The orogenic collisional period of the BP
                                             resulted in crustal thickening, especially along the West Gondwana O     ­ rogen27, however, Rayleigh wave tomogra-
                                             phy indicates that this orogen is marked by thinner lithosphere ­today77, thus suggesting that the orogenic roots
                                             have been removed after collision. Since this magmatism from 0.59 to 0.56 Ga is contemporaneous with transcur-
                                             rent deformation, post-dating the collisional event recorded by the West Gondwana Orogen, but synchronous
                                             to the collision along SBO, their origin can be attributed to the disturbance of the lithospheric mantle, including
                                             asthenospheric fluid ­influx76, unrelated to subduction ­processes78. This disturbance alongside collisional belts
                                             might be caused by lithospheric delamination, due to eclogitization at the base of the crust or by convective
                                             removal of a thickened thermal boundary l­ ayer79. The crustal structure of the province revealed by receiver func-
                                             tions and surface wave dispersion suggest that the delamination of the thickened lower crust might have occurred
                                                                   ­ eriod80 in combination with asthenospheric u
                                             after the collisional p                                              ­ pwelling81. Such delamination, caused renewed
                                             lithospheric weakening that facilitated lateral crustal ­flow82 contributing to the processes of escape tectonics and
                                             terrane dispersion in the Borborema ­Province19 (cross-section (c) to (d) in Fig. 6B). Finally, Ar–Ar cooling ages
                                             and U–Pb zircon emplacement ages (Fig. 2B) of poorly deformed to isotropic granitoids indicate slow cooling
                                             rates with continuous heat supply related to the delamination process until the Cambrian (0.52–0.50 Ga)19.
                                                 We conclude that the cratonic roots of the São Francisco Craton/Benino-Nigerian Shield, responsible for
                                             craton integrity, were weakened by Tonian-age extension. This created the conditions required for the involve-
                                             ment and dispersal of decratonized inliers within the Brasiliano orogen. We suggest that this may have been
                                             the general sequence of events for many of the Brasiliano/Pan-African orogens, where extension related to the
                                             break-up of cratonic masses and opening of oceanic realms with varying degrees of maturity were followed by
                                             convergence, wrench tectonics and late post-collisional magmatism during Gondwana amalgamation.

                                            Methods
                                            Nd isotopic maps. In order to discriminate terranes with similar signatures, a large compilation compris-
                                             ing 360 zircon U–Pb ages and 1331 Sm–Nd whole-rock isotope analyses of samples from the BP and the north-
                                             ern SFC were used (Fig. 1B). Sm–Nd isotope distribution of a large number of samples is suitable to identify
                                             correlated terranes due to resistance of the Sm–Nd isotope system to post-crystallization thermal ­disturbance83.
                                            Data were downloaded from the open sources ­DateView84 and the Geological Survey of Brazil database (http://
                                            geosg​b.cprm.gov.br). Due to the scarcity of U–Pb crystallization ages from the Benino-Nigerian Shield, for this
                                            region we used individual zircon ages for dated metaigneous rocks available in the global zircon compilation
                                            ­in85. The compilation was augmented with data from the literature. Most of the data are from (meta)igneous
                                            rocks with subordinated input from metasedimentary rocks.
                                                 In general, Nd model ages ­(TDM) do not correspond to a specific crust-formation event but instead reflect
                                            mixing of material derived from the mantle at different times and are determined by calculating the time when
                                            a sample had an isotopic composition identical to that of its ­source86, so they can be understood as minimum
                                            ages of crust formation. For Fig. 1B we used Sm–Nd T    ­ DM ages as originally reported by different authors (Sup-
                                            plementary Data 1). The Sm–Nd T    ­ DM ages were gridded in the ArcGIS software using Inverse-Distance-Weighted
                                            Interpolation (IDW). Since significant discontinuities modify the surrounding geological environment, we define
                                            the major shear zones as interpolation barriers. We also applied the Gaussian low-pass filter to attenuate high
                                            frequencies due to variable spacing between samples. We compared compiled U–Pb zircon crystallization ages
                                            and freely available geological maps (http://geosg​b.cprm.gov.br) of the Northern SFC and BP to cross-check
                                            terrane affinity based on the Sm–Nd ­TDM map, to identify old terranes within the Borborema orogenic province.
                                            From the compiled Sm–Nd dataset, only data with reported ppm contents of Nd, Sm and 143Nd/144Nd ratio
                                            associated with reliable reported U–Pb ages for the magmatic crystallization, were considered for building the
                                            εNd(t) maps and further crustal growth curve for the Borborema Province. This resulted in 889 data points with
                                            assigned geographical position and recalculated 147Sm/144Nd ratio. A new screening was applied to eliminate
                                            unreliable 147Sm/144Nd resulting in the 837 data points that were used for further gridding of εNd(t) values, as
                                             described above (Fig. 3). The gridded values resulted in a map with 46,426 pixels with a pixel size of 3.2 × 3.2 km
                                             (pixel area = 10.24 ­km2). These pixels where clipped to the igneous crystalline area of the Borborema Province,
                                            excluding the Neoproterozoic metasedimentary belts, Phanerozoic basins and more recent cover, resulting in a
                                            map with 17,235 pixels with a total area of 176,486 ­km2. Pixel values were further binned to 100 Myr intervals,
                                            from 0.5 to 3.4 Ga, along with the U–Pb crystallization ages from individual rock samples. The εNd(t) cut-off value
                                            of zero was used to distinguish juvenile addition of new crust from reworked crust in the Borborema Province
                                            through time. The pixels with εNd > 0 were binned to 100 Myr and the cumulative percentage area covered by
                                            juvenile rocks was used to infer the crustal growth of Borborema Province from 0.5 to 3.4 Ga. Supplementary
                                            Figure S1 shows the data distribution used to grid the Sm–Nd T    ­ DM ages.

                                            Magnetic maps.          The airborne magnetic database comprises data from seven surveys between 2001 and
                                            2010, with 500 m flight-line spacing in the N–S direction and flight height of 100 m (http://geosg​b.cprm.gov.
                                            br). The Total Magnetic Intensity map (TMI) was created by interpolating the magnetic data into a 125 m grid
                                            cell size using the bi-directional method and subsequently filtered by a Gaussian low-pass filter. To highlight the
                                            regional tectonic framework, we calculated the first vertical derivative of TMI (1VD).

          Scientific Reports |   (2021) 11:5719 |                https://doi.org/10.1038/s41598-021-84703-x                                                    10

Vol:.(1234567890)
www.nature.com/scientificreports/

                                  Data availability
                                  All isotopic data used during the study are available in Supplementary Information file (Supplementary Data 1
                                  and 2). Airborne geophysical data from the Geological Survey of Brazil used in Figs. 1B and 4A, and geological
                                  maps used in Fig. 1A are freely available at http://geosg​b.cprm.gov.br. P-wave tomography model for the Bor-
                                  borema Province in Fig. 1C is available in ­reference35.

                                  Received: 27 November 2020; Accepted: 10 February 2021

                                  References
                                   1. Jordan, T. H. The continental tectosphere. Rev. Geophys. 13, 1 (1975).
                                   2. Liégeois, J. P., Abdelsalam, M. G., Ennih, N. & Ouabadi, A. Metacraton: Nature, genesis and behavior. Gondwana Res. 23, 220–237
                                      (2013).
                                   3. Sonder, L. J. & England, P. Vertical averages of rheology of the continental lithosphere: Relation to thin sheet parameters. Earth
                                      Planet. Sci. Lett. 77, 81–90 (1986).
                                   4. England, P. C. Diffuse continental deformation: Length scales, rates and metamorphic evolution. Philos. Trans. R. Soc. Lond. Ser.
                                      A Math. Phys. Sci. 321, 3–22 (1987).
                                   5. Neil, E. A. & Houseman, G. A. Geodynamics of the Tarim basin and the Tian Shan in central Asia. Tectonics 16, 571–584 (1997).
                                   6. Ziegler, P. A., Cloetingh, S. & van Wees, J.-D. Dynamics of intra-plate compressional deformation: The Alpine foreland and other
                                      examples. Tectonophysics 252, 7–59 (1995).
                                   7. Butler, R. W. H., Holdsworth, R. E. & Lloyd, G. E. The role of basement reactivation in continental deformation. J. Geol. Soc. Lond.
                                      154, 69–71 (1997).
                                   8. Holdsworth, R. E., Butler, C. A. & Roberts, A. M. The recognition of reactivation during continental deformation. J. Geol. Soc.
                                      Lond. 154, 73–78 (1997).
                                   9. Yang, J. H., Wu, F. Y., Wilde, S. A., Belousova, E. & Griffin, W. L. Mesozoic decratonization of the North China block. Geology 36,
                                      467–470 (2008).
                                  10. Wang, C., Song, S., Niu, Y. & Su, L. Late Triassic adakitic plutons within the Archean terrane of the North China Craton: Melting
                                      of the ancient lower crust at the onset of the lithospheric destruction. Lithos 212–215, 353–367 (2015).
                                  11. Kusky, T. M., Windley, B. F. & Zhai, M.-G. Tectonic evolution of the North China Block: From orogen to craton to orogen. In
                                      Mesozoic Sub-Continental Lithospheric Thinning Under Eastern Asia Vol. 280 (eds Zhai, M.-G. et al.) 1–34 (Geological Society,
                                      2007).
                                  12. Griffin, W. L., O’Reilly, S. Y., Afonso, J. C. & Begg, G. C. The composition and evolution of lithospheric mantle: A re-evaluation
                                      and its tectonic implications. J. Petrol. 50, 1185–1204 (2009).
                                  13. Guimarães, I. P. et al. U-Pb zircon ages of orthogneisses and supracrustal rocks of the Cariris Velhos belt: Onset of Neoproterozoic
                                      rifting in the Borborema Province, NE Brazil. Precambrian Res. 192–195, 52–77 (2012).
                                  14. Salgado, S. S. et al. The Ni-Cu-PGE mineralized Brejo Seco mafic-ultramafic layered intrusion, Riacho do Pontal Orogen: Onset
                                      of Tonian (ca. 900 Ma) continental rifting in Northeast Brazil. J. S. Am. Earth Sci. 70, 324–339 (2016).
                                  15. Caxito, F. A. et al. A complete Wilson Cycle recorded within the Riacho do Pontal Orogen, NE Brazil: Implications for the Neo-
                                      proterozoic evolution of the Borborema Province at the heart of West Gondwana. Precambrian Res. 282, 97–120 (2016).
                                  16. Oliveira, E. P., Windley, B. F. & Araújo, M. N. C. The Neoproterozoic Sergipano orogenic belt, NE Brazil: A complete plate tectonic
                                      cycle in western Gondwana. Precambrian Res. 181, 64–84 (2010).
                                  17. Ganade, C. E. et al. Tightening-up NE Brazil and NW Africa connections: New U-Pb / Lu-Hf zircon data of a complete plate
                                      tectonic cycle in the Dahomey belt of the West Gondwana Orogen in Togo and Benin. Precambrian Res. 276, 24–42 (2016).
                                  18. Caxito, F. A. et al. Toward an integrated model of geological evolution for NE Brazil-NW Africa: The Borborema Province and its
                                      connections to the Trans-Saharan (Benino-Nigerian and Tuareg shields) and Central African orogens. Braz. J. Geol. https​://doi.
                                      org/10.1590/2317-48892​02020​19012​2 (2020).
                                  19. Ganade, C. E., Weinberg, R. F. & Cordani, U. G. Extruding the Borborema Province (NE-Brazil): A two-stage Neoproterozoic
                                      collision process. Terra Nova 26, 157–168 (2013).
                                  20. Vauchez, A. et al. The Borborema shear zone system, NE Brazil. J. S. Am. Earth Sci. 8, 247–266 (1995).
                                  21. Neves, S. P. Proterozoic history of the Borborema province (NE Brazil): Correlations with neighboring cratons and Pan-African
                                      belts and implications for the evolution of western Gondwana. Tectonics 22, 1–13 (2003).
                                  22. Neves, S. P. et al. The age distributions of detrital zircons in metasedimentary sequences in eastern Borborema Province (NE
                                      Brazil): Evidence for intracontinental sedimentation and orogenesis?. Precambrian Res. 175, 187–205 (2009).
                                  23. Fetter, A. H. et al. Evidence for Neoproterozoic continental arc magmatism in the Santa Quitéria Batholith of Ceará State, NW
                                      Borborema Province, NE Brazil: Implications for the assembly of West Gondwana. Gondwana Res. 6, 265–273 (2003).
                                  24. Ganade, C. E. et al. Tracing Neoproterozoic subduction in the Borborema Province (NE-Brazil): Clues from U-Pb geochronology
                                      and Sr-Nd-Hf-O isotopes on granitoids and migmatites. Lithos 202–203, 167–189 (2014).
                                  25. Santos, L. C. M. L. et al. Accretion tectonics in Western Gondwana deduced from Sm-Nd isotope mapping of Terranes in the
                                      Borborema Province, NE Brazil. Tectonics 37, 2727–2743 (2018).
                                  26. Ferreira, A. C. D., Dantas, E. L., Fuck, R. A. & Nedel, I. M. Arc accretion and crustal reworking from late Archean to Neoproterozoic
                                      in Northeast Brazil. Sci. Rep. 10, 7855 (2020).
                                  27. Ganade, C. E. et al. Ediacaran 2,500-km-long synchronous deep continental subduction in the West Gondwana Orogen. Nat.
                                      Commun. 5, 5198 (2014).
                                  28. Santos, T. J. S. et al. U-Pb age of the coesite-bearing eclogite from NW Borborema Province, NE Brazil: Implications for western
                                      Gondwana assembly. Gondwana Res. 28, 1183–1196 (2015).
                                  29. Caxito, F., Uhlein, A., Stevenson, R. & Uhlein, G. J. Neoproterozoic oceanic crust remnants in northeast Brazil. Geology 42, 387–390
                                      (2014).
                                  30. Pitombeira, J. P. A. et al. Vestiges of a continental margin ophiolite type in the Novo Oriente region, Borborema Province, NE
                                      Brazil. J. S. Am. Earth Sci. 73, 78–99 (2017).
                                  31. Silva, L. C., Pedrosa-Soares, A. C., Teixeira, L. R. & Armstrong, R. Tonian rift-related, A-type continental plutonism in the Araçuaí
                                      Orogen, eastern Brazil: New evidence for the breakup stage of the São Francisco-Congo Paleocontinent. Gondwana Res. 13,
                                      527–537 (2008).
                                  32. Brito Neves, B. B., dos Santos, E. J. & Van Schmus, W. R. Tectonic history of the Borborema Province. In Tectonic Evolution of South
                                      America Vol. 31 (eds Cordani, U. G. et al.) 151–182 (31st International Geological Congress, Rio de Janeiro, 2000).
                                  33. Brito Neves, B. B., Van Schmus, W. R., Santos, E. J., Campos Neto, M. C. & Kozuch, M. O evento Cariris Velhos na Província
                                      Borborema: Integração de dados, implicações e perspectivas. Rev. Bras. de Geoci. 25, 279–296 (1995).
                                  34. Van Schmus, W. R., Kozuch, M. & Brito Neves, B. B. Precambrian history of the Zona Transversal of the Borborema Province, NE
                                      Brazil: Insights from Sm-Nd and U-Pb geochronology. J. S. Am. Earth Sci. 31, 227–252 (2011).

Scientific Reports |   (2021) 11:5719 |                    https://doi.org/10.1038/s41598-021-84703-x                                                                   11

                                                                                                                                                                   Vol.:(0123456789)
www.nature.com/scientificreports/

                                            35. Simões-Neto, F. L., Julià, J. & Schimmel, M. Upper-mantle structure of the Borborema Province, NE Brazil, from P-wave tomog-
                                                raphy: Implications for rheology and volcanism. Geophys. J. Int. 216, 231–250 (2019).
                                            36. Rocha, M. P. et al. Delimiting the Neoproterozoic São Francisco Paleocontinental Block with P-wave traveltime tomography.
                                                Geophys. J. Int. 219, 633–644 (2019).
                                            37. Silva, T. R., Ferreira, V. P., de Lima, M. M. C., Sial, A. N. & da Silva, J. M. R. Synkinematic emplacement of the magmatic epidote
                                                bearing Major Isidoro tonalite-granite batholith: Relicts of an Ediacaran continental arc in the Pernambuco-Alagoas domain,
                                                Borborema Province, NE Brazil. J. S. Am. Earth Sci. 64, 1–13 (2015).
                                            38. Dantas, E. L. et al. Crustal growth in the 3.4–2.7Ga São José de Campestre Massif, Borborema Province, NE Brazil. Precambrian
                                                Res. 227, 120–156 (2013).
                                            39. Cruz, R. F., Pimentel, M. M., de Accioly, A. C. A. & Rodrigues, J. B. Geological and isotopic characteristics of granites from the
                                                Western Pernambuco-Alagoas Domain: Implications for the crustal evolution of the Neoproterozoic Borborema Province. Braz.
                                                J. Geol. 44, 627–652 (2014).
                                            40. Pitarello, M. Z., do Santos, T. J. S. & Ancelmi, M. F. Syn-to post-depositional processes related to high grade metamorphic BIFs:
                                                Geochemical and geochronological evidences from a Paleo to Neoarchean (3.5–2.6 Ga) terrane in NE Brazil. J. S. Am. Earth Sci.
                                                96, 102312 (2019).
                                            41. Guimarães, I. P. et al. Tonian granitic magmatism of the Borborema Province, NE Brazil: A review. J. S. Am. Earth Sci. 68, 97–112
                                                (2016).
                                            42. Cruz, R. F., Pimentel, M. M. & de Accioly, A. C. A. Provenance of metasedimentary rocks of the Western Pernambuco-Alagoas
                                                Domain: Contribution to understand the crustal evolution of southern Borborema Province. J. S. Am. Earth Sci. 56, 54–67 (2014).
                                            43. Caxito, F. A., Uhlein, A. & Dantas, E. L. The Afeição augen-gneiss Suite and the record of the Cariris Velhos Orogeny (1000–960
                                                Ma) within the Riacho do Pontal fold belt, NE Brazil. J. S. Am. Earth Sci. 51, 12–27 (2014).
                                            44. Neves, S. P. et al. From extension to shortening: Dating the onset of the Brasiliano Orogeny in eastern Borborema Province (NE
                                                Brazil). J. S. Am. Earth Sci. 58, 238–256 (2015).
                                            45. Chaves, A. O., Ernst, R. E., Söderlund, U., Wang, X. & Naeraa, T. The 920–900 Ma Bahia-Gangila LIP of the São Francisco and
                                                Congo cratons and link with Dashigou-Chulan LIP of North China craton: New insights from U-Pb geochronology and geochem-
                                                istry. Precambrian Res. 329, 124–137 (2019).
                                            46. Dilek, Y. & Furnes, H. Ophiolites and their origins. Elements 10, 93–100 (2014).
                                            47. Kotschoubey, B., Hieronymus, B. & de Albuquerque, C. A. R. Disrupted peridotites and basalts from the Neoproterozoic Araguaia
                                                belt (northern Brazil): Remnants of a poorly evolved oceanic crust?. J. S. Am. Earth. Sci. 20, 211–230 (2005).
                                            48. Carlson, R. W., Pearson, D. G. & James, D. E. Physical, chemical, and chronological characteristics of continental mantle. Rev.
                                                Geophys. 43, 1–24 (2005).
                                            49. Tappe, S. et al. Craton reactivation on the Labrador Sea margins: 40Ar/39Ar age and Sr–Nd–Hf–Pb isotope constraints from alkaline
                                                and carbonatite intrusives. Earth Planet. Sci. Lett. 256, 433–454 (2007).
                                            50. Williamson, P. A. & Pereira, R. S. The Salvador 01 Kimberlite, Bahia, Brazil: Its Regional and Local Geological Setting with Com-
                                                ments on the Sequence of Prospecting Activities Leading to Its Discovery. (Strategic Services Unit - De Beers Brasil, Internal
                                                Report, Brásilia, 1991).
                                            51. Pereira, R. S. & Fuck, R. A. Archean nucleii and the distribution of kimberlite and related rocks in the São Francisco Craton, Brazil.
                                                Rev. Bras. Geosci. 35, 93–104 (2005).
                                            52. Donatti-Filho, J. P., Tappe, S., Oliveira, E. P. & Heaman, L. M. Age and origin of the Neoproterozoic Brauna kimberlites: Melt
                                                generation within the metasomatized base of the São Francisco craton, Brazil. Chem. Geol. 353, 19–35 (2013).
                                            53. Tommasi, A. & Vauchez, A. Continental-scale rheological heterogeneities and complex intraplate tectono-metamorphic patterns:
                                                Insights from a case-study and numerical models. Tectonophysics 279, 327–350 (1997).
                                            54. Burov, E. B. Rheology and strength of the lithosphere. Mar. Pet. Geol. 28, 1402–1443 (2011).
                                            55. Grand, S. P. Mantle shear structure beneath the Americas and surrounding oceans. J. Geophys. Res. 99, 11591 (1994).
                                            56. VanDecar, J. C., James, D. E. & Assumpção, M. Seismic evidence for a fossil mantle plume beneath South America and implications
                                                for plate driving forces. Nature 378, 25–31 (1995).
                                            57. Caby, R. Precambrian terranes of Benin-Nigeria and northeast Brazil and the Late Proterozoic south Atlantic fit. In Terranes in
                                                the Circum-Atlantic Paleozoic Orogens (ed. Dallmeyer, R. D.) 145–158 (Geological Society of America, 1989).
                                            58. Castaing, C., Feybesse, J. L., Thiéblemont, D., Triboulet, C. & Chèvremont, P. Palaeogeographical reconstructions of the Pan-
                                                African/Brasiliano orogen: Closure of an oceanic domain or intracontinental convergence between major blocks?. Precambrian
                                                Res. 69, 327–344 (1994).
                                            59. Ganade, C. E. et al. Magmatic flare-up causes crustal thickening at the transition from subduction to continental collision. Com-
                                                mun. Earth Environ. 2, 41 (2021).
                                            60. Oliveira, E. P., McNaughton, N. J., Zincone, S. A. & Talavera, C. Birthplace of the São Francisco Craton, Brazil: Evidence from 3.60
                                                to 3.64 Ga Gneisses of the Mairi Gneiss Complex. Terra Nova 32, 281–289 (2020).
                                            61. Barbosa, J. S. F. & Sabaté, P. Archean and Paleoproterozoic crust of the São Francisco Craton, Bahia, Brazil: geodynamic features.
                                                Precambrian Res. 133, 1–27 (2004).
                                            62. Kröner, A., Ekwueme, B. N. & Pidgeon, R. T. The oldest rocks in West Africa: SHRIMP zircon age for early archean migmatitic
                                                orthogneiss at Kaduna, Northern Nigeria. J. Geol. 109, 399–406 (2001).
                                            63. Kröner, A. & Cordani, U. African, southern Indian and South American cratons were not part of the Rodinia supercontinent:
                                                evidence from field relationships and geochronology. Tectonophysics 375, 325–352 (2003).
                                            64. Heilbron, M. et al. The Ribeira belt. In São Francisco Craton, Eastern Brazil: Tectonic Genealogy of a Miniature Continent (eds
                                                Heilbron, M. et al.) 277–302 (Springer, 2017).
                                            65. Egydio-Silva, M., Vauchez, A., Fossen, H., Gonçalves Cavalcante, G. C. & Xavier, B. C. Connecting the Araçuaí and Ribeira belts
                                                (SE—Brazil): Progressive transition from contractional to transpressive strain regime during the Brasiliano orogeny. J. S. Am. Earth
                                                Sci. 86, 127–139 (2018).
                                            66. Pedrosa-Soares, A. C. & Alkmim, F. F. How many rifting events preceded the development of the Araçuaí-West Congo orogen?.
                                                Geonomos 19, 244–251 (2011).
                                            67. Basei, M. A. S. et al. The tectonic history of the southern Adamastor Ocean based on a correlation of the Kaoko and Dom Feliciano
                                                belts. In Geology of Southwest Gondwana (eds Siegesmund, S. et al.) 63–85 (Springer, New York, 2018).
                                            68. Johnson, S. P. et al. Geochronology of the Zambezi Supracrustal sequence, Southern Zambia: A record of neoproterozoic divergent
                                                processes along the southern Margin of the Congo Craton. J. Geol. 115, 355–374 (2007).
                                            69. Frimmel, H. E., Basei, M. S. & Gaucher, C. Neoproterozoic geodynamic evolution of SW-Gondwana: A southern African perspec-
                                                tive. Int. J. Earth Sci. 100, 323–354 (2011).
                                            70. Konopásek, J., Košler, J., Tajčmanová, L., Ulrich, S. & Kitt, S. L. Neoproterozoic igneous complex emplaced along major tectonic
                                                boundary in the Kaoko Belt (NW Namibia): Ion probe and LA-ICP-MS dating of magmatic and metamorphic zircons. J. Geol.
                                                Soc. Lond. 165, 153–165 (2008).
                                            71. Passarelli, C. R., Basei, M. A. S., Siga, O. & Harara, O. M. M. The Luis Alves and Curitiba Terranes: Continental fragments in the
                                                Adamastor Ocean. In Geology of Southwest Gondwana (eds Siegesmund, S. et al.) 189–215 (Springer, New York, 2018).

          Scientific Reports |   (2021) 11:5719 |                    https://doi.org/10.1038/s41598-021-84703-x                                                                    12

Vol:.(1234567890)
You can also read