Geophysical perspective on the structural interference zone along the Neoproterozoic Brasília and Ribeira fold belts in West Gondwana

Manuscript ID: 20160144. Received in: 07/22/2016. Approved in: 02/08/2017. ABSTRACT: The Brasília and Ribeira fold belts have been esta‐ blished in south‐southwestern São Francisco Craton during the Brasiliano‐Pan African orogeny (0.9–0.5 Ga – Tonian to Cam‐ brian), and played an important role in West Gondwana continent assembly. The region is given by a complex regional fold and thrust belt superposed by shearing during the orogeny late times, with superposing stress fields forming a structural interference zone. These thrust sheets encompasses assemblies from lower‐ to upper‐ ‐crust from different major tectonic blocks (Paranapanema, São Francisco), and newly created metamorphic rocks. Re‐evaluation of ground gravity datasets in a geologically constrained approach in‐ cluding seismology (CRUST1 model) and magnetic data (EMAG2 model) unveiled details on the deep‐ crust settings, and the overall geometry of the structural interference zone. The Simple Bouguer Anomaly map shows heterogeneous density distribution in the area, highlighting the presence of high‐density, high metamorphic gra‐ de rocks along the Alterosa suture zone in the Socorro‐Guaxupé Nappe, lying amid a series of metasedimentary thrust scales in a regional nappe system with important verticalization along regio‐ nal shear zones. Forward gravity modeling favors interpretations of structural interference up North into Guaxupé Nappe. Compari‐ son to geotectonic models shows similarities with modern accretio‐ nary belts, renewing the discussion.


INTRODUCTION
Geophysical modeling is an approach to get an insight into Earth's crust architecture based on indirect information by means of its physical properties, contributing to classical structural geology evaluations (Jessell & Valenta 1996, Jessell & Jessell 2001).Gravity modeling applications are widespread as data are relatively easy to access via public surveys, with well-established and powerful processing techniques becoming faster and reliable (Nabighian et al. 2005a), also in integrative interpretation in Global Information Systems and 3D-modeling suites.The use of geophysics presents a possibility to unveil crust geometry with patterns sometimes not envisioned just on field observation, such as effects caused by overburden or cryptic structures (Nabighian et al. 2005a,b, Thomas et al. 2008, Stewart & Betts 2010).
The South American platform main crustal growth episode took place in the Brasiliano-Pan African event, spanning from Rodinia breakup in Tonian through West Gondwana (Fig. 1) assembly in Cambrian (Almeida et al. 1981, 2000, Brito Neves & Fuck 2013).Neoproterozoic collision and accretion tectonics resulted in a collage of several elements ageing from Archean to Mesoproterozoic due to diachronic basin formation, its deformation and closure (Brito Neves et al. 1999, Cordani & Sato 1999, Strieder & Suita 1999, Campos Neto & Caby 2000).
This work provides renewed insight into the geophysical framework of the Tocantins and Mantiqueira orogenic systems (Fig. 1) around the Archean São Francisco Craton Southern margin.Its evolution is under discussion since the first widescale work by Almeida et al. (1981) being successively reappraised (Almeida et al. 2000, Brito Neves et al. 1999) and recently reviewed on Brito Neves andFuck (2013, 2014), highlighting its diverse geo-chronological inventory.In fact, most information derives from geochemistry and geochronology in a poorly developed scenario of integrative discussion of structural geology and geophysics.The later counts on contributions by Haralyi and Hasui (1982), Hasui et al. (1993), Mantovani and Brito Neves (2005) and Mantovani et al. (2005) dealing with tectonic elements outline such as cratonic blocks, regional suture zones and fold belts.Although present knowledge on tectonic and structural development during this orogeny is under construction, several questions remain unanswered, mainly on geochronology and crust geometry -especially in the linking point between the two fold belts, which constitutes a structural interference zone (Trouw et al. 2013).A point to stress is that the re-evaluation of previous data, discussed with new evidences from different data sources, presents a good opportunity to settle past questions.

