Body-waveforms and source parameters of some moderate-sized earthquakes near North Island, New Zealand

P-wave first motion and synthetic seismogram analysis of P- and SH-waveforms recorded at teleseismic distances on the WWSSN are used to estimate source parameters of seven of the largest earthquakes (6.1 ≤ m_b ≤ 6.3) that occurred in the vicinity of North Island, New Zealand since 1965. The source parameters of three other (m_b ≥ 6.1) events determined outside of this study are included and considered in the final analysis. Four of the earthquakes occurred at shallow depths (< 20 km), of which three were located within and to the north of North Island. Two of the shallow events show strike-slip and normal focal mechanisms with T-axes oriented in a manner consistent with their location in an area of known back-arc extension. One of the shallow events occurred in northern South Island and shows a reverse-type mechanism indicating horizontal contraction of the crust in an easterly azimuth. Six events occurred at intermediate depths (h = 39 to 195 km) of which five exhibit thrust mechanisms with T-axes consistently oriented near vertical. In the light of previously published plate tectonic models, the near vertical orientation of T-axes of the intermediate-depth events may be used to infer that the southern Kermadec plate boundary immediately north of North Island is not strongly coupled, and hence, not likely capable of producing great earthquakes. A similar inference cannot be made for the section of the Hikurangi Margin adjacent to North Island since the intermediate-depth events considered in this study lie to the north of this segment of the plate boundary.     

P-wave travel-time anomalies: Aleutian-Alaska region

The P-wave travel-time residuals in the Aleutian-Alaska region have been obtained by combining data from three underground nuclear tests so far carried out by the Atomic Energy Commission of the United States in Amchitka Island. The travel-time residuals show close correlation with the tectonism of the area. Attempts have been made to interpret qualitatively the pattern of the residuals in regard to the plate tectonics of the area. Areas of negative residuals appear interpretable as due to the underthrusted high-velocity lithospheric plate while areas with positive residuals seem to be associated with normal crust and upper mantle.     

Teleseismic P-wave residual investigations at Shillong,India

The northeast India region is seismically very active and it has experienced two large earthquakes of magnitude 8.7 during the last eight decades (1897 and 1950). We have analysed teleseismic P-wave residuals at Shillong, the only reliable seismic station operating in the region, to investigate a possible association of travel-time residual anomaly with earthquake occurrence. The period covered is from October 1964 through March 1976. The total number of events is 9479, including 1767 events with depth >/ 100 km. Six-monthly average residuals have been calculated. The standard deviations are less than 0.10 sec for these data sets. During the period of investigations, no major earthquake took place close to Shillong. The earthquake of June 1, 1969 with a magnitude (M_b) of 5.0, at an epicentral distance of 20 km from Shillong is the only significant event. This earthquake is found to be associated with a travel-time increase with a maximum amplitude of 0.4 sec. It appears that, in general, the P-wave velocity has decreased in the neighbourhood of Shillong since 1969. A quadrant-wise analysis of residuals indicates that the residual anomaly is most prominent in the SE quadrant from Shillong.     

Lithospheric structure of the Tornquist Zone resolved by nonlinear P and S teleseismic tomography along the TOR array

The main aim of the TOR project is to study the lithospheric–asthenospheric boundary structure under the Sorgenfrei–Tornquist Zone, across northern Germany, Denmark and southern Sweden. Relative arrival-time residuals of teleseismic P and S phases from 51 earthquakes, recorded by 150 seismic stations along the TOR array, were used to delineate the transition zone in the studied area. The effects of crustal structures were investigated by correcting the teleseismic residuals for travel-time variations in the crust based on a 3D crustal model derived from other data. The inversion was carried out for S phases. The results were then compared with the corresponding P-wave models. As expected, the derived models show that the relatively old and cold Baltic Shield has higher velocity at depth than the younger lithosphere farther South. The models show two sharp and distinct increases in depth to velocities which are low compared to our reference model, as we move from South to North. The location and sharpness of these boundaries suggests that the features resolved are, at least partially, compositional in origin, presumably related to mantle depletion. A sharp and steep subcrustal boundary is found roughly coincident with the southern edge of Sweden. This is below where the edge of the Baltic Shield is usually placed, based on surface geological evidence (the Sorgenfrei–Tornquist Zone). Another less significant transition is recognised more or less beneath the Elbe-lineament. Relatively high d(V_p / V_s) ratios under the central part of the profile (Denmark) indicate relatively low S-velocity in an area where a gravity high supports the hypothesis of extensive mafic intrusions.     

