Geochemistry of the Northern Paraná Continental Flood Basalt ( PCFB ) Province : implications for regional chemostratigraphy

2Departamento de Física da Terra e do Meio Ambiente, Instituto de Física, Universidade Federal da Bahia – Salvador (BA), Brazil. E-mail: eduardo.junior@ufba.br 3Instituto de Astronomia, Geofísica e Ciências Atmosféricas, Universidade de São Paulo – São Paulo (SP), Brazil. E-mail: leila.marques@iag.usp.br 4Instituto de Geociências e Ciências Exatas, Universidade Estadual Paulista “Júlio de Mesquita Filho” – Rio Claro (SP), Brazil. E-mail: nardy@rc.unesp.br 5Instituto Geociências, Universidade de São Paulo – São Paulo (SP), Brazil. E-mail: natasha.smarteleto@gmail.com


INTRODUCTION
Continental Flood Basalts (CFB) are a type of Large Igneous Province (LIP) that occur as extensive outpourings of juvenile material to the continental crust and represent major, short-lived igneous events on Earth.Despite the wealth of studies on CFB, as well as on the associated geodynamic processes, their genesis remain a controversial subject (Coffin & Eldholm 1994, 2005, Bryan & Ernst 2008, Rocha-Júnior et al. 2012, 2013, Ernst 2014, Marques et al. 2017).Usually, CFB provinces are associated with sedimentary basins, having an extremely complex internal architecture, preserving essential information on the onset of flood volcanism, as well as on paleo-environments (Jerram & Widdowson 2005).
Several chemostratigraphy studies have been carried out on volcanic piles of CFB provinces, such as the Siberian Traps, Karoo, Deccan, and Columbia River, in order to understand the architecture of magmatic dynamics (Ernst 2014).
On the other hand, few chemostratigraphy investigations of lava piles have been accomplished on the Paraná Continental Flood Basalt (PCFB) province.It is noteworthy that, until now, no chemostratigraphy research has been conducted in the northern part of the magmatic province, with all studies concentrated on the central and southern parts of the PCFB province (e.g., Nardy et al. 2002Nardy et al. , 2008)).This study describes for the first time samples from eight boreholes, as well as some surface samples arranged in an NEE-SWW profile located in the northern PCFB, to better understand the complex internal architecture of lava flows Brazilian Journal of Geology, 48(2): 177-199, June 2018 and reconstruct the eruptive processes associated with the different magma types.

