San José de Guaviare Syenite, Colombia: Repeated Ediacaran intrusions in the northwestern Amazonian Craton

The Neoproterozoic igneous rocks found in the municipality of San José del Guaviare include several isolated plutonic bodies that protrude from the Phanerozoic sedimentary cover in belts aligned NW-SE. Limited to the Guaviare department, these intrusions stretch from the La Lindosa mountain range to the corregimiento El Capricho. These plutonic bodies consist of nepheline syenites, nepheline monzosyenites, nepheline-bearing alkali-feldspar syenites, syenites, quartz-syenites, quartz-alkali-feldspar syenites, syenogranites, and quartz-rich granitoids, which have been grouped and termed the San José de Guaviare Syenite unit (SJGS). The intrusion of the unit occurred in the Ediacaran (604 ± 7 Ma and 620.5 ± 7.5 Ma) by mantle-derived alkaline magmas formed in anorogenic settings, most likely in rift-like stretching zones. The silica-subsaturated magma may have reacted with host rocks at the crust level, producing some silica-saturated igneous rocks, such as syenogranites and quartz-syenites, which are found in the El Capricho and Cerritos bodies.


INTRODucTION
The San José de Guaviare Syenite is located in eastern Colombia, in an area with exposed Precambrian plutonic igneous and metamorphic rocks, which belong to the crystalline basement of the NW Amazonian Craton. This unit forms several isolated bodies aligned NW, some of which are structurally controlled and related locally to Mesoproterozoic metamorphic rocks known as the Guaviare Complex (Complejo Guaviare), which shows contact metamorphism superimposed by the syenite intrusion (Maya et al., 2018).
Previous studies have focused on reporting Paleozoic or Precambrian syenite intrusive rocks located in the mountain range west of San José del Guaviare, Guaviare department, Colombia, belonging to the crystalline basement of the Amazonian Craton (Trumpy, 1943(Trumpy, , 1944Gansser, 1954;Pinson et al., 1962;Galvis et al., 1979). The most detailed early studies were performed by Vesga and Castillo (1972), followed by Arango et al. (2011Arango et al. ( , 2012. These authors focused on describing textures and structures and the petrography, geochemistry, and geochronology of 2.3 Total-rock isotope analysis Sm, Nd, and Sr isotope analyses were performed at the Geochronology and Isotope Geochemistry Laboratory of the University of Brasilia. The analytical procedures applied in this study to determine the 147 Sm/ 144 Nd and 143 Nd/ 144 Nd isotope ratios were those described by Gioia and Pimentel (2000). 149 Sm and 150 Nd spike solutions were added to the crushed and pulverized rock. Sm and Nd were separated on cation exchange columns. Two drops of 0.025 N H 3 PO 4 were added to the resulting fractions, which were then evaporated. The residue was dissolved in 1 µl of 5 % HNO 3 and mounted on a double Re filament on a Finnigan MAT 262 thermal ionization mass spectrometer with seven collectors in static mode. The uncertainties of the 147 Sm/ 144 Nd and 143 Nd/ 144 Nd isotope ratios were lower than 0.2 % and 0.0033 % (2σ), respectively, based on the analysis of the international standard BHVO-2. The 143 Nd/ 144 Nd ratio was normalized using a 146 Nd/ 144 Nd ratio = 0.7219 and a decay constant of 6.54 × 10 -12 a -1 (Lugmair and Marti, 1978).
The procedures followed for Sr isotope analysis were those presented by Gioia et al. (1999). The samples were pulverized and subjected to acid dissolution and separation on cation exchange columns. Subsequently, the Sr-containing fractions were deposited together with 1 µl of H 3 PO 4 on a Ta filament in the mass spectrometer described above. Based on the analysis of the international standard NBS-987, the uncertainties of the 87 Sr/ 86 Sr ratio were lower than 0.0069 % (2σ). The 87 Rb/ 86 Sr ratio was calculated based on the Rb and Sr concentrations of the samples assessed by ICP-MS (Annex 1), according to the procedure of Faure and Mensing (2005).

LA-ICP-MS U-Pb zircon geochronology
The samples were initially prepared at the EAFIT University (originally Escuela de Administración, Finanzas e Instituto Tecnológico), where zircon crystals were separated following standard crushing, sieving, and mineral concentration methods. Subsequently, at the University of Rochester, several crystals of each sample were mounted in epoxy resin and polished to expose an internal face, on which the textural and isotopic analyses were performed. Cathodoluminescence images, taken to reveal the internal structure of the study zircons and to guide the geochronological analyses, were acquired under a JEOL 7100FT field-emission scanning electron microscope with a Deben panchromatic cathodoluminescence detector at the Mackay Microbeam Laboratory, University of Nevada, Reno.
The U-Pb isotope analyses were performed at the University of Rochester on an Agilent 7900 ICP-MS coupled to a Photon Machine Analyte G2 laser ablation system. This contained a HelEx2 rapid purge cell and generated pulses lasting approximately 8 ns using an excited ArF excimer. During the analyses, the laser was operated at a repetition rate of 7 Hz and a spot size of 30-µm diameter, generating a constant energy density of approximately 7 J/cm 2 on the surface of the crystals. The ablations were conducted under an ultrahigh-purity He atmosphere, and the same gas was used to transport the ablation aerosol to the ICP-MS. Each analysis consisted of 15 s of an "analytical blank" measurement with the laser off, immediately followed by a 20-s measurement of the isotope composition of the crystals with the laser on. The isotopes 202 Hg, 204 (Pb+Hg), 206 Pb, 207 Pb, 208 Pb, 232 Th, and 238 U were measured on the mass spectrometer. The instrumental bias was corrected using the standard-sample bracketing method with fragments of a Sri Lankan zircon (SL2) with a known age of 563.6 ± 3.2 Ma, assessed by isotope dilution thermal ionization mass spectrometry (ID-TIMS), as the primary reference material (Gehrels et al., 2008).   U-Pb data were processed using the procedures and algorithms described by Pullen et al. (2018) to make the necessary linearity corrections of the detection system. In addition to the primary reference material and to zircons from the San José de Guaviare Syenite, during our analytical session, several Plešovice and FC-1 zircon fragments (Duluth Gabbro Complex, Minnesota) were also analyzed to assess the precision, accuracy, and reproducibility of the method. Said secondary materials yielded weighted average ages of 206 Pb/ 238 U = 335.0 ± 1.2/3.9 (2σ, n = 21, MSWD = 0.8) for Plešovice and 207 Pb/ 206 Pb = 1098 ± 8/11 (2σ, n = 23, MSWD = 0.3) for FC-1 zircons. The first level of uncertainty reported represents only internal analytical uncertainties, and the second level of uncertainty includes the propagation of systematic uncertainty and standardization sources. By including the systematic uncertainties, the values assessed in our analytical session become indistinguishable from the reference ages of 337.13 ± 0.37 Ma for Plešovice (Sláma et al., 2008) and 1099.96 ± 0.58 Ma for FC-1 (Ibáñez-Mejía and Tissot, 2019), determined by CA-ID-TIMS, which confirms that the geochronological results reported in this study are accurate and that the analytical uncertainties have been correctly assigned.

gEOLOgIcAL fRAMEwORk
The SJGS represents a period of tectonism-related intrusions after the Precambrian bodies of the Amazonian Craton basement (Pinson et al., 1962;Vesga and Castillo, 1972), and this magmatic event was the last relevant igneous activity that occurred in the western section of the Amazonian craton (Galvis et al., 1979).
According to the K-Ar and Rb-Sr biotite ages assessed by Pinson et al. (1962), these rocks would be Paleozoic, dated between the Early Ordovician and the Late Ordovician, with different episodes of intrusion, although the K-Ar and Rb-Sr biotite dating methods used actually indicate a cooling below approximately 300 °C. Toussaint (1993) indicates that the nepheline syenite bodies correspond to the last early Paleozoic magmatism of the Autochthonous Block (Bloque Autóctono), resulting from local partial melting of the base of the continental crust of the Eastern Range (Llanos Orientales), which would have been affected by a stretching movement that increased the thermal flow in the region.
Nevertheless, the U-Pb zircon age assessed by Arango et al. (2012) dated the unit to the Ediacaran (Neoproterozoic), with a proposed origin from anatexis processes of migmatitic metamorphic rocks, which are related to the Pan-African event during the Cambrian in the continent Gondwana. The U-Pb ages assessed in the study conducted by Franco et al. (2018) are similar, but the authors diverge in their interpretation of the origin of these rocks, which indicates that the intraplate alkaline magmatism of the Neoproterozoic-Cambrian, east of Colombia and absent from the Colombian Andes, could have derived from Pan-African orogeny during the accretion of mobile belts in the NW continental margin of Gondwana.
The SJGS rocks are overlaid by sedimentary formations deposited in unconformity over the unit, such as the San José Formation from the Cretaceous (Figure 1) (Arango et al., , 2012 and the Caja Formation from the Neogene (Maya et al., 2018). Furthermore, in the southern section of the Cerritos body, the SJGS is in contact with Proterozoic metamorphic rocks (Maya et al., 2018), where igneous rocks intrude into metamorphic rocks, creating a thermal effect, evidenced in the recrystallization to hornfels and in the formation of injection migmatites. Approximately 6 km west of the El Capricho body, a small syenite outcrop ( Figure 2) shows that the extent of the unit is mostly under the sedimentary cover.
Approximately 15 km southwest of the syenite bodies, basement rocks, named the Guaviare Complex by Maya et al.

Lithology
The San José de Guaviare Syenite, in the southern section, encompasses two bodies aligned NW-SE located in the El Capricho and Cerritos sectors ( Figure 2). In general, these bodies correspond to holocrystalline, phaneritic, hypidiomorphic rocks, ranging from equigranular ( Figure 3A) to inequigranular. They are leucocratic, with a medium to coarse grain size, up to pegmatitic (7 cm) ( Figure 3B), with a low to medium color index (7-20 %) and with white, gray, and pink hues mottled with black. The El Capricho body is the most homogeneous in both composition and texture. The Cerritos body tends to be heterogeneous, with locally dark-gray fine-grained facies, which preserve the mineralogical composition, as well as rocks oriented by a mafic mineral flow (Figure 3 A, C).
In the southern section of the Cerritos body, SJGS intrusion generated a contact aureole in the regional metamorphic rocks of the Guaviare Complex, forming hornfels covering an area of 0.1 km 2 . Although this type of rock only crops out in this zone, it provides considerable information about the origin of the syenites. The hornfels have a fine to medium grain, with black, gray, and red hues, some of which show banding and vestiges of the original orientation of the rock. At the same place, the rocks show silicification and form injection migmatites (Figure 3D), with stromatic textures. They are also characterized by the presence of syenite dikes, intersecting the injection migmatites and enveloping hornfels xenoliths.
The rocks with quartz correspond to syenogranite, syenite, quartz-syenite, and quartz-alkali feldspar syenite ( Figure 4). According to the quartz content, two facies are differentiated (   ). The zone of contact between the SJGS and metamorphic rocks of the Guaviare Complex shows silicification (Qz > 68 %), most likely through host rock assimilation; this sample was classified as a quartz-rich granitoid ( Figure 4). In general, these rocks have a holocrystalline, hypidiomorphic, inequigranular to equigranular texture, with a medium to coarse grain. Some samples show mineral orientation, defined by biotite, and others show recrystallization with a polygonal granoblastic texture ( Figure 5c). Mineralogically, they consist of microcline, plagioclase, nepheline or quartz, biotite, arfvedsonite, and aegirine.
Microcline presents as subhedral to anhedral, tabular crystals with gridiron and Carlsbad twinning, which implies formation from orthoclase. Sometimes the microcline forms an exsolution (perthite and mesoperthite) and myrmekitic fan texture, with plagioclase, nepheline, biotite, and calcite inclusions, and is slightly altered to clays and sericite. Albite to oligoclase plagioclase, determined using the Michel-Levy method, forms individual subhedral to euhedral, tabular, zoned crystals with albite, albite-Carlsbad, and Carlsbad twinning and with mesoperthitic and anti-perthitic textures, along with wedge-shaped lamellae and with calcite and biotite inclusions, sometimes altered to sericite and recrystallized along the grain boundaries, which causes the loss of extinction that is characteristic of this mineral.
Nepheline is found in subhedral to euhedral crystals with shapes ranging from rectangular to hexagonal, with simple twinned crystals and microcline, plagioclase, biotite, and calcite inclusions in poikilitic textures, showing slight alteration to   (Streckeisen, 1976) sericite, cancrinite, and natrolite. Quartz is found in colorless, clean, uniformly distributed anhedral crystals with undulose extinction, forming granular mosaics of complete contact with microcline and plagioclase.
Biotite is found in subhedral to euhedral, tabular crystals with strong pleochroism with X (α): pale yellowish brown, Y (β) = Z (γ): dark reddish or greenish brown, simple twinning, and bird's eye and undulose extinction. Sometimes the biotite shows a preferential orientation in undeformed clustered laths, which is therefore interpreted as the flow orientation. Some crystals have pleochroic halos around inclusions of zircon, while others are embayed, skeletal crystals, and they are replaced by muscovite, sericite, and chlorite along the grain boundaries.
Arfvedsonite amphibole forms subhedral to anhedral, columnar to hexagonal crystals, with simple twinning and strong pleochroism with X (α): bluish green, Y (β): greenish brown, and Z (γ): olive green, in columnar aggregates with biotite, titanite, and zircon. It is replaced by biotite towards the edges and fractures. Pyroxene-like aegirine is found in short, columnar, subhedral to anhedral crystals with strong pleochroism with X (α): dark green, Z (γ): brownish yellow, and Y (β): green. Aegirine is associated with arfvedsonite through possible transformation to amphibole. Cleavage planes are filled by opaque minerals possibly due to oxidation.
The primary accessory minerals correspond to granular opaque minerals, with sharp rhombic, square shapes or with skeletal textures, which correspond to pyrite and magnetite/ilmenite, in addition to titanite with simple twinning, some fractured and overgrown zircons, calcite-like carbonate, epidote, apatite, tourmaline, allanite, pyrochlore, and garnet, possibly melanite. Franco et al. (2018) mention that the syenites contain large magmatic zircons with a partly metamictic character.
Secondary accessory minerals, such as cancrinite, sericite, and natrolite, result from the alteration of nepheline in the cleavage planes, edges, and fractures of the host mineral. Sericite aggregates with carbonates are also an alteration product from the interior and edges of microcline, plagioclase, and biotite crystals.

Hornfels
The different types of hornfels identified by petrographic analysis highlight how the compositional variation of the host rock was affected by thermal metamorphism during the SJGS intrusion. Thus, pelitic and mafic protolith hornfels were initially subjected to regional metamorphism within the Guaviare Complex and later to thermal metamorphism in amphibolite facies (previously hornblende hornfels) and pyroxene hornfels. In addition, feldspar-quartz hornfels correspond to a syenite  Arango et al. (2011Arango et al. ( , 2012 material recrystallized through intrusive processes, which occurred at various stages during the period of magmatism that originated the SJGS, as also shown by the wide range of ages assessed by geochronology and the many multielement and REE patterns found in geochemistry. These rocks have a polygonal granoblastic texture and occasionally relict orientation derived from the original rock. The feldspar-quartz hornfels mineralogically consist of microcline, plagioclase, quartz, biotite, muscovite, and epidote (Annex 2), with zircon, calcite, apatite, titanite, fluorite, and opaque minerals as accessory minerals. The pelitic hornfels have a similar composition to that of feldspar-quartz hornfels but contain andalusite and sillimanite. The mafic hornfels (Figure 5d) show banding with mafic and quartz feldspar domains, compositionally differing in the occurrence of hornblende and clinopyroxene.
Quartz is detected as xenoblast crystals, with undulose extinction and granoblastic to polygonal texture. Microcline is identified as subidioblastic to xenoblastic, tabular crystals with gridiron twinning and slight sericitization. Albite and oligoclase plagioclase, as determined using the Michel-Levy method, ranges from subidioblastic to idioblastic and is tabular,  with albite twinning, apatite inclusions, and slight alteration to saussurite. Biotite ranges from subidioblastic to idioblastic and is laminar with strong pleochroism, X (α): very pale yellow to Y (β) = Z (γ): dark brown, some with mimetic orientation, whereas others are disordered. Occasionally, they exhibit irregular boundaries and moderate chloritization. Muscovite is identified as xenoblastic to subidioblastic, laminar crystals with irregular boundaries and without any orientation, in intergrowths of biotite. Hornblende is subidioblastic, columnar, with strong pleochroism X (α): pale yellow, Z (γ): dark green, Y (β): greenish yellow, forming aggregates with titanite. Clinopyroxene (diopside) is found associated with hornblende in xenoblasts with slight pleochroism X (α): pale bluish green, Z (γ): pale brownish green, Y (β): pale brownish yellow. Hornblende appears to replace clinopyroxene, most likely due to retrogression. Epidote is associated with biotite, suggesting epidotization, whereas zircon is found in inclusions, forming pleochroic halos in biotite. The opaque minerals are identified as pyrite, ilmenite, and magnetite.

Geochemistry
Eighteen representative samples of the petrographic groups were selected to analyze major, minor, and trace elements (Annex 1). The rocks from the San José de Guaviare Syenite are classified as nepheline syenites, syenites, and granites, according to the classification diagram of plutonic rocks by Cox et al. (1979) (Figure 6). They were subdivided into three groups according to their major and trace elements. In general, the rocks are peralkaline, according to the alkalinity index by Maniar and Piccoli (1989), but some rocks are placed in the metaluminous (group 3, mainly) or peraluminous (group 2) fields ( Figure 7).
The samples from group 2 are mostly classified as nepheline syenites, one of which is classified as syenite ( Figure 6). They show a higher SiO 2 content than those of group 1 (56.35-63.86 %) and are characterized by their low CaO (0.04-0.49 %), MgO (0.03-0.46 %), and P 2 O 5 (< 0.01-0.09 %) and high Al 2 O 3 (18.53-22.76 %) and Na 2 O contents (6.10-9.56 %). Some samples have abnormally high concentrations of Zr (4070 and 1970 ppm). These values may be related to the abundance of zircon. In this group, the high Nb content is also notable (up to 625 ppm in sample JDB-0011R), which can be explained by the presence of trace pyrochlore contents in the rock. The rocks of groups 1 and 2 are generally impoverished in TiO 2 , and the samples from group 2 have a very low modal proportion or do not contain titanite in their mineralogy, in contrast to those from groups 1 and 3.  The samples associated with group 3 include rocks saturated and subsaturated with SiO 2 (56.31-74.20 %). They are mainly classified as syenites and granites ( Figure 6). In contrast to the rocks from the two previous groups, these rocks are impoverished in Al 2 O 3 (12.63-18.53 %) and Na 2 O (3.30-6.63 %) and enriched in CaO (0.26-3.60 %), MgO (0.13-1.82 %), and P 2 O 5 (0.02-0.53 %).
Group 2 is heterogeneous and shows a pointed pattern, with negative Ba anomalies and positive U, K, Sr, and Zr anomalies in some samples. Group 3 shows a relatively coherent pattern, with negative Sr and Ti anomalies and with some variability in elements such as Ba, Th, U, and LREE.
In chondrite-normalized REE patterns (Figure 8), group 1 shows enrichment in LREE over middle (MREE) and heavy rare-earth elements (HREE), and the slope of the curve is relatively gentle (La/Yb) N = 19.2 -23.2. Group 2 displays a "reverse spoon" pattern, showing a relative impoverishment in MREE [(La/Sm) N = 4.0 -11.1 (Sm/Yb) N = 0.2 -2.9] compared to LREE and HREE; some samples show a slight (CAL-0023R and CAL-0017R) to strong positive Eu anomaly (ENA-0036R), which   indicates a possible plagioclase accumulation during the magmatic history or the inheritance of this pattern from the original source. Group 3 is less homogeneous than the other groups and shows an enrichment in LREE in relation to HREE, with a relatively gentle slope [(La/Yb) N = 8.9 -28.0]; these samples show a negative Eu anomaly, except two samples (CAL-0008R and CMR-0013R) had a slightly positive Eu anomaly. Most likely, the REE patterns reflect early differentiation of titanite and consequent loss of these magmatic elements.
The phase diagram of the NaAlSiO 4 -KAlSiO 4 -SiO 2 system (Hamilton and Mackenzie, 1965) (Figure 9) shows that the Si-subsaturated rocks of the San José de Guaviare Syenite are defined by a crystallization sequence in which, once the temperature decreases, two pathways occur: one in which nepheline crystallizes first, and another in which potassium feldspar crystallizes first, as shown by the textural relationships of euhedral nepheline inclusions in feldspar crystals, whereas other samples exhibit feldspar crystals in nepheline. In both cases, the residual liquid corresponds to nepheline syenite magma.
Si-saturated magmas (especially in group 3), which are found at the top of the Ab-Or line, when the temperature decreases, initially crystallize alkali feldspar or quartz. As the temperature decreases, the liquid increases or decreases in Si. The final outcome of the fractionation corresponds to granite magma. These paths are reflected in granitic rocks in which quartz is the last crystallization phase.
Rocks with intermediate compositions between both compositions, which correspond to syenites, should tend towards either of the previous two during the fractional crystallization process, meaning the abundance of these rocks results from magma mixing (Storey et al., 1989), assimilation of continental crust (Motoki et al., 2010), or absorption of previously precipitated cumulus feldspar material by new magma inside the magma chamber (Wolff, 2017).
The Silica Saturation Index (SSI) vs. (K + Na)/Al diagram (Figure 10a) shows that the rocks from groups 1 and 2 are classified as alkaline rocks, whereas the rocks from group 3 primarily belong to the nonalkaline field. In general, the samples define a linear pattern, especially those of group 3, in which (K + Na)/Al does not decrease when transitioning from a SiO 2 -subsaturated (SSI < 0) to a SiO 2 -saturated (SSI > 0) magma, which indicates crust assimilation, possibly alkaline. The patterns observed in the Rb/K vs. Nb/Y diagram (Figure 10b) suggest that the samples from group 1 are the least evolved and the closest to the original composition of the magma and that the samples from groups 2 and 3 were mainly formed by fractional crystallization and assimilation of continental crust, respectively.
The fractioning of titanite and zircon is noticeable in the diagrams of Figure 10c, d, e, in which the samples of the three groups are clearly separated. The anomalies observed in the REE patterns are also expressed in the diagrams of Figure 10f,   Figure 9. a) Rock composition in the Q-Ne-Ks system; b) Equilibrium phase diagram for the NaAlSiO 4 -KAlSiO 4 -SiO 2 system at 1 kbar Source: adapted from Hamilton and Mackenzie (1965).   g, h, which show the effect of assimilation processes of the continental crust in addition to fractional crystallization processes. In the latter, the samples of groups 1 and 2 are connected by trends of fractional crystallization, albeit opposite to those observed in Figure 10b. Field relationships showed the presence of metamorphic xenoliths partly absorbed in syenites in the southern section of the body. These xenoliths comprise rocks rich in SiO 2 , which include the Si-saturated rocks of group 3. However, the relationship between the rocks of groups 2 and 3 is not limited to a simple assimilation model. The patterns observed in the diagrams presented by Motoki et al. (2010Motoki et al. ( , 2015 indicate combined processes of fractional crystallization in groups 1 and 2 and assimilation in group 3.
In the Rb vs. (Y+Nb) tectonic discrimination diagram (Pearce et al., 1984), all rocks are placed in the within-plate granite (WPG) field (Figure 11), which corroborates the findings of other studies on the same unit, conducted north of the study area (Arango et al., , 2012. The comparison analysis of syenite rocks identified to the south and those to the north of SJGS reported in the study by Arango et al. (2011Arango et al. ( , 2012 shows that the compositional variation is higher in the southern SJGS rocks, which include nepheline syenites, syenites, and granites, whereas the northern SJGS rocks consist of only nepheline syenites ( Figure 6).
In the multielement and REE diagrams (Figure 12), the northern SJGS (Arango et al., , 2012 shows geochemical patterns similar to those of groups 2 and 3, defined for southern SJGS rocks; three samples show geochemical patterns similar to those of group 2, and the other three, to those of group 3. Sample 5000421 is deviates slightly from the pattern because it is impoverished in HREE relative to the group 3 field.
Although the composition of northern SJGS rocks is restricted to nepheline syenites, the main geochemical characteristics are similar to those of southern SJGS rocks, as are the concave-upward REE patterns, resulting from the impoverishment in MREE, slightly positive Eu anomaly, and positive Sr and Zr anomalies in group 2 and the relative enrichment in LREE in relation to the other groups, slightly negative Eu anomaly, and negative Sr and Zr anomalies in group 3.

Isotopes
The results from the Sm-Nd and Sr isotope analyses of two nepheline syenite (CAL-0014R, group 1, and JDB-0011R, group 2) and one syenogranite sample (JDB-0024R, Grupo 3) are outlined in Table 1. Based on geochronological data (Annex 3), a crystallization age of 604 Ma was used to calculate the ε Nd (T) and the initial 87 Sr/ 86 Sr ratio of the samples.
The samples have positive values of ε Nd (604), ranging from +0.8 to +1.8, with T DM model ages ranging from 920 to 1000 Ma (Figure 13a). The ε Nd values of the two nepheline syenite samples are very similar (+1.8 and +1.6) and slightly higher than that of the syenogranite sample (+0.6). In turn, the T DM model ages of the nepheline syenite samples are consistent Figure 11. Tectonic discrimination diagram (Pearce et al., 1984) of the SJGS The samples belong to the field of within-plate granites (WPG)  Wasserburg, 1980, 1984). T DM model ages according to the depleted mantle model of DePaolo (1981). The concentrations of Rb and Sr were determined by ICP-MS analysis (Annex 1). The 87 Rb/ 86 Sr ratio was calculated according to the procedure described by Faure and Mensing (2005).  (Arango et al., , 2012 and southern SJGS rocks The red and blue fields represent groups 2 and 3 defined in this study, respectively. The samples from the northern SJGS bodies are represented by lines.   (0.92 and 0.93 Ga), whereas the model age of the syenogranite sample is slightly older (1.0 Ga). With initial 87 Sr/ 86 Sr ratios ranging from 0.7048 to 0.7067, the isotope signatures of these samples range from those of typical mantle-derived magmas within the mantle array region until reaching decoupled signatures to the right of this field, characterized by high initial ratios of 87 Sr/ 86 Sr (Figure 13b). The latter is the case for the syenogranite sample, which has an initial 87 Sr/ 86 Sr ratio of 0.7067, whereas the nepheline syenite samples have lower ratios of 0.7048 and 0.7054.

Geochronology
Geochronological analyses were performed using the LA ICP-MS U-Pb zircon method in two syenite samples (Annex 3), one located in the El Capricho body (JDB-0001R) and the other located in the Cerrito body (OPP-0074R) (Figure 14). Authors such as Pinson et al. (1962), Arango et al. (2011Arango et al. ( , 2012, Franco et al. (2018), and Muñoz Rocha et al. (2019) have reported ages ranging from the Ordovician to the Neoproterozoic for syenite rocks, calculated using different methods (Table 2). By K-Ar and Rb-Sr isotope dating of biotites,   The ages calculated by Pinson et al. (1962) are corrected according to the constants of Steiger and Jager (1977) and to data published by Maya (1992). Pinson et al. (1962) assessed ages ranging from 445 to 495 Ma for the La Pizarra body, which were considered ages of cooling by uplift followed by erosion or metamorphism, processes placed between the early and the Late Ordovician. A similar age of 494.1 ± 5 Ma, assessed in the same La Pizarra body using the 40 Ar/ 39 Ar biotite method, was reported by Arango et al. (2011Arango et al. ( , 2012. This age corresponds to the late Cambrian (Furongian) and was interpreted as a cooling age. Another age reported by Arango et al. (2011Arango et al. ( , 2012, using the U-Pb zircon method, is 577.8 ± 6.3-9 Ma (Ediacaran), which would be considered the age of crystallization. Other, more recent studies have yielded crystallization ages, assessed using the U-Pb zircon method, of 602 ± 3 and 608 ± 1 Ma (Ediacaran) for El Jordán syenite (Franco et al., 2018;Muñoz Rocha et al., 2019), which is found the southern section of the Cerritos body ( Figure 14).

Nepheline syenite (JDB-0001R)
Zircons are in general anhedral and irregular, with a coarse concentric zoning of dark centers and lighter edges (Figure 15a, Annex 4). In dark zones, the uranium content is high, and the igneous fine oscillatory zoning is well defined. In contrast, in light zones, the uranium content is low, and the zoning is less perceptible or completely absent. The origin of this coarse zoning is not defined and could be associated with patchy zoning, attributed to local recrystallization along microfractures, according to Corfu et al. (2003). This characteristic is observed in one of the crystals with light veinlets that penetrate the dark zone   the age of the magmatic event ( Figure 16a) and corresponds to the Ediacaran period of the Neoproterozoic. The 206 Pb/ 238 U age that should be cited is 604 ± 7 Ma because this value was calculated considering all sources of uncertainty, both internal and systematic, as well as standardization. Two (inherited) zircon xenocrysts yielded individual 206 Pb/ 238 U ages of 647 ± 14 and 662 ± 19 Ma (Figure 16b, blue ellipses), and their ages support the assumption of a lengthy magmatism because the area shows no known igneous or late Neoproterozoic metamorphic events.  with good igneous fine zoning. The dark, uranium-rich areas may have begun to grow in a magma with a high content of this element, but the amount of uranium decreased due to assimilation of host rocks, as shown by the geochemical analysis. The weighted-average 206 Pb/ 238 U age of this rock is 604.3 ± 2.2/7.0 Ma (2σ, n = 45, MSWD = 1.6), where the first level of uncertainty represents only internal analytical uncertainties, and the second level of uncertainty includes propagation from systematic and standard sources. This dating is interpreted as

Syenite (OPP-0074R)
The zircons in this sample are euhedral and prismatic, with well-defined igneous oscillatory zoning along the grain boundaries, whereas cores with concentric, sectorial, and chaotic zoning, homogeneous light and dark textures, and inclusions are observed at the center in some crystals (Figure 15b, Annex 5). Other zircons have radial fractures, which are attributed to marked differences in uranium content within the zircon, according to Corfu et al. (2003). The weighted-average 206 Pb/ 238 U ages of this sample is 620.5 ± 3.6/7.5 Ma (2σ, n = 14, MSWD = 1.2), which is interpreted as the crystallization age of this intrusive body during the Ediacaran (Figure 17a). For sample JDB-0001R, the 206 Pb/ 238 U age that should be cited is 620.5 ± 7.5 Ma, which includes all sources of uncertainty. Besides the zircons from the magmatic event, two inherited zircons yielded older 206 Pb/ 238 U ages of 648 ± 10 and 685 ± 15 Ma ( Figure 17b). Four analyses with apparently younger 206 Pb/ 238 U ages are clearly discordant (gray ellipses in Figure 17b); these values are affected by Pb loss and were therefore disregarded when calculating the crystallization age of this sample.

DIScuSSION
The petrological evidence, the geochemical data, and the isotopic compositions of the nepheline syenite samples analyzed in this study (Figure 13b) indicate that the SJGS rocks derive from anorogenic and mantle-derived alkaline magmas subjected to magmatic differentiation and crustal assimilation.
Some of the magmas underwent fractional crystallization, resulting in Si-subsaturated compositions, which formed nepheline syenites. Others possibly reacted at the crustal level with the host rock, metamorphic rocks of the Guaviare Complex of varied compositions (Maya et al., 2018), producing some silica-saturated magmas, such as syenogranite. Such a contamination process is reflected in the isotope composition of a syenogranite sample, with an anomalously high initial 87 Sr/ 86 Sr ratio (0.7067, in contrast to the ratios of 0.7048 and 0.7054 of the nepheline syenite samples), which could have resulted from crustal assimilation of host rocks and/or rocks with low Rb/Sr ratios. Simultaneously, this contamination would also produce a lower ε Nd value, assessed in syenogranite (+0.6), and an apparently older T DM model age of 1000 Ma, in contrast to the model ages of 920 and 930 Ma assessed in the nepheline syenite samples. A similar origin has been suggested in other regions of the world where rocks are simultaneously subsaturated and saturated in silica (Foland et al., 1993;Zhu et al., 2016).
U-Pb zircon ages reported by different authors, such as Arango et al. (2012), Franco et al. (2018), and Muñoz Rocha et al. (2019), and those assessed in this study indicate that syenite was intruded in the Ediacaran period of the Neoproterozoic, in a magmatic age range of at least 30 Ma, between 621 − 7 and 577.8 + 6.3 Ma, which would indicate that the syenite magmas were not all intruded simultaneously but rather that multiple intrusions occurred during that period. The samples dated in this study contain inherited zircons with a minimum age of 662.3-19.2 Ma, that is, 643.1 Ma (Annex 3), and these zircons may be "antecrysts" formed in initial magma chambers of the syenites because no rocks of those ages are known in the area. Thus, the alkaline magmatism would have lasted approximately 59 million years. In the state of Bahia, Brazil, the alkaline magmatism associated with rifting lasted at least 58 million years (Rosa et al., 2007).
The Ar-Ar, K-Ar, and Rb-Sr biotite ages of intrusive rocks generally correspond to cooling ages. However, the large interval between the most recent U-Pb age (577.8 ± 6.3 Ma) and the oldest K-Ar age (508 ± 26 Ma) ( Table 2) indicates a cooling age, below 300 °C, of 37 million years. Considering that hornfels were formed in the contact between syenite and regional metamorphic rocks, with injection migmatites, and not regional anatectic migmatites, the intrusion would have occurred at a depth of a few kilometers. Accordingly, cooling to 300 °C unlikely lasted 37 million years. Moreover, data presented by Arango et al. (2012) show that stepwise-extracted Ar fractions do not plateau well, which may result from excessive recoil during the 39 K(n, p) 39 Ar reaction produced by irradiation, a complex thermal history, or the combination of both factors. The two possible events that may have generated a complex thermal history and that may have disturbed the K-Ar/ 40 Ar-39 Ar systems are: 1) the Cambro-Ordovician sedimentary cover in the area may have been quite thick, which would have resulted in a complex cooling history of the syenites or allowed the syenites to reheat during tectonic burial; and 2) an important thermal event may have affected the western margin of the Amazonian Craton during the Cretaceous, as suggested by the presence of basaltic dikes dated at 102 Ma in the Araracuara region (Ibáñez-Mejía and Cordani, 2020). Although these rocks have not been reported in the Guaviare region, currently available data prevent us from ruling out the presence of a recent thermal event (for example, Cretaceous) in the area. Thus, the thermal and post-Ediacaran tectonic burial history of the zone should be further studied to address these unanswered questions in the future.
The following model is proposed for the zone: At the end of the Neoproterozoic, a series of rifts enabled the disintegration of the Rodinia supercontinent into new continents. During this rifting period, at approximately 600 Ma, the northern zone of the Amazon was separated from Laurentia and Baltica, forming the Iapetus ocean (Condie, 2011;Cawood and Pisarevsky, 2017). Some of the rifts in the NW Amazon could have favored the intrusion of alkaline silica-subsaturated magmas, most likely generated by partial melting of mantle-derived igneous rocks during previous magmatic events, as indicated by the T DM model ages of 920 and 930 Ma assessed in the nepheline syenite samples. These ages suggest that the source rocks of the SJGS magmas were likely formed during the Putumayo tectono-magmatic event, whose thermal and lithological record is common in several crustal segments of the northern Andes (Kroonenberg, 1982;Restrepo-Pace et al., 1997;Ramos, 2010;Cardona et al., 2010;Ibáñez-Mejía et al., 2011. Based on the isotope characteristics described in this study, an alternative origin of nepheline syenites related to juvenile mantle magmas (separated from the mantle at ~600 Ma) contaminated with older crust can also be ruled out; if such an origin were true, this contamination should have produced Si-saturated magmas, and therefore nepheline syenites would not have been formed in the first place, and their isotope signature would also be altered, moving away from the mantle array, as in the case of syenogranite.
Other authors have proposed different origins, such as crustal anatexis (Arango et al., , 2012Toussaint, 1993) or intraplate magmatism associated with Pan-African orogeny (Franco et al., 2018), which are not supported by isotope or regional geological data. Furthermore, Si-subsaturated magmas most likely cannot be derived from partial melting of common crustal rocks, which are typically saturated in quartz. Most data on alkaline magmatism reported in the literature indicate mantle sources with isotope signatures similar to those found in this study (Fitton and Upton, 1987;Winter, 2001).

cONcLuSIONS
U/Pb zircon dating of SJGS rocks places the unit in the Ediacaran, between 604 ± 7 and 620.5 ± 7.5 Ma. The SJGS derives from anorogenic magmatism in an intraplate setting, most likely resulting from rift-like stretching during the fragmenta-tion of Rodinia at the end of the Neoproterozoic. The sources of the SJGS magmas would have been magmatic rocks directly formed by partial mantle melting during the Putumayo tectono-magmatic event.
The different syenite bodies consist of nepheline syenites, nepheline-bearing alkali-feldspar syenites, syenites, and syenogranites. The coexistence of nepheline-and quartz-bearing rocks could have resulted from the reaction between syenite magma and crustal rocks, which would have produced some silica-saturated magmas, such as the syenogranites that occur in southern SJGS bodies (El Capricho and Cerritos).
The southern SJGS syenite bodies tend to have broader compositions than the northern SJGS bodies, where the recorded rocks exclusively comprise nepheline syenites.
The SJGS is subdivided into three geochemical groups with contrasting characteristics that denote processes associated with the magmatic evolution of the igneous body, which include crustal differentiation and assimilation.
The intrusion of the SJGS generated a thermal metamorphism, in amphibolite and pyroxene-hornfels facies, in the host rocks of the Guaviare Complex, which allowed the recrystallization of minerals and the formation of hornfels and injection migmatites. A history of slow cooling, or a more recent thermal event not yet identified, affected the K-Ar/ 40 Ar-39 Ar systems of the syenites, yielding significantly younger ages than those assessed using the U-Pb zircon dating method.

AckNOwLEDgEMENTS
This research summarizes and interprets part of the results from the project Elaboración de la cartografía geológica en un conjunto de planchas a escala 1: 100 000 ubicadas en dos zonas del territorio colombiano, zona sur [Geologic mapping of a set of plates at a 1:100 000 scale located in two zones of southern Colombia] of the Dirección de Geociencias Básicas del Servicio Geológico Colombiano (SGC), conducted by the geology team of the Serviminas company under contract No. 508 of 2017. Finally, the authors thank the anonymous reviewers for their valuable comments and suggestions, which helped to improve this manuscript.