Genesis of high Mg# andesites and the continental crust

The continental crust has an andesitic composition with high Mg/(Mg+Fe) and Ni contents which may be too high to have formed by differentiation of basaltic magmas. Instead, mantle-derived, high Mg# andesites (HMA) may form a substantial component of the crust. HMA may be produced by partial melting of previously depleted, subsequently metasomatised mantle peridotite. However, they are more likely produced by reaction between ascending melts and mantle peridotite. HMA are less common than basalts among lavas in modern island arcs, but may have been more common in the past, may be produced in specific environments (such as ridge subduction), may be more common among plutonic rocks in the lower and middle crust than among lavas at the surface, and may be selectively preserved during later erosion and subduction processes.     

Baroclinic instabilities in Jupiter's zonal flow

The baroclinic stability of Jupiter's zonal flow is investigated using a model consisting of two continuously stratified fluid layers. The upper layer, containing a zonal shear flow and representing the Jovian cloudy regions above p ～ 5 bars, is the same as Eady's (1949) model for the Earth. The lower layer has a relatively large but finite depth with a quiescent basic state, representing the deep Jovian fluid bulk below p ～ 5 bars. Due to the presence of the lower layer, the linearized non-dimensional growth rates are drastically reduced from the O(1) growth rates of the original Early model. Only very long wavelengths relative to the upper fluid's radius of deformation L₁ are unstable. Eddy transports of heat are also reduced relative to estimates based on scaling arguments alone. Since the hydrostatic approximation for the lower-layer perturbation breaks down at great depths, a second model is presented in which energy propagates downward in an infinitely deep lower fluid obeying the full linearized fluid equations. In this model, the growth rates are again very small, but now all wavelengths are unstable with maximum growth rates occurring for wavelengths O(1) relative to L₁. These results illustrate the importance for the upper-layer meteorology of the interface boundary condition with the lower fluid, which is radically different from the rigid lower boundary of the Earth's troposphere.     

Calculations of the effects of angular momentum on the early evolution of Jupiter

The planet Jupiter is assumed to have formed as a subcondensation in the solar nebula. The initial phase of its evolution is one of hydrostatic contraction with radiative energy transport. Calculations of evolutionary sequences through this phase are presented, including the effects of angular momentum. The calculations are carried out in two space dimensions under the assumptions of axial symmetry, constancy of angular velocity on cylindrical surfaces about the rotation axis, a pressure-density relation given by the polytrope of index 3, conservation of angular momentum, and a homogeneous composition. The results show that under certain physically reasonable initial distributions of density and angular momentum the formation of a central planet and a rotating circumplanetary envelope is possible, while under assumptions a point of instability is reached that probably results in the breakup of the condensation by fission into two or more parts. The models are discussed with reference to the present angular momenta of Jupiter and its regular satellites.     

Ce and Nd isotope geochemistry on island arc volcanic rocks with negative Ce anomaly: existence of sources with concave REE patterns in the mantle beneath the Solomon and Bonin island arcs

¹³⁸Ce/¹⁴²Ce isotope ratios in Cenozoic island arc volcanic rocks are reported for the first time, together with isotope ratios of Nd and Sr and abundances of REE, Ba and Sr. The island arc volcanics studies here are boninites from Chichijima, the Bonin Islands, and basalts and andesites from the Solomon Islands. REE patterns of the island arc volcanic rocks from the Solmon Islands and the Bonin Islands are confirmed to have negative Ce anomalies. It is also disclosed that the majority of these island arc volcanic rocks show mainly positive values for both _Ce and _Nd. It is shown that these _Ce and _Ce values can hardly be interpreted by simple mixing between MORB and oceanic or continental crustal rocks; the former have positive _Nd and negative _Ce and the latter have negative _Ce and positive or negative _Nd. Existence of sources having positive _Ce and _Nd values is strongly suggested. If the sources are assumed to have been fractionated from CHUR (chondritic uniform reservoir) at the early or middle Precambrian era, the sources from which the volcanics were derived are concluded to have kept concave REE patterns with larger (La/Ce)_N and smaller (Nd/Sm)_N ratios than chondritic values over a substantial period of time, until the time of Cenozoic magmatism forming island arc volcanic rocks in question. During the periods of the Cenozoic magmatic activities and their related events, Ce anomalies are considered to have been created. From Ce and Nd isotope ratios, however, it is difficult to determine which of the following processes was responsible for the Ce anomaly; the incorporation process of subducted oceanic crust into magma at the mantle or the slab dehydration and metasomatism process. Nevertheless, so far as Ce and Nd isotopic ratios are concerned, incorporation of oceanic sediments did not take place to any clearly detectable degree.     

