Strong ground motion simulation of the 2003 Bam, Iran, earthquake using the empirical Green’s function method

The 2003 Bam, Iran, earthquake caused catastrophic damage to the city of Bam and neighboring villages. Given its magnitude (M _w ) of 6.5, the damage was remarkably large. Large-amplitude ground motions were recorded at the Bam accelerograph station in the center of Bam city by the Building and Housing Research Center (BHRC) of Iran. We simulated the Bam earthquake acceleration records at three BHRC strong-motion stations—Bam, Abaraq, and Mohammad-Abad—by the empirical Green’s function method. Three aftershocks were used as empirical Green’s functions. The frequency range of the empirical Green’s function simulations was 0.5–10 Hz. The size of the strong motion generation area of the mainshock was estimated to be 11 km in length by 7 km in width. To estimate the parameters of the strong motion generation area, we used 1D and 2D velocity structures across the fault and a combined source model. The empirical Green’s function method using a combination of aftershocks produced a source model that reproduced ground motions with the best fit to the observed waveforms. This may be attributed to the existence of two distinct rupture mechanisms in the strong motion generation area. We found that the rupture starting point for which the simulated waveforms best fit the observed ones was near the center of the strong motion generation area, which reproduced near-source ground motions in a broadband frequency range. The estimated strong motion generation area could explain the observed damaging ground motion at the Bam station. This suggests that estimating the source characteristics of the Bam earthquake is very important in understanding the causes of the earthquake damage.     

Degradation of Six Selected Ultraviolet Filters in Aquifer Materials Under Various Redox Conditions

Laboratory biodegradation batch studies were performed to investigate the degradation behavior of six selected UV filters, namely benzophenone‐3 (BP‐3), 3‐(4‐methylbenzylidene) camphor (4‐MBC), Octyl 4‐methoxycinnamate (OMC), Octocrylene (OC), 2‐(3‐t‐butyl‐2‐hydroxy‐5‐methylphenyl)‐5‐chloro benzotriazole (UV‐326), and 2‐(2’‐hydroxy‐5’‐octylphenyl)‐benzotriazole (UV‐329) in an aquifer microcosm (groundwater and aquifer sediment mixture) under aerobic and anaerobic (nitrate, sulfate, and Fe(III) reducing) conditions within 77 d. The results from the biodegradation experiments showed that the six UV filters were degraded well in the aquifer materials under different redox conditions. Rapid biodegradation was observed for BP‐3 and OMC in the aquifer materials, with their half‐lives of 1.5‐8.8 d and 1.3‐5.2 d, respectively. In most cases, aerobic conditions were more favorable for the degradation of the UV filters in aquifer materials. Relatively slow degradation of 4‐MBC, UV‐326, and UV‐329 under anaerobic conditions was noted with their half‐lives ranging between 47 d and 126 d, indicating potential persistence in anaerobic aquifers. The results showed that redox conditions could have significant effects on biodegradation of the UV filters in aquifers.     

Arbitrary amplitude ion-acoustic solitary waves in superthermal electron-positron-ion magnetoplasma

Arbitrary amplitude ion-acoustic solitary waves propagating in a magnetized plasma composed of positive ions, superthermal electrons and positrons are investigated. For this purpose, the ions are represented by the hydrodynamical fluid equations while the non-Maxwellian electrons and positrons densities are assumed to follow kappa (κ) distribution. The basic equations are reduced to a pseudoenergy-balance equation. Existence conditions for large amplitude solitary waves are presented. The analytical and numerical analysis of the latter show that the ion-acoustic solitary wave can propagate only in the subsonic region in our plasma system and it is significantly influenced by the plasma parameters. The present analysis could be helpful for understanding the nonlinear ion-acoustic solitary waves propagating in interstellar medium and pulsar wind, which contain an excess of superthermal particles.     

Stiffness pathologies in discrete granular systems: Bifurcation,neutral equilibrium,and instability in the presence of kinematic constraints

The paper develops the stiffness relationship between the movements and forces among a system of discrete interacting grains. The approach is similar to that used in structural analysis, but the stiffness matrix of granular material is inherently nonsymmetric because of the geometrics of particle interactions and of the frictional behavior of the contacts. Internal geometric constraints are imposed by the particles' shapes, in particular, by the surface curvatures of the particles at their points of contact. Moreover, the stiffness relationship is incrementally nonlinear, and even small assemblies require the analysis of multiple stiffness branches, with each branch region being a pointed convex cone in displacement space. These aspects of the particle-level stiffness relationship give rise to three types of microscale failure: neutral equilibrium, bifurcation and path instability, and instability of equilibrium. These three pathologies are defined in the context of four types of displacement constraints, which can be readily analyzed with certain generalized inverses. That is, instability and nonuniqueness are investigated in the presence of kinematic constraints. Bifurcation paths can be either stable or unstable, as determined with the Hill–Bažant–Petryk criterion. Examples of simple granular systems of three, 16, and 64 disks are analyzed. With each system, multiple contacts were assumed to be at the friction limit. Even with these small systems, microscale failure is expressed in many different forms, with some systems having hundreds of microscale failure modes. The examples suggest that microscale failure is pervasive within granular materials, with particle arrangements being in a nearly continual state of instability.     

