Site response using coda-wave techniques: Applications to short period data from central and southern Italy

The method of coda waves was applied to two different sets of data for the evaluation of the relative site response.The first set of data consists of low magnitude earthquakes with closely spaced locations, recorded at a small aperture array of velocimeters located in the Abruzzo region, central Italy. The second set of data is composed of events with epicentral distance ranging from 20 to 300 km, recorded at a seismological network with an aperture of about 100 km located in the Puglia region, southern Italy.Results show that the coda wave method furnishes stable estimates of the site effect. An amplification, relative to an arbitrary site, of a factor of about 2 occurs in the 1.7–6 Hz frequency band for two stations of the Abruzzo network, while an amplification factor of about 0.5 occurs in the whole frequency band (1–24 Hz) for one of the stations of the Puglia network. This station is located in an area which is correlated with a low macroseismic intensity anomaly.     

Earthquake activity during the 1983 Etna eruption

Seismic stations, with automatic P-picking and satellite retransmission were set up on Mount Etna following the eruption started on March 1983. Positions of the stations were chosen in order to complement the permanent telemetered network of Catania University.Comparison between locations obtained by both networks were made for earthquakes recorded by at least 5 ARGOS DCP (Data Collection Platform) stations. We observed a satisfactory agreement for events inside both networks.By merging data of both networks, it has been possible to locate more than 50 earthquakes for which separate computation was not possible due to the low number of arrivals.On 3rd-4th June a swarm of deep seismic events was observed. Hypocenters of these earthquakes are clearly located in a NNW-SSE-trending vertical zone of 5 km width at a depth of 7–36 km.Changes in the distribution of shallow seismic activity, before and after this swarm, have been observed.     

Combined discrete and continuous gravity observations at Mount Etna

Systematic investigation of discrete gravity measurements has continued at Mount Etna since 1986. The network now covers an area of 400 km² with about 70 stations 0.5–3 km apart. Mass redistributions occurring at depths ranging between about 8 km below sea level and a few hundred metres below the surface (magma level changes within the shallower parts of the feeding conduits) have been identified from these data. Conventional (discrete) microgravity monitoring on a network of stations furnishes only instantaneous states of the mass distribution at continuously active systems. In order to obtain information on the rate at which the volcanic processes (and thus mass transfers) occur, three stations for continuously recording gravity where installed on Mount Etna in 1998. A 16-month long sequence from one of the continuously running stations (PDN, located 2 km from the active northeast crater at the summit of Etna volcano) is presented. After removing the effects of Earth Tide and tilt, the correlation of the residual gravity sequence with simultaneous recordings of meteorological parameters acquired at the same station was analysed. Once the meteorological effects have also been removed, continuous gravity changes are within 10 μGal of gravity changes measured using conventional microgravity observations at sites very close to the continuous station. This example shows how discrete and continuous gravity observations can be used together at active volcanoes to get a fuller and more accurate picture of the spatial and temporal characteristics of volcanic processes.     

Magnetic measurements on La Soufriere Volcano, Guadeloupe (Lesser Antilles), 1976–1984: A re-examination of the volcanomagnetic effects observed during the volcanic crisis of 1976–1977

During the seismovolcanic crisis of 1976–1977 at La Soufrière on Guadeloupe, a magnetic network of 12 reference markers was set up. Measurements of the intensity of the earth's magnetic field, carried out up to once a day at each marker, showed volcano-magnetic variations of several nanoteslas (nT). The variations, at certain markers, were more or less concealed by transient magnetic variations due to anomalies in conductivity. As early as 1978, measurements were resumed and a telemetering network was coupled with the network of reference markers, to which 4 new markers were added. A detailed study of conductivity anomalies was carried out on the entire volcano. Contrasts in conductivity linked to the existence of a superficial conducting surface, on a SSW/NNE axis, located south of the volcano, caused a great lack of homogeneity in the field variations measured at the surface. Variations greater than about 10 nT appeared in the difference in intensity of the earth's magnetic field between two stations.No long-term magnetic variation was observed between 1978 and 1984. On the network of markers, the accuracy of measurements of volcanic effects was at best 2 nT. Measurements carried out on the telemetering network during the night refined these results, since their accuracy was 1 nT. The only significant volcanic crisis between 1978 and 1984 (5–7 January 1981) seems to be observed by telemetering stations. All the measurements carried out in periods of volcanic inactivity make it possible to re-examine the crisis of 1976–1977. Though volcanomagnetic effects over short periods cannot be accurately determined, variations with a time constant of several weeks were present over the entire volcano. These variations were as high as 7–8 nT in remote stations and they can be linked to the three major phases of eruptive activity at La Soufrière during the crisis of 1976–1977.     

