首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 15 毫秒
1.
Volcano seismology often deals with rather shallow seismic sources and seismic stations deployed in their near field. The complex stratigraphy on volcanoes and near‐field source effects have a strong impact on the seismic wavefield, complicating the interpretation techniques that are usually employed in earthquake seismology. In addition, as most volcanoes have a pronounced topography, the interference of the seismic wavefield with the stress‐free surface results in severe waveform perturbations that affect seismic interpretation methods. In this study we deal predominantly with the surface effects, but take into account the impact of a typical volcano stratigraphy as well as near‐field source effects. We derive a correction term for plane seismic waves and a plane‐free surface such that for smooth topographies the effect of the free surface can be totally removed. Seismo‐volcanic sources radiate energy in a broad frequency range with a correspondingly wide range of different Fresnel zones. A 2‐D boundary element method is employed to study how the size of the Fresnel zone is dependent on source depth, dominant wavelength and topography in order to estimate the limits of the plane wave approximation. This approximation remains valid if the dominant wavelength does not exceed twice the source depth. Further aspects of this study concern particle motion analysis to locate point sources and the influence of the stratigraphy on particle motions. Furthermore, the deployment strategy of seismic instruments on volcanoes, as well as the direct interpretation of the broad‐band waveforms in terms of pressure fluctuations in the volcanic plumbing system, are discussed.  相似文献   

2.
We have developed a new array method combining conventional migration with a slowness-backazimuth deviation weighting scheme. All seismic traces are shifted based on the theoretical traveltime of the scattered wave from specific gridpoints in a 3-D volume. Observed slowness and backazimuth are calculated for each raypath and compared with theoretical values in order to estimate slowness and backazimuth deviations. Subsequently, stacked energy calculated by a conventional migration method is weighted by the slowness and backazimuth deviations to suppress any arrival energy whose slowness and backazimuth are inconsistent with the expected theoretical values. This new method was applied to two P- wave data sets which comprise (1) underside reflections at the 410 and 660 km mantle discontinuities and (2) D" reflections as well as their corresponding synthetic data sets. The results show that the weighting scheme dramatically increases the resolution of the migrated images and enables us to obtain well-constrained, focused images, making upper-mantle discontinuities and D" reflections more distinct by reducing their surrounding energy.  相似文献   

3.
Generalized Born scattering of elastic waves in 3-D media   总被引：1，自引：0，他引：1  
It is well known that when a seismic wave propagates through an elastic medium with gradients in the parameters which describe it (e.g. slowness and density), energy is scattered from the incident wave generating low-frequency partial reflections. Many approximate solutions to the wave equation, e.g. geometrical ray theory (GRT), Maslov theory and Gaussian beams, do not model these signals. The problem of describing partial reflections in 1-D media has been extensively studied in the seismic literature and considerable progress has been made using iterative techniques based on WKBJ, Airy or Langer type ansätze. In this paper we derive a first-order scattering formalism to describe partial reflections in 3-D media. The correction term describing the scattered energy is developed as a volume integral over terms dependent upon the first spatial derivatives (gradients) of the parameters describing the medium and the solution. The relationship we derive could, in principle, be used as the basis for an iterative scheme but the computational expense, particularly for elastic media, will usually prohibit this approach. The result we obtain is closely related to the usual Born approximation, but differs in that the scattering term is not derived from a perturbation to a background model, but rather from the error in an approximate Green's function. We examine analytically the relationship between the results produced by the new formalism and the usual Born approximation for a medium which has no long-wavelength heterogeneities. We show that in such a case the two methods agree approximately as expected, but that in a media with heterogeneities of all wavelengths the new gradient scattering formalism is superior. We establish analytically the connection between the formalism developed here and the iterative approach based on the WKBJ solution which has been used previously in 1-D media. Numerical examples are shown to illustrate the examples discussed.  相似文献   

4.
We study properties of the energy-flux vector and other related energy quantities of homogeneous and inhomogeneous time-harmonic P and S plane waves, propagating in unbounded viscoelastic anisotropic media, both analytically and numerically. We propose an algorithm for the computation of the energy-flux vector, which can be used for media of unrestricted anisotropy and viscoelasticity, and for arbitrary homogeneous or inhomogeneous plane waves. Basic part of the algorithm is determination of the slowness vector of a homogeneous or inhomogeneous wave, which satisfies certain constraints following from the equation of motion. Approaches for determination of a slowness vector commonly used in viscoelastic isotropic media are usually difficult to use in viscoelastic anisotropic media. Sometimes they may even lead to non-physical solutions. To avoid these problems, we use the so-called mixed specification of the slowness vector, which requires, in a general case, solution of a complex-valued algebraic equation of the sixth degree. For simpler cases, as for SH waves propagating in symmetry planes, the algorithm yields simple analytic solutions. Once the slowness vector is known, determination of energy flux and of other energy quantities is easy. We present numerical examples illustrating the behaviour of the energy-flux vector and other energy quantities, for homogeneous and inhomogeneous plane P , SV and SH waves.  相似文献   

