Geophysical methods applied to fault characterization and earthquake potential assessment in the Lower Tagus Valley, Portugal

The study region is located in the Lower Tagus Valley, central Portugal, and includes a large portion of the densely populated area of Lisbon. It is characterized by a moderate seismicity with a diffuse pattern, with historical earthquakes causing many casualties, serious damage and economic losses. Occurrence of earthquakes in the area indicates the presence of seismogenic structures at depth that are deficiently known due to a thick Cenozoic sedimentary cover. The hidden character of many of the faults in the Lower Tagus Valley requires the use of indirect methodologies for their study. This paper focuses on the application of high-resolution seismic reflection method for the detection of near-surface faulting on two major tectonic structures that are hidden under the recent alluvial cover of the Tagus Valley, and that have been recognized on deep oil-industry seismic reflection profiles and/or inferred from the surface geology. These are a WNW–ESE-trending fault zone located within the Lower Tagus Cenozoic basin, across the Tagus River estuary (Porto Alto fault), and a NNE–SSW-trending reverse fault zone that borders the Cenozoic Basin at the W (Vila Franca de Xira–Lisbon fault). Vertical electrical soundings were also acquired over the seismic profiles and the refraction interpretation of the reflection data was carried out. According to the interpretation of the collected data, a complex fault pattern disrupts the near surface (first 400 m) at Porto Alto, affecting the Upper Neogene and (at least for one fault) the Quaternary, with a normal offset component. The consistency with the previous oil-industry profiles interpretation supports the location and geometry of this fault zone. Concerning the second structure, two major faults were detected north of Vila Franca de Xira, supporting the extension of the Vila Franca de Xira–Lisbon fault zone northwards. One of these faults presents a reverse geometry apparently displacing Holocene alluvium. Vertical offsets of the Holocene sediments detected in the studied geophysical data of Porto Alto and Vila Franca de Xira–Lisbon faults imply minimum slip rates of 0.15–0.30 mm/year, three times larger than previously inferred for active faults in the Lower Tagus Valley and maximum estimates of average return periods of 2000–5000 years for M 6.5–7 co-seismic ruptures.     

A way to synchronize models with seismic faults for earthquake forecasting: Insights from a simple stochastic model

Numerical models are starting to be used for determining the future behaviour of seismic faults and fault networks. Their final goal would be to forecast future large earthquakes. In order to use them for this task, it is necessary to synchronize each model with the current status of the actual fault or fault network it simulates (just as, for example, meteorologists synchronize their models with the atmosphere by incorporating current atmospheric data in them). However, lithospheric dynamics is largely unobservable: important parameters cannot (or can rarely) be measured in Nature. Earthquakes, though, provide indirect but measurable clues of the stress and strain status in the lithosphere, which should be helpful for the synchronization of the models.The rupture area is one of the measurable parameters of earthquakes. Here we explore how it can be used to at least synchronize fault models between themselves and forecast synthetic earthquakes. Our purpose here is to forecast synthetic earthquakes in a simple but stochastic (random) fault model. By imposing the rupture area of the synthetic earthquakes of this model on other models, the latter become partially synchronized with the first one. We use these partially synchronized models to successfully forecast most of the largest earthquakes generated by the first model. This forecasting strategy outperforms others that only take into account the earthquake series. Our results suggest that probably a good way to synchronize more detailed models with real faults is to force them to reproduce the sequence of previous earthquake ruptures on the faults. This hypothesis could be tested in the future with more detailed models and actual seismic data.     

地震快报系统

本文简述了开发地震快报系统的主要技术。     

三种高程传递方法的精度分析

通过试验对悬挂钢尺、GPS测高和全站仪垂直测高三种高程传递方法进行了理论分析与精度比较,并讨论了各种方法在实际应用中的优点与局限性。     

中国南极长城站1998年大气降水化学特征的初步研究

文中报告了 1998年 1～ 12月 ,在中国南极长城站 (6 2°13′S ,5 8°5 8′W ,海拔 10 .0m)采集的 115个有效降水样品的 pH值、电导率和化学组分分析结果。长城站地区降水的年平均 pH值为 5 .6 2 ,电导率年平均值为 85 .16 μS/cm。秋季期间降水的 pH值和电导率较高 ,其它季节较低。降水中最主要的离子为Cl-,Na+ ,其他离子按平均浓度值依次排列为SO42 -,Mg2 + ,Ca2 + ,K+ ,NO3 -,NH4+ 。除NO3 -,NH4+ 外 ,各主要离子浓度间呈显著正相关 ,表明其来源具有同一性。降水样品的离子组成比例与周围海区的海水接近 ,显示了海盐是降水中可溶性离子的主要来源 ,而其它源只对Ca2 + 有一定贡献。该地区的降水具有较典型的海洋性降水的化学特征     

