On the P-wave velocity and plate-tectonics implications for the Tyrrhenian deep-earthquake zone

P-wave velocities in the Tyrrhenian mantle have been determined for the 230–480 km depth range. Analysis of P-wave travel times for a set of Tyrrhenian deep earthquakes gives a velocity-distribution law which shows different behaviours in the 230–300 km and 300–480 km depth intervals. For the first interval the velocity gradient is 0.64 · 10⁻² sec⁻¹ and for the second one it is 0.59 · 10⁻² sec⁻¹. At a depth of 300 km the velocity decreases rapidly from 8.75 to 8.43 km/sec.The results have been analyzed in the framework of a Tyrrhenian structural model characterized by a lithospheric slab dipping 55–60° in the WNW direction.It is also pointed out that the analysis of some geodynamic features of the slabs of Pacific island arcs carried out by Oliver et al. (1973) and Sleep (1973) can be applied to the Tyrrhenian mantle geodynamic features.     

Recurrence of great earthquakes at subduction zones

Statistics of the recurrence times of great earthquakes at the Pacific subduction margins are made. The mean return period of great earthquakes is different from zone to zone, ranging from 27 to 117 years. The standard deviation of the return period proves to be very small, several years say, in some cases. The probabilities of a great earthquake recurring in each zone are estimated on the basis of Weibull distribution analysis.The mean return periods thus estimated are combined with the relative plate velocities at respective zones as obtained in the plate tectonics in order to estimate the ultimate displacement to rupture at the interface of the continental plate and the downgoing oceanic plate. It is presumed that great earthquakes at subduction zones occur as a result of a rebound of the continental plate at the time of rupture. The ultimate displacement thus estimated ranges from 2 to 8 m, and seems somewhat larger than that estimated on the basis of seismic observations, although the value of ultimate displacement seems to harmonize roughly with estimates based on geodetic observations on land. However, the ultimate displacement at the Aleutian—Alaska zone as estimated here seems much smaller than that estimated from actual observations.The ultimate strains, which are deduced from the displacements obtained on the assumption that the logarithmic extent of the deformed area is proportional to earthquake magnitude, are then calculated, and compared with those estimated for large inland earthquakes as revealed by repetition of geodetic surveys. The mean ultimate strain is estimated as 4.3 · 10⁻⁵ for subduction-zone earthquakes while that for inland earthquakes has been estimated as 4.7 · 10⁻⁵. As the agreement between both the ultimate strains is fairly good, it is tentatively concluded that the strength of the plate interface under the sea bottom is more or less the same as that in the crust on land.     

On the anomalous land uplift and aseismic creep prior to the 1976 Tangshan earthquake

A spatio-temporal analysis based on the data of eleven repeated levellings around the Tangshan region prior to the 1976 earthquake indicates that an uplift lasting for 2 years, from 1968 through 1969. with a magnitude of 50 mm, occurred in the epicentral area.Aseismic creep superimposed on the accumulated strain has been found in the vicinity of Tangshan and Baodi along both the Tangshan and the Jiyunhe faults.Assuming uniform strain accumulation and elastic dislocation, theoretical values of displacement at the various dislocation sites have been calculated and, using the least squares method, the optimal values of strain accumulation and the parameters of the creep faults in different years have been determined.The creep fault under Tangshan, a right-lateral normal fault, strikes N47°E and dips S87°E. and is 8 km long and 6 km wide. The upper boundary of the fault lies 2 km deep. The strike-slip and dip-slip offsets are, respectively, 104 cm and 8cm. The average rate of strain accumulation amounts to 0.9 × 10⁻⁷/yr. Creep at the fault amounted to 18.6 cm/yr and 1.4 cm/yr, respectively, in the strike and dip directions over the period 1969–1975. The Jiyunhe fault, although of smaller dimensions, has experienced a greater rate of creep than the Tangshan fault.A correlation of the above-mentioned uplift and creep with that of the Tangshan earthquake suggests that the uplift might have been a manifestation of the early development of the earthquake and that aseismic creep may be one of the precursory phenomena of shallow earthquakes. The sequence of processes preceding the Tangshan earthquake may be described as: strain accumulation-land upliftaseismic creep-inverse land deformation (or decrease in creep rate)-earthquake.     

