Evaluation of magmatic processes for the products of the Nevado del Ruiz Volcano, Colombia from geochemical and petrological data

Trace-element concentrations of the eruptive products from the Nevado del Ruiz volcano indicate that at least for the last 1 Ma, the major controlling factor in the evolution of the magma has been a simple fractional crystallization.The Zr/Hf ratios suggest that the source material is mantle-like in origin and whilst subduction-related contamination of the source material has taken place, there is no evidence that assimilation of continental crust has affected the series.No convincing geochemical evidence for magma mixing was found, and data on banded pumices of Pleistocene, Historic and the 1985 eruptions strongly support a crystal fractionation model to explain the compositional variation of their glasses. Therefore, if mixing has taken place it can only be a mixing of layers within a stratified or partially stratified magma chamber.Modeling of fractional crystallization using a stepwise program has been undertaken and least-squares approximations with small residuals are commensurate with analyzed samples. Crystal fractionation models with plagioclase dominating two-pyroxenes and iron oxides give good fits for all calculations. We conclude that simple crystal fractionation is the process most consistent with major- and trace-element variations in the Ruiz series. This closed-system model requires at least 77% crystal fractionation of a basaltic parent to generate the observed compositions of dacites.     

Geophysical investigations of the Platanares geothermal site, Honduras, Central America

Geophysical studies were carried out in Platanares in Western Honduras, as part of an investigation of the geothermal resource potential of Honduras by Los Alamos National Laboratory in cooperation with the Honduran government and the U.S. Geological Survey. Geological and geochemical studies later followed by exploratory bore-holes indicate this area of rugged topography has the best geothermal potential in Honduras. In order to better understand the geology of the area, gravity, audiomagnetotelluric and telluric surveys were carried out. The site is in a graben of Tertiary volcanics and Cretaceous red beds, with numerous hot springs along faults in the center of the graben. Gravity mapped the thickness of low-density tuffs and the position of associated normal faults. It also indicated the presence of persuasive, shallow, highly silicified, high-density tuffs and gravels. The electrical surveys mapped shallow, high-resistivity silicified zones and faults as low-resistivity zones. The red beds below the tuffs are considered to contain the geothermal reservoir, but their maximum thickness could not be determined by these surveys due to the lack of density contrast and the high surface electrical conductivities. Three geologic cross-sections constructed by others have been re-evaluated using gravity modeling and analysis of eletrical anomalies, showing general agreement but in detail significant differences.     

Simultaneous iterative time-domain sparse deconvolution to teleseismic receiver functions

In this paper, we present a new approach to estimate high-resolution teleseismic receiver functions using a simultaneous iterative time-domain sparse deconvolution. This technique improves the deconvolution by using reweighting strategies based on a Cauchy criterion. The resulting sparse receiver functions enhance the primary converted phases and its multiples. To test its functionality and reliability, we applied this approach to synthetic experiments and to seismic data recorded at station ABU, in Japan. Our results show Ps conversions at approximately 4.0 s after the primary P onset, which are consistent with other seismological studies in this area. We demonstrate that the sparse deconvolution is a simple, efficient technique in computing receiver functions with significantly greater resolution than conventional approaches.     

The Byers Basin: Jurassic-Cretaceous tectonic and depositional evolution of the forearc deposits of the South Shetland Islands and its implications for the northern Antarctic Peninsula

ABSTRACT

This paper addresses the Jurassic–Cretaceous stratigraphic evolution of fore-arc deposits exposed along the west coast of the northern Antarctic Peninsula. In the South Shetland Islands, Upper Jurassic deep-marine sediments are uncomformably overlain by a Lower Cretaceous volcaniclastic sequence that crops out on Livingston, Snow and Low islands. U-Pb zircon ages are presented for the upper Anchorage Formation (153.1 ± 1.7 Ma) and the Cape Wallace granodiorite of Low Island (137.1 ± 1.7 Ma) as well as ⁴⁰Ar/³⁹Ar ages of 136–139 Ma for Low Island andesites. Data are also presented for a U-Pb age of 109.0 ± 1.4 Ma for the upper volcanic succession of Snow Island. In combination with published stratigraphy, these data provide a refined chrono- and litho-stratigraphic framework for the deposits herein referred to as the Byers Basin. Tentative correlation is explored with previously described deposits on Adelaide and Alexander islands, which could suggest further continuation of the Byers Basin towards the south. We also discuss possible correlation of the Byers Basin with the Larsen Basin, a sequence that shows the evolution of foreland to back-arc deposits more or less contemporaneously with the fore-arc to intra-arc evolution of the Byers Basin.     

Non-Linear Gravitational Clustering Using the Zel'dovich Approximation

