Environmental Geology, (2005) 47:(6), 745-750
International Journal of Geosciences     

Original Article                                                        

Radon variations in an active landslide zone along the Pindar River, in Chamoli District, Garhwal Lesser Himalaya, India

V.M. ChoubeyContact Information, S.K. Bartarya1 and R.C. Ramola2

1 Wadia Institute of Himalayan Geology, 33 General Mahadeo Singh Road, Dehradun, 248001, India

2 Physics Department, HNB Garhwal University campus, New Tehri, 249001, India

Email: [email protected] Website: www.geocities.ws/vchoubey

Received: 22 May 2004  Accepted: 13 October 2004  

Abstract  Radon measurements were made in the soil and spring/seepage water in and around an active landslide located along the Pindar river in the Chamoli District of Uttaranchal in Garhwal Lesser Himalaya, to understand the application of radon in geological disasters. The landslide is a compound slide i.e. a slump in the crown portion, and debris slide and fall in the lower part. The bedrock consists of gneisses and schists of the Saryu Formation of the Almora Group of Precambrian age. The presence of several small slump scars and debris slide/fall scars along the length of the slide indicates continuous downward movement. The radon concentrations in the present study are much lower in comparison to values reported from other regions. However, the present radon data show relative variation in the slide zone. The concentration of radon measured in landslide zones varies from 3.1 Bq/l to 18.4 Bq/l in spring water and from 2.3 kBq/m3 to 12.2 kBq/m3 in the soil gas of the debris. Along the section of the slide, the radon values in water and soil are slightly higher in the upper slopes i.e. toward the crown portion of the landslide as compared to the distal portion. The relatively low concentration of radon both in soil gas and water in the toe portion of the landslide may be due to the high porosity of the debris, which does not allow radon to accumulate in the soil and water, whereas, towards the crown portion, the high frequency of fractures increases the surface area due to particle size reduction, and the near absence of debris enhances the radon emanation in soil.

Keywords  Radon - Landslides - Himalaya - Fracture/joints - Slide scar


Introduction

Radon gas is produced by the decay of naturally occurring uranium found in almost all rocks and soils. However, radon concentration varies with specific site and geological material. Especially at faults and active geodynamic zones, radon gas moves upward following the so-called ldquochannelsrdquo and then diffuses in groundwater and aerial zones (Hess et al. 1985; Ball et al. 1991; Choubey et al. 1999, 2001). This movement makes up abnormal fields of radon radiation in groundwater and soil gas (Kemsi et al. 2001). Determination of abnormal fields of radon radiation might enable assessment of the existence of active geodynamic or active seismic zones in an area. The enhanced radon concentration in subsoil in old mining areas may be indicative of subsidence due to fractioning of the underground rocks (Kies et.al. 2004). The profound impact of landslide-related slope instability is particularly observed in areas of active faults and thrusts of the Himalayan region. Rock types and their structural setting and degree of weathering immensely influence the stability of the hill slopes in all the four major litho tectonic belts of the Himalayas (Valdiya 1998). Even periodic seismic activity in the Himalayas is responsible for large-scale slope disturbances in the form of landslides and other mass movements (Sah and Bartarya 2002; Ravindran and Philip 2002). Deep dissection by streams and toe cutting by rivers have often led to local large landslides and other mass movements in the Himalayas (Bartarya et al. 1996).

The studies in the Himalayan region have demonstrated that the tectonic processes, types of rocks and geohydrological characteristics of rock mass control the concentration of radon in soil and water (Choubey et.al. 1998, 1999, 2001). Landslides, which result due to geodynamic processes, have been considered as a source of radon (Purtschller et. al. 1994, 1995). The indoor measurement in Koefels landslide (Tyrol, Austria) shows a high concentration of radon in houses located on old landslide deposits (Ennemoser et.al. 1994). However, enhanced emanation of radon in soil and water in the active landslide zone and enhanced radon emanation as a precursor of landslide are yet to be ascertained. This paper is aimed at analyzing the radon characteristics of soil and water in an active landslide zone and addressing their applicability in the study of geological disasters. Keeping this in mind, the radon measurements were made in the soil and spring water in and around an active landslide zone located along the Pindar River between the Hirmani and Mathera villages in the Chamoli District of Uttaranchal.


