Acceleration and evolution of faults: An example from the Hunter Mountain–Panamint Valley fault zone,Eastern California

We present new space geodetic data indicating that the present slip rate on the Hunter Mountain–Panamint Valley fault zone in Eastern California (5.0 ± 0.5 mm/yr) is significantly faster than geologic estimates based on fault total offset and inception time. We interpret this discrepancy as evidence for an accelerating fault and propose a new model for fault initiation and evolution. In this model, fault slip rate initially increases with time; hence geologic estimates averaged over the early stages of the fault's activity will tend to underestimate the present-day rate. The model is based on geologic data (total offset and fault initiation time) and geodetic data (present day slip rate). The model assumes a monotonic increase in slip rate with time as the fault matures and straightens. The rate increase follows a simple Rayleigh cumulative distribution. Integrating the rate-time path from fault inception to present-day gives the total fault offset.     

Active tectonics of the Northern Rif (Morocco) from geomorphic and geochronological data

We present results of a geomorphological and morphotectonic analysis of the northeastern part of the Rif. We show that the present day kinematics of the Rif is characterized by active deformation along the Trougout and Nekor faults in the North-East. Digital Elevation Models of offset drainage features (streams, fluvial terraces) allow determining a normal-left-lateral motion along the Trougout fault and a left-lateral strike-slip motion along the Nekor fault. Preliminary ³He cosmogenic dates of tectonic markers yield vertical and horizontal slip rates of ∼0.9 mm/yr and ∼0.5 mm/yr, respectively along the Trougout fault. The present-day localized transtension seen in the north-eastern Rif morphology (Ras Tarf) is coeval with uplifted marine terraces near the Al Hoceima Bay. U/Th dating of shells yield an average uplift rate of ∼0.2 mm/yr during the past 500 ka. These data show that active transtension in the northeastern Rif is also associated with uplift. These new morphotectonic constraints are consistent with the GPS measurements showing southwestward overall motion of most of the Rif belt with respect to stable Africa.     

GPS-derived velocity field of the Iznik-Mekece segment of the North Anatolian Fault Zone

Space-based tectonic studies on the western part of the North Anatolian Fault Zone (NAFZ) have been conducted over two decades. After the August 17, 1999, Izmit earthquake (M_w = 7.4), this region attracted greater scientific interest, and the collected data became more valuable. The Geodesy Department of the Kandilli Observatory and Earthquake Research Institute (KOERI) at Bogazici University established three micro-geodetic networks to the east of Akyazi, east of Iznik, and west of Lake Sapanca in the eastern part of the Marmara region; GPS data have been continually collected at these locations since 1994. The NAFZ branches out in the western part of the Marmara region and extends up to the Aegean Sea. Segments of the fault passing through the Marmara Sea are considered active, and this has increased concern regarding imminent earthquakes. Conventional geodetic measurements made between 1990 and 1994 are not sufficient for monitoring small movements. However, GPS has played a very important role in detecting such deformations in the area after 1994. The Iznik network, with 10 points, is bilaterally located on the Iznik-Mekece fault. Six years of GPS data for 2004–2010 collected for the monitoring of crustal deformation showed that the Iznik-Mekece fault segment moves westward at about 22 ± 1 mm/yr with respect to the Eurasia fixed reference frame. The GPS observations show that there is no strain accumulation in the area.     

Spatial-temporal characterization of the San Andreas Fault by fault-zone trapped waves at seismic experiment site,Parkfield, California

