Regolith mass balance inferred from combined mineralogical, geochemical and geophysical studies: Mule Hole gneissic watershed, South India

The aim of this study is to propose a method to assess the long-term chemical weathering mass balance for a regolith developed on a heterogeneous silicate substratum at the small experimental watershed scale by adopting a combined approach of geophysics, geochemistry and mineralogy. We initiated in 2003 a study of the steep climatic gradient and associated geomorphologic features of the edge of the rifted continental passive margin of the Karnataka Plateau, Peninsular India. In the transition sub-humid zone of this climatic gradient we have studied the pristine forested small watershed of Mule Hole (4.3 km(2)) mainly developed on gneissic substratum. Mineralogical, geochemical and geophysical investigations were carried out (i) in characteristic red soil profiles and (ii) in boreholes up to 60 m deep in order to take into account the effect of the weathering mantle roots. In addition, 12 Electrical Resistivity Tomography profiles (ERT), with an investigation depth of 30 m, were generated at the watershed scale to spatially characterize the information gathered in boreholes and soil profiles. The location of the ERT profiles is based on a previous electromagnetic survey, with an investigation depth of about 6 m. The soil cover thickness was inferred from the electromagnetic survey combined with a geological/pedological survey. Taking into account the parent rock heterogeneity, the degree of weathering of each of the regolith samples has been defined using both the mineralogical composition and the geochemical indices (Loss on Ignition, Weathering Index of Parker, Chemical Index of Alteration). Comparing these indices with electrical resistivity logs, it has been found that a value of 400 Ohm m delineates clearly the parent rocks and the weathered materials, Then the 12 inverted ERT profiles were constrained with this value after verifying the uncertainty due to the inversion procedure. Synthetic models based on the field data were used for this purpose. The estimated average regolith thickness at the watershed scale is 17.2 m, including 15.2 m of saprolite and 2 m of soil cover. Finally, using these estimations of the thicknesses, the long-term mass balance is calculated for the average gneiss-derived saprolite and red soil. In the saprolite, the open-system mass-transport function T indicates that all the major elements except Ca are depleted. The chlorite and biotite crystals, the chief sources for Mg (95%), Fe (84%), Mn (86%) and K (57%, biotite only), are the first to undergo weathering and the oligoclase crystals are relatively intact within the saprolite with a loss of only 18%. The Ca accumulation can be attributed to the precipitation of CaCO3 from the percolating solution due to the current and/or the paleoclimatic conditions. Overall, the most important losses occur for Si, Mg and Na with -286 x 10(6) mol/ha (62% of the total mass loss), -67 x 10(6) mol/ha (15% of the total mass loss) and -39 x 10(6) mol/ha (9% of the total mass loss), respectively. Al, Fe and K account for 7%, 4% and 3% of the total mass loss, respectively. In the red soil profiles, the open-system mass-transport functions point out that all major elements except Mn are depleted. Most of the oligoclase crystals have broken down with a loss of 90%. The most important losses occur for Si, Na and Mg with -55 x 10(6) mol/ha (47% of the total mass loss), -22 x 10(6) mol/ha (19% of the total mass loss) and -16 x 10(6) mol/ha (14% of the total mass loss), respectively. Ca, Al, K and Fe account for 8%, 6%, 4% and 2% of the total mass loss, respectively. Overall these findings confirm the immaturity of the saprolite at the watershed scale. The soil profiles are more evolved than saprolite but still contain primary minerals that can further undergo weathering and hence consume atmospheric CO2.

Non-linear regression model for spatial variation in precipitation chemistry for South India