THE PARANÁ CONTINENTAL FLOOD BASALTS AS A LIP
The attention to LIPs may be justified by their great relevance on Earth for testing magma generation models, providing temporal marks for stratigraphic correlations of lavas (mainly by geochemistry), paleocontinental reconstructions, and for the importance they have in economic geology (metals, diamonds, rare earth elements -REE, hydrocarbons and aquifers).
Another important aspect is that many LIPs that now lie partially on oceanic and continental lithosphere erupted where rifting led to continental breakup (e.g., Karoo, Ferrar, Central Atlantic Magmatic Province, and Deccan Traps), which may indicate that the locations of melting anomalies are controlled by stress (Anderson 2002), and consequently controlled by plate tectonics.Therefore, the PCFB volcanism agrees with this LIP criterion, since this magmatism took place approximately 10 Ma before the South Atlantic Ocean opening.
The term "Large Igneous Province" was initially proposed by Coffin & Eldholm (1991, 1992, 1993a, 1993b, 1994) to identify a variety of provinces composed by mafic igneous rocks that could reach > 0.1 Mkm 2 .These provinces are found in oceanic and continental areas and include both volcanic and associated intrusive rocks, isolated dyke swarms, sills and other forms of intrusions.Although many LIPs are restricted to the Cenozoic and Mesozoic (Coffin & Eldholm 2005), such as a large number of continental basaltic provinces (e.g., Paraná, Karoo, Afro-Arabia), rift volcanism (e.g., India-Western Australia, Bushveld), oceanic plateaux (e.g., Ontong Java, Shatsky Rise) and oceanic basaltic provinces (e.g., Nauru Basin, East Mariana), there are others with ages varying from the Archaean to Paleozoic (Ernst & Buchan 2001, 2003).
Studies (e.g., Sheth 2007, Bryan & Ernst 2008, Ernst 2014) have collected different types information on these provinces, such as the lithology diversity (presence of rare acid volcanic at the ending of the volcanism), geodynamic characteristics, origin (related or not with mantle plumes) and duration of the main phase of magmatism (~1-5 Ma), among others.
It is noteworthy that most of the LIPs are composed of SiO 2 < 56 wt.% rocks, represented by phenocrysts-poor tholeiitic basalts and basaltic andesites.However, a significant geochemical variation in the stratigraphy of flows in the most voluminous LIPs (CFB, Oceanic Flood Basalt -OFB and Silicic LIP -SLIP) is notable (Hooper 1997).Bryan and Ernst (2008) considered most of them, if not all, characterized by products with a bimodal magmatic activity with SiO 2 mostly clustered between 45-55 wt.% and 65-75 wt.% (e.g., .Acidic terms are always present at the end of the volcanic event, being usually less expressive (e.g., Paraná with 2.5 % of acidic rocks; Kirstein et al. 2000;Lustrino et al. 2005, 2010, Nardy et al. 2008).
Stratigraphy and correlation of LIP lavas, more precisely of CFB type, have demonstrated strong similarities (MacDougall 1988, Coffin & Eldholm 1994, Jerran & Widdowson 2005, Bryan & Ernst, 2008, Machado et al. 2014, 2015).Basic lava flows with a total thickness of several kilometers have shown intercalation with sediments in almost all volcanic sequences, and the geochemical stratigraphy is marked mainly by variation of TiO 2 and P 2 O 5 .
Geochemical stratigraphy has been used in classical models (e.g., Deccan, Paraná, Columbia River) for correlation of flows, mapping and temporal relation, as well as for recognizing formations, members and subgroups (Subbarao et al. 1988, Beane et al. 1986, Mahoney et al. 2000, Nardy et al. 2008, Riedel et al. 2013, and references therein).
As the sixth largest Meso-Cenozoic LIP on Earth, the PCFB (Fig. 1) and its extension in the African Plate (known as Etendeka) are preceded by Madagascar ( 1 (Bellieni et al. 1986, Nardy et al. 2008) and their continuation in the African Plate (Etendeka) during Rifting Subphase II of Gondwana (Almeida et al. 2012).Bertrand et al. 2014).The area occupies 75% of the entire Paraná Basin and the estimated volume is 0.8 Mkm 3 , but, according to Courtillot & Renne (2003), the original volume could have been > 2.3 Mkm 3 .In fact, post-Cretaceous denudation (Quintas et al. 1999, Franco-Magalhães et al. 2013) was deeper in the northern part than in the southern part of the PCFB.The uplift of marginal areas led to the occurrence of a sedimentary unit (Bauru Sub-Basin;Basilici et al. 2012) superposed to the volcanism.This situation considerably reduces the outcrops of basaltic rocks in the northern portion when compared to the southern one.

Brazilian
On the other hand, in the case of Deccan, there are 10 volcanic formations grouped into 3 subgroups, occupying 1.8 Mkm 2 and accounting for more than 1 Mkm 3 of lava (Beane et al. 1986, Mahoney et al. 2000).These formations are divided into more than 20 members that are separated by mineralogy, crystal size, and mainly by geochemistry (concentrations of SiO 2 , K 2 O, MgO, TiO 2 , P 2 O 5 , Ba, Sr, Zr and Sr, as well as 87 Sr/ 86 Sr isotope compositions).It should be noted that, according to Subbarao et al. (1988), some geological formations can exceed > 300 km of extension.
Similar to Deccan, Columbia River is divided in 5 formations and identified by geochemistry (Faure 2001, Camp & Hanan, 2008, Riedel et al. 2013) but also by magnetic polarity patterns.They spread over an area of ~0.24 Mkm 2 with an estimated volume of about 0.3 Mkm 3 , while flows are tabular and extend over for > 600 km.
They revealed the existence of two chemically distinct rock groups (and sub-provinces) based on the concentration of TiO 2 that were called high titanium (HTi; TiO 2 ≥ 2 wt.%) and low titanium (LTi, TiO 2 < 2 wt.%) rocks.
Compared to uncontaminated LTi rocks, the HTi ones are enriched in Ba, K 2 O, U, Sr, La, Ce, Ta, P 2 O 5 , Hf and Zr.In general, those of the HTi group are located in the northern PCFB, above the Piquiri River Lineament, while the LTi appears in the southern part, below the Uruguay River Lineament (Fig. 1).
In the southern part of the province, where the LTi rocks strongly prevail, the occurrence of HTi was identified, corresponding to only 7% of the basic igneous activity (Piccirillo et al. 1988, Marques et al. 1989, Peate et al. 1992).According to these references, scarce LTi basalts were emplaced in the HTi sub-province (especially in the northwest region), representing 6% of the basic rocks of the Paraná Basin.However, there are significant geochemical differences between the LTi basalts occurring in the northern (Ribeira type) and southern (Esmeralda and Gramado types) sub-provinces, as well as between the HTi rocks in these two regions, with Pitanga and Paranapanema occurring in northern PCFB and Urubici in the southern.After a detailed study of the rocks that outcrop at the extreme west of the province, Machado et al. (2009) proposed a new compositional variation for the Ribeira magma type.
In the African counterpart (Fig. 1; west of the African continent), basaltic rocks are equivalent to Gramado (Namibia) and Urubici (Angola), denoted as Tafelberg and Khumib, respectively (e.g., Marsh et al. 2001, Ewart et al. 2004).The differences between such magma types were mainly based on contents of TiO 2 and incompatible trace elements (Sr, Y and Zr), whose main geochemical characteristics are reported in Tab. 1.
As already mentioned, the distribution of these magma types in the PCFB is not random, but there seems to be a pattern of occurrence.An attempt to define a stratigraphy Table 1.Distinctive geochemical compositions of basic rocks of the so-called Serra Geral Formation within the Paraná sedimentary basin according to Piccirillo et al. (1988), Marques et al. (1989) and Peate et al. (1992).Ribeira magma type data are from Machado et al. (2009).for these magma types was accomplished by Turner et al. (1994) and Stewart et al. (1996).Through drilling, associated with geochemical data, it was verified that in the southern PCFB the Esmeralda basalts cover the Urubici and Gramado ones, whereas in the northern part the Paranapanema flows overlap the Pitanga lithotypes.However, the distribution and stratigraphy of northern HTi magma types were contested by Machado et al. (2007Machado et al. ( , 2014Machado et al. ( , 2015)), who reported the presence of some Urubici basalts in the extreme northeast of the province.These basic rocks, until then only found in southern PCFB, are present in flows and nearby sills in that region.Some Urubici dykes also intrude Archean basement rocks (Southern Espinhaço) located at the border of the São Francisco Craton (Marques et al. 2016).

LOCATION AND GEOLOGICAL ASPECTS
The boreholes are located at the center-northern PCFB in the Paraná Basin (Fig. 2).The maximum reach depths exceed 1,100 meters, crossing the complete sequence of basaltic rocks.Four boreholes cut up to 200 meters of sedimentary rocks of the Bauru Sub-Basin, while the other four cross soil, whose layers can be as thick as 40 m in the investigated region.
The eight boreholes are aligned in a roughly ENE-WSW direction and cover a distance of 671 km.Surface samples were collected near to FR, B, and IVI boreholes in the western limit of the basin, extending the section to about 877 km.
The oldest geological unit in the boreholes is the Triassic Piramboia Formation (Assine & Soares 1995).It consists of whitish, yellowish and reddish, medium-to very fine-grained sandstones with sub-angular and sub-rounded grains.There are also intercalations of thin layers of argillites and siltstones.
The Botucatu Formation (Early Jurassic-Cretaceous) is discordant and overlaps the previous unit (Scherer 2000).It is composed almost exclusively of medium to fine granules in sandstones, with good grain selection and high sphericity.
The PCFB province, which is named Serra Geral Formation in the stratigraphy of the Paraná Basin, lies on the top of the Botucatu Formation in concordant contact.Borehole data show that the thickness of the magmatic pile ranges from 80 m (ALT borehole) to nearly 1,020 m (IVI borehole).
Sometimes, upper and lower parts of the flows are marked by a few centimeters thick vesicular layers.Other boreholes (e.g., FR and GUA) show whole vesiculated flows as thick as 70-80 m.
It is worth mentioning the presence of flow intercalations with sandstones of the Botucatu Formation as thick as 14 m, 13 m and 22 m in IVI (at -350 m below sea level), PAI (at 470 m above sea level) and ALT (at 871 m above sea level) boreholes, respectively.Also, isolated sills were recognized between lava flows, with thicknesses similar to those of the sheet flows, despite being massive.They were identified in the FR (at 160 m) and B (at 210 m) boreholes.Sills were identified by their distinctive concentrations of minor and trace elements, as well as by their cross-cut relationships with overlain and underlain flows.As such, a single, massive (14 m-thick) sill occurs in the PAI borehole (at level 236 m) intruding sandstones of Piramboia Formation.

MATERIAL AND METHODS
The boreholes were sampled at 2 m intervals.Each sample, constituted by a set of rock chips, is representative of the respective entire interval (Fig. 3).(Zalán et al. 1987).Symbols mark the location of the studied boreholes and surface samples (KS) and will be used hereafter.Geology, 48(2): 177-199, June 2018 The samples were placed in beakers, washed in an ultrasonic bath with water and dried in a small fume hood within sieves.After drying, the least altered, massive chips were selected, resulting in nearly 300 samples of basalts.There was also the concern to maintain regular sampling intervals to represent the largest possible number of flows, which were primarily identified by their both upper and lower vesicular limits.

Brazilian Journal of
Major and trace elements were obtained by X-Ray fluorescence spectrometry in the Geosciences and Exact Sciences Institute of Universidade Estadual Paulista.Calibration parameters, equipment readings and analytical errors (< 1 % for major elements, and < 5% for the trace elements) are available in Nardy et al. (1997) and Machado (2010).
Samples with a loss on ignition (LOI) > 2 wt.% or that presented some anomalous concentration for any of the elements were discarded.To make the dataset internally consistent, the major elements were normalized to 100% on a volatile-free basis.
A subset of 13 samples was selected for REE and other trace element analysis by Inductively Coupled Plasma Mass Spectrometry (ICP-MS), at the Geosciences Institute of Universidade de São Paulo.The analytical error is < 5% and the methodological procedure adopted can be checked at Navarro et al. (2008).Also, a set of 7 surface samples was also selected to complete the NNE-SWW stratigraphic profile shown in Fig. 2 as A-A '.
All data of the analyzed samples (major and trace elements) can be accessed in the Supplementary Table A1.

PETROGRAPHIC ASPECTS
Selected chips from the eight boreholes were immersed in epoxy-type resin for the assembly and preparation of thin sections.After the resin was dried, the samples were cut and properly wiped down.
The rocks are characterized by the presence of five essential mineral phases: calcic (augite) and rarely calcium-poor (pigeonite) clinopyroxenes, plagioclase, opaque minerals (probably magnetite and ilmenite) and olivine.Plagioclase is the predominant mineralogical component, reaching up to 50% of the rock volume.It occurs as microphenocrysts and in the rock matrix, with grain sizes < 0.2 mm.Microphenocrysts are subhedral, lath-shaped, with Carlsbad and Albite-Carlsbad twinning.They are sometimes zoned or show corroded edges by the matrix.Occasionally the plagioclase crystals are saussuritized.
Clinopyroxenes constitute up to 45% of the rock volume (maximum of 2% are pigeonite).They are mostly anhedral to subhedral, occurring in the matrix with a granular habit or as inclusions in plagioclase.
Olivine is entirely or partially altered to a reddish brown mineral, probably iddingsite.It occurs in the matrix surrounded by clinopyroxenes.The olivine crystals comprise less than 3% of the volume of the rocks.
Opaque minerals (probably magnetite and ilmenite) constitute up to 20% of the total volume.They appear interstitially in the matrix, and as microphenocrysts, often subhedral.
Dark-brown to black glass occurs interstitially in the matrix in variable amounts, but may constitute as much as 60% of the volume of the rocks.
The variations of major and trace elements relative to MgO are shown in Figs. 5 and 6.The investigated samples define trends for all elements except for Fe 2 O 3(t) .In general, SiO 2 , TiO 2 , Na 2 O, K 2 O, and P 2 O 5 increase, whereas CaO and Al 2 O 3 decrease with decreasing MgO.Although some scattering is observed, the trends are compatible with fractional crystallization, and form a continuous, evolutionary trend, similar to those presented in previous studies (Bellieni et al. 1984, Piccirillo & Melfi 1988, Peate et al. 1992, Marques et al. 1999, Rocha-Júnior et al. 2013).Major element variations, for the same MgO content, are broadly similar for lava flows and intrusive rocks, and they do not vary with the stratigraphic height (for the lava piles) or sampling location.
Five samples from PIR are enriched in MgO (6.51-7.47wt.%) in comparison to those from the other boreholes located on the NEE side of the profile.These samples occur within a 10 m-thick sequence, close to level 400 m.
The ALT borehole (at the NEE end of profile) presents a sequence of 10 samples located at its bottom, close to level 900 m, corresponding to a thickness of at least 25 meters.The P 2 O 5 content of these samples is higher than in the other samples (up to 0.81 wt.%), and this sequence also appears to include two distinct groups: one with P 2 O 5 > 0.58 wt.% and the other with P 2 O 5 < 0.52 wt.%.Samples from ALT borehole with 0.73 to 0.81 wt.% of P 2 O 5 are also enriched in La (41-47 ppm), Ba (575-657 ppm) and Rb (30-39 ppm).
In general, the concentrations of these incompatible elements show positive correlations with each other (not shown).It should be noted that the samples from ALT borehole with 0.73 to 0.81 wt.% of P 2 O 5 are also enriched in La (41-47 ppm), Ba (575-657 ppm) and Rb (30-39 ppm).
All 13 samples have similar REE contents, and their patterns display a fractionation of light REE (LREE) to heavy (HREE), as indicated by (La/Yb) N ratios ranging from 2.3 to 8.8 (Fig. 7).This behavior is expected from melting of garnet-bearing sources, as also pointed out in previous studies (e.g., Piccirillo & Melfi 1988).For most of the samples the enrichment in LREE relative to HREE is directly correlated with REE and SiO 2 concentrations.The Pitanga (pit) sequences have (La/ Yb) N ratios ranging from 7.3 to 8.8; (La/Sm) N from 2.3 to 2.6; and (Sm/Yb) N from 3.1 to 3.6.Paranapanema (par) flow sequences have less enrichment in LREE to HREE, as indicated by their (La/Yb) N ratios (2.1-2.4).Also, there are no large variations in Eu anomalies (interpolation of the values of normalized Sm and Gd) in pit (Eu/Eu* = 0.94 -1.01) and par (Eu/Eu* = 0.92 -0.95) sequences.
A normalized multielement diagram for the studied samples is shown in Fig. 8.The patterns of these rocks exhibit significant enrichment in LILE, such as Th, Ba, Rb, and K, relative to some HFSE (such as Nb, Zr, Ti, and P).Compared to the Paranapanema sequences, the Pitanga rocks have slightly enriched patterns, as also shown by Rocha-Júnior et al. (2013).
As expected, the Urubici samples have a more enriched pattern while those from the Ribeira samples have the most depleted patterns of the set of rocks analyzed.In addition, the incompatible trace element patterns of the analyzed samples reveal distinct negative anomalies for Nb, weak negative anomalies for Th, Sr (except for Urubici), and P, implying the participation of subcontinental lithospheric mantle source (SCLM), as observed for most PCFB (Piccirillo & Melfi 1988, Hawkesworth et al. 1992, Marques et al. 1999).As observed by Rocha-Júnior et al. (2013) and Marques et al. (2016), these patterns resemble those found in subduction environments (island-arc basalts), as well as in continental crust.

CHEMOSTRATIGRAPHY OF LAVA FLOWS
The FR borehole was selected as a reference in order to put constraints on the chemostratigraphy of lava flows based on its following characteristics: • high variation in P 2 O 5 , TiO 2 and Sr concentrations, as well as in Zr/Y ratios, allowing to recognize the occurrence of Pitanga and Paranapanema magma types (Tab.1); • large thickness (more than 800 m) of basalts; • large number of flows (= 7); • large amount of samples (103 samples) with low LOI (< 1.5 wt.%).Acronyms were attributed to the FR borehole samples to distinguish the flow successions based on their geochemical characteristics, following the temporal order and the consequent correlation with the other boreholes.
In this correlation, it was possible to assign ranges in concentrations of TiO 2 , MgO, P 2 O 5 , Sr, and Zr/Y ratios for the flows to all eight boreholes shown in Tab. 3, including the surface samples.It is important to note that the distinction among different flows was carried out by using P 2 O 5 , while the discrimination between magma types was mainly based on Zr/Y ratios, along with TiO 2 and Sr concentrations.

ALT borehole
The ALT borehole (Fig. 2) top is at an altitude of 945 m with a 5 m-thick soil cover.The contact between the basalt and the Botucatu Formation sandstone occurs at ~ 80 m.The total thickness of basaltic flows is therefore ~ 75 m.
The geochemical data (Fig. 9A) show that pit1 type flows occur up to level 915 m.Pit5 type appears down to the contact with the sandstone, where the final 20 m correspond to an amygdaloidal facies, not sampled for analysis.A 2 m-thick inter-layered pit1 type basalt with contemporary sandstone was found at level 871 m.

PAI borehole
The PAI borehole is located northeastwards and 120 km away from the ALT borehole.It is at 686 m altitude, and the basalts are discordantly overlain by a 50 m-thick sequence of  Table 3. Compositional variation for selected elements of different sequences of basalt flows identified in the boreholes of the studied area.Caption: par (Paranapanema magma type); pit (Pitanga magma type); rib (Ribeira magma type); uru (Urubici magma type).Bauru sandstones, where a discordant contact occurs with basalt at 636 m.There is intercalation with the sandstone of the Botucatu Formation at 496 m (25 m thick), overlain by a vesicular 30 m-thick flow at the top, and another contact with contemporary sandstone at 466 m (Fig. 10A).Also, there is a sill in the Piramboia Formation (at least 4 m thick) close to the level 236 m.Geochemical data (Fig. 10A) show that the basaltic lava flows of pit1 type occur down to level 625 m.Pit3 overlies them up to level 600 m, and then pit2 type flow appears at 550 m.The pit1 sequence occurs again from level 550 m to 530 m whereas pit 4 appears downwards being underlain by sandstone at 496 m.Pit2 lava flows take place between the Botucatu sandstone sequences (471-466 m).A pit1 sill intrudes the Piramboia Formation at the bottom of the borehole.
All lava flows show vesiculated levels at their upper parts, with easy identification of the limits from each one.

PIR borehole
The PIR borehole is located 80 km south of PAI, with its top at level 585 m (altitude).The uppermost 10 m are composed of soil down to contact with basalt.Alternating sequences of Pitanga and Paranapanema appear, and the contact with sandstone occurs at 105 m.There are no intercalations with sandstones, but there are sheet flow zones between 269 m and 200 m, with multiple vesicular layers that had to be avoided for the geochemical analysis.In the upper part, there are pit1-type basalts up to level 443 m, preceded by pit2 at to 407 m.The par1 type occurs approximately between levels 407 m to 305 m.

GUM borehole
GUM is the northernmost borehole among the studied ones.It is located 80 km away from PIR borehole at level 495 m (altitude).There is 20 m of soil, and the basalt reaches a total thickness of 405 m.The contact with the Botucatu Formation sandstones occurs at 70 m (altitude).

FR borehole
The FR, the reference borehole, is located 300 km away from the GUM borehole, with its top at 540 m (altitude).As aforementioned, it constitutes a sequence with intercalations between Pitanga and Paranapanema magma types.The limits are marked by vesicular zones and a massive sill between flows.
The presence of such type of intrusion between flow limits is not unprecedented in the literature on volcanic fissural systems.It has already been described for flows in southern PCFB by Peate & Hawkesworth (1996) and in Iceland by Gudmundsson et al. (2014).
No sandstone intercalations occur and the sequence begins at level 500 m after 40 m of sandstone associated with the Bauru Sub-Basin.The borehole ends at -320 m in contact with the Botucatu Formation sandstone (Fig. 13A).
From top towards bottom, the sequence of lavas (level -196 m); par2 (level -260 m).The sill occurs at level 340 m between the last two sequences of flows, and is less than 1 m thick.
Two samples from the sill located between pit1 and par1, at the top of the sequence, have P 2 O 5 concentrations (0.40 and 0.41 wt.%; Fig. 13A) significantly different from those observed in pit1 and par1.The concentrations of TiO 2 (0.39 and 0.40 wt.%, Fig 13 .b)and Zr/Y ratios (6.2 and 6.3, Fig. 13D) are close to those of pit1 and pit2, while Sr is close only to pit1 (Fig. 13C).
Three surface samples were also shown in Figs.

B borehole
B borehole is 120 km west of FR borehole, with its top at level 370 m (altitude), and also features a wide variety of flows covered by 8 meters of Bauru Sub-Basin sandstones.The basalt occurs from level 350 m down to -430 m (780 m-thick), at a concordant contact with the Botucatu Formation (Fig. 14A).
There are no intercalations with contemporary sandstone.Flows are massive with vesicular levels at the upper part, and similarly to the FR borehole, there is also a massive sill, less than 4 m thick, close to the -236 m level between two flows (pit1 and par1).

ARA borehole
The ARA borehole is located at level 342 m (altitude), comprising the fewest samples when compared to others, and 27 km away from B borehole.The sequence of basaltic flows ends at level 308 m and is overlain by 34 m of Bauru Sub-Basin sandstone.It starts at the -613 m level in contact with the Botucatu Formation sandstone (Fig. 15A).
Brazilian Journal of Geology, 48( 2 This borehole has already been used by Torres et al. (2008) and is of great importance for the continuity of the flows related to the closest boreholes, despite the small number of suitable samples for geochemical analysis.
Thus, the sequence bears par1 flows up to approximately 110 m, preceded by a sample of pit1 (28 m), pit3 (-112 m), par1 (-368 m), pit1 (-562 m) and a pit4 sample close to the contact (-600 m).Because of the small number of samples, thicknesses of flow sequences are approximated and supported by the petrographic analyzes of the chips.

IVI borehole
The IVI borehole is located in the westernmost part of the area, close to the border of Paraná Basin, 360 km from the ARA borehole, crossing the Paraná River (Fig. 1).The contact with the basalt occurs at 190 m after 193 m of Bauru Sub-Basin sandstone, while the contact with the Botucatu Formation sandstone occurs at -563 m (Fig. 16A), adding 753 m of lava flows.
Intercalation with contemporary sandstone is observed at level -350 m (Fig. 16A), less than 2 m thick, separating two pit4 and pit1 flow sequences.All sequences were identified by vesicular layers.

DISCUSSION AND CONCLUSIONS
The different sequences of Paranapanema and Pitanga magma types (par1, par2, pit1-pit5), identified by the contents of P 2 O 5 , TiO 2 and Sr, and Zr/Y ratios, can be related to transitional compositions between the two magmas.This interpretation is also shown in Cox and Hawkesworth (1985), who described similar problems in defining the boundaries of different magma types in the Deccan.
The data obtained reinforce previous studies (e.g.Piccirillo & Melfi 1988, Marques et al. 1999, Rocha-Júnior et al. 2013), which demonstrated the cogeneticity (common mantle source) between the two high titanium magmas (Pitanga and Paranapanema).Distribution patterns for REE and other incompatible trace elements showed that there are not significant geochemical differences between Pitanga (pit) and Paranapanema (par) flows.Due to the strong depletion in incompatible elements, the Ribeira (rib) data are in perfect agreement to those reported by Machado et al. (2015) for this type of magma.
Urubici (uru) flows have a multi-element distribution pattern different from the other ones.Despite the small number (2) of samples, it is interesting to note that they are also different from the Urubici basalts that occur in southern PCFB province, where the incompatible element ratios are always higher than Pitanga type (Peate et al. 1999).This distinction could mean that the Urubici from the northern Paraná Continental Flood Basalt has a different origin and/or magmatic evolution when compared to southern Paraná, even though presenting similar contents of TiO 2 , P 2 O 5 and Sr, and Ti/Y ratios.Some of these geochemical differences can be associated with the mineralogy of such rocks, including the presence of primary sulfide (pyrite, arsenopyrite) up to 1%, also verified in the surface samples in Machado et al. (2007).Although the LILE enrichment of these rocks could be the result of crustal material involvement (Faure 2001), the Rb/Ba ratio (0.5-0.9), normalized to the primitive mantle, does not corroborate this possibility.
The data also show that, among the identified Pitanga and Paranapanema magma types, pit5 is the most enriched in P 2 O 5 , followed by pit4, pit1, pit2, pit3 and par1.It is important to emphasize that the thickest boreholes tend to present the largest diversity in P 2 O 5 , leading to the conclusion of possible multiple effusion points.Low-pressure magmatic differentiation processes may also have contributed to generating additional geochemical variability of the observed flows.
The maximum thicknesses calculated for the different flow sequences are:   It is remarkable that at the bottom of IVI, ARA and B boreholes there is a conspicuous presence of pit4 and pit1 flows for a distance as long as 387 km.Moreover, the B borehole has a sequence of par1 and pit1 that repeats in the FR borehole at the same level.This shows a good correlation between the ARA, B and FR boreholes in a distance of 147 km.Individual basaltic flows extending for more than 100 km (in length) in flat areas, such as those of Columbia River, are reported in the literature (Self et al. 1997).For the investigated area, the data indicate that, between the ARA and B boreholes, the flows run for at least 27 km, if single eruption sites are admitted for effusive points to the observed sequences.
The sequence of flows in the GUM borehole (pit1, pit4, and pit3), excluding the Urubici magma type, is repeated in the intermediate portion of the FR borehole between level 50 to -100 m for a distance of approximately 300 km.
Pit1 flows occur at the top of the ALT, PAI, PIR and FR boreholes, but have no equivalent in the lower sequences.On the other hand, pit5 is restricted to the ALT borehole, representing the basalts with the highest MgO investigated in the present study (Fig. 5).The Paranapanema (par1) flows would necessarily require multiple sources, possibly in the region around ARA, IVI and B boreholes, to reach distances greater than 100 km with decreasing thickening towards the east (FR borehole).The same happens with the pit3 flows, with great thickening in the ARA and B borehole and thinning towards GUM, located 400 km far.
Finally, considering only the magma types, it was possible to elaborate a schematic stratigraphic geochemical column.Fig. 18 shows a schematic chemostratigraphic column for the northern part of the Paraná Continental Flood Basalt province, considering only the magma types.Urubici and par2 are omitted in view of their restricted occurrence in the area.Pitanga are the most voluminous, reaching maximum thicknesses of 750 m, and are succeeded by up to 250 m of Paranapanema.Ribeira has much lower thicknesses (< 30 m) and is found at the basis of the sequence, directly over the Botucatu Formation (Machado et al. 2014) and, in the case of borehole IVI, intercalated with par1, above some Pitanga flows.

Figure 2 .
Figure 2. Isopach map of PCFB in Paraná Basin, with contours at 200 m intervals(Zalán et al. 1987).Symbols mark the location of the studied boreholes and surface samples (KS) and will be used hereafter.

Figure 3 .
Figure 3. Sample chips selection of PAI 238 for geochemical analysis.

Figure 4 .
Figure 4. Classification and nomenclature of the analyzed samples, according to the TAS scheme (Le Maitre 2002).The upper dashed line separates alkaline and sub-alkaline fields (from Miyashiro 1978), whereas the lower dashed line shows the dividing line of MacDonald and Katsura (1964) for Hawaiian basalts.Oxides in wt.% recalculated to 100% on a volatile-free basis.Symbols as shown in Fig. 2.

Figure 5 .
Figure 5. Major element variation versus MgO content for the investigated basalts.Values in wt.%.Total iron as Fe 2 O 3(t) .Symbols as shown in Fig. 2.

Figure 7 .
Figure 7. REE chondrite-normalized diagram with a range of concentration for the studied basalts associated with the Pitanga and Paranapanema rocks from the Paraná CFB.Normalization factors are from the primordial mantle of McDonough and Sun (1995).

Figure 8 .
Figure 8. Diagram of incompatible elements normalized to primordial mantle (McDonough & Sun 1995).The ratios presented in the diagram show the average of the representative samples calculated for the Pitanga and Paranapanema flows, while Urubici and Ribeira correspond to samples GUM123 and KS841 respectively.
■ pit1: 200 m; ■ pit2: 100 m; ■ pit3: 350 m; ■ pit4: 150 m; ■ pit5: 50 m; ■ par1: 250m; ■ par2: 170 m; ■ uru1: 30 m; ■ rib1: 150 m.The total pile of Pitanga magma type has a maximum thickness of 850 m, while the Paranapanema, Ribeira and Urubici have, 420, 150 and 30 m, respectively.The thicknesses increase towards the depocenter of the Paraná Basin (Fig. 17 ).The Paranapanema magma occurs at the upper part of the central portion of the Paraná Basin in the B, ARA and IVI boreholes.The last one also presents interbedded Ribeira type (Fig.17 ).The B and ARA boreholes have a similarity in the sequence of flows that suggest an erosion of at least 170 m of Paranapema type in the region nearby the B borehole.Still concerning the Paranapanema magma type, the par1 flows occur between -400 m to -200 m, with a similar a thickness in the ARA, B, and FR boreholes, whereas Par2 flows are only found at the bottom of the FR borehole.

Figure 17 .
Figure 17.Compilation of identified flow sequences from the 10 studied boreholes.The data show good correlations between FR, B, ARA and IVI boreholes.Colors are the same used in Tab. 2.

Table 2 .
Major and trace element data for representative samples of basalts in the study area.Total iron as Fe 2 O 3(t) .Oxides in wt.%.Elements in ppm.