In this article, we review our previous research for spatial and temporal characterizations of the San Andreas Fault (SAF) at Parkfield, using the fault-zone trapped wave (FZTW) since the middle 1980s. Parkfield, California has been taken as a scientific seismic experimental site in the USA since the 1970s, and the SAF is the target fault to investigate earthquake physics and forecasting. More than ten types of field experiments (including seismic, geophysical, geochemical, geodetic and so on) have been carried out at this experimental site since then. In the fall of 2003, a pair of scientific wells were drilled at the San Andreas Fault Observatory at Depth (SAFOD) site; the main-hole (MH) passed a ～200-m-wide low-velocity zone (LVZ) with highly fractured rocks of the SAF at a depth of ～3.2 km below the wellhead on the ground level (Hickman et al., 2005; Zoback, 2007; Lockner et al., 2011). Borehole seismographs were installed in the SAFOD MH in 2004, which were located within the LVZ of the fault at ～3-km depth to probe the internal structure and physical properties of the SAF. On September 282004, a M6 earthquake occurred ～15 km southeast of the town of Parkfield. The data recorded in the field experiments before and after the 2004 M6 earthquake provided a unique opportunity to monitor the co-mainshock damage and post-seismic heal of the SAF associated with this strong earthquake. This retrospective review of the results from a sequence of our previous experiments at the Parkfield SAF, California, will be valuable for other researchers who are carrying out seismic experiments at the active faults to develop the community seismic wave velocity models, the fault models and the earthquake forecasting models in global seismogenic regions.     

Velocity field and tectonic strain in Southern Spain and surrounding areas derived from GPS episodic measurements

The goal of this paper is to study the velocity field and deformation parameters in Southern Spain and surrounding areas (Ibero-Maghrebian region) using GPS episodic measurements. Results are compared to those previously published as well as deformation parameters derived from seismic data. For this purpose, a geodetic GPS network of 12 stations was observed during eight field campaigns from 1998 to 2005 by the San Fernando Naval Observatory (ROA), Spain. Relative GPS velocities in the Gulf of Cadiz with respect to the stable part of Eurasia are ～4.1 mm/yr in a NW–SE to NNW–SSE direction. In the Betics, Alboran Sea and North of Morocco, velocities are ～4.4 mm/yr in a NW–SE direction, and they are ～2.3 mm/yr in a N–S direction in the eastern part of the Iberian Peninsula. These results are in agreement with the anticlockwise rotation of the African plate. GPS strain tensors are determined from the velocity model, to obtain a more realistic crustal deformation model. The Gulf of Cadiz is subjected to uniform horizontal compression in a NNW–SSE direction, with a rotation to N–S in the Alboran Sea and Northern Morocco. An extensional regime in a NW–SE direction, which rotates to W–E, is present in the Internal Betics area. In the Betic, Alboran Sea and North of Morocco regions we compare seismic deformation rates from shallow earthquakes with the determined GPS deformation rates. The comparison indicates a seismic coupling of 27%, while the remaining 73% might be generated in aseismic processes. Deformations measured in the Ibero-Maghrebian region with GPS could be interpreted in terms of either elastic loading or ductile deformation.     

New insights in the geodynamics of the Lipari–Vulcano area (Aeolian Archipelago,southern Italy) from geological,geodetic and seismological data

Geological, geodetic and seismological data have been analyzed in order to frame the Lipari–Vulcano complex (Aeolian archipelago, southern Italy) into the geodynamic context of the southeastern Tyrrhenian Sea. It is located at the northern end of a major NNW–SSE trending right-lateral strike-slip fault system named “Aeolian–Tindari–Letojanni” which has been interpreted as a lithospheric discontinuity extending from the Aeolian Islands to the Ionian coast of Sicily and separating two different tectonic domains: a contractional one to the west and an extensional one to the north-east. Structural field data consist of structural measurements performed on well-exposed fault planes and fractures. The mesostructures are mostly represented by NW–SE striking normal faults with a dextral-oblique component of motion. Minor structures are represented by N–S oriented joints and tension gashes widespread over the whole analyzed area and particularly along fumarolized sectors. The analyzed seismological dataset (from 1994 to 2013) is based on earthquakes with magnitude ranging between 1.0 and 4.8. The hypocenter distribution depicts two major alignments corresponding to the NNW–SSE trending Aeolian–Tindari–Letojanni fault system and to the WNW–ESE oriented Sisifo–Alicudi fault system. GPS data analysis displays ∼3.0 mm/yr of active shortening between the two islands, with a maximum shortening rate of about 1.0 × 10⁻¹³ s⁻¹, between La Fossa Caldera and south of Vulcanello. This region is bounded to the north by an area where the maximum values of shear strain rates, of about 0.7 × 10⁻¹³ s⁻¹ are observed. This major change occurs in the area south of Vulcanello that is also characterized by a transition in the way of the vertical axis rotation. Moreover, both the islands show a clear subsidence process, as suggested by negative vertical velocities of all GPS stations which exhibit a decrease from about −15 to −7 mm/yr from north to south. New data suggest that the current kinematics of the Lipari–Vulcano complex can be framed in the tectonic context of the eastward migrating Sisifo–Alicudi fault system. This is dominated by transpressive tectonics in which contractional and minor extensional structures can coexist with strike-slip motion.     

