The WRF model performance for the simulation of heavy precipitating events over Ahmedabad during August 2006

The summer monsoon season of the year 2006 was highlighted by an unprecedented number of monsoon lows over the central and the western parts of India, particularly giving widespread rainfall over Gujarat and Rajasthan. Ahmedabad had received 540.2mm of rainfall in the month of August 2006 against the climatological mean of 219.8mm. The two spells of very heavy rainfall of 108.4mm and 97.7mm were recorded on 8 and 12 August 2006 respectively. Due to meteorological complexities involved in replicating the rainfall occurrences over a region, the Weather Research and Forecast (WRF-ARW version) modeling system with two different cumulus schemes in a nested configuration is chosen for simulating these events. The spatial distributions of large-scale circulation and moisture fields have been simulated reasonably well in this model, though there are some spatial biases in the simulated rainfall pattern. The rainfall amount over Ahmedabad has been underestimated by both the cumulus parameterization schemes. The quantitative validation of the simulated rainfall is done by calculating the categorical skill scores like frequency bias, threat scores (TS) and equitable threat scores (ETS). In this case the KF scheme has outperformed the GD scheme for the low precipitation threshold.     

Application of precise 142Nd/144Nd analysis of small samples to inclusions in diamonds (Finsch,South Africa) and Hadean Zircons (Jack Hills,Western Australia)

¹⁴⁶Sm–¹⁴²Nd and ¹⁴⁷Sm–¹⁴³Nd systematics were investigated in garnet inclusions in diamonds from Finsch (S. Africa) and Hadean zircons from Jack Hills (W. Australia) to assess the potential of these systems as recorders of early Earth evolution. The study of Finsch inclusions was conducted on a composite sample of 50 peridotitic pyropes with a Nd model age of 3.3 Ga. Analysis of the Jack Hills zircons was performed on 790 grains with ion microprobe ²⁰⁷Pb/²⁰⁶Pb spot ages from 3.95 to 4.19 Ga. Finsch pyropes yield 100 × ?¹⁴²Nd = ? 6 ± 12 ppm, ?¹⁴³Nd = ? 32.5, and ¹⁴⁷Sm/¹⁴⁴Nd = 0.1150. These results do not confirm previous claims for a 30 ppm ¹⁴²Nd excess in South African cratonic mantle. The lack of a ¹⁴²Nd anomaly in these inclusions suggests that isotopic heterogeneities created by early mantle differentiation were remixed at a very fine scale prior to isolation of the South African lithosphere. Alternatively, this result may indicate that only a fraction of the mantle experienced depletion during the first 400 Myr of its history. Analysis of the Jack Hills zircon composite yielded 100 × ?¹⁴²Nd = 8 ± 10 ppm, ?¹⁴³Nd = 45 ± 1, and ¹⁴⁷Sm/¹⁴⁴Nd = 0.5891. Back-calculation of this present-day ?¹⁴³Nd yields an unrealistic estimate for the initial ?¹⁴³Nd of ? 160 ?-units, clearly indicating post-crystallization disturbance of the ¹⁴⁷Sm–¹⁴³Nd system. Examination of ^146,147Sm–^142,143Nd data reveals that the Nd budget of the Jack Hills sample is dominated by non-radiogenic Nd, possibly contained in recrystallized zircon rims or secondary subsurface minerals. This secondary material is characterized by highly discordant U–Pb ages. Although the mass fraction of altered zircon is unlikely to exceed 5–10% of total sample, its high LREE content precludes a reliable evaluation of ¹⁴⁶Sm–¹⁴²Nd systematics in Jack Hills zircons.     

Controls on porphyroblast size along a regional metamorphic field gradient