Geological framework
The most of the structural framework in Brazilian crust dates back to the Neoproterozoic Brasiliano/Pan-African event (from 0.9 to 0.5 Ga) represented by Tocantins, Mantiqueira and Borborema orogens and pre-existing cratons, which undergone various degrees of reworking (Brito Neves & Fuck 2013, 2014).The episodes of sedimentation, reworking, accretion and magmatism throughout the complete geotectonic cycle were diachronic, spanning from Tonian to Early Cambrian with a complex structural evolution (Morales 1993, Brito Neves et al. 1999, Strieder & Suita 1999, Brito Neves & Fuck 2013, 2014).
The Brasília Fold Belt represents Tocantins province southern branch, and is given by metasedimentary and metavolcanic rocks thrusted forward São Francisco Craton western side, meanwhile Ribeira and Araçuaí fold belts are situated in its eastern flank (Almeida et al. 1981, 2000, Brito Neves & Fuck 2013, 2014).
The Southern Brasília belt (Fig. 2) presents a discontinuous stratigraphy and overlays the reworked São Francisco   Craton, and is covered itself by the younger Ribeira Fold Belt (590-550 Ma) metasedimentary thrust stacks (Embu and Paraíba do Sul terranes), although a precise limit between Brasília and Ribeira fold belts is still discussed (Heilbron et al. 2008, Trouw et al. 2013).The metasedimentary nappes along the intervening fold and thrust belts are interpreted as former oceanic systems that existed along São Francisco and Paranapanema cratons borders since its sedimentary provenience and age relationships are similar, despite debate on thrust stack stratigraphy (Ebert & Hasui 1998, Peternet et al. 2005, Heilbron et al. 2008, Hasui 2010, Campos Neto et al. 2011, Trouw et al. 2013, Westin & Campos Neto 2013).Previous structural continuity was obliterated by tight folding and regional branching shear belts from Late-Neoproterozoic to Cambrian in two regional scale near vertical systems (Ebert & Hasui 1998) -Campo do Meio (mostly E-W) and Paraíba do Sul (NE-SW) -due to late-orogenic lateral adjustments between São Francisco and Congo cratons.The products of this collision and lateral movement, the Socorro-Guaxupé and Andrelândia nappes, presents low-angle, southwest dipping nappe stacks superposed with sinistral shear deformation along Campo do Meio Shear Belt (Ebert & Hasui 1998).
Nappe systems tectono-stratigraphy is given by the rootless Andrelândia Nappe System in the bottom, Carrancas and Lima Duarte Nappe systems in middle levels, including preserved autochthonous sections, and Socorro-Guaxupé in the uppermost part (Campos Neto et al. 2011, Trouw et al. 2013).Andrelândia Nappe System comprises five thrust slices named in structural upward succession as Lima Duarte, undifferentiated/ indiscriminate, granulitic, Liberdade and Andrelândia-Carmo da Cachoeira.They comprise rocks varying from siliceous and chemical metasediments to ortho-derived metamorphic rocks on greenschist to high-pressure granulite facies (including eclogites) and migmatites (Alagoa and Rio Preto), juxtaposed by thrust surfaces (Campos Neto et al. 2007, 2011).Carrancas and Lima Duarte nappe systems represent the intermediate-level thrust sheets, with quartzites, graphite schists, marbles, ultramafic schists and paragneisses up to amphibolite facies, besides it is also found as autochthonous cover over São Francisco Craton (Campos Neto et al. 2011).Socorro-Guaxupé Nappe comprises an allochthonous lower crust slice, with high pressure/temperature facies metamorphism, given by metatexites, diatexites and granulites with magmatic arc affiliation (~630 Ma) which are older than the metasedimentary nappes (560 Ma) (Campos Neto & Caby 1999, 2000, Campos Neto et al. 2011).
Expressive sin-to post-orogenic, alkaline to calc-alkaline, acid magmatism took place on subduction-related settings in Paraíba do Sul Shear Belt and Socorro-Guaxupé Nappe from Criogenian to Ediacarian forming several suites and massifs, being emplaced mainly along shear zones with inheritance from Mesoproterozoic basement (Ebert et al. 1996, Sial et al. 1999, Mora et al. 2014).
The Alterosa Suture Zone gives the limit between São Francisco and Paranapanema blocks by the end of the Neoproterozoic, comprising the remainder of the subduction zone (Malagutti Filho et al. 1996, Ebert & Hasui 1998).The Alterosa Suture Zone trace is given by granulites from the bottom of Socorro-Guaxupé Nappe and has been modified by extensive NE-SW shear belts (Ebert & Hasui 1998, Campanha & Brito Neves 2004).Eclogites and other lower level crustal rocks are also found in the Andrelândia Nappes System indicating high pressure levels and crust subduction (and therefore, obduction) on the eastern collision front from the Congo Craton, involving middle-crust thrust slicing from 618-575 Ma (Campos Neto et al. 2011).
The older Amparo Complex and Rhyacian orthogneisses comprises high grade gneissic rocks with complex tectonic relationships that occur as slices between the lobes of Socorro-Guaxupé Nappe, mainly differing from the structural settings of the Neoproterozoic units (Heilbron et al. 2008, Valeriano et al. 2008).
The Bambuí Group overlies São Francisco Craton as an extensive low-grade to non-metamorphic mixed-siliciclastic and carbonatic Neoproterozoic platformal cover, with minor deformation in the Brasiliano/Pan-African event (Martins- Neto 2009).
In contrast to the Western part of the South American Platform, the Brazilian crust remained mostly absent of deformation throughout Phanerozoic, except for the formation of several sedimentary basins, and magmatic events developed by far-field stress deformation from the Andean chain (Stampfli et al. 2013, Granot & Dyment 2015).From Neo-Ordovician to Neo-Cretaceous the siliciclastic Paraná intra-cratonic basin is formed over Precambrian crust (mainly Paranapanema Block and adjacent fold belts), also hosting the expressive Serra Geral Formation basic magmatism (dykes, sills and basalts) from Lower Cretaceous (Zalan et al. 1990, Peate 1997).Alkaline magmatism took place during Upper Cretaceous over both metamorphic basement and Paraná Sedimentary Basin, with several volcanic and plutonic members given by alkaline massifs in the study area (Ulbrich and Gomes 1981) -Poços de Caldas (larger on west side, Fig. 2), Itatiaia (the larger in southeast) and Passa Quatro.Further evolution comprises South Atlantic Ocean opening from Cretaceous to the present, including complete break-up of Pangea and oceanic expansion (Granot & Dyment 2015).

