Seismotectonic deformation of the volcanic edifice prior to the 1998 lava eruption of Volcán de Colima, México

The andesitic stratovolcano Volcán de Colima is one of the most active volcanoes in Mexico. The recent eruption of Volcán de Colima began in November 1998 and was preceded by a 12-month period of seismic activity that included five earthquake swarms. About 600 events with magnitudes from -0.5 to 2.7 were located within a 50-km² area extending northward from the crater of Volcán de Colima to the Pleistocene volcano Nevado de Colima. The majority of hypocenters within this area did not exceed 5 km depth below sea level. We investigated earthquake focal mechanisms and seismotectonic deformations of the volcanic edifice. Focal mechanisms during four earthquake swarms indicated normal faulting associated with extensional processes, which is in agreement with the general stress regime near the volcano revealed by field measurements of fault slips. Earthquakes in the fifth swarm had focal mechanisms associated with inverse faulting, showing a significant change in the stress situation just before the beginning of the eruption. The calculated deformations varied from 1.3᎒^-11 to 2.7᎒^-9. The first swarm of November-December 1997 resulted in a N-S horizontal elongation that was two times greater than the E-W horizontal shortening. The volume was also subject to vertical shortening. The second and third swarms, observed in March and May 1998, showed uniform horizontal N-S and E-W elongations accompanied by a vertical elongation of the volume. In June-July 1998, the situation of November-December 1997 was repeated, with N-S horizontal elongation greater than the E-W horizontal shortening of the volume accompanied by intensive vertical shortening. During the last swarm of October-November 1998, slight E-W elongation of the seismic volume was accompanied by strong N-S shortening and very slight vertical shortening. We assume that the seismic activity prior to the 1998 eruption of Volcán de Colima developed along two intersecting tectonic structures, the N-S-trending Colima rift, and the E-W-trending system of faults associated with Tamazula fault. During the first stage (November 1997-July 1998) the passageway for magma was developed along the Tamazula fault system under horizontal extension without any surface manifestation. In October-November 1998, the seismic events began to cluster along the Colima rift structures under predominantly compressional stresses; this condition culminated with the extrusion of andesitic block lava from the summit crater.     

Seismic activity related to the 1991 eruption of Colima Volcano,Mexico

Ten years after the last effusive eruption and at least 15 years of seismic quiescence, volcanic seismic activity started at Colima volcano on 14 February 1991, with a seismic crisis which reached counts of more than 100 per day and showed a diversity of earthquake types. Four other distinct seismic crises followed, before a mild effusive eruption in April 1991. The second crisis preceded the extrusion of an andesitic scoriaceous lava lobe, first reported on 1 March; during this crisis an interesting temporary concentration of seismic foci below the crater was observed shortly before the extrusion was detected. The third crisis was constituted by shallow seismicity, featuring possible mild degassing explosion-induced activity in the form of hiccups (episodes of simple wavelets that repeat with diminishing amplitude), and accompanied by increased fumarolic activity. The growth of the new lava dome was accompanied by changing seismicity. On 16 April during the fifth crisis which consisted of some relatively large, shallow, volcanic earthquakes and numerous avalanches of older dome material, part of the newly extruded dome, which had grown towards the edge of the old dome, collapsed, producing the largest avalanches and ash flows. Afterwards, block lava began to flow slowly along the SW flank of the volcano, generating frequent small incandescent avalanches. The seismicity associated with the stages of this eruptive activity shows some interesting features: most earthquake foci were located north of the summit, some of them relatively deep (7–11 km below the summit level), underneath the saddle between the Colima and the older Nevado volcanoes. An apparently seismic quiet region appears between 4 and 7 km below the summit level. In June, harmonic tremors were detected for the first time, but no changes in the eruptive activity could be correlated with them. After June, the seismicity decreasing trend was established, and the effusive activity stopped on September 1991.     

The feeding of the eruptive activity of Etna volcano. The regional stress field as a constraint to magma uprising and eruption

