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.     

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.     

Residual deformations and seismic hazard of fault zones in underground workings

The results of long-term measurements of residual deformations in the area of the Severomuiskii tunnel during its construction are given. Comparison of spatial and temporal distributions between deformations and earthquakes shows that they are interrelated. The nonlinear behavior of disintegrated and waterlogged rocks within tectonic fault zones should be taken into account even for moderate intensity earthquakes, since they are accompanied by ground subsidence. Underground workings within such sites are subjected to loads that are greater than expected according to current hypotheses, which are the basis for calculations of rock pressure. It is established that the most probable mechanism of ground subsidence and residual deformation in underground workings could be the gravity load of the overlying strata of disintegrated and waterlogged rocks within a fault zone, which leads to instability of the rocks. Approximate stress estimations are made on the basis of the results of measurements of residual deformations of roof supports (linings).     

Impact of Sea-level Rise and Storm Surges on a Coastal Community

A technique to evaluate the risk of storm tides (the combination of a storm surge and tide) under present and enhanced greenhouse conditions has been applied to Cairns on the north-eastern Australian coast. The technique combines a statistical model for cyclone occurrence with a state-of-the-art storm surge inundation model and involves the random generation of a large number of storm tide simulations. The set of simulations constitutes a synthetic record of extreme sea-level events that can be analysed to produce storm tide return periods. The use of a dynamic storm surge model with overland flooding capability means that the spatial extent of flooding is also implicitly modelled. The technique has the advantage that it can readily be modified to include projected changes to cyclone behaviour due to the enhanced greenhouse effect. Sea-level heights in the current climate for return periods of 50, 100, 500 and 1000 years have been determined to be 2.0 m, 2.3 m, 3.0 m and 3.4 m respectively. In an enhanced greenhouse climate (around 2050), projected increases in cyclone intensity and mean sea-level see these heights increase to 2.4 m, 2.8 m, 3.8 m and 4.2 m respectively. The average area inundated by events with a return period greater than 100 years is found to more than double under enhanced greenhouse conditions.     

Peculiarities of Dynamic Slip Nucleation in a Thin Granular Layer

Izvestiya, Physics of the Solid Earth - Abstract—Granular media determine the dynamics of many natural systems including faults in the Earth’s crust. The paper addresses the laboratory...