Simultaneous inversion of the aftershock data of the 1993 Killari earthquake in Peninsular India and its seismotectonic implications

The aftershock sequence of the September 30th, 1993 Killari earthquake in the Latur district of Maharashtra state, India, recorded by 41 temporary seismograph stations are used for estimating 3-D velocity structure in the epicentral area. The local earthquake tomography (LET) method of Thurber (1983) is used. About 1500P and 1200S wave travel-times are inverted. TheP andS wave velocities as well asV _P/V_Sratio vary more rapidly in the vertical as well as in the horizontal directions in the source region compared to the adjacent areas. The main shock hypocentre is located at the junction of a high velocity and a low velocity zone, representing a fault zone at 6–7 km depth. The estimated average errors ofP velocity andV _P/V_Sratio are ±0.07 km/s and ±0.016, respectively. The best resolution ofP and S-wave velocities is obtained in the aftershock zone. The 3-D velocity structure and precise locations of the aftershocks suggest a ‘stationary concept’ of the Killari earthquake sequence.     

Compressional wave velocity of granite and amphibolite up to melting temperatures at 1 GPa

Compressional wave velocities (V_P) at above-solidus temperatures and at 1 GPa were obtained for a granite and amphibolite, which are considered to be major constituents of the continental crust. The temperature variation of velocities showed that the V_P values of granite decreased with rising temperature, but substantially increased beyond the melting temperature (850–900 °C). Such an increase may be caused by the α–β transition of quartz. The velocities of amphibolite decreased linearly with increasing temperature and dropped sharply at temperatures above the solidus (700 °C), indicating that partial melting of amphibolite acts to significantly lower the seismic velocities.     

Compressional and shear wave velocities of igneous rocks and volcanic glasses to 900° C and 20 kbar

Simultaneous measurements of compressional and shear wave velocities, V_p and V_s, in acidic and basic igneous rocks and volcanic glasses, were made up to 900°C and at 10–20 kbar.The effects of pressure and temperature on V_p and V_s in glasses and glassy rocks change at about 600°C, presumably the glass transition temperature. These effects are directly related to the silica content in the samples. and for obsidian are negative at room temperature and 245°C, but are positive at 655°C. The velocity—pressure relations for obsidian display an obvious hysteresis phenomena. for basalt glass is slightly negative, but is positive for usual substances at room temperature, and for obsidian and glassy andesite are positive up to about 600°C but are negative above that temperature. However, for basalt glass as well as other crystalline rocks, and are negative at all temperatures. Glass once heated above the glass transition temperature T_g under pressure P₁ retains the memory of pressure P₁ after it is cooled down below T_g and while subjected to another pressure P₂. An abrupt shift of the velocities correlating to pressure P₂ occurs when the glass is again heated to T_g. V_p−T and V_s−T relations for obsidian, glassy andesite, and basalt glass clearly exhibit this pressure memory.     

Iron-Titanium Oxide Geothermometry and Petrogenesis of Lava Flows and Dykes from Mandla Lobe of the Eastern Deccan Volcanic Province, India

Compositional studies on different forms of magnetite, ulvospinel, ilmenite and hematite mineral phases occurring in 37 lava flows and 6 dykes of the Mandla lobe are presented in this paper. Ilmenite (0001) in equilibrium with titanomanetite show high values of temperature of equilibration, ranging from 1172–974°C, for high alumina quartz normative tholeiitic lava flows of Chemical Type - A; 1129–1229°C for low alumina quartz normative tholeiitic lava flows of Chemical Type - B; 1283–1124°C for tholeiitic lava flows of Chemical Type - F and 1243°C and 99O°C for two diopside olivine normative tholeiite flows of Chemical Type D. High olivine normative flows of Chemical Type - G and H show 1095°C and 1092°C respectively. Whereas, high hypersthene normative tholeiite flow of Chemical me C shows temperature of 1187°C. Data plots disposition over iron-titanium oxide equilibration temperature vs – logfo₂, diagram for Mandla lava flows and other parts of the Deccan (Igatpuri, Mahabaleshwer, Nagpur and Sagar areas) revealed that tholeiitic (evolved) basalt of the eastern Deccan volcanic province formed at high temperatures whereas, picritic (primitive) lavas of Igatpuri and tholeiitic basalt of Mahabaleshwar areas were formed at low temperatures. Mahabaleshwer basalts follow FMQ (fayalite-magnetite-quartz) buffer curve but, plots of the Mandla basalts lie above this curve indicating higher temperatures of crystallisation of ilmenite-titanomagnetite than that of the lava flows from other parts of Deccan 'Raps. The eastern Deccan Traps are most evolved types of lava as characterised by its low Mg-number and Ni content whereas, Igatpuri lava flows are picritic (primitive), having high Mg-number and Ni contents. Temperature vs FeO + Fe₂O₃ / FeO + Fe₂O₃ + MgO ratio data plots for Mandla and other Deccan lava flows and liquidus data for Hawaiian tholeiites, indicated that Igatpuri basalts lie parallel to the liquidus line of Hawaiian tholeiite but at lower temperatures. Large data plots of Mandla lava flows lie along the liquidus line of the Hawaiian lava. The highly vesicular nature of compound lava flows having large amount of volatile is responsible for low temperature values whereas, lava flows represented by high temperatures show high modal values of glass and opaque minerals.     

