Managed water levels and the expansion of emergent vegetation along a lakeshore

Water-level fluctuations may be used to promote the expansion of emergent vegetation along lakeshores. We present the case of the lake Volkerak-Zoommeer in the Netherlands, a fresh-water lake created in 1987 after the enclosure of an estuary. Using an experimental area in which the water level could be manipulated, it was shown that partial summer drawdown of the shoreline created suitable conditions for germination and growth of tall emergent species (in particular Phragmites australis). Plant survival and growth depended on subsequent water-level fluctuations and grazing by waterbirds. Based on the experiment and empirical data, a model was developed to predict the effects of the water-level regime on potential reed bed development. The model was applied for four hydrological scenarios that have been considered for the water-level management of the lake.     

Mineral dissolution in the Cape Cod aquifer, Massachusetts, USA: I . Reaction stoichiometry and impact of accessory feldspar and glauconite on strontium isotopes, solute concentrations, and REY distribution

To compare relative reaction rates of mineral dissolution in a mineralogically simple groundwater aquifer, we studied the controls on solute concentrations, Sr isotopes, and rare earth element and yttrium (REY) systematics in the Cape Cod aquifer. This aquifer comprises mostly carbonate-free Pleistocene sediments that are about 90% quartz with minor K-feldspar, plagioclase, glauconite, and Fe-oxides. Silica concentrations and pH in the groundwater increase systematically with increasing depth, while Sr isotopic ratios decrease. No clear relationship between ⁸⁷Sr/⁸⁶Sr and Sr concentration is observed. At all depths, the ⁸⁷Sr/⁸⁶Sr ratio of the groundwater is considerably lower than the Sr isotopic ratio of the bulk sediment or its K-feldspar component, but similar to that of a plagioclase-rich accessory separate obtained from the sediment. The Si-⁸⁷Sr/⁸⁶Sr-depth relationships are consistent with dissolution of accessory plagioclase. In addition, solutes such as Sr, Ca, and particularly K show concentration spikes superimposed on their respective general trends. The K-Sr-⁸⁷Sr/⁸⁶Sr systematics suggests that accessory glauconite is another major solute source to Cape Cod groundwater. Although the authigenic glauconite in the Cape Cod sediment is rich in Rb, it is low in in-grown radiogenic ⁸⁷Sr because of its young Pleistocene age. The low ⁸⁷Sr/⁸⁶Sr ratios are consistent with equilibration of glauconite with seawater. The impact of glauconite is inferred to vary due to its variable abundance in the sediments. In the Cape Cod groundwater, the variation of REY concentrations with sampling depth resembles that of K and Rb, but differs from that of Ca and Sr. Shale-normalized REY patterns are light REY depleted, show negative Ce anomalies and super-chondritic Y/Ho ratios, but no Eu anomalies. REY input from feldspar, therefore, is insignificant compared to input from a K-Rb-bearing phase, inferred to be glauconite. These results emphasize that interpretation of groundwater chemistry, even in relatively simple aquifers, may be complicated by solute contributions from “exotic” accessory minerals such as glauconite. To detect such peculiarities, groundwater studies should combine the study of elemental concentration and isotopic composition of several solutes that show different geochemical behavior.     

The effect of dissolved barium on biogeochemical processes at cold seeps

A numerical model was applied to investigate and quantify the biogeochemical processes fueled by the expulsion of barium and methane-rich fluids in the sediments of a giant cold-seep area in the Derugin Basin (Sea of Okhotsk). Geochemical profiles of dissolved Ba²⁺, Sr²⁺, Ca²⁺, SO₄²⁻, HS⁻, DIC, I⁻ and of calcium carbonate (CaCO₃) were fitted numerically to constrain the transport processes and the kinetics of biogeochemical reactions. The model results indicate that the anaerobic oxidation of methane (AOM) is the major process proceeding at a depth-integrated rate of 4.9 μmol cm⁻² a⁻¹, followed by calcium carbonate and strontian barite precipitation/dissolution processes having a total depth-integrated rate of 2.1 μmol cm⁻² a⁻¹. At the low seepage rate prevailing at our study site (0.14 cm a⁻¹) all of the rising barium is consumed by precipitation of barite in the sedimentary column and no benthic barium flux is produced. Numerical experiments were run to investigate the response of this diagenetic environment to variations of hydrological and biogeochemical conditions. Our results show that relatively low rates of fluid flow (<∼5 cm a⁻¹) promote the dispersed precipitation of up to 26 wt% of barite and calcium carbonate throughout the uppermost few meters of the sedimentary column. Distinct and persistent events (several hundreds of years long) of more vigorous fluid flow (from 20-110 cm a⁻¹), instead, result in the formation of barite-carbonate crusts near the sediment surface. Competition between barium and methane for sulfate controls the mineralogy of these sediment precipitates such that at low dissolved methane/barium ratios (<4-11) barite precipitation dominates, while at higher methane/barium ratios sulfate availability is limited by AOM and calcium carbonate prevails. When seepage rates exceed 110 cm a⁻¹, barite precipitation occurs at the seafloor and is so rapid that barite chimneys form in the water column. In the Derugin Basin, spectacular barite constructions up to 20 m high, which cover an area of roughly 22 km² and contain in excess of 5 million tons of barite, are built through this process. In these conditions, our model calculates a flux of barium to the water column of at least 20 μmol cm⁻² a⁻¹. We estimate that a minimum of 0.44 × 10⁶ mol a⁻¹ are added to the bottom waters of the Derugin Basin by cold seep processes, likely affecting the barium cycle in the Sea of Okhotsk.     

