Results of locating the IASPEI GT(0-5) reference events using the standard ISC procedures

The International Seismological Centre (ISC) is charged with production of the definitive global bulletin of seismic events, based on the most comprehensive set of parametric data collected from all over the world. Almost every event in the bulletin retains the original hypocentral solutions reported to the ISC by contributing agencies. In addition, where possible, the ISC computes its own solution, which is intended to be the most accurate where the data from several networks are used. It is because of the requirement for consistency of the bulletin over the years that the procedures used at the centre to compute hypocentres have remained rather conservative despite considerable advances made in the field of earthquake location.The ISC has developed and put into operation a new data management system. As a result, it is now possible to review and subsequently introduce more up-to-date methods of locating seismic events into the operations. The ISC Governing Council called for a workshop dedicated to location procedures, which was held during the 2005 IASPEI General Assembly in Santiago, Chile.To compare the accuracy of different location algorithms, a list of 156 reference events (IWREF) was selected prior to the workshop. The list includes geographically well distributed earthquakes and explosions with positions known with an accuracy of up to 5 km. It covers the period of 1954-2001 and includes all station readings and hypocentral solutions of different agencies available for these events in the ISC bulletin. Although the original ISC solutions are included, these may be different from the solution obtainable now due to changes in the ISC procedures over the years. This paper presents the results of relocation of these events using standard ISC location procedures as of 2005. These new ISC locations and analysis of their shifts with respect to reference locations present a benchmark for further improvement.     

A Bayesian Approach to Seismic Hazard Estmation: Maximum Values of Magnitudes and Peak Ground Accelerations

Earthquake Researeh in Ch一na461 .METHODLet R be some value measured or estimated as a sequenee ina’‘Past”time interval(一丁,O)(I)万‘月,=(RI,…,R。),R,之R。,R=nlaX l二f匕11(RI,…,R,,) Values(l)eould have an arbitrary Physieal nature.BelowweshalleonsiderEq.(l)asearthquakemagnitudes in a given seismic aetive region or logarithms of seismie Peak ground aeeelerations at习given site.Ro isa而nimum eutoff value;it 15 defined by Possibilities of registration systems or wasehosen as the …     

海洋天然气水合物的地球物理研究(Ⅱ):地震方法

从野外观测手段、数据处理分析、振幅随炮检距变化分析、波形反演4个方面综述海洋天然气水合物地震研究的最新进展。     

Using an Artificial Intelligence Method to Recognize Seismic Precursory Anomaly (SPA)

In this work,we seek for a new way to recognize and analyze the Seismic Precursory Anomaly(SPA)quantitatively.For further study of the mechanism of SPA,an artificial intelligence method,the Machine Discovery(MD)method,is proposed to process the Seismic Precursor Observation(SPO)data and to fit an analytic function to them.The analytic function is used as a criterion to identify the SPA.The algorithm of MD is worked out based on the theory of Recursive Function(RF).The data of Amplitude Factor of Solid Tide(7 value)are processed by the proposed method.We find that the subsequent data fit to the extrapolation of the discovered function during the quiet period,whereas they deviate from it before a strong earthquake.Finally,the feasibility of quantitative analysis of the SPA is testified.     

High precision locations of multiplets on south-eastern flank of Mt. Etna (Italy): reconstruction of fault plane geometry

On 9 January 2001, a seismic swarm, located on the south-eastern flank of Mt. Etna and with nearly identical waveforms, caused some damage to Zafferana Etnea village, 3 km from the epicentral area.An analysis of the seismicity occurring in the last 8 years in this area has revealed other earthquakes with the same characteristics; some pre-empted and followed (up to a few months) the 2001 January swarm, others were recorded more than five years beforehand.Using similarity of waveforms, these earthquakes were classified into three families.The use of a multiplet-technique has allowed to obtain the spatial distribution of the events with higher precision (mean error of 10-20 m) with respect to traditional localization techniques.Mt. Etna earthquakes relocation clearly describes the geometry of the seismogenic tectonic structure; the hypocenters lie on a NE-SW oriented plane that is coincident with one of the focal planes obtained by first-arrival polarities. This alignment is also coherent with one of the main regional tectonic trends cutting the Mt. Etna area, and can be interpreted as a right-lateral strike seismic source on the south-eastern flank of Mt. Etna, distant from eruptive centres, which repeats from time to time and is able to produce strong energy releases.     

