A scheme to set preferred magnitudes in the ISC Bulletin

One of the main purposes of the International Seismological Centre (ISC) is to collect, integrate and reprocess seismic bulletins provided by agencies around the world in order to produce the ISC Bulletin. This is regarded as the most comprehensive bulletin of the Earth’s seismicity, and its production is based on a unique cooperation in the seismological community that allows the ISC to complement the work of seismological agencies operating at global and/or local-regional scale. In addition, by using the seismic wave measurements provided by reporting agencies, the ISC computes, where possible, its own event locations and magnitudes such as short-period body wave m _b and surface wave M _S . Therefore, the ISC Bulletin contains the results of the reporting agencies as well as the ISC own solutions. Among the most used seismic event parameters listed in seismological bulletins, the event magnitude is of particular importance for characterizing a seismic event. The selection of a magnitude value (or multiple ones) for various research purposes or practical applications is not always a straightforward task for users of the ISC Bulletin and related products since a multitude of magnitude types is currently computed by seismological agencies (sometimes using different standards for the same magnitude type). Here, we describe a scheme that we intend to implement in routine ISC operations to mark the preferred magnitudes in order to help ISC users in the selection of events with magnitudes of their interest.     

Event detectability of seismograph stations in Fennoscandia

The event detection probability as a function of ISC magnitude shows significant variations between seismic stations in Fennoscandia. In the teleseismic range, the 50% incremental detection thresholds extend from m_b = 4.5 to 5.5 among those stations that reported to the ISC during 1964–1969; the best stations being located in Finland. Most of this variation, though not all, can be explained from noise level fluctuation and regional signal attenuation effects. A common feature for all stations studied is a significantly better detectatbility for signals from eastern Asia than from western North America at a similar distance range. This is interpreted as resulting from higher attenuation in the mantle of high-frequency signals from the latter region.     

Extraction of soluble collagen and its feasibility in the palaeodietary research

In current palaeodietary research, gelatinization is the main method to extract insoluble collagen(ISC) from ancient bones. However, the degradation products of ISC, i.e., soluble collagen(SC), is often neglected and abandoned. In this work, we try to separate the extracts of ancient bones using gel chromatography and compare the contents of carbon and nitrogen, atomic C/N ratio, and stable carbon and nitrogen isotopic values of the extracts from three peaks to determine which peak can be attributed to SC. At last, the potential application of SC in palaeodietary research is discussed based on the comparison of stable isotopic values between ISC and SC. Among the three peaks, the second with the retention time between 17.5 min and 27.5 min had the most broad peak shape, indicating that the molecular weights of proteins collected were most variable. Besides, the contents of carbon and nitrogen and atomic C/N ratio of extracts in this peak were closest to the corresponding ISC. Based on the above, we conclude that the extract in second peak is SC. More important, the δ 13C and δ 15N values of ISC and SC are very similar. For ISC and SC with atomic C/N ratios within the normal range(2.9–3.6), the mean difference of δ 13C value was only(0.3±0.2)‰(n=2) while δ 15N value was(0.6±0.1)‰(n=2). Although the atomic C/N ratios of some SC are slightly beyond the normal range, the mean differences of δ 13C and δ 15N values were still only(0.4±0.1)‰ and(0.3±0)‰(n=2) respectively. These isotopic differences are quite below the isotope fractionation in one trophic level(δ 13C values of 1‰–1.5‰ and δ 15N values of 3‰–5‰), suggesting that SC had great application potentials in palaeodietary research.     

Grid-search event location with non-Gaussian error models

This study employs an event location algorithm based on grid search to investigate the possibility of improving seismic event location accuracy by using non-Gaussian error models. The primary departure from the Gaussian error model that is considered is the explicit use of non-Gaussian probability distributions in defining optimal estimates of location parameters. Specifically, the class of generalized Gaussian distributions is considered, which leads to the minimization of L_p norms of arrival time residuals for arbitrary p≥1. The generalized Gaussian error models are implemented both with fixed standard errors assigned to the data and with an iterative reweighting of the data on a station/phase basis. An implicit departure from a Gaussian model is also considered, namely, the use of a simple outlier rejection criterion for disassociating arrivals during the grid-search process. These various mechanisms were applied to the ISC phase picks for the IWREF reference events, and the resulting grid-search solutions were compared to the GT locations of the events as well as the ISC solutions. The results indicate that event mislocations resulting from the minimization of L_p residual norms, with p near 1, are generally better than those resulting from the conventional L₂ norm minimization (Gaussian error assumption). However, this result did not always hold for mislocations in event depth. Further, outlier rejection and iterative reweighting, applied with L₂ minimization, performed nearly as well as L₁ minimization in some cases. The results of this study suggest that ISC can potentially improve its location capability with the use of global search methods and non-Gaussian error models. However, given the limitations of this study, further research, including the investigation of other statistical and optimization techniques not addressed here, is needed to assess this potential more completely.     

