Deep electrical structure over the igneous arc of the Indo Burman Orogen in Sagaing province,Myanmar from magnetotelluric studies

Magnetotelluric studies over the igneous arc of the Indo Burman range in the Sagaing province of Myanmar have delineated the high resistivity Indian plate subducting westwards beneath the Burmese block to depths of 30 km and beyond. The thick moderately resistive (20–100 Ω m) layer overlying the subducting Indian plate may be due to the low resistivity sediments. The entire region is covered with prominent sedimentary layer with a conductance varying between 20 and 3000 S showing a general increase from the east to west, suggesting that their thickness increases toward the west. The large unsystematic variations in the conductance are indicative of the widely varying depositional environments and also possible vertical block movements during the course of their deposition. A west dipping low resistivity zone to the east of Burmese block seems to demarcate its eastern limit, suggesting the possibility of a hitherto unknown deep seated fault, which is also supported by the several earthquake foci located over this zone. The nature of the crustal movements over this fault is not immediately apparent. Possibility exists that the Sagaing fault is an en echelon fault and the present feature observed here is a part of this en echelon fault. The possibility of channel flows of the weakened rocks in the deep crust observed in the vicinity of the eastern Himalayan syntaxis may also cause such low resistivity zones.     

Geophysical constraints on geodynamic processes at convergent margins: A global perspective

Convergent margins, being the boundaries between colliding lithospheric plates, form the most disastrous areas in the world due to intensive, strong seismicity and volcanism. We review global geophysical data in order to illustrate the effects of the plate tectonic processes at convergent margins on the crustal and upper mantle structure, seismicity, and geometry of subducting slab. We present global maps of free-air and Bouguer gravity anomalies, heat flow, seismicity, seismic Vs anomalies in the upper mantle, and plate convergence rate, as well as 20 profiles across different convergent margins. A global analysis of these data for three types of convergent margins, formed by ocean–ocean, ocean–continent, and continent–continent collisions, allows us to recognize the following patterns. (1) Plate convergence rate depends on the type of convergent margins and it is significantly larger when, at least, one of the plates is oceanic. However, the oldest oceanic plate in the Pacific ocean has the smallest convergence rate. (2) The presence of an oceanic plate is, in general, required for generation of high-magnitude (M > 8.0) earthquakes and for generating intermediate and deep seismicity along the convergent margins. When oceanic slabs subduct beneath a continent, a gap in the seismogenic zone exists at depths between ca. 250 km and 500 km. Given that the seismogenic zone terminates at ca. 200 km depth in case of continent–continent collision, we propose oceanic origin of subducting slabs beneath the Zagros, the Pamir, and the Vrancea zone. (3) Dip angle of the subducting slab in continent–ocean collision does not correlate neither with the age of subducting oceanic slab, nor with the convergence rate. For ocean–ocean subduction, clear trends are recognized: steeply dipping slabs are characteristic of young subducting plates and of oceanic plates with high convergence rate, with slab rotation towards a near-vertical dip angle at depths below ca. 500 km at very high convergence rate. (4) Local isostasy is not satisfied at the convergent margins as evidenced by strong free air gravity anomalies of positive and negative signs. However, near-isostatic equilibrium may exist in broad zones of distributed deformation such as Tibet. (5) No systematic patterns are recognized in heat flow data due to strong heterogeneity of measured values which are strongly affected by hydrothermal circulation, magmatic activity, crustal faulting, horizontal heat transfer, and also due to low number of heat flow measurements across many margins. (6) Low upper mantle Vs seismic velocities beneath the convergent margins are restricted to the upper 150 km and may be related to mantle wedge melting which is confined to shallow mantle levels.     

P-wave tomography and dynamics of the crust and upper mantle beneath western Tibet

