The Chelyabinsk meteoroid: A seismologist’s view

This work determines the parameters of the meteoroid that exploded in Chelyabinsk oblast on February 15, 2013. Using video observation and automobile registration cameras, we identified the trajectory and velocity of the meteoric body from the beginning of its glow to the destruction. The precise explosion time (03:20:34UTC) was determined on the basis of records of the broadband seismic stations and arrivals of the surface wave. The estimated explosion energy with assumption of a point source similar to a high thermonuclear explosion was 100 kt in TNT equivalent based on the surface wave amplitude. The archive of records of seismic stations during the meteorite flight is composed.     

Crustal thickness and velocity structure beneath Singapore’s seismic network

We estimated the crustal thickness and velocity structure beneath the five stations comprising the Republic of Singapore’s seismic network. Our data set was composed of 697 teleseismic receiver functions and 7 months of broad-band data that was cross-correlated to produce inter-station Green’s functions. Surface wave group velocities were extracted from the Green’s functions to obtain dispersion data for a path from central Sumatra to Singapore in order to provide a complimentary data set to the receiver functions. Crustal thickness was estimated via an H − k stacking technique, and high-resolution 1D P-wave velocity profiles were generated beneath each station by jointly inverting receiver function stacks and the group velocity data using a linearised time-domain inversion scheme. Crustal thickness beneath four stations was found to be between 28.0 km and 32.0 km, while one station in the northeast of Singapore indicates 24.0 km thick crust. This implies a significant crustal thinning beneath Singapore over the lateral extent of 50.0 km. Inversion results exhibit several crustal features that are observable in the derived models at all five stations, indicating that they exist across Singapore as a whole. There appears to be an upper-crustal high-velocity zone beneath Singapore, underlain by a velocity inversion. Station NTU shows slower near-surface velocities than the other stations, consistent with its situation above the sedimentary Jurong formation. These results expand the available global velocity data set, as well as being useful for assessing the seismic hazard in Singapore.     

Distinct variations of crustal shear wave velocity structure and radial anisotropy beneath the North China Craton and tectonic implications

This study presents the crustal shear wave velocity structure and radial anisotropy along two linear seismic arrays across the North China Craton (NCC) from ambient noise tomography. About a half to one year long ambient noise data from 87 stations were used for obtaining the inter-station surface wave empirical Green's functions (EGFs) from cross-correlation. Rayleigh and Love dispersion curves within the period band 5–30 s were measured from the EGFs of the vertical and transverse components, respectively. These dispersion data were then used to determine the crustal shear wave velocity structure (V_SV and V_SH) and radial anisotropy (2(V_SH ? V_SV) / (V_SH + V_SV)) from point-wise linear inversion with constraints from receiver function analysis. Our results reveal substantial structural variations among different parts of the NCC. The Bohai Bay Basin in the eastern NCC is underlain by a thin crust (~ 30 km) with relatively low velocities (particularly V_SV) and large positive radial anisotropy in the middle to lower crust. Such a crustal structure is no longer of a cratonic type and may have resulted from the widespread tectonic extension and intensive magmatism in this region since late Mesozoic. Beneath the Ordos Basin in the western NCC, the crust is relatively thicker (≥ 40 km) and well stratified, and presents a large-scale low velocity zone in the middle to lower crust and overall weak radial anisotropy except for a localized lower crust anomaly. The overall structural features of this region resemble those of typical Precambrian shields, in agreement with the long-term stability of the region. The crustal structure under the Trans North China Orogen (TNCO, central NCC) is more complicated and characterized by smaller scale velocity variations, strong positive radial anisotropy in the middle crust and rapid change to weak-to-negative anisotropy in the lower crust. These features may reflect complex deformations and crust–mantle interactions, probably associated with tectonic extension and magmatic underplating during the Mesozoic to Cenozoic evolution of the region. Our structural images in combination with previous seismic, geological and geochemical observations suggest that the Phanerozoic lithospheric reactivation and destruction processes may have affected the crust (especially the middle and lower crust) of the eastern NCC, and the effect probably extended to the TNCO, but may have minor influence on the crust of the western part of the craton.     