2014年美国加州纳帕M_W6.1地震断层参数的Sentinel-1A InSAR反演

2014年8月24日,在美国加州旧金山海湾北部的纳帕地区发生了MW6.1地震.发震断层是西纳帕断裂系统中的一部分,但是该断层之前并未被足够重视.本文利用欧洲空间局最近发射成功并刚刚投入使用的Sentinel-1A卫星获取的第一对同震干涉像对(20140807-20140831),得到了该地震的地表同震形变场,结合震后24h内区域GPS同震形变资料作为约束条件,反演了纳帕地震的断层几何参数以及滑动分布.Sentinel-1A干涉结果表明,此次地震造成了明显的地面形变,视线向最大抬升和最大沉降量均达到了10cm.联合反演结果表明,该发震断层的走向为344°,倾角为80°.主要破裂以右旋走滑为主,平均倾滑角为-146.5°,最大倾滑量达到了1.1m,位于地表下约4km,存在明显的滑动亏损现象.此次地震,累计释放地震矩达1.5×1018 N·m,约合矩震级MW6.1.该结果略小于InSAR单独约束结果,可能与Sentinel-1A像对中包含的快速震后形变分量有关.     

Initiation of the San Jacinto Fault and its Interaction with the San Andreas Fault: Insights from Geodynamic Modeling

The San Andreas Fault (SAF) is the Pacific-North American plate boundary, yet in southern California a significant portion of the relative plate motion is accommodated by the San Jacinto Fault (SJF). Here we investigate the initiation of the SJF and its interaction with the SAF in a three-dimensional visco-elasto-plastic finite-element model. The model results show that the restraining bend of the southern SAF causes strain localization along the SJF, thus may have contributed to its initiation. Slip on the SJF tends to reduce slip rate on the SAF and enhance deformation in the Eastern California Shear Zone. The initiation of the SJF and its interaction with the SAF reflect the evolving plate boundary zone as it continuously seeks the most efficient way to accommodate the relative plate motion.     

Nanoscale porosity in SAFOD core samples (San Andreas Fault)

With transmission electron microscopy (TEM) we observed nanometer-sized pores in four ultracataclastic and fractured core samples recovered from different depths of the main bore hole of the San Andreas Fault Observatory at Depth (SAFOD). Cutting of foils with a focused ion beam technique (FIB) allowed identifying porosity down to the nm scale. Between 40 and 50% of all pores could be identified as in-situ pores without any damage related to sample preparation. The total porosity estimated from TEM micrographs (1–5%) is comparable to the connected fault rock porosity (2.8–6.7%) estimated by pressure-induced injection of mercury. Permeability estimates for cataclastic fault rocks are 10^? 21–10^? 19 m² and 10^? 17 m² for the fractured fault rock. Porosity and permeability are independent of sample depth. TEM images reveal that the porosity is intimately linked to fault rock composition and associated with deformation. The TEM-estimated porosity of the samples increases with increasing clay content. The highest porosity was estimated in the vicinity of an active fault trace. The largest pores with an equivalent radius > 200 nm occur around large quartz and feldspar grains or grain-fragments while the smallest pores (equivalent radius < 50 nm) are typically observed in the extremely fine-grained matrix (grain size < 1 μm). Based on pore morphology we distinguish different pore types varying with fault rock fabric and alteration. The pores were probably filled with formation water and/or hydrothermal fluids at elevated pore fluid pressure, preventing pore collapse. The pore geometry derived from TEM observations and BET (Brunauer, Emmett and Teller) gas adsorption/desorption hysteresis curves indicates pore blocking effects in the fine-grained matrix. Observations of isolated pores in TEM micrographs and high pore body to pore throat ratios inferred from mercury injection suggest elevated pore fluid pressure in the low permeability cataclasites, reducing shear strength of the fault.     

