Transient processes in Stromboli’s shallow basaltic system inferred from dolerite and magmatic breccia blocks erupted during the 5 April 2003 paroxysm

We describe the mineralogy, geochemistry, and mesomicrostructure of fresh subvolcanic blocks erupted during the 5 April 2003 paroxysm of Stromboli (Aeolian Islands, Italy). These blocks represent ∼50 vol.% of the total erupted ejecta and consist of fine- to medium-grained basaltic lithotypes ranging from relatively homogeneous dolerites to strongly or poorly welded magmatic breccias. The breccia components are represented by angular fragments of dolerites entrapped in a matrix of vesiculated (lava-like to scoriae) crystal-rich (CR) basalt. All of the studied blocks are cognates with the CR basalt of the normal Strombolian activity or lavas and they are often coated by a few-centimeter thick layer of crystal-poor (CP) basaltic pumice erupted during the paroxysm. We suggest that they result from the rapid increase of pressure and related subvolcanic rock failure that occurred shortly before the 5 April 2003 explosion, when the uppermost portion of the edifice inflated and suffered brecciation as the result of the sudden rise of the gas-rich CP basalt that triggered the eruption. Dolerites and magmatic matrix of the breccias show major and trace element compositions that match those of the CR basalts erupted during normal Strombolian activity and effusive events at Stromboli volcano. Dolerites consist of (a) phenocrysts normally found in the CR basalts and (b) late-stage magmatic minerals such as sanidine, An_60-28 plagioclase, Fe–Mn-rich olivines (Fo_68-48), phlogopite, apatite, and opaque mineral pairs (magnetite and ilmenite), most of which are never found both in lava flows and scoriae erupted during the persistent explosive activity that characterizes typical Strombolian behavior. Subvolcanic crystallization of the Stromboli CR magma, leading to slowly cooled equivalents of basalts, could result from transient drainage of the magma from the summit craters to lower levels. Fingering and engulfing of the material that collapsed from the summit crater floor into the shallow basaltic system during the late evening of 28 December 2002 coupled with the short break in the summit persistent explosions between December 2002 and March 2003 permitted the CR magma pockets to solidify as dolerites, which were confined to the uppermost portion of the system and thus not involved in the ongoing flank effusive activity. Crystal size distribution of the basaltic blocks and crystallization of the finer-grained (<0.1 mm) mafic minerals of the dolerites over a time interval of ∼100 days closely agrees with the above interpretation. Vesicle filling (miarolitic cavities) locally found in some dolerites, with minerals deposited as vapor-phase crystallization is a result of continuous gas percolation through the rocks of the uppermost portion of the volcanic system. Poorly welded magmatic breccias formed during syn-eruptive processes of 5 April 2003, when the paroxysm strongly shattered the shallow subvolcanic system and many dolerite fragments were entrapped in the CR magma. In contrast, the high degree of welding between the dolerite clasts and the CR basaltic matrix in the strongly welded magmatic breccias provides a snapshot of subvolcanic intrusions of the CR basalt into the dolerite when, after a 2-month break in activity, CR magmas started to rise again to the summit craters. Blocks similar to these subvolcanic ejecta of 5 April 2003 were also erupted during previous paroxysms (e.g., 1930) suggesting that changes in the usual Strombolian activity (e.g., short breaks in the persistent mild explosions and/or flank effusive activity) lead to transient crystallization of dolerites in the shallow plumbing system.     

Regularized tomographic inversion with geological constraints