H,O, Sr,Nd, and Pb isotope geochemistry of the Latir volcanic field and cogenetic intrusions,New Mexico,and relations between evolution of a continental magmatic center and modifications of the lithosphere

Over 200 H, O, Sr, Nd, and Pb isotope analyses, in addition to geologic and petrologic constraints, document the magmatic evolution of the 28.5–19 Ma Latir volcanic field and associated intrusive rocks, which includes multiple stages of crustal assimilation, magma mixing, protracted crystallization, and open- and closed-system evolution in the upper crust. In contrast to data from younger volcanic centers in northern New Mexico, relatively low and restricted primary ¹⁸O values (+6.4 to +7.4) rule out assimilation of supracrustal rocks enriched in ¹⁸O. Initial ⁸⁷Sr/⁸⁶Sr ratios (0.705 to 0.708), ¹⁸O values (-2 to-7), and ²⁰⁶Pb/²⁰⁴Pb ratios (17.5 to 18.4) of metaluminous precaldera volcanic rocks and postcaldera plutonic rocks suggest that most Latir rocks were generated by fractional crystallization of substantial volumes of mantle-derived basaltic magma that had near-chondritic Nd isotope ratios, accompanied by assimilation of crustal material in two main stages: 1) assimilation of non-radiogenic lower crust, followed by 2) assimilation of middle and upper crust by inter-mediate-composition magmas that had been contaminated during the first stage. Magmatic evolution in the upper crust peaked with eruption of the peralkaline Amalia Tuff (26 Ma), which evolved from metaluminous parental magmas. A third stage of late, roofward assimilation of Proterozoic rocks in the Amalia Tuff magma is indicated by trends in initial ⁸⁷Sr/⁸⁶Sr and ²⁰⁶Pb/²⁰⁴Pb ratios from 0.7057 to 0.7098 and 19.5 to 18.8, respectively, toward the top of the pre-eruptive magma chamber. Highly evolved postcaldera plutons are generally fine grained and are zoned in initial ⁸⁷Sr/⁸⁶Sr and ²⁰⁶Pb/²⁰⁴Pb ratios, varying from 0.705 to 0.709 and 17.8 to 18.6, respectively. In contrast, the coarser-grained Cabresto Lake (25 Ma) and Rio Hondo (21 Ma) plutons have relatively homogeneous initial ⁸⁷Sr/⁸⁶Sr and ²⁰⁶Pb/²⁰⁴Pb ratios of approximately 0.7053 and 17.94 and 17.55, respectively. ¹⁸O values for all the postcaldera plutons overlap those of the precaldera rocks and Amalia Tuff, except for those for two late-stage rhyolite dikes associated with the Rio Hondo pluton that have ¹⁸O values of-8.6 and-9.5; these dikes are the only Latir rocks which may be largely crustal melts.Chemical and isotopic data from the Latir field suggest that large fluxes of mantle-derived basaltic magma are necessary for developing and sustaining large-volume volcanic centers. Development of a detailed model suggests that 6–15 km of new crust may have been added beneath the volcanic center; such an addition may result in significant changes in the chemical and Sr and Nd isotopic compositions of the crust, although Pb isotope ratios will remain relatively unchanged. If accompanied by assimilation, crystallization of pooled basaltic magma near the MOHO may produce substantial cumulates beneath the MOHO that generate large changes in the isotopic composition of the upper mantle. The Latir field may be similar to other large-volume, long-lived intracratonal volcanic fields that fundamentally owe their origins to extensive injection of basaltic magma into the lower parts of their magmatic systems. Such fields may overlie areas of significant crustal growth and hybridization.     