中朝板块元古宙板内地震带与盆地格局

地史中发生的强地震事件在地层中留下固定的记录 (图 1～ 3) ,这些记录在区域上呈带状分布 ,代表地史中的地震带。中朝板块元古宙目前可识别出两个板内地震带 (图 5 )。中元古代板内地震带 (170 0～ 12 0 0Ma)西起太行山北段 ,经燕山山脉、辽宁西部、穿越辽河平原至辽宁北部的泛河流域分布 ,即燕山—泛河地震带 ,现今呈NEE向延伸。新元古代震旦纪地震带沿吉林南部、辽东半岛、山东中部及苏皖北部现今呈NNE走向分布 ,即古郯庐地震带 (6 5 0～ 6 0 0Ma)。上述两个板内地震带是元古宙不同时期超大陆裂解的响应。中元古代与新元古代两个不同方向的地震断裂带分别控制着两个时期的盆地边界。燕山泛河地震断裂带构成中元古代海盆南界 (指现在的位置 ) ,形成向北开放的海域。古郯庐地震断裂带将中朝板块裂解为华北块体与胶辽朝块体。古郯庐地震断裂带构成震旦纪海域的边界 ,震旦纪海盆通过朝鲜半岛与当时的外海相连接 ,华北块体则为陆源剥蚀区。文内四幅古地理图 (图 6～ 9)是以地震灾变思想为指导 ,以新的地层研究、对比为基础编制的 ,侧重反映了盆地的格局及其变化。根据地震、同沉积断裂新的思路 ,可提供地质学家重新认识与解释某些沉积矿床的成因 ,它们的成矿元素均来自地球深部而非地表风化作用。文中编制     

Structure and tectonic evolution of the Southern Eurasia Basin, Arctic Ocean

Multichannel seismic reflection data acquired by Marine Arctic Geological Expedition (MAGE) of Murmansk, Russia in 1990 provide the first view of the geological structure of the Arctic region between 77–80°N and 115–133°E, where the Eurasia Basin of the Arctic Ocean adjoins the passive-transform continental margin of the Laptev Sea. South of 80°N, the oceanic basement of the Eurasia Basin and continental basement of the Laptev Sea outer margin are covered by 1.5 to 8 km of sediments. Two structural sequences are distinguished in the sedimentary cover within the Laptev Sea outer margin and at the continent/ocean crust transition: the lower rift sequence, including mostly Upper Cretaceous to Lower Paleocene deposits, and the upper post-rift sequence, consisting of Cenozoic sediments. In the adjoining Eurasia Basin of the Arctic Ocean, the Cenozoic post-rift sequence consists of a few sedimentary successions deposited by several submarine fans. Based on the multichannel seismic reflection data, the structural pattern was determined and an isopach map of the sedimentary cover and tectonic zoning map were constructed. A location of the continent/ocean crust transition is tentatively defined. A buried continuation of the mid-ocean Gakkel Ridge is also detected. This study suggests that south of 78.5°N there was the cessation in the tectonic activity of the Gakkel Ridge Rift from 33–30 until 3–1 Ma and there was no sea-floor spreading in the southernmost part of the Eurasia Basin during the last 30–33 m.y. South of 78.5°N all oceanic crust of the Eurasia Basin near the continental margin of the Laptev Sea was formed from 56 to 33–30 Ma.     

Interpretation of seismic and potential field data from western New York State and Lake Ontario