Chemical composition of rainwater changes from sea to inland under the influence of several major factors - topographic location of area, its distance from sea, annual rainfall. A model is developed here to quantify the variation in precipitation chemistry under the influence of inland distance and rainfall amount. Various sites in India categorized as 'urban', 'suburban' and 'rural' have been considered for model development. pH, HCO3, NO3 and Mg do not change much from coast to inland while, SO4 and Ca change is subjected to local emissions. Cl and Na originate solely from sea salinity and are the chemistry parameters in the model. Non-linear multiple regressions performed for the various categories revealed that both rainfall amount and precipitation chemistry obeyed a power law reduction with distance from sea. Cl and Na decrease rapidly for the first 100 km distance from sea, then decrease marginally for the next 100 km, and later stabilize. Regression parameters estimated for different cases were found to be consistent (R-2 similar to 0.8). Variation in one of the parameters accounted for urbanization. Model was validated using data points from the southern peninsular region of the country. Estimates are found to be within 99.9% confidence interval. Finally, this relationship between the three parameters - rainfall amount, coastline distance, and concentration (in terms of Cl and Na) was validated with experiments conducted in a small experimental watershed in the south-west India. Chemistry estimated using the model was in good correlation with observed values with a relative error of similar to 5%. Monthly variation in the chemistry is predicted from a downscaling model and then compared with the observed data. Hence, the model developed for rain chemistry is useful in estimating the concentrations at different spatio-temporal scales and is especially applicable for south-west region of India. (C) 2008 Elsevier Ltd. All rights reserved.

ESTIMATION OF EVAPOTRANSPIRATION IN AN EXPERIMENTAL WATERSHED IN A DECIDUOUS FOREST OF SOUTH INDIA

Forested areas play a dominant role in the global hydrological cycle. Evapotranspiration is a dominant component most of the time catching up with the rainfall. Though there are sophisticated methods which are available for its estimation, a simple reliable tool is needed so that a good budgeting could be made. Studies have established that evapotranspiration in forested areas is much higher than in agricultural areas. Latitude, type of forests, climate and geological characteristics also add to the complexity of its estimation. Few studies have compared different methods of evapotranspiration on forested watersheds in semi arid tropical forests. In this paper a comparative study of different methods of estimation of evapotranspiration is made with reference to the actual measurements made using all parameter climatological station data of a small deciduous forested watershed of Mulehole (area of 4.5 km2 ), South India. Potential evapotranspiration (ETo) was calculated using ten physically based and empirical methods. Actual evapotranspiration (AET) has been calculated through computation of water balance through SWAT model. The Penman-Montieth method has been used as a benchmark to compare the estimates arrived at using various methods. The AET calculated shows good agreement with the curve for evapotranspiration for forests worldwide. Error estimates have been made with respect to Penman-Montieth method. This study could give an idea of the errors involved whenever methods with limited data are used and also show the use indirect methods in estimation of Evapotranspiration which is more suitable for regional scale studies.

Water balance modelling in a tropical watershed under deciduous forest (Mule Hole, India): Regolith matric storage buffers the groundwater recharge process

Accurate estimations of water balance are needed in semi-arid and sub-humid tropical regions, where water resources are scarce compared to water demand. Evapotranspiration plays a major role in this context, and the difficulty to quantify it precisely leads to major uncertainties in the groundwater recharge assessment, especially in forested catchments. In this paper, we propose to assess the importance of deep unsaturated regolith and water uptake by deep tree roots on the groundwater recharge process by using a lumped conceptual model (COMFORT). The model is calibrated using a 5 year hydrological monitoring of an experimental watershed under dry deciduous forest in South India (Mule Hole watershed). The model was able to simulate the stream discharge as well as the contrasted behaviour of groundwater table along the hillslope. Water balance simulated for a 32 year climatic time series displayed a large year-to-year variability, with alternance of dry and wet phases with a time period of approximately 14 years. On an average, input by the rainfall was 1090 mm year(-1) and the evapotranspiration was about 900 mm year(-1) out of which 100 mm year(-1) was uptake from the deep saprolite horizons. The stream flow was 100 mm year(-1) while the groundwater underflow was 80 mm year(-1). The simulation results suggest that (i) deciduous trees can uptake a significant amount of water from the deep regolith, (ii) this uptake, combined with the spatial variability of regolith depth, can account for the variable lag time between drainage events and groundwater rise observed for the different piezometers and (iii) water table response to recharge is buffered due to the long vertical travel time through the deep vadose zone, which constitutes a major water reservoir. This study stresses the importance of long term observations for the understanding of hydrological processes in tropical forested ecosystems. (C) 2009 Elsevier B.V. All rights reserved.

Using a structural approach to identify relationships between soil and erosion in a semi-humid forested area, South India

