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.     

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.     

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.     

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.     

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.     

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.     

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.     

The gravity field and crustal thickness of Venus

The gravity and topography of Venus obtained from observations of the Magellan mission, as well as the gravity and topography from our numerical mantle convection model, are discussed in this paper. We used the hypothesis that the geoid of degrees 2–40 is produced by sublithospheric mantle density anomalies that are associated with dynamical process within the mantle. We obtained the model dynamical admittance(the geoid topography ratio based on a convection model) by a numerical simulation of the Venusian mantle convection, and used it to correct the dynamical effect in the calculation of crustal thickness. After deducting the dynamical effect, the thickness of the Venusian crust is presented. The results show that the gravity and topography are strongly correlated with the Venusian mantle convection and the Venusian crust has a significant influence on the topography. The Venusian crustal thickness varies from 28 to 70 km. Ishtar Terra, and Ovda Regio and Thetis Regio in western Aphrodite Terra have the highest crustal thickness(larger than 50 km). The high topography of these areas is thought to be supported by crustal compensation and our results are consistent with the hypothesis that these areas are remnants of ancient continents. The crustal thickness in the Beta, Themis, Dione, Eistla, Bell, and Lada regiones is thin and shows less correlation with the topography, especially in the Atla and Imdr regiones in the eastern part of Aphrodite Terra. This is consistent with the hypothesis that these highlands are mainly supported by mantle plumes. Compared with the crustal thickness calculated with the dynamical effect, our results are more consistent with the crust evolution and internal dynamical process of Venus.     

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.     

Microseismic monitoring of columnar jointed basalt fracture activity: a trial at the Baihetan Hydropower Station,China

Severe stress release has occurred to the surrounding rocks of the typically columnar jointed basalt after excavation at the Baihetan Hydropower Station, Jinsha River, China, where cracking, collapse, and other types of failure may take place occasionally due to relaxation fracture. In order to understand the relaxation fracture characteristics of the columnar jointed basalt in the entire excavation process at the diversion tunnel of the Baihetan Hydropower Station, real-time microseismic monitoring tests were performed. First, the applicability of a geophone and accelerometer was analyzed in the columnar jointed basalt tunnel, and the results show that the accelerometer was more applicable to the cracking monitoring of the columnar jointed basalt. Next, the waveform characteristics of the microseismic signals were analyzed, and the microseismic signals were identified as follows: rock fracture signal, drilling signal, electrical signal, heavy vehicle passing signal, and blast signal. Then, the attenuation characteristics of the microseismic signals in the columnar jointed basalt tunnel were studied, as well as the types and characteristics of the columnar jointed basalt fracture. Finally, location analysis was conducted on the strong rock fracture events, in which four or more sensors were triggered, to obtain the temporal and spatial evolution characteristics and laws of the columnar jointed basalt relaxation fracture after excavation. The test results are not only of important reference value to the excavation and support of diversion tunnel at the Baihetan Hydropower Station, but also of great referential significance and value to the conduction of similar tests.     

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.     

A statistical model for vesuvius and its volcanological implications

The last 300 years of Vesuvius history are reconstructed as a chronological succession of 4 phenomenological states: i) repose, ii) persistent activity, iii) intermediate eruption and iv) final eruption. It turns out that the times of permanence in each state are distributed according to the same exponential law. Vesuvius activity is then described by a Markov chain of these 4 states, with transition probabilities determined from the previous phenomenological analysis. The model reproduces the Vesuvius activity between 1694 and 1872 and possibly also in the 1872–1944 period. It turns out that, at least between 1694 and 1872, the volcano was behaving like a quasistationary system with 4 equilibrium states, perturbed by a stochastic noise responsible for occasional transitions from an equilibrium state to another. Major physical or structural changes of the volcanic system around 1872 and possibly in the whole subsequent period, are clearly shown by the statistical analysis.     

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.     

Voluminous acid volcanism in the Busheveld Complex: A review of the Rooiberg Felsite

The 2.1 b.y. old Rooiberg Felsite roofs and is intruded by the mafic layered rocks and granites of the Bushveld Complex. The felsite unit, which locally exceeds 5 km in thickness and may represent an originally erupted volume of more than 300,000 km³, is dominated by rhyolitic to dacitic lavas with minor pyroclastic and sedimentary rock types. Volcanic rocks of more mafic composition occur towards the base of the sequence. The Rooiberg episode essentially terminated the volcanic activity in the Transvaal basin and heralded the emplacement of the Bushveld Complex. Despite the close spatial and temporal relationships between the Rooiberg Felsite and the Bushveld Complex, the precise nature of the petrogenetic link is obscure. Chemical analyses of felsite have been variously interpreted to suggest cyclic differentiation along a comagmatic trend or to demonstrate anomalous enrichment in SiO₂. Several characteristics delineate the Rooiberg Felsite as a possibly unique occurrence of rhyolitic magmatism, notably the immense volume of the unit, the marked preponderance of lavas over pyroclastic types, and the unusually great thickness and lateral extent of the flows. The thesis that the Rooiberg Felsite represents a shock-produced, meteorite-impact melt cannot be supported on the available evidence.     

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.     

