Expressions for the Global Gravimetric Moho Modeling in Spectral Domain

We apply a newly developed numerical method to improve the Moho geometry by the implementation of gravity data. This method utilizes expressions for the gravimetric forward and inverse modeling derived in a frequency domain. Methods for a spectral analysis and synthesis of the gravity field and crust density structures are applied in the gravimetric forward modeling of the consolidated crust-stripped gravity disturbances, which have a maximum correlation with the (a priori) Moho model. These gravity disturbances are obtained from the Earth’s gravity disturbances after applying the topographic and stripping gravity corrections of major known anomalous crust density structures; in the absence of a global mantle model, mantle density heterogeneities are disregarded. The isostatic scheme applied is based on a complete compensation of the crust relative to the upper mantle density. The functional relation is established between the (unknown) Moho depths and the complete crust-stripped isostatic gravity disturbances, which according to the adopted isostatic scheme have (theoretically) a minimum correlation with the Moho geometry. The system of observation equations, which describes the relation between spherical functions of the isostatic gravity field and the Moho geometry, is defined by means of a linearized Fredholm integral equation of the first kind. The Moho depths are determined based on solving the gravimetric inverse problem. The regularization is applied to stabilize the ill-posed solution. This numerical procedure is utilized to determine the Moho depths globally. The gravimetric result is presented and compared with the seismic Moho model. Our gravimetric result has a relatively good agreement with the CRUST2.0 Moho model by means of the RMS of differences (of 3.5 km). However, the gravimetric solution has a systematic bias. We explain this bias between the gravimetric and seismic Moho models by the unmodelled mantle heterogeneities and uncertainties in the CRUST2.0 global crustal model.     

Evidence for an Independent 26-s Microseismic Source near the Vanuatu Islands

The 26 s peak in the ambient seismic noise spectrum is persistently excited and observed at stations globally. Using noise cross-correlation functions (NCFs), the location suggests that the source could be situated in the Gulf of Guinea and Fiji Basin. However, the Fiji Basin was proposed to be the mirror site (near antipode) of the Gulf of Guinea source instead of an independent source, assuming that the surface waves more efficiently propagate along the major-arc paths of oceanic movements. To investigate the propagation of the Rayleigh waves along continental and oceanic paths, we analyzed the surface wave data recorded from an earthquake near the Gulf of Guinea and found that Rayleigh waves travel along continental minor-arc paths more efficiently than along oceanic major-arc paths. We then located the source in the western Pacific Ocean from group velocities measured with earthquake data by using the travel time misfit in NCFs after calibration and concluded that the source is in the Vanuatu Islands. Moreover, the temporal variation of the 26 s microseismic peak observed in the western Pacific seismic stations is very different from that in stations near the Gulf of Guinea, which suggests that they are excited by independent sources. Therefore, the Vanuatu source should be an independent microseismic source. As it is close to volcanoes in the Vanuatu islands, the Pacific 26 s microseismic source might be excited by magmatic processes, which are also responsible for very-long-period volcanic tremors.     

Land–Water Boundary Treatment for a Tsunami Model With Dimensional Splitting

The Method of Splitting Tsunamis (MOST) model adapted by National Oceanic and Atmospheric Administration (NOAA) for tsunami forecasting operations is praised for its computational efficiency, associated with the use of splitting technique. It will be shown, however, that splitting the computations between \(x\) and \(y\) directions results in specific sensitivity to the treatment of land–water boundary. Slight modification to the reflective boundary condition in MOST caused an appreciable difference in the results. This is demonstrated with simulations of the Tohoku-2011 tsunami from the source earthquake to Monterey Bay, California, and in southeast Alaska, followed by comparison with tide gage records. In the first case, the better representation of later waves (reflected from the coasts) by the modified model in a Pacific-wide simulation resulted in twice as long match between simulated and observed tsunami time histories at Monterey gage. In the second case, the modified model was able to propagate the tsunami wave and approach gage records at locations within narrow channels (Juneau, Ketchikan), to where MOST had difficulty propagating the wave. The modification was extended to include inundation computation. The resulting inundation algorithm (Cliffs) has been tested with the complete set of NOAA-recommended benchmark problems focused on inundation. The solutions are compared to the MOST solutions obtained with the version of the MOST model benchmarked for the National Tsunami Hazard Mitigation Program in 2011. In two tests, Cliffs and MOST results are very close, and in another two tests, the results are somewhat different. Very different regimes of generation/disposal of water by Cliffs and MOST inundation algorithms, which supposedly affected the benchmarking results, have been discussed.     

