Theory of the piezomagnetic effect in titanomagnetite-bearing rocks

Summary The stress sensitivitiesS _x andS _R of susceptibility and remanence for titanomagnetite-bearing rocks are calculated in terms of magnetostriction constants ₁₀₀ and ₁₁₁ and anisotropy constantsK ₁,K ₂ of the magnetic minerals.S _x andS _R are represented by quite different algebraic expressions but happen to have comparable numerical values over the whole range of titanomagnetite compositions. Both increase strongly with titanium content. This leads to more optimistic calculations of tectonomagnetic effects than with the previously assumed stress sensitivity for pure magnetite.     

(ϖμ/ϖT) P of the lower mantle

We estimate (/T) _P of the lower mantle at seismic frequencies using two distinct approaches by combining ambient laboratory measurements on lower mantle minerals with seismic data. In the first approach, an upper bound is estimated for |(/T) _P | by comparing the shear modulus () profile of PREM with laboratory room-temperature data of extrapolated to high pressures. The second approach employs a seismic tomography constraint ( lnV _S / lnV _P ) _P =1.8–2, which directly relates (/T) _P with (K _S /T) _P . An average (K _S /T) _P can be obtained by comparing the well-established room-temperature compression data for lower mantle minerals with theK _S profile of PREM along several possible adiabats. Both (K _S /T) and (/T) depend on silicon content [or (Mg+Fe)/Sil of the model. For various compositions, the two approaches predict rather distinct (/T) _P vs. (K _S /T) _P curves, which intersect at a composition similar to pyrolite with (/T) _P =–0.02 to –0.035 and (K _S /T) _P =–0.015 to –0.020 GPa/K. The pure perovskite model, on the other hand, yields grossly inconsistent results using the two approaches. We conclude that both vertical and lateral variations in seismic velocities are consistent with variation due to pressure, temperature, and phase transformations of a uniform composition. Additional physical properties of a pyrolite lower mantle are further predicted. Lateral temperature variations are predicted to be about 100–250 K, and the ratio of ( lnp/ lnV _S ) _P around 0.13 and 0.26. All of these parameters increase slightly with depth if the ratio of ( lnV _S / lnV _P ) _P remains constant throughout the lower mantle. These predicted values are in excellent agreement with geodynamic analyses, in which the ratios ( ln / lnV _S ) _P and ( / lnV _S ) _P are free parameters arbitrarily adjusted to fit the tomography and geoid data.     

A relationship between phonon conductivity and seismic parameter Φ for silicate minerals

Summary The relationship between the phonon conductivity at room temperature (K _N ) and the seismic parameter () for silicate minerals is suggested. The considerations are based on the Debye model of thermal energy transport phenomena in solids and on the seismic equation of state for silicates and oxides given byAnderson (1967). The semiempirical relationship is the formK _N = 0.43^0.82 where is in km²/s² andK _N in mcal/cm s K, and the empirical relationship isK _N =(0.528±0.006) –(8.18±2.11). The laboratory data on thermal and elastic properties for several silicates were taken fromHorai andSimmons (1970).     

An extended theory of equi-risk line for a general release rule

The expression of equi-risk line derived by the authors represents the relationship between discharge capacityy ₀ ^u and storage capacityz ₀ ^u to keep flood frequency under a certain risk level represented by the return periodT, i.e.,z ₀/z ₀ ^u ={(y ₀ ^u –y ₀ ^u )/y ₀ ^u } ^S , wherey ₀ ^u andz ₀ ^u areT-year probability peak discharge and total volume of a hydrograph. The shape parametersS is evaluated in this paper for various release rules of the storage facilities and correlations of durations and peaks of hydrographs. The expression forS is: , whereS ₀ andS are the values ofS forp=0 and , andp is the exponent of a general storage-release relation,q=az' ^p, wherea is the storage constant, andz' andq are the volume of stored water and the corresponding release. The values ofS ₀ andS are expressed in terms of the correlation coefficient of durations and peaks of inflow hydrographs.     

The elastic properties of single-crystal fayalite as determined by dynamical measurement techniques

We present new elasticity measurements on single-crystal fayalite and combine our results with other data from resonance, pulse superposition interferometry, and Brillouin scattering to provide a set of recommended values for the adiabatic elastic moduliC _ij and their temperature variations. We use a resonance method (RPR) with specimens that were previously investigated by pulse superposition experiments. The nineC _ij of fayalite are determined from three new sets of measurements. One set of our newC _ij data is over the range 300–500 K. We believe that the relatively large discrep ncies found in someC _ij are due in large part to specimen inhomogeneities (chemical and microstructural) coupled with differences in the way various techniques sample, rather than only systematic errors associated with experimental procedures or in the preparations of the specimens.Our recommendeaC _ij's (GPa) and (C _ij/T) _p (GPa/K) are: The resulting values for the isotropic bulk and shear moduli,K _s and , and their temperature derivatives are:K _s=134(4) GPa; =50.7(0.3) GPa; (K _s/T) _p =–0.024(0.005) GPa/K; and (/T) _p =–0.013(0.001) GPa/K. An important conclusion is thatK _s increases as the Fe/(Fe+Mg) ratio in olivine is increased.     

The effect of the number of charges carried by the argon ion on the darkening of ilford Q2 emulsion as ion detectors in mass-spectrometer

Summary The darkening (S) of Illford Q₂ photographic plates as ion detectors in mass spectrometer has been investigated. The dependence of the darkening (S) on the ion density (n=ions/mm²) i.e.S=S(n)E for constant energy (E)=z U ranging from 4U20 Kv of the impinging⁴⁰A⁺¹-,⁴⁰A⁺²- and⁴⁰A⁺³-ions whenS does not exceed the value 0.15 and the second relationn=n(z U) _S for darkening 0.05S0.15 constructed from the above relationS=S(n) _E has been determined. The darkening was found to increase with increasing ion-density which inturn decreases with the ionenergy. For⁴⁰A⁺¹-,⁴⁰A⁺²-, and⁴⁰A⁺³-ion of equal energy and ion-density the darkening effect was independent of the number of the charges carried by the argon ion.     

About digital alternative to theKp-index

For more than 40 yr the so-calledK-indices have been used to characterize irregular geomagnetic activity. Now that new methods of recording and digital computers are at our disposal it may be questioned, whether some alternative procedure for characteerizing of irregular geomagnetic activity has to be preferred.A proposal for a new procedure is presented. The indices such obtained are namedaK (meaning alternativeK), and some conformity with the oldK-indices is aimed at. It is not attempted to distinguish between wave radiation and corpuscular radiation from the Sun: consequently, the so-calledS _R variation is not removed.Elaboration of the procedure is realized by inserting a few extra statements in an extensive algorithm for general processing of digitally recorded data. A gradual modification of these statements, while working with the computer language PL/1, has lead to the final formulas.An appendix is added at the end of the paper. It is a condensed version of the procedure together with apossibly—clarifying example.Now retired; private address: Stigaardsvej 3, DK-2900 Hellerup     

Flow and transport through two-dimensional isotropic media of binary conductivity distribution. Part 2: NUMERICAL simulations and comparison with theoretical results

In the present part the results of numerical simulations of flow and transport in media made up from circular inclusions of conductivity K that are submerged in a matrix of conductivity K ₀, subjected to uniform mean velocity, are presented. This is achieved for a few values of =K/K ₀ (0.01, 0.1 and 10) and of the volume fraction n (0.05, 0.1 and 0.2). The numerical simulations (NS) are compared with the analytical approximate models presented in Part 1: the composite elements (CEA), the effective medium (EMA), the dilute system (DSA) and the first-order in the logconductivity variance (FOA). The comparison is made for the longitudinal velocity variance and for the longitudinal macrodispersivity. This is carried out for n<0.2, for which the theoretical and simulation models represent the same structure of random and independent inclusions distribution. The main result is that transport is quite accurately modeled by the EMA and CEA for low , for which _L is large, whereas in the case of =10, the EMA matches the NS for n<0.1. The first-order approximation is quite far apart from the NS for the values of examined . This material is based upon work supported by the National Science Foundation under Grant No. 0218914. Authors also wish to thank the Center of Computational Research, University at Buffalo for assistance in running numerical simulations.     

Case 25 application of the concentration parameter of seismoactive faults to Southern California

In this paper we evaluate the present state of the seismic regime in Southern California using the concentration parameter of seismogenic faults (K _sf ,Sobolev andZavyalov, 1981). The purpose of this work is to identify potential sites for large earthquakes during the next five or ten years. The data for this study derived from the California Institute of Technology's catalog of southern California earthquakes, and spanned the period between 1932 to June 1982. We examined events as small asM _L 1.8 but used a magnitude cutoff atM _L =3.3 for a detailed analysis. The size of the target earthquakes (M _M ) was chosen as 5.3 and 5.8.The algorithm for calculatingK _sf used here was improved over the algorithm described bySobolev andZavyalov (1981) in that it considered the seismic history of each elementary seismoactive volume. The dimensions of the elementary seismoactive volumes were 50 km×50 km and 20 km deep. We found that the mean value ofK _sf within 6 months prior to the target events was 6.1±2.0 for target events withM _L 5.3 and 5.41.8 for targets withM _L 5.8. Seventy-three percent of the targets withM _L 5.8 occurred in areas whereK _sf was less than 6.1. The variance of the time between the appearance of areas with lowK _sf values and the following main shocks was quite large (from a few months to ten years) so this parameter cannot be used here for accurate predictions of occurrence time.Regions where the value ofK _sf was below 6.1 at the end of our data set (June, 1982) are proposed as the sites of target earthquakes during the next five to ten years. The most dangerous area is the area east of San Bernardino whereK _sf values are presently between 2.9 and 3.7 and where there has been no earthquake withM _L 5.3 since 1948.     

Stress constraints and hydrofracturing stress data for the continental crust

Faulting and seismicity in the upper continental crust require considerable differential stresses. Application of experimentally developed friction, fracture and flow laws shows that high differential stresses can only exist in the uppermost crust. Direct hydraulic fracturing measurements in deep boreholes seem to support this rock mechanics conclusion. The experimental data base presently consists of approximately 500 hydrofrac tests conducted in about 100 boreholes at about 30 different geographical locations. To illustrate the variation of measured stresses with depth, the data are expressed as dimensionless horizontal stresses in the formS _H,h/S_v=(/z)+, whereS _v=gz Extrapolation of the experimental data to greater depth shows that the minor horizontal stress approaches the valueS _h/S_v=0.5 which limits friction on wet faults, and that the major horizontal stress approaches a value close toS _H/S_v=1 at rather shallow depth (5 to 10 km.). This limits faulting and seismicity in most of the upper crust to either strike-slip or normal faults. The lower boundary for seismicity is mainly dependent on tectonic strain accumulation and rock creep at the environmental conditions at depth.     

On estimating frequency-magnitude relations from heterogeneous catalogs

The paper studies the effect of magnitude errors on heterogeneous catalogs, by applying the apparent magnitude theory (seeTinti andMulargia, 1985a), which proves to be the most natural and rigorous approach to the problem. Heterogeneities in seismic catalogs are due to a number of various sources and affect both instrumental as well as noninstrumental earthquake compilations.The most frequent basis of heterogeneity is certainly that the recent instrumental records are to be combined with the historic and prehistoric event listings to secure a time coverage, considerably longer than the recurrence time of the major earthquakes. Therefore the case which attracts the greatest attention in the present analysis is that of a catalog consisting of a subset of higher quality data, generallyS ₁, spanning the interval T ₁ (the instrumental catalog), and of a second subset of more uncertain magnitude determination, generallyS ₂, covering a vastly longer interval T ₂ (the historic and/or the geologic catalog). The magnitude threshold of the subcatalogS ₁ is supposedly smaller than that ofS ₂, which, as we will see, is one of the major causes of discrepancy between the apparent magnitude and the true magnitude distributions. We will further suppose that true magnitude occurrences conform to theGutenberg-Richter (GR) law, because the assumption simplified the analysis without reducing the relevancy of our findings.The main results are: 1) the apparent occurrence rate exceeds the true occurrence rate from a certain magnitude onward, saym _GR; 2) the apparent occurrence rate shows two distinct GR regimes separated by an intermediate transition region. The offset between the two regimes is the essential outcome ofS ₁ being heterogeneous with respect toS ₂. The most important consequences of this study are that: 1) it provides a basis to infer the parameters of the true magnitude distribution, by correcting the bias deriving from heterogeneous magnitude errors; 2) it demonstrates that the double GR decay, that several authors have taken as the incontestable proof of the failure of the GR law and of the experimental evidence of the characteristic earthquake theory, is instead perfectly consistent with a GR-type seismicity.     

