Petrological Features of the Burlakski and Nizhne ‐ Derbinsk Mafic ‐ Ultramafic Plutons (East Sayan Mountains, Siberia, Russia)

: The Nizhne ‐ Derbinsk mafic ‐ ultramafic complex is located between the Central Asian Orogenic Belt and the Siberian Craton and, is associated with the Ballyk fault. The largest, spatially related to each other, plutons in the central part of the complex are the Burlakski and Nizhne ‐ Derbinsk. Rocks in the main units of these plutons are divided into three groups: peridotites (ultramafic), pyroxenites (sub ‐ ultramafic), and gabbroic rocks (mafic). The ultramafic and sub ‐ ultramafic cumulate series are devoid of plagioclase and contain <3 vol. % chromian spinel. The Fo content of olivine in the sub ‐ ultramafic cumulates from both plutons ranges from Fo79 to Fo86. The En content [= Mg/(Mg + Fe + Ca) × 100 atomic ratio] of clinopyroxenes and orthopyroxenes varies from 46–56, and 63–80, respectively. Plagioclase corresponds to labradorite with An contents between 55 and 57. Hornblende is compositionally similar to pargasite. The sequence of change of rock units corresponds to the paragenesis: olivine − olivine + clinopyroxene (orthopyroxene) − clinopyroxene + orthopyroxene – clinopyroxene + orthopyroxene + plagioclase – orthopyroxene. Petrographic, mineralogical, and mineral chemical features of the Burlakski and Nizhne ‐ Derbinsk plutons suggest that the diversity of the material composition of these plutons is due to the processes of magmatic differentiation in deep ‐ seated conditions. Estimates of crystallization pressures and temperatures of the Burlakski and Nizhne ‐ Derbinsk plutons suggest that they crystallized at high pressures ≥ 10kb and temperatures ranging from 1000–1400 °C. Mineralogical and petrological features suggest that the mafic ‐ ultramafic cumulates were derived from a high ‐ Mg basaltic magma. The presence of magmatic hornblende and hydrous mineral assemblages within the ultramafic cumulates indicates that the parental melts had been enriched in dissolved volatile constituents. Taking into account the age of the gabbronorites of the Burlakski pluton (~490 ± 11.8 Ma), the magmatism likely occurred during the Ordovician collision stage of the evolution of the Central Asian Fold Belt.


Introduction
Numerous researchers [1][2][3][4][5][6][7] have focused their efforts on characterizing the petrogenesis of mafic-ultramafic intrusions, as they are ideal environments for understanding magmatic differentiation. By characterizing the primitive magma composition of a magma series, it is possible to model the composition of the rock from which a melt was formed. Moreover, Cu-Ni-platinumgroup elements (PGE), Cr, and Ti ore deposits are often associated with mafic-ultramafic intrusions [8][9][10][11], so understanding how these intrusions form helps guide ore deposit exploration and extraction. Plutonic rocks of the dunite-clinopyroxenite series (dunites, wehrlites, olivine clinopyroxenites, and clinopyroxenites) are especially important as they are the main constituents of the Moho transition zone (MTZ) and their petrogenesis still remains unclear [12]. There have been many studies that have focused on dunites and clinopyroxenite xenoliths [13][14][15][16], but comparatively few studies on dunite-clinopyroxenite series rocks [14]. The main feature of the Nizhne-Derbinsk complex is the presence of ultramafic and subultramafic cumulates. The virtual absence of plagioclase in these cumulates containing both orthopyroxene and clinopyroxene, as early fractionating phases, is indicative of medium to high-pressure crystal fractionation of primary basaltic melts [17]. Our understanding of these plutons formation mechanisms, therefore, plays a critical role in understanding several petrogenetic issues such as parental and evolution of magma composition, metal and sulfur enrichment. These include the mechanism that determines the internal structure of intrusive bodies, the significance behind the various textures and structures exhibited by magmatic rocks, the mechanisms by which melts differentiate, and the development of migration pathways within intrusive bodies. Thus, according to [18][19][20], deep-seated mafic-ultramafic plutons relate to asthenosphere deep-lithosphere interaction processes prior to the melt-modifying processes that influence the parent melts during ascent toward the shallow crust. Moreover, these complexes have high volatile activity and they are at or close to sulfide saturation [20].
This study provides new petrological data about two poorly studied and poorly exposed maficultramafic plutons that contain dunite-clinopyroxenite series rocks. The Burlakski and Nizhne-Derbinsk plutons have long attracted the attention of researches as their mineralogical, textural, and compositional characteristics make them ideal systems to address petrological questions regarding the formation of such rarely described mafic-ultramafic systems. The petrogenesis of the Nizhne-Derbinsk mafic-ultramafic complex is still debated. Some authors [46] consider the plutons of the Nizhne-Derbinsk complex to be ophiolites related to Proterozoic magmatism, which manifested themselves in the Kuzeevsky greenstone belt. Other researchers [47] argue that the rocks of the complex are derivatives of the gabbroic-monzodioritic magmatism that occurred in the Altai-Sayan folded region. Serdyk S.S. and other geologists [48] do not exclude the fact that both Proterozoic peridotites and Ordovician gabbroic rocks are present in the Nizhne-Derbinsk complex. This contribution presents a detailed description of the rock units of the Nizhne-Derbinsk and Burlakski plutons, and combines this information with detailed mineralogy and mineral chemistry provide a better understanding of the genesis of mafic-ultramafic systems.

