Identification of parameters in unconfined aquifers

A method is presented for identification of parameters in unconfined aquifers from pumping tests, based on the optimisation of the objective function using the least squares approach. Four parameters are to be evaluated, namely: The hydraulic conductivity in the radial and the vertical directions, the storage coefficient and the specific yield. The sensitivity analysis technique is used for solving the optimisation problem. Besides eliminating the subjectivity involved in the graphical procedure, the method takes into account the field data at all time intervals without classifying them into small and large time intervals and does not use the approximation that the ratio of the storage coefficient to the specific yield tends to zero. Two illustrative examples are presented and it is found that the parameter estimates from the computational and graphical procedures differ fairly significantly.

A linear root water uptake model

A model of root water extraction is proposed, in which a linear variation of extraction rate with depth is assumed. Five crops are chosen for simulation studies of the model, and soil moisture depletion under optimal conditions from different layers for each crop is calculated. Similar calculations are also made using the constant extraction rate model. Rooting depth is assumed to vary linearly with potential evapotranspiration for each crop during the vegetative phase. The calculated depletion patterns are compared with measured mean depletion patterns for each crop. It is shown that the constant extraction rate model results in large errors in the prediction of soil moisture depletion, while the proposed linear extraction rate model gives satisfactory results. Hypothetical depletion patterns predicted by the model in combination with a moisture tension-dependent sink term developed elsewhere are indicated.

Stochastic modelling of monthly rainfall

A transformation is suggested which can transform a non-Gaussian monthly hydrological time series into a Gaussian one. The suggested approach is verified with data of ten Indian rainfall time series. Incidentally, it is observed that once the deterministic trends are removed, the transformation leads to an uncorrelated process for monthly rainfall. The procedure for normalization is general enough in that it should be also applicable to river discharges. This is verified to a limited extent by considering data of two Indian river discharges.

Depressed weir with two unequal sheetpiles in anisotropic soil

An analytical solution is presented, making use of the Schwartz-Christoffel transformation, for determining the seepage characteristics for the problem of flow under a weir having two unequal sheetpiles at the ends and embedded in an anisotropic porous medium of finite thickness. Results for several particular cases of simple hydraulic structures can be obtained from the general solution presented. Numerical results in nondimensional form have been given for quantity of seepage and exit gradient distribution for various conditions in the equivalent transformed isotropic section and, by making use of the physical parameters in the actual anisotropic plane and the set of transformation relations given, these quantities (seepage loss, exit gradient) can be interpreted in the actual anisotropic physical plane.

Classification of Indian meteorological stations using cluster and fuzzy cluster analysis, and Kohonen artificial neural networks

The present study deals with the application of cluster analysis, Fuzzy Cluster Analysis (FCA) and Kohonen Artificial Neural Networks (KANN) methods for classification of 159 meteorological stations in India into meteorologically homogeneous groups. Eight parameters, namely latitude, longitude, elevation, average temperature, humidity, wind speed, sunshine hours and solar radiation, are considered as the classification criteria for grouping. The optimal number of groups is determined as 14 based on the Davies-Bouldin index approach. It is observed that the FCA approach performed better than the other two methodologies for the present study.

Hybrid moving block bootstrap for stochastic simulation of multi-site multi-season streamflows

The Hybrid approach introduced by the authors for at-site modeling of annual and periodic streamflows in earlier works is extended to simulate multi-site multi-season streamflows. It bears significance in integrated river basin planning studies. This hybrid model involves: (i) partial pre-whitening of standardized multi-season streamflows at each site using a parsimonious linear periodic model; (ii) contemporaneous resampling of the resulting residuals with an appropriate block size, using moving block bootstrap (non-parametric, NP) technique; and (iii) post-blackening the bootstrapped innovation series at each site, by adding the corresponding parametric model component for the site, to obtain generated streamflows at each of the sites. It gains significantly by effectively utilizing the merits of both parametric and NP models. It is able to reproduce various statistics, including the dependence relationships at both spatial and temporal levels without using any normalizing transformations and/or adjustment procedures. The potential of the hybrid model in reproducing a wide variety of statistics including the run characteristics, is demonstrated through an application for multi-site streamflow generation in the Upper Cauvery river basin, Southern India. (C) 2004 Elsevier B.V. All rights reserved.

Mahanadi streamflow: climate change impact assessment and adaptive strategies

Impacts of climate change on hydrology are assessed by downscaling large scale general circulation model (GCM) outputs of climate variables to local scale hydrologic variables. This modelling approach is characterized by uncertainties resulting from the use of different models, different scenarios, etc. Modelling uncertainty in climate change impact assessment includes assigning weights to GCMs and scenarios, based on their performances, and providing weighted mean projection for the future. This projection is further used for water resources planning and adaptation to combat the adverse impacts of climate change. The present article summarizes the recent published work of the authors on uncertainty modelling and development of adaptation strategies to climate change for the Mahanadi river in India.

Detection of moisture stress by indicator plots — A simulation

The possibility of advanced indication of moisture stress in a crop by small prepared plots with compacted or partially sand-substituted soils is examined by an analytical simulation. A series of soils and three crops are considered for the simulation. The moisture characteristics of the soils are calculated with an available model. Using average potential evapotranspiration values and a simple actual evapotranspiration model, the onset of moisture stress in the natural and indicator plots is calculated for different degrees of sand substitution and compaction. Cases where sand substitution fails are determined. The effect of intervening rainfall and limited root depth on the beginning of moisture stress is investigated.

Ecology of the swampy relic forests of Kathalekan from central Western Ghats, India

