Evolution of non-linear α 2-dynamos and taylor's constraint

A key non-linear mechanism in a strong-field geodynamo is that a finite amplitude magnetic field drives a flow through the Lorentz force in the momentum equation and this flow feeds back on the field-generation process in the magnetic induction equation, equilibrating the field. We make use of a simpler non-linear?α?²-dynamo to investigate this mechanism in a rapidly rotating fluid spherical shell. Neglecting inertia, we use a pseudo-spectral time-stepping procedure to solve the induction equation and the momentum equation with no-slip velocity boundary conditions for a finitely conducting inner core and an insulating mantle. We present calculations for Ekman numbers (E) in the range 2.5× 10^?3 to 5.0× 10^?5, for?α?=α ₀cos?θ?sin?π?(r?r_i ) (which vanishes on both inner and outer boundaries). Solutions are steady except at lower E and higher values of?α?₀. Then they are periodic with a reversing field and a characteristic rapid increase then equally rapid decrease in magnetic energy. We have investigated the mechanism for this and shown the influence of Taylor's constraint. We comment on the application of our findings to numerical hydrodynamic dynamos.     

Space-Time scales of internal waves

Abstract

We have contrived a model E(αω) α μ^?1ω^?p+1(ω ²?ω _i ²)^?+ for the distribution of internal wave energy in horizontal wavenumber, frequency-space, with wavenumber α extending to some upper limit μ(ω) α ω ^r-1 (ω ²?ω _i ²)^½, and frequency ω extending from the inertial frequency ω _i to the local Väisälä frequency n(y). The spectrum is portrayed as an equivalent continuum to which the modal structure (if it exists) is not vital. We assume horizontal isotropy, E(α, ω) = 2παE(α₁, α₂, ω), with α₁, α₂ designating components of α. Certain moments of E(α₁, α₂, ω) can be derived from observations. (i) Moored (or freely floating) devices measuring horizontal current u(t), vertical displacement η(t),…, yield the frequency spectra F _(u,η,…)(ω) = ∫∫ (U ², Z ²,…)E(α₁, ∞₂, ω) dα₁ dα₂, where U, Z,… are the appropriate wave functions. (ii) Similarly towed measurements give the wavenumber spectrum F _(…)(α1) = ∫∫… dα₂ dω. (iii) Moored measurements horizontally separated by X yield the coherence spectrum R(X, ω) which is related to the horizontal cosine transform ∫∫ E(α1, α2 ω) cos α₁ Xdα₁ dα₁. (iv) Moored measurements vertically separated by Y yield R(Y, ω) and (v) towed measurements vertically separated yield R(Y, α₁), and these are related to similar vertical Fourier transforms. Away from inertial frequencies, our model E(α, ω) α ω ^?p-r for α ≦ μ ω ω ^r, yields F(ω) ∞ ω ^?p, F(α₁) ∞ α₁ ^?q, with q = (p + r ? 1)/r. The observed moored and towed spectra suggest p and q between 5/3 and 2, yielding r between 2/3 and 3/2, inconsistent with a value of r = 2 derived from Webster's measurements of moored vertical coherence. We ascribe Webster's result to the oceanic fine-structure. Our choice (p, q, r) = (2, 2, 1) is then not inconsistent with existing evidence. The spectrum is E(∞, ω) ∞ ω ^?1(ω ²?ω _i ^{2 ?1}, and the α-bandwith μ ∞ (ω ²?ω _i ²)⁺ is equivalent to about 20 modes. Finally, we consider the frequency-of-encounter spectra F([sgrave]) at any towing speed S, approaching F(ω) as S ≦ S _o, and F(α₁) for α₁ = [sgrave]/S as S ≧ S _o, where S _o = 0(1 km/h) is the relevant Doppler velocity scale.     

A laboratory study of baroclinic instability

Abstract

Experimental and theoretical results are presented for a simple system which exhibits baroclinic instability. We consider the motion of two immiscible fluids with densities ρ ₁ and ρ ₂ contained in a cylinder rotating with angular frequency ω. The motion is driven by a contact lid rotating with frequency ω + ω. In this paper ω, ω, 2(ρ ₂ – ρ ₁)/(ρ ₂ + ρ ₁), and the geometry are such that the interface does not intersect the “ground” (e.g. an almost horizontal boundary). The motions are described by two-layer quasi-geostrophic equations which are identical, except perhaps for the presence of interfacial friction and tension, with those used in meteorology and oceanography. For small enough internal Froude number F = 4ω² L ²/(g(Δρ/ρ)H) or small enough Rossby number ? = ω/2ω the flow is steady and axisymmetric, the velocity field in each layer being determined primarily by frictional effects in top, bottom, and interfacial Ekman layers. For certain (F, ?) the flow becomes non-axisymmetric. The transition points for the case where the basic potential vorticity gradient is due to interface slope alone have been carefully measured and are in very good agreement with a linear instability theory which neglects sidewall effects. Some preliminary observations of supercritical motion, which include repeatable amplitude and wavenumber vacillation, are reported.     

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.     

The relationship between vertical eddy diffusion and buoyancy gradient in the deep sea

Vertical eddy diffusivities (K_v's) have been estimated at fourteen widely separated locations from fourteen²²²Rn profiles and two²²⁸Ra profiles measured near the ocean floor as part of the Atlantic and Pacific GEOSECS programs. They show an inverse proportionality to the local buoyancy gradient [(g/?)(??_pot/?z)] calculated from hydrographic measurements. The negative of the constant of proportionality is the buoyancy flux [?K_v(g/?)(??_pot/?z)] which has a mean of ?4 × 10^?6 cm²/sec³. Our results suggest that the buoyancy flux varies very little near the ocean floor. K_v's for the interior of the deep Pacific calculated from the relationship K_v = (4 × 10^?6cm²/sec³)/[(g/?)(??_pot/?z)] agree well with published estimates. K_v's calculated for the pycnocline are one to two orders of magnitude smaller than upper limits estimated from tritium and⁷Be distributions.Heat fluxes calculated with the model K_v's obtained from the²²²Rn profiles average 31 μcal cm^?2 sec^?1 in the Atlantic Ocean and 8 μcal cm^?2 sec^?1 in the Pacific Ocean.     

First approximations of soil moisture retention curves using the filter‐paper method

Relationships between gravimetric soil moisture content (w) and matric potential (ϕ), and between volumetric soil moisture content (θ_v) and pressure head (h) were approximated for the unsaturated zone on Long Island, New York. Soil samples were collected from two sites using a hand auger. The soil moisture content was determined using the filter‐paper (w_f) and gravimetric (w) methods, respectively. The w_f was then used in an empirical equation to estimate ϕ_m. Each set of ϕ_m and w was combined with a straight‐line empirical model to obtain a w(ϕ_m) relationship. Soil ϕ_m was converted to h, and w to the volumetric moisture content θ_v, in order to produce a θ_v(h) curve. Graphical and statistical comparison showed that the resulting θ_v(h) curves fell within one order of magnitude of similar curves generated by a more sophisticated non‐linear model developed previously. The simplicity and low cost of the filter‐paper approach described in this study recommends it for preliminary studies of hydraulic properties in the unsaturated zone. Copyright © 2000 John Wiley & Sons, Ltd.     

Non-axisymmetric α2Ω-dynamo waves in thin stellar shells

Linear α²Ω-dynamo waves are investigated in a thin turbulent, differentially rotating convective stellar shell. A simplified one-dimensional model is considered and an asymptotic solution constructed based on the small aspect ratio of the shell. In a previous paper Griffiths et al. (Griffiths, G.L., Bassom, A.P., Soward, A.M. and Kuzanyan, K.M., Nonlinear α²Ω-dynamo waves in stellar shells, Geophys. Astrophys. Fluid Dynam., 2001, 94, 85–133) considered the modulation of dynamo waves, linked to a latitudinal-dependent local α-effect and radial gradient of the zonal shear flow. These effects are measured at latitude θ by the magnetic Reynolds numbers R _α f(θ) and R _Ω g(θ). The modulated Parker wave, which propagates towards the equator, is localised at some mid-latitude θ_p under a Gaussian envelope. In this article, we include the influence of a latitudinal-dependent zonal flow possessing angular velocity Ω_*(θ) and consider the possibility of non-axisymmetric dynamo waves with azimuthal wave number m. We find that the critical dynamo number D _c?=?R _α R _Ω is minimised by axisymmetric modes in the αΩ-limit (Rα→0). On the other hand, when Rα?≠?0 there may exist a band of wave numbers 0?m?m ^? for which the non-axisymmetric modes have a smaller D _c than in the axisymmetric case. Here m ^? is regarded as a continuous function of R _α with the property m?→0 as R _α→0 and the band is only non-empty when m??>1, which happens for sufficiently large R _α. The preference for non-axisymmetric modes is possible because the wind-up of the non-axisymmetric structures can be compensated by phase mixing inherent to the α²Ω-dynamo. For parameter values resembling solar conditions, the Parker wave of maximum dynamo activity at latitude θ_p not only propagates equatorwards but also westwards relative to the local angular velocity Ω_*(θ _p ). Since the critical dynamo number D _c?=?R _α R _Ω is O (1) for small R _α, the condition m ^??>?1 for non-axisymmetric mode preference imposes an upper limit on the size of |dΩ_*/dθ|.     

Effect of non‐uniformity of flow on velocity and turbulence intensities over a cobble‐bed

Non‐uniform flows encompassing both accelerating and decelerating flows over a cobble‐bed flume have been experimentally investigated in a flume at a scale of intermediate relative submergence. Measurements of mean longitudinal flow velocity u, and determinations of turbulence intensities u′, v′, w′, and Reynolds shear stress ?u_fw_f have been made. The longitudinal velocity distribution was divided into the inner zone close to the bed and the outer zone far from the bed. In the inner zone of the boundary layer (near the bed) the velocity profile closely followed the ‘Log Law’; however, in the outer zone the velocity distribution deviated from the Log Law consistently for both accelerating and decelerating flows and the changes in bed slopes ranging from ?2% to + 2% had no considerable effect on the outer zone. For a constant bed slope (S = ±0·015), the larger the flow rate, the smaller the turbulence intensities. However, no detectable pattern has been observed for u′, v′ and w′ distributions near the bed. Likewise, for a constant flow rate (Q = 0·040 m³/s), with variation in bed slope the longitudinal turbulent intensity profile in the longitudinal direction remained concave for both accelerating and decelerating flows; whereas vertical turbulent intensity (w′) profile presented no specific form. The results reveal that the positions of maximum values of turbulence intensities and the Reynolds shear stress depend not only on the flow structure (accelerating or decelerating) but also on the intermediate relative submergence scale. Copyright © 2009 John Wiley & Sons, Ltd.     

Experimental study of the Ekman layer instability in steady or oscillating flows

To investigate the instabilities of steady and oscillating Ekman layers, an 8 m×2 m horizontal plate was moved at controlled speed in homogeneous water at rest in solid body rotation in the “Coriolis” 13 m diameter rotating tank. For a steady Ekman layer two distinct wave types were found, in agreement with previous experimental or numerical studies. Type I was stationary, was oriented positively with respect to the flow direction and had a wavelength of about 10 times the Ekman layer thickness. Type II was oriented negatively with respect to the flow direction and had a wavelength which was more than 20 times the Ekman layer thickness and a phase-speed between 0.3 and 0.5 the forcing interior velocity. The growth rates of both type I and type II waves for various Reynolds numbers Re (computed with the Ekman layer thickness) were estimated and their Re-variations qualitatively agree with previous numerical results. For an oscillating Ekman layer, experimental results depended strongly on Ro_t, the temporal Rossby number: only when Ro_t<1 was it possible to observe either type I or type II instabilities. Moreover, for all Ro_t and average to high Re, there was a noticeable upward turbulent transport occurring during each cycle between the flow maximum and the flow reversal. Such an upward turbulent transport is consistent with observations in the English Channel where maximum upward benthic movements and maximum turbidity were recorded at the flow reversal, hence Ekman layer instabilities and transition to turbulence are likely to occur in shallow tidal seas where they may be relevant for sediment resuspension and transport as well as for some biological processes.     

Determining thermal diffusivity using near‐surface periodic temperature variations and its implications for tracing groundwater movement at the eastern margin of the Tibetan Plateau

The thermal diffusivity is the key parameter that controls near‐surface temperature where periodic temperature variation is progressively attenuated and delayed with depth. This article presents the results of apparent thermal diffusivity using temperatures recorded by a bedrock temperature measurement network in the fault zones of western Sichuan. High sensitivity temperature sensors (10^?4 K) were installed at a maximum depth reaching 30 m. The apparent thermal diffusivities were deduced from both amplitude damping and phase shifting of annual temperature variations between two different depths. Under pure conduction, the thermal diffusivity determined through the phase method (α_Φ) should be equivalent to that determined through the amplitude method (α_A), whereas effects of the upward (downward) water flow are evidently reflected in the amplitude decay to make α_Φ larger (lesser) than α_A. The discrepancy between α_Φ and α_A can thus be a tracer of water movement or convective heat transfer. The calculated α_Φ of the measurement stations varies from 1.22 × 10^?6 to 3.00 × 10^?6 m²/s, and the estimated α_A ranges from 0.93 × 10^?6 to 2.41 × 10^?6 m²/s. Two regimes of heat transfer underground were suggested from the results. Conductive heat transport prevails over the nonconductive processes at five stations, which is characterized by α_Φ coincident with α_A for the same depth pair. On the contrary, the values of α_Φ differ from α_A at six stations in the intersection area of the Y‐shaped fault system, implying that convective heat transfer also plays a comparably important role. This finding is consistent with the hot springs distribution of the area. The results also indicate that water moves upward with an average Darcy velocity of approximately ?1 × 10^?7 m/s in this region. Our research provides new evidence for the hydrothermal activity in the fault zones at the eastern margin of the Tibetan Plateau.     

The equatorial dynamics of a shallow,homogeneous ocean

Abstract

It is shown that the linear equatorial dynamics of a shallow ocean is characterized by two boundary layers of width γ^? L and γL (γ is the Ekman number of the flow, assumed small, and L is a horizontal dimension of the basin). In the γ^? layer stress in the bottom Ekman layer is comparable to that in the surface Ekman layer. In the γ layer vertical friction is important throughout the depth of the ocean. Should the Rossby number ? be so large as to invalidate a linear theory (? > γ^5/3), then inertial effects become important at a distance ?^2/5 L from the equator. The role played in the circulation of the basin by the non-linear equatorial current first studied by Charney (1960) is shown to be similar to that of the γ layer of the linear theory. Though lateral friction is unimportant in a linear model of the flow, shear layers at the equator are found to be a necessary feature of non-linear flow.     

Controls on effective settling velocities of mud flocs in the Changjiang Estuary,China

Theoretical and experimental studies were undertaken to gain insight into physical parameters controlling the flocculation and settling properties of mud flocs in the Changjiang Estuary, China. The Rouse equation is applied to vertical profiles of suspended sediment concentration to determine the bulk mean settling velocity (w_s) of sediment suspended in the Changjiang Estuary. Both in situ point‐sampled and acoustically measured profiles of suspended mud concentrations were fit selectively. The calculated settling velocities w_s mainly ranged from 0·4 to 4·1 mm s^?1 for the point‐sampled data set, and from 1·0 to 3·0 mm s^?1 for the acoustically measured data set. Furthermore, the settling velocities of mud flocs increased with mean concentration (C?) of mud flocs in suspension and were proportional to increasing bottom shear stress (τ_b) of tidal flow. The best equation for the field settling velocity of mud flocs in the Changjiang Estuary can be expressed by the power law: w_s = mC?ⁿ (m, 1·14–2·37; n, 0·84–1·03). It is suggested that C? and τ_b were the dominant physical parameters controlling the flocculation and w_s of mud flocs in suspension. Copyright © 2004 John Wiley & Sons, Ltd.     

Excess222Rn and the benthic boundary layer in the western and southern Indian Ocean

We present 9 bottom²²²Rn profiles measured from the western and southern Indian Ocean during the 1977–1978 GEOSECS expedition. These profiles can be grouped into three cypes: one-layer, two-layer, and irregular types. The one-layer profiles with quasi-exponential distributions allow one to estimate the apparent vertical eddy diffusivity,K_v, with a simple model. The two-layer profiles show that there is a benthic boundary layer of the order of 50–100 m in which the excess²²²Rn distribution shows a vertical gradient much smaller than that of the layer immediately above. Within the boundary layer, the STD potential temperature (θ) and density(σ₄) profiles are practically constant, and theK_v values are of the order of 1000 cm²/s. The STD profiles for the water column above the boundary layer show gradients of increasing stability, and theK_v values are of the order of 100 cm²/s. Modeling of the Rn data in the water column above the boundary layer indicates that there is a transition layer which effectively reduces the penetration of excess Rn from the benthic boundary layer into the upper layer.Sarmiento et al. [10] have shown that the buoyancy gradient or stability is inversely correlated with the apparent vertical eddy diffusivity, and the resulting buoyancy flux is fairly uniform, ranging from 1 to 14 × 10^?6 cm²/s³ in the Atlantic and Pacific Oceans. However, Sarmiento et al. [11] show that a much higher buoyancy flux is associated with an intensified flow of the bottom water through a passage. In the Indian Ocean basins, we have found that the buoyancy flux has a comparable range (3–14 × 10^?6 cm²/s³), except for a couple of stations where both stability and apparent vertical diffusivity are higher, resulting in a much higher buoyancy flux, probably indicative of rapid bottom water flow.     

α2-Dynamos and taylor's constraint

Abstract

An idealised α²ω-dynamo is considered in which the α-effect is prescribed. The additional ω-effect results from a geostrophic motion whose magnitude is determined indirectly by the Lorentz forces and Ekman suction at the boundary. As the strength of the α-effect is increased, a critical value α? _c is reached at which dynamo activity sets in; α? _c is determined by the solution of the kinematic α²-dynamo problem. In the neighbourhood of the critical value of α? the magnetic field is weak of order E ^1/4(μηρω)^½ due to the control of Ekman suction; E(?1) is the Ekman number. At certain values of α?, viscosity independent solutions are found satisfying Taylor's constraint. They are identified by the bifurcation of a nonlinear eigenvalue problem. Dimensional arguments indicate that following this second bifurcation the magnetic field is strong of order (μηρω)^½. The nature of the transition between the kinematic linear theory and the Taylor state is investigated for various distributions of the α-effect. The character of the transition is found to be strongly model dependent.     

Specific features of development of two-dimensional convection in a rotating compressible conducting fluid

An axisymmetric model of convection in a rotating cylinder in an external uniform magnetic field has been considered. In the considered model, the meridional circulation is created by a nonuniform rotation of the lower boundary relative to the other boundaries. In the considered model, the time of formation of the stationary regime in the magnetic field considerably increases if the vertical density (compressibility) inhomogeneity is taken into account for Ekman numbers of E = E _M = 3 × 10⁻³. This example shows that the compressibility of a medium should be taken into account in the convection and dynamics of the magnetic field when the magnetohydrodynamics of the Earth is analyzed.     

Ion microprobe studies of water in silicate melts. Concentration-dependent water diffusion in obsidian

The ion microprobe at Johnson Space Center has been calibrated for in situ water determinations on a 10-μm scale over the range 0.2 wt.% H₂O to 1.8, 6.8, and 3.7 wt.%, for basaltic, albitic, and rhyolitic glasses, respectively. The basalt glass calibration curve differs substantially from those of albite and rhyolite glasses, indicating a need to carefully match composition and/or melt structure between H₂O standards and unknowns.A value for the diffusivity of water as a function of concentration and time has been calculated from water diffusion profiles measured in rhyolite glasses prepared at 850°C and 700 barsP_t(H₂O) [1]. Transient diffusion into a semi-infinite medium is described by the equation:?(φ/2)?¸/?φ=?(D_w?¸/?φ)/?φ #x003B8;=1, φ=0, θ→ 0, θ→∞, wherex =distance from the cylinder edge,t =time,C₀ =initial concentration,C_s =concentration at the edge,C =concentration at x,θ = C ? C₀/C_s ? C₀,φ = x/t^1/2, andD_w =diffusivity of water. An iterative technique has been used to calculate solutions to the diffusion equation as a function ofD_w [2]. Comparison of these solutions with the ion probe data indicate that, for0.2wt.% ≤ C ≤ 3.7wt.%H₂O,D_w can be described by an exponential function of θ, of the formD_w = D₀exp(bθ), withD₀ (i.e.,D_w at 0.2%) = (0.8?2.2) × 10^?8 cm²/s and2 ≤ b ≤ 4.     

Measurements and estimation of flow velocity in mobile bed rills

Rills are primary sediment sources and hillslope water/sediment runoff transport channels. Water flow velocities in rills are easily affected by bed condition over eroding and non-eroding slopes, which is an important hydrodynamic process in soil erosion research. This research is done to demonstrate the poorly understood “feedback mechanism” related to slope independence of flow velocity to slope gradient. A series of experiments were done on silt loam soil slopes to measure water flow velocity...     

Source Characteristics of the 22 January 2003 M w?=?7.5 Tecom��n, Mexico, Earthquake: New Insights

Aftershock locations, source parameters and slip distribution in the coupling zone between the overriding North American and subducted Rivera and Cocos plates were calculated for the 22 January 2003 Tecomán earthquake. Aftershock locations lie north of the El Gordo Graben with a northwest-southeast trend along the coast and superimposed on the rupture areas of the 1932 (M _w?=?8.2) and 1995 (M _w?=?8.0) earthquakes. The Tecomán earthquake ruptured the northwest sector of the Colima gap, however, half of the gap remains unbroken. The aftershock area has a rectangular shape of 42?±?2 by 56?±?2?km with a shallow dip of roughly 12° of the Wadati-Benioff zone. Fault geometry calculated with the Náb??lek (1984) inversion procedure is: (strike, dip, rake)?=?(277°, 27°, 78°). From the teleseimic body wave spectra and assuming a circular fault model, we estimated source duration of 20?±?2?s, a stress drop of 5.4?±?2.5?MPa and a seismic moment of 2.7?±?.7?×?10²⁰?Nm. The spatial slip distribution on the fault plane was estimated using new additional near field strong motion data (54?km from the epicenter). We confirm their main conclusions, however we found four zones of seismic moment release clearly separated. One of them, not well defined before, is located toward the coast down dip. This observation is the result of adding new data in the inversion. We calculated a maximum slip of 3.2?m, a source duration of 30?s and a seismic moment of 1.88?×?10²⁰?Nm.     

Using a transient infiltrometric technique for intensively sampling field‐saturated hydraulic conductivity of a clay soil in two runoff plots

The point measurement of soil properties allows to explain and simulate plot scale hydrological processes. An intensive sampling was carried out at the surface of an unsaturated clay soil to measure, on two adjacent plots of 4 × 11 m² and two different dates (May 2007 and February–March 2008), dry soil bulk density, ρ_b, and antecedent soil water content, θ_i, at 88 points. Field‐saturated soil hydraulic conductivity, K_fs, was also measured at 176 points by the transient Simplified Falling Head technique to determine the soil water permeability characteristics at the beginning of a possible rainfall event yielding measurable runoff. The ρ_b values did not differ significantly between the two dates, but wetter soil conditions (by 31%) and lower conductivities (1.95 times) were detected on the second date as compared with the first one. Significantly higher (by a factor of 1.8) K_fs values were obtained with the 0.30‐m‐diameter ring compared with the 0.15‐m‐diameter ring. A high K_fs (> 100 mm h^?1) was generally obtained for low θ_i values (< 0.3 m³m^?3), whereas a high θ_i yielded an increased percentage of low K_fs data (1–100 mm h^?1). The median of K_fs for each plot/sampling date combination was not lower than 600 mm h^?1, and rainfall intensities rarely exceeded 100 mm h^?1 at the site. The occurrence of runoff at the base of the plot needs a substantial reduction of the surface soil permeability characteristics during the event, probably promoted by a higher water content than the one of this investigation (saturation degree = 0.44–0.62) and some soil compaction due to rainfall impact. An intensive soil sampling reduces the risk of an erroneous interpretation of hydrological processes. In an unstable clay soil, changes in K_fs during the event seem to have a noticeable effect on runoff generation, and they should be considered for modeling hydrological processes. Copyright © 2012 John Wiley & Sons, Ltd.     

SUR LES PROJECTS FUTURS DE L'AISH / ON THE FUTURE PLANS OF IAHS

Abstract

The baseflow characteristics of some of the numerous small basins in southeastern Nigeria have been analysed to estimate the developable groundwater in the basins. It is shown that from 5.62 × 10⁴ to 1.59 × 10⁶ m³ of groundwater can be developed per square kilometre of basin per annum. The relationship between the baseflow characteristics and other attributes of the basins, such as geology and stream density, were studied statistically, leading to the development of empirical equations for predicting the hydrological features of the several ungauged streams in the region. It is shown, for example, that the basin geology (represented as the percentage of sands), the drainage density, the basin area, the baseflow depletion rate and the total groundwater stored in the basin, Q_tp, are related by the equation:

Q_tp = ?1.85 × 10⁹?7.96 × 10⁸ dd+4.18 × 10⁷ gf?2.01 × 10⁶ df+6.25 × 10⁵ wa

where dd is drainage density; gf geological factor; df depletion factor; and wa basin area.