精确测量的地幔体波走时及地幔非均匀性的强度

用波形相关法精确地测定了在世界各地发生的87个6级以上地震的P波, PP波和Pdiff波的503个走时数据。记录这些地震波形的是新建于西太平洋地区的海洋半球地震观测网。我们利用这些高精度的走时数据研究了地幔体波的走时残差的范围及地幔非均匀性的强度。结果表明,P波、PP波和Pdiff波的走时残差最大分别为9 s ,11 s和15 s ,这为地幔层析成像反演中应该使用的体波走时残差数据的范围提供了重要信息。超出这一范围的走时残差数据不应该用于反演中,以免歪曲成像结果。我们发现,当震中距小于40°时,P波走时残差的范围为-6到+9 s。而对于40°到99°之间的震中距,P波走时残差的范围为-3到+5s。由于震中距越大,P波穿透地幔越深,我们这一结果提供了直接和确凿的证据,表明上地幔和地幔转换带中的横向非均匀性的强度要远胜于下地幔。我们精确测量的Pdiff波的走时数据表明,在地幔底部存在显著的低速异常,可能与地幔热柱或者超级地幔柱有关。我们使用了一个最新的三维全球层析成像模型来解释这些体波走时数据的空间变化。     

砂岩单轴压缩试验P波速度层析成像及声发射特性研究

深部矿井开采极易诱发矿柱型岩爆,矿柱型岩爆严重威胁着矿山的安全高效开采,为探究矿柱型岩爆破坏机制和前兆信息,选用自贡红砂岩进行单轴压缩试验,并采用主动超声和被动声发射对破裂过程进行监测,联合主动超声与被动声发射监测数据对波速进行层析成像反演,分析试样破裂过程中波速演化规律。研究结果表明：砂岩试样在加载过程中速度模型呈现高度非均质性,在加载过程中会出现低速区,且声发射事件集中出现在低速区内部;P波速度离散度可反映全局波速变化特征,在峰值附近变化剧烈,随荷载增加P波速度离散度持续增大;声发射事件在峰前、峰后定位结果相差较大,峰前阶段声发射事件随机分布,峰后阶段声发射事件定位结果集中。此外,研究发现选用均质速度模型进行声发射事件定位会增大定位误差,试样最终失稳前b值的降低表明试样大尺度裂纹活动加剧,导致岩石非均质性增大,也间接证明了采用非均质速度模型进行声发射震源定位的必要性。该研究结果可用于现场矿柱稳定性监测,定期反演获得矿柱波速变化为矿柱型岩爆提供先兆预警信息。     

Litho-asthenospheric structures of northern Italy as inferred from teleseismic P-wave tomography

Regional three-dimensional inversions of teleseismic P-wave travel time residuals recorded by high-frequency regional and local seismic networks operating along the Western Alps and surrounding regions were carried out and lithosphere and upper mantle P-wave velocity models down to 300 km were obtained.

Residuals of more than 500 teleseismic events, recorded by 98 fixed and temporary seismic stations, have been inverted.

The comparison between real residuals and the ones obtained from tomographic model indicates that the method is able to solve the feature of the regional heterogeneities.

Where the resolution is good, coherent lithospheric and upper mantle structures are imaged. In the shallower layers, high- and low-velocity anomalies follow the structural behaviour of the Alpine-Apenninic chains showing the existence of very strong velocity contrasts. In the deepest layers, velocity contrast decreases however two deep-seated high-velocity structures are observed. The most extended in depth and approximately trending NE-SW has been interpreted as a wreck of the oldest subduction responsible of the Alpine orogenesis. The second one, connected to the northwestern sector of the Apenninic chain, appears to vanish at depths greater than 180 km and is probably due to still active Apenninic roots.