Flow Estimation in a Basin by Using a Hydrological Model: Application to the Basin of Wadi Djendjen (Jijel,North-East Algeria)

The aim of this study was to transform the rain into flow in the basin of Wadi Djendjen by using a simulation model in absence of the continuous measurements of the flow. The obtained results indicated and confirmed the strong adaptability of this model with varied hydroclimatic situations. The validation process showed an interesting result, which lead to conclude that the model is well calibrated and has a good performance to be used for the basins with Mediterranean climate. The obtained results for the simulation by the GR2M model over the studied period (1973–1986) showed that there no significant difference between the obtained value for the yearly average flow (197.70 hm³) and that measured at Missa hydrometric station (200.80 hm³), which allow us to estimate quantitatively the flow in Missa hydrometric station. While, for the period of 2000–2012, the results showed that the yearly average flow value (171.90 hm³), is significantly different (reduction of 15% for 12 years), in comparison with that measured at Missa hydrometric station (200.80 hm³). This can be due to the dryness which struck the region since 2000, and the local degradation of the forest ecosystem, which has considerably affected the runoff surface.     

A simplified approach for vulnerability assessment in moderate-to-low seismic hazard regions: application to Grenoble (France)

Due to the moderate seismic risks in France, the building vulnerability assessment methods developed for high seismic risk countries could not easily be used here because of their cost and the low-risk perception among the public and officials. A light vulnerability assessment method is proposed and tested in Grenoble (France), based on classes and scores provided in the GNDT method but simplified in terms of visual screening and number of structural parameters used. Compared to the RiskUE method, the damage obtained by our approach shows that 90% of buildings have residuals smaller than 0.2, i.e. one grade of the EMS98 damage scale. A large scale survey is devised and conducted among the inhabitants of Grenoble in order to collect the main structural parameters. By comparing the results from the survey to the historical urbanization of Grenoble and to expert surveys performed in two urban districts, the information useful for the light method of vulnerability assessment can be rapidly collected by non-experts reducing substantially the estimate cost. The average damage is then computed using the GNDT formula considering the probable intensities which could be observed in Grenoble (VII and VIII). The average damage reaches 0.4 in the oldest part of Grenoble mainly made of masonry buildings and 0.2 in reinforced concrete suburbs where reinforced concrete predominates. The results are a relative vulnerability assessment that provides useful initial information for the urban zones of Grenoble where the vulnerability is higher. This method can be used to classify the seismic vulnerability in wide seismic-prone regions to a fair degree of accuracy and at low cost.     

Monte Carlo Inversion of DInSAR Data for Dislocation Modeling: Application to the 1997 Umbria-Marche Seismic Sequence (Central Italy)

— The study of surface deformation due to seismic activity is often made using dislocations with uniform slip and simple geometries. A better modeling of coseismic and postseismic surface displacements can be obtained by using dislocations with variable slip and nonregular shapes. This is consistent with the asperity model of fault surfaces, assuming a friction distribution on faults made of locked zones with much higher friction than surrounding zones. In this paper we consider the 1997–1998 Colfiorito seismic sequence. The coseismic surface displacements in the Colfiorito zone are used in order to infer the slip distribution on the fault surface at different stages of the sequence. The displacement field has been modeled varying the slip distribution on the fault, and comparing the deformation observed by SAR and GPS techniques with model results. The slip distribution is calculated by Monte Carlo simulations on a normal fault with the dip angle equal to 40°. A good approximation is obtained by using square asperity units of 1.5×1.5 km². In the first stage, we employed a simplified model with uniform slip, in which each asperity unit is allowed to slip a constant amount or not to slip at all, and in the second stage, we evaluate the slip distribution in the dislocation area determined by the Monte Carlo inversion: in this case we allow unit cells to undergo different values of slip in order to refine the initial dislocation model. The results show that the 1997 seismic events of the sequence can be modeled by irregular dislocations, obtaining a good fit to the DInSAR and GPS observations. The model also confirms the results of previous studies by a different methodology, defining the distribution of asperities on the fault plane using the fault geometry, the geodetic data and the seismic moment of the 1997–1998 Colfiorito seismic sequence. Furthermore, the analysis of 1997 aftershocks in the seismogenic region shows a strong correlation between most events and the asperity distribution, which can be considered as an independent test of the validity of the model.     

