Upward continuation of potential fields from a polyhedral surface1

An equivalent source procedure is derived for upward continuation of unevenly spaced gravity and magnetic data. The dipole layer is placed on a topographic relief approximated by a polyhedral surface, the stations being the vertices of the triangular faces. The dipoles have linear magnitudes, being directed along the normal vector over each triangle. The unknown values of the dipole magnitudes at each station are obtained by a suitable modification of the usual integral equation considering the discontinuity of the normal vector at each vertex of the dipole surface. Profile data processing is also studied. A numerical test outlines the accuracy and the limitations of the model for the case of a magnetic field significantly perturbed by a rough topographic relief.     

Line integrals of potential field data1

Information on the mass and the spatial location of an arbitrary source body can be obtained by performing suitable integrations of 3D gravity and magnetic data along an infinite straight line. No assumptions on the density/magnetization distribution or the shape and location of the source are required. For an oblique borehole, a relationship between the lower limit of the source mass and the distance to the body is obtained. The mass contrast and the magnetic moment of the source can also be estimated. For a vertical borehole, both gravity and vertical magnetic component anomalies have equal areas to the left and right of the depth axis. The particular case of a horizontal gallery not intersecting the body is also studied. If the source is intersected, a lower limit is estimated for the maximum thickness of the body along the gallery. Information on the vertical coordinate of the centre of mass of the source can also be obtained. Numerical tests with synthetic gravity data support the theoretical results.     

The photoelectron spectra of some Tl-Sb sulphosalts

The results of X-ray induced photoelectron spectroscopy (XPS) experiments on several phases of the ternary system Tl-Sb-S are reported. The binding energies of the inner S, Sb and Tl electrons increase with increasing quantities of Sb and decreasing amounts of Tl in these compounds. This is explained by the influences of the proportions of the bonded metals on the effective electron affinity of S. The higher proportions of the more electronegative element bonded to S cause the increase of its effective electron affinity. The results for Tl₂S (carlinite), Tl₃SbS₃, TlSbS₂ (weissbergite), TlAsS₂ (lorandite) and Sb₂S₃ (antimonite) can be interpreted in this way. The results for Tl₄S₃ suggest a predominantly covalent character of bonding for both Tl(III) and Tl(I), which are present in this sulfide. From comparison with Tl₃SbS₄ it could be supposed that Tl(III)-S bond has a more covalent character than Sb(V)-S bond. The results for Tl₃SbS₄ are in agreement with crystal structure data and the results of Moessbauer spectroscopy. For AsS (realgar) the binding energies of the inner electrons of As and S significantly increase, showing that the electrons in molecular orbitals are less strongly bonded to individual atoms, as compared to pure elements. The results for the amorphous TlSb₅S₈ (corresponding in composition to parapierrotite) suggest that in amorphous compound the Tl-S bonding is stronger and the coordination of Tl more regular than in a crystalline one.     

The alleged kimberlite-carbonatite relationship: evidence from ilmenite and spinel from Premier and Wesselton mines and the Benfontein sill,South Africa

Carbonate-rich, SiO₂-poor residua are developed in some kimberlites solidifying as ocelli, layers, or discrete dikes which satisfy petrographic definitions of carbonatite. Arguments that these rocks have mineralogies, antecedents, and comagmatic rocks differing from those of the carbonatites in alkaline rock complexes, including the specific observation that kimberlites and carbonatites contain ilmenites and spinels of different composition, have been used to refute the alleged kimberlite-carbonatite relationship. New microprobe analyses of ilmenites and spinels from carbonate-rich rocks associated with kimberlites in three South African localities correspond to spinels and ilmenites of carbonatites from alkalic complexes, or have characteristics intermediate between those of carbonatites and kimberlites. The ilmenites are distinguished from kimberlite ilmenites by higher MnO, FeTiO₃, and Nb₂O₅, and by negligible Cr₂O₃. The spinels are distinguished from kimberlite spinels by their Al₂O₃ and Cr₂O₃ contents. There is clearly a genetic relationship between the kimberlites and the carbonate-rich rocks, despite the observation that their ilmenites and spinels are distinctly different, which indicates that the same observation is not a valid argument against a petrogenetic relationship between kimberlites and carbonatites. These rocks are among the diverse products from mantle processes influenced by CO₂, and we believe that the petrogenetic links among them are forged in the upper mantle. We see insufficient justification to deny the name carbonatite to carbonate-rich rocks associated with kimberlites if they satisfy the petrographic definition in terms of major mineralogy.     

An experiment with joint mapping of the vertical velocity and its horizontal gradient in the study of recent tectonic movements