Assessment of genotoxic and humoral immune system alterations in silica exposed workers from pottery industries in South India

Chemicals, especially silica, have been suspected to cause genetic alterations in pottery industry workers. The present study aims to analyze the frequency of chromosomal aberrations (CA), micronucleus (MN) and DNA damage (comet assay) in the peripheral blood lymphocytes and the immunological alterations workers chronically exposed to silica and in control subjects. In our study, 50 silica exposed workers and 35 control subjects were recruited and silica level was measured by respirable dust and respirable quartz concentrations of badge dosimeter. The serum immunoglobulins (IgM, IgG, and IgA) of all the subjects were measured by using ELISA method. The individuals exposed to silica have a significant increase in the frequency of CA, MN and the total DNA damage (p < 0.05). Immunoglobulin elevation on silica exposed workers was statistically significant (p < 0.05) on comparison with their respective controls. Investigation of the smoking and alcohol habitats coupled with silica exposure in exposed and control subjects represents alcohol consumption and smoking as additional risk factors and must be avoided. Multiple linear regression analysis obtained for CA, MN and comet assay confirm these tests as biomarkers for silica exposed pottery workers. Some confounding factors also showed significant influence on exposed subjects. These results indicate the mutagenic risk in the working environment has a high probability of association with the silica dust exposure in pottery industries. Nevertheless, the present study will create awareness and public concern not only among the silica exposed workers but also to the welfare of their progeny.     

FEM × DEM modelling of cohesive granular materials: Numerical homogenisation and multi-scale simulations

The article presents a multi-scale modelling approach of cohesive granular materials, its numerical implementation and its results. At microscopic level, Discrete Element Method (DEM) is used to model dense grains packing. At the macroscopic level, the numerical solution is obtained by a Finite Element Method (FEM). In order to bridge the micro- and macro-scales, the concept of Representative Elementary Volume (REV) is applied, in which the average REV stress and the consistent tangent operators are obtained in each macroscopic integration point as the results of DEM’s simulation. In this way, the numerical constitutive law is determined through the detailed modelling of the microstructure, taking into account the nature of granular materials. We first elaborate the principle of the computation homogenisation (FEM × DEM), then demonstrate the features of our multiscale computation in terms of a biaxial compression test. Macroscopic strain location is observed and discussed.     

PALEOMAGIA: A PHP/MYSQL database of the Precambrian paleomagnetic data

Most paleomagnetic applications require a precise, rationally organized and up-todate catalogue or database of paleomagnetic results worldwide. These include reconstructions of continents, calculations of the Apparent Polar Wander Paths (APWPs) or paleolatitude drift curves, testing the Geocentric Axial Dipole (GAD) model, studies of geomagnetic paleosecular variation or reversal asymmetries, comparison of coeval results obtained from different types of rocks, estimation of inclination shallowing in sedimentary rocks and understanding the delay in remanence acquisition caused by slow cooling in large intrusions. For this purpose, various databases, such as the Global Paleomagnetic Database (GPMDB), and the Magnetics Information Consortium Database (MagIC) have been generated. This paper presents a new relational database (PALEOMAGIA) where 3278 entries of Precambrian data have been split geographically, sorted according to age and rock types and ranked using a revised version of the Van der Voo grading scheme. The latest geochronologic information is included wherever available. Significant effort has been put to the retrieval and archiving of data published in the last decade, which are virtually nonexistent in GPMDB. Here we present the database and its browser-based user interface from a scientific and a technical point of view.     

Recent developments of the Middle East catalog

