Abstract

The Asunción-Sapucai-Villarrica graben (ASV) in Eastern Paraguay at the westernmost part of the Paraná Basin was the site of intense magmatic activity in Mesozoic and Tertiary times. Geological, petrological, mineralogical, and geochemical results indicate that the following magmatic events are dominant in the area: (1) tholeiitic basalt and basaltic andesites, flows and sills of low- and high-titanium types; (2) K-alkaline magmatism, where two suites are distinguished, that is, basanite to phonolite and alkali basalt to trachyte and their intrusive analogues; (3) ankaratrite to phonolite with strong Na-alkaline affinity, where mantle xenoliths in ultramafic rocks are high- and low-potassium suites, respectively. The structural and geophysical data show extensional characteristics for ASV. On the whole, the geochemical features imply different mantle sources, consistently with Sr-Nd isotopes that are Rb-Nd enriched and depleted for the potassic and sodic rocks, respectively. Nd model ages suggest that some notional distinct “metasomatic events” may have occurred during Paleoproterozoic to Neoproterozoic times as precursor to the alkaline and tholeiitic magmas. It seems, therefore, that the genesis of the ASV magmatism is dominated by a lithospheric mantle, characterized by small-scale heterogeneity.

1. Introduction

Velázquez et al. [1] presented a structural analysis of the central segment of the “Asunción Rift,” mainly based on the previous papers related to the Eastern Paraguay magmatism in general, and to the Asunción-Sapucai-Villarrica graben (ASV) in particular, and also based on extensive field data collected earlier on the dyke swarms cropping out in the area. However, some aspects as, for example, the close association in space of potassic and sodic alkaline rock-types with tholeiitic dykes and flows (both of high- and low-Ti types; cf. [2]), have not been discussed in detail by the above authors.

This paper focuses on general aspects of the magmatism from the western margin of the Paraná-Angola-Etendeka system (PAE), where tholeiitic flows and dykes (Early Cretaceous, both of high-Ti and low-Ti types; cf. Figure 1(a) and [3]) are associated in time and space with a wide variety of alkaline rock-types (both potassic and sodic) and carbonatites.

The investigated alkaline and tholeiitic rocks from Eastern Paraguay span in age mainly from Early Cretaceous to Paleogene times. Therefore, they are germane to the magmatic and tectonic evolution of the PAE and of the Atlantic Ocean (cf. Figure 1(b) and [4, 5]).

This review aims at discussing (1) the relationships between tectonics and magmatic activity on the basis of field and geophysical data; (2) a detailed description of the geological characteristics of the Asunción-Sapucai-Villarrica (ASV) graben; (3) petrography and petrochemistry of the magmatic rock-types (and the associated mantle xenoliths) from ASV; (4) the most important geochemical and Sr-Nd isotopic features of the magmatism; (5) the petrogenesis and geodynamic implications.

2. Eastern Paraguay: Geological Outlines and Geophysical Aspects

Eastern Paraguay shows a complex block-faulted structure between the southern Precambrian tip of the Amazonian Craton (Apa block) and the northern one of the Rio de la Plata Craton (Caacupú high; Figure 2). This intercratonic area, including the westernmost fringe of the Paraná Basin (PB), represents an undeformed basin along the western Gondwana (Figure 1), where the sedimentation, started in Cambrian times, was topped by Early Cretaceous tholeiitic basalts of the Serra Geral Formation [8, 9] and followed by younger sedimentation (Figure 2).

Moreover, it should be also noted that Eastern Paraguay lies along the former western margin of the Gondwana, bounded by an anticlinal structure established since Early Paleozoic, the Asunción Arch, which separates the Paraná Basin (East) from the Gran Chaco Basin (West) [10]. PB shows a high-velocity upper mantle “lid” with a maximum S-wave velocity of 4.7 km/s (Moho 37 km depth), with an unresolved low-velocity zone to a depth of at least 200 km [11]. Between the two blocks, that is, Apa and Caacupú of Figure 2, Eastern Paraguay was subjected to NE-SW-trending crustal extension during Late Jurassic-Early Cretaceous, probably related to the western Gondwana breakup (Figure 1; [10, 12] and therein references). NW-SE trending faults, paralleling the dominant orientation of Mesozoic alkaline and tholeiitic dykes, reflect this type of structure [13, 14].

From the aeromagnetic survey [15, 16], the linear magnetic anomalies trend mainly N40-45W (Figure 3). These anomalies have been interpreted as resulting from Early Cretaceous tholeiitic dyke swarms [16], but field evidence does not support this hypothesis. Most magnetic anomalies correspond to older Precambrian tectonic lineament [17]. The Landsat lineaments, trending mainly NE and EW, may therefore reflect tectonic lineaments from the basement [18].

