Tropical Indian Ocean variability revealed by self-organizing maps

The tropical Indian Ocean climate variability is investigated using an artificial neural network analysis called self-organizing map (SOM) for both observational data and coupled model outputs. The SOM successfully captures the dipole sea surface temperature anomaly (SSTA) pattern associated with the Indian Ocean Dipole (IOD) and basin-wide warming/cooling associated with ENSO. The dipole SSTA pattern appears only in boreal summer and fall, whereas the basin-wide warming/cooling appears mostly in boreal winter and spring owing to the phase-locking nature of these phenomena. Their occurrence also undergoes significant decadal variation. Composite diagrams constructed for nodes in the SOM array based on the simulated SSTA reveal interesting features. For the nodes with the basin-wide warming, a strong positive SSTA in the eastern equatorial Pacific, a negative Southern Oscillation, and a negative precipitation anomaly in East Africa are found. The nodes with the positive IOD are associated with a weak positive SSTA in the central equatorial Pacific or positive SSTA in the eastern equatorial Pacific, a positive (negative) sea level pressure anomaly in the eastern (western) tropical Indian Ocean, and a positive precipitation anomaly over East Africa. The warming in the central equatorial Pacific appears to correspond to El Niño Modoki discussed recently. These results suggest usefulness of SOM in studying large-scale ocean–atmosphere coupled phenomena.     

新疆西准噶尔地区塔斯特岩体锆石LA-ICPMS年龄及其意义

西准噶尔萨吾尔地区地处新疆阿勒泰地区吉木乃县及塔城地区和丰县。酸性侵入岩在区内较发育,其研究对于准噶尔地区后碰撞构造-岩浆演化时限的厘定具有重要的意义。塔斯特岩体是萨吾尔地区重要的复式花岗岩体,与区内铜、金矿化关系密切。塔斯特岩体锆石 LA—ICPMS 年龄为337±4Ma(1σ),形成于早石炭世。结合已有的地球化学特征研究,这一年龄显示西准噶尔地区在早石炭世已处于后碰撞环境。本文获得的塔斯特岩体成岩年龄是目前在准噶尔地区获得的最老的后碰撞花岗岩锆石 U—Pb 年龄之一,为准噶尔地区后碰撞构造-岩浆演化时限的厘定提供了新的证据。     

Geological setting and timing of the Chah Zard breccia-hosted epithermal gold–silver deposit in the Tethyan belt of Iran

The breccia-hosted epithermal gold–silver deposit of Chah Zard is located within a high-K, calc-alkaline andesitic to rhyolitic volcanic complex in the central part of the Urumieh-Dokhtar Magmatic Arc (UDMA), west central Iran. The total measured resource for Chah Zard is ∼2.5 million tonnes of ore at 12.7 g/t Ag and 1.7 g/t Au (28.6 t Ag, 3.8 t Au), making it one of the largest epithermal gold deposits in Iran. Magmatic and hydrothermal activity was associated with local extensional tectonics in a strike-slip regime formed in transtensional structures of the Dehshir-Baft strike-slip fault system. The host rocks of the volcanic complex consist of Eocene sedimentary and volcanic rocks covered by Miocene sedimentary rocks. LA-ICP–MS U–Pb zircon geochronology yields a mean age of 6.2 ± 0.2 Ma for magmatic activity at Chah Zard. This age represents the maximum age of mineralization and may indicate a previously unrecognized mineralization event in the UDMA. Breccias and veins formed during and after the waning stages of explosive brecciation events due to shallow emplacement of rhyolite porphyry. Detailed systematic mapping leads to the recognition of three distinct breccia bodies: volcaniclastic breccia with a dominantly clastic matrix; gray polymict breccia with a greater proportion of hydrothermal cement; and mixed monomict to polymict breccia with clay matrix. The polymictic breccias generated bulk-mineable ore, whereas the volcaniclastic breccia is relatively impermeable and largely barren. Precious metals occur with sulfide and sulfosalt minerals as disseminations, as well as in the veins and breccia cements. There is a progression from pyrite-dominated (stage 1) to pyrite-base metal sulfide and sulfosalt-dominated (stages 2 and 3) to base metal sulfide-dominated (stage 4) breccias and veins. Hydrothermal alteration and deposition of gangue minerals progressed from illite-quartz to quartz-adularia, carbonate, and finally gypsum-dominated assemblages. Free gold occurs in stages 2 and 4, principally intergrown with pyrite, quartz, chalcopyrite, galena, sphalerite, and Ag-rich tennantite–tetrahedrite, and also as inclusions in pyrite. High Rb/Sr ratios in ore-grade zones are closely related to sericite and adularia alteration. Positive correlations of Au and Ag with Cu, As, Pb, Zn, Sb, and Cd in epithermal veins and breccias suggest that all these elements are related to the same mineralization event.     

