Shape and distribution analysis of Merensky Reef potholing, Northam Platinum Mine, western Bushveld Complex: implications for pothole formation and growth

Syn-magmatic removal of the cumulate pile during the formation of the Bushveld Complex resulted in “potholes”. Erosion progressed downward in the cumulate pile, resulting in a series of steep, transgressive contacts between locally conformable potholed reefs in the regional pothole sub-facies of the Swartklip Facies in the western limb of the Bushveld Complex. The deepest of these potholes, “third-order” or “FWP2” potholing, occurs where the base of the Merensky Cyclic Unit transgresses the Upper Pseudo-Reef Chromitite marker horizon. The base of a FWP2 pothole on Northam Platinum Mine consists of an unconformable stringer Merensky Chromitite overlain by a medium-grained, poikilitic orthopyroxenite and underlain by either a pegmatitic harzburgite or the medium-grained Lower Pseudo-Reef Anorthosite. Detailed shape and distribution analysis of FWP2 potholes reveals underlying patterns in their shape and distribution which, in turn, suggest a structural control. The ratio between pothole short vs long axes is 0.624 (N=1,385), although the ratio increases from 0.48 to 0.61 in the long axis range 10 to 60 m, then decreases from 0.61 to 0.57 from 61 to 100 m, increasing again from 0.57 to 0.61 from 101 to 400 m, suggesting that there is not a simple relationship between pothole shape and size. Shape (circularity, eccentricity, and dendricity) analysis of a subset of 638 potholes indicates that potholes with long axes <100 m have an elliptical, average normalized shape, elongate on a 120–150° orientation. Potholes with long axis lengths >100 m have an average normalized shape that is bilobate and elongate on a 120° orientation. The average aspect ratio (short axis length divided by long axis length) of potholes is highest for potholes with long axis lengths >100 m and lowest for potholes with long axis lengths between 35 and 60 m. The most common long axis orientation for potholes with long axis lengths <100 m is 150° but 120° for long axis lengths >100 m. Fractal analysis indicates that the distribution of pothole centers is controlled neither by a single nor several interacting fractal dimensions. Autocorrelation (Fry) analysis of the distribution of pothole centers shows recurring pothole distribution trends at 038, 070, and 110° for potholes over the full range of long axis lengths, while the trends of 008 and 152° occur in potholes with long axes lengths between 60 and 100 m. Chi-squared (X ²) analysis of the locations of pothole centers suggests that the distribution of small potholes is highly non-uniform but becomes exponentially more uniform with increasing pothole size. The model which best fits the observed shape and distribution analysis is a combination of protracted independent growth and “nearest neighbor” merging along specific orientations. For instance, the clustered distribution of original pothole centers resulted in merged potholes with long axes lengths of up to 60 m, exhibiting short vs long axes ratios of 0.61, preferred orientations of 150°, and alignment along 010 and 150° trends. Further independent growth allowed for merging of similar-sized (and smaller) neighboring potholes, generating potholes with long axes of up to 100 m in length, a preferred long axis orientation of 150°, and alignment along 010, 040, 075, and 150°. Subsequent preferential merging occurred along a 120° trend, thereby preserving a bilobate form. This implies that while pothole initiation and enlargement may be driven by a “top-down” (i.e., possibly thermomechanical) process, an underlying linear or structural catalyst/control is revealed in changes in pothole shape during enlargement and, furthermore, in the preferred trends along which potholes merged over a considerable period, possibly concomitant with adjustment of major structures in the footwall to the Bushveld Complex and pulses into the magma chamber.     

Mixing and differentiation in the Oruanui rhyolitic magma, Taupo, New Zealand: evidence from volatiles and trace elements in melt inclusions

Large pyroclastic rhyolites are snapshots of evolving magma bodies, and preserved in their eruptive pyroclasts is a record of evolution up to the time of eruption. Here we focus on the conditions and processes in the Oruanui magma that erupted at 26.5 ka from Taupo Volcano, New Zealand. The 530 km³ (void-free) of material erupted in the Oruanui event is comparable in size to the Bishop Tuff in California, but differs in that rhyolitic pumice and glass compositions, although variable, did not change systematically with eruption order. We measured the concentrations of H₂O, CO₂ and major and trace elements in zoned phenocrysts and melt inclusions from individual pumice clasts covering the range from early to late erupted units. We also used cathodoluminescence imaging to infer growth histories of quartz phenocrysts. For quartz-hosted inclusions, we studied both fully enclosed melt inclusions and reentrants (connecting to host melt through a small opening). The textures and compositions of inclusions and phenocrysts reflect complex pre-eruptive processes of incomplete assimilation/partial melting, crystallization differentiation, magma mixing and gas saturation. ‘Restitic’ quartz occurs in seven of eight pumice clasts studied. Variations in dissolved H₂O and CO₂ in quartz-hosted melt inclusions reflect gas saturation in the Oruanui magma and crystallization depths of ∼3.5–7 km. Based on variations of dissolved H₂O and CO₂ in reentrants, the amount of exsolved gas at the beginning of eruption increased with depth, corresponding to decreasing density with depth. Pre-eruptive mixing of magma with varying gas content implies variations in magma bulk density that would have driven convective mixing. Electronic Supplementary Material Supplementary material is available for this article at and is accessible for authorized users.     

