. Introduction

The collision of the Indian and the Eurasian plates resulted in north-south contraction of the lithosphere and development of southward propagating fold and thrust belt (FTB) (Molnar and Tapponnier, 1975; Molnar, 1984; Guha et al., 2007). Himalayas are one of the few longest mountain ranges in the world, comparable with Andes, Rocky Mountains, Great Dividing Range etc. The Himalayan FTB has been divided by northward dipping major boundary thrusts; Main Central Thrust (MCT) separates the crystalline Higher Himalayas from the metamorphosed Lesser Himalayas, Main Boundary Thrust (MBT) separates the Lesser Himalayas from the sedimentary rocks of the Siwaliks and Main Frontal Thrust (MFT) being the southernmost expression of FTB at mountain front separates the Siwaliks from the Quaternary deposits (Gansser, 1964, 1981). Himalayas provide a natural laboratory to the earth scientist trying to understand the balance between uplift, incision and erosion in an actively deforming orogen, resulting in topographic evolution on the 105–107 years timescales (Lavé and Avouac, 2000; Wobus et al., 2005). In present day scenario, the convergence between the Indian and the Eurasian plate is majorly accommodated by the movement on the MFT (Lavé and Avouac, 2000; Kumar et al., 2001; Malik et al., 2003; Malik and Nakata, 2003; Lavé et al., 2005). The MFT has mappable surface expressions in piedmont zones of Uttranchal, Darjeeling, and Arunachal Pradesh Himalayas (Raiverman et al., 1990; Valdiya et al., 1992; Srivastava and John, 1999; Sarkar and Som, 2001; Mukhopadhyay and Mishra, 2004). It is a commonly accepted hypothesis that either the thrust tip is blind in most of the areas or its trace is buried by rapid accumulation of sediments (Raiverman et al., 1990; Valdiya et al., 1992; Srivastava and John, 1999; Mukul 2000; Sarkar and Som, 2001; Mukhopadhyay and Mishra, 2004). In the Darjeeling foothill Himalayas, to the east of the Tista River, the MFT has discontinuous and nonlinear surface expressions (Nakata, 1989; Guha et al., 2007; Goswami et al., 2012).

Himalayan foreland basins are suitable to understand the role of climate and tectonics in alluvial fan sedimentation, as this region has been affected by both the factors (Nakata, 1972; Ruddiman et al., 1989; Valdiya et al., 1992, 1993; Singh et al., 2001; Goodbred, 2003; Srivastava et al., 2003; Gibling et al., 2005; Sinha et al., 2005). Singh (1996a, 1996b) suggested that the mega fans formed in the initial part of humid phases preceded by a long period of arid climate, which helped in generating sediments. Srivastava et al. (2003) attempted luminescence chronology of Yamuna Ganga mega fans with multiple aliquot additive dose method and reported its formation during 26–22 ka, but did not rule out the possibility of insufficient bleaching and hence an overestimation of the ages. Jaiswal, (2006, unpublished thesis) reassessed the dates of (Srivastava et al., 2003) using Single Aliquot Regeneration (SAR) method and suggested that the mega fan sedimentation was during the humid phase (~14–16 ka) preceded by Last Glacial maxima (LGM), which supports the hypothesis of (Singh 1996a, 1996b). Under wet conditions rivers not only flush out the produced sediments but also incise the bedrock forming various levels of terraces and intense uplifting in tectonically active areas is reflected by large number of rockcut terraces (Starkel, 2003). The study of quaternary deposits may provide important clues to the sedimentation style at long mountain chains spread all over the world which are similar to Himalayas. Thus the study of alluvial fans and fluvial terraces may provide information on the river dynamics on regional scale and effect of climate-tectonic interactions on the evolution of these landscapes.

