Reinterpreting the thermocline feedback in the western-central equatorial Pacific and its relationship with the ENSO modulation

Vertical stratification changes at low frequency over the last decades are the largest in the western-central Pacific and have the potential to modify the balance between ENSO feedback processes. Here we show evidence of an increase in thermocline feedback in the western-central equatorial Pacific over the last 50 years, and in particular after the climate shift of 1976. It is demonstrated that the thermocline feedback becomes more effective due to the increased stratification in the vicinity of the mean thermocline. This leads to an increase in vertical advection variability twice as large as the increase resulting from the stronger ENSO amplitude (positive asymmetry) in the eastern Pacific that connects to the thermocline in the western-central Pacific through the basin-scale ‘tilt’ mode. Although the zonal advective feedback is dominant over the western-central equatorial Pacific, the more effective thermocline feedback allows for counteracting its warming (cooling) effect during warm (cold) events, leading to the reduced covariability between SST and thermocline depth anomalies in the NINO4 (160°E–150°W; 5°S–5°N) region after the 1976 climate shift. This counter-intuitive relationship between thermocline feedback strength as derived from the linear relationship between SST and thermocline fluctuations and stratification changes is also investigated in a long-term general circulation coupled model simulation. It is suggested that an increase in ENSO amplitude may lead to the decoupling between eastern and central equatorial Pacific sea surface temperature anomalies through its effect on stratification and thermocline feedback in the central-western Pacific.     

LiDAR-based digital terrain analysis of an area exposed to the risk of lava flow invasion: the Zafferana Etnea territory, Mt. Etna (Italy)

The town of Zafferana Etnea, located on the southeastern slope of Mt. Etna volcano (Italy), has been repeatedly threatened by lava flows in recent centuries. The last serious threat occurred during the 1991–1993 eruption, when the lava front came to a halt only 1.7 km from the centre of town. Morphostructural data derived from light detection and ranging (LiDAR) surveys carried out on Etna in 2005 have enabled us to evaluate the risk of lava invasion in a section (16 km²) of the Zafferana Etnea territory. Qualitative and quantitative results are obtained combining the information derived from LiDAR analysis with geological, morphological and structural data using geographic information systems technology (GIS). The study quantifies in unprecedented detail the areal extent and volume of forested and urban areas and its degree of exposure to different levels of hazard from future lava invasion. Nearly 52% of the urban texture fall into areas of moderate to high risk from lava invasion. Future land use planning should take these findings into account and promote new development preferentially in areas of lower risk.     

Effect of moderate shock waves on magnetic susceptibility and microstructure of a magnetite‐bearing ore

This study demonstrates a relationship between changes of magnetic susceptibility and microstructure developing in minerals of a magnetite‐bearing ore, experimentally shocked to pressures of 5, 10, 20, and 30 GPa. Shock‐induced effects on magnetic properties were quantified by bulk magnetic susceptibility measurements while shock‐induced microstructures were studied by high‐resolution scanning electron microscopy. Microstructural changes were compared between magnetite, quartz, amphibole, and biotite grains. In the 5 GPa sample, a sharp drop of magnetic susceptibility correlates with distinct fragmentation as well as with formation of shear bands and twins in magnetite. At 10 GPa, shear bands and twins in magnetite are accompanied by droplet‐shaped nanograins. In this shock pressure regime, quartz and amphibole still show intensive grain fragmentation. Twins in quartz and foam‐shaped, highly porous amphibole are formed at 20 and 30 GPa. The formation of porous minerals suggests that shock heating of these mineral grains resulted in localized temperature spikes. The identified shock‐induced features in magnetite strongly advise that variations in the bulk magnetic susceptibility result from cooperative grain fragmentation, plastic deformation and/or localized amorphization, and probably postshock annealing. In particular, the increasing shock heating at high pressures is assumed to be responsible for a partial defect annealing which we suggest to be responsible for the almost constant values of magnetic susceptibility above 10 GPa.     

Experimental study of aluminum speciation in fluoride-rich supercritical fluids