An analysis of the structural lineaments, as observed on the middle and upper slopes of Etna volcano, was made with the aim of checking stress distribution within the volcanic structure. The observed features suggest that the deformation pattern of the volcanic edifice is compatible with a deviating stress field dominated by an E-W sinistral shear, with the maximum stress axes being oriented NE (σ¹) and NW (σ³) respectively. Such a framework appears to be consistent with the active regional stress field, as deduced from structural data and focal mechanism analyses. The role of the active stress field in the penetration, uprising and eruption of magma in the Etnean area is therefore discussed and some more general suggestions about the geodynamic evolution of eastern Sicily are made.     

Seismic signature of a phreatic explosion: hydrofracturing damage at Karthala volcano, Grande Comore Island, Indian Ocean

Karthala volcano is a basaltic shield volcano with an active hydrothermal system that forms the southern two-thirds of the Grande Comore Island, off the east coat of Africa, northwest of Madagascar. Since the start of volcano monitoring by the local volcano observatory in 1988, the July 11th, 1991 phreatic eruption was the first volcanic event seismically recorded on this volcano, and a rare example of a monitored basaltic shield. From 1991 to 1995 the VT locations, 0.5<M_l<4.3, show a crack shaped pattern (3 km long, 1 km wide) within the summit caldera extending at depth from –2 km to +2 km relative to sea level. This N-S elongated pattern coincides with the direction of the regional maximum horizontal stress as deduced from regional focal mechanism solutions. This brittle signature of the damage associated with the 1991 phreatic eruption is a typical pattern of the seismicity induced by controlled fluid injections such as those applied at geothermal fields, in oil and gas recovery, or for stress measurements. It suggests the 1991 phreatic eruption was driven by hydraulic fracturing induced by forced fluid flow. We propose that the extremely high LP and VT seismicity rates, relative to other effusive volcanoes, during the climax of the 1991 phreatic explosion, are due to the activation of the whole hydrothermal system, as roughly sized by the distribution of VT hypocenters. The seismicity rate in 1995 was still higher than the pre-eruption seismicity rate, and disagrees with the time pattern of thermo-elastic stress readjustment induced by single magma intrusions at basaltic volcanoes. We propose that it corresponds to the still ongoing relaxation of pressure heterogeneity within the hydrothermal system as suggested by the few LP events that still occurred in 1995.Editorial responsibility: H Shinohara     

Seismogenic stress field beneath Mt. Etna (South Italy) and possible relationships with volcano-tectonic features

Shallow shear-type seismic activity occurring beneath the Etna volcano during 1990–1995 has been analysed for hypocenter locations, focal mechanisms and stress tensor inversion. The results have been examined jointly with Electronic Distance Measurements and tiltmeter data collected in the same period and reported in the literature. Significant seismicity located in the upper 10 km was found to be confined to the time intervals in which ground deformation data indicated inflation of the volcano edifice (e.g., the periods preceding the December 1991–March 1993 and August 1995–March 1996 eruptive phases). The shocks mostly occurred in a sector approximately centered on the crater area and elongated in the East–West direction. The causative seismogenic stress shows a low-dip East–West orientation of σ₁. In agreement with existing knowledge on relationships between local fault systems and magma uprise processes, the shallow seismicity in question is tentatively explained as being due to lateral compression by magma inside a nearly North–South system. The volcano deflation phase revealed by Electronic Distance Measurements and tilt data during the 1991–1993 major eruption was not accompanied by any significant shear-type shallow event. Below the depth of 10 km, the North–South prevailing orientation of σ₁ reflects the dominant role of the regional stress.     

Seismicity and crustal deformation preceding the January 1996 eruptions at Karymsky Volcanic Center,Kamchatka