An experiment is described dealing with the detailed mapping of the vertical velocity V of recent crustal movements in the central Danubian riverside region of the Moesian platform in Northern Bulgaria. The error in V being significant here, the problem of statistical control on the reliability of the initial data and results of such studies is investigated; this question is also relevant for other platform regions in Eastern and Central Europe. With the aim of assessing the significance of a movement, the sampling distribution of V is used. The problem is reduced to a checking of a standard statistical hypothesis and to a probabilistic characterization of V. In order to estimate the statistical significance of the relative velocity δV between neighbouring benchmarks, the distribution of the module G of the horizontal gradient of V is studied. An agreement between the sampling distribution obtained and the theoretical Weibull distribution is established. A sufficiently high probability P>86% is determined for the significance of δV in this region and an optimal interval of the isobases is suggested. The experiment itself embraces three variant solutions: independent mapping of V, independent mapping of G, and mapping of V with respect to G. The analysis and comparison of the results obtained demonstrate the high efficiency of the third solution.     

Dynamics of lava flow fronts, Arenal Volcano, Costa Rica

Arenal Volcano has effused basaltic andesite lava flows nearly continuously since September, 1968. The two different kinds of material in flows, lava and lava debris, have different rheologic properties and dynamic behavior. Flow morphology depends on the relationship between the amount and distribution of the lava and the debris, and to a lesser extent the ground morphology.Two main units characterize the flows: the channel zone and the frontal zone. The channel zone consists of two different units, the levées and the channel proper. A velocity profile in the channel shows a maximum value at the plug where the rate of shear is zero, and a velocity gradient increasing outward until, at the levées, the velocity becomes zero. Cooling produces a marked temperature gradient in the flow, leading to the formation of debris by brittle fracture when a critical value of shear rate to viscosity is reached. When the lava supply ceases, much of this debris and part of the lava is left behind after the flow nucleus drains out, forming a collapsed channel.Processes at the frontal zone include levée formation, debris formation, the change in shape of the front, and the choice of the flow path. These processes are controlled primarily by the rheological properties of the lava.Frontal zone dynamics can be understood by fixing the flow front as the point of reference. The lava flows through the channel into the front where it flows out into the levées, thereby increasing the length of the channel and permitting the front to advance. The front shows a relationship of critical height to the yield strength (τ₀) surface tension, and slope; its continued movement is activated by the pressure of the advancing lava in the channel behind. For an ideal flow (isothermal, homogeneous, and isotropic) the ratio of the section of channel proper to the section of levées is calculated and the distance the front will have moved at any time t_x can be determined once the amount of lava available to the front is known. Assuming that the velocity function of the front {G(t)} during the collapsing stage is proportional to the entrance pressure of the lava at the channel-front boundary, an exponential decrease of velocity through time is predicted, which shows good agreement with actual frontal velocity measurements taken on two flows. Local variations in slope have a secondary effect on frontal velocities.Under conditions of constant volume the frontal zone can be considered as a machine that consumes energy brought in by the lava to perform work (front advancement). While the front will use its potential energy to run the process, the velocity at which it occurs is controlled by the activation energy that enters the system as the kinetic energy of the lava flowing into the front. A relation for the energy contribution due to frontal acceleration is also derived. Finally the entrance pressure, that permits the front to deform, is calculated. Its small value confirms that the lava behaves very much like a Bingham plastic.     

Instrumental Neutron Activation Analysis of NBS 1633a Fly Ash and 1632a Bituminous Coal Reference Samples With the Use of Short Irradiation

Instrumental neutron activation analysis has been applied to the analysis of recently issued NBS reference materials SRM 1633a Fly Ash and SRM 1632a Bituminous Coal. Only the nuclides with half-lives of less than 48 hours have been used. The results for the determination of 18 elements (Fly Ash) and 19 elements (Coal) are compared with the values given by NBS and with recent results obtained by other investigators.     

Numerical Modeling of the Ston-Slano (Croatia) Aftershock Sequence

The rate of aftershock occurrence after the M6 Ston-Slano (Croatia) earthquake is modeled as the Epidemic Type Aftershock Sequence (ETAS). Increase of the modeled cumulative number of aftershocks with time was fitted to observations by the least-squares criterion using the combined grid-search and Monte-Carlo approach. This enabled not only the estimation of the most probable ETAS parameters, but also the determination of their confidence limits, as well as the estimation of the bias between them. It has been found that the bias is significant for some of the parameter pairs, regardless of the threshold magnitude assumed. Residual analyses revealed that all strong aftershocks (M _L 4.5) occurred during the periods of normal to high aftershock activity. There were two periods of quiescence in the sequence, both of which were followed by a strong aftershock.     