Reflection tomography is the industry standard tool for velocity model building, but it is also an ill‐posed inverse problem as its solution is not unique. The usual way to obtain an acceptable result is to regularize tomography by feeding the inversion with some a priori information. The simplest regularization forces the solution to be smooth, implicitly assuming that seismic velocity exhibits some degree of spatial correlation. However, velocity is a rock property; thus, the geometry and structure of rock formations should drive correlation in velocity depth models. This observation calls for constraints driven by geological models. In this work, we present a set of structural constraints that feed reflection tomography with geometrical information. These constraints impose the desired characteristics (flatness, shape, position, etc.) on imaged reflectors but act on the velocity update. Failure to respect the constraints indicates either velocity inaccuracies or wrong assumptions concerning the constraints. Reflection tomography with structural constraints is a flexible framework that can be specialized in order to achieve different goals: among others, to flatten the base of salt bodies or detachment surfaces, to recover the horizontalness of oil–water contacts, or to impose the co‐location of the same imaged horizon between PP and PS images. The straightforward application of structural constraints is that of regularizing tomography through geological information, particularly at the latest stages of the depth imaging workflow, when the depth migration structural setting reached a consistent geological interpretation. Structural constraints are also useful in minimizing the well‐to‐seismic mis‐ties. Moreover, they can be used as a tool to check the consistency of interpreters' hypothesis with seismic data. Indeed, inversion with structural constraints will preserve image focusing only if the interpreters' insights are consistent with the data. Results from synthetic and real data demonstrate the effectiveness of reflection tomography with structural constraints.     

A neuro fuzzy approach to recognize rock magnetic and lithological patterns in a stratigraphic well from the Llanos Foreland Basin (Colombia)

The Neuro Fuzzy System (NFS) is a hybrid algorithm that combines fuzzy logic with neural networks. Since it can be used as a pattern recognition technique, we explore its potential to characterize the major lithological units encompassed by the first 512 m of the Colombian stratigraphic well Saltarin 1A (Guayabo and León Formations). Thus, we employ the NFS to infer the magnetic remanence S-ratio using bulk magnetic susceptibility (κ), κ-normalized saturation isothermal magnetization (SIRM κ) and/or volume of shale (V_sh) obtained from a gamma-ray log. The best results in terms of their corresponding Root Mean-Square Error (RMSE) values, throughout most of the upper Guayabo Formation, where magnetite seems to be an important magnetic phase, are attained with logκ and SIRM κ as input variables. Beyond 350 m downcore, the quality of the inference decreases over the León Formation, characterized by a significant presence of pyrrhotite. However, the extra input variable V_sh adjusts the inferred S-ratio to their experimental counterparts throughout this formation suggesting that the early diagenesis process that led to the formation of dispersed clay in these samples was also responsible for the formation of pyrrhotite. Hence the inclusion of manifold input data increases the ability of the net to predict S-ratio in complex geological settings with a sequence of changing lithologies, varying amounts and types of magnetic minerals, and different distributions of mineral grain sizes. In case these variables do not properly infer the actual S-ratio data, the extent of the different lithostratigraphic units would be still identifiable in some cases by the uneven quality of the correlation observed between inferred and experimental values.     

Asymptotic Analysis and Layer Decomposition for the Couplex Exercise

This document explores the possibilities of multiscale expansions and domain decomposition to solve part of the Couplex 1 exercise. We concentrate on the hydrostatic pressure and show that the slenderness of the domain and the large variations of the Darcy constants allows an analytical approximation which our test reveals to be true to relative errors smaller than 10^–3. The numerical tests are done in 2D with freefem+ and in 3D with freefem3D. Some considerations are also given for Iodine transport.     

Numerical relationships between magnetic parameters measured in Quaternary sediments and global paleoclimatic proxies

The complexity of most geological and geophysical problems prompts sometimes the use of non linear mathematical methods to handle them. An adaptive neuro fuzzy inference system (ANFIS) that combines fuzzy logic with neural networks, is applied here to study a paleoclimate section from the Quaternary sedimentary fill of the Lake Mucubají (western Venezuela). The purpose of this work is to find a set of numerical relationships that could predict the possible connections between oxygen isotope (δ¹⁸O) values from two different locations in the northern hemisphere (Ammersee in southern Germany and an ice core from the Greenland Ice Core Project — GRIP) and rock-magnetic parameters measured in Mucubají samples (i.e. mass-specific magnetic susceptibility — χ, magnetic remanence S-ratio, mass-specific saturation isothermal remanent magnetization — SIRM and anhysteretic remanent magnetization — ARM). The best inferences in terms of coefficient of determionation R² and the Root Mean-Square Error (RMSE) are obtained using those magnetic data as input that include information about magnetite grain size distributions, e.g., SIRM and ARM in FIS structures [1χ, 4ARM] and [4ARM, 1SIRM]. A comparison between Ammersee and GRIP actual data, as well as their corresponding inferences for the FIS structure [4ARM, 1SIRM], reveals a reasonable good inference of global trends for both records, overlooking the regional and/or local paleoclimate forcings that might have affected Ammersee. A better correlation between global isotope paleoclimate records and magnetic proxies, is perhaps prevented by the role played by local and regional paleoclimate and tectonism in Mucubají. We also argue that the ratio of ARM over SIRM appears to be related in a complex way to the onset and to the end of the Younger Dryas. Our novel approach to the assessment of a specific paleoclimate case study shows the potential of the ANFIS technique in solving problems where traditional univariate and multivariate linear regression methods could prove inadequate.     