Precursory Seismic Quiescence before the 1994 Kurile Earthquake (Mw = 8.3) Revealed by Three Independent Seismic Catalogs

—We have found that the M _w = 8.3 Kurile earthquake on October 4, 1994 followed an outstanding seismic quiescence starting 5–6 years before the mainshock near the ruptured area. We have analyzed three independent seismic catalogs Institute of Seismology and Volcanology, Hokkaido University (ISV), Japan Meteorological Agency (JMA) and International Seismology Center (ISC). In spite of selecting different magnitude bands and time windows all three catalogs presented the common feature of the seismic quiescence. This fact strongly suggests that the seismic quiescence should not be a man-made change but actually occurred. Moreover we have confirmed that the seismic quiescence was the most significant and the earthquake was the largest in the past twenty-five years in this region. Therefore we confidently interpret this seismic quiescence as an indication of a preparation process for the M _w = 8.3 Kurile earthquake.     

Effect of leaf distribution pattern on the interception storage capacity of leaf litter under simulated rainfall conditions

Rainwater interception by leaf litter is an important part of forest hydrological processes. The objective of this study was to investigate the interception storage capacity (ISC) of woodland leaf litter for three leaf distribution patterns, one flow path, two flow paths, and three flow paths, manually simulated via one-by-one leaf connection in the top leaf litter layer. A random pattern served as the control. Three different slopes (0°, 5° and 25°, representing flat, gentle and steep slopes, respectively) and two contrasting leaf litters (needle-leaf litter, represented by P. massoniana leaves, and broad-leaf litter, represented by C. camphora leaves) with a biomass of 0.5 kg/m² per unit area were applied, at a rainfall intensity of 50 mm/h. Results suggested that leaf distribution pattern greatly impacts litter drainage and hence affects leaf litter ISC. The delaying capacity of litter drainage initiation and ISC of broad-leaf litter were higher than those of needle-leaf litter under the same slope conditions. The maximum ISC (C_max) and minimum ISC (C_min) of leaf litter at flat and gentle slopes were higher than those at steep slope. C_min of the broad-leaf litter was two times higher than that of needle-leaf litter on average. When raindrops reached the litter layer, some were temporarily intercepted by the top litter layer while others infiltrated leaf litter sublayer along leaf edges, and in the process, some rainwater flowed through litter layer and contributed to lateral litter drainage along the potential flow path formed by leaves. The lateral litter drainage of broad-leaf litter was higher than that of needle-leaf litter, and the partitioning of rainwater into lateral litter drainage increased with increases in slope. The difference in leaf litter C_max among different slopes and leaf shapes decreased with flow path increasing. Therefore, leaf distribution pattern notably impact leaf litter ISC, which is similar to leaf shape and slope impacts. On inclined slopes, avoiding leaf accumulation to form flow path is helpful for improving ISC.     

Preliminary reference Earth model

A large data set consisting of about 1000 normal mode periods, 500 summary travel time observations, 100 normal mode Q values, mass and moment of inertia have been inverted to obtain the radial distribution of elastic properties, Q values and density in the Earth's interior. The data set was supplemented with a special study of 12 years of ISC phase data which yielded an additional 1.75 × 10⁶ travel time observations for P and S waves. In order to obtain satisfactory agreement with the entire data set we were required to take into account anelastic dispersion. The introduction of transverse isotropy into the outer 220 km of the mantle was required in order to satisfy the shorter period fundamental toroidal and spheroidal modes. This anisotropy also improved the fit of the larger data set. The horizontal and vertical velocities in the upper mantle differ by 2–4%, both for P and S waves. The mantle below 220 km is not required to be anisotropic. Mantle Rayleigh waves are surprisingly sensitive to compressional velocity in the upper mantle. High S_n velocities, low P_n velocities and a pronounced low-velocity zone are features of most global inversion models that are suppressed when anisotropy is allowed for in the inversion.The Preliminary Reference Earth Model, PREM, and auxiliary tables showing fits to the data are presented.     