Geology of the area

The study area lies in the north central part of the Kumaun Lesser Himalaya. The geology of this area has been studied in detail by various workers (Gansser 1964; Saklani 1972; Agarwal and Kumar 1973; Valdiya 1978, 1980). Geologically, the area belongs to the Almora Group (also known as Almora Crystalline) of Lesser Himalaya (Valdiya 1980). The rock consists of over-thrusted succession of a variety of schists, micaceous quartzite and gneisses belonging to the lower amphibolite facies of regional metamorphism and tectonically emplaced bodies of granite and granodiorites. The Almora Group is further divided into three units, namely Suryu formation (with Champawat granodiorite in the upper part), the Gumalikhet formation and the Munsari formation. The lithology of the landslide zone that was studied belongs to the Saryu formation (Fig. 1). The rocks of the slide zone consist of alternate bands of garnetiferous muscovite schists and gneisses. At places, mylonitized quartz porphyry and weathered micaceous quartzite are also present. Bands of augon gneisses are noticed further east near the Gwaldam area. The rocks of the Almora group are separated from underlying sedimentaries of the Berinag formation of the Jaunsar group by the Munsiari thrust. The dip of the bedrocks of the slide area varies from 20° to 30° in NE and NNE direction. The rocks are characterized by multiple deformations resulting in superimposed folding and repeated faulting and thrusting. The rocks are cut by numerous NW to SE- and NE to SW-trending faults and fractures.

254_2004_1196_f1.jpg

Fig. 1  Geological map of the study area showing the location of the slide


The bedrocks of the slide area have evidently been affected by severe tectonic movements, which have occurred along the Almora thrust. A number of shear planes are developed parallel, oblique and perpendicular to the bedding. The more competent granite gneisses are highly fractured, jointed and sheared, whereas, the less competent schists are extensively sheared and pulverized. The joints and fractures are generally open near the surface and their openness can be partially attributed to the de-compression, which accompanied the valley formation. Prominent sets of joints and fractures are noticed and are dipping in 60–65° towards the NE, 35° towards NW and 75° due ESE direction.


Description of the slide zone

Figure 2 shows the geomorphologic features of the slide zone. The slide is present on the northeast-facing slope of the Pindar River. The rainfall-triggered debris finds its way on the road cut and then goes to the Pindar River. A number of seepage zones and springs are present in the toe, main body and crown portion of the slide (Fig. 2). The seepages were generally confined to the shear plane and along the contact of gneisses and schists. Two seasonal streams are present on either side of the main body of the slide. Runoff, during the heavy rain, flows and seeps down into the slide mass and streambed. This water emerges as seepages and springs in the main body of the slide. The surface runoff and spring water led to the formation of 1-m-wide and 0.5-m-deep gullies in the main body of the slide.

254_2004_1196_f2.jpg

Fig. 2  Map of the landslide zone showing geomorphic features of the landslide zone and location of the samples for radon measurement


The slide measures about 150 m in length and 300 m in width. Several scars of slump and debris fall and slide are present in the crown portion and main body of the slide. The exposed slide scars dip at 60–65° towards the NE direction. There are at least three sets of colluvial benches present in old landslide deposits. Two 1st order streams (the eastern one is lineament controlled) cutting the slide on either side are present in the slide zone. The fall of debris scars the middle portion of the slide near the vertical dipping at 65°NE and about 10 m in length and 8 m in width. Scars caused by the fall of debris in the toe portion of the slide are present along the road and shown as road-cut scar (Fig. 3). Because the whole rock mass laden with debris from the head to toe region has slided down towards the Pindar river, lateral scars of 3 m in height, sloping 60° towards N, 30°E and 65°NW are present on either side of the slided mass.