Tectonics and crustal structures related to Bhuj earthquake of January 26, 2001: based on gravity and magnetic surveys constrained from seismic and seismological studies

Gravity and magnetic data of the Kachchh basin and surrounding regions have delineated major E–W and NW–SE oriented lineaments and faults, which are even extending up to plate boundaries in the north Arabian Sea and western boundary of the Indian plate, respectively. The epicentral zone of Bhuj earthquake and its aftershocks is located over the junction of Rann of Kachchh and median uplifts viz. Kachchh mainland and Wagad uplifts, which are separated by thrust faults. Gravity data with constraints from the results of the seismic studies along a profile suggest that the basement is uplifted towards the north along thrust faults dipping 40–60° south. Similarly gravity and magnetic modeling along a profile across Wagad uplift suggest south dipping (50–60°) basement contacts separating rocks of high susceptibility and density towards the north. One of these contacts coincides with the fault plane of the Bhuj earthquake as inferred from seismological studies and its projection on the surface coincides with the E–W oriented north Wagad thrust fault. A circular gravity high in contact with the fault in northern part of the Wagad uplift along with high amplitude magnetic anomaly suggests plug type mafic intrusive in this region. Several such gravity anomalies are observed over the island belt in the Rann of Kachchh indicating their association with mafic intrusions. The contact of these intrusives with the country rock demarcates shallow crustal inhomogeneities, which provides excellent sites for the accumulation of regional stress. A regional gravity anomaly map based on the concept of isostasy presents two centers of gravity lows of −11 to −13 mGal (10⁻⁵ m/s²) representing mass deficiency in the epicentral region. Their best-fit model constrained from the receiver function analysis and seismic refraction studies suggest crustal root of 7–8 km (deep crustal inhomogeneity) under them for a standard density contrast of −400 kg/m³. It is, therefore, suggested that significant amount of stress get concentrated in this region due to (a) buoyant crustal root, (b) regional stress due to plate tectonic forces, and (c) mafic intrusives as stress concentrators and the same might be responsible for the frequent and large magnitude earthquakes in this region including the Bhuj earthquake of January 26, 2001.     

Is the plastic flow uniformly distributed below the seismogenic region?

We compared the cutoff depth of seismicity in and around the Nojima fault broken by the 1995 Kobe earthquake occurring in intraplate Japan with the brittle–ductile transition depth of the widely accepted strength profile model of the crust. We successfully determined the temperature profile from borehole measurements, since almost the same geothermal gradients were observed at two boreholes located about 4 km apart from each other, and the thermal conductivity and heat production were also measured by taking numerous core samples. We found that the cutoff depth was much deeper than the transition depth under the assumption that wet granite is deformed at a strain rate of 3×10⁻¹⁵ s⁻¹. This small strain rate implies, however, that plastic flow is uniformly distributed below the seismogenic region. When the strain rate is assumed to be greater than 10⁻¹³ s⁻¹, the cutoff depth can be attributed to the transition depth. This suggests that deformation is localized in a narrow fault zone below the seismogenic region, even in the intraplate region.     

Thermal and mechanical structure of the central Iberian Peninsula lithosphere

The central Iberian Peninsula (Spain) is made up of three main tectonic units: a mountain range, the Spanish Central System and two Tertiary basins (those of the rivers Duero and Tajo). These units are the result of widespread foreland deformation of the Iberian plate interior in response to Alpine convergence of European and African plates. The present study was designed to investigate thermal structure and rheological stratification in this region of central Spain. Surface heat flow has been described to range from 80 to 60 mW m⁻². Highest surface heat flow values correspond to the Central System and northern part of the Tajo Basin. The relationship between elevation and thermal state was used to construct a one-dimensional thermal model. Mantle heat flow drops from 34 mW m⁻² (Duero Basin) to 27 mW m⁻² (Tajo Basin), and increases with diminishing surface heat flow. Strength predictions made by extrapolating experimental data indicate varying rheological stratification throughout the area. In general, in compression, ductile fields predominate in the middle and lower crusts and lithospheric mantle. Brittle behaviour is restricted to the first 8 km of the upper crust and to a thin layer at the top of the middle crust. In tension, brittle layers are slightly more extended, while the lower crust and lithospheric mantle remain ductile in the case of a wet peridotite composition. Discontinuities in brittle and ductile layer thickness determine lateral rheological anisotropy. Tectonic units roughly correspond to rheological domains. Brittle layers reach their maximum thickness beneath the Duero Basin and are of least thickness under the Tajo Basin, especially its northern area. Estimated total lithospheric strength shows a range from 2.5×10¹² to 8×10¹² N m⁻¹ in compression, and from 1.3×10¹² to 1.6×10¹² N m⁻¹ in tension. Highest values were estimated for the Duero Basin.Depth versus frequency of earthquakes correlates well with strength predictions. Earthquake foci concentrate mainly in the upper crust, showing a peak close to maximum strength depth. Most earthquakes occur in the southern margin of the Central System and southeast Tajo Basin. Seismicity is related to major faults, some bounding rheological domains. The Duero Basin is a relative quiescence zone characterised by higher total lithospheric strength than the remaining units.     

