Strong fast long-period waves in the Efpalio 2010 earthquake records: explanation in terms of leaking modes

The January 18, 2010, shallow earthquake in the Corinth Gulf, Greece (M _w  5.3) generated unusually strong long-period waves (periods 4–8 s) between the P and S wave arrival. These periods, being significantly longer than the source duration, indicated a structural effect. The waves were observed in epicentral distances 40–250 km and were significant on radial and vertical component. None of existing velocity models of the studied region provided explanation of the waves. By inverting complete waveforms, we obtained an 1-D crustal model explaining the observation. The most significant feature of the best-fitting model (as well as the whole suite of models almost equally well fitting the waveforms) is a strong velocity step at depth about 4 km. In the obtained velocity model, the fast long-period wave was modeled by modal summation and identified as a superposition of several leaking modes. In this sense, the wave is qualitatively similar to P long or Pnl waves, which however are usually reported in larger epicentral distances. The main innovation of this paper is emphasis to smaller epicentral distances. We studied properties of the wave using synthetic seismograms. The wave has a normal dispersion. Azimuthal and distance dependence of the wave partially explains its presence at 46 stations of 70 examined. Depth dependence shows that the studied earthquake was very efficient in the excitation of these waves just due to its shallow centroid depth (4.5 km).     

Seismic array evidence for a two-layer core transition zone

Measurements of traveltimes and traveltime gradients for PKP phases recorded at the Warramunga seismic array from distances between 113° and 176° provide evidence for a two-layer transition zone in the earth's core. The basic data consist of paper recordings played out at 40 mm/sec from analog magnetic tape with amplitude gain control and narrow bandpass filters. Measured values of $\text{d}\text{T}\text{dΔ}$ are perturbed by structure beneath the array and it is necessary to correct for this effect by an empirical approach. On the basis of the $\text{d}\text{T}\text{dΔ}$ values and traveltimes, two precursor branches to PKP_DF for Δ < 143° are identified consistently. These continue as branches intermediate to PKP_DF and PKP_AB for Δ > 145°. The corrected $\text{d}\text{T}\text{dΔ}$ values of all phases are smoothed and used to derive a velocity model for the core. The preferred model includes two velocity discontinuities of less than 0.1 km/sec each in the core transition zone to generate the two additional PKP branches. Together with the velocity jump at the outer-inner core boundary, these discontinuities define two layers each a few hundred kilometers thick. It is argued that recent proposals concerning inhomogeneities at the base of the mantle or within the transition zone as the source of PKP precursors do not agree with our observations nor those of other researchers. Inhomogeneities are not excluded but are considered a secondary effect for PKP phases, the primary effect being due to two discontinuous velocity increases in the transition zone.     

Automatic detection,location and quantification of earthquakes: Application to tsunami warning

We describe a fully automated seismic event detection and location system, providing for real-time estimates of the epicentral parameters of both local and distant earthquakes. The system uses 12 telemetered short-period stations, with a regional aperture of 350 km, as well as two 3-component broad-band stations. Detection and location of teleseismic events is achieved independently and concurrently on the short-period and long-period channels. The long-period data is then used to obtain an estimate of the seismic momentM ₀ of the earthquake through the mantle magnitudeM _m, as introduced byOkal andTalandier (1989). In turn, this estimate ofM ₀ is used to infer the expected tsunami amplitude at Papeete, within 15 minutes of the recording of Rayleigh waves. The performance of the method is discussed in terms of the accuracy of the epicentral parameters and seismic moment obtained in real time, as compared to the values later published by the reporting agencies. Our estimates are usually within 3 degrees of the reported epicenter, and the standard deviation on the seismic moment only 0.19 unit of magnitude for a population of 154 teleseismic events.     

Intrinsic Q and its frequency dependence

Various techniques for estimating $\text{t1}$ (travel time/quality factor Q) from short-period seismic-array records of body waves have been investigated. Spectral analysis in the frequency domain seems to be more appropriate for this purpose than time domain methods, because of the relative ease with which source and instrument effects can be removed. Of the techniques available, those based on maximum likelihood and homomorphic deconvolution give estimates of relative power versus frequency which best represent the power contained in a time-domain wavelet of short duration. The latter technique seemed to have better noise-eliminating properties than the former. Therefore, homomorphic deconvolution was used to obtain estimates of $\text{t1}$ values from P, PcP, ScP and S phases recorded at the Warramunga array in the Northern Territory of Australia. The source regions for the event studied were the Sunda, Mariana, New Hebrides, Kermadec and Tonga trench zones.The short-period $\text{t1}$ estimates obtained using the above method were much smaller than estimates from published free-oscillation Q models, indicating that the values of Q for compressional and shear waves are frequently-dependent. It was found that short-period $\text{t1}$ values and free-oscillation Q models could be made consistent with one another by assuming Q = Q₀(1+τω) where Q₀ and τ are constants. The results of this investigation suggest another approach to how the Q structure of the mantle can be investigated.     

Lithospheric structure of NW Iran from P and S receiver functions

