Comparison of soil solution chemistry sampled by centrifugation, two types of suction lysimeters and zero-tension lysimeters

The choice of sampling method for soil solution is of great importance. In this paper soil solution chemistry sampled by centrifugation, two types of suction lysimeters and zero-tension lysimeters have been studied with the purpose of investigating systematic differences between them. The samples were taken at 4 depths from an acidified forest soil as well as from adjacent lime and ash treated soils. A centrifugation drainage method was compared with two types of suction lysimeters (‘Rhizon’ and ‘Prenart’) and zero-tension lysimeters. About half of the 27 variables measured showed a significant difference between the sampling methods used. Typically the centrifuged samples had lower pH (4.0 vs. 4.4), Ca (21 μM vs. 30 μM) and Mg (25 μM vs. 34 μM) concentrations and higher Cl (330 μM vs. 230 μM) and DOC (4.4 mM vs. 3.2 mM) concentrations than the Rhizon lysimeters. Also the other lysimeters showed significant differences compared to the centrifuged samples for about half the number of analytes. Centrifuged samples had higher concentrations of all analytes except NO₃ and PO₄ compared to zero-tension lysimeters and also for all analytes except NO₃ and Al compared to Prenart lysimeters. Among the environmental factors considered depth showed an influence to some extent, while sampling occasion had a great significant impact on the difference between the centrifugation method and the Rhizon lysimeters. Factors like individual pits or soil treatment did not show any influence on the difference between the methods.     

Precambrian fluvial deposits: Enigmatic palaeohydrological data from the c. 2–1.9 Ga Waterberg Group, South Africa

Precambrian fluvial systems, lacking the influence of rooted vegetation, probably were characterised by flashy surface runoff, low bank stability, broad channels with abundant bedload, and faster rates of channel migration; consequently, a braided fluvial style is generally accepted. Pre-vegetational braided river systems, active under highly variable palaeoclimatic conditions, may have been more widespread than are modern, ephemeral dry-land braided systems. Aeolian deflation of fine fluvial detritus does not appear to have been prevalent. With the onset of large cratons by the Neoarchaean–Palaeoproterozoic, very large, perennial braided river systems became typical. The c. 2.06–1.88 Ga Waterberg Group, preserved within a Main and a smaller Middelburg basin on the Kaapvaal craton, was deposited largely by alluvial/braided-fluvial and subordinate palaeo-desert environments, within fault-bounded, possibly pull-apart type depositories.

Palaeohydrological data obtained from earlier work in the Middelburg basin (Wilgerivier Formation) are compared to such data derived from the correlated Blouberg Formation, situated along the NE margin of the Main basin. Within the preserved Blouberg depository, palaeohydrological parameters estimated from clast size and cross-bed set thickness data, exhibit rational changes in their values, either in a down-palaeocurrent direction, or from inferred basin margin to palaeo-basin centre. In both the Wilgerivier and Blouberg Formations, calculated palaeoslope values (derived from two separate formulae) plot within the gap separating typical alluvial fan gradients from those which characterise rivers (cf. [Blair, T.C., McPherson, J.G., 1994. Alluvial fans and their natural distinction from rivers based on morphology, hydraulic processes, sedimentary processes, and facies assemblages. J. Sediment. Res. A64, 450–489.]). Although it may be argued that such data support possibly unique fluvial styles within the Precambrian, perhaps related to a combination of major global-scale tectono-thermal and atmospheric–palaeoclimatic events, a simpler explanation of these apparently enigmatic palaeoslope values may be pertinent. Of the two possible palaeohydrological formulae for calculating palaeoslope, one provides results close to typical fluvial gradients; the other formula relies on preserved channel-width data. We suggest that the latter will not be reliable due to problematic preservation of original channel-widths within an active braided fluvial system. We thus find no unequivocal support for a unique fluvial style for the Precambrian, beyond that generally accepted for that period and discussed briefly in the first paragraph.     

Geochemistry of silicic magmas in the Macolod Corridor, SW Luzon, Philippines: evidence of distinct, mantle-derived, crustal sources for silicic magmas

