Abstract

The Indian monsoon activity, coinciding with the Inter-Tropical Convective Zone (ITCZ), progresses from the southern Indian Ocean during the boreal summer and withdraws towards the south in winter. Islands situated to the south of India receive, therefore, the first and last showers of the monsoon; speleothems in such islands have not yet been explored for their potential to reconstruct past monsoon rainfall. Here, we present the first measurements of stable carbon and oxygen isotopic compositions (δ13C and δ18O) of a stalagmite collected from the Baratang Island of Andamans, along with new data on δ18O of modern monsoon precipitation (May to July 2010). The aim was to detect (i) whether these samples are amenable to dating using 14C, (ii) whether their oxygen isotopes indicate precipitation under isotopic equilibrium, and (iii) if (i) and (ii) above are true, can we reconstruct monsoon activity during the past few millennia? Our results indicate that while δ18O of speleothem does show evidence for precipitation under isotopic equilibrium; dating by 14C shows inversions due to varying contributions from dead carbon. The present work highlights the problems and prospects of speleothem paleomonsoon research in these islands.

1. Introduction

The Indian monsoon is a mesoscale weather system affecting a very densely populated area of the globe. It is important to attempt high-quality paleoclimate reconstructions for a better understanding of the monsoon. The dynamics of Indian monsoon are not well understood due to lack of high-resolution and long-term monsoon records. The longest instrumental record of monsoon rainfall available spans the last ~150 years [1] and, therefore, too short to fully reveal monsoon variability on century and longer time scales, for which limited high resolution, geographically more widely distributed, paleomonsoon records based on speleothems [25], ice cores [6], tree-rings [7], and marine sediments [811] are currently available.

Speleothems are secondary carbonates precipitated in karst regions that record the paleoenvironmental changes in atmospheric circulation as the oxygen isotopic ratio (δ18O) of speleothems is primarily controlled by rainfall δ18O in the tropics [35], as the latter depends on the amount of rainfall (the amount effect); the δ18O in rain decreases with increasing amount of rainfall. On seasonal time scales, variations of δ18O in precipitation arise from variations in the source of rainfall and its characteristics such as the frequency, the duration, and the intensity of rain. On longer timescales (centennial to millennial), additional factors such as changes in the δ18O of the ocean (the ice volume effect), changes in the seasonality of precipitation (e.g., change in the relative proportions of winter and summer precipitation), and shifts in the source of moisture may be important (e.g., [12]). Despite these complications, δ18O of speleothem calcite is currently the most frequently used proxy as it can help reconstruct past precipitation on time scales ranging from centuries to millennia.

The carbon isotopic composition (δ13C) of speleothems is somewhat complicated to interpret. Nevertheless, changes in δ13C in most cases are the likely result of climate change [13], and a combined study with δ18O may reveal important aspects of climate and vegetation changes above the cave [14, 15]. δ13C in a speleothem is governed by the sources of dissolved carbon in the drip water, which are atmospheric CO2, soil CO2, and carbonate bedrock, the relative contributions of which depend on the mechanisms controlling the bedrock dissolution and carbonate precipitation in the cave [16, 17]; δ13C of speleothems can be linked to climate-driven vegetation change when the most carbon in the deposited carbonate derives from soil CO2 under isotopic equilibrium; the contribution from the bedrock limestone and atmospheric CO2 are relatively insignificant. This is because δ13C of soil CO2 is controlled by the plant communities exhibiting significant variations in δ13C. C3 plants like trees with δ13C value ranging from −19 to −32‰   (average −27‰) generally occur in cold and humid climates while C4 plants such as tropical grasses, with δ13C ranging from −9 to −19‰  (average −13‰) grow in warm and water-stressed conditions [18, 19]. The soil CO2 under C3 and C4 vegetation has typical δ13C values of −22‰  and −8‰ , respectively. Carbonate precipitated under isotopic equilibrium is enriched in 13C by ~10‰  with respect to soil CO2 at 25°C [20]. Therefore, speleothem precipitated under exclusively C3 and C4 vegetation should have δ13C values of ~−12‰   and +2, respectively. This large difference allows the reconstruction of shifts in vegetation type from forests to grasslands and vice versa [21]. However, one difficulty could be that the δ13C of speleothem is the integrated and weighted average contribution from the CO2 of the whole soil column; the deeper soil CO2 could be thousands of years old, and, therefore, it is difficult to decipher any short term climate change using the δ13C of speleothems.