The Bouguer gravity map (Figure 3) consists of prevailing NW-trending gravity “highs” and “lows” that represent shallow to exposed basement and sedimentary basins, respectively [19]. The boundaries between gravity highs and lows are generally marked by step gradients that may reflect abrupt basement offsets along faults, or basement dip changes by crustal warping. Notably, the gravity lows and highs parallel the dominant NW attitude of the magnetic lineaments, with post-Palaeozoic magmatism (both tholeiitic and alkaline) associated with gravity lows. Therefore, the ASV represents the most important geological structure of the region (see later).

Notably, the present seismic activity (earthquakes with hypocenters <70 km) indicates that the NW-trending fault systems continue up to present day, that is, the Pilcomayo lineament (inset of Figure 3), also aligned with the Piquirí lineament and the Torres syncline (Figure 1(b)).

In conclusion, the resulting structural pattern controlled the development of grabens or semigrabens as a response to NE-SW-directed extension and continued evolving into Cenozoic times [12, 24]. According to Tommasi and Vauchez [30], rift orientations seem to have been controlled by the preexisting lithospheric mantle fabric, as indicated by deep geophysical data.

The thermal history, using apatite fission track analyses (AFTA), reveals that at least two main episodes have been identified in sedimentary and igneous/metamorphic samples ranging in age from Late Ordovician to Early Cretaceous [31, 32]. AFTA data from ASV show evidence for rapid cooling beginning sometime between 90 and 80 Ma, similar to the results found along Brazilian and Uruguayan coasts. A Tertiary heating/cooling episode is also suggested by AFTA data (60–10 Ma), supporting early work in the area. The time of the first event is significantly younger than any rifting activity related to the Paraná flood basalts and to the opening of the South Atlantic. Late Cretaceous cooling may have involved several kilometers of differential uplift and erosion and would have played an important role in the control of the geomorphology and drainage patterns of the region, especially in the ASV system [31, 32].

Beginning from Mesozoic times (Figure 2), almost five main alkaline magmatic events have occurred in Eastern Paraguay, other than the Early Cretaceous tholeiitic magmatism (133-134 Ma; [4]). Three of these include rocks of sodic affinity, corresponding geographically to the provinces of Alto Paraguay (241.5 ± 1.3 Ma), Misiones (118.3 ± 1.6 Ma), and Asunción (58.7 ± 2.4 Ma), whereas two involve rocks of potassic affinity associated with the Apa and Amambay provinces, both of similar age (138.9 ± 0.7), and with the Central Province (ASV, 126.4 ± 0.4 Ma).

3. Asunción-Sapucai-Villarrica Graben (ASV)

The ASV graben (Figure 4) is the dominant rifting structure in Eastern Paraguay, linked to NE-SW extensional vectors [13], and it is characterized by a gravity anomaly trending N30-45W near Asunción townships (Figure 3). The graben, nearly symmetrical and defined by major faults along each margin [1, 12, 18, 37], may be subdivided into three main segments (Figure 4): (1) the northwestern part, between the Asunción and Paraguarí townships, indicated also by the intragraben Ypacaray rifting and by the presence of Na-alkaline mafic-ultramafic plugs carrying mantle xenoliths (Asunción Province; age between 66 and 39 Ma; cf. Figure 4); (2) a central area, defined between Paraguarí and La Colmena townships, characterized by K-alkaline complexes and dykes (average age 127 Ma; cf. Figures 3 and 4) and by subordinate Na-alkaline complexes and dykes (age around 66–60 Ma); (3) the eastern part, between La Colmena and the Ybytyruzú hills, characterized by K-alkaline complexes and dykes (age around 130–125 Ma) and by tholeiitic lavas and flows (age around 133–130 Ma; cf. [4]).

Geology and gravity results indicate that the ASV graben extends up to 100 km into the Chaco Basin (Figures 2 and 3) and turns to an N80W trend between Paraguarí and Villarrica townships, in the area marked by a gravity low anomaly that increases eastwards (Figures 3 and 4). Paleomagnetic data [38] revealed that the potassic Cretaceous rocks (75 samples) have reversal polarity remnants corresponding to a paleomagnetic pole at 62.3°E and 85.4°S, similar to that calculated for the Serra Geral Formation in the Paraná Basin. It is therefore inferred that ASV K-alkaline rocks were emplaced, while the tholeiitic magmatic activity was still taking place (Figure 1).

Tertiary sodic rocks, carrying mantle xenoliths, mainly occur in an area characterized by relatively gravity highs (Figure 3), whereas the potassic complexes and dykes occur in a gravity low belt characterized by block faulting (Figures 3 and 4). Finally, an NS-trending fault system bounds the western side of the Ybytyruzú hills, close to several alkaline intrusions [24] East of Villarrica.

A total of 527 samples from the ASV area, including 220 dykes, have been investigated in detail, including field relationships, age, orientation (dykes), petrography, and mineral chemistry of the rocks reported in [10, 13, 3336]. 55% of the dykes are preferentially oriented N20-60 W (Figure 5, (III)).

4. Classification, Petrography, and Petrochemistry of the Magmatic Rocks

The fine-grained texture of the ASV rocks, including the intrusive variants, is plotted on the De La Roche chemical classification [39] diagram (Figure 5), with an additional subdivision based on wt% K2O and Na2O on Le Maitre’s [40] diagram (inset (I)).

