Implications of Kali-Hindon inter-stream aquifer water balance for groundwater management in western Uttar Pradesh

The Kali-Hindon inter-stream region extends over an area of 395 km² within the Ganga-Yamuna interfluve. It is a fertile tract for sugarcane cultivation. Groundwater is a primary resource for irrigation and industrial purposes. In recent years, over-exploitation has resulted in an adverse impact on the groundwater regime. In this study, an attempt has been made to calculate a water balance for the Kali-Hindon inter-stream region. Various inflows and outflows to and from the aquifer have been calculated. The recharge due to rainfall and other recharge parameters such as horizontal inflow, irrigation return flow and canal seepage were also evaluated. Groundwater withdrawals, evaporation from the water table, discharge from the aquifer to rivers and horizontal subsurface outflows were also estimated. The results show that total recharge into the system is 148.72 million cubic metres (Mcum), whereas the total discharge is 161.06 Mcum, leaving a deficit balance of −12.34 Mcum. Similarly, the groundwater balance was evaluated for the successive four years. The result shows that the groundwater balance is highly sensitive to variation in rainfall followed by draft through pumpage. The depths to water level are shallow in the canal-irrigated northern part of the basin and deeper in the southern part. The pre-monsoon and post-monsoon water levels range from 4.6 to 17.7 m below ground level (bgl) and from 3.5 to 16.5 m bgl respectively. It is concluded that the groundwater may be pumped in the canal-irrigated northern part, while withdrawals may be restricted to the southern portion of the basin, where intense abstraction has led to rapidly falling water table levels.     

Use of sinkhole and specific capacity distributions to assess vertical gradients in a karst aquifer

The carbonate-rock aquifer in the Great Valley, West Virginia, USA, was evaluated using a database of 687 sinkholes and 350 specific capacity tests to assess structural, lithologic, and topographic influences on the groundwater flow system. The enhanced permeability of the aquifer is characterized in part by the many sinkholes, springs, and solutionally enlarged fractures throughout the valley. Yet, vertical components of subsurface flow in this highly heterogeneous aquifer are currently not well understood. To address this problem, this study examines the apparent relation between geologic features of the aquifer and two spatial indices of enhanced permeability attributed to aquifer karstification: (1) the distribution of sinkholes and (2) the occurrence of wells with relatively high specific capacity. Statistical results indicate that sinkholes (funnel and collapse) occur primarily along cleavage and bedding planes parallel to subparallel to strike where lateral or downward vertical gradients are highest. Conversely, high specific capacity values are common along prominent joints perpendicular or oblique to strike. The similarity of the latter distribution to that of springs suggests these fractures are areas of upward-convergent flow. These differences between sinkhole and high specific capacity distributions suggest vertical flow components are primarily controlled by the orientation of geologic structure and associated subsurface fracturing.     

The double solid reactant method: II. An application to the shallow groundwaters of the Porto Plain,Vulcano Island (Italy)

This paper presents an example of application of the double solid reactant method (DSRM) of Accornero and Marini (Environmental Geology, 2007a), an effective way for modeling the fate of several dissolved trace elements during water–rock interaction. The EQ3/6 software package was used for simulating the irreversible water–rock mass transfer accompanying the generation of the groundwaters of the Porto Plain shallow aquifer, starting from a degassed diluted crateric steam condensate. Reaction path modeling was performed in reaction progress mode and under closed-system conditions. The simulations assumed: (1) bulk dissolution (i.e., without any constraint on the kinetics of dissolution/precipitation reactions) of a single solid phase, a leucite-latitic glass, and (2) precipitation of amorphous silica, barite, alunite, jarosite, anhydrite, kaolinite, a solid mixture of smectites, fluorite, a solid mixture of hydroxides, illite-K, a solid mixture of saponites, a solid mixture of trigonal carbonates and a solid mixture of orthorhombic carbonates. Analytical concentrations of major chemical elements and several trace elements (Cr, Mn, Fe, Ni, Cu, Zn, As, Sr and Ba) in groundwaters were satisfactorily reproduced. In addition to these simulations, similar runs for a rhyolite, a latite and a trachyte permitted to calculate major oxide contents for the authigenic paragenesis which are comparable, to a first approximation, with the corresponding data measured for local altered rocks belonging to the silicic, advanced argillic and intermediate argillic alteration facies. The important role played by both the solid mixture of trigonal carbonates as sequestrator of Mn, Zn, Cu and Ni and the solid mixture of orthorhombic carbonates as scavenger of Sr and Ba is emphasized.

Evaluation of approaches to calculate debris-flow parameters for hazard assessment