The continental collision between the Indian and Asian plates plays a key role in the geologic and tectonic evolution of the Tibetan plateau. In this article we present high-resolution tomographic images of the crust and upper mantle derived from a large number of high-quality seismic data from the ANTILOPE project in western Tibet. Both local and distant earthquakes were used in this study and 35,115 P-wave arrival times were manually picked from the original seismograms. Geological and geochemical results suggested that the subducting Indian plate has reached northward to the Lhasa terrane, whereas our new tomography shows that the Indian plate is currently sub-horizontal and underthrusting to the Jinsha river suture at depths of ~ 100 to ~ 250 km, suggesting that the subduction process has evolved over time. The Asian plate is also imaged clearly from the surface to a depth of ~ 100 km by our tomography, and it is located under the Tarim Basin north of the Altyn Tagh Fault. There is no obvious evidence to show that the Asian plate has subducted beneath western Tibet. The Indian and Asian plates are separated by a prominent low-velocity zone under northern Tibet. We attribute the low-velocity zone to mantle upwelling, which may account for the warm crust and upper mantle beneath that region, and thus explain the different features of magmatism between southern and northern Tibet. But the upwelling may not penetrate through the whole crust. We propose a revised geodynamic model and suggest that the high-velocity zones under Lhasa terrane may reflect a cold crust which has interrupted the crustal flow under the westernmost Tibetan plateau.     

Rigid Indian plate: Constraints from GPS measurements

We analyze GPS data from 26 sites located on the Indian plate and along its boundary. The large spatial coverage of the Indian plate by these sites and longer data duration helped us in refining the earlier estimates of the Euler pole for the Indian plate rotation. Our analysis suggests that the internal deformation of the Indian plate is very low (< 1–2 mm/year) and the entire plate interior region largely behaves as a rigid plate. Specifically, we did not infer any significant difference in motion on sites located north and south of the Narmada Son failed rift region, the most prominent tectonic feature within the Indian plate and a major source of earthquakes. Our analysis also constrains the motion across the Indo-Burmese wedge, Himalayan arc, and Shillong Plateau and Kopili fault in the NE India.     

Evidence of Mid-ocean ridge and shallow subduction forearc magmatism in the Nagaland-Manipur ophiolites,northeast India: constraints from mineralogy and geochemistry of gabbros and associated mafic dykes

We discuss here the mineralogical and geochemical characteristics of mafic intrusive rocks from the Nagaland-Manipur Ophiolites (NMO) of Indo-Myanmar Orogenic Belt, northeast India to define their mantle source and tectonic environment. Mafic intrusive sequence in the NMO is characterized by hornblende-free (type-I) and hornblende-bearing (type-II) rocks. The type-I is further categorized as mafic dykes (type-Ia) of tholeiitic N-MORB composition, having TiO₂ (0.72–1.93 wt.%) and flat REE patterns (La_N/Yb_N = 0.76–1.51) and as massive gabbros (type-Ib) that show alkaline E-MORB affinity, having moderate to high Ti content (TiO₂ = 1.18 to 1.45 wt.%) with strong LREE-HREE fractionations (La_N/Yb_N = 4.54–7.47). Such geochemical enrichment from N-MORB to E-MORB composition indicates mixing of melts derived from a depleted mantle and a fertile mantle/plume source at the spreading center. On the other hand, type-II mafic intrusives are hornblende bearing gabbros of SSZ-type tholeiitic composition with low Ti content (TiO₂ = 0.54 wt.%–0.86 wt.%) and depleted LREE pattern with respect to HREE (La_N/Yb_N = 0.37–0.49). They also have high Ba/Zr (1.13–2.82), Ba/Nb (45.56–151.66) and Ba/Th (84.58–744.19) and U/Th ratios (0.37–0.67) relative to the primitive mantle, which strongly represents the melt composition generated by partial melting of depleted lithospheric mantle wedge contaminated by hydrous fluids derived from subducting oceanic lithosphere in a forearc setting. Their subduction related origin is also supported by presence of calcium-rich plagioclase (An_16.6–32.3). Geothermometry calculation shows that the hornblende bearing (type-II) mafic rocks crystallized at temperature in range of 565°–625 °C ± 50 (at 10 kbar). Based on these available mineralogical and geochemical evidences, we conclude that mid ocean ridge (MOR) type mafic intrusive rocks from the NMO represent the section of older oceanic crust which was generated during the divergent process of the Indian plate from the Australian plate during Cretaceous period. Conversely, the hornblende-bearing gabbros (type-II) represent the younger oceanic crust which was formed at the forearc region by partial melting of the depleted mantle wedge slightly modified by the hydrous fluids released from the subducting oceanic slab during the initial stage of subduction of Indian plate beneath the Myanmar plate.     