Silicic volcanic deposits (>65 wt% SiO₂), which occur as domes, lavas and pyroclastic deposits, are relatively abundant in the Macolod Corridor, SW Luzon, Philippines. At Makiling stratovolcano, silicic domes occur along the margins of the volcano and are chemically similar to the silicic lavas that comprise part of the volcano. Pyroclastic flows are associated with the Laguna de Bay Caldera and these are chemically distinct from the domes and lavas at Makiling stratovolcano. As a whole, samples from the Laguna de Bay Caldera contain lower concentrations of MgO and higher concentrations of Fe₂O_3(t) than the samples from domes and lavas. The Laguna de Bay samples are more enriched in incompatible trace elements. The silicic rocks from the domes, Makiling Volcano and Laguna de Bay Caldera all contain high alkalis and high K₂O/Na₂O ratios. Melting experiments of primitive basalts and andesites demonstrate that it is difficult to produce high K₂O/Na₂O silicic magmas by fractional crystallization or partial melting of a low K₂O/Na₂O source. However, recent melting experiments (Sisson et al., Contrib Mineral Petrol 148:635–661, 2005) demonstrate that extreme fractional crystallization or partial melting of K-rich basalts can produce these silicic magmas. Our model for the generation of the silicic magmas in the Macolod Corridor requires partial melting of mantle-derived, evolved, moderate to K-rich, crystallized calc-alkaline magmas that ponded and crystallized in the mid-crust. Major and trace element variations, along with oxygen isotopes and ages of the deposits, are consistent with this model. Electronic Supplementary Material Supplementary material is available for this article at     

Mineralogical and chemical evolution of the Ernest Henry Fe oxide–Cu–Au ore system, Cloncurry district, northwest Queensland, Australia

The Ernest Henry Cu–Au deposit was formed within a zoned, post-peak metamorphic hydrothermal system that overprinted metamorphosed dacite, andesite and diorite (ca 1740–1660 Ma). The Ernest Henry hydrothermal system was formed by two cycles of sodic and potassic alteration where biotite–magnetite alteration produced in the first cycle formed ca 1514±24 Ma, whereas paragenetically later Na–Ca veining formed ca 1529 +11/−8 Ma. These new U–Pb_titanite age dates support textural evidence for incursion of hydrothermal fluids after the metamorphic peak, and overlap with earlier estimates for the timing of Cu–Au mineralization (ca 1540–1500 Ma). A distal to proximal potassic alteration zone correlates with a large (up to 1.5 km) K–Fe–Mn–Ba enriched alteration zone that overprints earlier sodic alteration. Mass balance analysis indicates that K–Fe–Mn–Ba alteration—largely produced during pre-ore biotite- and magnetite-rich alteration—is associated with K–Rb–Cl–Ba–Fe–Mn and As enrichment and Na, Ca and Sr depletion. The aforementioned chemical exchange almost precisely counterbalances the mass changes associated with regional Na–Ca alteration. This initial transition from sodic to potassic alteration may have been formed during the evolution of a single fluid that evolved via alkali exchange during progressive fluid-rock interaction. Cu–Au ore, dominated by co-precipitated magnetite, minor specular hematite, and chalcopyrite as breccia matrix, forms a pipe-like body at the core of a proximal alteration zone dominated by K-feldspar alteration. Both the core and K-feldspar alteration overprint Na–Ca alteration and biotite–magnetite (K–Fe) alteration. Ore was associated with the concentration of a diverse range of elements (e.g. Cu, Au, Fe, Mo, U, Sb, W, Sn, Bi, Ag, F, REE, K, S, As, Co, Ba and Ca). Mineralization also involved the deposition of significant barite, K(–Ba)–feldspar, calcite, fluorite and complexly zoned pyrite. The complexly zoned pyrite and variable K–(Ba)–feldspar versus barite associations are interpreted to indicate fluctuating sulphur and/or barium supply. Together with the alteration zonation geochemistry and overprinting criteria, these data are interpreted to indicate that Cu–Au mineralization occurred as a result of fluid mixing during dilation and brecciation, in the location of the most intense initial potassic alteration. A link between early alteration (Na–Ca and K–Fe) and the later K-feldspathization and the Cu–Au ore is possible. However, the ore-related enrichments in particular elements (especially Ba, Mn, As, Mo, Ag, U, Sb and Bi) are so extreme compared with earlier alteration that another fluid, possibly magmatic in origin, contributed the diverse element suite geochemically independently of the earlier stages. Structural focussing of successive stages produced the distinctive alteration zoning, providing a basis both for exploration for similar deposits, and for an understanding of ore genesis.     