Geological Setting of Plutons
The Nizhne-Derbinsk mafic-ultramafic complex is located between the Central Asian Orogenic Belt and the Siberian Craton (Figures 1 and 2a). This area is confined to the northwest part of the Eastern Sayan Mountains, roughly 100 km south of Krasnoyarsk city (Siberia, Russia) ( Figure 2). The Nizhne-Derbinsk belt is roughly 40 km in length and has northeast strike. It consist of groups of spatially related plutons, including: Ashtatski, Azertakski, Burlakski, Nizhne-Derbinsk, Medvezi, Konzulski, and Tubilski ( Figure 2b) [49]. Volokhov et al. [50] considered these plutons as blocks of a large single pluton. Several separate plutons, however, were mapped during geological surveying at scale of 1:50,000 [51]. The plutons are sandwiched between series of tectonic plates broken by transverse faults. Along these faults, entire blocks of different rocks units are displaced. This phenomenon is manifested in various sizes of erosive section of plutons. As a result, the primary structure of the intrusive bodies is disturbed. The block structure of intrusions, as well as sometimes observed phase relationships, lead to complex relationships between individual rock types. The Nizhne-Derbinsk complex is poorly exposed with limited outcrop because of the forest cover, making it difficult to establish relationships between different rock series. According to gravity and magnetic anomalies [52], the Tubilski pluton has the form of a stock and is comprised of peridotite cumulates. It occupies an area of about 9 km 2 . The Azertakski pluton has an ellipsoidal-shape form and contains ultramafic and subultramafic (pyroxenites) cumulates. The Ashtatkki pluton has an area of about 3 km 2 and is located 5 km west of the Azertakski pluton. It occurs as an oval shaped body and intruded shales of the Urman Formation. The pluton has both intrusive and tectonic contacts [50]. It consists mainly of diallagites and websterites; gabbros and peridotites, represented by wehrlites. The pluton is located in the center of the residual local gravity anomaly [52]. The Medvezi pluton has an area of about 1 km 2 and is located 1 km northeast of the Burlakski pluton. It has an oval shape. The Konzhulski pluton is located 5 km northeast of the Burlakski pluton. The shape of the pluton is irregular with an area of 4.5 km 2 . The pluton is composed mainly of clinopyroxenites, with lesser gabbroic rocks.
This study focuses on the Burlakski and Nizhne-Derbinsk plutons located in the central part of the complex and associated with the Ballyk fault. These plutons contain several rock types of the Nizhne-Derbinsk complex and have the potential to host Cu-Ni-PGE mineralization [49,51,53]. Within the area, there are two varieties of metamorphic rocks represented mainly by Derbinsk and Urman Formation of Proterozoic age, and non-metamorphic intrusive complexes, including Derbinsk granitoids of Proterozoic age, Buejul gabbro-monzonite-syenite-granosyenite intrusive complex of Ordovician age, and subvolcanic syenite complex of the Devonian age.  The location of emplacement of many of the Nizhne-Derbinsk plutons was controlled by the northwest-and northeast-oriented Ballyk and Kuvay fault zones [49,50]. The Precambrian crystalline basement comprises rocks of the Derbinsk and Kuvay Formation. The overlying Paleozoic sedimentary rocks are largely composed of Cambrian marble, Devonian sandy-aleuritic rocks, limestones, and marls. Intensive tectonic activity and magmatism occurred during the Proterozoic, Ordovician, and Devonian. The area of the Nizhne-Derbinsk complex belongs to the junction zone of several structures with their own stratigraphic features: Derbinsk, Kuznetsk-East Sayan, East Sayan, Minusinsk-East Sayan, and Minusinsk structural facies areas. The plutons consist of a series of tectonic slices that were cut by transverse downthrows and over faults, along which the entire blocks were displaced [49]. The Niznhe-Derbinsk complex intruded the carbonate-rich shale of the Proterozoic Derbian and Urman formations (Figure 2b,c) and was intruded by Devonian subvolcanic alkaline syenitic rocks. According to the K-Ar data of gabbronorites from the Burlakski pluton [55] the Nizhne-Derbinsk complex is Ordovician in age (~490 ± 11.8 Ma) and was involved in the accretion-collisional development stage of the Central Asian fold [56].

