Geology of the Gobi Altai and Tseel terranes in the central part of the Sagsai River Watershed, SE Mongolian Altai

ABSTRACT A geological map is an indispensable tool for understanding the structure of the Earth’s crust but high-quality geological maps are usually lacking in remote areas of mountainous Central Asia covered by vast deserts. The progress in remote sensing and geographical information system (GIS), as well as the advancement in analytical methods, have generated new challenges in producing modern geological maps in such regions. The presented 1:50,000 geological map along the Sagsai River summarizes new and more accurate geological data from the geologically interesting region at the contact of the supracrustal and deep crustal Tugrug and Tseel units forming the metamorphosed accretionary wedge on the S and SW slopes of the Mongolian Altai. These geological units are formed by the lower Palaeozoic volcano-sedimentary sequences affected by Devonian Barrovian metamorphism ranging from lower greenschist to granulite facies. This metamorphic basement was subsequently intruded by the post-orogenic late Carboniferous Sagsai Pluton. The presented map shows complex relationships between different crustal levels of the metamorphosed accretionary wedge and post-tectonic intrusion evaluated using a GIS, numerical processing of remote sensing data as well as field documentation and laboratory studies.


Introduction
A geological map is a scaled-down interpretation of the structure of a selected area of the upermost layer of the Earth's crust, usually drawn on a topographic base. Geological maps can be used solely for academic purposes but are also essential for economic development of individual countries. This is true in Central Asia, and especially in Mongolia, where the economy is based on mining and raw material exploration.
By 2003, Mongolia has been covered by geological maps at a scale of 1:200,000, while geological mapping at the scale 1:50,000 has been under way since 1990s. Even though there are available geological maps produced by foreign geological surveys (e.g. Japan, France, South Korea and the Czech Republic) only 40% of Mongolian territory is covered by maps at this scale and the maps are of variable quality.
As southern Mongolia is poorly covered by vegetation, modern geographical information system (GIS) methods combining the interpretation of remote sensing data with field reconnaissance and up-to-date analytical data represent an effective approach for the compilation of new geological maps. We present here an example of such geological map , which covers the drainage area of the Sagsai River between the Tseel and Tsogt districts in the SE Mongolian Altai. This Main Map demonstrates the efficiency of combined satellite imagery and field investigations in producing reliable low-cost geological maps in regions of high economic interest.

Map construction
The topographic base of the map was derived from ASTER GDEM V2 data with a resolution of 30 m. The shaded relief was prepared using the 3D Analyst tool of ArcGIS in 4 layers with azimuths of 60°, 150°, 20°and 330°and a constant Sun angle of 30°. The course of the rivers was digitized from Landsat 8 imagery and roads were outlined using a global positioning system (GPS) receiver. Contour lines with 25 m interval were generated using the Contour tool in ArcGIS based on ASTER GDEM V2 data. False contours and apparent mistakes were manually deleted and the remaining contours digitally smoothed.
Geological boundaries were interpreted from the Landsat 8 scenes p137-138r28-29. These satellite images, free of cloud cover for the year 2015, were downloaded from the USGS website and processed using image enhancement techniques with GRASS GIS and Python. Principal component analysis and factor loading analysis with an emphasis on geological information showed that band combinations (6,7,4) had the highest contrast that could be further enhanced with a decorrelation stretch. Band ratio combinations (4/3, 6/2, 7/4) had good contrast for geological information and formed the input data for visual classification.
The geological interpretation of Landsat 8 imageries was supported by field documentation on 212 reference points complemented by nearly 60 thin sections and 21 whole-rock geochemical analyses.
Lithostratigraphic columns were constructed according to the observed field relationships. The chronological assignment of geological events was done according the published data cited in the text and supplemented by unpublished radiometric data of the Czech Geological Survey.
Main Map (map sheet No. L-47-85-G) covers the boundary between the volcano-sedimentary sequences of the Gobi Altai Terrane and metamorphic rocks of the Tseel Terrane ( Figure 1) according to regional nomenclature of Badarch et al. (2002). The boundary between these units is lithological in the West and East, and sealed in the South by intrusion of the late Carboniferous mafic Sagsai Pluton (Buriánek et al., 2016). Recent study by Hanžl et al. (2016) suggests that these two units represent varying degrees of metamorphism of an early Palaeozoic accretionary wedge.