Water structure at hematite-water interfaces

The atomic-level structure of water at mineral surfaces is an important controlling factor in interfacial reactions such as foreign ion incorporation, crystal growth and dissolution, and redox reactions. Molecular dynamics simulations with four different models based on interatomic potentials have been carried out to determine the atomic-level structure of three hematite-water interfaces. In addition, for each of the three surfaces, different terminations or protonation schemes were considered. The availability of surface X-ray scattering data for the surfaces considered here allowed for an extensive comparison with experimental data. Qualitatively, with the exception of one termination with one model, all models predict the correct arrangement of water molecules at the interface. Quantitatively, the agreement with experimental positions, distances, and layer occupancies is good to excellent, especially given the range of values reported in published experimental studies. Therefore, this study provides further evidence that interatomic potential models can be used to reliably predict the structure of mineral-water interfaces. In addition, molecular simulations are a valuable source of information to complement surface X-ray scattering experiments owing to their ability to directly determine the position of hydrogen atoms and to yield three-dimensional predicted structures at no added cost, as demonstrated in this work. Indeed, the molecular dynamics trajectories were analyzed to determine the surface structural controls on the interfacial water structure. Each of the three surface functional groups present at the surfaces considered in this work, namely, triply-coordinated oxo, doubly-coordinated hydroxo, and singly-coordinated aquo groups, was found to form similar hydrogen bond configurations with adsorbed water molecules at all surfaces. Oxo groups accept long-lasting and linear hydrogen bonds from adsorbed water molecules; hydroxo groups can form hydrogen bonds with other surface functional groups as well as with adsorbed water molecules; and aquo groups normally only donate hydrogen bonds to other surface groups or adsorbed water molecules. Additionally, the majority of adsorbed water molecules were found to adopt multiple configurations and orientations. This information was used to evaluate three-dimensional structural models of the interfaces, which were previously derived experimentally from one-dimensional electron density profiles and steric considerations.     

Key features of the IPSL ocean atmosphere model and its sensitivity to atmospheric resolution

This paper presents the major characteristics of the Institut Pierre Simon Laplace (IPSL) coupled ocean–atmosphere general circulation model. The model components and the coupling methodology are described, as well as the main characteristics of the climatology and interannual variability. The model results of the standard version used for IPCC climate projections, and for intercomparison projects like the Paleoclimate Modeling Intercomparison Project (PMIP 2) are compared to those with a higher resolution in the atmosphere. A focus on the North Atlantic and on the tropics is used to address the impact of the atmosphere resolution on processes and feedbacks. In the North Atlantic, the resolution change leads to an improved representation of the storm-tracks and the North Atlantic oscillation. The better representation of the wind structure increases the northward salt transports, the deep-water formation and the Atlantic meridional overturning circulation. In the tropics, the ocean–atmosphere dynamical coupling, or Bjerknes feedback, improves with the resolution. The amplitude of ENSO (El Niño-Southern oscillation) consequently increases, as the damping processes are left unchanged.     

Molecular simulation of the diffusion of uranyl carbonate species in aqueous solution

