A source-Sink laboratory model of the ocean circulation

Abstract

A two gyre circulation and inertial western boundary currents have been observed in a sloping bottom laboratory model of a barotropic ocean circulation. Water of viscosity v is contained in a rotating (angular velocity ω), square basin of side L (30 cm) with a flat top and a bottom slope (tan θ) such that the depth (H) varies from 12 to 15 cm. The flow is driven by a distributed source and sink at the upper surface, a plate drilled with 342 holes. The hole distribution and size is arranged so that the average imposed vertical velocity, w = w ₀ sin (2πy′/30), models the Ekman divergence from a two gyre zonal wind stress. Fluid flow is observed with the thymol blue technique over the ranges of Rossby numbers (w ₀/2ωL tan θ) from 1.44 × 10^?3 to 1.41 × 10^?2 and Ekman numbers (v/2ωH ²) from 2.13 × 10^?5 to 2.10 × 10^?3. At the largest Rossby numbers the flow pattern changes markedly, but the non-uniformity of the imposed vertical velocity also penetrates deep into the fluid in this regime.     

Stability of buoyancy opposed mixed convection in a vertical channel and its dependence on permeability

Non-Darcy mixed convective flow of water due to external pressure gradient and buoyancy opposed forces are considered in a vertical channel filled with porous medium, which can be either isotropic or anisotropic. The linear theory of stability analysis has been used to numerically investigate the dependence of the transition behavior of the fully developed basic flow on the permeability of the medium. Numerical experiments indicate that mainly two main instability modes appear: Rayleigh–Taylor (R–T) and buoyant instability. For Darcy numbers (Da) ?10⁻⁹, R–T instability dominates within the entire Reynolds number (Re) range considered here. It was also found that for the same Re, the fully developed base flow is highly unstable (stable) for porous media with high (low) permeability. Further, it was seen that the disturbance isotherm cells migrate from the channel walls toward the centerline when permeability is reduced. Reducing the permeability by one order of magnitude (corresponding to a decrease of Darcy number from 10⁻⁶ to 10⁻⁷) increases base flow stability approximately 20-fold. For higher Reynolds numbers, buoyant, mixed and shear instability of the basic flow were found when Da was increased from 10⁻⁷ to 10⁻³. However, for cases in which permeability and porosity behaved as suggested by Carman–Kozeny relation (CKR), buoyant stability was the only mode of instability. Critical values of the Rayleigh (Ra) and Darcy (Da) numbers in the R–T mode of instability were related to each other by the hyperbolic function RaDa = −2.465.     

Debris-flow-dominated sediment transport through a channel network after wildfire

Field studies that investigate sediment transport between debris-flow-producing headwaters and rivers are uncommon, particularly in forested settings, where debris flows are infrequent and opportunities for collecting data are limited. This study quantifies the volume and composition of sediment deposited in the arterial channel network of a 14-km² catchment (Washington Creek) that connects small, burned and debris-flow-producing headwaters (<1 km²) with the Ovens River in SE Australia. We construct a sediment budget by combining new data on deposition with a sediment delivery model for post-fire debris flows. Data on deposits were plotted alongside the slope–area curve to examine links between processes, catchment morphometry and geomorphic process domains. The results show that large deposits are concentrated in the proximity of three major channel junctions, which correspond to breaks in channel slope. Hyperconcentrated flows are more prominent towards the catchment outlet, where the slope–area curve indicates a transition from debris flow to fluvial domains. This shift corresponds to a change in efficiency of the flow, determined from the ratio of median grain size to channel slope. Our sediment budget suggests a total sediment efflux from Washington Creek catchment of 61 × 10³ m³. There are similar contributions from hillslopes (43 ± 14 × 10³ m³), first to third stream order channel (35 ± 12 × 10³ m³) and the arterial fourth to fifth stream order channel (31 ± 17 × 10³ m³) to the total volume of erosion. Deposition (39 ± 17 × 10³ m³) within the arterial channel was higher than erosion (31 ± 17 × 10³ m³), which means a net sediment gain of about 8 × 10³ m³ in the arterial channel. The ratio of total deposition to total erosion was 0.44. For fines <63 μm, this ratio was much smaller (0.11), which means that fines are preferentially exported. This has important implications for suspended sediment and water quality in downstream rivers. © 2019 John Wiley & Sons, Ltd.     

