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.     

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.     

Triggering Mechanisms of Slope Instability and their Relationship to Earthquakes and Tsunamis

— Submarine and shoreline slope failures that accompany large earthquakes and large tsunamis are triggered by several mechanisms. Triggering mechanisms range from direct effects, such as inertial forces from earthquake shaking, to indirect effects, such as rapid drawdown that occurs when an earthquake-generated tsunami first approaches a shoreline. Soil shear strength also plays an important role in earthquake-related slope failures. Earthquakes change the shear strength of the soil by inducing excess pore water pressures. These excess pore water pressures change with time after the earthquake, resulting in changes in shear strength and slope stability with time. This paper reviews earthquake-related triggering mechanisms for submarine and shoreline slope failures. The variation in shear strength with time following an earthquake is examined and it is shown that delayed slope failures after an earthquake can occur as a result of changes in earthquake-induced excess pore water pressures and shear strength with time.     

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.     

Pn arrivals and lateral variations of Moho geometry beneath the Kaapvaal craton

Pn arrivals from mining-induced earthquakes on the edge of the Witwatersrand basin show that the P wavespeeds in the uppermost mantle are almost constant throughout most of the Kaapvaal craton. The presence of only small wavespeed variations allows the use of a simple method of estimating crustal thicknesses below the stations of the Kaapvaal broad-band network using Pn times that has been compared with results from receiver functions. One thousand three hundred thirty-seven Pn arrivals were used to derive crustal thicknesses at 46 stations on the Kaapvaal craton. The average crustal thicknesses for 19 centrally located stations on each of the northern and southern regions of the craton that yielded well-constrained thicknesses were 50.52±0.88 km and 38.07±0.85 km, respectively. In contrast, the corresponding average thicknesses determined from receiver functions were 43.58±0.57 km and 37.58±0.70 km, respectively. The systematically lower values for receiver functions in the northern part of the Kaapvaal craton that was affected by the Bushveld magmatism at 2.05 Ga, suggest that the receiver functions do not enable the petrological crust mantle boundary to be reliably resolved due to variations in composition and metamorphic grade in a mafic lower crust. The Pn times also suggest pervasive azimuthal anisotropy with maximum wavespeeds of about 8.40 km/s at azimuths of about 15° and 217° in the northern and southern regions of the craton, respectively, and minimum wavespeeds of about 8.25 km/s.