The Mw = 6.0, 7 September 1999 Athens Earthquake

On 7 September 1999 at 11:56 GMT a destructive earthquake (M_w = 6.0) occurred close to Athens (Greece). The rupture process is examined using data from the Cornet local permanent network, as well as teleseismic recordings. Data recorded by a temporary seismological network were analyzed to study the aftershock sequence. The mainshock was relocated at 38.105°N, 23.565°E, about 20 km northwest of Athens. Four foreshocks were also relocated close to the mainshock. The modeling of teleseismic P and SH waves provides a well-constrained focal mechanism of the mainshock (strike = 105°, dip = 55° and rake = -80°) at a depth of 8 km and a seismic moment M₀ = 1.010²⁵ dyn·cm. The obtained fault plane solution represents normal faulting indicating an almost north-south extension. More than 3500 aftershocks were located, 1813 of which present RMS < 0.1 s and ERH, ERZ < 1.0 km. Two main clusters were distinguished, while the depth distribution is concentrated between 2 and 11 km. Over 1000 fault plane solutions of aftershocks were constrained, the majority of which also correspond to N–S extension. No surface breaks were observed but the fault plane solution of the mainshock is in agreement with the tectonics of the area and with the focal mechanisms obtained by aftershocks. The hypocenter of the mainshock is located on the deep western edge of the fault plane. The relocated epicenter coincides with the fringe that represents the highest deformation observed on the differential interferometric image. The calculated source duration is 5 sec, while the estimated dimensions of the fault are 15 km length and 10 km width. The source process is characterized by unilateral eastward rupture propagation, towards the city of Athens. An evident stop phase observed in the recordings of the Cornet local stations is interpreted as a barrier caused by the Aegaleo Mountain.     

Large-scale mass wasting in ancient volcanic materials

Giant landslides are significant hazards associated with many active volcanic edifices. We describe a similar feature on ancient (>4 Ma) volcanic deposits subject to active tectonism. The landslide is approximately 3 km long by 1 km wide, with an estimated depth of 400 m. Side margins are straight and parallel, mimicking regional structure; narrow valleys incised down these margins provide low-strength side-release surfaces. Between these is a giant slump consisting of at least four, largely intact, discrete blocks that have moved down-dip a distance of >500 m. A series of flows with areal extents ranging from 0.01 to 0.5 km² extends from the front of the failure. The materials represent an eroded sequence of andesite flows on the flanks of a stratovolcano. These have undergone two phases of hydrothermal alteration, and are deeply weathered to low-density (1040±80 kg m⁻³) silt (59%) and clay (35%) materials with strength properties typical of weathered silts (c=26±3 kN m⁻²; φ=42±8°). The size and location of this landslide preclude detailed geotechnical investigation of the failure. The worth of numerical stability analysis as an alternative technique in assessing the nature of the failure and hence the risk it poses to nearby communities is investigated. Sensitivity analysis identified likely conditions under which initial failure may have occurred: analyses for sensitivity to strength and earthquake acceleration needed conversion to critical combinations (F=1.0) of water table and strength/acceleration to remove the overriding influence of water table fluctuations. Failure was likely initiated either by a high water table level (83-84%), or some combination of intensity VII-IX earthquake waves together with water table heights of 40-80%. A general hazard assessment indicates that the risk associated with creep and catastrophic failure of the main mass is small, whereas the risk from flow failures near the toe of the landslide may be high. Important parameters (hydrological regime, flow failure morphology, age of initiation, and rates of movement) requiring closer investigation are identified. Development of a model is crucial to assessing the hazard associated with a feature such as that described here. With limited resources, a detailed stability analysis is a powerful tool as an initial stage in hazard analysis.     

内蒙古宝音图岩群变质基性火山岩锆石U-Pb年龄及意义

内蒙古达茂旗宝音图岩群是在达茂旗、查干呼绍地区 1:5万区调和白云鄂博幅 1:2 5万区调过程中从过去的巴特敖包群和白云鄂博群未分岩组解体出来的中级变质岩系。通过对宝音图岩群中阳起钠长片岩、二云片岩 (原岩均为基性火山岩 )进行单颗粒锆石U -Pb法年龄测定 ,获得了两条不一致线上交点年龄 (2 486± 42 )Ma、(2 496± 2 6 )Ma ,代表了宝音图岩群基性火山岩及宝音图岩群的形成时代 ;在不一致线之外不同成因的锆石的2 0 7Pb/ 2 0 6Pb表面年龄 (2 2 2 7± 15 )Ma ,寓示着 2 2 0 0Ma左右发生了一次区域性变质热事件 ;下交点年龄 (319±486 )Ma和 (40 6± 46 )Ma ,代表了 2 5 0～ 480Ma时期该区大面积岩浆侵位及宝音图岩群后期变形变质改造的时期 ,两个年龄样的获得为研究古元古代古陆的裂解及白云鄂博海槽的形成提供了直接依据。     

圆弧状沉积盆地与软土单覆盖层出平面地表运动对比

采用解析方法，通过圆弧状沉积盆地与同样参数的软土均匀单覆盖层的对比，研究沉积盆地对平面SH波二维散射效应和其出平面地表运动。数值对比结果表明：（1）波垂直入射时盆地中央地表位移幅值反应，其随频率变化形式与软土单覆盖层的地表运动有较强的对应性，且此时单覆盖层地表位移幅值为其下降；（2）盆地的二维散射效应特点显著波的汇聚作用使盆地中的地表运动显著放大，总体上较软土单覆盖层的地表位移幅值要大，其频域内的最大反应则明显大于单覆盖层情况，但最大反应在盆地表面出现的位置依赖于入射波角度、盆地深宽比和盆地的软硬程度；（3）沉积盆地中地表位移运动有强烈扭转效应，且随盆地变软和入射波频率增高而增强。     