Effect of high-low quartz transition on compressional and shear wave velocities in rocks under high pressure

The compressional and shear wave velocities in quarzite, granite, and granulite are determined at a fixed confining pressure of 2 kb as a function of temperature up to 720° C. The high-low quartz transition of the constituent quartz minerals is associated with a pronounced decrease in velocity of the compressional waves when approaching the transition and with a significant velocity increase after the transition. In contrast, the effect of the α-β transition on shear wave velocities is small. The drop of V _P is explained by the elastic softening of structure of the constituent quartz minerals near the α-β transition and the opening of grain-boundary cracks, caused by the very high volumetric thermal expansion of the quartz relative to the other component minerals. The velocity increase in the β-field may be attributed to an elastic hardening of the quartz structure. Poisson ratios computed from the velocity data are anomalous for a solid: they become negative within the transition regime. The transition temperature, as indicated by the minimum velocities, is higher in the polycristalline rocks than is expected on grounds of single crystal behavior, and the discrepancy is more marked in granite than in quartzite. The shift of the transition temperature to higher values is explained by internal stresses that arise from the anisotropy of the thermal expansion and compressibility of individual grains and the differences in thermal expansion and compressibility between different component minerals. The role of the α-β quartz transition as a possible cause of low-velocity layers is discussed.     

3-D crustal structure of the extensional Granada Basin in the convergent boundary between the Eurasian and African plates

Three-dimensional P and S wave velocity models of the crust under the Granada Basin in Southern Spain are obtained with a spatial resolution of 5 km in the horizontal direction and 2 to 4 km in depth. We used a total of 15407 P and 13704 S wave high-quality arrival times from 2889 local earthquakes recorded by both permanent seismic networks and portable stations deployed in the area. The computed P and S wave velocities were used to obtain three-dimensional distributions of Poisson's ratio (σ) and the porosity parameter (V_p×V_s). The 3-D velocity images show strong lateral heterogeneities in the region. Significant velocity variations up to ±7% in P and S velocities are revealed in the crust below the Granada Basin. At shallow depth, high-velocity anomalies are generally associated with Mesozoic basement, while the low-velocity anomalies are related to the neogene sedimentary rocks. The south–southeastern part of the Granada Basin exhibits high σ values in the shallowest layers, which may be associated with saturated and unconsolidated sediments. In the same area, V_p×V_s is high outside the basin, indicating low porosity of the mesozoic basement. A low-velocity zone at 18-km depth is found and interpreted as a weak–ductile crust transition that is related to the cut-off depth of the seismic activity. In the lower crust, at 34-km depth, a clear slow V_p and V_s anomalous zone may indicate variations in lithology and/or with the rigidity of the lower crust rocks.     

Primary basalts and magma genesis

