Smithsonian institution's global volcanism network

Information included in this summary is based on more detailed reports published in the Bulletin of the Global Volcanism Network, vol. 30, no. 2, February 2005 (on the Internet at ). Edited by scientists at the Smithsonian, this Bulletin includes reports provided by a worldwide network of correspondents. The reports contain the names and contact information for all sources. Please note that these reports are preliminary and subject to change as events are studied in more detail. The Global Volcanism Program welcomes further reports of current volcanism, seismic unrest, monitoring data, and field observations.     

Correlation of foreshocks and aftershocks and asperities

A close correlation in spatial distribution of local seismic activity and energy release patterns before and after the 1979 Petatlan, Mexico earthquake suggests heterogeneity within the fault plane of this major low-angle thrust event associated with subduction along the Middle America Trench. A simple two-asperity model is proposed to account for the complexity. Foreshocks and aftershocks of the neighboring 1981 Playa Azul earthquake showed a similar pattern. As both events occurred at the junction of the Orozco Fracture Zone and the Middle America Trench, we speculate that the observed complex fault plane is caused by subduction of the rugged ocean floor of the Orozco Fracture Zone. Short-term precursory seismicity prior to the Petatlan earthquake can be explained by using the asperity model and migration of a slip front from the south-east to the north-west across the main shock source region.     

Testing the modified Press–Schechter model against N-body simulations

A modified version of the extended Press–Schechter model for the growth of dark-matter haloes was introduced in two previous papers, with the aim of explaining the mass–density relation shown by haloes in high-resolution cosmological simulations. In this model, major mergers are well separated from accretion, thereby allowing a natural definition of halo formation and destruction. This makes it possible to derive analytic expressions for halo formation and destruction rates, the mass accretion rate and the probability distribution functions of halo formation times and progenitor masses. The stochastic merger histories of haloes can be readily derived and easily incorporated into semi-analytical models of galaxy formation, thus avoiding the usual problems encountered in the construction of Monte Carlo merger trees from the original extended Press–Schechter formalism. Here we show that the predictions of the modified Press–Schechter model are in good agreement with the results of N -body simulations for several scale-free cosmologies.     

Chromospheric Height and Density Measurements in a Solar Flare Observed with RHESSI – II. Data Analysis