This article summarizes a recent study in the framework of the Global Earth model (GEM) and the Earthquake Model of the Middle East (EMME) project to establish the new catalog of seismicity for the Middle East, using all historical (pre-1900), early and modern instrumental events up to 2006. According to different seismicity, which depends on geophysical, geological, tectonic, and seismicity data, this region is subdivided to nine subregions, consisting of Alborz–Azerbaijan, Afghanistan–Pakistan, Saudi Arabia, Caucasus, Central Iran, Kopeh–Dagh, Makran, Zagros, and Turkey (Eastern Anatolia; after 30° E). After omitting the duplicate events, aftershocks, and foreshocks by using the Gruenthal method, and uniform all magnitude to Mw scale, 28,244 main events remain for the new catalog of Middle East from 1250 B.C. through 2006. The magnitude of completeness (Mc) was determined as 4.9 for five out of nine subregions, where the least values of Mc were found to be 4.2. The threshold of Mc is around 5.5, 5.0, 4.5, and 4.0, for the time after 1950, 1963, 1975, and 2000, respectively. The average of teleseismic depths in all regions is less than 15 km. Totally, majority of depth for Kopeh–Dagh and Central Iran, Zagros, and Alborz–Azerbaijan, approximately, is 15, 13, and 11 km and for Afghanistan–Pakistan, Caucasus, Makran, Turkey (after 30° E), and Saudi Arabia is about 9 km.     

Geochemistry of the volcanics of central Mexico

The Tertiary and Recent volcanics of Mexico occur in two provinces. The Cordillera Province is made up of about 1700 m of ignimbrite sheets overlain and intercalated in the upper part by olivine basalt and basaltic andesite. The Rio Lerma Province extends transversely across Mexico and in the Valley of Mexico the lavas consist mainly of andesite and dacite, 68 % of those analysed having 62 | 4.7 % SiO₂. A total of 108 chemical analyses were made for the major elements, 90 of these including determinations of Cr, Ni, Cu, Zn, Rb, Sr, Ba, Th and Pb for two areas, the Valley of Mexico in the Rio Lerma Province and the Guadalajara region which lies at the intersection of the two provinces. Computer constructions of normative components in the basalt tetrahedron and other projections support an origin of partial melting of tholeiitic to pyrolitic material for the production of andesite. The Guadalajara lavas have consistently higher K/SiO₂ and K/Rb ratios and lower Mg/SiO₂ ratio than the Valley of Mexico rocks suggesting generation at greater depth.     

Design and testing of a global climate prediction system based on a coupled climate model

A global climate prediction system(PCCSM4) was developed based on the Community Climate System Model, version 4.0, developed by the National Center for Atmospheric Research(NCAR), and an initialization scheme was designed by our group. Thirty-year(1981–2010) one-month-lead retrospective summer climate ensemble predictions were carried out and analyzed. The results showed that PCCSM4 can efficiently capture the main characteristics of JJA mean sea surface temperature(SST), sea level pressure(SLP), and precipitation. The prediction skill for SST is high, especially over the central and eastern Pacific where the influence of El Ni?o-Southern Oscillation(ENSO) is dominant. Temporal correlation coefficients between the predicted Ni?o3.4 index and observed Ni?o3.4 index over the 30 years reach 0.7, exceeding the 99% statistical significance level. The prediction of 500-hPa geopotential height, 850-hPa zonal wind and SLP shows greater skill than for precipitation. Overall, the predictability in PCCSM4 is much higher in the tropics than in global terms, or over East Asia. Furthermore, PCCSM4 can simulate the summer climate in typical ENSO years and the interannual variability of the Asian summer monsoon well. These preliminary results suggest that PCCSM4 can be applied to real-time prediction after further testing and improvement.     

The Beni-Ilmane (North-Central Algeria) Earthquake Sequence of May 2010