Crustal and upper mantle S-wave velocity structures across the Taiwan Strait from ambient seismic noise and teleseismic Rayleigh wave analyses

Although orogeny tapers off in western Taiwan large and small earthquakes do occur in the Taiwan Strait, a region largely untouched in previous studies owing mostly to logistical reasons. But the overall crustal structure of this region is of particular interest as it may provide a hint of the proto-Taiwan before the orogeny.By combining time domain empirical Green’s function (TDEGF) from ambient seismic noise using station-pairs and traditional surface wave two-station method (TS) we are able to construct Rayleigh wave phase velocity dispersion curves between 5 and 120 s. Using Broadband Array in Taiwan for Seismology (BATS) stations in Taiwan and in and across the Strait we are able to derive average 1-D Vs structures in different parts of this region. The results show significant shear velocity differences in the upper 15 km crust as expected. In general, the highest Vs in the upper crust observed in the coastal area of Mainland China and the lowest Vs appears along the southwest offshore of the Taiwan Island; they differ by about 0.6–1.1 km/s. For different parts of the Strait, the upper crust Vs structures are lower in the middle by about 0.1–0.2 km/s relative to those in the northern and southern parts. The upper mantle Vs structure (Moho – 150 km) beneath the Taiwan Strait is about 0.1–0.3 km/s lower than the AK135 model. The overall crustal thickness is approximately 30 km, much thinner and less variable than under the Taiwan Island. The inversion of seismic velocity structures using shorter period band dispersion data in the sea areas with water depth deeper than 1000 m should take water layer into consideration except for the continental shelves.     

Seismic anisotropy in the wedge above the Philippine Sea slab beneath Kanto and southwest Japan derived from shear wave splitting

We conduct shear wave splitting measurements on waveform data from the Hi-net and the broadband F-net seismic stations in Kanto and SW Japan generated by shallow and intermediate-depth earthquakes occurring in the subducting Philippine Sea and Pacific slabs. We obtain 1115 shear wave splitting parameter pairs. The results are divided into those from the shallow (depth < 50 km) and the deep (depth > 50 km) events. The deep events beneath Kanto are further divided into PHS1 and PHS2 (upper and lower planes of the double seismic zone in the Philippine Sea slab, respectively), PAC1 and PAC2 (western and eastern Pacific slab, respectively) events. The results from the shallow events represent the crustal anisotropy, and their fast directions are more or less aligned in the σ_Hmax directions, implying that the anisotropy is produced by the alignment of the vertical cracks in the crust induced by the compressive stresses. In Kanto, Kii Peninsula and Kyushu regions, the results from the deep events suggest a contribution from the mantle wedge anisotropy. Events from all groups beneath Kanto show NW, NE and EW fast directions. This complex pattern seems to be produced by the corner flows induced by both the WNW PAC plate subduction and the oblique NNW PHS slab subduction with the associated olivine lattice-preferred orientations (LPOs), and the anisotropy frozen in the PHS slab. The deep events beneath Kii Peninsula show NE and NW fast directions and may be produced by the corner flow produced by the NNW PHS slab subduction with the associated olivine LPOs. The NE directions might also be produced by the segregated melts in the thin layers parallel to the PHS slab subduction. The deep events beneath N Kyushu show NNW fast directions, which may result from the southeastward flow in the upper mantle inferred from the stresses in the upper plate. Results from the deep events beneath middle-south Kyushu show dominantly E–W fast directions, in both the fore- and back-arcs. They may be produced by the corner flow of the westward PHS slab subduction with the olivine LPOs. Because the source regions with multiple fast directions are not resolved in this study, further detailed analyses of shear wave splitting are necessary for a better understanding of the stress state, the induced mantle flow, and the melt-segregation processes.     

