Abrupt onset of carbonate deposition in Lake Kivu during the 1960s: response to recent environmental changes

This study interprets the recent history of Lake Kivu, a tropical lake in the East African Rift Valley. The current gross sedimentation was characterized from a moored sediment trap array deployed over 2 years. The past net sedimentation was investigated with three short cores from two different basins. Diatom assemblages from cores were interpreted as reflecting changes in mixing depth, surface salinity and nutrient availability. The contemporary sediment trap data indicate seasonal variability, governed by diatom blooms during the annual mixing in the dry season, similar to Lakes Malawi and Tanganyika. The ratio of settling fluxes to net sediment accumulation rates implies mineralization rates of 80–90% at the sediment-water interface. The sediment cores revealed an abrupt change ~40 years ago, when carbonate precipitation started. Since the 1960s, deep-water methane concentrations, nutrient fluxes and soil mineral inputs have increased considerably and diatom assemblages have altered. These modifications probably resulted from a combination of three factors, commonly altering lake systems: the introduction of a non-native fish species, eutrophication, and hydrological changes inducing greater upwelling. Both the fish introduction and increased rainfall occurred at the time when the onset of carbonate precipitation was observed, whereas catchment population growth accompanied by intensified land use increased the flux of soil minerals already since the early twentieth century due to more intense erosion.     

Mapping of erosion rates in marly badlands based on a coupling of anatomical changes in exposed roots with slope maps derived from LiDAR data

Black marls form very extensive outcrops in the Alps and constitute some of the most eroded terrains, thus causing major problems of sedimentation in artificial storage systems (e.g. reservoirs) and river systems. In the experimental catchments near Draix (France), soil erosion rates have been measured in the past at the plot scale through a detailed monitoring of surface elevation changes and at the catchment scale through continuous monitoring of sediment yield in traps at basin outlets. More recently, erosion rates have been determined by means of dendrogeomorphic techniques in three monitored catchments of the Draix basin. A total of 48 exposed roots of Scots pine have been sampled and anatomical variations in annual growth rings resulting from denudation analysed. At the plot scale, average medium‐term soil erosion rates derived from exposed roots vary between 1·8 and 13·8 mm yr^?1 (average: 5·9 mm yr^?1) and values are significantly correlated with slope angle. The dendrogeomorphic record of point‐scale soil erosion rates matches very well with soil erosion rates measured in the Draix basins. Based on the point‐scale measurements and dendrogeomorphic results obtained at the point scale, a linear regression model involving slope angle was derived and coupled to high‐resolution slope maps obtained from a LiDAR‐generated digital elevation model so as to generate high‐resolution soil erosion maps. The resulting regression model is statistically significant and average soil erosion rates obtained from the areal erosion map (5·8, 5·2 and 6·2 mm yr^?1 for the Roubine, Moulin and Laval catchments, respectively) prove to be well in concert with average annual erosion rates measured in traps at the outlet of these catchments since 1985 (6·3, 4·1 and 6·4 mm yr^?1). This contribution demonstrates that dendrogeomorphic analyses of roots clearly have significant potential and that they are a powerful tool for the quantification and mapping of soil erosion rates in areas where measurements of past erosion is lacking. Copyright © 2011 John Wiley & Sons, Ltd.     

Coupled discrete element and finite volume solution of two classical soil mechanics problems

One dimensional solutions for the classic critical upward seepage gradient/quick condition and the time rate of consolidation problems are obtained using coupled routines for the finite volume method (FVM) and discrete element method (DEM), and the results compared with the analytical solutions. The two phase flow in a system composed of fluid and solid is simulated with the fluid phase modeled by solving the averaged Navier–Stokes equation using the FVM and the solid phase is modeled using the DEM. A framework is described for the coupling of two open source computer codes: YADE-OpenDEM for the discrete element method and OpenFOAM for the computational fluid dynamics. The particle–fluid interaction is quantified using a semi-empirical relationship proposed by Ergun [12]. The two classical verification problems are used to explore issues encountered when using coupled flow DEM codes, namely, the appropriate time step size for both the fluid and mechanical solution processes, the choice of the viscous damping coefficient, and the number of solid particles per finite fluid volume.     