Analysis of sustained long-period activity at Etna Volcano, Italy

Following the installation of a broadband network on Mt. Etna, sustained Long-Period (LP) activity was recorded accompanying a period of total quiescence and the subsequent onset of the 2004–2005 effusive episode. From about 56000 events detected by an automatic classification procedure, we analyse a subset of about 3000 signals spanning the December 17th, 2003–September 25th, 2004, time interval. LP spectra are characterised by several, unevenly-spaced narrow peaks spanning the 0.5–10 Hz frequency band. These peaks are common to all the recording sites of the network, and different from those associated with tremor signals. Throughout the analysed time interval, LP spectra and waveforms maintain significant similarity, thus indicating the involvement of a non-destructive source process that we interpret in terms of the resonance of a fluid-filled buried cavity. Polarisation analysis indicates radiation from a non-isotropic source involving large amounts of shear. Concurrently with LP signals, recordings from the summit station also depict Very-Long-Period (VLP) pulses whose rectilinear motion points to a region located beneath the summit craters at depths ranging between 800 and 1100 m beneath the surface. Based on a refined repicking of similar waveforms, we obtain robust locations for a selected subset of the most energetic LP events from probabilistic inversion of travel-times calculated for a 3D heterogenous structure. LP sources cluster in a narrow volume located beneath the summit craters, and extending to a maximum depth of ≈ 800 m beneath the surface. No causal relationships are observed between LP, VLP and tremor activities and the onset of the 2004–2005 lava effusions, thus indicating that magmatic overpressure played a limited role in triggering this eruption. These data represent the very first observation of LP and VLP activity at Etna during non-eruptive periods, and open the way to the quantitative modelling of the geometry and dynamics of the shallow plumbing system.     

Petrographic and geochemical evidence for the formation of primary, bacterially induced lacustrine dolomite: La Roda 'white earth' (Pliocene, central Spain)

Upper Pliocene dolomites (‘white earth’) from La Roda, Spain, offer a good opportunity to evaluate the process of dolomite formation in lakes. The relatively young nature of the deposits could allow a link between dolomites precipitated in modern lake systems and those present in older lacustrine formations. The La Roda Mg‐carbonates (dolomite unit) occur as a 3·5‐ to 4‐m‐thick package of poorly indurated, white, massive dolomite beds with interbedded thin deposits of porous carbonate displaying root and desiccation traces as well as local lenticular gypsum moulds. The massive dolomite beds consist mainly of loosely packed 1‐ to 2‐μm‐sized aggregates of dolomite crystals exhibiting poorly developed faces, which usually results in a subrounded morphology of the crystals. Minute rhombs of dolomite are sparse within the aggregates. Both knobbly textures and clumps of spherical bodies covering the crystal surfaces indicate that bacteria were involved in the formation of the dolomites. In addition, aggregates of euhedral dolomite crystals are usually present in some more clayey (sepiolite) interbeds. The thin porous carbonate (mostly dolomite) beds exhibit both euhedral and subrounded, bacterially induced dolomite crystals. The carbonate is mainly Ca‐dolomite (51–54 mol% CaCO₃), showing a low degree of ordering (degree of ordering ranges from 0·27 to 0·48). Calcite is present as a subordinate mineral in some samples. Sr, Mn and Fe contents show very low correlation coefficients with Mg/Ca ratios, whereas SiO₂ and K contents are highly correlated. δ¹⁸O‐ and δ¹³C‐values in dolomites range from ?3·07‰ to 5·40‰ PDB (mean=0·06, σ=1·75) and from ?6·34‰ to ?0·39‰ PDB (mean=?3·55, σ=1·33) respectively. Samples containing significant amounts of both dolomite and calcite do not in general show significant enrichment or depletion in ¹⁸O and ¹³C between the two minerals. The correlation coefficient between δ¹⁸O and δ¹³C for dolomite is extremely low and negative (r=?0·05), whereas it is higher and positive (r=0·47) for calcite. The lacustrine dolomite deposit from La Roda is interpreted mainly as a result of primary precipitation of dolomite in a shallow, hydrologically closed perennial lake. The lake was supplied by highly saturated HCO₃^?/CO₃^2? groundwater that leached dolomitic Mesozoic formations. Precipitation of dolomite from alkaline lake waters took place under a semi‐arid to arid climate. However, according to our isotopic data, strong evaporative conditions were not required for the formation of the La Roda dolomite. A significant contribution by bacteria to the formation of the dolomites is assumed in view of both petrographic and geochemical evidence.