Inverting ground motion from a seismometer array to obtain the vertical component of rotation: A test using data from explosions

To study factors that affect reliability of the derived vertical rotational ground motions (about vertical axis) from translational array waveforms, we have conducted a field test using explosion. The dataset was collected on March 4, 2008 from an active source experiment of the Taiwan Geodynamics Research project (TAIGER, 2008), in which 750 kg of explosives were set off about 500 m away from a specially designed seismic array in free field consisting of 8 tri-axial rotational sensors, 13 tri-axial accelerometers, and 12 six-channel, 24-bit dataloggers with GPS receivers. Our results show that the waveform fits were greatly improved when the explosion-induced crustal strains were less than 2 × 10⁻⁵ for our station configuration. As a result, our modeling for the large amplitude part of the mainshock was not entirely successful. However, we have successfully recovered the first order features of the later part of the vertical rotational ground motions using translational velocity waveforms in the 0.5–20 Hz bandwidth using the Spudich and Fletcher’s procedure (2009). Moreover, we used a modified Jaeger’s (1969) formula and derived similar rotational synthetics. We found that the best rotational rate waveform fits were derived using a station configuration of a dense, small-aperture, translational array combined with a few large-offset stations. A small-aperture array is needed to fulfill the uniform rotation assumption. However, inverting data from an array of small spatial dimension requires waveforms of high signal to noise ratio and high sampling rates. Surprisingly, large-offset translation stations provided additional constraints in time and space, helping to fit the low frequency waveforms. Using the inverted result, we forward predicted the vertical translational waveforms not used in the inversion and got excellent fits. This suggests that the inverted results are of good quality for predicting the complete translational ground motions, even though only horizontal components were used in the inversion. Array-derived or directly measured rotational ground motions can be used to predict translational ground motions in the near-by region. Lessons learned from this study could be helpful for future studies using translational ground motions to derive dynamic ground strains, tilts, and torsions.     

3D thermal structure of the continental lithosphere beneath China and adjacent regions

Based on the Crust2.0 model and the topography data of Chinese continent and its adjacent regions, a three-dimensional finite element model is constructed in terms of the spherical coordinate system. In our numerical model, the average annual ground temperature from 195 meteorological stations and temperature of upper mantle derived from the seismic velocities are adopted as the top and bottom boundary conditions, respectively. The observed thermal conductivity and heat production from P wave velocity based on empirical formula are employed in our numerical model as well. The comparison between the calculated and observed surface heat flow proved that our results are reliable. The temperature beneath the Precambrian cratons is lower than that of other areas for 100–300 °C also. The typical temperature rang at the Moho is estimated to be 800–1000 °C beneath the Tibetan plateau and 500–700 °C beneath the Precambrian cratons (such as Indian plate, Sichuan basin, South China, North China and Tarim), respectively. The thermal state in the eastern part of Sino-Korean craton at the depth deeper than 60 km indicates that it was destructed. The thermal structure in center of Tibetan plateau (especially beneath Qiangtang area) supports the proposed flow of lower crustal or upper mantle material to the east. Generally, the distribution of volcanoes in Chinese continent is consistent with the high temperature areas in the crust or upper mantle. There are many obvious thermal transition zones across the orogenic belts. The thermal transition zone between eastern and western parts in the crust of Chinese continent is consistent with the north–south seismic zone.     

P,S wave velocity model of the crust and upper most mantle of Albania region