Burlakski Pluton
The Bulakski pluton is located in the middle reaches of the Derbina river ( Figure 2c) and occupies an area of about 16 km 2 , making it the largest pluton of the Nizhne-Derbinsk complex. In planar view, the Burlakski pluton has a lens-like shape (5.5 × 2.5 km) with its major axis oriented roughly NE-SW. It has a vertical thickness of about 2-2.5 km according geophysical data (vertical-component magnetics, gravity survey, electroprospecting) [51][52][53] and represents a lopolithic body whose contact dips toward the center at the angles ranging from 65° to 5°. In contrast, the Burlakski pluton has a cup-shaped geometry. The pluton intrudes rocks of the Urman Formation and forms aureoles of cherts, quartz-plagioclase-epidote, plagioclase-pyroxene-cordierite, garnet, and rutile that are tens of meters wide. These aureoles are most readily visible on the eastern flanks of the pluton. The rocks and structural features of the Burlakski pluton make it an ideal environment to understand the origin of the Nizhne-Derbinsk complex, because it is the largest pluton with the widest range of compositional variations. Based on the field observation data of [49,50,55] we interpreted the Burlakski pluton as a layered mafic-ultramafic body comprising serpentinites, wehrlites, pyroxenites, peridotite cumulates, and gabbroic rocks ( Figure 3). The central portion of the pluton consists of the cumulate peridotite series, whereas the marginal areas comprise gabbros and gabbronorites ( Figure  3b). Dykes that crosscut the Burlakski pluton are not widely distributed and comprise both undifferentiated and differentiated rocks ( Figure 2). Undifferentiated rocks include porphyritic altered gabbro-porphyrites, micro-granodiorites, and granodiorites. Differentiated rocks include spessartites and vogesites. Displacement of pluton blocks by oblique-slip faulting caused various levels of the pluton to become exposed at the surface, with the basal portions being exposed in the south-eastern and eastern parts. The basal portion of the Burlakski pluton are characterized by rhythmic layering, with layers having a thickness of 10-15 to 25-30 m of serpentinized dunites, clinopyroxenites, and wehrlites [49], with wehrlites being the predominant rock type (section C-D in Figure 3b). Pyroxenites and serpentinites occur as lenticular bodies within the wehrlites.
Upwards in the stratigraphic sequence, pyroxenites (clinopyroxenites and websterites) that contain lens-shaped bodies of trachytic gabbronorites occur instead of wehrlites. Wehrlites and pyroxenites generally have a massive texture and are characterized by magnetic properties. Accordingly, these rocks have been well mapped using magnetic techniques [49,51]. The magnetic properties are due to the higher content of primary magnetite and secondary magnetire that formed as a result of serpentinization of olivine. According to field observation data of [57] the cross-cutting intrusive relationships of the ultramafic cumulates with gabbronorites was established. A sharp crosscutting contact was observed in the northern part of Burlak Mountain: a thin pyroxenite zone (up to 2 cm) occurs at this contact. A chilled margin was detected on the side of coarse-grained gabbronorites [57]. These observations, together with petrographic data, indicate that coarse-grained gabbronoritesa later injection of magma. This interpretation can also explain the absence of gabbroic rocks in all other plutons of the Nizhne-Derbinsk complex. Gabbronorites of the Burlakski pluton are characterized by a well-defined trachytic texture. The gabbronorite layers strike mainly to the northeast, and dip to the north and northwest at an angle of 45-80°, similar to the bedding of the host rocks. The orientation of the aligned plagioclase crystals typically coincides with the orientation of the layers. Izokh et al. [58] identified taxitic gabbroic rocks with large porphyritic plagioclase crystals, as well as small chert xenoliths at intrusive contacts near the right flank of the Fadeev stream (left intake of the Derbina river). Copper and Pd mineralization is confined to this contact zone [58].