Strain partitioning between the Suruga Trough and the Zenisu Thrust: Neogene to present tectonic evolution in the onland and offshore Tokai district, Japan

The Tokai district in central Japan is located close to the convergent boundary between the Philippine Sea and Eurasian plates, and has experienced not only repeated large interplate earthquakes but also intense aseismic movement. In this paper, the spatial and temporal tectonic evolution of the Tokai district, particularly around the Omaezaki area, is discussed to assess whether the district has been and will be active or inactive. According to a geological survey, the horizontal crustal shortening strain can imply the hypothetical tectonic model that the area has been getting less active and the strain rate since the Neogene can be calculated as 12% and 2×10⁻⁶%/year, respectively. The present interseismic horizontal crustal strain and strain rate around the Omaezaki area are approximately 4×10⁻⁷% and 4×10⁻⁹%/year. By comparing these rates, the decrease since Neogene can imply the hypothetical tectonic model that the area has been getting less active influenced by the strain partitioning between the Suruga Trough and the Zenisu Thrust.     

Hydrogeochemical interpretations of apparent anomalies in base metals and radium in groundwater near Lake Maurice in the Great Victoria Desert

Apparently anomalous levels of Cu, Pb, Zn (up to 6.1, 26.0 and 10.8 mg 1⁻¹ respectively) and Ra (2000 pg 1⁻¹) have been noted in groundwaters from 28 drill holes within a 20 km × 20 km zone centred about a 10 times background airborne radiometric anomaly near Lake Maurice in the Great Victoria Desert in South Australia. Within 6 km of the anomaly centre the water table depth is generally less than 10 m, increasing to approximately 30 m in the drill holes furthest from the anomaly centre. All waters are very acid (pH 3.6 to 5.8) and deficient in carbonate species (all <0.5 mg 1⁻¹) but saturated with respect to calcium sulphate minerals. XRD traces of drill hole cuttings show the presence of quartz and halite at every sample site, pyrite at 75% of sites, variable amounts of kaolinite and muscovite at all sites, and variable amounts of feldspar, jarosite, calcium sulphate minerals, hydrated iron oxides, siderite, chlorite and calcite at certain locations only. Salinity of waters is very high ranging from that approximating sea water (Ionic strength (I) = 0.93 and Cl⁻ = 19 g 1⁻¹) to approximately six times sea water salinity (I = 4.61 and Cl⁻ = 120 g 1⁻¹).     

Ray path interpretation of the crustal structure beneath Saudi Arabia

Interpretation of a long-range seismic refraction line in Saudi Arabia has shown that beneath the Arabian Shield velocity generally increases with depth, from about 6 km s⁻¹ at the surface to about 7 km s⁻¹ at the top of the crust-mantle transition zone. The base of this transition zone (Moho) occurs at 37–44 km in depth. Intracrustal discontinuities can also be recognized, the most important being in the 10–20 km-depth range and separating the upper from the lower crust. Laterally, the variations in the intracrustal discontinuities and the total crustal thickness can be correlated with previously defined tectonic regions. Beneath the Red Sea shelf and coastal plain the crust, including 4 km of sediments, is only 15–17.5 km thick. With the aid of both seismic and gravity data an abrupt, steeply dipping transition from the crust of the Red Sea shelf and coastal plain to that of the Arabian Shield has been derived. With a jump of more than 20 km in Moho depth, this appears to be the major discontinuity between the Red Sea depression and the Arabian continental shield.     