Production of organic molecules in the outer solar system by proton irradiation: Laboratory simulations

In the past few years considerable attention has been given to the determination of likely compounds that could account for the various colors observed in the outer solar system: and to possible formation mechanisms for these compounds. Many experiments have been done using electrical discharges (Chadha, M. S., et al., 1971, Icarus15, 39) and ultraviolet light (Khare, B. N., and Sagan, C., 1973, Icarus20, 311) on mixtures of CH₄, NH₃, and H₂S, which are most likely the dominant minor constituents of the atmospheres of Jupiter, Saturn, Titan, and possibly the other satellites early in their histories. Colored polymers, usually brownish-red, have been produced in these experiments. With the passage of Pioneer 10 around Jupiter, there is another source of energy worthy of consideration, energetic protons (and electrons). Preliminary experiments to investigate the formation of colored polymers and other interesting molecules by the irradiation of gas mixtures by protons are discussed. Two to four Mev protons were used, with corresponding beam fluxes (as measured at 6R_J from the planet) equivalent to approximately 80 Earth years at Jupiter per hour of exposure. As in the other types of experiments, colored polymers have been produced. An important feature of this work is the presence or absence of absorption at 5 μm in the different materials produced; Titan is quite dark at this wavelength and Io is fairly bright. Such features may provide criteria for accepting or rejecting various materials produced in these experiments as reasonable coloring agents for the outer solar system.     

Point-connected rigid bodies in a topological tree

Identical equations of motion are shown to emerge for a system ofn+1 rigid bodies all interconnected byn points, each of which is common to two bodies, by means of each of the following derivation procedures, all of which employ a kinematical identity developed by Hooker and Margulies: The Hooker-Margulies/Hooker equations; Kane's quasicoordinate formulation of D'Alembert's principle; the combination of Lagrange's generalized coordinate equations and Lagrange's quasicoordinate equations; and the combination of Lagrange's generalized coordinate equations and the vector rotational equationM=H applied to the total system and resolved into a vector basis fixed in a reference body of the system. Thus the previously published Hooker-Margulies/Hooker equations are shown to be the natural result of several derivation procedures other than the Newton-Euler method originally used, provided that the central kinematical identity of the original derivation of Hooker and Margulies is employed.     

Incorporating solid solutions in reactive transport equations using a kinetic discrete-composition approach

A new approach is proposed for incorporating solid solution reactions into mass conservation equations describing reaction paths in both closed and open systems. The method is applicable to problems involving advective, dispersive, and diffusive transport in a porous medium. By representing the continuously variable solid solution composition with a discrete set of stoichiometric solids that span composition space, combined with a kinetic formulation of their rates of reaction, a self-determining spatial and temporal evolution of the solid solution concentration and composition is obtained. It is demonstrated that equilibrium of an aqueous solution with a stoichiometric solid derived from a solid solution corresponds to equilibrium of the solid solution itself if and only if equilibrium of the stoichiometric solid is stable. One advantage of this approach is that it is unnecessary to introduce any additional compositional variables to represent the solid solution. Discretization may be over the entire range of composition space, or over some subset depending on the system. A major consequence of the kinetic discrete-composition solid solution representation is that modeling solid solutions is similar to modeling pure mineral phases with the exception of a weighting factor applied to reaction rates of stoichiometric solids corresponding to a common solid solution. With this approach, precipitation leads to a discrete zonation of the solid solution that approximates the continuous variation in composition expected for the actual solid solution. The approach is demonstrated for a hypothetical ideal and non-ideal binary solid solution A_xB_1−xC for a reaction path formulation and reactive transport involving advection and diffusion.     

Semi-analytical solutions of a contaminant transport equation with nonlinear sorption in 1D

A new method to determine semi-analytical solutions of one-dimensional contaminant transport problem with nonlinear sorption is described. This method is based on operator splitting approach where the convective transport is solved exactly and the diffusive transport by finite volume method. The exact solutions for all sorption isotherms of Freundlich and Langmuir type are presented for the case of piecewise constant initial profile and zero diffusion. Very precise numerical results for transport with small diffusion can be obtained even for larger time steps (e.g., when the Courant-Friedrichs-Lewy (CFL) condition failed).