Present-day crustal motion around the Pamir Plateau from GPS measurements

We have used geodetic techniques to improve constraints on the crustal motion of the Pamir Plateau. Three campaigns of Global Position System data acquisition between 2011 and 2015 demonstrate that, in association with the India–Asia collision, a complex pattern of crustal motion exists in the Pamir Plateau. In a north–south direction from the Indian Plate to the Hazak Block, the crust has absorbed ~ 35 mm/yr of shortening, of which ~ 35% is distributed around the Hindu Kush region (~ 12 mm/yr), and another ~ 35% is taken up around the Alai Valley (also ~ 12 mm/yr). Global Position System measurements also show ~ 5 mm/yr of shortening between the Pamir Plateau and the Tajik Basin, whereas between the Pamir and the Tarim Basin, an ~ 10 mm/yr extension rate is observed. With respect to the stable Eurasian Plate, the Pamir rotates counterclockwise at a rate of ~ 1.822°Myr^− 1, with an Euler pole positioned about the west end of the Tajik Basin (37.03 ± 0.74°N, 65.89 ± 0.12°E). The strain rate field calculated from Global Position System velocities reveals that the crustal motion is consistent with localized deformation around the Hindu Kush and the Alai Valley, the latter representing a zone with strong shallow seismic activity.     

The large aftershocks triggered by the 2011 Mw 9.0 Tohoku-Oki earthquake,Japan

The Mw 9.0 Tohoku-Oki earthquake that occurred off the Pacific coast of Japan on March 11, 2011, was followed by thousands of aftershocks, both near the plate interface and in the crust of inland eastern Japan. In this paper, we report on two large, shallow crustal earthquakes that occurred near the Ibaraki-Fukushima prefecture border, where the background seismicity was low prior to the 2011 Tohoku-Oki earthquake. Using densely spaced geodetic observations (GPS and InSAR datasets), we found that two large aftershocks in the Iwaki and Kita-Ibarake regions (hereafter referred to as the Iwaki earthquake and the Kita-Ibarake earthquake) produced 2.1 m and 0.44 m of motion in the line-of-sight (LOS), respectively. The azimuth-offset method was used to obtain the preliminary location of the fault traces. The InSAR-based maximum offset and trace of the faults that produced the Iwaki earthquake are consistent with field observations. The fault location and geometry of these two earthquakes are constrained by a rectangular dislocation model in a multilayered elastic half-space, which indicates that the maximum slips for the two earthquakes are 3.28 m and 0.98 m, respectively. The Coulomb stress changes were calculated for the faults following the 2011 Mw 9.0 Tohoku-Oki earthquake based on the modeled slip along the fault planes. The resulting Coulomb stress changes indicate that the stresses on the faults increased by up to 1.1 MPa and 0.7 MPa in the Iwaki and Kita-Ibarake regions, respectively, suggesting that the Tohoku-Oki earthquake triggered the two aftershocks, supporting the results of seismic tomography.     

Surface displacements in Japan before the 11 March 2011 M9.0 Tohoku-Oki earthquake

Daily resolution data retrieved from the 1243 ground-based Global Positioning System (GPS) stations in Japan are utilized to expose surface displacements before the destructive M9 Tohoku-Oki earthquake (March 11, 2011). Variations in the residual GPS data, in which effects of the long-term plate movements, short-term noise and frequency-dependent variations have been removed through a band-pass filter via the Hilbert–Huang transform, are compared with parameters of the focal mechanism associated with the Tohoku-Oki earthquake for validation. Analytical results show that the southward movements, which were deduced from the residual displacements and agree with the strike of the rupture fault, became evident on the 65th day before the Tohoku-Oki earthquake. This observation suggests that the shear stress played an important role in the seismic incubation period. The westward movements, which are consistent with the angle of the maximum horizontal compressive stress, covered entire Japan and formed an impeded area (142°E, 42°N) about 75 km away from the epicenter on the 47th day prior to the earthquake. The horizontal displacements integrated with the vertical movements from the residual GPS data are very useful to construct comprehensive images in diagnosing the surface deformation from destructive earthquakes along the subduction zone.     