A lot of work has been done in the western Himalayas to understand the mechanism shaping the morphology of the area along with the geochronology (Singh et al., 2001; Srivastava et al., 2003; Suresh et al., 2007; Kumar et al., 2007), whereas such work is less common in Darjeeling-Sikkim Himalayas (Mukul et al., 2007; Kar et al., 2014). Current study focuses on an alluvial fan in Darjeeling-Sikkim Himalayas, drained by the rivers like Neora, Kurti and Murti in the northern part of West Bengal, India (Fig. 1). A few studies have been done earlier in the region to understand the formation of this fanbut they lack sufficient geochronological data. The alluvial fan east of the Tista River in the Mal - Murti interfluves area is named as the Matiali fan (Fig. 1). At mountain fronts in Himalayas, the Indian monsoon system plays an important role in regulating the sediment supply and its transportation, thus controls the evolution of the fan system. In addition to this the Matiali fan has been affected by several tectonic movements along the Chalsa and Matiali faults. Thus, this fan provides an opportunity to explore climate-tectonics interactions.

Fig. 1

(a) The study region lies in the northern part of West Bengal, India. (b) The figure shows the tectonic division of the Himalayas and the present study area, (c) The figure shows the major geomorphic features in the study area. The sampling sites are marked with red star. The northern most fault known as Gorubathan is marked by black line just below the Samsing surface. The two other E-W scarps/faults viz. Matiali and Chalsa are also shown in the figure. The NNE-SSW and nearly N-S scarps along which the Murti and Neora rivers flow respectively are also shown in the figure (after Goswami et al., 2013). (d) This is the Google earth image of the region showing the major rivers in the area. A red star shows the OSL sample location of S-5 coinciding with the location of radiocarbon date provided by Guha et al. (2007) (e) The figure shows the height of the different geomorphic features from the river bed near the Chalsa scarp along the line AB in Fig. 1c (after Goswami et al., 2013).


Guha et al. (2007), based on their two radiocarbon dates, one near the Gorubathan thrust and another near the Chalsa thrust attempted to explain the evolution of fan. However, the C-14 age of 33,875 ± 550 Ybp near the Gorubathan thrust (Guha et al., 2007) gives only a lower bound and shows that the fan formation started before that. This date have been found much lower than other radiocarbon age, reported to be ~46.4 ka by Ghosh et al., (2015) from a terrace in the distal part of the fan. However these dates were taken from proximal and distal part of the fan therefore cannot be compared. Kar et al., (2014) provided 3 optical ages from the terraces on the Neora and Kurti rivers in the distal part of the fan which are young ages (all below 11 ka). This shows that there is an ambiguity in terms of chronology. Goswami et al., (2013) using total station and remote sensing suggested that this fan is affected majorly by tectonics. Kar et al., (2014) discussed the evolution of the fan and linked it with the past climate neglecting any role of tectonics, in general. As the fan is located in between MBT and MFT, the development of fan is controlled by both climate and tectonics. However, which is a dominant factor, is still a controversy. The detailed chronology of various fans located along the MFT may be helpful in solving these issues. The Matiali fan is a complex geomorphic entity which has been under continuous deformation and needs detailed chronological study to frame these tectonic movements. The present work provides a detailed chronology of the area using optical dating and attempts to remove the ambiguity in geochronological data and thus, helps in understanding the geomorphic evolution of the alluvial fan with time.

. Geomorphology of the area