We computed P and S receiver functions to investigate the lithospheric structure beneath the northwest Iran and compute the Vp/Vs ratio within the crust of this seismologically active area. Our results enabled us to map the lateral variations of the Moho as well as those of the lithosphere–asthenosphere boundary (LAB) beneath this region. We selected data from teleseismic events (M_b?>?5.5, epicentral distance between 30° and 95° for P receiver functions and M_b?>?5.7, epicentral distance between 60° and 85° for S receiver functions) recorded from 1995 to 2008 at 8 three-component short-period stations of Tabriz Telemetry Seismic Network. Our results obtained from P receiver functions indicate clear conversions at the Moho boundary. The Moho depth was firstly estimated from the delay time of the Moho converted phase relative to the direct P wave. Then we used the H-Vp/Vs stacking algorithm of Zhu and Kanamori to estimate the crustal thickness and Vp/Vs ratio underneath the stations with clear Moho multiples. We found an average Moho depth of 48 km, which varies between 38.5 and 53 km. The Moho boundary showed a significant deepening towards east and north. This may reveal a crustal thickening towards northeast possibly due to the collision between the Central Iran and South Caspian plates. The obtained average Vp/Vs ratio was estimated to be 1.76, which varies between 1.73 and 1.82. The crustal structure was also determined by modeling of P receiver functions. We obtained a three-layered model for the crust beneath this area. The thickness of the layers is estimated to be 6–11, 18–35, and 38–53 km, respectively. The average of the shear wave velocity was calculated to be 3.4 km/s in the crust and reaches 4.3 km/s below the Moho discontinuity. The crustal thickness values obtained from P receiver functions are in good agreement with those derived by S receiver functions. In addition, clear conversions with negative polarity were observed at ~8.7 s in S receiver functions, which could be related to the conversion at the LAB. This may show a relatively thin continental lithosphere of about 85 km implying that the lithosphere was influenced by various geodynamical reworking processes in the past.     

Creep response of the lunar crust in mare regions from an analysis of crater deformation

The settling trends of 318 lunar mare craters are compared with predictions of numerical finite-element models in order to determine the creep response of the upper lunar mare crust. No settling is evident in craters smaller than 5 km in diameter. Settling rates of larger craters increase as function of crater size in a manner suggesting a non-linear lunar creep response corresponding to the power law $\text{ε}\text{?}\text{= 8.3 · 10}{}^{\mathrm{?34}}\text{σ}{}^2$ where ε&#x0301; is the strain rate and σ is the differential stress. However, the observed nonlinearity is probably an apparent nonlinearity resulting from the temperature induced viscosity decrease with depth due to a lunar crustal temperature gradient of 3° C/km and a creep activation energy of 20 kcal/mole. It is concluded that creep in the lunar medium is essentially Newtonian, and that the effective viscosity of the upper lunar mare crust is (1.6 ± 0.3) · 10²⁵ poise.     

Modeling short-period crustal phases (P, Lg) for long-range refraction profiles

The short-period seismic phases known as $\text{P}$ and Lg are often recorded at distances of 200–1000 km on long-range refraction profiles and are usually the largest-amplitude features on record sections for this distance range. $\text{P}$ and Lg propagate as multiply reflected compressional and shear waves in a crustal waveguide whose principal boundaries are the Moho and the free surface. Equivalently, they can be interpreted as the interference pattern produced by a superposition of higher-mode P, SV and SH waves propagating in a leaky waveguide. For compressional waves, the waveguide efficiency is a strong function of frequency and depends on the presence or absence of low-velocity layers within a few kilometers of the surface, such as deep sedimentary sections commonly found in active tectonic areas. Such low-velocity surface layers create constructive interference effects for upcoming P waves incident at near grazing angles at the free surface and lead to efficient $\text{P}$ propagation. Several good examples of strong $\text{P}$ phases can be found on long-range refraction profiles for the tectonically active western United States; the 550 km profile eastward from SHOAL to Delta, UT is analyzed here. We have used a modified reflectivity-method computer program to model crustal phases for the SHOAL-Delta profile. The reflectivity technique accounts for all body and surface waves contributing to the short-period seismograms. It is found that the synthetic waveforms realistically model the observed $\text{P}$ characteristics. In this case, the decay of $\text{P}$ amplitudes with distance appears to be dominated by surface-reflection leakage from the waveguide rather than by anelastic attenuation due to Q of crustal rocks.     

Rupture process of the Miyagi-Oki,Japan, earthquake of June 12, 1978

The faulting mechanism and multiple rupture process of the M = 7.4 Miyagi-Oki earthquake are studied using surface and body wave data from local and worldwide stations. The main results are as follows. (1) P-wave first motion data and radiation patterns of long-period surface waves indicate a predominantly thrust mechanism with strike N10° E, dip 20°W, and slip angle 76°. The seismic moment is 3.1 × 10²⁷ dyne-cm. (2) Farfield SH waveforms and local seismograms suggest that the rupture occurred in two stages, being concordant with the two zones of aftershock activity revealed by the microearthquake network of Tohoku University. The upper and lower zones, located along the westward-dipping plate interface, are separated by a gap at a depth of 35 km and have dimensions of 37 × 34 and 24 × 34 km², respectively. Rupture initiated at the southern end of the upper aftershock zone and propagated at N20°W subparallel to the trench axis. About 11 s later, the second shock, which was located 30 km landward (westward) of the first, initiated at the upper corner of the lower aftershock zone and propagated down-dip N80°W. Using Haskell modelling for this rupture process, synthetic seismograms were computed for teleseismic SH waves and nearfield body waves. Other parameters determined are: seismic moment $\text{M}{}_0\text{= 1.7 × 10}{}^{27}\text{dyne-cm, slip dislocation}\text{u}\text{= 1.9}\text{m}\text{, Δσ = 95}\text{bar, rupture velocity}\text{ν = 3.2}\text{km s}{}^{\mathrm{?1}}\text{,}\text{rise time}\text{τ = 2}\text{s}$, for the first event; $\text{M}{}_0\text{= 1.4 × 10}{}^{27}\text{dyne-cm}\text{,}\text{u}\text{= 2.4}\text{m}\text{, Δσ = 145}\text{bar}$, for the second event; and time separation between the two shocks ΔT = 11 s. The above two-segment model does not explain well the sharp onsets of the body waves at near-source stations. An initial break of a small subsegment on the upper zone, which propagated down-dip, was hypothesized to explain the observed near-source seismograms. (3) The multiple rupture of the event and the absence of aftershocks between the two fault zones suggests that the frictional and/or sliding characteristics along the plate interface are not uniform. The rupture of the first event was arrested, presumably by a region of high fracture strength between the two zones. The fracture energy of the barrier was estimated to be 10¹⁰ erg cm^?2. (4) The possible occurrence of a large earthquake has been noted for the region adjacent to and seaward of the area that ruptured during the 1978 event. The 1978 event does not appear to reduce the likelihood of occurrence of this expected earthquake.     