Geophysical framework
Geophysical approaches to the Southeastern South American platform started in 1980's engaging on regional structural-geotectonic elements outline (Lesquer et al. 1981, Haralyi & Hasui 1982).From 1990's to 2000's efforts gradually focused to more specific elements in the crust (Hasui et al. 1993, Malagutti Filho et al. 1996, França & Assumpção 2004, Slavec et al. 2004, Mantovani et al. 2005), returning to regional scale works by seismology and potential field methods in the 2010's (Assumpção et al. 2013a, 2013b, Chulick et al. 2013, Marangoni & Mantovani 2013, van der Meijde et al. 2013, 2015).Lesquer et al. (1981) presents first the interpretations over Southern São Francisco Craton geometry, observing a NW-SE Neoproterozoic structural imprint on gravity data over its Archean grain, suggesting the existence of a WNW-ESE trending suture zone (Alfenas granulitic belt).Haralyi and Hasui (1982) provided gravity interpretations of the Archean-Proterozoic framework in eastern Brazil.They describe the presence of a (stable) crustal block marked by a very low Bouguer Anomaly, corresponding to São Francisco Craton with crustal thickening, which is surrounded by higher density belts relating to metamorphic belts with high-grade elements formed in the Brasiliano-Pan Africano event.The authors interpret these higher signal zones as suture zones (Alterosa suture zone) which corresponds to the Alfenas Granulitic Belt by Lesquer et al. (1981).Hasui et al. (1993) proceeded into a more detailed assessment of the crust established during Brasiliano/Pan-African, outlining different orogenic belts, stable entities and regional lineaments with more detail, including geochronological evidences, and developing geotectonic reconstructions further into the plate tectonics panorama.Further work by Malagutti Filho et al. (1996) lacks detail into crustal structure along the subduction zone between São Francisco Craton and the Paraná-Paranapanema Block (Alterosa Suture).Their crustal duplication model depicts the high-grade granulites being thrust up the São Francisco Craton, with crustal thickness varying from 40 to 30 km (thinning to southeast), and including the suture zone displacement along Campo do Meio shear belt.Evidence for the concealed Paraná/ Paranapanema block emerges from gravity processing for the Paraná Sedimentary Basin effect removal (Mantovani et al. 2005, Mantovani & Brito Neves 2005).Their results suggests the block outline as a NE-SW wedge-like feature that acted as an obducted block in relation to São Francisco Craton, with its exposed roots at Socorro-Guaxupé Nappe, in agreement to previous interpretations (Malagutti Filho et al. 1996, Brito Neves et al. 1999).
Estimates for Mohorovičić (Moho) discontinuity depth varies from 34 to 43 km with thicker crust along both São Francisco craton and Paraná Sedimentary Basin, and thinning over the Neoproterozoic belts and the Atlantic coast (França and Assumpçao 2004; CRUST1 model -1 degree global crustal model by Laske et al. 2012, Assumpção et al. 2013a, Chulick et al. 2013, van der Meijde et al. 2013).Several crustal thickness and mantle structure models are available proving consistent, besides differing input datasets, constraining parameters, physical property nature, sample density and data distribution (Assumpção et al. 2013b, van der Meijde et al. 2015).
Regarding the Cretaceous alkaline intrusions geophysical responses, Marangoni and Mantovani (2013) points out that its emplacement was strongly influenced by basement and basin structures interpreted on magnetic and gravity data over the N125 mantle-plume trace azimuth as previously envisaged (Ulbrich & Gomes 1981).Gravity modeling estimates for Poços de Caldas alkaline massif depth ranges from 7 to 8 km, with a maximum (17 km) in its center (Slavec et al. 2004).