Decoupling of As and Fe release to Bangladesh groundwater under reducing conditions. Part I: Evidence from sediment profiles

This study reexamines the notion that extensive As mobilization in anoxic groundwater of Bangladesh is intimately linked to the dissolution of Fe oxyhydroxides on the basis of analyses performed on a suite of freshly collected samples of aquifer material. Detailed sediment profiles extending to 40 to 70 m depth below the surface were obtained at six sites where local groundwater As concentrations were known to span a wide range. The sediment properties that were measured include (1) the proportion of Fe(II) in the Fe fraction leached in hot 1.2 N HCl, (2) diffuse spectral reflectance, and (3) magnetic susceptibility.In parallel with local concentrations of dissolved As ranging from <5 to 600 μg/L, Fe(II)/Fe ratios in shallow (gray) Holocene sands tended to gradually increase with depth from values of 0.3 to 0.5 to up to 0.9. In deeper (orange) aquifers of presumed Pleistocene age that were separated from shallow sands by a clay layer and contained <5 μg/L dissolved As, leachable Fe(II)/Fe ratios averaged ∼0.2. There was no consistent relation between sediment Fe(II)/Fe and dissolved Fe concentrations in groundwater in nearby wells. The reflectance measurements indicate a systematic linear relation (R² of 0.66; n = 151) between the first derivative transform of the reflectance at 520 nm and Fe(II)/Fe. The magnetic susceptibility of the shallow aquifer sands ranged from 200 to 3600 (x 10⁻⁹ m³/kg SI) and was linearly related (R² of 0.75; n = 29) to the concentrations of minerals that could be magnetically separated (0.03 to 0.79% dry weight). No systematic depth trends in magnetic susceptibility were observed within the shallow sands, although the susceptibility of deeper low-As aquifers was low (up to ∼200 × 10⁻⁹ m³/kg SI).This set of observations, complemented by incubation results described in a companion paper by van Geen et al. (this volume), suggests that the release of As is linked to the transformation of predominantly Fe (III) oxyhydroxide coatings on sand particles to Fe(II) or mixed Fe(II/III) solid phases with a flatter reflectance spectrum such as siderite, vivianite, or magnetite, without necessarily resulting in the release of Fe to groundwater. The very low As/Fe ratio of magnetically separated minerals compared to the As/Fe of bulk acid leachate (2 vs. 40 10⁻⁶, respectively) suggests that such a transformation could be accompanied by a significant redistribution of As to a mobilizable phase on the surface of aquifer particles.     

Three-dimensional forward and backward modelling of diapirism: numerical approach and its applicability to the evolution of salt structures in the Pricaspian basin

Numerical studies of ductile deformations induced by salt movements have, until now, been restricted to two-dimensional (2D) modelling of diapirism. This paper suggests a numerical approach to model the evolution of three-dimensional (3D) salt structures toward increasing maturity. This approach is also used here to restore the evolution of salt structures through successive earlier stages. The numerical methodology is applied to study several model examples. We analyse a model of salt diapirs that develop from an initial random perturbation of the interface between salt and its overburden and restore the evolved salt diapirs to their initial stages. We show that the average restoration errors are less than 1%. An evolutionary model of a 2D salt wall loaded by a 2D pile of sediments predicts a decomposition of the salt wall into 3D diapiric structures when the overburden of salt is supplied by 3D synkinematic wedge of sediments. We model salt extrusion feeding a gravity current over the depositional surface and estimate an average rate of extrusion and horizontal velocity of salt spreading. Faulting of the overburden to salt overhangs initiates new secondary diapirs, and we analyse the growth of these secondary diapirs. We also study how lateral flow effects the evolution of salt diapirs. The shape of a salt diapir can be very different if the rate of horizontal flow is much greater than the initial rate of diapiric growth solely due to gravity. We discuss the applicability of the results of the models to the evolution of Late Permian salt structures in the Pricaspian basin (Russia and Kazakhstan). These structures are distinguishable into a variety of styles representing different stages of growth: salt pillows, diapirs, giant salt massifs, 2D diapiric walls and 3D stocks complicated by large overhangs. The different sizes, shapes and maturities of salt structures in different parts of the Pricaspian basin reflect areal differences in salt thickness and loading history. Our results suggest that the numerical methodology can be employed to analyse the evolution of all salt structures that have upbuilt through younger ductile overburdens.     

Crustal thickening in an extensional regime: application to the mid-Norwegian Vøring margin

The structure of the mid-Norwegian volcanic Vøring margin at the onset of the Maastrichtian–Paleocene extension phase reflects the cumulative effect of earlier consecutive rifting events. Lateral structural differences present on the margin at that time are a consequence of migration of the location of maximum extension in time between Norway and Greenland. The most important imprints (Moho depth, thermal structure) of these events on the lithosphere are incorporated in a numerical simulation of the final extension phase. We focus on a possible mechanism of formation of the Vøring Marginal High and address the relationship between spatial and temporal evolution of crustal thinning and thickening, uplift of the surface and strength of the lithosphere.It is found that the Vøring Basin formed the strongest part of the margin which explains why the Maastrichtian–Paleocene rift axis was not located here but instead jumped westward with respect to the earlier rift axes locations. The modeling study predicts that local crustal thickening during extension can be expected when large lateral thermal variations are present in the lithosphere at the onset of extension. Negative buoyancy induced by lateral temperature differences increases downwelling adjacent to the rifting zone; convergence of material at the particular part of the margin is mainly taken up by the lower crust. The model shows that during the final phase of extension, the crust in the Vøring Marginal High area was thickened and the surface uplifted. It is likely that this dynamic process and the effects of magmatic intrusions both acted in concert to form the Marginal High.