Alteration patterns of chromian spinels from La Cabaña peridotite,south-central Chile

La Cabaña peridotite is part of a dismembered ophiolite complex located within the metamorphic basement of the Coastal Cordillera of south-central Chile, and is the only location in Chile were Cr-spinels have been described so far. The La Cabaña peridotite is part of the Western Series unit, which comprises meta-sedimentary rocks, metabasites, and serpentinized ultramafic rocks. This unit has been affected by greenschist-facies metamorphism with reported peak PT conditions of 7.0–9.3 kbar and 380°–420 °C. Within La Cabaña peridotite Cr-spinels are present in two localities: Lavanderos and Centinela Bajo. In Lavanderos, Cr-spinel occurs in small chromitite pods and as accessory/disseminated grains with a porous or spongy texture in serpentinite, whereas in Centinela Bajo Cr-spinel is present as accessory zoned grains in partly serpentinized dunites, and in chromitite blocks. All Cr-spinels display variable degrees of alteration to Fe²⁺-rich chromite with a variation trend of major elements from chromite to Fe²⁺-rich chromite similar to those observed in other locations, i.e., an increase in Fe₂O₃ and FeO, a decrease in Al₂O₃ and MgO. Cr₂O₃ content increases from chromite to Fe²⁺-rich chromite in chromitite pods from Lavanderos and chromitite blocks from Centinela Bajo, but decreases in ferrian chromite zones in accessory grains from Centinela Bajo. The minor element (Ti, V, Zn, Ni) content is mostly low and does not exceed 0.4 wt.%, with the exception of MnO (<0.9 wt.%), which shows a correspondence with increasing degree of alteration. Cr# (Cr/Cr?+?Al) versus Mg# (Mg/Mg?+?Fe²⁺) and Fe³⁺/Fe³⁺+Fe²⁺ versus Mg# plots are used to illustrate the Cr-spinel alteration process. Overall, the Cr-spinels from Lavanderos (chromitite pods and disseminated grains) exhibit Cr# values ranging from 0.6 to 1.0, Mg# (Mg/Mg?+?Fe²⁺) below 0.5, and (Fe³⁺/Fe³⁺+Fe²⁺) <0.4. Cr-spinels from chromitites in Centinela Bajo have Cr# and Mg# values that range from 0.65 to 1.0, and 0.7-0.3, respectively, and (Fe³⁺/Fe³⁺+Fe²⁺)?3+/Fe³⁺+Fe²⁺) ratio is less than 0.4 in chromite cores and Fe²⁺-rich chromite, and >0.5 in ferrian chromite and Cr-magnetite. Interpretation of the data obtained and Cr-spinel textures indicate that the alteration of Cr-spinel is a progressive process that involves in its initial stages the reaction of chromite with olivine under water-saturated conditions to produce clinochlore and Fe²⁺-rich chromite. During this stage the chromite can also incorporate Ni, Mn, and/or Zn from the serpentinization fluids. As alteration progresses, Fe²⁺-rich chromite loses mass resulting in the development of a spongy texture. In a later stage and under more oxidizing conditions Fe³⁺ is incorporated in chromite/Fe²⁺-rich chromite shifting its composition to an Fe³⁺-rich chromite (i.e., ferrian chromite). Depending on the fluid/rock and Cr-spinel/silicate ratios, Cr-magnetite can also form over Fe²⁺-rich chromite and/or ferrian chromite as a secondary overgrowth. The compositional changes observed in Cr-spinels from La Cabaña reflect the initial stages of alteration under serpentinization conditions. Results from this study show that the alteration of Cr-spinels is dependent on temperature. The degree and extent of alteration (formation of Fe²⁺-rich and/or ferrian chromite) are controlled by the redox nature of the fluids, the Cr-spinel/silicate and the fluid/rock ratios.     

Mineralogy,geochemistry and petrogenesis of the Khopoli mafic intrusion,Deccan Traps,India

