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.     

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.     

Revised identifications in the 5C2 survey

In the present paper revised distances of optical objects from radio sources of the 5C2 survey are given. A statistical investigation of the data is given for blue objects and galaxies by the statistical method of the “first neighbour”. The identification rate on blue plates for both the blue objects and galaxies amounts to about 40% out of the total number of 26 identifications. For every blue object and galaxy which are proposed as an identification the statistical reliability is given.     

Fe3+ in shocked olivine crystals of the ALHA 77005 meteorite

ALHA 77005, a shocked achondrite of the shergottite group, is unique in containing cumulus olivine crystals which display a brownish color in thin section. Mössbauer, EPR, and optical spectroscopic studies show the presence of Fe³⁺ in the olivine. Approximately 4.5 wt.% of the total iron in the olivine (24 wt.% as FeO) is in the trivalent state. The Fe³⁺ ions preferentially occupy M2 lattice positions. Charge transfer between Fe³⁺ and Fe²⁺ as well as between Fe³⁺ and oxygen ligands causes a broad absorption band in the optical spectrum at high energies which is the reason for the olivine color. It is concluded that the color of the ALHA 77005 olivine is preterrestrial and may be the result of shock-induced oxidation of the olivine crystals on the meteorite parent body.     

Laboratory study of deep-water breaking waves

In an attempt to elucidate the mechanics of deep-water wave breaking, a variety of breaking waves, including spilling and plunging waves, of different length scales and geometries was studied. The waves were generated through wave-wave interactions using wave packets with constant-steepness components, constant-amplitude components, and also components following the Pierson-Moskowitz distribution. Wave steepening prior to breaking were found to cause an increase in the high frequency spectral slope of the wave spectrum. The slopes were correlated to the type of breaking and the intensity of the breaking. The energy loss through breaking varied with the spectral characteristics of the wave packet. On the other hand, it was also noted that, irrespective of the wave packet, the losses were from the higher frequency end of the first harmonics.     

First Byurakan Spectral Sky Survey. Stars of Late Spectral Types. XI. Zone -7° ≤ δ ≤ -3°

The eleventh list of faint late M and carbon type stars detected on the plates of the First Byurakan Spectral Survey in zone -7° -3° and covering about 1000 square degrees is presented. From 126 detected stars, 88 are newly discovered objects: they are 6 carbon stars, 8 carbon star candidates, and 74 M-type stars; among the latter 38 (26 PSC + 12 FSC) are unclassified IRAS sources, and one object is an unclassified ROSAT source. Distances to the 6 newly discovered early-type carbon stars are estimated. Equatorial coordinates, red magnitudes, and spectral classes determined from the Palomar E-charts are provided. The lack of optical counterparts on Palomar O and E maps for two detected late M-type stars indicates a large variability in brightnesses of these objects (amplitude not smaller than 7.0 magnitude).     

Risk Management for Petroleum Reservoir Production: A Simulation-Based Study of Prediction

We consider numerical solutions of the Darcy and Buckley–Leverett equations for flow in porous media. These solutions depend on a realization of a random field that describes the reservoir permeability. The main content of this paper is to formulate and analyze a probability model for the numerical coarse grid solution error. We explore the extent to which the coarse grid oil production rate is sufficient to predict future oil production rates. We find that very early oil production data is sufficient to reduce the prediction error in oil production by about 30%, relative to the prior probability prediction.     

The role of prior model calibration on predictions with ensemble Kalman filter

This paper, based on a real world case study (Limmat aquifer, Switzerland), compares inverse groundwater flow models calibrated with specified numbers of monitoring head locations. These models are updated in real time with the ensemble Kalman filter (EnKF) and the prediction improvement is assessed in relation to the amount of monitoring locations used for calibration and updating. The prediction errors of the models calibrated in transient state are smaller if the amount of monitoring locations used for the calibration is larger. For highly dynamic groundwater flow systems a transient calibration is recommended as a model calibrated in steady state can lead to worse results than a noncalibrated model with a well-chosen uniform conductivity. The model predictions can be improved further with the assimilation of new measurement data from on-line sensors with the EnKF. Within all the studied models the reduction of 1-day hydraulic head prediction error (in terms of mean absolute error [MAE]) with EnKF lies between 31% (assimilation of head data from 5 locations) and 72% (assimilation of head data from 85 locations). The largest prediction improvements are expected for models that were calibrated with only a limited amount of historical information. It is worthwhile to update the model even with few monitoring locations as it seems that the error reduction with EnKF decreases exponentially with the amount of monitoring locations used. These results prove the feasibility of data assimilation with EnKF also for a real world case and show that improved predictions of groundwater levels can be obtained.