Abstract
High-resolution analysis of a 3.80 m sediment core recovered from Deoria Tal, a mid-elevation lake located at 2393 m a.s.l. in the Garhwal Himalaya, documents long-term and abrupt hydroclimate fluctuations in northern India during the mid- to late Holocene. The sediment chronology, based on ten 14C dates, indicates the core spans 5200 years. Non-destructive, radiological imaging approaches (X-ray fluorescence (XRF), X-ray imaging, and CT scans) were used to assess the response of the lake system to changing hydroclimatic conditions. Variations in elemental concentrations and sediment density evidenced notable hydroclimate change episodes centered at 4850, 4200, and 3100 cal yr BP. Elevated detrital input, greater sediment density, decreased lake ventilation, and lower autochthonous productivity reflects lake deepening between 4350 and 4200 cal yr BP. An abrupt shift in elemental concentrations and sediment density indicated the onset of lake drawdown at 4200 cal yr BP and a negative hydroclimate anomaly between 4200 and 4050 cal yr BP. Lower detrital flux, decreased sediment density, increased oxygenation, and higher autochthonous productivity, reflects a reduction in lake volume between 3200 and 3100 cal yr BP. The potential link between abrupt climate change at 4200 cal yr BP and the contraction of the Indus civilization is explored.
Similar content being viewed by others
Introduction
The 4200 cal yr Before Present (BP) megadrought, which extended from 4200 to 3900 cal yr BP, was characterized by a reduction in the global monsoon and associated circulation systems, leading to droughts and seasonal precipitation failures worldwide1. This event, which included a reduction in the strength of the Mediterranean Westerlies2 and a widespread drought over much of the interior of North America3, also impacted Asia4,5,6. The East Asian Summer Monsoon (EASM) and the Indian Summer Monsoon (ISM) experienced periods of weakening in eastern China, Inner Mongolia, and India, “synchronous” with the 4200 cal yr BP event1. Given that the EASM and the ISM are the main sources of precipitation for their respective regions, climate scientists and archaeologists have postulated that a dramatic reduction in the strength of the Asian summer monsoon at 4200 cal yr BP would have impacted the complex societies present in India and China during the mid- to late Holocene7,8.
At 4200 cal yr BP, northern India was occupied by the well-established Indus civilization9,10. The Indus cvilization, also referenced as the Indus civilization, was a socially complex, agrarian and highly urbanized civilization located along the Indus and Ghaggar-Hakra river valleys (Fig. 1). During the Late Mature phase15 (between ~ 4200 and ~ 3900 cal yr BP), the Indus civilization was characterized by well-developed infrastructure including large cities with extensive sanitation facilities16 and agricultural and water delivery systems heavily reliant on hydraulic engineering17. A distinguishing Late Mature phase characteristic is that it spans the onset of the Indus civilization’s “collapse” or “metamorphosis”, which involved demographic and spatial contraction comprised of de-urbanization coinciding with an eastward migration15,18,19,20. Scholars have proposed a number of causes to account for the decline and eventual disappearance of the Indus civilization, including hydroclimate anomalies, disease, warfare, tectonic-induced changes in the course of major river systems, socio-political instabilities, changes in crop patterns, and/or the disruption of trade networks21. Although the specific cause remains unknown, archaeological and paleoclimate studies have identified abrupt climate change and hydroclimate variability, specifically an increase in climate-induced drought episodes from ~ 4200 cal yr BP onwards, as a trigger for the decline of this highly complex civilization11,18,22,23,24.
The ISM is the main source of precipitation for Peninsular India, which includes the region formerly occupied by the Indus civilization. Variations in the strength and spatial extent of the ISM on decadal timescales would have greatly influenced this agrarian society23. Previous paleoenvironmental and paleoclimate studies assessing potential links between the Indus civilization’s re-organization and climate include lake sediment-based reconstructions from northern India25,26,27,28,29,30,31, marine sediment records from the Arabian Sea and Bay of Bengal11,18,20,32,33, and stalagmites4,13,34 from the region. Other approaches include relating climatologically driven fluvial dynamics to settlement locations19,35,36 and inferring climate based on agricultural requirements17. Modeling and paleoenvironmental studies suggest that the rapid, step-wise, and sustained decline in population that characterized the Indus civilization, beginning at ~ 4200 cal yr BP, can be attributed to increasing water scarcity associated with a weakening of the ISM and changes in regional hydroclimate23,24,29. However, existing records provide inconsistent and sometimes-conflicting evidence of the timing and nature of climate and environmental changes during the mid- to late Holocene, limiting our ability to explicitly link hydroclimatic variability with the archeological evidence of civilization and cultural change in the region4.
This study was undertaken to determine if the decline in the Indus civilization occurred during an interval characterized by a sustained decrease in precipitation and increasing hydroclimate variability. Here, we present a sub-decadal-scale reconstruction of hydroclimate variability between 5200 and 200 cal yr BP. We focus on the interval between 5200 and 3000 cal yr BP, which spans the time when the Indus civilization flourished and declined. We utilize non-destructive, radiological approaches, e.g. X-ray fluorescence (XRF) and Computerized Tomography (CT) scans of a well-dated lake sediment core recovered from Deoria Tal, a small, closed basin lake located in the Garhwal Himalaya to: (1) characterize regional hydroclimate variability during the mid- to late Holocene; (2) determine if evidence of an abrupt, short-lived climate event at 4200 cal yr BP is present at this site; and (3) assess the correspondence between hydroclimate variability and Indus civilization re-organization during the late Holocene.
Study site
Deoria Tal (30.5222° N, 79.1277° E) is located on a hilltop (2393 m a.s.l.) in the northern Indian state of Uttarakhand (Fig. 1). The lake is situated in a high-grade metamorphic terrain above the Main Central Thrust (MCT) zone in the Garhwal Himalaya and is underlain by three major lithotectonic groups37,38. The lithotectonic groups are delineated by three main thrusts—Bhulkund Thrust, Okhimath Thrust, and Banswara Thrust (from north to south). Deoria Tal is located between the Bhulkund Thrust and the Okhimath Thrust and is underlain by porphyritic gneiss with mica schist and granite also present39 (Fig. S1). Like quite a few other lakes in the Indian Himalaya40,41, Deoria Tal appears to have formed as a result of neo-tectonic activity.
The majority of the Indian subcontinent, which is characterized by warm, wet summers, is largely influenced by the ISM, with approximately 70% of the total annual precipitation falling between June and September42. The seasonal migration of the Intertropical Convergence Zone (ITCZ) and its movement southward during the late boreal summer limits the amount of precipitation that falls over most of India during the winter and spring, with the northeast coast and northwestern India being the exception43. Lying south of the ITCZ, the major source of precipitation for Deoria Tal is the ISM with average annual rainfall totaling around 200 cm44 and daily amounts of precipitation ranging between 8 and 10 mm/day45. Deoria Tal also receives lesser amounts of precipitation from extra-tropical cyclones, also known as Western Disturbances, between December and February46. Western Disturbances, which are frontal systems that originate over the North Atlantic, Mediterranean, and Black Sea regions during the boreal winter and spring, contribute a significant amount to total annual precipitation in the far north and northwest of India47. On inter-annual to inter-decadal timescales, variations in the North Atlantic Oscillation (NAO), the Atlantic Multidecadal Oscillation (AMO), the El Nino–Southern Oscillation (ENSO), and the Indian Ocean Dipole (IOD), can also modulate ISM behavior and strength23,30,48,49,50,51.
Deoria Tal is a small (surface area = 2.7 ha), closed basin lake, surrounded by a relatively small catchment (11 ha). Limnological data for the lake are available in the Supplemental Material (Table S1). Deoria Tal is moderately deep (depth = 16 m) (Fig. 1) and was strongly stratified when sampled on May 18, 2018. Surface water temperature was 19.1 °C, with a bottom water temperature of 6.4 °C (Table S1). A mat of water chestnut (genus: Trapa), which requires relatively shallow (1–4 m), nutrient rich water52, covers the littoral of Deoria Tal. Human impact on the lake system is currently relatively limited as it is protected by the Uttarakhand Forestry Department; however, human influence has been greater in the past. The expansion of the Indus civilization eastward through the Ganga Plain began at ~ 5000 cal yr BP53 and Demske et al.31 suggest, based on their pollen analysis of a lake sediment core from a nearby site, Badanital (located 20 km east of Deoria Tal), that migrants may have potentially arrived in the region at ~ 4000 cal yr BP. Other evidence of human influence includes the presence of the well-documented Painted Grey Ware (PGW) culture and large-scale cultivation between 2000 and ~ 1300 cal yr BP in the region31.
Results
Chronology and sedimentation rates
The age-depth model, developed in BACON54, indicates that the base of the sediment core extends back to ~ 5200 cal yr BP (Table 1; Fig. 2). An extant Trapa seed case returned a modern radiocarbon date with a percent Modern Carbon (pMC) equal to 101.9. This result indicates that the 14C dates obtained on the Trapa seed cases do not need to be corrected for a reservoir effect. A prominent increase in the rate of sedimentation occurs at ~ 1500 cal yr BP. The rate of sedimentation goes from ~ 17 years/cm at 1600 cal yr BP (264 cm) to 10 years/cm at 1500 cal yr BP (256 cm). The timing of this change in sedimentation, which aligns with the onset of zone DRL-3 (1500 cal yr BP; 256 cm) and a decrease in the flux of detrital elements, e.g. Rb, Sr,Ti and Zr, corresponds to the occurrence of large-scale cultivation in the region31.
Chemostratigraphy
The results of a Principal Components Analysis (PCA) of the pXRF-derived elemental data was used to divide the core into four distinct zones: DRL-1 (5200–3000 cal yr BP) is characterized by fluctuating and anti-phased PCA axes 1 and 2 scores; DRL-2 (3000–1500 cal yr BP) is characterized by decreasing PCA Axis 1 and increasing PCA Axis 2 scores; DRL-3 (1500–800 cal yr BP) is characterized by decreasing PCA axes scores; and DRL-4 (800 cal yr BP–present) is characterized by anti-phased PCA axes scores (Fig. S2). The pXRF provides returns for ten elements: silicon Kα, potassium Kα, calcium Kα, titanium Kα, manganese Kα, iron Kα, rubidium Kα, strontium Kα, zirconium Kα, and thorium Lα (from hereon, referenced by the respective elemental symbol) (Fig. S3). The Spearman’s correlation coefficients indicate the existence of a statistically significant correlation between all of the elements (Fig. S4). A PCA with the 10 measured elements and five ratios (Fe/Mn, Zr/Rb, Si/Ti, Ca/Ti and Mn/Ti) indicates that approximately 78% of the variance in the dataset could be explained by the first two principal component axes (Fig. S5). PCA Axis 1, which captures 61% of the variance in the elemental dataset, is strongly correlated with K, Ti, Fe, Rb, Sr, Zr, and Th (Table S2). PCA Axis 2, which captures an additional 16% of the variance, is strongly correlated with Mn, Fe/Mn and Mn/Ti. A plot of the PCA axes against time indicates notable deviations in the elemental data centered at 4850, 4200, and 3100 cal yr BP (Fig. 3).
DRL-1 (380–318 cm; 5200–3000 cal yr BP): Elemental returns in DRL-1 are largely driven by variations in the elements strongly correlated with PCA Axis 1 (K, Ti, Fe, Rb, Sr, Zr, and Th) while smaller magnitude changes are evident in Si, Ca, and Mn through the early portion of this zone (Figs. S2, S3). Prominent fluctuations in the elemental counts are centered at 4850, 4200, and 3100 cal yr BP with a less pronounced event evident at 3500 cal yr BP. The events centered at 4850, 4200 and 3500 cal yr BP are initially characterized by increasing PCA Axis 1 scores, HU number (measure of radiodensity calibrated according to the Hounsfield Scale), Th, and Fe/Mn and decreasing PCA Axis 2 scores, Mn/Ti, Ca/Ti, and Si/Ti between 5000 and 4850, 4350 and 4200 and 3600 and 3500 cal yr BP, respectively (Figs. 3, S6, S7). Each of these events is then followed by an interval characterized by decreasing PCA Axis 1 scores, HU number, Th, and Fe/Mn and increasing PCA Axis 2 scores, Mn/Ti, Ca/Ti, and Si/Ti from 4850 to 4700, 4200–4050 and 3500 to 3400 cal yr BP, respectively. The last event, between ~ 3200 and 3000 cal yr BP, is initially defined by decreasing PCA Axis 1 scores, HU number, Ti, and Fe/Mn and increasing PCA Axis 2 scores, Mn/Ti, Ca/Ti, and Si/Ti between ~ 3200 and 3100 cal yr BP; followed by increasing PCA Axis 1 scores, HU number, Th, and Fe/Mn and decreasing Mn/Ti, Ca/Ti, and Si/Ti between 3100 and 3000 cal yr BP.
DRL-2 (318–256 cm; 3000–1500 cal yr BP): Prominent fluctuations in the elemental counts are centered at 2000 and 1650 cal yr BP (Figs. S2, S3). The events between 2100 and 1900 cal yr BP and 1700 and 1600 cal yr BP are initially characterized by decreasing PCA Axis 1 scores, HU number, Th, and Fe/Mn with increasing PCA Axis 2 scores, Mn/Ti, Ca/Ti, and Si/Ti (Figs. S6, S7). These events are followed by an interval characterized by increasing PCA Axis 1 scores, HU number, Th, and Fe/Mn and decreasing PCA Axis 2 scores, Mn/Ti, Ca/Ti, and Si/Ti between 2000 and 1900 cal yr BP and 1650 and 1600 cal yr BP, respectively.
DRL-3 (256–169 cm; 1500–800 cal yr BP): The PCA results are more variable in DRL-3 relative to DRL-1 and DRL-2 with large magnitude shifts in the elemental returns centered on 1450 and 1250 cal yr BP (Figs. S2, S3). The event between 1500 and 1400 cal yr BP is initially characterized by lower PCA Axis 1 scores, HU number, Th, and Fe/Mn and increasing PCA Axis 2 scores, Mn/Ti, Ca/Ti, and Si/Ti between 1500 and 1450 cal yr BP; followed by increasing PCA Axis 1 scores, HU number, Th, and Fe/Mn and decreasing PCA Axis 2 scores, Mn/Ti, Ca/Ti, and Si/Ti between 1450 and 1400 cal yr BP (Figs. S6, S7). The event between 1350 and 1150 cal yr BP is initially characterized by lower PCA Axis 1 scores, HU number, Th, and Fe/Mn and increasing Ca/Ti and Si/Ti between 1350 and 1250 cal yr BP; followed by increasing PCA Axis 1 scores, HU number, Th, and Fe/Mn and decreasing Ca/Ti and Si/Ti between 1250 and 1150 cal yr BP.
DRL-4 (169–50 cm; 800–180 cal yr BP): This zone is characterized by the greatest amount of variance in the PCA and CT scan results. The PCA identifies multiple shifts in the elemental returns, the most prominent of which are centered at 500, 450, and 400 cal yr BP (Figs. S2, S3). These events are characterized by lower PCA Axis 1 scores, HU number, and Fe/Mn and elevated PCA Axis 2 scores, Mn/Ti, Ca/Ti, and Si/Ti. The interval between 400 and 180 cal yr BP is also characterized by highly variable elemental counts and sediment density, as inferred from HU number (Figs. S6, S7).
Discussion
Given our interest in assessing the relationship between hydroclimate variability and documented changes in the Indus civilization we focus our efforts on developing a holistic paleoclimate reconstruction spanning the interval captured by DRL-1 (5200–3000 cal yr BP), paying particular attention to the hydroclimate events evidenced at 4200 and 3100 cal yr BP.
Environmental interpretation of elemental and radiological proxies
Hydroclimatic regimes can influence the allogenic processes responsible for shaping the detrital flux to a lake. Typically, in subtropical environments during a wetter precipitation regime an elevated flux of coarse-grained sediment will be deposited within the lake57,58,59. As a means to distinguish anomalous climatic conditions in lacustrine records, general associations have been recognized between specific elements, as identified by XRF, and grain size at numerous sites in varying geologic and hydroclimatic conditions60. Existing grain size-element associations enable the development of specific elemental ratios that further aid environmental interpretation of hydroclimate variability. Fluctuations in Ti, which are associated with grain-size variations, have been used to reliably reconstruct rainfall and run-off61,62 and detrital input63,64 in varied settings. At Deoria Tal, PCA Axis 1 is inferred to primarily reflect variations in detrital influx. The elements most strongly correlated with PCA Axis 1 (K, Ti, Fe, Rb, Sr, Zr, and Th) have been identified has common constituents of detrital-sourced sediment14,58,65,66,67,68,69. The grouping of Ti with PCA Axis 1 and the strong Spearman’s correlation coefficients between Ti and Rb (ρ2 = 0.94), and Zr (ρ2 = 0.92), further supports the inference that variations in Ti are related to detrital input from the catchment70 (Figs. S4 and S5). Changes in sediment density, as captured by the HU number, are used to provide additional support for the PCA-based inferences of detrital flux71,72. A higher HU number is reflective of increasing sediment density and is inferred to occur during intervals characterized by elevated clastic/minerogenic input (i.e. increasing detrital flux). Decreasing HU numbers are inferred to reflecting lower sediment density and reduced detrital input.
PCA Axis 2 is inferred to reflect variations in hypolimnetic oxygen conditions. Hypolimnetic oxygen conditions are often reconstructed using iron28 and manganese54,62. The relative abundance of iron and manganese is influenced by redox potential in the benthic zone of aquatic systems53,54,57,58,64, which in turn is influenced by hypolimnetic oxygen conditions. The element most strongly correlated with PCA Axis 2 is Mn, which is strongly influenced by redox conditions66,73 (Fig. S5). During periods of stratification, ventilation of the bottom water is limited resulting in the reduction of hypolimnetic oxygen concentration. Higher lake levels reduce the ability of the wind to mix (ventilate) the entire water column of the lake and facilitates the onset of anaerobic conditions in the hypolimnion; whereas, lower lake levels would facilitate wind-mixing of the water column and increase oxygen supply to the bottom water. Fe/Mn and Mn/Ti ratios are commonly utilized as proxies for hypolimnetic oxygen conditions and redox conditions74,75. It is expected that Fe/Mn would decrease during periods characterized by abundant hypolimnetic oxygen concentration (lower lake level, weakened/no stratification) given that Mn is highly insoluble in oxic conditions and increasingly soluble, relative to Fe, in reducing (anoxic) environments76. It is important to note that redox potential in freshwater, is controlled by a suite of variables (e.g. dissolved oxygen, pH, temperature, organic matter, and nutrient availability), all of which can be influenced by changing limnological conditions including variations in lake productivity and lake level77.
In variable environments, such as lake systems, it is standard to normalize the variations of the pXRF-derived element counts for Ca and Si using Ti, an element that is not affected by biological or redox processes78. Ca/Ti and Si/Ti are used to decipher lake level and changes in autochthonous productivity, respectively (Fig. 3). Calcium in lake sediments may be related to both carbonate weathering in the catchment and in-lake precipitation of carbonate. Carbonate is precipitated when lake waters become saturated with respect to cCaCO3, which can occur when solutes become increasingly concentrated in response to a reduction in lake volume79. Numerous studies have documented that Ca/Ti can serve as a sensitive hydrological proxy with high Ca/Ti values indicative of lower lake level and evaporative enrichment and low Ca/Ti values reflecting increasing lake volume80,81. The PCA analysis indicates that Ca/Ti loads negatively on the axis capturing detrital flux (PCA Axis 1), providing support that the variations in Ca are primarily driven by within lake processes rather catchment processes (Fig. S5). This, together with the existence of a relatively weak correlation between Ca and detrital elements (i.e. Ti, Rb, Zr and Sr) (see Fig. S4) supports the claim that the Ca present in the lake sediment is related to the autochthonous precipitation of carbonate. Treatment of sediment samples with HCl for on-going stable isotope analysis indicates the presence of CaCO3 throughout the core. Si/Ti can be used as a proxy for biogenic silica with elevated Si values representing an increase in diatom productivity64,82,83 or as an indicator of clastic input75,84. The PCA analysis indicates that Si/Ti loads negatively on PCA Axis 1 (Fig. S5), suggesting that the variations in Si evidenced in the sediment are driven by autochthonous biogenic productivity rather than clastic input. We infer that elevated diatom productivity at Deoria Tal, as evidenced by high Si/Ti, occurred during intervals characterized by weakened stratification (due to the increasing effectiveness of wind mixing of a shallower lake), which increased nutrient availability in the surface waters. The anti-phased relationship between Si/Ti and Fe/Mn between 5000 and 3000 cal yr BP (Fig. 3) supports this inference with increased biogenic productivity (reflected by elevated Si/Ti) occurring during intervals characterized by lower lake levels, weakened stratification and abundant hypolimnetic oxygen concentration (reflected by reduced Fe/Mn)85.
Paleoclimate reconstruction
The event between 4350 and 4050 cal yr BP is initially characterized by elevated precipitation and a positive hydroclimate anomaly (Fig. 3). This event is larger in magnitude than the event centered at 4850 cal yr BP. This inference is supported by the elevated Fe/Mn ratio suggesting the occurrence of anoxic conditions in the bottom water between 4350 and 4200 cal yr BP. The existence of anoxia in the deep water, also supported by lower Mn/Ti, is likely a response to a deepening of the lake and strengthened thermal stratification, both of which would reduce water column mixing and thereby reduce oxygenation of the bottom water. As with the 4850 cal yr BP event, reduced diatom productivity, reflecting limited nutrient availability, is also evident. The increase in allogenic detrital input between 4350 and 4200 cal yr BP, inferred from the elevated counts for PCA Axis 1, Th and HU number provide further support for this inference, as increased precipitation would increase erosion in the catchment and the transport of detrital material to the lake.
Evidence of wet conditions in northern India is documented at Sahiya Cave, with variations in δ18O inferred to reflect an interval of elevated ISM between at 4300 and 4150 cal yr BP4 (Fig. 4). An abrupt increase in wind mixing and oxygenation of lake water between 4400 and 4200 cal yr BP at Lake Rara, in western Nepal, is associated with a strengthened ISM during this interval14 (Fig. 4). Marine records from the Arabian Sea suggest that monsoonal precipitation peaked between 4500 and 4300 cal yr BP18, with evidence that the increase in precipitation during this interval was driven primarily by an increase in the IWM11 (Fig. 4). At 4200 cal yr BP the parameters outlined above shift in the opposite direction at DRL suggesting the onset of drier conditions, leading to a negative lake mass balance and a lowering of lake level. An abrupt decrease in effective moisture between 4200 and 4000 cal yr BP is documented in multiple records from South Asia. Peak drying is evident at 4100 cal yr BP at paleolake Kotla Dahar, in northwest India24; between 4200 and 4000 cal yr BP at Lonar Lake, in central India25; and at 4200 cal yr BP at Tso Kar Lake, in northern India28. A shift from pluvial conditions to a variable ISM characterized by multiple multi-decadal scale droughts is inferred to have occurred between 4200 and 4050 cal yr BP in northern India86. The onset of dry conditions at this time is also documented at Lake Rara in western Nepal14 and at Mawmluh Cave in northeast India34. Marine records indicate a brief interval of increasing surface water salinity in the Bay of Bengal at 4200 cal yr BP32; the occurrence of a centennial-scale reduction in the discharge of the Indus River at ~ 4200 cal yr BP18; and a surface water salinity episode in the northern Arabian Sea, inferred to reflect a concurrent reduction in the IWM and ISM beginning at 4100 cal yr BP11.
The event between 3200 and 3000 cal yr BP is inferred to initially reflect the occurrence of a severe drought episode followed by a notable increase in effective moisture at Deoria Tal. This inference is supported by the existence of low Fe/Mn and high Mn/Ti between 3200 and 3100 cal yr BP, which suggest that the bottom waters were well oxygenated, likely resulting from the increased effectiveness of wind-driven water column mixing in a shallow lake. The negative excursion in Ti provides support for lake drawdown during this interval as well, as lower counts of Ti have been associated with drought conditions in subtropical environments57. The elevated Ca/Ti values, which peak at 3100 cal yr BP, provide additional support for the inference of lower lake levels at this time. An increase in Ca in lake sediment is expected to occur when lake water becomes saturated with respect to calcium, which can occur with reduced lake volume79. The increase in diatom productivity, resulting from an increase in nutrient availability due to wind mixing of the shallow lake, provides further support for the inferred reduction of lake levels between 3200 and 3100 cal yr BP. Together with the PCA Axis 1 data, these lines of evidence indicate that this event was characterized by a hydroclimate anomaly that was large enough to significantly lower lake levels. The presence of visibly coarser materials along with an abrupt increase in the PCA Axis 1 scores and HU number at 3100 cal yr BP, suggest that this drought episode was terminated abruptly at 3100 cal yr BP and was followed by a positive hydroclimate anomaly and increasing lake levels.
Evidence of a transition to increasingly arid conditions at ~ 3200 cal yr BP is evident in other regional records with a decadal-scale interval of lowered effective moisture29 leading to forest thinning44 and marked changes in vegetation due to fire and human activity31. This drought episode and the subsequent increase in aridity is very well expressed in a speleothem record from nearby Sainji Cave13 (Fig. 4) and a peat sequence from Kedarnath12 (Fig. 4). Nakamura et al.14 identify a brief decline in the ISM occurring at 3100 cal yr BP (Fig. 4), signaling a transition to a drier regime characterized by a weakened monsoon driven by a reduction in solar insolation87. In coastal areas, the dry pulse at 3200 cal yr BP favored the expansion of salt tolerant mangroves88. At approximately 3000 cal yr BP this arid interval was terminated17,44,89. Following this transition, our record supports the inference of a trend towards wetter conditions, as evidenced in other regional records89,90, between 3000 and ~ 2500 cal yr BP. Importantly, the episode of late Holocene peak aridity captured by the Sainji Cave speleothem at ~ 3150 cal yr BP and our record at Deoria Tal coincides with the timing of societal collapse and the disappearance of Indus culture in northwest India15.
Hydroclimate variability and Indus civilization re-organization at 4200 cal yr BP
Broadly, DRL-1 (5200–3000 cal yr BP) lies within the middle of a climatic regime characterized by a weakening ISM4,68. The trend of decreasing precipitation though the interval captured by DRL-1 was not gradual nor unidirectional. Notable hydroclimate anomalies characterized this interval with distinctive, discrete multi-decadal to centennial-scale episodes of fluctuating hydroclimate evident at Deoria Tal. The interval between 5000 and 3000 cal yr BP is of particular interest because it spans the interval that captures the rise, expansion, and eventual contraction of the Indus civilization8. The Indus civilization, which reached its apogee between 4300 and 3900 cal yr BP, consisted of several thousand communities along the Indus River and the Ghaggar-Hakra River by 4200 cal yr BP91,92. The Indus civilization’s agriculture system was dependent on flooding rather than precipitation directly19. The rise of the cultural complex associated with the Indus civilization, between 5200 and 4600 cal yr BP, occurs during an interval characterized by a trend of decreasing intensity of ISM and IWM precipitation. This “Goldilocks” scenario, which lasted from approximately 5200 to 4200 cal yr BP, allowed for crops to be planted in areas inundated frequently enough by floods of sufficient magnitude to provide needed irrigation but not so frequently or by such large magnitude floods that crops were damaged20.
According to Staubwasser et al.18, the 4200 cal yr BP event impacted the agricultural base of this highly centralized agrarian society, which likely contributed to its geographic and demographic contraction and may have ultimately led to demise of their civilization. Modeling and paleoenvironmental studies suggest that the rapid, step-wise, and sustained decline in population that characterized the Indus civilization, which began at ~ 4200 cal yr BP, can be attributed to increasing water scarcity associated with changes in regional hydroclimate and a weakening of the ISM13,18,23,24,29. This hypothesis is supported by meteorological observations (1961–1990), which indicate that ~ 80% of total mean annual precipitation that falls in the ISM dominated region falls between June and August23. However, recent research has also highlighted the potential contribution of winter precipitation to the overall precipitation regime in this region during the mid- to late Holocene11,20. The interval between 4500 and 4100 cal yr BP was characterized by fluctuations in the strength of the IWM11,20,28,88, with evidence of a strengthened IWM between 4500 and 4300 cal yr BP and a weakened IWM from 4300 to 4100 cal yr BP, followed by a step- shift at 4100 cal yr BP. Our record supports the inference of positive hydroclimate anomaly peaking at 4200 cal yr BP in the Garhwal Himalaya and the onset of lake drawdown and a severe centennial-scale drought beginning at 4200 cal yr BP. This study cannot distinguish whether the precipitation anomalies identified at Deoria Tal during this interval were generated by variations in the ISM and/or the IWM; however, it is notable that the positive hydroclimate anomalies evidenced in this record correspond to intervals of strengthened IWM11.
Conclusion
The multi-proxy record generated from Deoria Tal identifies major hydroclimate anomalies in the Garhwal Himalaya at 4850, 4200, and 3100 cal year BP, with the last two events being particularly noteworthy. A decrease in allogenic detrital input and Fe/Mn and an increase in Mn/Ti, Ca/Ti and Si/Ti is inferred to reflect the onset of drier conditions at 4200 cal yr BP. Shifts in Fe/Mn and Mn/Ti; a negative excursion in the detrital flux; and elevated Ca/Ti provides support for the inference of lower lake levels and the onset of a centennial-scale drought episode between 3200 and 3100 cal yr BP. The drought events documented at Deoria Tal correspond to existing records from the northern portion of the Indian subcontinent, providing additional insight into the spatial expression of mid- to late Holocene hydroclimate anomalies in this highly populated region. This study sheds light on the temporal relationship between hydroclimate anomalies and large scale deurbanization, population decrease and the spatial contraction of Indus civilization in northern India between 4500 and 3000 cal yr BP. This record suggests that the onset of de-urbanization and the deterioration of the Indus civilization at ~ 4200 cal yr BP coincided with largest hydroclimate anomaly evidenced in the Garhwal Himalaya during the last five millennia.
Materials and methods
Core recovery and sediment sampling
A sediment core was recovered from Deoria Tal in May 2018 using a modified Livingstone corer deployed from a platform anchored at a depth of 5.65 m (see Fig. 1 for coring location). The core consists of five overlapping drives with the total composite core length equal to 380 cm. The overlapping core drives were correlated using visual stratigraphy (Fig. S8), the results from the pXRF analysis and radiological imagery (Fig. S9). The uppermost sediment was recovered using a plastic tube to ensure that the flocculent surface sediment was captured with minimal disturbance while the stiffer lower sediment was recovered using a stainless-steel barrel. The flocculent, watery nature of the uppermost sediment necessitated that this sediment was sub-sampled in the field (at 0.5 cm intervals) at the time of core recovery to prevent mixing (these samples were not used for radiological analyses). Limnological measurements were made and water samples for water chemistry were also collected at this time (Table S1). Stratigraphical notes were made in the field and later at Kumaun University, Nainital, Uttarakhand, where each drive was longitudinally split, photographed, and described.
Chronology development
Chronological control for the sediment core recovered from Deoria Tal is based on ten AMS radiocarbon (14C) dates obtained on Trapa seed cases (Table 1). Additionally, a seed case from Trapa extant in Deoria Tal (marked as “modern” in Table 1) was also dated to determine if ‘old’ carbon or a reservoir effect was present. The samples used for radiocarbon dating were rinsed in deionized water, bagged, and labeled before being submitted for analysis at the Center for Applied Isotope Studies (CAIS) at the University of Georgia. The radiocarbon dates were converted to calendar years using OxCal 4.355 with the reported 2σ age ranges following Reimer et al.56. An age-depth model was developed using BACON54 (Bayesian Accumulation Model, version 2.2), an open-source R code package (Fig. 2).
Laboratory analyses
Preparation of the cores for elemental analysis using X-ray fluorescence (XRF) was completed at the Environmental Change Lab at the University of Georgia. The sediment cores were scanned using a Bruker Tracer 5i mounted portable XRF (pXRF) with a rhodium tube source at CAIS. A scanning resolution of 2 mm/step was utilized for the pXRF analysis. To detect elements with a low atomic number (Z) a 60 s exposure time of 20 keV, 35 amps, and no filter was employed. To detect high Z elements, a 60 s exposure time of 40 keV, 35 amps, and a Cu 100 µm: Ti 25 µm: Al 300 µm filter was employed93. The scan counts, based on a single scan, were initially post-processed using Bruker Artax software to determine the specific high and low Z elements captured by the analysis. All scan count data are presented. The heterogenous composition of lake sediment and the nature of the sediment matrix can affect pXRF readings by differentially influencing the resultant emission spectra. The sediment matrix effects were corrected following standard procedures78. Elemental data from the pXRF analysis were normalized using the rhodium Kα peak to minimize the noise variation of the Compton peak. To distinguish natural variation from anthropogenic influence the resultant time-series spectra were normalized using the aluminum Kα record14,94. No error estimates are available for the pXRF analysis; however, the measured response rates (Kα) of the elements range from 15,878 to 247 counts/second, which are well above the nominal sensitivity of the Bruker Tracer 5i95. Principal component analysis (PCA), of the XRF-derived elemental data (non-transformed, centered and standardized) was utilized to identify which of the elements could account for a statistically significant amount of variance along a reduced number of environmental axes65,96,97,98 (Fig. S5; Table S2). The statistical analyses were undertaken in R99 and CANOCO version 4.0100.
Computerized Tomography (CT) scans, also known as Computerized Coaxial Tomography (CAT) scans71, were completed at the University of Georgia Veterinary Teaching Hospital’s Radiology and Imaging Department on a Siemens Sensation 64 CT scanner utilizing Syngo WinNT 5.1 Service pack 3 imaging processing software. The output from CT scans, known as CT numbers, are calibrated according to the Hounsfield Scale. Hounsfield units (HU) are a dimensionless unit universally used to express CT numbers in a standardized form. HU are a measure of radiodensity based on the linear transformation of the measured attenuation coefficients which are calibrated using the arbitrarily-assigned densities of air (0 HU) and pure water (− 1000 HU). This calibration results in a scale ranging from − 1000 HU for air to + 3000 for metal101. HU numbers (or CT numbers) are influenced by bulk sediment density, mineralogy, and porosity71,102 and therefore CT scans can be used to identify internal structure within sediment and aid in paleoenvironmental interpretation72,103. In this study, an increase in HU number is interpreted to reflect an increase in sediment density resulting from higher clastic/minerogenic input (i.e. elevated detrital flux); a decrease in HU number is interpreted to reflect lower sediment density resulting from a reduction in clastic/minerogenic input (i.e. lower detrital flux). Each core drive was imaged axially with two different exposure settings: Soft (40 HU window level; 300 HU window width) and Bone (400 HU window level; 4000 HU window width) (Fig. S9). The tomograms were compiled using a MATLAB™ package104. Presentation of the CT scan HU data (z-scores and standard deviation) is limited to the Soft CT settings output (Fig. S6).
References
Weiss, H. Global megadroughts, societal collapse and resilience at 4.2–3.9 ka BP across the Mediterranean and west Asia. Pages Mag. 24, 62–63 (2016).
Wolff, C. et al. Precipitation evolution of Central Asia during the last 5000 years. The Holocene 27(1), 142–154 (2017).
Booth, R. K. et al. A severe centennial-scale drought in midcontinental North America 4200 years ago and apparent global linkages. The Holocene 15(3), 321–328 (2005).
Kathayat, G. et al. The Indian monsoon variability and civilization changes in the Indian subcontinent. Sci. Adv. 3(12), e1701296 (2017).
Perşoiu, A., Ionita, M. & Weiss, H. Atmospheric blocking induced by the strengthened Siberian High led to drying in west Asia during the 4200ka BP event–a hypothesis. Climate Past 15(2), 781–793 (2019).
Park, J. et al. Abrupt Holocene climate shifts in coastal East Asia, including the 8.2 ka, 4200ka, and 2.8 ka BP events, and societal responses on the Korean peninsula. Sci. Rep. 9(1), 1–16 (2019).
Ran, M. & Chen, L. The 4.2 ka BP climatic event and its cultural responses. Quat. Int. 521, 158–167 (2019).
Dutt, S. et al. Climate variability and evolution of the Indus civilization. Quat. Int. 507, 15–23 (2019).
Wheeler, M. The Indus Civilization (Cambridge University Press Archive, 1968).
Kenoyer, J. M. The Indus valley tradition of Pakistan and western India. J. World Prehist. 5(4), 331–385 (1991).
Giesche, A., Staubwasser, M., Petrie, C. A. & Hodell, D. A. Indian winter and summer monsoon strength over the 4.2 ka BP event in foraminifer isotope records from the Indus River delta in the Arabian Sea. Climate Past 15, 73–90 (2019).
Srivastava, P. et al. 8000-year monsoonal record from Himalaya revealing reinforcement of tropical and global climate systems since mid-Holocene. Sci. Rep. 7(1), 14515 (2017).
Kotlia, B. S., Singh, A. K., Joshi, L. M. & Dhaila, B. S. Precipitation variability in the Indian Central Himalaya during last ca. 4,000 years inferred from a speleothem record: impact of Indian Summer Monsoon (ISM) and Westerlies. Quat. Int. 371, 244–253 (2015).
Nakamura, A. et al. Weak monsoon event at 42 ka recorded in sediment from Lake Rara, Himalayas. Quat. Int. 397, 349–359 (2016).
Madella, M. & Fuller, D. Q. Palaeoecology and the Harappan Civilisation of South Asia: a reconsideration. Quatern. Sci. Rev. 25, 1283–1301 (2006).
Schug, G. R., Belvins, K. E., Cox, B., Gray, K. & Mushrif-Tripathy, V. Infection, disease, and biosocial processes at the end of the Indus civilization. PLoS ONE 8(12), e84814 (2013).
Wright, R. P., Bryson, R. A. & Schuldenrein, J. Water supply and history: Harappa and the Beas regional survey. Antiquity 82, 37–48 (2008).
Staubwasser, M., Sirocko, F., Grootes, P. M. & Segl, M. Climate change at the 4.2 ka BP termination of the Indus valley civilization and Holocene south Asian monsoon variability. Geophys. Res. Lett. 30(8), 1425 (2003).
Giosan, L. et al. Fluvial landscapes of the Harappan civilization. Proc. Natl. Acad. Sci. 109(26), E1688-1694 (2012).
Giosan, L. et al. Neoglacial climate anomalies and the Harappan metamorphosis. Climate Past. 14, 1669–1686 (2018).
Brooks, N. Cultural responses to aridity in the Middle Holocene and increased social complexity. Quat. Int. 151(1), 29–49 (2006).
Staubwasser, M. & Weiss, H. Holocene climate and cultural evolution in late prehistoric–early historic West Asia. Quat. Res. 66(3), 372–387 (2006).
MacDonald, G. Potential influence of the Pacific Ocean on the Indian summer monsoon and Harappan decline. Quat. Int. 229, 140–148 (2011).
Dixit, Y., Hodell, D. A. & Petrie, C. A. Abrupt weakening of the summer monsoon in northwest India ~4100 yr ago. Geology 42(4), 339–342 (2014).
Prasad, S., Kusumgar, S. & Gupta, S. K. A mid to late Holocene record of palaeoclimatic changes from Nal Sarovar: a palaeodesert margin lake in western India. J. Quat. Sci. 12(2), 153–159 (1997).
Sharma, C. & Gupta, A. Vegetation and climate in Garhwal Himalaya during EarlyHolocene: Deoria Tal. The Palaeobotanist. 46(3), 111–116 (1997).
Enzel, Y. et al. High-resolution holocene environmental changes in the Thar Desert, Northwestern India. Science 284(5411), 125–128 (1999).
Wünnemann, B. et al. Hydrological evolution during the last 15 kyr in the Tso Kar lake basin (Ladakh, India), derived from geomorphological, sedimentological and palynological records. Quat. Sci. Rev. 29, 1138–1155 (2010).
Leipe, C., Demske, D., Tarasov, P. E., HIMPAC Project Members. A Holocene pollen record from the northwestern Himalayan lake Tso Moriri: Implications for palaeoclimatic and archaeological research. Quat. Int. 348, 93–112 (2014).
Menzel, P. et al. Linking Holocene drying trends from Lonar Lake in monsoonal central India to North Atlantic cooling events. Palaeogeogr. Palaeoclimatol. Palaeoecol. 410, 164–178 (2014).
Demske, D. et al. Record of vegetation, climate change, human impact and retting of hemp in Garhwal Himalaya (India) during the past 4600 years. The Holocene 26(10), 1661–1675 (2016).
Ponton, C. et al. Holocene aridification in India. Geophys. Res. Lett. 39, L03704 (2012).
Bolton, C. T. et al. A 500,000 year record of Indian summer monsoon dynamics recorded by eastern equatorial Indian Ocean upper water-column structure. Quatern. Sci. Rev. 77, 167–180 (2013).
Berkelhammer, M. et al. An abrupt shift in the Indian monsoon 4000 years ago. Geophys. Monogr. Ser. 198, 75–87 (2012).
Rajani, M. B. & Rajawat, A. S. Potential of satellite based sensors for studying distribution of archaeological sites along palaeo channels: Harappan sites a case study. J. Archaeol. Sci. 38, 2010–2016 (2011).
Chatterjee, A. & Ray, J. S. Sources and depositional pathways of mid-Holocene sediments in the Great Rann of Kachchh, India: implications for fluvial scenario during the Harappan Culture. Quatern. Int. 443, 177–187 (2017).
Bist, K. S. & Sinha, A. K. Some observations on the geological and structural setup of Okhimath area in Garhwal Himalaya. Himalayan Geol. 10, 467–475 (1980).
Chaudhary, S., Gupta, V. & Sundriyal, Y. P. Surface and sub-surface characterization of Byung landslide in Mandakini valley, Garhwal Himalaya. Himalayan Geol. 31(2), 125–132 (2010).
Valdiya, K. S. Geology of Kumaun Lesser Himalaya. Dehradun (Wadia Institute of Himalayan Geology, 1980).
Joshi, L. M. & Kotlia, B. S. Neotectonically triggered instability around the palaeolake regime in Central Kumaun Himalaya, India. Quat. Int. 371, 219–231 (2015).
Kotlia, B. S., Hinz-Schallreuter, I., Schallreuter, R. & Schwarz, J. Evolution of Lamayuru palaeolake in the Trans Himalaya: palaeoecological implications. E&G Quat. Sci. J. 48(1), 177–191 (1998).
Goswami, B. N. & Chakravorty, S. Dynamics of the Indian summer monsoon climate in climate systems and climate dynamics. Oxf. Res. Encycl. Climate Sci. https://doi.org/10.1093/acrefore/9780190228620.013.613 (2017).
Dixit, Y. & Tandon, S. K. Hydroclimatic variability on the Indian subcontinent in the past millennium; review and assessment. Earth Sci. Rev. 161, 1–15 (2016).
Sharma, C., Chauhan, M. S. & Rajagopalan, G. Vegetation and climate in Garhwal Himalaya during last 4,000 years. The Palaeobotanist. 49(3), 501–507 (2000).
Patwardhan, S., Kulkarni, A. & Sabade, S. Projected changes in semi-permanent system of indian summer monsoon in CORDEX-SA framework. Am. J. Clim. Chang. 5, 133–146 (2016).
Dimri, A. P., Yasunari, T., Kotlia, B. S., Mohanty, U. C. & Sikka, D. R. Indian winter monsoon: present and past. Earth Sci. Rev. 163, 297–322 (2016).
Kotlia, B. S. et al. Stalagmite based high resolution precipitation variability for past four centuries in the Indian Central Himalaya: chulerasim cave re-visited and data re-interpretation. Quatern. Int. 444, 35–43 (2017).
Saji, N. H., Goswami, B. N., Vinayachandran, P. N. & Yamagata, T. A dipole mode in the tropical Indian Ocean. Nature 401(6751), 360–363 (1999).
Wang, Y. et al. The Holocene Asian monsoon: links to solar changes and North Atlantic climate. Science 308(5723), 854–857 (2005).
Ihara, C., Kushnir, Y., Cane, M. A. & Victor, H. Indian summer monsoon rainfall and its link with ENSO and Indian Ocean climate indices. Int. J. Climatol. 27(2), 179–187 (2007).
Rao, S. A., Chaudhari, H. S., Pokhrel, S. & Goswami, B. N. Unusual central Indian drought of summer monsoon 2008: role of southern tropical Indian Ocean warming. J. Clim. 23(19), 5163–5174 (2010).
Hummel, M. & Kiviat, E. Review of world literature on water chestnut with implications for management in North America. J. Aquat. Plant Manag. 42, 17–28 (2004).
Gupta, A. K., Anderson, D. M., Pandey, D. N. & Singhvi, A. K. Adaptation and human migration, and evidence of agriculture coincident with changes in the Indian summer monsoon during the Holocene. Curr. Sci. 90(8), 1082–1090 (2006).
Blaauw, M. & Christen, J. A. Flexible paleoclimate age-depth models using an autoregressive gamma process. Bayesian Anal. 6, 457–474 (2011).
Bronk Ramsey, C. OxCal Program, Version 4.3. Oxford Radiocarbon Accelerator Unit https://c14.arch.ox.ac.uk/oxcal.html (University of Oxford, 2017).
Reimer, P. J. et al. IntCal13 and Marine13 radiocarbon age calibration curves 0–50,000 years cal BP. Radiocarbon 55(4), 1869–1887 (2013).
Haug, G. H. et al. Climate and the collapse of Maya civilization. Science 5613(299), 1731–7135 (2003).
Cuven, S., Francus, P. & Lamoureux, S. Mid to Late Holocene hydroclimatic and geochemical records from the varved sediments of East Lake, Cape Bounty, Canadian High Arctic. Quat. Sci. Rev. 30, 2651–2665 (2011).
Bohra, A. & Kotlia, B. S. Tectono-climatic signatures during late quaternary in the Yunam basin, Baralacha Pass (upper Lahaul valley, India), derived from multi-proxy records. Quatern. Int. 371, 111–121 (2015).
Davies, S. J., Lamb, H. F. & Roberts, S. J. Micro-XRF Core Scanning in Palaeolimnology: Recent Developments in Micro-XRF Studies of Sediment Cores (eds. Croudace, I. W. & Rothwell, R. G.) 17, 189–226 (Springer, 2015).
Metcalfe, S. E., Jones, M. D., Davies, S. J., Noren, A. & MacKenzie, A. Climate variability over the last two millennia in the North American Monsoon region, recorded in laminated lake sediments from Laguna de Juanacatlán, Mexico. The Holocene 20(8), 1195–1206 (2010).
Corella, J. P. et al. The 1.5-ka varved record of lake Montcortès (southern Pyrenees, NE Spain). Quat. Res. 78(2), 323–332 (2012).
Balascio, N. L. et al. Investigating the use of scanning X-ray fluorescence to locate cryptotephra in minerogenic lacustrine sediment: experimental results in Micro-XRF Studies of Sediment Cores (eds. Croudace, I. W. & Rothwell, R. G.) 17, 305–324 (Springer, 2015).
Stansell, N. D., Rodbell, D. T., Abbott, M. B. & Mark, B. G. Proglacial lake sediment records of Holocene climate change in the western Cordillera of Peru. Quatern. Sci. Rev. 70, 1–14 (2013).
Chassiot, L. et al. A 7000-year environmental history and soil erosion record inferred from the deep sediments of Lake Pavin (Massif Central, France). Palaeogeogr. Palaeoclimatol. Palaeoecol. 497, 218–233 (2018).
Heymann, C. et al. Late Glacial to mid-Holocene palaeoclimate development of Southern Greece inferred from the sediment sequence of Lake Stymphalia (NE-Peloponnese). Quatern. Int. 302, 42–60 (2013).
Olsen, J., Anderson, J. N. & Leng, M. J. Limnological controls on stable isotope records of late-Holocene palaeoenvironment change in SW Greenland: a paired lake study. Quat. Sci. Rev. 66, 85–95 (2013).
Peng, J. et al. Indian summer monsoon variations and competing influences between hemispheres since ~35 ka recorded in Tengchongqinghai Lake, southwestern China. Palaeogeogr. Palaeoclimateol. Palaeoecol. 516, 113–125 (2019).
Schmidt, M. et al. A multi-proxy palaeolimnological record of the last 16,600 years from coastal Lake Kushu in northern Japan. Palaeogeogr. Palaeoclimatol. Palaeoecol. 514, 613–626 (2019).
Koinig, K. A., Shotyk, W., Lotter, A. F., Ohlendorf, C. & Sturm, M. 9000 years of geochemical evolution of lithogenic major and trace elements in the sediment of an alpine lake—the role of climate, vegetation, and land-use history. J. Paleolimnol. 30, 307–320 (2003).
Francus, P., Kanamaru, K., & Fortin, D. Standardization and Calibration of X-Radiographs Acquared with the ITRA Core Scanner in Micro-XRF Studies of Sediment Cores. (eds. Croudace, I. W. & Rothwell, R. G.) 17, 491–505 (Springer, 2015).
Flisch, A. & Becker, A. Industrial X-ray computed tomography studies of lake sediment drill cores in Applications of X-ray Computed Tomography in the. Geosciences 215, 205–212 (2003).
Haenssler, E., Nadeau, M. J., Vött, A. & Unkel, I. Natural and human induced environmental changes preserved in a Holocene sediment sequence from the Etliko Lagoon, Greece: new evidence from geochemical proxies. Quatern. Int. 308–309, 89–104 (2013).
Moreno, A. et al. A 14 kyr record of the tropical Andes: the Lago Chungará sequence (18 S, northern Chilean Altiplano). Quatern. Int. 161(1), 4–21 (2007).
Kylander, M. E., Ampel, L., Wohlfarth, B. & Veres, D. High-resolution X-ray fluorescence core scanning analysis of Les Echets (France) sedimentary sequence: new insights from chemical proxies. J. Quat. Sci. 26(1), 109–117 (2011).
Davison, W. Iron and manganese in lakes. Earth Sci. Rev. 34(2), 119–163 (1993).
Hamilton-Taylor, J. & Davison, W. Redox-Driven Cycling of Trace Elements in Lakes in Physics and Chemistry of Lakes 217–263 (Springer, 1995).
Löwemark, L. et al. Normalizing XRF-scanner data: a cautionary note on the interpretation of high-resolution records from organic-rich lakes. J. Asian Earth Sci. 40(6), 1250–1256 (2011).
Cohen, A.S. Paleolimnology: The History and Evolution of Lake Systems (Oxford University Press, 2003).
Haberzettl, T. et al. Hydrological variability in southeastern Patagonia and explosive volcanic activity in the southern Andean Cordillera during Oxygen Isotope Stage 3 and the Holocene inferred from lake sediments of Laguna Potrok Aike, Argentina. Palaeogeogr. Palaeoclimatol. Palaeoecol. 259(2–3), 213–229 (2008).
Guillaume, J. et al. Microsedimentological characterization using image analysis and μ-XRF as indicators of sedimentary processes and climate changes during Lateglacial at Laguna Potrok Aike, Santa Cruz, Argentina. Quat. Sci. Rev. 71, 191–204 (2013).
Liu, X., Colman, S. M., Brown, E. T., Minor, E. C. & Li, H. Estimation of carbonate, total organic carbon, and biogenic silica content by FTIR and XRF techniques in lacustrine sediments. J. Paleolimnol. 50, 387–398 (2013).
Brown, E. T., Johnson, T. C., Scholz, C. A., Cohen, A. S. & King, J. W. Abrupt change in tropical African climate linked to the bipolar seesaw over the past 55,000 years. Geophys. Res. Lett.. 34(20), 1–5 (2007).
Brown, E. T. Estimation of biogenic silica concentrations using scanning XRF: insights from studies of Lake Malawi sediments in Micro-XRF Studies of Sediment Cores (eds. Croudace, I. W. & Rothwell, R. G.) 17, 267–277 (Springer, 2015).
McGlue, M. M. et al. Solar irradiance and ENSO affect food security in Lake Tanganyika, a major African inland fishery. Sci. Adv. 6(41), p.eabb2191 (2020).
Kathayat, G. et al. Evaluating the timing and structure of the 4.2 ka event in the Indian summer monsoon domain from an annually resolved speleothem record from Northeast India. Climate Past. 14(12), 1869–1879 (2018).
Ankit, Y., Kumar, P., Anoop, A., Mishra, P. K. & Varghese, S. Mid-late Holocene climate variability in the Indian monsoon: evidence from continental shelf sediments adjacent to Rushikulya river, eastern India. Quatern. Int. 443, 155–163 (2017).
Kar, R. & Quamar, M. F. Pollen-based quaternary palaeoclimatic studies in India: an overview of recent advances. Palynology 43(1), 76–93 (2019).
Rawat, S., Gupta, A. K., Sangode, S. J., Srivastava, P. & Nainwal, H. C. Late Pleistocene-Holocene vegetation and Indian summer monsoon record from the Lahaul, northwest Himalaya, India. Quat. Sci. Rev. 114, 167–181 (2015).
Kotlia, B. S. & Joshi, L. M. Late Holocene climatic changes in Garhwal Himalaya. Curr. Sci. 104(7), 911–919 (2013).
Pant, G. B. Long-term climate variability and change over monsoon Asia. J. Indian Geophys. Union. 7(3), 125–134 (2003).
Puri, V. M. K. Vedic Saraswati: Scientific signatures on its origin from Himalaya in Vedic River Sarasvati and Hindu civilization (ed. Kalyanaraman, S.) 14–15 (Aryan Books International, 2008).
Hunt, A. M. W. & Speakman, R. J. Portable XRF analysis of archaeological sediments and ceramics. J. Archaeol. Sci. 53, 626–638 (2015).
Elder, J. F. Metal Biogeochemistry in Surface-Water Systems—A Review of Principles and Concepts. Department of the Interior, United States Geological Survey. 1013 (1988).
Bruker Corporation. Application Note #1901: TRACER 5 Sensitivity. (2021).
Olsen, J. et al. Lacustrine evidence of Holocene environmental change from three Faroese lakes: a multiproxy XRF and stable isotope study. Quatern. Sci. Rev. 29, 2764–2780 (2010).
Hahn, A. et al. Elemental composition of the Laguna Potrok Aike sediment sequence reveals paleoclimatic changes over the past 51 ka in southern Patagonia, Argentina. J. Paleolimnol. 52, 349–366 (2014).
Hendy, I. L., Napier, T. J. & Schimmelmann, A. From extreme rainfall to drought: 250 years of annually resolved sediment deposition in Santa Barbara Basin, California. Quat. Int. 387, 3–12 (2015).
R Core Team. R: A language and environment for statistical computing. Vienna, Austria. https://www.r-project.org/ (R Foundation for Statistical Computing, 2020).
Ter Braak, C. J. & Smilauer, P. CANOCO reference manual and CanoDraw for Windows user's guide: Software for canonical community ordination (version 4.0). https://www.canoco.com (2002).
Seeram, E. Computed Tomography: Physical Principles, Clinical Applications, and Quality Control (Elsevier Health Sciences, 2015).
St-Onge, G., Mulder, T., Francus, P. & Long, B. Continuous Physical Properties of Cored Marine Sediments in Proxies in Late Cenozic Paleoceanography (eds. Hillaire-Marcel, C. & De Vernal, A.) 1, 63–98 (Elsevier Science and Technology, 2007).
Fouinat, L. et al. A new CT scan methodology to characterize a small aggregation gravel clast contained in a soft sediment matrix. Earth Surf. Dyn. 5, 199–209 (2017).
Reilly, B. T., Stoner, J. S. & Wiest, J. SedCT: MATLAB™ tools for standardized and quantitative processing of sediment core computed tomography (CT) data collected using a medical CT scanner. Geochem. Geophys. Geosyst. 18(8), 3231–3240 (2017).
Acknowledgements
A National Science Foundation Human-Environment and Geographical Sciences Award (GSS-2026311) to D.F.P and B.S.K. and a U.S. Fulbright-Nehru Fellowship and Global Collaborative Research Grant from the University of Georgia (UGA) to D.F.P. supported this research. The authors thank the following individuals from Kumaun University: Anoop Singh, Harish Bisht, Shahdita Bakshi and Harish Singh for fieldwork assistance; Neha Kholia for assisting in the description of the core stratigraphy; and Kimani Bisht and Harish Bisht for providing Figure 1A,B, respectively. The authors recognize the following individuals from UGA: Hunter Kunzelmann and Andrew Knazek for laboratory assistance; Jeff Speakman from the Center of Applied Isotope Studies (CAIS) for assistance with pXRF analyses; Alex Cherkinsky and Carla Hadden from CAIS for assistance with 14C dating; Billie Rae King, Dana Duncan, Charleen Romaine, and Dr. Scott Secrest from the Radiology and Imaging Department at UGA’s Veterinary Teaching Hospital for the acquisition of the radiological data. Dr. Brendan Reilly from Oregon State University is thanked for his guidance using the SedCT software.
Funding
The National Science Foundation USA.
Author information
Authors and Affiliations
Contributions
This research was initially conceived by D.F.P. and B.S.K. Funding was obtained by D.F.P. Sediment core recovery was undertaken by D.F.P. and B.S.K. E.A.N. prepared, processed and interpreted the results of the radiological analyses of the sediment core with the assistance of D.F.P. Figures and tables were prepared by E.A.N. and D.F.P. The initial manuscript was written by D.F.P. in collaboration with E.A.N. and B.S.K. and edited and discussed by all three authors.
Corresponding author
Ethics declarations
Competing interests
The authors declare no competing interests.
Additional information
Publisher's note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary Information
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Niederman, E.A., Porinchu, D.F. & Kotlia, B.S. Hydroclimate change in the Garhwal Himalaya, India at 4200 yr BP coincident with the contraction of the Indus civilization. Sci Rep 11, 23082 (2021). https://doi.org/10.1038/s41598-021-02496-5
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/s41598-021-02496-5
Comments
By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.