Local versus regional processes: can soil characteristics overcome climate and fire regimes by modifying vegetation trajectories?

We analysed charcoal and pollen from sediments obtained from two lakes in the northwestern mixed‐wood Canadian boreal forest in order to reconstruct fire‐return intervals and vegetation dynamics over the last 8000 years. Sites were selected with contrasting soil properties (mesic versus dry‐sandy soils), allowing an estimation of the potential influence of soils on long‐term vegetation and fire dynamics. The sites likely experienced fewer fires during the period extending from 8000 to 4000 cal. a BP than over the last 4000 years. At both sites, eastern white pine (Pinus strobus) populations were most extensive shortly after deglaciation, with vegetation later shifting towards mixed woodlands with less P. strobus and more extensive Picea and Pinus banksiana populations. This gradual vegetation shift was probably induced by the establishment of colder and moister conditions along with a fire‐regime change. In spite of the parallel long‐term vegetation trajectories, vegetation composition differed between the two sites in both the past and present. Whereas Picea was more abundant at the mesic site, the fire‐adapted P. banksiana populations were more extensive at the sandy‐soil site. These differences in vegetation composition indicate that, in addition to climate changes and fire occurrence, soil properties also influenced vegetation dynamics. A likely increase in fire frequency in the Canadian boreal forest during the 21st century might therefore favour the expansion of these two disturbance‐adapted trees with spatial heterogeneity in the populations due to varying soil types. Copyright © 2012 John Wiley & Sons, Ltd.     

Mesozoic to Cenozoic tectono‐metamorphic history of the South Pamir–Hindu Kush (Chitral,NW Pakistan): Insights from phase equilibria modelling,and garnet–monazite petrochronology

The Karakoram–Hindu Kush–Pamir and adjacent Tibetan plateau belt comprise a series of Gondwana‐derived crustal fragments that successively accreted to the Eurasian margin in the Mesozoic as the result of the progressive Tethys ocean closure. These domains provide unique insights into the thermal and structural history of the Mesozoic to Cenozoic Eurasian plate margin, which are critical to inform the initial boundary conditions (e.g. crustal thickness, structure and thermo‐mechanical properties) for the subsequent development of the large and hot Tibetan–Himalaya orogen, and the associated crustal deformation processes. Using a combination of microstructural analyses, thermobarometry modelling and U–Th–Pb monazite and Lu–Hf garnet geochronology, the study reappraises the metamorphic history of exposed mid‐crustal metapelites in the Chitral region of the South Pamir–Hindu Kush (NW Pakistan). This study also demonstrates that trace elements in monazite (especially Y and Dy), combined with thermodynamical modelling and Lu–Hf garnet dating, provides a powerful integrated toolbox for constraining long‐lived and polyphased tectono‐metamorphic histories in all their spatial and temporal complexity. Rocks from the Chitral region were progressively deformed and metamorphosed at sub‐ and supra‐solidus conditions through at least four distinct episodes from the Mesozoic to the Cenozoic. Rocks were first metamorphosed at ~400–500°C and ~0.3 GPa in the Late Triassic–Early Jurassic (210–185 Ma), likely in response to the accretion of the Karakoram during the Cimmerian orogeny. Pressure and temperature subsequently increased by ~0.3 GPa and 100°C in the Early‐ to Mid Cretaceous (140–80 Ma), coinciding with the intrusion of calcalkaline granitic plutons across the Karakoram and Pamir regions. This event is interpreted as the record of crustal thickening and the development of a proto‐plateau within the Eurasian margin due to a long‐lived episode of slab flattening in an Andean‐type margin. Peak metamorphism was reached in the Late Eocene–Early Oligocene (40–30 Ma) at conditions of 580–600°C and ~0.6 GPa and 700–750°C and 0.7–0.8 GPa for the investigated staurolite schists and sillimanite migmatites respectively. This crustal heating up to moderate anatexis likely resulted in the underthrusting of the Indian plate after a NeoTethyan slab‐break off or to the Tethyan Himalaya–Lhasa microcontinent collision and subsequent oceanic slab flattening. Near‐isothermal decompression/exhumation followed in the Late Oligocene (28–23 Ma) as marked by a pressure decrease in excess of ~0.1 GPa. This event was coeval with the intrusion of the 24 Ma Garam Chasma leucogranite. This rapid exhumation is interpreted to be related to the reactivation of the South Pamir–Karakoram suture zone during the ongoing collision with India. The findings of this study confirm that significant crustal shortening and thickening of the south Eurasian margin occurred during the Mesozoic in an accretionary‐type tectonic setting through successive episodes of terrane accretions and probably slab flattening, transiently increasing the coupling at the plate interface. Moreover, they indicate that the south Eurasian margin was already hot and thickened prior to Cenozoic collision with India, which has important implications for orogen‐scale strain‐accommodation mechanisms.     