This paper describes the one-dimensional (1D) velocity model computed by VELEST in the SEISAN seismic analysis system, inverting re-picked P-wave and S-wave arrival times recorded during 2002–2006 by the Albanian, Montenegro, Thessalonica and Macedonia seismic networks. The re-picked data yield P-wave and S-wave velocities proved to be more suitable compared to bulletin data for this detailed inversion study. Seismic phases recorded by the Albania seismic network and integrated with data from the Montenegro, Thessalonica and Macedonia networks are used to prepare the Albanian seismic bulletin. Earthquake hypocenters from the Albanian bulletins have also location errors that are negligible for civil protection purposes, large scale seismotectonic analyses and more accurate hypocentral determinations which are necessary for detailed seismotectonic and geodynamic studies.It was noted that the smoothness of the velocity variation increased with depth. A velocity of 5.5 km/s was calculated for the upper crust, 6.1 km/s was calculated for the middle crust and 6.9 km/s was computed for the lower crust. P wave velocity was 7.85 km/s at depth of 50 km and for the upper mantle it is 8.28 km/s. Using the improved velocity model, the earthquakes which occurred in Albania in the past 5 years were able to be relocated, achieving constrained hypocentral determinations for events in Albania. The interpretation of the 1 D velocity models infers interesting features of the deep structure of Albania. These results represent an important step towards more detailed seismotectonic analyses.     

Crustal structure beneath two seismic broadband stations revealed from teleseismic P-wave receiver function analysis in the Virunga volcanic area,Western Rift Valley of Africa

The shear velocity structure beneath the Virunga volcanic area was estimated by using an average solution in the time domain inversion of stacked teleseismic receiver functions provided by two seismic broadband stations KUNENE (KNN) and KIBUMBA (KBB). These two stations are 29 km apart and located at the eastern and western escarpment of the Western Rift Valley of Africa in the Virunga area, respectively. The velocity model was presented as P-wave velocity models. From these models, the crust mantle transition zone beneath the area sampled by KNN and KBB in the Virunga area was determined at depth from about 36 to 39 km and 30 to 41 km, respectively. A low velocity zone was observed below stations KNN and KBB at depths between 20–30 km and 18–28 km, respectively, and with average velocity 5.9 km/s and 6.0 km/s. This low velocity zone may probably related to a magma chamber or a melt-rich sill. The models show also high velocity material (6.8–7.4 km/s) lying beneath stations KNN and KBB at depths 3–20 km and 3–10 km, respectively, which is indicative of magma cumulates within the volcanic edifice. The result obtained in this study was applied to the determination of epicentres during the period prior to the 27 November 2006 Nyamuragira eruption. This eruption was preceded by a swarm of hybrid volcanic earthquakes with clear P-waves onset. Using the receiver function model was found to improve the location of events. The located events correlate well with the location of the eruptive site and data provided by the INSAR observations of surface deformation associated with eruption.     

Geological variation in S-wave velocity structures in Northern Taiwan and implications for seismic hazards based on ambient noise analysis

Ambient noise analysis in Northern Taiwan revealed obvious lateral variations related to major geological units. The empirical Green’s functions extracted from interstation ambient noise were regarded as Rayleigh waves, from which we analyzed the group velocities for period from 3 to 6 s. According to geological features, we divided Northern Taiwan into seven subregions, for which regionalized group velocities were derived by using the pure-path method. On average, the group velocities in mountain areas were higher than those in the plain areas. We subsequently inverted the S-wave velocity structure for each subregion down to 6 km in depth. Following the analysis, we proposed the first models of geology-dependent shallow S-wave structures in Northern Taiwan. Overall, the velocity increased substantially from west to east; specifically, the mountain areas, composed of metamorphic rocks, exhibited higher velocities than did the coastal plain and basin, which consist of soft sediment. At a shallow depth, the Western Coastal Plain, Taipei Basin, and Ilan Plain displayed a larger velocity gradient than did other regions. At the top 3 km of the model, the average velocity gradient was 0.39 km/s per km for the Western Coastal Plain and 0.15 km/s per km for the Central Range. These S-wave velocity models with large velocity gradients caused the seismic waves to become trapped easily in strata and, thus, the ground motion was amplified. The regionalized S-wave velocity models derived from ambient noises can provide useful information regarding seismic wave propagation and for assessing seismic hazards in Northern Taiwan.     