Crustal neon: a striking uniformity

By combining data from a diverse suite of crustal fluid samples representing a broad geographical distribution, we have identified a well-defined nucleogenic (crustal) neon component. The neon is produced from (α, n) and (n, α) nuclear interactions involving nuclei of O, Mg, and F [1]. In the limiting case of ²⁰Ne/²²Ne = 0, the composition is: ²¹Ne/²²Ne = 0.47 ± 0.01 and ²¹Ne/⁴He = (0.46 ± 0.08) × 10⁻⁷. A crustal O/F ratio of 110 (atomic) calculated from the ²¹Ne/²²Ne ratio is 4–10 times less than the average crustal O/F ratio. The discrepancy can be accounted for by an enhanced O/F ratio within the 10–40 μm range of the U-Th-generated α-particles.     

PRINCIPES NOUVEAUX POUR DISTINGUER CE QUI EST Dû A L'ÉCOULEMENT SOUTERRAIN SUR L'HYDROGRAPHE DES FLEUVES

Abstract

Throughflow has been measured from three soil horizons on a 12 slope with impermeable, bedrock. Storm flow comes from the 10–45 cm horizon and is controlled by the upslope extent of saturated conditions. Base flow comes from the 45–75 cm horizon and is supplied by slow unsaturated flow from the whole soil mass to a small constant zone of saturation.

Differences between input and output stream hydrographs over 270 metres of channel are attributed to throughflow and correlate well with measured values providing a basis for separating throughflow components from the stream hydrograph. Observed stream flows contain no overland flow or ground water flow components. The main basin flood peak is not generated within this control section of channel but is produced in the headwater zone (0.1 km²) by the faster runoff characteristics of the soils in that area and by topographic factors which lead to rapid channel extension.     

Evolution of the spatial and temporal characteristics of the isotope hydrology of a montane river basin

Abstract

Precipitation and streamwater were analysed weekly for δ¹⁸O in seven tributaries and five main stem sites of a 2100 km² catchment; >60% of it is upland in character. Precipitation δ¹⁸O followed seasonal patterns ranging from –20‰ in winter to –4‰ in summer. Seasonality was also evident in stream waters, though much more damped. Mean transit times (MTTs) were estimated using δ¹⁸O input–output relationships in a convolution integral with a gamma distribution. The MTTs were relatively similar (528–830 days): tributaries exhibited a greater range, being shorter in catchments with montane topography and hydrologically responsive soils, and longer where catchments have significant water storage. Along the main stem, MTTs increased modestly from 621 to 741 days. This indicates that montane headwaters are the dominant sources of runoff along the main stem of the river system. Models suggest that around 10% of runoff has transit times of less than two weeks during higher flows whilst older (>10-year old) water sustains low flows contributing around 5% of runoff.

Citation Speed, M., Tetzlaff, D., Hrachowitz, M. & Soulsby, C. (2011) Evolution of the spatial and temporal characteristics of the isotope hydrology of a montane river basin. Hydrol. Sci. J. 56(3), 426–442     

The influence of flow discharge variations on the morphodynamics of a diffluence–confluence unit on a large river

