Geoscience Frontiers xxx (2016) 1e22

HOSTED BY Contents lists available at ScienceDirect

China University of Geosciences (Beijing) Geoscience Frontiers

journal homepage: www.elsevier.com/locate/gsf

Focus paper Mid-Miocene thermal impact on the lithosphere by sub-lithospheric convective mantle material: Transition from high- to moderate-Mg magmatism beneath Vitim Plateau, Siberia

Irina Chuvashova a,b,*, Sergei Rasskazov a,b, Tatyana Yasnygina a a Institute of the Earth’s Crust, Siberian Branch of RAS, Irkutsk, Russia b Irkutsk State University, Irkutsk, Russia article info abstract

Article history: High-Mg lavas are characteristic of the mid- volcanism in Inner Asia. In the Vitim Plateau, small Received 7 July 2015 volume high-Mg volcanics erupted at 16e14 Ma, and were followed with voluminous moderate-Mg lavas Received in revised form at 14e13 Ma. In the former unit, we have recorded a sequence of (1) initial basaltic melts, contaminated 10 May 2016 by crustal material, (2) uncontaminated high-Mg basanites and of transitional (KeNaeK) com- Accepted 21 May 2016 positions, and (3) picrobasalts and basalts of K series; in the latter unit a sequence of (1) initial basalts Available online xxx and basaltic andesites of transitional (NaeKeNa) compositions and (2) basalts and trachybasalts of KeNa series. From pressure estimation, we infer that the high-Mg melts were derived from the sub- Keywords: Volcanism lithospheric mantle as deep as 150 km, unlike the moderate-Mg melts that were produced at the e Geodynamics shallow mantle. The 14 13 Ma rock sequence shows that initial melts equilibrated in a garnet-free Cenozoic mantle source with subsequently reduced degree of melting garnet-bearing material. No melting of Asia relatively depleted lithospheric material, evidenced by mantle xenoliths, was involved in melting, Asthenosphere however. We suggest that the studied transition from high- to moderate-Mg magmatism was due to the Lithosphere mid-Miocene thermal impact on the lithosphere by hot sub-lithospheric mantle material from the Transbaikalian low-velocity (melting) domain that had a potential temperature as high as 1510 S. This thermal impact triggered rifting in the lithosphere of the . Ó 2016, China University of Geosciences (Beijing) and Peking University. Production and hosting by Elsevier B.V. This is an open access article under the CC BY-NC-ND license (http://creativecommons.org/ licenses/by-nc-nd/4.0/).

1. Introduction 2011; Wang and Gaetani, 2008; Zhao, 2009; Hanano et al., 2010; Herzberg, 2011). The origin of basalts from continents and other The explanation for the origin of within-plate volcanism from ocean islands, where no high-Mg liquids erupted, is attributed to deep mantle sources is currently based on global seismic tomography extension of the lithosphere and precursor metasomatism without images in combination with results of ultrahigh pressure experi- high temperature processes (Griffin et al., 1988; Barry et al., 2003). ments (Anderson, 2007; Ghosh et al., 2007, 2009; Maruyama et al., The Emeishan and Siberian large igneous provinces (LIPs) that 2007; Karato, 2012). The emphasis in the ongoing debate about hot formed in Asia in the late Permian and at the PermianeTriassic mantle plumes rising from the coreemantle boundary is placed on boundary were dominated by moderate-Mg (moderate-tempera- the crucial value of seismic velocity structure of the mantle ture) basaltic melts with minor high-Mg (high-temperature) mei- (Hofmann, 1997; Foulger, 2010). The structure of the mantle and mechites and, probably, komatiites (Vasiliev and Zolotukhin, 1995; components of volcanic rocks associated with the Hawaiian plume Fedorenko et al., 1996; Fedorenko and Czamanske, 1997; Kogarko are adequately interpreted by high-temperature processes respon- and Ryabchikov, 2000; Hanski et al., 2004; Carlson et al., 2006; sible for generating high-Mg magmatic liquids (Sobolev et al., 2005, Zhang et al., 2006; Sobolev et al., 2009; Rasskazov et al., 2010; Shellnutt, 2014). The Siberian LIP was simulated by a thermo- chemical plume in the upper mantle and transition zone with a * Corresponding author. Institute of the Earth’s Crust, Siberian Branch of RAS, diameter of nearly 1000 km (Sobolev et al., 2011). However, no low- Irkutsk, Russia. velocity region that might correspond to the Siberian LIP is recor- E-mail address: [email protected] (I. Chuvashova). ded (Koulakov and Bushenkova, 2010). Peer-review under responsibility of China University of Geosciences (Beijing). http://dx.doi.org/10.1016/j.gsf.2016.05.011 1674-9871/Ó 2016, China University of Geosciences (Beijing) and Peking University. Production and hosting by Elsevier B.V. This is an open access article under the CC BY-NC- ND license (http://creativecommons.org/licenses/by-nc-nd/4.0/).

Please cite this article in press as: Chuvashova, I., et al., Mid-Miocene thermal impact on the lithosphere by sub-lithospheric convective mantle material: Transition from high- to moderate-Mg magmatism beneath Vitim Plateau, Siberia, Geoscience Frontiers (2016), http://dx.doi.org/ 10.1016/j.gsf.2016.05.011 2 I. Chuvashova et al. / Geoscience Frontiers xxx (2016) 1e22

Cenozoic volcanic fields of Asia are often underlain by low-velocity highlighted by both a drastic lateral change of S-wave velocities at anomalies residing in the upper mantle (Rasskazov et al., 2003a; Lei the depth of 250 km and high- to moderate-Mg transition of lava and Zhao, 2005; Zhao et al., 2009; Chuvashova et al., 2012; Wei compositions (Fig. 1). et al., 2012). Decreasing velocities were firstly detected therein at a Most of the high-Mg low-volume lavas with 12e32 wt.% MgO shallow level of the upper mantle through a temporal delay of seismic (Le Bas, 2000; Kerr and Arndt, 2001) erupted in the outline of the waves from nuclear explosions in Nevada, USA (Rogozhina and Transbaikalian domain in the short time interval of 16e14 Ma. Kozhevnikov, 1979). At the beginning of the 1990s, some volcanic Voluminous moderate-Mg lavas flooded large areas afterward areas in Central and East Asia were referred to deep mantle plumes and lasted for millions of years. In the Kamar volcanic field, a (Nakamura et al., 1990; Yarmolyuk et al.,1990,1994; Rasskazov, 1991, single high-Mg lava flow was found within the 16e15 Ma unit 1994). There was also an attempt to locate mantle plumes through that was preceded and followed by moderate-Mg basalts dated at geophysical data (Zorin et al., 2003, 2006). It was suggested that the 18.1e17.6 and 13e12 Ma, respectively (Rasskazov et al., 2003b, dispersed within-plate Cenozoic volcanism of Asiawas derived due to 2013b). In the Udokan volcanic field, the initial low-volume a broad asthenospheric upwelling from a deeper mantle level (“hot (few km3) eruptions of high-Mg melaleucitites occurred region”)(Tatsumi et al., 1990) or convective flow upraised from the at ca. 14 Ma, whereas the subsequent high-volume (450 km3) coreemantle boundary (“hot mantle field”)(Zonenshain et al., 1991). alkali olivine basaltetrachyte and basanitesetephriphonolite The extensive asthenospheric upwelling was assumed to be marked volcanic series flooded the area of 3000 km2 in the last 8 Ma by the DUPAL isotopic anomaly (Flower et al., 1998). However, from (Rasskazov et al., 1997, 2000). Small outcrops of high-Mg rocks the measured Pb-isotope ratios in basalts, no anomaly was detected in and extensive moderate-Mg volcanic strata were identified also Cenozoic basalts erupted within the Caledonian terranes of Inner Asia, in the Vitim and Dariganga volcanic fields (Ashchepkov, 1991; although the anomaly was traced in lavas erupted through the Ashchepkov et al., 2003; Chuvashova et al., 2012, 2016). In this TuvaeMongolian and other Riphean micro-continents that separated paper, we present results of a comparative study of mid-Miocene from East Gondwana and drifted across Paleo-Asian Ocean from the rocks that erupted in the Vitim volcanic field at 16e14 and 14e13 Southern to Northern Hemisphere in the late Precambrian. It was Ma and exhibited high- and moderate-Mg liquids, respectively. inferred that the anomaly was brought by these tectonic units and could not belong to the deep “hot mantle region” (Rasskazov et al., 2. Geologic background and sampling 2002). The S-wave model of Asia by Yanovskaya and Kozhevnikov The Vitim Plateau is situated within the BaikaleVitim Fold (2003) was used for interpreting spatial-temporal distribution of System, the geological structure of which has long been the object Cenozoic basalts in terms of large low-velocity (melting) domains of lively debate. Comprehensive geological mapping of the area of two mantle levels: the Transbaikalian of 200e410 km and the showed sedimentary and -sedimentary units of late Pre- Sayan-Mongolian, Sea of Okhotsk, and Philippine Sea of cambrian and early designated as Baikalides and early 50e200km(Rasskazov et al., 2003a, 2004). The identities of the Caledonides, respectively (Belichenko, 1977). Later on, the whole SayaneMongolian and Transbaikalian melting domains were area of Western Transbaikal was referred to a single Barguzin

Figure 1. Spatial distribution of high- and moderate-Mg mid-Miocene volcanic rocks in Inner Asia relative to the S-wave velocity image at a depth of 250 km after Yanovskaya and Kozhevnikov (2003). Volcanic fields with mid-Miocene rocks: high-Mg e Vitim (VT), Udokan (UD), Kamar (KM), Dariganga (DR); moderate-Mg e Tunka (TN), Dzhida (DZ), Eastern Sayans (ES), Khubsugul (KH), Dzabkhan (DB), Valley of Lakes (VL), Eastern Hangay (EH), Ugii-Nur (UN), Hannuoba (HN). Next to each of mid-Miocene locations of high- and moderate-Mg rocks, a maximum content of MgO (wt.%) is shown in parentheses.

Please cite this article in press as: Chuvashova, I., et al., Mid-Miocene thermal impact on the lithosphere by sub-lithospheric convective mantle material: Transition from high- to moderate-Mg magmatism beneath Vitim Plateau, Siberia, Geoscience Frontiers (2016), http://dx.doi.org/ 10.1016/j.gsf.2016.05.011 I. Chuvashova et al. / Geoscience Frontiers xxx (2016) 1e22 3 super-terrane, much of which was cross-linked by the 310e280 Ma 2001, 2007; Chernyaeva et al., 2007; Usoltseva et al., 2010). The granites of the huge AngaraeVitim batholith (Melnikov et al.,1994). sequences were assigned as the Kularikta, Dzhilinda, Khoygot, and Due to thorough geological study, the area was subdivided, never- Bereya formations dated back to the , middleeupper theless, into early Paleozoic terranes of different origin (Belichenko Miocene, , and , respectively. et al., 2006). According to this interpretation, the Vitim Plateau is The Kularikta Formation accumulated in conditions of slightly located between the Ikat back-arc turbidite and West Stanovoy dissected relief. No reliable geological and geochronological evi- composite terranes within the Eravna island-arc tectonic unit that dence has so far been obtained on the Oligocene volcanic activity in was affected by volcanism in the and Cenozoic (Fig. 2). the area. The Dzhilinda Formation was deposited within erosional Stratigraphic and paleontological records in small fragments paleovalleys as deep as 500 m that were filled firstly with coarse disseminated in the AngaraeVitim batholith showed also Upper sediments, containing few volcanic layers, and then with lavas and SilurianeCarboniferous sedimentary rocks (Ruzhentsev et al., 2012) lacustrine sediments. From borehole records of the Dzhilinda For- that could originate on a continental margin of the closing Mon- mation, the major Northern, Southern, and Central paleovalleys goleOkhotsk Bay of the Paleo-Pacific(Sengör and Natal’in, 1996; have been traced over 120, 100, and 50 km, respectively. After Parfenov et al., 1999). partial erosion of lavas and sediments of the Dzhilinda Formation, The Cenozoic volcanism of the Vitim Plateau has become known the new erosional channels were buried beneath lavas and sedi- due to numerous publications on deep-seated inclusions from ments of the Khoygot, and Bereya Formations. basanites and picrobasalts in this area (Volyanyuk et al., 1976; The Vitim volcanic field occupies an area of about 5000 km2. The Kiselev et al., 1979; Rasskazov, 1985, 1993; Ashchepkov, 1991; total volume of volcanic rocks exceeds 1500 km3. Slag-lava edifices, Ionov et al., 1993, 2005; Litasov et al., 2000a,b; Aschepkov and eighty eight of which erected on the surface of the volcanic plane Andre, 2002; Litasov and Taniguchi, 2002; Ashchepkov et al., and another eighteen edifices buried beneath lavas, belong to seven 2011; Goncharov and Ionov, 2012). Cenozoic volcanic and large and six small volcanic centers (Fig. 2a). Each large volcano volcano-sedimentary sequences of the Vitim Plateau have been with diameter of 16e20 km consists of between seven to seventeen studied since the beginning of the 1980s using borehole cores for slag-lava edifices. A volcano of smaller size comprises two to four KeAr dating of volcanic layers and paleontological records of sed- slag-lava cones. A comparative study of volcanic sequences in the iments (Rasskazov and Batyrmurzaev, 1985; Rasskazov et al., 2000, eastern and western parts of the Vitim volcanic field demonstrated

Figure 2. Mid-Miocene river valleys buried beneath lavas of the Vitim volcanic field (a) and tectono-stratigraphic terranes of Western Transbaikal; (b) In the scheme (a), the Northern (N), Central (C) and Southern (S) paleovalleys and volcanic centers: Bereya (Br), Sharbahtay (Sh), Ingur (In), Salbulin (Sl), Antase (An), Kolichikan (Kl), Db (Dybryn), Nm (Namaru), Sirinikta (Sr), Muhal (Mh), Yaole (Ya), Khoigot (Kh), and Malyi Amalat (Ma) are shown after Rasskazov et al. (2000). In the scheme (b), the terrane boundaries of the island arc Eravna (Er), back-arc turbidite Ikat (Ik), Riphean BaikaleMuya (BM), and composite West Stanovoy (WS) are designated after Belichenko et al. (2006).

Please cite this article in press as: Chuvashova, I., et al., Mid-Miocene thermal impact on the lithosphere by sub-lithospheric convective mantle material: Transition from high- to moderate-Mg magmatism beneath Vitim Plateau, Siberia, Geoscience Frontiers (2016), http://dx.doi.org/ 10.1016/j.gsf.2016.05.011 4 I. Chuvashova et al. / Geoscience Frontiers xxx (2016) 1e22 similar lava compositions in the units formed between 13.0e9.5 Ma, center. The remaining 30% of the 14e13 Ma unit are basalts and significant difference between the earlier and later rocks trachybasalts that occur above the basaltebasaltic andesite package (Chuvashova et al., 2016). In the Bereya volcanic center, located in in the northern part of the volcanic center. This unit was sampled in the eastern part of the volcanic field, the early volcanic interval is both boreholes and outcrops. exhibited by the 16e14 and 14e13 Ma volcano-sedimentary units. The volcanic center is marked by seven volcanic cones that overtop 3. Sample preparation and analytical procedures a lava shield in an area of at least 250 km2. The representative data on rock compositions of the volcanic center were obtained from Collected samples were prepared for analytical work by crush- both outcrops and borehole cores (Figs. 3 and 4). ing in an agate mortar about 400e500 g of rock fragments. From The 16e14 Ma high-Mg rocks occur as fragments, pillow lavas, high-Mg rocks, xenocrystic material was removed under a binoc- and lenses in the layered tuffaceous-sedimentary strata, exposed in ular microscope. The effect of this removal was tested by analyzing a quarry at the 76th km of the Romanovka to Bagdarin road aliquots with and without xenocrysts for 5 samples selected from (Ashchepkov, 1991). This sequence was formed when the relief has the lower layer of the 16e14 Ma unit (Fig. 5). not yet been dissected. At present the tuffaceous-sedimentary Analytical work was performed in the laboratory for isotopic and strata occupies a relatively high altitude (>1000 m). Fragments of geochronological studies of the Institute of the Earth’s Crust ac- volcanic rocks and pillow lavas are composed primarily of volcanic cording to analytical procedures described by Rasskazov et al. (2012) glass. Some rocks are aphyric; others contain large olivine crystals and Yasnygina et al. (2015). Trace elements were determined by or fragments of disintegrated mantle and crustal rocks. inductively coupled plasma mass spectrometer (ICP-MS) technique The volcanic units were sampled from the tuffaceous- using an Agilent 7500 ce. Samples were prepared by microwave acid sedimentary sequence in the open quarry as it was deepening. digestion with a mixture of HNO3 and HF. For a more complete The top layer was opened in September of 1986, the middle one in decomposition, each sample was re-evaporated with HF and with the mid-1990s, the lower by 2013. Picrobasalts, basalts, and additions H2O2. In addition, samples of high-Mg rocks were fused numerous mantle xenoliths were found in the upper layer, high-Mg with lithium metaborate. Before measurement, In and Bi internal basanites and basalts as well as pillow lavas of the same rocks e in standards were introduced into the sample. Internal standard the middle layer, and high-Mg basalts with inclusions of crustal correction was done by interpolation. Standard reference solutions material e in the lower one (Fig. 5). of HPS (High-Purity Standards) production and geological reference The erosional paleovalleys, deepened into the to the materials of the USGS (BIR-1a, DNC-1a, BHVO-2, AGV-2, and RGM-1) altitude of 850 m, were filled with moderate-Mg lavas erupted were used for calibration and validation of analyses. Measurement between 14e13 Ma. A volume of the accumulated lava unit as thick of Sr and Nd isotope ratios were performed on a Finnigan MAT 262 3 as 150 m exceeded 15 km . About 70% of the stratum exhibited the thermal ionization mass spectrometer (TIMS). To control the quality group of basalts and basaltic andesites throughout the volcanic of measurements, standard reference samples JNd-1 (Japan) and

Figure 3. Location of the quarry and lines of drill holes that exposed, respectively, high- and moderate-Mg rocks in the Bereya volcanic center.

Please cite this article in press as: Chuvashova, I., et al., Mid-Miocene thermal impact on the lithosphere by sub-lithospheric convective mantle material: Transition from high- to moderate-Mg magmatism beneath Vitim Plateau, Siberia, Geoscience Frontiers (2016), http://dx.doi.org/ 10.1016/j.gsf.2016.05.011 I. Chuvashova et al. / Geoscience Frontiers xxx (2016) 1e22 5

Figure 4. Relations between the 14e13 Ma and younger (12.6e0.6 Ma) volcano-sedimentary units in the Bereya volcanic center. The drill-hole lines of the sections BA, BC, and DE are shown in Fig. 3. NBS SRM-987 (USA) were used. The obtained ratio 4. Major oxides 143Nd/144Nd ¼ 0.512104 4(2s) was close to the standard value 0.512103 (JNd-1), the ratio 87Sr/86Sr ¼ 0.710248 0.000013 (2s)was The volcanic rocks from the mid-Miocene units of the Bereya comparable to the standard value 0.71025 (NBS SRM-987) (Faure, volcanic center are well grouped on the classification diagrams K2O 2001). Major oxides were measured by classic “wet chemistry”. vs. K2O/Na2O, K2O vs. SiO2, (Na2O þ K2O) vs. SiO2, and MgO vs. SiO2 Data on volcanic rocks were included into the database (Rasskazov (Fig. 6). Representative compositions of the groups are placed in et al., 2013b). Additional data on volcanic rocks from the studied Table 1. units were used after Ashchepkov (1991), Harris (1998), Litasov and The KeNa ratio is usually used as the main criterion for Taniguchi (2002), Johnson et al. (2005). definition of potassic (K2O > Na2O), sodic (K2O < 0.25 Na2O),

Please cite this article in press as: Chuvashova, I., et al., Mid-Miocene thermal impact on the lithosphere by sub-lithospheric convective mantle material: Transition from high- to moderate-Mg magmatism beneath Vitim Plateau, Siberia, Geoscience Frontiers (2016), http://dx.doi.org/ 10.1016/j.gsf.2016.05.011 6 I. Chuvashova et al. / Geoscience Frontiers xxx (2016) 1e22

and potassicesodic (Na2O > K2O > 0.25 Na2O) series (in wt.%) (Foley et al., 1987; Rollinson, 1993). In the diagram of K2O vs. K2O/ Na2O, the mid-Miocene rocks from the Bereya volcanic center overlap a field of KeNa series and partly occupy potassic and sodic fields. Among the 16e14 Ma rocks, the high-Mg groups of picrobasaltsebasalts and contaminated basaltsetrachybasalts belong to K and KeNa series, respectively. The group of high-Mg basanitesebasalts shows intermediate position, i.e. transition from K to KeNa series. Unlike the UgiieNur volcanic field of Central Mongolia, where potassic affinity was provided by increasing K concentrations in moderate-Mg lavas, the potassic high-Mg series of the Bereya volcanic center was due to relatively decreasing Na content. Among the 14e13 Ma rocks, the group of basaltsetrachybasalts exhibits KeNa series, the group of basaltsebasaltic andesites demonstrates a transition from Na to KeNa series. Similarly, rocks of 16e14 and 14e13 Ma are well separated on Figure 5. Schematic cross-section that demonstrates relations between the 16e14 and the total alkalis vs. silica (TAS) diagram. High-Mg basanitesebasalts 14e13 Ma volcano-sedimentary units in the Bereya volcanic center. Location of the and picrobasaltsebasalts show low SiO2 contents (42.8e45.0 and quarry that exposed sediments with fragments of high-Mg volcanic rocks is shown in 44.2e45.6 wt.%), whereas high-Mg basaltsetrachybasalts demon- Fig. 3. strate the elevated SiO2 concentrations (46.2e50.2 wt.%). The most

Figure 6. Coviriation diagrams of K2O vs. K2O/Na2O(a), K2O vs. SiO2 (b), total alkalis vs. silica (TAS) (c), and MgO vs. SiO2 (d) diagrams for mid-Miocene rocks from the Bereya volcanic center. Among volcanic rocks from the 16e14 Ma layers (Fig. 5), we define three groups: (1) high-Mg basalts and trachybasalts of KeNa series, (2) basanites and basalts of transitional (KeNaeK) compositions, and (3) picrobasalts and basalts of K series. Among volcanic rocks of the 14e13 Ma unit (Fig. 4), we recognize two groups: (1) basalts and basaltic andesites of transitional (NaeKeNa) compositions and (2) basalts and trachybasalts of KeNa series. In the diagram (a), a field of K series of moderate-Mg lavas from the Ugii-Nur volcanic field, Central Mongolia is shown after Rasskazov et al. (2012). In the diagram (b), the dividing lines of the series in terms of the potassium contents are used after Rollinson (1993). Fields in the diagram (c): P e picrite, PB e picrobasalt, BSN, TPH e basanite, tephrite, PHT e phonotephrite, B e , TB e trachybasalt, BA e basaltic andesite, BTA e basaltic trachyandesite, TA e trachyandesite, F e foidite. The dividing lines in (c) are after Le Bas and Streckeisen (1991). Contents of major oxides are recalculated to 100% without loss on ignition.

Please cite this article in press as: Chuvashova, I., et al., Mid-Miocene thermal impact on the lithosphere by sub-lithospheric convective mantle material: Transition from high- to moderate-Mg magmatism beneath Vitim Plateau, Siberia, Geoscience Frontiers (2016), http://dx.doi.org/ 10.1016/j.gsf.2016.05.011 I. Chuvashova et al. / Geoscience Frontiers xxx (2016) 1e22 7

Table 1 Major oxides (wt.%), trace elements (ppm), and Sr, Nd isotope ratios in representative compositions of the rock groups from the Bereya volcanic center.

Age (Ma) 16e14 14e13

Group Contaminated rocks, High-Mg basanitesebasalts Picrobasaltsebasalts, Basaltsebasaltic Basaltsetrachybasalts, KeNa series of transitional (KeNaeK) K series andesites of transitional KeNa series composition (NaeKeNa) composition

Sample VT-13-9 VT-13-9 gm VT-13-20 VT-13-19 V9 V18

SiO2 56.20 47.90 42.55 42.51 52.89 50.71 TiO2 1.13 1.89 1.85 2.09 1.85 2.28 Al2O3 15.30 13.03 9.79 10.27 14.36 16.13 Fe2O3 3.95 5.84 4.20 5.46 11.18* 10.75* FeO 2.96 5.29 7.54 6.61 N.d. N.d. MnO 0.11 0.16 0.16 0.15 0.18 0.17 MgO 5.52 10.12 17.23 13.87 7.23 5.61 CaO 5.22 7.11 8.06 8.99 8.87 8.69

Na2O 3.80 2.90 2.22 1.28 3.08 3.88 K2O 3.31 2.11 1.60 1.61 0.63 1.70 P2O5 0.36 0.62 0.64 0.73 0.23 0.48 LOI 2.04 2.98 4.07 6.61 0.00 0.00 Total 99.90 99.95 99.91 100.18 100.50 100.40 Sc 15.9 22.2 12.3 11.5 22.7 19.0 V 105 163 147 156 169 183 Co 30 47 63 57 44 34 Cu 5 35 40 35 55 40 Ni 214 326 744 458 131 54 Cr 238 550 1301 651 193 85 Rb 73 45 26 22 9 24 Sr 632 769 791 4042 389 713 Y 22.4 26.1 18.9 20.9 23.0 24.0 Zr 88 156 201 227 123 198 Nb 21.9 38.2 49.1 57.0 14.5 32.7 Cs 1.26 0.62 0.19 0.25 N.d. N.d. Ba 1103 845 651 1048 194 470 La 48.3 56.7 44.4 54.3 10.8 25.8 Ce 95.6 114.5 89.6 107.0 23.4 52.1 Pr 10.9 13.1 10.9 12.7 3.2 6.6 Nd 42.2 50.9 41.3 49.4 14.2 26.5 Sm 7.23 9.87 8.64 9.71 4.08 6.16 Eu 2.37 2.92 2.85 3.34 1.52 2.09 Gd 6.12 9.68 6.59 7.71 3.89 5.83 Tb 0.77 1.24 0.98 1.13 0.73 0.87 Dy 4.69 5.73 4.54 5.08 4.15 4.59 Ho 0.91 0.97 0.72 0.79 0.78 0.84 Er 2.26 2.35 1.56 1.68 2.00 2.08 Yb 1.72 1.80 0.86 0.97 1.57 1.63 Lu 0.21 0.26 0.09 0.11 0.22 0.22 Hf 1.96 4.33 4.70 5.47 3.12 4.47 Ta 1.36 2.33 3.18 3.62 0.88 2.03 Pb 12.4 9.7 4.5 5.1 N.d. N.d. Th 5.39 5.77 5.51 6.51 1.18 2.88 U 0.81 1.11 1.17 2.44 0.32 0.81 87Sr/86Sr 0.706556 N.d. 0.704069 0.704243 0.70374 0.70431 2s 0.000036 0.000050 0.000075 0.00001 0.00001 144Nd/143Nd 0.512444 N.d. 0.512852 0.512809 0.512934 0.512788 2s 0.000049 0.000024 0.000044 0.000009 0.000009

The group of high-Mg contaminated rocks is exhibited with bulk and groundmass compositions VT-13-9 and VT-13-9 gm, respectively. *Defined only total iron content in the form of Fe2O3. N.d. e no data. Data on samples from the 16e14 Ma unit were obtained by authors, on those from the 14e13 Ma unit by Harris (1998). silica-rich sample VT-13-09 (SiO2 57.6 wt.%) falls into the tra- Mg basaltsetrachybasalts contain 10.3e12.6 wt.% MgO. The chyandesite field. We assign this rock to a trachyandesite-like 14e13 Ma lavas display a dramatic decrease of MgO contents to composition. While compared to bulk compositions, ground- 6e10 wt.%. masses separated from the high-Mg basaltsetrachybasalts reveal The rock groups of 16e14 and 14e13 Ma are clearly seen systematically lower SiO2 abundances in the range of in diagrams of normative minerals (Fig. 7). High- 45.9e49.6 wt.%. Lavas of 14e13 Ma have the whole range of SiO2 Mg basaltsetrachybasalts differ from high-Mg basanitesebasalts from 47.3 to 53.5 wt.%. KeNa basaltsetrachybasalts differ from the and picrobasaltsebasalts by lower contents of normative anorthite transitional (NaeKeNa) basaltsebasaltic andesites by lower SiO2 100 An/(An þ Ab). High-Mg basanitesebasalts show 3e9% Ne, contents. whereas picrobasaltsebasalts 0e8% Hy. In terms of these petro- MgO contents in high-Mg basanitesebasalts and pic- chemical characteristics, picrobasaltsebasalts are similar to absar- robasaltsebasalts range from 12.5 to 20.4 wt.%. With increasing okites from Southwest Japan (Tatsumi and Koyaguchi, 1989). The SiO2 contents in the bulk compositions of the high-Mg basaltetrachybasalt, Na basalt, and basaltebasaltic andesite groups basaltsetrachybasalts, MgO contents decrease slightly relative yield sub-parallel trends of slightly increasing normative anorthite to high-Mg basanitesebasalts and picrobasaltsebasalts, ranging consistent with increasing 100 An/(An þ Ab) and a general from 10.7 to 16.8 wt.%. The groundmass separates from the high- transition from Ae-normative to AtzeHy-normative compositions.

Please cite this article in press as: Chuvashova, I., et al., Mid-Miocene thermal impact on the lithosphere by sub-lithospheric convective mantle material: Transition from high- to moderate-Mg magmatism beneath Vitim Plateau, Siberia, Geoscience Frontiers (2016), http://dx.doi.org/ 10.1016/j.gsf.2016.05.011 8 I. Chuvashova et al. / Geoscience Frontiers xxx (2016) 1e22

Figure 7. Diagrams of CIPW normative Nee(Hy þ Qtz) vs. 100 An/(An þ Ab) (a) and Nee(Hy þ Qtz) vs. MgO (b). Symbols are as in Fig. 6.

Simultaneous reduction of MgO concentrations and normative incompatible elements from Ba to Sr, although peaks and troughs of anorthite from 16e14 to 14e13 Ma rocks demonstrates decreasing these two samples are consistent with those in moderate-Mg lavas. crystallization temperatures of magmatic liquids. Basalts of Na series reveal specific minimum of K relative to U and Nb. The trachyandesite-like rock VT-13-9 indicates an effect of crustal contamination by elevated Cs, Rb, K, Pb and reduced Nb, Ta, 5. Trace elements Sr, P, Zr, Hf, and Ti abundances.

The diagrams of rare earth elements (REE) show distinct separa- 6. Discussion tion of rocks from the 16e14 and 14e13 Ma units (Fig. 8). Light REE (LREE) decrease from the high-Mg basanite VT-13-14 and picrobasalt For understanding magmatic processes that led to successive VT-13-19 to moderate-Mg rocks, whereas heavy REE (HREE) increase eruptions of high- and moderate-Mg liquids, we associate the accordingly. As a result, REE patterns for high- and moderate-Mg mid-Miocene appearance of high-Mg with high temper- rocks intersect each other. Chondrite-normalized La/Yb ratios range atures in the Transbaikalian melting domain, define geochemical from 22 to 38 in the groups of high-Mg basanitesebasalts and pic- signatures of melts erupted from deep and shallow mantle sour- robasaltsebasalts and from 5 to 18 in the groups of moderate-Mg ces and, finally, present a model of a thermal impact of a hot sub- basalts. The trachyandesite-like rock VT-13-9 is comparable to the lithospheric material on the lithosphere beneath the Vitim high-Mg basanite VT-13-14 in terms of LREE and to moderate-Mg Plateau. lavas in terms of HREE. Four samples of the high- Mg basaltsetrachybasalts demonstrate elevated concentrations of REE in groundmass relative to bulk compositions (not shown in 6.1. Crustal contamination diagrams). The rock grouping is supported by a multi-element diagram. High-Mg volcanic rocks of the initial volcanic episode 16e14 Ma Compared to moderate-Mg lavas, the high-Mg basanite VT-13-14 are sufficiently contaminated by crustal material. In diagrams, data and picrobasalt VT-13-19 show elevated concentrations of points of high-Mg contaminated basaltsetrachybasalts and the

Please cite this article in press as: Chuvashova, I., et al., Mid-Miocene thermal impact on the lithosphere by sub-lithospheric convective mantle material: Transition from high- to moderate-Mg magmatism beneath Vitim Plateau, Siberia, Geoscience Frontiers (2016), http://dx.doi.org/ 10.1016/j.gsf.2016.05.011 I. Chuvashova et al. / Geoscience Frontiers xxx (2016) 1e22 9

Figure 8. Pyrolite-normalized concentrations of incompatible elements (a) and chondrite-normalized contents of REE (b). The pyrolite and chondrite compositions are after McDonough and Sun (1995). trachyandesite-like rock VT-13-9 are offset from the trends of high- moderate-Mg basaltsetrachybasalts to the basaltic trachyandesite Mg rocks. The contamination is expressed in increasing SiO2 and V-9. The shift is resulted from the generally increasing degree of decreasing 100 An/(An þ Ab) to values typical of moderate-Mg partial melting in the mantle sources. Meanwhile, Sr isotope ratios rocks (Fig. 6). may be overstated due to rock alteration. Although this effect can In the diagram Th/Yb vs. Ta/Yb, high-Mg contaminated be reduced by acid treatment of a sample (Song et al., 1990; basaltsetrachybasalts and trachyandesite-like rock VT-13-9 are Rasskazov et al., 2015), no criterion for a strontium isotope ratio shifted above the mantle array. This shift is a characteristic of in an unaltered rock is substantiated and no true value of this data points not only for bulk compositions of rocks, but also for parameter can be estimated. groundmasses (Fig. 9a). Thus, crustal material was partially melted and introduced into crystallized liquids. Bulk composition 6.2. Trends of high- and moderate-Mg rocks of the trachyandesite-like rock VT-13-9 reveals lower Ta/Yb, while its groundmass is comparable with those of other On plots of compatible (Ni) and incompatible (Th, La, K2O, Rb, contaminated species. This sample is the most contaminated by Sr, Ba) trace elements versus MgO, high- and moderate-Mg rocks crustal material. yield individual groups (Fig. 9aeg). On the Ni vs. MgO plot, data Crustal contamination of high-Mg basaltsetrachybasalts and points of bulk compositions for contaminated rocks show the trachyandesite-like rock VT-13-9 is reflected in high initial elevated Ni and MgO contents, as compared with groundmass 87 86 87 86 Sr/ Sr (Fig. 9b). In the diagram ( Sr/ Sr)i vs.1000/Sr, these rocks compositions, except strongly contaminated trachyandesites-like demonstrate a typical effect of crustal contamination (i.e. rock VT-13-9. Basaltsetrachybasalts from the lower 16e14 Ma 87 86 increasing both ( Sr/ Sr)i and 1000/Sr) relative to uncontami- unit resulted from partial assimilation of crust material accom- nated 16e14 Ma rocks. Various groups of uncontaminated rocks are panied with fractional crystallization of olivine. From shifted along the ordinate with increasing 1000/Sr from pic- an assimilationefractionalecrystallization (AFC) model (DePaolo, robasaltsebasalts through high-Mg basanitesebasalts and 1981) on the Ni vs. MgO plot (Fig. 10a), groundmasses of these

Please cite this article in press as: Chuvashova, I., et al., Mid-Miocene thermal impact on the lithosphere by sub-lithospheric convective mantle material: Transition from high- to moderate-Mg magmatism beneath Vitim Plateau, Siberia, Geoscience Frontiers (2016), http://dx.doi.org/ 10.1016/j.gsf.2016.05.011 10 I. Chuvashova et al. / Geoscience Frontiers xxx (2016) 1e22

Figure 9. Binary diagrams of Th/Yb vs. Ta/Yb (a) and initial 87Sr/86Sr vs. 1000/Sr (b). Due to young ages of rocks and low Rb/Sr, the initial Sr isotopic ratios do not differ from the measured values. E-MORB and OIB are from McDonough and Sun (1995). Symbols are as in Fig. 6.

rocks underwent fractional crystallization F ¼ 9e13%. For pic- region, whereas the trend of moderate-Mg melts results from robasalts this value decreases to a range of 3e10% and for melting in a shallow mantle source. Increasing MgO consistent with basanites to a range of 0e4%. decreasing concentrations of incompatible trace elements in the The AFC model curve assumes a possible origin of the mid- groups of high-Mg basanitesebasalts and picrobasaltsebasalts in- Miocene rocks from a common primary high-Mg melt at the pro- dicates increasing temperature and degree of partial melting. In the portion of initialecontaminated compositions 1:4. This proportion group of high-Mg contaminated basaltsetrachybasalts, the highest is adopted from concentrations of MgO (16.98 wt.%), Ni (756 ppm) concentrations of incompatible elements are defined in a ground- in the sample VT-13-8 and from concentrations of MgO (1.97 wt.%), mass of the basalt VT-13-21 (Th ¼ 7.3 ppm, La ¼ 68 ppm, Ni (78.5 ppm) in the sample VT-13-9, respectively. From trace- Sr ¼ 4912 ppm, Ba ¼ 1508 ppm). Conodes of these elements in bulk element modeling, a source is calculated to contain Ol compositions and groundmasses have steeper slopes than the (57%), Opx (25%), Cpx (5%), Grt (6%), and Amph (5%) (see details in trend of other high-Mg rocks. The conodes of uncontaminated Section 6.3). Concentrations of MgO and Ni in minerals and min- high-Mg basaltic samples lay sub-parallel to those of the basalt VT- eral/melt distribution coefficients are presented in Tables 2 and 3, 13-21 on the Th vs. MgO and La vs. MgO plots and more gentle on respectively. Bulk distribution coefficients Dmineral/melt(Ni) ¼ 7.39, the Sr vs. MgO and Ba vs. MgO ones. Minimal concentrations of La, Dmineral/melt(Mg) ¼ 4.48 are calculated from the mineral proportions Sr, and Ba in bulk compositions and groundmasses of high-Mg in the source with Ni abundance of 1750 ppm. The melting line of basalts overlap each other. As MgO decreases in the with the Ni concentration of 1900 ppm (Sobolev et al., trachyandesite-like rock VT-13-9 from the bulk composition to the 2005) demonstrates compositions with higher MgO contents. De- groundmass, La and Th increase and Ba decreases with minor viations of MgO from 20.4 to 12.5 wt.% in high-Mg rocks appear to variations of Sr (Fig. 10b,c,f and g). Unlike Th, La, Sr, and Ba, con- reflect relative variations mineral proportions in a source region centrations of K and Rb depend substantially on crustal contami- (see Fig. 13a in Section 6.3). nation. The group of high-Mg contaminated basaltsetrachybasalts Assumption on the origin of both moderate and high-Mg rocks and trachyandesite-like rock VT-13-9 show enrichment by both from a common primary high-Mg melt contradicts the distribution elements relative to other high-Mg groups (Fig. 10d and e). of incompatible elements. We suggest that the trend of high-Mg Grouping of volcanic rocks is expressed also in trace-element rocks reflects temperature variations in a sub-lithospheric source ratios. Zr/Nb is lower in high-Mg rocks than in moderate-Mg ones

Please cite this article in press as: Chuvashova, I., et al., Mid-Miocene thermal impact on the lithosphere by sub-lithospheric convective mantle material: Transition from high- to moderate-Mg magmatism beneath Vitim Plateau, Siberia, Geoscience Frontiers (2016), http://dx.doi.org/ 10.1016/j.gsf.2016.05.011 I. Chuvashova et al. / Geoscience Frontiers xxx (2016) 1e22 11

Figure 10. Ni, Th, La, K2O, Rb, Sr, and Ba vs. MgO (aeg) and Zr/Y vs. Zr/Nb (h). Symbols are as in Fig. 6.

Please cite this article in press as: Chuvashova, I., et al., Mid-Miocene thermal impact on the lithosphere by sub-lithospheric convective mantle material: Transition from high- to moderate-Mg magmatism beneath Vitim Plateau, Siberia, Geoscience Frontiers (2016), http://dx.doi.org/ 10.1016/j.gsf.2016.05.011 12 I. Chuvashova et al. / Geoscience Frontiers xxx (2016) 1e22

(2.3e4.3 and 5.0e12.2, respectively). Zr/Y is higher in pic- Table 3 robasaltsebasalts than in high-Mg basanitesebasalts. In the dia- The mineral/basaltic melt distribution coefficients used for the partial melting and AFC modeling. gram Zr/Y vs. Zr/Nb, the high-Mg groups of basalts-trachybasalts, basanitesebasalts, and the moderate-Mg basalt group of KeNa Element Ol Opx Cpx Grt Phl Amph Ap series show correlated trends of these ratios. Those of moderate- Rb 0.00002 0.0001 0.0004 0.00002 5.18 0.2 0.01 Mg transitional (NaeKeNa) basaltsebasaltic andesites demon- Ba 0.000005 0.00006 0.0007 0.00007 3.48 0.16 1.555 strate sufficient variations of Zr/Nb and limited ranges of Zr/Y Th 0.000007 0.0002 0.013 0.0021 0.0014 0.0039 2.131 U 0.000009 0.00004 0.006 0.011 0.0011 0.0041 2.131* (Fig. 10h). K 0.00002 0.0001 0.07 0.013 3.67 0.58 N.d. Nb 0.00005 0.0013 0.0027 0.02 0.0853 0.159 0.001 Ta 0.00005 0.0025 0.012 0.05 0.2255 0.16 0.001* 6.3. Components of mantle sources La 0.0002 0.003 0.054 0.0007 0.0004 0.055 7.39 Ce 0.00007 0.0021 0.086 0.0026 0.0006 0.096 5.965 In terms of Nd and Sr isotope ratios, sources of the 16e14 Ma Pr 0.0003 0.0022 0.154 0.01 0.0009 0.17 4.64* Sr 0.00004 0.0015 0.096 0.0007 0.183 0.298 3.055 high-Mg magmas have a common depleted isotopic component Zr 0.001 0.012 0.123 0.2 0.0232 0.18 0.012 (depleted mantle, DM) and other less depleted compositions. This Hf 0.0029 0.019 0.256 0.4 0.048 0.29 0.01 is indicative for convective isotopic homogenization of a sub- Sm 0.0009 0.0037 0.291 0.1 0.0008 0.33 4.603 lithospheric mantle material and its fractional differentiation, Ti 0.015 0.12 0.384 0.32 1.768 1.29 N.d. accompanied by incompatible element enrichment of resulting Y 0.0082 0.02 0.467 2.5 0.007 0.52 3.43 Yb 0.024 0.032 0.43 6.4 0.023 0.47 0.93 melts (Figs. 11 and 12). The separated melts, therefore, are derived Ni 7 3.1 1.4 0.19 5 3.4 N.d. due to high-temperature sub-lithospheric convective processes in Mg 5 4.7 5.4 N.d. N.d. 3.8 N.d. the mantle and decreasing melting temperatures connected with Data sources: olivine (Ol): Ni, Mg (Wang and Gaetani, 2008), other elements increasing incompatible element abundances in a source. This type (Halliday et al., 1995); orthopyroxene (Opx): Ta (Green et al., 2000), Ti (Girnis et al., of magma generation differs from the generation of komatiite 2006), Yb (Kennedy et al., 1993), Ni (O’Reilly et al., 1991), Mg (Henderson, 1982), magmas in a plume that shows depletion in both isotopic and other elements (Halliday et al., 1995); clinopyroxene (Cpx): Rb (Halliday et al., elemental signatures (Arndt et al., 2008). 1995), Ba, K, REE, Zr, Hf, Ti, Y (Hart and Dunn, 1993), Th and U (Hauri et al., 1994), Nb and Sr (Foley et al., 1996), Ta (Green et al., 1989), Ni (O’Reilly et al., 1991), Mg High-Mg liquids from the Vitim Plateau are enriched by LREE, (Villemant et al., 1981); garnet (Grt): Rb, Ca, Th, U, K, REE, Sr, Y (Halliday et al., 1995), Th, Zr, and other incompatible elements that are introduced into Nb, Ta, Zr, Hf, Ti (Girnis et al., 2006), Ni (O’Reilly, Griffin, 1995); phlogopite (Phl): Rb, melts due to prevailing partitioning from various mineral phases. Ba, Nb, Sr, Zr, Y (Foley et al., 1996), Th, U, K, Ti (LaTourrette et al., 1995), Ta, La ’ Peridotite xenoliths, extracted from the spinelegarnet transition (Gregoire et al., 2000), Ce, Pr, Hf, Sm, Yb (Ionov et al., 1997), Ni (O Reilly et al., 1991); amph (Amph): Ta (Zack et al., 1997), Zr, Hf, Sm (Ionov et al., 1997), Ni (O’Reilly et al., zone by ascending liquids of this area, show evidence on melt 1991), Mg (Henderson, 1982), other elements (LaTourrette et al., 1995); apatite infiltration in the upper mantle (Glaser et al.,1999). Various types of (Ap): Rb, Nb, Ta, Hf (Ionov et al., 1997), other elements (Chazot et al., 1996). *Values mantle xenoliths have mutual genetic relationships. It was sug- after interpolation. N.d. e no data. gested that the high-temperature garnet-bearing rocks result from repeated affection of basaltic melts on the ambient upper mantle material. From the AFC calculations, it was also inferred that metasomatized lherzolites and are examined as megacrysts and shallow xenoliths crystallized during the rise of possible sources for the mid-Miocene liquids. basaltic liquids to the surface (Ashchepkov and Andre, 2002; The group of high-Mg basanitesebasalts has initial Ashchepkov et al., 2011). From these observations, 143 144 87 86 ( Nd/ Nd)t ¼ 0.512867e0.512892 and ( Sr/ Sr)i ¼ 0.703831 e0.704183 in the range of the depleted mantle. The group of pic- robasaltsebasalts reveals less depleted isotopic signatures: initial 143 144 87 86 Table 2 ( Nd/ Nd)t ¼ 0.512730e0.512842, ( Sr/ Sr)i ¼ 0.703778 Trace element concentrations (ppm) in the pyrolite and mean concentration in e0.704406 with a relative shift of data points along the mantle array minerals from xenoliths used for calculations of the model magmatic sources towards the primitive mantle composition. Crustal contamination of beneath Vitim Plateau. high-Mg basaltsetrachybasalts and the trachyandesite-like rock VT- Element Pyrolite Cpx Grt Phl Amph Ap 87 86 143 144 13-9 is reflected both in high ( Sr/ Sr)i and low ( Nd/ Nd)t.The Rb 0.6 0.45 N.d. 263 6 15 group of moderate-Mg basaltsetrachybasalts overlaps both the high- Ba 6.6 11 0.03 3364 178 187 Mg basaniteebasalt and picrobasaltebasalt compositions. The Th 0.0795 0.07 N.d. 0.005 1.50 97 U 0.0203 0.01 N.d. N.d. 0.50 23 basaltic trachyandesite V-9 shows a more depleted signature: 143 144 87 86 K 240 83 N.d. 73390 15994 N.d. ( Nd/ Nd)t ¼ 0.512934, ( Sr/ Sr)i ¼ 0.703740 (Fig. 11a). Nb 0.658 0.79 0.14 19.1 42.8 0.66 The differences between the eruption age of 16e14 Ma and Ta 0.037 0.028 N.d. 0.38 3.23 o. e. SmeNd ages of 42e24 Ma obtained for garneteclinopyroxene pairs La 0.648 2.33 0.01 0.69 6.01 2300* of xenoliths from tuffaceous-sedimentary strata of the Bereya vol- Ce 1.675 9.42 0.07 0.9 21.3 4042 Pr 0.254 2* 0.02 0.27 4* 377 canic center was explained by incomplete isotopic equilibrium in Sr 19.9 104 0.18 212 467 13295 the xenoliths at the time of eruption due to insufficiently rapid Zr 10.5 95.5 24.3 18 119 24 diffusional exchange in coarse-grained rocks. SmeNd and RbeSr Hf 0.283 2.76 0.24 0.31 4.94 0.71* isotope data on xenoliths have been interpreted in terms of a Sm 0.406 3.25 0.48 0.14 5.5 170* w Ti 1260 7008 1112 47098 25120 N.d. 2 Gyr depletion event relative to the primitive mantle values after Y 4.3 8.29 27.7 0.78 11.4 173 Zindler and Hart (1986). It was suggested that the mantle perido- Yb 0.441 0.40 2.92 0.06 0.52 7.67 tites have been attached to the continental lithosphere since that Ni 1960 468 47 N.d. 868 N.d. time and survived major tectono-thermal events (Ionov et al., MgO (wt.%) 37.80 14.68 20.72 19.40 16.76 N.d. 2005). In the area of study, no geological units of w2 Gyr have Data sources: Ionov et al., 1997; Glaser et al., 1999; Litasov et al., 2000a,b; O’Reilly been recorded, though. fi and Grif n, 2000; Ashchepkov and Andre, 2002; Ashchepkov et al., 2003, 2011. From Pb isotope data on Cenozoic volcanic rocks, it was inferred Th and U in Amph is after Zheng et al. (2006). The pyrolite composition is after e McDonough and Sun (1995). *The values obtained after interpolation from that the mantle of the Baikal Mongolian region was isotopically normalized spectra. N.d. e no data. homogenized by vigorous convection at about 700e600 Ma in

Please cite this article in press as: Chuvashova, I., et al., Mid-Miocene thermal impact on the lithosphere by sub-lithospheric convective mantle material: Transition from high- to moderate-Mg magmatism beneath Vitim Plateau, Siberia, Geoscience Frontiers (2016), http://dx.doi.org/ 10.1016/j.gsf.2016.05.011 I. Chuvashova et al. / Geoscience Frontiers xxx (2016) 1e22 13

Figure 11. Initial 143Nd/144Nd vs. initial 87Sr/86Sr (a) and 143Nd/144Nd vs. 147Sm/144Nd (c). Insets (b) and (d) show data on xenoliths from the 16e14 Ma high-Mg rocks (Ionov et al., 2005). Symbols are as in Fig. 6. connection with disintegration of the supercontinent Rodinia Plateau. On Fig. 11c and d, data points of high-Mg basanites and (Chuvashova, 2013). This event seems to be displayed also in Nd mantle xenoliths with extremely high 143Nd/144Nd ratio are isotopic composition of the lithospheric mantle beneath the Vitim distributed along a line with a slope of 640 Ma. One can assume that the material of lithospheric mantle xenoliths is a comple- mentary product of mantle depletion of convective material responsible for disintegration of the supercontinent Rodinia, while the high-Mg basanites show time-integrated 143Nd/144Nd ratio of sub-lithospheric convective mantle. Offset of individual data points of xenoliths with increasing 147Sm/144Nd ratio might reflect minor transformation of the lithospheric mantle, associated with the formation of the island-arc Eravna terrane at about 440 Ma. Evolution of the convective mantle material at shallow depths is reflected in shifting data points parallel to abscissa-axis with increasing 147Sm/144Nd. Under generally depleted isotopic signa- tures of Nd in garnet and spinel from the 16e14 Ma tuffaceous-sedimentary stratum, the only composite xenolith of garnet-bearing xenolith 313-4 falls into a field of volcanic rocks. In terms of isotopic compositions of Nd, volcanic rocks are comparable also to minerals from veins of composite xenoliths. Consequently, the shallow-mantle volcanic rocks cannot be derived from a source of the depleted lithospheric mantle and might be related to the processes, exhibited by the formation of veins in mantle rocks. Figure 12. Initial 143Nd/144Nd vs. 100/Nd. Symbols are as in Fig. 6. The reported data on veins (Ionov et al., 2005) exhibit hetero- geneities of the lithospheric mantle on the scale of individual

Please cite this article in press as: Chuvashova, I., et al., Mid-Miocene thermal impact on the lithosphere by sub-lithospheric convective mantle material: Transition from high- to moderate-Mg magmatism beneath Vitim Plateau, Siberia, Geoscience Frontiers (2016), http://dx.doi.org/ 10.1016/j.gsf.2016.05.011 14 I. Chuvashova et al. / Geoscience Frontiers xxx (2016) 1e22

Figure 13. (a) (La/Yb)N vs. (Yb)N and (b) La/Yb vs. MgO diagrams. For simplification the diagram (a), data points of bulk and groundmass compositions for only two samples of the contaminated basalts that yield the highest (La/Yb)N are plotted. Abbreviations: Grt e garnet, Cpx e clinopyroxene, Amph e amphibole, Phl e phlogopite, Ap e apatite. Ratios of minerals in sources are based on trace-element modeling (see the text). Mineral/melt distribution coefficients for La and Yb are shown in Table 3. The primitive mantle composition is adopted after McDonough and Sun (1995). In diagram (b), data obtained are compared with data fields of Miocene lavas from the Hannuoba Province of China: the groups of evolved alkali olivine basalts (EVeAOB), primitive alkali olivine basalts (PReAOB), transitional basalts (TRB), and tholeiitic basalts (TLB) after Zhi et al. (1990).

Please cite this article in press as: Chuvashova, I., et al., Mid-Miocene thermal impact on the lithosphere by sub-lithospheric convective mantle material: Transition from high- to moderate-Mg magmatism beneath Vitim Plateau, Siberia, Geoscience Frontiers (2016), http://dx.doi.org/ 10.1016/j.gsf.2016.05.011 I. Chuvashova et al. / Geoscience Frontiers xxx (2016) 1e22 15 samples. On the one hand, the isotopic differences that might exist mantle source, which expresses the mantle temperature, projected among mineral grains of chemical layering on centimeter or meter along the adiabatic line to the atmospheric pressure at the Earth’s presumably, will not be seen except in xenolithic studies. On the surface. This parameter can vary from 1250 to 1400 C in mid-ocean other hand, references to lithospheric mantle components imply ridges and rise up to 1550 C in hot plumes such as Hawaii. Change volumes of rock large enough to represent the source of individual to the conductive transfer of geothermal heat flow in the litho- batches of magma, perhaps of the dimension of several cubic ki- sphere causes a large temperature gradient, in which the temper- lometers. For instance, volume of 15 km3 moderately-Mg magmas, ature increases with depth for 500 K GPa 1 (16 C km 1)or erupted in the Bereya volcanic center at 14e13 Ma, assumes a 20% more. As a result, melting of the mantle is dominated by decom- partial melting and segregation of magmas from a mantle batch as pression mechanisms that cause eruptions of basaltic magma large as 75 km3. through the extended lithosphere (Asimow, 2005; Herzberg et al., Mantle components of volcanic rocks are identified by cooper- 2007). 143 144 # ative use of the diagrams ( Nd/ Nd)t vs. 100/Nd and (La/Yb)N vs. Rocks enriched in Ni and Cr and showing Mg ¼ Mg/ 2þ 3þ YbN (Figs. 12 and 13). The high-Mg basanite 93VBS-4 with the (Mg þ Fe )>70 with corrections Fe ¼ 0.15 Fetotal are usually 1 143 144 lowest 100/Nd (1.7 ppm ) and relatively high ( Nd/ Nd)t related to the primary mantle liquids. Besides, the composition of (0.512892, εNd ¼þ5) is the starting point for trends of high-Mg primary melt should be in equilibrium with the most magnesian basanitesebasalts (BSN e high Grt), picrobasaltsebasalts (PB), olivine (Herzberg et al., 2007). High-Mg volcanic rocks from the and high-Mg basanites with increasing degree of partial melting Bereya volcanic center are characterized by Ni and Cr as high as 788 (BSN e high F). This composition is designated as the common and 1503 ppm, respectively, and by interval of Mg# ¼ 66.3e77.5 component of the depleted mantle (DM). The trend of high-Mg (predominated values 75e76). Some Mg# values decrease to 64 143 144 contaminated rocks yields low ( Nd/ Nd)t and 100/Nd because of crustal contamination and increase to 78 due to variable (elevated Nd concentrations). The data points of the 14e13 Ma mineral proportions in the mantle source (see Figs. 10a and 13a). moderate-Mg lavas offset from compositions of the 16e14 Ma high- Unlike high-Mg rocks with typical signatures of primary mantle Mg rocks with elevated 100/Nd (reduced Nd concentrations). The liquids, moderate-Mg rocks show relatively low Cr and Ni con- common trend of moderately-Mg lavas and minerals from veins centrations and interval of Mg# ¼ 53.9e63.2 that demonstrate (Fig. 12) is indicative for origin of the lavas through melting of the generation of such magmas from a modified mantle source. vein-bearing shallow mantle source. A maximum liquidus temperature of olivine crystallization from In the Bereya volcanic center, high-Mg liquids were followed an erupted melt or primary eruption temperature (Tpe) is estimated with moderate-Mg ones within a time interval of 1e2 million years. from an assumption on anhydrous composition of magma by There was drastic transition from high-Mg (high temperature) to various equations. Three of these by Kutolin (1966), Herzberg et al. moderate-Mg (moderate temperature) melts. In Miocene lavas (2007), and Arndt et al. (2008): Tpe ( C) ¼ 1056.6 þ 17.34 MgO, Tpe 2 from the Hannuoba province (North China), a similar contrast ( C) ¼ 935 þ 33 MgO e 0.37 MgO , Tpe ( C) ¼ 20 þ 1000 MgO existed between the groups of evolved alkali olivine basalts with yield approximately comparable results. For the high-Mg basanite high La/Yb and three other groups of basalts (tholeiitic, transitional VT-13-2 (MgO ¼ 17.1 wt.%), that belong to the predominated group and primitive alkaline) with low La/Yb. In the La/Yb vs. MgO dia- with Mg# ¼ 75e76, the maximum liquidus crystallization tem- gram, the three groups of Hannuoba basalts are compared with the perature obtained is in the range of 1340e1390 C. Calculation of a moderate-Mg lavas from the Bereya volcanic center (Zhi et al.,1990) potential temperature using equation Tp ( C) ¼ 1463 þ 12.7 MgO (Fig. 13b). e 2924/MgO (Herzberg et al., 2007) provides more satisfactory The group of evolved alkali olivine basalts from the Hannuoba approximation to the temperature conditions in the mantle. The province is shifted relative to the high-Mg rocks from the Bereya high-Mg basanite VT-13-2 yields an estimate Tp ¼ 1510 C. center with decreasing MgO content. The Hannuoba and Bereya The assessment of potential temperature for the mantle beneath high La/Yb trends are parallel to each other. They both can be the Vitim Plateau was at its maximum in the outline of the Trans- interpreted in connection with increasing degree of partial melting baikalian melting domain. The Tp for generation of high-Mg olivine at elevated temperatures. From model calculations by Chuvashova melaleucetites from the Udokan volcanic field does not exceed et al. (2012), it was inferred that the Hannuoba evolved alkali 1430 C. This estimate is comparable with the potential mantle olivine basalts were derived by 2e5% partial melting of a source temperature for rocks from Ontong Java Plateau and Iceland (early that had 4.5e6.0% of garnet. These parameters are close to those of Tertiary). The modern rocks from Iceland yield the lower Tp esti- a source for high-Mg rocks from the Bereya volcanic center. mates. Velocity models are indicative for association of Ontong Java Interestingly, that MgO contents strongly vary in the high-Mg Plateau and Iceland with upper mantle melting anomalies, not Bereya rocks only under elevated La/Yb and decrease with extending below the transition layer (Foulger, 2010). It is noteworthy decreasing La/Yb (i.e. at the elevated temperature and increasing that the calculated Tp ¼ 1550 C for Hawaii (Herzberg et al., 2007) degree of partial melting). The Hannuoba evolved alkali olivine exceeds the upper mantle Tp interval and is comparable with the basalts, affected with olivine fractionation (Zhi et al., 1990), had higher Tp interval for high-Mg rocks from Gorgona Island (Fig. 14). primary liquids compositionally similar to high-Mg rocks from the These high temperatures might reflect raising hot plumes from the Bereya volcanic center and other localities of the Transbaikalian deeper mantle according to the model suggested by Morgan (1971). melting domain shown on Fig. 1. The unfractionated high-Mg liq- The maximum temperatures in the mantle can be constrained uids from the Bereya volcanic center might be overheated above a also from temperature calculations of mineral equilibrium in liquidus. mantle xenoliths. The highest temperatures were obtained for xe- noliths of pegmatoid orthopyroxenites and websterites with la- 6.4. Temperature estimates mellas in . These xenoliths were found in basanites of the Udokan and Kamar volcanic fields. The temperature interval of The mantle melts at relatively high temperature and low pres- 1350e1450 C was calculated for initial Ca solubility in orthopyr- sure. Convective heat transfer in the asthenosphere provides ver- oxene solid solutions (Rasskazov, 1985; Rasskazov and Chuvashova, tical temperature gradient reflected in a material flux without 2013). In the Udokan volcanic field, the high-temperature xenoliths additional or lost heat. Asthenospheric adiabatic gradient of about occur in moderate-Mg basanite lavas erupted along the axial part of 1 10 K GPa corresponds to the potential temperature (Tp)ofa Udokan Range in the last 4 Ma. These basanites (as well as their

Please cite this article in press as: Chuvashova, I., et al., Mid-Miocene thermal impact on the lithosphere by sub-lithospheric convective mantle material: Transition from high- to moderate-Mg magmatism beneath Vitim Plateau, Siberia, Geoscience Frontiers (2016), http://dx.doi.org/ 10.1016/j.gsf.2016.05.011 16 I. Chuvashova et al. / Geoscience Frontiers xxx (2016) 1e22 inclusions) were separated in time and space from the 14 Ma high- sources for volcanic rocks in the SayaneMongolian melting domain Mg melaleucitites located in the northern part of the volcanic field. is also confirmed by 3He/4He ratios in from lavas of the The potential temperature of the latter Tp ¼ 1430 C is consistent Hangay and Dzhida areas comparable with the values of this with the initial temperature in the mantle pegmatoid orthopyrox- parameter in basalts from mid-ocean ridges (Barry et al., 2007). enites. In the Kamar volcanic field, the high-temperature xenoliths were sampled in basanites from the Sukhoi volcano erupted at 6.5. Depth estimates 13e12 Ma. The calculated potential temperature Tr ¼ 1470 C for these basanites is close to the initial temperature of the mantle- LREE enrichment of magmas relative to HREE is controlled derived pegmatoid orthopyroxenites. principally by the presence of garnet in a source region. From In the outline of the SayaneMongolian melting domain, calcu- definition of Gd/Yb ¼ 0.15 in a garnet, it was inferred that its lations for a basalt with the highest MgO content from the Dzhida occurrence in a residual solid maintains a higher Gd/Yb ratio in a volcanic field (10.1 wt.%) yield crystallization temperature interval partial melt as compared to a source (Walter, 1998). The value Gd/ Tpe ¼ 1200e1230 C and estimate of potential temperature Yb ¼ 1.6 was suggested as the lower limit for a liquid from a garnet- Tp ¼ 1300 C. Melt inclusions in olivines from basalts of this melting bearing source (Sobolev et al., 2011). On the one hand, from the Gd/ domain also show relatively low homogenization temperatures: in Yb criterion, all the mid-Miocene melts, erupted in the Bereya Eastern Sayans 1220e1250 C, in the Dzhida volcanic field volcanic center, should be assigned to garnet-bearing sources. Gd/ 1175e1220 C, and in Hangay 1210e1300 C(Simonov et al., 2012). Yb ranges from 2.5 to 3.4 in basaltsebasaltic andesites and in- The highest value of the homogenization temperature is compa- creases to the intervals 4.9e5.7, 3.6e6.5, 5.3e8.1, and 7.9e9.3 in Na rable with a maximum estimate of potential temperature in the basalts, basaltsetrachybasalts, high-Mg basanitesebasalts, and SayaneMongolian domain. picrobasaltsebasalts, respectively. On the other hand, the groups of Relatively low temperatures in the mantle liquids of the basaltsebasaltic andesites demonstrate the narrow (La/Yb)N, which SayaneMongolian melting domain might reflect shallow depths of is 7e8 times higher than that in the pyrolite, and Yb concentrations the latter, thermally provided by asthenospheric convection in the varying from 1.38 to 1.95 ppm. These are characteristics of a source upper mantle no deeper than 200 km, unlike the Transbaikalian enriched by LREE with no garnet or low garnet content. melting domain that demonstrates decreasing velocities and Deep magma-generating processes in the sub-lithospheric possible relation of convective processes to the lower part of the mantle at 16e14 Ma are indicated by the potassic affinity of pic- upper mantle and transition layer. The evaluated potential tem- robasaltsebasalts. The high activity of potassium in a deep part of perature Tp ¼ 1300 C for the SayaneMongolian melting domain is the upper mantle is reflected in potassium-bearing clinopyroxenes within the temperature range of basalts from mid-ocean ridges in xenoliths from kimberlites and inclusions in diamonds, in the Tp ¼ 1280e1400 C. From seismic tomography, it was inferred that diamond-bearing pyroxeneegarnetecarbonaceous rocks of the low-velocity anomalies beneath mid-ocean ridges do not extend Kokchetav metamorphic complex and confirmed by experiments at deeper than 150 km (Anderson et al., 1992). The shallow nature of high pressures and temperatures (Sobolev, 1974; Tsuruta and Takahashi, 1998; Perchuk et al., 2002; Shatsky et al., 2006; Foley et al., 2009). Natural assemblages with potassium-bearing clino- pyroxenes are limited to the pressure range 5e7 GPa (Safonov et al., 2005) that corresponds to depths of 160e220 km. From the empirical equation P (kbar) ¼ 213.6e4.05 SiO2 (wt.%) (Scarrow and Cox, 1995), high-Mg basanitesebasalts and pic- robasaltsebasalts yield pressure ranges of 36.4e48.6 and 41.4e44.3 kbar, respectively. These pressures correspond to gen- erations of high-Mg melts beneath the Bereya volcanic center in the depth intervals 115e150 and 135e140 km during accumulation of the middle and upper layers of the 16e14 Ma stratum (Fig. 5). The maximum depth up to 150 km is characteristic of the sub-group of high-Mg basanitesebasalts with an elevated degree of melting (BSN e high F), which is the best estimate for Tp. A depth of lithospheric mantle sources, involved in melting about 14e13 Ma, is limited by a garnetespinel transition. Experi- mental data of the CaOeMgOeAl2O3eSiO2 system (CMAS) show this transition below the solidus in the pressure ranges 1.8e2.0 GPa at 1200 C and 2.6e2.7 GPa at 1500 C (at depth 60e85 km). In Cr- bearing systems, spinel coexists with garnet. Increasing concen-

Figure 14. Primary eruption and mantle potential temperatures as a function of the trations of Cr shift a reaction of garnet introduction to the higher 2þ MgO contents in primary magmas. Primary eruption temperatures are anhydrous pressures. The presence of Fe has the opposite effect (Klemme, olivine liquidus temperatures at 1 bar calculated for the primary magmas using the 2004; Klemme and O’Neill, 2000). Similar estimates for the gar- fi ¼ þ e 2 simpli ed equation Tpe ( C) 935 33 MgO 0.37 MgO . The bracket shows the netespinel transition have been obtained through experiments on uncertainty at the 2s level. Potential temperatures are calculated from the equation enriched and depleted peridotites. Under lower pressure condi- Tr ( C) ¼ 1463 þ 12.7 MgO e 2924/MgO (Herzberg et al., 2007). Arrows display the effects of H2O content on liquidus temperature. Partial crystallization of primary tions, melting is assumed to begin in garnet veins magmas will yield actual eruption temperatures that are lower than primary eruption (Robinson and Wood, 1998). temperatures. Abbreviations: MORB e mid-ocean ridge basalts, OJP e Ontong Java Plateau, I (ET) e Island (early Tertiary), I (PD) e Island (present day). Highlighted in 6.6. Partial melting modeling three temperature ranges are: (1) Tr ¼ 1550 C and above e plumes rising from the deeper parts of the mantle (Gorgona, Hawaii), (2) Tr ¼ 1430e1510 C e upper-mantle melting anomalies extending from the mantle transition zone (Ontong Java Plateau, For modeling partial melting of enriched sources, we use mean Iceland, Transbaikalian domain), and (3) Tr ¼ 1280e1400 C e melting anomalies in compositions of minerals from Grt- and GrtePhl lherzolites and the shallow part of the upper mantle (mid-ocean ridges, SayaneMongolian domain). also pyroxenites with phlogopite, ilmenite, and kaersutite, as well

Please cite this article in press as: Chuvashova, I., et al., Mid-Miocene thermal impact on the lithosphere by sub-lithospheric convective mantle material: Transition from high- to moderate-Mg magmatism beneath Vitim Plateau, Siberia, Geoscience Frontiers (2016), http://dx.doi.org/ 10.1016/j.gsf.2016.05.011 I. Chuvashova et al. / Geoscience Frontiers xxx (2016) 1e22 17 as other deep-seated xenoliths that satisfy to the conditions: (1) a content of 6% towards melt fractions F up to 0.075. The maximal mineral assemblage Ol þ Cpx þ Grt þ Amph Phl Ap represents depth estimates up to 150 km from these data points are consistent a possible model source, (2) PeT parameters of xenoliths corre- with increasing MgO content due to high-temperature magma spond to the upper mantle conditions (i.e. T>1000 C), and (3) generation. A trend for picrobasaltsebasalts (PB) crosses the middle analytical data for trace-element concentrations in minerals are of part of the BSN e high Grt one. Unlike the latter source, garnet high quality. The trace element concentrations were defined in content of the PB trend remains constant under varying degree of xenoliths from picrobasalts and basanites of the Vitim Plateau melting (F ¼ 0.039e0.050). Contaminated basaltsetrachybasalts (Ionov et al., 1997; Glaser et al., 1999; Litasov et al., 2000a,b; with elevated (La/Yb)N fall at the line of sources with 12% of garnet. Ashchepkov and Andre, 2002; Ashchepkov et al., 2003, 2011). In the (La/Yb)N vs. YbN diagram, the data points of lavas erupted Additional trace-element data on Ta, Th, U concentrations in about 14e13 Ma are plotted between model lines that show garnet amphibole were taken for xenoliths from basalts of the SW Tian content 3.5% and no garnet in a source region. Firstly, melted ma- Shan (Zheng et al., 2006). The data on the type-A apatite were used terial from the shallow garnet-free source (or containing a small after O’Reilly and Griffin (2000) (Table 2). amount of garnet) provided voluminous eruptions (2/3 of the Equilibrium partial melting was modeled using the equation whole 14e13 Ma volume). The degree of partial melting was at a ð Þ ¼ Ci 0 e maximum (F ¼ 0.15e0.20). The mantle source for coeval small (Shaw, 1970): Ci D þFð1P Þ, where Ci concentration of element i i i portion of Na-basaltic lavas, derived from the garnet-bearing in a melt; C (0) e concentration of element i in the source; D e bulk i i source (1e2% of Grt), experienced relatively weak fusion distribution coefficient for the source material, F e melt fraction; P (F ¼ 0.05e0.07). The final lava portion, which amounted to 1/4 of P ¼ XjKj e bulk distribution coefficient for material involved in i i the volume of all the 14e13 Ma eruptions, was also derived from melting; Xj e weight fraction of phase j; Kj e mineral/melt distri- the garnet-bearing source (up to 3.5% Grt) with an elevated degree bution coefficient for phase j. Mineral/melt distribution coefficients of partial melting (F ¼ 0.05e0.10). are given in Table 3. The clinopyroxene fraction in the source ma- terial was increased about two times and the olivine fraction was decreased relative to their contents in the source according to the 6.7. Change of magmatism at 16e13 Ma experimental study of garnet-bearing rocks (Walter, 1998). From results of trace-element modeling, we infer that the 16e14 We suggest that the mid-Miocene temporal change from high- Ma high-Mg basanitesebasalts and picrobasaltsebasalts were to moderate-Mg liquids in the Bereya volcanic center was due to originated from a source that has: Ol (48e57%), Opx (25%), Grt a thermal impact on the lithosphere from the convective mantle of (6e12%), Cpx (5%), Amph (5e7%), Phl (2e3%), and Ap (<0.07%). the Transbaikalian melting domain (Fig. 16). Some discrepancies between rocks and model melts were found in In the time interval of 16e14 Ma, initial high-Mg magmas were terms of Rb, Th, and K. Besides, the maximum of Sr in the trace contaminated by disintegrated crustal material. The composition of element patterns are not explained by this model. Melt fractions (F) primary magmas that rise from the mantle source was distorted vary from 0.04 to 0.06 (Fig. 15a and b). due to the partial assimilation of this material. Perhaps, the primary In Fig. 13, some data points of high-Mg basanitesebasalts fall at mantle melts were similar in composition to the high-Mg basani- a line of the melt fraction F ¼ w0.045. Neodymium isotope ratios of teebasaltic melts erupted in the Bereya volcanic center after high- this sub-group decrease as the garnet content in the source in- Mg contaminated basaltsetrachybasalts. creases from 6 to 12%. Another sub-group of high-Mg basani- The high-Mg basaniteebasaltic and picrobasaltsebasaltic tesebasalts yields melt fractions increasing up to 0.075 in the magmas of the second and third episodes exhibited uncontami- source with 6e8% of garnet. The neodymium isotope ratios of this nated primary mantle liquids. In adiabatically ascending mantle sub-group range in narrow interval as Nd concentrations decrease material, melting depth depends on a heat content of a magmatic due to increasing melt fractions. The source of picrobasaltsebasalts system. When heat supply is high, melting is initiated at a deeper shows about 10% Grt and melt fraction F ¼ 0.04e0.05. level; when heat supply is low, melting is only provided in rela- The 14e13 Ma lavas were derived from sources with relatively tively shallow conditions (Fukuyama, 1985). low garnet (or garnet-free) and elevated clinopyroxene contents: The mid-Miocene eruptions of high-Mg magmatic liquids in the basalts and trachybasalts e Ol (49.5%), Opx (25%), Grt (<3.5%), Bereya volcanic center from the 115e150 km mantle source were Cpx (12%), Amph (8%), Phl (2%), Ap (<0.07%), F ¼ 0.05e0.1; due to adiabatic introduction of heat enough to melting at such a basaltsebasaltic andesites e Ol (54.5%), Opx (25%), Cpx (15%), deep level. Created buoyancy of partially molten material exceeded Amph (5%), Phl (<0.5%), Ap (<0.05%), F ¼ 0.15e0.20. that of crystalline phases in the ambient asthenospheric mantle. As In basaltsetrachybasalts of KeNa series, LREE are variable in a result, the melts, released from the matrix, moved upwards and transitional (NaeKeNa) rocks, while HREE remain constant erupted on the surface. (Fig. 8b). These REE patterns are typical for equilibrium partial The deepest high-Mg basanitesebasalts were produced through melting in a mantle source. Increasing La/Yb reflects decreasing relatively homogeneous melting by excessive heat that caused degree of melting (Rollinson, 1993). moderate melting in the magmatic column F ¼ 6%. An initial melt From the results of trace element modeling in the (La/Yb) vs. N region had the upper limit at the depth 115 km governed by release Yb coordinates, high-Mg basanitesebasalts are examined as de- N of buoyant high-Mg liquids from a fractionated matrix with the rivatives from a source with the garnet content from 6 to 12%. A garnet content varying from 6 to 12%. Shortly afterward melting high-garnet trend (BSN e high Grt) in Fig. 12 corresponds to 6e12% processes localized within a region that had more constant garnet of Grt at relatively low degree of melting (F ¼ 0.042e0.046) in content and yielded picrobasaltebasaltic liquids. Fig.13a. A similar source was identified for basanites from the Kamar At the initial stage of 16e14 Ma, igneous channels provided heat volcanic field located near the South Baikal basin (Rasskazov et al., supply to the shallow part of the lithospheric mantle. High-degree 2013a,b). In the source region, garnet concentrations might in- decompression melting of a garnet-free (or low-garnet) source crease due to depth-dependent increasing rock density and due to caused thermal impact on the lithosphere, uplift and erosional garnet fractionation relative to less dense and olivine. dissection of the area and resulted in voluminous eruptions of Some data points of the high-Mg basaniteebasalt group from the moderate-Mg basalts and basaltic andesites about 14e13 Ma. This Bereya volcanic center are shifted, however, along a line with garnet impulse of the shallow melting, however, was short-lived and

Please cite this article in press as: Chuvashova, I., et al., Mid-Miocene thermal impact on the lithosphere by sub-lithospheric convective mantle material: Transition from high- to moderate-Mg magmatism beneath Vitim Plateau, Siberia, Geoscience Frontiers (2016), http://dx.doi.org/ 10.1016/j.gsf.2016.05.011 18 I. Chuvashova et al. / Geoscience Frontiers xxx (2016) 1e22

Figure 15. Comparisons of trace element (a, s) and REE (b, d) patterns of rocks (lines with markers) and calculated model liquids (black lines without markers). Symbols are as in Fig. 6. (a, b) e high-Mg basanitesebasalts and picrobasaltsebasalts of 16e14 Ma, (c, d) e moderate-Mg basalts, basaltic andesites, and trachybasalts of 14e13 Ma. Normalization values are after McDonough and Sun (1995). The calculated degrees of melting (melt fractions) and garnet contents in sources are indicated (see discussion in the text).

changed to less voluminous melting in the deeper garnet-bearing reflect adiabatic heat transfer, whereas the less advanced melting at part of the lithospheric mantle. shallower depths (up to 115 km) was caused by relative decreasing temperature beneath the lithosphere. Variations of major oxides, trace elements, and isotopes in the 14e13 Ma rock sequence 7. Conclusion showed initial melts equilibrated at the shallow garnet-free mantle source with subsequently reduced degree of melting of garnet- A comparative study of high- and moderate-Mg volcanic rocks bearing material. The shallow melts revealed no melting of rela- erupted in the Bereya volcanic center of the Vitim Plateau, at 16e14 tively depleted lithospheric material, evidenced by mantle xeno- and 14e13 Ma, respectively, has revealed specific temporal varia- liths from the lithosphere. tions of rock compositions in these units. The former (low-volume) The main purpose of this paper is to connect the mid-Miocene unit shows the sequence of high-Mg rocks: (1) contaminated volcanic activity, initiated in deep and further developed in basaltsetrachybasalts of KeNa series, (2) basanitesebasalts of shallow mantle sources, with destabilization and rifting of the transitional (KeNaeK) compositions, and (3) picrobasaltsebasalts lithosphere in the Baikal Rift Zone. It is suggested that voluminous of K series. The latter (voluminous) unit demonstrates the sequence lava eruptions from the shallow lithospheric source at about of moderate-Mg rocks: (1) basaltsebasaltic andesites of transitional 14e13 Ma reflected a thermal impact on the lithosphere beneath (NaeKeNa) compositions and (2) basaltsetrachybasalts of KeNa the Vitim volcanic field by hot convective material of the Trans- series. baikalian melting domain that had a potential temperature T as It has been found that the initial portion of high-Mg liquids in p high as 1510 C, comparable with that parameter in the upper the 16e14 Ma unit was strongly affected by crustal contamination. mantle of the Ontong Java Plateau and Iceland hot spots. This Later on, the uncontaminated melt portions of this unit revealed thermal impact caused extension of the lithosphere, accompanied various temperature effects in a mantle anomaly; the elevated by subsidence of rift basins on the surface. In other volcanic fields of degree of melting in its source as deep as 150 km appeared to

Please cite this article in press as: Chuvashova, I., et al., Mid-Miocene thermal impact on the lithosphere by sub-lithospheric convective mantle material: Transition from high- to moderate-Mg magmatism beneath Vitim Plateau, Siberia, Geoscience Frontiers (2016), http://dx.doi.org/ 10.1016/j.gsf.2016.05.011 I. Chuvashova et al. / Geoscience Frontiers xxx (2016) 1e22 19

Figure 16. Change of sources for rocks of the Bereya volcanic center in terms of the mid-Miocene thermal impact on the lithosphere by a sub-lithospheric hot convective material from the Transbaikalian melting domain. Tp e potential temperature, F e degree of partial melting.

Please cite this article in press as: Chuvashova, I., et al., Mid-Miocene thermal impact on the lithosphere by sub-lithospheric convective mantle material: Transition from high- to moderate-Mg magmatism beneath Vitim Plateau, Siberia, Geoscience Frontiers (2016), http://dx.doi.org/ 10.1016/j.gsf.2016.05.011 20 I. Chuvashova et al. / Geoscience Frontiers xxx (2016) 1e22

e the Transbaikalian melting domain (Udokan, Dariganga, Kamar), Tp petrogenetic model of magmatism at the asthenosphere lithosphere bound- e estimates decreased to 1430 C, whereas in the shallower ary. Geodynamics & Tectonophysics 4, 385 407. http://dx.doi.org/10.5800/GT- 2012-3-4-0081. SayaneMongolian melting domain, Tp did not exceed 1300 C. Chuvashova, I.S., Rasskazov, S.V., Yasnygina, T.A., Rudneva, N.A., 2016. Activation and termination of late Cenozoic extension in the lithosphere of the marginal part of the Baikal Rift Zone: changing sources of volcanism in the Vitim Plateau. Acknowledgments Volcanology and Seismology (in press). DePaolo, D.J., 1981. Trace element and isotopic effects of combined wallrock This work was supported by the Russian Science Foundation for assimilation and fractional crystallization. Earth and Planetary Science Letters 53, 189e202. Basic Research (project 14-05-31328). Sample preparation for Faure, G., 2001. Origin of Igneous Rocks: the Isotopic Evidence. Springer, p. 496. trace-element and isotope analytical work was performed by M.E. Fedorenko, V.A., Czamanske, G.K., 1997. Results of new field and geochemical Markova and E.V. Saranina, isotope mass-spectrometric de- studies of the volcanic and intrusive rocks of the Maymecha-Kotuy area, Si- berian flood-basalt province, Russia. International Geology Review 39, terminations by N.N. Fefelov, and major-oxide determinations by 479e531. G.V. Bondareva and M.M. Samoilenko. Agilent 7500ce and Finnigan Fedorenko, V.A., Lightfoot, C., Naldrett, A.J., Czamanske, G., Hawkesworth, C., MAT 262 mass spectrometers are affiliated with the Common Use Wooden, J., Ebel, D., 1996. Petrogenesis of the Siberian flood-basalt sequence at ’ e Centers “Ultramicroanalysis” and “Geodynamics and Geochro- Noril sk, north central Siberia. International Geology Review 38, 99 135. Flower, M.F.J., Chung, Sun-Lin, Lo, Ching-Hua, Lee, Tung-Yi (Eds.), 1998. Mantle nology”, respectively. We thank Nick Roberts and two anonymous Dynamics and Plate Interactions in East Asia. Geodynamic Series Volume 27. reviewers for constructive discussion of the manuscript. Am. Geophysical Union, Washington D.C, p. 418. Foley, S.F., Venturelli, G., Green, D.H., Toscani, L., 1987. The ultrapotassic rocks: characteristics, classification, and constraints for petrogenetic models. Earth- Appendix A. Supplementary data Science Reviews 24, 81e134. Foley, S.F., Jackson, S.E., Fryer, B.J., Greenough, I.D., Jenner, G.A., 1996. Trace element partition coefficients for clinopyroxene and phlogopite in an alkaline lamp- Supplementary data related to this article can be found at http:// rophyre from Newfoundland by LAM-ICP-MS. Geochimica et Cosmochimica dx.doi.org/10.1016/j.gsf.2016.05.011. Acta 60, 629e638. Foley, S.F., Yaxley, G.M., Rosenthal, A., et al., 2009. The composition of near-solidus melts of peridotite in the presence of CO2 and H2O between 40 and 60 kbar. References Lithos 112S, 274e283. Foulger, G.R., 2010. Plates vs. Plumes: a Geological Controversy. WileyeBlackwell, Anderson, D.L., 2007. New Theory of the Earth. Cambridge University Press, Cam- p. 328. bridge, p. 384. Fukuyama, H., 1985. Heat fusion of basaltic magma. Earth and Planetary Science Anderson, D.L., Tanimoto, T., Zhang, Y.S., 1992. Plume heads, continental litho- Letters 73, 407e414. sphere, flood basalts and tomography. In: Storrey, B.C., Albaster, T., Ghosh, S., Ohtani, E., Litasov, K.D., Suzuki, A., Sakamaki, T., 2007. Stability of Pankhurst, R.J. (Eds.), Magmatism and the Causes of Continental Break-up. carbonated basaltic melt at the base of the Earth’s upper mantle. Geophysical Geological Society Special Publication, 68, pp. 99e124. Research Letters 34, L22312. http://dx.doi.org/10.1029/2007GL031349. Arndt, N., Lesher, C.M., Barnes, S.J., 2008. Komatiite. Cambridge University Press, Ghosh, S., Ohtani, E., Litasov, K.D., Teresaki, H., 2009. Solidus of carbonated peri- p. 458. dotite from 10 to 20 GPa and origin of magnesiocarbonated melt in the Earth’s Ashchepkov, I.V., 1991. Deep Xenoliths of the Baikal Rift. Nauka, Sib. Branch, deep mantle. Chemical Geology 262, 17e28. Novosibirsk, p. 160. Girnis, A.V., Bulatov, V.K., Lahaye, Y., Brey, G.P., 2006. Partitioning of trace elements Ashchepkov, I.V., Andre, L., 2002. Pyroxenite xenoliths in picrite basalts (Vitim between carbonate-silicate melts and mantle minerals: experiment and Plateau): origin and differentiation of mantle melts: example. Russian Geology petrological consequences. Petrology 14, 492e514. and Geophysics 43 (4), 343e363. Glaser, S.M., Foley, S.F., Günter, D., 1999. Trace element compositions of minerals in Ashchepkov, I.V., Travin, A.V., Saprykin, A.I., Andre, L., Gerasimov, P.A., garnet and spinel peridotite xenoliths from the Vitim volcanic field, Trans- Khmelnikova, O.S., 2003. Age of xenolith-bearing basalts and mantle evolu- baikalia, eastern Siberia. Lithos 48, 263e285. tion in the Baikal rift zone. Russian Geology and Geophysics 44 (11), Goncharov, A.G., Ionov, D.A., 2012. Redox state of deep off-craton lithospheric 1121e1149. mantle: new data from garnet and spinel peridotites from Vitim, southern Ashchepkov, I.V., Andre, L., Downes, H., Belyatsky, B.A., 2011. Pyroxenites and Siberia. Contributions to Mineralogy and Petrology 164, 731e745. http:// megacrysts from Vitim picrate-basalts (Russia): polybaric fractionation of rising dx.doi.org/10.1007/s00410-012-0767-z. melts in the mantle? Journal of Asian Earth Sciences 42, 14e37. Green, T.H., Sie, S.H., Ryan, C.G., Cousens, D.R., 1989. Proton microprobe-determined Asimow, P.D., 2005. Igneous Processes. In: Selley, R.C., Cocks, L.R.M., Plimer, I.R. partitioning of Nb, Ta, Zr, Sr and Y between garnet, clinopyroxene and basaltic (Eds.), Encyclopedia of Geology 2. Elsevier, Oxford, UK, pp. 209e217. magma at high pressure and temperature. Chemical Geology 74, 201e216. Barry, T.L., Saunders, A.D., Kempton, P.D., Windley, B.F., Pringle, M.S., Dorjnamjaa, D., Green, T.H., Blundy, J.D., Adam, J., Yaxley, G.M., 2000. SIMS determination of trace Saandar, S., 2003. Petrogenesis of Cenozoic basalts from Mongolia: evidence for element partition coefficients between garnet, clinopyroxene and hydrous the role of asthenospheric versus metasomatized lithospheric mantle sources. basaltic liquids at 2e7.5 GPa and 1080e1200 C. Lithos 53, 165e187. Journal of Petrology 44 (1), 55e91. Gregoire, M., Moine, B.N., O’Reilly, S.Y., Gottin, J.Y., Giret, A., 2000. Trace element Barry, T.L., Ivanov, A.V., Rasskazov, S.V., Demonterova, E.I., Dunai, T.J., Davies, G.R., residence and partitioning in mantle xenoliths metasomatized by highly alka- Harrison, D., 2007. Helium isotopes provide no evidence for deep mantle line, silicate- and carbonate-rich melts (Kerguelen Islands, Indian Ocean). involvement in widespread Cenozoic volcanism across Central Asia. Lithos 95, Journal of Petrology 41, 477e509. 415e424. Griffin, W.L., O’Reilly, S.Y., Stabel, C.G., 1988. Mantle metasomatism beneath western Belichenko, V.G., 1977. Caledonides of the Baikal Mountain Region. Nauka Publisher, Victoria, Australia: II. Isotopic geochemistry of Cr-diopside lherzolites and Al- Novosibirsk, p. 134. augite pyroxenites. Geochimica et Cosmochimica Acta 52, 449e459. Belichenko, V.G., Geletii, N.K., Barash, I.G., 2006. The Barguzin microcontinent Halliday, A.N., Lee, D.-C., Tommasini, S., Davies, G.R., Paslick, C.R., Fitton, J.G., (Baikal mountain area): the problem of outlining. Russian Geology and James, D.E., 1995. Incompatible trace elements in OIB and MORB and source Geophysics 47 (10), 1035e1045. enrichment in the sub-oceanic mantle. Earth and Planetary Science Letters 133, Carlson, R.W., Czamanske, G., Fedorenko, V., Ilupin, I., 2006. A comparison of Si- 379e395. berian meimechites and kimberlites: implications for the source of high-Mg Hanano, D., Weis, D., Scoates, J.S., Aciego, S., DePaolo, D.J., 2010. Horizontal and alkalic magmas and flood basalts. Geochemistry, Geophysics, Geosystems 7, vertical zoning of heterogeneities in the Hawaiian mantle plume from the Q11014. http://dx.doi.org/10.1029/2006GC001342. geochemistry of consecutive postshield volcano pairs: Kohala-Mahukona and Chazot, G., Menzies, M.A., Harte, B., 1996. Determinations of partition coefficients Mauna KeaeHualalai. Geochemistry, Geophysics, Geosystems 11 (1), Q01004. between apatite, clinopyroxene, amphibole, and melt in natural spinel lher- http://dx.doi.org/10.1029/2009GC002782. zolites from Yemen: implications for wet melting of the lithospheric mantle. Hanski, E., Walker, R.J., Huhma, H., Polyakov, G.V., Balykin, P.A., Hoa, T.T., Geochimica et Cosmochimica Acta 60, 423e437. Phuong, N.T., 2004. Origin of the Permo-Triassic komatiites, northwestern Chernyaeva, G.P., Lyamina, N.A., Rasskazov, S.V., Rezanov, I.N., Savinova, V.V., 2007. Vietnam. Contributions to Mineralogy and Petrology 147, 453e469. Biostratigraphy and deposition environment conditions of accumulation of the Harris, N.R., 1998. Isotopic, Geochemical, and Geochronological Constraints on the Middle-Late Miocene volcano-sedimentary strata in the Dzhilinda basin, Origin and Evolution of Cenozoic Volcanism, Baikal Rift Zone, Siberia. PH Western Transbaikalia. Russian Geology and Geophysics 48 (4), 361e370. Dissertation. Massachusetts Institute of Technology, USA, p. 440. Chuvashova, I.S., 2013. Convective upper mantle homogenization in Asia: phe- Hart, S.R., Dunn, T., 1993. Experimental cpx/melt partitioning of 24 trace elements. nomenology of sources for recent volcanic rocks in the Sayan-Mongolian low- Contributions to Mineralogy and Petrology 113, 1e8. velocity domain. Izvestia of the Irkutsk State University. A series of Earth Sci- Hauri, E.H., Wagner, T.P., Grove, T.L., 1994. Experimental and natural partitioning of ences 6 (2), 187e198. Th, U, Pb and other trace elements between garnet, clinopyroxene and basaltic Chuvashova, I.S., Rasskazov, S.V., Yasnygina, T.A., Mikheeva, E.A., 2012. High-Mg melts. Chemical Geology 117, 149e166. lavas from the Dariganga volcanic field in the South-Eastern Mongolia:

Please cite this article in press as: Chuvashova, I., et al., Mid-Miocene thermal impact on the lithosphere by sub-lithospheric convective mantle material: Transition from high- to moderate-Mg magmatism beneath Vitim Plateau, Siberia, Geoscience Frontiers (2016), http://dx.doi.org/ 10.1016/j.gsf.2016.05.011 I. Chuvashova et al. / Geoscience Frontiers xxx (2016) 1e22 21

Henderson, P., 1982. Inorganic Geochemistry. Pergamon Press, Oxford-New York- O’Reilly, S.Y., Griffin, W.L., Ryan, C.G., 1991. Residence of trace element in meta- Toronto-Sydney-Paris-Frankfurt, p. 310. somatised spinel lherzolite xenoliths: a proton-microprobe study. Contribu- Herzberg, C., 2011. Identification of source lithology in the Hawaiian and Canary tions to Mineralogy and Petrology 109, 98e113. islands: implications for origins. Journal of Petrology 52 (1), 113e146. Parfenov, L.M., Popeko, L.I., Tomurtogo, O., 1999. Problems of tectonics in the Herzberg, S., Asimov, N., Arndt, N., Niu, Y., Lesher, C.M., Fitton, J.G., Cheadle, M.J., Mongol-Okhotsk orogenic belt. Russian Pacific Geology 18 (5), 24e43. Saunders, A.D., 2007. Temperatures in ambient mantle and plumes: Constraints Perchuk, L.L., Safonov, O.G., Yapaskurt, V.O., Barton Jr., J.M., 2002. Crystalemelt from basalts, picrites, and komatiites. Geochemistry, Geophysics, Geosystems 8, equilibria involving potassium-bearing clinopyroxene as indicator of mantle- Q02006. http://dx.doi.org/10.1029/2006GC001390. derived ultra-high-potassic liquids: an analytical review. Lithos 60, 89e111. Hofmann, A.W., 1997. Mantle geochemistry: the message from oceanic volcanism. Rasskazov, S.V.,1985. The Udokan Basaltoids. Nauka, SB Ac. Sci. USSR, Novosibirsk, p.142. Nature 385, 219e229. Rasskazov, S.V., 1991. Hot spot volcanism and structure of the western part of the Ionov, D.A., Ashchepkov, I.V., Stosch, H.-G., Witt-Eickschen, G., Seck, H.A., 1993. Baikal Rift System. Geology and Geophysics 9, 72e81. Garnet peridotite xenoliths from the Vitim Volcanic Field, Baikal Region e the Rasskazov, S.V., 1993. Magmatism of the Baikal Rift System. Nauka, Siberian Pub- nature of the garnet spinel peridotite transition zone in the continental mantle. lishing Firm, Novosibirsk, p. 288. Journal of Petrology 34, 1141e1175. Rasskazov, S.V., 1994. Comparison of volcanism and hot-spot structures of Yel- Ionov, D.A., Griffin, W.L., O’Reilly, S.Y., 1997. Volatile-bearing minerals and lithophile lowstone and Eastern Sayans. Geology and Geophysics 10, 67e75. trace elements in the upper mantle. Chemical Geology 141, 153e184. Rasskazov, S.V., Batyrmurzaev, A.S., 1985. Cenozoic basalts of the Vitim Plateau and Ionov, D.A., Ashchepkov, I., Jagoutz, E., 2005. The provenance of fertile off-craton determination of their ages. Geology and Geophysics 5, 20e28. lithospheric mantle: SreNd isotope and chemical composition of garnet and Rasskazov, S.V., Chuvashova, I.S., 2013. The Latest Mantle Geodynamics of Central spinel peridotite xenoliths from Vitim, Siberia. Chemical Geology 217, 41e75. Asia. Irkutsk State University Press, Irkutsk, p. 308. Johnson, J.S., Gibson, S.A., Thompson, R.N., Nowell, G.M., 2005. Volcanism in the Rasskazov, S.V., Boven, A., Andre, L., Liegeois, J.P., Ivanov, A.V., Punzalan, L., 1997. Vitim volcanic field, Siberia: Geochemical evidence for a mantle plume beneath Evolution of magmatism in the Northeastern Baikal Rift System. Petrology 5 (2), the Baikal Rift Zone. Journal of Petrology 46 (7), 1309e1344. 101e120. Karato, S-i, 2012. On the origin of the asthenosphere. Earth and Planetary Science Rasskazov, S.V., Logatchev, N.A., Brandt, I.S., Brandt, S.B., Ivanov, A.V., 2000. Letters 321e322, 95e103. Geochronology and Geodynamics in the Late Cenozoic (South Siberia e South Kennedy, A.K., Lofgren, G.E., Wasserburg, G.J., 1993. An experimental study of trace and East Asia). Nauka, SB RAS, Novosibirsk, p. 288. element partitioning between olivine, orthopyroxene and melt in chondrulites: Rasskazov, S.V., Logatchev, N.A., Ivanov, A.V., Misharina, V.A., Chernyaeva, G.P., equilibrium values and kinetic effects. Earth and Planetary Science Letters 115, Brandt, I.S., Brandt, S.B., Skoblo, V.M., Lyamina, N.A., 2001. Palynological and 177e195. diatom analyses of sediments from the Late Cenozoic Paleo-Amalat valley Kerr, A.C., Arndt, N.T., 2001. A note on the IUGS reclassification on high-Mg and (Western Transbaikalia). Russian Geology and Geophysics 42 (5), 773e785. picritic volcanic rocks. Journal of Petrology 42 (11), 2169e2171. Rasskazov, S.V., Bowring, S.A., Housh, T., Demonterova, E.I., Logatchev, N.A., Kiselev, A.I., Medvedev, M.E., Golovko, G.A., 1979. Volcanism of the Baikal Rift Zone Ivanov, A.V., Saranina, E.V., Maslovskaya, M.N., 2002. The Pb, Nd, and Sr isotope and Problems of Deep Magma Formation. Nauka, Novosibirsk, p. 200. systematics in heterogeneous continental lithosphere above a convecting Klemme, S., 2004. The influence of Cr on the garnetespinel transition in the Earth’s mantle domain. Doklady Earth Sciences 387 (9), 1056e1059. mantle: experiments in the system MgOeCr2O3eSiO2 and thermodynamic Rasskazov, S.V., Logatchev, N.A., Kozhevnikov, V.M., Yanovskaya, T.B., 2003a. modeling. Lithos 77, 639e646. http://dx.doi.org/10.1016/j.lithos.2004.03.017. Multistage dynamics of the upper mantle in Eastern Asia: relationships be- Klemme, S., O’Neill, H.StC., 2000. The near-solidus transition from garnet lher- tween wandering volcanism and low-velocity anomalies. Doklady Earth Sci- zolite to spinel lherzolite. Contributions to Mineralogy and Petrology 138, ences 390 (4), 492e496. 237e248. Rasskazov, S.V., Luhr, J.F., Bowring, S.A., Ivanov, A.V., Brandt, I.S., Brandt, S.B., Kogarko, L.N., Ryabchikov, I.D., 2000. Geochemical evidence for meimechite magma Demonterova, E.I., Boven, A.A., Kunk, M., Housh, T., Dungan, M.A., 2003b. Late generation in the subcontinental lithosphere of Polar Siberia. Journal of Asian Cenozoic volcanism in the Baikal rift system: evidence for formation of the Baikal Earth Sciences 18, 195e203. and Khubsugul basins due to thermal impacts on the lithosphere and collision- Koulakov, I., Bushenkova, N., 2010. Upper mantle structure beneath the Siberian derived tectonic stress. Berliner palaobiologische abhandlungen 4, 33e48. craton and surrounding areas based on regional tomographic inversion of P and Rasskazov, S., Taniguchi, H., Goto, A., Litasov, K., 2004. Magmatic expression of plate PP travel times. Tectonophysics 486, 81e100. subduction beneath East Asia in the Mesozoic through Cenozoic. Northeast Kutolin, V.A., 1966. To the problem of the order and temperatures of mineral crys- Asian Studies 9, 179e219. tallization in rocks of basic compositions. Geology and Geophysics 11, 42e51. Rasskazov, S.V., Lyamina, N.A., Chernyaeva, G.P., Skoblo, V.M., Luzin, I.V., LaTourrette, T., Hervig, R.L., Holloway, J.R., 1995. Trace element partitioning between Rudnev, A.F., Rezanov, I.N., 2007. Cenozoic Stratigraphy of the Vitim Plateau: amphibole, phlogopite and basanite melt. Earth and Planetary Science Letters The Phenomenon of Long-term Rifting in the South of Eastern Siberia. Aca- 135, 13e30. demic Publishing House “Geo”, Novosibirsk, p. 193. Le Bas, M.J., 2000. IUGS reclassification on high-Mg and picritic volcanic rocks. Rasskazov, S.V., Brandt, S.B., Brandt, I.S., 2010. Radiogenic Isotopes in Geologic Journal of Petrology 41 (10), 1467e1470. Processes. Springer, p. 306. Le Bas, M.J., Streckeisen, A.L., 1991. The IUGS systematics of igneous rocks. Journal of Rasskazov, S.V., Chuvashova, I.S., Yasnygina, T.A., Fefelov, N.N., Saranina, E.V., 2012. the Geological Society London 148, 825e833. Potassic and Potassicesodic Volcanic Series in the Cenozoic of Asia. Academic Lei, J., Zhao, D., 2005. P-wave tomography and origin of the Changbai intraplate Publishing House “GEO”, Novosibirsk, p. 351. volcano in Northeast Asia. Tectonophysics 397, 281e295. Rasskazov, S.V., Yasnygina, T.A., Chuvashova, I.S., Markova, M.E., 2013a. Cenozoic Litasov, K., Taniguchi, H., 2002. Mantle Evolution beneath the Baikal Rift. CNEAS Volcanic Rocks of Asia: Database on Chemical and Isotopic Compositions. Cer- Monograph Series (5). Center for Northeast Asian Studies Tohoku University, tificate of state registration of the database No. 2013620655. p. 221. Rasskazov, S.V., Yasnygina, T.A., Chuvashova, I.S., Mikheeva, E.A., Snopkov, S.V., Litasov, K.D., Foley, S.F., Litasov, Yu.D., 2000a. Magmatic modification and meta- 2013b. The Kultuk Volcano: spatial-temporal change of magmatic sources at the somatism of the subcontinental mantle beneath the Vitim volcanic field (East western terminus of the South Baikal basin between 18 and 12 Ma. Geo- Siberia): evidence from trace element data on pyroxenite and peridotite xe- dynamics & Tectonophysics 4 (2), 135e168. http://dx.doi.org/10.5800/ noliths from Miocene picrobasalt. Lithos 54, 83e114. GT2013420095. Litasov, K.D., Litasov, Y.D., Mekhonoshin, A.S., Malkovets, V.G., 2000b. Mineralogy of Rasskazov, S.V., Chuvashova, I.S., Mikolaichuk, A.V., Sobel, E.R., Yasnygina, T.A., mantle xenoliths from Pliocene basanites of Dzhilinda river (Vitim volcanic Fefelov, N.N., Saranina, E.V., 2015. Lateral change of sources for the Cretaceous- field). Geology and Geophysics 41 (11), 1534e1556. Paleogene magmatism of the Tian Shan. Petrology 23 (3), 281e308. Maruyama, S., Santosh, M., Zhao, D., 2007. Superplume, supercontinent, and post- Robinson, J.A.C., Wood, B.J., 1998. The depth of the spinel to garnet transition at the perovskite: mantle dynamics and anti-plate tectonics on the coreemantle peridotite solidus. Earth and Planetary Science Letters 164, 277e284. boundary. Gondwana Research 11, 7e37. Rogozhina, V.A., Kozhevnikov, V.M., 1979. Region of Anomalous Mantle beneath the McDonough, W.F., Sun, S.-S., 1995. The composition of the Earth. Chemical Geology Baikal Rift. Nauka, SB RAS, Novosibirsk, p. 104. 120, 223e253. Rollinson, H.R., 1993. Using Geochemical Data: Evaluation, Presentation, Interpre- Melnikov, A.I., Mazukabzov, A.M., Sklyarov, E.V., Vasiljev, E.P., 1994. Baikal rift tation. Longman, Harlow, p. 310. basement: structure and tectonic evolution. Bulletin des Centres de Recherches Ruzhentsev, S.V., Minina, O.R., Nekrasov, G.E., Aristov, V.A., Golionko, B.G., Exploration e Production Elf Aquitaine 18 (1), 99e122. Doronina, N.A., Lykhin, D.A., 2012. The Baikal-Vitim Fold System: structure and Morgan, W.J., 1971. Convective plumes in the lower mantle. Nature 230, 42e43. geodynamic evolution. Geotectonics (2), 3e28. Nakamura, E., McCulloch, M.T., Campbell, I.H., 1990. Chemical geodynamics in the Safonov, O.G., Perchuk, L.L., Litvin, Yu.A., 2005. Equilibrium K-bearing back arc region of Japan based on the trace element and SreNd isotope com- clinopyroxene-melt as a model for barometry of mantle derived mineral as- positions. Tectonophysics 174, 207e233. semblages. Russian Geology and Geophysics 46 (12), 1300e1316. O’Reilly, S.Y., Griffin, W.L., 1995. Trace-element partitioning between garnet and Scarrow, J.H., Cox, K.G., 1995. Basalts generated by decompressive adiabatic melting clinopyroxene in mantle-derived pyroxenites and eclogites: PeTeX control. of a mantle plume: a case study from the Isle of Skye, NW Scotland. Journal of Chemical Geology 121, 105e130. Petrology 36, 3e22. O’Reilly, S.Y., Griffin, W.L., 2000. Apatite in the mantle: implications for meta- Sengör, A.M., Natal’in, B.A., 1996. In: Yin, A., Harrison, M. (Eds.), Paleotectonics of somatic processes and high heat production in Phanerozoic mantle. Lithos 53, Asia: Fragments of a Synthesis. The Tectonic Evolution of Asia, 21. Cambridge 217e232. University Press, Cambridge, pp. 486e640.

Please cite this article in press as: Chuvashova, I., et al., Mid-Miocene thermal impact on the lithosphere by sub-lithospheric convective mantle material: Transition from high- to moderate-Mg magmatism beneath Vitim Plateau, Siberia, Geoscience Frontiers (2016), http://dx.doi.org/ 10.1016/j.gsf.2016.05.011 22 I. Chuvashova et al. / Geoscience Frontiers xxx (2016) 1e22

Shatsky, V.S., Ragozin, A.L., Sobolev, N.V., 2006. Some aspects of the metamorphic Journal of Petrology 47 (10), 1997e2019. http://dx.doi.org/10.1093/petrology/ evolution ultrahigh-pressure calc-silicate rocks of the Kokchetav Massif. egl034. Russian Geology and Geophysics 47 (1), 105e119. Zhao, D., 2009. Multiscale seismic tomography and mantle dynamics. Gondwana Shaw, D.M., 1970. Trace element fractionation during anatexis. Geochimica et Research 15, 297e323. http://dx.doi.org/10.1016/j.gr.2008.07.003. Cosmochimica Acta 34, 237e243. Zhao, D., Tian, Y., Lei, J., Liu, L., Zheng, S., 2009. Seismic image and origin of the Shellnutt, J.G., 2014. The Emeishan large igneous province. A synthesis Geoscience Changbai intraplate volcano in East Asia: role of big mantle wedge above the Frontiers 5, 369e394. http://dx.doi.org/10.1016/j.gsf.2013.07.003. stagnant pacific slab. Physics of the Earth and Planetary Interiors 173 (3e4), Simonov, V.A., Kudryashova, E.A., Yarmolyuk, V.V., Kovyazin, S.V., Kotlyarov, A.V., 197e206. Stupakov, S.I., 2012. Physico-Chemical Conditions of Petrogenesis of Late Basalt Zheng, J., Griffin, W.L., O’Reilly, S.Y., Zhang, M., Pearson, N., Luo, Zh, 2006. The Complexes South Baikal and South Khangai Volcanic Regions, Central Asia (Data lithospheric mantle beneath the southwestern Tianshan area, northwest China. on Melt Inclusions). Geodynamic Evolution of the Lithosphere of the Central Contributions to Mineralogy and Petrology 151, 456e479. Asian Mobile Belt (from Ocean to Continent), 2. Institute of the Earth’s Crust SB Zhi, X., Song, Y., Frey, F.A., Feng, J., Zhai, M., 1990. Geochemistry of Hannuoba basalts, RAS, Irkutsk, pp. 75e77. eastern China: constraints on the origin of continental alkalic and tholeiitic Sobolev, N.V., 1974. Deep Inclusions in Kimberlites and the Problem of the basalt. Chemical Geology 88 (1/2), 1e33. Composition of the Upper Mantle. Nauka, Novosibirsk, p. 264. Zindler, A., Hart, S., 1986. Chemical geodynamics. Annual Review of Earth and Sobolev, A.V., Hofmann, A.W., Sobolev, S.V., Nikogosian, I.K., 2005. An olivine-free Planetary Sciences 14, 493e571. mantle source of Hawaiian shield basalts. Nature 434, 590e597. Zonenshain, L.P., Kuzmin, M.I., Bocharova, N.Y., 1991. Hot field tectonics. Tectono- Sobolev, A.V., Sobolev, S.V., Kuzmin, D.V., Malitch, K.N., Petrunin, A.G., 2009. Sibe- physics 199, 165e192. rian meimechites: origin and relation to flood basalts and kimberlites. Russian Zorin, Yu.A., Turutanov, E.Kh, Mordvinova, V.V., Kozhevnikov, V.M., Yanovskaya, T.B., Geology and Geophysics 50 (12), 999e1033. Treussov, A.V., 2003. The Baikal rift zone: the effect of mantle plumes on older Sobolev, S.V., Sobolev, A.V., Kuzmin, D.V., Krivolutskaya, N.A., Petrunin, A.G., structure. Tectonophysics 371, 153e173. Arndt, N.T., Radko, V.A., Vasiliev, Yu.R., 2011. Linking mantle plumes, large Zorin, Y.A., Turutanov, E.H., Kozhevnikov, V.M., Rasskazov, S.V., Ivanov, A.V., 2006. igneous provinces and environmental catastrophes. Nature 477, 312e318. The nature of Cenozoic upper mantle plumes in East Siberia (Russia) and Song, Y., Frey, F.A., Zhi, H., 1990. Isotopic characteristics of Hannuoba basalts, eastern Central Mongolia. Russian Geology and Geophysics 47 (10), 1046e1059. China: implications for their petrogenesis and the composition of subconti- nental mantle. Chemical Geology 85, 35e52. Tatsumi, Y., Koyaguchi, T., 1989. An absarokite from a phlogopite lherzolite source. e Contributions to Mineralogy and Petrology 102, 34 40. Chuvashova Irina: Senior Research Scientist at the Insti- Tatsumi, Y., Maruyama, S., Nohda, S., 1990. Mechanism of backarc opening in the tute of the Earth’s Crust SB RAS, Irkutsk, Russia (since e Japan Sea: role of asthenospheric injection. Tectonophysics 181, 299 306. 2004), assistant professor of Dynamic Geology Chair at the e Tsuruta, K., Takahashi, E., 1998. Melting study of an alkali basalt JB 1upto Irkutsk State University (since 2015), M.S. (2006) from the 12.5 GPa: behavior of potassium in the deep mantle. Physics of the Earth and Irkutsk State University, Ph.D. (2010) from the Institute of e Planetary Interiors 107, 119 130. the Earth’s Crust SB RAS. She worked in late Mesozoic and ’ Usoltseva, M., Khursevish, G., Rasskazov, S., Vorob eva, S., Chernyaeva, G., 2010. Cenozoic volcanic fields of Mongolia, China, Kyrgyzstan, Morphology of Actinocyclus and Lobodiscus species from the Miocene deposits Sakhalin, Amur region, and Siberia. Results of her research e of the Vitim Plateau, Russia. Plant Ecology and Evolution 143 (3), 352 364. in geochemistry and petrology of volcanic rocks have been http://dx.doi.org/10.5091/plecevo.2010.410. presented in 118 publications that include 6 monographs e fi Vasiliev, Y.R., Zolotukhin, V.V., 1995. The Maimecha Kotui alkaline-ultrama c and 18 peer-reviewed papers. province of the northern Siberian Platform, Russia. Episodes 18, 155e164. Villemant, B., Jaffrezic, H., Joron, J.L., Treuil, M., 1981. Distribution coefficients of major and trace elements; fractional crystallization in the alkali basalt series of Chaine des Puys (Massif Central, France). Geochimica et Cosmochimica Acta 45, 1997e2016. Volyanyuk, N.Y., Semenova, V.G., Elizarieva, T.I., 1976. Inclusions of pyrope and Rasskazov Sergei: Head of Laboratory for isotopic and pyropeespinel lherzolites from basanites of the Vitim plateau. Doklady of the geochronological studies at the Institute of the Earth’s USSR Academy of Sciences 228 (3), 693e696. Crust SB RAS, Irkutsk, Russia (since 1996), Head of Dy- Walter, M.J., 1998. Melting of garnet peridotite and the origin of komatiite and namic Geology Chair at the Irkutsk State University (since depleted lithosphere. Journal of Petrology 39 (1), 29e60. 2012), M.S. (1976) from the Irkutsk State University, Ph.D. Wang, Z., Gaetani, G.A., 2008. Partitioning of Ni between olivine and siliceous (1982), Dr. (1992), and Professor (2005) of Petrology and eclogite partial melt: experimental constraints on the mantle source of Ha- Volcanology from the Institute of the Earth’s Crust SB RAS. waiian basalts. Contributions to Mineralogy and Petrology 156 (5), 661e678. He has studied Cenozoic rifts in Central and East Asia, Wei, W., Xu, J., Zhao, D., Shi, Y., 2012. East Asia mantle tomography: new insight into plate North America, and East Africa. His main interest is origin subduction and intraplate volcanism. Journal of Asian Earth Sciences 60, 88e103. of volcanic rocks, erupted at the latest geodynamic stage of Yanovskaya, T.B., Kozhevnikov, V.M., 2003. 3D S-wave velocity pattern in the upper the evolving Earth. He has published more than 500 works mantle beneath the continent of Asia from Rayleigh wave data. Physics of the that include 17 monographs and 132 peer-reviewed Earth and Planetary Interiors 138, 263e278. papers. Yarmolyuk, V.V., Kovalenko, V.I., Bogatikov, O.A., 1990. South Baikal “hot spot” of the mantle and its role in the formation of the Baikal Rift. Doklady of the USSR Academy of Sciences 312 (1), 187e191. Yarmolyuk, V.V., Ivanov, V.G., Kovalenko, V.I., Samoilov, V.S., 1994. Dynamics of formation and magmatism in the Late Mesozoic-Cenozoic South Khangai Yasnygina Tatyana: Senior Research Scientist at the mantle hotspot (Mongolia). Geotektoniks 5, 28e45. Institute of the Earth’s Crust SB RAS, Irkutsk, Russia (since Yasnygina, T.A., Markova, M.E., Rasskazov, S.V., Pakhomova, N.N., 2015. Determi- 2004), M.S. (1993) from the Irkutsk State University, Ph.D. nation of rare earth elements, Y, Zr, Nb, Hf, Ta, and Th in geological reference (2004) from the Institute of the Earth’s Crust SB RAS. Her materials of the DV series by ICP-MS. Zavodskaya Laboratoriya. Diagnostika scientific interests are geochemistry of volcanic rocks and materialov 81 (2), 10e20. crude oils, trace-element melting modeling. The results of Zack, T., Foley, S.F., Jenner, G.A., 1997. A consistent partitioning coefficient set for her research have been presented in more than 110 publi- clinopyroxene, amphibole and garnet from laser ablation microprobe analyses cations that include 2 monographs and 29 peer-reviewed of garnet pyroxenite from Kakanui, New Zealand. Neues Jahrbuch fur Miner- papers. alogie. Abhandlungen 172 (1), 23e41. Zhang, Z., Mahoney, J.J., Mao, J., Wang, F., 2006. Geochemistry of picritic and associated basalt flows of the Western Emeishan Flood Basalt Province, China.

Please cite this article in press as: Chuvashova, I., et al., Mid-Miocene thermal impact on the lithosphere by sub-lithospheric convective mantle material: Transition from high- to moderate-Mg magmatism beneath Vitim Plateau, Siberia, Geoscience Frontiers (2016), http://dx.doi.org/ 10.1016/j.gsf.2016.05.011