High-pressure phase transformations and compressions of ilmenite and rutile,I. Experimental results

Natural ilmenite (Fe,Mg)TiO₃ has been found to transform to the perovskite structure and then to disproportionate into its component oxides, (Fe,Mg)O plus a cubic phase of TiO₂, at loading pressures of 140 and 250 kbar respectively, and at temperatures of 1,400 to 1,800°C. Samples were compressed in a diamond-anvil press and heated by irradiation with a YAG laser. The lattice parameters of the perovskite phase of (Fe,Mg)TiO₃ at room temperature and 1 bar are $\text{a}{}_0\text{= 4.471 ± 0.004, b}{}_0\text{= 5.753 ± 0.005,}\text{and}\text{c}{}_0\text{= 7.429 ± 0.006}\text{A}\text{?}$ with 4 molecules per cell. The zero-pressure volume change is 8.0% for the ilmenite-perovskite transition, 13.3% for the perovskite-mixed-oxides transition, and 20.2% for the ilmenite-mixed-oxides transition. The cubic phase of TiO₂ can be indexed on the basis of space group Fm3m with $\text{Z = 4}\text{and}\text{a}{}_0\text{= 4.455 ± 0.008}\text{A}\text{?}$ at room temperature and 1 bar, which corresponds to a decrease in zero-pressure volume of 29.2% for the rutile-cubic-phase transition. An isentropic bulk modulus at zero pressure of 5.75 ± 0.30 Mbar and a pressure derivative greater than 8 were calculated for the high-pressure cubic phase. The calculated bulk modulus for the mixture of (Fe,Mg)O and cubic TiO₂ is 2.48 ± 0.25 Mbar. All the phase transformations, the calculated lattice parameters, and the bulk moduli observed in this study are in good agreement with published shock-Hugoniot data for ilmenite and rutile.     

Elastic properties of garnet solid-solution series

The elastic constants of sixteen garnet specimens of wide variety in chemical composition are accurately determined by means of the rectangular parallelpiped resonance method. The dependence of the elastic properties on chemical composition is analyzed using the present data and those for seven garnets investigated by other authors. The property X_i of a garnet solid solution i is given by a linear addition law in terms of the mole fraction n_ij of component j; X_i = Σn_ijX_j where the X_j's are the properties of the end-members j (j = pyrope, almandine, spessartine, grossular and andradite). The X_j's are determined for density ρ, bulk modulus K, and shear moduli C_s = (C₁₁ ? C₁₂)/2 and C₄₄. No systematic deviation is observed from the linear addition law for the elastic moduli nor for other quantities such as the elastic wave velocities. The extrapolated elastic moduli (Mbar) of the end-members are:

     

11.

Source Parameters of the Deadly M<Subscript>w</Subscript> 7.6 Kashmir Earthquake of 8 October, 2005     

Prantik Mandal  R. K. Chadha  N. Kumar  I. P. Raju  C. Satyamurty 《Pure and Applied Geophysics》2007,164(10):1963-1983

During the last six years, National Geophysical Research Institute, Hyderabad has established a semi-permanent seismological network of 5–8 broadband seismographs and 10–20 accelerographs in the Kachchh seismic zone, Gujarat with a prime objective to monitor the continued aftershock activity of the 2001 M_w 7.7 Bhuj mainshock. The reliable and accurate broadband data for the 8 October M_w 7.6 2005 Kashmir earthquake and its aftershocks from this network as well as Hyderabad Geoscope station enabled us to estimate the group velocity dispersion characteristics and one-dimensional regional shear velocity structure of the Peninsular India. Firstly, we measure Rayleigh-and Love-wave group velocity dispersion curves in the period range of 8 to 35 sec and invert these curves to estimate the crustal and upper mantle structure below the western part of Peninsular India. Our best model suggests a two-layered crust: The upper crust is 13.8 km thick with a shear velocity (Vs) of 3.2 km/s; the corresponding values for the lower crust are 24.9 km and 3.7 km/sec. The shear velocity for the upper mantle is found to be 4.65 km/sec. Based on this structure, we perform a moment tensor (MT) inversion of the bandpass (0.05–0.02 Hz) filtered seismograms of the Kashmir earthquake. The best fit is obtained for a source located at a depth of 30 km, with a seismic moment, Mo, of 1.6 × 10²⁷ dyne-cm, and a focal mechanism with strike 19.5°, dip 42°, and rake 167°. The long-period magnitude (M_A ~ M_w) of this earthquake is estimated to be 7.31. An analysis of well-developed sPn and sSn regional crustal phases from the bandpassed (0.02–0.25 Hz) seismograms of this earthquake at four stations in Kachchh suggests a focal depth of 30.8 km.     

12.

Is the lower mantle homogeneous?     

A.L. Hales 《Journal of Geodynamics》1984,1(1):3-9

Comparison of some array dt/dΔ studies with the global travel times of Dziewonski and Anderson (1983) leads to the conclusion that a discontinuity in the P travel times between 80° and 85° is consistent with both sets of data. This discontinuity in dt/dΔ corresponds to an increase in velocity of about 0.1 km/sec between depths of 2400 and 2600 km. Models of the P velocity distribution which match the Dziewonski and Anderson travel times reasonable well have the shadow zone for short-period “diffracted” P beginning at about 110° arc distance.     

13.

Aftershock Activity and Frequency-dependent Low Coda Qc in the Epicentral Region of the 1999 Chamoli Earthquake of Mw 6.4     