Spatio-temporal distribution of seismic activity during the Umbria-Marche crisis, 1997

We present the spatio-temporal distribution of more than 2000 earthquakesthat occurred during the Umbria-Marche seismic crisis, between September 26and November 3, 1997. This distribution was obtained from recordings of atemporary network that was installed after the occurrence of the first two largest shocks (Mw =, 5.7, Mw = 6.0) of September 26. This network wascomposed of 27 digital 3-components stations densely distributed in theepicentral area. The aftershock distribution covers a region of about 40 km long and about2 km wide along the NW-SE central Apennines chain. The activity is shallow,mostly located at less than 9 km depth. We distinguished three main zonesof different seismic activity from NW to SE. The central zone, that containsthe hypocenter of four earthquakes of magnitude larger than 5, was the moreactive and the more complex one. Sections at depth identify 40–50^°dipping structures that agree well with the moment tensor focalmechanisms results. The clustering and the migration of seismicity from NW to SE and the generalfeatures are imaged by aftershock distribution both horizontally and at depth.     

Geological evidence for strong historical earthquakes in an “aseismic” region: The Pollino case (Southern Italy)

The Pollino Range is the southernmost segment of the Southern Apennines at the boundary with the Calabrian Arc. While several strong earthquakes (magnitude 6.5–7.0) have occurred in nearby regions, the Pollino area has no known historical record of seismic events of magnitude > 5. We carried out an aerial photograph interpretation and a field survey of the Pollino fault (the major Quaternary normal fault of the area) in order to characterize geologically the seismic potential of this structure. We dug two sets of trenches across fault scarps within the apecies of latest Pleistocene to Holocene alluvial fans at the Masseria Quercia Marina (MQM) and Grotta Carbone (GC) sites, in the central segment of the southern Pollino Range front. At both sites we identified two surface faulting events affecting the alluvial fan deposits and two overlying colluvial units of historical age. The penultimate event produced a vertical offset of 80–90 cm at GC and 50–60 cm at MQM; while the last event produced a vertical offset of 40–50 cm at GC and few centimeters of offset at MQM. Detailed geomorphological field observations suggest that the two historical earthquakes reactivated the entire length of the Masseria Marzano-Civita segment of the Pollino fault (rupture length about 18 km). For events in this range of rupture length and vertical displacement, comparison with surface faulting earthquakes in the Apennines (and abroad) indicates a magnitude of 6.5–7.0. Therefore, the maximum potential earthquake and the seismic hazard of the Pollino area are significantly larger than that suggested by the available historical seismic catalogue.     

A Multidisciplinary Study of the DPRK Nuclear Tests

The Democratic People Republic of Korea announced two underground nuclear tests carried out in their territory respectively on October 9th, 2006 and May 25th, 2009. The scarce information on the precise location and the size of those explosions has stimulated various kinds of studies, mostly based on seismological observations, by several national agencies concerned with the Nuclear Test Ban Treaty verification. We analysed the available seismological data collected through a global high-quality network for the two tests. After picking up the arrival times at the various stations, a standard location program has been applied to the observed data. If we use all the available data for each single event, due to the different magnitude and different number of available stations, the locations appear quite different. On the contrary, if we use only the common stations, they happen to be only few km apart from each other and within their respective error ellipses. A more accurate relative location has been carried out by the application of algorithms such as double difference joint hypocenter determination (DDJHD) and waveform alignment. The epicentral distance between the two events obtained by these methods is 2 km, with the 2006 event shifted to the ESE with respect to that of 2009. We then used a dataset of VHR TerraSAR-X satellite images to detect possible surface effects of the underground tests. This is the first ever case where these highly performing SAR data have been used to such aim. We applied InSAR processing technique to fully exploit the capabilities of SAR data to measure very short displacements over large areas. Two interferograms have been computed, one co-event and one post-event, to remove possible residual topographic signals. A clear displacement pattern has been highlighted over a mountainous area within the investigated region, measuring a maximum displacement of about 45 mm overall the relief. Hypothesizing that the 2009 nuclear test had been carried out close to the area where the displacement has been observed through the DInSAR technique, its relation with the epicenter location obtained through seismological processing has been discussed as a possible alternative hypothesis with respect to the preferred solutions reported by the nuclear explosion database (NEDB). The distance of about 10 km between the two places can be considered acceptable in light of the possible systematic location shifts commonly observed in the seismological practice over a global scale. The difference between the m _b magnitudes of the two tests could reflect differences in geological conditions of the two test sites, even if the yield of the two explosions had been the same.     

