The magma budget of Volcan Arenal, Costa Rica from 1968 to 1980

The volume of magma emitted by Volcan Arenal from July 1968 to March 1980 has been calculated to be 304 × 10⁶ m³ (dense rock equivalent). Most of this magma has been emplaced as block lava flows on the western flanks of the volcano following the initial explosive eruptions in 1968. From 1968 to 1973 the volumetric discharge rate of magma decreased from about 3-2 m³ s⁻¹ to about 1 m³ s⁻¹. During a break in activity in late 1973 the site of effusion moved from Crater A to Crater C about 400 m higher. Subsequent effusion was at a lower rate (0.3 m³ s⁻¹) which remained constant for the next six years. Comparison of dry-tilt measurements during this latter period of steady-state effusion with numerical finite-element models of Arenal's elastic response to the evacuation of magma from an underlying reservoir favor a very shallow reservoir (< 2 km depth) to explain the data. However, the constraints imposed by the measured volumes of magma are not compatible with such a reservoir. Instead, it is argued that the steady downward tilting of the volcano's summit was caused by the loading of the western side of the volcano by about 19 × 10⁶ m³ of lava. Surface loading by lava flows may be an important deformational effect at other volcanoes. A system of magma supply involving open conduits (pipes) for the uppermost one kilometer and transitory conduits (cracks) to a crustal reservoir is proposed. This crustal reservoir initially contained a compositionally graded magma which was evacuated from 1968 to 1973. The subsequent abrupt decrease in effusion rate is compatible with the increased magmatic head required to reach Crater C. The constancy of magma composition and effusion rate from 1974 to 1980 implies a homogeneous magma reservoir.     

Intermediate and deep earthquakes in Spain

Recent improvements in the seismological networks on the Ibero-Maghrebian region have permitted estimation of hypocentral location and focal mechanisms for earthquakes which occurred at South Spain, Alboran Sea and northern Morocco of deep and intermediate depth, with magnitudes between 3.5 and 4.5. Intermediate depth shocks, range from 60 to 100 km, with greater concentration located between Granada and Málaga. Fault-plane solutions of 5 intermediate shocks have been determined; they present a vertical plane in NE-SW or E-W direction. Seismic moments of about 10¹⁵ Nm and dimensions of about 1 km have been determined from digital records of Spanish stations.P-wave forms are complex. This may be explained by the crustal structure near the station, discontinuities in the upper mantle and inhomogeneities near the source. Deep activity at about 650 km has only 3 shocks since 1954 (1954, 1973, 1990). Shocks are located at a very small region. Fault-plane solutions show a consistent direction of the pressure axis dipping 45° in E direction. For the 1990 shock seismic moment is 10¹⁶ Nm and dimensions 2.6 km. TheP-waves are of simpler form with a single pulse. The intermediate and deep activities are not connected and no activity has been detected between 100 and 650 km. The intermediate shocks may be explained in terms of a recent subduction from Africa under Iberia in SE direction. The very deep activity must be related to a sunk detached block of lithospheric material still sufficiently cold and rigid to generate earthquakes.     

Nemuro-Oki earthquake of June 17, 1973: A lithospheric rebound at the upper half of the interface

P-wave first motions, radiation patterns and amplitudes of long-period surface waves, relocated aftershock distributions, leveling and tsunami data indicate that the 1973 Nemuro-Oki earthquake is caused by a low-angle thrust-faulting, representing a rebound at the upper 50 km of the interface between the continental and oceanic lithospheres. Rebound, most likely aseismic, at depths below 50 km, is suggested to take place in the near future from a comparison of recent geologic crustal deformation with pre-seismic and co-seismic data. The estimated seismic moment is about $\text{1}\text{3}\text{–}\text{1}\text{4}$ of that for the neighboring great earthquakes. The macro-seismic data suggest that the 1973 earthquake is smaller than the 1894 Nemuro-Oki earthquake, the last great earthquake in this region.The 1973 earthquake had been predicted on the basis of a seismic gap. Although the prediction was successful as to the location and nature of the faulting and partly as to the occurrence time, it is smaller than the predicted one. A part of the seismic gap may still remain. The difference between the observed seismic slip (1.6 m) and that predicted on the basis of the pre-seismic crustal deformation (3.0 m) indicates either (1) the 1973 earthquake relieved only a part of the strain accumulated in the upper 50 km, or (2) a significant amount of aseismic slip took place on the seismic fault and completely relieved the accumulated strain in the focal region of the 1973 earthquake. If the former is the case, the remaining strain, not only in the focal region, but also in the remaining seismic gap adjoining it, may be relieved in a larger earthquake in the future.The source parameters obtained are as follows: fault plane, dip direction = N40°W, dip angle = 27°; seismic moment = 6.7 · 10²⁷ dyn cm; average slip dislocation, 1.6 m in N63°W direction; stress drop = 35 bars. In these calculations, the fault dimension and the rigidity are assumed to be 100 · 60 km² and 7.0 · 10¹¹ dyn/cm², respectively.     

