Use of mode-converted waves in marine seismic data to investigate the lithology of the sub-bottom sediments in Isfjorden,Svalbard

The compression wavefield is efficiently converted to shear-wave energy at post-critical angles in areas of high impedance contrast at the sea floor. We have analysed mode-converted shear waves in a data set acquired with a hybrid marine/land geometry in Isfjorden, Svalbard. Through a kinematic 2D ray-tracing modellingV _p/V_s ratios for part of the uppermost 5km of the crust are obtained. Low values (V _p /V _s =1.65) are tentatively associated with the section of Devonian sandstones which appears to attain a minimum thickness of 1.5km below 3 km depth about 10km west of Kapp Thorden.     

The elastic thickness of the lithosphere in the Pacific Ocean

In this study, we present determinations of the effective elastic thicknessT_e of the oceanic lithosphere along Pacific chains or archipelagoes.T_e is determined by computing the deflection of a continuous elastic plate under the load of volcanoes, and constrained by geoid heights provided by SEASAT. In the South Central Pacific, estimates of 14 km for the Marquesas and 6 km or less for the Pitcairn-Mururoa-Gloucester chain are in good agreement with a previous work in this region (Cook-Austral and Society chains). Around the Line Islands chain, SEASAT data reveal that the bathymetry is poorly known, preventing fine analysis. Meanwhile,T_e looks globally very low ( 6 km), except for three volcanoes but these results may be unreliable. The Easter chain features lowT_e values ( 6 km), with no noticeable variation along the chain. Higher values are found for a Samoan island, Manuae (24 km), and along the Hawaiian-Emperor seamounts chain (from 32 km at the eastern end of the chain to 21.5 km for the Hawaiian volcanoes, and from 25.5 to 15 km for the Emperor seamounts). The large number ofT_e estimates obtained in this study points out a noticeable difference between North and South Pacific results. Those from the North Pacific agree with the general trend (increase with the square root of age plate at loading time), while those from the South Central Pacific are much lower, according to their plate age. These lowT_e results from the South Pacific are only partly explained by taking account of thermal perturbations using the rejuvenation model. Therefore, these results then point out a regional difference in oceanic lithosphere.     

A crust-mantle model for the NORSAR area

Summary Elastic waves from explosions were recorded at NORSAR and at a number of field stations, and the data were used for determining a crust-mantle model under the array. The number of explosions was eleven distributed on seven shot points. The total number of recording points was fifty-one, and the interpretation was based on 350 individual records.The velocities obtained for the crustal phases were 6.2, 6.6 and 8.2 km/sec for theP _g ,P _g andP _n waves respectively. A deep crustal phase with a velocity of about 7.4 km/sec was observed. The mean depths to the discontinuities within the crust were determined to be 17 and 26 km. The depth to Moho varied greatly across the array from 31.5 km in the central part to 38 km under the C-ring. The maximum dip observed for the Moho was 12^o.Contribution No. 57 to Norwegian Geotraverse Project.     

云南地区地壳速度模型研究及应用

选取2009~2014年发生在云南地区、每个地震均在10个以上台站有记录的7412个地震数据,作走时曲线。同时为提高精度,重点对其中每个地震均在80%以上台站有记录的、ML≥3.0的83个地震数据,再作线性分析、折合走时曲线和区间稳定性分析,结合前人研究成果得到了研究区的初始地壳速度模型。选取2010~2014年云南省内M≥3.0的200次地震,采用Hyposat批处理方法迭代初始速度模型,以及对S波作分层速度拟合,得到云南地区的地壳速度模型,即2015云南模型:v_(P1)=6.01km/s,v_(P2)=6.60km/s,v_(Pn)=7.89km/s,H_1=20km,H_2=21km,v_(S1)=3.52km/s,v_(S2)=3.86km/s,v_(sn)=4.43km/s。基于新模型的地震重定位分析表明,云南地区地震事件大多发生于20km内的上地壳;对2011年3月10日盈江M5.8和2014年8月3日鲁甸M6.5典型地震进行重定位,得出震源深度分别与精定位结果和震中强震台震源距接近,表明新的一维速度模型能更好地反映研究区平均速度结构。     

