Subsurface imaging,TAIGER experiments and tectonic models of Taiwan

The seismicity, deformation rates and associated erosion in the Taiwan region clearly demonstrate that plate tectonic and orogenic activities are at a high level. Major geologic units can be neatly placed in the plate tectonic context, albeit critical mapping in specific areas is still needed, but the key processes involved in the building of the island remain under discussion. Of the two plates in the vicinity of Taiwan, the Philippine Sea Plate (PSP) is oceanic in its origin while the Eurasian Plate (EUP) is comprised partly of the Asian continental lithosphere and partly of the transitional lithosphere of the South China Sea basin. It is unanimously agreed that the collision of PSP and EU is the cause of the Taiwan orogeny, but several models of the underlying geological processes have been proposed, each with its own evolutionary history and implied subsurface tectonics.TAIGER (TAiwan Integrated GEodynamics Research) crustal- and mantle-imaging experiments recently made possible a new round of testing and elucidation. The new seismic tomography resolved structures under and offshore of Taiwan to a depth of about 200 km. In the upper mantle, the steeply east-dipping high velocity anomalies from southern to central Taiwan are clear, but only the extreme southern part is associated with seismicity; toward the north the seismicity disappears. The crustal root under the Central Range is strongly asymmetrical; using 7.5 km/s as a guide, the steep west-dipping face on the east stands in sharp contrast to a gradual east-dipping face on the west. A smaller root exists under the Coastal Range or slightly to the east of it. Between these two roots lies a well delineated high velocity rise spanning the length from Hualien to Taitung. The 3-D variations in crustal and mantle structures parallel to the trend of the island are closely correlated with the plate tectonic framework of Taiwan. The crust is thickest in the central Taiwan collision zone, and although it thins toward the south, the crust is over 30 km thick over the subduction in the south; in northern Taiwan, the northward subducting PSP collides with Taiwan and the crust thins under northern Taiwan where the subducting indenter reaches 50 km in depth. The low Vp/Vs ratio of around 1.6 at a mid-crustal depth of 25 km in the Central Range indicates that current temperatures could exceed 700 °C. The remarkable thickening of the crust under the Central Range, its rapid uplift without significant seismicity, its deep exhumation and its thermal state contribute to make it the core of orogenic activities on Taiwan Island.The expanded network during the TAIGER deployment captured broadband seismic data yielding enhanced S-splitting results with mainly SKS/SKKS data. The polarization directions of the fast S-waves follow very closely the structural trends of the island, supporting the concept of a vertically coherent Taiwan orogeny in the outer few hundred kilometers of the Earth.     

Shear velocity and Vp/Vs ratio structure of the crust beneath the southern margin of South China continent

We here present the results of the inverse modeling of crustal S-phases recorded from a 400-km-long seismic profile, with azimuth nearly N30W, from Lianxian, near Hunan Province, to Gangkou Island, near Guangzhou City, Guangdong Province, in the southern margin of South China continent. The finding in this case is that many shot gathers provided by this wide-angle seismic experiment show relatively strong reflected and refracted S-phases, in particular some crustal refractions (Sg waves) and Moho reflections (SmS waves or simply Sm waves). The P-wave velocity structure of the crust and uppermost mantle was already obtained through the interpretation of vertical-component shot gathers. Now, with constraints introduced by the P-wave velocity architecture and after picking up S-wave traveltime data on the seismograms, we have obtained the S-velocity model of the crust by adjusting these traveltimes but keeping the geometry of the crustal reflectors. Our results demonstrate: (1) the average crustal S-velocity is about 3.64 km/s to the northwest of the Wuchuan-Sihui fault, and 3.62 km/s to the southeast of this fault; (2) relatively constant S-velocity of about 3.42 km/s for the upper crust, 3.55 km/s for the middle crust and laterally varying shear velocity around 3.82 km/s for the lower crust; (3) correspondingly, Vp/Vs ratio is 1.73 for the upper crust, 1.71 for the middle crust and 1.74 for the lower crust. Both shear velocities and Vp/Vs ratio correlate well with the major active faults that break the study area, and show significant changes especially in the upper crust. High Poisson’s ratio (1.8) is observed at shallow depth beneath the Minzhong depression to the southeast of the Wuchuan-Sihui fault and the Huiyuan depression in the southern margin of South China continent. In contrast, a very low Vp/Vs ratio (1.68) is observed between 8 and 14 km depth beneath Huiyuan. At deeper depth, a high Vp/Vs ratio (1.76) is observed in the lower crust beneath the Minzhong depression.     