中国南北地震带北段及其周缘地壳厚度与泊松比研究

本文利用中国地震科学探测台阵2013—2015年在南北地震带北段及其周缘架设的673个台站所记录到的远震波形所提取到的接收函数并应用H-κ扫描方法获取了南北地震带北段及其周缘的地壳厚度和泊松比,结果显示研究区地壳厚度从青藏高原东北缘向鄂尔多斯块体逐渐减小,从65 km逐渐减薄至40 km,不同块体之间地壳厚度存在明显差异.祁连造山带西部地壳厚度超过60 km,而东部地壳厚度仅为约50 km左右,表明祁连造山带东、西部地壳增厚变形存在着明显差异.西秦岭造山带地壳厚度从60 km减薄到40 km,其东部具有较薄的地壳厚度可能经历了拆沉.阿拉善块体作为华北克拉通西部块体的一部分,西部地壳厚度约50 km,而东部约45 km,表明阿拉善块体西部由于印度—欧亚板块碰撞也受到了活化改造,其克拉通性质只在其中东部残留.研究区泊松比变化范围为0.20~0.31,平均泊松比约0.25,表明地壳主要由长英质矿物组成,较高的泊松比主要分布在六盘山断裂带和银川—河套地堑.研究结果显示地壳厚度与高程之间具有较好的相关性,表明地壳整体上处于相对均衡的状态,而西秦岭造山带和祁连造山带东部的部分区域地壳可能处于不均衡状态.

     

A possible case of SES selectivity at Kozu-shima Island, Japan

Geoelectric potential difference monitoring has been conducted in Kozu-shima Island about 170 km south of Tokyo. During the monitoring period (14 May 1997–30 November 1999), there were 21 anomalous geoelectric changes that could have been VAN-type seismic electric signals (SES) and 26 M>3 earthquakes within about 20 km. Possible correspondence between the geoelectric changes and earthquakes has been examined. There seems to be some systematic tendency regarding the dependence of appearance/non-appearance of changes and their polarity on epicentral location. This observation, however, is preliminary and admittedly inconclusive because of small number of cases and non-uniqueness in signal-earthquake correlation.     

Features of deep crustal structure and the onshore-offshore transition at the Iberian flank of the Valencia Trough (Western Mediterranean)

First results are presented of a recent onshore seismic survey complementary to the Valsis-2 Cruise, which consisted of ESP, COP and CDP marine seismic profiles across the Valencia Trough (Western Mediterranean).The marine energy source used was an airgun array of 5800 cubic inch recorded at 2 land stations on the western flank of the Valencia Trough, at distances between 10–120 km.The experiment has resulted in an extended sampling of the deep crustal structure of the eastern Mediterranean flank of the Iberian peninsula, as well as the offshore-onshore transition.Three transverse NW-SE profiles have been interpreted. Local thinning of the sedimentary cover has been determined towards the centre of the basin which, together with the shallow high velocities observed on the southern profile, could be related to volcanic episodes.A seismic continental basement has been found at depths between 3 and 5 km. A thin lower crust (3–5 km) with velocities around 6.8 km/s has been identified in the northern part of the basin. Alternative crustal models considered for the 3 profiles have been tested, not only from arrival times but also from relative amplitude distributions. A first-order Moho discontinuity fits the data best. The welldefined Moho boundary results in energetic P_MP reflections, and a clear updoming is observed towards the interior of the basin, from depths about 20–21 km inshore of Barcelona to 15–17 km depths 60 km offshore. An anomalous upper mantle with low P_n velocities of about 7.7 km/s is confirmed in most of the sampled areas.     