Starting on 14 May 2010 and lasting several months, the village of Beni-Ilmane (Msila District, North-Central Algeria) and its surroundings were struck by an important seismic crisis marked by three successive moderate shocks (5.0 ≤ M _d ≤ 5.2). This sequence of events caused severe damage in the Beni-Ilmane village and in the epicentral area. The poor quality of masonry construction and the cumulative effects of the large number of aftershock events played a key role in the destruction. To follow this earthquake sequence, 11 temporary seismic stations, in addition to the permanent stations of the Algerian seismic network, were deployed in the region. A representative set of well located aftershocks in the period of maximum activity (lasting 18 days) were selected. The horizontal distribution of the aftershocks shows two main earthquake clusters located near Beni-Ilmane village, one cluster oriented E–W and the other oriented NNE–SSW, crossing the first cluster at its eastern tip. The aftershocks distribution suggests that the three main shocks ruptured two distinct and adjacent fault segments of about 8 km length. The focal mechanisms of the first and third events, located in the NNE–SSW cluster, show near-vertical left-lateral strike-slip fault planes. In the second cluster, oriented E–W, focal mechanisms show a high-angle reverse fault. A field survey, initiated immediately after the first main shock, identified surface fissures generated by the three largest events in the sequence. The fissures, concentrated in a narrow area at the western termination of the NE–SW Jebel Choukchot anticline (location of Beni-Ilmane village), showed several orientations which were mainly related to gravity instabilities. The 2010 Beni-Ilmane earthquake sequence, located in the Bibans–Hodna Mountains transition zone, demonstrates that the Tellian Atlas–High Plateaus border region is an active seismic zone marked by moderate and possibly strong earthquakes; thus, a reevaluation of the seismic hazard in the region is needed.     

A comparison of the recent eruptions of Mt. Pelée,Martinique and Soufrière,St. Vincent

The simultaneous eruption of Mt. Pelée, Martinique and Soufrière, St. Vincent are regarded as the first recognized examples of Pelean-type and St. Vincent-type pyroclastic eruptions. Both produced nuées ardentes, the former usually laterally directed because of the presence of a dome and the latter vertically directed from an open crater. Both volcanoes have subsequently erupted for a second time this century. The 1902–05 and 1929–32 eruptions of Mt. Pelée produced andesite lava of almost identical composition and mineralogy. Both contain two generations of plagioclase, orthopyroxene, Fe-Ti oxide, corroded brown amphibole and olivine rimmed by pyroxene. In contrast, the Soufrière material is more basic in composition varying from basaltic andesite to basalt in 1902–03 and basaltic andesite in 1971–72. The Soufrière material contains two generations of plagioclase (with those of 1971–72 having additional zones of labradorite), clinopyroxene, orthopyroxene, olivine and Fe-Ti oxide. The pyroclastic deposits are strikingly different, those from the Pelean-type eruption are termed «block and ash deposits» being characterised by poorly vesicular lava blocks up to 7 m in diameter, while the St. Vincent-type eruption produced «scoria and ash deposits» containing vesicular ropey blocks or bombs no larger than 1 m in diameter. The differences in styles of eruption are attributed to differences in viscosity and mechanism of eruption of the magmas. Stratigraphic studies of Mt. Pelée reveal that the volcano has produced basaltic andesite scoria and ash deposits from St. Vincent-type eruptions. It is concluded that the recent eruptions of Pelée tapped a deep level magma during both eruptions releasing magma of similar composition, while the 1971 Soufrière magma is thought to be a remnant of the 1903 basaltic magma which remained at a high level within the volcano where it underwent enrichment in plagioclase and loss of olivine and oxide.     

Determination of Dips and Depths of Near Surface Layers by Radon Transform

The refracted arrivals on seismic shot records have long been recognized as an efficient tool for the computation of detailed near-surface information. In this paper, a new concept of refraction static, which is based on the Radon transform and avoids the tedious process of picking first arrival times, is proposed. This method is particularly suitable when a rough near-surface problem necessitates the utilization of numerous shallow refraction data for the one reflector case. Quasi-linearity of refractors and a constant velocity medium are assumed within the shooting range. Synthetic and real cases have been tested to evaluate the performance of the method. The result is revealed to be satisfactory. Comparison of the synthetic model with the results obtained through the Radon transform reveals a very good accuracy for the proposed method.     

Petrology of filicudi,aeolian archipelago

Volcanic rocks on the island of Filicudi are typical of the cale-alkaline association of the Aeolian archipelago, and range in composition from basalt to high-K andesite. All specimens contain abundant phenocrystic plagioclase (cores of An_82–97, rims becoming progressively more sodic with increasing Differentiation Index), clinopyroxene, and magnetite; the most basic samples contain olivine, which is replaced by orthopyroxene at around D.I.=44. SiO₂ generally increases with time, and this in conjunction with the continuity in mineralogy suggests that the lavas are the differentiation products of a single batch of magma. Least-squares mixing calculations are consistent with a shallow fractional model, as is the variation in trace element abundances.     

Eruption of Piip Crater (Kamchatka)