Effects of fractional wettability on capillary pressure–saturation–relative permeability relations of two-fluid systems

Capillary pressure (P_c)–saturation (S)–relative permeability (k_r) relationships must be quantified to accurately predict non-aqueous phase liquid (NAPL) distribution in the subsurface. Several experimental techniques are presented here for two-fluid P_c–S–k_r relationships for various saturation paths to better define the effect of fractional wettability on these relationships. During the primary drainage path of the P_c–S curves, the air–water system showed no distinct trend as a function of the fraction of sand treated by organosilane (S) to render it non-water wetting. In a NAPL–water system, however, a consistent decrease of capillary pressure with increase of the fraction of non-water wetting sands was observed. The much lower contact angle for air–water (a–w) system may result in the observed insensitivity of the a–w P_c–S curves to fractional wettability, at least for the PD pathway. For the main imbibition path of NAPL–water system, capillary pressure decreased as the fraction of the S component increased, requiring forced imbibition (negative capillary pressures) for a certain range of saturations. Systems with an increasing percentage of the S component also exhibited a higher water k_r and lower NAPL or air k_r at a given saturation for the primary drainage and main imbibition paths in both air–water and NAPL–water systems. The increase of water k_r with increase of the fraction of the S component can be explained by the ability of water to occupy larger and highly conductive pores in such a system. Experimental k_r–S data for the primary drainage path of NAPL–water system presented here were used to test the Bradford et al. [Bradford SA, Abriola LM, Leij FJ. Wettability effects on two- and three-fluid relative permeabilities. J Contam Hydrol 1997;28:171–91] model and the modified Mualem model for estimating the k_r–S curves from measured P_c–S data as a function of fractional wettability. Both models predicted significantly less variation in the k_r–S curves than measured indicating that they did not adequately represent the system under investigation.     

Study of calibration function for surface wave magnitude of DK1 seismographs

ＳｔｕｄｙｏｆｃａｌｉｂｒａｔｉｏｎｆｕｎｃｔｉｏｎｆｏｒｓｕｒｆａｃｅｗａｖｅｍａｇｎｉｔｕｄｅｏｆＤＫ１ｓｅｉｓｍｏｇｒａｐｈｓＦＥＮＧＸＵＥ（薛峰）ＹＯＮＧＺＨＡＯ（赵永）ＣｅｎｔｅｒｆｏｒＡｎａｌｙｓｉｓａｎｄＰｒｅｄｉｃｔｉｏｎ，Ｓｔａｔ...     

Moment Magnitude (M W) and Local Magnitude (M L) Relationship for Earthquakes in Northeast India

An attempt has been made to examine an empirical relationship between moment magnitude (M _W) and local magnitude (M _L) for the earthquakes in the northeast Indian region. Some 364 earthquakes that were recorded during 1950–2009 are used in this study. Focal mechanism solutions of these earthquakes include 189 Harvard-CMT solutions (M _W?≥?4.0) for the period 1976–2009, 61 published solutions and 114 solutions obtained for the local earthquakes (2.0?≤?M _L?≤?5.0) recorded by a 27-station permanent broadband network during 2001–2009 in the region. The M _W–M _L relationships in seven selected zones of the region are determined by linear regression analysis. A significant variation in the M _W–M _L relationship and its zone specific dependence are reported here. It is found that M _W is equivalent to M _L with an average uncertainty of about 0.13 magnitude units. A single relationship is, however, not adequate to scale the entire northeast Indian region because of heterogeneous geologic and geotectonic environments where earthquakes occur due to collisions, subduction and complex intra-plate tectonics.     

A worldwide comparison of predicted S-wave delays from a three-dimensional upper mantle model with P-wave station corrections