Two explosive eruptions occurred on 2 January 1996 at Karymsky Volcanic Center (KVC) in Kamchatka, Russia: the first, dacitic, from the central vent of Karymsky volcano, and the second, several hours later, from Karymskoye lake in the caldera of Akademia Nauk volcano. The main significance of the 1996 volcanic events in KVC was the phreatomagmatic eruption in Karymskoye lake, which was the first eruption in this lake in historical time, and was a basaltic eruption at the acidic volcanic center. The volcanic events were associated with the 1 January Ms 6.7 (Mw 7.1) earthquake that occurred at a distance of about 9–17 km southeast from the volcanoes just before the eruptions. We study the long-term (1972–1995) and short-term (1–2 January 1996) characteristics of crustal deformations and seismicity before the double eruptive event in KVC. The 1972–1995 crustal deformation was homogeneous and characterized by a gradual extension with a steady velocity. The seismic activity in 1972–1995 developed at the depth interval from 0 to 20 km below the Akademia Nauk volcano and spread to the southeast along a regional fault. The seismic activity in January 1996 began with a short sequence of very shallow microearthquakes (M ~0) beneath Karymsky volcano. Then seismic events sharply increased in magnitude (up to mb 4.9) and moved along the regional fault to the southeast, culminating in the Ms 6.7 earthquake. Its aftershocks were located to the southeast and northwest from the main shock, filling the space between the two active volcanoes and the ancient basaltic volcano of Zhupanovsky Vostryaki. The eruption in Karymskoye lake began during the aftershock sequence. We consider that the Ms 6.7 earthquake opened the passageway for basic magma located below Zhupanovsky Vostryaki volcano that fed the eruption in Karymskoye lake.     

The July 1994 episode of seismic activity at Colima Volcano, Mexico

A new period of seismic activity that culminated in a small phreatic explosion took place in Colima Volcano (Western Mexico) during the month of July 1994. In this note, we present our analysis of this seismicity based upon information from RESCO, the seismic network of the University of Colima. The activity began with a seismic swarm of type A (tectonic-like) earthquakes with epicenters towards the SSW of the summit, followed by shallow low-frequency events underneath the volcanic edifice. The activity was accompanied by landslides and culminated with an explosion that produced small ash falls on the surrounding area. The seismic activity ceased after this episode.     

Chemistry and mineralogy of high-temperature gas discharges from Colima volcano, Mexico. Implications for magmatic gas–atmosphere interaction

Gases, condensates and silica tube precipitates were collected from 400°C (Z2) and 800°C (Z3) fumaroles at Colima volcano, Mexico, in 1996–1998. Volcanic gases at Colima were very oxidized and contain up to 98% air due to mixing with air inside the dome interior, close to the hot magmatic body. An alkaline trap method was used to collect gas samples, therefore only acidic species were analysed. Colima volcanic gases are water-rich (95–98 mol%) and have typical S/C/Cl/F ratios for a subduction type volcano. δD-values for the high-temperature Z3 fumarolic vapour vary from −26 to −57‰. A negative δD–Cl correlation for the Z3 high-temperature fumarole may result from magma degassing: enrichment in D and decrease in the Cl concentration in condensates are likely a consequence of input of “fresh” batches of magma and an increasing of volcanic activity, respectively.The trace element composition of Colima condensates generally does not differ from that of other volcanoes (e.g. Merapi, Kudryavy) except for some enrichment in V, Cu and Zn. Variations in chemical composition of precipitates along the silica tube from the high-temperature fumarole (Colima 1, fumarole Z3), in contrast to other volcanoes, are characterized by high concentrations of Ca and V, low concentration of Mo and a lack of Cd. Mineralogy of precipitates differs significantly from that described for silica tube experiments at other volcanoes with reduced volcanic gas. Thermochemical modelling was used to explain why very oxidized gas at Colima does not precipitate halite, sylvite, and Mo- and Cd-minerals, but does precipitate V-minerals and native gold, which have not been observed before in mineral precipitates from reduced volcanic gases.     

Random patterns of occurrence of explosive eruptions at Colima Volcano, Mexico