Lithoprobe and industry seismic profiles have furnished evidence of major zones of easterly dipping Grenville deformed crust extending southwest from exposed Grenville rocks north of Lake Ontario. Additional constraints on subsurface structure limited to the postulated Clarendon–Linden fault system south of Lake Ontario are provided by five east–west reflection lines recorded in 1976. Spatial correlations between seismic structure and magnetic anomalies are described from both Lake Ontario and the newly reprocessed New York lines.In the Paleozoic to Precambrian upper crust, the New York seismic sections show: (1) An easterly thickening wedge of subhorizontal Paleozoic strata unconformably overlying a Precambrian basement whose surface has an apparent regional easterly dip of 1–2°. Minor apparent normal offsets, possibly on the order of tens of meters, occur within the Paleozoic section. The generally poorly reflective unconformity may be locally characterized by topographic relief on the order of 100 m; (2) Apparent local displacement on the order of 90 m at the level of the Black River Group diminishes upward to little or no apparent offset of Queenston Shale; (3) Within the limited seismic sections, there appears to be no evidence that the complete upper crustal section is vertically or subvertically offset; (4) Dipping structure in the Paleozoic strata (15° to 35°) resembles some underlying Precambrian basement elements; (5) The surface continuity of inferred faults constituting the Clarendon–Linden system is not strongly supported by the seismic data.Beneath the Paleozoic strata, the seismic sections show both linear and arcuate reflector geometry with easterly apparent dips of 15° to 35° similar to the deep structures imaged on seismic lines from nearby Lake Ontario and on Lithoprobe lines to the north. The similarity supports an extension of easterly dipping Central Metasedimentary Belt structures of the Grenville orogen from southern Ontario to beneath western New York State.From a comparison of the magnetic and gravity fields with the New York seismic sections, we suggest: (1) The largely nonmagnetic Paleozoic strata appear to contribute negligibly to magnetic anomalies. Seismically imaged fractures in the New York Paleozoic strata appear to lie mainly west of a positive gravity anomaly. The relationship between magnetic and gravity anomalies and the changes in the geometry of interpreted Precambrian structures remains enigmatic; (2) North to northeast trending curvilinear magnetic and gravity anomalies parallel, but are not restricted to the principal trend of the postulated Clarendon–Linden fault system. Paleozoic fractures of the Clarendon–Linden system may partly overlie a southward extension of the Composite Arc Belt boundary zone.     

Physical properties of ultrahigh-pressure metamorphic rocks from the Sulu terrain, eastern central China: implications for the seismic structure at the Donghai (CCSD) drilling site

About 30 samples representing major lithologies of Sulu ultrahigh-pressure (UHP) metamorphic rocks were collected from surface exposures and exploration wells, and compressional (Vp) and shear wave (Vs) velocities and their directional dependence (anisotropy) were determined over a range of constant confining pressures up to 600 MPa and temperatures ranging from 20 to 600 °C. Samples range in composition from acidic to ultramafic. P- and S-wave velocities measured at 600 MPa vary from 5.08 to 8.64 km/s and 2.34 to 4.93 km/s, respectively. Densities are in the range from 2.60 to 3.68 g/cm³. To make a direct tie between seismic measurements (refraction and reflection) and subsurface lithologies, the experimental velocity data (corresponding to shallow depths) were used to calculate velocity profiles for the different lithologies and profiles of reflection coefficients at possible lithologic interfaces across the projected 5000-m Chinese Continental Scientific Drilling Program (CCSD) crustal segment. Comparison of calculated in situ velocities with respective intrinsic velocities suggests that the in situ velocities at shallow depths are lowered by an increased abundance of open microcracks. The strongly reflective zone beneath the Donghai drill site can be explained by the impedance contrasts between the different lithologies. Contacts between eclogite/peridotite and felsic rocks (gt-gneiss, granitic gneiss), in particular, may give rise to strong seismic reflections. In addition, shear-induced (lattice preferred orientation (LPO)-related) seismic anisotropy can increase reflectivity. For the explanation of the high velocity bodies (>6.4 km/s) around 1000 m and below 3200-m depth, large proportions of eclogite/peridotite (about 40 and 30 vol.%, respectively) are needed.     

Stability of earthquake moment tensor inversions: effect of the double-couple constraint

We show that spurious large non-double-couple components can be obtained in inversions for the full deviatoric moment tensor for shallow crustal earthquakes due to inaccurate Earth models. The traditional “best double-couple” solution does not in general provide an optimal estimate of a double-couple mechanism, and is only reliable when the non-double-couple component of the full deviatoric solution is small. The inverse problem for the moment tensors of the 1998 Antarctic Plate and 2000 Wharton Basin strike-slip earthquakes is shown in each case to have two well-fitting minima in the misfit function of pure double-couple solutions. Such pairs of solutions are most likely to exist for earthquakes which are close either to vertical strike-slip or to dip-slip on a fault plane dipping at 45°. It is shown theoretically that these pairs of solutions arise from the combination of the pure double-couple constraint and the instability of two elements of the moment tensor. No significant non-double-couple component is found for the shallow thrusting 1996 Biak, Indonesia earthquake.