A Model of Magmatic Crystallization

A computer model simulating fractional crystallization at oneatmosphere pressure incorporates nine broadly-defined minerals—magnetite,olivine, hypersthene, augite, quartz, plagioclase, orthoclase,leucite, and nepheline. The crystallization temperature of eachmineral is considered to be a smooth function of the compositionof the magmatic liquid. These mineral temperature equationsare obtained by multiple linear regression analysis of informationfrom published silicate systems and rock melting experiments.The nine equations are solved for any primary liquid, withinthe broad range of common magma types, to select the crystallizingmineral or minerals. Partition ratios from published experimentsand analyses of lavas and phenocrysts permit calculation ofthe composition of the crystallizing mineral assemblage. Subtractionof a small amount of that composition from the primary liquidyields a new liquid, which may be recycled to yield a sequenceof liquids during fractional crystallization. The crystallizationmodel handles assemblages of co-precipitating minerals, andcan trace progressive saturation in new minerals, substitutionof a new mineral for an old mineral, and cessation of crystallizationof a mineral. The sequences of minerals and liquids derivedfrom a broad set of primary liquids are geologically realistic,so the model is useful in predicting phenocrysts in volcanicrocks and events during crystallization of shallow intrusions.     

Use of Rhodamine water tracer in the marshland upwelling system

Rhodamine water tracer (RWT) was used to characterize the migration of waste water within the saline subsurface of a marshland upwelling system (MUS), which is an alternative on-site waste water treatment system for coastal areas. Field tracer studies were performed to investigate the fresh waste water plume movement within the saline ground water. Pore velocities were calculated using first detection times and ranged from 0.68 to 10.7 x 10(-4) cm/sec for the loamy sandy soil matrix present at the site. Use of RWT in the field also allowed determination of main and preferential flowpaths. One- and two-dimensional laboratory experiments were performed using silica sand to investigate the interactions of the organically rich waste water with RWT within the zone surrounding the point of injection (one-dimensional) and the impact of background salinity on plume movement (two-dimensional). The results from these studies were used to help explain the field data. One-dimensional breakthrough curves revealed retardation factors for the RWT in the waste water mixture of 1.73 to 1.90. These results were similar to other researchers, indicating little interaction between the waste water and RWT. Variations in pore water salinity (5, 15, 25, and 35 ppt) were found to have a significant effect on pore water velocities of the fresh water plume (two-dimensional), indicating the need to incorporate background salinities into the design process for MUS.     

Precession and nutation in the η Carinae binary system: evidence from the X-ray light curve

It is believed that η Carinae is actually a massive binary system, with the wind–wind interaction responsible for the strong X-ray emission. Although the overall shape of the X-ray light curve can be explained by the high eccentricity of the binary orbit, other features like the asymmetry near periastron passage and the short quasi-periodic oscillations seen at those epochs have not yet been accounted for. In this paper we explain these features assuming that the rotation axis of η Carinae is not perpendicular to the orbital plane of the binary system. As a consequence, the companion star will face η Carinae on the orbital plane at different latitudes for different orbital phases and, since both the mass-loss rate and the wind velocity are latitude dependent, they would produce the observed asymmetries in the X-ray flux. We were able to reproduce the main features of the X-ray light curve assuming that the rotation axis of η Carinae forms an angle of  29°± 4°  with the axis of the binary orbit. We also explained the short quasi-periodic oscillations by assuming nutation of the rotation axis, with an amplitude of about  5°  and a period of about 22 days. The nutation parameters, as well as the precession of the apsis, with a period of about 274 years, are consistent with what is expected from the torques induced by the companion star.     

Metadolerite age relationships in the Dalradian of northwest Donegal,Ireland and their orogenic significance

The structural relationships of several intrusive sheets belonging to the extensive metadolerite suite of the Appin (lower) Dalradian rocks of the Creeslough area, northwest Donegal, Ireland have been examined. It is concluded that, relative to the regional chronology of deformation of the enclosing metasediments, sheet emplacement has occurred: pre-D2; syn-D2; between D2 and D3; and syn-D3. The height of activity came in the post-D2, pre-D3 interval and coincided with the peak of regional metamorphism in the area. These relationships in Donegal are compared with those of the pretectonic basic rocks of the southwest Highlands and models are discussed.     