DATA AND METHODS
Gravity is a quite simple physical property to interpret as it present straightforward relationships in geological contexts relating to density of earth materials, nevertheless, caution is advised for subjectivity and ambiguity issues (La Fehr & Nabighian 2012).Gravity data processing is a well-established procedure with plenty resourceful descriptions of mathematical aspects, method revisions and good practices (Blakely 1996, Milligan & Gunn 1997, Nabighian et al. 2005a, La Fehr & Nabighian 2012).Processing and interpretation of magnetics and gravity data follows the potential field theory and and are described by relatively similar mathematical operations (Blakely 1996).
Gravity data used in this contribution are given by 1214 ground stations from Universidade de São Paulo (USP) Instituto de Astronomia, Geofísica e Ciências Atmosféricas (São Paulo -Brazil) database acquired with varying station spacing from 8 to 20 km (Fig. 3).The dataset comprises surveys back to 1960's and includes data from Malagutti Filho et al. (1996).These data are stored in tables with location for each station, gravity measurements and its reduced form by International Union of Geodesy and Geophysics 1967 gravity formula and values for Simple Bouguer anomaly (2670 kg.m -3 density) calculated over altimeter elevation data.Gravity station data interpolation to a surface took place in 4 km cell-size array with the minimum curvature algorithm (Briggs 1974).Geophysical data processing took place on Geosoft Oasis Montaj platform, with information layers integration on ESRI Arc Map™ GIS system for comparison to geological data.
Earth Magnetic Anomaly 2-arc min grid (EMAG2, Maus et al. 2009) accounts for a magnetic data compilation collected around the globe by different platforms, mainly ground, satellite and airborne, integrated into a 2 arc-minute resolution surface referenced 4 km above the geoid, and containing long wavelength information (>330 km) from satellite data.The interpretation of the magnetic geophysical method has been reviewed in Blakely (1996) and Nabighian et al. (2005b).These data are useful for deep crust magnetic structure observation, displayed for comparison to the features observed on geology, gravity and crustal thickness data.EMAG2 data are hosted in GMT (Wessel & Smith, 2013) software-compatible gridded surfaces for visualization.CRUST1 model (Laske et al. 2012), data comprises point information on crustal thickness (Moho-depth) determinations by seismology studies for a world-wide model from several sources (active, passive, long and short-term stations).These data were interpolated into 1-degree cell size surfaces with the minimum curvature algorithm (Briggs 1974).The discussion of the development of both EMAG2 and CRUST1 models is not intended in this contribution and reader should refer to the mentioned publications.
Preliminary observation of the gravity dataset showed presence of high-frequency noise over some data stations.Its analysis with the Power spectrum method (Spector & Grant 1970) in a gridded surface allowed identification of spurious high-frequency noise.The spurious data have been suppressed through Nyquist frequency removal by low-pass filtering with 0.06-cycle frequency cut-off, resulting in the Simple Bouguer Anomaly map refered from here onwards as Bouguer Anomaly (Fig. 4).This map presents smooth contours, being absent of local high amplitude artifacts.
Two and a half dimensional modeling took place in Geosoft's Oasis Montaj GM-SYS platform over information on the Simple Bouguer anomaly map (Fig. 4) being constrained by geology information based on a geological map by Trouw et al. (2013) (Fig. 2) and interpretations as previously exposed in the geological framework section.Data from CRUST1 model were used to constrain Moho discontinuity depth.Three geological sections (Fig. 4) have been modeled with reference density information from compilations by Carmichael (1989), Telford et al. (1990) and Turcotte and Schubert (2002) and a paper by Slavec et al. (2004) -see details in Table 1.GM-SYS forward modeling platform is an iterative environment for potential field data modeling which complies with Talwani et al. (1959) and Won and Bevis (1987) proposals for anomaly calculation.

GRAVITY RESULTS
Results comprise qualitative interpretations over Bouguer Anomaly maps and forward models in the following sections unveiling the geometry of the structural interference zone.