The Pliocene-Holocene lavas of the Snake River Plain, Idaho, U.S.A., have a bimodal composition range, consisting predominantly of basalts (olivine-tholeiites), with subordinate intercalated tholeiitic andesites but with very few analyses falling between these groups. The more-magnesian of the tholeiitic andesites contain more total Fe, alkalis, TiO₂ and P₂O₅ but less SiO₂ than the less-magnesian basalts. Derivation of the tholeiitic andesites from the basalts by low-pressure fractional crystallization or by major-element crustal contamination does not seem possible, although some minor-element exchange with ancient crust apparently has occurred. Two lavas, representative of the least-magnesian basalts and the most-magnesian tholeiitic andesites, respectively, have been subjected to anhydrous experimental studies within their melting ranges at pressures up to 35kb. Both appear to show four-phase points on their liquidi at about 8kb and these are thought to have genetic significance. Microprobe analyses of the interstitial glasses in partially-crystalline runs on the basalt between 8 and 12kb show that these reproduce all the characteristic features of the Snake River Plain most-magnesian tholeiitic andesites, notably their reduced Si-saturation. The compositions of the most Mg-rich Snake River Plain basalts are such that they may perhaps be primary magmas, produced by partial fusion of a relatively Fe-rich spinel-lherzolite upper mantle at 50 to 60km depth; a proposal which accords well with the geophysics of this currently-active region. Partial crystallization of batches of this magma, delayed during ascent within the crust at depths of about 30 km, is thought to have given rise to the tholeiitic andesites.     

The effect of high temperature and high confining pressure on compressional wave velocities in quartz-bearing and quartz-free igneous and metamorphic rocks

Compressional wave velocities have been measured in granite, granulite, amphibolite and peridotite specimens under conditions of high temperature up to 700°C and confining pressures up to 6 kbar. In general, velocity increases with pressure and decreases with temperature.Quartz-bearing rocks show an anomalous behavior of their compressional wave velocities. The velocity—temperature relations exhibit a velocity-“deep” due to the high—low inversion of the constituent quartz crystals. The intrinsic effect of temperature on velocities is hard to determine due to thermal expansion and consequent loosening of the structure. The opening of new cracks and the widening of old cracks causes a large decrease in compressional wave velocities. The minimum pressure to prevent damage at a given temperature should, therefore, be about 1 kbar/100°C.The values obtained at these conditions are considered to be most nearly correct as intrinsic properties of the compact aggregates. Velocity anisotropies at high confining pressures and high temperatures correlate with preferred lattice orientation of the constituent minerals. The effect of dimensional orientation and microcracks on seismic anisotropy seems to be of minor importance in dry rocks. It is the more eliminated the higher the confining pressure. The data do not support the concept of a velocity maximum in depth of 10–20 km.     

Lower crustal fluid distribution in the northeastern Japan arc revealed by high-resolution 3D seismic tomography

The Ou Backbone Range strikes northwards through the central northeastern Japan arc and is bounded on both sides by the active reverse Uwandaira and Sen'ya faults. We have applied a traveltime inversion method (seismic tomography) with spatial velocity correlation to active and passive seismic data in order to investigate a three-dimensional (3-D) velocity structure. The data set contains 33,993 P- and 18,483 S-wave arrivals from 706 natural sources and 40 blasts, as well as 2803 P-wave traveltime data from 10 explosions detonated during the 1997 controlled source experiment. The traveltime inversion reveals a zone beneath the Ou Backbone Range in which P-wave velocities (V_P) are approximately 6–8% lower than the average velocity at equivalent depths. The low V_P and a low V_P to S-wave velocity (V_S) ratio (V_P/V_S) of about 1.65 suggest the presence of aqueous fluids in the middle crust.     

S-waves profiles from noise cross correlation at small scale

Cross correlations of ambient seismic noise recordings are studied to infer shear-wave velocities (V_S) profiles versus depth. Experiments are performed, over inter-station distances ranging from 60 m to about 400 m, in the crowded and noisy historical centre of Napoli, where the conditions for the use of active (controlled source) seismic spreadings are prohibitive, even for just one receiver. From the noise cross correlation, group velocity dispersion curves are extracted, with FTAN method, and their non-linear inversion supplies V_S profiles versus depth. The information of nearby stratigraphies and the range of V_S variability for samples of Neapolitan soils and rocks confirms the validity and reliability of the results obtained with our expeditious procedure. If successfully tested in other geological settings, the proposed approach is a low cost methodology to get reliable V_S velocity profiles versus depth.     

Directional properties of small strain shear stiffness in soils

A compiled database of shear wave velocity measurements in a variety of clays, silts and sands shows directional hierarchies between downhole (V_sVH), standard crosshole (V_sHV), and rotary crosshole (V_sHH) tests. The special in situ database has been collected from 33 well-documented geotechnical test sites. Expressions relating the small-strain shear modulus in terms of effective confining stress level, stress history and void ratio are explored for each of these three modes of directional shear wave velocity. The relationships are separated initially into soil groups (intact clays, fissured clays, sands and silts), and then generalised to consider all soil types together.     