Near-Wall Flow Development After A Step Change In Surface Roughness

An experimental study of the initial flow field downstream of a step change in surface roughness is presented. The roughness length of the downstream surface was approximately tenfold that of the upstream roughness and, unlike all previous studies, attention was concentrated on the roughness sublayer region beneath the inertial (log-law) region. The experiments were conducted at a boundary layer Reynolds number of about 6 × 10⁴ (based on layer thickness andfree-stream velocity) and around a longitudinal location where the (downstream) roughness length, z_o2, was about 1% of the boundary-layer thickness atthe roughness change point.The thickness of the roughness sublayer was found for the two roughness. It was observed that the vertical profiles of mean velocity and turbulence characteristics started to show similarity after about 160z₀₂ downstream of the roughness change. The presence of a shear stress overshoot is shown to depend strongly on the precise location (with respect to the roughness elements) at which the measurements are made and the thickness of the equilibrium layer is shown to be very sensitive to the way it is defined. It is demonstrated that the growing equilibrium layer has first to encompass the roughness sublayer before any thickness of inertial sublayer can be developed. It follows that, in somepractical cases, like flows across some urban environments, the latter(log-law) region may never exist at all.     

In situ measurement of Re-Os isotopes in mantle sulfides by laser ablation multicollector-inductively coupled plasma mass spectrometry: analytical methods and preliminary results

A method has been developed for the in situ determination of Re-Os isotopes in single grains of sulfides in mantle-derived peridotites using a laser ablation microprobe attached to a multicollector-induced coupled plasma mass spectrometer (MC-ICPMS). High-precision Os isotope analysis by MC-ICPMS is demonstrated by the measurement of interlaboratory Os standards. Evaluation of mass bias correction procedures shows that the exponential law provides the best fit to the Os isotope data and that the ratio of the mass fractionation coefficients for Re and Os remains constant for the range of typical instrument operating conditions. This relationship enables the accurate and precise correction of the isobaric interference of ¹⁸⁷Re on ¹⁸⁷Os for ¹⁸⁷Re/¹⁸⁸Os values up to 1.6.Results are presented for single sulfide inclusions in olivine macrocryts from kimberlites in the Siberian and Slave Cratons, and sulfides enclosed in silicates and interstitial to silicates in peridotite xenoliths from the Slave Craton and Massif Central, France. Enclosed sulfides larger than 50 μm in diameter and with Os contents ≥40 ppm give ¹⁸⁷Os/¹⁸⁸Os ratios with a precision of 0.1% (2 SE), which is equivalent to N-TIMS whole-rock data. Interstitial sulfides typically have lower Os (10 to 30 ppm) and give analyses with lower precision (∼1 to 2%) but still provide valuable information on the movement of Os within the lithosphere. The sulfide inclusions in silicates preserve significantly less radiogenic Os isotopic compositions than interstitial sulfides and accordingly produce significantly older and more realistic Re-Os age information. Whole-rock Os isotope compositions reflect the proportions of different generations of enclosed and interstitial sulfides; this calls into question the significance of many published “depletion ages.”     

Implications of Nb/U, Th/U and Sm/Nd in plume magmas for the relationship between continental and oceanic crust formation and the development of the depleted mantle