Solar activity dependence of ionospheric electron content and slab thickness using different solar indices

Ionospheric electron content (IEC) and slab thickness () data for the period 1977 to 1980 from Lunping (23.03°N; 121.90°E subionospheric) have been examined for their solar activity dependence. Local noontime monthly means as well as values for the 5 QQ days in a month have been examined separately with different solar indices, namely: solar EUV flux (170–190 Å),S _{10.7 cm} flux and sun spot number (SSN) on a seasonal basis. Both IEC and parameters exhibit better correlation with solar EUV andS _{10.7 cm} fluxes than with SSN for all seasons. IEC increases linearly with both EUV andS _{10.7 cm} flux whereas with SSN it shows a distinct nonlinear relationship during all seasons in both monthly mean and 5 QQ days' values. This study indicates that for correlating and predicting the variations (especially the medium term) in the ionospheric parameters, both EUV andS _{10.7 cm} fluxes have an advantage over SSN.     

Towards realistic flow modelling. Creation and evaluation of two-dimensional simulated porous media: An image analysis approach

This work is part of an attempt to quantify the relationship between the permeability tensor (K) and the micro-structure of natural porous media. A brief account is first provided of popular theories used to relate the micro-structure toK. Reasons for the lack of predictive power and restricted generality of current models are discussed. An alternative is an empirically based implicit model whereinK is expressed as a consequence of a few pore-types arising from the dynamics of depositional processes. The analytical form of that implicit model arises from evidence of universal association between pore-type and throat size in sandstones and carbonates. An explicit model, relying on the local change of scale technique is then addressed. That explicit model allows, from knowledge of the three-dimensional micro-geometry to calculateK explicitly without having recourse to any constitutive assumptions. The predictive and general character of the explicit model is underlined. The relevance of the change of scale technique is recalled to be contingent on the availability of rock-like three-dimensional synthetic media. A random stationary ergodic process is developed, that allows us to generate three-dimensional synthetic media from a two-dimensional autocorrelation functionr( _x , _y ) and associated probability density function measured on a single binary image. The focus of this work is to ensure the rock-like character of those synthetic media. This is done first through a direct approach:n two-dimensional synthetic media, derived from single set ( ,r( _x , _y )) yieldn permeability tensorsK _i-1,n ⁱ (calculated by the local change of scale) of the same order. This is a necessary condition to ensure thatr( _x , _y ) and carry all structural information relevant toK. The limits of this direct approach, in terms of required Central Process Unit time and Memory is underlined, raising the need for an alternative. This is done by comparing the pore-type content of a sandstone sample andn synthetic media derived fromr( _x , _y ) and measured on that sandstone-sample. Achievement of a good match ensures that the synthetic media comprise the fundamental structural level of all natural sandstones, that is a domainal structure of well-packed clusters of grains bounded by loose-packed pores.     