Bifurcations are key geomorphological nodes in anabranching and braided fluvial channels, controlling local bed morphology, the routing of sediment and water, and ultimately defining the stability of their associated diffluence–confluence unit. Recently, numerical modelling of bifurcations has focused on the relationship between flow conditions and the partitioning of sediment between the bifurcate channels. Herein, we report on field observations spanning September 2013 to July 2014 of the three‐dimensional flow structure, bed morphological change and partitioning of both flow discharge and suspended sediment through a large diffluence–confluence unit on the Mekong River, Cambodia, across a range of flow stages (from 13 500 to 27 000 m³ s^?1). Analysis of discharge and sediment load throughout the diffluence–confluence unit reveals that during the highest flows (Q = 27 000 m³ s^?1), the downstream island complex is a net sink of sediment (losing 2600 ± 2000 kg s^?1 between the diffluence and confluence), whereas during the rising limb (Q = 19 500 m³ s^?1) and falling limb flows (Q = 13 500 m³ s^?1) the sediment balance is in quasi‐equilibrium. We show that the discharge asymmetry of the bifurcation varies with discharge and highlight that the influence of upstream curvature‐induced water surface slope and bed morphological change may be first‐order controls on bifurcation configuration. Comparison of our field data to existing bifurcation stability diagrams reveals that during lower (rising and falling limb) flow the bifurcation may be classified as unstable, yet transitions to a stable condition at high flows. However, over the long term (1959–2013) aerial imagery reveals the diffluence–confluence unit to be fairly stable. We propose, therefore, that the long‐term stability of the bifurcation, as well as the larger channel planform and morphology of the diffluence–confluence unit, may be controlled by the dominant sediment transport regime of the system. © 2017 The Authors. Earth Surface Processes and Landforms published by John Wiley & Sons Ltd.     

Flow field-induced drag forces and swimming behavior of three benthic fish species

Modern ethohydraulics is the study of the behavioral responses of swimming fish to flow fields. However, the exact drag forces experienced by fish remain poorly studied; this information is required to obtain a better understanding of the behavioral responses of fish and their current resistance strategies. We measured near-ground frontal drag forces on preserved individuals of three benthic fish species, round goby (Neogobius melanstomus), gudgeon (Gobio gobio) and bullhead (Cottus gobio), in a flow channel. The forces were compared to acoustic Doppler velocity (ADV) measurements and fish tracking data based on video observations of live fish in the flow channel. Overall, we observed drag coefficients (C_D) of ∼10⁻³ at Reynolds numbers ∼10⁵. The frontal drag forces acting on preserved fish with non-spread fins ranged from -1.96 mN*g^-1 (force per fish wet weight, velocity 0.55 m*s^-1) to 11.01 mN*g^-1 (velocity 0.85 m*s^-1). Spreading the fins strongly increased the drag forces for bullhead and round goby. In contrast, the drag forces were similar for gudgeon with spread fins and all fish with non-spread fins. Video tracking revealed no clear relationship between the position of the fish in the flow field and the forces experienced by the preserved fish at these positions. Collectively, these results suggest that i) the differences in frontal drag forces between species are small in homogenous flow, ii) individuals chose their position in the flow field based on factors other than the drag forces experienced, and iii) whether fins are spread or non-spread is an essential quality that modulates species-specific differences. The methodology and results of this study will enable integration of flow measurements, fish behavior and force measurements and inform ethohydraulics research. More advanced force measurements will lead to a detailed understanding of the current resistance strategies of benthic fish and improve the design of fish passes.     

Construction dynamics of a lava channel