Analysis of the patterns of eruption occurrences may improve our understanding of volcanic processes. In this paper, the available historical data of an individual volcano, Colima in México, are used to classify its eruptions by size using the Volcanic Explosivity Index (VEI). The data shows that, if eruptions are only taken into account above a certain VEI level, the stochastic process associated with the explosive volcanic events can be represented by a non-stationary Poisson point process, which can be reduced to a homogeneous Poisson process through a transformation of the time axis. When eruptions are separated by VEI values, the occurrence patterns of each magnitude category can also be represented by a Poisson distribution. Analysis of the rate of occurrence of all eruptions with VEI greater than 1 permits the recognition of three distinct regimes or rates of volcanic activity during the last 430 years. A double stochastic Poisson model is suggested to describe this non-stationary eruptive pattern of Colima volcano and a Bayesian approach permits an estimation the present hazard.     

Earthquake activity related to the 1991 eruption of the Hekla volcano,Iceland

The 1991 eruption of the Hekla volcano started unexpectedly on 17 January. No long-term precursory seismicity was observed. The first related activity was a swarm of small earthquakes that began approximately half an hour before the eruption. Intensive seismicity, both earthquakes and volcanic tremor, accompanied the violent onset of the eruption. Almost 400 events up to M_L magnitude 2.5 were recorded during the first few hours. During the later phases of the eruption, the earthquake activity was modest and the main volcano-related seismic signal was the persistent volcanic tremor. The tremor died away, together with the eruption on 11 March, and Hekla was seismically quiet until the beginning of June 1991, when a sudden swarm of numerous small shallow earthquakes occurred. This activity is atypical for Hekla and is interpreted to be a failed attempt to resume the eruption.     

Volcanic Tremor at Mt. Etna,Italy, Preceding and Accompanying the Eruption of July – August, 2001

The July 17 – August 9, 2001 flank eruption of Mt. Etna was preceded and accompanied by remarkable changes in volcanic tremor. Based on the records of stations belonging to the permanent seismic network deployed on the volcano, we analyze amplitude and frequency content of the seismic signal. We find considerable changes in the volcanic tremor which mark the transition to different styles of eruptive activity, e.g., lava fountains, phreatomagmatic activity, Strombolian explosions. In particular, the frequency content of the signal decreases from 5 Hz to 3 Hz at our reference station ETF during episodes of lava fountains, and further decreases at about 2 Hz throughout phases of intense lava emission. The frequency content and the ratios of the signal amplitude allow us to distinguish three seismic sources, i.e., the peripheral dike which fed the eruption, the reservoir which fed the lava fountains, and the central conduit. Based on the analysis of the amplitude decay of the signal, we highlight the migration of the dike from a depth of ca. 5 km to about 1 km between July 10 and 12. After the onset of the effusive phase, the distribution of the amplitude decay at our stations can be interpreted as the overall result of sources located within the first half kilometer from the surface. Although on a qualitative basis, our findings shed some light on the complex feeding system of Mt. Etna, and integrate other volcanological and geophysical studies which tackle the problem of magma replenishment for the July–August, 2001 flank eruption. We conclude that volcanic tremor is fundamental in monitoring Mt. Etna, not only as a marker of the different sources which act within the volcano edifice, but also of the diverse styles of eruptive activity. An erratum to this article is available at .     

Dynamical parameter analysis of continuous seismic signals of Popocatépetl volcano (Central Mexico): A case of tectonic earthquakes influencing volcanic activity

The continuous background seismic activity contains information on the internal state of a volcanic system. Here, we report the influence of major regional tectonic earthquakes (M > 5 in most cases) on such state, reflected as changes in the spectral and dynamical parameters of the volcano continuous seismic data. Although changes do not always occur, analysis of five cases of earthquake-induced variations in the signals recorded at Popocatépetl volcano in central México reveal significant fluctuations following the tectonic earthquakes. External visible volcanic activity, such as small to moderate explosions and ash emissions, were related to those fluctuations. We briefly discuss possible causes of the variations. We conclude that recognition of fluctuations in the dynamical parameters in volcano monitoring seismic signals after tectonic earthquakes, even those located in the far field, hundreds of kilometers away, may provide an additional criterion for eruption forecasting, and for decision making in the definition of volcanic alert levels.     

