High-Energy 3He-Rich Solar Particle Events

Energetic particle observations of the ERNE instrument on board SOHO enable measurements of ³He and ⁴He fluxes beyond 15 MeV nucleon^–1 with good statistical resolution. We report results of a survey of the ERNE observations covering the period from 8 February 1999 to 6 December 2000. We find 10 and 5 days during which the ³He-to-⁴He ratio exceeds the levels of 20% and 50%, respectively. Those periods include, in particular, four ³He-rich events that are sufficiently strong for a reasonably accurate estimate of the time-intensity profiles. We analyze the history of solar and interplanetary phenomena associated with these high-energy ³He-rich events. Basic properties of such events and significant solar and interplanetary factors are formulated. The significant factors comprise, in particular, a strong, impulsive flare, typically observed about day before the ³He onset, and an interplanetary shock wave or magnetic field enhancement arriving at 1 AU about frac43 day after the ³He onset. The high-energy ³He-rich events make up a new kind of hybrid events, possessing the impulsive-type composition and the gradual-type time-profiles. We emphasize a dependence of the resultant particle event on the history of the particular solar eruption comprising coronal mass ejection (CME) and the flare associated with the CME.     

Plane wave propagation in nonlinear-elastic anisotropic media

Summary. Kelvin-Christoffel equations describing plane wave propagation in anisotropic media are generalized to account for the effects of nonlinear elasticity. The polarization and waveform of nonlinear distortions of a transient plane wave are investigated by means of perturbation theory. Detailed analysis for an anisotropic medium with hexagonal symmetry shows that for "pure" shear-waves the polarization vector of the nonlinear component is always perpendicular to that of the linear wave. In the case of a high-amplitude excitation (for instance, in the vicinity of large earthquakes) the influence of nonlinearity may cause distortions of shear-wave polarization, which contains the most reliable information on the presence and characteristics of anisotropy. The solutions presented in this paper make it possible to solve reflection-transmission problems in nonlinear-elastic anisotropic media.     

Partitioning of lanthanides and Y between immiscible silicate and fluoride melts, fluorite and cryolite and the origin of the lanthanide tetrad effect in igneous rocks

Some F-rich granitic rocks show anomalous, nonchondritic ratios of Y/Ho, extreme negative Eu anomalies, and unusual, discontinuous, segmented chondrite-normalised plots of rare earth elements (REE). The effects of F-rich fluids have been proposed as one of the explanations for the geochemical anomalies in the evolved granitic systems, as the stability of nonsilicate complexes of individual rare earths may affect the fluid-melt element partitioning. The lanthanide tetrad effect, related to different configurations of 4f-electron subshells of the lanthanide elements, is one of the factors affecting such complexing behaviour. We present the first experimental demonstration of the decoupling of Y and Ho, and the tetrad effect in the partitioning of rare earths between immiscible silicate and fluoride melts. Two types of experiments were performed: dry runs at atmospheric pressure in a high-temperature centrifuge at 1100 to 1200°C, and experiments with the addition of H₂O at 700 to 800°C and 100 MPa in rapid-quench cold-seal pressure vessels. Run products were analysed by electron microprobe (major components), solution-based inductively coupled plasma mass spectrometry (ICP-MS) (REE in the centrifuged runs), and laser ablation ICP-MS (REE and Li in the products of rapid-quench runs). All the dry centrifuge runs were performed at super-liquidus, two-phase conditions. In the experiments with water-bearing mixtures, minor amounts of aqueous vapour were present in addition to the melts. We found that lanthanides and Y concentrated strongly in the fluoride liquids, with two-melt partition coefficients reaching values as high as 100-220 in water-bearing compositions. In all the experimental samples, two-melt partition coefficients of lanthanides show subtle periodicity consistent with the tetrad effect, and the partition coefficient of Y is greater than that of Ho. One of the mixtures also produced abundant fluorite (CaF₂) and cryolite (Na₃AlF₆) crystals, which enabled us to study fluorite-melt and cryolite-melt REE partitioning. REE concentrations in fluorite are high and comparable to those in the fluoride melt. However, fluorite-melt partition coefficients appear to depend mostly on ionic radii and show neither significant tetrad anomalies, nor differences in Y and Ho partitioning. In contrast, REE concentrations in cryolite are low (∼5-10 times lower than in the silicate melt), and cryolite-melt REE partitioning shows very strong tetrad and Y-Ho anomalies. Our results imply that Y-Ho and lanthanide tetrad anomalies are likely to be caused mainly by aluminofluoride complexes, and the tetrad REE patterns in natural igneous rocks can result from fractionation of F-rich magmatic fluids.     

Element partitioning between immiscible borosilicate liquids: A high-temperature centrifuge study