The study area is an alluvial fan bound by major rivers, the Mal in the west and the Murti in the east, which originate in the Higher Himalayas. The alluvial fan is dissected by many smaller order tributaries of these major rivers, e.g. the Kurti and the Juranti which originate in the piedmont zone itself (Fig. 1c). Two East-West trending scarps, the Matiali scarp in north and the Chalsa scarp in the south considered as the splay of the MBT and the MFT respectively owe their origin to the movement on blind thrust below these scarps (Nakata 1989; Guha et al., 2007) (Fig. 1). One more scarp named as the Gorubathan scarp exist to the north of Matiali scarp but could not be studied due to lack of suitable exposure, also represents MBT (Fig. 1) (Guha et al., 2007). Nakata (1989) identified three geomorphic surfaces in the Quaternary deposits, viz., Matiali, Rangmati and Samsing in order of their time of formation, Matiali being the oldest and Samsing as the youngest. He also categorized terraces in three groups: (i) Lower River Terrace or Younger River Terrace (LRT or YRT), (ii) Middle river terraces (MRT) and (iii) Higher river terrace (HRT). However, based on pedological studies, Das and Chattopadhyay (1993a, 1993b) considered Matiali and Rangamati to be a single geomorphic surface and Samsing to be the oldest surface. Similar classification was followed by later workers (Guha et al., 2007; Goswami et al., 2012, 2013; Kar et al., 2014). Goswami et al., (2012, 2013) classified regional highest surface of the fan as T4 and subsequent lower surfaces as T3, T2 and T1. The regional slope of the fan is southerly (2.1°), at the head, height is 550 m above mean sea level (MSL) which decreases to 200 m above MSL at toe (Goswami et al., 2012, 2013). Goswami et al., (2013) reported that the Matiali and the Chalsa scarps are manifestation of antiforms over two north dipping blind thrusts. The regional fan surface T4 has been dissected and incised by many rivers like Neora, Murti, Kurti and Juranti forming lower terraces. T3 occurs as raised terrace north of the Matiali scarp on both the banks of the Murti, Kurti, Juranti and only on the east bank of the Neora river. South of the Matiali scarp, raised terraces on the banks of Neora, Kurti and Murti are designated as T2 and T1. Terrace T2 has an average height of 200 m above MSL. T1 is above recent flood plains and has an average height of 180 m above MSL. Matiali scarp is highest (60 m) on the eastern bank of the Neora River and its height decreases eastward and the scarp dies out on the western bank of the Murti River after taking a north-east bend (Fig. 1). Chalsa scarp is 90 m in the Neora-Murti interfluve area, traceable from the Chel River to Mal River and in the Murti-Jaldhaka interfluves area also. The higher surfaces show signature of deformation whereas the lower small terraces are undeformed. In addition to E-W trending scarps, there are transverse faults trending NNE-SSW and nearly N-S along which the present day Murti and Neora rivers flow (Goswami et al., 2012, 2013; Ghosh et al., 2015). The radial drainage pattern, north of the Matiali scarp has been modified by later tectonism (Goswami et al., 2012).

. Luminescence chronology-methods and sampling

Optically stimulated Luminescence (OSL) dating provides the date of the event of last exposure of sediments to sunlight, thus giving burial ages (Huntley et al., 1985; Aitken, 1998). The method is well suited for aeolian deposits as these samples are supposedly well bleached in daylight before burial, however in case of fluvial deposits due to various factors e.g. sediment depth, attenuation of daylight in water column and turbidity of the water along with transport depth, some part of the sediments may or may not be well bleached, thus giving overestimated ages. With the application of Single Aliquot Regenerative (SAR) protocol of OSL dating (Murray and Wintle, 2000) one can have a better insight into the bleaching of sediments and thus can provide better estimate of ages. Luminescence dating technique using SAR has been widely used in establishing the chronology of fluvial deposits in India and in particular to Himalayas (Srivastava et al., 2003; Suresh et al., 2007; Jaiswal et al., 2008, 2009; Juyal et al., 2010; Ray and Srivastava, 2010).

The samples were collected from six locations; these sites were chosen on the basis of relative stratigraphy of different terraces (Figs. 2, 3). Among these the location of the S-5 is outside of the Matiali fan near the Taljhora scarp, close to Taljhora stream (Fig. 1d) was collected to check with independent age control provided by Guha et al., (2007). Various tests were performed on S-5 for optical dating (Fig. 4). The location of sampling sites follows and is also summarised in Table 1.

Fig. 2