We use a kinematic GPS and laser range finder survey of a 200 m-long section of the Muliwai a Pele lava channel (Mauna Ulu, Kilauea) to examine the construction processes and flow dynamics responsible for the channel–levee structure. The levees comprise three packages. The basal package comprises an 80–150 m wide ′a′a flow in which a ∼2 m deep and ∼11 m wide channel became centred. This is capped by a second package of thin (<45 cm thick) sheets of pahoehoe extending no more than 50 m from the channel. The upper-most package comprises localised ′a′a overflows. The channel itself contains two blockages located 130 m apart and composed of levee chunks veneered with overflow lava. The channel was emplaced over 50 h, spanning 30 May–2 June, 1974, with the flow front arriving at our section (4.4 km from the vent) 8 h after the eruption began. The basal ′a′a flow thickness yields effusion rates of 35 m³ s⁻¹ for the opening phase, with the initial flow advancing across the mapped section at ∼10 m/min. Short-lived overflows of fluid pahoehoe then built the levee cap, increasing the apparent channel depth to 4.8 m. There were at least six pulses at 90–420 m³ s⁻¹, causing overflow of limited extent lasting no more than 5 min. Brim-full flow conditions were thus extremely short-lived. During a dominant period of below-bank flow, flow depth was ∼2 m with an effusion rate of ∼35 m³ s⁻¹, consistent with the mean output rate (obtained from the total flow bulk volume) of 23–54 m³ s⁻¹. During pulses, levee chunks were plucked and floated down channel to form blockages. In a final low effusion rate phase, lava ponded behind the lower blockage to form a syn-channel pond that fed ′a′a overflow. After the end of the eruption the roofed-over pond continued to drain through the lower blockage, causing the roof to founder. Drainage emplaced inflated flows on the channel floor below the lower blockage for a further ∼10 h. The complex processes involved in levee–channel construction of this short-lived case show that care must be taken when using channel dimensions to infer flow dynamics. In our case, the full channel depth is not exposed. Instead the channel floor morphology reflects late stage pond filling and drainage rather than true channel-contained flow. Components of the compound levee relate to different flow regimes operating at different times during the eruption and associated with different effusion rates, flow dynamics and time scales. For example, although high effusion rate, brim-full flow was maintained for a small fraction of the channel lifetime, it emplaced a pile of pahoehoe overflow units that account for 60% of the total levee height. We show how time-varying volume flux is an important parameter in controlling channel construction dynamics. Because the complex history of lava delivery to a channel system is recorded by the final channel morphology, time-varying flow dynamics can be determined from the channel morphology. Developing methods for quantifying detailed flux histories for effusive events from the evidence in outcrop is therefore highly valuable. We here achieve this by using high-resolution spatial data for a channel system at Kilauea. This study not only indicates those physical and dynamic characteristics that are typical for basaltic lava flows on Hawaiian volcanoes, but also a methodology that can be widely applied to effusive basaltic eruptions.     

Convection in a variable-viscosity fluid: Newtonian versus power-law rheology

A number of finite-element calculations of convection in a variable-viscosity fluid have been carried out to investigate the effects of non-Newtonian flow when rheology is also subject to a strong temperature and pressure influence. A variety of cases has been studied in the range of effective Rayleigh numbers between 10⁴ and 10⁶, including different modes of heating and a range of values for activation energy and activation volume. Power-law creep with a stress exponent of 3 turns out to lead to considerably different flow pattern and heat transfer properties than Newtonian rheology. In general, the effect is to reduce viscosity contrasts imposed by p,T dependence, which can lead in some circumstances to the mobilisation of otherwise stagnant regions within the cell. The properties of non-Newtonian flow can be closely imitated by a Newtonian fluid with a reduced value of the activation enthalpy bH^* with b?0.3–0.5. It appears possible that non-Newtonian rheology plays a key role in determining the convective style in a planetary mantle.     

Quantification of the hydraulic diffusivity of a bentonite-sand mixture using the water head decrease measured upon sudden flow interruption