Bouguer anomaly maps
The Bouguer Anomaly (Fig. 4) varies from -125.8 to -59.6 mGal in the area, showing complex outlines due to different density assemblages and the superposition of the structural grains.Gravity domains trends NE-SW, NW-SE and N-S similar to mapped geology, despite an immediate relationship is not straightforward, with regions overlapping different tectonic elements.The results are discussed in terms of its possible geological counterpart according to the geological background map (Fig. 2) by Trouw et al. (2013).
Higher Bouguer anomaly values in the area are up to -59 mGal and occurs in wide, E-W to roughly NE-SW trending features in the central part part of the study area (domains E and F, Fig. 4) interpreted as high-to mediumgrade metamorphics such as the Socorro-Guaxupé Nappe granulites amidst relatively less dense areas (domain C, Fig. 4).Other higher Bouguer anomaly regions along Paraná Sedimentary Basin border are considered as basic rocks (domain B, Fig. 4), or high grade nappes over São Francisco Craton or greenstone belt associations (domain H, Fig. 4).
Gravity Domain C is given by values lower than -90 mGal trending NE-SW to NNE-SSW and E-W directions corresponding to flat-lying do medium-angle, lower grade meta-sedimentary rocks from Andrelândia and Carrancas nappe systems and Lima Duarte Nappe, also ortho-derived lower crust from Socorro-Guaxupé nappe, including syn-to post-tectonic granitic bodies.Amphibolite facies Rhyacian orthogneiss-migmatite complexes and Archean gneisses from Amparo Complex and Mantiqueira Gneiss comprises this domain, with low Bouguer anomaly values along Ribeira do Sul shear belt further southeast.
Bouguer anomaly domains with values less than -95 mGal, and circular to oval outlines (domains D and G, Fig. 4) relates to Cretaceous alkaline massifs, respectively Poços de Caldas and Itatiaia-Passa Quatro.These gravity domains reaches the minimum values observed (-125.87mGal) and its response is given by its lower density relative to a denser basement.
The domains with higher values in Bouguer Anomaly (E and F) shows gentle signal attenuation according to the expected dip of the thrust sheets, with the signal attenuating to the SW in E and to the SE in F. In domain C there is not clear evidence for the dip of the rock packages, with high gradients in NE-SW directions, suggesting that they have near-vertical geometry, which is evidenced by geological information regarding its position centered in the shear belt system.Two gravity domains with Simple Bouguer Anomaly ranging from -95.0 to -80 mGal occur in the western flank of the area relates to Paraná Sedimentary Basin (A and B, Fig. 4).Domain A represents a region of thicker sedimentary section (from basal Itararé Group to the Pirambóia Formation) underlain by basic volcanics, meanwhile Domain B complies with a thinner sedimentary pile and preponderant response from the basic rocks in Serra Geral Formation.An alternative for the higher gravity values in domain B is the presence of relatively shallow high-grade basement near to the basin border.The latter hypothesis seems feasible as the B domain extends from the basin to the North-Northeast into the Neoproterozoic metamorphic basement.
Considering Bouguer Anomaly along the structural interference zone, domains C, E and F gives information as follows: ■ domain E (fig.4) presents anomaly values up to -59.6 mGal in a E-W trend representing granulitic rocks from Guaxupé Nappe, juxtaposed to Andrelândia Nappe System granulite-facies metamorphic rocks along Campo do Meio shear belt.Similar to E domain, F domain presents higher values; however, it presents a NE-SW trend, interpreted as granulite-facies metamorphic rocks from Socorro Nappe; ■ a set of NE-SW gravity lows occurs along C domain (Fig. 4) trending parallel to Paraíba do Sul shear belt and the structural interference zone limit proposed by Ebert andHasui (1998) andTrouw et al. (2013), being considered as expression of highly-deformed/stretched metasedimentary rock packages; ■ connecting E and F domains exists a relatively higher density zone trending N-S, varying from -90 to -80 mGal without surface expression of expectedly relatively denser rock types, being interpreted as rocks akin to Socorro and Guaxupé nappe lobes amid the intense shearing.

Forward models
Forward modeling took place on three sections along the nappe system main trace (section 1, Fig. 5) and also orthogonal to regional shear belts (sections 2 and 3, Figs.6 and  7) to depict the two structural patterns.
Despite that gravity data processing did not focused Moho depth determination, the gravity long wavelength signal/forms (regional field) complies with a 37 to 43 km seismic-modeled Moho, (CRUST1 and other models).This interface is deeper to the West beneath Paraná Basin and thins to East and Southeast in a sharp step close to São Francisco Craton.Models also shows a Moho surface given by wide arcuate segments with local pinch-up points beneath the Neoproterozoic fold and shear belt, contrasting with a linear outline bellow the Archean craton (Figs. 5 and 6).The first order structural trace is given by discontinuous, medium-to low-angle thrust sheets over São Francisco Craton (sections 1 and 3, Figs. 5 and 7).These thrust sheets presents flat-lying geometry with internal stack truncation, as well as rootless and folded sheets, like Andrelândia and Carrancas nappe systems (profiles 2 and 3, Figs.6 and 7).Steepening of the thrust stacks along the Late-Neoproterozoic shear zones affects the tectonic pile from basement to the uppermost nappes, including stack truncation and affecting the Rhyacian orthogneisses and migmatites, without significant imprint in São Francisco Craton (Figs. 6 and 7).
Granulite-diatexite-metatexite lower crust associations occurs as thick low-angle nappe stacks (profiles 1 to 3, Figs. 5 to 7) representing the higher density domains (E and F) on Bouguer anomaly maps, with a discontinuous outline.These lower crust stacks are part of Alterosa Suture zone (Figs. 5 to 7), which is juxtaposed to upper level meta-sedimentary and meta-volcanic thrust sheets from Andrelândia nappe system, and above the Guaxupé Nappe granulite slice.The granulitic nappes in profiles 1 and 3 (Figs.5 and 7) are assigned density 3310 kg.m -3 for being described as containing eclogite levels as well.
Deformation along the shear belt incorporates older gneisses amid the nappe system, causing displacement and block isolation in an intricate array along Paraíba do Sul shear belt central displacement zone (e.g.Amparo Complex, Alagoa migmatite and Rhyacian gneisses, Figs. 6 and 7).
São Francisco craton was modeled as a straight to slightly curved block bellow the Neoproterozoic fold and thrust belts (Figs. 5 and 7) in a flat ramp-like geometry, with crust duplication of overriding lower crust nappes from Socorro Guaxupé nappe on high-grade gneisses from its infrastructure (Figs. 5 to 7).
Syn-to post-tectonic granites occur amid thrust sheets with bounding surfaces following thrust stacks limits, secondary thrust surfaces and shear zones, including rootless bodies (Figs. 5 and 6), hosted along past higher permeability and stress relief zones.
Paraná sedimentary basin occurs as a thin cover (Fig. 5, west flank), being modeled with a 2500 kg.m -3 density, with increasing thickness to the West according to overall basin geology (Zalan et al. 1990) emplaced over high-grade associations from Guaxupé Nappe.
Poços de Caldas Cretaceous alkaline massifs appears on section 2, shaped as an voluminous irregular circular-like body and its geometry is dubious in absence of constraining evidence.