Possible kinetic isotope fractionation must be considered before interpreting δ13C variation in terms of paleoenvironmental changes, because this can override changes in organic productivity associated with the soil zone [22]. Evaporative conditions may lead to nonequilibrium fractionation causing a large variation in δ18O and a strong linear correlation between δ13C and δ18O [17, 23]. CO2 degassing in the unsaturated zone above the cave may lead to enrichment of 13C in the deposited carbonate [24]. The kinetic fractionation, if happens during calcite precipitation, can be checked by the Hendy’s test: both δ13C and δ18O would increase along a particular growth layer with distance from the central axis [17]. However, this is not a foolproof technique [25] as it is difficult to pick samples from a particular layer. Also kinetic effect may cause more enrichment at the farther distances from the central growth axis.

Unlike the speleothems growing in mainland India, which indicate past monsoon overland, those from the caves of Andaman and Nicobar Islands (Figure 1) can be a very good source to study the oceanic monsoon variability. Monsoon activity begins in the boreal summer from the southern hemisphere, and the Andamans are the first to get rainfall during the monsoon season. Likewise during the retreat of the monsoon in winter, they are the last to get rain. The monthly temperature and rainfall are shown in Figure 2. Further, the industrial pollution effects are smaller in the Andamans relative to mainland India.

We present here the first stable isotope records of a stalagmite from the Baratang cave of Andaman Islands, dated using radiocarbon. The aim was to detect (i) whether these samples are amenable to dating using 14C, (ii) Whether their oxygen isotopes indicate precipitation under isotopic equilibrium, and (iii) if (i) and (ii) above are true, can we reconstruct monsoon activity during the past few millennia?

2. Dating Speleothems

The best way to date is to count the annual layers [5] of actively growing speleothems. However, it is difficult to judge whether a speleothem is growing actively or not, especially if the collection is being made in the dry season. It is also more difficult to spot a speleothem with annual layers; for this, a rapidly growing speleothem in a region with high seasonal precipitation (e.g., monsoon) is necessary. The most widely used dating method for speleothems is using U-series, specifically “daughter deficiency” in 230Th/U. A majority of speleothems behaves as a closed system with respect to uranium and its decay products. Dating is achieved using a Thermal Ionisation Mass Spectrometer (TIMS) [26, 27] or a Multi-Collector Inductively Coupled Plasma Mass Spectrometer (MC-ICP-MS) [28, 29]. But for samples that are only a few millennia old, this could give results with large uncertainties because of contribution from detrital Th. Radiocarbon dating of speleothem is also useful for modern samples when the dilution due to dead carbon can be estimated [4]. One way to quantify the dead carbon is by using both 230Th/U and 14C chronologies simultaneously [16].

3. Materials and Methods

3.1. Study Area and Present-Day Climate

The Baratang Island is located around 12°05′N, 92°45′E in the junction of middle and south Andaman, about 100 km north of Port Blair (Figure 3). A limestone cave of size ~50 m with an elongated opening (height ~10 m, width ~1 m) housed speleothems. A stalagmite of length ~20 cm was collected from the cave in 2005. The present-day vegetation above the cave is dense tropical rain forest (C3 type) with a thin soil layer. The height of the cave above sea level is ~18 m. Meteorological data from Port Blair (the nearest available station to the cave site) shows that the present-day climate is humid with a mean annual precipitation of 3168 mm (average of 30 years from 1951–1980). The mean annual temperature and humidity are 26.6°C and 82%, respectively [30].

3.2. Sampling and Analytical Methods

For stable isotope study, subsamples of about 1 mg were collected along the growth direction of the stalagmite at an interval of about 2 mm using stainless steel drill bits of size 0.5 mm. Each powdered sample was reacted with 100% H3PO4 for ~8 hours at 25°C. The CO2 produced was extracted, purified cryogenically, and analysed for δ13C and δ18O using dual-inlet stable isotope ratio mass spectrometer (GEO 20-20). A laboratory marble standard was run after every 10 measurements to check the reproducibility of measurements. The precision for the δ13C measurements, for the procedure followed, is better than ±0.1‰  [31]. The accuracy of the measurements is checked occasionally by running an international standard (NBS-19 supplied by IAEA), for which the long-term (~5 years) average δ13C = 1.95 ± 0.08‰   and δ18O = −2.3 ± 0.1‰  (relative to PDB) (error at 2 sigma level). All the carbonate isotopic values are reported relative to the VPDB standard [32].