Many different runout prediction methods can be applied to estimate the mobility of future debris flows during hazard assessment. The present article reviews the empirical, analytical, simple flow routing and numerical techniques. All these techniques were applied to back-calculate a debris flow, which occurred in 1982 at La Guingueta catchment, in the Eastern Pyrenees. A sensitivity analysis of input parameters was carried out, while special attention was paid to the influence of rheological parameters. We used the Voellmy fluid rheology for our analytical and numerical modelling, since this flow resistance law coincided best with field observations. The simulation results indicated that the “basal” friction coefficients rather affect the runout distance, while the “turbulence” terms mainly influence flow velocity. A comparison of the velocity computed on the fan showed that the analytical model calculated values similar to the numerical ones. The values of our rheological parameters calibrated at La Guingueta agree with data back-calculated for other debris flows. Empirical relationships represent another method to estimate total runout distance. The results confirmed that they contain an important uncertainty and they are strictly valid only for the conditions, which were the basis for their development. With regards to the simple flow routing algorithm, this methods could satisfactorily simulate the total area affected by the 1982 debris flow, but it was not able to directly calculate total runout distance and velocity. Finally, a suggestion on how different runout prediction methods can be applied to generate debris-flow hazard maps is presented. Taking into account the definition of hazard and intensity, the best choice would be to divide the resulting hazard maps into two types: “final hazard maps” and “preliminary hazard maps”. Only the use of numerical models provided final hazard maps, because they could incorporate different event magnitudes and they supplied output-values for intensity calculation. In contrast, empirical relationships and flow routing algorithms, or a combination of both, could be applied to create preliminary hazard maps. The present study only focussed on runout prediction methods. Other necessary tasks to complete the hazard assessment can be looked up in the “Guidelines for landslide susceptibility, hazard and risk zoning” included in this Special Issue.     

Assessment of disaster resilience capacity of hillslope communities with high risk for geological hazards

This study presents a novel preparedness assessment method for assessing hazard mitigation and environmental planning of hillslope communities. A professional questionnaire was utilized to weight each indicator. Communities in Hsinchu, Taichung and Nantou counties with debris flow hazards were taken as study samples. Debris flow risk and landslide susceptibility for each community were determined using Geographic Information System (GIS) technology and logistic regression analysis. Thus, a novel risk assessment method for evaluating disaster resilience capacity of hillslope communities was established. This method was then applied to assess casualties caused by Typhoon Herb in 1996 and Typhoon Mindulle in 2004. Additionally, the analytical results generated by this assessment method were discussed with the aim of developing references for implementation of risk analysis, increasing the effectiveness of disaster mitigation, and reducing future loss of life and property.     

Isotopic variation of fishes in freshwater and estuarine zones of a large subtropical coastal lagoon

We used stable C and N isotope ratios of tissues from 29 fish species from a large subtropical lagoon in southern Brazil to examine spatial variability in isotopic composition and vertical trophic structure across freshwater and estuarine habitats. Nitrogen isotope ratios indicated a smooth gradation in trophic positions among species, with most fishes occupying the secondary and tertiary consumer level. Fish assemblages showed a significant shift in their carbon isotopic signatures between freshwater and estuarine sites. Depleted carbon signatures (from −24.7‰ to −17.8‰) were found in freshwater, whereas more enriched signatures (from −19.1‰ to −12.3‰) were obtained within the estuarine zone downstream. Based on our survey of the C₃ and C₄ plants and isotopic values for phytoplankton and benthic microalgae reported for ecosystems elsewhere, we hypothesized that the observed δ¹³C differences in the fish assemblage between freshwater and estuarine sites is due to a shift from assimilating organic matter ultimately derived from C₃ freshwater marsh vegetation and phytoplankton at the freshwater site (δ¹³C ranging from −25‰ to −19‰), to C₄ salt-marsh (e.g. Spartina) and widgeon grass (Ruppia maritima), benthic microalgae and marine phytoplankton at the estuarine sites (from −18‰ to −12‰). Our results suggested that fish assemblages are generally supported by autochthonous primary production. Freshwater fishes that likely were displaced downstream into the estuary during periods of high freshwater discharge had depleted δ¹³C values that were characteristic of the upper lagoon. These results suggest that spatial foodweb subsidies can occur within the lagoon.     

Numerical Experiment on Two-Dimensional Line Thermal

The time evolution of a two-dimensional line thermal-a turbulent flow produced by an initial element with signifi-cant buoyancy released in a large water body, is numerically studied with the two-equation k - e model for turbulence closure. The numerical results show that the thermal is characterized by a vortex pair flow and a kidney shaped concentra-tion structure with double peak maxima; the computed flow details and scalar mixing characteristics can be described by self-similar relations beyond a dimensionless time around 10. There are two regions in the flow field of a line thermal: a mixing region where the concentration of tracer fluid is high and the flow is turbulent and rotational with a pair of vortex eyes, and an ambient region where the concentration is zero and the flow is potential and well-described by a model of doublet with strength very close to those given by early experimental and analytical studies. The added virtual mass coeffi-cient of the thermal motion is found to be approximat     