5.
Decay times inferred from relative sea‐level (RSL) histories of previously glaciated regions provide a potentially important constraint on mantle rheology. We present a new compilation of RSL data from Richmond Gulf and James Bay, Canada. This recompilation reveals errors in previous compilations that led to inaccurate estimates for the Richmond Gulf decay time in a series of recently published articles. We derive updated estimates for the decay time at Richmond Gulf and James Bay using a methodology that incorporates errors in both the age and the height of the sea‐level markers. This exercise is guided by a series of synthetic RSL calculations that show that decay time estimates in the region can be significantly biased if the RSL time‐series are not corrected for global eustatic sea‐level trends, or if the estimates are based on composite RSL histories derived by combining data from both the Richmond Gulf and the James Bay regions. Our decay time analysis for Richmond Gulf applies the pioneering approach of Walcott (1980) to a large database and we derive a value of 4.0–6.6 kyr, where the range is defined by a misfit tolerance 10 per cent higher than the minimum. Our analysis for James Bay is based on the uplift curve derived by Hardy (1976) , and we estimate a decay time of about 2.0–2.8 kyr. The difference between our estimates for Richmond Gulf and James Bay may be due to errors in the observational record from these regions, but could also be influenced by lateral variations in lithospheric structure associated with the assembly of Laurentia.  相似文献   

6.
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.  相似文献   

7.
Summary . Synthetic seismograms represented by integrals generally display signals associated with the limits of integration. Sometimes these 'end-point' contributions are spurious (e.g. in the WKBJ seismogram) and sometimes they are the main physical interest (e.g. the Kirchhoff integral for an edge). The end-point contributions may be asymptotically approximated using integration by parts or Laplace's method and it may then be possible to reduce them if desired. We describe examples in the WKBJ seism ogram for reflected or transmitted waves in homogeneous layers and for turning waves. We also study signals due to discontinuities in reflection coefficients, by partitioning the real slowness integral so that the discontinuities lie at end points. Examples are the head wave, which is a physically correct signal, and spurious diffractions caused by using plane-wave coefficients for grazing rays in the WKBJ seismogram.  相似文献   

8.
ABSTRACT

In the mid-latitudes, snow plays a critical role in regional hydroclimate, with snow ablation variability in ephemeral regions representing an area of essential research. Due to a lack of historical snow-water-equivalent data in the eastern United States, recent research has substituted daily snow depth changes for ablation. These studies, however, do not explicitly examine if such a substitution yields a snowmelt hydrological signal, an important component of water resource management. As such, this study evaluates if ablation events, as defined as a daily snow depth decrease, subsequently result in increased river discharge within two similarly sized watersheds in the eastern United States: the Wabash and Susquehanna River basins. For both basins, >75% of snow ablation events resulted in a positive river discharge response (increase in discharge) at a 3-day lag. Furthermore, results show a significant and positive relationship between ablation event frequency and seasonal discharge response, such that an increase (decrease) in seasonal snow ablation event frequency yields an increase (decrease) in associated seasonal river discharge at a 3-day lag. These relationships indicate that inter-diurnal decreases in snow depth do carry hydrological implications, adding confidence that such a definition of ablation is appropriate for climatological applications.  相似文献   

9.
Summary. The asymptotic properties of spheroidal mode dispersion at high frequency for fixed phase velocity are related to the intercept times τβ( p ) for P and S waves. If the mode eigenfrequency and the ratio of horizontal to vertical displacement at the surface for the mode are known τα( p ) and τβ( p ) may be separately estimated. If discontinuities exist in the velocity model then 'solotone' effects occur, in frequency at fixed slowness, and in τα( p ), τβ( p ) estimated from the mode dispersion as a function of slowness. The coupling of P and S waves in the spheroidal modes means that the interaction of P waves with upper-mantle discontinuities affects also the estimates of the S wave τβ( p ) values for which the corresponding turning points lie in the lower mantle. The asymptotic formalism also shows that sharp pulses formed by superposition of spheroidal modes correspond to multiple PS reflections.
A study of τα( p ), τβ( p ) estimates derived from spheroidal modes with periods from 45–50s, calculated for model 1066B, shows that even in the presence of strong upper-mantle discontinuities the errors in intercept time are only about one-tenth of a period. The asymptotic properties may there-for provide a useful means of estimating intercept times from modes with a few seconds period as a supplement to travel-time methods.  相似文献   