The figure shows different sampling location along with the field photographs from different locations. (a) The figure shows the JRM section from which the sample JRM-1 was collected (from T4) on the banks of Juranti river near Matiali scarp. The 42° south dipping beds are marked by dashed black lines (Scale = 15 cm). (b) This figure shows a part of KRC section from which samples KRC-1 and KRC-2 (both from T2) were collected. KRC-3 was collected from T4 and was behind this section, KRC-4 was collected from T1 to the south of this photograph. The dashed black line shows the deformed beds and the orange line shows the unconformity in the T2 terrace. The height of section is about 20 m from river bed. (c) The photograph of the JRC section from which the samples JRC-1 and JRC-2 were collected, in the foreground Juranti river could be seen. Sampling location of JRC-1 is marked and the sample JRC-2 was collected from the west facing scarp of the same section (along Neora river). A scale of ~2m is marked in the lower right corner of the image. (d) NRM-1 was collected from this section, it represents the T3 terrace and parallel horizontal beds of sediments could be seen. NRM-2 was collected from the opposite bank. The black and white scale in the figure is 50 cm. (e) The photograph shows the NRC section, from which the samples NRC-1 and NRC-2 were collected (both from T4). The whole section has a very gentle slope towards north. For scale two students could be seen (~170 cm) in the right side of the photograph. The sample NRC-3 was collected from the opposite bank of the Neora river.

Fig. 3

Log showing stratigraphic position of different sections from which the samples were collected.

Fig. 4

(a) The BLSL shine down curve for the quartz from the JRM-1 sample. (b) The laboratory generated dose growth for the sample JRM-1. (c) The dose recovery test was performed on the S-5, laboratory dose of 151.4 Gy was administered and it was recovered, the red dots represent the recovered doses. It was found that all the five aliquots showed very good recovery with less than 10% of deviation from the given lab dose. (d) Histogram showing the dose distribution for the sample JRM-1.

Table 1

Details of sampling location along with coordinates and relative stratigraphic positions.

S. No.Sample NameHeight above mean sea level of River bed (m)LocationThickness of sand bed (m)Height from river bed (m)
287Western bank of Juranti river, near Matiali scarp, basement of T40.51.5
251Eastern bank of Neroa, near Matiali scarp, from T30.53
-do-Western bank of Neroa near Matiali scarp from the T40.31
358Uplifted scarp material near Taljhora0.54
178Western bank of Kurti river near Chalsa scarp, from T20.53
-do-Western bank of Kurti river near Chalsa scarp, from T40.2615
-do-Western bank of Kurti river near Chalsa scarp, from T10.283.5
188Eastern bank of Neora river near Chalsa scarp from T40.612
-do-Eastern bank of Neora river near Chalsa scarp from T40.330
166Western bank of Neora river near Chalsa scarp from T20.33
189Near confluence of Neora and Juranti rivers from T40.51.5
-do-Near confluence of Neora and Juranti from T1a0.32

Juranti River section near Matiali (JRM section)

Sample JRM-1 was collected from the western bank of Juranti River near Matiali scarp from the T4 surface. The beds here are dipping at 42° towards south (Fig. 2a). These dipping beds are of T4 fan surface material. This dipping of strata is due to the deformation on the Matiali thrust beneath (Goswami et al., 2012, 2013). The material is indurated arkosic pebbly sand (Fig. 2a). The exposed part of the section shows overall coarsening upward trend with intermittent sand layers, which is typical to alluvial fan sedimentation. The sand also varies in size from coarse to medium sized along with pebbles and cobbles.

Neora River section near Matiali (NRM section)

Two samples were collected from this site (i) NRM-1 was collected from the eastern bank of Neora river, from T3 terrace as defined by Goswami et al., (2013) and (ii) NRM-2 was collected from the western bank of Neora River near Matiali. NRM-2 was from a hard indurated polymictic conglomerate comprising sub-angular pebbles. Here the Neora river flows over this bed which represents the surface of the alluvial fan. T3 terrace shows both coarsening upward as well as fining upward trends with well imbricated pebbles and cobbles along with few boulders (Fig. 2d). Intermittent coarse sand beds also occur, indicating changes in flow energy conditions.