A numerical study of the heat transfer through a fluid layer with recirculating flow between concentric and eccentric spheres

A numerical study has been made of the heat transfer through a fluid layer with recirculating flow. The outer fluid surface was assumed to be spherical, while the inner surface consisted of a sphere concentrically or eccentrically located with respect to the outer spherical surface. The recirculating flow was assumed to be driven by a gas flow creating stress on the fluid's outer surface so that creeping (low Reynolds number) flow developed in its interior. The present study solves the Stokes equation of motion and the convective diffusion equation in bispherical coordinates and presents the streamline and isotherm patterns.Nomenclature a _i inner sphere radius - a _d outer sphere radius - A ₁ defined by equation (5) - A ₂ defined by equation (6) - B ₁ defined by equation (7) - B ₂ defined by equation (8) - c dimensional factor for bispherical coordinates - C constant in equation (4) - d narrowest distance between the two eccentric spheres - E ² operator defined by equation (1) in spherical coordinates and by equation (21) in bispherical coordinates - G modified vorticity, defined in equation (22) - G ^* non-dimensional modified vorticity, defined in equation (28) - h metric coefficient of bispherical coordinate system, defined in equation (18) - k _w thermal conductivity of water - K ₁ defined by equation (9) - K ₂ defined by equation (10) - N _Re Reynolds number=2a _dU/gn - N _Pe,h Peclet number=2a _dU/ - n integer counter - q heat flux - r radius - r ^* non-dimensional radius=r/a _d - S surface area - t time - t ^* non-dimensional time=t/a _d ² - T temperature - T _o temperature at inner sphere surface - T _a temperature at outer sphere surface - T ^* non-dimensional temperature;=(T–T _o)/(T_a–T_o) - u velocity - u _r radial velocity in spherical coordinates - u angular velocity in spherical coordinates - u radial velocity in bispherical coordinates - u angular velocity in bispherical coordinates - U free stream velocity - u _r ^* =u _r/U - u ^* =u /U - u ^* =u /U - u ^* =u /U Greek symbols a ₁ small displacement - vorticity, defined in equation (17) - ^* non-dimensional vorticity, defined in equation (27) - radial bispherical coordinates - _o bispherical coordinate of inner sphere - _a bispherical coordinate of outer sphere - angular coordinate in spherical coordinates - thermal diffusivity - _w thermal diffusivity of water - kinematic viscosity - angular bispherical coordinate - spherical coordinate - streamfunction - non-dimensional streamfunction for spherical coordinates, = /(U a _d ² ) - ^* non-dimensional streamfunction for bispherical coordinates, defined in equation (26)     

