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.     

Réflexions sur la contingence des lois de la Nature

There are two fundamental conceptions concerning the laws of nature — these laws are eternal or they may vary in the course of time. We suggest an experiment in order to find an intermediary thesis — i.e. to discover that law, if it really exists, according to which a given law of nature varies, or, at least, to find out the direction in which this variation is performed.     

A Spectrophotometric Study of the Be Star HD 164447

Spectra of the Be star HD 164447 have been studied. The equivalent widths and the intensities of the measured lines are given. No emission is apparent on the plates. Applying the thin layer theory some parameters of the stars atmosphere were derived.     

Measures of Parameter Uncertainty in Geostatistical Estimation and Geostatistical Optimal Design

Studies of site exploration, data assimilation, or geostatistical inversion measure parameter uncertainty in order to assess the optimality of a suggested scheme. This study reviews and discusses measures for parameter uncertainty in spatial estimation. Most measures originate from alphabetic criteria in optimal design and were transferred to geostatistical estimation. Further rather intuitive measures can be found in the geostatistical literature, and some new measures will be suggested in this study. It is shown how these measures relate to the optimality alphabet and to relative entropy. Issues of physical and statistical significance are addressed whenever they arise. Computational feasibility and efficient ways to evaluate the above measures are discussed in this paper, and an illustrative synthetic case study is provided. A major conclusion is that the mean estimation variance and the averaged conditional integral scale are a powerful duo for characterizing conditional parameter uncertainty, with direct correspondence to the well-understood optimality alphabet. This study is based on cokriging generalized to uncertain mean and trends because it is the most general representative of linear spatial estimation within the Bayesian framework. Generalization to kriging and quasi-linear schemes is straightforward. Options for application to non-Gaussian and non-linear problems are discussed.     

Directional dependence of P- and S-wave propagation and polarization in foliated rocks from the Kola superdeep well: Evidence from laboratory measurements and calculations based on TOF neutron diffraction

We have measured P- and S-wave velocities on two amphibolite and two gneiss samples from the Kola superdeep borehole as a function of pressure (up to 600 MPa) and temperature (up to 600 °C). The velocity measurements include compressional (Vp) and shear wave velocities (Vs1, Vs2) propagating in three orthogonal directions which were in general not parallel to inherent rock symmetry axes or planes. The measurements are accompanied by 3D-velocities calculations based on lattice preferred orientation (LPO) obtained by TOF (Time Of Flight) neutron diffraction analysis which allows the investigation of bulk volumes up to several cubic centimetres due to the high penetration depth of neutrons. The LPO-based numerical velocity calculations give important information on the different contribution of the various rock-forming minerals to bulk elastic anisotropy and on the relations of seismic anisotropy, shear wave splitting, and shear wave polarization to the structural reference frame (foliation and lineation). Comparison with measured velocities obtained for the three propagation directions that were not in accordance with the structural frame of the rocks (foliation and lineation) demonstrate that for shear waves propagating through anisotropic rocks the vibration directions are as important as the propagation directions. The study demonstrates that proper measurement of shear wave splitting by means of two orthogonal polarized sending and receiving shear wave transducers is only possible when their propagation and polarization directions are parallel and normal to foliation and lineation, respectively.     

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.     

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.     

High-precision repeat-groundtrack orbit design and maintenance for Earth observation missions

The focus of this paper is the design and station keeping of repeat-groundtrack orbits for Sun-synchronous satellites. A method to compute the semimajor axis of the orbit is presented together with a station-keeping strategy to compensate for the perturbation due to the atmospheric drag. The results show that the nodal period converges gradually with the increase of the order used in the zonal perturbations up to \(J_{15}\). A differential correction algorithm is performed to obtain the nominal semimajor axis of the reference orbit from the inputs of the desired nodal period, eccentricity, inclination and argument of perigee. To keep the satellite in the proximity of the repeat-groundtrack condition, a practical orbit maintenance strategy is proposed in the presence of errors in the orbital measurements and control, as well as in the estimation of the semimajor axis decay rate. The performance of the maintenance strategy is assessed via the Monte Carlo simulation and the validation in a high fidelity model. Numerical simulations substantiate the validity of proposed mean-elements-based orbit maintenance strategy for repeat-groundtrack orbits.     

Formation Mechanism of Mud Bank Along the Southwest Coast of India

Mud bank formation during the southwest monsoon along the southwest coast of India remains an enigma to the researchers and coastal community in spite of several earlier studies. The present study attempts to unravel the mystery through a high-frequency, season-long time-series observation at Alappuzha, located at the southern part of the west coast of India, a region of frequent occurrence of mud bank. Using 7-month-long weekly time-series observation, we identified strong winds and high waves associated with onset of the southwest monsoon and subsequent three episodic atmospheric low-pressure events (LPEs).With the help of in situ time-series data, we show that the strong winds and high waves associated with southwest monsoon pre-conditions the near shore bottom sediment to bring it into suspension. The high amplitude waves associated with the southwest monsoon, while propagating from the deep water to shallow water region, interact with the bottom initiating bottom-sediment movement and its suspension due to wave refraction and shoaling. The sporadic occurrence of the atmospheric LPEs enhances the process of suspension of bottom sediment in the near shore region leading to the formation of fluid mud. Simulations with a cohesive sediment transport model yielded realistic estimates of sediment transport, in the presence of an onshore current, a pre-requisite for transporting the fluid mud toward the coast. The prevailing onshore upwelling current during the southwest monsoon provides the favorable pre-requisite conditions for transporting the fluid mud through depression channel network towards the coast. Once sufficient quantity and thickness of fluid mud is accumulated in the near shore region, it acts as a wave damper for subsequent high monsoon waves, as indicated by the time-series wave data, leading to the formation of tranquil mud bank region. Depression channel networks extending from the shelf to the coast off Alappuzha, Kochi, Ponnani, Beypore, and Ullal were found in the bathymetric charts, thus explaining why mud banks occur only at few locations in spite of the prevalence of similar monsoon conditions.