This paper is a comparison, on a worldwide scale, between P-station corrections deduced from ISC residuals (Dziewonski and Anderson) and synthetic S-station corrections computed using a three-dimensional upper mantle model obtained from mantle wave data (Woodhouse and Dziewonski). The upper mantle S-velocity model is described by a spherical harmonic expansion up to degree 8; the P-station corrections are smoothed using a similar expansion, in order that the two data sets can be compared.Correlations between P-station corrections, δt_P, and synthetic S-station corrections, δt_S relevant to various depths of integration indicate that a station correction contains information about structures down to at least 670 km. For this depth of integration, the correlation coefficient of the two data sets is 0.59; the slope ‘a’ of the relation δt_S = aδt_P + b, obtained for the worldwide distribution of stations, is in good agreement with results of previous regional studies using direct readings of P and S arrival times (a = 3.61 ± 0.13).An analysis of regional variations of the relation δt_S = aδt_P + b is carried out on the basis of two published global tectonic patterns (Okal; Jodan). Results for oceanic regions are not reliable, due to the lack of data. On continental areas, a significant difference appears between mountains (a = 2.7 ± 0.3) and shields (a = 4.5 ± 0.4 for Okal's pattern, a = 5.7 ± 1.5 for Jordan's pattern). The largest a-value for shields rules out an explanation by partial melting, as proposed in previous studies. Thermal heterogeneities lead to low a-values; undulations of the lithosphere-asthenosphere boundary appear to be the most feasible explanation of the high slope beneath shields; they are also able to explain the range of variation of the station corrections; the lowest values of the station corrections correspond to a total vanishing of the low velocity zone beneath the oldest shields. For mountains, the mean values of the station corrections as well as the low a-value can be accounted for by a slight increase of Poisson's ratio together with a significant density increase.     

Revision of magnitudes of large shallow earthquakes, 1897–1912

In previous research, trace amplitudes of surface wave maxima recorded by undamped Milne seismographs were used to determine the surface-wave magnitudes M_s of large shallow earthquakes which occurred prior to 1912. For this purpose, the effective gain of these instruments was calibrated by using the surface-wave magnitudes M_s(GR) which were calculated from the unpublished worksheets for Seismicity of the Earth of Gutenberg and Richter. In this paper, the real quality of M_s(GR) is critically re-evaluated by using independent sets of data. It is found that M_s(GR) for the period 1904–1909 is considerably overestimated. The average excess from the real M_s is 0.5 units for 1904–1906, 0.4 for 1907, 0.3 for 1908–1909 and 0.0 for 1910–1912. This overestimation is so systematic and large that the previous results are all redetermined. The average effective gain of Milne instruments is revised to be 21.9; previously, the gain depended on M_s. This revision results in systematic reduction in the previously assigned magnitudes. The revised values of M_s for 264 shallow earthquakes, with M_s=6.8 and over in the period 1897–1912 inclusive, are listed. The present revision is large enough to preclude the possibility of the high activity of large shallow earthquakes around the turn of the century. The present results have a direct effect on all the magnitude catalogues of shallow earthquakes which occurred prior to 1909.     

Temporal variations of seismicity parameters in the Central Alborz,Iran

Temporal changes of b-value, fractal (correlation) dimensions of epicenters (D _e ) and occurrence time of earthquakes (D _t ) and relations between these parameters were calculated to investigate precursory changes before 28 May 2004, Baladeh-Kojour earthquake (M _w = 6.3) of Central Alborz, Iran. 2086 events with M _N ≥ 1.7 were selected for our analyses. A wide range of variation was seen in these parameters: b-value ～ 0.6–1.11, D _e ～ 0.97–1.64, and D _t ～ 0.13–0.93. The results showed decreases in all fractal parameters several months before the main shock. This decrease, which might have arisen due to clusters of events occurred between 2002–2003, was followed by a systematic increase, corresponding to the increased level of low-magnitude seismicity. It seems that changes in fractal parameters may be precursors of Baladeh-Kojour earthquake which was caused by seismic activation and quiescence. Furthermore, a positive correlation between b-value and D _e was detected before the main shock (D _e = 0.87 + 0.7b) and during aftershock sequences (D _e = 2b ± 0.09), which was further on changed to a negative one (D _e = 2.56–1.32b).     

Undrained behaviour of two silica sands and practical implications for modelling SSI in liquefiable soils

The aim of the present study is twofold. Firstly, the paper investigates the undrained cyclic and post-cyclic behaviour of two silica sands by means of multi-stage cyclic triaxial tests. Secondly, based on the post-cyclic response observed in the element test, the authors formulate a simplified stress–strain relationship that can be conveniently used for the construction of p–y curves for liquefiable soils. The multi-stage loading condition consists of an initial cyclic loading applied to cause liquefaction, followed by undrained monotonic loading that aimed to investigate the post-cyclic response of the liquefied sample. It was found that due to the tendency of the liquefied soil to dilate upon undrained shearing, the post-liquefaction strain–stress response was characterised by a distinct strain–hardening behaviour. The latter is idealized by means of a bi-linear stress–strain model, which can be conveniently formulated in terms of three parameters, i.e.: (i) take-off shear strain, γ_to, i.e. shear strain required to mobilize 1 kPa of shear strength; (b) initial secant shear modulus, G₁, defined as 1/γ_to; (c) post-liquefied shear modulus at large strain, G₂ (γ⪢γ_to). Based on the experimental results, it is concluded that these parameters are strongly influenced by the initial relative density of the sample, whereby γ_to decreases with increasing relative density. Differently both shear moduli (G₁ and G₂) increases with increasing relative density. Lastly, the construction of new p–y curves for liquefiable soils based on the idealized bi-linear model is described.     