Fine S wave velocity structure beneath Iwate volcano, northeastern Japan, as derived from receiver functions and travel times

A genetic algorithm inversion of receiver functions derived from a dense seismic network around Iwate volcano, northeastern Japan, provides the fine S wave velocity structure of the crust and uppermost mantle. Since receiver functions are insensitive to an absolute velocity, travel times of P and S waves propagating vertically from earthquakes in the subducting slab beneath the volcano are involved in the inversion. The distribution of velocity perturbations in relation to the hypocenters of the low-frequency (LF) earthquakes helps our understanding of deep magmatism beneath Iwate volcano. A high-velocity region (dV_S/V_S=10%) exists around the volcano at depths of 2–15 km, with the bottom depth decreasing to 11 km beneath the volcano’s summit. Just beneath the thinning high-velocity region, a low-velocity region (dV_S/V_S=−10%) exists at depths of 11–20 km. Intermediate-depth LF (ILF) events are distributed vertically in the high-velocity region down to the top of the low-velocity region. This distribution suggests that a magma reservoir situated in the low-velocity region supplies magma to a narrow conduit that is detectable by the hypocenters of LF earthquakes. Another broad low-velocity region (dV_S/V_S=−5 to −10%) occurs at depths of 17–35 km. Additional clusters of deep LF (DLF) events exist at depths of 32–37 km in the broad low-velocity zone. The DLF and ILF events are the manifestations of magma movement near the Moho discontinuity and in the conduit just beneath the volcano, respectively.     

Features of deep crustal structure and the onshore-offshore transition at the Iberian flank of the Valencia Trough (Western Mediterranean)

First results are presented of a recent onshore seismic survey complementary to the Valsis-2 Cruise, which consisted of ESP, COP and CDP marine seismic profiles across the Valencia Trough (Western Mediterranean).The marine energy source used was an airgun array of 5800 cubic inch recorded at 2 land stations on the western flank of the Valencia Trough, at distances between 10–120 km.The experiment has resulted in an extended sampling of the deep crustal structure of the eastern Mediterranean flank of the Iberian peninsula, as well as the offshore-onshore transition.Three transverse NW-SE profiles have been interpreted. Local thinning of the sedimentary cover has been determined towards the centre of the basin which, together with the shallow high velocities observed on the southern profile, could be related to volcanic episodes.A seismic continental basement has been found at depths between 3 and 5 km. A thin lower crust (3–5 km) with velocities around 6.8 km/s has been identified in the northern part of the basin. Alternative crustal models considered for the 3 profiles have been tested, not only from arrival times but also from relative amplitude distributions. A first-order Moho discontinuity fits the data best. The welldefined Moho boundary results in energetic P_MP reflections, and a clear updoming is observed towards the interior of the basin, from depths about 20–21 km inshore of Barcelona to 15–17 km depths 60 km offshore. An anomalous upper mantle with low P_n velocities of about 7.7 km/s is confirmed in most of the sampled areas.     

新疆地区一维地壳速度模型研究

利用新疆区域地震台网观测到的2009年1月—2014年7月Pn、Sn、Pg和Sg震相数据,综合使用线性拟合、折合走时、PTD定深方法和HypoSAT定位方法确定该地区Pg、Pb和Pn的平均传播速度(v_(Pg),v_(Pb),v_(Pn))、康拉德界面和莫霍面的深度(H_(conr)和H_(moho))范围,以速度和深度步长分别为0.1km/s、1km精度迭代计算样本数据,通过对比分析计算结果与全国地震统一编目和3400模型下样本数据的定位结果后,确定RMS平均值最小的一维速度模型。在新模型中v_(Pg)、v_(Pb)和v_(Pn)分别为6.10km/s、6.70km/s和8.20km/s,H_(conr)和H_(moho)分别为26km、54km。通过检验对比,认为本文获得的新模型优于新疆地区现有的3400模型。     

Tidal circulation and energy dissipation in a shallow, sinuous estuary