Nizhne-Derbinsk Pluton
The Nizhne-Derbinsk pluton is located on the left flank of the Derbina River and has an area of about 4 km 2 . The vertical thickness of the pluton does not exceed 600 m according to geophysical data [51][52][53]. In plain view, it has an oval shape and is elongated in the NE-SW direction (~2.7 × 2 km). The contacts of the pluton dip under it at an angle of roughly 60°, similarly to contacts in the Burlakski pluton. The pluton, which comprises wehrlites, amphibole pyroxenites, and minor hornblendites, intruded the carbonate-rich shales of the Urman Formation. Gabbroic rocks are absent.
Using previous field observations data about the internal structure of the plutons from Volokhov et al. [50] and Ekhanin et al. [49], Cherkasova and Chernyshov [59] created a schematic cross section of the Nizhne-Derbinsk pluton (Figure 2b, Section A-B). According to the cross section, we interpreted that the pluton has a layered internal structure, with wehrlites occurring at the base of the pluton and hornblende pyroxenites (clinopyroxenites and websterites) occurring in the upper portions of the pluton. Pyroxene in these rocks has typically been replaced by amphibole (uralitization and tremolitization). Small lenticular bodies (less than 100 m thick) of wehrlites were identified by [49] within the hornblende pyroxenites. The hornblendites in the southwestern part of the pluton likely represents the top of the Nizhne-Derbinsk pluton.
At the contact between the Nizhne-Derbinsk pluton and the carbonate-rich shales of the Urman Formation, the shales were metamorphosed to hornfels. Additionally, proximal to this contact, the shales and hornfels are crosscut by small dykes of granitoid rock. The granitoid melts were likely generated because of the high temperatures of contact metamorphism since granitoids are completely absent away from the contact. Similar contact melting are observed at the contact of the Chineysky massif in Transbaikalia [60], as well as the Bayankol and Bashkymugur plutons of Western Sangilen [55]. In the southwestern part of the Nizhne-Derbinsk pluton, Izokh et al. [47] identified hornfels that were crosscut by veins of fine-grained, melanocratic amphibole gabbro and granite, as well as dykes of porpyritic gabbro.

Materials and Methods
Thirty-seven samples of serpentinites, wehrlites, olivine clinopyroxenites, clinopyroxenites, websterites, olivine gabbro, gabbronorites, and hornblendites from the Burlakski and Nizhne-Derbinsk plutons were collected and characterized. Given that outcrop is limited, sampling was carried out partially from outcrops and from deluvial deposits at sampling intervals of 200 m or more using a 1:25,000 topographic base map. The locations where the samples were collected are illustrated in Figure 3a. Using oriented samples, cuts were made in three mutually perpendicular directions to determine the foliation of the rocks. Polished thin sections were prepared at the Optical Laboratory of Tomsk Polytechnic University, Tomsk, Russia. The textural relationships of the minerals were characterized using a Zeiss Axio Scope A1 microscope and a TESCAN VEGA 3 SBU scanning electron microscope operated at 20 kV. Modal mineralogy (in vol. %) was estimated by point counting [61]. Quantitative chemical analyses of olivine, pyroxene, plagioclase, hornblende, and chromian spinels were collected using a Camebax-micro electron microprobe (EPMA) equipped with wavelength-and energy-dispersive spectrometers (WDS and EDS, respectively) at the Analytical Center for Multi-Elemental and Isotope Research of the Sobolev Institute of Geology and Mineralogy in Novosibirsk, Russia. Silicate minerals were analyzed using a focused beam with a spot size of 2 μm, an accelerating voltage of 20 kV, counting times of 10 s, and a beam current of 40 nA. Instrument calibration was achieved using diopside (CaKα), Mn-bearing garnet (FeKα, MnKα), Sr-bearing glass (SrLα), and Babearing glass (BaLα) standards. The limit of detection for SiO2, Al2O3, Na2O, K2O, and CaO was 0.01 wt %, and 0.02 wt % for FeO, MnO, BaO, and SrO. Chromian spinels were analyzed using a focused beam with a spot size of 2 μm, an accelerating voltage of 20 kV, counting times of 10 s, and a beam current of 50 nA. Spectral interferences were corrected for in obtaining net peak intensities. The detection limits for Ti, Al, Cr an Fe, Mn, Mg, Ni, V, and Zn were 0.014 wt %, 0.013 wt %, 0.019 wt %, 0.020 wt %, 0.010 wt %, 0.022 wt % 0.016 wt %, and 0.039 wt %, respectively. The Fe 2+ and Fe 3+ contents of the spinel were estimated following the method of [62]. Cr#, Mg#, and Fe# were calculated using the atomic ratios Cr/(Cr + Al), Mg/(Mg + Fe 2+ ), and Fe 2+ /(Mg + Fe 2+ ) × 100, respectively. The two-pyroxene thermometer of [63][64][65] was used to characterize the crystallization temperature of the plutons. The Al IV and Al VI contents of pyroxene was calculated following the relationship Al IV = 2 − Si and Al VI = Al − Al IV .