We present an analysis of hard X-ray imaging observations from one of the first solar flares observed with the Reuven Ramaty High-Energy Solar Spectroscopic Imager (RHESSI) spacecraft, launched on 5 February 2002. The data were obtained from the 22 February 2002, 11:06 UT flare, which occurred close to the northwest limb. Thanks to the high energy resolution of the germanium-cooled hard X-ray detectors on RHESSI we can measure the flare source positions with a high accuracy as a function of energy. Using a forward-fitting algorithm for image reconstruction, we find a systematic decrease in the altitudes of the source centroids z(ε) as a function of increasing hard X-ray energy ε, as expected in the thick-target bremsstrahlung model of Brown. The altitude of hard X-ray emission as a function of photon energy ε can be characterized by a power-law function in the ε=15–50 keV energy range, viz., z(ε)≈2.3(ε/20 keV)^−1.3 Mm. Based on a purely collisional 1-D thick-target model, this height dependence can be inverted into a chromospheric density model n(z), as derived in Paper I, which follows the power-law function n _e(z)=1.25×10¹³(z/1 Mm)^−2.5 cm⁻³. This density is comparable with models based on optical/UV spectrometry in the chromospheric height range of h≲1000 km, suggesting that the collisional thick-target model is a reasonable first approximation to hard X-ray footpoint sources. At h≈1000–2500 km, the hard X-ray based density model, however, is more consistent with the `spicular extended-chromosphere model' inferred from radio sub-mm observations, than with standard models based on hydrostatic equilibrium. At coronal heights, h≈2.5–12.4 Mm, the average flare loop density inferred from RHESSI is comparable with values from hydrodynamic simulations of flare chromospheric evaporation, soft X-ray, and radio-based measurements, but below the upper limits set by filling-factor insensitive iron line pairs.     

Regularized Energy-Dependent Solar Flare Hard X-Ray Spectral Index

The deduction from solar flare X-ray photon spectroscopic data of the energy-dependent model-independent spectral index is considered as an inverse problem. Using the well-developed regularization approach we analyze the energy dependency of spectral index for a high-resolution energy spectrum provided by Ramaty High Energy Solar Spectroscopic Imager (RHESSI). The regularization technique produces much smoother derivatives while avoiding additional errors typical of finite differences. It is shown that observations imply a spectral index varying significantly with energy, in a way that also varies with time as the flare progresses. The implications of these findings are discussed in the solar flare context.     

Die deutschen Bodentypen nach dem heutigen Stande der Bodentypenlehre

Ohne Zusammenfassung     

Numerical Model of Two-Phase Flow in Dissolvable Porous Media and Simulation of Reservoir Acidizing

Mathematical model of porous media dissolution coupled with two-phase flow is proposed. The model is based on the conception of dissolvable porous medium with deformable mass-variable porous skeleton. Model can be used for simulation of coupled chemo- and hydrogeomechanical processes which are difficult to examine experimentally. Acidizing of calcite oil reservoir is used as an example of the process. Water solution of hydrochloric acid and oil are two fluid phases of the model with several components. Dissolvable porous media is treated as deformable mass-variable solid phase. Change in mass of the solid phase is caused by hydrochloric acid dissolving the calcite part of the solid phase. Dissolution is supposed to be congruent; kinetics is governed by the Nernst law. Software for numerical solution of the model is developed. It uses AmgCL parallel library for high-performance computing in order to deal with large algebraic systems on the each time step of calculations. The library uses algebraic multigrid methods for preconditioning and parallel iterative solvers. NVidia CUDA framework is used as a backend to perform GPGPU calculations, because it proved to be faster than OpenCL framework on this problem. Numerical experiments on the basis of data set from real reservoirs are conducted with the developed software. Good correlation between field and calculated data is achieved. Numerical experiments for different configurations of heterogeneous layer are performed. Acidizing of layers with highly permeable conduit and with random distribution of permeability is modeled.     

Effect of hydrocarbon-contaminated fluctuating groundwater on magnetic properties of shallow sediments

We investigate magnetic phase (trans)formation in the presence of petroleum hydrocarbons and its relation to bacterial activity, in particular in the zone of remediation driven fluctuating water levels at a former military air base in the Czech Republic. In a previous study an increase of magnetite concentration from the groundwater table towards the top of the groundwater fluctuation zone (GWFZ) was reported, however with limited reliability as there was no control on small-scale effects. To recognize statistically significant magnetic signatures versus depth, we obtained multiple sediment cores from three locations in January 2011 and April 2012, penetrating the unsaturated zone, the GWFZ and the uppermost one meter below the groundwater level (～2.3 m depth at the time of sampling). Magnetic concentration variation versus depth was determined by measuring magnetic susceptibility (MS) and remanence parameters. Small-scale features were identified and eliminated by statistical processing of multiple cores. A trend of increasing MS values from the lowermost position of the groundwater table upward was verified and highest magnetic concentration was found at the top of the GWFZ. Magnetic mineralogy indicates that newly formed fine-grained magnetite in the single domain to small pseudo-single domain range is responsible for the MS enhancement confirming previous results. There is no correlation with the depth variation of hydrocarbon (HC) concentrations; however, total organic carbon is linked to MS and may represent a degradation product of HC that is bioavailable for microorganisms. Bacterial activity is likely responsible for magnetite formation as indicated by most probable number (MPN) results of iron-metabolizing bacteria. The comparison of our results with an earlier study conducted at the same site revealed that magnetic concentration clearly decreased since remediation was terminated in 2008, possibly due to dissolution of magnetite.     

Magnetic particles in atmospheric particulate matter collected at sites with different level of air pollution

Magnetic measurements of deposited atmospehric dust can serve as an additional parameter in assessing environmental pollution. This method is based on the assumption that atmospherically deposited particles contain significant portion of ferrimagnetic iron oxides of anthropogenic origin, which can be easily detected. Aim of this paper is to identify clearly magnetic fraction of daily samples of particulate matter less than 10 μm (PM10), routinely used for air quality assessment and monitoring. We used combination of thermomagnetic analyses and other physical and chemical methods, including scanning electron microscopy (SEM) and Mössbauer spectroscopy. Our results show that daily samples of PM10, collected at sites with different degree of atmospheric pollution, contain magnetite of spherical shape, which is presumably of industrial origin. Thus, magnetic methods can be applied directly to the same substances, which are used routinely in air quality assessment and monitoring.     

Petrophysical and geochemical constraints on alteration processes in granites

The hydrothermal alteration of granites has large influence on their petrophysical properties. To reveal the impact of alteration on magnetic and porosity properties of granites we have conducted a complex study of effects of two largely independent alteration processes, related to chemically different fluids, in granites of the Vysoký Kámen stock (the Krudum granite body, Czech Republic). It includes the whole-rock geochemical, magnetic and pore-space characterization. The alkali feldspathization resulted in decomposition of Li-mica, quartz removal, depletion in mafic cations and growth of new alkali feldspars (albite, K-feldspar), decreasing the overall magnetic susceptibility and disrupting the pore space by its discontinuation. The preservation of the orientation of the principal susceptibility axes is likely related to insignificant influence of the feldspathization process on the paramagnetic and diamagnetic phases orientation acquired during the magma emplacement. The greisenization, on the other hand had considerably more significant effects on microstructure and physical properties of the granite. The microstructure was modified by the growth of large amounts of new phases (Li-mica, quartz and topaz). This changed the mineral density of the rock, the porosity, size and character of pores to larger, flatter and probably more connected. This led also to the complete reworking of the original anisotropy of magnetic susceptibility during the greisenization.