Dynamic subsidence of Eastern Australia during the Cretaceous

During the Early Cretaceous Australia's eastward passage over sinking subducted slabs induced widespread dynamic subsidence and formation of a large epeiric sea in the eastern interior. Despite evidence for convergence between Australia and the paleo-Pacific, the subduction zone location has been poorly constrained. Using coupled plate tectonic–mantle convection models, we test two end-member scenarios, one with subduction directly east of Australia's reconstructed continental margin, and a second with subduction translated ~ 1000 km east, implying the existence of a back-arc basin. Our models incorporate a rheological model for the mantle and lithosphere, plate motions since 140 Ma and evolving plate boundaries. While mantle rheology affects the magnitude of surface vertical motions, timing of uplift and subsidence depends on plate boundary geometries and kinematics. Computations with a proximal subduction zone result in accelerated basin subsidence occurring 20 Myr too early compared with tectonic subsidence calculated from well data. This timing offset is reconciled when subduction is shifted eastward. Comparisons between seismic tomography and model temperature cross-sections, and an absence of subduction zone volcanism in eastern Australia in the Early Cretaceous provide support for the back-arc basin scenario.     

Neotectonics at the Arabian plate margins

Opening of the Red Sea is accompanied by convergence between the Arabian plate and Eurasia. Regional topography and structure favour gravity glide as the main driving force of plate translation. At the leading edge of the plate, the Zagros Mountains undergo coseismic serial folding which is equivalent to Holocene shortening by ∼20 mm/year and which has led to major episodes of coastal uplift of which the last was ∼1700 years BP. At the Jordan Rift transform, which bounds the Arabian plate on the west, a recurrence interval of ∼1600 years is reported for events of M_L≥5.5. The palaeomagnetic record for the last 3.2 Ma indicates an average spreading rate for the Red Sea of ∼20 mm/year; there is some evidence that hydrothermal activity in the Red Sea is pulsatory, with a period of ∼2000 year, and that it reflects discontinuous spreading. The Holocene neotectonic records of the Zagros, the Jordan Rift and the Red Sea are the product of complex plate interactions and of the accumulation and release of strain in the crust along the plate margins. But they also reflect elastic strain energy storage and release within the Arabian plate, whence parallels in the period of major deformation episodes in the three deforming zones and the apparent discrepancy between the seismic moment predicted by plate kinematics and that recorded in the Zagros. Any associated intraplate deformation, if detected geodetically, would thus help the assessment of seismic hazard.     

Salt marshes as archives of recent relative sea level change in West Greenland

Temperate latitude salt marshes are a proven environment from which high precision (±0.10–0.20 m) relative sea level (RSL) records can be developed over recent decades and centuries. Such records provide an important link between millennial histories of RSL change and instrumental records. The high latitude salt marshes in Greenland have not previously been explored as potential archives of recent RSL change. Here we develop four diatom-based transfer functions using contemporary diatom data collected from three salt marshes located 40 km south of the coastal town of Sisimiut, West Greenland. Our preferred model has a good fit between observed and predicted elevations (r² = 0.94) and a root mean square error prediction of ±0.19 m. We apply the four models to a short sediment profile collected from one of the marshes that formed between c. 600 cal. year BP and the present. Three of the four models predict the same trend in which RSL rose from ?0.55 ±0.19 m mean tide level (MTL) to ?0.05 ± 0.19 m MTL between c. 600 and 400 cal. year BP at a rate of c. 2.7 mm year^?1. After 400 cal. year BP RSL slowed and remained stable until the present day. The results of this study demonstrate that Greenland salt marshes are potentially valuable archives of data for developing quantitative estimates of RSL change during the last few centuries, thereby bridging the gap between existing millennial-scale approaches and more recent direct observations of ice sheet behaviour and associated vertical land motions.     