P. Mandal  S. Padhy  B. K. Rastogi  H. V. S. Satyanarayana  M. Kousalya  R. Vijayraghavan  A. Srinivasan 《Pure and Applied Geophysics》2001,158(9-10):1719-1735

— On 28 March, 1999 (19:05:10.09, UT) a significant earthquake of M _w 6.4 occurred in the Garhwal Himalaya (30.555°N, 79.424°E). One hundred and ten well-recorded aftershocks show a WNW-ESE trending northeasterly dipping seismic zone extending from a depth of 2 to 20?km. As the main shock hypocenter occurred at the northern end of this seismic zone and aftershocks extended updip, it is inferred that the main-shock rupture nucleated on the detachment plane at a depth of 15?km and then propagated updip along a NE-dipping thrust plane. Further, the epicentral distribution of aftershocks defines a marked concentration near a zone where main central thrust (MCT) takes a significant turn towards the north, which might be acting as an asperity in response to the NNE compression due to the underthrusting of Himalayan orogenic process prevalent in the entire region. Presence of high seismicity including five earthquakes of magnitude exceeding 6 and twelve earthquakes of magnitude exceeding 5 in the 20th century is presumed to have caused a higher level of shallow crustal heterogeneity in the Garhwal Himalaya, a site lying in the central gap zone of the Himalayan frontal arc. Attenuation property of the medium around the epicentral area of the 1999 Chamoli earthquake, covering a circular area of 61,500?km² with a radius of 140?km, is studied by estimating the coda Q _c from 48 local earthquakes of magnitudes varying from 2.5–4.8. These earthquakes were recorded at nine 24-bit REFTEK digital stations; two of which were equipped with three-component CMG40T broadband seismometers and others with three-component L4-3D short-period seismometers. The estimated Q _o values at different stations suggest on average a low value of the order of (30?±?0.8), indicating an attenuating crust beneath the entire region. The frequency-dependent relation indicates a relatively low Q _c at lower frequencies (1–3?Hz) that can be attributed to the loss of energy due to scattering on heterogeneities and/or the presence of faults and cracks. The large Q _c at higher frequencies may be related to the propagation of backscattered body waves through deeper parts of the lithosphere where less heterogeneities are expected. An important observation is that the region north of MCT (more rigid highly metamorphosed crystalline rocks) is less attenuative in comparison to the region south of MCT (less rigid slightly metamorphosed rocks (sedimentary wedge)). The acceleration decays to 50% at 20?km distance and to 7% at 100?km. Hence, even 1g acceleration at the source may not cause significant damage beyond 100?km in this region.     

14.

Global electromagnetic induction in the Moon and planets     

Palmer Dyal  Curtis W Parkin 《Physics of the Earth and Planetary Interiors》1973,7(3):251-265

A summary of experiments and analyses concerning electromagnetic induction in the Moon and other extraterrestrial bodies is presented. Magnetic step-transient measurements made on the lunar dark side show the eddy current response to be the dominant induction mode of the Moon. Analysis of the poloidal field decay of the eddy currents has yielded a range of monotonic conductivity profiles for the lunar interior: the conductivity rises from 3·10^?4 mho/m at a depth of 170 km to 10^?2 mho/m at 1000 km depth. The static magnetization field induction has been measured and the whole-Moon relative magnetic permeability has been calculated to be $\text{μ}\text{μ}{}_0\text{= 1.01 ± 0.06}$. The remanent magnetic fields, measured at Apollo landing sites, range from 3 to 327 γ. Simultaneous magnetometer and solar wind spectrometer measurements show that the 38-γ remanent field at the Apollo 12 site is compressed to 54 γ by a solar wind pressure increase of 7·10^?8 dyn/cm². The solar wind confines the induced lunar poloidal field; the field is compressed to the surface on the lunar subsolar side and extends out into a cylindrical cavity on the lunar antisolar side. This solar wind confinement is modeled in the laboratory by a magnetic dipole enclosed in a superconducting lead cylinder; results show that the induced poloidal field geometry is modified in a manner similar to that measured on the Moon. Induction concepts developed for the Moon are extended to estimate the electromagnetic response of other bodies in the solar system.     

15.

Nemuro-Oki earthquake of June 17, 1973: A lithospheric rebound at the upper half of the interface     

Kunihiko Shimazaki 《Physics of the Earth and Planetary Interiors》1974,9(4):314-327

P-wave first motions, radiation patterns and amplitudes of long-period surface waves, relocated aftershock distributions, leveling and tsunami data indicate that the 1973 Nemuro-Oki earthquake is caused by a low-angle thrust-faulting, representing a rebound at the upper 50 km of the interface between the continental and oceanic lithospheres. Rebound, most likely aseismic, at depths below 50 km, is suggested to take place in the near future from a comparison of recent geologic crustal deformation with pre-seismic and co-seismic data. The estimated seismic moment is about $\text{1}\text{3}\text{–}\text{1}\text{4}$ of that for the neighboring great earthquakes. The macro-seismic data suggest that the 1973 earthquake is smaller than the 1894 Nemuro-Oki earthquake, the last great earthquake in this region.The 1973 earthquake had been predicted on the basis of a seismic gap. Although the prediction was successful as to the location and nature of the faulting and partly as to the occurrence time, it is smaller than the predicted one. A part of the seismic gap may still remain. The difference between the observed seismic slip (1.6 m) and that predicted on the basis of the pre-seismic crustal deformation (3.0 m) indicates either (1) the 1973 earthquake relieved only a part of the strain accumulated in the upper 50 km, or (2) a significant amount of aseismic slip took place on the seismic fault and completely relieved the accumulated strain in the focal region of the 1973 earthquake. If the former is the case, the remaining strain, not only in the focal region, but also in the remaining seismic gap adjoining it, may be relieved in a larger earthquake in the future.The source parameters obtained are as follows: fault plane, dip direction = N40°W, dip angle = 27°; seismic moment = 6.7 · 10²⁷ dyn cm; average slip dislocation, 1.6 m in N63°W direction; stress drop = 35 bars. In these calculations, the fault dimension and the rigidity are assumed to be 100 · 60 km² and 7.0 · 10¹¹ dyn/cm², respectively.     