The tidal dynamics in a pristine, mesotidal (>2 m range), marsh-dominated estuary are examined using moored and moving vessel field observations. Analysis focuses on the structure of the M ₂ tide that accounts for approximately 80% of the observed tidal energy, and indicates a transition in character from a near standing wave on the continental shelf to a more progressive wave within the estuary. A slight maximum in water level (WL) occurs in the estuary 10–20 km from the mouth. M ₂ WL amplitude decreases at 0.015 m/km landward of this point, implying head of tide approximately 75 km from the mouth. In contrast, tidal currents in the main channel 25 km inland are twice those at the estuary mouth. Analysis suggests the tidal character is consistent with a strongly convergent estuarine geometry controlling the tidal response in the estuary. First harmonic (M ₄) current amplitude follows the M ₂ WL distribution, peaking at mid-estuary, whereas M ₄ WL is greatest farther inland. The major axis current amplitude is strongly influenced by local bathymetry and topography. On most bends a momentum core shifts from the inside to outside of the bend moving seaward, similar to that seen in unidirectional river flow but with point bars shifted seaward of the bends. Dissipation rate estimates, based on changes in energy flux, are 0.18–1.65 W m⁻² or 40–175 μW kg^–1. A strong (0.1 m/s), depth-averaged residual flow is produced at the bends, which resembles flow around headlands, forming counter-rotating eddies that meet at the apex of the bends. A large sub-basin in the estuary exhibits remarkably different tidal characteristics and may be resonant at a harmonic of the M ₂ tide.     

Receiver function study in northern Sumatra and the Malaysian peninsula

In this receiver function study, we investigate the structure of the crust beneath six seismic broadband stations close to the Sunda Arc formed by subduction of the Indo-Australian under the Sunda plate. We apply three different methods to analyse receiver functions at single stations. A recently developed algorithm determines absolute shear-wave velocities from observed frequency-dependent apparent incidence angles of P waves. Using waveform inversion of receiver functions and a modified Zhu and Kanamori algorithm, properties of discontinuities such as depth, velocity contrast, and sharpness are determined. The combination of the methods leads to robust results. The approach is validated by synthetic tests. Stations located on Malaysia show high-shear-wave velocities (V _S) near the surface in the range of 3.4–3.6 km s^− 1 attributed to crystalline rocks and 3.6–4.0 km s^− 1 in the lower crust. Upper and lower crust are clearly separated, the Moho is found at normal depths of 30–34 km where it forms a sharp discontinuity at station KUM or a gradient at stations IPM and KOM. For stations close to the subduction zone (BSI, GSI and PSI) complexity within the crust is high. Near the surface low V _S of 2.6–2.9 km s^− 1 indicate sediment layers. High V _S of 4.2 km s^− 1 are found at depth greater than 6 and 2 km at BSI and PSI, respectively. There, the Moho is located at 37 and 40 km depth. At station GSI, situated closest to the trench, the subducting slab is imaged as a north-east dipping structure separated from the sediment layer by a 10 km wide gradient in V _S between 10 and 20 km depth. Within the subducting slab V _S ≈ 4.7 km s^− 1. At station BSI, the subducting slab is found at depth between 90 and 110 km dipping 20° ± 8° in approximately N 60° E. A velocity increase in similar depth is indicated at station PSI, however no evidence for a dipping layer is found.     

Vertical distribution features of atmospheric water vapour in the Po valley area

Two sets of radiosounding measurements, taken at different hours from some stations in the Po Valley area, were examined in order to obtain the values of precipitable waterw and determine the shapes of the vertical distribution curves of absolute humidity. From these curves, we calculated the bestfit values of scale heightsH ₁ andH ₂ in the 0÷5 km and 5÷10 km altitude intervals, respectively. The analysis of the radiosounding data clearly shows that parametersw andH ₁ are closely related to the evolutionary features of the meteorological conditions on synoptic scale and are influenced, to a lesser extent, by the diurnal variations in the atmospheric ground layer. Seasonal average curves of temperature and absolute humidity, as functions of altitude, were also defined from the radiosounding measurements taken at various hours of the day. Moreover, interpolation methods in time (applied to a 12-hour range) and in space (range of about 300 km) were proposed for determining parametersw andH ₁ from the radiosounding measurements taken at different hours from the same station or from different stations at the same hour. Reliability tests, made by comparing the same station or from different stations at the same hour. Reliability tests, made by comparing the values given by the interpolation methods with those directly obtained from the radiosounding measurements, show that these evaluations ofw andH ₁ are affected by standard errors of estimate, which are comparable to the errors usually made in analyzing the radiosounding data.     