DISCUSSION
Discussion on the crustal framework and structural interference response on geophysical data follows gravity interpretations in the light of current geological models, including comparison to previous information, data from CRUST1 and EMAG2 models, and geotectonic models.

Gravity and crustal framework
Power spectrum analysis results are in the range of crustal thickness calculated via other methods, being shallower (24 km) than estimates in CRUST1 model ranging from 37 to 42 km, similar to other independent studies (França & Assumpção 2004, Assumpção et al. 2013a, Chulick et al. 2013, van der Meijde et al. 2013).Concerning the crust stabilization age (Durrheim andMooney 1991, 1994) the Moho depth for the crust formed in Archean-and Proterozoic-ranges from 32 to 45 km on the first, and 42 to 55 km in the second.Overall lower Moho depth values are observed in the area in comparison to these reference values.São Francisco Craton crustal thickness on forward models is up to 30 km, which is inferior to mean values for crust stabilized from 2.5 to 2.6 Ga (Abbott et al., 2013).The overall trace of the geophysical domains is compatible to the structural trace in shear zones and axial zones of major folds (Fig. 8).Gravity domains E and F relates to Alterosa Suture (Fig. 8) as an discontinuous high-density crust fragment taken as the limit for Brasília fold belt and São Francisco Craton (Alkmim et al. 1993, Malagutti Filho et al. 1996), and displaced by the sinistral-motion of Paraíba do Sul shear belt by the end of the Neoproterozoic (Ebert & Hasui 1998).The link between the domains E and F, and, moreover, the Alterosa Suture are missing amid the shear system.The models constructed here indicate that Alterosa suture zone traces needs re-observation due to its structural position just below the granulite sets (bellow domains E and F).For suggesting the new outline of the Alterosa suture zone we followed the boundaries of domains E (Guaxupé Nappe) and F (Socorro Nappe).We are aware that the scale of the gravity survey and lithology maps used in our interpretation, also the density contrast of the lithologies, can differ from mapped surface geology.This suture trace is given by geophysical information and has been constrained by geological facts as long as possible.
A crustal duplex structure along the lower crust nappes and upper nappes has been proposed previously by Haralyi and Hasui (1982) Neto & Caby 1999, Iyer et al. 1996, Lama et al. 1998).Complex geometries for both thrust-folding and shear systems observed in forward models, structural and metamorphism features are expected in a continuum crustal model for orogen evolution (Boyer & Elliott 1982, Sanderson 1984, Sylvester 1988).Syn-to late-orogenic granite bodies and massifs presents shallow roots along thrust and shear-related surfaces, confirming its relation to local magma generation and emplacement during orogeny late-stages (Ebert et al. 1996, Sial et al. 1999, Mora et al. 2014).Relationships to the different crustal levels observed and tentatively modeled in the gravity datasets are further discussed in the following section.
Survey spacing, potential field ambiguity, source body shape and superposition, along with poor density contrast are probably involved in absence of clear geophysical signatures for tectonic elements, except for high-grade metamorphic rocks and Cretaceous alkaline massifs.