Crustal structure in the Mizuho Plateau, East Antarctica, by a two-dimensional ray approximation

Studies on the crustal and upper-mantle structure in Antarctica have been one of the major contributions to Antarctic research since the International Geophysical Year of 1957–1958. Many refraction surveys with small charge size have been conducted in Antarctica, but long-range experiments were also made in 5 regions on the margin of the Antarctic continent.In 1979–1981, the scientific program of the Japanese Antarctic Research Expedition was focused on the earth sciences, and in particular, an explosion seismic experiment along a long survey line was the major item during these years. An experiment along a 300 km-long line with three shots and 27 observation stations was successfully made in the northern Mizuho Plateau, East Antarctica.From the analysis of travel times and the amplitude study of synthetic seismograms, the crustal structure of the northern Mizuho Plateau was determined. The depths of the Conrad and the Moho discontinuties were determined as 31 km and 42 km, respectively. The P-wave velocity and depth relation was determined as 6.0 km/s on the surface of the bedrock, 6.3 km/s at a depth of 2 km, 6.8 km/s at the Conrad and 7.9 km/s at the Moho. The velocity in the crust increases gradually. The crustal structure obtained is representative of East Antarctica.     

Seismic investigations of crustal structure in east-central Minnesota

Relative arrival times from 120 mine blasts recorded on the Central Minnesota Seismic Array (CMSA) out to distances of 350 km have been statistically analysed and interpreted to yield a crustal velocity-depth function for the Minnesota area. Velocities increase continuously from 5.9 km s^?1 at the surface to 7.4 km s^?1 at 40 km depth, with a steepening of the velocity gradient at about 20 km. A long-range refraction profile extending from the Mesabi Iron Range in Northern Minnesota to the seismic array, has also been recorded. A striking feature of the profile is the absence of direct P waves in the distance range 60–120 km. Our preferred explanation for this shadow zone is a fault or fracture zone near x = 60 km. Reflections within the shadow zone have been analysed to give a Moho depth of 42 km under northeastern Minnesota. Teleseismic residuals have been computed for 85 earthquakes recorded on the CMSA. No systematic dependence of the residuals upon azimuth or distance was found. The average residuals agree quantitatively with time delays predicted from earlier upper crustal refraction studies (Mooney et al., 1970).     

腾冲火山区的地壳厚度和平均泊松比研究

腾冲是青藏高原东南缘重要的第四纪火山活动区域,全新世以来的火山主要集中在腾冲盆地的中央,由北向南形成一个串珠状的火山链.为了深入探索这一火山区的深部结构和岩浆活动特征,我们在腾冲北部开展了为期一年的流动地震观测,利用接收函数方法计算了台站下方的地壳厚度、平均波速比和泊松比,研究结果揭示出测线下方地壳结构与岩浆活动及火山分布的对应关系.测线北部7个台站的地壳厚度在35.4~37.6 km之间,平均波速比为1.82~1.92、泊松比为0.28~0.31,其中马站附近莫霍面抬升幅度最大,与相邻地区莫霍面深度相差1~2 km,平均波速比和泊松比也达到最大值.相比之下,测线南端两个台站的地壳厚度接近40 km,平均波速比和泊松比仅为1.61~1.64和0.18~0.20,与测线北部7个台站的地壳结构相差甚大.分析表明地幔上涌对火山区莫霍面的局部抬升产生了一定影响,火山湖、黑空山、大-小空山和打鹰山下方应该存在一个相互联通的壳内岩浆囊.该岩浆囊在南北方向上的尺度约为20 km,热流活动以及幔源物质的侵入是地壳平均波速比和泊松比偏高的主要原因,它与热海附近的地温异常区分属两个不同的壳内岩浆存储系统.

     

云南思茅—中甸地震剖面的地壳结构

云南思茅—中甸宽角反射/折射地震剖面切割松潘—甘孜、扬子和华南三个构造单元的部分区域. 我们利用初至波和壳内反射波走时层析成像获得地壳纵波速度结构. 在获得新的地壳速度结构模型基础上,利用地震散射成像思想和低叠加次数的叠前深度偏移方法重建了研究区的地壳、上地幔反射结构. 综合分析研究区地壳P波速度模型和壳内地震反射剖面发现：沿测线从北至南地壳厚度从约50 km减薄至35 km左右,地壳厚度的减薄量主要体现在下地壳,剖面北段下地壳厚度约为30 km,剖面南段下地壳厚度仅为15 km左右;上地幔顶部局部位置P波速度值偏低,一般为7.6～7.8 km/s,反映出云南地区是典型的构造活动区的特点.剖面沿线地壳内地震反射发育,其中莫霍强反射出现在景云桥下方;在景云桥弧形断裂带8～10 km深处出现宽约50 km的强反射带.     