A new scenario for the last magmatic eruption of La Soufrière of Guadeloupe (Lesser Antilles) in 1530 A.D. Evidence from stratigraphy radiocarbon dating and magmatic evolution of erupted products

La Soufrière of Guadeloupe is a dangerous volcano characterized over the last decade by moderate seismic and fumarolic unrest. In the last 15,000 years it has experienced phreatic and magmatic eruptions and unusually numerous flank collapse events sometimes associated with a magmatic eruption. We propose a new age of 1530 A.D. and a new eruptive scenario for the last magmatic eruption on the basis of a novel statistical analysis of radiocarbon age dates, and new field and geochemical data. This eruption is the only magmatic eruption likely to have occurred in Guadeloupe during the last 1400 years. The eruption mainly involved an andesitic magma which, in the first phase of the eruption, partially mixed with a slightly more differentiated magma stored in a small and shallow magma chamber. Ascent of magma to the surface generated a partial collapse of the hydrothermally altered edifice that increased the magma discharge and led to a sub-plinian phase with scoria fallout and column-collapse pyroclastic flows followed by near-vent pyroclastic scoria fountains. The eruption ended with growth of a lava dome. Our revised interpretation of the last magmatic eruption of La Soufrière constitutes the most likely key to a future magmatic eruption scenario for this volcano which displays strong evidence of unrest since 1992.     

Chromite crystallization in a multicellular magma flow: Evidence from a chromitite dike in the Oman ophiolite

The investigated chromitite dike is located at the top of an upwelling mantle structure of the Oman ophiolite (Maqsad diapir), in undeformed dunites displaying evidence for magma impregnation and circulation, just below the paleo-ridge axis. The chromitite dike is undeformed, its shape is that of an upward widening tube. It exhibits an internal layering which is roughly perpendicular to the cavity axis and comprises a vertical succession of four main layers showing a graded-bedding. Chromitite magmatic structures are beautifully preserved and result from a progressive crystallization from small euhedral crystals to wide octahedron-shaped nodules; dissolution textures provide evidence for late magmatic desequilibrium; sedimentation structures include flattening of the largest nodules. The silicate matrix comprises poikilitic forsterite and a locally abundant association of primary pargasite and plagioclase and alteration minerals (vesuvianite-chlorite-dolomite); pargasite inclusions are very abundant in the chromite. Chromite composition changes from one layer to the other and from core to rim in the chromite nodules (chromium decreases and titanium increases); Ti contents are generally high (0.4 to 0.8 wt.% TiO₂) with respect to podiform chromites. Platinum-group elements are not abundant but they show a strong fractionation at the scale of the orebody and of the main graded-bedded layers (Pd/Ir ratio varies from 0.5 to 11.5). REE patterns of chromitite parallel to those of gabbros and furthermore display a sea water related hydrothermal alteration (Ce negative anomaly).

The chromitite dike of Maqsad provides evidence for the crystallization of chromitite bodies in subvertical magma conduits below oceanic ridges; it corroborates the model of Cassard et al. (1981) and Lago et al. (1982) concerning the formation of chromitite pods in ophiolites which were later deformed and transposed into the horizontal plane due to the plastic flow prevailing away from the paleo-axial zone. Layering and chromite compositional variations are ascribed to a multicellular convective system segregating various stocks of chromite particles either in the upwelling flow of fresh magma or in the convective cells of fractionated residual magma in the confined part of the cavity. The estimated life-time for the magma influx is very short (<2 months). The parent-magma was probably of MORB-type and already fractionated (Ti-rich and PGE-poor), which is consistent with the strong evidence of magma-peridotite interactions in the core of the Maqsad diapir. Hydrous fluids were present during chromite crystallization (pargasite inclusions) suggesting that fluid-rich melts occur in the upper mantle.     