Cross-sections depict the spatial trend of perturbations and in particular outline the sub-vertical character of the Alpine and Apenninic anomalies. Under the Ligurian Sea, the 3-D inversion confirms the uplift of the asthenosphere in agreement with the tectonic evolution of the basin.     

Crustal heterogeneity and seismotectonics of the region around Beijing, China

A detailed three-dimensional (3-D) P-wave velocity model of the crust and uppermost mantle under the Chinese capital (Beijing) region is determined with a spatial resolution of 25 km in the horizontal direction and 4–17 km in depth. We used 48,750 precise P-wave arrival times from 2973 events of local crustal earthquakes, controlled seismic explosions and quarry blasts. These events were recorded by a new digital seismic network consisting of 101 seismic stations equipped with high-sensitivity seismometers. The data are analyzed by using a 3-D seismic tomography method. Our tomographic model provides new insights into the geological structure and tectonics of the region, such as the lithological variations and large fault zones across the major geological terranes like the North China Basin, the Taihangshan and the Yanshan mountainous areas. The velocity images of the upper crust reflect well the surface geological and topographic features. In the North China Basin, the depression and uplift areas are imaged as slow and fast velocities, respectively. The Taihangshan and Yanshan mountainous regions are generally imaged as broad high-velocity zones, while the Quaternary intermountain basins show up as small low-velocity anomalies. Velocity changes are visible across some of the large fault zones. Large crustal earthquakes, such as the 1976 Tangshan earthquake (M=7.8) and the 1679 Sanhe earthquake (M=8.0), generally occurred in high-velocity areas in the upper to middle crust. In the lower crust to the uppermost mantle under the source zones of the large earthquakes, however, low-velocity and high-conductivity anomalies exist, which are considered to be associated with fluids. The fluids in the lower crust may cause the weakening of the seismogenic layer in the upper and middle crust and thus contribute to the initiation of the large crustal earthquakes.     

Critical evaluation and geological correlation of teleseismic phase data from Indian seismological stations

We evaluated the quality of seismic phase data from Indian seismological stations through the analysis of teleseismic travel times reported during 1976–83 and infer that only WWSSN stations (NDI, SHL, POO, KOD) apart from GBA and HYB can be rated satisfactory while the majority of stations (more than 40) produce very poor quality data sets. Detailed analysis of teleseismic P-wave travel time residuals shows that while the average structure of the upper mantle beneath India has high velocity (negative residuals) there are marked lateral variations. In particular, three zones of anomalous positive residuals (low velocity) are observed: one beneath the north western part of the Deccan trap, the second covering the southernmost peninsula (granulite terrain) and a third rather localized one, to the north of Delhi coinciding with Delhi-Haridwar ridge. New Delhi exhibits strong negative residuals in the E-SE quadrant along with negative station anomaly, implying that it is underlain by an anomalous high velocity crust/upper mantle. The negative residuals observed over India, continue beneath the Himalaya till the south of Lhasa but change sign further northward, suggesting the northern limit of the Indian upper mantle structure.     

Regional patterns in the paragenesis and age of inclusions in diamond, diamond composition, and the lithospheric seismic structure of Southern Africa