Rainwater samples were collected from the cave site during summer 2010 (from May to July). δ18O of water samples were measured using CO2 equilibration method [31]. All values are reported relative to VSMOW, and the analytical uncertainty was better than 0.1‰. For reproducibility check, some samples were repeated, and the results of the replicates were well within the uncertainty of measurements.

For radiocarbon dating, about 5–7 g of carbonate were taken from different parts of the stalagmite along its growth axis (shown as patches in Figure 4), reacted with 85% H3PO4 to produce CO2 which was then quantitatively converted to benzene. The residual activity of the benzene was measured in a liquid scintillation counter (Quantulus LKB-1220). Further details of experimental procedures are discussed elsewhere [31].

4. Results

Four radiocarbon ages obtained for the stalagmite are shown in Figure 4 and Table 1. The ages are anomalous at the tip and the base; the ages are the same within 1σ uncertainty of radiocarbon ages. Therefore, the speleothem must have grown during the past few millennia. More precise dating could be possible by U-Th mass spectrometry, but only for older samples.

To check if any kinetic fractionation had occurred during the calcite precipitation, we performed the Hendy’s test (Figure 5). There are no significant trends obtained in the δ18O indicating that the oxygen isotopic equilibrium between calcite and drip water was maintained during stalagmite precipitation. A little enrichment in δ13C (Figure 5(a)) on both sides of the central axis shows a little kinetic fractionation during the degassing of CO2 leading to an enrichment in 13C in the residual carbonate. Also, the plot between δ13C and δ18O (Figure 5(b)) does not show a strong correlation (not significant at 0.01 level) indicating limited evaporation during stalagmite deposition. This is not surprising as the island is close to the sea and the cave is likely humid for most of the year.

A total of 91 stable isotope measurements were carried out on the stalagmite. The δ18O values of rain water collected between May to July, 2010 are shown in Figure 6. They vary between 0.7 to −5.2‰  with respect to VSMOW. The relationship between fractionation factor () and temperature at which calcite precipitates is given by the following equation [33]: where , is the temperature in Kelvin. At 26.6°C, the fractionation factor = 1.02796 and hence enrichment factor (α – 1)103 = ε = 27.96‰. Thus, carbonate precipitated under isotopic equilibrium with rain water should be enriched by ~28‰. Considering δ18O of rain water varying between 0.7 to −5.2‰  with respect to VSMOW, the stalagmite δ18O should vary between 28.7 to 22.5‰ , relative to the same standard, VSMOW. The relation between δ18O of the VSMOW and PDB standards are given by [34, 35]

It is observed that the δ18O of stalagmite (Figure 7) mostly varies between −4 and −6‰  relative to PDB. Using the above relations, we get stalagmite δ18O relative to VSMOW varies between 26.8 and 24.7‰  which is well within the predicted range for the present-day calcite precipitating under isotopic equilibrium. Therefore, the southwest monsoon in the past few millennia was quite similar to that at present. Small variations in δ18O shows little change in the mean annual precipitation throughout the deposition history.

δ13C also shows very little variation throughout the deposition history except at the tip. The δ13C value lies between −10 and −12‰  for most of the samples while it is somewhat higher at the tip.

5. Discussion

