About the specification of erosion flux for soft stratified cohesive sediments

The present paper deals with the specification of bed erosion flux that accounts for the effects of sediment-induced stratification in the water column. Owing to difficulties in measuring the bed shear stress _b and the erosive shear strength _s, we suggest a series of methods that combine laboratory and numerical experiments. A simplified turbulent transport model that includes these effects helps to quantify _b and _s. Focusing on soft stratified beds, the present study considers erosion rate formulas of the form =_f exp {[T_b-T_s]} where is a model constant (=1 for Gularte's (1978) formula and =1/2 for Parchure's (1984) formula). First, the bed erosive strength profile _s(Z) is adjusted by forcing the turbulent transport model with measured erosion rates. Second, three procedures are suggested to determine the erosion rate formula coefficients _f and : a global procedure and two different layer-by-layer procedures. Each procedure is applied to an erosion experiment conducted in a rotating annular flume by Villaret and Paulic (1986). The use of the layer-by-layer procedure based on a least squares fitting technique provides a closer fit than the global procedure. The present study points out the complementarity of experimental and numerical approaches and also suggests possible improvements in laboratory test procedures.     

Timing of the late-glacial climate reversal in the Southern Hemisphere using high-resolution radiocarbon chronology for Kaipo bog, New Zealand

The pattern of climate change in the Southern Hemisphere during the Younger Dryas (YD) chronozone provides essential constraint on mechanisms of abrupt climate change only if accurate, high-precision chronologies are obtained. A climate reversal reported previously at Kaipo bog, New Zealand, had been dated between 13,600 and 12,600 cal yr B.P. and appeared to asynchronously overlap the YD chron, but the chronology, based on conventionally radiocarbon-dated bulk sediment samples, left the precise timing questionable. We report a new high-resolution AMS ¹⁴C chronology for the Kaipo record that confirms the original chronology and provides further evidence for a mid-latitude Southern Ocean cooling event dated between 13,800 and 12,400 cal yr B.P. (2σ range), roughly equivalent to the Antarctic Cold Reversal.     

Planation, bauxites and epeirogeny: One or two paleosurfaces on the West African margin?

Mapping of lateritic bauxites over the West African rifted margin and analysis of the geomorphic properties of these bauxites, combined with available geological data lead to a discussion of the presence of either two Meso-Cenozoic planation surfaces or a single Eocene surface to account for the morphotectonic and paleoclimatic evolution of the Guinean landforms. At large scale, two stepped bauxitic levels are documented. Ongoing or episodic uplift following Gondwana breakup and Meso-Cenozoic climate change are proposed to have allowed the formation and abandonment of an Early mid-Cretaceous surface today preserved as the higher bauxitic level, and the setting of an Eocene planation surface bearing a second generation of bauxites, making the lower bauxitic level. The single Eocene surface hypothesis requires that Paleogene bauxitization preserved large pre-existing relief to explain two stepped bauxitic levels of the same age. The two-surface hypothesis is favored because it would explain rebauxitization of alluvial pebbles of bauxites under the lower lateritic level.     

3D numerical modelling of fluid flow in the Val-d’Or orogenic gold district: major crustal shear zones drain fluids from overpressured vein fields

Fluid flow patterns have been determined using oxygen isotope isopleths in the Val-d’Or orogenic gold district. 3D numerical modelling of fluid flow and oxygen isotope exchange in the vein field shows that the fluid flow patterns can be reproduced if the lower boundary of the model is permeable, which represents middle or lower crustal rocks that are infiltrated by a metamorphic fluid generated at deeper levels. This boundary condition implies that the major crustal faults so conspicuous in vein fields do not act as the only major channel for upward fluid flow. The upper model boundary is impermeable except along the trace of major crustal faults where fluids are allowed to drain out of the vein field. This upper impermeable boundary condition represents a low-permeability layer in the crust that separates the overpressured fluid from the overlying hydrostatic fluid pressure regime. We propose that the role of major crustal faults in overpressured vein fields, independent of tectonic setting, is to drain hydrothermal fluids out of the vein field along a breach across an impermeable layer higher in the crust and above the vein field. This breach is crucial to allow flow out of the vein field and accumulation of metals in the fractures, and this breach has major implications for exploration for mineral resources. We propose that tectonic events that cause episodic metamorphic dehydration create a short-lived pulse of metamorphic fluid to rise along zones of transient permeability. This results in a fluid wave that propagates upward carrying metals to the mineralized area. Earthquakes along crustal shear zones cause dilation near jogs that draw fluids and deposit metals in an interconnected network of subsidiary shear zones. Fluid flow is arrested by an impermeable barrier separating the hydrostatic and lithostatic fluid pressure regimes. Fluids flow through the evolving and interconnected network of shear zones and by advection through the rock matrix. Episodic breaches in the impermeable barrier along the crustal shear zones allow fluid flow out of the vein field.     

Avoiding pitfalls of least‐squares inversion by using the energy‐flux error criterion

