Petrographical and Mineralogical Characteristics of Magmatic Rocks in the Northwestern Siberian Traps Province, Kulyumber River Valley. Part I: Rocks of the Khalil and Kaya Sites

The origin of the Siberian Traps province has been under discussion for the last three decades. Up to now, there is no real model of its formation in a good agreement with geological data on the magmatic evolution at P–T boundary in Eastern Siberia. Modern geochemical data on magmatic rocks around the province is a key to reconstructing magmatic development in time and space. Such data have been obtained for the Norilsk and Meimecha–Kotuy and not for other parts of the Siberian province. For the first time, we studied the geochemistry and mineralogy of magmatic rocks at the Kulyumber river valley, located in the intersection of the Tunguska syneclise and Norilsk–Igarka zone in the NW Siberian platform. In this article, we present data from the Khalil and Kaya sites of this area belonging to the Syverminsky, Gudchikhinsky, Khakanchansky and Nadezhdinsky formations. Their mineralogical and geochemical features (including Sr, Nd and Pb isotope data) are similar to the same formations in the Norilsk area, while the rocks belonging to the Gudchikhinsky formation show differences. The Syverminsky tuffs are also described for the first time. The intrusive rocks are attributed to four intrusive complexes, i.e., Ergalakhsky, Kureysky, Katangsky and Norilsk. The Ergalakhsky complex comprises trachydolerites similar to the trachydolerites of the Norilsk area. The rocks of the Norilsk complex at the Khalil site differ from the rocks of the same complex at the Norilsk area by the (U/Nb)n = 1.8, (La/Yb)n = 2.1 in comparison with 3.7 and 2.3 of the rocks of the Norilsk 1 intrusion. The intrusions of the Kureysky complex are more differentiated than those of the Katangsky intrusions but show comparable TiO2 and trace elements distribution. Thus, the magmatism of the Kulyumber area is characterized by features matching those of the Tunguska syneclise and Norilsk area, i.e., suggesting rift and platform regimes.

