Feasibility of AFC models for the petrogenesis of calc-alkaline magma series

 One of the most popular petrogenetic concepts is that of simultaneous assimilation and fractional crystallisation (AFC). The equations governing these processes are well known and widely applied. The power and attraction of the model are that, with the minimum of constraints, petrogenetic schemes can be erected for essentially any group of rocks. These models are usually based on the observed variations in a small number of trace elements and/or isotopic systems. In this paper we have generated well-constrained AFC models for the Variscan, calc-alkaline Quérigut complex in the French Pyrénées. Our models are based on geological constraints on the identities of the parent magma and candidate assimilants, and the Sm-Nd and Rb-Sr isotope systematics of the rocks. The modelling is successful in mimicking the isotope array for the actual rocks. However, the results are massively inconsistent with those produced using the major- and trace-element concentrations. Our major conclusion is that AFC modelling should be undertaken only with robust initial constraints on the compositions of the parent magma and the proposed assimilant, with due respect for crystal-liquid partition coefficients, in so far as they are known. The results should then be accepted only if the isotope-based model parameters are consistent with those derived from major- and trace-element variations. Received: 19 August 1994/Accepted: 23 March 1995     

Equilibrium coexistence of three amphiboles

Electron probe and wet chemical analyses of amphibole pairs from the sillimanite zone of central Massachusetts and adjacent New Hampshire indicated that for a particular metamorphic grade there should be a restricted composition range in which three amphiboles can coexist stably. An unequivocal example of such an equilibrium three amphibole rock has been found in the sillimanite-orthoclase zone. It contains a colorless primitive clinoamphibole, space group P2₁/m, optically and chemically like cummingtonite with blue-green hornblende exsolution lamellae on (100) and (¯101) of the host; blue-green hornblende, space group C2/m, with primitive cummingtonite exsolution lamellae on (100) and (¯101) of the host; and pale pinkish tan anthophyllite, space group Pnma, that is free of visible exsolution lamellae but is a submicroscopic intergrowth of two orthorhombic amphiboles. Mutual contacts and coarse, oriented intergrowths of two and three host amphiboles indicate the three grew as an equilibrium assemblage prior to exsolution. Electron probe analyses at mutual three-amphibole contacts showed little variation in the composition of each amphibole. Analyses believed to represent most closely the primary amphibole compositions gave atomic proportions on the basis of 23 oxygens per formula unit as follows: for primitive cummingtonite (Na_0.02Ca_0.21 ^– Mn_0.06Fe²⁺ _2.28Mg_4.12Al_0.28) (Al_0.17Si_7.83), for hornblende (Na_0.35Ca_1.56Mn_0.02Fe_1.71Mg_2.85Al_0.92) (Al_1.37Si_6.63), and for anthophyllite (Na_0.10Ca_0.06Mn_0.06Fe_2.25Mg_4.11Al_0.47) (Al_0.47Si_7.53). The reflections violating C-symmetry, on X-ray single crystal photographs of the primitive cummingtonite, are weak and diffuse, and suggest a partial inversion from a C-centered to a primitive clinoamphibole. Single crystal photographs of the anthophyllite show split reflections indicating it is an intergrowth of about 80% anthophyllite and about 20% gedrite which differ in their b crystallographic dimensions. Split reflections are characteristic of all analyzed orthorhombic amphiboles so far examined from Massachusetts and New Hampshire except the most aluminous gedrites, and the relative intensity of the gedrite reflections is roughly proportional to the degree of Na and Al substitution. Thin sections of a few of these anthophyllite specimens show lamellae parallel to (010) that are just resolved with a high power objective.Publication approved by the Director, U.S. Geological Survey.     

Stepwise mass extinctions: the case for the Late Cenomanian event

The various extinctions of marine taxa in the Late Cenomanian can be shown as conforming to the 'step-wise extinction' model. The order of such extinctions is, however, compatible with the presence of an expanded oxygen minimum zone in the world ocean. The occurrence of iridium anomalies at this level in the succession may, however, indicate a volcanogenic control on the global environment. Detailed time-analysis of the Late Cenomanian event has been used to produce a correlation of some key successions.     

A three-dimensional numerical model of particulate transport for coastal waters

A Lagrangian marker particle in Eulerian finite difference cell solution to the three-dimensional incompressible mass transport equation was developed for predicting particulate transport in coastal and estuarine waters. Special features of the solution procedure include a finite difference grid network which translates horizontally and vertically with the mean particle motion and expands with the dispersive growth of the marker particle cloud. The cartesian vertical coordinate of the three-dimensional mass transport equation has been transformed, using instantaneous water column depth to allow adaptation to flow situations with a temporally and spatially varying bottom topography and free surface. Results from this model for turbulent diffusion and advection of a uniform plug flow of sediment in an unbounded uniform flow field with various sediment settling velocities were in excellent agreement with the corresponding analytic solutions. Using current information from a two-dimensional vertically averaged hydrodynamic's model, the model was utilized to predict the long term diffusion and advection of dilute neutrally and negatively buoyant suspended sediment clouds resulting from a hypothetical instantaneous release of dredge spoil waste at Brown's Ledge in Rhode Island Sound.     

A Model of Conductive Heat Flow in Forest Edges and Fragmented Landscapes

Although the creation of edges during forest fragmentation can have important abiotic and biotic impacts, especially under conditions of future climate change, mechanistic models of edge effects have not been forthcoming. A simple numerical model of two-dimensional heat flow is developed and applied to a vertical forest/clearcut edge profile and to simulated fragmented landscapes. Height-specific thermal diffusivity and conductivity in the forest were assumed to vary in proportion to foliage densities measured in the central Amazon. In the edge profile, the clearcut that abutted the edge served as a heat source and its temperature was maintained at a constant value higher than in the initially cooler forest. In the fragmented landscapes, simulated treefall gaps were heat sources whose temperature varied with sun movements during the day. Gap frequency was varied so as to approximate the gap coverage observed in selectively logged forests. In one set of simulations, temperature in the openings was systematically varied; in another, thermal diffusivity of the forest was varied. Along the edge profile, high temperatures in the clearcut were rapidly transmitted into the upper canopy due to additive edge effects. Temperatures in the forest understory were also very sensitive to clearcut temperatures due to relatively sparse understory foliage. An overall increase in forest diffusivity led to markedly higher temperatures close to the edge and a more even temperature distribution among height strata. In fragmented landscapes, total gap coverage and additivity from neighboring gaps strongly influenced forest temperatures. At low conductivities, heat flowed only into the forest close to the gaps and hence forest temperature increased almost linearly with gap area. However, at high conductivities, heat flowed far into the forest and forest temperature varied as a function of gap density in the surrounding neighborhood. Because of these additive effects, slight increases in total gap area led to disproportionate changes in the thermal profile of the landscape. These results have important implications for the conservation of forest ecosystems.     

Some methods for statistical analysis of multimodal distributions and their application to grain-size data

Some of the methods used in the resolution of mixed normal distributions are discussed under three headings: analytical, graphical, and numerical methods. Attention is given to their applicability in the analysis of grain-size data as derived from sieving. Comparisons are made by applying several methods to published data. It is concluded that the numerical methods offer most scope, especially the method of nonlinear least squares. Some analyses of beach sediments, using this method, are presented. The adoption of a convention for the number of individuals in the sample increases ease of interpretation.