Seismic velocities of granulite-facies xenoliths from Central Ireland: Implications for lower crustal composition and anisotropy

V_p and V_s values have been measured experimentally and calculated for granulite-facies lower crustal xenoliths from central Ireland close to the Caledonian Iapetus suture zone. The xenoliths are predominantly foliated and lineated metapelitic (garnet–sillimanite–K-feldspar) granulites. Their metapelitic composition is unusual compared with the mostly mafic composition of lower crustal xenoliths world-wide. Based on thermobarometry, the metapelitic xenoliths were entrained from depths of c. 20–25 ± 3.5 km and rare mafic granulites from depths of 31–33 ± 3.4 km. The xenoliths were emplaced during Lower Carboniferous volcanism and are considered to represent samples of the present day lower crust.V_p values for the metapelitic granulites range between 6.26 and 7.99 km s^{− 1} with a mean value of 7.09 ± 0.4 km s^{− 1}. Psammite and granitic orthogneiss samples have calculated V_p values of 6.51 and 6.23 km s^{− 1}, respectively. V_s values for the metapelites are between 3.86 and 4.34 km s^{− 1}, with a mean value of 4.1 ± 0.15 km s^{− 1}. The psammite and orthogneiss have calculated V_s values of 3.95 and 3.97 km s^{− 1}, respectively.The measured seismic velocities correlate with density and with modal mineralogy, especially the high content of sillimanite and garnet. V_p anisotropy is between 0.15% and 13.97%, and a clear compositional control is evident, mainly in relation to sillimanite abundance. Overall V_s anisotropy ranges from 1% to 11%. Poisson's ratio (σ) lies between 0.25 and 0.35 for the metapelitic granulites, mainly reflecting a high V_p value due to abundant sillimanite in the sample with the highest σ. Anisotropy is probably a function of deformation associated with the closure of the Iapetus ocean in the Silurian as well as later extension in the Devonian. The orientation of the bulk strain ellipsoid in the lower crust is difficult to constrain, but lineation is likely to be NE–SW, given the strike-slip nature of the late Caledonian and subsequent Acadian deformation.When corrected for present-day lower crustal temperature, the experimentally determined V_p values correspond well with velocities from the ICSSP, COOLE I and VARNET seismic refraction lines. Near the xenolith localities, the COOLE I line displays two lower crustal layers with in situ V_p values of 6.85–6.9 and 6.9–8.0 km s^{− 1}, respectively. The upper (lower velocity) layer corresponds well with the metapelitic granulite xenoliths while the lower (higher velocity) layer matches that of the basic granulite xenoliths, though their metamorphic pressures suggest derivation from depths corresponding to the present-day upper mantle.     

Geochemistry, petrofabrics and seismic properties of eclogites from the Chinese Continental Scientific Drilling boreholes in the Sulu UHP terrane, eastern China

We present an integrated study of geochemistry, petrofabrics and seismic properties of strongly sheared eclogites from the Chinese Continental Scientific Drilling (CCSD) project in the Sulu ultrahigh-pressure (UHP) metamorphic terrane, eastern China. First, geochemical data characterize diverse protoliths of the studied eclogites. The positive Eu- and Sr-anomalies, negative Nb anomaly and flat portion of heavy rare earth elements in coarse-grained rutile eclogites (samples B270 and B295) suggest a cumulate origin in the continental crust, whereas the negative Nb anomaly and enrichment of light rare earth elements in retrograde eclogites (samples B504, B15 and B19) imply an origin of continental basalts or island arc basalts. Second, P-wave velocities (V_p) of three typical eclogite samples were measured under confining pressures up to 500 MPa and temperatures to 700 °C. At 500 MPa and room temperature, the mean V_p reaches 8.50-8.53 km/s in samples B270 and B295 but drops to 7.86 km/s in sample B504, and the P-wave anisotropy changes from 1.7-2.7% to 5.5%, respectively. The pressure and temperature derivatives of V_p are larger in the retrograde eclogite than in fresh ones. Third, the electron backscatter diffraction (EBSD) measurements of the eclogites reveal random crystal preferred orientation (CPO) of garnet and pronounced CPO of omphacite, which is characterized by a strong concentration of [001]-axes sub-parallel to the lineation and of (010)-poles perpendicular to the foliation. The asymmetric CPO of omphacite in sample B270 recorded a top-to-the-south shear event during subduction of the Yangtze plate. The calculated fastest V_p is generally sub-parallel to the lineation, but a different deformation environment during exhumation could form second-order variations in omphacite CPO and affect the V_p distribution in eclogites (e.g., the fastest V_p is at ~ 35° from the foliation in sample B295). Comparison between measured and calculated seismic properties indicates that the CPO of omphacite controls the seismic anisotropy of eclogites at high pressure, and compositional layering and retrograde minerals will increase the anisotropy. Calculated P-wave velocities agree well with velocities measured at 500 MPa and room temperature for fresh eclogites, but much higher than those of retrograde eclogite. As a case study, the laboratory-derived V_p-P and V_p-T relationships were used to estimate P-wave velocities of eclogites and peridotites beneath the Western Superior Province, Canada. The results indicate that besides the fabric-induced anisotropy, the direction dependence of pressure and temperature derivatives of V_p can significantly increase seismic anisotropy of eclogites with depth, which results in eclogites being an important candidate for the seismic anisotropy in the upper mantle. Due to their very high density and velocity, garnet-rich eclogites within peridotite could be detected in seismic reflections in subduction zones.     