Using SPIRAL (Single Pollen Isotope Ratio AnaLysis) to estimate C3- and C4-grass abundance in the paleorecord

C₃ and C₄ grasses differ greatly in their responses to environmental controls and influences on biogeochemical processes (e.g. water, carbon, and nutrient cycling). Difficulties in distinguishing between these two functional groups of grasses have hindered paleoecological studies of grass-dominated ecosystems. Stable carbon isotopic analysis of individual grains of grass pollen using a spooling-wire microcombustion device interfaced with an isotope-ratio mass spectrometer holds promise for improving C₃ and C₄ grass reconstructions. This technique, SPIRAL (Single Pollen Isotope Ratio AnaLysis), has only been evaluated using pollen of known C₃ and C₄ grasses. To test the ability of SPIRAL to reproduce the abundance of C₃ and C₄ grasses on the landscape, we measured δ¹³C values of > 1500 individual grains of grass pollen isolated from the surface sediments of ten lakes in areas that span a large gradient of C₃- and C₄-grass abundance, as determined from vegetation surveys. Results indicate a strong positive correlation between the δ¹³C-based estimates of % C₄-grass pollen and the abundance of C₄ grasses on the landscape. The % C₄-grass pollen slightly underestimates the actual abundance of C₄ grasses at sites with high proportions of C₄ grasses, which can be corrected using regression analysis. Comparison of the % C₄-grass pollen with C/N and δ¹³C measurements of bulk organic matter illustrates the distinct advantages of grass-pollen δ¹³C as a proxy for distinguishing C₃ and C₄ shifts within the grass family. Thus SPIRAL promises to advance our understanding of grassland ecology and evolution.     

Time Evolution of Deformation Using Time Series of Differential Interferograms: Application to La Palma Island (Canary Islands)

Differential interferometry is a very powerful tool for detecting changes in the Earth’s crust where coherence conditions are good, but is difficult to employ in some volcanic areas due to dense vegetation. We apply two differential InSAR methods using the time series associated with the interferograms to perform a phase analysis on a data set for La Palma island (Canary Islands) from the ERS-1 and ERS-2 European Space Agency (ESA) satellites for the time period 1992 to 2000. Both methods involve choosing a master image from the database and creating a series of interferograms with respect to this image. The “Coherent Pixel Time Series” (CPTS) technique chooses pixels with good average coherence, aligns the unwrapped interferograms with a stable area and then performs an inversion to calculate the linear velocity to quantify the deformation. The Coherent Target Modeling (CTM) method calculates the temporal coherence of each pixel to identify stable targets and then determines the best velocity for each pixel by using a linear fit that maximizes the temporal coherence. Using these two methods we have been able to detect deformation on La Palma Island that has been previously undetectable by conventional InSAR methods. There is a roughly circular region on the Southern part of the island that is actively deforming at ~ −4 to −8 mm/yr. This region is located near the Teneguia valcano, the host of the last known eruption on La Palma in 1971. A thorough investigation of the possible sources for this deformation revealed that it was most likely created by a subsurface thermal source.     

The transition from Alboran to Algerian basins (Western Mediterranean Sea): Chronostratigraphy,deep crustal structure and tectonic evolution at the rear of a narrow slab rollback system