10.
The main goal of this study is to improve the modelling of the source mechanism associated with the generation of long period (LP) signals in volcanic areas. Our intent is to evaluate the effects that detailed structural features of the volcanic models play in the generation of LP signal and the consequent retrieval of LP source characteristics. In particular, effects associated with the presence of topography and crustal heterogeneities are here studied in detail. We focus our study on a LP event observed at Kilauea volcano, Hawaii, in 2001 May. A detailed analysis of this event and its source modelling is accompanied by a set of synthetic tests, which aim to evaluate the effects of topography and the presence of low velocity shallow layers in the source region. The forward problem of Green's function generation is solved numerically following a pseudo-spectral approach, assuming different 3-D models. The inversion is done in the frequency domain and the resulting source mechanism is represented by the sum of two time-dependent terms: a full moment tensor and a single force. Synthetic tests show how characteristic velocity structures, associated with shallow sources, may be partially responsible for the generation of the observed long-lasting ringing waveforms. When applying the inversion technique to Kilauea LP data set, inversions carried out for different crustal models led to very similar source geometries, indicating a subhorizontal cracks. On the other hand, the source time function and its duration are significantly different for different models. These results support the indication of a strong influence of crustal layering on the generation of the LP signal, while the assumption of homogeneous velocity model may bring to misleading results.  相似文献   

11.
Simulation of seismicity due to fluid migration in a fault zone   总被引：5，自引：0，他引：5  
Spatio-temporal variation of rupture activity is modelled assuming fluid migration in a narrow, porous fault zone formed along a vertical strike-slip fault in a semi-infinite elastic medium. The principle of the effective stress coupled to the Coulomb failure criterion introduces mechanical coupling between fault slip and the pore fluid. The fluid is assumed to flow out of a localized high-pressure fluid compartment in the fault at the onset of earthquake rupture. The duration of the earthquake sequence is assumed to be much shorter than the recurrence period of characteristic events on the fault. Both an earthquake swarm and a foreshock–main-shock sequence can be simulated by changing the relative magnitudes of the initial tectonic stress, pore fluid pressure, fracture strength and so on. When an inhomogeneity is introduced into the spatial distribution of fracture strength, high complexity is observed in the spatio-temporal variation of rupture activity. For example, the time interval between two successive events is highly irregular, and a relatively long quiescence of activity is sometimes observed in a foreshock–main-shock sequence. The quiescence is caused by the temporary arresting of rupture extension, due to an encounter with fault segments having locally high strengths. The frequency–magnitude statistics of intermediate-size events obey the Gutenberg–Richter relation. The calculations show the temporal variation of the b value during some foreshock sequences, and the degree of the change seems to depend on the statistical distribution of the fracture strength.  相似文献   

12.
ABSTRACT Zircon is the most widely used mineral in detrital dating studies because it is common to multiple rock types, is chemically and physically resistant, and can endure successive cycles of burial, metamorphism and erosion. Zircon also has the advantage that single grains may be dated by either the fission‐track (FT) or U–Pb method, which, because of their contrasting thermal sensitivities (total resetting occurs at temperatures > 320 °C for FT and > 700–900 °C for U–Pb), can provide unique information about both the age structure and the thermal evolution of a sediment source. However, single method‐based bias and difficulties associated with interpreting measured ages can influence both the quality and the level of useful provenance information. For example, the zircon FT system is sensitive to metamorphic overprinting and hence measured ages alone cannot be interpreted as unambiguously dating formation age of the source rock. In contrast, U–Pb zircon data have high resistance (700–900 °C) to thermal overprinting and therefore recorded formation ages may not relate to an immediate source but may instead reflect a polycyclical history. The focus of this paper is to examine, from a practical standpoint, the provenance potential of detrital zircon fission track data and to investigate the method's complementary role as an aide to the interpretation of high‐temperature detrital U–Pb zircon data by combining U–Pb and FT methods in a single study.  相似文献   

13.
Surface mass redistribution within the Earth system, especially in the atmosphere, oceans, continents and ice sheets, causes the position of the centre of mass to vary in a reference frame attached to the solid Earth. Space techniques are now precise enough to measure the centre of mass motion. Here we present a determination of the centre of mass coordinates at regular monthly intervals using DORIS data on SPOT‐2, SPOT‐3 and Topex–Poseidon (1993–1997) and laser data on Lageos‐1 and Lageos‐2 (1993–1996). The amplitude and phase of the space‐geodesy‐derived annual cycle for each coordinate are further compared to estimates based on surface mass redistribution at the Earth surface derived from various climatic data sources: surface pressure, soil moisture, snow depth and ocean mass variations.  相似文献   