This paper presents a new modelling approach to quantify the hydraulic diffusivity of low-permeability unconsolidated porous media under confined saturated-flow conditions in the laboratory. The derived analytical solution for the transient variation of the hydraulic head after flow interruption was applied to experimental data obtained from continuous measurements of the water pressure at two locations in the soil column. Three soil samples made of a mixture of natural bentonite (at different mass fractions) and medium sand were studied during a series of stepwise constant flow rates of water. The numerical results well fit the experimentally measured decrease of the dimensionless hydraulic head. The study shows that the increase of the mass fraction of bentonite in the soil sample from 10 to 30% is accompanied by a strong decrease of the hydraulic diffusivity from 2.4 × 10⁻² to 1.1 × 10⁻³ m² s⁻¹, which is clearly due to the decrease of the hydraulic conductivity of the soil sample. The specific storages obtained for each of the three samples are in the same order of magnitude and seem to decrease with the increase of mass fraction of bentonite. However, they clearly reflect the predominant portion of the compressibility of the porous medium compared with that of water. Compared with reported literature values for clayey soils, the specific storage values in this study are slightly higher, varying within the range of 2 × 10⁻³ to 8.1 × 10⁻³ m⁻¹._. The experimental results also give insight into the limitations of the modelling approach. In the case of low-permeability soils (K < 2 × 10⁻⁶ ms⁻¹) and steady-flow conditions with low Reynolds numbers, for example, Re < 0.003, it is recommended to choose a time duration for flow interruption between subsequent flow rate steps of longer than 5 s. For high-permeability porous media, to increase the precision of the quantified hydraulic diffusivity, it might be useful to select a measuring frequency significantly higher than 1 Hz.     

Estimation of water velocities in boreholes using a new passive flowmeter

Abstract

Preferential flow pathways in a fractured aquifer may yield abrupt reductions of the water velocity in a well. We propose a new device for measuring low (5–13 cm d^-1) velocities in wells originating from fractures at different depths. The presented flowmeter has been applied in a well in the Bari (southern Italy) fractured aquifer. In the same well, the horizontal flowmeter velocity (9.6 cm d^-1) at 0.5 m depth was compared with velocity (8 cm d^-1) derived from a field tracer test, providing a value 16.5% higher. Moreover, the flowmeter measurements at 1.5 m depth gave a horizontal velocity of 7.2 cm d^-1, which is 11% less than water flow velocity estimated from the field test. The new flowmeter implements the tracer point-dilution method in a plastic (PVC) pipe by causing the water flow to pass through an artificial filter. Laboratory calibration tests have confirmed the good performance of the proposed flowmeter technique, even for water flow up to 300 cm d^-1. The flowmeter was sensitive to 0.1 cm d^-1, with a detection limit of 1.5 cm d^-1, i.e. half the measurable flow velocity of existing flowmeters in wells.

Editor D. Koutsoyiannis; Associate editor S. Grimaldi     

Groundwater potential of the Egbe-Mopa basement area,central Nigeria

Abstract

An investigation on the groundwater potentials of the Egbe-Mopa area in central Nigeria, underlain by the Basement Complex, is presented. The investigation involved mapping of the subsurface by use of vertical electrical soundings; measurement of depth to groundwater; and evaluation of hydraulic conductivity, transmissivity and yield by means of pumping test interpretation. The results indicate subsurface units that range from three to five resistivity layers; depth to groundwater of 0–10 m; overburden thickness of 3–16 m; hydraulic conductivity of 6.2?×?10^?6 to 3.4?×?10^?4 m/s; transmissivity of 4.3?×?10^?7 to 2?×?10^?3 m²/s; and groundwater yield of 0.2–2.5 L/s. The hydraulic head assessments revealed a general northward groundwater flow direction. The study identified three aquifer potential types, of high, medium and low productivity, respectively. Based on the longitudinal conductance of the overburden units, four distinct Aquifer Protective Capacity zones were delineated, namely, poor, weak, moderate and good.

Citation Okogbue, C.O. and Omonona, O.V., 2013. Groundwater potential of Egbe-Mopa basement area, central Nigeria. Hydrological Sciences Journal, 58 (4), 826–840.     

A numerical study of detached shear layers in a rotating sliced cylinder

Abstract

The low Rossby number flow in a rotating cylinder with an inclined bottom, of small slope, is examined when part of the lid of the container is rotating at a slightly different rate. The resulting flow is calculated numerically by solving the governing equations for the two-dimensional geostrophic motion which approximates the flow in most of the fluid including the inertially-modified E ^¼ -layers. The presence of ageostrophic regions, on the container walls and beneath the velocity discontinuity on the lid, is accounted for in the governing equations and their boundary conditions. This study supplements previous work on this configuration, in which the zero Rossby number flow was calculated and experimental results were presented, by enabling a direct comparison to be made between the results of the low Rossby number theory and the experiments. The numerical results for a range of Rossby and Ekman numbers compare well with those from the experiments despite a severe limitation on the size of the Rossby number arising from the analysis in the ageostrophic part of the detached shear layer.     

