共查询到20条相似文献，搜索用时 15 毫秒
1.
Explosive volcanic eruptions - VI. Ejecta dispersal in plinian eruptions: the control of eruption conditions and atmospheric properties 总被引：5，自引：0，他引：5
Summary. A simple model is developed to relate the maximum down-wind and cross-wind ranges of pyroclasts forming a plinian airfall deposit to the dynamic processes in the eruption cloud from which they fall and the atmospheric wind conditions in the area. The eruption cloud dynamics are in turn related to the eruptive conditions in the vent (vent radius, exsolved magmatic volatile weight fraction, velocity with which material passes through the vent, and mass eruption rate), some or all of which can be deduced if the appropriate field measurements can be made. Some aspects of the stability of convecting volcanic eruption clouds are investigated, and the effects on eruption cloud height of the local atmospheric temperature profile and the value adopted for the entrainment constant (which relates the horizontal flow speed of atmospheric air entering the column to the vertical rise speed of the column material) are explored. It is confirmed that eruption-cloud rise height and pyroclast dispersal are mainly controlled by the mass eruption rate (per unit length of active fissure in the case of linear vents) and, hence, the heat input rate to the cloud; but a significant subsidiary dependence on the amount of exsolved magma volatiles is also found. The eruption cloud model is validated by application to observed historic eruptions, and its use in the analysis of palaeo-eruptions is discussed. 相似文献
2.
The commonest eruption styles of basaltic volcanoes involve Hawaiian lava fountaining or intermittent Strombolian explosions. We investigate the ways in which magma rise speed at depth, magma volatile content and magma viscosity control which of these eruption styles takes place. We develop a model of the degree of coalescence between gas bubbles in the magma which allows us to simulate the transition between the two extreme styles of activity. We find that magma rise speed is the most important factor causing the transition, with gas content and viscosity also influencing the rise speed at which the transition occurs. Counter to intuitive expectations, a decrease in gas content does not cause a transition from Hawaiian to Strombolian activity, but instead causes a transition to passive effusion of vesicular lava. Rather, a change from Hawaiian to Strombolian style requires a significant reduction in magma rise speed. 相似文献
3.
Massif Central (France): new constraints on the geodynamical evolution from teleseismic tomography 总被引：2，自引：0，他引：2
M. Granet G. Stoll J. Dorel U. Achauer G. Poupinet K. Fuchs 《Geophysical Journal International》1995,121(1):33-48
The Massif Central, the most significant geomorphological unit of the Hercynian belt in France, is characterized by graben structures which are part of the European Cenozoic Rift System (ECRIS) and also by distinct volcanic episodes, the most recent dated at 20 Ma to 4000 years BP. In order to study the lithosphere-asthenosphere system beneath this volcanic area, we performed a teleseismic field experiment.
During a six-month period, a joint French-German team operated a network of 79 mobile short-period seismic stations in addition to the 14 permanent stations. Inversion of P -wave traveltime residuals of teleseismic events recorded by this dense array yielded a detailed image of the 3-D velocity structure beneath the Massif Central down to 180 km depth. The upper 60 km of the lithosphere displays strong lateral heterogeneities and shows a remarkable correlation between the volcanic provinces and the negative velocity perturbations. The 3-D model reveals two channels of low velocities, interpreted as the remaining thermal signature of magma ascent following large lithospheric fractures inherited from Hercynian time and reactivated during Oligocene times. The teleseismic inversion model yields no indication of a low-velocity zone in the mantle associated with the graben structures proper. The observation of smaller velocity perturbations and a change in the shape of the velocity pattern in the 60–100 km depth range indicates a smooth transition from the lithosphere to the asthenosphere, thus giving an idea of the lithosphere thickness. A broad volume of low velocities having a diameter of about 200 km from 100 km depth to the bottom of the model is present beneath the Massif Central. This body is likely to be the source responsible for the volcanism. It could be interpreted as the top of a plume-type structure which is now in its cooling phase. 相似文献
During a six-month period, a joint French-German team operated a network of 79 mobile short-period seismic stations in addition to the 14 permanent stations. Inversion of P -wave traveltime residuals of teleseismic events recorded by this dense array yielded a detailed image of the 3-D velocity structure beneath the Massif Central down to 180 km depth. The upper 60 km of the lithosphere displays strong lateral heterogeneities and shows a remarkable correlation between the volcanic provinces and the negative velocity perturbations. The 3-D model reveals two channels of low velocities, interpreted as the remaining thermal signature of magma ascent following large lithospheric fractures inherited from Hercynian time and reactivated during Oligocene times. The teleseismic inversion model yields no indication of a low-velocity zone in the mantle associated with the graben structures proper. The observation of smaller velocity perturbations and a change in the shape of the velocity pattern in the 60–100 km depth range indicates a smooth transition from the lithosphere to the asthenosphere, thus giving an idea of the lithosphere thickness. A broad volume of low velocities having a diameter of about 200 km from 100 km depth to the bottom of the model is present beneath the Massif Central. This body is likely to be the source responsible for the volcanism. It could be interpreted as the top of a plume-type structure which is now in its cooling phase. 相似文献
4.
Magnetic history of a dyke on Mount Etna (Sicily) 总被引：1，自引：0，他引：1
During the 1989 eruption of Mount Etna, two fracture systems, trending c. N45°E and N150°E, opened at the foot of its 3000 m high SE Crater and propagated quickly downslope to distances of ≈3 and 7 km, respectively. The northeastern fracture fed a flank eruption, whereas the southeastern fracture remained dry and offered contrasting volcanological and geophysical evidence of the presence of magma at a shallow depth. During the opening of this non-eruptive fracture system, a differential magnetic network was set up on a short profile across its distant extremity. Initially, the magnetic field did not display any change along the profile between frequent surveys. However, repeated measurements at intervals of about 3 months for two years revealed the slow build-up of a 130 nT anomaly. The anomaly vanishes laterally within 0.2 km of the surface expression of the fracture system. This exceptional set of observations constrains the location and time of cooling of a shallow dyke. The increase in magnetization of the dyke inferred by the rate of growth of the anomaly leads to the interpretation that the dyke was emplaced near the end of the eruption. 相似文献
5.
Vyacheslav M. Zobin Carlos Navarro Gabriel Reyes-Dávila Justo Orozco Mauricio Bretón Armando Tellez Gabriel Reyes-Alfaro Homero Vázquez 《Geophysical Journal International》2006,167(1):467-478
A methodology is proposed for the quantification of volcanic explosions based on three parameters derived from broad-band seismic signals: the counter force of the eruption F , the power of the explosion P and the duration of the upward movement of the gas slug in the conduit to the free surface of magma, D . This methodology was applied to the 2004–2005 sequence of explosions at Volcán de Colima, Mexico. The broad-band records of more than 100 explosive events were obtained at a distance of 4 km from the crater. We determined the counter force of the eruption by modelling the low-frequency impulse of the seismic records of 66 volcanic explosions and estimated the power of 116 explosions from the spectra of the high-frequency impulse. The power of Colima explosions spans five orders of magnitude; the counter force spans four orders of magnitude. We show that the power of a volcanic explosion is proportional to the counter force of the eruption. These parameters may be used for the elaboration of a scale of volcanic explosions. 相似文献
6.
Magnetic anisotropy produced by magma flow: theoretical model and experimental data from Ferrar dolerite sills (Antarctica) 总被引：3，自引：0，他引：3
Volcanic rocks forming sills, dykes or lava flows may display a magnetic anisotropy derived from the viscous flow during their emplacement. We model a sill as a steady-state flow of a Bingham fluid, driven by a pressure gradient in a horizontal conduit. The magma velocity as a function of depth is calculated from the motion and constitutive equations. Vorticity and strain rate are determined for a reference system moving with the fluid. The angular velocity and the orientation of an ellipsoidal magnetic grain immersed in the fluid are calculated as functions of time or strain. Magnetic susceptibility is then calculated for a large number of grains with a uniform distribution of initial orientations. It is shown that the magnetic lineation oscillates in the vertical plane through the magma flow direction, and that the magnetic foliation plane changes periodically from horizontal to vertical. The results are compared with the magnetic fabric of Ferrar dolerite sills (Victoria Land, East Antarctica) derived from low-field susceptibility measurements. 相似文献
7.
8.
Lionel Wilson 《Geophysical Journal International》1999,136(3):609-619
Previous treatments of the relationship between the mass fraction of released magma volatiles and the eruption speeds of gas and pyroclasts in steady explosive eruptions have not taken detailed account of the dynamic effects associated with the finite size distribution of the pyroclasts. When this is done, it is found that previously published estimates of exsolved magma volatile contents obtained from the analysis of pyroclast size distributions in near-vent deposits overestimate the volatile content by approximately 20 per cent in the case of Plinian eruptions. The discrepancy is much worse for pyroclast size distributions skewed towards coarse clasts, as is common in basaltic lava fountains; in such cases pyroclast dispersal studies may overestimate the exsolved magma volatile content by at least 200 per cent. An analogous problem arises if released magma volatile amounts deduced from studies of fluid inclusions in pyroclasts are inserted into most current computer models of eruption dynamics, but the gas eruption speeds deduced have an even larger error, being underestimated by up to 300 per cent in the case of coarse-grained eruptions. The more sophisticated of the currently available numerical models of eruption dynamics can in principle avoid this problem, but practical implementation limitations have so far prevented such models being run with a sufficiently wide range of grain sizes for the importance of these effects to be fully appreciated. 相似文献
9.
戈壁风蚀面与植被覆盖面地表性质粗糙度长度的确定 总被引：9，自引：5，他引：4
以Monin-Obukhov相似性理论为基础, 利用量纲分析法分别推导出不同层结稳定度下确定戈壁风蚀面与植被覆盖面空气动力学参数的物理模型, 并利用该模型研究了粗糙度长度与粗糙元性质, 流经近地层流体特征以及大气层结稳定度之间的关系。得出以下结论: 戈壁风蚀面上空气动力学粗糙度长度与砾石粒径、高度、覆盖度、自由风速、摩擦速度以及大气层结稳定度有关; 植被覆盖面空气动力学粗糙度长度取决于植被类型、植被高度、覆盖度、零平面位移高度、自由风速、摩擦速度以及大气层结稳定度。 相似文献
10.
Summary During the 1979 eruption of the Soufrière of St Vincent, West Indies, a 14 min period of explosive activity occurred on April 22. The first three minutes of this activity were filmed. Five individual explosions formed distinct plumes which fed an eruption column which eventually ascended to a height of over 18 km. The ascent velocities of the fronts of these plumes ranged from 8.5 to 61.7 ms-1 . The half-angle of spreading of the plume fronts ranged from 21.5 to 24°. One of the plumes was observed to 8 km height and decelerated steadily from 51.5 to 23 ms-1 . The main column fed by individual explosions was, however, observed to accelerate from 8.5 to 58 m s-1 in the same height interval.
A theoretical treatment of volcanic plume motion is presented. Measurements of dimensions, velocities and spreading rates from the film are used to estimate plume parameters such as temperature, particle content and volume discharge rate of magma from the theoretical relationships. These calculations show that the individual plumes became hotter and richer in juvenile ash with time. The acceleration of the main eruption column was the result of being fed by increasingly hotter and more ash-rich explosions. An average volume discharge rate of 12600m3 s-1 is estimated from the analysis of the plume motions. This value agrees closely with an estimate of discharge rate based on the heat flux required to form an 18 km high column. This agreement suggests that theories of convective motions in plumes can be successfully adapted to the volcanic case, as long as the effect of hot entrained ash particles is considered. 相似文献
A theoretical treatment of volcanic plume motion is presented. Measurements of dimensions, velocities and spreading rates from the film are used to estimate plume parameters such as temperature, particle content and volume discharge rate of magma from the theoretical relationships. These calculations show that the individual plumes became hotter and richer in juvenile ash with time. The acceleration of the main eruption column was the result of being fed by increasingly hotter and more ash-rich explosions. An average volume discharge rate of 12600m
11.
Roger G. Newton 《Geophysical Journal International》1982,69(2):571-572
Summary It is pointed out that the data on which the Marchenko inversion of the second-order equation for fluid media reconstructs the impedance are not all independent. In particular, the derivative Z' 0 of the surface impedance is implied by the displacement velocity response at the surface. We also note that unless Z' 0 = 0 the zero-frequency reflection amplitude must always equal — 1. 相似文献
12.
Effects of fractures on seismic-wave velocity and attenuation 总被引：1，自引：0，他引：1
The effects of fractures on the seismic velocity and attenuation of a rock are investigated using theoretical results and experimental data. Fractures in a rock mass influence the traveltimes and amplitudes of seismic waves that have propagated through them. The displacement discontinuity model, recently employed in fracture investigations, is modified to describe the effect of fractures on seismic-wave velocity and attenuation. This new model, the modified displacement discontinuity model (MDD), is formulated in a way analogous to transmission-line analysis. The fractures are treated as transmission lines for the passage of seismic waves. The MDD takes into consideration realistic fracture parameters which include the fracture length, the fractional area of a fracture surface in contact, and the nature of the infilling material. A single fracture of varying geometric and material properties is shown to affect dramatically the transmission properties of a propagating waveform, and hence the seismic velocity and attenuation. These effects have been shown to result in a frequency-dependent velocity and attenuation. The sensitivity of the fracture parameters to seismic-wave velocity and attenuation was investigated and interesting results were obtained. Fracture parameters used in designing experimental models consisting of synthetically manufactured cracks were fed into the MDD and a well-known crack model, Hudson's model, for comparison. Velocities as a function of the incident-wave angle were obtained from both numerical models and were compared with the results from the experimental modelling. For P waves, the MDD model results show better agreement with those of the experimental model for all crack densities investigated than those from Hudson's model. 相似文献
13.
Anomalous seismic crustal structure of oceanic fracture zones 总被引：2，自引：0，他引：2
Robert S. White Robert S. Detrick Martin C. Sinha Marie H. Cormier 《Geophysical Journal International》1984,79(3):779-798
Summary. The seismic structure of crust found within fracture zones falls outside the range of velocity structures observed for normal oceanic crust in the North Atlantic. The crust in fracture zones is frequently very thin and is characterized by low crustal velocities and by the conspicuous absence of a refractor with a velocity typical of oceanic layer 3. Anomalous crust is present in both large- and small-offset fracture zones. Since they are among the most common tectonic features in the ocean basins, and are particularly closely spaced on slow-spreading ridges, fracture zones represent a major source of seismic crustal heterogeneity. We interpret the anomalous crust as a thin, intensely fractured, faulted and hydrothermally altered basaltic and gabbroic section overlying ultramafics that, in places, are extensively serpentinized. The unusually thin crust found within fracture zones and the gradual crustal thinning over a distance of several tens of kilometres on either side of the fracture zones can be explained by two main processes; firstly the cold lithosphere edge opposite the spreading centre at the ridgetransform intersection modifies the normal intrusive and extrusive processes of the spreading centre leading to the accretion of an anomalous and thin igneous section; and secondly each spreading ridge segment is fed from a separate subcrustal magma supply point, so as the magma flows laterally down the spreading centre it generates a crustal section of decreasing thickness, culminating in the very thin crust of the fracture zones at either end of the ridge segment. 相似文献
14.
B. Di Lieto G. Saccorotti L. Zuccarello M. La Rocca R. Scarpa 《Geophysical Journal International》2007,169(2):699-705
Using data from two dense array of short period seismometers, we analyse the kinematic properties of volcanic tremor preceding and accompanying the 2004–2005 eruption of Etna Volcano, Italy. Results from slowness analyses indicate the action of at least two distinct sources. The first dominates the pre-eruptive period, and is likely associated with the main plumbing system feeding the Summit and southeast craters. Following the onset of the eruption, secondary directions of wave-arrival encompass the eruptive fissures, located on the lower eastern flank of the southeast crater. Nonetheless, significant energy radiation from this latter source was also occurring prior to the onset of the lava effusion, likely suggesting the presence of a resident magma batch, in agreement with independent petrologic and geochemical data. 相似文献
15.
Significant faulting and deformation of the ground surface has been rarely known during volcanic eruptions. Usu Volcano, Hokkaido, Japan, is a unique example of deformation due to felsic magma intrusion. Usu Volcano has a history of such types of eruptions as phreatic, pumice eruption (Plinian type), pyroclastic flowing and lava doming since 1663. On March 31, 2000, phreatomagmatic to phreatic eruptions took place after 23 years of dormancy in the western piedmont, followed by explosions on the western flank of Usu Volcano. They were associated with significant deformation including faulting and uplift. The eruptions and deformation were continuing up to the end of May 2000.We identified the faulting using total nine sets of aerial photographs taken from before the eruption (March 31, 2000) to more than 1 year (April 27, 2001) after the end of the activity, and traced deformation processes through image processing using aerial photographs. We found that some of the new faults and the associated phreatic eruptions were related to old faults formed during the 1977–1981 eruptive episode.The image processing has revealed that the surface deformation is coincident with the area of faulting forming small grabens and the phreatic explosion vents. However, the faulting and main explosive eruptions did not take place in the highest uplift area, but along the margin. This suggests that the faulting and explosive activities were affected by small feeder channels diverging from the main magma body which caused the highest uplift. 相似文献
16.
Open fractures provide high-permeability pathways for fluid flow in sedimentary basins. The potential for flow along permeable or open fractures and faults depends on the continuity of flow all the way to the surface except in the case of convective flow. Upward flowing fluid cools and may cause cementation due to the prograde solubility of quartz, but in the case of carbonates such flow may cause dissolution. The rate and duration of these processes depend on the mechanisms for sustaining fluid flow into the fracture, the geometries of fracture and sedimentary beds intersected, permeability, pressure and temperature gradients. Heat loss to the adjacent sediments causes sloping isotherms which can induce non-Rayleigh convection. To analyse these problems we have used a simple model in which a single fracture acts as a pathway for vertically moving fluid and there is no fluid transport across the walls of the fracture except near its inlet and outlet. Four mechanisms for fluid flow into the lower part of the fracture are considered: decompression of pore water; compaction of intersected overpressared sediments; focusing of compaction water derived from sediments beneath the fracture; and finally focusing of pore water moving through an aquifer. Water derived from the basement is not considered here. We find that sustained flow is unlikely to have velocities much higher than 1–100 m/yr, and the flow is laminar. The temperature of the fluid expelled at the top of the fracture increases by less than 1% and the vertical temperature gradient in the fracture remains close to the geothermal gradient. Where hot water is introduced from basement fractures (hydrothermal water) during tectonic deformation, much higher velocities may be sustained in the overlying sediments, but here also this depends on the permeability near the surface. Most of the cooling of water with (ore) mineral precipitation will then occur near the surface. In most cases, pore water decompression and sediment compaction will yield only very limited pore water flux with no significant potential for cementation or heating of the sediments adjacent to the fracture. Focusing of compaction water from sediments beneath the fracture or from an intersected aquifer can yield fluxes high enough to cement an open fracture significantly but the flow must be sustained for a very long time. For velocities of 1–100 m/yr, it takes typically 0.3–30 Myr to cement a fracture by 50%. The highest velocities may be obtained when a fracture extends all the way to the surface or sea floor. When a fracture does not reach the sediment surface, the flow velocity is reduced by the displacement of water in the sediments near the top of the fracture. The flow into the fracture from the sediments may often be rate limiting rather than the flow on the fracture. Sedimentary rocks only a few metres from the fracture will receive a much lower flux than the fracture. The fracture will therefore close due to cementation before significant amounts of silica can be introduced into adjacent sandstones. The isotherm slope in the adjacent sediments will in most cases be less than 10–20°. Non-Rayleigh convection velocities in the sediments adjacent to the fracture are too small to cause any significant diagenetic reactions such as quartz cementation. These quantifications of fluid flow in fractures in sedimentary basins are important in terms of constraining models for diagenesis, heat transport and formation of ore minerals in a compaction-driven system. 相似文献
17.
Explosive eruptions of mafic magmas produce lava fountains whose heights are a function of the exsolved volatile content of the magma, its erupted mass flux, and the geometry of the vent (which may be an elongate fissure or a localized, near-circular conduit). The geometry of the initial vent (and the eruptive behaviour) can be distinctly modified by lava drainback and accumulating ejecta. Hot pyroclasts landing near the vent may coalesce to form rootless flows, some of which may drain back into the vent to be recycled into the eruption products. Rootless flows may be at least partially confined by pre-existing topographic features, or by spatter or cinder ramparts being built up by the eruption itself, so that they accumulate into a lava pond over and around the vent. The erupting jet of magmatic gas and pyroclasts must force its way through such a pond and will entrain some of the pond lava as it does so. The energy expended in entraining and accelerating previously erupted materials will reduce the eruption velocity and the lava fountain height by an amount which can be calculated as a function of the eruption conditions and the lava pond depth (or lava drainback rate). The results of such calculations are presented, and are used to assess the influence of this process on attempts to infer magma volatile contents from field observations of lava fountain heights. 相似文献
18.
19.
J. Julià C. J. Ammon R. B. Herrmann A. M. Correig 《Geophysical Journal International》2000,143(1):99-112
We implement a method to invert jointly teleseismic P wave receiver functions and surface wave group and phase velocities for a mutually consistent estimate of earth structure. Receiver functions are primarily sensitive to shear wave velocity contrasts and vertical traveltimes, and surface wave dispersion measurements are sensitive to vertical shear wave velocity averages. Their combination may bridge resolution gaps associated with each individual data set. We formulate a linearized shear velocity inversion that is solved using a damped leastsquares scheme that incorporates a priori smoothness constraints for velocities in adjacent layers. The data sets are equalized for the number of data points and physical units in the inversion process. The combination of information produces a relatively simple model with a minimal number of sharp velocity contrasts. We illustrate the approach using noisefree and realistic noise simulations and conclude with an inversion of observations from the Saudi Arabian Shield. Inversion results for station SODA, located in the Arabian Shield, include a crust with a sharp gradient near the surface (shear velocity changing from 1.8 to 3.5 km s1 in 3 km) underlain by a 5kmthick layer with a shear velocity of 3.5 km s1 and a 27kmthick layer with a shear velocity of 3.8 km s1 , and an upper mantle with an average shear velocity of 4.7 km s1 . The crustmantle transition has a significant gradient, with velocity values varying from 3.8 to 4.7 km s1 between 35 and 40 km depth. Our results are compatible with independent inversions for crustal structure using refraction data. 相似文献
20.
Rosaria Tondi Roberto de Franco & Riccardo Barzaghi' 《Geophysical Journal International》2000,141(3):679-698
A new algorithm is presented for the integrated 2-D inversion of seismic traveltime and gravity data. The algorithm adopts the 'maximum likelihood' regularization scheme. We construct a 'probability density function' which includes three kinds of information: information derived from gravity measurements; information derived from the seismic traveltime inversion procedure applied to the model; and information on the physical correlation among the density and the velocity parameters. We assume a linear relation between density and velocity, which can be node-dependent; that is, we can choose different relationships for different parts of the velocity–density grid. In addition, our procedure allows us to consider a covariance matrix related to the error propagation in linking density to velocity. We use seismic data to estimate starting velocity values and the position of boundary nodes. Subsequently, the sequential integrated inversion (SII) optimizes the layer velocities and densities for our models. The procedure is applicable, as an additional step, to any type of seismic tomographic inversion.
We illustrate the method by comparing the velocity models recovered from a standard seismic traveltime inversion with those retrieved using our algorithm. The inversion of synthetic data calculated for a 2-D isotropic, laterally inhomogeneous model shows the stability and accuracy of this procedure, demonstrates the improvements to the recovery of true velocity anomalies, and proves that this technique can efficiently overcome some of the limitations of both gravity and seismic traveltime inversions, when they are used independently.
An interpretation of field data from the 1994 Vesuvius test experiment is also presented. At depths down to 4.5 km, the model retrieved after a SII shows a more detailed structure than the model obtained from an interpretation of seismic traveltime only, and yields additional information for a further study of the area. 相似文献
We illustrate the method by comparing the velocity models recovered from a standard seismic traveltime inversion with those retrieved using our algorithm. The inversion of synthetic data calculated for a 2-D isotropic, laterally inhomogeneous model shows the stability and accuracy of this procedure, demonstrates the improvements to the recovery of true velocity anomalies, and proves that this technique can efficiently overcome some of the limitations of both gravity and seismic traveltime inversions, when they are used independently.
An interpretation of field data from the 1994 Vesuvius test experiment is also presented. At depths down to 4.5 km, the model retrieved after a SII shows a more detailed structure than the model obtained from an interpretation of seismic traveltime only, and yields additional information for a further study of the area. 相似文献