Seismic structure of the crust and uppermost mantle of South America and surrounding oceanic basins

We present a new set of contour maps of the seismic structure of South America and the surrounding ocean basins. These maps include new data, helping to constrain crustal thickness, whole-crustal average P-wave and S-wave velocity, and the seismic velocity of the uppermost mantle (P_n and S_n). We find that: (1) The weighted average thickness of the crust under South America is 38.17 km (standard deviation, s.d. ±8.7 km), which is ∼1 km thinner than the global average of 39.2 km (s.d. ±8.5 km) for continental crust. (2) Histograms of whole-crustal P-wave velocities for the South American crust are bi-modal, with the lower peak occurring for crust that appears to be missing a high-velocity (6.9–7.3 km/s) lower crustal layer. (3) The average P-wave velocity of the crystalline crust (P_cc) is 6.47 km/s (s.d. ±0.25 km/s). This is essentially identical to the global average of 6.45 km/s. (4) The average P_n velocity beneath South America is 8.00 km/s (s.d. ±0.23 km/s), slightly lower than the global average of 8.07 km/s. (5) A region across northern Chile and northeast Argentina has anomalously low P- and S-wave velocities in the crust. Geographically, this corresponds to the shallowly-subducted portion of the Nazca plate (the Pampean flat slab first described by Isacks et al., 1968), which is also a region of crustal extension. (6) The thick crust of the Brazilian craton appears to extend into Venezuela and Colombia. (7) The crust in the Amazon basin and along the western edge of the Brazilian craton may be thinned by extension. (8) The average crustal P-wave velocity under the eastern Pacific seafloor is higher than under the western Atlantic seafloor, most likely due to the thicker sediment layer on the older Atlantic seafloor.     

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.     

Genetic waveform modeling for the crustal structure in Northeast Japan

We propose a genetic algorithm (GA) search procedure for waveform modeling of local crustal earthquakes for optimal one-dimensional (1-D) crustal velocity model. Both waveforms and travel-time data are used for the structure determination. The use of travel times in model evaluation improves the waveform modeling performance in the sense of computation speed and accuracy. We applied this method to broadband waveforms of a local crustal earthquake (M 4.2) in Northeast Japan. P-wave velocities of the crustal model are found to be 4.95 ± 0.30, 5.9 ± 0.02, and 6.51 ± 0.20 km/s for a surface layer, upper crust and lower crust, respectively. The surface layer thickness and the Conrad and Moho depths are found to be 3.01 ± 0.8, 17.77 ± 0.4 and 34.59 ± 1.0 km, respectively. For epicentral distances <200 km, our synthetic waveforms match the observed ones generally well. Early arrivals are mainly observed at stations near the Pacific coast in the forearc area having a thinner crust. In contrast, delayed arrivals appear at stations near the volcanic front and back-arc areas where low-velocity anomalies exist due to the effect of the Pacific slab dehydration and the hot upwelling flows in the mantle wedge. In general, our results agree well with the main tectonic setting of the study area, which confirms the reliability of the proposed approach. Despite a 1-D velocity model is too simple to represent the complex crustal structure, it is still required for the conventional routine analysis of seismology, such as earthquake location and source parameter studies. The current approach is considered as a step toward the genetic full waveform modeling for the 3-D velocity model estimation.     

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.     

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 lithospheric thinning beneath the seismogenic Kachchh rift zone,Gujarat (India): Its implications toward the generation of the 2001 Bhuj earthquake sequence