River–floodplain hydrology of an embanked lowland Chalk river and initial response to embankment removal

Abstract

Rivers have been channelized, deepened and constrained by embankments for centuries to increase agricultural productivity and improve flood defences. This has decreased the hydrological connectivity between rivers and their floodplains. We quantified the hydrological regime of a wet grassland meadow prior to and after the removal of river embankments. River and groundwater chemistry were also monitored to examine hydrological controls on floodplain nutrient status. Prior to restoration, the highest river flows (～2 m³ s^?1) were retained by the embankments. Under these flow conditions the usual hydraulic gradient from the floodplain to the river was reversed so that subsurface flows were directed towards the floodplain. Groundwater was depleted in dissolved oxygen (mean: 0.6 mg O₂ L^?1) and nitrate (mean: 0.5 mg NO₃ ^?-N L^?1) relative to river water (mean: 10.8 mg O₂ L^?1 and 6.2 mg NO₃ ^?-N L^?1, respectively). Removal of the embankments has reduced the channel capacity by an average of 60%. This has facilitated over-bank flow which is likely to favour conditions for improved flood storage and removal of river nutrients by floodplain sediments.

Editor Z.W. Kundzewicz; Associate editor K. Heal

Citation Clilverd, H.M., Thompson, J.R., Heppell, C.M., Sayer, C.D., and Axmacher, J.C., 2013. River–floodplain hydrology of an embanked lowland Chalk river and initial response to embankment removal. Hydrological Sciences Journal, 58 (3), 627–650.     

Macroturbulent coherent structures in an ice‐covered river flow using a pulse‐coherent acoustic Doppler profiler

The aim of this work is to compare macroturbulent coherent structures (MCS) geometry and organization between ice covered and open channel flow conditions. Velocity profiles were obtained using a Pulse‐Coherent Acoustic Doppler Profiler in both open channel and ice‐covered conditions. The friction imposed by the ice cover results in parabolic shaped velocity profiles. Reynolds stresses in the streamwise (u) and vertical (v) components of the flow show positive values near the channel bed and negative values near the ice cover, with two distinctive boundary layers with specific turbulent signatures. Vertically aligned stripes of coherent flow motions were revealed from statistics applied to space‐time matrices of flow velocities. In open channel conditions, the macroturbulent structures extended over the entire depth of the flow whereas they were discontinued and nested close to the boundary walls in ice‐covered conditions. The size of MCS is consequently reduced in scale under an ice cover. The average streamwise length scale is reduced from 2.5 to 0.4Y (u) and from 1.5 to 0.4Y (v) where Y is the flow depth. In open channel conditions, the vertical extent of MCS covers the entire flow depth, whereas the vertical extent was in the range 0.58Y–1Y (u) and 0.81Y–1Y (v) in ice‐covered conditions. Under an ice cover, each boundary wall generates its own set of MCS that compete with each other in the outer region of the flow, enhancing mixing and promoting the dissipation of coherent structures. Copyright © 2012 John Wiley & Sons, Ltd.     

Effect of bed topography on soil aggregates transport by rill flow