4% of the dykes have tholeiitic affinity, sills, lava flows, and dykes, both high- and low-Ti variants, according to [2, 43] and straddle the fields of the olivine tholeiitic basalts and andesite-basalts, whereas 17% are sodic and 79% potassic in composition. The distribution of the potassic rocks shows that the variation of the dykes is consistent with that of the associated intrusive and volcanic rocks (inset (I)). Two main lineages are apparent for the potassic rock-types (Figure 5 and inset (II)): (1) a silica undersatured lineage (B-P) ranging from basanite to phonolite and peralkaline phonolite and (2) a silica-saturated lineage (AB-T) ranging from alkali basalt to trachyphonolite and trachyte [13, 42, 44].

In summary, tholeiitic rock-types (both high-titanium, H-Ti, and low-titanium, L-Ti, variants), K-alkaline, and Na-alkaline rocks are widespread in a relatively narrow area represented by ASV graben (ca. 35 × 200 km; cf. Figure 3). In Tables 1 to 4, representative and average analyses of these rocks are reported.

4.1. Tholeiites

The tholeiitic rocks are mainly lava flows (Ybytyruzú hills) and gabbro sills usually occurring near the Ybycuí township (Figure 4). The textures are porphyritic to aphanitic (lava flows and dykes) and equi- to subequigranular (sills), characterized by the presence of two clinopyroxenes (augite and pigeonite) and by pronounced variations of TiO2 (1.4 to 3.2 wt%; cf. Table 1) and incompatible elements (IE: e.g., Ba ca. 150 to 400 ppm, Zr ca. 100 to 200 ppm). These variations are similar to those common to tholeiitic basalts of the Paraná Basin that are dominated by low and high contents of TiO2 (L-Ti and H-Ti variants: southern and northern Paraná Basin with TiO2 <2 and >2.5 wt%, resp.; cf. Figure 1) and IE [2, 45, 46]. The mineral assemblages (augite, pigeonite, olivine, plagioclase, magnetite, and ilmenite) show crystallization temperatures around 1200–1000°C and fO2 conditions between NNO and QMF buffers [3].

Mass balance calculations and the remarkable differences in the ratios of the IE between L-Ti and H-Ti tholeiites suggest that such differences cannot be explained by fractional crystallization, or by means of melting and zone refining processes of a homogeneous mantle source. A basalt genesis from chemically different lithospheric mantle materials is therefore indicated. Mixing processes seem to have played a minor role [2, 47].

4.2. Potassic Suites

The different lithotypes of the potassic rocks (cf. Figure 5 and Tables 2 and 3) maybe subdivided into B-P and AB-T lineages, a distinction consistent with crystal fractionation models [13].

4.2.1. B-P Suite

The effusive rock-types and dykes (basanites, tephrites, and phonotephrites) typically show porphyritic textures with phenocrysts of clinopyroxene , olivine , and leucite pseudomorphed by sanidine + nepheline in a glassy groundmass showing microlites of clinopyroxene ± olivine, Ti-magnetite ± ilmenite, Ti-phlogopite-biotite, alkali feldspar , and nepheline-analcime . Phenocrystal plagioclase (up to An74) is also present. Accessory phases are amphibole (pargasite-kaersutite) and apatite ± zircon. Phonolites are characterized by phenocrysts of leucite-pseudomorphs, alkali feldspar , clinopyroxene , Fe-pargasite, and nepheline ± biotite ± titanite ± melanite (Ti-andradite up to 68 wt%) ± magnetite or haematite. Glassy groundmass contains microlites of alkali feldspar, nepheline, and clinopyroxene ± Ti-andradite ± magnetite. Haematite is found replacing magnetite or as a common groundmass constituent in peralkaline phonolites.

The intrusive analogues (theralites, essexitic gabbros, ijolites and essexites) show a holocrystalline inequigranular texture, with diopsidic pyroxene , olivine to , mica (Ti-phlogopite to Ti-biotite), Ti-magnetite, alkali feldspar, and nepheline   ± leucite ± amphibole. “Intrusive” leucite and leucite pseudomorphed by analcime and plagioclase are common, testifying to subvolcanic conditions [24]. Cumulitic textures are represented by clinopyroxene with alkali feldspar + nepheline + carbonate as intercumulus phases.

On the whole, the temperatures calculated for the mineral assemblages span from 1250 to 1200°C and from 1130 to 1080°C, pheno- and microphenocrysts, respectively, under about 1 kb H2O pressure and around 1000°C under atmospheric pressure; lower temperatures, for example, 900–700°C, suggest possible subsolidus exchange reactions (e.g., albitization). Magnetite-ilmenite pairs define temperatures between 1000 and 1100°C, along the NNO buffer [12].

