Deep Structure of the Lithosphere in the Central Tien Shan along the Son-Kul Magnetotelluric Sounding Profile

New results of a detailed study of the deep structure of the Central Tien Shan along the Son-Kul magnetotelluric (MT) profile crossing the Son-Kul Lake are reported. Based on the results of magnetotelluric data modeling, the regional and local geoelectric anomalies in the lithosphere are studied and their quantitative characteristics are given. Geological interpretation of the geoelectric cross-section was carried out, which supported the existing ideas about the block–hierarchical structure of the upper part of the Earth’s crust. This corresponds to the tectonophysical concepts of the sequential inserted subordination of large and smaller elements of the zone–block structure consisting of stable blocks and limiting mobile zones, which are distinguished by the high dislocation of the geological substrate. The integral pattern of the distribution and morphology of zones of high electrical conductivity in this segment of the Central Tien Shan crust may reflect discretely localized palm tree–type structures associated with the evolution of transgressive suture zones of localized deformation during the Hercynian and Alpine tectogenesis.


INTRODUCTION
The study of the distribution of deep structuralmaterial heterogeneities, as well as their complex geological and geophysical characterization, is one of the most important problems of modern tectonics and geodynamics. The fact that magnetotelluric sounding (MTS) is most sensitive to conductive structures in the Earth's crust and upper mantle, most likely due to the presence of fluids, makes it one of the most successful deep geophysical methods in the study of fault structures in geodynamically active regions [1][2][3]. The Son-Kul Basin is a flat area with dimensions of ~29 × 57 km (Son-Kul Lake has a size of 19 × 25 km), its shape is close to an ellipse elongated in the latitudinal direction. It is located at an altitude of just over 3000 m and separated by ridges, the crests of which are raised to 3800 m and more. These ridges are the part of the Moldotau-Terskey system of uplifts. Geographically, the Son-Kul Basin is a link in an extremely extended zone of highly elevated and relatively narrow (valleylike) intermontane basins [4]. In the west, it is sepa-rated by sublatitudinal faults and linear grabens of the Minkush-Kokomeren Zone; in the east, the Ottuk-Kokturpak and further Archalu-Dzhinalach basins are observed in the relief (Fig. 1). Noteworthy is the anomalously large width of the Son-Kul Basin (~20 km), which gives it an almost isometric outline. This isometric shape, which sharply distinguishes the Son-Kul Basin from the grabens of the Minkush-Kokomeren Zone, is controlled by the absence of large faults on the sides of the basin, as well as by the younger age of the Son-Kul Basin, overprinted on the more ancient troughs of the Minkush-Kokomeren Zone. There are no Jurassic and Oligocene-Neogene sediments (Kyrgyz and Naryn formations) within the Son-Kul Basin, and only low-thickness Quaternary sediments occur there. Structurally and historically, the Son-Kul Basin is not an element of the Minkush-Kokomeren Zone of linear grabens; it is the younger overprinted basin. It was previously believed that a vast and deep granitoid batholith was the base of the basin. However, in light of later reconstructions, the Son-Kul granitoid massif is not considered as a batholith any longer, but is a rather thin plate, which almost does not differ from the host formations in its rheological properties. The influence of deep tectonic dislocations of the regional order in this area was considered as a likely reason for the presence of a large granitoid massif with anomalous values of apparent resistivity as  [3,5]. The magnetotelluric sounding performed in this study can provide the prerequisites for clarifying the answer to the following question: are the geoelectric anomalies discovered associated with deep tectonic faults or with fracturing of the Son-Kul granitoid massif itself. The most important tectonic structure of the Central Tien Shan, known as the Nikolaev Line, stretches in the immediate vicinity. It is very important in many structural, paleotectonic, geophysical, and other reconstructions, being considered as one of the fundamental deep-seated elements of the Tien Shan structure. At the same time, this largest structural boundary between the uneven-aged formations of the Northern and Central Tien Shan practically does not manifest itself in the neotectonic structure, while being well reflected in the results of deep electromagnetic studies. Another structural element that may appear in the geoelectric model is the so-called Kyrgyz-Terskey oceanic suture zone [6], which is observed from the western part of the Kyrgyz Range to the southeast, to the area of Lake Son-Kul and then to the east, to the upper reaches of the Naryn River. On the other hand, the deep structure of the studied region as a part of the Issyk-Kul microcontinent is of particular interest to many scientists as well [7][8][9]. It may be expected that the detailed magnetotelluric sounding performed will provide more definite data for the solution of these problems, which are of not only regional, but general scientific interest as well. Therefore, the MTS study was carried out along the Son-Kul profile, which was laid across the strike of the major tectonic structures of the Central Tien Shan on the territory of the Son-Kul Basin, by the Research Station of the Russian Academy of Sciences in Bishkek, with the support of the Russian Science Foundation, project no. 17-16-10059 ( Fig. 1).