Regional Variation of Lg Coda Q in the Continental United States and its Relation to Crustal Structure and Evolution

—Records from broadband digital stations have allowed us to map regional variations of Lg coda Q across almost the entire United States. Using a stacked ratio method we obtained estimates of Q ₀ (Lg coda Q at 1 Hz) and its frequency dependence, <eta>, for 218 event-station pairs. Those sets of estimates were inverted using a back-projection method to obtain tomographic images showing regional variations of Q ₀ and <eta>. Q ₀ is lowest (250–300) in the California coastal regions and the western part of the Basin and Range province, and highest (650–750) in the northern Appalachians and a portion of the Central Lowlands. Intermediate values occur in the Colorado Plateau (300–500), the Columbia Plateau (300–400), the Rocky Mountains (450–550), the Great Plains (500–650), the Gulf Coastal Plain and the southern portion of Atlantic Coastal Plain (400–500), and the portions of the Central Lowlands surrounding the high-Q region (500–550). The pattern of Q ₀ variations suggests that the United States can be divided into two large Q provinces. One province spans the area from the Rocky Mountains to the Atlantic coast, is tectonically stable, and exhibits relatively high Q ₀?. The other extends westward from the approximate western margin of the Rocky Mountains to the Pacific coast, is tectonically active, and exhibits low Q ₀?. The transition from high to low Lg coda Q in the western United States lies further to the west than does an upper mantle transition for Q and electrical resistivity found in earlier studies. The difference in Q ₀ between the western and eastern United States can be attributed to a greater amount of interstitial crustal fluids in the west. Regions of moderately reduced Q within the stable platform often occur where there are accumulations of Mesozoic and younger sediments. Reduced Q ₀ in the southeastern United States may not be due to anelasticity but may rather be explained by a gradational velocity increase at the crust-mantle boundary that causes shear energy to leak into the mantle.     

Coda-Q and its lapse time dependence analysis in the interaction zone of the Dinarides,the Alps and the Pannonian basin

The single backscattering model was used to estimate total attenuation of coda waves (Q_c) of local earthquakes recorded on eight seismological stations in the complex area of the western continental Croatia. We estimated Q₀ and n, parameters of the frequency dependent coda-Q using the relation Q_c = Q₀fⁿ. Lapse time dependence of these parameters was studied using a constant 30 s long time window that was slid along the coda of seismograms. Obtained Q_c were distributed into classes according to their lapse time, t_L. For t_L = 20–50 s we estimated Q₀ = 45–184 and n = 0.49–0.94, and for t_L = 60–100 s we obtained Q₀ = 119–316 and n = 0.37–0.82. There is a tendency of decrease of parameter n with increasing Q₀, and vice versa. The rates of change of both Q₀ and n seem to decrease for lapse times larger than 50–80 s, indicating an alteration in rock properties controlling coda attenuation at depths of about 100–160 km. A very good correlation was found between the frequency dependence parameter n and the Moho depths for lapse times of 50, 60 and 70 s.     

Changes in Earthquake Source Properties across a Shallow Subduction Zone: Kamchatka Peninsula

—This paper studies the source properties of earthquakes originating within the shallow subduction zone near Kamchatka Peninsula. We use the regional catalog of 1962–1993 Kamchatkan earthquakes completed by the Institute of Volcanology, Russia. Our previous investigations (Zobin, 1990, 1996a) and this study allow us to show a gradual change in source properties of earthquakes from trench to coast.¶It was demonstrated that the swarm sequences change to the mainshock–aftershock sequences from trench to coast. The source area of aftershock sequences is generally smaller than the swarm areas for the same magnitude M _s of the mainshock or clue event of the swarm. Study of the M _s –K _s relation, where K _s is the energy class for Kamchatka earthquakes, reveals that the events radiate relatively higher frequencies from trench to coast.     

Long-term earthquake prediction in the Marmara region based on the regional time- and magnitude-predictable model