Temporal and spatial variation of infiltration in urban green infrastructure

Infiltration is the primary mechanism in green stormwater infrastructure (GSI) systems to reduce the runoff volume from urbanized areas. Soil hydraulic conductivity is most important in influencing GSI infiltration rates. Saturated hydraulic conductivity (Ksat) is a critical parameter for GSI design and post-construction performance. However, Ksat measurement in the field is problematic due to temporal and spatial variability and measurement errors. This review paper focuses on a comparison of methods for in-situ Ksat measurement and the causes of temporal and spatial variations of Ksat within GSI systems. Automated infiltration testing methods, such as the Modified Philip–Dunne (MPD) and SATURO infiltrometers, show promise for efficient Ksat measurements. Soil Ksat values can change over time and substantially vary throughout a GSI, which can be attributed to multiple factors, including but not limited to temperature changes, soil composition and properties, soil compaction level, plant root morphology and distribution, biological and macrofauna activities in the soil, inflow sediment characteristics, quality of infiltrating water, and measurement errors. There is evidence that infiltration rates in vegetated urban GSI systems are sustained given an appropriate GSI design, reasonable concentration of suspended sediments in the inflow runoff, and routine maintenance procedures. These observations indicate that clogging can be counteracted by processes that tend to increase the soil hydraulic conductivity (e.g., plant root and biological activities). This self-sustainability underlines that infiltration-based GSI systems are a reliable long-term stormwater management solution. Recommendations on how to incorporate the temporal changes of Ksat in GSI design and on obtaining a spatially-representative Ksat for the GSI design are presented.     

Probability of Peak Ground Horizontal and Peak Ground Vertical Accelerations at Tehran and Surrounding Areas

The active faults near Tehran are capable of producing Mw magnitudes of 6.62?C7.23; at epicentral distances of 25?km from the active faults, and Mw?=?7.23, the peak ground horizontal acceleration, PGH, is between 386 and 730?cm/s² and peak ground vertical acceleration, PGV, is between 192 and 261?cm/s²; the historic earthquake of the fourth century BC, Mw?=?7.16 produced the highest estimated PGH acceleration in Tehran, between 254 and 479?cm/s² and PGV acceleration between 127 and 173?cm/s². Earthquakes from 1909 to 2008, within 300?km from Tehran, are used for calculation of magnitude frequency relation, and results are applied to estimate PGH and PGV accelerations as a function of magnitudes for a set of fixed epicentral distance and site conditions; also as a function of epicentral distances for a set of fixed magnitudes and site conditions. Poisson??s distribution is used for probability calculation of PGH and PGV accelerations for several exposure times, site conditions and epicentral distances; accelerations with very high probability, near 1, are presented. At an epicentral distance of 10?km and exposure time of 450?years, in the northern part of Tehran, close to Mosha and the Northern Tehran faults, PGH acceleration is 800?C420?cm/s² and PGV is 400?C220?cm/s² with high probability. At an epicentral distance of 25?km and 1,000?years exposure time, PGH is 610?C320?cm/s² and PGV is 310?C160?cm/s² with high probability, where larger values are for soft soil and smaller values are for hard rock.     

塔里木盆地大气降尘变化特征及影响因素分析

利用哈密、塔中与和田沙尘暴观测站2005—2008年大气降尘观测资料,结合铁干里克、民丰等18个大气降尘监测点2007年4月以来的大气降尘、沙尘天气等观测资料,分析了塔里木盆地大气降尘变化特征,同时分析了影响大气降尘浓度变化的主要因素。结果表明：①2005—2008年3站中哈密大气降尘量最少,其次为和田,塔中的大气降尘量最多; 2006年是4 a中降尘量最多年份,其次为2005年,2008年相对较少。②春夏季是大气降尘集中季节,秋季略大于冬季。塔里木盆地中部最高,南部及西南部的站点降尘明显高于盆地的东部和西部。③在沙尘暴季节,无论是5月还是6月,2007年大气降尘量基本上都高于2008年。④2007年5月和6月,除盆地东面的哈密、铁干里克以及西缘的阿拉尔3站外,其他18个站沙尘天气较多。沙尘天气是塔里木盆地大气降尘量高的主要因素,沙尘天气越多,大气降尘量越高。