Estimation of shearwave velocity in drifts using multichannel analysis of surface wave (MASW) technique — A case study from Jammu & Kashmir,India

Multichannel analysis of surface waves (MASW) is a non-destructive seismic prospecting method utilizing Rayleigh waves for imaging and characterizing shallow sub-surface structure. Multichannel analysis of surface waves (MASW) studies were conducted in drift areas of two bridge sites in the hilly terrain of J&K for imaging and characterizing shallow sub-surface structure. The purpose of the present study is to estimate the shear wave velocity (V_S) and subsurface structure in four drifts made in a hilly terrain for construction of two bridges. Rayleigh waves are having dispersive properties, travelling along or near the ground surface and are usually characterized by relatively low velocity, low frequency, and high amplitude. The study area comprises of Tertiary group of rocks which are underlain by Siwalik group. The main rock type in the study area is dolomite which has undergone various geological processes like weathering, jointing, fracturing and shearing. MASW data was collected inside four drifts in the mountainous terrain of J&K state which are located on either sides of Chenab river. The data was analyzed by relevant processing software using dispersion and inversion technique. Shear wave velocities were estimated up to 30 m depth. Average shear wave velocity (V_S ³⁰) up to top 30m was also computed. It is observed that, V_S in the range 400–800 m/s upto 10–15 m corresponding to weathered rock, followed by compact dolomite rock up to the depth of about 30 m with V_S in the range 1200–1600 m/s. Some low velocity zones are also identified from these sections which represent shear zones.     

Surface wave studies for shear wave velocity and bedrock depth estimation over basalts

Shear wave velocity (V _S) estimation is of paramount importance in earthquake hazard assessment and other geotechnical/geo engineering studies. In our study, the shear wave velocity was estimated from ground roll using multichannel analysis of surface wave (MASW) technique making use of dispersive characteristics of Rayleigh type surface waves followed by imaging the shallow subsurface basaltic layers in an earthquake-prone region near Jabalpur, India. The reliability of MASW depends on the accurate determination of phase velocities for horizontally traveling fundamental mode Rayleigh waves. Inversion of data from surface waves resulted in a shear wave velocity (V _S) in the range of 200–1,200 m/s covering the top soil to weathering and up to bedrock corresponding to a depth of 10–30 m. The P-wave velocity (V _P) obtained from refraction seismic studies at these locations found to be comparable with V _S at an assumed specific Poisson’s ratio. A pair of selected set of V _S profiles over basalt which did not result in a hazardous situation in an earthquake of moderate magnitude are presented here as a case study; in other words, the shear wave velocity range of more than 200 m/s indicate that the area is highly unlikely prone to liquefaction during a moderate or strong earthquake. The estimated depth to basalt is found to be 10–12 m in both the cases which is also supported by refraction studies.     

A note on the petrology of some lava types from East New Guinea

Quaternary and upper Tertiary lavas from the volcanic islands along the north coast of east New Guinea consist of pyroxene andesites, andesitic basalts, and basalts; many contain accessory olivine. The lavas of Mts Lamington, Yelia, and Victory, in southeast Papua and New Guinea, are rather more acid hornblende and lamprobolite andesites. Peralkaline obsidian and rhyolite occur at Fergusson and Dobu Islands, east of the Papuan mainland. Chemical analyses show that most of the lavas are of calc‐alkaline type, although some from Karkar and Long Islands appear to have tholeiitic tendencies.     

