The nature and significance of the feedback of changes in terrestrial vegetation on global atmospheric and climatic change

The potential feedback on global atmospheric and climate change of climate-driven changes in terrestrial vegetation is examined by systematically relating the surface exchanges of energy, mass and momentum to two dimensions of vegetation, structure and taxonomy, such that the significance of climate driven changes in these characteristics can be assessed. A detailed quantitative understanding of this feedback is an important prerequisite to realistic and dynamic representations of the Earth's surface within general circulation and biological models (GCMs and GBMs). Without realistic representations of terrestrial vegetation within these models, any forecasts of future climates by these models must be suspect.Several general conclusions are drawn. The first is that the indirect feedbacks, those associated with the clouds and aerosols of the planetary boundary layer, appear to be very powerful but as yet their behaviour and connections with the underlying surface are both poorly understood and captured within GCMs.The physical structure of vegetation, the disposition of biomass in 3-D, is the characteristic that most strongly influences the exchange of momentum (via aerodynamic roughness) and solar radiation (via albedo). Vegetation structure and species composition determine the most important of the mass exchanges, evapotranspiration. Of all of the surface exchanges, the parameterization of evapotranspiration (E ) and the simulation of the water balance over time is the most critical.Lastly, the problems of scaling and spatial heterogeneity, the sub-grid variability of the modellers, looms as a difficult, but not insoluble, problem. It remains a critical problem however, and the detailed parameterization of the various big leaf models stands in absurd contrast to the simplistic generalization of the spatial heterogeneity of terrestrial landscapes.Plant ecologists can contribute to the task of improving the representation of vegetated landscapes within GCMs. There is need to simply and unify the way in which vegetation can be grouped at landscape scales. A classification that is based on function rather than phylogeny is required. The definition of Vegetation Functional Types (VFTs) would expedite research on both the impact of, and feedback on, climate change.     

Diffuse radiation,twilight, and photochemistry — I

A photochemical scheme which includes a detailed treatment of multiple scattering up to solar zenith angles of 96° (developed for use in a GCM) has been used to study partitioning within chemical families. Attention is drawn to the different zenith angle dependence of diffuse radiation for the two spectral regions <310 nm and >310 nm. The effect that this has on the so-called 40 km ozone problem is discussed. The importance of correctly including multiple scattering for polar ozone studies is emphasised.     

Flow directions in ash-flow tuffs: a comparison of geological and magnetic susceptibility measurements, Tshirege member (upper Bandelier Tuff), Valles caldera, New Mexico, USA

This study concludes that the elongation axis (K ₁) of the ellipsoid of anisotropic magnetic susceptibility (AMS) is a suitable proxy for flow axis in ashflow tuffs. 153 oriented samples (176 specimens) were studied from 18 sites in the 1.1 Ma Tshirege member of the Bandelier Tuff. These sites are distributed around the Valles caldera at distances of 5–25 km outside of the rim.K ₁ axes correlate well with postulated radial flow axes at 13 sites.K ₁ also agrees with measured geological flow indicators, mainly imbricated larger clasts, at 7 sites. At 2 of the 5 sites where significant disagreement is seen between theoretical radial flow directions and measuredK ₁ axes, theK ₁ axes correspond well with geological flow indicators, indicating that the divergence of flow from the predicted radial flow pattern is real. Two major topographic buttresses are suggested as the cause of flow divergence for the Tshirege ash flows: the San Pedro buttress northwest of the caldera, and the San Miguel buttress in the southeast. In situK ₁ axes plunge about 7° toward the source at two-thirds of the sites; therefore the plunge ofK ₁ is a plausible in situ indicator for thedirection of flow. Multiple flow zones in sections of several meters thickness indicate changes of flow direction that are both rapid and large during ash-flow emplacement. These observations raisre the question of how best to represent mean flow directions in ash-flow sheets: by eigenvector methods, by vector-sum methods, or by modes. A method for measuring imbrication of larger clasts using apparent dips in vertical joints is outlined. Imbrication, determined in this way at one-third of the sites, dips toward the source, i.e., up-flow. The minimum (K ₃) axis of the AMS ellipsoid correlates with the flow foliation rather than with the larger clast imbrication. The flow axes of ash flows correspond with theK ₁ axes, not with the declination ofK ₃ axes as suggested by some authors. Initial dip of the sampled ash flows is not large and does not affect the paleomagnetic remanence direction, which is reversed with a mean ofD=173.5°,I=-38.4°, ₉₅=3.4°N=18. This mean is not different at the 95% confidence level from that of earlier workers. The mean pole, at 098.0°E, 74.8°N,A ₉₅=3.3°,N=18, is about 15° far-sided relative to the expected time-averaged geomagnetic pole, suggesting a history of emplacement too short to adequately average secular variation.     

Anelasticity of the crust and upper mantle beneath north and central India from the inversion of observed Love and Rayleigh wave attenuation data