2-D crustal structure of Changbaishan-Tianchi volcanic region determined by seismic traveltime inversion

IntroductionTrial-and-error forward modeling of wide-angle seismic reflection/refraction traveltimes for 2-D velocity structure is extremely time-consuming, even for experienced data interpreters. For wide-angle seismic reflection/refraction experiments that consist of numerous shots along a single line, it is quite difficult through repeated trial-and-error forward modeling to construct a 2-D model that fits the data within acceptable limits (Cerveny, et al, 1977; ZHANG, et al, 200 . In ad…     

Transient foundation solution of saturated soil to impulsive concentrated loading

The transient dynamic response of saturated soil under suddenly applied normal and horizontal concentrated loading is studied in this paper. The behavior of saturated soil is governed by Biot's consolidation theory. The general solutions for Biot equations of equilibrium are derived in terms of displacements and variations of fluid volume, using Laplace–Hankel integral transforms. The solutions in the time domain can be evaluated by numerical inverse Laplace–Hankel transforms. Selected numerical results for displacements, stresses, and pore pressures are presented. Comparisons with existing closed-form solutions for the elastic half-space are made to confirm the accuracy of the present solutions. The solutions can be used to study a variety of transient wave propagation problems and dynamical interactions between saturated soil and structures.     

A methodology for the evaluation of long-term volcanic risk from pyroclastic flows in Campi Flegrei (Italy)

The volcanological history of Campi Flegrei suggests that the most frequent eruptions are characterized by the emplacement of pyroclastic flow and surge deposits erupted from different vents scattered over a 150-km² caldera. The evaluation of volcanic risk in volcanic fields is complex because of the lack of a central vent. To approach this problem, we subdivided the entire area of Campi Flegrei into a regular grid and evaluated the relative spatial probability of opening of vents based on geological, geophysical and geochemical data. We evaluated the volcanic risk caused by pyroclastic flows based on the formula proposed by UNESCO (1972), R=H×V×Va, where H is the hazard, V is the vulnerability and Va is the value of the elements at risk. The product H×V was obtained by performing simulations of type eruptions centered in each cell of the grid. The simulation is based on the energy cone scheme proposed by Sheridan and Malin [J. Volcanol. Geotherm. Res. 17 (1983) 187–202], hypothesizing a column collapse height of 100 m for eruptions of VEI=3 and 300 m for eruptions of VEI=4 with a slope angle of 6°. Each simulation has been given the relative probability value associated with the corresponding cell. We made use of the GIS software ArcView 3.2 to evaluate the intersection between the energy cone and the topography. The superposition of the areas invaded by pyroclastic flows (124 simulations for VEI=3 and 37 for VEI=4) was used to obtain the relative hazard map of the area. The relative volcanic risk map is obtained by superimposing the urbanization maps.     

Gravity and magnetic signatures of volcanic plugs related to Deccan volcanism in Saurashtra, India and their physical and geochemical properties

The Bouguer anomaly and the total intensity magnetic maps of Saurashtra have delineated six circular gravity highs and magnetic anomalies of 40-60 mGal (10⁻⁵m/s²) and 800-1000 nT, respectively. Three of them in western Saurashtra coincide with known volcanic plugs associated with Deccan Volcanic Province (DVP), while the other three in SE Saurashtra coincide with rather concealed plugs exposed partially. The DVP represents different phases of eruption during 65.5±2.5 Ma from the Reunion plume. The geochemical data of the exposed rock samples from these plugs exhibit a wide variation in source composition, which varies from ultramafic/mafic to felsic composition of volcanic plugs in western Saurashtra and an alkaline composition for those in SE Saurashtra. Detailed studies of granophyres and alkaline rocks from these volcanic plugs reveal a calc-alkaline differentiation trend and a continental tectonic setting of emplacement. The alkaline plugs of SE Saurashtra are associated with NE-SW oriented structural trends, related to the Gulf of Cambay and the Cambay rift basin along the track of the Reunion plume. This indicates a deeper source for these plugs compared to those in the western part and may represent the primary source magma. The Junagadh plug with well differentiated ring complexes in western Saurashtra shows well defined centers of magnetic anomaly while the magnetic anomalies due to other plugs are diffused though of the same amplitude. This implies that other plugs are also associated with mafic/ultramafic components, which may not be differentiated and may be present at subsurface levels. Paleomagnetic measurements on surface rock samples from DVP in Saurashtra suggest a susceptibility of 5.5×10⁻² SI units with an average Koenigsberger ratio (Q_n) of almost one and average direction of remanent magnetization of D=147.4° and I=+56.1°. The virtual geomagnetic pole (VGP) position computed from the mean direction of magnetization for the volcanic plugs and Deccan basalt of Saurashtra is 30°N and 74°W, which is close to the VGP position corresponding to the early phases of Deccan eruption. Modeling of gravity and magnetic anomalies along two representative profiles across Junagadh and Barda volcanic plugs suggest a bulk density of 2900 and 2880 kg/m³, respectively and susceptibility of 3.14×10⁻² SI units with a Q_n ratio of 0.56 which are within the range of their values obtained from laboratory measurements on exposed rock samples. The same order of gravity and magnetic anomalies observed over the volcanic plugs of Saurashtra indicates almost similar bulk physical properties for them. The inferred directions of magnetization from magnetic anomalies, however, are D=337° and 340° and I=−38° and −50° which represent the bulk direction of magnetization and also indicate a reversal of the magnetic field during the eruption of these plugs. Some of these plugs are associated with seismic activities of magnitude ≤4 at their contacts. Based on this analysis, other circular/semi-circular gravity highs of NW India can be qualitatively attributed to similar subsurface volcanic plugs.