Crustal shear-wave velocity structure beneath northeast India from teleseismic receiver function analysis

We investigated the seismic shear-wave velocity structure of the crust beneath nine broadband seismological stations of the Shillong–Mikir plateau and its adjoining region using teleseismic P-wave receiver function analysis. The inverted shear wave velocity models show ∼34–38 km thick crust beneath the Shillong Plateau which increases to ∼37–38 km beneath the Brahmaputra valley and ∼46–48 km beneath the Himalayan foredeep region. The gradual increase of crustal thickness from the Shillong Plateau to Himalayan foredeep region is consistent with the underthrusting of Indian Plate beyond the surface collision boundary. A strong azimuthal variation is observed beneath SHL station. The modeling of receiver functions of teleseismic earthquakes arriving the SHL station from NE backazimuth (BAZ) shows a high velocity zone within depth range 2–8 km along with a low velocity zone within ∼8–13 km. In contrast, inversion of receiver functions from SE BAZ shows high velocity zone in the upper crust within depth range ∼10–18 km and low velocity zone within ∼18–36 km. The critical examination of ray piercing points at the depth of Moho shows that the rays from SE BAZ pierce mostly the southeast part of the plateau near Dauki fault zone. This observation suggests the effect of underthrusting Bengal sediments and the underlying oceanic crust in the south of the plateau facilitated by the EW-NE striking Dauki fault dipping 30⁰ toward northwest.     

Crustal structure beneath Kottamiya broadband station,northern Egypt,from analysis of teleseismic receiver functions

We analyzed a total of 206 receiver functions beneath Kottamiya broadband station in northern Egypt to study the crustal structure and any azimuthal variations in the crustal thickness. The computed receiver functions are subdivided according to their azimuth into eight subgroups and analyzed separately using a genetic algorithm. The genetic algorithm is more appropriate than conventional linearized inversion schemes in regions where there is little a priori information about local crustal structures such as northern Egypt because it does not strongly depend on an initial model. The study region is located on the unstable shelf of Egypt in the northeastern corner of Africa. Little information about the deep structure of the crust beneath this region is available. For this reason, we have adopted the genetic algorithm to seismic waveform data recorded by Kottamiya broadband station. The crustal thickness varies slightly from 32 to 34 km with an average of 32.25 km, which is consistent with previous studies in the region. The crustal thickness shows a tendency of decrease toward the east and northeast being consistent with the general tectonic setting of the region including the opening of the Red Sea in the Tertiary times. Nonetheless, more teleseismic receiver functions from earthquakes recorded at denser seismic stations in northern Egypt and the southeastern Mediterranean combined with surface wave dispersion data as well as other geophysical investigations are necessary for more detailed imaging of the crustal structure which will deepen our understanding of the current tectonic and seismic activities of the region.     

Chemical composition of crustal xenoliths from southwestern Syria: Characterization of the upper part of the lower crust beneath the Arabian plate

The chemical bulk rock composition of 37 xenoliths, brought from depths of 25–30 km to the surface by penetrating Cenozoic alkali basaltic magma, from the Shamah Harrat, southwestern Syria, was determined by XRF spectroscopy. The geochemical character of these xenoliths points to original marls and within-plate igneous rocks. To obtain the mean chemical composition of the corresponding upper portion of the lower crust, the compositions of the 37 xenoliths were averaged and a leucogranitic and upper crustal component was added to account for assimilation by the Cenozoic magmas. This mean is more basic (SiO₂—50.5 wt%) and richer in HFSE, LREE, and LILE compared to compositions of the lower crust given by Taylor and McLennan [1985. The Continental Crust: Its Composition and Evolution. Blackwell, Oxford, 312pp.] and Rudnick and Gao [2005. Composition of the continental crust. In: Rudnick, R.L. (Ed.), The Crust. Treatise on Geochemistry, vol. 3. Elsevier, Amsterdam, pp. 1–64]. Calculations of the seismic compressional-wave velocity from our compositional mean, using the PERPLE_X computer software, yielded values around 6.85 km/s, which are in accordance with reported seismic studies for the corresponding depth levels (6.7–7.1 km/s).     