Potential-based molecular dynamics simulations of aqueous uranyl carbonate species (M_xUO₂(CO₃)_y^2+2x−2y with M = Mg, Ca, or Sr) were carried out to gain molecular-level insight into the hydration properties of these species. The simulation results were used to estimate the self-diffusion coefficients of these uranyl carbonate species, which often dominate uranyl speciation in groundwater systems. The diffusion coefficients obtained for the monoatomic alkaline-earth cations and polyatomic ions (uranyl, carbonate, and uranyl tri-carbonate) were compared with those calculated from the Stokes-Einstein (SE) equation and its variant formulation by Impey et al. (1983). Our results show that the equation of Impey et al. (1983), originally formulated for monovalent monoatomic ions, can be extended to divalent monoatomic ions, with some success in reproducing the absolute values and the overall trend determined from the molecular dynamics simulations, but not to polyatomic ions, for which the hydration shell is not spherically symmetrical. Despite the quantitative failure of both SE formulations, a plot of the diffusion coefficients of the uranyl carbonate complexes as a function of the inverse of the equivalent spherical radius showed that a general linear dependence is observed for these complexes as expected from the SE equation. The nature of the alkaline-earth cation in the uranyl carbonate complexes was not found to have a significant effect on the ion’s diffusion coefficient, which suggests that the use of a single diffusion coefficient for different alkaline-earth uranyl carbonate complexes in microscopic diffusion models is appropriate.The potential model reproduced well published quantum mechanical and experimental data of and of the individual constituent ions, and therefore is expected to offer reliable predictions of the structure of magnesium and strontium uranyl carbonate aqueous species, for which there is no structural data available to date. In addition, the interatomic distances reported for could help with the refinement of the interpretation of EXAFS data of these species, which is made difficult by the similar uranium-distant carbonate oxygen and uranium-calcium distances.An analysis of the dynamics of water exchange around the alkaline-earth cations revealed that the presence of the uranyl tri-carbonate molecule has a strong influence on the geometry of the cation’s first hydration shell, which, in turn, can considerably affect the water exchange kinetics depending on whether the imposed geometry matches that around the isolated alkaline-earth cation. This result shows that the alkaline-earth uranyl carbonate complexes have distinct water exchange dynamics, which may lead to different reactivities. Finally, significant changes in water residence time were also predicted when replacing carbonate for water ligands in the uranyl coordination shell.     

Late Holocene high precipitation events recorded in lake sediments and catchment geomorphology,Lake Vuoksjávrátje,NW Sweden

In this paper, we show the potential of combining multi‐proxy analysis of lake sediments with catchment geomorphology to better understand palaeoenvironmental changes. Previous studies have suggested that alpine lakes in N Sweden may be highly sensitive to variations in catchment erosion and precipitation, and that this sensitivity may influence ecologically based reconstructions of past temperature changes. We analysed lake sediments covering the last 5100 years from the alpine Lake Vuoksjávrátje in NW Sweden in order to identify different erosional regimes in the lake catchment, sediment sources and lake sedimentary processes, which ultimately affect the palaeoecological record. The measured proxies include elemental geochemistry from XRF core scanning, grain size, sediment accumulation rates, fraction of terrestrial organic carbon and geomorphological mapping, supported by previously published chironomid and total organic carbon data. From the integrated results we identified time intervals when increased intensity of precipitation altered sedimentation and lake catchment erosional processes. The most prominent event occurred c. 2900 cal. a BP and is interpreted to be the result of excessive precipitation, possibly related to the climatic shift towards cooler and wetter conditions referred to as the 2.8 ka event.     

Interpreting the Helioseismic and Magnetic Imager (HMI) Multi-Height Velocity Measurements

The Solar Dynamics Observatory/Helioseismic and Magnetic Imager (SDO/HMI) filtergrams, taken at six wavelengths around the Fe i 6173.3 Å line, contain information about the line-of-sight velocity over a range of heights in the solar atmosphere. Multi-height velocity inferences from these observations can be exploited to study wave motions and energy transport in the atmosphere. Using realistic convection-simulation datasets provided by the STAGGER and MURaM codes, we generate synthetic filtergrams and explore several methods for estimating Dopplergrams. We investigate at which height each synthetic Dopplergram correlates most strongly with the vertical velocity in the model atmospheres. On the basis of the investigation, we propose two Dopplergrams other than the standard HMI-algorithm Dopplergram produced from HMI filtergrams: a line-center Dopplergram and an average-wing Dopplergram. These two Dopplergrams correlate most strongly with vertical velocities at the heights of 30?–?40 km above (line center) and 30?–?40 km below (average wing) the effective height of the HMI-algorithm Dopplergram. Therefore, we can obtain velocity information from two layers separated by about a half of a scale height in the atmosphere, at best. The phase shifts between these multi-height Dopplergrams from observational data as well as those from the simulated data are also consistent with the height-difference estimates in the frequency range above the photospheric acoustic-cutoff frequency.