254_2004_1196_f3.jpg

Fig. 3  Photograph of the landslide showing deris slide scar along the Ist order stream on the western flank and along the road in the toe portion of the slide


Due to the slide, several 2–3-m long and 0.5-m-deep E–W and NW to SE-oriented transverse cracks and 0.3–0.5-m wide and about 1-m-deep N to S-oriented longitudinal cracks have developed in the upper and middle portion of the slide mass during the distress movements. Besides, several semi-circular radial cracks are also present in the crown portion of the slide. These cracks were open near the surface. At places, particularly towards the toe portion, secondary scars of 0.5–1 m height developed subsequently, indicating continued slow downward movement (creeping) of the sliding mass.


Causes and mechanism of the slide

It is a compound slide i.e. a slump in the crown portion and the fall of debris, and slides in the middle and lower portion. The rocks consist of highly fractured and jointed schists and gneisses and old landslide debris. The intense fracturing and joints converted the rocks into loose blocks. Water has played a major role in triggering the movement. Exposed backward-rotated scars and radial cracks in the crown portion in gneisses and schists suggest slumping. The disposition of fractures in joint planes (trending NE-SW, NW-SE, ESE-WNW) and schistosity (30° toward NE) create conditions for sliding, apart from the presence of old landslide debris and slumped mass. The presence of highly sheared, fractured and jointed rocks due to the close proximity of major tectonic plains (north and south Almora thrust) have also contributed to the slope failure. The presence of drainage on either side of the slide and spring and seepages in the main body of the slide indicates continuous percolation of rainwater. The schist rocks acted like a viscous body after absorbing water and facilitating mass movement by lubricating the slip surface. The Pindar River flowing at the base of the slide under cuts the toe of the slide. High rainfall and accumulated pore water pressure (evident from the presence of spring/seepages) over a long period have reduced the shear strength of the highly fractured and jointed rock mass along the dip slope (Fig. 4). The loose land mass slided down due to toe cutting by the Pindar River and continuous road cutting.

254_2004_1196_f4.jpg

Fig. 4  Cross section of the landslide showing the mechanism of the slide


Materials and methods

The radon measurements were made in water (spring and seepage) and soil gas of the debris in the landslide zone. About 750 ml of the water sample was taken in a radon-tight reagent bottle of 1-litre capacity connected in a close circuit with a ZnS-coated detection chamber through a hand-operated rubber pump and a glass bulb containing CaCl2 to absorb the moisture. Air was then circulated in close circuit for a period of ten min till the radon formed a uniform mixture with the air and the resulting alpha activity was recorded.

In soil-radon emanometry, auger holes, each 50 cm in depth and 6 cm in diameter, were left covered for 24 h so that the amount of radon and thoron became stable. The soil-gas probe was fixed in the auger hole and formed an airtight compartment (Ghosh and Bhalla 1981). The rubber pump, soil-gas probe and alpha detector were connected in a close circuit. Air was circulated through a ZnS-coated chamber for a period of 15 min till the radon formed a uniform mixture with the air.

The alpha particles produce scintillations on the sample chamber, which are sensed by the photomultiplier tube and converted to electronic pulses counted on a digital display. Clamping both the ends then isolated the detector and observations were taken after 4 h. A time gap of 4 h is necessary to allow alpha-emitting daughters of radon to come into equilibrium with their own daughters. The resulting number of the alpha counts is converted into Bq/m3 by using the calibration factor 1 count/min=0.11 KBq/m3, determined by Singh et al. (1986) under similar conditions.


Results and discussions