The Nb/U and Th/U of the primitive mantle are 34 and 4.04 respectively, which compare with 9.7 and 3.96 for the continental crust. Extraction of continental crust from the mantle therefore has a profound influence on its Nb/U but little influence on its Th/U. Conversely, extraction of midocean ridge-type basalts lowers the Th/U of the mantle residue but has little influence on its Nb/U. As a consequence, variations in Th/U and Nb/U with Sm/Nd can be used to evaluate the relative importance of continental and basaltic crust extraction in the formation of the depleted (Sm/Nd enriched) mantle reservoir.This study evaluates Nb/U, Th/U, and Sm/Nd variations in suites of komatiites, picrites, and their associated basalts, of various ages, to determine whether basalt and/or continental crust have been extracted from their source region. Emphasis is placed on komatiites and picrites because they formed at high degrees of partial melting and are expected to have Nb/U, Th/U, and Sm/Nd that are essentially the same as the mantle that melted to produce them. The results show that all of the studied suites, with the exception of the Barberton, have had both continental crust and basaltic crust extracted from their mantle source region. The high Sm/Nd of the Gorgona and Munro komatiites require the elevated ratios seen in these suites to be due primarily to extraction of basaltic crust from their source regions, whereas basaltic and continental crust extraction are of subequal importance in the source regions of the Yilgarn and Belingwe komatiites. The Sm/Nd of modern midocean ridge basalts lies above the crustal extraction curve on a plot of Sm/Nd against Nb/U, which requires the upper mantle to have had both basaltic and continental crust extracted from it.It is suggested that the extraction of the basaltic reservoir from the mantle occurs at midocean ridges and that the basaltic crust, together with its complementary depleted mantle residue, is subducted to the core-mantle boundary. When the two components reach thermal equilibrium with their surroundings, the lighter depleted component separates from the denser basaltic component. Both are eventually returned to the upper mantle, but the lighter depleted component has a shorter residence time in the lower mantle than the denser basaltic component. If the difference in the recycling times for the basaltic and depleted components is ∼1.0 to 1.5 Ga, a basaltic reservoir is created in the lower mantle, equivalent to the amount of basalt that is subducted in 1.0 to 1.5 Ga, and that reservoir is isolated from the upper mantle. It is this reservoir that is responsible for the Sm/Nd ratio of the upper mantle lying above the trend predicted by extraction of continental crust on the plot of Sm/Nd against Nb/U.     

Calcrete profile development in Quaternary alluvial sequences,southeast Spain: implications for using calcretes as a basis for landform chronologies

A detailed study of the morphology and micro‐morphology of Quaternary alluvial calcrete profiles from the Sorbas Basin shows that calcretes may be morphologically simple or complex. The ‘simple’ profiles reflect pedogenesis occurring after alluvial terrace formation and consist of a single pedogenic horizon near the land surface. The ‘complex’ profiles reflect the occurrence of multiple calcrete events during terrace sediment aggradation and further periods of pedogenesis after terrace formation. These ‘complex’ calcrete profiles are consequently described as composite profiles. The exact morphology of the composite profiles depends upon: (1) the number of calcrete‐forming events occurring during terrace sediment aggradation; (2) the amount of sediment accretion that occurs between each period of calcrete formation; and (3) the degree of pedogenesis after terrace formation. Simple calcrete profiles are most useful in establishing landform chronologies because they represent a single phase of pedogenesis after terrace formation. Composite profiles are more problematic. Pedogenic calcretes that form within them may inherit carbonate from calcrete horizons occurring lower down in the terrace sediments. In addition erosion may lead to the exhumation of older calcretes within the terrace sediment. Calcrete ‘inheritance’ may make pedogenic horizons appear more mature than they actually are and produce horizons containing carbonate embracing a range of ages. Calcrete exhumation exposes calcrete horizons whose morphology and radiometric ages are wholly unrelated to terrace surface age. Composite profiles are, therefore, only suitable for chronological studies if the pedogenic horizon capping the terrace sequence can be clearly distinguished from earlier calcrete‐forming events. Thus, a detailed morphological/micro‐morphological study is required before any chronological study is undertaken. This is the only way to establish whether particular calcrete profiles are suitable for dating purposes. Copyright © 2003 John Wiley & Sons, Ltd.     

Simulation and measurement of surface shear stress over isolated and closely spaced transverse dunes in a wind tunnel