In order to estimate the recurrence intervals for large earthquakes that occurred in the Marmara region, this region, limited with the coordinates of 39°–42°N, 25°–32°E, has been separated into seven seismogenic sources on the basis of certain seismological criteria, and regional time- and magnitude-predictable model has been applied for these sources. Considering the interevent time between successive mainshocks, the following two predictive relations were computed: log T _t = 0.26 M _min + 0.06 M _p –0.56 log M ₀ + 13.79 and M _f = 0.63 M _min ? 0.07 M _p + 0.43 log M ₀ ? 7.56. Multiple correlation coefficient and standard deviation have been computed as 0.53 and 0.35 for the first relation and 0.66 and 0.39 for the second relation, respectively. On the basis of these relations and using the occurrence time and magnitude of the last mainshocks in each seismogenic source, the probabilities of occurrence P(Δt) of the next mainshocks during the next five decades and the magnitude of the expected mainshocks were determined.     

Solar wind structure sources and periodicities of auroral electron power over three solar cycles

We assess the contributions of various types of solar wind structures (transients, coronal hole high-speed streams (HSS), and slow-speed wind) to hourly average auroral electron power (P_e). The time variation of the solar wind velocity (V_sw) and P_e are determined by HSS, which contribute ～47% to P_e and V_sw. Transients contribute ～42% of P_e in solar maxima, and ～6% in solar minimum. Cross-correlations of P_e with V_sw|B| for negative B_z are significant. P_e exhibits solar rotational periodicities similar to those for V_sw, with strong 7- and 9-day periodicities in 2005–2008 and equinox semiannual periodicities in 1995–1999.     

Spatial relation between source regions of strong earthquakes (M ≥ 7) on Southern Kamchatka and the distribution of velocities and elastic moduli in the Benioff zone

The presence of a phenomenological relationship between high velocity regions in the Benioff zone and sources of relatively strong earthquakes (M ≥ 6) was established for the first time from the comparison of such earthquakes with the velocity structure of central Kamchatka in the early 1970s. It was found that, in the region with P wave velocities of 8.1–8.5 km/s, the number of M ≥ 6 earthquakes over 1926–1965 was 2.5 times greater than their number in the region with velocities of 7.5–8.0 km/s. Later (in 1979), within the southern Kurile area, Sakhalin seismologists established that regions with V _P = 7.3–7.7 km/s are associated with source zones of M = 7.0–7.6 earthquakes and regions with V _P = 8.1–8.4 km/s are associated with M = 7.9–8.4 earthquakes. In light of these facts, we compared the positions of M = 7.0–7.4 earthquake sources in the Benioff zone of southern Kamchatka over the period 1907–1993 with the distribution of regions of high P velocities (8.0–8.5 to 8.5–9.0 km/s) derived from the interpretation of arrival time residuals at the Shipunskii station from numerous weak earthquakes in this zone (more than 2200 events of M = 2.3–4.9) over the period 1983–1995. This comparison is possible only in the case of long-term stability of the velocity field within the Benioff zone. This stability is confirmed by the relationship between velocity parameters and tectonics in the southern part of the Kurile arc, where island blocks are confined to high velocity regions in the Benioff zone and the straits between islands are confined to low velocity regions. The sources of southern Kamchatka earthquakes with M = 7.0–7.4, which are not the strongest events, are located predominantly within high velocity regions and at their boundaries with low velocity regions; i.e., the tendency previously established for the strongest earthquakes of the southern Kuriles and central Kamchatka is confirmed. However, to demonstrate more definitely their association with regions of high P wave velocities, a larger statistics of such earthquakes is required. On the basis of a direct correlation between P wave velocities and densities, the distributions of density, bulk modulus K, and shear modulus μ in the upper mantle of the Benioff zone of southern Kamchatka are obtained for the first time. Estimated densities vary from 3.6–3.9 g/cm³ in regions of high V _P values to 3.0–3.2 g/cm³ for regions of low V _P values. The bulk modulus K in the same velocity regions varies from (1.4–1.8) × 10¹² to (0.8–1.1) × 10¹² dyn/cm², respectively, and the shear modulus μ varies from (0.8–1.0) × 10¹² to (0.5–0.7) × 10¹² dyn/cm², respectively. Examination of the spatial correlation of the source areas of southern Kamchatka M = 7.0–7.4 earthquakes with the distribution of elastic moduli in the Benioff zone failed to reveal any relationship between their magnitudes and the moduli because of the insufficient statistics of the earthquakes used.