秦山核电站邻近水域的流况分析

本文分析了秦山核电站邻近水域各测站的实测流、潮流和佘流分布及季节变化。结果表明:测区内潮流为非正规丰日浅海潮流,其运动形式呈往复流。核电站二期工程取水口附近(L_1,18站)为一大流速区,而一期工程出水口周围(L_2,L_2′站)有一相对小的流速区。大流速区落潮流历时大于涨潮流历时,而相对小流速区则相反;大流速区的余流较小,而小流速区的余流大。在四个季节中,涨潮流速秋季最大,春季次之,夏季最小;落潮流速秋季最大,夏季次之,春季最小。余流冬、秋两季较大。在实测流中以潮流为主,径流和季风的作用也相当显著。     

珠江口沉积物酸挥发性硫化物与重金属生物毒性的研究

研究了2003年10月采于珠江口5个站位的沉积物剖面酸挥发性硫化物(AVS)和同时提取的重金属(SEM:Pb,Zn,Cu,Cd,Ni)。其中,站位1、2位于中滩,其沉积物的AVS含量变化范围较小,为0.25—4.06μmol.g-1;站位3、4位于西滩,其沉积物的AVS含量变化范围较大,为0—26.09μmol.g-1。中滩和西滩沉积物的AVS含量均随深度增加。西滩表层沉积物的AVS含量接近于零,这可能与该水域较强的底层流和沉积物的再悬浮作用有关。站位5位于珠江口外侧,其表层沉积物的AVS含量相对较高,且垂向变化较小,可能是还原性沉积物间歇性再悬浮后重新沉积的结果。站位1、2和5沉积物中同时提取的重金属含量大体在0.95±0.2μmol.g-1范围内,随深度增加略呈下降趋势;而西滩沉积物重金属含量相对较高,为1.43—2.42μmol.g-1,且在一定深度范围内随深度增加呈明显下降的趋势,表明珠江口西滩沉积物中的重金属污染有加重的趋势。对AVS/SEM摩尔比值和单个金属的毒性效应研究显示,珠江口内尤其是西滩的表层沉积物存在重金属污染,对其中生活的底栖生物具有潜在的毒性效应。     

How tides and river flows determine estuarine bathymetries

For strongly tidal, funnel-shaped estuaries, we examine how tides and river flows determine size and shape. We also consider how long it takes for bathymetric adjustment, both to determine whether present-day bathymetry reflects prevailing forcing and how rapidly changes might occur under future forcing scenarios.Starting with the assumption of a 'synchronous' estuary (i.e., where the sea surface slope resulting from the axial gradient in phase of tidal elevation significantly exceeds the gradient in tidal amplitude ), an expression is derived for the slope of the sea bed. Thence, by integration we derive expressions for the axial depth profile and estuarine length, L, as a function of and D, the prescribed depth at the mouth. Calculated values of L are broadly consistent with observations. The synchronous estuary approach enables a number of dynamical parameters to be directly calculated and conveniently illustrated as functions of and D, namely: current amplitude Û, ratio of friction to inertia terms, estuarine length, stratification, saline intrusion length, flushing time, mean suspended sediment concentration and sediment in-fill times.Four separate derivations for the length of saline intrusion, L_I, all indicate a dependency on (U_o is the residual river flow velocity and f is the bed friction coefficient). Likely bathymetries for `mixed' estuaries can be delineated by mapping, against and D, the conditions L<sub>I</sub>/L<1,E<sub>X</sub>/L<1 (E<sub>X</sub> is the tidal excursion) alongside the Simpson-Hunter criteria D/U<sup>3</sup><50 m<sup>−2</sup> s<sup>3</sup>. This zone encompasses 24 out of 25 `randomly' selected UK estuaries.However, the length of saline intrusion in a funnel-shaped estuary is also sensitive to axial location. Observations suggest that this location corresponds to a minimum in landward intrusion of salt. By combining the derived expressions for L and L_I with this latter criterion, an expression is derived relating D_i, the depth at the centre of the intrusion, to the corresponding value of U_o. This expression indicates U_o is always close to 1 cm s⁻¹, as commonly observed. Converting from U_o to river flow, Q, provides a morphological expression linking estuarine depth to Q (with a small dependence on side slope gradients).These dynamical solutions are coupled with further generalised theory related to depth and time-mean, suspended sediment concentrations (as functions of and D). Then, by assuming the transport of fine marine sediments approximates that of a dissolved tracer, the rate of estuarine supply can be determined by combining these derived mean concentrations with estimates of flushing time, F_T, based on L_I. By further assuming that all such sediments are deposited, minimum times for these deposition rates to in-fill estuaries are determined. These times range from a decade for the shortest, shallowest estuaries to upwards of millennia in longer, deeper estuaries with smaller tidal ranges.