Late Quaternary slip rates and slip partitioning on the Alpine Fault,New Zealand

Published geological data on Late Quaternary offsets on the Alpine Fault, New Zealand, have been assembled into a common format and analysed with respect to uncertainties. Uncertainties arise mainly from measurement of offset features, relating apparent offsets to actual fault slip, and dating the offset features. Despite the considerable uncertainties, the data form a coherent set consistent with a relatively constant rate of strike-slip of 27±5 mm/year between Milford Sound and Hokitika. This rate represents 70–75% of the fault-parallel interplate motion. North of the confluence with the Hope Fault, the rate drops substantially. Dip-slip rates, on the other hand, show considerable variation along strike, rising to a maximum of more than 10 mm/year in the central section and decreasing to zero at the southern end. Partitioning of c. 25% of the interplate slip on to structures east of the Alpine Fault occurs in the central section, consistent with predictions from critical wedge models. The partitioning of all the fault-normal component of displacement on to other structures in the south may be related, in part, to a doubling in width of the deforming wedge to the east. Most probably, however, the fault-normal displacement is mainly accommodated by underthrusting of the Australian plate offshore, due to a change in the nature of the crust from continental to oceanic.     

Field evidences from northern Dead Sea Fault Zone (South Turkey): New findings for the initiation age and slip rate

The left-lateral strike–slip Dead Sea Fault Zone (DSFZ) extends from the Red Sea in the south to the East Anatolian Fault Zone (EAFZ) in the north. This study examines the northern part of the DSFZ around Amik Basin and presents surface and subsurface geological evidence for the Quaternary activity and initiation age of the northernmost DSFZ. The DSFZ extends N–S in the south of the Amik Basin where clear geological and morphological evidence exists for faulting. Geological observations around Amik Basin, analyses of borehole data and electrical resistivity profiles within the Amik Basin indicate that the activity of the northern DSFZ started after Pliocene in the Amik Basin. Subsurface data in the basin suggest that the DSFZ offsets a pre-Quaternary basin sinistrally by about 7.9 km. The offset pre-Quaternary basin suggests at least 4.94 ± 0.13 mm/year slip rate for the northern part of the DSFZ. The Karasu Fault Zone (KFZ) extends in an en-echelon pattern along the western margin of the Karasu Valley and it transfers the significant amount of slip from DSFZ to the EAFZ.     

Holocene sea-level changes along the North Carolina Coastline and their implications for glacial isostatic adjustment models

We have synthesized new and existing relative sea-level (RSL) data to produce a quality-controlled, spatially comprehensive database from the North Carolina coastline. The RSL database consists of 54 sea-level index points that are quantitatively related to an appropriate tide level and assigned an error estimate, and a further 33 limiting dates that confine the maximum and minimum elevations of RSL. The temporal distribution of the index points is very uneven with only five index points older than 4000 cal a BP, but the form of the Holocene sea-level trend is constrained by both terrestrial and marine limiting dates. The data illustrate RSL rapidly rising during the early and mid Holocene from an observed elevation of ?35.7 ± 1.1 m MSL at 11062–10576 cal a BP to ?4.2 m ± 0.4 m MSL at 4240–3592 cal a BP.We restricted comparisons between observations and predictions from the ICE-5G(VM2) with rotational feedback Glacial Isostatic Adjustment (GIA) model to the Late Holocene RSL (last 4000 cal a BP) because of the wealth of sea-level data during this time interval. The ICE-5G(VM2) model predicts significant spatial variations in RSL across North Carolina, thus we subdivided the observations into two regions. The model forecasts an increase in the rate of sea-level rise in Region 1 (Albemarle, Currituck, Roanoke, Croatan, and northern Pamlico sounds) compared to Region 2 (southern Pamlico, Core and Bogue sounds, and farther south to Wilmington). The observations show Late Holocene sea-level rising at 1.14 ± 0.03 mm year^?1 and 0.82 ± 0.02 mm year^?1 in Regions 1 and 2, respectively. The ICE-5G(VM2) predictions capture the general temporal trend of the observations, although there is an apparent misfit for index points older than 2000 cal a BP. It is presently unknown whether these misfits are caused by possible tectonic uplift associated with the mid-Carolina Platform High or a flaw in the GIA model. A comparison of local tide gauge data with the Late Holocene RSL trends from Regions 1 and 2 support the spatial variation in RSL across North Carolina, and imply an additional increase of mean sea level of greater than 2 mm year^?1 during the latter half of the 20th century; this is in general agreement with historical tide gauge and satellite altimetry data.     