The Archean lithospheric mantle beneath the Kaapvaal–Zimbabwe craton of Southern Africa shows ±1% variations in seismic P-wave velocity at depths within the diamond stability field (150–250 km) that correlate regionally with differences in the composition of diamonds and their syngenetic inclusions. Seismically slower mantle trends from the mantle below Swaziland to that below southeastern Botswana, roughly following the surface outcrop pattern of the Bushveld-Molopo Farms Complex. Seismically slower mantle also is evident under the southwestern side of the Zimbabwe craton below crust metamorphosed around 2 Ga. Individual eclogitic sulfide inclusions in diamonds from the Kimberley area kimberlites, Koffiefontein, Orapa, and Jwaneng have Re–Os isotopic ages that range from circa 2.9 Ga to the Proterozoic and show little correspondence with these lithospheric variations. However, silicate inclusions in diamonds and their host diamond compositions for the above kimberlites, Finsch, Jagersfontein, Roberts Victor, Premier, Venetia, and Letlhakane do show some regional relationship to the seismic velocity of the lithosphere. Mantle lithosphere with slower P-wave velocity correlates with a greater proportion of eclogitic versus peridotitic silicate inclusions in diamond, a greater incidence of younger Sm–Nd ages of silicate inclusions, a greater proportion of diamonds with lighter C isotopic composition, and a lower percentage of low-N diamonds whereas the converse is true for diamonds from higher velocity mantle. The oldest formation ages of diamonds indicate that the mantle keels which became continental nuclei were created by middle Archean (3.2–3.3 Ga) mantle depletion events with high degrees of melting and early harzburgite formation. The predominance of sulfide inclusions that are eclogitic in the 2.9 Ga age population links late Archean (2.9 Ga) subduction-accretion events involving an oceanic lithosphere component to craton stabilization. These events resulted in a widely distributed younger Archean generation of eclogitic diamonds in the lithospheric mantle. Subsequent Proterozoic tectonic and magmatic events altered the composition of the continental lithosphere and added new lherzolitic and eclogitic diamonds to the already extensive Archean diamond suite.     

Structure of the uppermost mantle from long-range seismic observations in southern Germany and the Rhinegraben area

The crustal structure has been determined in the area between the Lorraine, the Bohemian massif and the northern Alps with considerable detail in recent years. But up to now little has been known about the velocity—depth structure of the uppermost mantle in this area. The situation changed recently when two recent seismic events near the northern and southern end of the Rhinegraben rift system were recorded to distances of 400 km. The explosions at the westernmost shotpoint of the international alpine refraction profile in 1975 were also observed in the Rhinegraben area up to the same distance. Earlier refraction seismic experiments between Steinbrunn near Basle and Boehmischbruck at the western border of the Bohemian massif also reach distances of 400 km. All these data lead to the rather high P-wave velocities of 8.5–8.6 km/s at depths between 40 and 50 km. These velocities are considerably higher than the average velocities of 8.2 km/s under other areas of western and central Europe, as for example the Bretagne in northwestern France and the North German Plain. There are indications of a minor velocity inversion in the uppermost mantle between Steinbrunn and Boehmischbruck. From the dispersion of surface waves there is good evidence that the regional high P-wave velocities are limited to a certain depth range only. This indicates rather pronounced lateral variations of the velocity—depth structure in the uppermost mantle of central Europe.     

Apparent velocity and azimuths of short period P-Waves recorded at indian WWSSN stations network

Over 200 earthquakes in the distance range 30°–90° and azimuthal range 0°–360°, recorded at Indian WWSSN stations, have been used in the present study. We have treated the four WWSSN recording stations i.e. New Delhi, Poona, Shillong and Kodai-Kanal, as Super Large Aperture Seismic Array (SLASA) Network with New Delhi being its cross-over point. Short period P-wave data as obtained from these stations have been analysed using a least square technique. Slowness and azimuthal anomalies have been computed for all these events. Relative time residuals have also been calculated. A velocity model has been derived on the basis of the slowness and travel-time data. The results do not indicate presence of any triplication in the travel-time curve. Variations in the relative residuals refer to the tectonic features beneath the recording stations. The P-wave velocity increases continuously in the lower mantle region and there is no indication for the presence of any appreciable velocity gradient.     

Land/ocean correlations during the last interglacial/glacial transition, Baffin Bay, northwestern North Atlantic: A review

