Abstract
The Kukkal basin, Tamil Nadu, India, receives most of its rain from the southwest monsoon (SWM). A sediment core from Kukkal Lake preserves a continuous sediment record from the early-Holocene to present (9000 yr BP to present). The present lake is situated at an elevation of ~1887 m a.s.l., in a small basin that appears to have alternated between a and wetland depositional environment. Climate proxies, including sediment texture, total organic carbon (TOC), total nitrogen (TN), C/N, pollen and geochemical composition indicate a steady progression to wetter conditions, with two stepwise changes at about 8000, and between 3200 and 1800 yr BP. The change at 8000 yr BP appears to correspond to a brief (100–150 years) dry spell recorded elsewhere in India. The change at 3200–1800 yr BP consisted in a rapid intensification of the SWM, and may correlate with the initiation of the ‘Roman Warm Period’. There is no clear evidence of changes at the times of the ‘Medieval Warm Period’ (‘MWP’) and the ‘Little Ice Age’ (‘LIA’). The C/N ratio of the sediments ranges from 14.02 to 8.31, indicating that the organic matter originated from a mixture of lacustrine algae, vascular and terrestrial plants. Chemical weathering indices (Chemical Index of Alteration (CIA), Chemical Index of Weathering (CIW), and Plagioclase Index of Alteration (PIA)) are consistent with extreme silicate weathering. Pollen data show a development from savanna vegetation prior to about 8000 yr BP, followed by grassland with palms, the appearance of ferns just prior to 3200 yr BP and the establishment of the tropical humid forest between 3200 and about 1800 yr BP.
Keywords
Introduction
The Indian subcontinent is one of the largest monsoon-dominated regions of the world, and the climate is controlled primarily by the southwest monsoon (SWM). The term ‘monsoon’ refers to seasonal moisture-rich winds, and it is a unique climatic factor that differentiates south Asian countries – particularly India and parts of Asia, Australia and Africa – from the rest of the world. The Indian subcontinent is predominantly affected by the changes in these winds, which have a southerly and a westerly component. Combined together, they constitute the SWM, the sign of life and the much needed rains. The influence of the SWM is felt as far north as China’s Xinjiang region. India gets most of its rainfall by the SWM from June to September, and only the coastal tract of the Tamil Nadu state receives the Northeast rains or the return monsoon during mid-September to mid-December. Around September, as the sun retreats south, the northern land mass of the Indian subcontinent cools fast, and with this, the air pressure rises over northern India, while the surrounding atmosphere of the Indian Ocean holds the heat. This causes the cold wind to move down from the Himalayas and Indo-Gangetic Plain towards the vast spans of the Indian Ocean south of the Deccan peninsula. This is known as the northeast monsoon (NEM) or Retreating Monsoon. However, the amount and intensity of both the monsoonal rains have fluctuated since the late Quaternary period (Clift and Plumb, 2008; Gupta et al., 2006; Singhvi et al., 2010, 2012). During the last two decades, scientific interest in the reconstruction of the Asian monsoon variations, particularly the SWM, since the late Quaternary to Holocene period has increased significantly. In particular, in the Indian subcontinent, many types of natural archives such as lake and marine sediment cores, palynology, tree rings and corals (Mayewski et al., 2004) have been used to restructure late Quaternary–Holocene palaeoclimatic changes. Such records include fluvial archives of the Gujarat region (Chamyal et al., 2002), deep marine sediment cores studied from the Arabian Sea and the Bay of Bengal (Achyuthan et al., 2014; Chauhan et al., 2010; Gupta et al., 2003, 2005; Nagasundaram, 2014; Nagasundaram et al., 2014), lacustral and aeolian sediments from the Thar desert (Achyuthan et al., 2007a, 2007b; Enzel et al., 1999), lacustrine sediments from Kerala (Farooqui et al., 2010; Farooqui and Sekhar, 2011; Kumaran et al., 2008; Veena et al., 2014), peat deposits from the foreland basin of the Himalayas (Rawat et al., 2015; Sharma et al., 2006), glacial moraines and proglacial deposits of the Ladakh region (Sant et al., 2011), speleothems (Yadava and Ramesh, 2005), stalagmites (Sinha et al., 2007) and tree rings (Esper et al., 2002a, 2002b, 2003; Shah and Bhattacharya, 2009; Yadav and Singh, 2002).
Lakes and wetlands in this region function as repositories of sediment deposited in a closed system and act as major sinks for organic carbon. Therefore, sediment cores from such sites in the area affected by the SWM can provide a continuous, high resolution of palaeoenvironmental record of the evolution of the SWM (Verschuren, 2003; Verschuren et al., 2000). It is important to reconstruct past climate and to study palaeoclimate records with high resolution especially since the early-Holocene period because the past performance of the monsoon may provide vital input for projecting future climate change and predictions. Recent high-resolution studies from the northern Andes show that vegetation has responded to climate change within decades to centuries (Gonzalez-Carranza et al., 2012). The growth and impacts of the human species worldwide had profound implications for human development, advance of major river valley civilizations took place because of rain-fed agriculture and were established during this time interval (Allchin, 1997; Enzel et al., 1998; Gupta, 2004; Ponton et al., 2012; Possehl, 1996) and that from the past records the early-Holocene period was wetter than today (Allchin, 1997; Enzel et al., 1998; Gupta, 2004; Kaufman et al., 2004; Possehl, 1996). Hence, it is important to reconstruct past climate changes with higher resolution and greater degree of accuracy.
The Holocene Megathermal, also known as Hypsithermal, Altithermal or Holocene Climatic Optimum (HCO) or Holocene Optimum or Holocene Thermal Maximum, was a warmer phase during roughly the interval 9000–5000 yr BP in areas that received East Asian monsoon (EAM) (An et al., 2000; Sun and Li, 2012; Wang et al., 2005). In the Asian region, palaeoclimate records based on lake sediments indicate that both southwest summer and northeast winter rains were stronger in the early-Holocene period (10,000–7000 yr BP) than today, resulting in a wetter and cooler climate over the Asian continent. Warming started around 8000 yr BP, and the Holocene Megathermal Maximum (HMM) that occurred between 7200 and 6000 yr BP observed particularly change in vegetation zones, increase in temperatures, precipitation and sea-level fluctuations. Annual mean temperature during this period was estimated at 1°C higher than today for South China; 2°C higher for the Changjiang River Valley; about 3°C higher for North China, Northeast China and Northwest China; and up to 4–5°C higher for southern Qinghai-Xizang Plateau (Shi et al., 1993). After 3000 yr BP, climate in China area returned to a cold and dry state (Sun and Li, 2012). Peat deposits from Sanjiang Plain, China, revealed a dry and warm climate from 1900 to 1000 yr BP, a warm and humid phase from 1000 to 500 yr BP, and a cool and dry climate since 500 yr BP (Xia and Wang, 2001). A period of warming in Oman coeval with that in China was observed from the index fossil Globigerina bulloides, which indicates wind-induced upwelling of ocean water, during the interval 1200–600 yr BP (Gupta et al., 2003). These conclusions are based on a wide variety of proxies, including pollen and diatoms (Prell and Kutzbach, 1987; Van Campo et al., 1996), oxygen isotopes from marine sediments (Sirocko et al., 1993, 1996) and Arabian Sea upwelling intensities (Overpeck et al., 1996; Prell and Van Campo, 1986). Staubwasser et al. (2002, 2003) provided some insight into the climate variability in the SWM region since the mid-Holocene.
A history of flow volume of the Indus River over the past 5000 years has been reconstructed by Von Rad et al. (1999a, 1999b). The Indus has a large catchment in a region that is now partly semi-arid. Being fed by melting Himalayan ice as well as precipitation in the semi-arid hinterland, its discharge is governed by a complex interplay of phenomena. Information obtained to date shows a stronger SWM since about 1500 yr BP. Temperature indicators are less clear, except around the ‘Medieval Warm Period’ (‘MWP’) (about 1000–600 yr BP) and the ‘Little Ice Age’ (‘LIA’) (Cronin et al., 2003; Crowley and Lowery, 2000). Chauhan and Vogelsang (2006) reported that the ‘LIA’ had a large temporal and spatial variability because of the variations in the intensity of SWM and had its impact across the globe. However, in the northern hemisphere, the ‘LIA’ was a mild cooling event (Anderson et al., 2002; Cronin et al., 2003; Crowley, 2000; Crowley and Lowery, 2000; Singhvi and Kale, 2008). Summarizing these studies, it is evident that a significant variation of the SWM, encompassing several cycles of fluctuations, occurred during the Holocene.
To date, the marine record from the Arabian Sea has proved to be the most consistent and continuous record of the SWM fluctuations (Kudrass et al., 2001; Lückge et al., 2001; Sirocko et al., 1993, 1996, 2000; Thamban et al., 2002; Tiwari et al., 2005, 2006), whereas limited studies have been carried out in the Bay of Bengal (Achyuthan et al., 2014; Chauhan et al., 2004; Chauhan and Vogelsang, 2006; Duplessy, 1982). However, lakes of peninsular India preserve a great wealth of palaeoenvironmental and palaeoclimate data spanning the last 6000 years as described by Limaye et al. (2010) and Veena et al. (2014). The lake and wetland archives of Tamil Nadu have received so far little attention. Kukkal Lake, in the Palani Hills, Dindigul District, Tamil Nadu State, lies in a basin affected by both the SWM and NEM rains. Hence, the main objectives of this study are to reconstruct the Holocene palaeomonsoon and palaeoenvironmental shifts using the Kukkal Lake sediments.
Study area
Kukkal Lake (10°16′N: 77°22′E) is a small, irregular, fresh water body covering an area of 18 ha, and situated at an altitude of 1887 m a.s.l. (Figure 1) about 33 km from Kodaikanal, Dindigul district, Tamil Nadu, India. The present lake is impounded by a low dam at the northern end, and floods an area that appears to have been a wetland similar to one about 300 m south of the lake. Both features may have been lakes at times in the past. Presently, the lake has a catchment area of 2–3 km2 that is forested and relatively undisturbed by human activities; some abandoned paddy terraces on the east side of the lake are the only evidence of agriculture in the catchment. Sediment in the basin beneath the lake is therefore suitable for carrying out palaeoenvironment research. Bedrock underlying the catchment is Precambrian charnockite associated with hornblende-biotite gneiss, granite and quartzite, and is covered by a thin veneer of red vertisol and lateritic soil (Bagyaraj and Gurugnanam, 2011). The geology and geomorphic observations in and around the lakes reveal that the Kukkal Lake basin is structurally controlled by faults and lineaments, as are the main tributary watercourses within the basin. The catchment receives rainfall from both SWM and NEM, the NEM being dominant (Meher-Homji and Gupta, 1999). The mean annual rainfall for Kodaikanal over a 4-year period (2001–2004) was 1690 mm (Bagyaraj and Gurugnanam, 2011). The dry season is from December through March, and there is also a short dry period in July. The climate is tropical humid for most of the year, with summer temperatures reaching 19.8°C (max.) and 11.3°C (min.) and winter temperatures between 17.3°C (max.) and 8.3°C (min.). Meher-Homji and Gupta (1999) state that Kodaikanal experiences a bixeric climate regime receiving high rainfall from October to December (i.e. NEM). The lake is surrounded by the ‘Shola’ forests (the word ‘Shola’ derived from the Tamil language word ‘cÕlai’ meaning grove or thick vegetation cover (Meher-Homji, 1967). The Sholas are patches of stunted tropical montane forest found in the valleys in the midst of rolling grassland in the higher montane regions of the Nilgiri Hills, South India. The soil texture varies from clay to silty clay; soils are acidic in nature and contain a high percentage of iron oxide and alumina. In the ‘Shola’ forests, accumulation of humus in the top layers of the soil renders it a black colour (Bera and Farooqui, 2000; Meher-Homji, 1967). The dominant trees in the Shola forest are represented by Michelia nilagirica, Bischofia javanica (bishop wood), Calophyllum tomentosum, Cedrela toona (Indian mahogany), Eugenia (myrtle) spp., Ficus glomerata and Mallotus spp.; and small trees such as Pygeum gardneri, Schefflera racemosa, Linociera ramiflora, Syzygium spp., Rhododendron nilgiricum, Mahonia nepalensis, Elaeocarpus recurvatus, Ilex denticulata, Michelia nilagirica, Actinodaphne bourdillonii and Litsea wightiana. These forests are interspersed with grasslands, characterized by frost and fire-resistant grass species like the Chrysopogon zeylanicus, Cymbopogon flexuosus, Arundinella ciliata, Arundinella mesophylla, Arundinella tuberculata, Themeda tremula and Sehima nervosum.

Map showing the location of the Kukkal Lake and the sediment coring site.
Materials and methods
Field methods
The field map was prepared using Survey of India toposheets, a Google Earth map and Satellite imagery data products IRS LISS III for delineating faults, lineaments and channels draining into the lake. Faults and the rock types exposed around the lake surroundings were also noted. A piston core of 1 m in length was used to retrieve the sediment core from lake floor (Figure 1). The cored sampled location was marked using a Garmin III Global Positioning System (GPS) (Figure 1). The sediment core was then subsampled at 2-cm intervals, collected in clean zip-locked polythene bag, marked, described (Table 1) and stored for further analyses.
Lithology of the Kukkal Lake sediment core.
Radiocarbon dating
Four organic carbon rich samples from the core were subjected to acid-base-acid pre-treatment, and dated by liquid scintillation counting on benzene at the Department of Geosciences, University of Arizona, Tucson, USA. The radiocarbon ages were then calibrated using Calib 5.2.0 (Stuiver et al., 2010; Stuiver and Reimer, 1993) (Figure 2, Table 2).

Sediment texture, radiocarbon ages, sedimentation rate (mm/a), CaO/MgO and CWI variation with depth.
Radiocarbon dates of the lake sediment core collected from the Kukkal Lake.
Chemical and physical analyses
For the sediment texture analyses, about of 5 g of the dried sediment was washed following Carver (1971) through a 230 American Standard Testing Machine (ASTM) sieve mesh opening of 0.063 mm until the water passing through was clear (Tables 1 and 2), and the washed water was collected in a large porcelain basin. This sieving separated the sand particles from the fine fractions of silt and clay particles. The sediments collected in the 230 sieve mesh were calculated as the sand content. Furthermore, the silt and clay percentages were determined using the measuring jar and settling velocity procedure put forward by Carver (1971). Procedures of Loring and Rantala (1992) and Piper (1947) titration methods were adopted for the determination of calcium carbonate content. The organic matter (OM) in the sediment sample was determined following the methods of Gaudette et al. (1974). The total organic carbon (TOC) and total nitrogen (TN) were determined using about 1 g of sediment sample, decarbonated with a 1M solution of hydrochloric acid, washed three times with deionized water, freeze-dried and then analysed in a Thermo Scientific FLASH 2000 CHNS-O Organic Elemental Analyser (OEA) at the Department of Geology, Anna University, Chennai (Figure 3, Table 3). Major oxide and trace element analyses were carried out at the Central Instrumentation Facility (CIF) at Pondicherry University, Pondicherry, using a wavelength dispersive x-ray fluorescence (WD-XRF) spectrometer with samples in powder form. The samples were analysed for SiO2, Al2O3, Na2O, CaO, Fe2O3, K2O, TiO2, P2O5, MnO, MgO, Sr, Cr, Cu, Ni, Co, Pb, Zn, V, Ba and Zr (Tables 4 and 5). Results were normalized to Al (Figure 4).

Down core variations of texture, OM, TOC, TN percentage and C/N ratio showing increase in sand from 7640 to 3215 yr BP because of sediment flux from the catchment area.
Sediment texture, OM, TOC, TN content, C/N ratio and age of Kukkal sediment core.
OM: organic matter; TOC: total organic carbon; TN: total nitrogen.
Ages: measured interpolated. Compositions: maximum in bold, minimum in bold italics.
Major oxides (%) and weathering index of the sediments from the Kukkal sediment core.
CIA: Chemical Index of Alteration; CIW: Chemical Index of Weathering; PIA: Plagioclase Index of Alteration; CWI: Chemical Weathering Intensity.
Maximum in bold and minimum in italic.
Trace element content (ppm) of the sediments from the Kukkal sediment core.
Maximum in bold and minimum in bold italic.

Down core variations of major elements normalized by Al and CWI values showing distinct variation during short extreme events such as the ‘RWP’ and ‘MWP’.
Calculation of chemical indices
The degree of rock alteration and weathering can be estimated quantitatively using several geochemical indices (e.g. Duzgoren-Aydin et al., 2002; Maslov et al., 2003; Price and Velbel, 2003; Yudovich and Ketris, 2011). However, for the charnockite and gneisses exposed around the study area, weathering indices to date have not been calculated. Most of these indices are expressed in terms of molecular or weight percentages of various oxides or groups of oxides. During chemical weathering, more labile minerals, such as feldspar and plagioclase, are depleted in Ca2+, K+ and Na+ and transformed into minerals that are more stable under surface conditions. As the intensity of weathering increases, rocks and sediments become enriched in Al, Ti, Fe and Mn (Minyuk et al., 2014).
Chemical Index of Alteration (CIA): The CIA was defined as follows by Nesbitt and Young (1982) using molecular proportions:
with CaO being the amount of CaO incorporated in the silicate fraction of the rock. The CIA measures the proportion of Al2O3 versus more labile oxides, and reflects the relative amount of feldspars and clay minerals in a sample. A CIA value of 50 corresponds to unaltered albite, anorthite and potassic feldspar. Typical values of the CIA are 30–40 for basalt, 45–55 for granites and granodiorites, 75–85 for illite, 75 for muscovite and close to 100 for kaolinite and chlorite (Nesbitt and Young, 1982).
Chemical Index of Weathering (CIW). This index was proposed by Harnois (1988) and is calculated using the equation:
with Al2O3 treated as an immobile component and CaO and Na2O as mobile elements. K, in contrast, is not included in this index because it may be leached and/or accumulated in the residual weathering products. Fedo et al. (1995) suggested that Al in this index be used without any correction for its occurrence in K-feldspar, and therefore, generally, K-feldspar-rich rocks yield very high CIW values. Potassic granite, for example, has a CIW value of 80, and fresh K-feldspar has an index of 100, values that are similar to those of the residual products of chemical weathering, for example, 80 for smectite and 100 for kaolinite, illite and gibbsite.
Plagioclase Index of Alteration (PIA). This index (Fedo et al., 1995) estimates the weathering intensity of plagioclase feldspar and is calculated as follows:
Chemical Weathering Intensity (CWI). Oxides such as CaO, MgO and Na2O are considered more soluble and mobile, while Al2O3, SiO2 and TiO2 are considered to be more insoluble and resistant (Engstrom and Wright, 1984; Mackereth, 1966). Sun et al. (2010) proposed the CWI as an alternative to the CIA in order to understand the intensity of weathering and apparently wet or dry conditions. The CWI was calculated as follows, using molecular proportions:
Palynology
Forty eight sediment samples (taken at 2-cm intervals) were analysed. Ten grams of each sample was treated with warm 10% potassium hydroxide for removing the OM coagulation from the sediments and later sieved through 150 ASTM mesh (105 µm) in order to remove roots/rootlets or any other macro organic/inorganic parts in the sediment. The sieved filtrate fraction was allowed to settle overnight and the supernatant was drained. The residue was treated with 40% hydrofluoric acid (HF) for 2–3 days depending on the sand content in sediments. The HF was decanted after centrifuging the samples at 2000 r/min for 10 min.
The sample was then acetolysed, following Fægri and Iversen (1989), using glacial acetic acid, and later after centrifuging and decanting, the samples were treated with warm acetolysis mixture (anhydrous acetic acid and concentrated H2SO4 in 9:1 ratio). Decanting of all the acid was done after centrifuging. The residue was washed with distilled water and decanting was done by centrifuging every time at 2000 r/min for 10 min. After decantation, the samples were passed through 650 mesh size (~10 µm). The residue was then made up to 5 mL in 50% glycerine and water. Whole of the 5-mL sample was studied, and the total count of each ecological group presented in Figures 6 and 7 includes pollen sum. Arboreals and non-arboreals along with aquatic pollen, monocots, pteridophytic spores and the ubiquitous taxa were included in pollen sum because these perhaps grew in the near vicinity of the studied site and their relative percentages provide insight into its equilibrium with the prevailing climatic conditions in the past (Figures 6 and 7). The 650 mesh (10 µm) fraction was collected and mounted on a glass slide in glycerine jelly for palynological study under a light microscope (Olympus BX-52). Pollen atlases of Guinet (1962), Vasanthy (1976), Nayar (1990), Tissot et al. (1994), Farooqui et al. (2010) and reference slides from the Herbarium collection of the Birbal Sahni Institute of Palaeobotany were consulted for identification of the palynomorphs. The phytolith cells were observed in the palynological slides and only qualitative report is presented in the present work. The quantitative analysis through standard methodology is in progress and will be published later.
Results
Lithology
Based on the colour and sediment texture, three lithological units were recognized in the core. Lithological unit I (100–47 cm) from the lower end of the sediment core is represented by black coloured (10 Y 2/1) (Munsell colour scale codes), fine sandy silt, clayey silt and silty sand, with peaty layers, overlain by lithological unit II (47–28 cm) which is brownish black (10 Y 2/3) in colour, moderately sorted, sandy silt and silty sand, with peat laminae within sub-angular to angular grains that are poorly sorted. Lithological unit III (28–0 cm) is brownish black (10 Y 3/1) in colour and consists of sandy silt alternating with silty sand laminae, with minute specks of peat with well rounded sand grains. This lithological unit contains abundant roots in varying stages of decay and root pores (Table 1, Figures 2 and 3).
Chronology
Radiocarbon dates (Figure 2, Table 2) range from 9105 yr BP (99 cm depth, ~9255 BC) to near modern (99.7 pMC at 14 cm depth, ~AD 1650–1955 or 300 yr BP). The plot of depth versus age indicates that the sedimentation rate from the early-Holocene to recent fluctuated, with greater sedimentation rate during 9100–7600 and less after 7600. Sedimentation rate from 7640 to 300 yr BP was fairly constant. The small variations are as follows: 0.19 mm/a from 9105 to 7640 yr BP; 0.06 mm/a from 7640 to 3215 yr BP (71–43 cm depth), 0.08 mm/a from 3215 to about 300 yr BP and about 0.4 mm/a since the last 300 years. Above 14 cm, the deposition rate is higher, with a minimum value, calculated by attributing an age of 300 yr BP at 14 cm, of 0.4 mm/a.
Stable carbon isotopes
Values of δ13C of organic carbon in the sediment show a general upward decrease from near −16‰ at 98–100 cm (9105 yr BP) to near −23‰ at present (Table 2). The transition can best be explained as the result of a decreasing proportion of OM derived from C4 plants over time. Both C3 and C4 sedges grow at present in a peat bog in a similar tropical montane environment in the Nilgiri Hills, Tamil Nadu State, India (Rajagopalan et al., 1999).
Sediment texture
Textural analysis of clastic fraction of the Kukkal sediment core is shown in Figure 3. The sand percentage ranges from 12.5% (depth 70 cm) to 92% (depth 50 cm). The silt percentage ranges from 8% (depth 50 cm) to 75% (depth 70 cm), while the clay ranges from 0.5% (depth 50 cm) to 17% (depth 85 cm). Three main zones are apparent: below about 65 cm, with silt > sand and clay > 10%; from about 60 to 40 cm, sand predominant, clay < 5%; above 40 cm, sand and silt about equal, and clay < 5% (Figure 3). Most of the sediment core contains >25% OM (peaty in nature) consistent with a lacustrine or wetland environment. Recurring peaty sediment suggests that a wetland environment may have prevailed over much or whole of the interval represented by the sediment core, prior to construction of the dam.
Sediment chemistry
TOC content increases steadily up the section from about 3.5% at 85–100 cm to 15.8% at 1 cm unit I to unit III (Figure 3, Table 3). Similarly, TN content increases gradually from unit I 0.4% (depth 85–100 cm) to unit III 1.0–1.3% (depth 0–30 cm), and the C/N ratio increases, mimicking the TN content (unit I to unit III) gradually from about 8% (depth 85–100 cm) to 14% (depth 1 cm). A C/N ratio 13–14 suggests a sub-equal mixture of algal and vascular plant contributions (Meyers and Lallier-Vergès, 1999). Intervals of low C/N ratio (<9) occur at 12–14, 22–24, 30–34, 40–44 and 76–100 cm core depth (Figure 3). The sample at 2–4 cm probably represents the present artificial lacustrine environment and has a C/N ratio of 14, indicating that vascular and algal plant material is deposited in shallow standing water in this environment. At the depth of 40–35 cm, a concurrent decrease in OM, TOC and TN immediately follows a sharp decrease in sand deposition (Figure 3).
The vertical profile of the major elements normalized by Al versus depth (Figure 4) shows a distinct variation with depth, and the major oxide content (Table 4) shows that the lake sediments have SiO2 varying from 31% (depth 70 cm) to 47% (depth 35 cm), with Al2O3 from 13% (depth 20 cm) to 22% (depth 80 cm) and Fe2O3 from 5% (depth 50 cm) to 9% (depth 75 cm). TiO2 varies from 0.9% (depth 10 cm) to 1.3% (depth 40 cm). K2O, CaO, MgO, MnO and Na2O are present in small (<2%) amounts. The relative concentrations of the major oxides in the sediment core can be summarized as SiO2 > Al2O3 > Fe2O3 > TiO2 > MgO > CaO > K2O > P2O5 > MnO > Na2O. Most of the oxide trends show two prominent features: decreases in SiO2, Al2O3, Fe2O3 and TiO2 at 30–35 cm, accompanied by increases in P2O5, MnO and CaO, and at 70–80 cm, decreases in SiO2, Al2O3, P2O5 and K2O (Tables 3 and 4). There is an overall significant decrease of Al2O3 (about 10%) since 9105 yr BP (Figure 4).
In the sediment sequence of the Kukkal Lake core, the values of the CWI increase from 4–6 (characteristic of depths below 40 cm) to 8–10 (characteristic of depths above 30 cm) (Table 4). Generally, higher values of CWI indicate stronger weathering and presumably wetter conditions and a strong monsoon (Sun et al., 2010). Concurrent decreases occur in the values of the CIA, PIA and CIW between 40 and 30 cm (Table 4). The shifts in all four indices reflect an increase in CaO concurrent with a decrease in Al2O3 between 40 and 30 cm. Scatter plots of the CIA versus CIW and PIA values reflect extreme silicate weathering of the bedrock exposed around the lake catchment and the deposition of extremely altered sediments into the lake (Figure 5).

Scatter plots of CIA versus CIW and CIA versus PIA values of the core.
Palynology
Three palynological Phases I–III were identified on the basis of relative percentages of ecological groups identified during the study. A cluster analysis of the results was carried out and the display of palynological spectrum was prepared following Grimm (1987).
Phase I (90–98 cm)
This phase shows low arboreal pollen (7.1%) and high non-arboreal pollen (16.8%). The average relative percentage of Poaceae and Cyperaceae accounts to 35.1%. Highest percentage (~40%) of pteridophytic spores (Monolete and Triletes) was recorded. The percentage of Arecaceae remained low. The arboreals constitute Acacia (1.7%), Bignoniaceae (3.6%), Elaeocarpus (2.4%), Eugenia (1.2%), Eurya (1.2%), Shorea/Hopea (2.4%); and non-arboreals comprise Malvaceae (2.1%), Apiaceae (2.2%), Artemisia (3.3%), Asteraceae (9.3%), Chenopodiaceae (0.9%). Cyperaceae (7.0%), Poaceae (28.1%), Arecaceae (8.1%), Monolete spore (16.5%) and Trilete spore (16.2%) were recorded.
Phase II (30–90 cm)
Poor count was observed in this phase with no arboreal pollen. The average relative percentages of these are Asteraceae (15.8%), Cyperaceae (4.9%), Poaceae (49.2%), Arecaceae (26.7%), Monolete spore (7.2%) and Trilete spore (5.4%).
Phase III (30–2 cm)
This phase is comparatively good in pollen/spore preservation. Total 29.1% of arboreal and 18.3% of non-arboreal pollen were recorded. Aquatic pollen account to about 11.6%. The ubiquitous taxa such as Poaceae and Cyperaceae account to 20.9%, along with pteridophytic spore showing 20.1% of the total count. Arecaceae accounts to about 4.9% of the total count. The arboreals constitute Acacia (1.7%), Bignoniaceae (3.6%), Elaeocarpus (2.3%), Eugenia (3.4%), Euphorbiaceae (4.3%), Eurya (0.5%), Shorea/Hopea (4.2%), Myrtaceae (5.8%), Syzygium (4.3%) and Terminalia (3.3%). The non-arboreals constitute Lamiaceae (2.3%), Malvaceae (1.2%), Apiaceae (2.9%), Artemisia (0.8%), Asteraceae (6.0%) and Chenopodiaceae (0.9%). The aquatic pollen present were Myriophyllum (0.7%) and Lemna (10.9%). Cyperaceae (11.9%) and Poaceae (9.0%) along with Arecaceae (4.9%), Monolete spore (8.1%) and Trilete spore (7.1%) were recorded (Figures 6 and 7).

Pollen spectra in sediments of Kukkal Lake.

A comparison of various palaeoclimate proxies: (a) Globigerina bulloides percent from the ODP Hole 723A, Arabian Sea 39; (b) haematite % in core MC52 in the North Atlantic; (c) Oxygen isotope record from Dongee Cave, southern China; (d) Hongyuan peat deposit; and (e) Qunf cave δ18O and (f to h) Sand, C/N and CWI of Kukkal Lake (present study).
Discussion
Palynological environment
Increased non-arboreal pollen along with highest percentage of Monolete/Trilete fern spores with low arboreal pollen in Phase I indicate moist, fluctuating and non-static climatic conditions with sparse vegetation, suggesting more open land with herbaceous taxa. The Phase II exhibits a complete absence of arboreal pollen followed by scarce non-arboreals with occurrence of high pollen of Poaceae and Arecaceae, suggesting increased open land and drier climatic conditions. Fern spores (Monoletes and Triletes) too decrease significantly, suggesting low moisture conditions. The beginning of Phase II shows total absence of Cyperaceae and pteridophytic spores, pointing towards cooler and dry 8200 yr BP global event which is also evident from different parts of Indian subcontinent (Dixit et al., 2014). The Phase III shows good percentage of arboreals and comparatively quite low non-arboreals. Poaceae increased along with Arecaceae and pteridophytic spores, suggesting the shift in climatic conditions from comparatively drier to moist conditions (strengthened precipitation) favouring the vegetation in the catchment vicinity. Overall palynological results in Phase III indicate the lacustrine ecosystem as the silt percentage increased along with high OM, TN and carbon/nitrogen ratio. Major elements too show high percentage in this phase. Phases I and II show fluctuating values with respect to all the studied proxies, suggesting fluctuating climatic conditions, but Phase III indicates static environmental conditions and the existence of lake in the studied area perhaps since ~3000 yr BP (Figure 6).
The increase in fern spores and their ecology, however, remains poorly understood because of the sustained lack of information on the gametophyte generation (Limm and Dawson, 2010). It is unsurprising that we know little about the factors influencing the decline/abundance of fern species (Bucharová et al., 2010; Kelsall et al., 2004; Testo and Watkins, 2013). By studying the development, reproductive biology and stress ecophysiology of fern gametophyte generation, it is noted that at both 20°C and 25°C temperatures, the sporophyte formation was affected (Limm and Dawson, 2010). It generally adapts itself from well to cooler (temperate) moist climate and also has the capacity to transfer moisture from the fog within the forest ecosystem even when the precipitation is comparatively low. In this study, it is observed that high percentage of fern spores (both Trilete and Monolete spores) correspond to the time period when relatively high arboreal/non-arboreals are recorded. Thus, it is inferred that high percentage of fern spores recorded here indicate an increase in moisture level with relatively cooler climate enhancing the stabilization of forest cover around the lake. An integration of all the proxies studied is presented below.
Sedimentary environment
The chronology of the core, and inferred measurements of sedimentation rate, indicates a sedimentary environment that has been continuous and close to stable in this small basin for the past 9000 years. However, high proportion of sand deposited between 9100 and 7600 (60–45 cm) interval indicates a sediment flux from the nearby catchment area and/or faster accumulation rates and less after 7600 yr BP, and a decrease in sand flux after 1700 yr BP probably because of drier conditions. Peat-like deposits of OM are consistent with a wetland environment like those present elsewhere in montane valleys of the tropical Western Ghats. Variations in the C/N ratio of OM may indicate changes in the abundance of algae in the basin and, by implication, changes in the amount of water depth in the lake (Figure 3). The basin may have been lacustrine for relatively brief intervals corresponding to 12–14, 22–24, 30–34, and 40–44 cm, in addition to a longer interval corresponding to 76–100 cm in the core (Figure 4). The present-day lacustrine environment is because of damming. Superimposed on this sedimentary system are changes revealed by our climate proxy measurements – both unidirectional steady changes and abrupt events.
Unidirectional changes
The chemistry of the OM in the core has evolved steadily over the past 9000 years. TOC, TN and C/N all show an upward increase in the core (Figure 3). The values of δ13C also suggest an increase. A progression from mainly C4 to mainly C3 species in the basin is indicated, and this is in support with an increasingly wet climate. Increasing rainfall may also be responsible for the increasing amount of OM over time, preservation being more likely when the basin sediments remain wet for long periods. The CaO/MgO ratio in the silicate fraction also increases towards the top of the core section, and Al2O3 decreases.
Stepwise changes
Several proxies show stepwise changes in the range of 30–40 cm. Major oxide measurements change in such a way that all four chemical weathering indices calculated for this study shift to new values corresponding to more intense weathering over that depth interval. Fe2O3, MnO and P2O5 also show shifts. The changes in chemical indices are concurrent with increases in Ca/Al and Si/Al (Figure 4), but there is no change in the relative abundance of clay (Figure 3). Proportions of Asteraceae, Poaceae and Palmae or Arecaceae decrease markedly in the record over this depth interval (Figure 6). These changes immediately follow a change in the character of clastic sedimentation from sand-dominant to sand and silt in approximately equal proportions. They are succeeded rapidly by a further change in the pollen record, namely, enhanced preservation of arboreal pollen. Ferns increase in abundance earlier, at depths of 46–40 cm in the core. This change is at ca. 3000 yr BP (at 40 cm) to 1800 yr BP (at 30 cm). A stepwise change similar to that reflected in enhanced preservation of pollen and the appearance of diatoms just above 30 cm in the Kukkal core and occurred at 1500 yr BP in Pookode Lake, Kerala, at an altitude of 770 m.a.s.l (Veena et al., 2014).
As noted before, the shifts in the chemical indices appear to relate to increased deposition of CaO as silicate minerals, and concurrent increase in Si/Al. A change in clay mineralogy, with deposition of Ca-smectite above 30 cm, is a possible way of accounting for the change. The increase in CWI indicates a greater intensity of weathering, but deposition of sediment bearing this signature is likely to have lagged behind the climate stimulus responsible for the change. Together, most of the observations suggest a period of increasing rainfall. The changes may even have begun as early as 6000 yr BP, when sand began to dominate the clastic sediment, and have led to more abundant fern vegetation by 3500 yr BP, and this forest vegetation largely replacing grasslands with palms by 1700 yr BP. This change might account for the decrease in sand deposition at about 1700 yr BP, if the change to forest vegetation led to less energetic runoff.
The increase in sand deposition near 6000 yr BP could be a separate stepwise event, independent of the events beginning around 3500 yr BP. A further stepwise event is indicated at 90 cm in the core (about 8700 yr BP), prior to which some arboreal pollen and fern spores were preserved. Above this horizon, the relative abundances of sand, silt and clay change abruptly towards the top in the sediment core. Major oxide trends show inflections at about 80 cm, possibly resulting from the same event, with a lag reflecting the time required for weathering chemistry to respond.
Human effects
The higher sediment deposition rate corresponds to events of the last 350 years (17–18 cm in depth to the core top) and is the only evidence consistent with human-caused changes in the lake catchment. As noted above, abandoned paddy terraces are present along the east side of the catchment. Agriculture may have become established in this montane site since the early 17th century. Cutting of trees in the forested part of the catchment may also have contributed to the increased sedimentation rate.
Regional and global correlations
The Kukkal lacustrine/wetland basin has responded to hydrological changes, principally changes in the amount of precipitation from the SWM, over the 9000-year record presented here. Several other proxy records of similar length are available for comparison in the region of the SWM and the EAM. Oxygen isotope records from radiometrically dated speleothems (Cheng et al., 2009; Fleitmann et al., 2003, 2004, 2007; Liu et al., 2013; Morrill et al., 2013a, 2013b) show a general unidirectional decrease in δ18O since about 8000 yr BP (Figure 7). The trend is interpreted by the authors, on the basis of the isotope amount effect, as evidence for a gradual increase in intensity of the SWM and EAM, and is consistent with the unidirectional changes at Kukkal, also interpreted as evidence of an increase in monsoon intensity. Similar unidirectional changes occur in another proxy, for example, the abundance of Globigerina bulloides in seafloor sediment in the Arabian Gulf near Oman (Gupta et al., 2003). The general increase in SWM rainfall is consistent with the results of Yadava and Ramesh (2005), who interpreted a speleothem record as indicating the highest rainfall around AD 1666 (‘LIA’).
The speleothem records also show a trough in the δ18O between 8500 and 8000 yr BP, interpreted in part as a decrease in monsoon intensity beginning at 8200 yr BP and lasting 100–150 years. An event close to this time is present in the Kukkal sediments at an interpolated age of 8700 yr BP, where it is marked by the replacement of savanna grassland and inflections in the curves of major oxide abundances. Given the uncertainties of the age interpolation, the proxy records may all be recording the same event.
At Kukkal, the more pronounced stepwise event between 3000 and 1800 yr BP has no clear counterpart in the speleothem or Globigerina bulloides records (Figure 7). Possible counterpart events in climate records from the Indian subcontinent include the warmer and humid event at 1500 yr BP at Pookode Lake, Kerala (Veena et al., 2014) and speleothem δ18O record from Orissa, India (Yadava and Ramesh, 2005). In the case of the latter, the authors interpret the record as indicating increasing rainfall up until 1500 yr BP, followed by a period of decreasing rainfall that is not evident at Kukkal. A high-rainfall phase at 3000–2300 yr BP is recorded in the discharge of the Indus River (Von Rad et al., 1999a, 1999b). However, this is at variance with a long drought in the Thar Desert, India, between 3600 and 2000 yr BP (Bryson and Bryson, 1996; Enzel et al., 1999; Swain et al., 1983) or aridity and weakening of the summer monsoon between 4500 and 2000 yr BP in the continental records (Swain et al., 1983). The lake records show a maximum Holocene wet period from 6000 to 4000 yr BP and a change to a drier climate after ~4000 yr BP. It is possible that some of the moisture at these sites in the Thar desert received the winter precipitation, which might explain as to why they are not in synch with other monsoon sites that respond to summer precipitation.
Climatic events of global scale during the last 3000 years include the ‘Roman Warm Period’ (‘RWP’) (2500–1600 yr BP (Wang et al., 2012, 2013), the ‘MWP’ (790–620 yr BP) and the ‘LIA’ (330–80 yr BP) (Bianchi and McCave, 1999). These phenomena have been proposed on the basis of multi-proxy studies in records from higher latitudes (e.g. Denniston et al., 2000; Lehmkuhl, 1997; Liu et al., 2013; Sukumar et al., 1993; Von Rad et al., 1999a, 1999b; Zhou et al., 1991, 2011). Crowley and Lowery (2000) described the ‘LIA’ as a mild cooling event in the Northern Hemisphere. Hu et al. (2001, 2008) determined that the ‘LIA’ led to ~1.7°C of cooling around Alaska (Keigwin, 1996). The ‘LIA’, therefore, has a large temporal and spatial variability across the globe (Intergovenrmental Panel on Climate Change (IPCC), 2001). The stepwise event at 3000–1800 yr BP in the Kukkal proxies appear to anticipate the ‘RWP’ by a few hundred years, but the interpolated dates are only approximate; the inflections in the curves may in fact correspond with the initiation of the ‘RWP’. There are no events at Kukkal that correspond convincingly to the ‘MWP’ and the ‘LIA’. Even though Zonneveld et al. (1997) suggested a rapid teleconnection between the climates of high latitude and the SWM, the only evidence of such relationships in the Kukkal core is the stepwise event that may correspond to the beginning of the ‘RWP’.
A possible scenario for Kukkal
When comparing and contrasting climatic proxy records from different environments, it must be borne in mind that each environment will have its own set of responses to a particular climatic forcing. Terrestrial records involving plant associations, in particular, may have much longer response times than the cave speleothems and marine sediments. In the case of the Kukkal record, it is possible to conceive of the evolution of the basin as the protracted response to a century or two of dry climate at about 8000 yr BP, followed by increasing rainfall. Savanna is initially replaced by grassland with palms, sedges and ferns appear later, and all are eventually succeeded by forest in a process that lasts several thousand years. Evolving vegetation and the amount of rainfall control the weathering and the character and availability of clastic sediments. Eventual establishment of mature Shola forest limits the energy of runoff and the availability of coarser clastics. Pollen preservation is poor following the period of dry climate, and recovers only when the SWM intensity increases more rapidly, augmenting the amount of standing water in the wetland/lake basin.
Conclusion
In summary, the following points emerge from this study:
The Kukkal sediment core reflects lacustrine or wetland depositional environments of early-Holocene age (9000 yr BP) to modern age. Four radiocarbon ages were used for producing an age model showing a slow average sedimentation rate of 0.11 mm/yr, except in the upper 14 cm of core, where the rate was 0.4 mm/yr or greater. The increased rate may reflect human activity because of deforestation, urbanization and agricultural practices around the lake.
Rapid climate change is indicated by the degree of preservation of OM and changes in its composition. Superimposed on this trend are two stepwise changes, at about 8000 yr BP and between about 3000 and 1800 yr BP. The former change appears to correspond to a dry spell of 100–150 years recorded elsewhere in India. The latter corresponds to a rapid intensification of the SWM, and may correlate with the initiation of the ‘RWP’, the two phenomena coinciding in time. The ‘MWP’ and the ‘LIA’ are not clearly represented in the record.
The C/N ratio of Kukkal Lake sediments ranges from 14.02 to 8.31, indicating that the organic carbon originates from a mixture of lacustrine algae, vascular plants and terrestrial plants. Stable carbon isotope data suggest early dominance of C4 plants such as sedges, and a progressive increase in the ratio of C3 vegetation.
Chemical weathering indices indicate extreme silicate weathering over the past 9000 years, with a stepwise increase in intensity beginning about 3000 yr BP. Palynological results reveal a succession of vegetation types, beginning with savanna prior to about 8000 yr BP, succeeded by grassland with palms and Asteraceae species. A vegetational change is observed just prior to 3200 yr BP, with ferns being abundant and the Shola forest established between 3200 and 1800 yr BP.
Footnotes
Funding
The author(s) received no financial support for the research, authorship and/or publication of this article.