As an alternative to computationally expensive N-body simulations for gravitional clustering, the Zel'dovich approximation (ZA) was studied in 3D, 2D and 1D. Plots of the density contrast were compared against linear theory and the exact solution. The ZA was found to perform very well in the linear regime, better than linear theory, and to give a good approximation well into the non-linear regime. This revised version was published online in July 2006 with corrections to the Cover Date.     

Enhanced groundwater vulnerability assessment in geological homogeneous areas: a case study from the Argentine Pampas

The southeast area of the Argentine Pampas is characterized by the presence of an unconfined aquifer in a wide plain. A methodology is proposed that deals with the aquifer vulnerability where the homogeneity of the hydrogeological variables used by traditional methods (in this case, DRASTIC-P) causes vulnerability maps to show more than 80% of the territory under the same class. This absence of discrimination renders vulnerability maps of little use to decision-makers. In addition, the proposed methodology avoids the traditional vague classification (high, low, and moderate vulnerability) which is highly dependent on subjectivity in its association of each class with hydrogeological considerations. That traditional vulnerability assessment methodology was adapted using a geographic information system to reclassify classes, based on the Natural Breaks (Jenks) method. The pixel-to-pixel comparison between the result obtained by the DRASTIC-P and the reclassified classes generates the so-called operational vulnerability index (OVI), which shows four classes, associating each with different hydrogeological requirements to make decisions.     

Catastrophic,rainfall-induced debris flows in Andean villages of Tarapacá, Atacama Desert,northern Chile

In March 2012, during the rainy season in the Altiplano plateau, a >100-year return period rainfall event affected the deeply incised valleys of the Precordillera of the Tarapacá Region, northern Chile. This extreme event in a very arid region triggered a number of debris and mud flows that caused severe damage and destruction in several small villages along the Camiña and Tarapacá valleys. The highly vulnerable location of the villages on top of alluvial fans due to socioeconomic and cultural reasons is a key factor to explain the level of destruction in most villages. In this paper, this unusual, remarkable landslide event is described, and the hazard faced by these settlements for future rainfall episodes and possible mitigation measures are discussed.     

Spectral Line and Continuum Radiation Propagation in a Clumpy Medium

We discuss the propagation of spectral line and continuum radiation in a clumpy medium and give general expressions for the observed absorption or emission from a cloud population. We show that the affect of the medium clumpiness can usually be characterised by a single number multiplying the mean column opacity. Our result provides a simpler proof and generalization of the result of Martin et al. (1984). The formalism provides a simple way to understand the effects of clumping on molecular line profiles and ratios; for example, how clumping effects the interpretation of ¹³CO(1–0) to ¹²CO(1–0) line ratios. It also can be used as a propagation operator in physical models of clumpy media where the incident radiation effects the spectral line emissivity. We are working to extend the formalism to the propagation of masers in a clumpy medium, but in this case, there are special difficulties because formal expectation values are not characteristic of observations because they are biased by rare events.     

A GPS-Supported Gravity Survey in the Amazon of Ecuador

The growing importance and effectiveness of the Global Positioning System (GPS) in positioning geophysical surveys is now well established. In Ecuador's Amazon basin the use of GPS was the determining factor in the successful completion of a gravity survey carried out in 1994. Over 500 gravity stations, spaced from 1 to 5 km apart and positioned by GPS, were collected in the jungle environment in the country's southeast, where the preexisting data, of variable and unknown quality, were mostly confined to areas along the major rivers. This station densification, which includes the reoccupation of selected points, contributed to increase the gravity coverage; it also allowed us to evaluate and correct 1100 preexisting data points in the area, so that the different surveys can be merged. The composite data set is used for a geologic analysis of the study area, where previously unknown structures are discovered. ? 1999 John Wiley & Sons, Inc.     

Three‐term inversion of prestack seismic data using a weighted l2, 1 mixed norm

We present a new inversion method to estimate, from prestack seismic data, blocky P‐ and S‐wave velocity and density images and the associated sparse reflectivity levels. The method uses the three‐term Aki and Richards approximation to linearise the seismic inversion problem. To this end, we adopt a weighted mixed l_{2, 1}‐norm that promotes structured forms of sparsity, thus leading to blocky solutions in time. In addition, our algorithm incorporates a covariance or scale matrix to simultaneously constrain P‐ and S‐wave velocities and density. This a priori information is obtained by nearby well‐log data. We also include a term containing a low‐frequency background model. The l_{2, 1} mixed norm leads to a convex objective function that can be minimised using proximal algorithms. In particular, we use the fast iterative shrinkage‐thresholding algorithm. A key advantage of this algorithm is that it only requires matrix–vector multiplications and no direct matrix inversion. The latter makes our algorithm numerically stable, easy to apply, and economical in terms of computational cost. Tests on synthetic and field data show that the proposed method, contrarily to conventional l₂‐ or l₁‐norm regularised solutions, is able to provide consistent blocky and/or sparse estimators of P‐ and S‐wave velocities and density from a noisy and limited number of observations.