A high-resolution passive seismic experiment in the Kachchh rift zone of the western India has produced an excellent dataset of several thousands teleseismic events. From this network, 500 good teleseismic events recorded at 14 mobile broadband sites are used to estimate receiver functions (for the 30–310° back-azimuth ranges), which show a positive phase at 4.5–6.1 s delay time and a strong negative phase at 8.0–11.0 s. These phases have been modeled by a velocity increase at Moho (i.e. 34–43 km) and a velocity decrease at 62–92 km depth. The estimation of crustal and lithospheric thicknesses using the inversion of stacked radial receiver functions led to the delineation of a marked thinning of 3–7 km in crustal thickness and 6–14 km in lithospheric thickness beneath the central rift zone relative to the surrounding un-rifted parts of the Kachchh rift zone. On an average, the Kachchh region is characterized by a thin lithosphere of 75.9 ± 5.9 km. The marked velocity decrease associated with the lithosphere–asthenoshere boundary (LAB), observed over an area of 120 km × 80 km, and the isotropic study of xenoliths from Kachchh provides evidence for local asthenospheric updoming with pockets of partial melts of CO₂ rich lherzolite beneath the Kachchh seismic zone that might have caused by rifting episode (at 88 Ma) and the associated Deccan thermal-plume interaction (at 65 Ma) episodes. Thus, the coincidence of the area of the major aftershock activity and the Moho as well as asthenospheric upwarping beneath the central Kachchh rift zone suggests that these pockets of CO₂-rich lherzolite partial melts could perhaps provide a high input of volatiles containing CO₂ into the lower crust, which might contribute significantly in the seismo-genesis of continued aftershock activity in the region. It is also inferred that large stresses in the denser and stronger lower crust (at 14–34 km depths) induced by ongoing Banni upliftment, crustal intrusive, marked lateral variation in crustal thickness and related sub-crustal thermal anomaly play a key role in nucleating the lower crustal earthquakes beneath the Kachchh seismic zone.     

Delineation of intra crustal horizon in Eastern Dharwar Craton – An aeromagnetic evidence

A correlative study of two mutually independent geophysical properties like magnetic susceptibility variations and shear wave velocity structure of the crust has been carried out in a part of the Eastern Dharwar Craton of Indian peninsular shield. Analysis of the aeromagnetic anomaly field over an area of 35,000 km² comprising the peninsular gneissic basement complex and a part of Cuddapah Basin has resulted in identification of two distinct magnetic horizons: one at a depth of 2 km and the other at a depth of 12 km. Correlation of these results with the inferences made by the inversion of Rayleigh wave phase velocity and other geophysical studies has confirmed the presence of a crustal layer at a depth of 12 km. This horizon has been inferred to be the depth to the lower boundary of the upper crust in this region.     

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.     

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).     

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.     

Geophysical,magmatic, and metallogenic manifestations of a mantle plume in the upper reaches of the Aldan and Amur Rivers

Gravity models of the crust and upper mantle to a depth of 100 km are analyzed to study structural relationships of tectonic and tectonophysical media of different rigidities with the distribution of shallow ore deposits above the Aldan-Zeya plume. The spatial correlation of ore clusters and districts with high crustal viscosity inhomoheneities at depths of 10, 20, and 35 km shows distinct stepwise behavior. On the other hand, media of decreased viscosity are observed in the lower crust (at depths of 25–30 km), subcrustal (40–50 km) layers, and asthenosphere (at a depth below 70 km). They are related to chambers of the complete or partial melting (heat sources) of magmatic and ore occurrences near the Earth’s surface. Lateral metallogenic zoning in the spatial distribution of the ore deposits is due to the spread and redistribution of magmas and ore-forming fluids, shielded by rigid plates in the lower crust. A naturally determined series of ore parageneses is observed from center to flanks of the plume: Au, Mo → Au, Ag, Pb, Zn → Au, Pb, Zn → Au, W → Au, Sb → W, Sn → Sn. The mutual position of the tectonomagmatic structures of different ranks within the plume head obeys hierarchical and fractal laws.     

Seismogenesis of clustered seismicity beneath the Kangra–Chamba sector of northwest Himalaya: Constraints from 3D local earthquake tomography