北京市北运河沉积物对氮、磷的吸附/解吸动力学特征

选取北运河土沟、榆林庄及和合站3个典型河道断面,进行河道0~60 cm深度沉积物对氨氮、磷酸盐的吸附/解吸批平衡静态试验研究,采用多种动力学模型分析北运河沉积物对氮、磷的吸附/解吸动力学特性,提出适用于北运河沉积物吸附/解吸特性的动力学模型,并对模型参数的主要影响因素进行探讨.结果表明:1)北运河典型断面各深度的沉积物对氨氮和磷酸盐的吸附/解吸过程总体呈现3个阶段:快反应阶段—慢反应阶段—平衡阶段;在0~0.5 h的快反应阶段可完成吸附或解吸总量的60%,且对氮、磷的吸附速率大于解吸速率.2)北运河各断面对氨氮与磷酸盐的平衡吸附量表现为S(榆林庄)S(土沟)S(和合站);对氨氮的平衡解吸量表现为S(土沟)S(榆林庄)S(和合站),对磷酸盐的平衡解吸量表现为S(和合站)S(土沟)S(榆林庄),沉积物对氨氮和磷酸盐的吸附以化学吸附为主,平衡吸附与解吸量随断面深度的增加而减小,0~20 cm表层沉积物对氮、磷的吸附能力较强.3)Lagergren二级动力学模型对各深度沉积物的吸附/解吸过程拟合最优,模型参数公式为k2=S-0.369max+0.163;qe=0.022 Smax+18.077 Kf+41.947.通过模拟得出在400 mg/L氮、磷浓度下吸附于沉积物中的氮、磷污染物会随着解吸过程释放52%~80%的氨氮和6%~42%的磷酸盐,可能不仅二次污染上覆水体,还随河水下渗从而对地下水质形成潜在污染风险.     