Rocks of the Burlakski and Niznhe-Derbinsk Plutons
Rocks of the Burlakski and Nizhne-Derbinsk plutons are compositionally variable and divided into three units: peridotites, pyroxenites, and gabbroic rocks (Figure 4) [56,59].The gabbroic rocks, however, only occur in the Burlakski pluton. The pyroxenites are the most prevalent rock type based on the current level of exposure, comprising roughly half of the outcrop area. Gabbroic rocks account for about a third of the outcrop area, with the peridotites occupying the remaining area. Although peridotites from the Burlakski and Nizhne-Derbinsk plutons both comprise wehrlites and serpentinites, they differ in that the latter contains significant amounts of hornblende.
Wehrlites are the predominant rock type of the peridotite unit in the Burlakski pluton and occur at the base of the differentiated series (Figure 2b). They comprise euhedral-subhedral olivine (~40-55 vol. %) and clinopyroxene (~40-60 vol. %) (Figure 5b). Olivine is represented by subisometric or anhedral that are 1-1.5 mm in size, sometimes up to 3 mm. Olivine are commonly serpentinized, the degree of which varies from 25% to 100%. Fragments of fresh olivine relict grains occur and are 0.2-0.4 mm in size. Clinopyroxene occurs as tabular, elongated or subisometric grains 1-3 mm in size, sometimes up to 5-10 mm. Rarely, olivine and clinopyroxene can exhibit a pokilitic texture. Sulfides include pentlandite-and pyrrhotite with grain size ranging from 0.01-0.3 mm.
The size of clinopyroxene grains are the same as those without hornblende. Hornblende grains have a size of about 2 mm. Biotite also occurs in the hornblende clinpyroxenites as isolated crystals up to 2.5 mm in size. Chlorite occasionally occurs as an overgrowth on biotite. Websterite is widely distributed and account for approximately 15-20 vol. % of the outcrop area of both plutons. They are composed mainly of clinopyroxene (60-90 vol. %) and orthopyroxene (up to 40 vol.%) (Figure 5c,d), but can contain olivine in abundances up to 10 vol. %. Clinopyroxene occurs as elongate or subisometric grains ranging in size from 0.5-5.0 mm, with most being 1-2 mm. Orthopyroxene grains also exhibit a range in sizes from 0.5-4.0 mm. Orthopyroxene in the websterites contain lamellae of clinopyroxene. The websterite in the Nizhne-Derbinsk pluton is different from that in the Burlakski pluton as it contains hornblend, which occurs as anhedral grains ranging in size from 0.5-2.0 mm. Hornblendites only occur in the Nizhne-Derbinsk pluton.    [66].
Mafic rocks in the Burlakski pluton comprise gabbro and gabbronorites (Figure 5f,g) and make up roughly 30% of the outcrop area. Anorthosites are not common in the Burlakski pluton, but do occur as separate thin layers in the leucocratic gabbronorites. In general, the mafic rocks consist of plagioclase (~50-65 vol. %), clinopyroxene (~20-30 vol. %), and orthopyroxene (~5-15 vol. %). They can be divided into i) gabbroic rocks of the layered series that exhibit a gabbroic structure (Figure 4f) and ii) lath-shaped gabbronorites that exhibit an ophitic texture (Figure 4g). Leucocratic olivine gabbro (Figure 4f (Table S1). These minerals are, however, irregularly distributed throughout the olivine gabbro, with some portions lacking melanocratic minerals altogether. Plagioclase in the layered series olivine gabbros occur as short, prismatic, and unzoned crystals. There are two generations of plagioclase in the gabbronorites. The first generation of plagioclase occurs as elongate crystals that typically range in size from 2-5 mm along their longest axis, with some crystals being as long as 25 mm in the lath-shaped gabbronorites. These elongate plagioclase crystals generally exhibit a preferred orientation.
The second generation of plagioclase, as well as clinopyroxene, and orthopyroxene are generally isometric-subisometric, characterized by a smaller grain size than first generation of plagioclase, and occur in the intervals between bands of first generation of plagioclase (Figure 5g). The first generation of plagioclase in the lath-shaped gabbronorites are unique in that they exhibit smooth, curved crystals as a result of plastic deformation. Although pyroxene is less prone to deformation, deformation twins are observed in some crystals.

Mineral description
Olivine is irregularly distributed throughout the rocks of the Burlakski and Nizhne-Derbinsk plutons, and occurs as euhedral to subhedreal crystals that typically range in size from 1-1.5 mm, but can be up to 3 mm in size (Figure 5b). Most olivine crystals have been serpentinized to variable degrees ranging from 25% to 90% (Figure 5a). Fragments of relic olivine do occur, but they are small (0.2-0.4 mm in size). Fractures in olivine crystals are common and filled with lizardite (Figure 5a). Cummingtonite occasionally occurs intergrown with lizardite in these fracture fills.
Clinopyroxene is one of the dominant mineral in the ultramafic rocks of both plutons. It occurs as tabular to subisometric crystals that typically range in size from 1-3 mm ( Figure 5), but can be up to 5-10 mm in size. Clinopyroxene are characterized by well-defined prismatic cleavage and twinning. The boundaries between crystals are rounded and smooth. Small crystals of clinopyroxene (<1.0 mm) are largely subisometric and can either occur interstitial to or as poikilitic inclusions within larger clinopyroxene crystals (Figure 5d). Amphibole commonly occurs as a pseudomorph after clinopyroxene (uralitization).
Orthopyroxene is irregularly distributed throughout the rocks of both plutons, ranging from as low as 3 vol. % in the clinopyroxenites to as high as 30 vol. % in websterites. In both plutons orthopyroxene crystals typically contain fine lamellae of clinopyroxene (Figure 5d). Some crystals of clinopyroxene and orthopyroxene were partially recrystallized. Large pyroxene porphyroclasts occur in a fine-grained (≤0.5 mm) matrix of pyroxene crystals. Brown-yellow iron hydroxides and aggregates of greenish-yellow chlorite are commonly developed along fractures in pyroxenes.
Plagioclase in gabbronorites vary in size, with porphyritic crystals ranging from 4-8 mm in length and rarely up to 25 mm. The porphyritic crystals typically exhibit a preferred orientation (i.e., flow banding) and are plastically deformed, as illustrated by their undulatory extinction and folded twin planes (see Figure 5g). In these deformed crystals, numerous simple shears have divided the porphyritic crystals into separate subgrains.
Hornblende occurs as anhedral crystals that range in size from 0.5-2.0 mm (Figure 5h). They occur as xenomorphic crystals among idiomorphic crystals of clinopyroxene or as rims on clinopyroxene (i.e., corona texture).
Chromian spinel from the Burlakski and Nizhne-Derbinsk plutons are confined to the peridotites and make up 1-3 vol. % of the rock. They occur as small (0.1-0.4 mm) euhedral, and anhedral crystals set in a matrix of serpentine ( Figure 7) and pyroxene. Chromian spinel from the Nizhne-Derbinsk pluton are reddish brown to deep red in color (Figure 7a,b). Some chromian spinels contain inclusions of olivine, pyroxene, and sulfides (Figure 7c,d).
Oxides from the Burlakski and Niznne-Derbinsk plutons are represented by magnetite I, II (Figure 6), and ilmenite.
Pentlandite is less common than pyrrhotite and either occurs as flames in pyrrhotite, or as grains at the rims of the pyrrhotite (Figure 8).
Chalcopyrite is rare. Metasomatic association within serpentinites of the Niznhe-Derbinsk complex was first described by [67]. According to the authors, metasomatic processes that actively took part in the redistribution of iron and nickel should also be involved in the redistribution of PGE + Au as confirmed by the works of, e.g., R. Larsen, M. Fiorentini, etc., [18][19][20] in similarly deep-seated magmatic systems.

Mineral Chemistry
Relic olivines in the Burlakski pluton are Mg-rich (Mg# ~85) and do not differ from relic olivine in the Nizhne-Derbinsk pluton (Mg# ~84). The main difference between olivine in these plutons is that these from the Burlakski pluton are characterized by higher nickel contents (NiO ~0.15 wt %) (Table S2).
Clinopyroxenes from both plutons have similar chemical composition (Table S3, Figure S1). Clinopyroxene in peridotites from both plutons range from augite to diopside (Table S3, Figure S1). In the Burlakski pluton, clinopyroxene in peridotites are compositionally different from those in the mafic units, with the former having lower Mg, Ni, and Cr, and higher Fe, Ti, and Na (Table S3).
Orthopyroxene in both plutons is compositionally represented by bronzite and hypersthene. They are characterized by Fe# that ranges from 18-21 wt % in the subultramafic cumulates of both plutons to 32-33 wt % in the mafic rocks of the Burlakski pluton.
Plagioclase is one of the dominant minerals in the gabbronorites. It corresponds to labradorite with an An content of 55-57.
The chemical composition of chromian spinels from the Burlakski pluton falls within the range of compositions of solid solution between chromite, magnesiochromite, and hercynite [69] (Table S5).
Chromian spinels in the Burlakski and Nizhne-Derbinsk plutons can be divided into three groups based on their Cr# [= Cr/(Cr + Al) × 100 atomic ratio]: type I spinels have Cr# between ~32-38, and type II spinels are characterized by Cr# between~9-10 occur in the Burlakski intrusion, type III spinels are characterized by Cr# between ~65-75 and only occur in the Nizhne-Derbinsk intrusion. Chromian spinels in the peridotites of the Burlakski pluton above the section of the cumulate serieshave lower Cr# than those in the Nizhne-Derbinsk pluton.

Pressure and Temperature Estimates
Estimation of the equilibrium temperature based on the Mg-Fe partition coefficients between olivine and clinopyroxene [70] revealed a temperature of formation of ~1100 °C for the Burlakski pluton and ~900 °C for the Nizhne-Derbinsk pluton. The same temperature estimates for the Nizhne-Derbinsk pluton were obtained using the olivine-spinel geothermometer [71]. Mg-Fe distribution between ortho-and clinopyroxenes may be also used for geothermometry [72], however, orthopyroxene-clinopyroxe pairs correlated with the Cr content, especially in clinopyroxenes, thus KD (= (Xcp × Mg/Xcp × Fe 2+ ) × (Xop × Fe 2+ /Xop × Mg)) within the Niznhe-Derbinsk pluton has significant variation in values and ranges from 0.9 to 1.4 ( Figure S3). The large scattering of KD is improbable according to [73] and can be explained by temperature dependence of the solubility of Cr in pyroxenes coexisting with chromian spinel [74]. According to the Ca/(Ca + Mg) of orthopyroxenes and clinopyroxenes [75] temperatures range from 1000-1382 °C. The temperatures of crystallization obtained using different pyroxene geobarometers and geothermometers [63][64][65] were broadly similar and, on average, range from ~1100-1400 °C (Table S6). The lower crystallization temperatures of the hornblende websterite from the Nizhne-Derbinsk pluton may be the result of amphibolization.
The virtual absence of plagioclase in ultramafic and sub-ultramafic rocks containing both orthopyroxene and clinopyroxene as early fractionating phases is indicative of medium to highpressure crystal fractionation from primary basaltic melts [17]. According to the Al VI /Al IV ratio of orthopyroxenes [76], the rocks of the Burlakski and Nizhne-Derbinsk plutons crystallized at pressures ranging from 14 to 17 kb. Based on data from [77], the Al2O3 and Mg# number of orthopyroxenes and clinopyroxenes fall near or within the high pressure field on the Al2O3 vs. Mg# diagram of [77]. The pressures obtained from these geobarometers may reflect the conditions under which pyroxene crystallized, rather than the conditions under which the rocks crystallized. Pressure estimates using Al IV and Altot content of amphiboles ( Figure S4a) from [78] also indicate high pressure emplacement of the plutons, which is consistent with interpretations made using data from [79,80] ( Figure S4b).
The presence of pargasite in the ultramafic rocks of the Nizhne-Derbinsk pluton and hydrous metasomatic mineral assemblages within the ultramafic cumulates of the both plutons are evidences for the enrichment of the parental melts in dissolved volatile constituents [81].

Petrological Features
According to our observations and the data of previous researchers [49,50,58], dunites and wehrlites occur at the base of the macro rhythm of the Burlakski pluton. Clinopyroxenites and websterites are recorded upsection. The central part of the macro rhythm, with a thickness of about 200 m, is composed of stratified gabbronorites (see Figure 3). The sequence of change of rock units corresponds to the paragenesis: olivine − olivine + clinopyroxene (orthopyroxene) -clinopyroxene + orthopyroxene -clinopyroxene + orthopyroxene + plagioclase -orthopyroxene. In this sequence from ultramafic cumulates to gabbronorites the Mg# of olivine, clinopyroxene, and orthopyroxene ranges from 86-79, 91-72, 82-66, respectively. In the same sequence, the chemical composition of pyroxenes, as the most common minerals in both plutons there is an increase in Ti, Na, Mn. Variation in minor oxides of the pyroxenes is shown in Figure S5, where the Cr2O3 content of clinopyroxene increases, and Na2O, MnO, and TiO2 contents decrease with increasing Mg# [= Mg/(Fe 2+ + Mg) × 100% atomic ratio]. Such trends are also exhibited by MnO and TiO2 in the case of orthopyroxene. An exception is aluminum oxide Al2O3.
Chromian spinels in both plutons have relatively high ZnO contents (up to 1.1 wt %) (Table S5), which may reflect secondary post-magmatic processes (i.e., metamorphism and serpentinization) [38,39,44]. Johan and Ohnenstetter [86] suggested that the enrichment of Zn in chromian spinels may be related to their alteration by high-temperature metasomatic fluids. Such fluid involvement in the peridotite rhythms of the Burlakski and Nizhne-Derbinsk plutons is corroborated by the presence of metasomatic minerals, such as awaruite, native brass, native iron, etc., (Figure 10) [67]. Sulfide liquation is another possible cause of the higher Zn concentrations. Along with low Ni concentrations, higher Zn concentrations can be considered as an indicator of the sulfide segregation [8,87,88]. For example, chromian spinels in the upper zone of the sulfide-bearing Bushveld complex contain about 0.8 wt % ZnO, whereas those in the lower zone contain about 0.5 wt % ZnO [88]. Chromian spinels with elevated ZnO occur in the sulfide ores of Kambalda (0.53-2.92 wt %) [89] and in the sulfide deposits of the Vammala nickel belt (ZnO > 0.8 wt %) [88]. The ZnO contents of spinel from Jinchuan range from 0.1-1.4 wt % ZnO [10]. Chromian spinel from barren plutons in the same region as Jinchuan (e.g., the Zangbutai and Yejili plutons) typically contain less than 0.5 wt % ZnO [10]. The presence of sulfide inclusions in chromian spinels (see Figure 7d) allows assuming the potential equilibrium between the sulfide liquid and the silicate melt. The occurrence of both sulfide and silicate inclusions in chromian spinel (see Figure 7c,d) suggests that sulfide liquated early relative to the crystallization of silicates. The presence of sulfides suggests that the conditions under which these plutons formed may be favorable for the formation of Ni-Cu-PGE mineralization.

Conclusions
Rocks of the main units in the Burlakski and Nizne-Derbinsk plutons were divided into three groups: peridotites (ultramafic), pyroxenites (subultramafic), and gabbroic rocks (mafic). There are two main differences between the rock units in the Burlakski and Nizhne-Derbinsk plutons. First, the rocks from the latter contain significant amounts of hornblende and the pluton itself contains a hornblendite unit. Second, the gabbro and gabbronorite units only occur in the Burlakski pluton.
Formation of the plutons suggested by the predominance of role of the ultramafic and subultramafic cumulates and a lack of more evolved rock types. Pressure and temperature estimates suggest that these plutons formed at pressures ≥10 kb and temperatures ranging from 1000-1400 °C. The presence of magmatic hornblende and hydrous mineral assemblages within the ultramafic cumulates indicates that the parental melts had been enriched in dissolved volatile constituents. The presence of sulfides suggests that the conditions under which these plutons formed may be favorable for the formation of Ni-Cu-PGE mineralization. Petrological and mineralogical features suggest that the rocks are mafic-ultramafic cumulates derived from a high-Mg basalt magma. Taking into account the age of the gabbronorites of the Burlakski pluton (~490 ± 11.8 Ma) the magmatism likely occurred during the Ordovician collision stage of the evolution of the Central Asian Fold Belt.
Author Contributions: T.Y. conceived and designed the study, and wrote the initial draft of the manuscript. T.Y. and A.C. collected the samples and conducted sample preparation for analytical work. T.Y., T.T., O.S., and A.M. analyzed the data, conducted processing graphic and analytical work. M.B. and G.G. revised and corrected the final version. All authors have read and agreed to the published version of the manuscript.