An attempt to detect temporal variations of crustal structure in the source area of the 2006 Wen-An earthquake in North China

Lei et al. (2008) revealed a low-velocity (low-V) anomaly in the lower crust under the source area of the 4 July 2006 Wen-An earthquake (M 5.1) in North China. In this work we tried to investigate the temporal variations of the crustal structure by using a number of P and PmP (Moho reflected) wave arrival times recorded by 107 digital seismic stations from earthquakes that occurred separately in 2002, 2003, 2004, and 2005–2006 to determine P-wave velocity structures in and around the source area of the Wen-An earthquake in different periods. Our results show that tomographic images inferred from the data sets in different years are all dominated by a low-V anomaly in the lower crust under the Wen-An source area. However, there exist some differences in the P-wave velocity image in the Wen-An source area inferred from the P and P + PmP data sets, suggesting that the PmP data have improved the tomographic images in the middle and lower crust, but the results from the 2005–2006 P and P + PmP data sets all show a relatively lager increase of the low-V anomaly under the Wen-An source area in the amplitude and extent as compared with those from the 2003 and 2004 P and P + PmP data sets. Incorporating the previous results, if this low-V anomaly may indicate the existence of fluids, then our results suggest that the occurrence of the Wen-An earthquake is not only related to the long-term influence of fluids that decrease the effective normal stress on the fault plane, but also closely associated with the drastic increase of such influence. However, this study is just an experimental work and the results are still preliminary because the resolution scale of the present tomographic model is much larger than the Wen-An source area and our extensive tests show that different samplings of seismic rays from different data sets have affected the details of the tomographic images, suggesting that the present sparse data coverage can hardly detect reliably any temporal variations of the velocity anomalies in the Wen-An source area. In future studies it is necessary to improve the resolution of crustal tomography to the size of the rupture zone and utilize identical seismic ray paths from the same pairs of sources and receivers in order to detect any temporal variations of the crust structure in the source area of a large earthquake.     

S-wave velocity images of the Dead Sea Basin provided by ambient seismic noise

Based on passive seismic interferometry applied to ambient seismic noise recordings between station pairs belonging to a small-scale array, we have obtained shear wave velocity images of the uppermost materials that make up the Dead Sea Basin. We extracted empirical Green’s functions from cross-correlations of long-term recordings of continuous data, and measured inter-station Rayleigh wave group velocities from the daily correlation functions for positive and negative correlation time lags in the 0.1–0.5 Hz bandwidth. A tomographic inversion of the travel times estimated for each frequency is performed, allowing the laterally varying 3-D surface wave velocity structure below the array to be retrieved. Subsequently, the velocity-frequency curves are inverted to obtain S-wave velocity images of the study area as horizontal depth sections and longitude- and latitude-depth sections. The results, which are consistent with other previous ones, provide clear images of the local seismic velocity structure of the basin. Low shear velocities are dominant at shallow depths above 3.5 km, but even so a spit of land with a depth that does not exceed 4 km is identified as a salt diapir separating the low velocities associated with sedimentary infill on both sides of the Lisan Peninsula. The lack of low speeds at the sampling depth of 11.5 km implies that there are no sediments and therefore that the basement is near 10–11 km depth, but gradually decreasing from south to north. The results also highlight the bowl-shaped basin with poorly consolidated sedimentary materials accumulated in the central part of the basin. The structure of the western margin of the basin evidences a certain asymmetry both whether it is compared to the eastern margin and it is observed in north–south direction. Infill materials down to ∼8 km depth are observed in the hollow of the basin, unlike what happens in the north and south where they are spread beyond the western Dead Sea shore.     