16.

Source Parameters of the 8 October, 2005 M<Subscript>w</Subscript>7.6 Kashmir Earthquake     

Prantik Mandal  R. K. Chadha  N. Kumar  I. P. Raju  C. Satyamurty 《Pure and Applied Geophysics》2007,164(11):2235-2254

During the last six years, the National Geophysical Research Institute, Hyderabad has established a semi-permanent seismological network of 5 broadband seismographs and 10 accelerographs in the Kachchh seismic zone, Gujarat, with the prime objective to monitor the continued aftershock activity of the 2001 M_w7.7 Bhuj mainshock. The reliable and accurate broadband data for the M_w 7.6 (8 Oct., 2005) Kashmir earthquake and its aftershocks from this network, as well as from the Hyderabad Geoscope station, enabled us to estimate the group velocity dispersion characteristics and the one-dimensional regional shear-velocity structure of peninsular India. Firstly, we measure Rayleigh- and Love-wave group velocity dispersion curves in the range of 8 to 35 sec and invert these curves to estimate the crustal and upper mantle structure below the western part of peninsular India. Our best model suggests a two-layered crust: The upper crust is 13.8-km thick with a shear velocity (Vs) of 3.2 km/s; the corresponding values for the lower crust are 24.9 km and 3.7 km/sec. The shear velocity for the upper mantle is found to be 4.65 km/sec. Based on this structure, we perform a moment tensor (MT) inversion of the bandpass (0.05–0.02 Hz) filtered seismograms of the Kashmir earthquake. The best fit is obtained for a source located at a depth of 30 km, with a seismic moment, Mo, of 1.6 × 10²⁷ dyne-cm, and a focal mechanism with strike 19.5°, dip 42°, and rake 167°. The long-period magnitude (M_A ~ M_w) of this earthquake is estimated to be 7.31. An analysis of well-developed sPn and sSn regional crustal phases from the bandpassed (0.02–0.25 Hz) seismograms of this earthquake at four stations in Kachchh suggests a focal depth of 30.8 km.     

17.

Synthesis of a new high-pressure ohase of tin dioxide and some geophysical implications     

Lin-Gun Liu 《Physics of the Earth and Planetary Interiors》1974,9(4):338-343

Tin dioxide (SnO₂) in the rutile structure as starting material has been found to transform to the orthorhombic α-PbO₂ structure (S.G. Pbcn) at about 155 kbar and 1000–1400°C when compressed in a diamond-anvil cell and heated by irradiation with a YAG laser. The lattice parameters at room temperature and 1 bar are $\text{a}{}_{\mathrm{<mtext>o</mtext>}}\text{= 4.719 ± 0.002, b}{}_{\mathrm{<mtext>o</mtext>}}\text{= 5.714 ± 0.002,}\text{and}\text{c}{}_{\mathrm{<mtext>o</mtext>}}\text{= 5.228 ± 0.002}\text{A}\text{?}\text{with}\text{Z = 4}$ for the orthorhombic form of SnO₂, which is 1.5% more dense than the rutile form. Crystal-chemical arguments suggest that stishovite (SiO₂) may also transform to the α-PbO₂ structure at elevated pressure and temperature with an increase in zero-pressure density of about 2–3%. Mineral assemblages containing the orthorhombic SiO₂ are unstable relative to those containing the perovskite MgSiO₃ under lower-mantle conditions.     

18.

The Campania-Lucania (southern Italy) earthquake of 23 November 1980     

A. Deschamps  G.C.P. King 《Earth and Planetary Science Letters》1983,62(2):296-304

We present the main seismological results of our study of the Campania-Lucania earthquake of 23 November 1980. A complete set of far field and local data has been analysed. From long-period body waves data we determine the fault plane solution (φ₁ = 140°,δ₁ = 60°,φ₂ = 75°,δ₂ = 54°), a depth of 15 km and calculate a seismic moment of 6 × 10²⁵ dyne cm and a source duration of 6 s. From data of a local network deployed immediately after the event we determine aftershock locations: they are aligned in a direction NW-SE that fit extremely well with the focal solution determined above. We can choose as fault plane the plane striking 140° and dipping at 60° and the event is a normal event with a large component of left-lateral strike slip. The source area evaluated from this aftershock distribution 14 km × 40 km is quite suitable for an earthquake of a seismic moment of 6 × 10²⁵ dyne cm.     

19.

An experimental study of magnetite-titanomagnetite interdiffusion     

R. Freer  Z. Hauptman 《Physics of the Earth and Planetary Interiors》1978,16(3):223-231

Interdiffusion experiments were performed between Fe₃O₄ (single crystal) and Fe_2.8Ti_0.2O₄ (powder), under self-buffering conditions (temperature range 600–1034°C), and for various oxygen potentials at 1400°C. Profiles of Fe and Ti were obtained by electronprobe microanalysis, and the interdiffusion coefficient $\text{D}$ was calculated by the Boltzmann-Matano method. Low-temperature data at 3 mole% Ti could be described by $\text{D}\text{= (3.85}{}_{\mathrm{?1.11}}{}^{\mathrm{+1.68}}\text{) × 10}{}^{\mathrm{?3}}\text{exp}\text{(2.23 ± 0.04}\text{eV}\text{/kT)}\text{cm}{}^2\text{/}\text{s}$. An estimate is given for the time to interdiffuse 2μm at various temperatures, and the results compared with recent experiments.     

20.