Evidence from terrestrial sections, ice cores, and marine cores are reviewed and used to develop a scenario for environmental change in the area of the extreme northwest North Atlantic during marine isotope stages 5 and 4. The critical physical link between the landbased glacial chronology and marine events in Baffin Bay is the presence of carbonate rich drift along the Baffin Bay coast of Bylot Island and a detrital carbonate facies (Facies B) in Baffin Bay sediments. Cores from Baffin Bay/Labrador Sea can be dated by means of oxygen isotope variations and by peaks in the abundance of volcanic glass shards. One occurrence of Facies B is dated between late stage 5 and stage 4 and we correlate this event with the Eclipse Glaciation of Bylot Island and the Ayr Lake stade of the Foxe Glaciation of Baffin Island (= Kogalu aminozone). In contrast on West Greenland, amino acid racemization evidence suggests that the Greenland Ice Sheet developed throughout stage 4 and reached a maximum in stage 3 (Svartenhuk advance >40 ka). The oxygen isotope record in the Devon Island Ice Cap (northwest Baffin Bay) indicates that Baffin Bay was largely open during marine isotope stage 5. Analyses of shallow water molluscan and foraminiferal assemblages, deep-water foraminifera, pollen from Iand sections and deep-sea cores, and dinoflagellates from marine cores indicate that interglacial conditions prevailed during much of the stage glaciation.     

Teleseismic p-wave traveltime residuals and deep structure of the Aegean region

Teleseismic P-wave traveltime residuals have been measured at the Greek seismic stations with respect to the Herrin 68 tables. In spite of the large scatter, some insight into crustal and upper-mantle structure of the Aegean region can be gained. The average absolute residuals (observed minus Herrin traveltimes) are of the order of + 2 s. The most plausible interpretation is an efficient low-velocity zone in the upper mantle. Simple estimates of densities and subsequently gravity with the aid of Birch's law suggest that the Aegean region is underlain by hot expanded upper mantle, perhaps involving partial melting. The relative P residuals (observed minus Herrin traveltime differences between a station and Athens) are generally positive and can be interpreted with lateral variations of the LVZ or of the crust. The latter interpretation is supported in some cases by seismic refraction data. The azimuthal variation of the relative residuals at stations on the non-volcanic arc bears a distinct relation with the arc orientation. At Archangelos (Rhodes) where we “see” through the Benioff zone, the residuals from N to W are between -1 and -2 s and indicate a high-velocity slab sinking below the Aegean sea. At Vamos (Crete), Valsamata (Kephallenia), and Joanina (Pindus Range) the largest (smallest) residuals are along directions parallel (perpendicular) to the arc. This can be interpreted by crustal thickening under the sedimentary arc and/or by velocity anisotropy with the maximum perpendicular to the arc. On the whole, our study supports the hypothesis that the Aegean region is a trench—island-arc—marginal-sea system.     

Neocene and Quaternary extent and geometry of the subducted Pacific Plate beneath North Island, New Zealand: Implications for Kaikoura tectonics

The extent and geometry of the obliquely subduced oceanic Pacific Plate beneath North Island, New Zealand, for five million year intervals through the mid-Miocene to Quaternary, are presented in a series of maps and cross-sections. These show that the subducted plate progressively increased its extent from NE to SW beneath the North Island, and in the more northern regions where it was first emplaced, concomitantly increased its dip from 10° to 50°.The changing extent and geometry of the subducted slab has been established from the age pattern of orogenic andesites and from the geochemical K₂O-h parameter of depth of magma generation. The radiometric dates show a migration of the volcanic front back towards the trench at an average rate of 20 km/My. The trenchward migration is explained by a model of increasing slab dip which is corroborated by the K₂O data calibrated against the presently active arc (Taupo Volcanic Zone). With the exception of northern Coromandel Peninsula, the andesitic magmas were generated at 85–100 km depth. The interpretation of the dates adopted here indicates that the subducted slab originated at the NE-SW trending Kermadec-Hikurangi Trench, and implies a different and much simpler evolution of the Australia-Pacific plate boundary in the vicinity of North Island than other recent models.Subduction geometry has been found elsewhere to be a principal influence upon the state of stress and deformational style in an over-riding plate. The possibility is explored that the timing, nature and pattern of the Neogene to Quaternary Kaikoura Orogeny in North Island is due to this influence. Apart from the effect of oblique subduction in eastern North Island, there is an accord between the onset of deformation and the emplacement sequence of the shallow slab beneath North Island, and between the change in subduction geometry and a progressive north to south change in northern North Island from compression to extension.     