Phase relations of a high-Mg basalt from the Aleutian Island Arc: Implications for primary island arc basalts and high-Al basalts

Many volcanic centers in the Aleutian Islands have erupted lavas that range in composition from high-Mg basalt (MgO>9 wt%) to more fractionated and voluminous high-Al basalts and basaltic andesites. The petrogenetic relationships between these rock types and the composition of primary magmas has been vigorously debated. The phase relations of a typical high-Mg basalt from the Makushin volcanic field on Unalaska Island provide important constraints on petrogenetic models. Results of one-atmosphere and moderate-to high-pressure (5–20 kb) anhydrous experiments are similar to results obtained from primitive MORB. At low pressures olivine is the liquidus phase joined by plagioclase and clinopyroxene at progressively lower temperatures. Clinopyroxene is the second phase to crystallize at pressures greater than 5 kb and replaces olivine on the liquidus at approximately 10 kb. Above 10 kb the liquidus pyroxene is aluminous augite and orthopyroxene is the second phase to crystallize. Glasses in equilibrium with olivine and clinopyroxene at intermediate-pressure (5 to 10 kb) are similar in composition to high-Al basalt. Plagioclase is not involved and most likely does not become a liquidus phase until the liquid has evolved significantly. Although our studies do not confirm the primary nature of high-Mg basalts they do support a model in which high-Al basalts are generated by moderate amounts of crystal fractionation from more primitive (high Mg/Mg+Fe, lower Al₂O₃) basaltic magmas near the arc crust-mantle boundary.Abbreviations Ol olivine - Cpx Clinopyroxene - Pl plagioclase - L liquid - Sp spinel - Pig pigeonite - Opx Orthopyroxene     

Liquidus phase relations on the join diopside-forsterite-anorthite from 1 atm to 20 kbar: Their bearing on the generation and crystallization of basaltic magma

In the system CaO-MgO-Al₂O₃-SiO₂, the tetrahedron CaMgSi₂O₆(di)-Mg₂SiO₄(fo)-SiO₂-CaAl₂ SiO₆(CaTs) forms a simplified basalt tetrahedron, and within this tetrahedron, the plane di-fo-CaAl₂Si₂O₈(an) separates simplified tholeiitic from alkalic basalts. Liquidus phase relations on this join have been studied at 1 atm and at 7, 10, 15, and 20 kbar. The temperature maximum on the 1 atm isobaric quaternary univariant line along which forsterite, diopside, anorthite, and liquid are in equilibrium lies to the SiO₂-rich side of the join di-fo-an. The isobaric quaternary invariant point at which forsterite, diopside, anorthite, spinel, and liquid are in equilibrium passes, with increasing pressure, from the silica-poor to the silica-rich side of the join di-fo-an, which causes the piercing points on this join to change from forsterite+diopside+anorthite+liquid and forsterite +spinel+anorthite+liquid below 5 kbar to forsterite +diopside+spinel+liquid and diopside +spinel+anorthite+liquid above 5 kbar. As pressure increases, the forsterite and anorthite fields contract and the diopside and corundum fields expand. The anorthite primary phase field disappears entirely from the join di-fo-an between 15 and 20 kbar. Below about 4 kbar, the join di-fo-an represents, in simplified form, a thermal divide between alkalic and tholeiitic basalts. From about 4 to at least 12 kbar, alkalic basalts can produce tholeiitic basalts by fractional crystallization, and at pressures above about 12 kbar, it is possible for alkalic basalt to be produced from oceanite by crystallization of both olivine and orthopyroxene. If alkalic basalts are primary melts from a lherzolite mantle, they must be produced at high pressures, probably greater than about 12 kbar.Department of Geosciences, University of Texas at Dallas Contribution No. 327. Hawaii Institute of Geophysics Contribution No. 814.     

Three-dimensional P- and S-wave velocity structures beneath the Japan Islands obtained by high-density seismic stations by seismic tomography