Water movement within lac du bonnet batholith as revealed by detailed thermal studies of three closely-spaced boreholes

Successive temperature logs have been obtained over a period of two years in three closely-spaced boreholes in the Lac du Bonnet batholith of the Superior Province of the Canadian Shield. Two of the boreholes, of depth 450 m and 830 m, intersect a dipping fracture zone at 435–450 m. In both holes water is flowing from near the surface to the fracture zone at approximately 1.5–1.9·10⁻⁵ m³ s⁻¹, the flow being inferred from analysis of the temperature logs. Below 25 m, temperatures in these two holes are 0.22–0.28 K lower than those in the third, 145 m, hole.The temperature data have been combined with over 200 thermal conductivity measurements on core samples to produce heat flow values. In the deepest hole heat flow above the fracture zone is 16% higher than that below the zone. This indicates that water is flowing up the fracture zone. The flow rate is approximately 0.3 g s⁻¹ m⁻¹, and the flow has existed for thousands of years.Observation of thermal effects of water flow in massive, relatively unfractured plutons in a region having little topographic relief causes one to be concerned about the reliability of heat flow values measured in similar environments.The regional heat flow is taken to be 50 mW m⁻² after correction for glaciation effects. The average value of 24 determinations of radioactive heat generation in granitic core samples is 5.23 ± 1.11 μW m⁻³, which is more than three times higher than expected for such a heat flow in the Superior Province. This implies that the layer of high radioactive heat generation is thin, being not more than 4 km and probably about 1.3 km thick.     

Strain measurements and tectonics of New Zealand

Measurements of shear strain from triangulation data have been made at 30 locations in New Zealand. The standard error of measurement in terms of strain rate is about ±1 · 10⁻⁷ y⁻¹ and values of up to 7 · 10⁻⁷ y⁻¹ are observed. Together with 22 fault-plane solutions for crustal earthquakes the measurements indicate broad-scale patterns of deformation. Between the Hikurangi and Flordland active margins is a 100-km-wide belt, the axial tectonic belt, with shear strain rate averaging 5 ± 1 · 10⁻⁷y⁻¹ and an azimuth of the principal axis of compression of 114 ± 8°. The rate of movement (45 mm y⁻¹) and direction (085°) between the Pacific and Indian plates from the Minster et al. pole can be accounted for by the measured strain in the axial tectonic belt through simple shear parallel to, and compression normal to, the belt. The similarity in the rates determined from triangulation data averaged over 20–100 years and from plate movement averaged over 5 m.y. indicates plate movement to be uniform in time. West of the axial tectonic belt in Nelson and Fiordland are two zones in which movement is highly oblique to plate movement, and can be explained by slip line deformation analogous to the deformation of Asia. The azimuth of the principal axis of compression in the Taupo rift and East Cape region is NE—SW, perpendicular to its direction in the axial tectonic belt, suggesting extension in the rift and East Cape region normal to the subduction zone.     

Generalized geophysical model and dynamic properties of the continental crust

Detailed seismic investigations of the continental crust have produced evidence of definite regularities in the general layering of the consolidated crust despite its high degree of inhomogeneity. Three main layers may be resolved in the inner part of a continent: an upper layer with velocities of 5.8–6.4 km/s and a velocity gradient about 0.04–0.05 s⁻¹, an intermediate layer with velocities of 6.2–6.6 km/s and velocity gradient about zero, and a lower layer with velocities of 6.8–7.2 km/s and a high-velocity gradient of 0.05–0.1 s⁻¹. The intermediate layer is characteristically different not only because of its low average velocity gradient, but also because of its more pronounced horizontal layering, inversion zones, and its higher “transparency” and V_p/V_s ratio. The gravity and magnetic data have shown that basement inhomogeneities disappear at the top of the intermediate layer. Also there are few earthquakes in this layer. These pecularities may be interpreted as the result of partial melting (weakening) of rocks and their possible horizontal mobility inside this layer.Thus, dynamic models of tectonic processes must take into consideration the possible existence of a weak zone in the crust.     