The eastern Alboran basin and its transition to the Algerian basin is a key area in the Mediterranean realm where controversial kinematic and geodynamical models are proposed. Models imply striking differences regarding the nature of the crust, the prevalence of brittle faulting and ductile shear, the origin of magmatism, the style of Miocene deformation and the driving mechanisms of the Alboran plate kinematics. Combining a new chronostratigraphic chart of the Alboran and Algerian basins based on the Habibas (HBB-1) core drill, deep seismic sections striking WSW-ENE and SSE-NNW, and potential field data, we re-assess the tectonic evolution that controlled the sedimentation and basement deformation of the westernmost limit of the Algerian basin and its transition with the Alboran domain. A WSW-directed extensional tectonic phase has shaped a stretched continental crust with typical tilted blocks along ∼100 km from Burdigalian to Tortonian times, which is assumed to result from the WSW-directed migration of the Alboran block driven by a narrow slab rollback. In the Algerian basin, this event was followed by the emplacement of an oceanic-type crust. Potential field signatures of the deep basin as well as geometrical correlations with onland outcrops of inner zones suggest a minimum WSW-directed displacement of the Alboran terrane of ∼200 km. At the southern foot of the Algerian basin, the continent-ocean transition is sharp and may result from the westward propagation of a slab tear at depth, forming two segments of STEP (Subduction-Transform Edge Propagator) margins. Our results support models of intense shear tractions at the base of an overriding plate governed by slab rollback-induced mantle flow. Finally, Messinian salt tectonics affected overlying deposits until today. A late Tortonian to Quaternary dominantly transpressive tectonic episode linked to the Africa-Iberia convergence post-dates previous events, deforming the whole margin.     

Biomonitoring of Heavy Metals Using the Bivalve Molluscs in Sunderban Mangrove Wetland,Northeast Coast of Bay of Bengal (India): Possible Risks to Human Health

The suitability of using four bivalve molluscs (Sanguinolaria acuminata, Anadara granosa, Meretrix meretrix, and Pelecyora trigona) in biomonitoring of heavy metals (Cu, Pb, Cd, Zn, and Hg) collected from intertidal regions of the Sunderban mangrove wetland, northeastern part of the Bay of Bengal, were evaluated. Both speciesdependent variability and temporal variations were pronounced. A high degree of organ specificity was evident in the bivalves where gill and mantle exhibited higher metal accumulation due to ion exchange property of the mucous layer covering these organs while shells represent very poor accumulation. Elevated values of Zn and Cu reflect high potential for biomagnification through marine food chain. Metal concentrations in different body size groups of the bivalves do not follow uniform trend. Correlation coefficient between different metal couplings as tested statistically revealed significant coupling for Pb–Zn, Pb–Cu, Zn–Cu, and Hg–Cu. Concentrations of all the metals in specific organs (visceral mass, mantle and gill) of the bivalves exceeded the safe levels according to the international standards for metals compiled by Food and Agricultural Organization of the United Nations and would be of great risk for human consumption. It is concluded that the mussel and clam are suitable biomonitors to employ in programs designed to assess changes in metal pollution in the Sunderban mangrove wetland.     

Using hydrologic and hydraulically derived geometric parameters of perennial rivers to determine minimum water requirements of ecological habitats (case study: Mazandaran Sea Basin—Iran)

This paper studies one of the most important problems of dry countries that are confronted with water deficit and the competition of rivals to allocate water. Some common methods have been investigated for computing the minimum water requirement to save a river's biological activity. After a discussion of the currently used method in Iran (the Tenant Method), the application of some other methods, which are known as Hydraulic and Hydrological Methods, is illustrated. The case study is a river in the northern part of Iran and this research addresses the critical situation of this river in near future regarding the planned anthropogenic alteration and its consequences. It has been shown that the application of environmental water allocation methods that have no background in a region could be misleading. The first proposed method is the Texas Method, in which flexibility in water allocation helps to develop an integrated river management paradigm in the study area. The second preferred method is a Hydraulic Method, by which the implementation of morphological parameters or flow geometrical properties could sustain physical habitat within an acceptable range in terms of depth, width, velocity, and bed shear stress. In the case study, the Maximum Curvature Method was superior to the Slope Method. The investigation revealed that using a widely recommended slope of 1 for the discharge‐wetted perimeter function can lead to an overestimated and unacceptable discharge. The Tenant Method in respect to minimum environmental flow requirement yielded the weakest result, and it has been illustrated that its application might impose irrecoverable shock to the ecosystem. The Flow Duration Curve Method (the Q₉₅ Method), in spite of its subjectivity, showed more compatibility with the river's condition in comparison with the Tenant Method. Copyright © 2011 John Wiley & Sons, Ltd.