14.
Multivariate numerical analyses (DCA, CCA) were used to study the distribution of chironomids from surface sediments of 100 lakes spanning broad ecoclimatic conditions in northern Swedish Lapland. The study sites range from boreal forest to alpine tundra and are located in a region of relatively low human impact. Of the 19 environmental variables measured, ordination by CCA identified mean July air temperature as one of the most significant variables explaining the distribution and the abundance of chironomids. Lossonignition (LOI), maximum lake depth and mean January air temperature also accounted for significant variation in chironomid assemblages. A quantitative transfer function was created to estimate mean July air temperature from sedimentary chironomid assemblages using weightedaveraging partial least squares regression (WAPLS). The coefficient of determination was relatively high (r2 = 0.65) with root mean squared error of prediction (RMSEP, based on jack-knifing) of 1.13 °C and maximum bias of 2.1 °C, indicating that chironomids can provide useful quantitative estimates of past changes in mean July air temperature. The paper focuses mainly on the relationship between chironomid composition and July air temperature, but the relationship to LOI and depth are also discussed.  相似文献   

15.
Summary. In Part I of this paper we modelled shear-wave splitting observed in crystalline rock bordering an active, normal fault-zone at Oroville, California, with Červený's ray-tracing system applied to anisotropic heterogeneous media using Hudson's formulation of elastic constants for a medium containing aligned cracks. In Part II we use the ray-tracing results of Part I to quantitatively interpret P -wave polarization anomalies observed in the three-component seismograms recorded in the Oroville fault zone. We show that the eigenvectors of the first-order Christoffel tensor defined by the ray-tracing slowness vector and Hudson's first-order anisotropic corrections to the isotropic elastic tensor correctly account for P -wave first motion that deviates from the ray vector.  相似文献   

16.
Summary. The transformation of a set of seismograms to the delay time-slowness, τ—p, domain is presented as a sequence of Fourier and Bessel transforms, For a horizontally layered medium, this sequence gives an exact cylindrical wave decomposition of the response to a point source; correctly compensating for the phase shifting and geometrical spreading associated with transmission through the Earth. The resultant τ—p map or 'slant stack' contains true amplitude and phase information. The spatial aliasing properties of the transformation, when applied to a dataset, are greatly improved by the use of only outgoing waves in the Bessel transform. This is equivalent to using Hankel functions rather than Bessel functions, and is justified by the absence of incoming waves from most datasets. The WKBJ approximation to the medium response enables predictions to be made about the shape and amplitude variation with slowness of truncation effects. Theoretically the τ—p transformation is reversible, thus the τ—p domain is a suitable one in which to perform filtering operations before seismogram reconstruction.  相似文献   

17.
Katherine Donovan 《Area》2010,42(1):117-126
Social volcanology refers to the integration of social science research methods into the traditionally physical domain of volcanology. This emerging multi-methodological research area draws from many disciplines in order to examine hazard-mitigation strategies that are community focused. A key facet of social volcanology is the role of culture and this paper explores the influence of traditional cultural values in relation to the 2006 volcanic crisis at Mt Merapi (Java). This paper describes the complex amalgam of cultural and socio-economic factors that influence community reactions to volcanic hazards and demonstrates the need for interdisciplinary hazard research.  相似文献   

18.
高超超 《极地研究》2014,26(4):460-468
火山气溶胶通过反射太阳辐射导致地表降温。重建历史火山活动序列,是评估其气候效应,从而区分人为和自然因子的相对贡献、预测温室气体对全球气候影响的基础之一。此综述简要介绍了火山活动影响气候的作用机理。通过对冰芯研究历史及利用两极冰芯记录重建火山活动序列的研究进展的回顾,阐述了极地冰芯在火山活动序列重建中的作用与优势,以及近期研究在冰芯钻取和数据分析手段、极地火山信号向大气辐射强迫转化方法等方面所取得的成果;同时指出在火山沉积物的空间分布、冰芯定年及高低纬度火山喷发源的确定方面仍存在较大不确定性。因此应结合火山碎屑研究及非质量同位素分馏分析等方法,并钻取更多时间尺度更长且现有冰芯钻取较少的地区的深冰芯,以更精确地重建火山活动序列,从而提高对历史火山活动及其对气候影响的认识。  相似文献   

19.
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.  相似文献   

20.
Several years of broad-band teleseismic data from the GRSN stations have been analysed for crustal structure using P -to- S converted waves at the crustal discontinuities. An inversion technique was developed which applies the Thomson-Haskell formalism for plane waves without slowness integration. The main phases observed are Moho conversions, their multiples in the crust, and conversions at the base of the sediments. The crustal thickness derived from these data is in good agreement with results from other studies. For the Gräfenberg stations, we have made a more detailed comparison of our model with a previously published model obtained from refraction seismic experiments. The refraction seismic model contains boundaries with strong velocity contrasts and a significant low-velocity zone, resulting in teleseismic waveforms that are too complicated as compared to the observed simple waveforms. The comparison suggests that a significant low-velocity zone is not required and that internal crustal boundaries are rather smooth.  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号