Decarbonation of subducting slabs: Insight from petrological–thermomechanical modeling

Subduction of heterogeneous lithologies (sediments and altered basalts) carries a mixture of volatile components (H₂O ± CO₂) into the mantle, which are later mobilized during episodes of devolatilization and flux melting. Several petrologic and thermodynamic studies investigated CO₂ decarbonation to better understand carbon cycling at convergent margins. A paradox arose when investigations showed little to no decarbonation along present day subduction geotherms at subarc depths despite field based observations. Sediment diapirism is invoked as one of several methods for carbon transfer from the subducting slab. We employ high-resolution 2D petrological–thermomechanical modeling to elucidate the role subduction dynamics has with respect to slab decarbonation and the sediment diapirism hypothesis. Our thermodynamic database is modified to account for H₂O–CO₂ binary fluids via the following lithologies: GLOSS average sediments (H₂O: 7.29 wt.% & CO₂: 3.01 wt.%), carbonated altered basalts (H₂O: 2.63 wt.% & CO₂: 2.90 wt.%), and carbonated peridotites (H₂O: 1.98 wt.% & CO₂: 1.50 wt.%). We include a CO₂ solubility P–x[H₂O wt.%] parameterization for sediment melts. We parameterize our model by varying two components: slab age (20, 40, 60, 80 Ma) and convergence velocity (1, 2, 3, 4, 5, 6 cm year^− 1). 59 numerical models were run and show excellent agreement with the original code base. Three geodynamic regimes showed significant decarbonation. 1) Sedimentary diapirism acts as an efficient physical mechanism for CO₂ removal from the slab as it advects into the hotter mantle wedge. 2) If subduction rates are slow, frictional coupling between the subducting and overriding plate occurs. Mafic crust is mechanically incorporated into a section of the lower crust and undergoes decarbonation. 3) During extension and slab rollback, interaction between hot asthenosphere and sediments at shallow depths result in a small window (~ 12.5 Ma) of high integrated CO₂ fluxes (205 kg m^− 3 Ma^− 1).     

The hydraulic conductivity field and groundwater flow in the unconfined aquifer system of the Keta Strip,Ghana

This study investigates the hydraulic conductivity field and the groundwater flow pattern as predicted by a calibrated steady state groundwater flow model for the Keta Strip, southeastern Ghana. The hydraulic conductivity field is an important parameter in evaluating aquifer properties in space, and in general basin-wide groundwater resources evaluation and management. This study finds that the general hydraulic conductivity of the unconsolidated unconfined aquifer system of the Keta Strip ranges between 2 m/d and 20 m/d, with an average of 15 m/d. The spatial variation in horizontal hydraulic conductivity appears to take the trend in the variations in the nature of the material in space. Calibrated groundwater recharge suggests that 6.9–34% of annual precipitation recharges the shallow aquifer system. This amount of recharge is significant and suggests high fortunes in terms of groundwater resources development for agriculture and industrial activities in the area. A spatial distribution of groundwater recharge from precipitation is presented in this study. The spatial pattern appears to take the form of the distribution in horizontal hydraulic conductivity, and suggests that the vertical hydraulic conductivity takes the same pattern of spatial variation as the horizontal hydraulic conductivity. This is consistent with observations in other areas. The resulting groundwater flow is dominated by local flow systems as the unconfined system is quite shallow. A general northeast – southwest flow pattern has been observed in the study area.     

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.     

ESR dating of the Wangsan fault,South Korea