On reproducibility of relative locations of recurrent earthquakes recovered from similarity of their waveforms at a single station

Relative locations of recurrent earthquakes can be recovered from values of cross-correlation coefficients between their waveforms at a single station, provided that an idea of Menke (1999) that the values depend inversely exponentially from separation between the sources is true. An algorithm developed in Gnyp (2013) based on this idea was applied here for recovering relative locations of the two groups of the 2005–2006 recurrent earthquakes in Mukacheve, Ukraine, from the sets of cross-correlation coefficients between their records at the four local network stations to find out whether the locations were in principle reproducible independently at each of the stations, at epicentral distances from about 6 to 64 km. The resulting sets of locations occurred in good agreement with each other, which may be considered as an experimental evidence for validity of the Menke’s dependence over the corresponding range of distances. The reproducibility of relative locations has been analyzed in terms of existence and uniqueness of the inverse problem’s solution, and in connection with data completeness and the influence of noise on it.     

Real-time seismology for the 05/12/2008 Wenchuan earthquake of China: A retrospective view

The devastating 05/12/2008 Wenchuan earthquake (M_w7.9) in Sichuan Province of China showed very few precursory phenomena and occurred on a fault system once assigned to be of moderate long term seismic risk. Given the existing coverage of seismograph stations in Sichuan Province, real-time seismology could have been effective in avoiding some earthquake damage and helping post-earthquake emergency response. In a retrospective view, we demonstrated that the epicenter can be located with 20 km accuracy using just two broadband stations with three-component, which takes only about 10 s after the onset of the earthquake. Initial magnitude is estimated to be M7 with the Tc measurement over first 4 seconds of P waves. Better magnitude estimate can be obtained within 2 min by modeling Pnl waves for stations about 500 km away where the S waveforms are clipped. The rupture area is well revealed by teleseismically-recorded >M5 early aftershocks within two hours after the mainshock. Within a few hours, teleseismic body waves were inverted to derive a more detailed rupture process and the finite fault model can be readily used to calculate ground motions, thus providing vital information for rescue efforts in the case where no real-time strong motion records are available. Supported by the Knowledge Innovation Program of the Chinese Academy of Sciences (Grant No. KZCX2-YW-116-1) and National Key Technology Research and Development Program of China (Grant No. 2006BAC03B00)     

The 1997 Umbria-Marche (Italy) earthquake sequence: analysis of the data recorded by the local and temporary networks

We present some preliminary results obtained from thejoined analysis of the data collected by the permanentand the temporary networks operating in the area ofthe earthquake sequence that followed (andanticipated) the 26 September, Central Italy, mainshocks. In particular, these earthquake data haveallowed us to determine a well constrainedwave-velocity model (both P and S) with stationcorrections which demonstrated to produce robusthypocentral locations. These velocity modelswith station corrections have been used forre-locating the whole September 1997–July 1998subset of data of the permanent network, and theprevious background seismicity, starting from May1996. The focal mechanisms of the largest events werealso obtained from an analysis of the first-motionpolarities.Our results indicate that 1) the seismic activityaligns on a SE-NW trend for a total length of about50 km of extension; 2) the focal depth of theseevents is restricted to the range 0–9 km; 3) mostevents can be related to sub-parallel SW dipping faultplanes; 4) focal mechanisms of the largest shocks(M_L > 4) show a coherent behaviour, withnormal fault solution on SSE-NNW striking, SW dippingplanes; 5) the space-time evolution of the activitydisplays a discontinuous mode of energy release, withdifferent episodes of activation and an apparentclustering of aftershocks at the edges of the areaswhich presumably ruptured in the main shocks.     