Behaviour of REE and mass balance calculations in a lateritic profile over chlorite schists in South Cameroon

An in situ weathering profile overlying chlorite schists in the Mbalmayo-Bengbis formations (South Cameroon) was chosen for the study of the behaviour of REE and the evaluation of geochemical mass balance. After physical and mineralogical studies, the chlorite schists and the undisturbed weathered materials were chemically analyzed for major elements (X-ray fluorescence and titrimetry) and REE (ICP-MS). The behaviour of the REE in the Mbalmayo weathering system was established in comparison with the REE of the reference parent rock. Mass balance calculations were applied to both major elements and REE. The mineralogy of the materials was determined with the aid of a Philips 1720, diffractometer. The chlorite schists of the Mbalmayo sector show low REE contents (Σ=153.44 ppm). These rocks are relatively rich in LREE (about 125 times the chondritic value) and relatively poor in HREE (about 20 times the chondritic value). The REE diagram normalized to chondrites shows a slightly split graph ((La/Yb)_N=6.18) with marked enrichment in LREE (LREE/HREE=9.50) in relation to HREE. Moreover, these spectra do not present any Ce anomaly, but a slightly positive Eu anomaly. The imperfectly evolved profile, whose materials are genetically linked, shows an atypical behaviour of REE. In effect, the LREE are more mobile than the HREE during weathering ((La/Yb)NASC<1) with weak Ce anomalies. This has been rarely reported in lateritic profiles characterized by higher HREE mobility than LREE during weathering processes with high Ce anomalies. This is either due to the difference in the stability of REE-bearing minerals, or to the weak acidic to basic pH conditions (6.70<pH<7.80), or even due to the average evolution of the weathering materials. The pathway of the REE along the profile is as follows: (1) leaching in the saprolites and summit of the profile, except for Ce, which precipitates very weakly in the nodular materials and the coarse saprolite materials, (2) at the base of the profile, solutions come in contact with chlorite schist formations, at this level, the pH increases (pH=7.79), HREE and a part of LREE partially void of Ce precipitate and (3) the other part of LREE precipitates further up in the profile. The geochemical mass balance calculations reveal that these elements are leached in the same phases as the relatively high Si, Al, K and Fe²⁺ contents.     

Arsenite sorption at the magnetite-water interface during aqueous precipitation of magnetite: EXAFS evidence for a new arsenite surface complex

The interaction of aqueous As(III) with magnetite during its precipitation from aqueous solution at neutral pH has been studied as a function of initial As/Fe ratio. Arsenite is sequestered via surface adsorption and surface precipitation reactions, which in turn influence the crystal growth of magnetite. Sorption samples were characterized using EXAFS spectroscopy at the As K-edge in combination with HRTEM observations, energy dispersive X-ray analysis at the nanoscale, electron energy loss spectroscopy at the Fe L₃-edge, and XRD-Rietveld analyses of reaction products. Our results show that As(III) forms predominantly tridentate hexanuclear As(III)O₃ complexes (³C), where the As(III)O₃ pyramids occupy vacant tetrahedral sites on {1 1 1} surfaces of magnetite particles. This is the first time such a tridentate surface complex has been observed for arsenic. This complex, with a dominant As-Fe distance of 3.53 ± 0.02 Å, occurs in all samples examined except the one with the highest As/Fe ratio (0.33). In addition, at the two highest As/Fe ratios (0.133 and 0.333) arsenite tends to form mononuclear edge-sharing As(III)O₃ species (²E) within a highly soluble amorphous As(III)-Fe(III,II)-containing precipitate. At the two lowest As/Fe ratios (0.007 and 0.033), our results indicate the presence of additional As(III) species with a dominant As-Fe distance of 3.30 ± 0.02 Å, for which a possible structural model is proposed. The tridentate ³C As(III)O₃ complexes on the {1 1 1} magnetite surface, together with this additional As(III) species, dramatically lower the solubility of arsenite in the anoxic model systems studied. They may thus play an important role in lowering arsenite solubility in putative magnetite-based water treatment processes, as well as in natural iron-rich anoxic media, especially during the reductive dissolution-precipitation of iron minerals in anoxic environments.     