Geodetic constraints on active tectonics of the Western Mediterranean: Implications for the kinematics and dynamics of the Nubia-Eurasia plate boundary zone

We present GPS observations in Morocco and adjacent areas of Spain from 15 continuous (CGPS) and 31 survey-mode (SGPS) sites extending from the stable part of the Nubian plate to central Spain. We determine a robust velocity field for the W Mediterranean that we use to constrain models for the Iberia-Nubia plate boundary. South of the High Atlas Mountain system, GPS motions are consistent with Nubia plate motions from prior geodetic studies. We constrain shortening in the Atlas system to <1.5 mm/yr, 95% confidence level. North of the Atlas Mountains, the GPS velocities indicate Nubia motion with respect to Eurasia, but also a component of motion normal to the direction of Nubia-Eurasia motion, consisting of southward translation of the Rif Mountains in N Morocco at rates exceeding 5 mm/yr. This southward motion appears to be directly related to Miocene opening of the Alboran Sea. The Betic Mountain system north of the Alboran Sea is characterized by WNW motion with respect to Eurasia at ～1–2 mm/yr, paralleling Nubia-Eurasia relative motion. In addition, sites located in the Betics north of the southerly moving Rif Mountains also indicate a component of southerly motion with respect to Eurasia. We interpret this as indicating that deformation associated with Nubia-Eurasia plate motion extends into the southern Betics, but also that the Betic system may be affected by the same processes that are causing southward motion of the Rif Mountains south of the Alboran Sea. Kinematic modeling indicates that plate boundary geometries that include a boundary through the Straits of Gibraltar are most compatible with the component of motion in the direction of relative plate motion, but that two additional blocks (Alboran-Rif block, Betic Mountain block), independent of both Nubia and Eurasia are needed to account for the motions of the Rif and Betic Mountains normal to the direction of relative plate motion. We speculate that the southward motions of the Alboran-Rif and Betic blocks may be related to mantle flow, possibly induced by southward rollback of the subducted Nubian plate beneath the Alboran Sea and Rif Mountains.     

Deriving fault-slip histories to test for secular variation in slip,with examples from the Kunlun and Awatere faults

Although offset and age data from displaced landforms are essential for identifying earthquake clusters and thus testing whether faults slip at uniform or secularly varying rates, it is not clear how the uncertainties in such measurements should be propagated so as to yield a robust fault-slip history (i.e., record of fault displacement over time). Here we develop a Monte Carlo approach for estimating the distribution of geologically reasonable fault-slip histories that fit the offset and age data from a population of dated and displaced landforms. The model assumes that the landforms share common faulting histories, the offset and age constraints are correct, and the fault has not reversed shear sense. Analysis of the model results yields both a precise average slip rate, in the case where a linear fit is applied to the data, and a best-fit fault-slip history, in the case where the linear constraint is removed. The method can be used to test for secular variation in slip because the uncertainty on this best-fit history is quantified. By applying the method to previously published morphochronologic data from faulted late Quaternary terrace risers along the Kunlun fault in China and the Awatere fault in New Zealand, we have assessed the extent to which our modeled average slip rates match previously reported values and the data support previous interpretations of uniform slip rate. The Kunlun data set yields average slip rates of 8.7 + 3.6/?2.1 mm/yr and 5.1 + 1.6/?1.2 mm/yr (68.27% confidence), for the central and eastern reaches of the fault, respectively, both of which match previously published slip rates. Our analysis further indicates that these fault reaches have both slipped uniformly over the latest Quaternary. In contrast, analysis of data from the Saxton River site along the Awatere fault reveals a mid-Holocene deceleration in slip rate from 6.2 + 1.6/?1.4 mm/yr to 2.8 + 1.0/?0.6 mm/yr. This result contradicts previous interpretations of uniform slip along the Awatere fault. The Monte Carlo method we present here for quantifying fault-slip histories using the offset and age data from a population of faulted landforms provides an important tool for distinguishing temporally uniform from secularly varying fault slip.     