Monitoring the freely dissolved concentrations of polycyclic aromatic hydrocarbons (PAH) and alkylphenols (AP) around a Norwegian oil platform by holistic passive sampling

In order to assess the environmental impact of aquatic discharges from the offshore oil industry, polar organic chemical integrative samplers (POCIS) and semipermeable membrane devices (SPMDs) were deployed around an oil platform and at reference locations in the North Sea. Exposure to polycyclic aromatic hydrocarbons (PAH) and alkylated phenols (AP) was determined from passive sampler accumulations using an empirical uptake model, the dissipation of performance reference compounds and adjusted laboratory derived sampling rates. Exposure was relatively similar within 1–2 km of the discharge point, with levels dominated by short chained C1–C3 AP isomers (19–51 ng L⁻¹) and alkylated naphthalenes, phenanthrenes and dibenzothiophenes (NPD, 29–45 ng L⁻¹). Exposure stations showed significant differences to reference sites for NPD, but not always for more hydrophobic PAH. These concentrations are several orders of magnitude lower than those reported to give both acute and sub-lethal effects, although their long term consequences are unknown.     

Locating pyroclastic flows on Soufriere Hills Volcano, Montserrat, West Indies, using amplitude signals from high dynamic range instruments

Pyroclastic flows are located using amplitude signals from a seven-station high dynamic range seismograph array located 1.9–6.1 km from Soufriere Hills Volcano in Montserrat, West Indies. Locations are determined by measuring the seismograph signal amplitude for an event recorded at several stations in a moving time window analysis. For a given window, the measured amplitudes are corrected to a trial source location by removing the effect of the surface wave geometric spreading, instrument gain, and the attenuation at calculated travel-times. The trial source location is then compared to other trial locations via an iterative localised grid search where the root-mean-squared amplitude residual (δA) is minimised. The process is repeated for subsequent time steps resulting in a best-fit event location and size through time. The method has been tested on four small events occurring on April 8, 1999, August 12, 1999, February 25, 2001, and July 4, 2001, when visual observations of pyroclastic flows coincided with good seismograph station coverage (number of stations 5, azimuthal gap <160°). Based on the location results the four events propagated 0.5, 1.4, 1.3 and 1.0 km from the dome, and had maximum attenuation-corrected reduced displacements (D_RQ) of 9.0, 2.8, 6.9 and 2.3 cm² and maximum pyroclastic flow velocities of 7, 30, 20 and 8 ms⁻¹, respectively. A time-lapse video of the event of August 12, 1999, shows that amplitude-based location through time closely matches the observed run-out distance and velocity. In contrast, amplitude-based locations for the events of April 8, 1999, and July 4, 2001, underestimated the actual flow run-out by 1.5 km. Underestimation of the true run-out distance is probably due to both the increased distribution of sources as coherent dome material disaggregates into many blocks, and signal contamination from other sources. Results indicate that pyroclastic flows and rockfalls can be located using amplitude signals from high dynamic range seismograph stations yielding estimates of size and trajectory, regardless of visibility conditions on the volcano. This new method is being tested as a hazard mitigation and research tool on Montserrat.     

2013年4月20日四川芦山M7.0级地震与余震精确定位及发震构造初探

使用汇集在四川台网中心的固定台站、震后架设的流动台站、周边水库台站等震中距150 km以内的震相数据,选用分层速度模型,对芦山7.0级地震及震后9天内的余震利用双差定位法进行了重新定位.给出了芦山7.0级地震的发震时刻为2013-04-20 08:02:46.8,震中位置30.278°N,102.989°E,震源深度16.67 km,给出了3324次余震的双差定位结果,并对发震构造进行了分析.结果表明:芦山地震主破裂长度约40 km,下倾宽度约20 km,破裂视面积约800 km²,主破裂沿南西走向,倾角约40°.余震震源优势深度为10~22 km.余震沿南西走向,主要集中于大邑-名山断裂上盘.     