Insight into the Crustal Structure of the Eastern Marmara Region, NW Turkey

In order to investigate crustal structure beneath the eastern Marmara region, a seismic refraction survey was conducted across the North Anatolian Fault (NAF) zone in north west Turkey. Two reversed profiles across two strands of the NAF zone were recorded in the Armutlu Highland where a tectonically active region was formed by different continents. We used land explosions in boreholes and quarry blasts as seismic sources. A reliable crustal velocity and depth model is obtained from the inversion of first arrival travel times. The velocity-depth model will improve the positioning of the earthquake activities in this active portion of the NAF. A high velocity anomaly (5.6–5.8 km s⁻¹) in the central highland of Armutlu block and the low velocity (4.90 km s⁻¹) pattern north of Iznik Lake are the two dominant features. The crustal thickness is about 26 ± 2 km in the north and increases to about 32 ± 2 km beneath the central Armutlu block in the south. P-wave velocities are about 3.95 km s⁻¹ to 4.70 km s⁻¹ for the depth range between about 1 km and 5 km in the upper crust. The eastern Marmara region has different units of upper crust with velocities varying with depth to almost 8 km. The high upper crust velocities are associated with Armutlu metamorphic rocks, while the low velocity anomalies are due to unconsolidated sedimentary sequences. The western side of Armutlu block has complex tectonics and is well known for geothermal sources. If these sources are continuous throughout the portions of the crust, it may be associated with a granitic intrusion and deformation along the NAF zone. That is, the geothermal sources associated with the low velocity may be due to the occurrence of widespread shear heating, even shear melting. The presence of shear melting may indicate the presence of crustal fluid imposed by two blocks of the NAF system.     

汶川MS8.0地震：地壳上地幔S波速度结构的初步研究

2008年5月12日我国四川省汶川地区发生了震惊世界的M_S8.0地震.历史上,同类地震在大陆内部极为罕见.该地震深部构造背景的研究对理解其成因极为重要.本文利用中国地震局地质研究所地震动力学国家重点实验室在川西地区布设的大规模密集流动宽频带地震台阵记录的远震P波波形数据和接收函数非线性反演方法,得到了沿北纬31°线的19个台站下方120 km深度范围内的S波速度结构及台站下方地壳的平均泊松比.该观测剖面穿越了主震区,总长度约为420 km. 我们的结果揭示了川滇地块、松潘-甘孜地块和四川盆地三个不同地块构造差异.上述三个地块的地壳结构特征可以概括为：（1）四川盆地前陆壳幔界面向西侧倾斜并有较为明显的横向变形,地壳厚度存在46～52 km的横向变化,中下地壳S波速度存在横向变化,地壳平均泊松比值较高（0.28～0.31）,但在龙门山断裂带附近,显示了坚硬地壳的特征,地壳平均泊松比仅为0.2;（2）松潘-甘孜地块地壳厚度由西侧靠近鲜水河断裂的60 km,向东减薄为52 km,在14～50 km深度范围内存在S波速度2.75～3.15 km/s的楔状低速区,其厚度由西侧的～30 km向东逐渐减薄为～15 km,相应区域的地壳平均泊松比高达0.29～0.31; （3）鲜水河断裂西侧,川滇地块地壳结构相对简单,地壳厚度为58 km,并在26 km深度存在约10 km厚度的高速层,地壳内平均泊松比约为0.25;（4）汶川大震区在12～23 km深度上具有近乎4.0 km/s的S波高速结构,而其下方的地壳为低速结构,地壳平均泊松比0.31～0.32,汶川大震的余震序列主要分布在高速介质区域内. 本文的结果表明松潘-甘孜地块的地壳相对软弱;而且并不存在四川盆地向西侧的俯冲.我们认为在青藏高原东向挤压的长期作用下,四川盆地强硬地壳的阻挡作用可导致松潘-甘孜地块内部蓄积很大的应变能量以及上、下地壳在壳内低速层顶部边界的解耦,在龙门山断裂带附近形成上地壳的铲形逆冲推覆.汶川大地震及其邻近区域所具有的坚硬上地壳和四川盆地的阻挡作用为低应变率下的高强度应力积累创造了必要条件,而松潘-甘孜地块长期变形积累的高应变能构成了孕育汶川大地震的动力来源.     