Structural heterogeneity of an ultra-low-velocity zone beneath the Philippine Islands: Implications for core–mantle chemical interactions induced by massive partial melting at the bottom of the mantle

An ScP phase reflected and converted at the core–mantle boundary (CMB) beneath the region east of the Philippine Islands shows clear pre- and postcursors, recorded on short-period seismic networks in Japan. These waveform variations can be explained by interaction of the ScP wavefield with thin layers at the CMB. The results of forward modeling of double-array stacks reveal two different structural heterogeneities in the lowermost mantle beneath the region east of the Philippine Islands. One of the structures represents a decreased velocity, and increased density across the reflector at the lowermost ～10 km of the mantle, with P- and S-wave velocity reductions of 5–10% and ～30%, respectively, and an increase in density of 5–10%. Another structure consists of a pair of reflectors at ～10 km and ～5 km above the CMB, both of which are characterized by reduced P- and S-wave velocities. The upper reflector is the interface of a low-velocity zone in which P- and S-wave velocities decrease of 10% and 30%, respectively, accompanied by an extremely large increase in density (20–25%). The lower reflector is characterized by a 25% reduction in S-wave velocity relative to the above low-velocity layer, as well as a 5% decrease in P-wave velocity and no change in density. The nature of the low-velocity zone detected locally at the CMB is comparable with that of ultra-low-velocity zones (ULVZs) observed by various seismic probes in the South Pacific and Central America. Extensive observations of the ULVZ beneath the region east of the Philippine Islands indicate massive partial melting at the bottom of the mantle. Low-S-velocity basal layer partly detected within the ULVZ may be resulting from core–mantle chemical interactions, driven by massive partial melting.     

Intra-plate deformation in west-central Europe

The central part of Europe north of the alpine orogenic belt is generally seen as a relatively stable area of the western tip of the Eurasian plate. Indeed, up to now, no geodetically significant motions have been detected although an active rift system running roughly in SSE–NNW direction along the Rhine valley could have some effect on the stability of this region. Presently, the increasing accuracy of geodetic point motions should allow the study of small motions at levels down to nearly 0.1 mm/yr. We start our investigation with a closer look at the ‘true’ accuracy and significance of GPS derived point velocities of permanent stations. We compare and discuss the different levels of formal errors obtained by the three analysis centers considered in this study (EPN, JPL and SOPAC) and present additional ways of assessing the accuracy using the redundancy offered by different independent analyses and multiple systems operating at one site. On the average, all results indicate that a one-sigma level of ±0.3 mm/yr can be seen as a conservative estimate for the horizontal accuracy of point motions in central Europe. On the basis of this assumption we find that at present, the actual velocity field as determined by different analysis groups and centers does not show any significant east–west extensional deformation. We do however see a prominent north–south compressional velocity gradient of about 1 mm/yr/1000 km (1 nanostrain/yr) which could be associated with the Alpine thrust in conjunction with a south-directed horizontal component of the Fennoscandian Glacial Isostatic Adjustment.     

Hypoxia in Manila Bay,Philippines during the northeast monsoon

Herein we present results from one of the first extensive bay-wide oceanographic surveys of Manila Bay, wherein 31 stations were sampled during the northeast monsoon (cold and dry season). A band of hypoxic bottom water (dissolved oxygen <2.8 mg/L) spanned the midsection of the bay from east to west. Bottom nitrate concentrations (5.7–16.8 μM; avg. 11.1 μM) and total organic carbon values in sediments (1.7–3.1%; avg. 2.4%) were high in the midsection, which coincided with the band of hypoxic bottom water. Physical processes and site-specific accumulation of organic material likely lead to hypoxic conditions in Manila Bay, even during the northeast monsoon period when the water column is relatively well mixed. The results of this study complement the previously reported widespread hypoxia that occurs during the rainy season. Thus, hypoxia may be pervasive in the bay throughout the year, although it varies in intensity and spatial extent.     

Processing and preliminary interpretation of noisy high-resolution seismic reflection/refraction data across the active Ostler Fault zone,South Island,New Zealand