Location of the pilot borehole for investigations of reservoir triggered seismicity at Koyna,India

Artificial water reservoir triggered earthquakes are now known to have occurred at over 120 sites globally. The part played by the reservoirs in triggering is not exactly known due to lack of near field observations of triggered earthquakes. Koyna, located near the west coast of India, where triggered earthquakes have been occurring since 1962 provides an excellent site for near field observations of the target M ≥ 2 earthquakes. A 6 borehole seismic network has been deployed recently in the Koyna region at depths of 981–1522 m to improve the hypocenter locations. During May–December 2015, a total of 1039 earthquakes of M_L ≥ 0.5 were located using the borehole seismic network. The region is also monitored through a dense network of 23 surface broad-band stations. Our analysis indicates a significant improvement in the estimation of absolute locations of earthquakes with errors of the order of ± 300 m, combining both the networks. Based on seismicity, and logistics, a block of 2 × 2 km² area has been chosen for drilling the first pilot borehole of ~ 3 km depth, where M ≥ 2 earthquakes have been occurring frequently since 2005.     

Amount of Asian lithospheric mantle subducted during the India/Asia collision

Body wave seismic tomography is a successful technique for mapping lithospheric material sinking into the mantle. Focusing on the India/Asia collision zone, we postulate the existence of several Asian continental slabs, based on seismic global tomography. We observe a lower mantle positive anomaly between 1100 and 900 km depths, that we interpret as the signature of a past subduction process of Asian lithosphere, based on the anomaly position relative to positive anomalies related to Indian continental slab. We propose that this anomaly provides evidence for south dipping subduction of North Tibet lithospheric mantle, occurring along 3000 km parallel to the Southern Asian margin, and beginning soon after the 45 Ma break-off that detached the Tethys oceanic slab from the Indian continent. We estimate the maximum length of the slab related to the anomaly to be 400 km. Adding 200 km of presently Asian subducting slab beneath Central Tibet, the amount of Asian lithospheric mantle absorbed by continental subduction during the collision is at most 600 km. Using global seismic tomography to resolve the geometry of Asian continent at the onset of collision, we estimate that the convergence absorbed by Asia during the indentation process is ~ 1300 km. We conclude that Asian continental subduction could accommodate at most 45% of the Asian convergence. The rest of the convergence could have been accommodated by a combination of extrusion and shallow subduction/underthrusting processes. Continental subduction is therefore a major lithospheric process involved in intraplate tectonics of a supercontinent like Eurasia.     

Crustal structure in an onshore–offshore transitional zone near Hong Kong,northern South China Sea

To study the crustal structure beneath the onshore–offshore transitional zone, a wide-angle onshore–offshore seismic experiment was carried out in northern South China Sea near Hong Kong, using large volume airgun sources at sea and seismic stations on land. The crustal velocity model constructed from traveltime fitting shows that the sedimentary thickness abruptly increases seaward of the Dangan Islands based on the characteristics of Pg and Multiple Pg, and the crustal structure beneath the sedimentary layer is relatively simple. The Moho depth is about 25–28 km along the profile and the P-wave velocity increases gradually with depth. The velocities in the upper crust range from 5.5 to 6.4 km/s, while that in the lower crust is 6.4–6.9 km/s. It also reveals a low velocity zone with a width of more than 10 km crossing the crust at about 75–90 km distance, which suggests that the Littoral Fault Zone (LFZ) exists beneath the onshore–offshore transitional zone. The magnetism anomalies, bouguer gravity anomalies and active seismic zone along the coastline imply the LFZ is a main tectonic fault in the onshore–offshore area. Combined with two previously published profiles in the continental South China (L–G profile) and in the northern margin of South China Sea (OBS1993) respectively, we constructed a land-sea super cross-section about 1000 km long. The results show the onshore–offshore transitional zone is a border separating the unstretched and the stretched continental crust. The low velocity layer (LVL) in the middle crust was imaged along L–G profile. However, the high velocity layer (HVL) in the lower crust was detected along OBS1993. By analyzing the mechanisms of the LVL in the middle crust and HVL in the base of crust, we believe the crustal structures had distinctly different attributes in the continental South China and in the northern SCS, which indicates that the LFZ could be the boundary fault between them.     

