共查询到20条相似文献,搜索用时 10 毫秒
1.
The 1783–1784 Laki tholeiitic basalt fissure eruption in Iceland was one of the greatest atmospheric pollution events of
the past 250 years, with widespread effects in the northern hemisphere. The degassing history and volatile budget of this
event are determined by measurements of pre-eruption and residual contents of sulfur, chlorine, and fluorine in the products
of all phases of the eruption. In fissure eruptions such as Laki, degassing occurs in two stages: by explosive activity or
lava fountaining at the vents, and from the lava as it flows away from the vents. Using the measured sulfur concentrations
in glass inclusions in phenocrysts and in groundmass glasses of quenched eruption products, we calculate that the total accumulative
atmospheric mass loading of sulfur dioxide was 122 Mt over a period of 8 months. This volatile release is sufficient to have
generated ∼250 Mt of H2SO4 aerosols, an amount which agrees with an independent estimate of the Laki aerosol yield based on atmospheric turbidity measurements.
Most of this volatile mass (∼60 wt.%) was released during the first 1.5 months of activity. The measured chlorine and fluorine
concentrations in the samples indicate that the atmospheric loading of hydrochloric acid and hydrofluoric acid was ∼7.0 and
15.0 Mt, respectively. Furthermore, ∼75% of the volatile mass dissolved by the Laki magma was released at the vents and carried
by eruption columns to altitudes between 6 and 13 km. The high degree of degassing at the vents is attributed to development
of a separated two-phase flow in the upper magma conduit, and implies that high-discharge basaltic eruptions such as Laki
are able to loft huge quantities of gas to altitudes where the resulting aerosols can reside for months or even 1–2 years.
The atmospheric volatile contribution due to subsequent degassing of the Laki lava flow is only 18 wt.% of the total dissolved
in the magma, and these emissions were confined to the lowest regions of the troposphere and therefore important only over
Iceland. This study indicates that determination of the amount of sulfur degassed from the Laki magma batch by measurements
of sulfur in the volcanic products (the petrologic method) yields a result which is sufficient to account for the mass of
aerosols estimated by other methods.
Received: 30 May 1995 / Accepted: 19 April 1996 相似文献
2.
Two magnetotelluric (MT) surveys were carried out on the Mt. Etna volcano after two of the most intense eruptions of the last 30 years which took place in summer 2001 and winter 2002–2003. Surveying was pursued for two main reasons. First, we sought to contribute to the definition of the first-order structure and physico-chemical state (temperature, fluids, melts) of a volcano that has been extensively explored and monitored by means of various geophysical methods, but where only few electrical and electromagnetic surveys have been performed. Secondly, we acquired MT data in the same sites in the two different surveys with the aim of monitoring the possible changes of the first-order structure, since conditions are expected to vary on an active volcano such as Etna, and are supposed to be linked to the eruptive events. Soundings have been acquired in an E-W 10 km-long profile across the southern flank of Mt. Etna, at a distance of almost 6 km south from the Central Crater. The first survey was carried out three months after the 2001 eruption. Inverse models define a pronounced (4 km thickness) low resistivity section at a depth of about 1 km b.s.l. to the west. To the east, a low resistivity section is still present, but appears deeper, thinner and more resistive, and a shallow low resistivity anomaly also exists. The shallow anomaly to the east is tentatively correlated with altered and clayey volcanic units and/or temporary groundwater storage. The deep anomalies are interpreted as being due to melt storage at shallow depths which was not exhausted during the eruption. This would be confirmed by the abundance of lava erupted within one year from the end of the survey. The few good sites retrieved in the second survey, carried out a few weeks after the eruption of 2002–2003, confirm the picture defined in the first survey, and provide a better definition of the bottom of the deep anomaly located in the sedimentary basement. 相似文献
3.
The complexity of the processes responsible for volcanic eruptions makes a theoretical approach to forecasting the evolution of volcanic unrest rather difficult. A feasible strategy for this purpose appears to be the identification of possible repetitive schemes (patterns) in the pre-eruptive unrest of volcanoes. Nevertheless, the limited availability and the heterogeneity of pre-eruptive data, and the objective difficulty in quantitatively recognizing complex pre-eruptive patterns, make this task very difficult. In this work we address this issue by using a pattern recognition approach applied to the seismicity recorded during 217 volcanic episodes of unrest around the world. In particular, we use two non-parametric algorithms that have proven to give satisfactory results in dealing with a small amount of data, even if not normally distributed and/or characterized by discrete or categorical values. The results show evidence of a longer period of instability in the unrest preceding an eruption, compared to isolated unrest. This might indicate, even if not necessarily, a difference in the energy of processes responsible for the two types of unrest. However, if the unrest is followed by an eruption, it seems that the seismic energy released during the unrest (parameterized by the duration of the swarm and the maximum magnitude recorded) is not indicative of the magnitude of the impending eruption. We also found that, in general, unrest followed by the largest explosive eruptions have a longer repose time than those related to moderate eruptions. This evidence supports the fact that the occurrence of a large eruption needs a sufficient amount of time after the last event in order to re-charge the feeding system and to achieve a closed-conduit regime so that a sufficiently large amount of gas can be accumulated.Editorial responsibility: T. Druitt 相似文献
4.
Sediment transport models require appropriate representation of near-bed processes. We aim here to explore the parameterizations of bed shear stress, bed load transport rate and near-bed sediment erosion rate under the sheet flow regime. To that end, we employ a one-dimensional two-phase sheet flow model which is able to resolve the intrawave boundary layer and sediment dynamics at a length scale on the order of the sediment grain. We have conducted 79 numerical simulations to cover a range of collinear wave and current conditions and sediment diameters in the range 210–460 μm. The numerical results confirm that the intrawave bed shear stress leads the free stream velocity, and we assess an explicit expression relating the phase lead to the maximum velocity, wave period and bed roughness. The numerical sheet flow model is also used to provide estimates for the bed load transport rate and to inspect the near-bed sediment erosion. A common bed load transport rate formulation and two typical reference concentration approaches are assessed. A dependence of the bed load transport rate on the sediment grain diameter is observed and parameterized. Finally, the intrawave near-bed vertical sediment flux is further investigated and related to the time derivative of the bed shear stress. 相似文献
5.
The role of carbon dioxide in the dynamics of magma ascent in explosive eruptions is investigated by means of numerical modeling.
The model is steady, one-dimensional, and isothermal; it calculates the separated flow of gas and a homogeneous mixture of
liquid magma and crystals. The magma properties are calculated on the basis of magma composition and crystal content and are
allowed to change along the conduit due to pressure decrease and gas exsolution. The effect of the presence of a two-component
(water + carbon dioxide) exsolving gas phase is investigated by performing a parametric study on the CO2/(H2O+CO2) ratio, which is allowed to vary from 0 to 0.5 at either constant total volatile or constant water content. The relatively
insoluble carbon dioxide component plays an important role in the location of the volatile-saturation and magma-fragmentation
levels and in the distribution of the flow variables in the volcanic conduit. In detail, the results show that an increase
of the proportion of carbon dioxide produces a decrease of the mass flow rate, pressure, and exit mixture density, and an
increase of the exit gas volume fraction and depth of the fragmentation level. A relevant result is the different role played
by water and carbon dioxide in the eruption dynamics; an increasing amount of water produces an increase of the mass flow
rate, and an increasing amount of carbon dioxide produces a decrease. Even small amounts of carbon dioxide have major consequences
on the eruption dynamics, implying that the multicomponent nature of the volcanic gas must be taken into account in the prediction
of the eruption scenario and the forecasting of volcanic hazard.
Received: 6 March 1998 / Accepted: 28 October 1998 相似文献
6.
Huang Qingtuan 《中国地震研究》2005,19(2):214-223
INTRODUCTION TheLonghai ZhangpucoastalareaofFujianProvinceliesonthesouthernsideoftheoutletofthe JiulongjiangRiver.Tectonically,itislocatedonthesouthernsegmentoftheChangle Zhao’anfault zone.Previously,alotofseismogeologicresearchworkhasbeencarriedoutinthi… 相似文献
7.
David J. Kratzmann Steven Carey Roberto Scasso Jose-Antonio Naranjo 《Bulletin of Volcanology》2009,71(4):477-439
The August 1991 eruptions of Hudson volcano produced ~2.7 km3 (dense rock equivalent, DRE) of basaltic to trachyandesitic pyroclastic deposits, making it one of the largest historical
eruptions in South America. Phase 1 of the eruption (P1, April 8) involved both lava flows and a phreatomagmatic eruption
from a fissure located in the NW corner of the caldera. The paroxysmal phase (P2) began several days later (April 12) with
a Plinian-style eruption from a different vent 4 km to the south-southeast. Tephra from the 1991 eruption ranges in composition
from basalt (phase 1) to trachyandesite (phase 2), with a distinct gap between the two erupted phases from 54–60 wt% SiO2. A trend of decreasing SiO2 is evident from the earliest part of the phase 2 eruption (unit A, 63–65 wt% SiO2) to the end (unit D, 60–63 wt% SiO2). Melt inclusion data and textures suggest that mixing occurred in magmas from both eruptive phases. The basaltic and trachyandesitic
magmas can be genetically related through both magma mixing and fractional crystallization processes. A combination of observed
phase assemblages, inferred water content, crystallinity, and geothermometry estimates suggest pre-eruptive storage of the
phase 2 trachyandesite at pressures between ~50–100 megapascal (MPa) at 972 ± 26°C under water-saturated conditions (log fO2 –10.33 (±0.2)). It is proposed that rising P1 basaltic magma intersected the lower part of the P2 magma storage region between
2 and 3 km depth. Subsequent mixing between the two magmas preferentially hybridized the lower part of the chamber. Basaltic
magma continued advancing towards the surface as a dyke to eventually be erupted in the northwestern part of the Hudson caldera.
The presence of tachylite in the P1 products suggests that some of the magma was stalled close to the surface (<0.5 km) prior
to eruption. Seismicity related to magma movement and the P1 eruption, combined with chamber overpressure associated with
basalt injection, may have created a pathway to the surface for the trachyandesite magma and subsequent P2 eruption at a different
vent 4 km to the south-southeast.
Electronic supplementary material The online version of this article (doi:) contains supplementary material, which is available to authorized users. 相似文献
8.
Following an intersection of rising magma with drifts of the potential Yucca Mountain nuclear waste repository, a pathway
is likely to be established to the surface with magma flowing for days to weeks and affecting the performance of engineered
structures located along or near the flow path. In particular, convective circulation could occur within magma-filled drifts
due to the exsolution and segregation of magmatic gas. We investigate gas segregation in a magma-filled drift intersected
by a vertical dyke by means of analogue experiments, focusing on the conditions of sustained magma flow. Degassing is simulated
by electrolysis, producing micrometric bubbles in viscous mixtures of water and golden syrup, or by aerating golden syrup,
producing polydisperse bubbly mixtures with 40% of gas by volume. The presence of exsolved bubbles induces a buoyancy-driven
exchange flow between the dyke and the drift that leads to gas segregation. Bubbles segregate from the magma by rising and
accumulating as a foam at the top of the drift, coupled with the accumulation of denser degassed magma at the base of the
drift. Steady-state influx of bubbly magma from the dyke into the drift is balanced by outward flux of lighter foam and denser
degassed magma. The length and time scales of this gas segregation are controlled by the rise of bubbles in the horizontal
drift. Steady-state gas segregation would be accomplished within hours to hundreds of years depending on the viscosity of
the degassed magma and the average size of exsolved gas bubbles, and the resulting foam would only be a few cm thick. The
exchange flux of bubbly magma between the dyke and the drift that is induced by gas segregation ranges from 1 m3 s−1, for the less viscous magmas, to 10−8 m3 s−1, for the most viscous degassed magmas, with associated velocities ranging from 10−1 to 10−9 m s−1 for the same viscosity range. This model of gas segregation also predicts that the relative proportion of erupted degassed
magma, that could potentially carry and entrain nuclear waste material towards the surface, would depend on the value of the
dyke magma supply rate relative to the value of the gas segregation flux, with violent eruption of gassy as well as degassed
magmas at relatively high magma supply rates, and eruption of mainly degassed magma by milder episodic Strombolian explosions
at relatively lower supply rates. 相似文献
9.
M. Nakagawa K. Wada T. Thordarson C. P. Wood J. A. Gamble 《Bulletin of Volcanology》1999,61(1-2):15-31
Ruapehu volcano erupted intermittently between September and November 1995, and June and July 1996, producing juvenile andesitic
scoria and bombs. The volcanic activity was characterized by small, sequential phreatomagmatic and strombolian eruptions.
The petrography and geochemistry of dated samples from 1995 (initial magmatic eruption of 18 September 1995, and two larger
events on 23 September and 11 October), and from 1996 (initial and larger eruptions on 17–18 June) suggest that episodes of
magma mixing occurred in separate magma pockets within the upper part of the magma plumbing system, producing juvenile andesitic
magma by mixing between relatively high (1000–1200 °C)- and low (∼1000 °C)- temperature (T) end members. Oscillatory zoning
in pyroxene phenocrysts suggests that repeated mixing events occurred prior to and during the 1995 and 1996 eruptions. Although
the 1995 and 1996 andesitic magmas are products of similar mixing processes, they display chronological variations in phenocryst
clinopyroxene, matrix glass, and whole-rock compositions. A comparison of the chemistry of magnesian clinopyroxene in the
four tephras indicates that, from 18 September through June 1996, the tephras were derived from at least two discrete high-temperature
(high-T) batches of magma. Crystals of magnesian clinopyroxene in the 23 September and 11 October tephras appear to be derived
from different high-T magma batches. Whole-rock and matrix-glass compositions of all tephras are consistent with their derivation
from distinct mixed melts. We propose that, prior to 1995 there was a shallow low-temperature (low-T) magma storage system
comprising crystal-rich mush and remnant magma from preceding eruptive episodes. Crystal clots and gabbroic inclusions in
the tephras attest to the existence of relict crystal mush. At least two discrete high-T magmas were then repeatedly injected
into the mush zone, forming discrete and mixed magma pockets within the shallow system. The intermittent 1995 and 1996 eruptions
sequentially tapped these magma pockets.
Received: 1 April 1998 / Accepted: 22 December 1998 相似文献
10.
11.
Results are presented from 11 microgravity surveys on Mt. Etna between 1987 and 1993, a period including the major 1989 and
1991–1993 flank eruptions and subordinate 1990 activity. Measurements were made with LaCoste and Romberg D-62 and D-157 gravity
meters along a network around the volcano between 1000 and 1900 m a.s.l. and, since 1992, a N–S summit profile. Gravity changes
of as much as 200 μGal were observed at scales from the size of the summit region to that of the volcano. None was associated
with significant changes in ground elevation. The data show an increase in gravity for 2 years before the 1989 eruption. The
increase is attributed to the accumulation of magma (0.25–1.7×109 m3) in an elongate zone, oriented NNW–SSE, between 2.5 and 6 km below sea level. Part of this magma was injected into the volcanic
pile to supply the 1989 and 1990 eruptions. It also probably fed the start of the 1991–1993 eruption, since this event was
not preceded by significant gravity changes. A large gravity increase (up to 140 μGal) detected across the volcano between
June and September 1992 is consistent with the arrival in the accumulation zone of 0.32–2.2×109 m3 of new magma, thus favoring continued flank effusion until 1993. A large gravity decrease (200 μGal) in the summit region
marked the closing stages of the 1991–1993 event and is associated with magma drainage from the upper levels of Etna's central
feeding system.
Received: 15 July 1995 / Accepted: 27 October 1997 相似文献
12.
New and detailed petrographic observations, mineral compositional data, and whole-rock vs glass compositional trends document
magma mixing in lavas erupted from Kilauea's lower east rift zone in 1960. Evidence includes the occurrence of heterogeneous
phenocryst assemblages, including resorbed and reversely zoned minerals in the lavas inferred to be hybrids. Calculations
suggest that this mixing, which is shown to have taken place within magma reservoirs recharged at the end of the 1955 eruption,
involved introduction of four different magmas. These magmas originated beneath Kilauea's summit and moved into the rift reservoirs
beginning 10 days after the eruption began. We used microprobe analyses of glass to calculate temperatures of liquids erupted
in 1955 and 1960. We then used the calculated proportions of stored and recharge components to estimate the temperature of
the recharge components, and found those temperatures to be consistent with the temperature of the same magmas as they appeared
at Kilauea's summit. Our studies reinforce conclusions reached in previous studies of Kilauea's magmatic plumbing. We infer
that magma enters shallow storage beneath Kilauea's summit and also moves laterally into the fluid core of the East rift zone.
During this process, if magmas of distinctive chemistry are present, they retain their chemical identity and the amount of
cooling is comparable for magma transported either upward or laterally to eruption sites. Intrusions within a few kilometers
of the surface cool and crystallize to produce fractionated magma. Magma mixing occurs both within bodies of previously fractionated
magma and when new magma intersects a preexisting reservoir. Magma is otherwise prevented from mixing, either by wall-rock
septa or by differing thermal and density characteristics of the successive magma batches.
Received: July 10, 1995 / Accepted: October 10, 1995 相似文献
13.
Sheet flows occur widely in natural free-surface flows including rivers in flood, tidal estuaries and coastal waters in storm conditions when bed shear stress becomes sufficiently high. Particle volumetric concentration in sheet flows normally follows a linear distribution with the Rouse [Rouse H. Modern conceptions of the mechanics of fluid turbulence. Trans ASCE, 1937;102:463–543] distribution applicable in the dilute water column above the sheet-flow layer. However, a well-verified formula for determining particle velocity distribution in sheet flows is still lacking. This paper presents formulas to describe the particle velocity profile in steady or oscillatory sheet flows. They compare well with measured data. In particular, the novel formula for determining the particle velocity at the top of bedload–sediment-dominated sublayer in sheet flows is also well verified with measured data. 相似文献
14.
The rates of passive degassing from volcanoes are investigated by modelling the convective overturn of dense degassed and
less dense gas-rich magmas in a vertical conduit linking a shallow degassing zone with a deep magma chamber. Laboratory experiments
are used to constrain our theoretical model of the overturn rate and to elaborate on the model of this process presented by
Kazahaya et al. (1994). We also introduce the effects of a CO2–saturated deep chamber and adiabatic cooling of ascending magma. We find that overturn occurs by concentric flow of the magmas
along the conduit, although the details of the flow depend on the magmas' viscosity ratio. Where convective overturn limits
the supply of gas-rich magma, then the gas emission rate is proportional to the flow rate of the overturning magmas (proportional
to the density difference driving convection, the conduit radius to the fourth power, and inversely proportional to the degassed
magma viscosity) and the mass fraction of water that is degassed. Efficient degassing enhances the density difference but
increases the magma viscosity, and this dampens convection. Two degassing volcanoes were modelled. At Stromboli, assuming
a 2 km deep, 30% crystalline basaltic chamber, containing 0.5 wt.% dissolved water, the ∼700 kg s–1 magmatic water flux can be modelled with a 4–10 m radius conduit, degassing 20–100% of the available water and all of the
1 to 4 vol.% CO2 chamber gas. At Mount St. Helens in June 1980, assuming a 7 km deep, 39% crystalline dacitic chamber, containing 4.6 wt.%
dissolved water, the ∼500 kg s–1 magmatic water flux can be modelled with a 22–60 m radius conduit, degassing ∼2–90% of the available water and all of the
0.1 to 3 vol.% CO2 chamber gas. The range of these results is consistent with previous models and observations. Convection driven by degassing
provides a plausible mechanism for transferring volatiles from deep magma chambers to the atmosphere, and it can explain the
gas fluxes measured at many persistently active volcanoes.
Received: 26 September 1997 / Accepted: 11 July 1998 相似文献
15.
Effect of Mesozoic volcanic eruptions in the western Liaoning Province,China on paleoclimate and paleoenvironment 总被引:4,自引:0,他引:4
Volatiles emitted during volcanic eruptions have an significant effect on the climate and environmental changes[1]. Different compositions of volatiles released may result in different trends of climate and environmental changes. Amounts of the erupted … 相似文献
16.
Franck Lavigne 《地球表面变化过程与地形》2004,29(8):1045-1058
Sediment yields were calculated on the ?anks of Merapi and Semeru volcanoes in Java, Indonesia, using two different methods. During the ?rst year following the 22 November 1994 eruption of Merapi, a sediment yield in excess of 1·5 × 105 m3 km?2 yr?1 was calculated in the Boyong River drainage basin, based on the volumes of sediment that were trapped by ?ve check dams. At Semeru, sediment discharges were assessed in the Curah Lengkong River from direct measurements on the lahars in motion and on the most signi?cant stream?ows. The calculated rate of sediment yield during one year of data in 2000 was 2·7 × 105 m3 km?2 yr?1. Sediment yields are dominated by rain‐triggered lahars, which occur every rainy season in several drainage basins of Merapi and Semeru volcanoes, mostly during the rainy season extending from October to April. The return period of lahars carrying sediment in excess of 5 × 105 m3 is about one year in the Curah Lengkong River at Semeru. At Merapi, the volume of sediments transported by a lahar did not exceed 2·8 × 105 m3 in the Boyong River during the rainy season 1994–95. On both volcanoes, the sediments are derived from similar sources: pyroclastic‐?ow/surges deposits, rockfalls from the lava domes, and old material from the riverbed and banks. However, daily explosions of vulcanian type at Semeru provide a more continuous sediment supply than at Merapi. Therefore, sediment yields are larger at Semeru. Copyright © 2004 John Wiley & Sons, Ltd. 相似文献
17.
18.
利用宽频带流动地震台阵和首都圈固定台网记录到的近震波形数据,研究了首都圈地区(38.6°N~41.0°N,〖JP〗115.0°E~119.7°E)的横波分裂,给出了快波偏振优势方向的场分布,讨论了首都圈地区的应力场特征. 在此基础上,采用二维线弹性有限元数值模拟方法,探讨了断层不均匀滑动对区域构造应力场的影响. 结果表明:(1)首都圈地区的应力场整体特征表现为NE向的背景应力场和受张家口-蓬莱断裂带控制的NW向的局部应力场;(2)在研究区域的西半部分和中部,最大主压应力方向为NE60°~70°,在唐山大震区及其东部区域,最大主压应力方向近WE向;(3)首都圈地区的局部应力场最大主压应力方向比较一致,基本上都与张家口-蓬莱断裂带走向平行,为120°~130°;〖JP2〗(4)首都圈区域内断层的存在及其郴均匀滑动是研究区内出现大量局部应力场的一个重要原因,张家口-〖JP〗蓬莱断裂带对首都圈局部应力场起着重要的控制作用. 相似文献
19.
对2017年九寨沟MS7.0地震序列的横波分裂的时空变化特征进行了分析.通过横波窗内S波质点运动图的分析,从九寨沟地震震源区各个地震台站的近震横波记录中提取了横波分裂的快波偏振方向和慢波延迟时间.观测结果显示,震源区各台站的上地壳各向异性在空间上存在分区特征,时间上有随时间的趋势性变化特征.空间上,位于震源区北部余震区内的3个台站中,发震断层东面的台站L5112和L5111只有一个突出的快波偏振优势方向(NNE向),而西面的L6202台有两个快波偏振优势方向(除了NNE向,还有一个近EW向),体现了余震区剧烈调整的地壳应力和构造复杂断裂的综合作用;余震区外的3个台站中,震源区东部靠近塔藏断裂(东)附近的JZG台的快波偏振优势方向为NW向,与塔藏断裂(东)的走向一致,南部的台站L5110和L5113的快波偏振优势方向为近EW向,与区域主压应力方向一致;余震区内各台站的平均慢波延迟时间大于余震区外各台站,反映了九寨沟地震孕育过程中余震区的应力积累强于其周边区域.时间上,快波偏振方向在主震后前期离散度较大,随着时间的推移,离散度在后期有逐渐变小的趋势;慢波延迟时间在主震后较大,但随着时间的推移,也表现出逐渐减小,趋于稳定.横波分裂随时间逐渐减小和趋于稳定的变化特征反映了九寨沟地震在孕震中积累的应力,随着主震和余震的发生而导致的应力释放和调整,应力大小和调整幅度逐渐减小,后期趋于稳定.
相似文献20.
During explosive eruptions the deposition of fine-grained volcanic ash fallout reduces soil permeability, favouring runoff
of meteoric water and thus increasing the occurrence of catastrophic floods. A fully dynamic, two-dimensional model was used
to simulate flooding scenarios in the Vesuvian area following an explosive volcanic eruption. The highest risk occurs in the
catchment area of the Acerra-Nola Plain N and NE of Vesuvius. This plain has a population of 70,000 living in low-lying areas.
This catchment area is vulnerable to ash fall because it lies downwind of the dominant synoptic circulation and it lacks a
natural outflow toward the sea. Our numerical simulations predict dangerous scenarios, even in quiescent periods, during extreme
rain events (return periods of 200 years have been considered), and a significant increase in the extent of the flooded areas
due to renewed volcanic activity. Based on these simulations a hazard zonation has been proposed.
Editorial responsibility: A Woods 相似文献