Abstract
Paleoclimate investigations of the peat deposits in the Nilgiri Plateau, an important paleoclimate archive of India, are mainly restricted to the carbon isotope composition ((δ13C) of plant-derived materials and pollen studies. However, it is unclear whether these proxies reflect past variability in temperature or hydrology. Here, we report the hydrogen and carbon isotopic variability of n-alkanoic acid of chain length 28 (δDC28 and δ13CC28, respectively) and demonstrate that the peatland leaf wax hydrogen isotopes provide a sensitive record of past hydrology. The decoupling of δ13CC28 and δD of vegetation-corrected rain during the Holocene indicate that δD of the leaf wax compounds mainly respond to past hydrological variability whereas δ13C variations might reflect the temperature-controlled variability of C3 and C4 vegetation. Conforming with the other paleoclimate records from the region, the δDC28 variations showed a reducing precipitation trend since the early Holocene. However, a large amplitude of reconstructed δD of rain (~44‰) during the Holocene indicated changes in the moisture source and trajectory could be an additional factor contributing to the orbital-scale δD variability of proxies from the Indian region.
Introduction
The Holocene paleomonsoon records provide an opportunity to assess the causal link between well-established climate forcing and monsoon (Banerji et al., 2020; Kaushal et al., 2018; Mayewski et al., 2004; Misra et al., 2019; Morrill et al., 2003; Zhang et al., 2011, 2019). For understanding continental rainfall variability and underlying forcing, it is preferable to use terrestrial rainfall proxy archives such as lake sediments, peat bogs and speleothems. The lake sediments (Misra et al., 2019) and speleothem (Kaushal et al., 2018) records have been used to reconstruct paleomonsoon variability in the Indian region. However, due to the limited availability of appropriate sampling sites, the terrestrial records that offer sampling resolution and cover a period enough to reveal a trend in the Holocene precipitation are conspicuously absent from India.
The Nilgiri Plateau, southern India, hosts paleoclimate archives such as peat deposits that are ideal for the preservation of pollen (Sutra, 1997; Sutra et al., 1997; Vasanthy, 1988) and other organic matter (Ramya Bala et al., 2022; Sukumar et al., 1993), soil profiles (Caner et al., 2007), and lake sediments (Raja et al., 2018). In contrast to the other archives in the Indian region, peat deposits offer a continuous record of >40 kyr (Rajagopalan et al., 1997; Ramya Bala et al., 2022) and are important to reveal the orbital-scale climate variability and the factors influencing it. The pollen records from the Nilgiri Plateau have established not only the long-term stability of the Shola-Grassland mosaic but also the heterogeneity of landcover change, going back >30,000 years (Sutra, 1997). The carbon isotopic composition (δ13C) of peat (Rajagopalan et al., 1997; Ramya Bala et al., 2022; Sukumar et al., 1993) and soil profiles (Caner et al., 2007) from the Nilgiri Plateau, reflecting the relative proportion of C3 and C4 plants in the catchment area, has been used to reconstruct the past climate variability. However, as the trees to grass ratio in the Nilgiri Plateau is controlled by the winter frost (Joshi et al., 2020), δ13C and pollen-based proxies might not reflect past precipitation variability alone. Further, short-term disturbances such as fire and local attributes led to divergence in the δ13C records of peat cores (Ramya Bala et al., 2022). Thus the potential of peat deposits from the Nilgiri Plateau in reconstructing past rainfall variability is yet to be fully realized.
The deuterium to hydrogen ratio (δD) of the leaf waxes is known to preserve δD of concurrent environmental water or precipitation (Hou et al., 2008; Huang et al., 2004; Sachse et al., 2012). δD of n-alkanoic acid has been extensively used to reconstruct past precipitation isotope variability (Konecky et al., 2011; Niedermeyer et al., 2010; Pagani et al., 2006; Shuman et al., 2006; Tierney et al., 2008). Although several studies have used δD studies of leaf wax compounds in the peat bogs from higher latitudes (e.g. Daniels et al., 2017; Seki et al., 2011), its potential in reconstructing past rainfall variability in the tropics is not adequately demonstrated. Here, we report the results of δD and δ13C-based investigations of n-alkanoic acid fractions of the leaf wax in the peat deposits from the Nilgiri Plateau. The influence of change in the vegetation composition and seasonality in rainfall on the δD of reconstructed rain was assessed. The amplitude of the Holocene-scale changes in δD of reconstructed rain was used to infer the past changes in precipitation amount and associated atmospheric circulation.
Materials and methods
Sampling site
The Nilgiri Plateau in southern India (Figure 1) receives rainfall during the southwest monsoon (SWM) and northeast monsoon (NEM) from June to September and October to December, respectively. The Western and Central Indian Oceans are the major sources of moisture during the SWM (Pathak et al., 2017) (Figure 1a). The rainfall and its seasonality are spatially variable (inset in Figure 1b): the western part (Region I) is dominated by the SWM and receives 2500–5000 mm of rain; the central region (Region II) receives both the monsoons but with a lower annual total (900–1200 mm) and the eastern and southern parts (Region III) receive a higher contribution from the NEM and a lower annual total (1500–2000 mm) (Caner et al., 2007). The sampling site Avalanche–B (76°35′58.40″E, 11°18′06.03″N, ~2000 m amsl, Figure 1) is located in Region I. The contribution of the SWM rainfall to the annual total at Terrace Factory (Region I), Avalanche (Region I), Ootacamund (Region II) and Coonoor (Region III) are 75%, 71%, 42%, and 21%, respectively. The mean monthly temperature variations over the plateau are shown in Figure 1b. Regions I and II are more prone to frost formation which can occur up to 50 days from December to March with its frequency the highest during January (Von Lengerke, 1977). The most important factors controlling vegetation dynamics between trees and grasses here are the daily thermal differences and the extreme temperatures (Joshi et al., 2020; Legris and Blasco, 1969; Legris and Meher-Homji, 1977).

(a) Major sources of moisture to Indian summer monsoon rainfall and vertically integrated moisture flux divergence over the region (after Pathak et al., 2017). WIO: Western Indian Ocean; UIO: Upper Indian Ocean; CIO: Central Indian Ocean; SIO: Southern Indian Ocean. 1: Lonar Lake (Sarkar et al., 2015), 2: site U1446 (McGrath et al., 2021). Thick black and red arrows schematically indicate proximal and distal transport of moisture during the SWM, respectively. The thin blue arrow indicates the direction of the NEM. The study site is indicated by a star. (b) The Nilgiri Plateau. The three climatic regions (separated by the dashed lines) are as in Caner et al. (2007). Data for the plots in the inset (Von Lengerke, 1977) indicate mean monthly precipitation and temperature in the region. The sampling location Avalanche-B is shown by a star.
The vegetation in the Nilgiri Plateau is characterized by a mosaic of grassland, sedgeland and shola (i.e. stunted evergreen forest) (Blasco, 1971, 1987; Blasco and Tassy, 1975; Meher-Homji, 1967; Ramya Bala et al., 2022; Ranganathan, 1938). The grassland occurs on the hilltops; shola, in the sheltered valleys; and sedgeland, on the waterlogged valley floors (Ranganathan, 1938). While shola vegetation is mainly C3, grassland and sedgeland comprise both C3 and C4 vegetation (Ramya Bala et al., 2022). At the sampling site, the shola vegetation is in proximity in the south, east and west; and the grassland, in the north and northwest (Supplemental Figure S1, available online).
The isotopic characteristics of the modern-day rainfall
Seasonal isotopic variability
As the study area receives SWM and NEM, it is important to characterize their isotopic composition. Four Global Network of Isotopes in Precipitation (GNIP) (IAEA/WMO, 2022) stations and a few year-long investigations (Supplemental Table S1, available online; Deshpande et al., 2003; Lekshmy et al., 2015; Warrier et al., 2010; Yadava et al., 2007) reveal that the NEM is relatively more depleted in the heavy isotopes than the SWM. The differences in the δD of SWM and NEM at GNIP stations varied from 9‰ to 26‰. The δD of stream water at Ootacamund during July (−32‰) and October (−59‰) (Parthasarathy et al., 2021), gave an estimate of the δD of rainfall during the SWM and NEM, respectively.
Isotopic variability in the Nilgiri Plateau
Monthly or seasonal rainwater isotopic data are not available for the Nilgiri Plateau. The δD of surface and groundwater samples collected in the Nilgiri Plateau (Tripti et al., 2019) suggested −32‰, −45‰, and −53‰ as the mean annual δDrain for Region I, II, and III, respectively (See Supplemental Discussion, available online). The expected δD of rain at the sampling location (~2000 m a.s.l.) considering an elevation effect of −0.8‰/100 m and data given by Resmi et al. (2016) (Supplemental Table S1 available online) is −38‰.
The amount effect in the rainfall and the response of the proxy records
An inverse relationship between the amount of precipitation and its isotopic composition, the “amount effect” (Araguás-Araguás et al., 2000; Dansgaard, 1964), forms a basis for interpreting proxies, such as δDC28, that inherit the isotopic composition of rains. Based on the monthly rainwater isotopic data Lekshmy et al. (2015) showed that the amount effect is apparent only in the region where the amount of rain received during the NEM and SWM were comparable. Further, the daily rain samples revealed that the lower δ18O values reflect large-scale organized convection rather than the amount of local rainfall (Lekshmy et al., 2014).
Although the daily and monthly isotopic data are useful in understanding the “amount effect” in terms of atmospheric processes, their utility in interpreting the isotopic records of proxies that are more likely to reflect the isotopic signal of rainfall integrated over a year (e.g. in tree rings) or longer (e.g. in speleothem) is limited. Inter-annual variability in the isotopic composition of rain, more useful while interpreting the isotope-based proxies, at GNIP station Kozhikode showed an inverse relationship between the weighted δ18O of the annual precipitation and its amount (r = −0.62, N = 11, p < 0.05). Further, the inter-annual δ18O record of a teak tree (1943–1988) from southern India showed a negative correlation with the amount of concurrent rainfall (Managave et al., 2011) suggesting the utility of the isotopic composition of plant-derived material as a proxy for rainfall.
Sample and radiometric dates
The peat core analyzed was taken using a modified Russian Peat corer (semi-circular) during a field campaign carried out in 1995 by French Institute of Pondicherry (IFP) in the framework of the PhD thesis of Sutra (1997). The core of 7.5 m depth was sub-sampled at 1 cm resolution for pollen analyses and a part of the remaining samples stored at the Laboratory of Palynology and Paleoecology (IFP) was used in this study. Here we report the results of the analysis of a part of the core to a depth of 529 cm (corresponding to 9429 BP). The dates associated with samples from various depths and an age-depth chronology built using Clam (version 2.2) (Blaauw, 2010) were shown in Supplemental Table S2 and Figure S3, available online. The peat accumulation rate was 116 cm/kyr prior to 6.7 ka and 26 cm/kyr after this time.
The samples were selected to get a uniform temporal sampling during the Holocene; an average of ~270 years between the two samples. Further, as the thrust of the present study was hydroclimatic reconstruction using δD, the number of samples analyzed for δD (N = 35) was more than those analyzed for δ13C (N = 15).
Lipid extraction and the isotopic measurements
Peat samples were freeze-dried and labile organic fractions from them were extracted using DCM:MeOH in a 9:1 ratio in an Accelerated Solvent Extractor (Dionex ASE 200). Neutral and acid fractions were separated from the total lipid extract using aminopropylsilyl as a solid phase extractor. The neutral fraction was extracted first using DCM and Isopropanol (2:1, v/v). This was followed by the extraction of acid fraction using acetic acid in ether (4%, v/v). The acid fraction was later methylated using MeOH with acetyl chloride to produce fatty acid methyl esters (FAME) following Huang et al. (2004).
δ13C and δD of FAME were analyzed in Huang’s lab at Brown University by gas chromatography stable isotope mass spectrometry (GC-IRMS). Peaks associated with different chain lengths of FAMEs were identified from their retention times. For δ13C and δD, each of the samples was analyzed twice and thrice, respectively. The precision (1-sigma) of these analyses was better than 0.15 ‰ (for δ13C) and 2‰ (for δD). Laboratory isotopic standards of alkane and FAME were routinely measured to check the accuracy of the analysis which showed similar standard deviations for δ13C (±0.15‰) and
where
Similarly,
where
The effect of change in vegetation on δD of n-alkanoic acid of chain length 28 (δDC28)
An approach adopted by Konecky et al. (2016) was used to assess the effect of changes in the relative contribution of trees and grasses on δDC28. We first estimated the fraction of C4 plants using
where

(a) δDC28 variability in the peat core. Two horizontal gray lines indicate the mean of δDC28 during the early Holocene (shown in a blue rectangle) and the Late-Holocene (shown in a pink rectangle); associated dotted lines indicate 1-σ standard deviations. (b) δ13CC28 variability (the left axis) and estimated ƒ
The vegetation corrected apparent fractionation between
where
The vegetation corrected δD of rainfall (
To assess the effect of the choice end-member
Results and discussion
δD and δ13C variations of n-alkanoic acid and factors controlling them
δD of different homologs of n-alkanoic acid revealed an overall increasing trend during the Holocene (Supplemental Figure S4, available online). δDC28 variations showed a range of 44‰: the heaviest (−166‰) in the most recent times (0.6 kyr BP) and the lightest (−210‰) at ~8.6 ky BP (Figure 2a). δ13C of homologs of higher chain lengths showed the lowest δ13C values for a period from 7.7 to 9.4 ky BP. It showed the highest δ13C values for a period between 4.4 and 6.6 ky (Supplemental Figure S4, available online and Figure 2b).
δDC28 and δ13CC28 variability suggest some similarities (Figure 2). For example, around ~ 3 ka BP, and from 7.7 to 9.4 ky BP lower values of δDC28 and δ13CC28 can be seen. However, the Holocene-scale variations of both are different. The highest δ13CC28 values were seen during 4.4–6.1 ky BP whereas the highest δDC28 values were associated with the most recent sample. The decreasing trend in δDC28 observed during the Holocene is absent in the δ13CC28 variations.
Effect of vegetation changes
The
Figure 2c shows
Reconstructing δD of rainfall
The absolute
The
Effect of varying contributions from the NEM
The relative contributions from the SWM and NEM could influence δD of the annual rainfall. The pre-Monsoon, SWM and NEM contribute about 13%, 71%, and 16%, respectively, to the annual rainfall at Avalanche (the data from Von Lengerke, 1977). If Avalanche was to receive the same amount of NEM rainfall as at Coonoor, these contributions would change to 10%, 59%, and 31%, respectively. As long as the differences among δD of the pre-Monsoon, SWM and NEM remain the same, the increase in the contribution of the NEM from 16 to 31% decreases δD of the annual rainfall by ~4‰ (Supplemental Table S4, available online), which is much less than 44‰ range in δD28 observed in the present study. It appears that the variation in the contribution of the NEM might not have significantly influenced δDC28.
Comparison with other rainfall reconstructions
δDC28 trend revealed in this study follows the solar insolation (Figure 3a and b). The coherence of the Holocene-scale monsoonal reconstruction of the present study with that from other reconstructions from monsoonal Asia (Figure 3) validates the utility of δD of n-alkanoic acid of peat from Nilgiri Plateau as a proxy for rainfall.

The Holocene precipitation proxy records from the Asian region. (a) Solar insolation at 30⁰N. (b) δDC28 variability in the peat deposit of Nilgiri plateau (present study). (c) Hydroclimatic changes inferred from a multi-proxy study at Lonar Lake (Prasad et al., 2014). In (c), darker shades of blue and brown represent, respectively, wetter and drier climates. (d) δ13C variations of n-alkanoic acids from a sediment core that have a contribution from the Godavari Basin (Ponton et al., 2012). (e) δ18O variability speleothem from Qunf cave, southern Oman (Fleitmann et al., 2003). (f) δ18O record of a speleothem from Sahiya cave, India (Kathayat et al., 2017). Precipitation reconstruction based on the core retrieved from (g) Tianchi lake and (h) Gonghai lake (Chen et al., 2017). (i) δ13C record of Mayinghai lake sediments (Cheng et al., 2020). (j) δD record of Qinghai lake sediments (Thomas et al., 2016). δ18O record of a speleothem from (k) Dongge cave (Dykoski et al., 2005), (l) Heshang cave (Hu et al., 2008), (m) Jinfo cave (Yang et al., 2019), and (n) Sanbo cave (Wang et al., 2008). (o) Humification index in Hongyuan bog (Yu et al., 2006). (p) δ18O of bulk carbonate in Seling Co. (Gu et al., 1993). (q) a composite δ18O record of a speleothem from Kesang cave (Cheng et al., 2012). Note that the values of δD and δ18O are expressed in permil (‰) and are reversed. The blue rectangle indicates a wetter phase observed around 3 ka.
Synthesis of paleomonsoon variability reconstructed using various terrestrial and oceanic proxies from the Indian region show an overall drying trend during the Holocene (Misra et al., 2019; Banerji et al., 2020). The pollen records from the Nilgiri Plateau have suggested a particularly very wet early Holocene as compared to recent times (Sutra, 1997). The Holocene rainfall trend and the wetter and drier phases revealed in a sediment core from Lonar Lake (Prasad et al., 2014) match well with that reported here (Figure 3b and c). Increasing aridification in the Godavari basin during the Holocene (Ponton et al., 2012) (Figure 3d) corroborates the rainfall variability inferred here. The oxygen isotope variability of speleothem from Qunf cave (Fleitmann et al., 2003) (Figure 3e) also revealed a drying trend during the Holocene. The monsoon variability inferred from various ocean sediment cores from the Arabian Sea and Bay of Bengal also suggested decreasing monsoonal rainfall during the Holocene. These studies are based on (i) percent variations in Globigerina bulloides abundance (Gupta et al., 2003) in the western Arabian sea, (ii) high-resolution terrigenous proxy studies (Thamban et al., 2007), and oxygen isotopic variations (δ18O) of G. ruber (Saraswat et al., 2013; Tiwari et al., 2015) at the eastern Arabian sea, (iii) sea surface salinity records based on δ18O variability of G. ruber from the western Bay of Bengal (Rashid et al., 2011). The locations of the ocean sediment cores from the eastern Arabian sea in Saraswat et al. (2013) and Tiwari et al. (2015) are significantly influenced by the discharge of rivers originating in the Western Ghats. As the present study site also lies in the Western Ghats, the decrease in the rainfall during the Holocene inferred at these sites mutually corroborates each other.
The decreasing trend in precipitation during the Holocene has also been reported for various sites in China (Zhang et al., 2017, 2019; Zhao et al., 2021). δ18O records of the speleothems from Dongee (Figure 3k) (Dykoski et al., 2005), Heshang (Figure 3l) (Hu et al., 2008) and Sanbao caves (Figure 3n) (Wang et al., 2008) demonstrate the decreasing trend in rainfall during the Holocene. δ13C studies of peat bogs at Dahu (Zhong et al., 2010), Dajiuhu peat (Ma et al., 2008), Hongyuan peat (Hong et al., 2003) and Hani Peat (Hong et al., 2005) also exhibit weakening of the monsoon. A compilation of eight lake records from various parts of China also suggests the weakening of the East Asian Summer Monsoon during the Holocene (Zhang et al., 2011).
Our δDC28 record indicated a wetter phase during ~2.4–3.6 kyr BP; the degree of wetness, however, was less than that during the early Holocene (Figure 2a). Interestingly, a coeval wet phase was observed in multiple paleoclimate reconstructions from the Asian region (Figure 3c, f–j and m–q). Additionally, the abundance of G. bulloides in the eastern Arabain Sea (Saravanan et al., 2019) and δ18O records from Xianglong (Tan et al., 2018) and Dongshiya (Zhang et al., 2018) caves also exhibit a concurrent wetter phase.
Factors influencing δD of the reconstructed rainfall
Even though the extreme values of δDrain during the Holocene suggested a change of ~44 ‰, the mean isotopic change was of lesser extent. When the mean δDrain during the early (older than 8.2 ky) and the late (younger than 4.2 ky) Holocene are considered, this range reduces to 16 ‰ (Figure 2a). If the amount effect used by McGrath et al. (2021) (−28.19 mm/1‰ δD) was considered, 16‰ translated to ~22% more rainfall during the early Holocene than during the Late-Holocene. The relationship between annual rainfall and δ18O of tree-ring cellulose from Kerala (Managave et al., 2011) yielded a −24.7 mm/1‰ increase in δD which, when used, suggested ~20% higher rainfall during the early Holocene.
A modeling study suggested that annual δ8Ο of precipitation varied in the precession band (i.e. during high and low summer insolation scenarios) by 2.4–2.8‰ in the Nilgiri region (from Figure 6 in Tabor et al., 2018). This amplitude translated to a 19–22‰ increase in δDrain (i.e. 8 times δ18Ο signal). The modeled difference in annual δ8Ο of precipitation during high insolation (such as during the early Holocene) and modern conditions suggested ~16‰ increase in δD of precipitation (i.e. 8 times 2‰) in southern India (from Figure 7 in Battisti et al., 2014). It thus appears that the average change in δDrain during the Holocene is consistent with the model predictions.
The range of extreme values in δDrain during the Holocene, however, was ~44 ‰ which translated to ~60% and ~55% higher rainfall during the wettest phase of the early Holocene, considering the amount effect of −28.19 mm/1‰ and −24.7 mm/1‰, respectively. Further, the range of ~44‰ was 2 and 2.75 times higher than that reported by Tabor et al. (2018) and Battisti et al. (2014), respectively, for the region and was comparable to the glacial-interglacial range of 50‰ in δDrain reported by McGrath et al. (2021). It is also much higher than the range of 12‰ observed in the Holocene-level δD variability of n-alkane in a sediment core from the upper Bay of Bengal (Contreras-Rosales et al., 2014). The large amplitude of δDrain observed in the present study thus indicated a role played by factors other than the amount of precipitation in controlling δDrain during the wettest phase of the early Holocene.
The emerging view (Battisti et al., 2014; Hu et al., 2019; McGrath et al., 2021; Sarkar et al., 2015; Tabor et al., 2018) suggests that δ8Οrain (or δDrain) on an orbital scale in Asia is influenced by factors such as precipitation amount, moisture source and variable rainout during transport of vapor. The insolation-only isotope-enabled climate model indicated changes in the wind strength and directions could affect the relative contributions of vapor from different vapor source regions (Tabor et al., 2018). During precession minima, a higher contribution from a more distal moisture source, associated with more rainout during longer transport, led to the lighter isotopic composition of rainfall in the Indian region (Tabor et al., 2018). It is speculated that δDrain during the wettest part of the early Holocene reported here plausibly indicates changes in the circulation pattern. WIO is a major source of moisture for rainfall over the Western Ghats followed by the CIO (Pathak et al., 2017) (Figure 1a). The latter constitutes a more distal source than the former (Figure 1a). A stronger monsoon during the higher summer insolation was associated with moisture from more distal sources (Tabor et al., 2018). Therefore, during the wettest phase of the early Holocene the study area might have received a higher contribution of moisture from CIO (Figure 1a), and the rainout during the longer transport path led to a lower δD of the associated rain.
The implication to the peat-based past climate variability in the Nilgiri Plateau
The C3 and C4 plants on the Nilgiri Plateau have δ13C values of −26.6 ± 2.3‰ and −10.9 ± 0.97‰, respectively (Rajagopalan et al., 1999). However, considerable ambiguity exists regarding the interpretation of the relative abundance of C3-C4 vegetation inferred from δ13C-based studies. Initial investigations of δ13C records in the peat (Sukumar et al., 1993) were interpreted in terms of precipitation variability whereas Caner et al. (2007) stressed the role of temperature in controlling δ13C in the soil profiles. The recent experimental evidence (Joshi et al., 2020) suggests that winter frosting in the open grasslands hinders the establishment of native trees in the grassland. This suggests that the ratio of grass to trees in the region is mainly controlled by temperature. The role of local disturbance (fire) and topography-controlled hydrology was also invoked to explain the vegetation composition in the area (Ramya Bala et al., 2022).
Sukumar et al. (1993) and Ramya Bala et al. (2022) have carried out δ13C-based investigations in peat cores from Sandynallah, from Region II. In the latter, the δ13C variability varied depending upon whether the core is from an ecotone between shola and sedgeland (in Core 1, Figure 4c) or the ecoregion of stable sedgeland (in Core 2, Figure 4d).

Peat-based paleoclimate investigation in the Nilgiri. (a) δDC28 and (b) δ13CC28 variability in the peat core from Avalanche (this study). δ13C variations in (c) Core 1 and (d) Core 2 from Sandynallah (Ramya Bala et al., 2022).
In contrast to the distinct seasonality in rainfall, the seasonality in temperature between Region I and II is not pronounced (the inset in Figure 1). The climate records from Region I and II show overlapping average monthly temperature and the monthly average of daily minimum temperature (Supplementary Figure S5). If the frost-mediated temperature effect controls the C3 and C4 composition in the region, then the similar ecotones from Region I and II are likely to show coherent variability in δ13C. A period from ~7.5 to ~4 ka was characterized by the complete dominance of C4-type sedgeland vegetation (shown by a horizontal arrow) followed by the appearance of shola (shown by an inclined arrow) (Figure 4c). These trends can be observed in our δ13CC28 record as well (Figure 4c). The similarity in the δ13C variability in our study (Figure 4c) and in Core 1 (Figure 4c), both representing variations in the ecotone between shola and sedgeland, perhaps attests to the frost-mediated temperature control of δ13C.
At this stage, it is not clear whether the Holocene scale drying in Region I reported in this study occurred in Region II as well. If it had occurred, then this study can help in resolving the uncertain interpretations in Ramya Bala et al. (2022). The δ13C variability in the Core 1 was attributed either to local disturbance (such as a fire) or changes in the climate-initiated hydrological conditions. The hydrological reasons invoked by Ramya Bala et al. (2022) for explaining the changes in the vegetation composition in the Core 1 are at odds with the rainfall variability reconstructed based on δDC28 in this study (Figure 4a). The higher δ13C values during a period from 7.4 to 4 ka in the Core 1 were interpreted to reflect wetter conditions but our δDC28 record showed a relatively drier phase for the same period. An arid phase was suggested for explaining the lowering of δ13C values post 4 ka (shown by an inclined arrow in Figure 4c) but the concurrent phase was wetter in our study. It thus appears that C3-C4 variability in Core 1 may not represent the climate-initiated hydrological changes. Further, the role played by temperature and precipitation in controlling the relative abundances of C3 and C4 sedges in Core 2 was also not fully resolved. Interestingly, the δ13C trend in Core 2 and δDC28 variability show similar variations during the Holocene (Figure 4) plausibly suggesting the role of rainfall in controlling δ13C of Core 2. At this stage, it appears that the δ13C variability in the ecotone between trees and grassland is controlled by temperature whereas precipitation controls δ13C of stable sedgeland ecoregion.
To optimally utilize the peat-based proxies for reconstructing past climate variability, characterization of the isotope-based proxies from the Nilgiri Plateau is necessary. The end-member characterization of δD and δ13C of various leaf wax compounds in grasses, shola and sedges needs to be estimated. Based on the δD of the reconstructed rain of the latest sample, it appears that the plant community-averaged
Conclusion
The Nilgiri Plateau hosts important paleoclimate archives having the potential for reconstructing climate since the last glacial maxima. Its location in a monsoonal environment and its higher altitude, however, result in both rainfall and temperature controlling proxy response in the region. This work demonstrates that the leaf wax δD record of peat deposit provides sensitive and continuous records of past hydrology. It appears to be a better proxy for the reconstruction of past hydrological variability whereas δ13C-based proxies reflect the temperature-controlled variability in the composition of vegetation. This is the first report to show the precipitation over the Nilgiri Plateau followed the solar insolation during the Holocene and also the first to show the utility of δD studies of tropical peat in reconstructing past rainfall. However, the large amplitude of the reconstructed δDrain suggested rainout associated with a longer monsoon transport path could have also contributed to lower δDrain during the early Holocene.
Supplemental Material
sj-docx-1-hol-10.1177_09596836231183110 – Supplemental material for Holocene precipitation hydrogen isotopic values on Nilgiri Plateau (southern India) suggest a combined effect of precipitation amount and transport paths
Supplemental material, sj-docx-1-hol-10.1177_09596836231183110 for Holocene precipitation hydrogen isotopic values on Nilgiri Plateau (southern India) suggest a combined effect of precipitation amount and transport paths by Shreyas Managave, Yongsong Huang, Jean-Pierre Sutra, Krishnamurthy Anupama and Srinivasan Prasad in The Holocene
Footnotes
Acknowledgements
We thank Lekshmy PR for providing rainwater δD data. Sarah M. McGrath is acknowledged for her critical suggestions. Devi Maheshwori’s help during the preparation of Figure 1 is acknowledged. We thank the Tamil Nadu Forest Department for permission to collect the sediment core analyzed in this study. T. Gopal (IFP) took an active role in extracting the core on the field.
Funding
The author(s) disclosed receipt of the following financial support for the research, authorship, and/or publication of this article: This work is supported by funding from UGC, New Delhi under Raman Fellowship (Grant No. F. 5-104/2014(IC)).
Supplemental material
Supplemental material for this article is available online.
References
Supplementary Material
Please find the following supplemental material available below.
For Open Access articles published under a Creative Commons License, all supplemental material carries the same license as the article it is associated with.
For non-Open Access articles published, all supplemental material carries a non-exclusive license, and permission requests for re-use of supplemental material or any part of supplemental material shall be sent directly to the copyright owner as specified in the copyright notice associated with the article.