Seismic crustal structure of the Limpopo mobile belt,Zimbabwe

A 145 km N–S seismic traverse was deployed to determine the crustal structure of the Limpopo mobile belt in southern Zimbabwe and the nature of its northern boundary with the Zimbabwean craton. Rockbursts from South African gold mines to the south and regional seismicity from the Kariba-South Zambia belt to the north were used as seismic sources. P-wave relative teleseismic residuals were also measured to assess whether any velocity contrast between the craton and the mobile belt extended into the upper mantle.Interpretation of reduced travel times from the local Buchwa iron-ore mine blasts, which were broadside to the traverse, revealed an upper crustal interface in the Limpopo mobile belt at a depth of 5.8 ± 0.6 km, dividing material with a velocity of about 5.8 km/s from that of about 6.4 km/s. On the craton, arrivals from the same source showed a 4.4 ± 0.5 km thick 5.5 km/s layer overlying crust of about velocity 6.5 km/s. P-wave arrivals from the regional seismicity were used to construct a crustal cross-section. Absolute crustal thickness was tentatively estimated from the identification of a Moho reflection on the mine blast recordings. To the south of Rutenga, the crust thins from around 34 km to 29 km in association with a positive gravity anomaly centred over the late-Karoo Nuanetsi Igneous Province and Karoo Tuli Syncline. North of Rutenga to the boundary with the Zimbabwean craton, the crust is about 34 km thick. The craton boundary was found to be a steeply southerly dipping zone associated with high-velocity material, which could either be deep-seated greenstones or mafic material associated with the margin in the region studied. This zone divides cratonic crust, which was found to be about 40 km thick, from that typical of the mobile belt and implies a step in the Moho of around 6 km.Analysis of relative teleseismic residuals showed that the velocity contrasts are not confined to the crust but extend into the uppermost upper mantle with the cratonic lithosphere being about 4% faster than that of the Limpopo mobile belt. The resolution of the technique is such that it is difficult to ascertain whether these differences are features of Precambrian evolution or are due to reactivation of the upper mantle during Karoo igneous and tectonic activity.     

Crustal structure of the Eastern Alps along the TRANSALP profile from wide-angle seismic tomography

The objective of the TRANSALP project is an investigation of the Eastern Alps with regard to their deep structure and dynamic evolution. The core of the project is a 340-km-long seismic profile at 12°E between Munich and Venice. This paper deals with the P-wave velocity distribution as derived from active source travel time tomography. Our database consists of Vibroseis and explosion seismic travel times recorded at up to 100 seismological stations distributed in a 30-km-wide corridor along the profile. In order to derive a velocity and reflector model, we simultaneously inverted refractions and reflections using a derivative of a damped least squares approach for local earthquake tomography. 8000 travel time picks from dense Vibroseis recordings provide the basis for high resolution in the upper crust. Explosion seismic wide-angle reflection travel times constrain both deeper crustal velocities and structure of the crust–mantle boundary with low resolution. In the resulting model, the Adriatic crust shows significantly higher P-wave velocities than the European crust. The European Moho is dipping south at an angle of 7°. The Adriatic Moho dips north with a gentle inclination at shallower depths. This geometry suggests S-directed subduction. Azimuthal variations of the first-break velocities as well as observations of shear wave splitting reveal strong anisotropy in the Tauern Window. We explain this finding by foliations and laminations generated by lateral extrusion. Based on the P-wave model we also localized almost 100 local earthquakes recorded during the 2-month acquisition campaign in 1999. Seismicity patterns in the North seem related to the Inn valley shear zone, and to thrusting of Austroalpine units over European basement. The alignment of deep seismicity in the Trento-Vicenza region with the top of the Adriatic lower crust corroborates the suggestion of a deep thrust fault in the Southern Alps.     

Frequency-magnitude,depth, and time relations for earthquakes in an island arc: North Island,New Zealand