Deep geophysical expression of the structural interference zone
Here we discuss the proposed outline of the Alterosa Suture zone and modeled gravity sections in the light of the Moho discontinuity (CRUST1) and magnetic anomaly (EMAG2) maps.The discontinuous Alterosa Suture outline previously defined by Malagutti Filho et al. (1996) was redefined, with the new outline presented being situated near a crustal thinning region in crustal thickness maps, following the trace of the Moho depth contours (Fig. 9).The relationship to crustal thickness suggests that late-Neoproterozoic evolution included deeper crust re-organization or that such level of reconstruction was already established.Cenozoic alkaline magmatism could account for this discrete thinning in São Francisco Craton border but needs supporting information.The discontinuity in the trace of the suture zone along with the interference zone confirms that the lateral structural system is post-kinematic in relation to the suturing event.
EMAG2 data magnetic grain shows high-magnetization areas according to regional low-angle thrust surfaces direction and to the high-angle shear systems main trace as presented by Ebert and Hasui (1998) (Fig. 9).The wedge-like continuous high magnetization zone follows the nappe system outline, with decreasing signal to the southwest (into the Paraná Sedimentary Basin) and extending further northeast into the São Francisco Craton with subdued signal.The magnetic response in the Paraná Basin is given by extreme lower values and into the Neoproterozoic regions in a cuspid shape, also beneath the Poços de Caldas alkaline massif.
Interestingly, the magnetic grain near the proposed suture zone trace (Fig. 9) is given by gentle signal decay to the SW in the Guaxupé Nappe lobe, suggesting a dip to that direction, as also argued from structural evidences, the gravity anomaly and the gravity models.For the Socorro lobe things are not that straightforward, with a general resemblance from the magnetic outline to the suture line, with more pronounced gradients in the magnetic anomaly, suggesting a dip angle not so gentle as in Guaxupé, probably due to the major influence of the Paraíba do Sul Shear Belt and the structural interference zone.
It is noteworthy that the structural interference zone limit proposed by Trouw et al. (2013) fairly centers on the high magnetization wedge-like zone, suggesting that the developing the structural interference zone did not altered its magnetic properties by giving a magnetic signature by itself.This relies on the fact that the EMAG2 model in the area is given mainly by satellite information, representing long-wavelength information -so it depicts deep crust.
Discontinuity from the high-grade nappes and São Francisco craton further north is not evident on magnetic data.This possibly relates to poor magnetic contrast between the two high-grade terrains, as they comprise exposed lower crust sections commonly with high magnetization and long-wavelength signal (Shive et al. 1992).These long wavelength characteristics are products of deep-seated continental roots with its common mafic nature, rich in magnetic minerals, also Moho, upper-lower crust and Curie interface undulations, or higher metamorphic grade of deeply-eroded crust (Shive et al. 1992).
The new trace for the Alterosa suture is more conformable to the gravity, magnetic and Moho depth information, and has been conceived with geologic constraints, differing from the outline by Malagutti Filho et al. (1996).It also has a clear oblique geometry to the structural interference zone limit defined by Trouw et al.
(2013) (Fig. 9).Regarding EMAG2 information a total match from the trace of the features mapped in gravity and surface geology is reasonably unexpectable given the satellite nature of the magnetic information, and presents deep crust information which implies a possible mismatch regarding dipping structures.
To define a precise structural interference zone limit by geophysical methods is challenging as it concerns a complex structural framework with weak density contrast between entities.Surface information suggests that structural interference goes further North over Socorro-Guaxupé Nappe (Peternel et al. 2005, Trouw et al. 2013), which is supported by gravity data as Alterosa suture discontinuity is fairly coincident with the proposed structural interference zone limit marked by Três Corações shear zone on Paraíba do Sul shear belt (Malagutti Filho et al. 1996, Ebert & Hasui 1998).

Insights into geotectonic models
Debate over crust-mantle interactions and its linkage to orogenesis since Archean has poor evidence for lower crust development (Windley & Tarney 1986, Percival et al. 1992).However, present models on mantle thermal evolution and mass transfer points out to a dynamic scenario from early Earth to nowadays tectonics with different evolution of collisional systems along time, markedly post-Proterozoic orogens (Cawood et al. 2009, Abbott et al. 2013, Gerya 2014, Jamieson & Beaumont 2013, Kamber 2015).Lower-and upper-crust sections juxtaposition by thrust systems supports its evolution by thick-skinned tectonics style (Coward 1983).This contrasts to the observation that thrust systems are common in orogens involving shallower crust conditions (thin-skinned tectonics) as observed on Brasilia Fold Belt further north from study area (Simões 1995, Strieder & Suita 1999, Valeriano et al. 2008).Collision tectonics models supports late-colisional continental block adjustment in large-scale shear belts, overprinting prior structural grain (Cawood et al. 2009, Gerya 2014, Jamieson & Beaumont 2013).
Comparing the study area based on a geodynamic model for accretionary orogens by Jamieson and Beaumont (2013) shows similarity with large-hot orogen systems, with distinct periods of crustal material accretion involving sutures, alternating with stages of thrusting of deep crust plutonic and high grade metamorphic rocks.According to the author's Prowedge/Retrowedge/Uplifted plug/subduction Conduit (P-U-R-C) systematics, Brasília and Ribeira fold belts lies between São Francisco, and Paranapanema cratons and related African cratons (Congo) -similar to a PUR system, as follows.The closure of the former ocean took place with the fragment containing the Prowedge (P) overriding plate forming thrust stacks (Ediacaran to lower crust ultra-high pressure Cryogenian nappe systems) and an internal and less deformed Retrowedge (R) with its subduction Conduit (C) far from the collision front, probably not outcropping today.In-between P and R segments lies an extruded bivergent Uplifted plug (U) of highly deformed, upper and lower crust nappes with the meta-sediments from the former oceanic(?)basins and interlayered granitoid bodies and extruded granulites.Late-stage Early Cambrian tectonics took place with PUR assemblage deformation by Paraíba do Sul and Campo do Meio shear belts including granitic magma intrusion along structural surfaces, with ongoing sin-to post tectonic granitic magmatism.Nowadays deep erosion levels (probably) precludes entire large-hot orogen system recognition by suppression of uppermost orogen architecture, as well as most of Uplifted plug and Retrowedge being mostly buried bellow Paraná Sedimentary Basin and demands further investigation.