We construct fine-scale 3D P- and S-wave velocity structures of the crust and upper mantle beneath the whole Japan Islands with a unified resolution, where the Pacific (PAC) and Philippine Sea (PHS) plates subduct beneath the Eurasian (EUR) plate. We can detect the low-velocity (low-V) oceanic crust of the PAC and PHS plates at their uppermost part beneath almost all the Japan Islands. The depth limit of the imaged oceanic crust varies with the regions. High-V_P/V_S zones are widely distributed in the lower crust especially beneath the volcanic front, and the high strain rate zones are located at the edge of the extremely high-V_P/V_S zone; however, V_P/V_S at the top of the mantle wedge is not so high. Beneath northern Japan, we can image the high-V subducting PAC plate using the tomographic method without any assumption of velocity discontinuities. We also imaged the heterogeneous structure in the PAC plate, such as the low-V zone considered as the old seamount or the highly seismic zone within the double seismic zone where the seismic fault ruptured by the earthquake connects the upper and lower layer of the double seismic zone. Beneath central Japan, thrust-type small repeating earthquakes occur at the boundary between the EUR and PHS plates and are located at the upper part of the low-V layer that is considered to be the oceanic crust of the PHS plate. In addition to the low-V oceanic crust, the subducting high-V PAC plate is clearly imaged to depths of approximately 250 km and the subducting high-V PHS zone to depths of approximately 180 km is considered to be the PHS plate. Beneath southwestern Japan, the iso-depth lines of the Moho discontinuity in the PHS plate derived by the receiver function method divide the upper low-V layer and lower high-V layer of our model at depths of 30–50 km. Beneath Kyushu, the steeply subducting PHS plate is clearly imaged to depths of approximately 250 km with high velocities. The high-V_P/V_S zone is considered as the lower crust of the EUR plate or the oceanic crust of the PHS plate at depths of 25–35 km and the partially serpentinized mantle wedge of the EUR plate at depths of 30–45 km beneath southwestern Japan. The deep low-frequency nonvolcanic tremors occur at all parts of the high-V_P/V_S zone—within the zone, the seaward side, and the landward side where the PHS plate encounters the mantle wedge of the EUR plate. We prove that we can objectively obtain the fine-scale 3D structure with simple constraints such as only 1D initial velocity model with no velocity discontinuity.     

An Experimental Characterization of Shear Wave Velocity (V<Subscript>s</Subscript>) in Clean and Hydrocarbon-Contaminated Sand

The characteristics of hydrocarbon-contaminated soils have been among major concerns of geotechnical engineers due to its significant frequency of event and also its influential consequences on our surroundings from various environmental and engineering viewpoints. Heretofore, the effects of diverse kinds of hydrocarbon contaminants on majority of geotechnical properties of fine- and coarse-grained soils such as grain size, hydraulic conductivity, plasticity, compressibility, internal friction, cohesion, and shear strength have been investigated. However, there has not been a concentrated research study examining shear wave velocity (\({\text{V}}_{\text{s}}\)) of hydrocarbon-contaminated soils as an important geotechnical property of soil due to this fact that, in small/very small strain levels, the maximum shear modulus of soils (\({\text{G}}_{ \hbox{max} }\)) can be determined using shear wave velocity (\({\text{G}}_{ \hbox{max} } =\uprho{\text{V}}_{\text{s}}^{2}\)). This paper aims to investigate effects of hydrocarbon contamination on shear wave velocity of sandy soils by comparing shear wave velocities in identically prepared clean and contaminated samples. To this aim, an Iranian light crude oil, a standard type of silica sand (Ottawa sand), and a bender element apparatus were used to minutely measure shear wave velocity of clean and crude oil contaminated sand samples. Moreover, dry and quasi-moist tamping methods were employed in order to provide comparable clean and contaminated specimens (containing 4, 6, 8, 10, and 12 wt% of crude oil), respectively. Firstly, a comprehensive bender element (BE) and resonant column tests were conducted on the identically prepared clean sand samples at various amounts of frequency (2–20 kHz) and under various confining pressure (50–500 kPa) to find the best methods of accurately determining shear wave travel time in BE tests. Thereafter, BE tests were conducted to examine shear wave velocity in contaminated specimens. Based on the results, it was found that there was a critical value for crude oil content with the maximum shear wave velocity so that shear wave velocity of 4 wt% contaminated sand (V_s-4 wt%) was about 1.2 times higher than clean one (V_s-clean), and contrastingly adding further crude oil up to 6 wt% made a significant reduction in value of shear wave velocity to some extent that V_s-6 wt% was slightly lower than V_s-clean (V_s-6 wt% = 0.95–0.97V_s-clean). Moreover, adding more contaminant (8–12 wt%) into sand had negligible influences on shear wave velocity. In this paper, the effects of crude oil contamination on sand microstructure were also evaluated using scanning electron microscopy.