Multiple scattering of seismic waves in fractured media: Velocity and effective attenuation of the coherent components ofP waves andS waves

We show that the multiple scattering by small fractures of seismic waves with wavelengths long compared to the fracture size and fracture spacing is indistinguishable from multiple-scattering effects produced by regular porosity, except for an orientation factor due to fracture alignment. The fractures reduce theP-wave andS-wave velocities and produce an effective attenuation of the coherent component of the seismic waves. The attenuation corresponds to 1000/Q of about unity for a Gaussian spectrum of fractures, and it varies with frequencyf asf ³. For a Kolmogorov spectrum of fractures of spectral index the attenuation is an order of magnitude or so larger and varies with frequency asf ^3-v The precise degree of attenuation depends upon the matrix properties, the fracture porosity, the degree of fracture anisotropy, the type of fluid filling the fractures, and the incidence angle of the wave.For fracture porosities less than about 15% theP-wave andS-wave velocities are decreased by the order of 5–10% with a lesser dependence on the type of fluid filling the fractures (gas, oil, or brine) and with a dependence on both the degree of anisotropy and the incident angle made by the wave. The tendency of fractures to occur perpendicularly to bedding suggests that the best way to measure seismically fractured rock behavior in situ is by using the travel-time delay and reflection amplitude. As both the offset and the azimuth of receivers vary from a shot, the travel-time delay and reflection amplitude should both show an elliptical pattern of behavior—the travel-time delay in response to the varying seismic speed, and the reflection amplitude in response to angular variations in the multiple scattering. Observations of attenuation at several frequencies should permit (a) determination of the spectrum of fractures (Gaussian versus Kolmogorovian) and (b) determination of the contribution of viscous damping to the effective attenuation.     