The earthquake of March 4, 1977, its recording peculiarities and source parameters

A detailed analysis of recording peculiarities at seismic stations of the Uniform System of Seismic Observations (USSO) is presented a complicated nature of the source being shown. Consideration is given to parameters of the earthquake source, including the seismic moment and the length of the rupture.Comparison of magnitudes M_LH and M_PV indicates an anomalous attenuation in surface waves, itis is 3–4 times weaker than it had been noticed in case of other intermediate-depth Carpathian earthquakes.On the basis of comparison of the logarithm of the ratio of P-wave spectra at different epicentral distances (30° –70° ), the fac tor characterizing the absorption of P wave is found to remain practically unchanged.Average value of the seismic moment is estimated to be 2.6 × 10²⁷ dyne × cm, the most reasonable length of the rupture 58 km, and its focus 100 –130 km. The source parameters of the earthquake in question are compared with those of the earthquake of November 10, 1940.     

Induced earthquakes in the Izu peninsula by the Izu-Hanto-Oki earthquake of 1974, Japan

Recently small earthquakes in the Izu Peninsula, central Japan, occurred in a region where differential strain, or shear strain on the nodal planes, may have been enhanced by the Izu-Hanto-oki earthquake of 1974 (M = 6.9 after JMA). It is suggested that the seismic ctivity was induced by the redistribution of strain accompanying the Izu-Hanto-oki earthquake. The activity from August, 1975, may have also been affected by an abnormal uplift in the northeastern part of the peninsula. Based on plausible models, the uplift caused the accumulation of differential strain in the focal region of the subsequent earthquakes. Quantitatively, this change of crustal strain was of the order of 10⁻⁶; it is ten times as much as the average annual accumulation. Consequently, the sudden or rapid change of strain was likely to have played an essential role in the subsequent seismic activity. This effect could be one of the factors which trigger a shallow intra-plate earthquake.     

Source parameters of the Burma-India border earthquake of July 29, 1970, from body waves

The source parameters are determined for the Burma-India border earthquake of July 29, 1970, from body-wave spectra. We obtain seismic moment [ , ] · 10²⁶ dyne cm, source dimension [ ] km, radiated energy [ , −E_R (S) = 1.35] · 10²⁰ ergs and the stress drop = 11 bars.     

Earthquake forerunner as probable precursor — an example from north burma subduction zone

The Burmese Arc seismic activity is not uniform for its ∼ 1100 km length; only the Northern Burmese Arc (NBA) is intensely active. Six large earthquakes in the magnitude range 6.1–7.4 have originated from the NBA Benioff zone between 1954–2011, within an area of 200 × 300 km² where the Indian plate subducts eastward to depths beyond 200 km below the Burma plate. An analysis on seismogenesis of this interplate region suggests that while the subducting lithosphere is characterized by profuse seismicity, seismicity in the overriding plate is rather few. Large earthquakes occurring in the overriding plate are associated with the backarc Shan-Sagaing Fault (SSF) further east. The forecasting performance of the Benioff zone earthquakes in NBA as forerunner is analysed here by: (i) spatial earthquake clustering, (ii) seismic cycles and their temporal quiescence and (iii) the characteristic temporal b-value changes. Three such clusters (C1–C3) are identified from NBA Benioff Zones I & II that are capable of generating earthquakes in the magnitude ranges of 7.38 to 7.93. Seismic cycles evidenced for the Zone I displayed distinct quiescence (Q₁, Q₂ and Q₃) prior to the 6^th August 1988 (M 6.6) earthquake. Similar cycles were used to forecast an earthquake (Dasgupta et al. 2010) to come from the Zone I (cluster C1); which, actually struck on 4 February 2011 (M 6.3). The preparatory activity for an event has already been set in the Zone II and we speculate its occurrence as a large event (M > 6.0) possibly within the year 2012, somewhere close to cluster C3. Temporal analysis of b-value indicates a rise before an ensuing large earthquake.     

Secular change of permeability in the fracture zone near the Nojima fault estimated using strain changes due to water injection experiments