After its formation, a rill may remain in the field for months, often receiving lower flow rates than the formative discharge. The objective of this work was to evaluate the rill flow transport capacity of soil aggregates at discharges unable to erode the rill, and to analyse the influence of the rill macro‐roughness on this transport process. A non‐erodible rill was built in which roughness was reproduced in detail. In order to assess only the rill macro‐roughness, a flat channel with a similar micro‐roughness to that in the rill replica was built. Rill and channel experiments were carried out at a slope of 8 and at six discharges (8·3 × 10^?5 to 5·2 × 10^?4 m³ s^?1) in the rill, and eight discharges (1·6 × 10^?5 to 5·2 × 10^?4 m³ s^?1) in the channel. Non‐erodible aggregates of three sizes (1–2, 3–5 and 5–10 mm) were released at the inlet of the rill/channel. The number of aggregates received at the outlet was registered. The number and position of the remaining aggregates along the rill/channel were also determined. The rill flow was a major sediment transport mechanism only during the formation of the rill, as during that period the power of the flow was great enough to overcome the influence of the macro‐roughness of the rill bed. At lower discharges the transport capacity in the previously formed rill was significantly less than that in the flat channel under similar slope and discharge. This was determined to be due to local slowing of flow velocities at the exit of rill pools. Copyright © 2006 John Wiley & Sons, Ltd.     

The destabilising nature of differential rotation

Abstract

In a rapidly rotating, electrically conducting fluid we investigate the thermal stability of the fluid in the presence of an imposed toroidal magnetic field and an imposed toroidal differential rotation. We choose a magnetic field profile that is stable. The familiar role of differential rotation is a stabilising one. We wish to examine the less well known destabilising effect that it can have. In a plane layer model (for which we are restricted to Roberts number q = 0) with differential rotation, U = sΩ(z)1 _?, no choice of Ω(z) led to a destabilising effect. However, in a cylindrical geometry (for which our model permits all values of q) we found that differential rotations U = sΩ(s)1 _? which include a substantial proportion of negative gradient (dΩ/ds ≤ 0) give a destabilising effect which is largest when the magnetic Reynolds number R _m = O(10); the critical Rayleigh number, Ra _c, is about 7% smaller at minimum than at R_m = 0 for q = 10⁶. We also find that as q is reduced, the destabilising effect is diminished and at q = 10^?6, which may be more appropriate to the Earth's core, the effect causes a dip in the critical Rayleigh number of only about 0.001%. This suggests that we see no dip in the plane layer results because of the q = 0 condition. In the above results, the Elsasser number A = 1 but the effect of differential rotation is also dependent on A. Earlier work has shown a smooth transition from thermal to differential rotation driven instability at high A [A = O(100)]. We find, at intermediate A [A = O(10)], a dip in the Ra_c vs. R_m curve similar to the A = 1 case. However, it has Ra_c ≤ 0 at its minimum and unlike the results for high A, larger values of R_m result in a restabilisation.     

Spatial variability of hydraulic conductivity and bulk density along a blanket peatland hillslope

This article presents the results of a field investigation of saturated hydraulic conductivity K_sat and bulk density (ρ_bd) in an Atlantic blanket bog in the southwest of Ireland. Starting at a peatland stream and moving along an uphill transect toward the peatland interior, ρ_bd and K_sat were examined at regular intervals. Saturated horizontal hydraulic conductivity (Kh_sat) and vertical (Kv_sat) was estimated at two depths: 10–20 and 30–40 cm below the peat surface, whereas ρ_bd was estimated for the full profile. We consider two separate zones, one a riparian zone extending 10 m from the stream and a second zone in the bog interior. We found that the K_sat was higher (~10^–5 m s^–1) in the bog interior than that in the riparian zone (~10^–6 m s^–1), whereas the converse applied to bulk density, with lowest density (~0.055 g cm^–3) at the interior and highest (~0.11 g cm^–3) at the riparian zone. In general, we found Kh_sat to be approximately twice the Kv_sat. These results support the idea that the lower K_sat at the margins control the hydrology of blanket peatlands. It is therefore important that the spatial variation of these two key properties be accommodated in hydrological models if the correct rainfall runoff characteristics are to be correctly modelled. Stream flow analysis over 3 years at the peatland catchment outlet showed that the stream runoff was composed of 8% base flow and 92% flood flow, suggesting that this blanket peatland is a source rather than a sink for floodwaters. Copyright © 2011 John Wiley & Sons, Ltd.