METHODS OF STUDY
Our work was aimed at study of the deep geoelectric structure and its comparison with the geological structure in order to substantiate the structure of the lithosphere along the Son-Kul MTS profile. For this purpose, new detailed MTS was performed along a profile parallel to the 75th meridian ( Fig. 1) and crossing the major tectonic elements of the Northern Tien Shan: the Kyrgyz, Dzhumgal-Too, Son-Kul-Too, and Moldo-Too ranges, Son-Kul Basin, fault zones of different ranks, and the ophiolite zone confined to Son-Kul Lake. The total length of the profile was 90 km, and it included 30 sounding points. The study along the profile was carried out with different steps (from 1 to 15 km), depending on the field conditions. We used a Phoenix MTU-5A equipment in the range of periods from 0.001 to 1600 s. The measurements of the three components of the magnetic field were carried out using standard MTC-50 induction coils oriented to the north, east, and vertically. The components of the electric field were measured by dipoles with a length of 50 m, grounded by nonpolarizable electrodes with

Chuy Basin Chuy Basin Chuy Basin
T a l a s -F e r g a n a F a u l t T a l a s -F e r g a n a F a u l t T a l a s -F e r g a n a F a u l t Issyk-Kul Lake Issyk-Kul Lake Issyk-Kul Lake C e n t r a l T i e n S h a n C e n t r a l T i e n S h a n C e n t r a l T i e n S h a n S o u t h e r n T i e n S h a n S o u t h e r n T i e n S h a n S o u t h e r n T i e n S h a n manganese-carbon powder, designed and patented by the Research Station of the Russian Academy of Sciences. The data processing was performed in the onepoint mode only. Mathematical processing of the data of field soundings obtained along the Son-Kul profile was carried out using the standard software package SSMT-2000 (Phoenix). As a result, the frequency dependences of the components of the impedance tensor and the Wiese-Parkinson matrix were obtained in the range of periods from 0.001 to 1600 s.
Based on the results of MTS data processing, estimates of the impedance tensor, tippers, and horizontal magnetic tensor were obtained using the MT-Corrector software (worked out by the Nord-West Company, Russia), according to the standard research scheme for the Tien Shan region [10][11][12][13].

RESULTS OF THE STUDY AND ANALYSIS
The qualitative and quantitative interpretation of the data obtained was carried out by the methodology considered in [10][11][12][13]. The construction of the geoelectric model along the Son-Kul profile included the parameterization of the impedance, i.e., the determination of the inhomogeneity parameters N and the amplitude parameter of asymmetry Skew (Fig. 2). These constructions are necessary for regionalization of the area studied by to the "electrical dimension" of the environment and allow us to assess the degree of horizontal heterogeneity of the environment and localize the structures identified.
The inhomogeneity parameter N is determined by the formula where and are the principal values of the impedance tensor determined by the Eggers orthogonalization method [14]. It characterizes the degree of horizontal heterogeneity of the environment. If the environment is horizontally homogeneous (1D), then N = 0.
The amplitude parameter of asymmetry [14] skew is a measure of the asymmetry of the environment as well. In the case of a two-dimensional model, skew = 0.
The pseudo-sections skew and N characterize the degree of horizontal inhomogeneity of the environment; the skew value is 0.1 up to significant depths, which indicates the probability of using 2D inversion of MTS data. This is evident from the pseudo-section of the parameter N; this parameter does not have anomalous behavior (as is often the case). A two-dimensional smoothing inversion of MTS data along the Son-Kul profile (Fig. 3) was performed using the Rodi-Mackie software [15].
This program uses a nonlinear conjugate gradient method that attempts to minimize the objective function, which is the total of the normalized data residuals and model smoothness. The relationship between data discrepancies and model smoothness is controlled by the regularization parameter τ. The input data for the inversion are the longitudinal and transverse curves of the apparent resistivity and impedance phases (TEand TM-modes in the designations of this program), and tippers at 30 observation points in the period range from 0.01 to 1000 s. The mesh of the approximation model consists of 179 cells horizontally, sparse at the edges of the model and more frequent near sounding points, and 151 cells vertically, increasing in size with depth. The terrain relief is taken into account for detailed construction of the mesh (with a step of 50 m). The resistivity of the starting model is 100 Ω m. The error limits of the inverted data were set as follows: apparent resistivity modulus (TE mode), 100%; apparent resistivity modulus (TM mode), 10%; impedance phase (TE mode), 5%; impedance phase (TM mode), 5%; tipper, 0.05%. Based on several test calculations of the inversion, the value of the regularization parameter τ = 3 was selected. The model shown in Fig. 3 was obtained after 200 iterations. The value of the RMS residual (inversion convergence) was 2.54.