Water injection experiments were performed in 1997, 2000 and 2003 at the 1800 m borehole near the fracture zone of the 1995 Hyogo-ken Nanbu earthquake. During these experiments, a contraction of about 10^− 8–10^− 7 was observed with three-component strainmeters at a bottom of the 800 m borehole, 70 m southwest of the 1800 m borehole. We estimated hydraulic properties of the fracture zone near the Nojima fault by using the strain data to investigate a healing of the fault during the postseismic stage. We calculated pore pressure changes due to the water injection using Darcy's equation and obtained strain changes due to the pore pressure changes as elastic deformations of the crust. The calculated strain changes have a nearly agreement with the observed strain changes. Hydraulic conductivity in 1997, 2000 and 2003 was determined to be 0.9 ± 0.2 × 10^− 6, 0.8 ± 0.2 × 10^− 6 and 0.4 ± 0.1 × 10^− 6 m/s, respectively. The reduced hydraulic conductivities in 2000 and 2003 suggest that the fractures had been healing.     

3D modeling of earthquake cycles of the Xianshuihe fault,southwestern China

We perform 3D modeling of earthquake generation of the Xianshuihe fault, southwestern China, which is a highly active strike-slip fault with a length of about 350 km, in order to understand earthquake cycles and segmentations for a long-term forecasting and earthquake nucleation process for a short-term forecasting. Historical earthquake data over the last 300 years indicates repeated periods of seismic activity, and migration of large earthquake along the fault during active seismic periods. To develop the 3D model of earthquake cycles along the Xianshuihe fault, we use a rate- and state-dependent friction law. After analyzing the result, we find that the earthquakes occur in the reoccurrence intervals of 400–500 years. Simulation result of slip velocity distribution along the fault at the depth of 10 km during 2694 years along the Xianshuihe fault indicates that since the third earthquake cycle, the fault has been divided into 3 parts. Some earthquake ruptures terminate at the bending part of the fault line, which may means the shape of the fault line controls how earthquake ruptures. The change of slip velocity and displacement at 10 km depth is more tremendous than the change of the shallow and deep part of the fault and the largest slip velocity occurs at the depth of 10 km which is the exact depth of the seismic zone where fast rupture occurs.     

Electro-magnetic earthquake bursts and critical rupture of peroxy bond networks in rocks

We propose a mechanism for the low frequency electromagnetic emissions and other electromagnetic and electric phenomena which have been associated with earthquakes. The mechanism combines the critical earthquake concept and the concept of crust acting as a charging electric battery under increasing stress. The electric charges are released by activation of dormant charge carriers in the oxygen anion sublattice, called peroxy bonds or positive hole pairs (PHP), where a PHP represents an O₃X/^OO YO₃ with X,Y = Si⁴⁺,Al³⁺,…, i.e. O⁻ in a matrix of O²⁻ of silicates. We propose that PHP are activated by plastic deformations during the slow cooperative build-up of stress and the increasingly correlated damage culminating in a large “critical” earthquake. Recent laboratory experiments indeed show that stressed rocks form electric batteries which can release their charge when a conducting path closes the equivalent electric circuit. We conjecture that the intermittent and erratic occurrences of EM signals are a consequence of the progressive build-up of the battery charges in the Earth crust and of their release when crack networks percolate through the stressed rock volumes, providing a conductive pathway for the battery currents to discharge. EM signals are thus expected close to the rupture, either slightly before or after, that is, when percolation is most favored. The proposed mechanism should be relevant for the broader understanding of fractoemissions.     

The seismicity of Kohistan, Pakistan: Source studies of the Hamran (1972.9.3), Darel (1981.9.12) and Patan (1974.12.28) earthquakes

Long period body waves are examined to show that the Hamran (1972.9.3), Darel (1981.9.12) and Patan (1974.12.28) earthquakes in Kohistan had focal depths of about 8–10 km. All involved high angle reverse faulting (thrusting) and had seismic moments of about 2.2 to 2.7·10²⁵ dyne cm. These shallow depths contrast with the deeper hypocentres found in the Hindu Kush and northeast Karakoram to the north and in Hazara to the south. The Hamran and Patan shocks were assigned depths of 45 km by the ISC, indicating that even well-recorded events in this region may have focal depths in error by 30 km