Fast growing speleothems are suitable for the high-resolution paleomonsoon reconstruction. Such speleothems could be useful for the reconstruction of annual and even seasonal rainfall variations in the past. Unpredictable dead carbon contribution makes it difficult to date a speleothem using the 14C method. Dead carbon is the carbon present in the speleothem that originated from the old (typically millions of years) limestone bedrock having no radiocarbon. Therefore, the 14C age of speleothem is overestimated (the reservoir effect [36]). In the case of an actively growing speleothem, the dead carbon contribution can be estimated using the tip age, because of the presence of bomb carbon (e.g., [5]). Assuming constant dead carbon contribution throughout the calcite deposition, a chronology can be established [5, 37]. However, the assumption of constant dead carbon contribution needs to be checked. The ages of the stalagmite listed in Table 1 show an anomalous behaviour, varying from 2500 to 3200 yrs BP (uncalibrated ages). This could be due to the changes in the dead carbon contribution during stalagmite deposition. Such material could also cause problems in U-Th dating because the increase in detrital matter introduces large errors. Yadava et al. [5] estimated −880 and −1230 yr as age corrections for the addition of dead carbon in two stalagmites from Central India. Hendy et al. [17] estimated corrections of 530–715 yrs in a stalagmite and 775–880 yrs in another from Villars cave, southwest France and Han-sur-Lesse cave, Belgium, respectively. Yanjun et al. [38] showed a dead carbon fraction of 10–15 % in a stalagmite from Shihua Cave, Beijing, leading to an age overestimation of 928–1500 yrs. From all these comparisons we conclude that the dead carbon contribution is not more than 15% in most cases, for samples of recent geological age. Therefore, the dead carbon could have overestimated radiocarbon ages by 1000–1500 years in the present sample.

δ18O for the present stalagmite mostly ranges between −4 and −6‰  excluding a few sharp peaks (Figure 7). By and large, we can conclude that the rainfall fluctuations were not more than 200 mm during the past few millennia, assuming an amount effect of −1.5‰  per 100 mm increase of rain [5, 36].

A comparison with other paleomonsoon proxy records from regions affected by ISM is shown in Figure 8. The comparison is made only for the period 3200–2500 yrs BP, the approximate growth period of the stalagmite and for more recent time to compare with the present-day monsoon. Yadava and Ramesh [4] reconstructed the southwest monsoon rainfall variation for the last ~3400 yrs based on a stalactite δ18O from Orissa (location shown in Figure 1), and this is compared with the present results (Figure 8(A) and 8(B)). The monsoon prior to 2500 yrs BP was similar to that of the last 0-400 yrs, a similar finding with the present stalagmite. However, there was a large variation in monsoon between 2100-1500 yrs (not shown in Figure 8(A)). Von Rad et al. [39] reconstructed a high-resolution (~7 yr) paleomonsoon record using thickness variations in the varved sediments collected from the northern Arabian Sea (for location, see Figure 1). The varve thickness is controlled by precipitation and river runoff, and, hence, the variability in varve thickness was interpreted as a proxy for past rainfall variations [40]. Lower thickness of the varved sediments is associated with the lower discharge from the Indus river and vice versa. A comparison of varve thickness variation with stalagmite δ18O from the Baratang Island is shown in (Figure 8(C) and 8(B)). Though there were variations in monsoon recorded in the varved sediments between 2500 and 400 yrs BP, but the monsoon after 400 yrs BP is similar to that present between 3200–2500 yrs BP (Figure 8(C)).

δ13C in the stalagmite mostly varies between −11 and −13‰, a dominant C3 signature since carbonate precipitated under equilibrium with soil CO2 under C3 vegetation should have δ13C value around −12‰  at ~26°C (the mean annual temperature). The bedrock carbonate (δ13C~0‰) contribution shifts the calcite δ13C values towards more positive values. Therefore, the bedrock carbon contribution is not very high in the stalagmite except at the tip. A gradual enrichment in 13C towards the tip shows an apparently larger bedrock contribution. The possibility of shift in vegetation from C3 to C4 at the later phase of the stalagmite growth cannot be ruled out. Probably tropical grasses (C4) grew at the surface causing enriched 13C in the stalagmite carbonate. But this is less probable as there is no major change in the precipitation (inference based on δ18O) observed during this time to cause a water-stressed condition. Therefore, we speculate that the 13C enrichment at the tip was caused by apparent higher bedrock contribution because of lesser soil CO2. This needs further work and verification.

6. Conclusions

This is the first isotopic study of speleothems reported from Andamans. The summer monsoon over the Andaman Islands was similar to the present a few millennia ago; the southwest monsoon was the main source of annual precipitation. The radiocarbon dating of speleothem has highlighted problems associated with incorporation of an unknown proportion of dead carbon during calcite precipitation.

Acknowledgments

The authors thank ISRO-GBP for funding and Mr. Felix for collecting samples of rainwater at Andamans. The authors are grateful to Professor M. S. Srinivasan for kindling our interest by sending us excellent photographs of the Andaman caves and Dr. J. S. Ray for help during field work.