Biogeochemical and hydrological cycles are currently studied on a small experimental forested watershed (4.5 km(2)) in the semi-humid South India. This paper presents one of the first data referring to the distribution and dynamics of a widespread red soil (Ferralsols and Chromic Luvisols) and black soil (Vertisols and Vertic intergrades) cover, and its possible relationship with the recent development of the erosion process. The soil map was established from the observation of isolated soil profiles and toposequences, and surveys of soil electromagnetic conductivity (EM31, Geonics Ltd), lithology and vegetation. The distribution of the different parts of the soil cover in relation to each other was used to establish the dynamics and chronological order of formation. Results indicate that both topography and lithology (gneiss and amphibolite) have influenced the distribution of the soils. At the downslope, the following parts of the soil covers were distinguished: i) red soil system, ii) black soil system, iii) bleached horizon at the top of the black soil and iv) bleached sandy saprolite at the base of the black soil. The red soil is currently transforming into black soil and the transformation front is moving upslope. In the bottom part of the slope, the chronology appears to be the following: black soil > bleached horizon at the top of the black soil > streambed > bleached horizon below the black soil. It appears that the development of the drainage network is a recent process, which was guided by the presence of thin black soil with a vertic horizon less than 2 in deep. Three distinctive types of erosional landforms have been identified: 1. rotational slips (Type 1); 2. a seepage erosion (Type 2) at the top of the black soil profile; 3. A combination of earthflow and sliding in the non-cohesive saprolite of the gneiss occurs at midslope (Type 3). Types 1 and 2 erosion are mainly occurring downslope and are always located at the intersection between the streambed and the red soil-black soil contact. Neutron probe monitoring, along an area vulnerable to erosion types 1 and 2, indicates that rotational slips are caused by a temporary watertable at the base of the black soil and within the sandy bleached saprolite, which behaves as a plane of weakness. The watertable is induced by the ephemeral watercourse. Erosion type 2 is caused by seepage of a perched watertable, which occurs after swelling and closing of the cracks of the vertic clay horizon and within a light textured and bleached horizon at the top of black soil. Type 3 erosion is not related to the red soil-black soil system but is caused by the seasonal seepage of saturated throughflow in the sandy saprolite of the gneiss occurring at midslope. (c) 2006 Elsevier B.V. All rights reserved.

Ferrolysis induced soil transformation by natural drainage in Vertisols of sub-humid South India

In sub-humid South India, recent studies have shown that black soil areas (Vertisols and vertic Intergrades), located on flat valley bottoms, have been rejuvenated through the incision of streambeds, inducing changes in the pedoclimate and soil transformation. Joint pedological, geochemical and geophysical investigations were performed in order to better understand the ongoing processes and their contribution to the chemistry of local rivers. The seasonal rainfall causes cycles of oxidation and reduction in a perched watertable at the base of the black soil, while the reduced solutions are exported through a loamy sand network. This framework favours a ferrolysis process, which causes low base saturation and protonation of clay, leading to the weathering of 2:1 then 1:1 clay minerals. Maximum weathering conditions occur at the very end of the wet season, just before disappearance of the perched watertable. Therefore, the by-products of soil transformation are partially drained off and calcareous nodules, then further downslope, amorphous silica precipitate upon soil dehydration. The ferrolysed area is fringing the drainage system indicating that its development has been induced by the streambed incision. The distribution of C-14 ages of CaCO3 nodules suggests that the ferrolysis process started during the late Holocene, only about 2 kyr B.P. at the studied site and about 5 kyr B.P. at the watershed outlet. The results of this study are applied to an assessment of the physical erosion rate (4.8x10(-3) m/kyr) since the recent reactivation of the erosion process. (C) 2010 Elsevier B.V. All rights reserved.

Modelling the chemical weathering fluxes at the watershed scale in the Tropics (Mule Hole, South India): Relative contribution of the smectite/kaolinite assemblage versus primary minerals