Aspetto elettrogeosmotico bidimensionale variabile col tempo

Riassunto Fino ad oggi sono apparse solo teorie elettrogeosmotiche transitorie unidimensionali. Quì risolviamo un problema del genere bidimensionale, con elettrodi verticali cilindrici. Le curve rappresentative delle portate d'acqua catodicaq _K a cui perveniamo, denotano diminuzioni abbastanza rapide fino ai tempi dell' ordine dei «tempi elettrodici », (t=), portate che poi diminuiscono in modo molto lento. Si riscontrano quì più estesi «pianerottoli»q _K che non quelli delle monodimensionalità (t=/9). Inoltre, col crescere dei rapporti distanze eteropolari-raggi catodici, o col diminuire dei raggi catodici, tali portate (che si compongono dei contributi dei singoli anodi) s'incrementano nel modo diagrammato.

---

Summary Till now only mono-dimensional transient electro-geoosmotical theories have been published. In this paper a solution is given for a problem of bi-dimensional type, with vertical cylindrical electrodes. The plots which represent the amounts of cathodic waterq _K show quite rapid lowerings as far as to reach time durations of the order of «electrodic times » (t=) volumes which, after, diminish very slowly. In this case more extendedq _K «lobbies» are found than those in mono-dimensional case (t=/9). Furthermore increasing the ratios heteropolar spacings—cathodic radii, or diminishing the cathodic radii, water discharges (which result from the contribution of the single anods) increase as plotted in the diagrams.

     

Flow directions in ash-flow tuffs: a comparison of geological and magnetic susceptibility measurements, Tshirege member (upper Bandelier Tuff), Valles caldera, New Mexico, USA