Tseel Terrane
The Tseel Terrane (Bibikova et al., 1992) forms an E-W trending narrow belt along the southern slopes of the Mongolian Altai. It is composed of lower Palaeozoic volcano-sedimentary rocks with Devonian to late Carboniferous polymetamorphic history (e.g. Burenjargal et al., 2014;Jiang et al., 2012;Kozakov et al., 2007) intruded by the Devonian to Triassic plutons (Burenjargal et al., 2016;Buriánek et al., 2016;Cai et al., 2015). Two parts of the Tseel Terrane are exposed on the map: the Tseel Metamorphic Complex in the SW and the Tsogt Metamorphic Complex in the SE.

Tseel Metamorphic Complex
The Tseel Metamorphic Complex is composed of phyllites to garnet-bearing micaschists which alternate with layers of metarhyolites and greenschists to amphibolites, and two orthogneiss belts of Devonian protolith age . The southern orthogneiss belt is formed by lithologically homogeneous coarsegrained orthogneiss, while the northern one consists of fine-to medium-grained orthogneisses containing layers of micaschists and amphibolites. Dominant structure is a NW-SE trending metamorphic foliation (S 1 ) steeply dipping to the NE ( Figure 2). Two generations of folds and lineations are preserved; the earlier linear elements are steeply plunging to the N, while the younger ones are sub-horizontal and trending NW-SE.

Tsogt Metamorphic Complex
The Tsogt Metamorphic Complex is characterized by the presence of high-grade metamorphic rocks represented by amphibole-biotite migmatites, migmatitized amphibolites and layered metagabbros. These rocks contain lenses of orthopyroxene granulite yielding a peak metamorphic age of 384±2 Ma (Kozakov, Glebovitsky, Bibikova, Azimov, & Kirnozova, 2002). The metamorphic sequence is intruded by foliationparallel sheets of deformed granite and tonalite. The main metamorphic fabric S 1 defines a dome-like structure as shown by poles of foliation concentrated along a NNE-SSW trending girdle in the stereonet (Figure 2).

Tugrug Formation
In the study area, the Gobi Altai Terrane is represented by the Tugrug Formation, which was defined by Markova (1975) on the southern slopes of the Mongolian Altai. In the northern part of the region, the Tugrug Formation is formed by slates to sericite-chlorite phyllites with carbonaceous admixture and layers of tuffitic sandstones. Towards the south prevail epidote-chlorite greenschists alternating with tuffitic sandstone layers with locally preserved flattened pillow lavas. The southernmost part of the Sagsai River section is characterized by metagreywackes and metasiltstones.
Volcanic rocks of the Tugrug Formation are represented by subalkaline series forming a basaltbasaltic andesitedaciterhyolite compositional trend ( Figure 3). The Silurian stratigraphic age of the Tugrug Formation was assumed by analogy to similar sequences in the Shine-Jinst district to the east (Markova, 1975). However, Cambro-Ordovician age was proposed for this slightly metamorphosed volcanosedimentary sequence in the geological maps of Borzakovskii (1985) and Rauzer et al. (1987). As there are no palaeontological data from the study area we assume Cambrian to Ordovician age for the Tugrug Formation based on detrital zircon age spectra from the paragneisses of the Mongolian Altai further west (Jiang et al., 2012) and the lithologically identical Habahe Group rocks from the Chinese Altai (Jiang et al., 2011). This proposition is further supported by the unpublished clastic zircon U-Pb age spectra from two samples collected from the Tugrug Formation north of Tsogt village. The youngest peaks at 506 and 475 Ma are interpreted to reflect the sedimentation of the Tugrug Formation volcano-sedimentary detritus (Y. Jiang, pers. communication).
The volcano-sedimentary rocks of the Tugrug Formation in the Sagsai River valley are cross-cut by a dacite dyke, which yields a U-Pb zircon age of 380 Ma (LA ICP-MSunpublished data on zircon of the Czech Geological Survey). A similar Devonian age (∼380 Ma) was reported for the slightly deformed rhyolite body NE of the Sagsai Pluton intruding the Tugrug Formation (Y. Jiang, pers. communication). Both dykes and rhyolites are chronologically and geochemically identical to syn-tectonic felsic orthogneisses of the structurally deeper Tseel Metamorphic Complex implying genetic link between these magmatic rocks ).
The regional metamorphic grade of the Tugrug Formation increases to the SW where metagreywackes gradually pass into garnet-biotite schists at the boundary with the underlying Tseel complex. The boundary between the Tseel Metamorphic Complex and the Tugrug Formation is thus arbitrary and may be defined by lithological differences and metamorphic grade. In the present map it is placed between the variegated staurolite bearing schists and orthogneisses of the Tseel Metamorphic Complex and greenschist-facies metagreywackes of the Tugrug Formation W of the Sagsai Pluton. In contrast, the NE-SW trending boundary of the low-grade Tugrug Formation with the underlying high-grade Tsogt Metamorphic Complex is defined by the sharp metamorphic and lithological gradient.
Approaching the Sagsai Pluton, the low-grade greywackes and shales are transformed into biotite-chlorite and cordierite-biotite schists indicating increasing thermal metamorphic effects of the pluton. Septa inside the pluton are metamorphosed to cordierite-biotite schists and garnet-bearing sillimanite-biotite migmatitized paragneisses.
The Tugrug Formation rocks reveal well preserved bedding (S 0 ) of variable orientation (Figure 2). The bedding is transposed in many places to composite S 1 foliation trending dominantly WNW-ESE. This pervasive fabric is refolded by open to closed, upright folds with steep NW-SE trending axial planar cleavage (S 2 ).