Seismic analysis of non-linearly base-isolated soil-structure interacting reactor building by way of the hybrid frequency-time-domain procedure

The seismic analysis of a reactor building on a sliding-type base isolation is performed by way of the hybrid frequency-time-domain procedure. The frequency dependence of the foundation stiffness coefficients is duly accounted for in the analysis, although the problem is non-linear. The response results are shown to be reliable by way of comparison with the results of a time-stepping algorithm in the specialized case of constant foundation stiffness coefficients. The fact that such an analysis can be performed is an outstanding result, particularly when considering that no difficulties are encountered in the implementation and in the calculation. The flow chart used to implement the procedure is presented in the paper. The reactor building is analysed for three different sites with each site having three sets of material properties. It appears that the structural response is not strongly affected by the frequency dependence of the foundation stiffness coefficients. This also holds true as far as the response of equipment attached to the superstructure is concerned. For all practical purposes this frequency dependence may be disregarded in favour of constant spring and damping coefficients.     

Study of two alteration systems as natural analogues for radionuclide release and migration

U-deposit hosted in hydrothermally altered tuffs in Mexico, together with weathering profiles from Cameroon were studied as natural analogues of radionuclide release and migration. Using petrological and spectroscopic methods (infrared and electron paramagnetic resonance), we have distinguished successive secondary mineral parageneses and the behaviour of radionuclides.

In the U-deposit, the mineral parageneses show that uranium migration is mainly controlled by the redox potential and silica activity of the altering solutions. The high silica content of the solutions is caused by the intense alteration of volcanic rocks. Two types of secondary clay mineral parageneses are evidenced: a kaolinization, intense where uranium is accumulated in the welded tuffs, and a smectitization mainly developed in the underlying weakly welded tuffs.

Several types of kaolinite have been defined according to their genesis (fillings in fissures and feldspar pseudomorphs), their location relative to a breccia pipe where uranium has accumulated (core and rim of the pipe; surrounding rhyolitic tuffs), and particle morphology, structural order and substitutional Fecontent. It is shown that the variations of the concentration of paramagnetic defect centres, always more than ten times as important than those measured in weathering kaolinites, are only correlated to the location of the kaolinites. The highest values correspond to the breccia pipe kaolinites, e.g. kaolinites located in the uranium accumulation zones. Moreover, one or two main defects centres are detected depending on the intimate association of kaolinites with uranium-bearing minerals. Besides, in weathering kaolinites from U-depleted laterites, defect centre concentrations are correlated to the total Fe203 content in bulk samples. This means that the defect centre acts as a memory of the travel of uranium when this element was sorbed onto iron gels in the first stage of weathering.

It is concluded that paramagnetic defect centres in kaolinites might allow an efficient fingerprint of successive irradiations in the natural analogues under study and could be an useful tool to control radionuclides migration through kaolinite-containing clayey materials such as those used for waste repository.

A better understanding of radiation efficiency as well as accurate dose-ratekaolinite-containing clayey materials such as those used for waste reposit estimation are needed for a quantitative tracing of the migration ofA better understanding of radiation efficiency as well as accurate dose-ratekaolin radionuclide elements. With this aim, a simulation has been undertaken withestimation are needed for a quantitative tracing of the migration of various radiations sources. We have determined for each irradiation the parameters of the paramagnetic centres created in order to understand the way they are forming. The knowledge of the parameters governing the formation and the stability of the radiation centres in kaolinites allow to use this mineral as a natural dosimeter.     

Electron density profiles above the F-layer

A method for the determination of electron densities up to 1000 km by using data from vertical sounding and Faraday rotation techniques is suggested. Electron density profiles and contours plotted following this method are compared with results obtained by Evans using the incoherent scatter technique.