A tectonic earthquake sequence preceding the April–May 1999 eruption of Shishaldin Volcano, Alaska

On 4 March 1999, a shallow M_L 5.2 earthquake occurred beneath Unimak Island in the Aleutian Arc. This earthquake was located 10-15 km west of Shishaldin Volcano, a large, frequently active basaltic-andesite stratovolcano. A Strombolian eruption began at Shishaldin roughly 1 month after the mainshock, culminating in a large explosive eruption on 19 April. We address the question of whether or not the eruption caused the mainshock by computing the Coulomb stress change caused by an inflating dike on fault planes oriented parallel to the mainshock focal mechanism. We found Coulomb stress increases of ~0.1 MPa in the region of the mainshock, suggesting that magma intrusion prior to the eruption could have caused the mainshock. Satellite and seismic data indicate that magma was moving upwards beneath Shishaldin well before the mainshock, indicating that, in an overall sense, the mainshock cannot be said to have caused the eruption. However, observations of changes at the volcano following the mainshock and several large aftershocks suggest that the earthquakes may, in turn, have influenced the course of the eruption.     

The 1997–1998 Activity of Volcán de Colima,Western Mexico: Some Aspects of the Associated Seismic Activity

Volcán de Colima, the most active volcano in Mexico, had a climactic episode on 20 November, 1998. On this date, a dome formed on the small summit crater during the previous few days, collapsed generating block-and-ash flows. The event was preceded by almost twelve months of seismic activity, which continued afterwards for several more months. We analyzed the main seismic activity, which occurred from 20 March, 1998 to 31 March, 1999. The seismicity was dominated by volcano-tectonic earthquakes before the climax, and subsequently by hybrid and long-period earthquakes. We determined the frequency of events for the entire period, and located most of the volcano-tectonic events. To assess the possibility that these earthquakes were generated by the same source, they were tested for their similitude through cross correlation in the time domain. Six groups of similar events, or earthquake families, were generated. The members of these families appeared before the 20 November event, apparently ceasing afterwards. We examined the location of the families' events with respect to an existing gravity model in which an anomalous body of negative density contrast suggests the presence of the magma chamber. Most of the family events occur on top of the anomalous body, which suggests they were associated with the passage of magma through the feeding conduits of the volcano.     

Mechanism of volcanic earthquakes of the Sheveluch volcano,Kamchatka

Determinations of the local mechanisms of three volcanic earthquakes are given connected with the eruption of the Sheveluch volcano (November, 1964). As initial material the data on first arrivals of P-waves are used. The local mechanism of all three earthquakes is close to a strike-slip type of faulting and similar to the focal mechanism of tectonic earthquakes of Kamchatka. One nodal surface of all the volcanic earthquakes strikes in the same direction as the outbursts of the directed volcano explosions.     

The geological evolution of Merapi volcano, Central Java, Indonesia