A low magnitude seismic sequence near Isernia (Molise,Central Italy) in January 1986

Following the increase in seismic activity which occurred near Isernia (Molise, Central Italy) in January 1986, a digital seismic network of four stations with three-component, short-period seismometers, was installed in the area by the Osservatorio Vesuviano. The temporary network had an average station spacing of about 8–10 km and, in combination with permanent local seismic stations, allowed the accurate determination of earthquake locations during an operating period of about one month. Among the 1517 detected earthquakes, 170 events were selected with standard errors on epicentre and depth not greater than respectively 0.5 and 1.5 km. The most frequent focal depths ranged between 4 and 8 km, while the epicentres distribution covered a small area NE of Isernia of about 10 km². A main rupture zone could not be clearly identified from the spatial distribution of the earthquakes, suggesting a rupture mechanism involving heterogeneous materials. The activity was characterized by low energy levels, the largest earthquake, on January 18, 1986, havingM _D =4.0. The time sequence of events and pattern of seismic energy release revealed a strong temporal clustering of events, similar to the behaviour commonly associated with seismic swarms.     

Exploration of geothermal structure in Puga geothermal field, Ladakh Himalayas, India by magnetotelluric studies

To understand the crustal electric structure of the Puga geothermal field located in the Ladakh Himalayas, wide band (1000 Hz–0.001 Hz) magnetotelluric (MT) study have been carried out in the Puga area. Thirty-five MT sites were occupied with site spacing varying from 0.4 to 1 km. The measurements were carried out along three profiles oriented in east–west direction. After the preliminary analysis, the MT data were subjected to decomposition techniques. The one-dimensional inversion of the effective impedance data and the two-dimensional inversion of the TE (transverse electric) and TM (transverse magnetic) data confirm the presence of low resistive (5–25 Ω m) near surface region of 200–300 m thick in the anomalous geothermal part of the area related to the shallow geothermal reservoir. Additionally, the present study delineated an anomalous conductive zone (resistivity less than 10 Ω m) at a depth of about 2 km which is possibly related to the geothermal source in the area. A highly resistive basement layer separates the surface low resistive region and anomalous conductive part. The estimated minimum temperature at the top of conductive part is about 250 °C. The significance of the deeper conductive zone and its relation to the geothermal anomaly in the area is discussed.     

Surface heat flow and igneous intrusion in the Cambay Basin, India

Heat flow values from some additional locations in the Cenozoic Cambay Basin have been determined. Together with the previously published data, they show that the heat flow is moderate (55–67 mW/m′) in the southern part of the basin towards Broach and Ankleswar, and that there is a clear trend of high heat flow (75–93 mW/m²; range of average values for six different, widely separated, locations) in a part of the basin located north of the Mahisagan river between Cambay and Mehsana along a stretch of about 140 km. Conductive steady state geotherms, calculated using observed high surface heat flow values and appropriate models show, beneath the Cambay-Mehsana area, a large degree of melting in the lower crust and upper mantle, which is not suggested by the existing geodata. Considering this aspect and taking into account the existence of a normal crust about 37 km thick below the Cambay-Tarapur and Ahmedabad-Mehsana blocks (as obtained from deep seismic soundings), it has been inferred that the heat flow anomaly is due to transient thermal perturbations introduced from tectonic activity in the form of magmatic intrusions. A careful analysis of heat flow, gravity and other related geodata point out and support the possibility of a Miocene/Pliocene basic intrusive body at a depth of around 10 km under the Cambay-Mehsana area. Further, the consistent trend of the thermal and gravity fields indicates thinning of the postulated intrusive body from Cambay towards Mehsana.     

Evidence for fault-related directionality and localized site effects from strong motion recordings of the 2003 Boumerdes (Algeria) earthquake: Consequences on damage distribution and the Algerian seismic code