Interpretation of S-wave data from Tai’an-Xinzhou DSS profile and its relationship with Xingtai earthquakes

On the basis of S wave information from Tai’an-Xinzhou DSS profile and with reference to the results from P-wave interpretation, the 2-D structures, including S-wave velocity V _s, ratio γ between V _p and V _s; and Poisson’ s ratio σ, are calculated; the structural configuration of the profile is presented and the relevant inferences are drawn from the above results. Upwarping mantle districts (V _s≈4.30 km/s) and sloping mantle districts (V _s≈4.50 km/s) of the profile with velocity difference about −4% at the top of upper mantle are divided according to the differences of V _s, γ and σ in different media and structures, also with reference to the information of their neighbouring regions; the existence of Niujiaqiao-Dongwang high-angle ultra-crustal fault zone is reaffirmed; the properties of low and high velocity blocks (zones) including the crust-mantle transitionalzone and the boudary indicators of North China rift valley are discussed. A comprehensive study is conducted on the relation of the interpretation results with earthquakes. It is concluded that the mantle upwarps, thermal material upwells through the high-angle fault, the primary hypocenter was located at the crust-mantle juncture 30.0∼33.0 km deep, and additional stress excited the M _S=6.8 and M _S=7.2 earthquakes at specific locations around 9.0 km below Niujiaqiao-Dongwang, the earthquakes took place around the high-angle ultra-crustal fault and centered in the brittle media and rock strata with low γ and low σ values. This subject is part of the 85-907-02 key project during the “8th Five-Year Plan” from the State Science and Technology Commission.     

Moho depth of the European Plate from teleseismic receiver functions

Crustal structure and the Moho depth are exceptionally well known beneath Europe. The first digital, high-resolution map of the Moho depth for the whole European Plate was compiled in 2007 and recently published in Geophysical Journal International. In the past few years, considerable developments have taken place in the receiver function techniques. Different receiver function techniques provide new, independent information, in particular on the S-wave velocity distribution in the crust and on the Moho depth. This gives an opportunity to compare the Moho depth from the Moho depth map of the European Plate (H _MM) and the Moho depth from receiver function studies (H _RF). Herein, we also compile and analyze the uncertainty of the crustal thickness determinations data obtained with receiver function analysis. The uncertainty is found to be ±2 km for 20-km-thick crust and about ±4 km for 60-km-thick crust. Comparison of the Moho depths shows an approximately linear trend between H _RF and H _MM. For the Moho depth of 30–40 km, the values are approximately equal, while for thin crust, H _RF is about 5 km shallower than H _MM, and for thick crust, it is about 5 km deeper than H _MM. Possible reasons for this, the observed discrepancy between the Moho depths H_MM and H_RF, are discussed.     

Attenuation of Coda Waves at the Changbaishan Tianchi Volcanic Area in Northeast China

21 earthquakes recorded by a temporary seismic network in the Changbaishan Tianchi volcanic area in Northeast China operated during the summer of 2002 and 2003 were analyzed to estimate the S coda attenuation. The attenuation quality factor Q_c was estimated using the single scattering attenuation model of Sato (1977) in the frequency band from 4 to 24 Hz. All the events studied in this paper occurred at depths from 2 to 6 km with ML of 1.4–2.8. The epicentral distances are less than 25 km. For all events which occurred near the Tianchi Lake (caldera), the Q_c patterns obtained at the stations near the lake are similar, and the Q_c values are relatively small. At the stations located about 15 km east of the Tianchi Lake, however, the average Q_c is significantly higher. For an event which occurred 25km from the lake to the west, Q_c patterns derived at the stations near the lake are quite similar to the above mentioned Q_c for stations located in the east. Further study shows that Q_c value in the north and central areas of the volcano is relatively lower than that in the surrounding area. Compared to other volcanic areas in the world, the average Q_c of the Changbaishan Tianchi volcanic area is obviously lower. The deep seismic sounding and teleseismic receiver function studies indicated more than one lower velocity layer in the crust. The MT studies suggested the presence of high conductive bodies beneath the area. We interpret the strong attenuation of coda waves near the Changbaishan Tianchi volcano as being possibly related to high temperature medium caused by shallow magma chambers.     