Introduction of agriculture three millennia ago in Peninsular India’s Western Ghats altered substantially ancient tropical forests. Early agricultural communities, nevertheless, strived to attain symbiotic harmony with nature as evident from prevalence of numerous sacred groves, patches of primeval forests sheltering biodiversity and hydrology. Groves enhanced heterogeneity of landscapes involving elements of successional forests and savannas favouring rich wildlife. A 2.25 km2 area of relic forest was studied at Kathalekan in Central Western Ghats. Interspersed with streams studded with Myristica swamps and blended sparingly with shifting cultivation fallows, Kathalekan is a prominent northernmost relic of southern Western Ghat vegetation. Trees like Syzygium travancoricum (Critically Endangered), Myristica magnifica (Endangered) and Gymnacranthera canarica (Vulnerable) and recently reported Semecarpus kathalekanensis, are exclusive to stream/swamp forest (SSF). SSF and non-stream/swamp forest (NSSF) were studied using 18 transects covering 3.6 ha. Dipterocarpaceae, its members seldom transgressing tropical rain forests, dominate SSF (21% of trees) and NSSF (27%). The ancient Myristicaceae ranks high in tree population (19% in SSF and 8% in NSSF). Shannon-Weiner diversity for trees is higher (>3) in six NSSF transects compared to SSF (<3). Higher tree endemism (45%), total endemic tree population (71%) and significantly higher above ground biomass (349 t/ha) cum carbon sequestration potential (131 t/ha) characterizes SSF. Faunal richness is evident from amphibians (35 species - 26 endemics, 11 in IUCN Red List). This study emphasizes the need for bringing to light more of relic forests for their biodiversity, carbon sequestration and hydrology. The lives of marginal farmers and forest tribes can be uplifted through partnership in carbon credits, by involving them in mitigating global climatic change through conservation and restoration of high biomass watershed forests.

Estimation of variables at ungaged locations by empirical orthogonal functions

A technique based on empirical orthogonal functions is used to estimate hydrologic time-series variables at ungaged locations. The technique is applied to estimate daily and monthly rainfall, temperature and runoff values. The accuracy of the method is tested by application to locations where data are available. The second-order characteristics of the estimated data are compared with those of the observed data. The results indicate that the method is quick and accurate.

Analysis of an Aquifer-water table aquitard system

The problem of pumping an aquifer in an aquifer-water table aquitard system is considered, accounting for the elastic properties of both the aquifer and the aquitard, the gravity drainage in the aquitard and treating the water table as an unknown boundary. The coupled partial differential equations are nondimensionalised, yielding three principal parameters governing the problem. The numerical solution of these equations is obtained for a wide range of parameter values. Type curves are generated and their use is illustrated through a field application.

Analysis of flow near a dug well in an unconfined aquifer

A numerical analysis of flow to a dug well in an unconfined aquifer is made, taking into account well storage, elastic storage release, gravity drainage, anisotropy, partial penetration, vertical flow and seepage surface at the well face, and treating the water table in the aquifer and water level in the well as unknown boundaries. The pumped discharge is maintained constant. The solution is obtained by a two-level iterative scheme. The effects of governing parameters on the drawdown, development of seepage surface and contribution from aquifer flow to the total discharge are discussed. The degree of anisotropy and partial penetration are found to be the parameters which affect the flow characteristics most significantly. The effect of anisotropy on the development of seepage surface is very pronounced.

Meridional propagation of large-scale monsoon convective zones

Observational studies indicate that the convective activity of the monsoon systems undergo intraseasonal variations with multi-week time scales. The zone of maximum monsoon convection exhibits substantial transient behavior with successive propagating from the North Indian Ocean to the heated continent. Over South Asia the zone achieves its maximum intensity. These propagations may extend over 3000 km in latitude and perhaps twice the distance in longitude and remain as coherent entities for periods greater than 2-3 weeks. Attempts to explain this phenomena using simple ocean-atmosphere models of the monsoon system had concluded that the interactive ground hydrology so modifies the total heating of the atmosphere that a steady state solution is not possible, thus promoting lateral propagation. That is, the ground hydrology forces the total heating of the atmosphere and the vertical velocity to be slightly out of phase, causing a migration of the convection towards the region of maximum heating. Whereas the lateral scale of the variations produced by the Webster (1983) model were essentially correct, they occurred at twice the frequency of the observed events and were formed near the coastal margin, rather than over the ocean. Webster's (1983) model used to pose the theories was deficient in a number of aspects. Particularly, both the ground moisture content and the thermal inertia of the model were severely underestimated. At the same time, the sea surface temperatures produced by the model between the equator and the model's land-sea boundary were far too cool. Both the atmosphere and the ocean model were modified to include a better hydrological cycle and ocean structure. The convective events produced by the modified model possessed the observed frequency and were generated well south of the coastline. The improved simulation of monsoon variability allowed the hydrological cycle feedback to be generalized. It was found that monsoon variability was constrained to lie within the bounds of a positive gradient of a convective intensity potential (I). The function depends primarily on the surface temperature, the availability of moisture and the stability of the lower atmosphere which varies very slowly on the time scale of months. The oscillations of the monsoon perturb the mean convective intensity potential causing local enhancements of the gradient. These perturbations are caused by the hydrological feedbacks, discussed above, or by the modification of the air-sea fluxes caused by variations of the low level wind during convective events. The final result is the slow northward propagation of convection within an even slower convective regime. The ECMWF analyses show very similar behavior of the convective intensity potential. Although it is considered premature to use the model to conduct simulations of the African monsoon system, the ECMWF analysis indicates similar behavior in the convective intensity potential suggesting, at least, that the same processes control the low frequency structure of the African monsoon. The implications of the hypotheses on numerical weather prediction of monsoon phenomenon are discussed.