Seismic inversion by modelling and data fitting depends on the criterion chosen to measure the misfit between observed and modelled data. The popular least‐squares error criterion has an important drawback: it is sensitive both to the shape of the recording surface and to velocity variations along this surface. Tests on synthetic seismic reflection data show that least‐squares inversion may work surprisingly poorly in situations where (i) the range of angles between reflected rays and the acquisition surface is large, (ii) the velocity varies significantly along this surface, or (iii) a compensation for the effects of dissipation is applied to the gradients. In these situations, the gradients may contain important artefacts and have incorrect amplitudes. The outgoing flux of energy of the residual wavefield across the acquisition surface provides an alternative measure of the data misfit which is independent of the recording surface, provided this surface is closed, and which is only sensitive to the aperture in the practical situation of an open surface or line of receivers. Energy‐flux inversion presents a strong resemblance to reverse‐time migration, but with the additional possibility of iteratively improving the images. In all the tests, energy‐flux inversion provided better images than least‐squares inversion.     

Investigating the Halogen Chemistry From High-Latitude Nighttime Stratospheric Measurements of OClO and NO2

Two cases of simultaneous nighttime measurements of NO₂ and OClO in the winter polar stratosphere are analyzed in order to test our present knowledge of halogen chemistry in the presence of high amount of NO₂ at low temperature. Comparisons with Lagrangian model calculations using several hypotheses are performed. First simulations, using the admitted constant rates of chemical reaction, strongly underestimate the measured OClO while the NO₂ profiles are correctly reproduced. If uncertainties in actinic fluxes calculations are taken into account, simulation results do not show a significant reduction of the underestimation. A better agreement can be achieved if the formation of unstable isomers of ClONO₂ and of BrONO₂ occurs in the cold conditions of the polar stratosphere. An approximate value of the branching ratios of the channels leading to ClONO₂ and ClOONO, and to BrONO₂ and BrOONO, necessary to reproduce both OClO and NO₂ is given and discussed.     

Weathering: Toward a Fractal Quantifying

Weathering occurs over a wide range of scales. To link features through these scales is a major challenge for interdisciplinary weathering studies. Fractal approach seems to be specially useful for this purpose. We introduce a multistep fractal weathering assessment scheme devoted to extract fractal weathering classifiers from texture analysis of the mineral's image. Our scheme enables to quantitatively estimate the global and local information about the geometry of the weathering pattern. This information is basic to develop geometrical indices of weathering, which can significantly enrich the common qualitative and semiquantitative weathering assessment schemes. To justify the fractal approach, a strong statistical self-similarity has been documented for both the weathering and fresh features of two common silica minerals: quartz and biogenic A-opal (phytolith) over four orders of length scales. The procedure is fast, drastically reduces thresholding bias, promises to be universal, it is valid for genetically different minerals and rock types, scale independent, and specially useful for monitoring the changes in the mineral's roughness during the alteration. Two of the proposed classifiers seem to be potentially useful for direct application in the field and be used by nonspecialist.     

Weights of Evidence Modeling of Mineral Potential: A Case Study Using Small Number of Prospects,Abra, Philippines

Weights of evidence (WofE) modeling usually is applied to map mineral potential in areas with large number of deposits/prospects. In this paper, WofE modeling is applied to a case study area measuring about 920 km² with 12 known porphyry copper prospects. A pixel size of 100 m × 100 m was used in the spatial data analyses to represent in a raster-based GIS lateral extents of prospects and of geological features considered as spatial evidence. Predictor maps were created based on (a) estimates of studentized values of positive spatial association between prospects and spatial evidence; (b) proportion of number of prospects in zones where spatial evidence is present; and (c) geological interpretations of positive spatial association between prospects and spatial evidence. Uncertainty because of missing geochemical evidence is shown to have an influence on tests of assumption of conditional independence (CI) among predictor maps with respect to prospects. For the final predictive model, assumption of CI is rejected based on omnibus test but is accepted based on a new omnibus test. The final predictive model, which delineates 30% of study area as zones with potential for porphyry copper, has 83% success rate and 73% prediction rate. The results demonstrate plausibility of WofE modeling of mineral potential in large areas with small number of mineral prospects.     

Rainfall thresholds for landsliding in the Himalayas of Nepal

Landsliding of the hillslope regolith is an important source of sediment to the fluvial network in the unglaciated portions of the Himalayas of Nepal. These landslides can produce abrupt increases of up to three orders of magnitude in the fluvial sediment load in less than a day. An analysis of 3 years of daily sediment load and daily rainfall data defines a relationship between monsoonal rainfall and the triggering of landslides in the Annapurna region of Nepal. Two distinct rainfall thresholds, a seasonal accumulation and a daily total, must be overcome before landslides are initiated. To explore the geomorphological controls on these thresholds, we develop a slope stability model, driven by daily rainfall data, which accounts for changes in regolith moisture. The pattern of rainfall thresholds predicted by the model is similar to the field data, including the decrease in the daily rainfall threshold as the seasonal rainfall accumulation increases. Results from the model suggest that, for a given hillslope, regolith thickness determines the seasonal rainfall necessary for failure, whereas slope angle controls the daily rainfall required for failure.