The pairwise velocity dispersion of galaxies: effects of non-radial motions

I discuss the effect of non-radial motions on the small-scale peculiar pairwise velocity dispersions (PVD) of galaxies in a cold dark matter (CDM) model and calculate the PVD for the SCDM model by means of the refined cosmic virial theorem (CVT), taking account of non-radial motions by means of the Del Popolo & Gambera model. I compare the results of the present model with the data from Davis & Peebles, the IRAS value at 1  h ⁻¹ Mpc of Fisher et al. and Marzke et al. I show that while the SCDM model disagrees with the observed values, as pointed out by several authors, taking account of non-radial motions produces smaller values for the PVD. At r ≤1  h ⁻¹ Mpc the result is in agreement with Bartlett & Blanchard. In the light of this last paper, the result may be also read as a strong dependence of the CVT prediction on the model chosen to describe the mass distribution around galaxies, suggesting that the CVT cannot be taken as a direct evidence for a low-density Universe. Similarly to what is shown by Del Popolo & Gambera and Del Popolo et al., the agreement of our model to the observational data is because of a scale-dependent bias induced by the presence of non-radial motions. As the assumptions on which CVT is based have been questioned by several authors, I also calculated the PVD using the redshift distortion in the redshift-space correlation function, ξ _z( r _p, π), and compared it with the PVD measured from the Las Campanas Redshift Survey by Jing et al. The result confirms that non-radial motions influence the PVD making them agree better with the observed data.     

Granulation deformation near and in sunspot regions

High-resolution white-light pictures are analyzed to study the differences between the granular size near sunspot penumbrae and in light-bridges presenting granular structure and that of the quiet photosphere. No difference is found between the mean granular diameter in light-bridges and the quiet photosphere. The dispersion found in the results corresponding to different zones around the sunspots indicates that the size of the granulation may vary from place to place near the sunspots, its mean value not differing significantly from that of the quiet photosphere. A possible systematic bias in the selection of the granules by Macris (1979) is found.     

A focal plane detector design for a wide-band Laue-lens telescope

The energy range above 60 keV is important for the study of many open problems in high energy astrophysics such as the role of Inverse Compton with respect to synchrotron or thermal processes in GRBs, non thermal mechanisms in SNR, the study of the high energy cut-offs in AGN spectra, and the detection of nuclear and annihilation lines. Recently the development of high energy Laue lenses with broad energy bandpasses from 60 to 600keV have been proposed for a Hard X ray focusing Telescope (HAXTEL) in order to study the X-ray continuum of celestial sources. The required focal plane detector should have high detection efficiency over the entire operative range, a spatial resolution of about 1mm, an energy resolution of a few keV at 500keV and a sensitivity to linear polarization. We describe a possible configuration of the focal plane detector based on several CdTe/CZT pixelated layers stacked together to achieve the required detection efficiency at high energy. Each layer can operate both as a separate position sensitive detector and polarimeter or work with other layers to increase the overall photopeak efficiency. Each layer has a hexagonal shape in order to minimize the detector surface required to cover the lens field of view. The pixels would have the same geometry so as to provide the best coupling with the lens point spread function and to increase the symmetry for polarimetric studies.