CONCLUDING REMARKS
Previous geophysical crustal models (Haralyi & Hasui 1982, Malagutti Filho et al. 1996) stands confirmed for crustal duplication and counts now with revised and detailed assessment, although Alterosa suture outline in southeast-to east regions still demands further scrutiny with field information and alternative geophysical approaches.Models proposed here also favors previously envisaged late-stage shearing tectonic models (Morales 1993, Ebert & Hasui 1998), seemingly extendable to deeper levels (Boyer & Elliott 1982, Sanderson 1984, Sylvester 1988).Shear-related structures superposes the low-lying to medium-inclined thrust stacks and controls syn-to late-tectonic granitic bodies emplacement (Ebert et al. 1996, Ebert & Hasui 1998, Mora et al. 2014).The extension of the structural interference zone proposed by Trouw et al. (2013) stands confirmed by geophysical means, also observations by Campanha and Brito Neves (2004) are somewhat in agreement regarding the thrust system formation, and the oblique kinematics for its formation.
Thrust stacks are extended to deep levels and external/ independent information is necessary to reduce subjectivity and ambiguity inherent to the gravity interpretation (Skeels 1947, Jessell & Jessell 2001).Gravity itself does not address entirely the structural complexity and joint interpretation.The CRUST1 model presents evidence for crustal-scale structure along the fold and thrust system limit along São Francisco Craton south/south-western edge with similar interpretation on EMAG2 data.Since complete understanding is far from complete, our proposed analogue model for accretion-related evolution takes the syntax zone as the product of a two-fold event, accounting for late-shearing structural imprint.This analogue model takes basis in Jamieson and Beaumont (2013) PURC systematic for accretion since Proterozoic times (at least) in disagreement to pure fold and thrust belt proposals (Campos Neto et al. 2011, Trouw et al. 2013).
From the geophysical perspective the outline of the Alterosa Suture zone still demands further scrutiny allied by field work and, maybe, by other geophysical surveys.This work presents an alternative to current views and stands to be proven by information other than geophysics -relying in the supportive nature of the geophysical datasets -absent of trueness by itself.

Figure 2 .
Figure 2. Tectono-stratigraphic units along southern Brasília thrust and fold belt with study area outline (dashed line).Note that 'alk' stands for Cretaceous alkaline massifs.Modified from Trouw et al. (2013).Black dots for cities, black dashed line for the structural interference zone according to Trouw et al. (2013), long dashed black line for the suture zone outline by Malagutti Filho et al. (1996) and black continuous outline for the study area.

Figure 3 .
Figure 3. Gravity data station location over SRTM digital elevation model.Black dashed line for the structural interference zone according to Trouw et al. (2013).

Figure 4 .Figure 6 .
Figure4.Bouguer Anomaly map and gravity domains in this study (dashed black lines) and modeled sections trace (dashed white lines, from 1 to 3).

Figure 5 .
Figure 5. Forward model 1-1' along the thrust system over São Francisco Craton (see figure 4 for model location).Gravity domain limits are marked in lower level.Vertical exaggeration is 0.5.

Figure 7 .Figure 8 .
Figure 7. Forward model 3-3' along the thrust system over São Francisco Craton (see figure 4 for model location).Gravity domain limits are marked in lower level.Vertical exaggeration is 0.8.

Figure 9 .
Figure 9. Structural and geophysical elements around Alterosa suture zone region.
suture outline (this work) Structural interference zone limit by Trouw et al. (2013) Moho depth estimates (in km) on

Table 1 .
Density values for the rock unities in the modeled sections.