The great scale of the province and many remote areas in Eastern Siberia are the main reasons for incomplete information on this great volume of magmatic rocks. Modern geochemical studies have so far been conducted only for the Norilsk and Meimecha-Kotuy regions, which comprise PGE-Cu-Ni magmatic deposits and high alkaline rocks, respectively [27,28]. In these areas, major, trace and isotopic data are a valuable for lavas [29][30][31][32][33][34][35][36][37]. Data on intrusive rocks are more restricted in comparison with basalts and refer mainly to the Norilsk ore-bearing intrusions [14,[38][39][40][41][42][43][44][45]. Other parts of the province have been studied only in a reconnaissance way [46,47]-or not examined at all.
This study is focused on the Kulyumber area located in the middle part of the Kulyumber river valley which represents one of the least studied areas of the province ( Figure 1). It is located between the Tunguska syncline and Norilsk-Igarka zone [48][49][50] and thus it provides information on the magmatic events within both structures at the Permian-Triassic boundary (250 ± 0.3 Ma; [51]). Within this territory effusive rocks and numerous intrusive bodies were found during the geological mapping. The volcanic rocks were subdivided into five formations according with the formations of the lower part of the tuff-lava sequence in the Norilsk area, i.e., Syverminsky, Gudchikhinsky, Khakanchansky, Tuklonsky and Nadezhdinsky formations [14]. Intrusions have been attributed to eight intrusive complexes located both in the Norilsk district and in Tunguska syneclise [14]. This division is doubtful, since it is based on the inner structure of intrusive bodies and structural and textural features of rocks.
The purpose of this work is to study the mineral composition and geochemical characteristics of igneous rocks and compare them with the rocks of the Norilsk area and Tunguska syneclise. As has been shown in [14,17,[31][32][33], both OIB (Oceanic Island Basalts) and WPB (Within-Plate Basalts) occur in the Norilsk area and were formed in a rift system and during a platform magmatism. For the Tunguska syneclise only platform magmatism is described. Magmatism of the Norilsk area occurs mainly within the Yenisey-Khatangsky trough (YKT). The Kulyumber river valley is not considered as a part of the Norilsk and YKT but it is involved in the rift zone. There are many questions that are not answered about the formation of this region. What types of igneous rocks occur in the Kulumbe River valley? Are there any variations in the compositions of rocks along the strike of the tectonic structure from the north to the south? Or are they a constant homogeneous in composition? What types of primary magmas formed the igneous rocks of the Kulyumber river valley? How does a transition from rift magmatism to a platform one happen-abruptly or gradually? To answer these questions, we must collect information first of all on mineralogy and geochemistry of igneous rocks of this area.
Furthermore, the Kulyumber area is considered as a priority region for the discovery of new deposits because the geologists of the company, "Norilskgeology" have found the disseminated and massive sulfide ore in gabbro within this area, in contrast to other regions of the Siberian platform where only poor mineralization was detected. Does the rift zone from the Norilsk area and the Kulyumber area show a consistent trend of mineralization?

Tectonic Structure
The studied area is located in NW Eastern Siberia, 150 km to the south of the Norilsk ore district. The tectonic structure of the territories is complex due to their occurrence within two different structural blocks. The eastern part belongs to the Siberian craton consisting of a crystalline basement and a Minerals 2020, 10, 409 3 of 39 platform cover. The western part has a more complex structure due to the depth of the basement and a huge volume of sediments and intrusions within it.
First of all, we have performed a new geophysical map of the crystalline basement of north-western part of the Siberian platform where Kulyumber area is located (Figure 2). The initial data were the results of 1:200,000 scale gravimetric and 1:100,000 scale aeromagnetic mapping performed at constant barometric altitudes (700 m and 2400 m), as well as a digital landscape model GTOP030. The geophysical fields were converted to a horizontal plane with a height of 2400 m by approximation by equivalent sources [52]. This allowed to reduce the influence of the relief of the earth's surface and local inhomogeneities in the upper part of the geological section. The measured values of the geomagnetic field intensity T were converted into the vertical component of the anomalous magnetic field ∆Z, which has a more distinct relationship with the anomaly forming objects. The network of points for digital models of fields and terrain used in processing is 2 × 2 km Minerals 2020, 10, x FOR PEER REVIEW 3 of 39 basement and a platform cover. The western part has a more complex structure due to the depth of the basement and a huge volume of sediments and intrusions within it. First of all, we have performed a new geophysical map of the crystalline basement of north-western part of the Siberian platform where Kulyumber area is located (Figure 2). The initial data were the results of 1:200,000 scale gravimetric and 1: 100,000 scale aeromagnetic mapping performed at constant barometric altitudes (700 m and 2400 m), as well as a digital landscape model GTOP030. The geophysical fields were converted to a horizontal plane with a height of 2400 m by approximation by equivalent sources [52]. This allowed to reduce the influence of the relief of the earth's surface and local inhomogeneities in the upper part of the geological section. The measured values of the geomagnetic field intensity T were converted into the vertical component of the anomalous magnetic field ΔZ, which has a more distinct relationship with the anomaly forming objects. The network of points for digital models of fields and terrain used in processing is 2 × 2 km.  River, about 100 km wide, bordering the ancient craton and extending beyond the southern frame of the area. This zone is probably a fragment of an ancient rift structure with increased mobility and permeability for magma, saturated with basic-ultrabasic magmatic rocks. The gravity anomaly has the form of a flattened funnel, elongated in the meridional (N-S) direction, with a "vent" (or "main supply channel") going into the upper mantle, which is located beneath the in area of Igarka city. The zone is characterized by low values of the magnetic field due to the presence of iron only in silicate form.
Of particular interest are the sub-latitudinal zone of low magnetic field values and slightly elevated values of the gravity field, covering the basins of the Kulyumber river and stretching to the upper flows of the Kureyka river (to the east from Snezhnogorsk city, Figure 2). According to  Figure 2 shows a submeridional band of positive values of the gravitational field stretched along the Yenisey River, about 100 km wide, bordering the ancient craton and extending beyond the southern frame of the area. This zone is probably a fragment of an ancient rift structure with increased mobility and permeability for magma, saturated with basic-ultrabasic magmatic rocks. The gravity anomaly has the form of a flattened funnel, elongated in the meridional (N-S) direction, with a "vent" (or "main supply channel") going into the upper mantle, which is located beneath the in area of Igarka city. The zone is characterized by low values of the magnetic field due to the presence of iron only in silicate form.
Of particular interest are the sub-latitudinal zone of low magnetic field values and slightly elevated values of the gravity field, covering the basins of the Kulyumber river and stretching to the upper flows of the Kureyka river (to the east from Snezhnogorsk city, Figure 2). According to geophysical data this zone is close to the Norilsk-Igarka zone and corresponds to the sub-latitude Dyupkun branch Minerals 2020, 10, x FOR PEER REVIEW 6 of 39 Figure 3. Geological map of the Kulyumber river valley (after Ltd. Norilskgeology data with the authors' corrections). Position of the studied area is shown in Figure 1.

Volcanic Rocks
Volcanic rocks (mainly lavas and tuffs) occur in the eastern part of the area, within the Nirungdinsky trough ( Figure 3). Five formations were recognized: Syverminsky, Gudchikhinsky, Khakanchansky, Tuklonsky and Nadezhdinsky. We studied in detail the structure and composition of the lower part of the lavas sequence in the Khalil site, Section 1 ( Figure 4) and compare it with the lower part of the volcanic section in the eastern Norilsk area, Section 2, Lake Lama. Several samples were also analyzed at the Kaya site ( Figure 5).

Structure of Tuff-Lava Sequence and Petrography of Volcanic Rocks
The volcanic Section 1 ( Figure 4, points X-20-X-27) in the Khalil site ( Figure 6) begins with the Syverminsky formation which consists of eight lava flows (total thickness of 140 m). The flows have an average thickness of 10-15 m. They are overlain by olivine basalts of the Gudchikhinsky formation. Two sills cut these basalts.

Volcanic Rocks
Volcanic rocks (mainly lavas and tuffs) occur in the eastern part of the area, within the Nirungdinsky trough ( Figure 3). Five formations were recognized: Syverminsky, Gudchikhinsky, Khakanchansky, Tuklonsky and Nadezhdinsky. We studied in detail the structure and composition of the lower part of the lavas sequence in the Khalil site, Section 1 ( Figure 4) and compare it with the lower part of the volcanic section in the eastern Norilsk area, Section 2, Lake Lama. Several samples were also analyzed at the Kaya site ( Figure 5).

Structure of Tuff-Lava Sequence and Petrography of Volcanic Rocks
The volcanic Section 1 ( Figure 4, points X-20-X-27) in the Khalil site ( Figure 6) begins with the Syverminsky formation which consists of eight lava flows (total thickness of 140 m). The flows have an average thickness of 10-15 m. They are overlain by olivine basalts of the Gudchikhinsky formation. Two sills cut these basalts.
Minerals 2020, 10, x FOR PEER REVIEW 9 of 39 Figure 6. Sections of volcanic rocks in the Khalil area (a, Section 1) and the eastern Norilsk area (b, Section 2 shown in Figure 1 and in [58]).  Figure 1 and in [58].
The Syverminsky formation comprises middle-size grained aphyric lavas and tuffs. Poikilophytic varieties dominate in the lower part of the Syverminsky formation while its upper part mostly consists of tholeiitic basalts (Figure 7a,b). Poikilophytic structure characterizes the three lower flows of the section where large clinopyroxene oikocrysts (0.5 mm) comprise elongated plagioclase chadacrysts. Tholeiitic basalts consist of plagioclase (40-55 vol.%), clinopyroxene (15-40 vol.%), volcanic glass (5-15 vol.%), Fe-Ti oxides (1-3 vol.%). Sometimes they are enriched in titanite (Figure 8a,d). Glass is devitrified (Figure 8a,b) and altered (opacitated) and pallagonite occurs in groundmass (Figures 7a  and 8a). Secondary minerals are chlorite, amphibole. Margins of flows are represented by fine-grained amygdaloidal basalts (Figure 7b) The Syverminsky formation comprises middle-size grained aphyric lavas and tuffs. Poikilophytic varieties dominate in the lower part of the Syverminsky formation while its upper part mostly consists of tholeiitic basalts (Figure 7a,bFigure 7a; Figure 7b). Poikilophytic structure characterizes the three lower flows of the section where large clinopyroxene oikocrysts (0.5 mm) comprise elongated plagioclase chadacrysts. Tholeiitic basalts consist of plagioclase (40-55 vol.%), clinopyroxene (15-40 vol.%), volcanic glass (5-15 vol.%), Fe-Ti oxides (1-3 vol.%). Sometimes they are enriched in titanite (Figure 8a,d). Glass is devitrified (Figure 8a,b) and altered (opacitated) and pallagonite occurs in groundmass (Figures 7a and 8a). Secondary minerals are chlorite, amphibole. Margins of flows are represented by fine-grained amygdaloidal basalts (Figure 7b). Plagioclase comprises 52-70 mol% An; clinopyroxene composition varies in narrow range (Wo36-38En46-47Fs15-16); this mineral contains up to, in wt.%, TiO2: 0.52, MnO: 0.31, Al2O3: 3.63, Cr: 0.25 (Table 1).   The maximum thickness of pyroclastic rocks was found on the eastern side of the Khalil river valley. They occur as tuff horizons (first meters thick) or form extrusions consisting of tuff-lavas breccias.   Table 1. We compared the tuff-lava sequence of Section 1 ( Figure 6a) with the lower sequence of volcanic rocks of the north coast of Lake Lama located in the eastern Norilsk area (Section 2). Early we payed particular attention to the middle part of the volcanic sequence, i.e., formations overlapping the Ivakinsky formation in the Khalil area should be underlined (Figure 6b). The Syverminsky basalts lie on the sedimentary rocks of the Tunguska Group (C 2 -P 2 ) and are overlain by basalts of the Gudchikhinsky formation everywhere around the Kulyumber area. The thickness of the Syverminsky formation in both sections is similar (140 and 160 m, respectively). The Lake Lama section consists of 9 flows of Minerals 2020, 10, 409 13 of 39 tholeiitic basalts overlapped by two flows of poikilophytic basalts. The Gudchikhinsky formation in the Lama Lake area consists of picrites [36] or picritic basalts attributed to the middle subformation while in the Khalil area it is represented by olivine basalts (lower subformation).
The studied volcanic rocks of the Kaya site coming from the bottom of the tuff-lava section. They belong to the Syverminsky formation and are represented by tholeiitic and poikilophytic basalts similar to the rocks of Section 1.

Geochemistry of Volcanic Rocks
Major elements. Chemical compositions of the volcanic rocks from the Khalil and Kaya sites of the Kulyumber area and Norilsk area are given in Table 2. Figure 9 shows diagram (SiO 2 -Na 2 O+K 2 O) of the analyzed rocks.   We compared the tuff-lava sequence of Section 1 (Figure 6a) with the lower sequence of volcanic rocks of the north coast of Lake Lama located in the eastern Norilsk area (Section 2). Early we payed particular attention to the middle part of the volcanic sequence, i.e., formations overlapping the Ivakinsky formation in the Khalil area should be underlined (Figure 6b). The Syverminsky basalts lie on the sedimentary rocks of the Tunguska Group (C2-P2) and are overlain by basalts of the Gudchikhinsky formation everywhere around the Kulyumber area. The thickness of the Syverminsky formation in both sections is similar (140 and 160 m, respectively). The Lake Lama section consists of 9 flows of tholeiitic basalts overlapped by two flows of poikilophytic basalts. The Gudchikhinsky formation in the Lama Lake area consists of picrites [36] or picritic basalts attributed to the middle subformation while in the Khalil area it is represented by olivine basalts (lower subformation).
The studied volcanic rocks of the Kaya site coming from the bottom of the tuff-lava section. They belong to the Syverminsky formation and are represented by tholeiitic and poikilophytic basalts similar to the rocks of Section 1.

Geochemistry of Volcanic Rocks
Major elements. Chemical compositions of the volcanic rocks from the Khalil and Kaya sites of the Kulyumber area and Norilsk area are given in Table 2. Figure 9 shows diagram (SiO2-Na2O+K2O) of the analyzed rocks. Yellow, light brown and brown circles are rocks of the Syverminsky formation from different areas (Kylyumber-Khalil and Kaya sites, Norilsk-Lake Lama, respectively). They show a wide range of compositions belonging to basalt, basaltic andesite, trachybasalt, basaltic trachyandesite and trachyandesite ( Figure 9). Most of the samples are concentrated in the basaltic andesite and basaltic trachyandesite fields: these rocks have similar SiO2 content (Figures 9 and 10a) but differ in alkalis concentrations. Rocks compositions of the Syverminsky formation of the Khalil area plot separately due to their enrichment in alkalis as compared with the rocks from the Kaya site and the Norilsk area (Section 2) and even with the Ivakinsky formation from Section 2. Figure 10 demonstrates the difference in alkalis between the Syverminsky rocks of the Khalil, Kaya and Norilsk area ( Figure  10g,h) and the Khalil area rocks being higher in Na2O whereas K2O and P2O5 contents are similar for Yellow, light brown and brown circles are rocks of the Syverminsky formation from different areas (Kylyumber-Khalil and Kaya sites, Norilsk-Lake Lama, respectively). They show a wide range of compositions belonging to basalt, basaltic andesite, trachybasalt, basaltic trachyandesite and trachyandesite ( Figure 9). Most of the samples are concentrated in the basaltic andesite and basaltic trachyandesite fields: these rocks have similar SiO 2 content (Figures 9 and 10a) but differ in alkalis concentrations. Rocks compositions of the Syverminsky formation of the Khalil area plot separately due to their enrichment in alkalis as compared with the rocks from the Kaya site and the Norilsk area (Section 2) and even with the Ivakinsky formation from Section 2. Figure 10 demonstrates the difference in alkalis between the Syverminsky rocks of the Khalil, Kaya and Norilsk area (Figure 10g,h) and the Khalil area rocks being higher in Na 2 O whereas K 2 O and P 2 O 5 contents are similar for all three areas (Figure 10f-h). Concentrations of other oxides in the Syverminsky rocks from different areas are very similar (Figure 10b-e).  We have studied only one sample from the Gudchikhinsky formation in the Khalil site (Section 1), which was compared with two samples from the same formation in the Norilsk area (green triangles in Figure 9). This sample of the Gudchikhinsky formation correlates with the rocks of the lower part of this formation in the Norilsk area (lower subformation, T 1 gd 1 ) represented by basalts. The basalts of the Khalil area show MgO = 7.96 wt.%) and Na 2 O and K 2 O (2.80 wt.% and 1.33 wt.%, respectively) (Figure 10g,h). These alkaline contents are not typical of the Gudchikhinsky formation in the Norilsk area (Table 2, Nos 23 and 24) where the total alkali content do not exceed 1.5 wt.%.
Tuffs compositions fall in the field of andesite and basaltic andesite (Figure 9, red symbols). They are characterized by low TiO 2 and high Al 2 O 3 . We attributed pyroclastic rocks to different formations based on trace element distribution (see later).
Trace elements. The patterns of trace elements for the Syverminsky rocks of both sections and of rocks outcropping in the Kaya river valley (Figure 11a; Table 2 Here and in Figures 15 and 16 elements are normalized to primitive mantle after [60]. We have studied only one sample from the Gudchikhinsky formation in the Khalil site (Section 1), which was compared with two samples from the same formation in the Norilsk area (green triangles in Figure 9). This sample of the Gudchikhinsky formation correlates with the rocks of the lower part of this formation in the Norilsk area (lower subformation, T1gd1) represented by basalts. The basalts of the Khalil area show MgO = 7.96 wt.%) and Na2O and K2O (2.80 wt.% and 1.33 wt.%, respectively) (Figure 10g,h). These alkaline contents are not typical of the Gudchikhinsky formation There is a difference between the rocks of the Gudchikhinsky formation of the Khalil and Norilsk areas in trace elements. The Khalil rocks contain less Th and U in comparison with the similar rocks of the Norilsk area where there is essential variations in the content of many elements [27], but the Th and U content do not almost change. Thus, the low Th and U concentrations are a specific feature of the Gudchikhinsky basalts of the Khalil area (Figure 11b).
We studied two tuff samples, one is located in the tuff horizon (X-0) and the second one was taken from the volcanic structure outcropping in the Khalil site (X-37). They have similar compositions in term of major components (high SiO 2 content) but different rare elements distributions that allowed us to attribute them to different formations (Figure 11c). Sample X-37 is characterized by a deep slope of the right part of the pattern with (Gd/Yb)n = 1.9 and was attributed to the Syverminsky formation. Sample X-0 has a different pattern (with low (Gd/Yb)n = 1.5) and belongs to the Khakanchansky formation.

Isotope Composition of Volcanic Rocks
We have analyzed eight samples from the Khalil and Norilsk areas, mostly from the Syverminsky formation. The Gudchikhinsky picrites and Khakanchansky tuffs and Tuklonsky basalts have also been analyzed analyzed (Table 3, Figure 12).
Data for flood basalts demonstrate a wide range of Nd-, Sr-isotopic compositions, these variations are typical of the Siberian trap province [30,31,61]. formation. Sample X-0 has a different pattern (with low (Gd/Yb)n = 1.5) and belongs to the Khakanchansky formation.

Isotope Composition of Volcanic Rocks
We have analyzed eight samples from the Khalil and Norilsk areas, mostly from the Syverminsky formation. The Gudchikhinsky picrites and Khakanchansky tuffs and Tuklonsky basalts have also been analyzed analyzed (Table 3, Figure 12).

General Characteristic of Intrusive Rocks
Intrusions represent almost the 10% of the studied area ( Figure 3). They are represented by sills or sill-like bodies located in sedimentary rocks (Devonian and Tunguska Group, C 2 -P 2 ), i.e., they extend in submeridional (N-S) direction at 15-20 km and fall to the east at 10 • -12 • of the Nirungdinsky trough's center. Several intrusive bodies form dikes, which cut volcanic rocks. The thickness of intrusions varies from 5-10 m to 150 m. The attribution of the intrusive rocks of the studied area to the intrusive plutonic complexes of the literature is a problem that is not easy to solve. The location of the territory on the border of the two tectonic structures (Figure 1) has led to emplacement of many intrusive bodies that are typical of both environments. According to the Legend to the Geological Map [62] this area comprises 3 complexes: Ergalakhsky (typical of the Norilsk area), Katangsky (typical of the Tunguska syneclise) and Kureysky (located between Kulyumber and Kureyka rivers). The intrusive bodies of these complexes have similar morphology, texture and composition. It is very difficult to distinguish between them without a detailed geochemical study. Traditionally coarse-grained and olivine gabbro-dolerites (with elevated MgO = 8-9 wt.%) were combined into the Kureysky complex while middle-grained olivine-bearing and olivine-free gabbro-dolerites were attributed to the Katangsky complex.
These three types of intrusive rocks, Ergalakhsky, Kureysky and Katangsky, were found at the Khalil site during the geological mapping ( Figure 4). They are roughly equal in percentage. Bodies of the Ergalakhsky complex are localized along the valley of the stream Khalil (Intrusion 1 in Figure 4). A dike-like body of submeridional (N-S) direction located in volcanic rocks was named as the Khalil intrusion and attributed to the Kureysky complex (Intrusion 2 in Figure 4). It is 100-500-m thick, dips to the east (angles 45 • -75 • ) and consists of gabbro-dolerites and leucogabbro-dolerites in the hanging side. One intrusive body, earlier attributed to the Katangsky complex, was attributed by us to the Norilsk complex based on its geochemical composition, which is close to the Norilsk complex (intrusion 3 in Figure 4). Several small subconcordant and nonconcordant intrusive bodies of the Katangsky complex were recognized in the Khalil and Kaya sites (Intrusion 4 in Figure 4; in Section 1, sample X-21, Figure 6).
Kureysky complex. Here we consider a large dike-like massif in the Khalil site which was previously attributed to the Kureysky complex. It was named the Khalil intrusion during geological mapping of this area [62] (Figure 4, intrusion 2). It consists of middle-, large-grained olivine gabbro-dolerites and gabbro-dolerites. The main minerals are olivine, clinopyroxene, plagioclase; rare minerals include biotite, magnetite and apatite. The composition of olivine ( Figure 9, Table 3 Table 4).
The reference object of the Kureysky complex is the Dzhaltulsky massif located several kilometers south of the Kulyumber area. We have analyzed samples from the borehole OKG-13. Olivine from one sample corresponds to the average value of the olivine compositions in the Kylyumber massif: it is characterized by very low concentrations of calcium and nickel (0.07-0.17 wt.% CaO and 0.03-0.06 wt.% NiO) at the same Fo 72-73 of the mineral. Pyroxene is more magnesian-rich (Mg# = 77-78).
Intrusions of the Katangsky complex are usually composed of fine-medium-grained homogeneous gabbro-dolerites and olivine-bearing gabbro-dolerites. They have usually a dolerite, poikilophytic structure where plagioclase laths (0.1-0.3 microns) are enclosed in large grains of pyroxene (1-1.5 mm). The composition of olivine in olivine gabbro-dolerites (sample X-15) is close to the olivine composition from the Kureysky complex but differs in higher nickel contents (up to 0.19 for Fo 73 ; it is close in composition to olivine from the Norilsk 1 intrusion. Plagioclase changes from An 56 to An 70 . The Mg# number of clinopyroxene varies from 53 to 73 and orthopyroxene from 58 to 69. Magnetite contains 12-14 wt.% TiO 2 . The mineral composition of the Khalil intrusion was studied in the sample X-19 and X-35 (Table 4). The compositions of rock-forming minerals from the Katangsky complex will be described in the second part of the article.   Table 4.

Chemical Composition of Intrusive Rocks
The compositions of the different intrusive complexes are given in Table 5 and illustrated in Figure 14. Trachydolerites of the Ergalakhsky complex (Intrusion 1) differ from the other rocks in higher SiO 2 , P 2 O 5 , Na 2 O and especially in TiO 2 and K 2 O contents (up to 2.5 wt.%. Figure 12e-g) whereas rocks from the other complexes contain less than 1 wt.% K 2 O, and low SiO 2 ., MgO and CaO concentrations.   Low MgO contents (<10 wt.%) characterize all studied samples excluding one sample of picritic gabbro-dolerite from the Norilsk 1 intrusion (24.5 wt.% MgO). This sample contains sulfide mineralization and, due to high sulfide content, has high Fe and low CaO and alkalis' concentrations. The Norilsk 1 intrusion is characterized by wide compositional range because it consists of different rocks that form separate layers within the intrusive body. These rocks are enriched in Al 2 O 3 and depleted in MnO, TiO 2 and Fe 2 O 3 in comparison with the other intrusions of the Khalil site excepting the intrusive body 3. There is a good correlation between rocks from the Norilsk 1 and Intrusion 3 ( Figure 4) for all major elements-TiO 2, SiO 2 , alkalis and P 2 O 5 . Only CaO content in gabbro-dolerites in the Intrusion 3 is slightly lower than that in the Norilsk 1 intrusion. We attributed this Intrusion 3 to the Norilsk intrusive complex (Norilsk-type). In both cases the lower contents of TiO 2 in the Norilsk 1 intrusion and Intrusion 3 are typical of the Norilsk complex. According to numerous data, the TiO 2 concentration from the Norilsk 1 intrusion is 0.87 wt.%, i.e., <1% [8,40,66].
Most interesting is the comparison of the compositions of the Kureysky and Katangsky complexes; these complexes are the most common in this area. They consist of gabbro-dolerites belonging to rocks of normal alkalinity (less than for the Ergalakhsky rocks). Figure 14 shows their identical composition in major components and their difference with the Ergalakhsky intrusion. Only Na 2 O contents are higher in respect to the rocks of the Katangsky complex. Probably, due to alteration of the Katangsky rocks.
The large Dzhaltulsky massif is considered as a reference intrusion of the Kureysky complex. We used three samples from this massif to compare these rocks with other intrusions of the Kulyumber area. Our data show that the Dzhaltulsky massif differs from both the Kureysky and Katangsky intrusions by low TiO 2, alkalis, Fe 2 O 3 and high MgO, Al 2 O 3 , CaO. Thus, we believe that intrusions attributed to the Kureysky complex in the Khalil site belong to the Katangsky complex because they differ of the Dzhaltulsky massif in all major elements.
The distribution of trace elements in rocks of different complexes is shown in a series of spider diagrams ( Figure 15) and on binary diagrams ( Figure 16).       In Figure 15 the intrusive rocks show two distinct compositional patterns: (1) enrichment in LILE with small Ta-Nb negative anomaly and depletion in other HFSE and REE elements (Ergalakhsky complex; Figure 15a); (2) strong negative Ta-Nb and Sr anomalies and minor depletion in other HFSE and REE elements (Figure 15b-f). The second pattern has some differences between intrusions. The greatest variation in the content of trace elements characterizes the Norilsk 1 massif (Figure 15b) Minerals 2020, 10, 409 28 of 39 due to different olivine content of the rock: the most depleted patterns are typical of the picritic gabbro-dolerites. The patterns of the Katangsky and Kureysky (Khalil intrusion) complexes are very close to each other (Figure 15c,d). Rocks of the Dzhaltulsky massif are the most depleted in all incompatible elements. Intrusive body of the Norilsk-type (intrusion 3 in Figure 4) approaches to the composition of the Dzhaltulsky massif (Figure 15e,f).
Ratios of indicator elements detect differences between intrusions of the second type. The (Gd/Yb)n and (La/Sm) n ratios reflect the slope steepness of the right and central parts of the patterns (Figure 16a), whereas the (La/Yb)n ratio characterizes the pattern inclination. To use these indicator element ratios see we did not consider samples of the Ergalakhsky complex in Figure 16a because much higher values of these ratios characterize them. This diagram shows that the (U/Nb)n ratio is only significant difference between them (Figure 16b). The highest values are typical of the Khalil intrusion presumably assigned to the Kureysky complex. However, the Dzhaltulsky massif, which is considered the rock-type of this complex, is characterized by the lowest values of the (U/Nb)n ratio. Intrusion 3 defined preliminarily as Norilsk-type intrusion is close to the Dzhaltulsky massif according to this parameter. Therefore, the (U/Nb)n ratio allows separating intrusions of the different complexes.
Comparison of rocks from different massifs in term of ore elements (excluding rocks with rich sulfide mineralization) indicates their same content of metals such as vanadium and chromium (Figure 17a), nickel and copper (Figure 17b), excepting the Ergalakhsky complex, which is depleted in these metals. Intrusion 3 defined preliminarily as Norilsk-type intrusion is close to the Dzhaltulsky massif according to this parameter. Therefore, the (U/Nb)n ratio allows separating intrusions of the different complexes. Comparison of rocks from different massifs in term of ore elements (excluding rocks with rich sulfide mineralization) indicates their same content of metals such as vanadium and chromium (Figure 17a), nickel and copper (Figure 17b), excepting the Ergalakhsky complex, which is depleted in these metals.

Volcanic Rocks
For the first time we have studied in detail the igneous rocks outside the Norilsk district using the modern geochemical methods. The volcanic rocks of the Kulyumber river valley have been previously subdivided into several formations similar to the formations of the Norilsk region (Table 5, [55]). The Ivakinsky formation is absent in the Nirungdinsky trough. Its composition studied in the Lama Lake area (Section 2) corresponds to the composition of the Ivakinsky rocks in the Kharaelakh trough [29,32]. The Syverminsky formation has been studied also in two sites of the Kulyumber river valley, the Khail and Kaya. Our data allowed to conclude that this formation varies along its stretch. Although the textural and structural features of rocks vary slightly, their chemical composition change significantly. In the Norilsk region (Lama Lake, Section 2), the rocks show a constant homogeneous composition, with small variations in the concentrations of the major elements such as MgO-5-6 wt.%, TiO2-1.6-2 wt.%, alkalis-2.6-3.7. Similar variations of these elements have been found in the Kharaelakh trough [29] and Sunduk mountain [32]. Within the Kulyumber river valley, the contents of major oxides in rocks change more significantly. Trachydolerites of Section 1 contain (wt.%): MgO = 3.4-6.9, TiO2 = 1.2-1.8, alkalis = 5.2-6.3, whereas those (samples 3-2-2 and 3-3-3) at the Khalil significantly differ showing low magnesium (3.4 wt.% MgO) and high alkalis (up to 7.1 Na2O + K2O wt.%) content. Probably, these differences are due to secondary alterations of the rocks as suggested by the presence of many pyrite crystals. The homogenous compositions are typical of the Syverminsky rocks with respect to the La/Sm, Gd/Yb and U/Nb ratios, especially of the rocks at the Nirungdinsky trough.
The Gudchikhinsky formation in the Norilsk region is subdivided into three subformations [14,55] that differ in composition. They are represented by porphyry basalts, picrites and glomeroporphyritic basalts, respectively. Picrites of the middle subformation are the most stable along the stretch, while the basalts of the lower and upper subformations pinch out from the west to east. The studied sample from the Kulyumber river valley corresponds to the rocks of the lower

Volcanic Rocks
For the first time we have studied in detail the igneous rocks outside the Norilsk district using the modern geochemical methods. The volcanic rocks of the Kulyumber river valley have been previously subdivided into several formations similar to the formations of the Norilsk region (Table 5, [55]). The Ivakinsky formation is absent in the Nirungdinsky trough. Its composition studied in the Lama Lake area (Section 2) corresponds to the composition of the Ivakinsky rocks in the Kharaelakh trough [29,32]. The Syverminsky formation has been studied also in two sites of the Kulyumber river valley, the Khail and Kaya. Our data allowed to conclude that this formation varies along its stretch. Although the textural and structural features of rocks vary slightly, their chemical composition change significantly. In the Norilsk region (Lama Lake, Section 2), the rocks show a constant homogeneous composition, with small variations in the concentrations of the major elements such as MgO-5-6 wt.%, TiO 2-1.6-2 wt.%, alkalis-2.6-3.7. Similar variations of these elements have been found in the Kharaelakh trough [29] and Sunduk mountain [32]. Within the Kulyumber river valley, the contents of major oxides in rocks change more significantly. Trachydolerites of Section 1 contain (wt.%): MgO = 3.4-6.9, TiO 2 = 1.2-1.8, alkalis = 5.2-6.3, whereas those (samples 3-2-2 and 3-3-3) at the Khalil significantly differ showing low magnesium (3.4 wt.% MgO) and high alkalis (up to 7.1 Na 2 O + K 2 O wt.%) content. Probably, these differences are due to secondary alterations of the rocks as suggested by the presence of many pyrite crystals. The homogenous compositions are typical of the Syverminsky rocks with respect to the La/Sm, Gd/Yb and U/Nb ratios, especially of the rocks at the Nirungdinsky trough.
The Gudchikhinsky formation in the Norilsk region is subdivided into three subformations [14,55] that differ in composition. They are represented by porphyry basalts, picrites and glomeroporphyritic basalts, respectively. Picrites of the middle subformation are the most stable along the stretch, while the basalts of the lower and upper subformations pinch out from the west to east. The studied sample from the Kulyumber river valley corresponds to the rocks of the lower subformation, although it differs from them in its aphyric texture. The Gudchikhinsky rocks (lower and middle subformations) were studied within the Kharaelakh [29] and Norilsk [67] troughs. The studied basalts of the Khalil site contain more MgO (8 wt.%) in comparison with similar basalts of the Norilsk region, where the Gudchikhinsky basalts contain 6.3 ± 0.5 wt.% in average (67 analyses, [68]). Similar data were obtained for basalts in the Kharaelakh trough, where MgO varies in narrow range of 6-7 wt.%, The basalts of the lower subformation of the Gudchikhinsky formation in the Norilsk trough (boreholes OM-6 and OM-25, [67]) show less magnesium content, with MgO ranges from 4.3 wt.% to 4.8 wt.%.
The trace element ratios show the differences between the rocks of the studied areas more clearly. The (Gd/Yb)n ratio reflecting the inclination of the right part of patterns changes from 2.0 to 2.3 in rocks of the Norilsk area and it is close to the ratio established for the Gudchikhinsky rocks in the Khalil site (2.0). The La/Sm ratio is more variable, and it demonstrates the differences between rocks of the Gudchikhinsky formation located in the different tectonic structures of the Norilsk area. The picrites of the Kharaelkh trough are characterized by values of 1.3-1.5, while basalts have 1.5-1.7. The highest difference between basalts of the lower subformation and picrites of the middle subformation occurs in the rocks of the Norilsk rough, where they are characterized by (La/Sm)n ratio of 2.2 and 1.3-1.7, respectively. The U/Nb ratio has a wide range in the Gudchikhinsky rocks. The basalts of the Norilsk trough have low ratio (1.4), while the picrites have (U/Nb)n = 2.0-4.0. The difference between basalts and picrites within the Kharaelakh trough is also significant, 0.9 for basalts and 1.7-2.2 for picrites. These values are very close to the values that we obtained for basalts of the Nirungdinsky trough (0.7) and for picrites of the eastern part of the Norilsk area (1.8 and 2.4).
The Khakanchansky formation consists of tuffs. The data on the Khakanchansky rocks in literature are restricted in comparison with the data on other formations. Previously we described these rocks from the upper part of Section 2 [58], Lama Lake area and in the Norilsk trough [67]. The contents of major elements in tuffs vary significantly but the distribution of trace elements and their ratios are quite homogenous (La/Sm)n = 2.3-3.3, (Gd/Yb)n = 1.5-1.6, (U/Nb)n = 4.7-6. The composition of tuff from the Khalil area is close to the compositions of tuffs in the Norilsk area (excepting U/Nb). Published data on the composition of the Khakanchansky rocks at the Khantayskoe Lake [69] significantly differ in the La/Sm ratio, which varies from 1.5 to 1.8. According to our data, similar values are typical of the Tuklonsky formation, so it is likely that these tuffs belong to this formation.
In summary, the obtained geochemical data on the volcanic rocks in the Kulyumber river valley allows us to attribute them to the known formations despite some variations in their compositions.

Intrusive Rocks
The diagnostic of intrusive rocks in the Kulyumber river valley is inconveniently due to their similar texture and structure and close chemical composition. Preliminary attribution to intrusive complexes is based on the morphology and inner structure of intrusive bodies. Geological map [53] indicates that three intrusive complexes occur in this area, i.e., Ergalakhsky, Katangsky and Kureysky. We have studied four main intrusions in the Khalil site ( Figure 4) and two intrusions from the Norilsk area (the Norilsk 1) and Kureyka river valley (the Dzhaltulsky massif) representing the Norilsk and the Kureysky complexes. To compare these intrusions between them and with standard Norilsk and Kureysky complexes we used the same parameters as we used for the volcanic rocks (Tables 6 and 7). Fabric of the intrusive rocks belonging to different complexes are similar and are typical of the igneous rocks crystallized near the surface (subvolcanic bodies). Poikilophytic texture dominates and dolerite texture is widespread. Compositions of rock-forming minerals reflect the composition of primary magmas that formed the intrusive bodies. The most magnesium minerals (olivines, pyroxens) are typical of the intrusions with high MgO contents, i.e., Norilsk 1 and Dzhaltulsky massifs. Plagioclases demonstrate a huge range of compositions, so the comparison of different intrusive bodies in the respect of this mineral needs a large number of analyses. We believe that chemical compositions of rocks allow better to distinguish intrusions.
The composition of major elements in the intrusions demonstrates the essential difference between the Ergalakhsky complex and the other rocks. The intrusion of the Ergalakhsky complex contains very low MgO (3.1wt.%) and high TiO 2 (2.4 wt.%) and alkalis (5.4 wt.%). Three intrusions have similar compositions in respect of these elements (high MgO, low TiO 2 and Na 2 O + K 2 O), they are the Norilsk 1 intrusion, the Dzhaltulsky massif. Although the Intrusion 3 in the Khalil site is not differentiated, we referred it to the Norilsk complex based on major elements contents. The Intrusion 2 preliminary attributed to the Kureysky complex differs from the Dzhaltulsky massif (standard object of this complex) in lower MgO and higher TiO 2 and it is close to the Katangsky complex (Intrusion 4) in respect of these elements. The comparison of the studied intrusive bodies in trace elements shows the principle difference between the Ergalakhsky complex and the other complexes as well. The rocks of the Ergalakhsky complex are characterized by very high Gd/Yb and especially La/Yb ratios reflecting the patterns inclination to X axis. The rocks of the other complexes have no such dramatical differences and form one range of compositions with varying characteristics (trace element ratios). There is no intrusive bodies in the Kulyumber river valley which is in a good correlation with the Dzhaltulsky massif. We concluded that the intrusions preliminary attributed to the Kureysky complex are similar to the intrusions of the Katangsky complex.
This aspect of the subdivision of magmatic rocks in the Kulyumber river valley will be regarded in the next part of the article where many intrusions of the Katangsky complex are widespread.

Magmatic Evolution
The geochemistry of igneous rocks located within the Kulyumber area demonstrates their origin from two principally different primary magmas. The first type of magmas has a mantle source which produced the volcanic rocks of the Gudchikhinsky formation. These rocks are characterized by high MgO content in comparison with the other volcanic rocks and absence of Ta-Nb negative and Pb positive anomalies, high εNd (+4) and low 87 Sr/ 86 Sr (0.704). These isotope characteristics allow to attribute the Gudchikhinsky rocks to products of mantle magmas. We previously studied the compositions of these magmas in the Gudchikhinsky picrites of the Norilsk area on the basis of data on melt inclusions in olivines [27]. We highlighted that the magma composition changes from the east to the west in the Norilsk area varying from a primitive mantle composition to a mantle enriched in the crustal material (up to 10%). The parental magma contains 12-14 wt.% MgO and was undersaturated in water. It crystallized near surface at T−1200 • C. Based on high Ni in olivine it was suggested that pyroxenite was a source for this magma [36]. The Gudchikhinsky basalts of the Kulymber valley river show similar geochemical features of the rocks of the Norilsk area but they enriched in LREE (La, Sm especially) and contain very low U and Th. Thus, there is an evolution of magma composition not only in sub latitudinal direction but also in sub meridional from the Norilsk area to the Kulyumber river. This type of magma is an analog of the Hawaiian tholeiites and represent an OIB source (Ocean Island Basalts) [27]. There was not found intrusive analogs of the Gudchikhinsky volcanic rocks within the Khalil and Kaya sites while they were found in the Norilsk area (the Fokinsky intrusive complex [31]).
The second type of magmas is substantially different from that which formed the Gudchikhinsky rocks. These type of magmas combines magmas of many volcanic and intrusive rocks in the Kulyumber area. It is characterized by negative Ta-Nb, positive Pb anomalies and very low εNd (up to −4) and high 87 Sr/ 86 Sr (reaching up 0.707) that reflect a crustal origin (WPB). It was suggested [29][30][31][32] that a sublithospheric mantle contaminated by the crust material (at 25-30 wt.%) is the source magma for these rocks. We recognized two different rock groups formed by this way. The first group combine the Syverminsky trakhybasalts and Ergalakhsky intrusive rocks while the second group comprises the volcanic rocks of the Khakanchansky, Tuklonsky and Nadezhdinsky formations and intrusions of the Katangsky, Kureysky and Norilsk complexes. The rocks of the second group were crystallized from tholeiitic magmas of normal or with relevantly high magnesium (MgO = 6−8 wt.%) at T-1200-1250 • C and contain <1% H 2 O [70,71]. The magmas that formed the intrusions in the Kulyumber river valley are characterized by high TiO 2 (1.3-2.0) in comparison with the magmas of the Norilsk area where only the Daldykansky rocks have similar TiO 2 contents. During the magmatic evolution the TiO 2 content increased with time in the Norilsk area (as shown by the volcanic rock sequence, [72] and, partially by the intrusive rocks). On the basis of these data we suggest that the intrusive bodies of the Kulyumber area were formed at the end of the magmatic evolution within the Siberian platform while volcanic rocks are in good agreement with the volcanic rocks of the Norilsk are and coeval with them.
Thus, we suggest the occurrence of both type of magmatism-rift (OIB) and platform (WPB) in the Kulyumber area.

1.
Geochemical and mineralogical data on the volcanic rocks in the Khalil and Kaya sites of the Kulyumber area allow their attribution to the Syverminsky, Gudchikhinsky, Khakanchansky, Tuklonsky and Nadezhdinsky formations. No rocks of the Ivakinsky formation were found in the studied area, which pinches out to the south from the Norilsk area. Thus, volcanic rocks of the Syverminsky formation overlap sub-concordantly the terrigenous deposits of the Tunguska Group. 2.
The composition of the Guchikhinskiy formation changes from the Norilsk district to the Kulyumber river valley: rocks are enriched in alkalis (up to 4.3 wt.% in total) in the Khalil area (in comparison with 1.15% in the Norilsk area) and depleted in the Th and U contents, respectively. These data reflect an evolution of the formation within the Norilsk-Igarka paleo rift zone.

3.
Four intrusive complexes were recognized within the Khalil area, including Ergalakhsky, Kureysky, Katangsky and Norilsk. The latter is represented by one dyke-like body of gabbro-dolerites attributed to the Norilsk complex on the basis of it low TiO 2 and high Cr contents although the (U/Nb)n ratio is low.
A detailed analysis of intrusive rocks and discussion on their origin is given in the second part of the article devoted to the study of rocks of the southern part of the territory (the Kulyumber site), where igneouse products dominate over volcanic. The petrographic and geochemical study allowed their better attribution to the various intrusive complexes recognized in this area of the Siberian province.