Mapping New IOCG Mineral Systems in Brazil: The Vale do Curaç á and Riacho do Pontal Copper Districts

: The Vale do Curaç á and Riacho do Pontal copper districts are located within the northern part of the Archaean S ã o Francisco Craton and represent two pulses of mineralization. The copper districts have been identiﬁed as Iron-Oxide-Copper-Gold (IOCG) classes of deposits. An older metallogenic event associated with the Cara í ba copper deposit, which is located in the Vale do Curaç á district, is related to Palaeoproterozoic (ca. 2 to 2.2 Ga) hydrothermal processes. A younger Neoproterozoic (ca. 750 to 570 Ma) episode of volcanism and associated plutonism is represented by the Riacho do Pontal mineral district. Seismic tomography data from across east-central Brazil show that the multiage Caraj á s province and Vale do Curaç á and Riacho do Pontal copper districts sit along either side of a prominent NW-trending upper lithospheric high-velocity zone. The edges of the high-velocity zone point to long-lived subparallel transcrustal structures that have been the focus of multiple reactivations and copper mineralization events. Regional gravity and magnetic maps show that the Vale do Curaç á copper district extends over an area greater than 110 km by 22 km. The magnetic and gravity values show signiﬁcant variations correlated with this area. The district includes high gravity values associated with the Cara í ba copper mine ( > − 35 mGal), which has a greater density (3.13 g / cm 3 ) than the nonmineralized host rock density (2.98 g / cm 3 ). The gravity anomaly signature over the Riacho do Pontal copper district is characterized by a 40-km long NW–SE trending Bouguer gravity low. The Ria4 occurrences of the Riacho do Pontal copper district are situated in these regional low-gravity domains. Data from regional airborne magnetic and ground gravity surveys were inverted to obtain a 3D magnetic susceptibility and density model, respectively, for the known districts. The results show that the Cara í ba deposit is characterized by a both dense and magnetic source showing structural control by thrust shear zones. The 2D and 3D geological models show two main NNW prospective trends. Trends I and II have a sigmoidal shear shape and are positioned in the contact zone between domains with high magnetic susceptibility (SI > 0.005) and density > 0 g / cm 3 ). Trend I is 40 km × 10 km in size and hosts the Cara í ba, Surubim, and Vermelho copper mines and other minor deposits. The results obtained from the 3D magnetic inversion model for the region of the Riacho do Pontal district show weak magnetic anomaly highs extending along a NW–SE magnetic gradient trend. The gradient is related to mapped shear zones that overprint older and deeper NE–SW features of the S ã o Francisco cratonic root. The area includes high gravity values associated with the Cara í ba copper deposit, which has a greater density (3.13 g / cm 3 ) than the nonmineralized host rock density (2.4

magnetic and gravity data are inverted to obtain 3D models of magnetic susceptibility and density, respectively, to understand the concealed geology related to physical property distributions and geometries underlying the districts. Third, seismic tomography data are analyzed to place the Vale do Curaçá and Riacho do Pontal IOCG districts within a tectonic and cratonic framework and to characterize potential lithospheric controls of the mineral provinces. Our objective is to use these approaches to obtain an advanced understanding of the mineral system(s) that controlled the emplacement of the NE Brazil IOCG mineral districts.  [20,27,31].  [20,27,31].

Purpose and Scope
The purpose of our study is to synthesize available regional magnetic, gravity, and seismic tomography data with geologic information to evaluate the 2D and 3D geologic and tectonic setting of the Vale do Curaçá and Riacho do Pontal IOCG districts. An overarching goal of our study is to develop a better understanding of the deep crustal and lithospheric structural controls for Palaeoproterozoic and Neoproterozoic mineralization events.

Classification as IOCG Deposits
Early studies of the Caraíba copper mine within the Vale do Curaçá district emphasized orthomagmatic processes as responsible for the genesis of copper mineralization [13,14,30]. noted that the unusual aspects of the Caraíba deposit shared similarities with other copper deposits of magmatic origin [33], including (a) the presence of primary sulfides such as bornite and chalcopyrite; (b) the presence of large amounts (>50 wt %) of magnetite; (c) a high Cu/Ni ratio (~40); and (d) orthopyroxenites with abundant biotite related to shear zones. The importance of metasomatic processes in the genesis of the Caraíba deposit was initially suggested by Rocha et al. [34] and more recently confirmed by Teixeira et al. [20], who proposed that copper deposits in the Vale do Curaçá copper district are analogous to other IOCG deposits worldwide.
Early studies fulfilled by several authors [35,36] on the copper deposits in the Riacho do Pontal district show that they were characterized as the volcanogenic massive sulphide (VMS) type. In the region of Riacho Seco. The following criteria have been used to classify copper mineralization in the Riacho do Pontal district as IOCG deposits [24,28]: (a) strong association with potassic and albite alteration; (b) copper ore that is hydrothermal and epigenetic, which is closely associated with shear zones; (c) the presence of hydrothermal iron oxides spatially and temporally related to copper orebodies; (d) mineralization that is not spatially related to proximal granites; and (e) magnetite or hematite (hydrothermal iron > 10%; [7,8]).

Palaeoproterozoic Vale do Curaçá Copper District
Since the beginning of the last century, the Vale do Curaçá copper district (B in Figure 1) has been the focus of intense geologic mapping and exploration activity. During the 1970s, a mineral exploration boom occurred in this region. The focus of the exploration program was on copper in volcanic-hosted massive sulfide (VMS) mineralization [37,38] in greenstone belts and copper mineralization hosted in mafic-ultramafic layered complexes [14,31,39]. The exploration activities led to the discovery of the Caraíba mine, which is the largest mine in the district. In addition, a series of copper deposits and occurrences in the Vale do Curaçá district were discovered. Approximately 300,000 m of drilling has taken place in the mining district [40].
The geology of the Vale do Curaçá copper district (B in Figure 1) consists mainly of rocks belonging to the Itabuna-Salvador-Curaçá Orogen (ISCO). The ISCO is located east of the Gavião and Jequié blocks and west of the Serrinha block ( Figure 2) and represents an 800 km long belt of Archaean to Palaeoproterozoic high metamorphic grade granulitic rocks [20,41]. These rocks were deformed and metamorphosed due to the collision of the Gavião, Jequié, and Serrinha blocks [20,[41][42][43]. The ISCO is a Rhyacian to Orosirian orogenic belt resulting from the collision of Archaean blocks in the range of ca. 2 to 2.2 Ga [26,44]. This was followed by deformation with high metamorphic grade and syntectonic granites [45]. A SHRIMP U/Pb age of 2695 ± 12 Ma was obtained for granulitic rocks of the Caraíba complex [20] and a U/Pb age of 2072 Ma was found for the metamorphic event [45].
The mafic-ultramafic complex should have mafic-ultramafic intrusions of two different ages (2695 Ma and 2580 Ma). The specific geology of the Curaçá copper district comprises the Tanque Novo and Caraíba complexes. The Tanque Novo complex consists of Bom Despacho, Arapuá, Surubim, and Bogó gneissic rocks [13] composed mainly of feldspathic quartz gneiss, cordierite-garnet-biotite gneiss, amphibolite, and magnetite (quartzite)-rich rocks. The Caraíba complex contains tonalite to granodiorite bodies. A series of mafic-ultramafic bodies individually intruded into the supracrustal rocks. These rocks have undergone episodes of deformation, metamorphism, and related granitic intrusion. The U/Pb SHRIMP ages of ultramafic rocks within the Caraíba mine show that the rocks formed at ca 2580 ± 10 Ma [45]. Field relationships show that the Itiúba syenite is younger than the Medrado gabbro. Importantly, Itiúba syenite is a syn tectonic to shear zone of 2084 ± 9 Ma. This age is also related to copper mineralization situated along shear zones [26,46]. Shear zones have been mapped along the Vale do Curaçá district and are often marked by intense biotite-rich rocks related to mylonitic rocks. These zones are characterized by two phases of progressive deformation [34,46]. (a) The first phase (D1) involved transpressive thrust structures showing convergence to the west. The associated metamorphism occurred in the amphibolite to granulite facies. (b) The D2 phase was related to strike-slip shear zones with N-S isoclinal folds and penetrative vertical shearing structures. Kinematic studies indicate that D1 was overprinted by sinistral strike-slip shear zones [26,47]. These shear zones were responsible for the positioning of mylonitized granites and the intense fluid flux. The shearing process destroyed the preexisting fabric. The foliation of shear zones has a strike of 350° and dips towards SE at 45° to 75°. The main mineralization event in the Vale do Curaçá district occurred between ca 2050 and 2010 Ma [26,29]. The specific geology of the Curaçá copper district comprises the Tanque Novo and Caraíba complexes. The Tanque Novo complex consists of Bom Despacho, Arapuá, Surubim, and Bogó gneissic rocks [13] composed mainly of feldspathic quartz gneiss, cordierite-garnet-biotite gneiss, amphibolite, and magnetite (quartzite)-rich rocks. The Caraíba complex contains tonalite to granodiorite bodies. A series of mafic-ultramafic bodies individually intruded into the supracrustal rocks. These rocks have undergone episodes of deformation, metamorphism, and related granitic intrusion. The U/Pb SHRIMP ages of ultramafic rocks within the Caraíba mine show that the rocks formed at ca 2580 ± 10 Ma [45]. Field relationships show that the Itiúba syenite is younger than the Medrado gabbro. Importantly, Itiúba syenite is a syn tectonic to shear zone of 2084 ± 9 Ma. This age is also related to copper mineralization situated along shear zones [26,46]. Shear zones have been mapped along the Vale do Curaçá district and are often marked by intense biotite-rich rocks related to mylonitic rocks. These zones are characterized by two phases of progressive deformation [34,46]. (a) The first phase (D1) involved transpressive thrust structures showing convergence to the west. The associated metamorphism occurred in the amphibolite to granulite facies. (b) The D2 phase was related to strike-slip shear zones with N-S isoclinal folds and penetrative vertical shearing structures. Kinematic studies indicate that D1 was overprinted by sinistral strike-slip shear zones [26,47]. These shear zones were responsible for the positioning of mylonitized granites and the intense fluid flux. The shearing process destroyed the preexisting fabric. The foliation of shear zones has a strike of 350 • and dips towards SE at 45 • to 75 • . The main mineralization event in the Vale do Curaçá district occurred between ca 2050 and 2010 Ma [26,29].

Geology of the Caraíba Mine
The geology of the Caraíba mine is well marked by the contact between orthogneiss (G2) and mafic-ultramafic rocks. Frequently, biotite has been mapped along N-S trending geologic contacts. The supracrustal rocks include banded gneiss, diopsidite, forsterite marbles, garnetiferous mafic granulite, biotite-schist, and iron formations [14].
Ore minerals at Caraíba include chalcopyrite, bornite, and chalcocite, which are usually associated with magnetite. Intercumulus crystals of disseminated chalcopyrite, bornite, and magnetite grown in the interstices of pyroxenes and amphiboles constitute the primary (magmatic) mineralization. Minerals associated with hydrothermal alteration include biotite, microcline, epidote, chlorite, magnetite, and quartz [36]. Ore often occurs as veining or magnetite hydrothermal breccia. Magnetite sometimes occurs in equilibrium with spinel. Frequently, ilmenite exhibits replacement by chalcopyrite, bornite, and magnetite, especially when the host rock is biotite mylonite. Magmatic ilmenite is replaced by ilmenite hydrothermal. The latter IOCG system is richer in copper, highly magnetic, and dense.

Geology of the Neoproterozoic Riacho do Pontal Copper District
More than twenty copper occurrences have been found in the Riacho do Pontal region (A in Figure 1; [26]). Geologically, these occurrences are mostly located in the contact zone between the Borborema province and the San Francisco craton (Figures 2 and 3). The area includes a small slab of the ca. 750 Ma to 570 Ma Brasiliano collisional zone, which developed during the convergence of Neo-to Mesoproterozoic terrains with the São Francisco craton. Copper occurrences are structurally controlled by shear zones (Figure 1). The Ria4 occurrences and Riacho Seco deposit ( Figure 1) are the main copper occurrences in the district. The drill hole intercepted a mineralized zone (0-26 m) related to biotite-malachite (after cpy; 1-3%) parallel to mylonitic foliation.
Several alteration types related to hydrothermal processes have been identified in the Riacho do Pontal district [28]. Hydrothermal alteration has led to a pervasive calcic-potassic and calcic-ferric overprint of the host gneisses and migmatites. Early-stage sodic hydrothermal alteration is distal in relation to areas with calcic-potassic and potassic-ferric alteration with a higher strain rate. Hydrothermal alteration processes led to the replacement of metamorphic minerals within banded gneiss with a suite of minerals, including albite, biotite, hematite, and amphibole. These replacement minerals are spatially connected to penetrative shear deformation, and ore zones plunge parallel to stretching lineations along shear zones. Pyrite, chalcopyrite, and chalcocite are related to the main sulfide minerals found in ore zones. Magnetite and/or hematite typically occurs in association with copper ore minerals and constitutes less than 1% of the paragenesis [29].
The Riacho do Pontal belt was overprinted by two main shearing events: (a) D1-thrust shear zones with convergence to the north representing a series of northwest-trending parallel structures related to hydrothermal alteration and mineralization and (b) D2-dextral strike slip shear zones related to multiple reactivation responsible for the strong shearing of paragneissic and orthogneissic rocks. The foliation of shear zones is significantly tilted with a strike of 320 • and a dip to the NE at 60 • to 75 • . Both events were related to transpressive structural regimes.
The Riacho do Pontal Copper District is located at the apex of three orogenic belts (the Riacho Preto, Riacho do Pontal, and Sergipana belts), which are located along the northern border of the São Francisco Craton ( Figure 2). The area has been located within the Brasiliano collisional zone (from ca. 750 Ma to 570 Ma), which developed during the convergence of the Neo-to Mesoproterozoic terrains of Borborema Province (Figure 3) with the São Francisco Craton [49]. The suture is marked by a gravity anomaly heavily flanked by anomalous lows (Figure 4b). The positive anomaly corresponds to the lower crustal uplift of the Borborema province, and the flanking negative anomalies correspond to lower density nappes pushed towards the craton [50,51]. According to [50], supracrustal rocks of the Riacho do Pontal and Sergipana belts were placed over the São Francisco Craton with displacements to the south of the order of 30-60 km [52]. This overlapping structure is observed in both magnetic and gravimetric data with the truncation and displacement of old N-S structures occurring along the border between the São Francisco Craton and Borborema terrain. According to [52], the collisional event was overprinted by transcurrent shear zones.
Minerals 2020, 10, x FOR PEER REVIEW 7 of 30 supracrustal rocks of the Riacho do Pontal and Sergipana belts were placed over the São Francisco Craton with displacements to the south of the order of 30-60 km [52]. This overlapping structure is observed in both magnetic and gravimetric data with the truncation and displacement of old N-S structures occurring along the border between the São Francisco Craton and Borborema terrain. According to [52], the collisional event was overprinted by transcurrent shear zones.

Tectonic Setting
IOCG provinces and districts in Brazil and Africa are mainly situated along the borders of major crustal blocks and in Neoproterozoic orogenic belts ( Figure 3). There are at least three IOCG metallogenic provinces with different ages in Brazil, which include the Carajás province (ca 2.72 Ga and ca 2.5 Ga), the Vale do Curaçá district (ca 2 Ga), and the Neoproterozoic (ca 750 Ma-570 Ma) districts of Riacho do Pontal.

Tectonic Setting
IOCG provinces and districts in Brazil and Africa are mainly situated along the borders of major crustal blocks and in Neoproterozoic orogenic belts ( Figure 3). There are at least three IOCG metallogenic provinces with different ages in Brazil, which include the Carajás province (ca 2.72 Ga and ca 2.5 Ga), the Vale do Curaçá district (ca 2 Ga), and the Neoproterozoic (ca 750 Ma-570 Ma) districts of Riacho do Pontal.

Gravity and Magnetic Data
Gravity and magnetic anomaly data over an area of 7500 km 2 were extracted from public databases available from the Brazilian Geological Survey (CPRM) and the gravity data reduced to Bouguer gravity anomalies using procedures described in [50]. A total of 60 gravity points spaced an average of 8 to 15 km distance were extracted from the gravity database and used for this study (Figure 4a). Given the coarse spacing of the gravity data, we expect crustal density sources that occur at depths of 5 km and beyond to be effectively resolved. The gravity data were gridded by applying a minimum curvature algorithm to a 500 m grid (Figure 4b). The gravimetric signature is quite distinct in both districts while the Vale do Curaça Copper District is associated with a higher gravimetric signature, and the Riacho do Pontal Copper District is associated with a lower valued gravimetric domain (Figure 4b).
Airborne magnetic anomaly data were acquired in 2001 as part of the Riacho Seco and Andorinhas Project executed by the Brazilian Geological Survey [58]. Data were collected at a nominal height above ground of 100 m along N-S flight lines spaced 500 m with E-W tie lines every 10,000 m. A total of 48,641 km of magnetic data were collected. The total magnetic field (TMI) flight line data were gridded onto a 250 m grid ( Figure 5). Given the flight-line spacing of the magnetic surveys and height at which the magnetic field measurements were observed, we expect to resolve crustal magnetic sources that lie at depths of 400 m and deeper [59]. The horizontal resolution is 200 m.

Gravity and Magnetic Data
Gravity and magnetic anomaly data over an area of 7500 km² were extracted from public databases available from the Brazilian Geological Survey (CPRM) and the gravity data reduced to Bouguer gravity anomalies using procedures described in [50]. A total of 60 gravity points spaced an average of 8 to 15 km distance were extracted from the gravity database and used for this study (Figure 4a). Given the coarse spacing of the gravity data, we expect crustal density sources that occur at depths of 5 km and beyond to be effectively resolved. The gravity data were gridded by applying a minimum curvature algorithm to a 500 m grid (Figure 4b). The gravimetric signature is quite distinct in both districts while the Vale do Curaça Copper District is associated with a higher gravimetric signature, and the Riacho do Pontal Copper District is associated with a lower valued gravimetric domain (Figure 4b).
Airborne magnetic anomaly data were acquired in 2001 as part of the Riacho Seco and Andorinhas Project executed by the Brazilian Geological Survey [58]. Data were collected at a nominal height above ground of 100 m along N-S flight lines spaced 500 m with E-W tie lines every 10,000 m. A total of 48,641 km of magnetic data were collected. The total magnetic field (TMI) flight line data were gridded onto a 250 m grid ( Figure 5). Given the flight-line spacing of the magnetic surveys and height at which the magnetic field measurements were observed, we expect to resolve crustal magnetic sources that lie at depths of 400 m and deeper [59]. The horizontal resolution is 200 m.

Gradient Functionals Calculated from the Magnetic Data
Several gradient quantities of the magnetic field were calculated from the grid of the TMI to enhance magnetic variations and lineations related to geologic structures, map magnetic signatures of altered and unaltered rocks, and identify magnetic signatures related to deep magnetic sources. The TMI data were reduced to the magnetic north pole (RTP) using an inclination of −20° and declination of −22.31° [60]. The RTP grid was continued upward to 1000 m above the terrain to reduce

Gradient Functionals Calculated from the Magnetic Data
Several gradient quantities of the magnetic field were calculated from the grid of the TMI to enhance magnetic variations and lineations related to geologic structures, map magnetic signatures of altered and unaltered rocks, and identify magnetic signatures related to deep magnetic sources. The TMI data were reduced to the magnetic north pole (RTP) using an inclination of −20 • and declination of −22.31 • [60]. The RTP grid was continued upward to 1000 m above the terrain to reduce noise caused by the low geomagnetic latitude and to enhance the magnetic signatures of deeper sources ( Figure 5).
The grid of upward-continued RTP data was used to calculate the tilt angle of the total horizontal gradient (TAHG; [60]). The TAHG transform clearly identifies tracks of high anomaly gradients that map magnetic contrasts related to geologic structures at a range of depths ( Figure 6a). Tracks of the TAHG ridgelines, marking major geologic structures, are shown relative to the gravity anomaly map (Figure 6b). The TAHG effectively identifies changes in the prevailing geologic structures as the northern limit of the SFC, where the main strike of the gradient changes from sub-NS orientations to SSW-NNE directions on the SFC and to SE-NW oriented strikes north of the SFC along the Sergipano foldbelt. The SE-NW strike farther north is replaced by more chaotic strike directions when moving northwards out of the foldbelt. Major copper mines Ria 1 to Ria 5 are located in the well-defined SE-NW striking magnetic field anomalies. The copper mineralizations are tied to shear movements and to magnetite and hematite occurrences. It is to be expected that the magnetic anomalies are correlated with the copper mines and have shearing strike directions.
Euler deconvolution solutions [61] were calculated from the grid of the magnetic anomaly data, and depths to magnetic source (structural index = 0) are shown along major shear zones. The structural index = 0 is adequate for detecting the top of a thick slab [62].

Petrophysical Data
Magnetic susceptibility and density data were acquired from drill cores and surface rock samples. Measurements are provided as susceptibility in SI (susceptibility intensity) × 10 −3 and density in g/cm 3 . Magnetic susceptibility was measured using a KT10 PLUS S/C susceptibility meter, and rock density was measured by vacuum saturating the samples using a hydrostatic balance. During the field seasons, 33 drill core samples were investigated to determine density and susceptibility from samples collected in the Riacho do Pontal copper district, and 11 rock samples were collected from the outcrops and mines of the Vale do Curaça copper district (Figure 7). The densities found for the Riacho do Pontal copper district are representative of the rock types found in the boreholes, which are gneissic rocks (Table 1).  The densities found for the Riacho do Pontal copper district are representative of the rock types found in the boreholes, which are gneissic rocks (Table 1). The characteristic density and susceptibility values in the two districts are very different, as is the relation between the two quantities. The Riacho do Pontal copper district has quite standard crustal-rock-type densities, with values between 2.6 and 3.1 g/cm 3 . For all of the rock samples the magnetic susceptibility is proportional to density, but with two groups of proportionality. A first group is only mildly proportional, with susceptibility varying between 0 and 1.8·10 −3 SI, and the other group is highly proportional, with susceptibility varying between 0 and 13·10 −3 SI. For the second group the hydrothermal process was responsible for increase of magnetite. The Vale do Curaça copper district has much higher densities between 3.5 and 4.4 g/cm 3 , with much higher proportionally varying susceptibilities values, ranging between 0 and 800·10 −3 SI. The densities of the Vale do Curaça copper district are typical of mafic and ultramafic rocks, rather than crustal type rocks. In the area the metamorphism has overprinted the ultramafic rocks, transforming them into granulites. The IOCG hydrothermal process in the shear zones increases the magnetization through deposition of magnetite in veins associated with copper mineralization.

Geologic Data
A compilation of data from the Brazil Geological Survey (CPRM) was used as a geological base map [20]. Supplemental geologic information, including structural and image interpretation data (Landsat 8), was collected for this study along secondary roads, resulting in an additional 300 geological observations [28]. Geologic sections along selected profiles were developed to better understand the structural setting of the area.

Seismic Tomography
Seismic tomography data have been collected for an area covering east-central Brazil [63,64]. The seismic velocity model which emerged from that study were provided to us by Minerals Targeting International Pty Ltd. and were extracted from a database produced as part of the Global Lithospheric Architecture Mapping (GLAM) Project [65] and refer to the methodology used for the model of [63]. The tomographic image shows the S-wave velocity (Vs) of the Carajás IOCG province, Riacho do Pontal and Vale do Curaçá copper districts in Brazil [63]. The data were used to determine seismic properties related to the lithosphere underlying major IOCG provinces in Brazil ( Figure 8). The data represent a model of average seismic velocities at the depth slice of 0-100 km [63]. Our studied area is inboard of a relatively high velocity region, bordered eastwards by a low velocity region. The high-velocity region presumably indicates the cratonic area, with the lower velocity region demarcating the craton limit.

3D Magnetic and Gravity Inversion Models
The magnetic potential method is sensitive to the magnetic properties of the rocks in terms of residual and induced magnetization. However, applications of this practice have been limited by

3D Magnetic and Gravity Inversion Models
The magnetic potential method is sensitive to the magnetic properties of the rocks in terms of residual and induced magnetization. However, applications of this practice have been limited by conventional assumptions that rock magnetization is dominated by induced magnetization and that the magnetization direction is aligned with the geomagnetic field [66]. Modeling the magnetization vector direction and strength using magnetization vector inversion (MVI) has challenged these conventions, and MVI modeling has become an important exploration tool in the mining industry [66]. VOXI Earth Modelling is an Oasis Montaj cloud and clustered computing package that allows the inversion of geophysical data in 3D [66,67]. It uses a Cartesian Cut Cell (CCC) and the algorithm has been simplified by Ellis and MacLeod [67] in order to represent geological surfaces with better accuracy. VOXI Earth Modelling is an Oasis Montaj cloud and clustered computing package that allows the inversion of geophysical data in 3D [66,67]. It uses a Cartesian Cut Cell (CCC) and the algorithm has been simplified by Ellis and MacLeod [67] in order to represent geological surfaces with better accuracy.
Let us begin with the very general assumption that the magnetic properties of the earth can be represented by volume magnetization [68]. We make no assumptions about whether the source of this magnetization is induced or remanent [66,68].
Magnetic and gravity data were inverted to create 3D models of magnetic susceptibility and density, respectively, of selected areas of the two studied copper districts to better understand the anatomy of the crust underlying them. The three-dimensional models were integrated with other project data using a multicore cloud-based computing platform.

3D Magnetic Inversion
The grid of total field magnetic anomaly data was inverted to create a 3D magnetic susceptibility model using the MVI method developed by [67] and employing software available from the Geosoft Oasis software package [66,68]. The MVI method is a particularly effective inversion methodology to use in areas of low geomagnetic latitude and in geologic environments where little is known about the remanent magnetization of the rocks under study [66].
No palaeomagnetic studies have been conducted on the rocks within these copper districts. The method assumes that the magnetic signal observed at the surface is associated both with varying susceptibility contrasts and with the varying direction of magnetic field magnetization, which allows one to verify the presence of remanence. In the end, we conclude that remanence is not present.
Voxels for the magnetic inversion were set to 250 m (x) × 200 m (y) × 250 m (z). Model inputs included petrophysical properties, surface geology, and structural information as qualitative constraints.

3D Density Inversion
VOXI MVI modeling is a Geosoft Oasis Montaj cloud and clustered computing module used for the inversion of ground gravity data. The Geosoft package uses a cartesian cut cell (CCC) and an iterative reweighting inversion algorithm [67,68]. The algorithm has been simplified by [65,66] and represents geological surfaces with good accuracy. For the gravity models, voxels sized 250 m × 250 m × 250 m were used.
The model used topography to represent the top of the model and a 5000-m depth to constrain the bottom surface of the model. The Z voxel cell size follows a logarithmic progression with depth. The acceptable absolute error level for all models was set to 0.02 mGal. A linear trend background was removed from the input gravity grid to facilitate the modeling process and avoid generating erroneous results and undesirable edge effects. We use the Bouguer gravity field instead of the isostatic anomalies, because in the specific inversion areas, the topography is flat.

Vale do Curaçá District
Gravity anomaly signatures over the Vale do Curaçá copper district (Figure 6b) show that the district is underlain by rocks with high densities in contrast with those surrounding it. The district is located within an elongate 110 km × 22 km long N-S trending gravity high that occurs over rocks related to the Bom Despacho gneiss, Surubim gneiss, Banguê gneiss, mafic and ultramafic complex and tonalitic to granodioritic orthogneiss (Figure 1). In this study, density measurements were made from these units, including gneissic-granite, mylonitic, and pyroxenite rocks, which have average densities of 2.78, 3.1, and 3.3 g/cm 3 , respectively (Figure 7b). The lower gravity anomalies surrounding the gravity high are associated with the following units: Riacho da Onça augen gneiss and the Acauã Group. The densities measured from these units are generally lower and range from 2.65 to 3.0 g/cm 3 . The gravity low is hardly to be expected to be due to an isostatic compensation of topography, as there are little variations of topography.
The mines of Caraíba, Surubim, and Vermelhos and other copper occurrences are located in a region where rocks are widely affected by granulite-facies metamorphism, and strong hydrothermal processes which show a more positive gravimetric region situated on the eastern portion of the Vale do Curaçá.
The mines of Caraíba, Surubim, and Vermelhos and other copper occurrences are located in the eastern region of the Vale do Curaça district, which includes a relatively higher gravimetric region. Density measurements taken from rocks sampled from ore zones within the Caraíba and Surubin mines show very high average densities of 3.5 to 4.35 g/cm 3 . Given the coarse gravity station spacing (Figure 4a), we do not expect anomalies associated with the high-density ore to be resolved in the gravity anomaly map. The ore is comprised of high-density rocks. The ore includes chalcopyrite, magnetite, and bornite associated with hornblende, biotite, apatite, and zircon. Chalcocite and ilmenite are rare [13].
The Vale do Curaçá district sits within a NS magnetic high positioned alongside sets of elongate magnetic lows that record major shear structural episodes and juxtaposed rocks with different magnetic susceptibilities. The magnetic anomaly lows are mapped as gneissic rocks and mafic and felsic dikes. These rocks exhibit weak levels of hydrothermal alteration as shown by outcrop and petrographic studies. Copper mineralization in the Vale do Curaçá district is associated with linear NS-trending magnetic anomalies located along the west side of the Itabuna-Salvador-Curaçá Belt ( Figure 1). Late NW magnetic trends are related to hydrothermal alteration ( Figures 5 and 6a). Late mafic and felsic dikes show the same trend.
The north-to northeast-trending magnetic anomalies of the Vale do Curaçá copper district correspond to anastomosed structures that have been geologically mapped throughout the district. The NS-trending structures have been overprinted by at least two cycles of shearing [34,46]. These later shear faults are recorded on the TAHG map (Figure 6a).
The copper occurrences and mines in this district, including areas of strong hydrothermal alteration, are associated with magnetic anomalies (TMIs) of~25.98 to~100.59 nT ( Figure 5). Magnetic susceptibility ranges from 100 to 800 × 10 −3 (Figure 7b). The magnetic susceptibility of ore zones of the Vale do Curaça copper district is ten times higher than that of the Riacho do Pontal district (Figure 7a,b).

Riacho Do Pontal Copper District
The Riacho do Pontal district sits on a gravity gradient that separates a gravity high to the north from a NW-trending gravity low to the south (Figure 6b). The Bouguer anomaly low is likely related to low-density gneiss that partly outcrops in the Sobradinho Remanso and the southern portion of the Riacho Seco complex (Figure 1). Ria4 copper occurrences within the Riacho do Pontal district are situated along the edge of the gravity gradient. Densities measured from rocks within the Ria4 mine show that the less altered rocks have densities of 2.7 to 2.78 g/cm 3 . Rocks that have been more heavily altered have densities that overlap with those of less altered rocks and range from 2.75 to 3.1 g/cm 3 (Figure 7b).
The northwest trending Bouguer gravity low underlying the Riacho do Pontal district geologically corresponds to the collisional suture zone located between the higher density cratonic block and the lower density crust forming the Riacho do Pontal mobile belt [50,68]. Ophiolite bodies have been mapped in this region [52].
The northeast portion of the Riacho do Pontal copper district is characterized by a more positive gravimetric anomaly with values >−35 mgal. The gravity high corresponds to foliated orthogneissic rocks of the Riacho Seco complex. The Riacho Seco copper occurrences are situated along the gravity gradient, which likely marks a structure that has juxtaposed the denser crust of the craton against the lower density crust underlying the Riacho do Pontal mobile belt. The area of the gravity gradient includes biotite-garnet rich rocks, which occur along shear zones and host a small number of copper occurrences.
In the Riacho do Pontal Copper District, the copper occurrences (Ria4 and Riacho Seco) are not related to high magnetic anomalies. The Ria4 and Riacho Seco occurrences are hosted, respectively, in nonmagnetic gneiss of the Sobradinho Remanso and Riacho Seco complexes. The gneissic rocks are pervasively hydrothermally altered and have a range of magnetic susceptibility of 0.2 × 10 −3 SI in less altered rocks to 12.2 × 10 −3 SI in more hydrothermally altered rocks. On a regional scale, the NW-SE trending magnetic anomalies overprint longer wavelength NS trending anomalies related to older and deeper parts of the São Francisco Craton [69].
Copper deposits located in the belt of the Riacho do Pontal copper district are mainly situated in zones with or without weak magnetic anomalies (TMI) with intensities below 40 nT ( Figure 5).

Vale do Curaçá Copper District
The 3D magnetic susceptibility model of the Vale do Curaçá copper district reveals magnetic features corresponding to numerous geologic structures that reflect episodic tectonic activity. The 3D interpretation was built using field geological and structural data. More than 150 geological points were visited, more than 50 thin sections were described, and at least 20 drill holes were observed. Structural data were systematically collected. A series of parallel magnetic trends aligned N-S and dipping 75-85 • to the east are overprinted by at least two episodes of NS-trending shear faults (D1 and D2; [20,34,47]). The Caraíba mine is positioned along a NS-trending subvertical shear zone, showing an anastomosing shape in both the horizontal and vertical directions. A series of subparallel N-S shear structures have been active and responsible for multiple reactivation events. N-S shearing magnetic zones are related to thrust shear zones (D1) with convergence to the west and are overprinted by strike slip faults. The D1 and D2 events controlled the distribution of hydrothermal alteration and mineralization. The NE-SW late shear zones are subvertical [20,34,47].
The Caraíba mine bodies are strongly related to magnetic and gravimetric NS trends which are evident when analyzing the inversion results. In Figure 9, the inverted model is represented through the delineation of the body with high magnetic susceptibility, and high-density contrast. We have chosen a discriminant level that allows to represent the increased density and magnetization through positive density contrast values (>0 g/cm 3 ) and susceptibility values above 0.005 SI. In a 3D rendering the model of these relatively high values is illustrated in Figure 9, where the known shear faults and the position of the Caraiba open pit mine has been included. A characteristic large-scale densification and increased susceptibility is found below the mine, as is expected considering the results concerning the rock properties we have described above, regarding the rocks that encase the IOCG mineralization ( Figure 9). We find that the Caraíba mine and Paredão occurrences show a strong relationship to magnetic anomalies. The other occurrences along the Vale do Curaçá district are not directly related to high gravimetric and magnetic anomalies. The host rocks of the Caraíba mine are granulitic orthogneiss and ultramafic bodies. The mine is situated approximately 1 km from the N-S shear zone that was displaced by NW-SE shear zones. Sodic and potassic hydrothermal alterations are mapped along shear zones [20,29]. To demonstrate the inversion results in more detail, we show a section of the model focused on the mine. In a 3D tectonic setting, the mine is intrinsically related to high magnetic susceptibility associated with a positive density anomaly ( Figure 10). The orebodies are oriented in a subvertical direction and are controlled by stretching lineation [34]. The model resulting from the inversion of the potential field data is consistent and correlates with the outcropping geology and structural features mapped in the Caraíba mine. In Figure 10, we show a cross section through the model, centered on the mine, and extending 500 m on each side of the mine. The depth extent of the model reaches 8000 m. The blue model shows the relatively high susceptibility, with values being greater than 0.009 SI, whereas the mixed blue and red graph shows the superposition of susceptibility and positive density contrast. The black circle shows the location of the open pit of the Caraíba mine. Shown in grey is the topography, which is completely flat.
The magnetic trends (magnetic susceptibilities > 0.009 SI) are juxtaposed with shear structures and faults with different magnetism intensities ( Figure 10). The Caraíba mine correlates with magnetic signatures higher than those of occurrences and deposits situated near the mine.
In Figure 11, the results from the gravimetric and magnetic inversions along an EW profile are shown, with the Caraiba mine in its eastern portion. Highly magnetic (>0.005 SI) and dense (>3.20 g/cm 3 ) areas are related to rock properties at superficial depths. In the same profiles, at 3000 m depth, magnetization anomalies show continuity up to deeper zones (−8000 m), in the case of density anomalies, in deeper zones (−3000 m), the anomalies do not show continuity: this is interpreted geologically as graffitiated rocks that emerge at shallower depths and probably metasedimentary rocks at greater depths. We find that hydrothermal zones are related to magnetic zones in shear zones, opening the possibility for the continuity of mineralization in deeper zones. model shows the relatively high susceptibility, with values being greater than 0.009 SI, whereas the mixed blue and red graph shows the superposition of susceptibility and positive density contrast. The black circle shows the location of the open pit of the Caraíba mine. Shown in grey is the topography, which is completely flat. Figure 10. 3D representation of the results from magnetic (susceptibility in terms of >0.009 SI) and gravimetric (density > 0 g/cm3) inversion integrated for the Caraíba mine. The Caraíba mine is strongly related to susceptibility and density anomalies. It shows potential for continuity in depth in ore zones related to susceptibility and density anomalies along shear traps.
The magnetic trends (magnetic susceptibilities > 0.009 SI) are juxtaposed with shear structures and faults with different magnetism intensities ( Figure 10). The Caraíba mine correlates with magnetic signatures higher than those of occurrences and deposits situated near the mine.
In Figure 11, the results from the gravimetric and magnetic inversions along an EW profile are shown, with the Caraiba mine in its eastern portion. Highly magnetic (>0.005 SI) and dense (>3.20 g/cm³) areas are related to rock properties at superficial depths. In the same profiles, at 3000 m depth, magnetization anomalies show continuity up to deeper zones (−8000 m), in the case of density anomalies, in deeper zones (−3000 m), the anomalies do not show continuity: this is interpreted geologically as graffitiated rocks that emerge at shallower depths and probably metasedimentary Figure 10. 3D representation of the results from magnetic (susceptibility in terms of >0.009 SI) and gravimetric (density > 0 g/cm3) inversion integrated for the Caraíba mine. The Caraíba mine is strongly related to susceptibility and density anomalies. It shows potential for continuity in depth in ore zones related to susceptibility and density anomalies along shear traps.  . 3D representation of the results from magnetic (susceptibility in terms of >0.009 SI) and gravimetric (density > 0 g/cm 3 ) inversion integrated for the Caraíba mine. The Caraíba mine is strongly related to susceptibility and density anomalies. It shows potential for continuity in depth in ore zones related to susceptibility and density anomalies along shear traps. Geologically, the regional magnetic trends are related to more intense metamorphic granulitic processes overprinting geological units: the Bom Despacho gneiss and the Surubim, Banguê, and tonalitic to granodioritic orthogneiss. Late NE-SW shear zones are quite common, crosscutting N-S early magnetic trends [20,29,34]. Ore zones from the Caraíba and Surubim mines have a range of magnetic susceptibilities of 400 × 10 −3 SI to 750 × 10 −3 SI.
In Figure 12 we show the integrated model, with the gravimetric Bouguer values, overlaid on the inverted magnetic and density model. In the 3D view, the red domains are the increased magnetic susceptibility areas, and in grey the positive density contrast areas. Moreover, the faults are shown. The Caraíba mine is associated with a shear structure, and a high density and magnetic susceptibility zone. It can be seen that the copper occurrences and mines are prevalently associated to increased susceptibility and increased density values. In fact, the occurrences align along the trend I in which we find structures and high density and susceptibility. Moreover, it is found that in the area Central North-East, these conditions are similar, but up to now no occurrences have been explored. In the future, this area should be studied in more detail.
5 Figure 12. The IOCG mineral system in the region of the Vale do Curaçá Copper District, in regional terms, shows two main prospective trends. The trends are positioned in the contact zone between domains with high magnetic susceptibility (>0.005 SI) and density (>0 g/cm 3 ) and show contrasts in density and magnetite content.
Three-dimensional magnetic and gravity inversions show two main NNW prospective trends, both with a sigmoidal shear shape. These trends are positioned in the zone of contact between domains with high magnetic (susceptibility > 0.005 SI) and density anomalies (density > 0 g/cm 3 ) and show contrasts in dense rocks and magnetite content. The trends envelope granulitic rocks but are also the most hydrothermal rocks in the Vale do Curaça district. Trend I is 40 km in length and 10 km in width, hosts the Caraíba, Surubim, and Vermelho mines and other deposits and is the most important trend for host IOCG deposits in the Vale do Curaçá Copper District. Trend II is 15 km in length and 5 km in width. Geologically, these trends are hosted mainly along points of contact between Surubim gneiss, tonalitic to granodioritic orthogneiss, and gabbronorite and/or norite bodies ( Figure 12).
The alteration at the Vale do Curaça copper district has a hydrothermal signature characterized by intense potassic-iron alteration indicative of a high temperature system [1,8,10].
Alteration and associated mineralization at the Vale do Curaça Copper District are associated with magnetic and gravity anomalies as observed in other IOCG deposits from Carajás Mineral Province, Brazil; the Gawler Craton, Australia; Kiruna, Sweden; and the Great Bear Magmatic Zone, Canada [1,5,8,10].

Riacho do Pontal District
The 3D magnetic inversion over the Vale do Curaçá copper district shows that the Ria4 prospect and Riacho Seco deposit do not show a correlation with magnetic susceptibilies. Occurrences of copper in the region are positioned at magnetic lows. This observation is explained by the fact that the mineralized zone is related to hematite zones and iron oxide content is low (<2%). The exploratory drilling program (Ria4-DH01) crosscut 32 m @ 1% copper [24,28].
In the Riacho do Pontal district, the copper occurrences are not associated to evident positive density and magnetic susceptibility values. This is coherent with the fact that the rock samples showed gneiss characteristics, with no evident density and susceptibility anomalies. The copper occurrences are associated to the shear faults, but not to increased density and magnetite content. Probably the mineralization and magnetization in the shear zone is much weaker compared to the Val do Curacao district.
The results obtained from magnetic inversions conducted in the region of the Riacho do Pontal copper district show weak magnetic anomalies along a NW-SE trend. The discrete linear magnetic features are related to shear zones. These NW-SE features overlap older and deeper NE-SW features. However, an analysis of Riacho do Pontal district magnetic anomalies in 3D in the Riacho Seco deposit shows a low magnetic anomaly close to the surface and a change in depth (>400 m) to an area with a high magnetic signature (>0.005 SI) ( Figure 13). These results allow for the possible development of a continuous model for follow-up research on these deposits as related to deeper magnetic susceptibility anomalies.
More recently, a new generation of ore system exploration models has been developed to take advantage of the massive volume of data available to the exploration industry [76]. Improving exploration competitiveness will be necessary with the more sophisticated integration of different datasets. The challenge remains to decrease the size of exploration areas and increase the success rate of drilling programs [76,77]. Thus, for the Vale do Curaçá and Riacho do Pontal copper districts, the main footprint that shows the ore system signature was identified ( Figure 13). The more important geological, structural, geochemical, mineralogical, and petrophysical parameters of these copper districts are summarized below and illustrated in Figure 13. At a regional scale, large major lithospheric structures control the traps of IOCG mineralizing systems or heat magmas, or fluid flow is frequently evident in gravity and magnetic databases [4]. The outline of the inversion model for Riacho do Pontal obtained using these methods fits well with geological and structural features mapped in the field.

Lithospheric Setting Expressed by the Seismic Tomography Potential
The high zones (in red) on the tomographic map represent high-velocity domains that denote microcontinental blocks (mantle lithosphere domain). They are flanked by low velocity domains interpreted as the refertilization of the subcontinental lithospheric mantle (SCLM) and/or higher temperatures along fault zones. Frequently, deposits are concentrated along prominent translithospheric structures, particularly in lower-velocity regions (blue) or on the flanks of velocity highs (courtesy of [63,64]).
A world map of the lithosphere produced through the integration of geophysical, geological, and geochemical data shows that approximately 70% of the existing SCLM may be of Archaean origin [63,64]. The implication for an Archaean SCLM is that most preserved Proterozoic crust overlies the older Archaean SCLM variably refertilized and metasomatized by mantle melts associated with convergent margin and mantle plume processes [63,78].
In general, seismic tomography images can show boundaries of cratons and discontinuities related to the SCLM. These discontinuities are favorable for the repeated ascent of magmas [76][77][78]. Preserved cratonic regions (SCLMs) can be mapped as areas of high seismic velocity. At a global scale, it has been shown that smaller volumes of magmas directly originating from the mantle are positioned at the craton edge (Figure 4).
Domains of high seismic velocity can represent ancient blocks of metasomatized SLCM within a discrete mantle lithosphere domain flanked by major faults [63,64]. The IOCG Carajás province and Vale do Curaçá and Riacho do Pontal copper districts are located in the transition zone between domains of high and low velocities ( Figure 4) and may reflect the refertilization of the SCLM and/or zones of higher lithospheric and mantle temperatures. Ore deposits that originate from the mantle or from transition mantle/crust magmas such as diamond, platinum-group elements (PGEs), Ni-Cu-(PGE), copper porphyry and IOCG deposits are situated along cratonic edges and can be related to the tectonic evolution of the SCLM [63,64,78]. In this context, the extensive gradients flanking either side of the high velocity zone through east-central Brazil should be evaluated in more detail (Figure 4). Craton edges are characterized by a thinning of lithosphere and are strongly tectonized zones creating transtensional or transpressional sites along translithospheric faults [62]. These major shear or fault zones are frequently used for melt introduction into the crust. For example, IOCG districts in Australia, including the world class Olympic Dam in the Gawler craton and within the Eastern Succession of Mount Isa Inlier in Australia, are related to the presence of large igneous provinces (LIPs) that may contain a considerable amount of intrusive felsic rocks [8,10,71,[79][80][81][82][83].
Fertile provinces that lead to the formation of mantle-related world-class deposits require melting of the SCLM on the margins of cratons, which is responsible for magmatic ascent associated high heat flow and volatile flow [25,65]. The deposits occur in long-lived subparallel regional structures that have been the focus of multiple reactivation events [24,26,28,30,34,47]. In the Riacho do Pontal and Vale do Curaçá Copper Districts, structural interpretations show that early thrust zones are reused several times and are overprinted by strike-slip shear zones [47].

Structural Control
Sigmoidal or "SC-shaped" structures are indicative of major shear systems [18,23,25,74,75]. Fractal patterns that mimic "SC" structures are systematically repeated from macro-to microscales. These structures serve as critical pathways for the upward flux of fluids promoting interactions between fluid and rock. Duplex, shear band, and asymmetrical structures are quite common along the shear system. The structural controls and traps used for orogenic gold in shear zones are the same as those used for IOCG deposits [24]. The more important traps for host orebodies include (a) interconnection zones between "SC" structures; (b) fractal second-or third-order "C" structures; (c) bends along shear zones in both horizontal and vertical planes where high-grade breccia are typically hosted; (d) fault intersections highly favorable for host breccia ores; and (e) ductile shear zones reactivated by late brittle faults that enrich older zones and that can generate new mineralized bodies. The permeability and degree of fluid overpressure (where fluid pressure > lithostatic pressure) greatly influence the formation of brecciated zones along shear zones.
Frequently, lithological boundaries between units with different competencies are crucial for hosting high deformation corridors. In this way, zones of contact with competent rocks such as paragneiss, mafic granulite, granulitic orthogneiss, and banded iron formations (BIFs) are usually trapped for the percolation of hydrothermal fluids related to shear zones. In the Vale do Curaçá copper district, considering only the epigenetic IOCG mineralization process, the host rocks are predominantly mafic and felsic granulitic rocks.
The orebodies often assume a sigmoidal shape, as they are often hosted in shear zones [8,11,70]. Often, early, sin and late-tectonic breccia pipe zones are found along hydrothermal IOCG systems. The evolution of IOCG-rich systems entails development by the direct melting of the SLCM, positioning in shear zones initially as pipe breccia, and continuous reworking by shear deformation. Variation in the direction and dipping along shear zones, both horizontally and vertically, acts as a structural trap for both hydrothermal and pipe breccias. Typically, breccia-rich ores are situated in a central hydrothermal system zone in transtensional structures.
Hydrothermal alteration is related to intermediate-to high-temperature fluids in the Riacho do Pontal and Vale do Curaçá districts. The alteration halo extends from 1 to 5 km in width, following the major shear zones. Sodic alteration is more distal and potassic, and calcic and ferric alterations are more proximal in both districts. The southern deposits and Sussuarana, Caraíba, and Surubim mines are related to deeper and high-temperature sodic/potassic hydrothermal alteration in contrast with the Vermelho deposit, in which hydrolytic alteration occurred in shallower zones in the IOCG hydrothermal system as initially described by [5,10]. Typically, a series of magnetic dipoles are linked to shear zones and hydrothermal halos.
The general hydrothermal zoning pattern for IOCG deposits is usually vertical from magnetite-dominant at depth to hematite-dominant at shallower levels [5,10,[84][85][86][87][88][89][90]. Using this model, it is possible to interpret the Riacho do Pontal District as positioned in the shallowest portion of the IOCG system ( Figure 14). Thus, additional work must be performed to identify deeper (or more eroded) areas of the IOCG system related to more magnetic and density targets, such as those identified in the Riacho Seco project.
Hydrothermal alteration is related to intermediate-to high-temperature fluids in the Riacho do Pontal and Vale do Curaçá districts. The alteration halo extends from 1 to 5 km in width, following the major shear zones. Sodic alteration is more distal and potassic, and calcic and ferric alterations are more proximal in both districts. The southern deposits and Sussuarana, Caraíba, and Surubim mines are related to deeper and high-temperature sodic/potassic hydrothermal alteration in contrast with the Vermelho deposit, in which hydrolytic alteration occurred in shallower zones in the IOCG hydrothermal system as initially described by [5,10]. Typically, a series of magnetic dipoles are linked to shear zones and hydrothermal halos.
The general hydrothermal zoning pattern for IOCG deposits is usually vertical from magnetitedominant at depth to hematite-dominant at shallower levels [5,10,[84][85][86][87][88][89][90]. Using this model, it is possible to interpret the Riacho do Pontal District as positioned in the shallowest portion of the IOCG system ( Figure 14). Thus, additional work must be performed to identify deeper (or more eroded) areas of the IOCG system related to more magnetic and density targets, such as those identified in the Riacho Seco project.

Exploration Potential
Through 2D and 3D mapping, it is possible to verify variation in these parameters with distance in the direction of mines, deposits, or occurrences. An increase or decrease in the intensity of these processes represents vectors towards the IOCG systems studied (Figure 14).
The main footprints of the two IOCG districts are described below:

Exploration Potential
Through 2D and 3D mapping, it is possible to verify variation in these parameters with distance in the direction of mines, deposits, or occurrences. An increase or decrease in the intensity of these processes represents vectors towards the IOCG systems studied (Figure 14).
The main footprints of the two IOCG districts are described below: Riacho do Pontal District: (a) Structural: increase in shearing intensity and in the density of structures in the direction of secondary structures. Foliation tends to be vertical. Regardless of the structural regime, most deposits in both districts are controlled by bends in the host rocks. (b) Stream sediment and soil: anomalies of Cu, Ni, Fe, and La in stream sediment and Cu (Ni, Ce and La) in soil. (c) Low magnetic anomalies related to Kd (Kd-anomalous potassium) and shear zones [24]. (d) Low magnetic anomalies (abundant hematite related to ore zones).
Vale do Curaçá District: (a) The district of Vale do Curaçá is associated regionally with high gravimetric and magnetic anomalies corresponding with the portion west of the ISCO. (b) Our studied area is inboard of a relatively high velocity region, bordered eastwards by a low velocity region. The high-velocity region presumably indicates the cratonic area, with the lower velocity region demarcating the craton limit. (c) The Caraíba mine bodies are strongly related to magnetic and gravimetric NS trends interpreted to have susceptibility with >0.09 SI and density >0 g/cm 3 .
Hydrothermal alteration: an increase in the anomalies is related to hydrothermal processes: (i) distal sodic alteration and (ii) the mineralized zone intensely related to calcic-ferric alteration.
However, an atypical paradigm remains in the Vale do Curaçá district: the selection of new targets not necessarily hosted in ultramafic rocks. The IOCG system can develop in any host rock at varying depths. A full understanding of the IOCG model for mining companies can greatly expand exploration in the districts via data integration and the processing of structural, geological, and geophysical databases [11,65].

Conclusions
Knowledge of the signature of IOCG deposits in Brazil shows that these deposits can be hosted in any type of rock with specific styles of hydrothermal alteration and at different depths. Three-dimensional geological modeling, when combined with geophysical and structural interpretation, can facilitate the spatial and temporal understanding of the hydrothermal systems of the Vale do Curaça and Riacho do Pontal districts.
The main conclusions obtained from the present study of the Vale do Curaçá and Riacho do Pontal copper districts are as follows: 1.
The Vale do Curaçá and Riacho do Pontal copper districts are hosted in distinct geological and tectonic settings, and an older IOCG mineralization event occurring within the São Francisco Craton is related to Palaeoproterozoic (ca 2 to 2.2 Ga) hydrothermal processes. A younger phase of IOCG mineralization in the Neoproterozoic (ca 750 to 570 Ma) has been recognized in the Brazilian mobile belt.

2.
Both IOCG districts developed in an arc tectonic setting [50] overprinted by strike slip shear zones.

3.
The interpretation of gravimetric and magnetic data shows different magnetic and gravimetric signatures for these two copper districts. 4.
The Vale do Curaçá copper district is in a high regional gravimetric domain. It is 110 km in length and 22 km in width and shows a strong relationship with magnetic and gravimetric anomalies. The Caraíba, Surubim, and Vermelhos mines and some copper occurrences are positioned at this site. 5.
The regional gravimetric signature of the Riacho do Pontal District is mainly related to a negative Bouguer NW-SE anomaly. This negative Bouguer anomaly represents the collisional suture between the cratonic block and mobile belt Riacho do Pontal. The IOCG occurrences of the Ria4 prospects are found in these low gravimetric signature domains. The rocks, slightly altered, have a density of 2.7-2.78 g/cm 3 , and the most heavily altered rocks show densities of 2.75-3.1 g/cm 3 . 6.
The magnetic signatures are quite distinct. Copper corridors in the Vale do Curaçá copper district are associated with magnetic anomalies located in the western portion of the Itabuna-Salvador-Curaçá orogen. On the other hand, the copper deposits located in the belt of the Riacho do Pontal district are situated in zones with or without very weak magnetic anomalies. 7.
MVI shows two main NNW prospective trends in the Vale do Curaça district. Trends I and II have a sigmoidal shear shape. These trends are positioned in the contact zone between domains with high magnetization (susceptibility > 0.005 SI) and density anomalies (>0 g/cm 3 ) and show a contrast in density and magnetite content. Trend I is 40 km in length and 10 km in width. It hosts the Caraíba, Surubim, and Vermelho mines and other deposits and is the most important trend for hosting IOCG deposits in the district. 8.
The general hydrothermal zoning pattern for IOCG deposits is largely vertical from magnetite-dominant at depth to hematite-dominant in the upper levels. It is possible to interpret the Neoproterozoic Riacho do Pontal district as positioned in the shallowest portion of the IOCG system. Thus, additional work must be performed by searching for deeper (or more eroded) areas of the IOCG system related to more magnetic and dense areas (targets in the northern portion of the Riacho Seco target). 9.
Large areas flanking the high-velocity seismic tomography image that already control, on a broad scale, the Carajás IOCG province and Riacho do Pontal and Vale do Curaçá copper districts in Brazil should be evaluated in more detail in a further exploratory program.