Analysis of over 1400 earthquakes in the North Island of New Zealand from 1955 to 1969, comprising all shocks with m_l ? 4.3 for shallow, and M_L ? 4.0 for deep events, reveals several empirical relationships between the depth and the equivalent radius of the area occupied by shocks, the number and density of the shocks, and the coefficient b and the maximum magnitude. The coefficient b increases linearly with depth from 1.0 for shallow earthquakes to 1.4 for those at a depth of 120 km, and then decreases to 0.75 at 300—350 km. The variation with depth shows clear inverse correlation with the distribution of maximum stress along the downgoing slab, calculated for several slab models by Smith and Toksöz. Similarly, the maximum magnitude at different depths correlates distinctly with the distribution of the principal stress. Time variations of the coefficient b and the rate of earthquake occurrence, for both shallow and deep earthquakes, have an oscillatory character, with a period of 7–8 years. These variations also imply that shallow and deep seismicity are mutually dependent.     

Residual gravity variations in volcanic areas: Constraints to the interpretation of uplift episodes at Campi Flegrei,Italy

An elastically homogeneous, compressible half space model with vertical density stratification is employed to compute the displacement field and the gravity variations produced by the inflation of a spherically symmetric deformation source. Contributions to gravity variations are produced by (1) the displacements of the free surface and of subsurface layers, (2) dilation/contraction of the medium, (3) the displacement of source boundaries and, possibly, by (4) new mass input from remote distances into the source volume. Two extreme cases were examined in detail, in which the magma chamber is identified as the deformation source: in the first model the volume and pressure increase is due to input of magma with the same density of the magma already resident in the magma chamber (deformation source with constant density), in the second model it is due to magma differentiation (deformation source with constant mass). In recent years (1970–72 and 1982–84) two inflation episodes took place in the Campi Flegrei caldera (Italy), characterized by significant ground uplift and strongly correlated gravity variations. From the comparison between measured and computed gravity residuals (free-air-corrected gravity variations) we can assess that an inflation source with constant density would predict gravity residuals consistent with those measured during the phase of uplift (within experimental errors), while an expansion at constant mass would predict gravity residuals much lower than observed. However, during the subsidence phase, which followed after the maximum uplift in 1984, gravity residuals at most sites were completely different from those measured during the uplift phase, suggesting that more mass left the system during deflation than entered during inflation. An alternative model is then proposed, in which the deformation source is ascribed to fluid pressure increase within a geothermal system, close to the critical point, shallower than the magma chamber.     

Structure of the upper mantle in the Northeastern Atlantic close to the azores-gibraltar ridge from surface-wave and body-wave observations

Phase velocities of teleseismic Rayleigh waves have been measured in the central North Atlantic on both sides of the Azores-Gibraltar Ridge (AGR) by means of a specially designed long-period station network. The dispersion data obtained were regionalized and then subjected to a “hedgehog” inversion, which gives a set of upper mantle models compatible with the observational data within specified error bounds.Reasonable model solutions were selected by using regional body-wave observations, such as P_n- and S_n-wave velocities determined from earthquakes along the AGR. The S_(itn) velocities measured indicate that the shear-wave velocity in the mantle part of the lithosphere is much higher on the northern side of the AGR. Strongly negative P-wave residuals in this area indicate faster seismic propagation than implied by the Jeffreys-Bullen travel-time tables, while propagation is much slower in the Gulf of Cadiz area. Furthermore the residuals show a clear difference for paths through oceanic and continental domains and suggest that the transition between these two domains extends much further into the ocean on the southern side of the AGR than on the northern side.The proposed model for the structure of the upper mantle in that region shows that there exists a pronounced velocity contrast across the AGR. Thickening of the lithospheric plate with increasing plate age is indicated to the south of the ridge. The greatest thickness is reached close to the continental margin within a zone about 500 km wide, whose velocity close to the Canary Islands and Madeira is significantly lower, probably due to the well-known volcanic activity there. These observations together with the travel time residuals reveal that this zone seems to be of a transitional nature somewhere between a continental and oceanic structure.