The solubility of the albite-paragonite-quartz mineral assemblage was measured as a function of NaCl and fluorine concentration at 400°C, 500 bars and at 450°C, 500 and 1000 bars. Decreasing Al concentrations with increasing NaCl molality in F-free fluids of low salinity (mNaCl < 0.01) demonstrates that Al(OH)₄⁻ dominates Al speciation and is formed according to the reaction 0.5 NaAl₃Si₃O₁₂H_2(cr)+2 H₂O = 0.5 NaAlSi₃O_8(cr)+Al(OH)₄⁻+H⁺. Log K results for this reaction are −11.28 ± 0.10 and −10.59 ± 0.10 at 400°C, 500 bars and 450°C, 1000 bars, respectively. Upon further salinity increase, Al concentration becomes constant (at 400°C, 500 bars) or even rises (at 450°C, 1000 bars). The observed Al behavior can be explained by the formation of NaAl(OH)₄⁰_(aq) or NaAl(OH)₃Cl_(aq)⁰. The calculated constant for the reaction Al(OH)₄⁻+Na⁺=NaAl(OH)₄⁰_(aq) expressed in log units is equal to 2.46 and 2.04 at 400°C, 500 bars and 450°C, 1000 bars, respectively. These values are in good agreement with the predictions given in Diakonov et al. (1996). Addition of fluoride at m(NaCl) = const = 0.5 caused a sharp increase in Al concentration in equilibrium with the albite-paragonite-quartz mineral assemblage. As fluid pH was also constant, this solubility increase indicates strong aluminum-fluoride complexation with the formation of NaAl(OH)₃F_(aq)⁰ and NaAl(OH)₂F₂⁰_(aq), according to 0.5 NaAl₃Si₃O₁₂H_2(cr)+Na⁺+HF_(aq)⁰+H₂O = 0.5 NaAlSi₃O_8(cr)+ NaAl(OH)₃F_(aq)⁰+H⁺, log K = −5.17 and −5.23 at 400°C and 450°C, 500 bars, respectively, and 0.5 NaAl₃Si₃O₁₂H_2(cr)+Na⁺+2 HF_(aq)⁰ = 0.5 NaAlSi₃O_8(cr)+NaAl(OH)₂F₂⁰_(aq)+H⁺, log K = −2.19 and −1.64 at the same P-T conditions. It was found that temperature increase and pressure decrease promote the formation of Na-Al-OH-F species. Stability of NaAl(OH)₂F₂⁰_(aq) in low-density fluids also increases relative to NaAl(OH)₃F_(aq)⁰. These complexes, together with Al(OH)₂F_(aq)⁰ and AlOHF₂⁰_(aq), whose stability constants were calculated from the corundum solubility measured by Soboleva and Zaraisky (1990) and Zaraisky (1994), are likely to dominate Al speciation in metamorphic fluids containing several ppm of fluorine.     

On Motions Near the Lagrange Equilibrium Point L 4 in the Case of Routh's Critical Mass Ratio

We deal with the planar restricted circular problem of three bodies. We study trajectories in a small neighborhood of the Lagrange equilibrium point L ₄ when mass ratio is close to Routh's value. In particular, we show that the case of proper degeneracy takes place and for most initial conditions trajectories are conditionally-periodic. We obtain an approximate representation of families of periodic solution emanating from the equilibrium point L ₄. We also show that in the case of instability of L ₄ the trajectories starting in a vicinity of L ₄ remain in a finite domain forever. We give an upper bound of this domain. To carry out our investigation, we analyze the dynamics of a general Hamiltonian system with two degrees of freedom in the case of 1 : 1 resonance in detail.This revised version was published online in October 2005 with corrections to the Cover Date.     

Finite deformation analysis of geomaterials

The mathematical structure and numerical analysis of classical small deformation elasto–plasticity is generally well established. However, development of large deformation elastic–plastic numerical formulation for dilatant, pressure sensitive material models is still a research area. In this paper we present development of the finite element formulation and implementation for large deformation, elastic–plastic analysis of geomaterials. Our developments are based on the multiplicative decomposition of the deformation gradient into elastic and plastic parts. A consistent linearization of the right deformation tensor together with the Newton method at the constitutive and global levels leads toward an efficient and robust numerical algorithm. The presented numerical formulation is capable of accurately modelling dilatant, pressure sensitive isotropic and anisotropic geomaterials subjected to large deformations. In particular, the formulation is capable of simulating the behaviour of geomaterials in which eigentriads of stress and strain do not coincide during the loading process. The algorithm is tested in conjunction with the novel hyperelasto–plastic model termed the B material model, which is a single surface (single yield surface, affine single ultimate surface and affine single potential surface) model for dilatant, pressure sensitive, hardening and softening geomaterials. It is specifically developed to model large deformation hyperelasto–plastic problems in geomechanics. We present an application of this formulation to numerical analysis of low confinement tests on cohesionless granular soil specimens recently performed in a SPACEHAB module aboard the Space Shuttle during the STS‐89 mission. We compare numerical modelling with test results and show the significance of added confinement by the thin hyperelastic latex membrane undergoing large stretching. Copyright © 2001 John Wiley & Sons, Ltd.     