Kurti River section near Chalsa (KRC section)

Four samples were collected from this site; all the samples were from the western bank of the Kurti River. KRC-1 was collected from the base of T2, around 3 m from the river bed. The sand beds were dipping 17° towards south and moving further south the dip of beds changes abruptly to 25° within a distance of 10–12 meters (Fig. 2b). KRC-2 was collected from the top of T2 at a height of about 18 meters from the river bed and the sand beds were horizontal. KRC-4 was collected from the T1 surface, at a height of 3.5 m from the river bed. KRC-3 was collected from the T4 surface. The sand beds of T4 are highly oxidised and have reddish brown colour while sand beds of T2 are grey coloured. All the samples are coarse to medium grained (Fig. 2b).

Neora River section near Chalsa (NRC section)

Three samples were collected from this site. NRC-1 and NRC-2 were from the eastern bank and NRC-3 was from the western bank of the Neora river near the Chalsa scarp. NRC-1 and NRC-2 were at a height of about 12 m and 30 m from the river bed respectively. The beds were showing a gentle slope towards north and the sand was oxidised having reddish brown colour, both of these samples were from the T4 surface (Fig. 2e). NRC-3 was collected from the T2 terrace on the western bank of the Neora river.

Juranti River section near confluence of Neora and Juranti Rivers (JRC section)

JRC-1 was collected from a section near the confluence of Neora and Juranti rivers, from the sediments of T4. The section was south facing, the course of Juranti river is more or less N-S and when it reaches near Chalsa scarp it starts turning westward and finally meets with Neora river within the alluvial fan (Fig. 2c). JRC-2 was collected from a terrace designated as T1a. It is the youngest terrace and is an additional modern terrace of Neora river. A litholog of the studied area and the relative stratigraphic location of the samples collected can be seen in Fig. 3.

Samples for OSL dating were collected in galvanized iron (GI) pipes to avoid any accidental exposure to daylight, a routine method of collecting sediment sample for OSL dating. Summary of the section studied has been given in Table 1. All the samples were processed in sub-dued red light conditions under dark room environment. The sediments were treated 3–4 times with 1N HCl and 30% H2O2 to remove the carbonates and organic materials respectively. The sediments were then sieved to obtain size fraction of 90–150 gm. This was followed by heavy liquid density separation using Sodium polytungstate (density = 2.58 gm/cm3) to separate quartz and feldspar. Obtained quartz grain were etched with 40% HF for 80 minutes to remove alpha dose affected layer and feldspar contaminations present (if any after heavy liquid separation) followed by 12 N HCl treatment to dissolve the fluorides formed during etching. All samples were checked for feldspar contamination using Infrared stimulated luminescence (IRSL) at room temperature.