Garnet-bearing schists from the Waterville Formation of south-central Maine provide an opportunity to examine the factors governing porphyroblast size over a range of metamorphic grade. Three-dimensional sizes and locations for all garnet porphyroblasts were determined for three samples along the metamorphic field gradient spanning lowest garnet through sillimanite grade, using high-resolution X-ray computed tomography. Comparison of crystal size distributions to previous data sets obtained by stereological methods for the same samples reveals significant differences in mode, mean, and shape of the distributions. Quantitative textural analysis shows that the garnets in each rock crystallized in a diffusion-controlled nucleation and growth regime. In contrast to the typical observation of a correlation between porphyroblast size and position along a metamorphic field gradient, porphyroblast size of the lowest-grade specimen is intermediate between the high- and middle-grade specimens’ sizes. Mean porphyroblast size does not correlate with peak temperatures from garnet-biotite Fe-Mg exchange thermometry, nor is post-crystallization annealing (Ostwald Ripening) required to produce the observed textures, as was previously proposed for these rocks. Robust pseudosection calculations fail to reproduce the observed garnet core compositions for two specimens, suggesting that these calc-pelites experienced metasomatism. For each of these two specimens, Monte Carlo calculations suggest potential pre-metasomatism bulk compositions that replicate garnet core compositions. Pseudosection analyses allow the estimation of the critical temperatures for garnet growth: ∼481, ∼477, and ∼485°C for the lowest-garnet-zone, middle-garnet-zone, and sillimanite-zone specimens, respectively. Porphyroblast size appears to be determined in this case by a combination of the heating rate during garnet crystallization, the critical temperature for the garnet-forming reaction and the kinetics of nucleation. Numerical simulations of thermally accelerated, diffusion-controlled nucleation, and growth for the three samples closely match measured crystal size distributions. These observations and simulations suggest that previous hypotheses linking the garnet size primarily to the temperature at the onset of porphyroblast nucleation can only partially explain the observed textures. Also important in determining porphyroblast size are the heating rate and the distribution of favorable nucleation sites.     

The iron-isotope fractionation dictated by the carboxylic functional: An ab-initio investigation

The ground-state geometries, electronic energies and vibrational properties of carboxylic complexes of iron were investigated both in vacuo and under the effect of a reaction field, to determine thermodynamic properties of iron-acetates and the role of the carboxylic functional on the isotopic imprinting of this metal in metalorganic complexation. The electronic energy, zero point corrections and thermal corrections of these substances at variational state were investigated at the DFT/B3LYP level of theory with different basis set expansions and the effect of the reaction field on the variational structures was investigated through the Polarized Continuun Model. Thermochemical cycle calculations, combined with solvation energy calculations and appropriate scaling from absolute to conventional properties allowed to compute the Gibbs free energy of formation from the elements of the investigated aqueous species and to select the best procedure to be applied in the successive vibrational analysis. The best compliance with the few existing thermodynamic data for these substances was obtained by coupling the gas phase calculations at DFT/B3LYP level with the [6-31G(d,p)]-[6-31G+(d,p)] (for cations and neutral molecules - anions; respectively) with solvation calculations adopting atomic radii optimized for the HF/6-31G(d) level of theory (UAHF). A vibrational analysis conducted on ⁵⁴Fe, ⁵⁶Fe, ⁵⁷Fe and ⁵⁸Fe gaseous isotopomers yielded reduced partition function ratios which increased not only with the nominal valence of the central cation, as expected, but, more importantly, with the extent of the complexation operated by the organic functional. Coupling thermodynamic data with separative effects it was shown that this last is controlled, as expected, by the relative bond strength of the complex in both aggregation states. Through the Integral Equation Formalism of the Polarized Continuum Model (IEFPCM) the effect of the ionic strength of the solution and of a T-dependent permittivity on the energy and separative effects of the solvated metalorganic complexes were analyzed in detail. The solvent effect in the standard state (hypothetical one-molal solution referred to infinite dilution; T = 298.15 K, P = 1 bar) is a limited reduction of the separative effects of all the isotopomeric couples. With an increase in T (and the concomitant decrease in the dielectric constant of the solvent) this effect diminishes progressively.     

Thermodynamic model for arsenic speciation in sulfidic waters: A novel use of ab initio computations

Recent discoveries demonstrate that the chemistry of arsenic in sulfidic waters is much more complex that previously believed. One implication is that all earlier thermodynamic data on stabilities of As thioanions require revision. Previously used experimental approaches for determining As thioanion stabilities may be inadequate to deal with the full range of complexity. Here we use computational as well as empirical information to construct a provisional model for equilibrium As thioanion distributions in sulfidic waters. Whereas previous authors have argued for either As(III) or As(V) thioanions, the new model predicts that both are important and can occur simultaneously under commonly encountered pH and ΣS^−II conditions. At the order of magnitude level, the model reasonably predicts the solubility of As₂S₃ in sulfidic solutions, provides tentative peak assignments for published Raman spectroscopic data and plausibly accounts for how sulfide modifies the bacterial toxicity of As. The model yields a thermodynamic justification for how sulfide, which is usually regarded as a reducing agent, can counter-intuitively drive oxidation of As(III) to As(V), as has been observed both in the laboratory and in the field. Despite its uncertain accuracy, the model serves as a useful source of new, testable hypotheses about As geochemistry and highlights crucial experimental data needs.     