The Wangsan fault is exposed at Kyeongju, Korea. The andesite is unconformably covered by Quaternary alternating conglomerate and sandstone deposits. The unconformity is cut by the thrust fault which displaces a hanging wall block of about 30 m. Exposed at the surface of this fault is a light gray and brown fault gouge, about 40–120 cm thick. In order to test the consistency of ESR ages from a single fault gouge zone, we collected six gouge samples systematically along and across the Wangsan fault. We found that six samples collected from the same gouge zone show consistent ESR date estimates (average 550 ka). Because fault rock is rare along faults in unconsolidated sediments due to low confining stress near the surface, we consider that this fault gouge has been moved up with the hanging wall block along the fault. The estimated average uplift rate of the hanging wall block is about 0.04 cm/year based on the age of the displaced Quaternary deposits (vertical separation; about 20 m) dated by OSL dating methods. The depth of the fault gouge at the time of reactivation, which was estimated from uplift rate and the ESR ages, is about 220 m. Therefore, we conclude that the results of ESR age estimates represent the time of reactivation of the fault gouge at a depth of 220 m in the past during fault movement, because later movements which occurred during uplift near the surface may not have zeroed ESR signals significantly for ESR dating of fault movements.     

Effects of citrate on the dissolution and transformation of biotite,analyzed by chemical and atomic force microscopy

The aim of this study was to investigate the dissolution and transformation characteristics of phyllosilicate under low molecular weight organic acids in the farmland environment (pH 4.0–8.0). Changes of dissolution and morphology of biotite were evaluated using chemical extraction experiments and in situ/ex situ atomic force microscopy (AFM) with fluids of citric acid (CA) solution at pH 4.0, 6.0, and 8.0. Results of extracting experiments show that CA solutions contributed to the release rate of potassium (K), silicon (Si), and aluminum (Al) from biotite relative to a control aqueous solution. In situ AFM observations indicate that the dissolution of biotite from the biotite (0 0 1) surface occurred on the terrace, segment, and fringe of pits, while new etch pits did not readily form on biotite (0 0 1) surfaces in aqueous solutions. However, dissolution rates of terraces can be greatly accelerated with the help of citrate. In pH 4.0 CA solution, 70 min dissolution reactions of biotite (0 0 1) surfaces result in more etch pits than in pH 6.0 and 8.0 solutions. In addition, the transformation of biotite occurred simultaneously with the dissolution process. Secondary coating was observed on the biotite (0 0 1) surface after 140 h of immersion in a weak acid environment. Thus, the protons have a dominant role in the dissolution process of biotite with organic (carboxyl) acting as a catalyst under acidic condition. Based on the theory of interactions on a water–mineral interface in a weak acid environment, dissolution of biotite starts from defect/kink sites on the surface, one layer by one layer, and develops along the [h k 0] direction. A secondary coating that forms on the biotite (0 0 1) surface may restrain the formation and growth of etch pits, whereas this process may have a positive role on the stability of soil structure during long-term soil management.     

A study on the uplift mechanism of Tongjiezi dam using a coupled hydro-mechanical model

During filling in 1992 of the Tongjiezi reservoir, in Southwest China, it was noticed that the dam body and the rock masses on the right side of the valley were uplifted up to 22.2 mm and 24.3 mm respectively. After reservoir filling in 1993, the uplift continued, but the rate decreased, and the measured maximum uplift increased from 22.2 mm to 27.5 mm in the dam body and from 24.3 mm to 28.9 mm in the rock masses from 1993 to 2004. Based on the geological features of the dam site and observed geomechanics data, a numerical method is used to evaluate the representative elementary volume (REV) and to identify the parameters related to the mechanical and hydraulic properties of the rock mass. Furthermore, a coupled hydro-mechanical model is presented to describe the uplift process. The simulated results agree well with those measured. In addition, the time-dependent deformation is studied in the laboratory and numerically. It is concluded that, with the specific hydrogeological conditions at the Tongjiezi dam site, hydro-mechanical coupling during and after the reservoir impoundment is the main factor contributing to the uplift, and the time-dependent deformation is due to the rheological behavior of rock masses under seepage pressure.