We investigate the chemical weathering processes and fluxes in a small experimental watershed (SEW) through a modelling approach. The study site is the Mule Hole SEW developed on a gneissic basement located in the climatic gradient of the Western Ghats, South India. The model couples a lumped hydrological model simulating the water budget at the watershed scale to the WITCH model estimating the dissolution/precipitation rates of minerals using laboratory kinetic laws. Forcing functions and parameters of the simulation are defined by the field data. The coupled model is calibrated with stream and groundwater compositions through the testing of a large range of smectite solubility and abundance in the soil horizons. We found that, despite the low abundance of smectite in the dominant soil type of the watershed (4 vol.%), their net dissolution provides 75% of the export of dissolved silica, while primary silicate mineral dissolution releases only 15% of this flux. Overall, smectites (modelled as montmorillonites) are not stable under the present day climatic conditions. Furthermore, the dissolution of trace carbonates in the saprolitic horizon provides 50% of the calcium export at the watershed scale. Modelling results show the contrasted behavior of the two main soil types of the watershed: red soils (88% of the surface) are provider of calcium, while black soils (smectite-rich and characterized by a lower drainage) consumes calcium through overall carbonate precipitation. Our model results stress the key role played by minor/accessory minerals and by the thermodynamic properties of smectite minerals, and by the drainage of the weathering profiles on the weathering budget of a tropical watershed. (C) 2010 Elsevier B.V. All rights reserved.

Formation and preservation of pedogenic carbonates in South India, links with paleo-monsoon and pedological conditions: Clues from Sr isotopes, U-Th series and REEs

The influence of the pedogenic and climatic contexts on the formation and preservation of pedogenic carbonates in a climosequence in the Western Ghats (Karnataka Plateau, South West India) has been studied. Along the climosequence, the current mean annual rainfall (MAR) varies within a 80 km transect from 6000 mm at the edge of the Plateau to 500 mm inland. Pedogenic carbonates occur in the MAR range of 500-1200 mm. In the semi-arid zone (MAR: 500-900 mm), carbonates occur (i) as rhick hardpan calcretes on pediment slopes and (ii) as nodular horizons in polygenic black soils (i.e. vertisols). In the sub-humid zone (MAR: 900-1500 mm), pedogenic carbonates are disseminated in the black soil matrices either as loose, irregular and friable nodules of millimetric size or as indurated botryoidal nodules of centimetric to pluricentimetric size. They also occur at the top layers of the saprolite either as disseminated pluricentimetric indurated nodules or carbonate-cemented lumps of centimetric to decimetric size. Chemical and isotopic (Sr-87/Sr-86) compositions of the carbonate fraction were determined after leaching with 0.25 N HCl. The corresponding residual fractions containing both primary minerals and authigenic clays were digested separately and analyzed. The trend defined by the Sr-87/Sr-86 signatures of both labile carbonate fractions and corresponding residual fractions indicates that a part of the labile carbonate fraction is genetically linked to the local soil composition. Considering the residual fraction of each sample as the most likely lithogenic source of Ca in carbonates, it is estimated that from 24% to 82% (55% on average) of Ca is derived from local bedrock weathering, leading to a consumption of an equivalent proportion of atmospheric CO2. These values indicate that climatic conditions were humid enough to allow silicate weathering: MAR at the time of carbonate formation likely ranged from 400 to 700 mm, which is 2- to 3-fold less than the current MAR at these locations. The Sr, U and Mg contents and the (U-234/U-238) activity ratio in the labile carbonate fraction help to understand the conditions of carbonate formation. The relatively high concentrations of Sr, U and Mg in black soil carbonates may indicate fast growth and accumulation compared to carbonates in saprolite, possibly due to a better confinement of the pore waters which is supported by their high (U-234/U-238) signatures, and/or to higher content of dissolved carbonates in the pore waters. The occurrence of Ce, Mn and Fe oxides in the cracks of carbonate reflects the existence of relatively humid periods after carbonate formation. The carbonate ages determined by the U-Th method range from 1.33 +/- 0.84 kyr to 7.5 +/- 2.7 kyr and to a cluster of five ages around 20 kyr, i.e. the Last Glacial Maximum period. The young occurrences are only located in the black soils, which therefore constitute sensitive environments for trapping and retaining atmospheric CO2 even on short time scales. The maximum age of carbonates depends on their location in the climatic gradient: from about 20 kyr for centimetric nodules at Mule Hole (MAR = 1100 mm/yr) to 200 kyr for the calcrete at Gundlupet (MAR = 700 mm/yr, Durand et al., 2007). The intensity of rainfall during wet periods would indeed control the lifetime of pedogenic carbonates and thus the duration of inorganic carbon storage in soils. (C) 2010 Elsevier Ltd. All rights reserved.