Crustal structure of Northwest Zagros (Kermanshah) and Central Iran (Yazd and Isfahan) using teleseismic Ps converted phases     

Narges Afsari  Forogh Sodoudi  Fataneh Taghizadeh Farahmand  Mohammad Reza Ghassemi 《Journal of Seismology》2011,15(2):341-353

Receiver functions are widely employed to detect P-to-S converted waves and are especially useful to image seismic discontinuities in the crust. In this study we used the P receiver function technique to investigate the velocity structure of the crust beneath the Northwest Zagros and Central Iran and map out the lateral variation of the Moho boundary within this area. Our dataset includes teleseismic data (M _b ≥ 5.5, epicentral distance from 30° to 95°) recorded at 12 three-component short-period stations of Kermanshah, Isfahan and Yazd telemetry seismic networks. Our results obtained from P receiver functions indicate clear Ps conversions at the Moho boundary. The Moho depths were firstly estimated from the delay time of the Moho converted phase relative to the direct P wave beneath each network. Then, we used the P receiver function inversion to find the properties of the Moho discontinuity such as depth and velocity contrast. Our results obtained from PRF are in good agreement with those obtained from the P receiver function modeling. We found an average Moho depth of about 42 km beneath the Northwest Zagros increasing toward the Sanandaj-Sirjan Metamorphic Zone and reaches 51 km, where two crusts (Zagros and Central Iran) are assumed to be superposed. The Moho depth decreases toward the Urmieh-Dokhtar Cenozoic volcanic belt and reaches 43 km beneath this area. We found a relatively flat Moho beneath the Central Iran where, the average crustal thickness is about 42 km. Our P receiver function modeling revealed a shear wave velocity of 3.6 km/s in the crust of Northwest Zagros and Central Iran increasing to 4.5 km/s beneath the Moho boundary. The average shear wave velocity in the crust of UDMA as SSZ is 3.6 km/s, which reaches to 4.0 km/s while in SSZ increases to 4.3 km/s beneath the Moho.     

	Almandine	Pyrope	Spessartine	Grossular	Andradite
$\text{K}$	1.779 ± 0.008	1.730 ± 0.009	1.742 ± 0.009	1.691 ± 0.008	1.379 ± 0.017
$\text{C}{}_{\mathrm{<mtext>s</mtext>}}$	0.981 ± 0.004	0.925 ± 0.004	0.964 ± 0.004	1.106 ± 0.004	0.979 ± 0.007
$\text{C}{}_{\mathrm{<mtext>44</mtext>}}$	0.958 ± 0.005	0.919 ± 0.005	0.937 ± 0.005	1.017 ± 0.006	0.827 ± 0.010

Source Parameters of the Deadly M<Subscript>w</Subscript> 7.6 Kashmir Earthquake of 8 October, 2005

During the last six years, National Geophysical Research Institute, Hyderabad has established a semi-permanent seismological network of 5–8 broadband seismographs and 10–20 accelerographs in the Kachchh seismic zone, Gujarat with a prime objective to monitor the continued aftershock activity of the 2001 M_w 7.7 Bhuj mainshock. The reliable and accurate broadband data for the 8 October M_w 7.6 2005 Kashmir earthquake and its aftershocks from this network as well as Hyderabad Geoscope station enabled us to estimate the group velocity dispersion characteristics and one-dimensional regional shear velocity structure of the Peninsular India. Firstly, we measure Rayleigh-and Love-wave group velocity dispersion curves in the period range of 8 to 35 sec and invert these curves to estimate the crustal and upper mantle structure below the western part of Peninsular India. Our best model suggests a two-layered crust: The upper crust is 13.8 km thick with a shear velocity (Vs) of 3.2 km/s; the corresponding values for the lower crust are 24.9 km and 3.7 km/sec. The shear velocity for the upper mantle is found to be 4.65 km/sec. Based on this structure, we perform a moment tensor (MT) inversion of the bandpass (0.05–0.02 Hz) filtered seismograms of the Kashmir earthquake. The best fit is obtained for a source located at a depth of 30 km, with a seismic moment, Mo, of 1.6 × 10²⁷ dyne-cm, and a focal mechanism with strike 19.5°, dip 42°, and rake 167°. The long-period magnitude (M_A ~ M_w) of this earthquake is estimated to be 7.31. An analysis of well-developed sPn and sSn regional crustal phases from the bandpassed (0.02–0.25 Hz) seismograms of this earthquake at four stations in Kachchh suggests a focal depth of 30.8 km.     

Is the lower mantle homogeneous?

Comparison of some array dt/dΔ studies with the global travel times of Dziewonski and Anderson (1983) leads to the conclusion that a discontinuity in the P travel times between 80° and 85° is consistent with both sets of data. This discontinuity in dt/dΔ corresponds to an increase in velocity of about 0.1 km/sec between depths of 2400 and 2600 km. Models of the P velocity distribution which match the Dziewonski and Anderson travel times reasonable well have the shadow zone for short-period “diffracted” P beginning at about 110° arc distance.     

Aftershock Activity and Frequency-dependent Low Coda Qc in the Epicentral Region of the 1999 Chamoli Earthquake of Mw 6.4