Notably, the B-P suite is also characterized by carbonatitic magmatism, although it is very limited in volume. A silicobeforsitic flow is present in the Sapucai lava sequences and as “ocelli” at the Canãda and Cerro E. Santa Helena K-alkaline complexes [48]. These primary carbonates are relevant as they reveal a CO2 imprinting after the tholeiitic magmatism.

4.2.2. AB-T Suite

Alkali gabbros, syenogabbros, and syenodiorites are usually porphyritic and seriate in texture. They contain clinopyroxene , olivine , Ti-biotite, Ti-magnetite ± ilmenite, plagioclase , alkali feldspar , interstitial nepheline-analcime , and alkali feldspar . Accessories are apatite ± amphibole ± titanite ± zircon. Nepheline syenites and syenites are texturally equi- to subequigranular and seriate. The rock-types are characterized by alkali feldspar , clinopyroxene , and nepheline (Ne85) and hastingsite. Common accessories include titanite and apatite ± carbonate ± zircon.

Alkali basalts, trachybasalts, and trachyandesites are porphyritic rocks exhibiting phenocrysts and/or microphenocrysts of clinopyroxene , olivine , plagioclase , magnetite, and biotite set in a glassy groundmasses consisting of microlites of clinopyroxene , magnetite, ilmenite, biotite, plagioclase , alkali feldspar , nepheline-analcime , amphibole, and apatite ± titanite ± zircon.

Trachyphonolites and trachytes are porphyritic to aphyric. The phenocrysts are alkali feldspar   ± clinopyroxene   ± plagioclase , pseudomorphosed leucite, amphibole, and biotite in a hypocrystalline to glassy groundmass containing microlites of alkali feldspar, biotite, and clinopyroxene ± biotite ± amphibole ± magnetite ± Ti-andradite ± haematite.

Notably, the equilibration temperatures of the mineral assemblages are very similar to those suggested for the B-P suite [10, 13, 49].

The differentiation history of the potassic ASV magmas remained largely into the stability fields of the ferromagnesian phases [12], that is, olivine + clinopyroxene, in the system Mg2SiO4-KAlSiO4 [10, 13, 51], suggesting that relatively high temperature subvolcanic conditions were prevalent. Moreover, the extreme differentiates from both B-P and AB-T suites approaches to the composition of peralkaline residua, pointing to an extensive subvolcanic crystallization of K-rich, aluminous phases, likely phlogopite and leucite, at least for the B-P suite [52].

4.3. The Sodic Rocks

These consist mainly of ankaratrites + (mela)nephelinites (45%) and phonolites (42%) (Comin-Chiaramonti et al. [10, 13, 44]). Representative chemical analyses are reported in Table 4.