Land-surface scheme validation using the Oklahoma Atmospheric Surface-layer Instrumentation System (OASIS) and Oklahoma Mesonet data: Preliminary results

Summary The Oklahoma Atmospheric Surface-layer Instrumentation System (OASIS) is a recently-developed observational system that collects, archives, and quality controls atmospheric, surface, and soil data in real-time from 90 stations across Oklahoma. Ten of the 90 sites, termed “super sites”, are equipped with additional sonic anemometry and four-component net radiometers to provide complete observations of the surface energy balance. Oklahoma Mesonet and OASIS data are used in this study to validate the sensitivity and accuracy of a land-surface scheme within a numerical prediction model. The Advanced Regional Prediction System (ARPS) is a three-dimensional, nonhydrostatic mesoscale model developed by the Center for Analysis and Prediction of Storms (CAPS) at the University of Oklahoma. The land-surface model (LSM) used within ARPS is the Interactions Soil Biosphere Atmosphere (ISBA) scheme. Mesonet and OASIS data collected from the super site located in Norman, Oklahoma, are used as verification for the ISBA. Research presented in this study outlines the challenges in developing, maintaining, and using in-situ data for model validation. Such problems as instrument error, surface heterogeneity, and non-closure of the surface energy budget limit data accuracy. Preliminary results of model validation focus on the sensitivity of the soil physics within the ISBA scheme. Model sensitivity to vegetation cover, surface roughness, and soil type are investigated. Furthermore, several recent improvements to ISBA are evaluated and compared to observations. This study concludes that the sensitivity of the ISBA to a priori soil and vegetation type is detrimental for this scheme to be used in a mesoscale model without improved treatment of surface heterogeneity. Received November 18, 2001 Revised December 28, 2001     

Capture-zone design in an aquifer influenced by cyclic fluctuations in hydraulic gradients

Design of a groundwater pumping and treatment system for a wood-treatment facility adjacent to the tidally influenced Fraser River estuary required the development of methodologies to account for cyclic variations in hydraulic gradients. Design of such systems must consider the effects of these cyclic fluctuations on the capture of dissolved-phase contaminants. When the period of the cyclic fluctuation is much less than the travel time of the dissolved contaminant from the source to the discharge point, the hydraulic-gradient variations resulting from these cycles can be ignored. Capture zones are then designed based on the average hydraulic gradient determined using filter techniques on continuous groundwater-level measurements. When the period of cyclic fluctuation in hydraulic gradient is near to or greater than the contaminant travel time, the resulting hydraulic-gradient variations cannot be ignored. In these instances, procedures are developed to account for these fluctuations in the capture-zone design. These include proper characterization of the groundwater regime, assessment of the average travel time and period of the cyclic fluctuations, and numerical techniques which allow accounting for the cyclic fluctuations in the design of the capture zone. Electronic Publication     

Late Cretaceous reactivation of major crustal shear zones in northern Namibia: constraints from apatite fission track analysis

Namibia's passive continental margin records a long history of tectonic activity since the Proterozoic. The orogenic belt produced during the collision of the Congo and Kalahari Cratons in the Early Proterozoic led to a zone of crustal weakness, which became the preferred location for tectonism during the Phanerozoic. The Pan-African Damara mobile belt forms this intraplate boundary in Namibia and its tectonostratigraphic zones are defined by ductile shear zones, where the most prominent is described as the Omaruru Lineament–Waterberg Thrust (OML–WT). The prominance of the continental margin escarpment is diminished in the area of the Central and Northern Zone of the Damara belt where the shear zones are located. This area has been targeted with a set of 66 outcrop samples over a 550-km-long, 60-km-broad coast-parallel transect from the top of the escarpment in the south across the Damara sector to the Kamanjab Inlier in the north. Apatite fission track age and length data from all samples reveal a regionally consistent cooling event. Thermal histories derived by forward modelling bracket this phase of accelerated cooling in the Late Cretaceous. Maximum palaeotemperatures immediately prior to the onset of cooling range from ca. 120 to ca. 60 °C with the maximum occurring directly south of the Omaruru Lineament. Because different palaeotemperatures indicate different burial depth at a given time, the amount of denudation can be estimated and used to constrain vertical displacements of the continental crust. We interpret this cooling pattern as the geomorphic response to reactivation of basement structures caused by a change in spreading geometry in the South Atlantic and South West Indian Oceans.