Topographic interactions generate multidirectional and unsteady air?ow that limits the application of velocity pro?le approaches for estimating sediment transport over dunes. Results are presented from a series of wind tunnel simulations using Irwin‐type surface‐mounted pressure sensors to measure shear stress variability directly at the surface over both isolated and closely spaced sharp‐crested model dunes. Findings complement existing theories on secondary air?ow effects on stoss transport dynamics and provide new information on the in?uence of lee‐side air?ow patterns on dune morphodynamics. For all speeds investigated, turbulent unsteadiness at the dune toe indicates a greater, more variable surface shear, despite a signi?cant drop in time‐averaged measurements of streamwise shear stress at this location. This effect is believed suf?cient to inhibit sediment deposition at the toe and may be responsible for documented intermittency in sand transport in the toe region. On the stoss slope, streamline compression and ?ow acceleration cause an increase in ?ow steadiness and shear stress to a maximum at the crest that is double that at the toe of the isolated dune and 60–70 per cent greater than at ?ow reattachment on the lower stoss of closely spaced dunes. Streamwise ?ow accelerations, rather than turbulence, have greater in?uence on stress generation on the stoss and this effect increases with stoss slope distance and with incident wind speed. Reversed ?ow within the separation cell generates signi?cant surface shear (30–40 per cent of maximum values) for both spacings. This supports ?eld studies that suggest reversed ?ow is competent enough to return sediment to the dune directly or in a de?ected direction. High variability in shear at reattachment indicates impact of a turbulent shear layer that, despite low values of time‐averaged streamwise stress in this region, would inhibit sediment accumulation. Downwind of reattachment, shear stress and ?ow steadiness increase within 6 h (h = dune height) of reattachment and approach upwind values by 25 h. A distance of at least 30 h is suggested for full boundary layer recovery, which is comparable to ?uvial estimates. The Irwin sensor used in this study provides a reliable means to measure skin friction force responsible for sand transport and its robust, simple, and cost‐effective design shows promise for validating these ?ndings in natural dune settings. Copyright © 2003 John Wiley & Sons, Ltd.     

The stability of a remediated bed in Hamilton Harbour,Lake Ontario,Canada

ABSTRACT In situ measurements of lakebed sediment erodibility were made on three sites in Hamilton Harbour, Lake Ontario, using the benthic flume Sea Carousel. Three methods of estimating the surface erosion threshold (τ_c(0)) from a Carousel time series were evaluated: the first method fits measures of bed strength to eroded depth (the failure envelope) and evaluates threshold as the surface intercept; the second method regresses mean erosion rate (E_m) with bed shear stress and solves for the floc erosion rate (E_f) to derive the threshold for E_m = E_f = 1 × 10^?5 kg m^?2 s^?1; the third method extrapolates a regression of suspended sediment concentration (S) and fluid transmitted bed shear stress (τ₀) to ambient concentrations. The first field site was undisturbed (C) and acted as a control; the second (W) was disturbed through ploughing and water injection as part of lakebed treatment, whereas the third site (OIP) was disturbed and injected with an oxidant used for remediation of contaminated sediment. The main objectives of this study were: (1) to evaluate the three different methods of deriving erosion threshold; (2) to compare the physical behaviour of lacustrine sediments with their marine estuarine counterparts; and (3) to examine the effects of ploughing and chemical treatment of contaminated sediment on bed stability. Five deployments of Sea Carousel were carried out at the control site. Mean erosion thresholds for the three methods were: τ_c(0) = 0·5 (±0·06), 0·27 (±0·01) and 0·34 (±0·03) Pa respectively. Method 1 overpredicted bed strength as it was insensitive to effects in the surface 1–2 mm, and the fit of the failure envelope was also highly subjective. Method 2 exhibited a wide scatter in the data (low correlation coefficients), and definition of the baseline erosion rate (E_f) is largely arbitrary in the literature. Method 3 yielded stable (high correlation coefficients), reproducible and objective results and is thus recommended for evaluation of the erosion threshold. The results of this method correlated well with sediment bulk density and followed the same trend as marine counterparts from widely varying sites. Mass settling rates, expressed as a decay constant, k, of S(t), were strongly related to the maximum turbidity at the onset of settling (S_max) and were also in continuity with marine counterparts. Thus, it appears that differences in salinity had little effect on mass settling rates in the examples presented, and that biological activity dominated any effects normally attributable to changes in salinity. Bedload transport of eroded aggregates (2–4 mm in diameter) took place by rolling below a mean tangential flow velocity (U_y) of 0·32 ms^?1 and by saltation at higher velocities. Mass transport as bedload was a maximum at U_y = 0·4 ms^?1, although bedload never exceeded 1% of the suspended load. The proportion of material moving as bedload was greatest at the onset of erosion but decreased as flow competence increased. Given the low bulk density and strength of the lakebed sediment, the presence of a bedload component is notable. Bedload transport over eroding cohesive substrates should be greater in estuaries, where both sediment density and strength are usually higher. Significant differences between the ploughed and control sites were apparent in both the erosion rate and the friction coefficient (φ), and suggest that bed recovery after disruption is rapid (< 24 h). τ_c(0) increased linearly with time after ploughing and recovered to the control mean value within 3 days. The friction coefficient was reduced to zero by ploughing (diagnostic of fluidization), but increased linearly with time, regaining control values within 6 days. No long‐term reduction in bed strength due to remediation was apparent.