In an attempt to understand the structure of active faults as they emerge from bedrock into shallow semi-consolidated and unconsolidated sediments, we have recorded a comprehensive high-resolution seismic reflection/refraction data set across the Ostler Fault zone on the central South Island of New Zealand. This fault zone, which absorbs 1–2 mm/yr of compression associated with oblique convergence of the Pacific and Australian tectonic plates, consists of a series of surface-rupturing N–S trending, west-dipping reverse faults that offset a thick sequence of Quaternary glacial outwash and late Neogene fluvio-lacustrine sediments of the Mackenzie Basin. Our study focuses on a region of the basin where two non-overlapping fault segments are separated by a transfer zone. Deformation in this area is accommodated by offsets on multiple small faults and by folding in their hanging walls. The seismic data with source and receiver spacing of 6 and 3 m and nominal CMP fold of 60 was acquired along twelve 1.2 km long lines orthogonal to fault strike and an additional 1.6 km long tie-line parallel to fault strike. The combination of active deformation and shallow glacial outwash sediments results in particularly complicated seismic data, such that application of relatively standard processing schemes yields only poor quality images. We have designed a pre- and post-stack reflection/refraction processing scheme that focuses on minimising random and source-generated noise, determining appropriate static corrections and resolving contrasting reflection dips. Application of this processing scheme to the Ostler Fault data provides critical information on fault geometry and offset and on sedimentary structures from the surface to ～ 800 m depth. Our preliminary interpretation of one of the lines includes complex deformation structures with folding and multiple subsidiary fault splays on either side of a ～ 50° west-dipping primary fault plane.     

Solar activity imprints in tree ring width from Chile (1610–1991)

We have investigated the solar activity signal in tree ring data from two locations in Chile. The tree ring time series extended over a period of ∼400 yr. Spectral and wavelet analysis techniques were employed. We have found evidence for the presence of the solar activity Schwabe (∼11 yr), Hale (∼22 yr), fourth-harmonic of the 208-yr Suess cycle (∼52 yr) and Gleissberg (∼80 yr) cycles. The Gleissberg cycle of tree ring data is in anti-phase with solar activity. Wavelet and cross-wavelet techniques revealed that the periods found are intermittent, possibly because solar activity signals observed in tree rings are mostly due to solar influence on local climate (rainfall, temperature, and cloud cover) where trees grow up. Further, cross-wavelet analysis between sunspot and tree ring time series showed that the cross power around the 11 yr solar cycle is more significant during periods of high solar activity (grand maximum) than during periods of low solar activity (grand minimum). As Glaciar Pio XI is practically at the Pacific Ocean level, the tree-ring response may be stronger due to the heating of the Pacific Ocean water following an increase of the solar radiation incidence rather than at the higher altitudes of Osorno region.     

DORIS as a potential part of a Global Geodetic Observing System

We have processed all available DORIS data from all available satellites, except Jason-1 over the past 10 years (from January 1993 to April 2003). Weekly solutions have been produced for stations positions coordinates, geocenter motion and scale factor stability. We present here accuracy presently achievable for all types of potential geodetic products. Typically weekly stations positions can be derived with a repeatability of 1.0–1.5 cm using data from 5 satellites simultaneously, showing the significant improvement in precision that has been gained recently using the additional new DORIS satellites. As an example, we show how such new results can detect displacement from large magnitude earthquakes, such as the 2003 Denali fault earthquake in Alaska. Displacements of −5 cm in latitude and +2 cm in longitude were easily detected using the DORIS data and are confirmed by recent GPS determination. The terrestrial reference frame was also well be monitored with DORIS during this 10-year period. Other geodetic products, such as tropospheric corrections for atmospheric studies are also analyzed. Finally, we discuss here the possible advantages and weaknesses of the DORIS system as additional geodetic tool, in conjunction with the already existing GPS, VLBI and SLR services, to participate in an Global Geodetic Observing System (GGOS).     

Drainage migration and out of sequence thrusting in Bhalukpong,Western Arunachal Himalaya,India