Fine upper crustal structure of Jiashi strong earthquake swarm region in Xinjiang inferred from high resolution seismic refraction profile data

Introduction Jiashi-Artux area in southwest Xinjiang is one of the most active earthquake provinces at pre-sent in Chinese mainland. In the last century, about 3/4 strong earthquakes in Chinese mainland hit this area, and especially from January 21 to April 16 in 1997, 7 earthquakes with the magnituderanging from 6.0 to 6.9 occurred in a very small area of 9 km×18 km near Jiashi (ZHU et al, 1998). It has never taken place before in Chinese mainland that a series of strong earthquakes shoo…     

大姚地震前连续形变短期前兆特征

总结了2003年大姚6．2级地震云南模拟与数字化形变变化的短期异常特征。由模拟形变观测资料显示倾斜单分向加速变化显著,震前均出现转折变化;地应力阶变及大幅度波动为震前异常形态,而重力没有异常出现。异常展布距震中97～346km,异常台站比(距震中400km范围)为35．7％。数字化形变的9个测项中,楚雄倾斜和应变最为显著,尽管与降雨影响有关,但仍认为具有短期前兆意义。     

Structure of crust and upper mantle beneath the Ordos Block and the Yinshan Mountains revealed by receiver function analysis

A temporary seismological network of broadband three-component stations has been deployed N–S to investigate the crust and upper mantle structure across the Ordos Block and the Yinshan Mountains. P wave receiver functions reveal the Moho depth to be about 41 km beneath the central Ordos Block and down to 45 km beneath the northern Ordos Block, a slight uplifting to 42–43 km beneath the Hetao Graben, increasing to 47–48 km beneath the Yinshan Mountains and then decreasing to 44 km beneath the northern Yinshan Mountains along the profile. In the Ordos Block, the crustal Vp/Vs ratio (about 1.80) south to the Hetao Graben differs from that (about 1.75) beneath the center Ordos Block. The crustal Vp/Vs ratio is significantly lower (about 1.65–1.70) beneath the Yinshan Mountains. The P wave receiver function migration imaging suggests relatively flat discontinuities at 410 and 660 km, indicating the lack of a strong thermal anomaly beneath this profile at these depths, and a low S wave velocity anomaly in the upper mantle beneath the Hetao Graben. We suggest that the low S wave velocity anomaly may be attributable to heat and that the thermal softening advances the evolution of the Hetao Graben, while the lower-crustal ductile flows transfer from the Hetao Graben to the northern Ordos Block, resulting in crustal thickening.     

雷琼拗陷地壳结构特征

根据1984～1985年在雷琼地区完成、长约210 km的白延—江洪深地震测深资料,应用有限差分法、时间项法反演和二维非均匀动力学射线追踪等方法,得到白延—江洪地壳剖面的二维非均匀构造图像.结果显示：琼北及雷州半岛地壳厚约25～26 km,地壳平均速度在琼州海峡南段琼北凸起约6.16 km·s-1,海峡北段雷琼拗陷区约5.70 km·s-1;地壳速度偏低,下地壳尤为明显,仅约6.3～6.5 km·s-1;结晶基底显示为厚约2～4 km、速度5.5～6.0 km·s-1的速度梯度带,Moho界面上、下界面速度分别是6.5 km·s-1和8.05 km·s-1的一级速度间断面,没有发现明显的壳幔过渡带.白延—江洪剖面位于雷琼拗陷更新世玄武岩覆盖区,推测地壳厚度减薄以及低速度构造主要反映了该区域自新生代以来一直延续至第四系晚期的多期火山岩浆喷发热地幔对地壳底部的底浸、拆层与地壳仍处于相对温热状态相关.     

兴蒙造山带及邻区地壳速度结构特征

