Intrusion and differentiation of granitic magma at a high level in the crust: the Puscao pluton, Lima Province, Peru

Regression surface techniques have been used on partial whole-rockchemical analyses from the Puscao pluton to establish the three-dimensionalnature of the chemical variations within this body, which wasemplaced at a depth of between 4 km and 8 km within the crust.Block diagrams produced by the computer program are consistentwith the hypothesis that during the period preceding the totalconsolidation of the magma, the pluton has undergone in situdifferentiation to give a layered structure. Statistical analysis confirms that there is a large scale variationwithin the pluton and that local variation, such as might beexpected to result from the contamination of the Puscao magmaby assimilated volcanic material, did not contribute significantlyto the overall variation. Evidence is presented which suggests that the layering withinthe Puscao pluton was disturbed as a result of the emplacementof the much later Cañhas pluton, and that this disturbancehas the form of a simple doming centred upon the site of thepresent Cañas intrusion.     

Coronal disturbances

Time lapse movies acquired with the Sacramento Peak Observatory 6-in filter coronagraph very infrequently show dramatic transient events. Many of these events (at 1.04 < r/r < 1.20) are correlated with H activity at the limb, especially with large surges and eruptive prominences. Metric or centrimetric radio bursts are found associated with certain types of 5303 coronal events. The heights to which such disturbances will propagate can often be estimated from the green line movie data. We report the results of a study of all the examples of such transients observed during the last 16 years.     

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.     

A regional model for studies of atmosphere-ice-ocean interaction in the western Arctic

Summary Asa step in the development of a fully coupled regional model of the atmosphere-ice-ocean system, atmospheric and sea ice models have been adapted to a western Arctic domain centered on the Bering Strait. Lateral boundary conditions derived from operational analyses drive the models through simulations on grids having horizontal resolutions of 21 km and 7 km. Sensitivities to the presence of sea ice are large after only 48 hours, by which time the surface temperatures in the Bering and Chukchi Seas are 10–15°C higher without sea ice than with sea ice. The temperatures, in turn, modify the fields of sea level pressure, surface wind and precipitation. By influencing the surface wind stress through the static static stability, the surface state feeds back to the surface momentum exchange, ice/ocean transport, and the rate of formation of new ice. The results also show a resolution-dependence of the surface winds, precipitation rates and new ice formation rates, particularly in areas in which the coastal configuration and topography are spatially complex. The experiments will be augmented by the implementation of an ocean model on the same grids.With 12 Figures     

Solubility of alumina in orthopyroxene plus spinel as a geobarometer in complex systems. Applications to spinel-bearing alpine-type peridotites

The addition of Fe and Cr to the simple system MgO-SiO₂-Al₂O₃ markedly affects the activities of phases involved in the equilibrium

Refinement of the modified time-to-failure method for intermediate-term earthquake prediction

The modified time-to-failure method for intermediate-term earthquake prediction utilizes empirical relationships to reduce the number of unknown parameters providing a stable and unique solution set. The only unknown parameters in the modified time-to-failure method are the time and size of the impending main shock. The modified time-to-failure equation is used to model the precursory events and a prediction contour diagram is constructed with the magnitude and time-of-failure as the axes of the diagram. The root-mean-square (rms) is calculated for each set of time and magnitude on the prediction diagram representing the difference between the model (calculated) acceleration and the actual accelerated energy release of the precursory events. A small region, corresponding to the low rms region on the diagram, defines the prediction. The prediction has been shown to consistently under-estimate the magnitude and over-estimate the time-of-failure. These shortcomings are caused by an underestimation in energy release of the modified time-to-failure equation at the very end of the sequence. An empirical correction can be applied to the predicted results to minimize this problem. A main shock location search technique has been developed for use with the modified time-to-failure method. The location technique is used to systematically search an earthquake catalog and identify locations corresponding to precursory sequences that display accelerated energy releases. It has shown good results when applied in retrospective predictions, and is essential for the practical application of the modified time-to-failure method. In addition, an observed linear characteristic in long-term energy release can be used to minimize false predictions. The refined empirical relationships that eliminate or constrain unknown constants used in the modified time-to-failure method and the main shock location search technique are used in a practical application in the New Madrid Seismic Zone (NMSZ). The NMSZ, which is over due for a magnitude 6 event according to recurrence rates (Johnston and Nava, 1985), makes this region ideal for testing the method. One location was identified in the NMSZ as a high risk area for an event in the magnitude 4.5 range. The prediction, if accurate, is of scientific interest only because of the relatively small size of the main shock.