Ground deformations on Mt. Vesuvius from 1977 to 1981

The results of a ground deformation study (1) of Vesuvius are reported. The study is a part of the geophysical program for the surveillance of this extremely dangerous active volcano. The results obtained in the period 1977–1981 show no significant changes in the dynamic state of the volcano. Corresponding to the southern Aapennines Nov. 23, 1980 earthquake a modification of the zenithal distance between the rim and bottom of the Vesuvius crater was observed which was likely connected with fairly shallow effects.     

Volcan ceboruco: A major composite volcano of the mexican volcanic belt

Ceboruco is a major composite volcano at the western end of the Mexican Volcanic Belt, near the junction between the North American and Pacific plates. The volcano is built from successive eruptions of andesite lavas and pyroclastic rocks, and major eruptions during its history have resulted in the formation of two concentric calderas. The youngest volcanic activity has included the extrusion of dacites within the inner caldera and a voluminous flank eruption of andesite during 1870–72. Fumarolic activity persists to the present day. Chemical analyses show that the lavas are of cale-alkaline type and rangs from andesite (SiO₂=58–61%) to acid dacite (SiO₂=68%) in composition. The rate of increase of K₂O relative to SiO is greater than that in volcanic rocks from the Mexican Volcanic Belt as a whole. This indicates that simple models based on the application of such relationships may not be adequate to explain the petrogenesis of calc-alkaline lavas.     

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.     

Impact of Land Surface Processes on the South Asian Monsoon Simulations Using WRF Modeling System

The Weather Research and Forecasting model has been used to examine the role of land surface processes on Indian summer monsoon simulations. Isolated experiments have been carried out with physical parameterization schemes (land surface and planetary boundary layer) and data assimilation to examine their relative roles in the representation of regional hydroclimate in model simulations. The impact of vegetation green fraction on the model simulations has been extensively studied by replacing the default United States Geological Survey (USGS) vegetation cover data with that of Indian Space Research Organisation (ISRO) data. Results indicate that differences in the treatment of surface processes in the model lead to large differences in precipitation simulation over the Indian domain. Several hydroclimate parameters from the simulations using ISRO and USGS vegetation green fractions were examined. It is seen that the role of vegetation green fraction in these experiments has been to increase latent heat flux to the atmosphere. Two sets of data assimilation experiments were also carried out for an entire year using the same set of observed data but with different land surface parameterization schemes. It is found that evenwhen using the same observed data, the differences in land surface schemes reduce the impact and contribution of observed data being assimilated into the model. The hydroclimate over the region becomes a function of the land surface scheme. This study highlights the importance of vegetation green fraction and land surface schemes in the context of the regional hydroclimate over South Asia.     

Major-element petrochemistry of some extrusive rocks from the volcanically active Mariana Islands

Volcanic rocks from six of the currently or recently active volcances of the Mariana Island are show little variation in major element abundances. SiO₂ content averages 51.5 wt.%. The flows are high in Al₂O (mean 17.7 wt.%) and Fe oxides (mean 10.1 wt.% calculated as FeO only), and moderate in MgO content (mean 4.7 wt.%), Na₂O (mean 2.7 wt.%), and K₂O (mean 0.7 wt.%). Only the rocks from Farallon de Pajaros, the northernmost of the Mariana Islands, deviate slightly from the average of the analyses. Three analyses from this island are slightly higher in SiO₂ (about 54 wt.%) and Al₂O₃, and are lower in total Fe oxides and MgO. According to preferred classification, the lavas of the Mariana Islands can be termed mela-andesites, high-alumina basalts, or calc-alkaline (orogenic) basalts. The K₂O values (mean 0.7 wt.%) obtained from lavas of the Mariana Islands are significantly higher than the K₂O values (about 0.33 wt.%) from volcanics of the Izu chain to the north. Inasmuch as the substantial scatter in location of earthquake foci beneath both arcs prevents accurate delineation of the upper boundary of the Benioff zone, it presently cannot be determined whether this discrepancy in K₂O values reflects a difference in depth from the volcanic are to the dipping seismic zone or relates to other phenomena. The older volcanic islands within the Mariana-Bonin island chain apparently defined an island arc system during Eocene to Miocene time. This indicates that the present plane of convergence between the Pacific plate and the Philippine Sea plate has defined the convergence between these plates since Eocene time.