Seismic structure and rheology of the crust under mainland China

The crust and upper mantle in mainland China were relatively densely probed with wide-angle seismic profiling since 1958, and the data have provided constraints on the amalgamation and lithosphere deformation of the continent. Based on the collection and digitization of crustal P-wave velocity models along related wide-angle seismic profiles, we construct several crustal transects across major tectonic units in mainland China. In our study, we analyzed the seismic activity, and seismic energy releases during 1970 and 2010 along them. We present seismogenic layer distribution and calculate the yield stress envelopes of the lithosphere along the transects, yielding a better understanding of the lithosphere rheology strength beneath mainland China. Our results demonstrate that the crustal thicknesses of different tectonic provinces are distinctively different in mainland China. The average crustal thickness is greater than 65 km beneath the Tibetan Plateau, about 35 km beneath South China, and about 36–38 km beneath North China and Northeastern China. For the basins, the thickness is ~ 55 km beneath Qaidam, ~ 50 km beneath Tarim, ~ 40 km beneath Sichuan and ~ 35 km beneath Songliao. Our study also shows that the average seismic P-wave velocity is usually slower than the global average, equivalent with a more felsic composition of crust beneath the four tectonic blocks of mainland China resulting from the complex process of lithospheric evolution during Triassic and Cenozoic continent–continent and Mesozoic ocean–continent collisions. We identify characteristically different patterns of seismic activity distribution in different tectonic blocks, with bi-, or even tri-peak distribution of seismic concentration in South Tibet, which may suggest that crustal architecture and composition exert important control role in lithosphere deformation. The calculated yield stress envelopes of lithosphere in mainland China can be divided into three groups. The results indicate that the lithosphere rheology structure can be described by jelly sandwich model in eastern China, and crème brulee models with weak and strong lower crust corresponding to lithosphere beneath the western China and Kunlun orogenic belts, respectively. The spatial distribution of lithospheric rheology structure may provide important constraints on understanding of intra- or inter-plate deformation mechanism, and more studies are needed to further understand the tectonic process(es) accompanying different lithosphere rheology structures.     

The Nebula Winter: The united view of the snowball Earth,mass extinctions,and explosive evolution in the late Neoproterozoic and Cambrian periods

Encounters with nebulae, such as supernova remnants and dark clouds in the galaxy, can lead to an environmental catastrophe on the Earth through the negative climate forcings and destruction of the ozone layer by enhanced fluxes of cosmic rays and cosmic dust particles. A resultant reduction in primary productivity leads to mass extinctions through depletion of oxygen and food starvations as well as anoxia in the ocean. The model shows three levels of hierarchical time variations caused by supernova encounters (1–10 kyrs), dark cloud encounters (0.1–10 Myrs), and starbursts (~ 100 Myrs), respectively. This “Nebula Winter” model can explain the catastrophic phenomena such as snowball Earth events, repeated mass extinctions, and Cambrian explosion of biodiversities which took place in the late Proterozoic era through the Cambrian period. The Late Neoproterozoic snowball Earth event covers a time range of ca. 200 Myrs long spanning from 770 Ma to the end of Cambrian period (488 Ma) with two snowball states called Sturtian and Marinoan events. Mass extinctions occurred at least eight times in this period, synchronized with large fluctuations in δ¹³C of carbonates in the sediment. Each event is likely to correspond to each nebula encounter. In other words, the late Neoproterozoic snowball Earth and Cambrian explosion are possibly driven by a starburst, which took place around 0.6 Ga in the Milky Way Galaxy. The evidences for a Nebula Winter can be obtained from geological records in sediment in the deep oceans at those times.