The quartz grains obtained from sediments in Himalayas suffers from low and changing sensitivity, in general. This limits the possibility of reducing the aliquot size, although small aliquot increases the chance of detecting poorly bleached samples (Wallinga, 2002). Even with low photon counts and bigger aliquots, this quartz can be used meaningfully in age calculation (Jaiswal et al., 2008; Ray and Srivastava, 2010). Aliquots of ~8 mm were prepared on stainless steel disks of 10 mm diameter, the grains were mounted on the disks using Silko-Spray silicone oil. Standard SAR protocol (Murray and Wintle, 2000) was applied to obtain equivalent dose for each sample. The optical luminescence measurements were made on Lexsygsmart TL-OSL reader (Freiberg Instruments, Germany) at luminescence dating laboratory in IISER Kolkata. The instrument beta dose rate was calibrated using calibration quartz supplied from the Risoe National Laboratory, Denmark Technical University, Denmark. The instrument is equipped with blue light emitting LEDs (λ = 458 nm), the detection window consists of a combination of optical filters Hoya U340 and Delta BP 365/50 EX mounted on a solid state photomultiplier tube (PMT) designed to have very low background counts. A 90Sr/90Y beta source is used for irradiation. The pre heat plateau test was performed on S-5 by heating aliquots to different temperatures for different time periods and it was found that the 260°C for 20 seconds is the most suitable temperature and time for measuring natural signal and regeneration doses. Back ground subtracted integral of 0.1 second of initial part of shine down (rapid bleaching part) curve was used for equivalent doses calculations (Fig. 4a). The dose recovery test was performed on the laboratory bleached aliquots of sample S-5. A laboratory beta dose of 151.4 Gy was administered to a set of 5 aliquots and treating them natural, recovered using SAR to evaluate equivalent dose (Fig. 4c). The results have shown excellent recovery using SAR and were recovered within 10% of the given laboratory doses. Some samples had very low signals due to which the recuperation varied from 0–10%. The disks showing best fit of data and recycling ratio within 0.9–1.1 were chosen for equivalent dose estimation. Since all the ages were well above 1 ka, average ED was preferred for age calculation (Jain et al., 2004). The growth curve was plotted for the JRM-1 to check the suitability of the sample for OSL (Fig. 4b). The histograms showing dose distribution were plotted for sample JRM-1 (Fig. 4d). ICPMS was used to measure the concentration of radioactive elements (U, Th and K). Estimation of moisture content is a very difficult task, especially in areas like North Bengal where there is a lot of variation in rainfall during different seasons therefore; estimation of 15 ± 5% in samples close to river bed and 10 ± 5% in areas at greater height from the river bed has been made.

. Results

OSL ages obtained from the sediments of alluvial fan and river terraces constrain the age of different depositional and incision phases. The results of OSL dates have been summarized in Table 2. The quartz obtained from the sample JRM-1, near Matiali scarp from the T4 surface, has the oldest depositional age of 171 ± 18 ka. Further west of this location from NRM section, sample NRM-1 was collected from the terrace T3 on the eastern bank of Neora river. It gave a depositional age of 59.7 ± 16 ka. Sample NRM-2, from the western bank of Neora river (T4) has a depositional age of 84.6 ± 8.6 ka.

Table 2

OSL age table of analyzed samples showing various parameters for age calculation.

S.No.Sample NameTerraceU (ppm)Th (ppm)K (%)Moisture (%)ED (Gy)Dose Rate (Gy/ka)Age (ka)No. of discs

In KRC section on the banks of Kurti river, KRC-1 from the base of T2 terrace has a depositional age of 47.7 ± 10.7 ka; KRC-2 from the top of this terrace has an age of 41.4 ± 10 ka, bracketing the event of sedimentation of terrace T2. At the same location, KRC-3 has depositional age of 80.8 ± 4.7 ka, representing the material from the alluvial fan (T4 surface). KRC-4, located at the southern margin of Chalsa fault, shows horizontal bedding is from the T1 terrace and has depositional age of 28.6 ± 8.5 ka. These beds are unaffected from the deformational events which modified the T2 terrace. NRC-1, NRC-2, and NRC-3 were collected from the NRC section. NRC-1 from the T4 surface has depositional age of 89.1 ± 17.3 ka, NRC-2 lies stratigraphically above NRC-1 and gives an age of 72.6 ± 19.9 ka. NRC-3 from the T2 terrace on the western bank of the river is 59.4 ± 18.1 ka old. JRC-1 collected from the JRC section near the confluence of Neora and Juranti Rivers, is the material from the T4 and has OSL age of 99.03 ± 25.3 ka. JRC-2 from the youngest terrace has deposition age of 5.9 ± 1.3 ka. The OSL age of quartz obtained from the S-5 sample from the adjacent fan located further east of the Matiali and Chalsa scarps near Taljhora is 27.8 ± 7 ka and is in close match with C-14 age of 27,210 ± 240 Ybp provided by Guha et al., (2007).

. Discussion