利用CAP方法和瑞利面波振幅谱联合反演宁强5.3级地震震源深度

2018年9月12日19时6分,陕西省汉中市宁强县发生5.3级地震,不同机构给出的震源深度结果相差较大。为进一步确定宁强5.3级地震的震源深度,基于区域速度模型,首先利用CAP方法反演得到该地震的震源机制解,然后采用瑞利面波振幅谱和CAP深度误差函数联合反演,进一步测定了此次地震的矩心深度。结果显示:CAP方法得到的陕西宁强5.3级地震矩心深度约为12km,瑞利面波振幅谱测定的矩心深度为13km,结合引入的误差函数联合反演,最终确定陕西宁强5.3级地震的矩心深度为13km左右,表明此次地震仍属于发生于上地壳的地震。     

Inversion of seismic and gravity data for the composition and core sizes of the Moon

We model the internal structure of the Moon, initially homogeneous and later differentiated due to partial melting. The chemical composition and the internal structure of the Moon are retrieved by the Monte-Carlo inversion of the gravity (the mass and the moment of inertia), seismic (compressional and shear velocities), and petrological (balance equations) data. For the computation of phase equilibrium relations and physical properties, we have used a method of minimization of the Gibbs free energy combined with a Mie-Gr@uneisen equation of state within the CaO-FeO-MgO-Al₂O₃-SiO₂ system. The lunar models with a different degree of constraints on the solution are considered. For all models, the geophysically and geochemically permissible ranges of seismic velocities and concentrations in three mantle zones and the sizes of Fe-10%S core are estimated. The lunar mantle is chemically stratified; different mantle zones, where orthopyroxene is the dominant phase, have different concentrations of FeO, Al₂O₃, and CaO. The silicate portion of the Moon (crust + mantle) may contain 3.5–5.5% Al₂O₃ and 10.5–12.5% FeO. The chemical boundary between the middle and the lower mantle lies at a depth of 620–750 km. The lunar models with and without a chemical boundary at a depth of 250–300 km are both possible. The main parameters of the crust, the mantle, and the core of the Moon are estimated. At the depths of the lower mantle, the P and S velocities range from 7.88 to 8.10 km/s and from 4.40 to 4.55 km/s, respectively. The radius of a Fe-10%S core is 340 ± 30 km.     

The effective elastic lithosphere under the Cook-Austral and Society islands

We have determined the elastic thicknessT_e of the oceanic lithosphere along two volcanic chains of the South Central Pacific: Cook-Austral and Society islands. We used a three-dimensional spatial method to model the lithospheric flexure assuming a continuous elastic plate. The model was constrained by geoid height data from the SEASAT satellite.Along the Cook-Austral chain the elastic thickness increases westward, from 2–4 km at McDonald hot spot to 14 km at Rarotonga. At McDonald seamount, however, the data are better explained by a local compensation model. The observed trend shows an increase ofT_e with age of plate at loading time. However, the elastic layer under the Cook-Austral appears systematically thinner by several kilometers than expected for “normal” seafloor, suggesting that substantial thermal thinning has taken place in this region. Considering the apparent thermal age of the plate instead of crustal age improves noticeably the results. Along the Society chainT_e varies from 20 km under Tahiti to 13 km under Maupiti which is located 500 km westward. When plotting together the Society and Cook-AustralT_e results versus age of load, we notice that within the first five million years after loading,T_e decreases significantly while tending rapidly to an equilibrium value. This may be interpreted as the effect of initial stress relaxation which occurs just after loading inside the lower lithosphere and suggests that the presently measured elastic thickness under the very young Tahiti load ( 0.8 Ma) is not yet the equilibrium thickness.     

Experimental evidence for flexibility of a building foundation supported by concrete friction piles