The Khopoli intrusion, exposed at the base of the Thakurvadi Formation of the Deccan Traps in the Western Ghats, India, is composed of olivine gabbro with 50–55 % modal olivine, 20–25 % plagioclase, 10–15 % clinopyroxene, 5–10 % low-Ca pyroxene, and <5 % Fe-Ti oxides. It represents a cumulate rock from which trapped interstitial liquid was almost completely expelled. The Khopoli olivine gabbros have high MgO (23.5–26.9 wt.%), Ni (733–883 ppm) and Cr (1,432–1,048 ppm), and low concentrations of incompatible elements including the rare earth elements (REE). The compositions of the most primitive cumulus olivine and clinopyroxene indicate that the parental magma of the Khopoli intrusion was an evolved basaltic melt (Mg# 49–58). Calculated parental melt compositions in equilibrium with clinopyroxene are moderately enriched in the light REE and show many similarities with Deccan tholeiitic basalts of the Bushe, Khandala and Thakurvadi Formations. Nd-Sr isotopic compositions of Khopoli olivine gabbros (εNd_t?=??9.0 to ?12.7; ⁸⁷Sr/⁸⁶Sr?=?0.7088–0.7285) indicate crustal contamination. AFC modelling suggests that the Khopoli olivine gabbros were derived from a Thakurvadi or Khandala-like basaltic melt with variable degrees of crustal contamination. Unlike the commonly alkalic, pre- and post-volcanic intrusions known in the Deccan Traps, the Khopoli intrusion provides a window to the shallow subvolcanic architecture and magmatic processes associated with the main tholeiitic flood basalt sequence. Measured true density values of the Khopoli olivine gabbros are as high as 3.06 g/cm³, and such high-level olivine-rich intrusions in flood basalt provinces can also explain geophysical observations such as high gravity anomalies and high seismic velocity crustal horizons.     

A reduced-order model for Monte Carlo simulations of stochastic groundwater flow

We explore the ability of the greedy algorithm to serve as an effective tool for the construction of reduced-order models for the solution of fully saturated groundwater flow in the presence of randomly distributed transmissivities. The use of a reduced model is particularly appealing in the context of numerical Monte Carlo (MC) simulations that are typically performed, e.g., within environmental risk assessment protocols. In this context, model order reduction techniques enable one to construct a surrogate model to reduce the computational burden associated with the solution of the partial differential equation governing the evolution of the system. These techniques approximate the model solution with a linear combination of spatially distributed basis functions calculated from a small set of full model simulations. The number and the spatial behavior of these basis functions determine the computational efficiency of the reduced model and the accuracy of the approximated solution. The greedy algorithm provides a deterministic procedure to select the basis functions and build the reduced-order model. Starting from a single basis function, the algorithm enriches the set of basis functions until the largest error between the full and the reduced model solutions is lower than a predefined tolerance. The comparison between the standard MC and the reduced-order approach is performed through a two-dimensional steady-state groundwater flow scenario in the presence of a uniform (in the mean) hydraulic head gradient. The natural logarithm of the aquifer transmissivity is modeled as a second-order stationary Gaussian random field. The accuracy of the reduced basis model is assessed as a function of the correlation scale and variance of the log-transmissivity. We explore the performance of the reduced model in terms of the number of iterations of the greedy algorithm and selected metrics quantifying the discrepancy between the sample distributions of hydraulic heads computed with the full and the reduced model. Our results show that the reduced model is accurate and is highly efficient in the presence of a small variance and/or a large correlation length of the log-transmissivity field. The flow scenarios associated with large variances and small correlation lengths require an increased number of basis functions to accurately describe the collection of the MC solutions, thus reducing significantly the computational advantages associated with the reduced model.     

The caranto paleosol and its role in the early urbanization of Venice

As rising sea level threatens Venice, there is a need to construct a historical framework for interpreting modern environmental changes. Environmental conditions that would later help support Venice's urbanization were established during the Late Glacial period when calcic soils began to develop in the Venetian alluvial paleoplain. A calcic paleosol, buried by Middle to Late Holocene marine transgressive deposits, represents a subsurface layer long known in the Venice area as “caranto.” Referenced in the ancient chronicles of architects and builders, the caranto exhibits relatively high compressive and shear strength, making it an important substrate for supporting building foundations, some dating back to the Gothic era (12th–15th centuries A.D.). Hence, the caranto paleosol documents local post‐glacial environmental changes while playing an important role in Venetian building construction and human settlement. Here we provide geochemical, sedimentological, paleoecological, and chronological analyses of the caranto paleosol and related deposits based on recent coring of the Venetian Lagoon. © 2011 Wiley Periodicals, Inc.     

Magnetic field modulated dust streams from Jupiter in interplanetary space

High speed dust streams emanating from near Jupiter were first discovered by the Ulysses spacecraft in 1992. Since then the phenomenon has been re-observed by Galileo in 1995, Cassini in 2000, and Ulysses in 2004. The dust grains are expected to be charged to a potential of , which is sufficient to allow the planet's magnetic field to accelerate them away from the planet, where they are subsequently influenced by the interplanetary magnetic field (IMF). A similar phenomenon was observed near Saturn by Cassini. Here, we report and analyze simultaneous dust, IMF and solar wind data for all dust streams from the two Ulysses Jupiter flybys. We find that compression regions (CRs) in the IMF – regions of enhanced magnetic field – precede most dust streams. Furthermore, the duration of a dust stream is roughly comparable with that of the precedent CR, and the occurrence of a dust stream and the occurrence of the previous CR are separated by a time interval that depends on the distance to Jupiter. The intensity of the dust streams and their precedent CRs are also correlated, but this correlation is only evident at distances from the planet no greater than 2 AU. Combining these observations, we argue that CRs strongly affect dust streams, probably by deflecting dust grain trajectories, so that they can reach the spacecraft and be detected by its dust sensor.     

Analysis of Titan CH4 3.3 μm upper atmospheric emission as measured by Cassini/VIMS

After molecular nitrogen, methane is the most abundant species in Titan’s atmosphere and plays a major role in its energy budget and its chemistry. Methane has strong bands at 3.3 μm emitting mainly at daytime after absorption of solar radiation. This emission is strongly affected by non-local thermodynamic equilibrium (non-LTE) in Titan’s upper atmosphere and, hence, an accurate modeling of the non-LTE populations of the emitting vibrational levels is necessary for its analysis. We present a sophisticated and extensive non-LTE model which considers 22 CH₄ levels and takes into account all known excitation mechanisms in which they take part. Solar absorption is the major excitation process controlling the population of the v₃-quanta levels above 1000 km whereas the distribution of the vibrational energy within levels of similar energy through collisions with N₂ is also of importance below that altitude. CH₄-CH₄ vibrational exchange of v₄-quanta affects their population below 500 km. We found that the ν₃ → ground band dominates Titan’s 3.3 μm daytime limb radiance above 750 km whereas the ν₃ + ν₄ → ν₄ band does below that altitude and down to 300 km. The ν₃ + ν₂ → ν₂, the 2ν₃ → ν₃, and the ¹³CH₄ν₃ → ground bands each contribute from 5% to 8% at regions below 800 km. The ν₃ + 2ν₄ → 2ν₄and ν₂ + ν₃ + ν₄ → ν₂ + ν₄ bands each contribute from 2% to 5% below 650 km. Contributions from other CH₄ bands are negligible. We have used the non-LTE model to retrieve the CH₄ abundance from 500 to 1100 km in the southern hemisphere from Cassini-VIMS daytime measurements near 3.3 μm. Our retrievals show good agreement with previous measurements and model results, supporting a weak deviation from well mixed values from the lower atmosphere up to 1000 km.     

Wind Tunnel And Field Measurements Of Turbulent Flow In Forests. Part I: Uniformly Thinned Stands

Many forest management methods alterstand density uniformly. The effectsof such a change on the wind andturbulence regimes in the forest arecritical to a number of processes governingthe stability of the stand and itsmicroclimate. We measured wind speed andturbulence statistics with a Dantec tri-axialhot-film probe in model forests of variousdensities (31–333 trees m^-2), created byremoving whole trees in a regular pattern in awind tunnel, and compared them with similarmeasurements made with propeller anemometers insimilarly thinned plots (156–625 trees ha^-1)within a Sitka spruce stand in Scotland. The results agree well, in general, with measurements made inother such studies with diverse canopy types.The systematic variations with density and verticalleaf-area distribution (which differed betweenwind-tunnel and field trees) in our work can explainmuch of the variability shown in scaled profiles ofbasic turbulence statistics reported in theliterature. The wind tunnel and field results are shown to be in good agreement overalldespite the difference in vertical leaf-areadistribution. Within-canopy and isolated-treedrag coefficients in the wind tunnel showthat tree-scale shelter effects increase astree density increases. The measurements indicatethat turbulence in the canopy is dominated bylarge-scale structures with dimensions of the sameorder as the height of the canopy as found inother studies but suggest that inter-tree spacing also modulates the size of these structures. These structures are associated with the sweeps that dominatemomentum exchange in the canopy and it is thisfact that allows the tri-axial probe to operate sowell despite the relatively narrow range of anglesin which the wind vector is correctly measured. Theratio of streamwise periodicity of these structuresto vorticity thickness varies systematically withtree density in the range 2.7–5.1, which spans theexpected range of 3.5–5 found in a laboratorymixing-layer, suggesting that tree spacing imposes another relevant length scale. This test andothers show that the results are in agreement withthe idea that canopy turbulence resembles that of a mixing layer even though they disagree with, and challenge the linear relationship between, streamwise periodicity andshear length scale presented recently in theliterature. The measurements are also in goodoverall agreement with simple drag models presented recently by other researchers.