GEOLOGICAL INTERPRETATION OF THE GEOELECTRIC SECTION
The resulting geoelectric model is shown in Fig. 3, where the key heterogeneities of the cross-section are numbered for easier interpretation and the numbers in the text are indicated by numerals in parentheses. In general, the model is characterized by a contrasting distribution of objects with different resistivity. For the upper part of the cross-section, there is a good correlation between the high-resistivity blocks of the model and mountain ranges, such as Moldo-Too (4), Dzhumgal-Too, and Kyrgyz (8), with the exception of the Son-Kul-Too Range. The Nikolaev Line Zone (1) is the most contrasting subvertical geoelectric structure, which is characterized by an anomalously low resistivity (2-10 Ω m) with a depth of 45-50 km, and it is the border between the Northern and Central Tien Shan. The maximum width of the anomalous zone of high electrical conductivity reaches 5-7 km. The physical nature of the high conductivity can be associated both with the general fluid saturation of adjacent faults and with the syn-deformational processes of illitization of clays and hydration of feldspars in the clays of friction over Jurassic rocks in the neighboring   N   30  29  28  27  26  25  24  23  21  20  19  18  17  16  15  14  13  12  11  10  9  8  5  4  3  2  1  7  6  22   1000  630  400  205  160  100  60  40  25  16  10  6  4 2.5 90 km Minkush-Kokomeren suture zone detected in the upper part of the section [16].
The high-resistivity block (5) forming the central part of the cross-section may be composed of felsic igneous rocks (?) and is the axial part of the Issyk-Kul microcontinent [7,8,17]. Based on the results of observations of rock outcrops, the conducting horseshoe-shaped structure surrounding the high-resistivity block (2) is most likely composed of Precambrian mica-bearing granite-gneiss and metamorphic rocks, which are characterized by intense fracturing and, as a result, high fluid saturation. Such rocks compose the lower part of the Issyk-Kul microcontinent, the boundary of which passes into the lower crustal conductive layer. In addition, the Central Terskey fault zone appears in the geoelectric section near the area of point 30 (Fig. 3). The mosaic distribution of areas (objects) with high electrical conductivity (points 29 and 30) in the conjunction zone of the Kyrgyz Range and the Dzhumgal Basin may be explained by the fragmentary structure of the deep fault consisting of a number of separate backstage faults with a length up to 15 km. The zone of the high gradient of electrical resistivity at sounding points 22 and 23 on the geoelectric cross-section corresponds to the ophiolite suture [6], which includes fragments of the ancient oceanic crust.
Noteworthy are the two high-resistivity bodies located directly beneath Son-Kul Lake (6a) and at depths of 6-12 km (6b), with the values of apparent resistivity (AR) 1000-1600 Ω m and separated by a horizontal electrically conductive zone (AR = 160-200 Ω m) at a depth of ~5 km (6c), which may be due to intense fracturing of rocks. The high-resistivity objects (6a) and (6b) may have previously composed a single body, a large intrusive massif comparable in size to a batholith, the depth of which reaches 15 km and the thickness of which is ~6 km.
The subhorizontal zone of high electrical conductivity (7) may be due to tectonic disintegration of rocks and their high fluid saturation; this is in a good agreement with the position of the upper crustal conductive layer, sporadically distributed throughout the entire Central Tien Shan [18].
In general, the integral pattern of the distribution and morphology of zones of high electrical conductivity in this crustal segment of the Central and Northern Tien Shan may reflect, in our opinion, a discretely localized manifestation of palm tree-type structures associated with the evolution of transpressive suture zones of local deformation during the Hercynian and Alpine tectogenesis [19]. As a rule, the latter are combined with compensatory structures of transtension in the form of pull-apart intermontane basins spatially and hierarchically [5,20].

CONCLUSIONS
Thus, a comparative analysis of geological and geophysical data on key elements of the geoelectric model developed along the Son-Kul MTS profile was carried out. Due to the significantly more detailed and accurate MTS data, the boundaries of the Issyk-Kul microcontinent corresponding to modern tectonic concepts were obtained in the conductive structures of the constructed geoelectric cross-section at the stage of qualitative interpretation and the first quantitative estimates. The fault zones are clearly identified in the geoelectric cross-section in the form of subvertical conductive heterogeneities with different slopes of the interpreted deep structures. According to the currently available geological and geophysical data (the results of study on earthquake sources, the results of GPS measurements, structural-geological studies) for the Central Tien Shan, the zone considered occurs under the conditions of general horizontal submeridional compression. Thrust faults formed in horizontal compression settings can subside both northward and southward. This conclusion is supported by the crosssection constructed. Most likely, geoelectric anomalies (highly conductive structures) have different vergence and shape, similar to those of palm tree-type transpressive suture structures. The distribution of high-resistivity objects in the upper part of the crosssection (to depths of ~25 km) supports the existing ideas about the block-hierarchical structure of the upper part of the Earth's crust in the Central Tien Shan, and the size of the blocks reflects the degree of discreteness of its tectonic processing. Such representation for the upper part of the Earth's crust is in agreement with the tectonophysical concepts of a system of successively enclosed smaller elements in larger stable ones (a zone-block structure), as well as the presence of mobile zones limiting them. It is accepted that such mobile zones differ from stable blocks by the high dislocation of the geological substrate.