Incorporating variable source area hydrology into a curve-number-based watershed model

Many water quality models use some form of the curve number (CN) equation developed by the Soil Conservation Service (SCS; U.S. Depart of Agriculture) to predict storm runoff from watersheds based on an infiltration-excess response to rainfall. However, in humid, well-vegetated areas with shallow soils, such as in the northeastern USA, the predominant runoff generating mechanism is saturation-excess on variable source areas (VSAs). We reconceptualized the SCS–CN equation for VSAs, and incorporated it into the General Watershed Loading Function (GWLF) model. The new version of GWLF, named the Variable Source Loading Function (VSLF) model, simulates the watershed runoff response to rainfall using the standard SCS–CN equation, but spatially distributes the runoff response according to a soil wetness index. We spatially validated VSLF runoff predictions and compared VSLF to GWLF for a subwatershed of the New York City Water Supply System. The spatial distribution of runoff from VSLF is more physically realistic than the estimates from GWLF. This has important consequences for water quality modeling, and for the use of models to evaluate and guide watershed management, because correctly predicting the coincidence of runoff generation and pollutant sources is critical to simulating non-point source (NPS) pollution transported by runoff. Copyright © 2007 John Wiley & Sons, Ltd.     

Comment on the effect of anisotropy on the onset of convection in a porous medium

The effect of anisotropy on the onset of convection in a saturated porous medium is discussed. In particular, the case of time-dependent density-driven convection is examined. The applicability of the value of an equivalent Rayleigh number as the criterion for the onset of convection is discussed.     

人工增雨火箭发射仰角及用弹量的估算

根据南阳气象观测站探空资料及AgI在云中成核率与温度的关系,确定了豫南BL-1型火箭人工增雨可播云体的判别标准、火箭发射仰角及火箭用弹量。     

Remote sensing and GIS for artificial recharge study, runoff estimation and planning in Ayyar basin, Tamil Nadu, India

This paper focuses on artificial groundwater recharge study in Ayyar basin, Tamil Nadu, India. The basin is covered by hard crystalline rock and overall has poor groundwater conditions. Hence, an artificial recharge study was carried out in this region through a project sponsored by Tamil Nadu State Council for Science and Technology. The Indian Remote Sensing satellite 1A Linear Imaging Self Scanning Sensor II (IRS 1A LISS II) satellite imagery, aerial photographs and geophysical resistivity data were used to prioritize suitable sites for artificial recharge and to estimate the volume of aquifer dimension available to recharge. The runoff water available for artificial recharge in the basin is estimated through Soil Conservation Service curve number method. The land use/land cover, hydrological soil group and storm rainfall data in different watershed areas were used to calculate the runoff in the watersheds. The weighted curve number for each watershed is obtained through spatial intersection of land use/land cover and hydrological soil group through GeoMedia 3.0 Professional GIS software. Artificial recharge planning was derived on the basis of availability of runoff, aquifer dimension, priority areas and water table conditions in different watersheds in the basin.     

Drag on non-spherical particles in power law non-Newtonian media

The free settling velocity of cylinders and disks falling in quiescent Newtonian and power law liquids has been measured over wide ranges of experimental conditions of the particle Reynolds number (10^− 5–∼300), power law flow behaviour index (0.31–1) and the length-to-diameter ratio, ∼0.4–∼14. The corresponding range of sphericity is 0.62 to 0.86. An existing drag expression which has been tested extensively for spherical particles falling in Newtonian and in power law fluids has been slightly modified here for non-spherical particles. In particular, the use of this drag expression necessitates a knowledge of an equal volume sphere diameter (to evaluate the Reynolds number and drag coefficient) and the ratio of the surface area to the projected area of a non-spherical particle. With these modifications, the approach outlined here reproduces the present and the literature data for a wide range of non-spherical particles including cones, prisms, needles, cylinders settling in both Newtonian and power law fluids with reasonable levels of accuracy.     

The Effect of Relative Sunspot Numbers,Solar Flare Numbers and Variable Component of 10.7 cm Solar Flux on The Seasonal Variation of 6300 Å Line Intensity at Calcutta

The paper presents the seasonal variation of 6300 Å line intensity at Calcutta with relative sunspot number, solar flare number and variable component of 10.7 cm solar flux. A study has been made and important results have been obtained which are as follows. (i) Intensity of 6300 Å line shows periodic variation with relative sunspot number, solar flare number and variable component of 10.7 cm solar flux during the period 1984–1986 which is the secondary peak of the descending phase of 21st solar cycle. (ii) 6300 Å line intensity at Cachoeira Paulista station, taken by Sahai et al. (1988), also shows periodic variation with solar parameters during the period 1978–1980 which is the peak phase of the solar cycle. (iii) A possible explanation of such a type of variation is also presented.     

Phosphorus in granitic rocks of Colombia

One hundred and three samples of granitic rock taken systematically from the Andean part of Colombia reveal that the observed variation in phosphorus content reflects the variation observed for the wholerock chemical composition. This relationship was established using a form of trend analysis termed COMTRENA and information-theory statistics.     

Classification of Vertical Wind Speed Profiles Observed Above a Sloping Forest at Nighttime Using the Bulk Richardson Number

Wind speed profiles above a forest canopy relate to scalar exchange between the forest canopy and the atmosphere. Many studies have reported that vertical wind speed profiles above a relatively flat forest can be classified by a stability index developed assuming wind flow above a flat plane. However, can such a stability index be used to classify vertical wind speed profiles observed above a sloping forest at nighttime, where drainage flow often occurs? This paper examines the use of the bulk Richardson number to classify wind speed profiles observed above a sloping forest at nighttime. Wind speed profiles above a sloping forest were classified by the bulk Richardson number Ri _B . Use of Ri _B distinguished between drainage flow, shear flow, and transitional flow from drainage flow to shear flow. These results suggest that Ri _B is useful to interpret nighttime CO₂ and energy fluxes above a sloping forest. Through clear observational evidence, we also show that the reference height should be high enough to avoid drainage-flow effects when calculating Ri _B .     

An Approximate Analytical Solution Of Flux-Profile Relationships For The Atmospheric Surface Layer With Different Momentum And Heat Roughness Lengths

An approximate method for calculating the relationship between z/L(z = reference height, L = Obukhov length) and the bulk Richardsonnumber is presented. If this relationship is known, the momentum andheat fluxes can be computed easily without any iteration. The avoidance of iteration can speed up computationsin large-scale models considerably (up to 10 times) and cases which do not converge or converge very slowly cannot occur. The proposed formulae take into account the difference between momentum (z_0M) and heat roughnesslengths (z_0H). Because the roughness lengths are not neglected at any step of the derivation, the resulting analytical formulae can be used not only between the surface and the reference height but also between two finite levels z₁ andz₂ (by replacing z_0M and z_0H by z₁ and z by z₂). Theequations remain correct even in the limit z₁ z₂.The formulae are based upon the (partially modified) Businger–Dyer flux–profile relationships and,consequently, they are restricted to predominantly homogeneous terrain.These new approximations are an improvement over the existing solutions because they are simpler than most of the formulae in the literature and are able to match the numerical exact solution for different parameter sets (Businger, Dyer, Högström) with an maximum error of about 2% for a wide range of z/L, z/z_0M and z_0M/z_0H.Furthermore, in stable conditions, schemes with and without a finitecritical bulk Richardson number can be approximated. The possibleambiguity of the exact solution =f(RI_B) in (moderately) stable conditions is discussed briefly. The performance of the new formulae is compared to the exact numerical solution and to different formulae proposed in the literature.