Guha et al., (2007) suggested that the oldest deposits in the region date 33,875 ± 550 Ybp. Goswami et al., (2013); Kar et al., (2014) took this date as a bench mark and restricted the formation of the alluvial fan to this time period. However the oldest sample dated by OSL in our study gives an age of 171 ± 18 ka, suggesting that the formation of this alluvial fan started much before than 33,875 ± 550 years.

The JRM-1 from the section along the bank of the Juranti river gives a depositional age of 171 ka suggesting that the formation of the fan started before 171 ka. This formation was accompanied with a strong fluvial activity as boulders of several meters (2–3 m on an average) diameter and a few up to 7–8 m were found in the proximal part of the fan suggesting strong monsoon. The sediments in the proximal part of the fan are disorganized; lack any stratification, suggesting mass-flow deposits. The Neora river section just below the Matiali scarp from the western bank is polymictic conglomerate with angular pebbles having age of 84.6 ± 5 ka (NRM-2). There is a large gap of ages after 171 ka. This could be due to the tectonic movement on the Matiali thrust the sediments near the Juranti section (JRM-1 dated 171 Ka) have been pushed up and are exposed for sampling, while the samples collected from the Neora river section (NRM-2 dated 84.6 ka) near Matiali (to the south of the scarp) are not affected much by this tectonic movement as shown by horizontal beds. However, sediments between these dates may be present but not accessible. Also, there is a possibility that alluvial fan deposits are not continuous deposit at one place but spread over the whole fan area, which could be one of the reasons for gap.

The samples collected from the T4 in the distal part of the fan near the Chala scarp give ages of 99 ka (JRC-1), 89 ka (NRC-1), 81 ka (KRC-3) and 72 ka (NRC-2) showing that the deposition of the fan continued during this time and the Matiali surface and the Rangmati surface (to the south of Matiali) are the facies variation of the alluvial fan and are not different surfaces as defined by Nakata (1989). The deposition on the fan continued till 72 ka or a little beyond as suggested by the date of the sample NRC-2. The parallel stratified beds of cobble, pebble and sand are seen in the distal part of the fan suggesting sheet-flood deposits. The monsoon weakened around 72 ka limiting the supply of sediments from the Himalayas to the foreland basin, also supported by the work of Rahman et al., (2009). Thus the sediment supply to the basin stopped. Since Kurti and Juranti Rivers originated in the alluvial fan, it is possible that they are synchronous with the fan formation and started incision after the fan formation and have exposed the fan sediments to its present level. OSL ages, thus suggest that alluvial fan formation started before 171 ka and continued till 72 ka at least. After that due to reduced sediment/water ratio, the river system incised the fan. After the deglaciation phase during 63 ka (Ray and Srivastava, 2010), the river got enhanced sediment supply due to increased sediment/water ratio. This helped to make fill terraces in the valley (Pazzaglia, 2013).

The undeformed thick sand beds with coarse sand below the T3 terrace near the Matiali scarp is dated 59.7 ± 16 ka (NRM-1). Benn and Owen, (1998) reported enhanced summer monsoon during the last interstadial (58–24 ka) in the Himalaya. The formation of this terrace could be attributed to the initial phase of enhanced monsoon from a transition of arid phase to humid phase. This terrace shows no signature of deformation to the north of Matiali suggesting that the formation of this terrace post dated the reactivation on the Matiali scarp. In the north of Matiali, it has an average height of 35 m above the river bed; to the south of Matiali it has an average height of 40 m above the Neora river bed (Goswami et al., 2013). Goswami et al., (2012, 2013) designated the term T3 to the terrace to the north of Matiali and T2 to the terrace to the south of Matiali, similar classification of HRT and MRT by Nakata, (1989) and Guha et al., (2007) is not well explained. In our study we took two samples, one from the T3 terrace to the north of Matiali which gave an age of 59.7 ± 16 ka (NRM-1) and another sample to the south of Matiali near Chalsa from T2 which gave an age of 59.4 ± 18.1 (NRC-3) ka. Thus we deduce that the time of formation of this terrace is around 59 ka. The sediment/water ratio was high during the beginning of monsoon as there was lot of weathered product (produced during glaciation period) available for transportation so there was an aggradational phase which led to the formation of T3 and T2 terraces along the banks of the Neora River.