The Algiers–Boumerdes region has been struck by a destructive magnitude 6.8 (M_w) earthquake on May 21, 2003. The study presented in this paper is based on main shock strong motions from 13 stations of the Algerian accelerograph network. A maximum 0.58g peak ground acceleration (PGA) has been recorded at 20 km from the epicenter, only about 150 m away from a PGA of 0.34g, with both a central frequency around 5 Hz, explained by a strong very localized site effect, confirmed by receiver function technique results showing peaks at 5 Hz with amplitudes changing by a factor of 2. Soil amplifications are also evidenced at stations located in the quaternary Mitidja basin, explaining the higher PGA values recorded at these stations than at stations located on firm soil at similar distances from the epicenter. A fault-related directionality effect observed on the strong motion records and confirmed by the study of the seismic movement anisotropy, in agreement with the N65 fault plan direction, explains the SW–NE orientation of the main damage zone. In the near field, strong motions present a high-frequency content starting at 3 Hz with a central frequency around 8 Hz, while in the far field their central frequency is around 3 Hz, explaining the high level of damage in the 3- to 4-story buildings in the epicentral zone. The design spectra overestimate the recorded mean response spectra, and its high corner frequency is less than the recorded one, leading to a re-examination of the seismic design code that should definitively integrate site-related coefficient, to account for the up to now neglected site amplification, as well as a re-modeling of the actual design spectra. Finally, both the proposed Algerian attenuation law and the worldwide laws usually used in Algeria underestimate the recorded accelerations of the 6.8 (M_w) Boumerdes earthquake, clearly showing that it is not possible to extrapolate the proposed Algerian law to major earthquakes.     

The establishment of the national gravity standard base net of Egypt (N.G.S.B.N.-77)

The earliest gravity measurements in Egypt were carried out by the Anglo-Egyptian Oil Company using Holweck-Lejay pendulum in the period 1937–1940.The measurement comprised few stations (about 41) in restricted areas along both the Western and Eastern coasts of the Gulf of Suez.During the period 1950–1951 Prof. G. Woolard made 21 gravimetric measurements in Egypt as a part of the world-wide gravity base net. Thi was refeered to Potsdam absolute value which was considered as the World Gravity Datum.In 1974, a new world-wide gravity net was published (I.G.S.N.71) which introduced a new concept for the world gravity reference datum, and consequently has replaced the Potsdam System. Since then, it was recommended that a national gravimetric network of base stations in Egypt on the basis of the I.G.S.N.-71 be established within the frame of establishing a gravity map of Egypt.In the period from January 1975, The General Petroleum Company carried out 624 gravity observations throughout the country. These measurements have been combined with the I.G.S.N.-71 values to form 65 gravity base points which in turn from the N.G.S.B.N. (Fig. 1).     

区域地震台网地震定位能力分析

本文提出了分析区域地震台网地震定位能力的一种通用方法.在分析过程中考虑了如下事实:当改变地震震级和震源在台网区域内的位置时,能读出 P 震相和 S 震相到时的台站组合也随之改变.作为一个实例,具体分析了北京遥测地震台网的定位能力.北京台网现有19个台站,覆盖面积约为300公里400公里.根据观测资料,建立了确定每个台站检测能力(作为台站仪器放大倍数和地震震级的函数)的经验公式.然后,对每个给定震级和震源位置的可能事件,确定其子台网的构成.在这些子台网的基础上,利用奇异值分解技术,对事先设定的到时数据标准误差计算出震源座标的误差分布和线性化条件方程组的条件数分布.对给定震级(例如,ML=1.0,2.0或3.0),将计算结果绘成等值线图.为了进行比较,还分析了一个包括61个台站的北京台网扩建方案,从而可以预先估计未来台网的地震定位能力.      

Low-latitude total ozone measurements in the Brazilian sector

Regular measurements of the atmospheric ozone in the Brazilian sector were started at Cachoeira Paulista (22.7°S, 45.0°W), and Natal (5.8°S, 35.2°W) in May 1974 and November 1978, respectively. The results of the total ozone measurements carried out at these two stations up to 1981 are presented in this communication and compared with other low-and mid-latitude stations. Although Natal is an equatorial station, it presents a prominent annual variation, and the average total ozone content is high compared to satellite measurements. During 1977–78, abnormally low values of total ozone were observed at Cachoeira Paulista. Some preliminary results about the QBO 9quasi-biennial oscillation) during 1974–81 are also presented.