The fundamental mode Love and Rayleigh waves generated by earthquakes occurring in Kashmir, Nepal Himalaya, northeast India and Burma and recorded at Hyderabad, New Delhi and Kodaikanal seismic stations are analysed. Love and Rayleigh wave attenuation coefficients are obtained at time periods of 15–100 seconds, using the spectral amplitude of these waves for 23 different paths along northern (across Burma to New Delhi) and central (across Kashmir, Nepal Himalaya and northeast India to Hyderabad and Kodaikanal) India. Love wave attenuation coefficients are found to vary from 0.0003 to 0.0022 km^–1 for northern India and 0.00003 km^–1 to 0.00016 km^–1 for central India. Similarly, Rayleigh wave attenuation coefficients vary from 0.0002 km^–1 to 0.0016 km^–1 for northern India and 0.00001 km^–1 to 0.0009 km^–1 for central India. Backus and Gilbert inversion theory is applied to these surface wave attenuation data to obtainQ ^–1 models for the crust and uppermost mantle beneath northern and central India. Inversion of Love and Rayleigh wave attenuation data shows a highly attenuating zone centred at a depth of 20–80 km with lowQ for northern India. Similarly, inversion of Love and Rayleigh wave attenuation data shows a high attenuation zone below a depth of 100 km. The inferred lowQ value at mid-crustal depth (high attenuating zone) in the model for northern India can be by underthrusting of the Indian plate beneath the Eurasian plate which has caused a low velocity zone at this shallow depth. The gradual increase ofQ ^–1 from shallow to deeper depth shows that the lithosphere-asthenosphere boundary is not sharply defined beneath central India, but rather it represents a gradual transformation, which starts beneath the uppermost mantle. The lithospheric thickness is 100 km beneath central India and below that the asthenosphere shows higher attenuation, a factor of about two greater than that in the lithosphere. The very lowQ can be explained by changes in the chemical constitution taking place in the uppermost mantle.     

Na-rich carbonate inclusions in perovskite and calzirtite from the Guli intrusive Ca-carbonatite,polar Siberia

Carbonate phases, some rich in Na₂O and comparatively rich in SrO and BaO, occur as inclusions in perovskite and calzirtite (Ca₂Zr₅Ti₂O₁₆) in the carbonatite of the Guli complex, Siberia. This is the first record of alkali carbonates, akin to nyerereite [Na₂Ca(CO₃)₂], in plutonic igneous rocks. The inclusion populations suggest that the parental magma of the complex was Ca-rich but developed Na-rich differentiates during the latest stages. This points to the dominant calcic carbonatites of the complex not being derivatives of alkali-rich parental carbonatites. These alkali-rich carbonate inclusions (and rare inclusions of djerfisherite) have been preserved due to the resistance of perovskite and calzirtite to processes of leaching, hydrothermal alteration and weathering.     

High pressure experimental calibration of the olivine-orthopyroxene-spinel oxygen geobarometer: implications for the oxidation state of the upper mantle

Synthetic spinel harzburgite and lherzolite assemblages were equilibrated between 1040 and 1300° C and 0.3 to 2.7 GPa, under controlled oxygen fugacity (f _{O 2}). f _{O 2} was buffered with conventional and open double-capsule techniques, using the Fe−FeO, WC-WO₂-C, Ni−NiO, and Fe₃O₄−Fe₂O₃ buffers, and graphite, olivine, and PdAg alloys as sample containers. Experiments were carried out in a piston-cylinder apparatus under fluid-excess conditions. Within the P-T-X range of the experiments, the redox ratio Fe³⁺/ΣFe in spinel is a linear function of f _{O 2} (0.02 at IW, 0.1 at WCO, 0.25 at NNO, and 0.75 at MH). It is independent of temperature at given Δlog(f _{O 2}), but decreases slightly with increasing Cr content in spinel. The Fe³⁺/ΣFe ratio falls with increasing pressure at given Δlog(f _{O 2}), consistent with a pressure correction based on partial molar volume data. At a specific temperature, degree of melting and bulk composition, the Cr/(Cr+Al) ratio of a spinel rises with increasing f _{O 2}. A linear least-squares fit to the experimental data gives the semi-empirical oxygen barometer in terms of divergence from the fayalite-magnetite-quartz (FMQ) buffer:

An Example of Island-Arc Petrogenesis: Geochemistry and Petrology of the Southern Luzon Arc, Philippines