Zalaa-Oroin Formation
Lower Carboniferous turbidites with prevailing siltstones of the Zalaa-Oroin Formation (Bezzubtsev, Zaporotets, & Lysov, 1959) unconformably cover paragneisses and schists at the NW corner of the Main Map. Carboniferous greywackes and sandstones are unmetamorphosed to weakly metamorphosed. The S 0 bedding is folded by open upright NW-SE F 2 folds and locally reworked by pronounced fracture to slaty S 2 cleavage.

Sagsai Pluton
The Sagsai Pluton intruded the southern part of the Tugrug Formation at its contact with the Tseel Terrane. The body consists of monzogabbros to quartz diorites and biotite to amphibole-biotite syenites to quartz monzonites (Figure 4). Mafic rocks with common mingling textures were intruded by dykes of two-mica granites arranged in a ring like structure along the rim of the pluton. U-Pb zircon dating of monzonite yielded ages of 307±2 and 309±2 Ma,  (Pearce, 1996). Unpublished data of the Czech Geological Survey. while the gabbro revealed an older age of 322±2 Ma (Buriánek et al., 2016).
Orientation of primary magmatic layering is scattered all around the pluton and is clearly discordant to NE-SW trending steep S 1 and S 2 foliations inside and outside the metamorphic aureole.

Cenozoic cover
The area has a denuded nature with insignificant Cenozoic cover. Recent fluvial sediments have evolved in the valley of the Sagsai River and its tributaries with intermittent streams. They are composed of gravels and sandy gravels with clasts up to 50 cm in size and also of blocks occurring close to the nearby slopes. Colluvial sediments are of minor areal extent. They are usually characteristic of debris fans under steeper slopes or silty deposits in flat depressions. Two small gold placers were found in these sediments and weathered rocks of the Tugrug Formation in the NW quadrant of the map. Sandy to silty aeolian deposits were only seen in one dune below the western rocky slope of the lower part of the Sagsai River valley.

Conclusions
The presented geological map sheet No. L-47-85-G summarizes recent Landsat 8 imagery, geological field and analytical data from the region at the southern slopes of the Mongolian Altai. In comparison with previous geological maps it shows a: . more precise subdivision of lithological and lithostratigraphic units; . chronostratigraphic calibration of the main lithostratigraphic units; . increase of the metamorphic grade in the Tugrug Formation from NE to SW towards the contact with the Tseel Metamorphic Complex; . similar pattern of regional metamorphic foliation trends of all units forming the Tseel and Gobi Altai terranes; . efficiency of using free Landsat imagery in combination with precise field lithological and structural data as well as modern geochronological and whole-rock geochemical data.

Software
The geological map was constructed using ArcGIS 10.2.1 and improved in Adobe Illustrator CS 11. Landsat 8 data were processed using Grass GIS and Python. The geochemical data were interpreted using the freeware GCDkit (Janoušek, Farrow, & Erban, 2006) and structural data by the SpheriStat software.