This study concludes that the elongation axis (K ₁) of the ellipsoid of anisotropic magnetic susceptibility (AMS) is a suitable proxy for flow axis in ashflow tuffs. 153 oriented samples (176 specimens) were studied from 18 sites in the 1.1 Ma Tshirege member of the Bandelier Tuff. These sites are distributed around the Valles caldera at distances of 5–25 km outside of the rim.K ₁ axes correlate well with postulated radial flow axes at 13 sites.K ₁ also agrees with measured geological flow indicators, mainly imbricated larger clasts, at 7 sites. At 2 of the 5 sites where significant disagreement is seen between theoretical radial flow directions and measuredK ₁ axes, theK ₁ axes correspond well with geological flow indicators, indicating that the divergence of flow from the predicted radial flow pattern is real. Two major topographic buttresses are suggested as the cause of flow divergence for the Tshirege ash flows: the San Pedro buttress northwest of the caldera, and the San Miguel buttress in the southeast. In situK ₁ axes plunge about 7° toward the source at two-thirds of the sites; therefore the plunge ofK ₁ is a plausible in situ indicator for thedirection of flow. Multiple flow zones in sections of several meters thickness indicate changes of flow direction that are both rapid and large during ash-flow emplacement. These observations raisre the question of how best to represent mean flow directions in ash-flow sheets: by eigenvector methods, by vector-sum methods, or by modes. A method for measuring imbrication of larger clasts using apparent dips in vertical joints is outlined. Imbrication, determined in this way at one-third of the sites, dips toward the source, i.e., up-flow. The minimum (K ₃) axis of the AMS ellipsoid correlates with the flow foliation rather than with the larger clast imbrication. The flow axes of ash flows correspond with theK ₁ axes, not with the declination ofK ₃ axes as suggested by some authors. Initial dip of the sampled ash flows is not large and does not affect the paleomagnetic remanence direction, which is reversed with a mean ofD=173.5°,I=-38.4°, ₉₅=3.4°N=18. This mean is not different at the 95% confidence level from that of earlier workers. The mean pole, at 098.0°E, 74.8°N,A ₉₅=3.3°,N=18, is about 15° far-sided relative to the expected time-averaged geomagnetic pole, suggesting a history of emplacement too short to adequately average secular variation.     

The ratio between corner frequencies of source spectra of <Emphasis Type="Italic">P</Emphasis>- and <Emphasis Type="Italic">S</Emphasis>-waves—a new discriminant between earthquakes and quarry blasts

This study analyzes and compares the P- and S-wave displacement spectra from local earthquakes and explosions of similar magnitudes. We propose a new approach to discrimination between low-magnitude shallow earthquakes and explosions by using ratios of P- to S-wave corner frequencies as a criterion. We have explored 2430 digital records of the Israeli Seismic Network (ISN) from 456 local events (226 earthquakes, 230 quarry blasts, and a few underwater explosions) of magnitudes Md?=?1.4–3.4, which occurred at distances up to 250 km during 2001–2013 years. P-wave and S-wave displacement spectra were computed for all events following Brune’s source model of earthquakes (1970, 1971) and applying the distance correction coefficients (Shapira and Hofstetter, Teconophysics 217:217–226, 1993; Ataeva G, Shapira A, Hofstetter A, J Seismol 19:389-401, 2015), The corner frequencies and moment magnitudes were determined using multiple stations for each event, and then the comparative analysis was performed.The analysis showed that both P-wave and especially S-wave displacement spectra of quarry blasts demonstrate the corner frequencies lower than those obtained from earthquakes of similar magnitudes. A clear separation between earthquake and explosion populations was obtained for ratios of P- to S-wave corner frequency f ₀(P)/f ₀(S). The ratios were computed for each event with corner frequencies f ₀ of P- and S-wave, which were obtained from the measured f ₀ ^I at individual stations, then corrected for distance and finally averaged. We obtained empirically the average estimation of f ₀(P)/f ₀(S)?=?1.23 for all used earthquakes, and 1.86 for all explosions. We found that the difference in the ratios can be an effective discrimination parameter which does not depend on estimated moment magnitude M _w .The new multi-station Corner Frequency Discriminant (CFD) for earthquakes and explosions in Israel was developed based on ratios P- to S-wave corner frequencies f ₀(P)/f ₀(S), with the empirical threshold value of the ratio for Israel as 1.48.