In autumn of 1966 on the northern slope of Kliuchevskoy volcano a chain of new adventive craters broke out at the height of about 2200 m. Eighty-four hours before the beginning of the eruption a swarm of preliminary volcanic earthquakes had appeared. The number of preliminary shocks was 457 with total energy of 4 × 10¹⁷ erg. With the beginning of the lava flow the earthquakes stopped and a continuous volcanic tremor appeared. The total energy of volcanic tremor amounts to 10¹⁶ erg. During the eruption numerous explosive earthquakes with the energy of 10¹⁵–10¹⁶ erg were recorded and besides the microbarograph of the Volcanostation recorded 393 explosions with an energy more than 10¹³ erg and their total energy was equal to 10¹⁷ erg. All together it has been formed 8 explosive craters and the lowest 9th crater was effusive. The slag cone was formed round this effusive crater, the lava effusion of basaltic-andesite composition (52,5% SiO₂) tooke place from the lava boccas at the cone base and from the crater. The lava flow covered a distance of 10 km along the valley of the Sopochnoy river and descended to a height of about 800 m. The lava flow velocity at the outflow reached 800 m/hr, the lava temperature was 1050°C. The effused lava volume amounts to 0.1 km³. The eruption stopped on December 25–26, 1966.     

Calc-alkaline volcanic rocks from NW Sardinia: Evaluation of a fractional crystallization model

The volcanic centre of Monte Seda Oro, N. W. Sardinia, representative of a Cenozoic calc-alkaline andesitic suite of rocks is composed of a variety of rocks ranging from high alumina basalts to dacites. The minerals of basaltic, andesitic and dacitic rocks show only limited variation in chemical composition. The geochemical data suggest that the various rock-types are related by a crystal-liquid fractionation. Least-squate numerical calculations, using major element data, support the derivation of andesites with SiO₃ content ranging from 53.8 to 59.0% from basalts having about 48.7% of SiO₂ by low pressure crystal fractionation of the phenocryst phases present in these rocks. However, the origin of dacites cannot be readily explained by this mechanism.     

Volcanic ash hazard in the Central Mediterranean assessed from geological data

Volcanic ash produced during explosive eruptions can have very severe impacts on modern technological societies. Here, we use reconstructed patterns of fine ash dispersal recorded in terrestrial and marine geological archives to assess volcanic ash hazards. The ash-dispersal maps from nine Holocene explosive eruptions of Italian volcanoes have been used to construct frequency maps of distal ash deposition over a wide area, which encompasses central and southern Italy, the Adriatic and Tyrrhenian seas and the Balkans. The maps are presented as two cumulative-thickness isopach maps, one for nine eruptions from different volcanoes and one for six eruptions from Somma-Vesuvius. These maps represent the first use of distal ash layers to construct volcanic hazard maps, and the proposed methodology is easily applicable to other volcanic areas worldwide.     

Erosion and aggradation on persistently active volcanoes—a case study from Semeru Volcano,Indonesia

Erosion processes on active volcanoes in humid climates result in some of the highest sediment yields on Earth. Episodic sediment yields after large eruptions have been evaluated, but not the long-term and continuous patterns on persistently active volcanoes. We have used high-spatial resolution satellite imagery and DEMs/DSMs along with field-based geologic mapping to assess accurately sediment budgets for the active Semeru Volcano in Java, Indonesia. Patterns of aggradation and degradation on Semeru differ from that of other active volcanoes because (1) both episodic pyroclastic density currents (PDC) and continuous supplies of tephra generate pulses of sediment, (2) sediment is transferred via cycles of aggradation and degradation that continue for >15 years in river channels after each PDC-producing eruption, and (3) rain-triggered lahars remove much greater material than fluvial transport during long, intense rainfall events. The geomorphic response of two of Semeru’s rivers to volcanic sediment migration indicates that (1) each river experiences alternating aggradation and degradation cycles following PDC-producing eruptions and (2) spatial patterns of sediment transfer are governed by geomorphic characteristics of the river reaches. Usually high degradation in the steep source reach is followed by a long bypassing middle reach. Aggradation predominates in the depositional reaches further down valley on the ring plain. Average sediment yields (10³–10⁵ t/km²/year) at persistently active volcanoes are two to three orders of magnitude lower than sediment yields after large and infrequent eruptions, but the continuous and steady sediment transfer in rivers removes more sediment on a mid-term (10 years) to long-term (30 years) basis. In contrast to the trend observed on composite cones after large and infrequent eruptions, decay of sediment yields is not exponential and river channels do not fully recover at steadily active volcanoes as episodic inputs from BAF eruptions, superimposed on the background remobilization of daily tephra, have a greater cumulative effect.     