Experimental study of gold-hydrosulphide complexing in aqueous solutions at 350-500°C, 500 and 1000 bars using mineral buffers

The solubility of gold was measured in KCl solutions (0.001-0.1 m) at near-neutral to weakly acidic pH in the presence of the K-feldspar-muscovite-quartz, andalusite-muscovite-quartz, and pyrite-pyrrhotite-magnetite buffers at temperatures 350 to 500°C and pressures 0.5 and 1 kbar. These mineral buffers were used to simultaneously constrain pH, f(S₂), and f(H₂). The experiments were performed using a CORETEST flexible Ti-cell rocking hydrothermal reactor enabling solution sampling at experimental conditions. Measured log m(Au) (mol/kg H₂O) ranges from −7.5 at weakly acid pH to −5.9 in near-neutral solutions, and increases slightly with temperature. Gold solubility in weakly basic and near-neutral solutions decreases with decreasing pH at all temperatures, which implies that Au(HS)₂ is the dominant Au species in solution. In more acidic solutions, solubility is independent of pH. Comparison of the experimentally measured solubilities with literature values for Au hydrolysis constants demonstrates that at 350°C dominates Au aqueous speciation at the weakly acidic pH and f(S₂)/f(H₂) conditions imposed by the pyrite-pyrrhotite-magnetite buffer. In contrast, at temperatures >400°C becomes less important and predominates in weakly acid solutions. Solubility data collected in this study were used to calculate the following equilibrium reaction constants:

     

Carbonatite diversity in the Central Andes: the Ayopaya alkaline province,Bolivia

The Ayopaya province in the eastern Andes of Bolivia, 100 km NW of Cochabamba, hosts a Cretaceous alkaline rock series within a Palaeozoic sedimentary sequence. The alkaline rock association comprises nepheline-syenitic/foyaitic to ijolitic intrusions, carbonatite, kimberlite, melilititic, nephelinitic to basanitic dykes and diatremes, and a variety of alkaline dykes. The carbonatites display a wide petrographic and geochemical spectrum. The Cerro Sapo area hosts a small calciocarbonatite intrusion and a multitude of ferrocarbonatitic dykes and lenses in association with a nepheline-syenitic stock. The stock is crosscut by a spectacular REE-Sr-Th-rich sodalite-ankerite-baryte dyke system. The nearby Chiaracke complex represents a magnesiocarbonatite intrusion with no evidence for a relationship to igneous silicate rocks. The magnesiocarbonatite ( REE up to 1.3 wt%) shows strong HREE depletion, i.e. unusually high La/Yb ratios (520–1,500). Calciocarbonatites ( REE up to 0.5 wt%) have a flatter REE distribution pattern (La/Yb 95–160) and higher Nb and Zr contents. The sodalite-ankerite-baryte dyke system shows geochemical enrichment features, particularly in Na, Ba, Cl, Sr, REE, which are similar to the unusual natrocarbonatitic lavas of the recent volcano of Oldoinyo Lengai, Tanzania. The Cerro Sapo complex may be regarded as an intrusive equivalent of natrocarbonatitic volcanism, and provides an example for carbonatite genesis by late-stage crystal fractionation and liquid immiscibility. The magnesiocarbonatite intrusion of Chiaracke, on the other hand, appears to result from a primary carbonatitic mantle melt. Deep seated mantle magmatism/metasomatism is also expressed by the occurrence of a kimberlite dyke. Neodymium and strontium isotope data (_Nd 1.4–5.4, ⁸⁷Sr/⁸⁶ Sr<Bulk Earth) indicate a depleted mantle source for the alkaline magmatism. The magmatism of the Ayopaya region is attributed to failed rifting of western South America during the Mesozoic and represents the only occurrence of carbonatite and kimberlite rocks in the Andes.