In the foreland regions of the Western Arunachal Himalaya (WAH), geological studies along the Kameng river (between Tipi village and the Himalayan Frontal Thrust (HFT)) reveal four levels of unpaired terraces and a paired terrace. In WAH, wrench deformation of HFT zone resulted in a SE propagation of the Balipara anticline and it is suggested that the Mikir high basement controls its orientation. Ages of terrace surfaces from Siwaliks suggest that since the Late Pleistocene, Kameng River migrated at a rate varying between ∼7.5 cm/yr in upper reaches and ∼13.5 cm/yr towards northeast due to HFT related uplift. In the Brahmaputra plains, luminescence ages of abandoned paleochannel deposits suggest eastward shifting of the Kameng river at an average rate of ∼1 m/yr. Field evidences between Bhalukpong and Tipi villages show Pliocene strath and Quaternary terrace surfaces, displaced by faults that do not correspond to the mapped faults in the foreland region. We interpret them as out-of-sequence thrusts (OOSTs). This is the first such report of OOST in the NE Himalaya. Presence of active OOST is inferred by similar age (∼1 ka) and differing incision rates of the surface of same terrace (T_2b) in adjacent locations. This suggests that OOSTs in the western Arunachal Siwalik are <1 ka. Average slip rate and horizontal shortening rate on OOST during the Holocene, are calculated as ∼12 mm/yr and 7 mm/yr respectively. Thus any estimation of Holocene shortening in the Siwalik therefore, needs to incorporate slip along the OOSTs given that it accommodates a significant amount of N-S compression of the Himalayan fold-and-thrust belt. The reason for OOST in the WAH Siwalik foreland is discussed in terms of the critical wedge dynamics arising from erosion via tectonics-climate interaction. We estimate a minimum slip rate of Siwalik as ∼27 mm/yr during the Holocene and suggest acceleration in shortening rates east of Bhutan.     

Occurrence and distribution of polycyclic aromatic hydrocarbons in surface sediments of San Diego Bay marinas

Polycyclic aromatic hydrocarbons (PAHs) have garnered much attention due to their bioaccumulation, carcinogenic properties, and persistence in the environment. Investigation of the spatial distribution, composition, and sources of PAHs in sediments of three recreational marinas in San Diego Bay, California revealed significant differences among marinas, with concentrations in one site exceeding 16,000 ng g^? 1. ‘Hotspots’ of PAH concentration suggest an association with stormwater outfalls draining into the basins. High-molecular weight PAHs (4–6 rings) were dominant (> 86%); the average percentage of potentially carcinogenic PAHs was high in all sites (61.4–70%) but ecotoxicological risks varied among marinas. Highly toxic benzo(a)pyrene (BaP) was the main contributor (> 90%) to the total toxic equivalent quantity (TEQ) in marinas. PAHs in San Diego Bay marina sediments appear to be derived largely from pyrogenic sources, potentially from combustion products that reach the basins by aerial deposition and stormwater drainage from nearby streets and parking lots.     

Crustal deformation before the 2008 Wenchuan MS8.0 earthquake studied using GPS data

We used GPS velocities from approximately 700 stations in western China to study the crustal deformation before the Wenchuan M_S8.0 earthquake. The processing methods included analyses of the strain rate field, inversion of fault locking and the GPS velocity profiles. The GPS strain rate in the E-W direction in the Qinghai-Tibet block shows that extensional deformation was dominant in the western region of the block (west of 92.5° E), while compressive deformation predominated in the eastern region of the block (from 92.5° E to 100° E). On a regional scale, the hypocentral region of the Wenchuan earthquake was located at the edge of an intense compression deformation zone of about 1.9 × 10⁻⁸/a in an east-west direction. The characteristic deformation in the seismogenic fault was compressive with a dextral component. The compression deformation rate was greater in the fault's western region than in its eastern region, and the strain accumulation was very slow on the fault scale. The results of a fault locking inversion show that the locking fraction and slip deficit was greater in the middle-northern section of the seismogenic fault than in the southern section. The GPS velocity profile before the Wenchuan earthquake shows that the compression deformation was smaller than the dextral deformation, which is asymmetrical with respect to the distribution of co-seismic displacement. These deformation characteristics should provide some clues to the Wenchuan earthquake which occurred in the later period of the earthquake cycle.