To investigate subsurface structure and seismogenic layers, 3D velocity inversion was carried out in the source zone of 1905 Kangra earthquake (M8.0) in the northwestern Himalaya. P-wave and S-wave phase data of 159 earthquakes recorded by a network of 21 stations were used for this purpose. Inverted velocity tomograms up to a depth range of 18 km show significant variations of 14% in Vp and Vs and 6% in the Vp/Vs across the major tectonic zones in the region. Synthesis of seismicity pattern, velocity structure, distinctive focal mechanisms coupled with nature of stress distribution allows mapping of three different source regions that control regional seismotectonics. Accumulating strains are partly consumed by sliding of Chamba Nappe to the southwest through reverse-fault movements along Chamba/Panjal/Main Boundary Thrusts. This coupled with normal-fault type displacements along Chenab Normal Fault in the north account for low magnitude widespread seismicity in upper 8–10 km of the crust. At intermediate depths from 8 to 15 km, adjusting to residual compressive stresses, the detachment or lower end of the MBT slips to produce thrust dominated seismicity. Nucleation of secondary stresses in local NE–SW oriented structure interacts in complex manner with regional stresses to generate normal type earthquakes below the plane of detachment and therefore three seismic regimes at different depths produce intense seismicity in a block of 30 × 30 km² centered NE to the epicenter of Kangra earthquake.     

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.     

Sea-level change and paleogeographic reconstructions,southern Vancouver Island,British Columbia,Canada

Forty-eight new and previously published radiocarbon ages constrain deglacial and postglacial sea levels on southern Vancouver Island, British Columbia. Sea level fell rapidly from its high stand of about +75 m elevation just before 14 000 cal BP (12 000 radiocarbon yrs BP) to below the present shoreline by 13 200 cal BP (11 400 radiocarbon years BP). The sea fell below its present level 1000 years later in the central Strait of Georgia and 2000 years later in the northern Strait of Georgia, reflecting regional differences in ice sheet retreat and downwasting. Direct observations only constrain the low stand to be below ?11 m and above ?40 m. Analysis of the crustal isostatic depression with equations utilizing exponential decay functions appropriate to the Cascadia subduction zone, however, places the low stand at ?30 ± 5 m at about 11 200 cal BP (9800 BP). The inferred low stand for southern Vancouver Island, when compared to the sea-level curve previously derived for the central Strait of Georgia to the northwest, generates differential isostatic depression that is consistent with the expected crustal response between the two regions. Morphologic and sub-bottom features previously interpreted to indicate a low stand of ?50 to ?65 m are re-evaluated and found to be consistent with a low stand of ?30 ± 5 m. Submarine banks in eastern Juan de Fuca Strait were emergent at the time of the low stand, but marine passages persisted between southern Vancouver Island and the mainland. The crustal uplift presently occurring in response to the Late Pleistocene collapse of the southwestern sector of the Cordilleran Ice Sheet amounts to about 0.1 mm/yr. The small glacial isostatic adjustment rate is a consequence of low-viscosity mantle in this tectonically active region.     

S-wave velocity structure of the North China from inversion of Rayleigh wave phase velocity

We constructed the S-wave velocity structure of the crust and uppermost mantle (10–100 km) beneath the North China based on the teleseismic data recorded by 187 portable broadband stations deployed in this region. The traditional two-step inversion scheme was adopted. Firstly, we measured the interstation fundamental Rayleigh wave phase velocity of 10–60 s and imaged the phase velocity distributions using the Tarantola inversion method. Secondly, we inverted the 1-D S-wave velocity structure with a grid spacing of 0.25° × 0.25° and constructed the 3-D S-wave velocity structure of the North China. The 3-D S-wave velocity model provides valuable information about the destruction mechanism and geodynamics of the North China Craton (NCC). The S-wave velocity structures in the northwestern and southwestern sides of the North–South Gravity Lineament (NSGL) are obviously different. The southeastern side is high velocity (high-V) while the northeastern side is low velocity (low-V) at the depth of 60–80 km. The upwelling asthenosphere above the stagnated Pacific plate may cause the destruction of the Eastern Block and form the NSGL. A prominent low-V anomaly exists around Datong from 50 to 100 km, which may due to the upwelling asthenosphere originating from the mantle transition zone beneath the Western Block. The upwelling asthenosphere beneath the Datong may also contribute to the destruction of the Eastern Block. The Zhangjiakou-Penglai fault zone (ZPFZ) may cut through the lithosphere and act as a channel of the upwelling asthenosphere. A noticeable low-V zone also exists in the lower crust and upper mantle lid (30–50 km) beneath the Beijing–Tianjin–Tangshan (BTT) region, which may be caused by the upwelling asthenosphere through the ZPFZ.     