Explosive felsic volcanism on El Hierro (Canary Islands)

The Canary Islands consist of seven basaltic shield volcanoes whose submerged portion is much more voluminous than the subaerial part of each island. Like so many other volcanic oceanic islands, the indicative deposits of explosive felsic volcanism are not a common feature on the Canary archipelago. Hitherto, they have only been documented from the central islands of Gran Canaria and Tenerife, which are the largest volcanic complexes of the islands. On the other Canary Islands, the presence of felsic rocks is mostly restricted to intrusions and a few lava flows, generally within the succession in the oldest parts of individual islands. In this paper, we present a detailed stratigraphic, lithological and sedimentological study of a significant felsic pumice deposit on the island of El Hierro, referred here as the Malpaso Member, which represents the only explosive episode of felsic volcanism found on the Canary Islands (outside of Gran Canaria and Tenerife). The products of the eruption indicate a single eruptive event and cover an area of about 15 km². This work provides a detailed stratigraphic and chronological framework for El Hierro, and four subunits are identified within the member on the basis of lithological and granulometric characteristics. The results of this study demonstrate the importance of an explosive eruption in a setting where the activity is typified by effusive basaltic events. Given the style and the spatial distribution of the Malpaso eruption and its products, a future event with similar characteristics could have a serious impact on the population, infrastructure and economy of the island of El Hierro.     

Petrology of tholeiitic lavas from Tanna Island (New Hebrides): Importance of cumulative processes in Island arc magmatism

Tanna, one of the southernmost islands of the New Hebrides volcanic arc, is made of Late Pliocene to Recent island arc tholeiitic basalts and andesites, with SiO₂ contents ranging from 45 to 57%. These lavas are highly porphyritic (30–50% in volume): phenocrysts of plagioclase are the most abundant, together with olivine and clinopyroxene. The groundmass contain plagioclase, augite, olivine, magnetite and glass; pigeonite, tridymite, sanidine and, rarely, biotite may also occur. The olivines and clinopyroxenes show an iron enrichment from the cores of phenocrysts to their rims and the groundmass crystals, but their compositional variations are not correlated with the Mg/Fe ratio of bulk host rocks, the most Fe-rich compositions being found in Mg-rich lavas. Plagioclase compositions range from An₉₅ to An₆₀ in the basalts and An₆₀ to An₅₀ in the andesites, but, within each group, they are not correlated with SiO₂ or Na₂O contents of host lavas. Consequently, the bulk major element compositions of Tanna volcanic rocks cannot be considered as primarily controlled by crystal separation from successive liquids. The oxyde-SiO₂ variations diagrams, and the modal compositions and mineral chemistry show that crystal accumulation is the predominant mechanism accounting for bulk rock compositions. However, this does not exclude fractional crystallization: the variation of the calculated groundmass mineralogy strongly suggest the occurrence of crystal removal mainly clinopyroxene and magnetite.     

The determination of deep temperatures by means of the CO-CO2-H2-H2O geothermometer: an example using fumaroles in the Campi Flegrei,Italy. A comment

Three points raised in the paper by Tedesco and Sabroux (1987) are dealt with. (1) The inconsistency between the water partial pressure calculated by Tedesco and Sabroux (1987) and saturation pressure is due to the improper use of the water-gas-shift reaction as a geothermometer. In fact Tedesco and Sabroux (1987) do not take into account the distribution of gas species between the coexisting vapour and liquid phases. (2) The depth of the “steam reservoir” is evaluated by Tedesco and Sabroux (1987) in too simplistic a way. This matter should be treated with greater care owing to the high social impact of any consideration on the Phlegraean Fields system. (3) The reliability of carbon monoxide determination at the concentration level encountered at Solfatara depends on the collection method rather than on the gas-chromatographic technique.