— On 28 March, 1999 (19:05:10.09, UT) a significant earthquake of M _w 6.4 occurred in the Garhwal Himalaya (30.555°N, 79.424°E). One hundred and ten well-recorded aftershocks show a WNW-ESE trending northeasterly dipping seismic zone extending from a depth of 2 to 20?km. As the main shock hypocenter occurred at the northern end of this seismic zone and aftershocks extended updip, it is inferred that the main-shock rupture nucleated on the detachment plane at a depth of 15?km and then propagated updip along a NE-dipping thrust plane. Further, the epicentral distribution of aftershocks defines a marked concentration near a zone where main central thrust (MCT) takes a significant turn towards the north, which might be acting as an asperity in response to the NNE compression due to the underthrusting of Himalayan orogenic process prevalent in the entire region. Presence of high seismicity including five earthquakes of magnitude exceeding 6 and twelve earthquakes of magnitude exceeding 5 in the 20th century is presumed to have caused a higher level of shallow crustal heterogeneity in the Garhwal Himalaya, a site lying in the central gap zone of the Himalayan frontal arc. Attenuation property of the medium around the epicentral area of the 1999 Chamoli earthquake, covering a circular area of 61,500?km² with a radius of 140?km, is studied by estimating the coda Q _c from 48 local earthquakes of magnitudes varying from 2.5–4.8. These earthquakes were recorded at nine 24-bit REFTEK digital stations; two of which were equipped with three-component CMG40T broadband seismometers and others with three-component L4-3D short-period seismometers. The estimated Q _o values at different stations suggest on average a low value of the order of (30?±?0.8), indicating an attenuating crust beneath the entire region. The frequency-dependent relation indicates a relatively low Q _c at lower frequencies (1–3?Hz) that can be attributed to the loss of energy due to scattering on heterogeneities and/or the presence of faults and cracks. The large Q _c at higher frequencies may be related to the propagation of backscattered body waves through deeper parts of the lithosphere where less heterogeneities are expected. An important observation is that the region north of MCT (more rigid highly metamorphosed crystalline rocks) is less attenuative in comparison to the region south of MCT (less rigid slightly metamorphosed rocks (sedimentary wedge)). The acceleration decays to 50% at 20?km distance and to 7% at 100?km. Hence, even 1g acceleration at the source may not cause significant damage beyond 100?km in this region.     

Global electromagnetic induction in the Moon and planets

A summary of experiments and analyses concerning electromagnetic induction in the Moon and other extraterrestrial bodies is presented. Magnetic step-transient measurements made on the lunar dark side show the eddy current response to be the dominant induction mode of the Moon. Analysis of the poloidal field decay of the eddy currents has yielded a range of monotonic conductivity profiles for the lunar interior: the conductivity rises from 3·10^?4 mho/m at a depth of 170 km to 10^?2 mho/m at 1000 km depth. The static magnetization field induction has been measured and the whole-Moon relative magnetic permeability has been calculated to be $\text{μ}\text{μ}{}_0\text{= 1.01 ± 0.06}$. The remanent magnetic fields, measured at Apollo landing sites, range from 3 to 327 γ. Simultaneous magnetometer and solar wind spectrometer measurements show that the 38-γ remanent field at the Apollo 12 site is compressed to 54 γ by a solar wind pressure increase of 7·10^?8 dyn/cm². The solar wind confines the induced lunar poloidal field; the field is compressed to the surface on the lunar subsolar side and extends out into a cylindrical cavity on the lunar antisolar side. This solar wind confinement is modeled in the laboratory by a magnetic dipole enclosed in a superconducting lead cylinder; results show that the induced poloidal field geometry is modified in a manner similar to that measured on the Moon. Induction concepts developed for the Moon are extended to estimate the electromagnetic response of other bodies in the solar system.     

Nemuro-Oki earthquake of June 17, 1973: A lithospheric rebound at the upper half of the interface

P-wave first motions, radiation patterns and amplitudes of long-period surface waves, relocated aftershock distributions, leveling and tsunami data indicate that the 1973 Nemuro-Oki earthquake is caused by a low-angle thrust-faulting, representing a rebound at the upper 50 km of the interface between the continental and oceanic lithospheres. Rebound, most likely aseismic, at depths below 50 km, is suggested to take place in the near future from a comparison of recent geologic crustal deformation with pre-seismic and co-seismic data. The estimated seismic moment is about $\text{1}\text{3}\text{–}\text{1}\text{4}$ of that for the neighboring great earthquakes. The macro-seismic data suggest that the 1973 earthquake is smaller than the 1894 Nemuro-Oki earthquake, the last great earthquake in this region.The 1973 earthquake had been predicted on the basis of a seismic gap. Although the prediction was successful as to the location and nature of the faulting and partly as to the occurrence time, it is smaller than the predicted one. A part of the seismic gap may still remain. The difference between the observed seismic slip (1.6 m) and that predicted on the basis of the pre-seismic crustal deformation (3.0 m) indicates either (1) the 1973 earthquake relieved only a part of the strain accumulated in the upper 50 km, or (2) a significant amount of aseismic slip took place on the seismic fault and completely relieved the accumulated strain in the focal region of the 1973 earthquake. If the former is the case, the remaining strain, not only in the focal region, but also in the remaining seismic gap adjoining it, may be relieved in a larger earthquake in the future.The source parameters obtained are as follows: fault plane, dip direction = N40°W, dip angle = 27°; seismic moment = 6.7 · 10²⁷ dyn cm; average slip dislocation, 1.6 m in N63°W direction; stress drop = 35 bars. In these calculations, the fault dimension and the rigidity are assumed to be 100 · 60 km² and 7.0 · 10¹¹ dyn/cm², respectively.     

Source Parameters of the 8 October, 2005 M<Subscript>w</Subscript>7.6 Kashmir Earthquake