Merapi is an almost persistently active basalt to basaltic andesite volcanic complex in Central Java (Indonesia) and often referred to as the type volcano for small-volume pyroclastic flows generated by gravitational lava dome failures (Merapi-type nuées ardentes). Stratigraphic field data, published and new radiocarbon ages in conjunction with a new set of ⁴⁰K–⁴⁰Ar and ⁴⁰Ar–³⁹Ar ages, and whole-rock geochemical data allow a reassessment of the geological and geochemical evolution of the volcanic complex. An adapted version of the published geological map of Merapi [(Wirakusumah et al. 1989), Peta Geologi Gunungapi Merapi, Jawa Tengah (Geologic map of Merapi volcano, Central Java), 1:50,000] is presented, in which eight main volcano stratigraphic units are distinguished, linked to three main evolutionary stages of the volcanic complex—Proto-Merapi, Old Merapi and New Merapi. Construction of the Merapi volcanic complex began after 170?ka. The two earliest (Proto-Merapi) volcanic edifices, Gunung Bibi (109?±?60?ka), a small basaltic andesite volcanic structure on Merapi’s north-east flank, and Gunung Turgo and Gunung Plawangan (138?±?3?ka; 135?±?3?ka), two basaltic hills in the southern sector of the volcano, predate the Merapi cone sensu stricto. Old Merapi started to grow at ~30?ka, building a stratovolcano of basaltic andesite lavas and intercalated pyroclastic rocks. This older Merapi edifice was destroyed by one or, possibly, several flank failures, the latest of which occurred after 4.8?±?1.5?ka and marks the end of the Old Merapi stage. The construction of the recent Merapi cone (New Merapi) began afterwards. Mostly basaltic andesite pyroclastic and epiclastic deposits of both Old and New Merapi (<11,792?±?90 ¹⁴C years BP) cover the lower flanks of the edifice. A shift from medium-K to high-K character of the eruptive products occurred at ~1,900 ¹⁴C years BP, with all younger products having high-K affinity. The radiocarbon record points towards an almost continuous activity of Merapi since this time, with periods of high eruption frequency interrupted by shorter intervals of apparently lower eruption rates, which is reflected in the geochemical composition of the eruptive products. The Holocene stratigraphic record reveals that fountain collapse pyroclastic flows are a common phenomenon at Merapi. The distribution and run-out distances of these flows have frequently exceeded those of the classic Merapi-type nuées ardentes of the recent activity. Widespread pumiceous fallout deposits testify the occurrence of moderate to large (subplinian) eruptions (VEI 3–4) during the mid to late Holocene. VEI 4 eruptions, as identified in the stratigraphic record, are an order of magnitude larger than any recorded historical eruption of Merapi, except for the 1872?AD and, possibly, the October–November 2010 events. Both types of eruptive and volcanic phenomena require careful consideration in long-term hazard assessment at Merapi.     

The November 2002 eruption at Piton de la Fournaise volcano,La Réunion Island: ground deformation,seismicity, and pit crater collapse

An eruption on the eastern flank of Piton de la Fournaise volcano started on 16 November, 2002 after 10 months of quiescence. After a relatively constant level of activity during the first 13 days of the eruption, lava discharge, volcanic tremor and seismicity increased from 29 November to 3 December. Lava effusion suddenly ceased on 3 December while shallow earthquakes beneath the Dolomieu summit crater were still recorded at a rate of about one per minute. This unusual activity continued and increased in intensity over the next three weeks, ending with the formation of a pit crater within Dolomieu. Based on ground deformation, measured by rapid-static and continuous GPS and an extensometer, seismic data, and lava effusion patterns, the eruptive period is divided into five stages: 1) slow summit inflation and sporadic seismicity; 2) rapid summit inflation and a short seismic crisis; 3) rapid flank inflation, onset of summit deflation, sporadic seismicity, accompanied by stable effusion; 4) flank inflation, coupled with summit deflation, intense seismicity, and increased lava effusion; and finally 5) little deflation, intense shallow seismicity, and the end of lava effusion. We propose a model in which the pre-intrusive inflation of Stage 1 in the months preceding the eruption was caused by a magma body located near sea level. The magma reservoir was the source of an intrusion rising under the summit during Stage 2. In Stage 3, the magma ponded at a shallow level in the edifice while the lateral injection of a radial dike reached the surface on the eastern flank of the basaltic volcano, causing lava effusion. Pressure decrease in the magmatic plumbing system followed, resulting in upward migration of a collapse front, forming a subterranean column of debris by faulting and stoping. This caused intense shallow seismicity, increase in discharge of lava and volcanic tremor at the lateral vent in Stage 4 and, eventually the formation of a pit crater in Stage 5.     

长白山天池火山减灾对策初探