Radio spectra of some markarian galaxies

Radio spectra of 16 Markarian galaxies are studied between 2.7 and 15.5 GHz. The Seyfert-type galaxies have power-law spectra with an average spectral index of -0.82, except for MRK 348, which is an active galaxy with a complex spectrum. The galaxies with featureless optical spectra have rather flat radio spectra, similar to BL Lac-type objects.     

Discrepancy of masses in double galaxies

The estimations of integral masses of double galaxies are obtained taking as a basis recent observational data on 150 pairs of galaxies with measured radial velocities for both components. It is shown that the mean orbital mass-to-luminosity ratio essentially depends on a morphological type of double galaxies and on a kind of interaction between pair components. Minimal values of mass-to-luminosity ratio take place for double galaxies having linear features of interaction, namely: bridges, tails; as well for pairs with compact and Markarian components, too. The empirical relation between mass-to-luminosity ratio and linear separation of pair components do not indicate on an existence of a hidden stellar coronae around double galaxies.     

Numerical Modeling of the Observed Wiechert Seismograph Magnification

—The paper presents results of experiments designed to measure the actual dynamic magnification of the Wiechert 1000 kg horizontal seismometer when excited by seismic waves. This is accomplished by comparing 51 digital records of seismic events recorded by the Wiechert and a well calibrated reference seismometer. The results obtained indicate that the magnification of the Wiechert seismometer is influenced by the interaction of its mass and frame, especially for high frequencies. This interaction has been modeled by considering a system of two coupled pendulums, yielding a theoretical dynamic magnification curve which exhibits main features of the observed magnification. The dis crepancy between the nominal and the actual response of the Wiechert seismograph may lead to errors in studies involving spectral analyses of recorded seismograms, and to overestimation of local earth quake magnitudes.     

Metasomatism and the formation of greisen in Grainsgill,Cumbria, England

A section of hornfelsed Skiddaw Slate adjacent to the margin of the Grainsgill Greisen is described and chemical analyses of the various rock types presented. The results confirm that hydrothermal fluids liberated K and Na from the granite to produce the greisen. This gave rise to K metasomatism of the adjacent hornfels and to retrogressive metamorphism over a radius of 200m from the intrusion. Na leached from the hornfels near the intrusion and that liberated during the formation of the greisen are probably the sources of the Na in the brines recorded in fluid inclusions within the nearby quartz-tungsten veins.     

Seismic Ground Motion Expected for the Eastern District of Naples

The seismic ground motion of a test area in the eastern district of Naples is computed with a hybrid technique based on the mode summation and the finite difference methods. This technique allows us the realistic modelling of source and propagation effects, including local soil conditions. In the modelling, we consider the 1980 Irpinia earthquake, a good example of strong shaking for the area of Naples, which is located about 90 km from the epicenter.The detailed geological setting is reconstructed from a large number of drillings. The sub-soil is mainly formed by alluvial (ash, stratified sand and peat) and pyroclastic materials overlying a pyroclastic rock (yellow neapolitan tuff), representing the neapolitan bedrock. The detailed information available on mechanical properties of the sub-soil and its geometry warrants the application of the sophisticated hybrid technique.As expected, the sedimentary cover causes an increase of the signal's amplitudes and duration. If thin peat layers are present, the amplification effects are reduced, and the peak ground accelerations are similar to those observed for the bedrock model. This can be explained by the backscattering of wave energy at such layers, that tend to seismically decouple the upper from the lower part of the structure.For SH-waves, the influence of the variations of the S-wave velocities on the spectral amplification is studied, by considering locally measured velocities and values determined from near-by down-hole measurements. The comparison between the computed spectral amplifications confirms the key role of an accurate determination of the seismic velocities of the different layers.The comparison performed between a realistic 2-D seismic response and a standard 1-D response, based on the vertical propagation of waves in a plane layered structure, shows considerable difference, from which it is evident that serious caution must be taken in the modelling of expected ground motion at a specific site.