兴蒙造山带位于中亚造山带东段, 作为古亚洲构造域的重要组成部分, 由西伯利亚板块与华北板块碰撞拼合而成, 其经历了大陆裂解、洋盆扩张、洋壳俯冲消减和碰撞拼合造山等复杂的构造演化过程.为了利用壳幔结构约束造山带演化的深部过程, 跨越华北地块北缘、松辽—锡林浩特地块、兴安地块以及索伦—西拉木伦缝合带和二连—贺根山缝合带, 实施了一条520 km长的深地震测深剖面, 获得了高质量的人工源大当量的宽角反射和折射地震资料, 并采用地震动力学射线方法获得地壳速度结构.结果显示: (1)研究区地壳平均速度为6.15~6.3 km·s^-1, Pn波速度为7.8~8.2 km·s^-1; (2)地壳厚度约为36.1~42.2 km, 最厚位置(~42.2 km)对应地表大兴安岭主峰, 说明大兴安岭在此位置存在山根; (3)地壳速度在1.5~6.8 km·s^-1范围内, 认为在该区地壳内不存在洋壳物质; (4)主要断裂带或缝合带位于速度等值线变化剧烈的梯度带上; (5)速度结构显示研究区具有明显的横向分区和纵向分层的特点.地壳内速度剧烈变化特征表明兴蒙造山带的地壳物质组成不均匀, 尤其中下地壳, 速度等值线起伏剧烈.这种复杂的地壳速度结构应该与中生代以来多板块汇聚引发的多期区域性伸展和挤压作用有关.

     

Crustal structure in the Falcón Basin area, northwestern Venezuela, from seismic and gravimetric evidence

The Falcón Basin in northwestern Venezuela has a complex geological history driven by the interactions between the South American and Caribbean plates. Igneous intrusive bodies that outcrop along the axis of the basin have been associated with crustal thinning, and gravity modeling has shown evidence for a significantly thinned crust beneath the basin. In this study, crustal scale seismic refraction/wide-angle reflection data derived from onshore/offshore active seismic experiments are interpreted and forward-modeled to generate a P-wave velocity model for a 450 km long profile. The final model shows thinning of the crust beneath the Falcón Basin where depth to Moho decreases to 27 km from a value of 40 km about 100 km to the south. A deeper reflected phase on the offshore section is interpreted to be derived from the downgoing Caribbean slab. Velocity values were converted to density and the resulting gravimetric response was shown to be consistent with the regional gravity anomaly. The crustal thinning proposed here supports a rift origin for the Falcón Basin.     

A study of the crustal structure beneath the Himalayas from body waves

Summary The crustal structure beneath the Himalayas has been investigated using body wave data from near earthquakes having epicentres over the Himalayas and recorded by the observatories situated over, or very near, the foothills of the mountains. A three-layered crustal model, without the top sedimentary layer, with velocities for theP wave group in Granite I, Granite II and the Basaltic layer as 5.48, 6.00 and 6.45 and for theS wave group as 3.33, 3.56 and 3.90 km/sec respectively, has been interpreted. The upper mantle velocity for theP wave has been observed to be 8.07 km/sec and for theS wave as 4.57 km/sec. Average thickness for the Granite I layer has been computed as 22.7 km, for the Granite II layer as 16.3 km and for the Basaltic layer as 18.7 km. Crustal and sub-crustal velocities indicate a lower trend under the mountain. A thicker crust has been obtained beneath the Himalayas.     

The lithosphere in the central-eastern Mediterranean area

The lithosphere beneath the central-eastern Mediterranean area has been investigated by the inversion of the regional dispersion relations derived from analysis of surface waves. It is possible to distinguish several types of crust with average S-wave velocities in the range 3.0–3.8 km/sec, and thicknesses varying from a minimum of about 30 km, which corresponds to the Apennines, Crete and Otranto Channel regions, to a maximum of about 51 km beneath the Ionian Sea, which can be considered as a submerged continent. Associated with these crustal features, large lateral variations have been detected in the lithosphere thickness, which varies from a minimum of about 30 km corresponding to the Tyrrhenian Sea and south of Crete to a maximum of about 130 km corresponding to south-eastern Alps and north-central Greece, while the sub-Moho S-wave velocity varies in the range 4.2–4.8 km/sec. The constraint furnished by our results to the geological-tectonic setting of the investigated area, characterized by the continent continent collision between Africa and Europe, is pointed out.Publication No. 405, P. F. Geodinamica, CNR, Roma, Italy.