This article explores the possibility to measure deformations of building foundations from measurements of ambient noise and strong motion recordings. The case under study is a seven-storey hotel building in Van Nuys, California. It has been instrumented by strong motion accelerographs, and has recorded several earthquakes, including the 1971 San Fernando (M_L=6.6, R=22 km), 1987 Whittier–Narrows (M_L=5.9, R=41 km), 1992 Landers (M_L=7.5, R=186 km), 1992 Big Bear (M_L=6.5, R=149 km), and 1994 Northridge (M_L=6.4, R=1.5 km) earthquake and its aftershocks (20 March: M_L=5.2, R=1.2 km; 6 December, 1994: M_L=4.3, R=11 km). It suffered minor structural damage in 1971 earthquake and extensive damage in 1994. Two detailed ambient vibration tests were performed following the Northridge earthquake, one before and the other one after the 20 March aftershock. These included measurements at a grid of points on the ground floor and in the parking lot surrounding the building, presented and analyzed in this article. The analysis shows that the foundation system, consisting of grade beams on friction piles, does not act as a “rigid body” but deforms during the passage of microtremor and therefore earthquake waves. For this geometrically and by design essentially symmetric building, the center of stiffness of the foundation system appears to have large eccentricity (this is seen both from the microtremor measurements and from the earthquake recordings). This eccentricity may have contributed to strong coupling of transverse and torsional responses, and to larger than expected torsional response, contributing to damage during the 1994 Northridge, earthquake.     

Auroral E-region electron density gradients measured with EISCAT

In the theory of E-region plasma instabilities, the ambient electric field and electron density gradient are both included in the same dispersion relation as the key parameters that provide the energy for the generation and growth of electrostatic plasma waves. While there exist numerous measurements of ionospheric electric fields, there are very few measurements and limited knowledge about the ambient electron density gradients, N_e, in the E-region plasma. In this work, we took advantage of the EISCAT CP1 data base and studied statistically the vertical electron density gradient length, L_z = N_e/(dN_e/d_z), at auroral E-region heights during both eastward and westward electrojet conditions and different ambient electric field levels. Overall, the prevailing electron density gradients, with L_z ranging from 4 to 7 km, are found to be located below 100 km, but to move steadily up in altitude as the electric field level increases. The steepest density gradients, with L_z possibly less than 3 km, occur near 110 km mostly in the eastward electrojet during times of strong electric fields. The results and their implications are examined and discussed in the frame of the linear gradient drift instability theory. Finally, it would be interesting to test the implications of the present results with a vertical radar interferometer.     

Geothermal profile and crust-mantle transition beneath east-central Queensland: Volcanology, xenolith petrology and seismic data

Geothermobarometry of garnet granulite and garnet websterite xenoliths in basalts from numerous localities in east-central Queensland gives P-T points that fall along the geotherm previously defined for southeastern Australia. This elevated geotherm is ascribed to the advective transport of heat by Tertiary-Recent magmas ponded at the crust-mantle boundary. The lower crust in this region consists dominantly of mafic granulites, representing frozen basaltic melts and cumulates. Spinel lherzolite becomes a dominant rock type at depths of ca. 30 km, and persists, interlayered with pyroxenites, to depths of ca. 55 km. Seismic reflection profiles show a “layered lower crust” between depths of 20 and 36 km depth. The lithologically defined crust-mantle boundary lies within this zone, at least 6 km above the seismically defined Moho. This interpretation is consistent with the observed velocity (V_p) gradient downward through the layered zone. The constructed geotherm implies that the bottom of the lithosphere beneath eastern Australia is shallower than ca. 100 km. This makes it unlikely that the diamonds of eastern Australia are derived from local intrusions, unless these are > 200 Ma old.     

Earthquake ground motion in Mexico City: An analysis of data recorded at Roma array

We analysed in detail three earthquakes recorded in a small-aperture accelerometric array in Mexico City, using the correlation of the records as a function of time along the accelerogram and frequency. Ground response is strongly conditioned by the fundamental period of the soft soils at the site of the array (T₀). Energy at periods longer than 2T₀ is guided by the crustal structure (with a thickness of 45 km). The wave field at periods between T₀ and 2T₀ also consists of surface waves but guided by the upper 2–3 km of volcanic sediments in central Mexico. For periods smaller than T₀, ground motion is uncorrelated among the stations. Our results indicate that seismic response of Mexico City, including its very long duration, results from deeply guided surface waves (between 2 and 45 km depth) interacting with the very local response of the soft surficial clay layer.