Normative quartz as an indicator of the mass transfer intensity during the postmagmatic alteration of the Botuobinskaya pipe kimberlites (Yakutia)

The intensity of postmagmatic processes in the Botuobinskaya pipe kimberlites was estimated from the calculated content of normative secondary quartz (Q). Several simple algorithms are proposed to calculate the Q content from chemical analyses of kimberlites. Ten groups of altered kimberlites have been recognized from the Q contents. The contents of MgO, some trace elements, and LREE in the groups, the contents of Cr and Ca of crimson garnets, the diamond contents of kimberlites, and the average weight of diamonds decrease as the Q content increases. It is shown that the negative SiO₂–MgO correlation is the most effective indicator of the postmagmatic alteration of kimberlites. As the degree of their secondary alteration increases, the kimberlites transform into an assemblage of quartz and clay minerals enriched in some trace elements and almost completely lacking REE and diamonds.     

Integration of local–global upscaling and grid adaptivity for simulation of subsurface flow in heterogeneous formations

We propose a methodology, called multilevel local–global (MLLG) upscaling, for generating accurate upscaled models of permeabilities or transmissibilities for flow simulation on adapted grids in heterogeneous subsurface formations. The method generates an initial adapted grid based on the given fine-scale reservoir heterogeneity and potential flow paths. It then applies local–global (LG) upscaling for permeability or transmissibility [7], along with adaptivity, in an iterative manner. In each iteration of MLLG, the grid can be adapted where needed to reduce flow solver and upscaling errors. The adaptivity is controlled with a flow-based indicator. The iterative process is continued until consistency between the global solve on the adapted grid and the local solves is obtained. While each application of LG upscaling is also an iterative process, this inner iteration generally takes only one or two iterations to converge. Furthermore, the number of outer iterations is bounded above, and hence, the computational costs of this approach are low. We design a new flow-based weighting of transmissibility values in LG upscaling that significantly improves the accuracy of LG and MLLG over traditional local transmissibility calculations. For highly heterogeneous (e.g., channelized) systems, the integration of grid adaptivity and LG upscaling is shown to consistently provide more accurate coarse-scale models for global flow, relative to reference fine-scale results, than do existing upscaling techniques applied to uniform grids of similar densities. Another attractive property of the integration of upscaling and adaptivity is that process dependency is strongly reduced, that is, the approach computes accurate global flow results also for flows driven by boundary conditions different from the generic boundary conditions used to compute the upscaled parameters. The method is demonstrated on Cartesian cell-based anisotropic refinement (CCAR) grids, but it can be applied to other adaptation strategies for structured grids and extended to unstructured grids.     

Accurate local upscaling with variable compact multipoint transmissibility calculations

We propose a new single-phase local upscaling method that uses spatially varying multipoint transmissibility calculations. The method is demonstrated on two-dimensional Cartesian and adaptive Cartesian grids. For each cell face in the coarse upscaled grid, we create a local fine grid region surrounding the face on which we solve two generic local flow problems. The multipoint stencils used to calculate the fluxes across coarse grid cell faces involve the six neighboring pressure values. They are required to honor the two generic flow problems. The remaining degrees of freedom are used to maximize compactness and to ensure that the flux approximation is as close as possible to being two-point. The resulting multipoint flux approximations are spatially varying (a subset of the six neighbors is adaptively chosen) and reduce to two-point expressions in cases without full-tensor anisotropy. Numerical tests show that the method significantly improves upscaling accuracy as compared to commonly used local methods and also compares favorably with a local–global upscaling method.     

Fractional crystallization models and B–Be–Li systematics at Mt Somma-Vesuvius volcano (Southern Italy)

A compilation of B–Be–Li data on rocks that cover the entire eruptive history of Somma-Vesuvius is presented and interpreted in the light of evolution models for the Somma-Vesuvius rocks. Using major and trace element data, fractional crystalllization models are presented for different geochemical units. These data were used to constrain the source mineralogy of the Somma-Vesuvius rocks (ol-opx-cpx-gar-amp of 0.4-0.3-0.1-0.1-0.1), the amount of sediment added (5–10%) and the melt fraction from batch partial melting computations (0.05–0.1). From the B–Li data it is inferred that the main process responsible for the B isotopic signature is sediment recycling. However, the B–Li data show a major variation in Li abundances respect to B which is explained with Li dehydration before the fluid enriched the mantle wedge that produced the arc magmas. The Somma-Vesuvius B isotope composition is intermediate between that of the Campi Flegrei and the broad field of the Eolian Island arc. A low Be isotopes in the recent volcanic rocks can be explained as: (a) the top 1–22 m of the incoming sediment is accreted, (b) large amounts of sediment erosion, (c) a slow rate of subduction which have provoked a long magmatic history for the Vesuvius magma, (d) the sediment component takes several Myr longer than the subducting plate to reach the magma source region beneath Italy.     

Mt. Melbourne volcano,antarctica: Evidence of seismicity related to volcanic activity