The Kurti and Juranti Rivers originate within the piedmont zone near the Matiali scarp; these are small rivers and did not respond in the similar manner as the Neora and Murti rivers to climate or tectonic changes, which have their origin in the Higher Himalayas. The Kurti and Juranti rivers probably developed as a product of gully erosion and later modified in a channel as the tectonic activity and monsoonal discharge increased along the Matiali fault. From satellite image, it is clear that Juranti river has changed its course and turned west to meet the Neora River due to the activity along the Chalsa fault (Fig. 1c). A sample at height of 3 m from the Kurti river bed near Chalsa gave an age of 47.7 ± 10.7 ka (KRC-1); these beds are deformed and dip towards south. The dip is not constant, as we move north beds tend to become horizontal, pointing that this tilting is because of tectonic activity on the Chalsa fault. In the same section at a height of about 18 m sediments are dated to be 41.4 ± 10 ka (KRC-2). These beds are not deformed and lie unconformably over the lower deformed beds. Thus we can infer that the movement on the Chalsa scarp started somewhere around 48 ka and stopped by 41 ka. This movement helped in aggradation of the system and formation of terrace. Gully erosion along the T2 terrace exposed the T4 surface in contact with T2 terrace. A sample from this surface gave an age of 80.8 ± 4.7 ka (KRC-3), this is indicative of a fill terrace, i.e. T2 formed adjacent to T4.

The sample dated from T1 terrace along the Kurti River gives an age of 28.6 ± 8.5 ka (KRC-4), this terrace shows no signature of deformation suggesting that these terraces are manifestation of climatic changes. Youngest terrace along the Neora river near its confluence with Juranti river is named as T1a gives an age of 5.9 ± 1.3 ka (JRC-2) showing that the last phase of aggradation experienced by the fan was around 6 ka and after that the rivers have been continuously incising their valleys.

. Conclusion

The following conclusion could be drawn from the study:

  1. The formation of the alluvial fan started before 171 ka and continued till 72 ka covering a time span of nearly 100 ky.

  2. The luminescence dates suggest weakening of monsoon around 72 ka which reduced the sediment supply and thus the sedimentation on the fan stopped.

  3. The fan has been deformed by later tectonism along Matiali and Chalsa faults. There is no direct evidence to bracket the activation in Matiali fault however; deformed section at Juranti suggests it to be younger than 171 ka.

  4. OSL dates suggest activations in the Chalsa fault during 48–41 ka.

  5. There have been several phase of aggradations which led to formation of the various level of terraces from the regional surface T4, (i) first phase was around 59 ka which led to the formation of T3 surface, (ii) second aggradational phase was around 48 ka which led to the formation of the T2 terrace, (iii) third aggradational phase was observed around 28 ka which led to the formation of T terrace and (iv) the last aggradational phase was around 6 ka which led to the formation of T1a terrace. This suggests that climatic fluctuation during the period shaped the morphology of the alluvial fan and this was also modified by tectonic activities along the two faults.

  6. The OSL age of quartz from sample S-5 (27.8 ± 7 ka) obtained from adjacent fan located east of the Matiali fan near Taljhora scarp is in close match with C-14 age (27,210 ± 240 Ybp) provided by Guha et al., (2007). Similarly Ghosh et al., (2015) shows that the terraces on the Neora river are ~50 ka in the distal part of the fan are comparable within error to OSL dates of 57 ± 18.6 ka (Table 2). It shows that the OSL dates are in good correlation with independent studies and well suited for further detailed study in this and adjacent regions of the north Bengal parts of Himalayas.