Volcanism throughout the Luzon arc is associated with eastwardsubduction of the South China Sea floor along the Manila Trench.The southern section of the arc, the focus of this study, extendsfrom the Lingayen-Dingalan fault to the small islands just southof Luzon. Two segments appear to exist along this section ofthe arc the northern Bataan and southern Mindoro segments whichare separated by the Macolod Corridor. The volcanic rocks have typical arc phenocryst mineralogies:olivine, clinopyroxene, plagioclase, and titanomagnetite inthe most mafic rocks and clinopyroxene, plagioclase, orthopyroxene,titanomagnetite, ? amphibole in the more felsic samples. Complexzoning, sieve textures, and decoupling of incompatible traceelements suggest that processes such as assimilation have takenplace. The rocks from the study area range from basalts to rhyolitesand show typical calc-alkaline features. The rocks of the MacolodCorridor and Mindoro segment are particularly enriched in largeion lithophile elements (LILE), light rare earth elements (LREE),and radiogenic Sr compared with the Bataan segment. The datafall within the mantle array on Sr-Nd isotopic diagrams andgrade toward higher Sr and lower Nd isotopic values from northto south. A likely source for the volcanics of this study is either amid-ocean ridge basalt (MORB)-type mantle that undergoes higherdegrees of partial melting than regions involved in MORB generationor a previously depleted source. We suggest that the high fieldstrength element (HFSE) anomalies have been derived throughdifferential element partitioning during fluid transport fromthe subducted lithosphere to the mantle wedge. Continental crustal material seems to play a significant roleparticularly in the Macolod Corridor and the Mindoro segment,based on the high LILE, La/Sm ratios, radiogenic Sr isotopes,and ¹⁸O values. The Macolod Corridor and the Mindoro segmenthave undergone source contamination by crustal material fromthe North Palawan-Mindoro crustal block either during the collisionof this block with the Manila Trench or by subduction of sedimentsrich in this crustal material. A similar component has alsobeen detected in the Bataan segment but in minor amounts. Thetrace element and isotopic differences between the northernand southern sections of the arc are interpreted in terms ofvariable composition (i. e., variable amounts of a crustal componentintroduced from the Palawan-Mindoro crustal terrain) of themetasomatic fluids released into the source.     

The effect of bulk rock composition on the stability of amphibole in the upper mantle: Implications for solidus positions and mantle metasomatism

Summary The stability of pargasitic amphibole in the upper mantle is a function of water content and bulk rock composition, and under water-undersaturated conditions, the stability of amphibole controls the solidus position. Experiments in the system Tinaquillo peridotite +0.2% H₂O, a refractory peridotite under water-undersaturated conditions, show that amphibole is stable to 1030°C and 26 kb. In contrast, pargasitic amphibole is stable to 1150°C and 30 kb in Hawaiian pyrolite, a more fertile peridotite composition. This indicates that under water-undersaturated conditions, the most fertile part of a crystallizing mantle diapir with an inhomogeneous composition will solidify first while a more refractory component will contain an alkali-rich melt which will have the ability to metasomatize adjacent regions. The relative stabilities of amphibole in refractory and fertile bulk compositions may result in increasing rather than diminishing chemical contrasts in high temperature lherzolite, i.e. a process of metamorphic differentiation. Ti, Fe, Al and Na metasomatism can therefore be considered a normal occurrence associated with the upward migration and solidification of an H₂O-bearing mantle diapir.

---

Der Einfluß der Gesamtgesteins-Zusammensetzung auf die Stabilität von Amphibol im oberen Mantel: Bedeutung für Solidus-Positionen und Mantel-Metasomatose

Zusammenfassung Die Stabilität pargasitischer Amphibole im oberen Mantel ist eine Funktion von Wassergehalt und Gesamtgesteins-Zusammensetzung. Unter wasser-untersättigten Bedingungen, kontrolliert die Stabilität von Amphibol die Solidus-Position. Experimente in dem System Tinaquillo Peridotit +0,2% H₂O, einem refraktären Peridotit unter wasser-untersättigten Bedingungen, zeigen daß Amphibol bis 1030°C und 26 Kb stabil ist. Im Gegensatz dazu ist pargasitische Hornblende in einem Hawaii-Pyrolit, von mehr fertiler Peridotit-Zuammensetzung, bis 1150°C und 30 Kb stabil. Das zeigt, daß bei wasser-untersättigten Bedingungen der am meisten produktive Teil eines kristallisierenden Mantel-Diapirs mit inhomogener Zusammensetzung sich zuerst verfestigen wird, während eine mehr refraktäre Komponente eine alkali-reiche Schmelze enthalten wird, die wiederum die Fähigkeit hat, umliegende Bereiche metasomatisch zu beeinflussen. Die relativen Stabilitäten von Amphibol in refraktären und fertilen Gesamtzusammensetzungen können dazu führen, daß die chemischen Gegensätze in Hochtemperaturlherzoliten eher zunehmen als abnehmen, d. h. ein Prozeß metamorpher Differentiation. Ti, Fe, Al und Na Metasomatose können deshalb als ein verbreiteter Vorgang, der mit der Aufwärtsbewegung und Verfestigung eines H₂O-führenden Mantel-Diapirs assoziiert ist, betrachtet werden.

---

---

With 4 Figures     

An internal variable model for the creep of rocksalt

Summary The creep strain rate of rocksalt, like that of other ductile crystalline materials, can be described by a power law equation of the type ( ) ⁿ , where the active stress is the difference between the total deviatoric applied stress and an internal stress _i . In this paper, the origin and the nature of this internal stress, which develops during inelastic deformation of the material, are discussed. It is shown that this internal stress can serve as an internal (or state) variable in the constitutive model of rocksalt, which reflects the microstructure evolution of the material under the competitive action of hardening and recovery mechanisms.An analysis of experimental data, both our own and those taken from the literature, demonstrates that such a law is able to correctly reproduce rocksalt creep test results in the steady-state domain. The proposed model is in accordance with the macroscopic and microscopic behavior of salt, and with direct measurements of the internal stresses made by others on this material.