Characterization of the microseismic activity (M _L <2.0) has been performed at Mt. Melbourne since 1990. We recorded a group of low frequency events with common morphological characteristics, i.e., an emerging onset, an unclear second phase and a sharply dropping coda. Spectral analysis of events recorded at more than one station indicates that the seismogram characteristics and spectral content are largely due to source effects. A polarization filter applied to a set of three component data revealed a first phase made up ofP waves followed (after about 0.9–1.4 sec.) by a second phase probably composed ofSH-type waves. Particle motion analysis detected a seismic ray angle direction mainly between N70°E and N110°E and apparent angle of incidence between 35° and 48° for the first phase. The studied seismicity was localized in an area on the eastern slope of Mt. Melbourne Volcano which presents a surface temperature anomaly (Mazzarini andSalvini, 1994). We formulate two hypotheses for the type of earthquakes recorded: 1) long-period events involving active presence of magmatic fluids in the source processes; 2) or the result of fracturing processes (shear?) in a medium characterized by transition between brittle and plastic behaviors. In the latter hypothesis the superficial thermal anomaly may be a symptom of this behavior at depth and is confirmed by the lown values observed for the exponential fit in the codaQ analysis.     

Fe-Cu-Zn-Mo同位素示踪氧化还原过程

非传统稳定同位素（Fe-Cu-Zn-Mo）理论与数据相结合提高了科研工作者对地质体系氧化还原过程的理解。本文对这一相对较新的领域进行了综述,包括与氧化还原过程相关的同位素分馏理论和实验约束、时空尺度下的氧逸度以及同位素示踪氧化还原过程。稳定同位素理论预测,Fe-Cu-Zn-Mo同位素应该对氧化还原状态的变化能够做出响应。结果表明,Fe同位素作为岩浆过程、表生过程、俯冲带流体性质"氧逸度计"应用前景广阔;Cu同位素在岩浆、热液、陆地系统可以很好地示踪氧化还原过程;Zn同位素由于络合过程分馏已经被用在许多不同环境中作为含硫/碳流体迁移的敏感示踪剂;Mo同位素作为古氧逸度计可有效重建古海洋-大气氧化还原状态。     

Aragonite cycles: diagenesis caught in the act

Organic biomarker and nannofossil assemblages are used in combination with sedimentary petrology to identify the role of diagenesis for the formation of a rhythmic calcareous succession. A record from the Early Pliocene of the Maldives (Ocean Drilling Program, ODP Leg 115) is presented containing precession and eccentricity cycles expressed as variable aragonite content. Additional sub‐Milankovitch cycles are caused by rhythmic precipitation of calcite cement in the lower part of the interval. Comparison with palaeo‐productivity indicators (nannofossils, chlorin, total organic carbon) suggests that cementation occurs preferentially in intervals characterized by increasing or decreasing productivity. The coupled variability in productivity and carbonate diagenesis is attributed to the effect of organic matter degradation in the sediment. The observed combination of primary and diagenetic factors hampers the interpretation of the cyclicity on Milankovitch and sub‐Milankovitch scales. Diagenetically stable proxies for palaeo‐productivity, such as nannofossil assemblage data, were used to distinguish between palaeoclimate and diagenetic influences.     

Mais comment s'écoule donc un glacier ? Aperçu historique

Ice and snow have often helped physicists understand the world. On the contrary it has taken them a very long time to understand the flow of the glaciers. Naturalists only began to take an interest in glaciers at the beginning of the 19th century during the last phase of glacier advances. When the glacier flow from the upslope direction became obvious, it was then necessary to understand how it flowed. It was only in 1840, the year of the Antarctica ice sheet discovery by Dumont d'Urville, that two books laid the basis for the future field of glaciology: one by Agassiz on the ice age and glaciers, the other one by canon Rendu on glacier theory. During the 19th century, ice flow theories, adopted by most of the leading scientists, were based on melting/refreezing processes. Even though the word ‘fluid’ was first used in 1773 to describe ice, more the 130 years would have to go by before the laws of fluid mechanics were applied to ice. Even now, the parameter of Glen's law, which is used by glaciologists to model ice deformation, can take a very wide range of values, so that no unique ice flow law has yet been defined. To cite this article: F. Rémy, L. Testut, C. R. Geoscience 338 (2006).     

Flow and form in rehabilitation of large-river ecosystems: An example from the Lower Missouri River

On large, intensively engineered rivers like the Lower Missouri, the template of the physical habitat is determined by the nearly independent interaction of channel form and flow regime. We evaluated the interaction between flow and form by modeling four combinations of modern and historical channel form and modern and historical flow regimes. The analysis used shallow, slow water (shallow-water habitat, SWH, defined as depths between 0 and 1.5 m, and current velocities between 0 and 0.75 m/s) as an indicator of habitat that has been lost on many intensively engineered rivers and one that is thought to be especially important in rearing of young fishes. Two-dimensional hydrodynamic models for modern and historical channels of the Lower Missouri River at Hermann, Missouri, indicate substantial differences between the two channels in total availability and spatial characteristics of SWH. In the modern channel, SWH is maximized at extremely low flows and in overbank flows, whereas the historical channel had substantially more SWH at all discharges and SWH increased with increasing discharge. The historical channel form produced 3–7 times the SWH area of the modern channel regardless of flow regime. The effect of flow regime is evident in increased within-year SWH variability with the natural flow regime, including significant seasonal peaks of SWH associated with spring flooding. Comparison with other reaches along the Lower Missouri River indicates that a) channel form is the dominant control of the availability of habitat even in reaches where the hydrograph is more intensively altered, and b) rehabilitation projects that move toward the historical condition can be successful in increasing topographic diversity and thereby decreasing sensitivity of the availability of habitat to flow regime. The relative efficacy of managing flow and form in creating SWH is useful information toward achieving socially acceptable rehabilitation of the ecosystem in large river systems.