The results of radon measurements in spring water and soil are given in Table 1 and their location in Fig. 2. The radon concentration measured varies from 3.1 Bq/l to 18.4 Bq/l in spring water, whereas in soil-gas samples, it varies from 2.3 kBq/m3 to 12.2 kBq/m3. The average (background) value in soil and water in the surrounding area from similar lithology varies from 0.2 kBq/m3 to 0.9 kBq/m3 and 0.7 Bq/l to 1.2 Bq/l, respectively (Choubey et al. 2004) in preparation). The radon concentration across a cross section of the landslide (Figs. 2, 4) shows that the radon values in soil are slightly higher (4.7–12.2 kBq/m3) in the upper slopes i.e. towards the crown portion of the landslide as compared to the toe portion (2.3 kBq/m3 to 6.3 kBq/m3). The water samples also show a similar trend (15.3 Bq/l to18.3 Bq/l in crown portion and 3.1 Bq/l to 4.2 Bq/l in toe portion, except sample# SW 1). The relatively high concentration of radon (18.4 Bq/l) in sample# SW1 in the toe portion of the slide is possibly because the spring is present close to a lineament-controlled stream on the side of the slide mass. Several fractures and joints are present in the freshly exposed rocks along the stream and may also contribute to the enhanced radon in this sample.

Table 1  Radon concentration in soil gas and water in the land slide area

Sample no.

Location of sample in slide area

Radon concentration

In soil

 SL1

Crown area

10.3+0.83 kBq/m3

 SL2

Crown area

9.1+0.78 kBq/m3

 SL3

Crown area

12.2±0.90 kBq/m3

 SL4

Crown area

11.7±0.88 kBq/m3

 SL4a

Crown area

4.7±0.56 kBq/m3

 SL5

Crown area

12.1±0.89 kBq/m3

 SL6

Toe area

2.8±0.43 kBq/m3

 SL7

Toe area

2.3±0.39 kBq/m3

 SL8

Toe area

6.3±0.65 kBq/m3

 SL9

Toe area

5.7±0.61 kBq/m3

 SL11

Away from the slide

4.1±0.52 kBq/m3

 SL12

Away from the slide

2.3±0.39 kBq/m3

In spring water

 Sw1

Toe area

18.4±1.10 Bq/l

 Sw2

Toe area

3.1±0.45 Bq/l

 Sw3

Toe area

4.2±0.53 Bq/l

 Sw4

Main body of slide

8.1±0.73 Bq/l

 Sw5

Main body of slide

9.8±0.81 Bq/l

 Sw6

Main body of slide

8.4±0.75 Bq/l

 Sw7

Crown area

15.3±1.01 Bq/l

 Sw8

Crown area

18.3±1.10 Bq/l

The radon concentrations (8.1–9.8 Bq/l) in spring water of the middle portion of the slide are in between the radon values of the crown and toe portion of the slide. The radon in soil/debris of the middle portion of the slide could not be measured due to inaccessibility to the site. Towards the crown portion, the increased intensity of fractures and near absence of debris enhanced the radon emanation in soil. The low concentration of radon both in soil and water in the toe portion of the landslide may be due to the high porosity and permeability of the debris, which does not allow radon to accumulate in the soil and water. Schery et al. 1982 observed that the permeability/porosity of the sub-surface material plays an important role in the escape of radon to the atmosphere.

The radon concentration measured in the present landslide is relatively low when compared to other parts of the world (Kemsi et al. 2001), and also, there are not many references related to landslides as possible sources of radon; a comparison with the reported studies of Ennemoser et.al. (1994) and Partscheller et al. (1995) of the Koefels (Tyroll, Austria) and Himang (Nepal) landslide will be a worthy attempt. Lithologically, the present landslides of Chamoli District are also present in the gneisses and schists of the Saryu formation of the Almora Group as that of the Koefels (Tyroll, Austria) and Himang (Nepal) landslides. The higher brittle deformation and shearing of rocks occur near the slide scar. As a result, the rocks in the slide area show fragmentation into blocks due to the high frequency of fracture and joints. The intensity of fracturing in the landslide differs from place to place. Steeply dipping, crosscutting and partly open joints reveal good permeabilities and circulation pathways for soil air. Therefore, the highest radon emanations in the landslide zone occur within the landslide scar area where strongly fractured gneisses and schists are exposed.

Measurements of radon in the landslide studied have shown that higher radon levels are in accordance with the frequency of fractures in augen gneisses and schists exposed in the landslide zone. In the basal part of the sliding plane, the highest degree of fracturing occurs. As a result, the radon concentration in soil gas and water show higher values near the sliding plane. The grain-size reduction due to high shear stress near the slide scar and further weathering of the slided mass expose the uranium-rich mineral (apatite, zircon), thus giving rise to high values of radon concentration in this portion of the slide.

However, Surbeck (1992) considered rockslides as an important factor controlling radon distribution but mainly related to permeabilities and not as a possible source. However, as pointed out by Ennemoser et al. (1993), the emanation rate of radon from soil depends not only on the concentration of uranium and radium content but to some extent on emanating power and the diffusion coefficient for radon in soil. The freshly fallen landslide debris, compared to solid rocks, has an increased emanating power and diffusion coefficient due to intense fracturing, and crushing down to the grain and sub-grain size, resulting in an increased surface area and porosity (Purtschller et al. 1995). Semkov (1990), Sutherland (1994) and Pirchl et al. (1994) have also suggested correlation between specific surface area and emanating powers of radon. Thus, particle-size fractions due to fracturing and crushing close to the landslide scar in the crown area in Garhwal Himalaya have possibly played a major role in increased radon emanation. Between the slide scar and toe of the slide, the distribution of radon takes place through a circulation pathway that exists between slided masses. But the concentration decreases towards the toe portion as radon possibly escapes into the air through highly porous, permeable and pulverized loose colluvial debris


Inference

In our study, radon concentrations are not high when compared to Koefels (Tyroll, Austria) and Himal (Nepal). However, elevated radon concentrations possibly indicate that the landslide scar and failure plane of an active landslide may be the source of enhanced emanation (when compared to mean radon concentration in the surroundings) in an active landslide zone. The enhanced radon emanation is because of the highly active surface area and circulation pathways due to fracturing and crushing of rocks near the landslide scar. However, in the toe portion, increased porosity of the colluvial debris increased the diffusion of radon in the atmosphere, thus resulting in low radon concentrations. More measurements from other active landslides would confirm the above inference.

Acknowledgements  We are thankful to the Director, Wadia Institute of Himalaya Geology, Dehradun, for providing necessary facilities to carry out the present study. The authors (V.M. Choubey and R.C. Ramola) are also thankful to DST, Government of India, for providing the funds in the form of a research project.


References

Agarwal NC, Kumar G (1973) Geology of upper Bhagirathi and Yamuna valleys, Uttarkashi Distt. Kumaun Himalaya. Him Geol 3:1–23

 

Ball TK, Cameron DG, Colman TB, Robert PD (1991) Behaviour of radon in the geological environment: a review. Q J Eng Geol Lond 24:169–182

 

Bartarya SK, Virdi NS, Sah MP (1996) Landslide Hazards: Some case studies from the Satluj valley, Himachal Pradesh. J Him Geol 17:193–207

 

Choubey VM (1998) Radon measurements in soil and water and its relation with geology, Garhwal Himalaya. In: Katase A, Shimo M (eds) Preceeding volume of 7th Tohwa international symposium ldquoRadon and Thoron in the Human Environmentrdquo held at Fukuoka, Japan October 1997. World Scientific Publ, Singapore, pp 193–198

 

Choubey VM, Bist KS, Saini NK, Ramola RC (1999) Relation between soil-gas radon variation and the different lithotectonic units of Garhwal Himalaya, India. J Appl Radiat Isot 51(5):587–592

 

Choubey VM, Bartarya SK, Saini NK, Ramola RC (2001) Impact of geohydrology and neotectonic activity on radon concentration in groundwater of intermontane Doon Valley, Outer Himalaya, India. Environ Geol 40(3):257–266.

 

Choubey VM, Bartarya SK, Ramola RC (2004) Radon variations in soil and water in eastern Uttaranchal, Kumaun Himalaya (communicated).

 

Ennemoser O, Ambach W, Brunner P, Schneider P, Oberaigner W, Purtschller R, Stingl V (1993) Unusually high indoor radon concentrations. Atmos Environ 27A(14):2169–2172

 

Ennemoser O, Ambach W, Brunner P, Schneider P, Oberaigner W, Purtschller F, Stingl V, Keller G (1994) Unusaly high indoor radon concentrations from a giant rockslide. Sci Total Environ 151:235–240.

 

Gansser A (1964) Geology of the Himalaya. Inter Science Publication. Wiley, New York, p 289

 

Ghosh PC, Bhalla NS (1981) On the behaviour and measurement of thoron (220Rn) in soil. Ind J Earth Sci 8:1–9

 

Hess CT, Michel J, Horton TR, Prichard HM, Coniglio WA (1985) The occurrence of radio-activity in public water supplies in the United States. Health Phys 48(5):553–586.

 

Kemsi J, Siehl A, Stegemann R, Valdivia-manchego M (2001) Mapping the geogenic radon potential in Germany. Sci Total Environ 272:217–230.

 

Kies A, Storoni A, Tosheva Z, Hofmann H (2004) Radon measurements as a monitoring possibility for mining subsidence occurrences, Norm IV, Szczyrk, pp 625–632

 

Pirchl T, Sieder G, Brunner P, Purtschller F, Stingl V, Tessadri R (1994) Geochemische, mineralogische und geologische Aspekte zur Radon-Anomalie in Umhausen/Ötztal. Mitt.Österr. Min Ges 139:360–361, Wien

 

Purtschller F, Stingl V, Brunner P, Ennemoser O (1994) The Tsergo Ri Landslide (Langtang Himal)—a case study of radon emanation from giant landslides. J Nepal Geol Soc 10:102–104

 

Purtschller F, Pirchl P, Sieder G, Stingl V, Tessadri T, Brunner P, Ennemoser O, Schneider P (1995) Radon emanation from giant landslides of Koefels (Tyrol, Austria) and Langtang Himal (Nepal). Environ Geol 26:32–38.

 

Ravindran KV and Philip G (2002) Mapping of 29th March 1999 Chamoli earthquake induced landslides using IRS-IC-ID data. Him Geol 23(1&2):69–76

 

Sah MP and Bartarya SK (2002) The impact of March 29th, 1999 Chamoli earthquake on slope stability and spring discharge in Chamoli and Rudraprayag districts of Garhwal Himalaya. Him Geol 23(1&2):121–133

 

Saklani PS (1972) Lithostratigraphy and structure of the area between the Bhagirathi and Bhilangana rivers, Garhwal Himalaya. Him Geol 2:342–355

 

Schery SD, Gaeddert DH, Wilkening MH (1982) Transport of radon from fractured rock. J Geophys Res 87(B4):2969–2976

 

Semkov TM (1990) Recoil-emanation theory applied to radon release from mineral grain. Geochim. Cosmochim Acta. 54:425–440.

 

Singh M, Singh NP, Singh S and Virk HS (1986) Calibration of radon detectors. Nucl Tracks-Radiat Meas 12:739–742.

 

Surbeck H (1992) Das nationale Schweizer Radonprogramm. In: Tagungsber 14th Jahrestagg. OSRAD, Vienna, pp 6–14

 

Sutherland DS (1994) Radon workshop- geology, environment, techniques. Geoscientist 4(2):27–29

 

Valdiya KS (1978) Extension and analague of the Chail nappe in Kumaun. Ind. J Earth Sci 5:1–19

 

Valdiya KS (1980) Geology of Kumaun Lesser Himalaya. WIHG Publ., Dehra Dun, p 291

 

Valdiya KS (1998) Catastrophic landslides in Uttaranchal, Central Himalaya. J Geol Soc Ind 52:483–486