During the last six years, the National Geophysical Research Institute, Hyderabad has established a semi-permanent seismological network of 5 broadband seismographs and 10 accelerographs in the Kachchh seismic zone, Gujarat, with the prime objective to monitor the continued aftershock activity of the 2001 M_w7.7 Bhuj mainshock. The reliable and accurate broadband data for the M_w 7.6 (8 Oct., 2005) Kashmir earthquake and its aftershocks from this network, as well as from the Hyderabad Geoscope station, enabled us to estimate the group velocity dispersion characteristics and the one-dimensional regional shear-velocity structure of peninsular India. Firstly, we measure Rayleigh- and Love-wave group velocity dispersion curves in the range of 8 to 35 sec and invert these curves to estimate the crustal and upper mantle structure below the western part of peninsular India. Our best model suggests a two-layered crust: The upper crust is 13.8-km thick with a shear velocity (Vs) of 3.2 km/s; the corresponding values for the lower crust are 24.9 km and 3.7 km/sec. The shear velocity for the upper mantle is found to be 4.65 km/sec. Based on this structure, we perform a moment tensor (MT) inversion of the bandpass (0.05–0.02 Hz) filtered seismograms of the Kashmir earthquake. The best fit is obtained for a source located at a depth of 30 km, with a seismic moment, Mo, of 1.6 × 10²⁷ dyne-cm, and a focal mechanism with strike 19.5°, dip 42°, and rake 167°. The long-period magnitude (M_A ~ M_w) of this earthquake is estimated to be 7.31. An analysis of well-developed sPn and sSn regional crustal phases from the bandpassed (0.02–0.25 Hz) seismograms of this earthquake at four stations in Kachchh suggests a focal depth of 30.8 km.     

Synthesis of a new high-pressure ohase of tin dioxide and some geophysical implications

Tin dioxide (SnO₂) in the rutile structure as starting material has been found to transform to the orthorhombic α-PbO₂ structure (S.G. Pbcn) at about 155 kbar and 1000–1400°C when compressed in a diamond-anvil cell and heated by irradiation with a YAG laser. The lattice parameters at room temperature and 1 bar are $\text{a}{}_{\mathrm{<mtext>o</mtext>}}\text{= 4.719 ± 0.002, b}{}_{\mathrm{<mtext>o</mtext>}}\text{= 5.714 ± 0.002,}\text{and}\text{c}{}_{\mathrm{<mtext>o</mtext>}}\text{= 5.228 ± 0.002}\text{A}\text{?}\text{with}\text{Z = 4}$ for the orthorhombic form of SnO₂, which is 1.5% more dense than the rutile form. Crystal-chemical arguments suggest that stishovite (SiO₂) may also transform to the α-PbO₂ structure at elevated pressure and temperature with an increase in zero-pressure density of about 2–3%. Mineral assemblages containing the orthorhombic SiO₂ are unstable relative to those containing the perovskite MgSiO₃ under lower-mantle conditions.     

The Campania-Lucania (southern Italy) earthquake of 23 November 1980

We present the main seismological results of our study of the Campania-Lucania earthquake of 23 November 1980. A complete set of far field and local data has been analysed. From long-period body waves data we determine the fault plane solution (φ₁ = 140°,δ₁ = 60°,φ₂ = 75°,δ₂ = 54°), a depth of 15 km and calculate a seismic moment of 6 × 10²⁵ dyne cm and a source duration of 6 s. From data of a local network deployed immediately after the event we determine aftershock locations: they are aligned in a direction NW-SE that fit extremely well with the focal solution determined above. We can choose as fault plane the plane striking 140° and dipping at 60° and the event is a normal event with a large component of left-lateral strike slip. The source area evaluated from this aftershock distribution 14 km × 40 km is quite suitable for an earthquake of a seismic moment of 6 × 10²⁵ dyne cm.     

An experimental study of magnetite-titanomagnetite interdiffusion

Interdiffusion experiments were performed between Fe₃O₄ (single crystal) and Fe_2.8Ti_0.2O₄ (powder), under self-buffering conditions (temperature range 600–1034°C), and for various oxygen potentials at 1400°C. Profiles of Fe and Ti were obtained by electronprobe microanalysis, and the interdiffusion coefficient $\text{D}$ was calculated by the Boltzmann-Matano method. Low-temperature data at 3 mole% Ti could be described by $\text{D}\text{= (3.85}{}_{\mathrm{?1.11}}{}^{\mathrm{+1.68}}\text{) × 10}{}^{\mathrm{?3}}\text{exp}\text{(2.23 ± 0.04}\text{eV}\text{/kT)}\text{cm}{}^2\text{/}\text{s}$. An estimate is given for the time to interdiffuse 2μm at various temperatures, and the results compared with recent experiments.     

Crustal structure of Northwest Zagros (Kermanshah) and Central Iran (Yazd and Isfahan) using teleseismic Ps converted phases

Receiver functions are widely employed to detect P-to-S converted waves and are especially useful to image seismic discontinuities in the crust. In this study we used the P receiver function technique to investigate the velocity structure of the crust beneath the Northwest Zagros and Central Iran and map out the lateral variation of the Moho boundary within this area. Our dataset includes teleseismic data (M _b ≥ 5.5, epicentral distance from 30° to 95°) recorded at 12 three-component short-period stations of Kermanshah, Isfahan and Yazd telemetry seismic networks. Our results obtained from P receiver functions indicate clear Ps conversions at the Moho boundary. The Moho depths were firstly estimated from the delay time of the Moho converted phase relative to the direct P wave beneath each network. Then, we used the P receiver function inversion to find the properties of the Moho discontinuity such as depth and velocity contrast. Our results obtained from PRF are in good agreement with those obtained from the P receiver function modeling. We found an average Moho depth of about 42 km beneath the Northwest Zagros increasing toward the Sanandaj-Sirjan Metamorphic Zone and reaches 51 km, where two crusts (Zagros and Central Iran) are assumed to be superposed. The Moho depth decreases toward the Urmieh-Dokhtar Cenozoic volcanic belt and reaches 43 km beneath this area. We found a relatively flat Moho beneath the Central Iran where, the average crustal thickness is about 42 km. Our P receiver function modeling revealed a shear wave velocity of 3.6 km/s in the crust of Northwest Zagros and Central Iran increasing to 4.5 km/s beneath the Moho boundary. The average shear wave velocity in the crust of UDMA as SSZ is 3.6 km/s, which reaches to 4.0 km/s while in SSZ increases to 4.3 km/s beneath the Moho.