国内外专家学者认为，长白山天池火山是一座具潜在灾害性喷发危险的活火山，因此制定火山减灾对策理应提到议事日程。针对天池火山研究现状和火山灾害特点，制定了火山活动各阶段的减灾对策。中长期阶段应加强火山监测与研究和火山知识宣传工作，采取必要的工程防护措施，重大工程进行火山安全性评价，制定火山喷发应急预案；短期阶段请求国际火山流动监测台网给予支援；临近喷发阶段重点是有组织的撤离；喷发及其后阶段应及时救灾抢险，对火山喷发趋势进行科学判定，合理地重建家园。     

Seismicity and magma supply rate of the 1998 failed eruption at Iwate volcano,Japan

Iwate volcano, Japan, showed significant volcanic activity including earthquake swarms and volcano inflation from the beginning of 1998. A large earthquake of magnitude 6.1 hit the south-west of the volcano on September 3. Although a 1 km² fumarole field formed, blighting plants on the ridge in the western part of the volcano in the spring of 1999, no magmatic eruptions occurred. We reconcile the spatio-temporal distributions of volcanic pressure sources determined by previously reported studies in which GPS, strain and tilt data from dense geodetic station networks are analyzed (Miura et al. Earth Planet Space 52:1003–1008, 2000; Sato and Hamaguchi J Volcanol Geotherm Res 155:244–262, 2006). We calculate the magma supply rates from their results and compare them with the occurrence rates of volcanic earthquakes. The results show that the magma supply rates are almost constant or even decrease with time while the earthquake occurrence rate increases with time. This contrast in their temporal changes is interpreted to result from stress accumulation in the volcanic edifice caused by constant magma supply without effusion of magma to the surface. We further show that data showing slight acceleration in strain can be best explained by magma ascent at a constant velocity, and that there is no evidence for increased magma buoyancy resulting from gas bubble growth. This consideration supports the interpretation that the magma stayed at 2 km depth and horizontally migrated. These findings relating magma supply rate and seismicity to magma ascent process are clues to understanding why no magmatic eruption occurred at Iwate volcano in 1998.     

INVERSION RESEARCH OF PALEOSTRESS FIELD IN TENGCHONG AREA BASED ON VENT ALIGNMENT SURVEY METHOD

The acquisition of paleostress data has always been a key and difficult problem in the field of tectonics and crustal dynamics. Previous findings show that the distribution of dykes and the alignment of volcano vents can indicate the direction of major principal stresses during the volcanic eruption period. The extension direction of underground dyke is generally parallel to the linear arrangement direction of vent and both of them give the direction of maximum horizontal principal stress in this period. This law has similar mechanical principle to hydraulic fracturing. Vent alignment survey method is a kind of geostress inversion method based on volcanic linear mode. In this paper theoretical basis of vent alignment survey method and the linear model of volcanic array are introduced. Based on the panchromatic high-resolution remote sensing image and digital elevation model(DEM)data in the Tengchong area, combining the geological exploration data about lithologic distribution of igneous rocks and the past research results of volcanic activities in the area, the research divides the volcanoes in this area, according to the age of volcanic eruption, into Holocene, late Pleistocene, early Pleistocene and Pliocene. Nakamura model was established based on the morphological characteristics and spatial distribution of the crater, and volcanic eruption period to evaluate the linear degree of the volcano in different age, and to carry out inversion of the paleostress field of the four volcanic eruption periods in Tengchong area. The reliability of the inversion results is evaluated according to Paulsen's data reliability assessment system. The data from inverting the paleostress of each show that maximum horizontal principal stress directions in the early Pleistocene period and late Pleistocene period are 39° and 37° respectively and reliability grades of inverted paleostress results in the two eras all belong to grade B; Inverted paleostress results in Holocene and Pliocene are not reliable because of the poor linearity of volcanic arrangement. The geostress direction determined by the vent alignment survey method is fairly consistent with the current geostress field obtained by the focal mechanism solution and the anelastic strain recovery(ASR method)and matches the tectonic history of Tengchong area since Neogene. The research findings demonstrate that the vent alignment survey method is an effective paleostress inversion method. The application of this method in the paleostress inversion of the Tengchong volcanic area is of important reference significance for the study of paleostress field in other volcanic areas in China.