The main objectives of this study were to investigate conditions for stable and metastable liquid immiscibility in dry borosilicate synthetic systems and to evaluate effects of temperature and bulk melt composition on two-liquid element partitioning and boron speciation. To distinguish between the stable immiscibility and quench heterogeneity, we used high-temperature centrifuge phase separation. For the case of stable liquid immiscibility, silica-rich (LS) and borate-rich (LB) conjugate liquids formed two distinct layers separated by a sharp meniscus. The liquids were quenched into glasses, which were analysed by electron microprobe. Some of the glasses were also studied by Raman spectroscopy. We used several synthetic mixtures along the danburite-anorthite (CaB₂Si₂O₈-CaAl₂Si₂O₈) and danburite-reedmergnerite (CaB₂Si₂O₈-NaBSi₃O₈) joins. In addition, we studied four complex, six-component, Mg-bearing compositions with variable Na₂O and Al₂O₃ contents. The experiments show that the width of the LS-LB miscibility gap decreases more rapidly with the B-Al substitution (in the danburite-anorthite join) than with the Ca-Na substitution, implying that interactions between network-forming elements have a greater effect on borate-silicate unmixing than the nature of network-modifying cations. Ca and Mg partition strongly to the depolymerised borate-rich liquid with LB-LS partition coefficients of ∼40 and higher. On the other hand, two-liquid partition coefficients of Na and Al in most cases are close to 1 and show complex variations with temperature and bulk melt composition. Raman spectra of LB glasses quenched at different temperatures suggest that the proportion of trigonal boron in bulk boron content decreases with decreasing temperature. The change in boron speciation appears to affect Al and Na two-liquid partitioning in such a way that at low temperatures, the latter element becomes more compatible with LS.     

Nonconforming Finite Volume Methods

We extend and generalize recently proposed finite volume methods using the framework of mixed finite element methods. Proposed discretizations are defined for tensor permeabilities and naturally produce a generalization of harmonic averaging, and are therefore well suited for heterogeneous and anisotropic media. They are locally mass conservative and work on extremely flexible distorted meshes. Flux variables can be excluded locally and the resulting discretizations for the pressures has the same stencils for Voronoi/PEBI grids as 2-point finite volume discretizations currently used in many simulators.     

Moveout analysis for tranversely isotropic media with a tilted symmetry axis

The transversely isotropic (TI) model with a tilted axis of symmetry may be typical, for instance, for sediments near the flanks of salt domes. This work is devoted to an analysis of reflection moveout from horizontal and dipping reflectors in the symmetry plane of TI media that contains the symmetry axis. While for vertical and horizontal transverse isotropy zero-offset reflections exist for the full range of dips up to 90°, this is no longer the case for intermediate axis orientations. For typical homogeneous models with a symmetry axis tilted towards the reflector, wavefront distortions make it impossible to generate specular zero-offset reflected rays from steep interfaces. The ‘missing’ dipping planes can be imaged only in vertically inhomogeneous media by using turning waves. These unusual phenomena may have serious implications in salt imaging. In non-elliptical TI media, the tilt of the symmetry axis may have a drastic influence on normal-moveout (NMO) velocity from horizontal reflectors, as well as on the dependence of NMO velocity on the ray parameter p (the ‘dip-moveout (DMO) signature’). The DMO signature retains the same character as for vertical transverse isotropy only for near-vertical and near-horizontal orientation of the symmetry axis. The behaviour of NMO velocity rapidly changes if the symmetry axis is tilted away from the vertical, with a tilt of ±20° being almost sufficient to eliminate the influence of the anisotropy on the DMO signature. For larger tilt angles and typical positive values of the difference between the anisotropic parameters ε and δ, the NMO velocity increases with p more slowly than in homogeneous isotropic media; a dependence usually caused by a vertical velocity gradient. Dip-moveout processing for a wide range of tilt angles requires application of anisotropic DMO algorithms. The strong influence of the tilt angle on P-wave moveout can be used to constrain the tilt using P-wave NMO velocity in the plane that includes the symmetry axis. However, if the azimuth of the axis is unknown, the inversion for the axis orientation cannot be performed without a 3D analysis of reflection traveltimes on lines with different azimuthal directions.     

Liquid Immiscibility and the Evolution of Basaltic Magma

This experimental study examines relationships between alternativeevolution paths of basaltic liquids (the so-called Bowen andFenner trends), and silicate liquid immiscibility. Syntheticanalogues of natural immiscible systems exhibited in volcanicglasses and melt inclusions were used as starting mixtures.Conventional quench experiments in 1 atm gas mixing furnacesproved unable to reproduce unmixing of ferrobasaltic melts,yielding instead either turbid, opalescent glasses, or crystallizationof tridymite and pyroxenes. In contrast, experiments involvingin situ high-temperature centrifugation at 1000g (g = 9·8m/s²) did yield macroscopic unmixing and phase separation. Centrifugationfor 3–4 h was insufficient to complete phase segregation,and resulted in sub-micron immiscible emulsions in quenchedglasses. For a model liquid composition of the Middle Zone ofthe Skaergaard intrusion at super-liquidus temperatures of 1110–1120°C,centrifugation produced a thin, silicic layer (64·5 wt%SiO₂ and 7·4 wt% FeO) at the top of the main Fe-richglass (46 wt% SiO₂ and 21 wt% FeO). The divergent compositionsat the top and bottom were shown in a series of static runsto crystallize very similar crystal assemblages of plagioclase,pyroxene, olivine, and Fe–Ti oxides. We infer from theseresults that unmixing of complex aluminosilicate liquids maybe seriously kinetically hampered (presumably by a nucleationbarrier), and thus conventional static experiments may not correctlyreproduce it. In the light of our centrifuge experiments, immiscibilityin the Skaergaard intrusion could have started already at thetransition from the Lower to the Middle Zone. Thus, magma unmixingmight be an important factor in the development of the Fe-enrichmenttrend documented in the cumulates of the Skaergaard LayeredSeries. KEY WORDS: liquid immiscibility; Skaergaard; layered intrusions; experimental petrology     

Correction to: Multilevel Graph Partitioning for Three-Dimensional Discrete Fracture Network Flow Simulations

Mathematical Geosciences - The publication of this article unfortunately contained a mistake. The assignment to the affiliations of author Satish Karra was not correct