Comparison of electrical resistivity by geophysical method and neutron probe logging for soil moisture monitoring in a forested watershed

Geophysical methods are becoming more popular nowadays in the field of hydrology due to their time and space efficiency. So an attempt has been made here to relate electrical resistivity with soil moisture content in the field. The experiments were carried out in an experimental watershed `Mulehole' in southern India, which is a forested watershed with approximately 80% red soil. Five auger holes were drilled to perform the soil moisture and electrical resistivity measurements in a toposequence having red and black soils, with sandy weathered soil at the bottom. Soil moisture was measured using neutron probe and electrical resistivity was measured using electrical logging tool. The results indicate that electrical resistivity measurements can be used to measure soil moisture content for red soils only.

Inverse models to analyze the spatiotemporal variations of chemical weathering fluxes in a granito-gneissic watershed: Mule Hole, South India

Water-rock reactions are driven by the influx of water, which are out of equilibrium with the mineral assemblage in the rock. Here a mass balance approach is adopted to quantify these reactions. Based on field experiments carried out in a granito-gneissic small experimental watershed (SEW), Mule Hole SEW (similar to 4.5 km(2)), quartz, oligoclase, sericite, epidote and chlorite are identified as the basic primary minerals while kaolinite, goethite and smectite are identified as the secondary minerals. Observed groundwater chemistry is used to determine the weathering rates, in terms of `Mass Transfer Coefficients' (MTCs), of both primary and secondary minerals. Weathering rates for primary and secondary minerals are quantified in two steps. In the first step, top red soil is analyzed considering precipitation chemistry as initial phase and water chemistry of seepage flow as final phase. In the second step, minerals present in the saprolite layer are analyzed considering groundwater chemistry as the output phase. Weathering rates thus obtained are converted into weathering fluxes (Q(weathering)) using the recharge quantity. Spatial variability in the mineralogy observed among the thirteen wells of Mule Hole SEW is observed to be reflected in the MTC results and thus in the weathering fluxes. Weathering rates of the minerals in this silicate system varied from few 10 mu mol/L (in case of biotite) to 1000 s of micromoles per liter (calcite). Similarly, fluxes of biotite are observed to be least (7 +/- 5 mol/ha/yr) while those of calcite are highest (1265 791 mol/ha/yr). Further, the fluxes determined annually for all the minerals are observed to be within the bandwidth of the standard deviation of these fluxes. Variations in these annual fluxes are indicating the variations in the precipitation. Hence, the standard deviation indicated the temporal variations in the fluxes, which might be due to the variations in the annual rainfall. Thus, the methodology adopted defines an inverse way of determining weathering fluxes, which mainly contribute to the groundwater concentration. (C) 2011 Elsevier B.V. All rights reserved.

Controls on intense silicate weathering in a tropical river, southwestern India

The Silicate Weathering Rate (SWR) and associated Carbon dioxide Consumption Rate (CCR) in tropical silicate terrain is assessed through a study of the major ion chemistry in a small west flowing river of Peninsular India, the Nethravati River. The specific features of the river basin are high mean annual rainfall and temperature, high runoff and a Precambrian basement composed of granitic-gneiss, charnockite and minor metasediments. The water samples (n = 56) were collected from three locations along the Nethravati River and from two of its tributaries over a period of twelve months. Chemical Weathering Rate (CWR) for the entire watershed is calculated by applying rainwater correction using river chloride as a tracer. Chemical Weathering Rate in the Nethravati watershed is estimated to 44 t.km(-2).y(-1) encompassing a SWR of 42 t.km(-2).y(-1) and a maximum carbonate contribution of 2 t.km(-2).y(-1). This SWR is among the highest reported for granito-gneissic terrains. The assessed CCR is 2.9 . 10(5) mol.km(-2).y(-1). The weathering index (Re). calculated from molecular ratios of dissolved cations and silica in the river, suggests an intense silicate weathering leading to kaolinite-gibbsite precipitation in the weathering covers. The intense SWR and CCR could be due to the combination of high runoff and temperature along with the thickness and nature of the weathering cover. The comparison of silicate weathering fluxes with other watersheds reveals that under similar morpho-climatic settings basalt weathering would be 2.5 times higher than the granite-gneissic rocks. (C) 2012 Elsevier B.V. All rights reserved.