Phenocrysts-microphenocrysts in ankaratrites and nephelinites are characterized by olivine (phenocrysts 1–7 vol%, mole%; 1–6 vol% microphenocrysts, mole%), clinopyroxene (phenocrysts 1–6 vol%, Mg# ~0.8), titanomagnetite (0.3–0.7 vol%, up to 38 ulv. mole%), and occasionally phlogopite microphenocrysts. The hypocrystalline groundmass contains clinopyroxene (39–46 vol%, Mg# ~0.75), olivine (3–6 vol%, mole%), titanomagnetite (4–7 vol%, up to 43% ulv. mole%), nepheline (16–21 vol%), and glass (11–25 vol%) [53].

Phonolites are typically microphyric to hypocrystalline with alkali feldspar phenocrysts or microphenocrysts , nepheline , occasionally altered to cancrinite, acmitic clinopyroxene (acmite up to 63 wt%), and ferroedenite-ferropargasite amphibole. Haüyne, mica, haematite, and zircon are typical accessory minerals. Analcime occurs as a nepheline pseudomorph. Ti-andradite and titanite phenomicrophenocrysts are present in some dykes [10].

Lherzolite, harzburgite, dunite mantle xenoliths, and xenocrystic debris are common and abundant in the ankaratrites and nephelinites [54] from eastern ASV (Figure 3).

The generation of the ASV nephelinites may be modelled for 4–6% degrees of partial melting, with crystallization temperatures between 1200 and 1000°C from garnet peridotite sources [10, 55].

4.4. Mantle Xenoliths

The xenoliths are mainly spinel-lherzolites, harzburgites, and subordinate dunites. The dominant texture is protogranular, rarely tabular or porphyroclastic [54]. The Paraguay mantle xenoliths contain variable amounts of glassy patches (blebs) and glassy drops in clinopyroxenes [56]. The latter show an overprinted spongy texture [54]. The blebs (whole composition: Mg# 0.88–0.91) consist mainly of a glassy matrix having microlites of olivine (Mg# 0.91-0.92), clinopyroxene (Mg# 0.91–0.93), Cr-spinel, and (rarely) phlogopite (Mg# 0.86–0.92). According to [56, 57], they were formed by decompression melting of amphibole and phlogopite.

The Paraguayan xenoliths are characterized by a large range of K2O (0.02 to 0.51 wt%). Some xenoliths have K2O abundances comparable or even higher than those reported for metasomatized mantle peridotites [58], resembling in some cases to amphibole-mica-apatite-bearing mantle-xenolith suites [59]. K2O contents and the abundance of blebs and glassy drops allow to group the Paraguay xenoliths into two main suites, that is, a low-K suite (LK, K2O <0.15 wt%) and a high-K suite, the latter with abundant glassy drops and/or variable amounts of blebs and spongy clinopyroxenes (HK, K2O ≥0.2 wt%). Representative analyses are listed in Table 5, and complete sets of chemical analyses are in [54].

Based on Spera’s [60] approach, Comin-Chiaramonti et al. [55] suggested that the ascent of the xenoliths to the surface took place in a very short time, for example, less than nine days (assuming a diameter of 45 cm, corresponding to the largest xenolith size, a density of 3.3 g/cm3, and an origin of the hosting liquids from depth ~70–75 km, i.e., ~boundary between garnet and spinel peridotite). These mantle xenoliths provide a direct sampling of the subcontinental mantle in ASV along with the metasomatic processes [61, 62].

Coherent variations of major element in both suites follow a dunite-lherzolite sequence trending to the mantle composition [54]. The population is represented mainly by lherzolitic-harzburgitic (dunitic) compositions, mostly with a 0.55–0.63 range of (SiO2 + Al2O3)/(MgO + FeOt) molar ratio. The residual character of the harzburgitic mantle xenoliths (believed to be consistent with melting and basalt-component removal) is also indicated by the decrease in the cpx/opx modal ratio with decreasing modal cpx, which fits the model variation trend induced by partial melting of lherzolite [56, 63].

Equilibration temperatures for orthopyroxene-clinopyroxene pairs [64] and for olivine-spinel pairs [65] vary between 862° and 1075°C and between 748° and 968°C, respectively. Intracrystalline temperatures [6668] for clinopyroxene and orthopyroxene pairs vary between 936° and 1033°C and between 920 and 1120°C (LK and HK suites, resp., with pressures in the range of 1.1–2.3 GPa for both suites, with the higher values related to the more depleted xenoliths from the LK suite [53, 56]).

The oxygen isotope compositions on separates of clinopyroxene and coexisting olivine (δ18O‰) vary from 5.5 to 6.0‰ and from 5.0 to 6.1‰, respectively [53]. These isotopic ratios are in the range for worldwide mantle phases: olivine, 4.4 to 7.5‰, and clinopyroxene, 4.8 to 6.7‰  [69, 70], and for South America mantle xenoliths: olivine, 4.9 to 6.4‰, and clinopyroxene, 5.0 to 6.0‰  [71]. The calculated Cpx-Ol isotopic temperatures [53, 71] are around 970–1070°C, 1030–1130°C (LK), and 1100–1180°C (HK), suggesting equilibration temperatures higher in the HK suite than in LK suite.

The highest temperature of the HK suite is substantially due to higher δ18O‰ in olivine than δ18O‰ in clinopyroxene. According to [72], the relatively low δ18O in the pyroxenes reflects metasomatism by a silicate melt from subducted altered oceanic crust. Therefore, the Paraguayan HK suite would be for some geologists the best candidate for a subduction-related environment.

Notably, the ASV clinopyroxenes [73] have and sites intermediate between those of plagioclase- and garnet-bearing mantle peridotites, that is, in a pressure range between 1.2 and 2.2 GPa. Thus, isotopic and crystallographic results are internally consistent.

On the whole, the ASV xenoliths define a geotherm [53] which, starting at about 830°C and 1.0 GPa, at the transition mantle-crust, intersects the hydrous peridotite solidus at about 1140°C and 2.1-2.2 GPa in the spinel peridotite facies, near to the transition with the garnet peridotite facies, believed to be the source of the sodic alkaline magmatism.

5. Geochemistry

Incompatible elements (i.e., large ion lithophile elements, IE, LILE, and high field strength elements (HFSE)), considered with the Sr-Nd isotopic composition, indicate that the ASV magmatic was generated from geochemically distinct (enriched versus depleted) mantle sources (cf. [75]).

5.1. Incompatible Elements

Mantle-normalized IE patterns for the various and different magmatic rock groups are represented in Figure 6.

The L-Ti and H-Ti ASV tholeiites are distinct in terms of their relatively low elemental abundances and high LILE/HFSE ratios. In particular, their marked Ta-Nb negative spikes are similar to those of the potassic alkaline magmas from ASV, but clearly different in comparison to the Cenozoic sodic alkaline rocks from the same area.

The suites are similarly enriched in REE and exhibit steep, subparallel LREE trends ( = 26–161, 17–62 and 11–46 for B-P, AB-T and Na rocks, resp.), which tend to flatten out for HREE ( = 1.24–1.96, 1.09–2.00 and 0.56–2.05 for B-P, AB-T and Na rocks, resp.); HK-dykes are excepted. REE profiles with LREE enrichment and flat HREE suggest mantle sources previously depleted by melt extraction and subsequently enriched [76].

Multielemental diagrams, normalized to a primordial mantle composition (Figure 6) show a substantial overlap of B-P and AB-T compositions, negative Nb, Ta, P, Ti, and Y spikes, and positive U, K, and Sr anomalies. In general, B-P compositions are higher in Rb, K, Zr, Hf, Ti, Y, LREE, and MREE than AB-T compositions. On the contrary, the Na rocks yielded La/Nb and La/Ta ratios close to unity, respectively, and Nb/K and Ta/K ratios greater than 1.0, respectively. Notably, the variations of incompatible elements of the AB-T compositions mimic to some extent those of the ASV H-Ti (and L-Ti) tholeiitic basalts, which approach the lower (Rb to Nd) and higher (Sm to Lu) elemental concentrations of that suite (Figure 6).

In summary, substantial overlap in bulk-rock chemistry exists between the investigated B-P and AB-T rocks, both characterized by variable K/Na ratio, the K types being dominant. REE and other incompatible elements show similar concentration levels and variation trends in the two suites. The mantle normalized incompatible element patterns of both ASV suites show strong affinities, including negative “Ta-Nb-Ti anomalies,” with the Paraná tholeiites [10].

On the whole, the geochemical features suggest that the enrichment processes were related to small-volume melts in a lithospheric mantle [4, 10, 74].

5.2. Sr-Nd Isotopes

The investigated rocks from ASV show a large distribution of Sr-Nd isotopic compositions (Figure 7), delineating a trend similar to the “low Nd” array of [77] (cf. also “Paraguay array” of [4]). Due to the high Sr and Nd of the most “primitive” alkaline rocks and associated carbonatites from Eastern Paraguay, Comin-Chiaramonti et al. [12] suggested that initial Sr-Nd isotopic ratios of these rocks can be considered uncontaminated by the crust and, as a result, representative of the isotopic composition of their mantle source(s). The ASV potassic rocks have the highest (time integrated) Sri and the lowest Ndi. Including the carbonatites of the BP-suite [24, 78], the Sri and Ndi range from 0.70645 to 0.70716 and from 0.51151 to 0.51179, respectively (Tables 2 and 3). These values are quite distinct from those of the ASV Paleocene sodic rocks (ca. 60 Ma), which plot within the depleted quadrant (Sri = 0.70364–0.70391, Ndi = 0.51264–0.51273; cf. Table 4), towards the HIMU-DMM depleted mantle components (Figure 7). Notably, Sri and Ndi of the tholeiites, believed to be uncontaminated [10], both H- and L-Ti, are intermediate between the potassic and sodic rocks: 0.70456–0.70576 and 0.51225–0.51256, respectively (cf. Table 4). These values approach to the range of the Early Cretaceous uncontaminated tholeiites from the Paraná Basin [47], that is, Sri = 0.70527 ± 0.00051 and Ndi = 0.51264 ± 0.00011. To be stressed that the genesis of these tholeiites requires lithospheric mantle components, as indicated by K-alkaline and carbonatitic rocks from ASV [4, 5].

Considering the whole Paraná-Agola-Etendeka system (PAE), and that the ASV is located at the central westernmost side of the PAE [5, 6], the different geochemical behaviour in the different PAE sectors implies also different sources. Utilizing the (Nd) model ages [78, 79], it should be noted that (1) the H-Ti flood tholeiites and dyke swarms from the Paraná Basin and the Early Cretaceous potassic rocks and carbonatites from Eastern Paraguay range mainly from 0.9 to 2.1 Ga, whereas in Angola and Namibia the Early Cretaceous K-alkaline rocks vary from 0.4 to 0.9 Ga; (2) the L-Ti tholeiites display a major variation, from 0.7 to 2.4 Ga (mean 1.6 ± 0.3) with an increase of the model ages from North to South; (3) Late Cretaceous alkaline rocks show model ages ranging from 0.6 to 1 Ga, similar to the age span shown by the Triassic to Paleocene sodic alkaline rock-types lying along the Paraguay river [4].

Thus, model ages suggest that some distinct “metasomatic events” may have occurred during Paleoproterozoic to Neoproterozoic times as precursor to the alkaline and tholeiitic magmas in the PAE [78].

The significance of model ages is supported by (1) the isotopic overlapping of different igneous rocks (e.g., H-Ti and L-Ti tholeiites or K-alkaline rocks and carbonatites), which cannot be accidental and points to sampling of ancient reservoirs formed at different times from the same subcontinental upper mantle [79]; (2) whatever the implication, that is, heterogeneity induced by recycled crust in the mantle [80, 81], or occurrence of variably veined material in the subcontinental upper mantle, or both [82, 83], Pb isotope data indicate a mantle source of ca. 1.8 Ga for the Paraná H-Ti tholeiites. Since much of the crust in southern Brazil appears to have been formed at ca. 2 Ga ago [84], it follows that magma genesis involved ancient lithospheric mantle reset at well-defined isotopic ranges. A veined lithospheric mantle (amphibole/phlogopite-carbonate-lherzolite and amphibole-lherzolite + CO2-fluid type III and IV veins of Meen et al. [85] of Proterozoic age) may well account for the magmatism both of the PAE and ASV (Figure 8).

5.3. Mantle Sources

The origin of the ASV magmas is closely related to, and probably constrained by, the geodynamic processes which promoted the generation of the adjacent and coeval magmatism in Brazil [10]. The origin and emplacement of the ASV magmatism occurred before, during, and after the opening of the South Atlantic and appears to be related to early stages of lithospheric extension, associated with an anomalously hot mantle [3, 90, 91]. The K/Na and trace element variations and Sr-Nd isotope characteristics of the ASV suites support the view that the lithospheric mantle played an important role in their genesis as well as that of the Paraná flood basalts in Brazil [47, 92, 93].

It should be noted that the most mafic tholeiitic and potassic ASV rock-types are relatively evolved compared with expected primary compositions (e.g., Mg# >0.60–0.65, Ni >235 ppm), ankaratrite is excepted. Possible ASV primary melts (e.g., Mg# ~0.74) in equilibrium with are expected to have fractionated olivine (Ol) and clinopyroxene (Cpx) at or near the mantle source and certainly during the ascent to the surface [13]. Neglecting the effects of polythermal-polybaric fractionation on the chemistry of the fractionates, crystal fractions have been calculated to restore the selected ASV parental compositions to possible near-primary melts in equilibrium with their mantle sources (cf. Table 3 of [10]).

Mass balance calculations, starting from different garnet-phlogopite peridotites (cf. [94]; whole-rock and mineral compositions in Table 3 of [10] and Table 10 of [43]), indicate that the ASV composition of primary melts can be derived from high and relatively high melting degrees of anhydrous garnet or phlogopite-bearing peridotite, that is, 12 and 30% for tholeiitic basalts (H-Ti and L-Ti, respectively) and 4–11% for the alkaline rock-types (Table 3 of [10]). However, the presence of a relative K enrichment in the ASV potassic and tholeiitic rocks (Figure 6) suggests that a K-bearing phase (e.g., phlogopite) did not represent a residual phase during partial melting. Phlogopite, instead, was probably a residual phase in the mantle source(s) of the ASV sodic rocks. Notably, the parent melts (Figure 6) show high abundances of IE and high LREE/HREE ratios, which require mantle sources enriched in IE before or subcontemporaneously with melting process [94]. Nd model ages approximately indicate that the IE enrichment in the source mantle of the ASV tholeiitic and potassic rocks probably occurred during Paleoproterozoic times (mean 2.03 ± 0.08 Ga), whereas those relative to the sodic rocks would be related to Late Neoproterozoic events (0.57 ± 0.03 Ga). On the other hand, the mantle xenoliths display Nd model ages of 444 ± 11 Ma (cf. Tables 1 to 5).

Concluding, the patterns of the mantle sources of the ASV potassic rocks (along with H-Ti and L-Ti tholeiites) are characterized by negative “Ta-Nb-Ti” and positive Ba and Sm spikes. On the contrary, the patterns of the mantle sources of the ASV sodic rocks (ankaratrites-melanephelinites) show positive Ta-Nb and Zr and negative K and Sm spikes. It seems, therefore, that the genesis of the ASV alkaline magmatism is dominated by a lithospheric mantle, characterized by small-scale heterogeneity, and documented by the occurrence of bleb-like glass in spinel peridotite nodules from the sodic ultramafic rocks [54, 55].

Finally, the isotopic signature of the tholeiitic and K-Na-alkaline-magmatism from the PAE in general, and from ASV in particular, may reflect ancient heterogeneities preserved in the subcontinental lithospheric mantle. As matter of fact, all the geochemical data indicate that the PAE magmatism requires heterogeneous mantle sources, also in terms of radiogenic isotopes [6, 10]: such heterogeneities probably represent metasomatic processes which have occurred mainly at 0.5–1.0 Ga and at 1.5–2.0 Ga (cf. of Tables 14).

6. Geodynamic Implications

Some general considerations on the geodynamics are discussed in [46, 8789]. The ASV geodynamic evolution is related to that of the Western Gondwana (Figure 1), which reflects the regional amalgamation processes and requires an overall understanding of the pre-Gondwanic setting, at least at the time of the Brasiliano cycle, both at the Atlantic and Pacific systems. The Brasiliano cycle was developed diachronically between about 890 and 480 Ma, until the final basement of the South American Platform was attained [101]. This occurred during the Lower Ordovician, where Unrug [102] suggested a mosaic of lithospheric fragments linked by several (accretionary, collisional) Neoproterozoic mobile belts. After amalgamation, the Gondwana supercontinent accumulated Paleozoic and Triassic-Jurassic sediments. Concomitantly, it was continuously and laterally accreted to the western borders by successive orogenic belts, until the Pangea Formation [103, 104]. The main cratonic fragments descending from ancestors of the Pangea, for example, the Amazonic and Rio de La Plata cratons, were revoked, and smaller ancient crustal blocks formed at the Paraguay boundaries [105]. The magmatism was likely driven by the relative extensional regimes derived from the relative movements of the ancient blocks (Figure 9), probably induced by counterclockwise and clockwise movements (North and South, resp.) hinged at about 20° and 24°S [106].

The earthquake distribution and mechanisms [29] with hypocenters >500 km and <70 km, respectively (Figure 9), coincide with the inferred Nazca plate subduction at the Paraguay latitude (Figure 10), where the depth of the lithospheric earthquakes, together with the paleomagnetic results [38, 95, 107], delineates different rotational paths at about 18–24°S (i.e., Chaco-Pantanal Basin). These data indicate extensional subplate tectonics along the Andean system [95], coupled with the lineaments of lithospheric earthquakes versus the Atlantic system that parallels the Ponta Grossa Arch and the lineaments of Rio Paraguay and Rio Pilcomayo, also clearly related to extensional environments [24, 29, 78].

Crucial to the genesis of the ASV magma types is the link with the whole PAE system and with the geodynamic processes responsible for the opening of the South Atlantic. According to Nürberg and Müller [108], the sea-floor spreading in the South Atlantic at ASV latitude started at about 125–127 Ma (Chron M4), but North of the Walvis-Rio Grande ridge (Figure 1(b)) the onset of the oceanic crust would be younger, for example, 113 Ma [109]. Therefore, the Early Cretaceous alkaline and alkaline-carbonatite complexes from ASV (and PAE) are subcoeval with respect to the main flood tholeiites of the Paraná Basin and occurred during the early stages of rifting, before the main continental separation.

The origin of alkaline-carbonatitic magmatism in terms of plate tectonics is currently debated [77, 110] and often considered in terms of models only. In the present context, the different westward angular velocity of the lithospheric fragments in the South American plate, as defined by the “second order plate boundaries” (e.g., the Pylcomaio lineament [18]) as well as the different rotational trends at S19–24° Latitude, would favour decompression and melting at different times of variously metasomatized (wet spot) portion of the lithospheric mantle with variable isotopic signatures [24]. Small-scale metasomatism is indicated by the different suites of mantle xenoliths and significant H-O-C and F variations, expected in the mantle source from the occurrence of the related carbonatites. The latter imply lowering of the solidus and, along with an extensional tectonic, favour the melting of the mantle and magma ascent. This scenario also accounts for the presence of Late Cretaceous to Tertiary sodic magmatism in the PAE and in the Eastern Paraguay, where rifting structures are evident (cf. Figures 2 and 4).

7. Concluding Remarks

Geological and geophysical results show that the Mesozoic-Tertiary block-faulting tectonics in Eastern Paraguay is responsible for NW-trending grabens (ASV and Amambay), fault systems, and fault-controlled basins, for example, Jejui-Aguaray-Guazu, Asunción-Encarnación, San Pedro (cf. Figure 2).

Distinct magmatic events are widespread in the Asunción-Sapucai-Villarrica graben, the major extensional structure in Eastern Paraguay, that is, Early Cretaceous tholeiites and K-alkaline-carbonatitic complexes and dykes, along with Tertiary sodic magmatism. It is expected that beneath a rifted continental area, as ASV, the lithospheric mantle may have been thinned and subjected to high heat flow. This allowed the geotherm to intersect the peridotite solidus in the presence of H2O-CO2-rich fluids and favour the melt formations. Under these conditions, the derived melts may form, dissolving most of the “volatile” phases and, percolating through a porous and deformed rock matrix, invaded the base of the lithosphere and favoured metasomatic processes at different levels in the peridotitic lithosphere with the formation of amphibole and/or phlogopite (“veined mantle” of [82, 83]). Notably, the different tholeiitic (high- and low-Ti), K-alkaline suites (B-P and Ab-T), and sodic magmatism with HK and LK xenoliths are consistent with variously depleted lithospheric mantle at different times, pervasively and locally invaded by metasomatizing fluids and/or melts. Thereafter, the newly formed veins (“enriched component”) and peridotitic matrix (“depleted component”) undervent different isotopic compositions with time, depending on their parent/daughter ratio (cf. Figure 8), testifying heterogeneous mantle sources beneath ASV.

Acknowledgments

This paper is dedicated to the memory of Enzo Michele Piccirillo who was an inspiration to the authors over many years of geological work in the region. P. Comin-Chiara-monti gratefully acknowledges CNR and MURST (Italian Agencies) for the field financial support since 1981 and the Department of Geosciences of the Trieste University for the use of chemical laboratories where all the analyses were performed. Thanks are due to Exxon Co., USDMA, ITAIPU, and YACIRETA (Paraguayan Agencies) for their collaborations and to engineer E. Debernardi (deceased) and geologists D. Orué, Luis Lúcia, and L. A. Martinez. This research was also supported by grants (Procs. 08/3807-4 and 10/50887-3) to C. B. Gomes from the Brazilian agency Fapesp. The Ms benefit of the accurate revision of Keith Bell and Lalou Gwalani.