Synchronism in compression of the lithosphere in the central Atlantic and Western Tethys at the Tortonian-Messinian transition

The structure of anomalously uplifted areas in transverse ridges of the Vema, S o Paulo, and Romanche fracture zones is considered. It is concluded that their formation and eventual development in the present-day structure of the central Atlantic bottom proceeded during two stages. The first stage that corresponds to a short period at the Tortonian-Messinian transition (10 Ma ago) was marked by transportation of deep-seated rocks into the upper part of the lithosphere along thrust faults with mass motion in the meridional direction along the axis of the Mid-Atlantic Ridge. The second stage was characterized by contrasting highamplitude vertical movements from 10 to 3 Ma ago. It is suggested that near-meridional compression in the domains surrounding the Western Tethys in the Tortonian-Messinian resulted in deformation of the upper lithosphere within large transform fracture zones of the central Atlantic. The deformation that occurred 10 Ma ago was a manifestation of the global neotectonic epoch of the Earth.     

Large and giant hydrocarbon accumulations in the transitional continent-ocean zone

The petroleum resource potential is considered for the Atlantic, West Pacific, and East Pacific types of deepwater continental margins. The most considerable energy resources are concentrated at the Atlantic-type passive margins in the zone transitional to the ocean. The less studied continental slope of backarc seas of the generally active margins of the West Pacific type is currently not so rich in discoveries as the Atlantic-type margin, but is not devoid of certain expectations. In some of their parameters, the margins bounded by continental slopes may be regarded as analogs of classical passive margins. At the margins of the East Pacific type, the petroleum potential is solely confined to transform segments. In the shelf-continental-slope basins of the rift and pull-apart nature, petroleum fields occur largely in the upper fan complex, and to a lesser extent in the lower graben (rift) complex. In light of world experience, the shelf-continental-slope basins of the Arctic and Pacific margins of Russia are evaluated as highly promising.     

Biofacies analysis of Upper Cretaceous deposits in the Ust-Yenisei region: Implications of palynomorphs

The results of palynomorph biofacies analysis in the Upper Cretaceous deposits of the Ust-Yenisei region are presented. The established facies confinement and indicative features of separate palynomorph groups are used, along with identified dinocyst morphotypes and taxa, in paleogeographic reconstructions. Seven palynomorph associations characterizing continental, coastal-marine, shallow-and deep-water facies are distinguished based on quantitative proportions between morphological groupings and individual taxa. As boundaries between distinguishable biostratigraphic and facies subdivisions do not coincide, dinocysts were likely insignificantly dependent in distribution on facies in the West Siberian epicontinental basin at least. On the other hand, distribution trends of particular dinocyst morphotypes and other microphytofossils are correlative with transgressive-regressive cycles and can be used for reconstruction of paleoenvironments.     

An Attempt to Predict the Failure Time of Abandoned Mine Pillars

Summary  Left in place pillars of abandoned mines are subject to weathering (e.g action of water, bacteria) that degrades their mechanical strength and eventually leads to collapse. A simple weathering model is proposed, that is governed by two parameters: the rate of progression of weathered front and the rate of degradation of the compressive strength with time. Both plane strain and axisymmetric analyses are performed and closed form solutions of the variation with time of the bearing capacity of the pillar are given. Experimental data of the tests conducted on gypsum and anhydrite specimens attacked by water are presented. It is shown that in order to fit the experimental data a third parameter must be introduced. New closed form solutions are given and the data are used for estimating the time to failure of abandoned gypsum mines in Northern Italy. Authors’ address: Dr. Riccardo Castellanza, Research Assistant, Department of Structural Engineering, Milan University of Technology (Politecnico), Piazza Leonardo da Vinci 32, 20133 Milan, Italy