Crustal structure and nature of emplacement of the 85°E Ridge in the Mahanadi offshore based on constrained potential field modeling: Implications for intraplate plume emplaced volcanism

An integrated interpretation of multi-channel seismic reflection, gravity and magnetic datasets belonging to northern most part of the 85°E Ridge in the Mahanadi offshore is carried out to study the crustal structure and mode of its emplacement. The basement structure map of the ridge reveals that it is 130–150 km wide and is composed of an eastern high which appears as a continuous, broad and smooth topographyand the western high characterized by several steep isolated highs. The seismic velocities reported for the first time over the ridge indicate several sedimentary sequences ranging in velocities between 1.6 and 4.0 km/s above the acoustic basement top. The salient aspects of the sedimentary velocities are; a low velocity layer (2.6–3.2 km/s) within the Cretaceous sequence in the intervening depressions encompassing the flank region, and a regionally widespread higher velocity layer (3.5–3.8 km/s) belonging to the Eocene–Oligocene section overlying the ridge. A layer having a velocity of 4.2–4.7 km/s probably made of volcanoclastic rocks is observed immediately below the acoustic basement. The sediment isopach maps presented here for three major horizons are used to compute the 3-D sediment gravity effect to obtain a crustal Bouguer anomaly map of the region. Detailed analysis of the gravity and magnetic anomaly maps clearly demonstrates the continuity of ridge up to the Mahanadi coast at Chilka Lake. Seismically constrained gravity and magnetic models indicate that the ridge is composed of volcanic material that was emplaced on continental crust in the shelf-slope areas and over the oceanic crust in the deep offshore areas. The modeled crustal structure below the ridge further indicates volcanic emplacement of the ridge on a relatively younger lithosphere. We propose two alternative models for the emplacement of the ridge.     

Crustal structure of the Gujarat region,India: New constraints from the analysis of teleseismic receiver functions

In this study, receiver function analysis is carried out at 32 broadband stations spread all over the Gujarat region, located in the western part of India to image the sedimentary structure and investigate the crustal composition for the entire region. The powerful Genetic Algorithm technique is applied to the receiver functions to derive S-velocity structure beneath each site. A detail image in terms of basement depths and Moho thickness for the entire Gujarat region is obtained for the first time. Gujarat comprises of three distinct regions: Kachchh, Saurashtra and Mainland. In Kachchh region, depth of the basement varies from around 1.5 km in the eastern part to 6 km in the western part and around 2–3 km in the northern part to 4–5 km in the southern part. In the Saurashtra region, there is not much variation in the depth of the basement and is between 3 km and 4 km. In Gujarat mainland part, the basement depth is 5–8 km in the Cambay basin and western edge of Narmada basin. In other parts of the mainland, it is 3–4 km. The depth of Moho beneath each site is obtained using stacking algorithm approach. The Moho is at shallower depth (26–30 km) in the western part of Kachchh region. In the eastern part and epicentral zone of the 2001 Bhuj earthquake, large variation in the Moho depths is noticed (36–46 km). In the Saurashtra region, the crust is more thick in the northern part. It varies from 36–38 km in the southern part to 42–44 km in the northern part. In the mainland region, the crust is more thick (40–44 km) in the northern and southern part and is shallow in Cambay and Narmada basins (32–36 km). The large variations of Poisson’s ratio across Gujarat region may be interpreted as heterogeneity in crustal composition. High values of σ (∼0.30) at many sites in Kachchh and few sites in Saurashtra and Mainland regions may be related to the existence of high-velocity lower crust with a mafic/ultramafic composition and, locally, to the presence of partial melt. The existing tectono-sedimentary models proposed by various researchers were also examined.