. Introduction
The evolution of paleoseismogeological studies clearly demonstrates that in order to properly understand the seismic potential of a region, and to assess the associated seismic hazard, extensive studies are necessary to take full advantage from the geological evidence of past earthquakes. The period of instrumental seismological observations is insignificant in comparison with the recurrence interval of strong earthquakes. Thus to achieve these goals the historical data are being involved. Paleoseismogeology supplements historical and instrumental records of seismicity by characterizing strong prehistoric earthquakes. It is focused on studying ground effects (both primary — ruptures, and secondary — gravitational deformations (Solonenko, 1973 and McCalpin, 2009)) from past earthquakes preserved in the geologic and geomorphic environment. In this context the recurrence interval of strong paleoearthquakes is a critical characteristic. It assumes applying numerical dating techniques for absolute age determination of seismically induced land-forms, associated sediments or other evidences of strong past seismic events.
The most direct methods (fault scarp formation applying scarp degradation modeling, cosmogenic surface exposure dating, analysis of scarp soils etc.) return relatively large age uncertainties (McCalpin, 2009). Dating landforms is typically less precise than dating sediments (Burbank and Anderson, 2001 and McCalpin, 2009), although relatively precise indirect dating of associated deposits also may not provide close constraints of ancient seismic events. Today radiocarbon analysis is the most exploitable and widely applicable technique in paleoseismogeology. It allows determining the age of soils, tree fragments and other organic materials which were deformed and/or buried during the earthquake.
Generally, the application of the radiocarbon method is limited by ~50 ka but some problems occurred while using this technique for dating recent (about 2 ka) seismic events due to high relative methodological error and presence of several “plateaus” at the calibration curve which make it difficult the age correlation with the calendar scale (Wagner, 1998). In case of spreading forest vegetation in seismically active areas tree-ring analysis can be used as an additional or alternative approach. It has a great potential due to utmost precision of dating — annual and sometimes even seasonal resolution, which is an obvious advantage in comparison with the radiocarbon analysis.
Seismically induced surface processes can affect vital activity of trees and cause their reaction on environmental changes. Thus rings of trees — witnesses of seismic events, could be potentially valuable nature archive containing both the evidence for the earthquake and its age. This forms the basis for precise dating seismically induced geomorphic processes, landforms and sediments produced or deformed by earthquakes (Jacoby, 1997). Tree rings can record evidences of earthquakes associated with tectonic deformations or ruptures (primary effects), seismogravitational deformations or mass movement processes (secondary effects), and various other seismically induced environmental changes (side effects) (Table 1). Some most exploitable in paleoseismogeological investigations aspects of tree-ring analysis are discussed in literature (Jacoby, 1997; McCalpin, 2009 and others).
Table 1
For dating ruptures and tectonic scarps formation analysis of splitting tree-trunks and tilting stems could be applied (Page, 1970; Ruzhich et al., 1982 and Bekker, 2004).
Regarding the earthquake induced mass movement processes the germination age of forest vegetation growing on the bare surface (estimated by the age of the eldest tree) returns the upper possible time of landslide formation. The direct dates of landslide events can be obtained by analyzing buried dead trees and/or reaction wood formation from the living ones as a response on the stem inclination. Subsequent growth of the tilted tree will aim to restore the trunk’s vertical position. Thus eccentric tree-ring growth after a tilting event can be dated (Sheppard and Jacoby, 1989 and Stoffel and Corona, 2014). Wood penetrating injuries and scars are also a common feature in trees affected by earthquake induced rock-falls. Callus tissue formation is commonly regarded as a reliable indicator of past geomorphic process activity (Lundström et al., 2009 and Trappmann and Stoffel, 2013). Partial stem burial with debris leads to a reduced activity of the roots, and to additional mechanical pressure of slope material (Hupp et al., 1987 and Friedman et al., 2005). The resulted growth suppression is indicated by reduction in the annual ring widths and therefore can be dated (Stoffel and Corona, 2014). Exceptionally the growth increase is possible in case of rich nutrition and water supply (Strunk, 1997). Roots exposure induced by slope mass movements leads to shortage of water and nutrient supply and, therefore, to suppressed tree growth and the formation of narrow rings in the stem (McAuliffe et al., 2006). In case of decapitation or the removal of branches, as a result of seismically triggered rock-falls, trees react with distinct radial growth suppression.
Various side effects of strong earthquakes such as ground shaking, hydrological changes etc. can also affect the forest vegetation. The negative impact on their living conditions can stimulate growth suppression, missed rings, and even mortality of trees. Sometimes the elimination of neighboring trees can result in a new less competitive environment for survivor ones. The survival trees can respond on these favorable changes with growth increase and wider tree rings (Strunk, 1997).
This paper presents the results of paleoseismogeological investigations in the pleistoseist zone of the 2003 Chuya earthquake (MS=7.3) that stroke southeastern part of Russian Altai. During this strongest historical seismic event the southern fault boundary between Kurai-Chuya system of intermountain depressions and framing ridges, as well as fault boundaries of the Chagan-Uzun massif (an independent neotectonic block, which separates these depressions) was reactivated. The studied sites are located in the northern and southern bounders of this massif, which is one of the major seismogenic structure in the region. We report new radiocarbon dates related to the geological evidences, which define the chronology of strong prehistoric earthquakes in the region. These data together with the results of tree-ring analysis and previously published radiocarbon dates allow specifying the recurrence interval of strong earthquakes for the last 4000 years.
. Study area
Geological settings and evidences of high regional seismicity
The Altai Mountains are the northern part of the Central Asia collision belt (Fig. 1). They stretch northwest more than 1500 km across the borders of Mongolia, China, Kazakhstan and Russia, and form a wedge shape narrowest in the southeast and widest in the northwest. Elevation increases in the opposite direction from 400 to more than 4000 m a.s.l. The neotectonic structure of the Altai Cenozoic uplift is formed by four NW trending dextral strike-slip fault zones. These major regional faults, together with feathering ones, divide Earth crust into separate tectonic blocks — ridges and intermountain depressions (Novikov, 2004).
The high-mountain southeastern part of the Russian Altai (SE Altai) is the northern extension of Mongolian Altai which is known by its high seismicity. Within the neotectonic structure of the Altai Mountains, SE Altai is a transpressional zone formed due to oblique thrusting. Significant lithosphere destruction and block movement along the sub-latitudinal reverse faults and thrusts take place here.
The SE Altai is the highest part of Russian Altai. It includes the Chuya and Kurai intermountain depressions and framing ridges. Bottoms of the Chuya and Kurai intermountain depressions are located at about 1750–1900 and 1500–1600 m a.s.l. correspondingly, with altitudes of framing ridges up to 3500–4200 m a.s.l. These depressions are separated by the Chagan-Uzun massif with the highest Sukor peak (2926 m a.s.l.) in its north-western part. Developing presently as a single structure the Chagan-Uzun massif together with the western part of the Chuya intermountain depression is characterized by the largest rates of vertical movement 7.84–8.00 mm a–1 within the whole Chuya river basin (Rogozhin and Platonova, 2002). Numerous large seismically induced paleolandslides and rock-falls, that mark fault boundaries of the Chagan-Uzun massif, as well as smaller blocks of its inner structure, provide evidence for its high neotectonic activity and make it one of the major seismogenic structures in the region (Rogozhin and Platonova, 2002; Emanov et al., 2009 and Agatova et al., 2014b).
The high seismicity of neighboring Mongolian Altai (Zhalkovsky and Muchnaya, 1975 and Rogozhin et al., 2013) was one of the reasons to mark the SE Altai as the most seismically active part of the Russian Altai. Almost its whole territory is characterized by seismic intensity of VIII point according to MSK-64 (Medvedev-Sponheuer-Karnik macroseismic intensity scale). The intensity IX zone was restricted to a small area including the south-western part of the Chuya depression (Zhalkovsky et al., 1995). It was outlined proceeding from the seismic origin of numerous landslides in the Chagan-Uzun and Dzhazater river basins and from the high seismicity of neighboring areas of Mongolia, shocked by several great seismic events in the 20th century (Devyatkin, 1965).
High seismicity of the SE Altai is evidenced by paleoruptures, numerous large Holocene earthquake induced landslides, and seismic convolutions in soft sediments (Fig. 2) (Butvilovsky, 1993; Rogozhin and Platonova, 2002; Rogozhin et al., 2007; Deev et al., 2009 and Agatova et al., 2006, 2014a). Recently it was supported by the 2003 Chuya earthquake (Ms=7.3). This strongest historical seismic event generated a complicated NW trending system of primary seismotectonic deformations of about 70 km in length. Main rupture is a dextral strike slip fault with a horizontal displacement up to 2 m, and subsidiary branch faults of reverse and dextral strike slip kinematics (Rogozhin et al., 2007). Landslide on the right Taldura river bank is the largest seismogravitational deformation triggered by this earthquake (Fig. 3).
Widely distributed unconsolidated Cenozoic sediments provide favorable conditions to produce earthquake induced landslides. Rockslides are mainly distributed within mountain ranges, framing Chuya and Kurai intermountain depressions. In contrast, the vast extent of the Neogene lacustrine and Pleistocene moraine deposits in the most active areas of the depression-range transition within the Chuya basin provides geological conditions for developing deep seated landslides.
It should be mentioned that generally all types of landslides triggered by earthquakes may also occur without seismic triggering (Keefer, 2002). Strong earthquakes are the leading factor in the generation of giant landslides in the SE Altai. Our field observations of the impact of strong modern and prehistoric earthquakes reveal several criteria that indicate seismic origin of observed ground failures within studied area (Nepop and Agatova, 2008): i) the great extent of ground failures upon strongly cemented Neogene and Pleistocene permafrost sediments, as well as high concentration of numerous seismically induced landforms (both huge and small ones) within the same area; ii) location of the landslides in the most active areas at the basin-range transition along fault boundaries of landforms and fault crossings; iii) absence of correlation between location of landslides and exposure of slopes to insolation, their age and rock composition; and iv) conjunction (relative position) of landslide locations with the ruptures on the watersheds and valley slopes. It should be also noted that location of all giant landslides in an area of extremely dry climate rules out their rainfall or a snowmelt trigger.
Regional climate and potential of applying tree-ring analysis in the SE Altai
Contemporary climate in the region is determined by its intracontinental position, with the main moisture transfer from the west (Atlantic Ocean), and to a lesser degree from the north (Arctic), with a dominant influence of the Mongolian anticyclone, giving rise to increasing aridity southeastwards and complicated latitude-longitude orographic climatic zoning. The SE Altai is characterized by extreme ultra continental cryoarid permafrost affected environment. The mean annual precipitation within the SE Altai is less than 200 mm in the floor of intermountain depressions, increasing with height. The mean annual precipitation near the snow line decreases along the W-E axis from 2000 mm (in the western part of the Katun range) to less than 500 mm in the Chikhachev range near the Mongolian border (Narozny and Osipov, 1999).
In the western part the stony steppe on the floor of the Kurai depression grades into taiga vegetation on the range’s slopes. In the western mountains flanking the Chuya depression, forests have an insular distribution, vanishing completely in the southeastern part. The plateau-toppled highlands are covered by alpine meadows. The higher topographic belt is represented by alpine land-scapes with mountain tundra and tundra-steppe vegetation changing with height into a glacial zone.
Under arid conditions, local orographic climate factors, first of all, the slope exposure, mainly control tree ring characteristics and tree growth near the upper timber limit in the region. There is a good correlation among growth indexes of all obtained regional tree-ring chronologies (TRC), located above 2200 m a.s.l., as well as a good correlation with tree-ring chronologies for neighboring areas of Tuva (TRC “Mongun”) and Mongolia (TRC “Khalzam Khamar”) (Tainik et al., 2015). Radial growth of trees in the region is affected by early summer (June or July) temperatures (Tainik et al., 2015; Barinov et al., 2016). Today a number of TRC are developed for the SE Altai: 1105-year TRC “Aktru” (Ovchinnikov et al., 2002), 1896-year TRC “Jelo” (Myglan et al., 2009) (both on Larix sibirica Ledeb); 1241-year TRC on Pinus sibirica Du Tour (Nazarov and Myglan, 2012). All these chronologies characterize environmental conditions of tree growth within the northern slope of the North Chuya range.
Generally, distribution of forest vegetation and the age of trees are the main limitations of applying tree-ring analysis in paleoseismogeological investigations. This technique is a promising tool for dating seismically induced landforms within the SE Altai due to several reasons: i) seismically induced slope processes occur in the forest stand zone in immediate vicinity with the modern upper timber limit (at this altitude the tree-ring growth is mainly controlled by climate which is favorable for building long TRC); ii) arid climate promotes good preservation of wood (up to two thousand years on stone surfaces); iii) a number of long absolute TRC has been developed for the SE Altai and adjusted Tuva region; iv) the 2003 Chuya earthquake gives unique opportunity to analyze the distribution patterns of wood penetrating injuries and other growth anomalies of specifically seismic origin.
History of paleoseismogeological investigations
First historic accounts of earthquakes within the Russian Altai are dated back to the 18th century and the first Earthquake Catalog of the Russian Empire, which mentioned Altai seismic events, was issued in the 19th century. The earliest instrumental seismic data had come from remote stations that provided reliable records of quite strong (M >5) earthquakes (Kirnos et al., 1961) before a regional seismological network was set up in 1962 (Zhalkovsky et al., 1995). Numerous instrumental records of K=7–14 seismic events (corresponding M~1.66–5.5) collected by the regional network allowed calculating location of earthquake epicenters, estimating their magnitudes, and making seismic zoning (Zhalkovsky and Muchnaya, 1978). Applying these results together with geological and paleoseismological data the recurrence interval 1–3 ka of strong Holocene earthquakes for the SE Altai was suggested (Zhalkovsky and Muchnaya, 1978). The idea of high regional seismicity of the SE Altai later was supported by the results of seismotectonic investigations (Reisner and Ioganson, 1996). The first extensive paleoseismogeological researches including numerous radiocarbon dating within the SE Altai were carried out before the 2003 Chuya earthquake (Rogozhin and Platonova, 2002), when a number of previously unknown strong past seismic events was reported.
For the first time in paleoseismogeological researches within the SE Altai Rogozhin introduced routine application of radiocarbon analysis for numerical age estimations of seismically deformed/cut paleosoils, organic material buried in colluvium wedges and pressure ridges, and for dating other effects of past earthquakes (Rogozhin and Platonova, 2002 and Rogozhin et al., 2007, 2008). A broad set of obtained radiocarbon dates allowed specifying the recurrence interval of strong earthquakes in the region — 500–900 years for the last 9 ka (Rogozhin et al., 2007, 2008).
In our paleoseismogeological investigations within the SE Altai we complete existing radiocarbon dataset and for the first time utilize dendroseismological approach (Agatova et al., 2014b).
. Methods
Detailed geomorphological and paleoseismogeological investigations and process analyses were based on interpretation of aerial photographs, Landsat-TM images, topographic maps (scale of 1:25000), and field investigations including mapping of landforms and deposits of different genesis.
Radiocarbon dating was done in the Cenozoic Geochronology Center SB RAS, Novosibirsk, and in the Institute of Geography RAS, Moscow. Determination of carbon residual activity in both laboratories was done with the QUANTULUS-1220. The conventional radiocarbon ages were calibrated applying OxCal v.4.3.2 program (Bronk Ramsey, 2009) with the IntCal13 calibration data set (Reimer et al., 2013). In the text calibrated results are reported as intervals at 95.4% significance levels.
For making dendroseismological investigations local TRC was developed. Sampling was carried out with the main goal of building the longest possible TRC and obtaining the representative sample collection of wood penetrating injuries. Samples were taken from trees (both living and dead ones) located on the surfaces of the talus fans and landslides near the scarps or most active talus channels. Trees with traumas on the slope facing and lateral sides of the trunks (affected by frontal and tangent impact of falling stones) were selected among the standing ones. Cores were taken from living trees and discs – from dead ones. In order to provide reliable dating additional discs were collected from the uninjured parts of tree trunks as well as cores from undamaged living trees.
Preliminary sample preparation and dendrochronological analysis was carried out in the Laboratory of natural science methods, Siberian Federal University, Krasnoyarsk. The standard procedure of building the tree ring chronology (Shiyatov et al., 2000) on Pinus Sibirica Du Tour was applied. Measurements of radial growth for each sample were made on two and more radii applying LINTAB 5 equipment with an accuracy of 0.01 mm. Tree ring series were dated by combining cross-correlation analysis (in the special DPL software (Holmes, 1983)) and graphical cross dating in the TSAP system V3.5 software (Rinn, 1996) with a visual control of variability curves. The age trend was removed using standard procedure (Fritts, 1976) with a negative exponent and spline (2/3 of the length of each series) in the ARSTAN program (Cook and Krusic, 2008). Tree ring width measurements of several samples collected from the same tree were brought into generalized series. Thus finally every tree was represented by a single tree-ring series.
The germination ages of trees growing on the bare surfaces of earthquake triggered landslides were calculated, applying age correction for colonization time gap and time of surface stabilization. Besides estimating the germination ages and assessing the minimal possible age of landforms, tree-ring analyses was also applied for estimating the age of earthquakes occurred already after forest regeneration on the surface of these landslides. We suppose that for trees (both dead and living ones) located near the scarps and talus fans, wood penetrating injures could be caused by earthquake induced rock falls as it was during the 2003 Chuya earthquake. Developing local tree ring chronology gives the opportunity to date such events.
The suggested approach implies analysis of wood penetrating injuries in the individual tree ring series for identification and dating of earthquake triggered paleorock-falls. As well as the number of such anomalies, the simultaneity of injuries sustained by several trees growing on different earthquake induced landslides was taken as a criterion of their seismic origin. The accuracy of this approach was supported by analysis of wood penetrating injuries caused by the 2003 Chuya earthquake induced rock-falls. Growth anomalies (disturbance of the tree ring structure) were visually determined with the stereomicroscope and dated.
It should be noted that besides earthquake triggered rock-falls there are climatically driven ones when heavy rains or snow melting act as a trigger. Thus determination the seismic origin of rock-falls plays a crucial role. We suggest that, similarly to the 2003 Chuya earthquake, the number and simultaneity of wood penetrating injuries sustained by many trees grown in different parts of an area affected by seismic shock could be such criteria. It should be emphasized that even after that, the dendrochronologically obtained dates of abrupt intensifications of rock-falls are just supposed dates of earthquakes, which should be verified by alternative proxy data.
. Results
Our paleoseismogeological researches were focused on the pleistoseist zone of the 2003 Chuya earthquake. Newly arisen and renewed scarps and ruptures gave a unique opportunity to study previously unknown evidences of strong prehistoric seismic events. Witnesses of several events were found in sediments of the largest landslide triggered by the 2003 Chuya earthquake and in Taldura river valley (southward of the Chagan-Uzun massif) (Fig. 4 – T1).
Dendroseismological investigations were conducted within two largest landslides (located at 2250–2200 and 2170–2130 m a.s.l. in the northern part of the Chagan-Uzun massif (Fig. 4 – A)) from the complex of earthquake triggered landslides in the Arydjan valley (Fig. 5) (Agatova et al., 2014b). A reactivated ancient fault, which laid the foundation of the head of Arydjan valley, formed main scarps of these landslides.
High concentration of numerous seismically induced landforms (both large and small ones) within the same area; absence of correlation between location of landslides and exposure of slopes to insolation, as well as between age of deposits and their rock composition; correspondence of landslides to fault boundaries of land-forms and fault crossings; conjunction of landslide locations with the ruptures on the watersheds and valley slopes argue for the seismic origin of the studied landslides.
Siberian pine (Pinus sibirica Du Tour) is the main forest forming tree species at this altitude in the Arydjan valley. Single trees settled the surface of taluses and rock-streams leaned on landslides. The upper timber limit in this valley extends slightly above the landslides but the scarp of the upstream landslide with the height of about 150–170 m is not settled by trees due to slope steepness and active taluses. More gentle main scarp of the older downstream landslide is settled by single trees.
As a result of precedes geomorphological study, collecting and culling the samples, altogether 61 discs (49 from the downstream landslide and 12 from the upstream one) from the 36 dead trees (30 from the downstream landslide and 6 from the upstream one) with the traumatic traces were analyzed. Utilizing discs and cores collected from trees, which were not affected by seismic shocks, the local 1154-years (AD 856–2009) TRC on Pinus sibirica Du Tour was built (Fig. 6), and verified by the “Jelo” tree-ring chronology (Agatova et al., 2014b). This local TRC has provided the dendrochronological dating of revealed wood penetrating injuries.
The year of birth of the eldest examined trees which colonized the surface of the talus fans covering landslide bodies is AD 1069 and AD 856 for the upstream and downstream landslides respectively. In order to estimate the time of creation of earthquake triggered landslides in addition to the germination ages of trees growing on their surfaces, the colonization time gap (no less than 50–100 years as it was evidenced by germination ages of trees that settled the LIA moraines in the SE Altai (Okishev, 1982)) and the formation time of talus fans covered landslide bodies, as well as surface stabilization period, should be taken into consideration. Due to these reasons the applied age correction can reach two centuries. Thus, it could be asserted that by AD 600–700 the studied seismically triggered landslides in the Arydjan valley had been existed. The currently available data do not allow clearly distinguish the time of landslides formation.
From 120 traumatic injuries of tree rings, three and more coincide into the years 1316, 1422, 1532 and 2003 (Agatova et al., 2014b). At the same time numerous wood penetrating injuries dated by AD 1532 are indicated for trees grown at both landslides and are displayed at various heights of tree trunks. This fact can imply a high magnitude of a medieval seismic event or/and its nearby epicenter location. Same patterns of wood penetrating injuries caused by rock-falls, triggered by the 2003 Chuya earthquake with its epicenter located on the southern border of the Chagan-Uzun massif, allow us to assume that a strong earthquake struck SE Altai in 1532. Most likely, as in 2003, in 1532 the fault boundaries of the Chagan-Uzun massif were reactivated. The expected date of this medieval seismic event was later verified by radiocarbon dating.
In addition to dendroseismological investigations the radiocarbon dating of earthquake induced landforms and seismically affected sediments was applied. Twenty-five radiocarbon ages were obtained for previously unknown effects of prehistoric earthquakes located along fault boundary between Chuya depression and North Chuya range, southward of the Chagan-Uzun massif (Table 2).
Table 2
Calibrated age (cal BP) | ||||||
---|---|---|---|---|---|---|
Sample Lab. cod | Location Fig. 4 | Sample type | 14C age* (BP) | Interpretation | ||
Significance level 68.2% | Significance level 95.4% | |||||
IGAN 4090 | T1 | seismically cut fossil soil exhumed in rupture of the 2003 Chuya earthquake | ha 2240±60 | 2340–2290 (17.3%) 2260–2150 (50.9%) | 2360–2110 (95.4%) | low possible date of strong earthquake that cut and buried soil layer |
IGAN 4092 | T1 | char coal from seismically deformed fossil soil | cc 1460±70 | 1420–1290 (68.2%) | 1530–1280 (95.4%) | low possible date of strong earthquake that deformed and buried soil layer |
IGAN 4103 | T2 | paleosoils overlaid by colluvium sediments | ha 510±60 | 630–600 (13.3%) 560–500 (54.9%) | 660–460 (95.4%) | period of tectonic respite and soil formation at the foot of the steep slope |
IGAN 4104 | T2 | paleosoils overlaid by colluvium sediments | ha 600±60 | 650–580 (50.8%) 570–540 (17.4%) | 670–520 (95.4%) | period of tectonic respite and soil formation at the foot of the steep slope |
IGAN 4105 | T1 | seismically cut fossil soil exhumed in rupture of the 2003 Chuya earthquake | ha 950±60 | 930–790 (68.2%) | 970–730 (95.4%) | low possible date of strong earthquake that cut and buried soil layer |
IGAN 4106 | T2 | paleosoils overlaid by colluvium sediments | ha 710±70 | 730–630 (50.6%) 600–560 (17.6%) | 770–540 (95.4%) | period of tectonic respite and soil formation at the foot of the steep slope |
SOAN 8416 | T1 | seismically cut peat layer covered lacustrine sediments with seismic convolutions | ha 2040±55 | 2100–2080 (3.9%) 2070–1920 (64.3%) | 2150–1880 (95.4%) | low possible date of earthquake that formed convolution structures |
SOAN 8417 | T1 | charcoal from paleosoil layer | cc 2755±40 | 2880–2780 (68.2%) | 2950–2770 (95.4%) | low possible date of strong earthquake that deformed and buried soil layer |
SOAN 8418 | T1 | fossil soil with charcoal overlaid by colluvium sediments | ha 2535±40 | 2750–2690 (26.8%) 2640–2610 (9.0%) 2590–2500 (32.5%) | 2750–2480 (95.4%) | period of tectonic respite and soil formation |
SOAN 8419 | T1 | fossil soil with charcoal overlaid by colluvium sediments | ha 3335±40 | 3640–3550 (48.6%) 3540–3490 (19.6%) | 3690–3660 (3.9%) 3650–3460 (91.5%) | period of tectonic respite and soil formation |
SOAN 8425 | K | wood fragment from humified loam layer overlaid by colluvium sediments | w 845±35 | 790–700 (68.2%) | 900–860 (6.4%) 830–810 (1.1%) 800–680 (87.9%) | stabilization period of the alluvial-colluvial fan at the foot of the rocky slope |
SOAN 8549 | K | humified loam layer overlaid by colluvium sediments | ha 3275±80 | 3590–3400 (68.2%) | 3700–3340 (95.4%) | stabilization period of the alluvial–colluvial fan at the foot of the rocky slope |
SOAN 8658 | T1 | wood fragments from fossil peat layer in the proximal part of the giant landslide triggered by the 2003 Chuya earthquake | w 650±65 | 670–620 (31.9%) 610–550 (36.3%) | 690–530 (95.4%) | upper possible date of strong earthquake that coursed peat bog formation |
SOAN 8659 | T1 | contemporary soil overlapping seismically cut paleosoils | ha 250±30 | 320–280 (50.2%) 170–150 (18.0%) | 430–360 (15.0%) 330–270 (55.2%) 190–140 (21.3%) 20–0 (4.0%) | upper possible date of strong paleoearthquake |
SOAN 8663 | T1 | seismically cut fossil soil with charcoal | ha 2380±65 | 2680–2640 (9.5%) 2610–2600 (1.6%) 2500–2340 (57.0%) | 2720–2310 (95.4%) | low possible date of strong earthquake that cut and buried soil layer |
SOAN 8664 | T1 | fossil peat layer with charcoal in the proximal part of the giant landslide triggered by the 2003 Chuya earthquake | ha 1035±45 | 1050–1030 (4.6%) 990–910 (63.6%) | 1060–900 (87.2%) 870–800 (8.2%) | upper possible date of strong earthquake that coursed peat bog formation |
SOAN 8665 | T1 | seismically cut fossil soil with charcoal | ha 2295±40 | 2360–2300 (55.0%) 2230–2200 (13.2%) | 2360–2290 (59.4%) 2270–2150 (36.0%) | low possible date of strong earthquake that cut and buried soil layer |
SOAN 8666 | T1 | wood fragments from fossil soil layer | w 1935±35 | 1930–1860 (54.5%) 1850–1820 (13.7%) | 1990–1810 (95.4%) | low possible date of strong earthquake that deformed and buried soil layer |
SOAN 8667 | T1 | wood fragments from fossil soil layer | w 2390±65 | 2680–2640 (10.5%) 2610–2590 (2.3%) 2500–2340 (55.4%) | 2720–2320 (95.4%) | low possible date of strong earthquake that deformed and buried soil layer |
SOAN 8668 | T1 | wood fragments from fossil soil layer | w 2490±35 | 2720–2680 (13.1%) 2640–2490 (55.1%) | 2740–2430 (95.4%) | low possible date of strong earthquake that deformed and buried soil layer |
SOAN 8669 | T1 | fossil soil with charcoal | ha 2540±50 | 2750–2690 (25.5%) 2640–2610 (8.4%) 2600–2500 (34.3%) | 2760–2460 (95.4%) | low possible date of strong earthquake that interrupt soil formation |
SOAN 8670 | T1 | charcoal in sandy loams | cc 3785±45 | 4240–4090 (68.2%) | 4380–4370 (0.5%) 4360–4320 (1.6%) 4300–4060 (85.4%) 4050–3980 (7.9%) | intensification of slope processes |
SOAN 8671 | T1 | charcoal in lacustrine sediments | cc 3765±45 | 4230–4080 (62.5%) 4030–4010 (5.7%) | 4290–4270 (1.1%) 4260–3980 (94.3%) | upper possible date of strong earthquake that coursed lake formation |
SOAN 8672 | T1 | buried peat | ha 3590±45 | 3970–3940 (11.2%) 3930–3830 (57.0%) | 4080–4040 (3.3%) 4000–3810 (84.4%) 3800–3720 (7.7%) | upper possible date of strong earthquake that coursed peat bog formation |
SOAN 8673 | T1 | buried peat layer in the proximal part of the giant landslide triggered by the 2003 Chuya earthquake | ha 770±40 | 730–670 (68.2%) | 770–660 (95.4%) | upper possible date of strong earthquake that coursed peat bog formation |
Most samples were collected in newly arisen outcrops produced within the largest landslide triggered by the 2003 Chuya earthquake in the Taldura valley. These are evidences of past seismic events: cut paleosoils and peat layers, redeposited fragments of charcoals resulted in inversion of radiocarbon dates in outcrops, seismic convolutions in sediments of ancient lakes, which was developed in the same manner as modern lake in the back area of recently triggered landslide.
The ages of seismically cut fossil soils with wood/charcoal fragments return the low possible date of seismic event while overlaid undistorted layers give information about its upper possible date. Thus obtained radiocarbon chronology from the newly arisen outcrop in the proximal part of the earthquake triggered landslide in the Taldura valley (Fig. 7; Fig. 3, location 1) supports the expected date of the 1532 seismic event. The ages of seismically deformed fossil soil and buried peat layer are 2360–2110 cal BP (IGAN 4090) and 970–730 cal BP (IGAN 4105) and the age of undistorted fossil soil covering them is 430–140 cal BP (SOAN 8659). The dendrochronologically obtained earthquake age lies inside the time limits determined by radiocarbon dating.
Radiocarbon dates of buried peat layers as well as fragment of tree trunk in lacustrine deposits returns the upper possible age of earthquakes that coursed peat bog or lake formation in the proximal part of previously existed at the same place paleolandslide (690–530, 770–660, 1060–800, 4080–3720, 4290–3980 cal BP (SOAN 8658, 8673, 8664, 8672, 8671)). In section (Fig. 8; Fig. 3, location 1) the seismically cut buried peat layer covers lacustrine sediments with seismic convolutions and underlie deposits of younger lake. Marking the period of tectonic quiescence and peat bog formation the age of peat (2150–1880 cal BP (SOAN 8416)) also indicate time period between two possible prehistoric seismic events. One of them stipulated developing of seismic convolutions and another caused cutting of peat layer and lake formation.
Deposits exhumed in an outcrop in one of the scarps of the giant landslide triggered by the 2003 Chuya earthquake in the Taldura valley indicate a complicated slope sedimentation patterns. A section of redeposited slope material with broken and deformed thin soil layers stuffed with charcoals and charred tree fragments overlays brownish sandy loams above moraines and layer of buried soil of 2760–2460 cal BP (SOAN 8669) (Fig. 3, location 2). Radiocarbon dates of organic material from deformed layers of fossil soils evidence for at least two periods (about 1400 and 2600 cal BP) of abrupt intensification of slope activity which can be related to seismic events. Trees with the age about 1400 cal BP are buried by older deposits with the charcoals and tree fragments (SOAN 8417, 8666, 8667). Intensification of slope activity about 2600 cal BP is evidenced by the age of tree fragment in colluvium sediments (SOAN 8668). Seismically deformed paleosoils from another outcrop in the same part of the landslide (2750–2480, 3690–3460 cal BP (SOAN 8418, 8419)) support the earthquake occurred about 2600 cal BP and allows suggesting presence of another more old seismic event about 3500 cal BP.
Fossil soils buried by colluvium at the foot of steep slopes mark periods of tectonic quiescence and soil formation which were altered by intensive slope mass movements. Ages of these fossil soils could be also interpreted as the upper possible ages of earthquakes that triggered intensive displacement of slope strata. Several such estimates were revealed from radiocarbon dates of fossil soils in outcrops in the upper part of the Taldura valley (660–460, 670–520, 770–540 cal BP (IGAN 4103, 4104, 4106) (Fig. 4 – T2)).
Radiocarbon dates of fossil soils buried by talus fans at the foot of the Kurai range in Kurai intermountain depression (900–680, 3700–3340 cal BP (SOAN 8425, 8549) (Fig. 4 – K)) also return the upper ages of possible past seismic events that caused soils inhumation.
. Discussion
Successful implementation of tree-ring analysis in paleoseismology implies basing seismic origin of studied tree-ring anomalies. Despite of seismic reasons, changes in growth rate can be also caused by many environmental factors and nonseismic geomorphic disturbances. To conduct effective dendrochronological study Jacoby (McCalpin, 2009) suggests the following requirements: 1) damaged trees should be confined to the rupture zone; 2) diverse evidences of paleoearthquake must exist; and 3) muster chronology should be developed from trees unaffected by the paleoearthquake.
Seismic genesis of studied landforms is proved by applying a number of criteria (Nepop and Agatova, 2008). Location of forest stand in an area affected by earthquake induced slope processes is also evidenced by analyzing ground effects of the 2003 Chuya earthquake. Local muster TRC was developed from undamaged trees and from parts of tree trunks unaffected by rock-falls, thus it gives an opportunity to calculate the precise age of detected growth anomalies. Further analysis of spatial distribution of injured trees, a number of traumatic disturbances of tree-rings and their location at different height within the tree trunk, as well as comparison these results with the distribution patterns of traumatic injuries caused by the 2003 Chuya earthquake, allowed to suggest the date of seismic event, which was a trigger of rock-falls. It should be mentioned that at first, this dendrochronologically obtained date of abrupt intensification of rock-falls was just supposed date of strong earthquake. Later it was verified by paleoseismogeological data and radiocarbon dating.
Twenty-five radiocarbon dates and results of tree-ring analysis most completely cover the period from about 4000 cal BP till present time. We have obtained just few radiocarbon dates of ground effects related to strong earthquakes occurred in the first half of the Holocene. This can be explained by worse preservation of surface effects and other features induced by ancient shocks. So in this paper we present dating results only for the last 4000 years.
It should be also noted that accuracy of the radiocarbon method doesn’t allow estimating the precise earthquake ages, but it is acceptable for reliable distinguishing the time of formation of ground effects associated with different earthquakes. This is important for marking time intervals for single seismic events and then for calculating the recurrence interval of strong earthquakes.
As it was shown above, dendrochronologically obtained date (AD 1532) of strong medieval earthquake was verified by radiocarbon dating. Moreover, radiocarbon age of paleosoils buried in colluvium wedge within the modern landslide in the Taldura valley is interpreted as a date of strong prehistoric earthquake (370 ± 30 BP IGAN 3007 (Rogozhin et al., 2007)) and also supports this thesis. After calibration it falls within the time interval AD 1440–1640 (1440–1530 (55.0%) or 1550–1640 (40.4%)). Another evidence of this medieval seismic event was detected within the southern flank of the Central Kurai neotectonic foreberg. 1.5 m scarp of an active fault was dissected by trenches (Rogozhin and Platonova, 2002). A cross section of ancient fault area represents a set of three-four steep and more flat reverse faults. Radiocarbon ages (324 ± 30 BP (IGAN 1692) and 345 ± 30 BP (IGAN 1702)) of wooden splinters and charcoals from peat layers, which were cut by these reverse faults, indicate the formation of the ancient fault area after AD 1460.
Thus geological evidences of the strong medieval earthquake are presented i) at the cross points of the fault borders of the Chagan-Uzun massif with the fault borders of the Chuya intermountain depression and South Chuya range; ii) within the inner structure of the Chagan-Uzun massif along its northern border (in the Arydjan valley); and iii) at the southern flank of the Central Kurai neotectonic foreberg. The traced length of rapture zone about 70–75 km allows to assume the magnitude of this seismic event as about 7.5 (Rogozhin and Platonova, 2002).
The obtained ages of the oldest trees and estimated minimal ages of the largest seismically induced landslides in the Arydjan valley don’t allow clearly correlate them with the known earthquakes identified by region paleoseismogeological data. Most likely they are related to the paleoearthquakes that caused ruptures and fossil soils deformations in the Kuskunur valley and in the Uzuk stow (Kuskunur – Taldura watershed). Paleosoils in colluvium wedges and pressure ridges (compressional moletracks) buried by coarse-grained material have radiocarbon ages 1280 ± 30, 1160 ± 30, 1100 ± 30 BP (IGAN 2818, 2980, 2983) (Rogozhin et al., 2007). Fossil soils in colluvium wedge along tectonic scrap at the foothill of the Central Kurai neotectonic foreberg has similar radiocarbon age 1040 ± 80 BP (GIN 9082) and indicate the date of strong prehistoric earthquake (Rogozhin and Platonova, 2002). Calibrated ages indicate that strong past earthquakes, which caused these surface disturbances, struck the SE Altai between AD 660 and AD 1170.
Our own paleoseismogeological investigations also support this thesis. The radiocarbon ages of fossil soil with char coals and buried peat layer in the proximal part of the largest landslide triggered by the 2003 Chuya earthquake in the Taldura valley (1035 ± 45, 770 ± 40 BP (SOAN 8664, 8673)) return the minimal possible age of prehistoric earthquake which coursed the lake and peat bog formation. The radiocarbon age of char coal from seismically deformed fossil soil within this landslide is 1460 ± 70 BP (IGAN 4092). Thus we can assert that strong earthquakes struck the SE Altai between AD 420 and AD 1290. Probably one of these earthquakes triggered studied landslides in the Arydjan valley.
Our analysis of unknown ground effects of past seismic events indicates that strong earthquakes occurred about 600–700, 1300–1500, 2400–2700, 3400–3700 and 3800–4200 cal BP. This conclusion supports previously published data on high regional seismicity (Rogozhin and Platonova, 2002 and Rogozhin et al., 2007, 2008) and argues for the repeated reactivation of the same earthquake source zones during the Holocene. Obtained radiocarbon dates together with previously published ones clarify the chronology of seismic events within the SE Altai (Fig. 9). It should be mentioned that due to specific problems of applying radiocarbon technique for dating young events (Wagner, 1998) several too young radiocarbon dates (about 200 years) presented in some publications (Rogozhin and Platonova, 2002 and Rogozhin et al., 2007, 2008) couldn’t be utilized for reliable time distinguishing of ground effects for the seismic events of the recent past. Thus reconstructed dates of prehistoric seismic events together with the dendrochronologically obtained age of the medieval earthquake allow estimating the recurrence interval of strong earthquakes in the SE Altai as about 400 years for the last 4000 years.
Paleoseismogeological investigations in adjusted seismically active areas of the Mongolian Altai reveal traces of several more strong (M≥8) earthquakes. In the epicenter zone of the 1931 Fuyun earthquake (M=8.0) at the western part of the Mongolian Altai, the historical and paleoseismological data evidence for presence of five such earthquakes over the last 10000 years (Shumo et al., 1996 and Rogozhin et al., 2012). Five strongest seismic events with M≥8 over the last 7500 years, including the 1761 Great Mongolian earthquake, were also reconstructed in the eastern part of the Mongolian Altai (Rogozhin et al., 2013). The specified recurrence interval of strong earthquakes for the SE Altai correlates with these data published for neighboring areas of Mongolian Altai. It is shorter than intervals for adjusted seismically active areas of the Mongolian Altai due to smaller magnitudes of strongest shocks.
. Conclusion
Paleoseismogeological researches within the pleistoseist zone of the 2003 Chuya earthquake, SE Altai, reviled unknown evidences of strong past seismic events. Seismically deformed fossil soils overlaid by undistorted ones, soils and peat layers buried by colluvium at the foot of the slopes, organic material from colluvium wedges and overlaid by pressure ridges, paleopeats affected by seismic convolutions were studied in newly arisen and renewed outcrops at the southern fault boundary of the Chagan-Uzun massif in the Taldura valley. Twenty-five radiocarbon dates indicate six earthquakes that stroke SE Altai during the last 4000 years and mark time limits of these seismic events: 300–500, 600–700, 1300–1500, 2400–2700, 3400–3700, 3800–4200 cal BP.
Perspectives of applying tree-ring analysis in paleoseismogeological researches within the SE Altai were demonstrated by studying seismically triggered landslides in the northern part of the Chagan-Uzun massif in the Arydjan valley. In addition to estimating the germination ages of trees growing on the bare surfaces of these landslides, analysis of wood penetrating injuries in the individual tree ring series was suggested for dating earthquake induced paleorock-falls. Spatial distribution of injured trees, a number of traumatic disturbances of tree-rings and their location at different height within the tree trunk, as well as comparison of these results with the distribution patterns of traumatic injuries caused by the 2003 Chuya earthquake, was used as criteria of seismic origin of past rock-falls. Dendrochronologically obtained date of abrupt intensifications of rock-falls was considered as supposed date of strong medieval earthquake, which was verified by radiocarbon dating and alternative proxy data. It allows estimating the precise date (AD 1532) for seismic event which was limited 300–500 cal BP by radiocarbon analysis.
Studied effects of strong prehistoric earthquakes preserved in the northern part of the Chagan-Uzun massif, and along fault boundary between Chuya depression and North Chuya range supports thesis about repeated reactivations of the same earthquake source zones within the SE Altai.
Obtained dendrochronological and radiocarbon dates together with the previously published data confirm the high occurrence of strong (M≥7) earthquakes in the region. Generally, eight such seismic events were revealed in the region during last 4000 years. Specified recurrence interval of strong earthquakes for this time interval is about 400 years. This data argues for the high Holocene seismicity of the SE Altai and defines the Chagan-Uzun massif, which separates Chuya and Kurai intermountain depressions, as one of the major regional seismogenic structure.