Abstract
The heat transfer from the low latitudes to high latitudes is responsible for maintaining the earth’s climate dynamics. Thus, deciphering the possible mechanism driving the variability of the Indian summer monsoon (ISM) during the Holocene Epoch has been critical to understand the hydroclimatic changes of the low latitudes. Despite several efforts, the teleconnection of ISM with the global climate dynamics remains under-represented and poorly understood. The present study aims to delineate the ISM variability and its possible forcing mechanism from western India (Gujarat). In this study, a sediment core (~65 cm long) was raised from the Jaffrabad mudflat (MIT) in western Gujarat. The sediment samples were subjected to geochemical analysis to investigate paleomonsoon, paleo-sediment source and paleoweathering changes. The results show that, with the addition of intermediate sources, the sediments were principally derived from the hinterland’s Deccan basalts. Further, the study suggested a warm and wet climate due to strong ISM during 10,650−5500 cal yr BP associated with the solar as well as orbital forcings. The weak monsoon during 5500−2700 cal yr BP has been linked with southward migration of the Intertropical Convergence Zone (ITCZ) along with the increased El Niño-like conditions. Further, the wavelet analysis revealed that a combined influence of solar, orbital and North Atlantic forcings led to monsoon variability along western India, during the Holocene Epoch. By reconciling the geochemical proxies, the present study has implications in the reconstruction of paleomonsoon and establishing the possible teleconnection with the global climate system.
Keywords
Introduction
The dynamics of the Earth’s climate system rely on the oceanic and atmospheric heat transport from low latitudes to high latitudes. The current interglacial period, the Holocene Epoch (~11.8 ka to present) was initially assumed to be of consistent climate. However, the globally distributed paleoclimate reconstruction provided mounting evidence to demonstrate that the Holocene has evolved as a fluctuating climate modulated by the internal and external forcing mechanism (Banerji et al., 2020; Mayewski et al., 2004; Misra et al., 2019; Wanner et al., 2008). The solar inputs have been the prime factor which can be attributed to long-term changes caused by Earth’s orbital variability while short-term changes have been led by perturbations in the sun itself during the Holocene epoch (Mayewski et al., 2004). In higher latitudes, the enhanced solar insolation of the early Holocene driven by orbital forcing resulted in the stepwise withdrawal of the largest ice sheets of the northern hemisphere (Laurentide Ice Sheet) that prevailed during the last glaciation and finally disappeared ~6.8 ka (Carlson et al., 2008; Clark et al., 1999). While at low latitudes, the increased solar insolation during the early Holocene triggered the northward migration of the mean latitudinal position of the Intertropical Convergence Zone (ITCZ) promoting the boreal summer monsoon rainfall thereby leading to prominent changes in the hydrological cycle (DeMenocal et al., 2000; Fleitmann et al., 2007; Haug et al., 2001). Concisely, it can be underscored that the high latitude climate variability is demonstrated through the growth and decay of ice sheets while the monsoon intensities predominantly govern the climate system of low-latitude regions (DeMenocal and Rind, 1993). Despite the colossal evidence from both high and low latitudes, the response of the low-latitude climate system towards climatic signals emanating from the high latitudes remains severely unaddressed.
The Indian Summer Monsoon (ISM), being a part of the low-latitude Asian Monsoon System exhibits a major control on the global hydrological cycle (Clift and Plumb, 2008). The ISM causes >80 % of rainfall over the Indian subcontinent and nearby regions and generally prevails during June-August as a result of northward migration of the ITCZ. The southward shift in the ITCZ leads to the winter monsoon [western disturbances (WD) from December to February and northeast monsoon (NEM) from October to December] in the selected regions of India (Dimri et al., 2016; Dixit and Tandon, 2016; Resmi and Achyuthan, 2018).
The variations in the ISM strength have been linked with the rejuvenation and collapse of civilisation over thousands of years (Enzel et al., 1999; Pokharia et al., 2017; Prasad et al., 2014; Staubwasser and Weiss, 2006). Moreover, in the last few decades, the conspicuous occurrence of torrential rainfall events or severe droughts (Krishnan et al., 2009; Shaw and Nguyen, 2011) have profoundly impacted agricultural production, economic status, and a few billion of human lives. Such ISM perturbations have invoked the stark need to provide robust predictions of the seasonal monsoon rains with better accuracy that has canonical implications in the socio-economy of the country. However, the direct observational data of ISM rainfall seldom goes beyond the last few decades to centuries. Hence, a comprehensive view of the large-scale mechanism operative on the local, regional and global scale climate system is needed with improved spatiotemporal scale for advanced and accurate prediction of the Indian monsoon and for developing mitigation strategies for the ongoing global warming scenarios (Krishna Kumar et al., 2005).
Gujarat, the westernmost state of India, witnesses varied climatic patterns ranging from humid, subhumid, semi-arid, and arid to the extremely arid zones (Merh, 1995; Prasad et al., 2014) which have been endorsed due to variability in ISM rainfall pattern (Kshetrimayum, 2007). The prevalence of climatic variability led to substantial attempts to address the ISM perturbations during the Holocene Epoch. The multiproxy attempt on the Wadhwana Lake deposit, mainland Gujarat demonstrated a cool and moist climate during ~7500−5560 cal yr BP and stronger ISM after ~3500 cal yr BP, interrupted by a dry climate during 5560−4250 cal yr BP (Prasad et al., 2014). A multiproxy approach on the sediment core from Pariyaj Lake, Vatrak River basin, Gujarat revealed a wet climate and high lake stand during ~11,000 cal yr BP, 7630 cal yr BP and 4680−3500 cal yr BP interspersed by dry climate during 5864−4680 cal yr BP and 8000−9000 cal yr BP. The wet climate during 7630 cal yr BP and 4680−3500 cal yr BP has been attributed to winter rainfall in the Pariyaj Lake region (Raj et al., 2015). The lacustrine sequence of Timbi lake, Dhadhar river basin underscored three large flash floods from 4830 to 2730 cal yr BP and during ~1770 cal yr BP while enhanced monsoon and warm climate prevailed during 1730−880 cal yr BP and 880−360 cal yr BP (Sridhar et al., 2020). The Nal Lake, central Gujarat indicated a transition phase from arid to humid during 6.6−6 ka following which a dry climate interrupted by wet spells was observed during 6−4.8 ka. From 4.8 to 3 ka wettest climate was observed which later trended towards an arid phase till 2 ka and the present climate was attained after 2 ka (Prasad et al., 1997). Likewise, a wet phase before 2185 cal yr BP followed by an arid phase till 1809 cal yr BP has been revealed by the Narmada valley sediments (Sridhar et al., 2015). A multiproxy approach on the three sediment units at the Khari, Saraswati and Sabarmati Rivers demonstrated marine influence and high sea stand during 6070 ± 150 cal yr BP which was followed by high terrestrial flux till 2650 ± 180 cal yr BP interspersed with reduced rainfall and dry climate during 2650 ± 180 and 1530 ± 350 cal yr BP (Sridhar et al., 2018). The palynological attempt in the Mahi Estuary, mainland Gujarat revealed a coupled summer and winter monsoon rainfall during 3400−3000 cal yr BP followed by an arid climate from ~2850 cal yr BP onwards (Prasad et al., 2007). Similarly, two dry events at ~2.1 ka and 1.3 ka preceded by a subhumid condition at ~3 ka has been demonstrated from the lower Narmada valley, mainland Gujarat (Laskar et al., 2013). Further, four episodes of enhanced monsoon such as >5 ka, 4.6 ka and 4.6−1.7 ka were reported from slack water deposits of Mahi River, mainland Gujarat (Sridhar, 2007) while larger flood events were demonstrated during 1400−1440 CE (Sridhar et al., 2014). The mainland Gujarat has been significantly studied for the Holocene monsoon variability.
The Saurashtra peninsula of western Gujarat has provided enormous evidence of being the focal point for maritime activities during the Harappan civilisation (Gaur et al., 2009; Gaur and Bhatt, 2008; Pramanik, 2004). Other studies have combined archaeology and palynology to demonstrate the presence of rich forest cover due to the moist climate during the mid-Holocene which also resulted in the settlement and growth of Harappans along the Saurashtra coast as they significantly sustained based on agriculture, domestication and maritime activity (Farooqui et al., 2013). Considering the present and past socioeconomic significance of the Saurashtra region, there have been attempts to demonstrate the monsoon variability during the mid-Late-Holocene (Banerji et al., 2021, and references therein) using relict and active mudflats of the southern Saurashtra coast. However, a continuous reconstruction of monsoon variability throughout the Holocene epoch with special emphasis on the early Holocene and the response of the Saurashtra region towards global climate scenarios during the early Holocene remains poorly explored.
In general, the sheltered coastal environments such as mudflats, lagoons, estuaries, bays, rias, inlets and isolated basins are the prime realms that are capable of preserving thick Holocene sediment sequences and thus provide an opportunity to reconstruct a continuous paleoclimate and paleoenvironmental reconstruction, while occasionally it may also offer to investigate past sea level and land level changes (Lamb et al., 2006). The southern Saurashtra coast is marked by 40–50 m vertical cliffs of miliolite limestone and extensive tidal mudflats. These active mudflats receives terrestrial contribution only during the monsoon rainfall which leads to the activation of seasonal rivers and during high tides when they receive marine sediments, thereby demonstrating its inimitable characteristics. Considering the above, the present study broadly aims to reconstruct the paleomonsoon variability and decipher the global teleconnection during the Holocene Epoch from the mudflats of the Southern Saurashtra coast, Western Gujarat. The objectives for the present study are: (i) To understand the paleoweathering and paleo-sediment source in the mudflats of southern Saurashtra coast; (ii) To reconstruct the past ISM rainfall variability from the southern Saurashtra region that predominantly witnesses ISM rainfall; (iii) To deconvolute the possible teleconnection of the ISM with regional and global climate scenarios.
Regional settings
Geologically, the Saurashtra peninsula of western Gujarat, India comprises Deccan basalt with Tertiary and Quaternary limestone along with alluvium fringed around the coastal region (Figure 1). The extensive thickness of tholeiitic flood basalts (Najafi et al., 1981), several volcano-plutonic complexes (Naushad et al., 2019), the prevalence of rhyolite and granophyre and extensive compositional diversity (Sheth et al., 2011, 2012) are the prime features that distinguish Saurashtra Deccan basalts from the rest of the Deccan plateau of west-central India. The Deccan trap is unconformably overlain by miliolites, Gaj beds, fluvial gravels and clays (Figure 1). Three age ranges viz. 200–140, 115–75 and 70–50 ka have been identified for the miliolites based on U-Th dating techniques, nevertheless, the age and its origin are still contentious (Baskaran et al., 1986).

Geology of southern Saurashtra along with the seasonal rivers draining into the mudflats of Jafrabad and Rohisa. The MIT core location has been marked by a red star. The Jafrabad mudflat is drained by Raidy River which is active during ISM.
The mean annual rainfall of the Saurashtra region is ~600 mm with a maximum precipitation during the ISM (June to September) (Farooqui et al., 2013). The mean maximum and minimum temperatures vary between 34°C and 19°C respectively (Gundalia and Dholakia, 2013). Based on the forest type classification (Champion and Seth, 1968), the study area falls under 5A/C-1a (very dry teak forest). The dry deciduous teak forest constitutes Terminalia, Diospyros melanoxylon, Tectona grandis, Wrightia tinctoria. Savannah includes Acacia grasses such as Dicanthium annulatum, Sehima nervosum, Apludamutica, Chrysopogon, Heteropogon, Cymbopogon, Aristida. Thorn scrub forest consists of A. catechu, A. marmalos, A. nilotica, B. aegyptica, A. leucophloea, etc. (Farooqui et al., 2013). The previous study at the Vasoj relict mudflat, NE Diu Island, southern Saurashtra coast revealed the prevalence of extensive mangrove forest including core mangrove species (Sonneratia sp., Bruguiera gymnorrhiza, Rhizophora mucronata, and Excoecaria agallocha) during 4.7 ka which gradually diminished as a result of sea regression and monsoon weakening (Banerji et al., 2015). The present study has been conducted on the exposed mudflats of Jaffrabad which is located ~40 km NE of Vasoj village along the southern Saurashtra coast, Gujarat.
Material and methods
The exposed mudflats Jaffrabad, southern Saurashtra coast, Gujarat were selected for the sediment core collection. Geologically, the mudflat is surrounded by biogenic carbonate rocks (miliolite) dated to ~178 ka (Baskaran et al., 1987). Presently the region seldom receives any marine influence, but the seasonal rivers get activated and deposit sediment during the ISM. The modern tidal flat sediments are dominated by clay, whereas the beach sediments are dominated by biogenic carbonate (Hardas and Patel, 1982).
A sediment core of 65 cm in length was raised from the exposed mudflat of Jaffrabad (MIT; Latitude: 20.889°N; Longitude: 71.382°E) with the help of a PVC hand-held corer (Figure 1). The sediment core was sub-sampled at 1 cm intervals with a Perspex knife to avoid metallic contamination. During the sampling of the sediment core, small-sized gastropod shells, broken shell pieces and bivalves were encountered from 26 cm to 65 cm. While the section from surface to 26 cm depth the sediment core was barren in terms of gastropod shells. The subsamples were preserved in zip lock polythene bags in a deep freezer till further analysis.
Chronology
The chronology of the core was established using AMS radiocarbon (14C) dating on the decarbonated sediment samples from the depths of 10.5, 22.5, 46.5 and 62.5 cm. The homogenised crushed bulk sediment samples were decarbonated with 0.5 N HCl to remove the carbonate fraction following which they were repeatedly washed with deionised water. The AMS radiocarbon dating for the organic fraction was conducted at the AURiS-AMS Facility at the Physical Research Laboratory, Ahmedabad, Gujarat (Bhushan et al., 2019a, 2019b).
Total organic carbon (TOC) and total nitrogen (TN)
The total organic carbon (TOC) indicates the in-situ productivity while the ratio between TOC and total nitrogen (TN) is frequently used for organic carbon source (Lamb et al., 2006; Tribovillard et al., 2006). In general, the TOC/TN >12 represents C3 plants with high cellulose and lignin (Prahl et al., 1980) while TOC/TN <10 indicates C4 plant of aquatic origin having high proteins and nucleic acids (Blackburn, 1983).
The total organic carbon (TOC) and total nitrogen (TN) was estimated on ~65 bulk homogenised sediment samples. The sediment samples were decarbonated with 0.1N HCl and thoroughly washed with the deionised water. Later, the dried decarbonated sediment samples were crushed and packed in tin capsules for the measurement of TOC and TN in NC analyzers (FISONS model NA 1500). The NC analyzer was calibrated using Low Organic Soil Sample (LOSS; Batch no. 647582814) as a reference material with 1.65% and 0.14% of carbon and nitrogen contents respectively (Bhushan et al., 2001).
Textural and geochemical analysis
The combined sieving and pipette analysis techniques of Lewis (1984) were used for the estimation of sand, silt and clay contents in the sediments. Nearly 20 samples were treated with 30% H2O2 and 10% HCl followed by wet sieving (diameter: 63 µm) to separate the sand fraction. The remaining fraction consisting of silt and clay was determined using pipette analysis.
Nearly, 33 samples from the MIT sediment core were selected, oven-dried at ~70°C, crushed, homogenised and washed with deionised water to eliminate the salt content in the sediment. Approximately 2 gm of the dried and combusted sediment samples were utilised in the preparation of pressed pellets to estimate the major oxides (SiO2, Al2O3, Fe2O3, MgO, CaO and K2O) and Zr using X-ray Fluorescence (XRF) technique at National Centre for Earth Science Studies (NCESS), Thiruvananthapuram, India. The finely powdered sediment samples were sprinkled over a boric acid binder filled in aluminium cups and pressed with a 40-ton hydraulic press for the 30 s. The measurement was performed on a Bruker S8 Tiger wavelength dispersive (WD) XRF instrument. The reproducibility of the instrument was monitored using international reference material SCO-1 and MAG.
Nearly 0.3 g of bulk sediment samples were digested by treating them with concentrated acids (HCl, HF and HNO3) using a microwave digestion system. The digested samples were analysed for various trace elements (Ba, Cr, Co, Ni, Mn, Ti, Cu) by aspirating the sample solutions in the ICP-MS (Thermo-X series2) and rare earth elements (REE), Th and Sc were measured in HR-ICP-MS (Element XR HR-ICP-MS; Thermo Fisher Scientific) at Geosciences Division, Physical Research laboratory, Ahmedabad, India. The reproducibility of the geochemical data was monitored using NOVA (Agnihotri, 2001; Amin et al., 1972) and MAG. The precision and accuracy for the major and trace elements through XRF and ICP-MS were better than 5% (Banerji et al., 2017; Shaji et al., 2022) while for REE, Th and Sc measured through HR-ICP-MS were better than 7%.
Results
In the present study, the sediment core from Jaffrabad mudflats (MIT) was analysed with geochemical proxies supported by radiocarbon ages and textural variability to investigate the past climate and hydrological variability based on detrital flux, weathering changes, in-situ productivity, and possible influence of the sea level variability.
Chronology
The MIT sediment core was chronologically constrained using four AMS radiocarbon ages. The organic fractions at the depths of 10.5, 22.5, 46.5 and 62.5 cm yielded an age of 3644 ± 73, 4236 ± 61, 9279 ± 61 and 9786 ± 61 years BP respectively. The obtained radiocarbon ages were calibrated using Marine20 using CALIB Rev 8.2 (Heaton et al., 2020) with reservoir age (ΔR) correction of −4.0 ± 57.0 (Dutta et al., 2001; Southon et al., 2002). The calibrated ages obtained were 3391 cal yr BP, 4150 cal yr BP, 9893 cal yr BP and 10,567 cal yr BP for the depths 10.5, 22.5, 46.5 and 62.5 cm respectively (Figure 2a; Table 1). The age-depth model was prepared using the help of the RBACON package in R software (Figure 2b) which utilises the gamma autoregressive semiparametric method and yields the age estimations with the most reliable uncertainty estimation (Blaauw and Christen, 2011; Wang et al., 2019). The sedimentation rates varied from 0.042 to 2.4 mm/year with the lowest sedimentation rate observed during 22−47 cm (4150−9893 cal yr BP) and the MIT represented a sediment depositional history between 2790 and 10,650 cal yr BP. In order to investigate the past monsoon and hydroclimate variability from the MIT core, the respective sample depths have been transformed into ages based on the sedimentation rates and following linear extrapolation using the calibrated radiocarbon ages.

(a) The lithology of the sediment core MIT indicates that most of the core is composed of silty Clay and silty Mud. (b) the age-depth plot for the MIT sediment core generated using RBACON package [Note. zM – silty Mud; cM – clayey Mud; and zC – silty Clay].
Sample details along with the AMS radiocarbon ages and their respective calibrated ages for the MIT sediment core collected from the southern Saurashtra coast.
Textural variation
The textural variation in the sediment core provides cues on the hydraulic activity (Pettijohn, 1957) thereby demonstrating the depositional environment. The textural variation for the MIT sediment core indicated sand, silt and clay ranging from 1.1% to 32.4%, 35.0% to 48.8% and 21.0% to 62.4 % respectively (Figure 3; Table 2). The MIT sediment core is predominantly composed of silty clay from 26 to 60 cm (~5130−10,440 cal yr BP) while silty mud prevails from surface to 21 cm (2790−4100 cal yr BP) and 60−65 cm (10,440−10,650 cal yr BP). However, a conspicuous layer of clayey mud is found at a depth of ~24.5 cm (~4655 cal yr BP). The textural variation and its classification is based on Picard (1971) while the classification of hydrodynamic conditions through the textural variation has been adopted from Pejrup (1988). Thus, the sediment from 24 to 60 cm consisting of silty clay (zC) deposited under calm conditions (Figure 3b and c) while the deposition of silty mud (zM) deposited under relatively higher hydrodynamic conditions (Pejrup, 1988).

(a) The downcore variation of sand, silt, clay and mud indicates that the MIT core is primarily composed of mud. The ternary plots of the MIT core depicting (b) textural variation broadly falling within silty mud, clayey mud and silty clay (Picard, 1971) and (c) energy regimes of the sediment depositional environment varying from calm to marginally violent (Pejrup, 1988). [Note. C – Clay; S – Sand; Z – Silt; zC – Silty Clay; cZ – Clayey Silt; sZ – Sandy Silt; cM – Clayey Mud; zM – Silty Mud; sC – Sandy Clay; sM – Sandy Mud; cS – Clayey Sand and zS – Silty Sand.]. The obtained AMS ages are marked as red star on the depth profile.
Textural variation for the MIT sediment core.
Major and trace elemental variation
In the MIT core, a marginally high TOC content with a gradual lowering has been observed during 10,650−9900 cal yr BP followed by a consistently low TOC until 6000 cal yr BP (Figure 4). From 6000 to 2700 cal yr BP, the TOC content has significantly fluctuated with a noticeable increase between 3300 and 3900 cal yr BP. However, unlike TOC, the TN content consistently remained low from 10,650−5000 cal yr BP while like TOC it also indicated prominent enhancement between 3300 and 4000 cal yr BP. In case of TOC/TN, higher values are observed during ~10,650 cal yr BP with a gradual declining trend in two steps viz during 7000−8000 cal yr BP and ~3900 cal yr BP (Figure 4). Until 8000 cal yr BP, the TOC/TN >12 demonstrating significant contribution of terrestrially derived TOC while after 8000 cal yr BP the region plausibly witnessed more of aquatic environment till 2700 cal yr BP.

Downcore variation of CaCO3, TOC, TN and TOC/TN for the sediment core MIT. The TOC/TN >12 indicated enhanced contribution of terrestrially derived TOC. The obtained AMS ages are marked as red stars on the depth profile.
The SiO2 variation remained nearly consistent with marginal enhancement during 10,350−10,650 cal yr BP and reduction during 5000−6000 cal yr BP. The overall trend of Al2O3 and TiO2 displayed a similar trend (r = 0.71) with higher values between 10,650 and 5000 cal yr BP followed by significant lowering from 5000 to 2790 cal yr BP (Figure 5). Additionally, other major elements such as Na2O, K2O, Fe2O3 and trace elements viz Mn, Cu, Co, and Sr concentration also followed a trend comparable with Al2O3 and TiO2. In the MIT core, the SiO2/Al2O3, TiO2/Al2O3 and Zr/Al2O3 remained nearly consistent during 10,650−5000 cal yr BP with a conspicuous spike in case of TiO2/Al2O3 between 10,350 and 10,600 cal yr BP while a significant enhancement in SiO2/Al2O3, TiO2/Al2O3 and Zr/Al2O3 is inferred during 5000−2700 cal yr BP (Figure 6).

Downcore variation of elemental concentration indicating its mutual association [Note. The concentration of SiO2, Al2O3, TiO2, Na2O, K2O and Fe2O3 are in %; and Mn, Cu and Co are in ppm]. The obtained AMS ages are marked as red stars on the depth profile.

Downcore variation of textural and chemical weathering proxies. The high values of textural proxies corroborated with the presence of increased sand and silt content while chemical weathering agrees with the clay fraction of the sediment core MIT. The obtained AMS ages are marked as red stars on the depth profile.
The chemical weathering is an imperative process in controlling the global hydrogeochemical cycle of elements resulting in the elemental partition between continental rocks and water (Taylor and McLennan, 1985; Zhao and Zheng, 2015). Sediment weathering can result due to microbes, climate, the nature of source rocks, and relief (Joshi, 2014; Madhavaraju et al., 2016; Xi et al., 2014). The chemical weathering has been addressed using Chemical Index of Alteration (CIA) that measures the degree and intensity of chemical weathering (Fedo et al., 1996; Nesbitt and Young, 1982) and the Plagioclase Index of Alteration (PIA) that access the plagioclase weathering (Fedo et al., 1995). The CIA and PIA values in the MIT sediment core ranged from 47 to 61 and 46 to 63 respectively with a gradual increase during 10,650−10,200 cal yr BP and remained consistently high till 5000 cal yr BP following which abruptly declined and remained consistent till 2790 cal yr BP (Figure 6). Additionally, the index of compositional variability (ICV) is also being estimated which deciphers the abundance of alumina with respect to the other major cations in minerals or rocks, and it is also frequently used to study mineralogical maturity and weathering intensity (Cox et al., 1995). In addition to this, the ICV values were also estimated wherein it ranged from 2.7 to 8 and remained consistently low during 1065−5000 cal yr BP followed by enhancement till 2700 cal yr BP. Overall the ICV demonstrated an anticorrelation pattern with CIA and PIA variation. The ICV values suggested immature sediments as the rock forming minerals shows ICV >0.84 while clay minerals such as kaolinite, illite, etc. show <0.84 (Cox et al., 1995). Usually, the Ti/Ca successfully provided cues on Ti-bearing clastic sedimentary flux compared to the in situ carbonates (Hu et al., 2012).
The organic carbon productivity was estimated for the MIT core based on TOC, however, the TOC content is usually prone to degradation and therefore, the Cu/Al2O3, Zn/Al2O3, and Ni/Al2O3 account for paleoproductivity changes. The Cu/Al2O3, Zn/Al2O3 and Ni/Al2O3 demonstrated high values during 10,650−10,050 cal yr BP following which they consistently remained low except a marginal enhancement observed during 3500−5500 cal yr BP reminiscing the TOC trend signifying the in-situ productivity (Figure 7). Usually, the TOC is being annihilated by the aerobic organism that devours the available ambient dissolved oxygen. The continued consumption of oxygen can lead to sub-oxic, anoxic or euxinic conditions depending on the concentration of the oxygen. Under such oxygen-depleted environment, the migration and enrichment of the selected elements prevail between sediment and porewater or between water and sediment-water interface (Banerji et al., 2016; Tribovillard et al., 2006) thereby leading to depletion (Mn and Fe) and enrichment (Cr, Mo, V, etc.) of selected elements, providing cues on paleoredox conditions (Tribovillard et al., 2006). The oxygen-deficient conditions can be estimated using redox-sensitive elements (Banerji et al., 2016; Tribovillard et al., 2006). In the MIT core, the enriched Mo/Al2O3, V/Al2O3, Co/Al2O3 and Cr/Al2O3 while reduced Mn/Al2O3 observed during 10,650−10,200 cal yr BP demonstrated anoxic conditions while an intermittent occurrence of euxinic prevailed during 10,100−10,350 cal yr BP as indicated from the conspicuous spike of Mo/Al2O3.

Comparison of redox-sensitive proxies with the TOC content. A reducing condition has been depicted during 10,440−10,650 cal yr BP as indicated by the Al2O3 normalised V, Mo, Zn, Cr and Mn. The obtained AMS ages are marked as red stars on the depth profile.
Major oxides estimated from MIT core raised from the mudflats of southern Saurashtra plotted on the Harker variation diagram. A broad negative correlation between SiO2 versus Al2O3, K2O, MgO and Fe2O3 (Figure 8) while Al2O3 (Figure 8) versus TiO2, K2O, Fe2O3 and Na2O showed a positive correlation. The binary plot of log (Fe2O3/K2O) versus log (SiO2/Al2O3) indicates that most of the samples fall within the Fe-sand and Fe-shale field (Figure 9). The Al2O3/TiO2, SiO2/Al2O3, Fe2O3/K2O and K2O/Al2O3 ratios in the studied sediments range from 5 to 11, 3.34 to 8.63, 4.84 to 8.21 and 0.09 to 0.20, respectively. To decipher the enrichment and depletion of major and trace elements with respect to the upper continental crust (UCC), the elemental concentrations of the sediment samples from the MIT sediment core of southern Saurashtra have been normalised with the UCC (Figure 9b). It is noted that most of the major elements (TiO2, FeO, CaO and MgO) show enrichment as compared to the UCC.

Binary variation diagram for the selected major oxides of the MIT samples from the Jaffrabad mudflats, southern Saurashtra.

(a) Geochemical classification diagram for the studied sample (Herron, 1988); (b and c) UCC normalised major and trace element concentrations wherein the UCC values taken from elsewhere (Rudnick and Gao, 2013); (d) Chondrite normalised REE diagram for the Jaffrabad mudflats; (e) PAAS normalised REE diagram for the Jaffrabad mudflats.
The relative concentrations of trace elements with respect to UCC are illustrated in Figure 9c. The low field strength elements (LFSEs) such as Rb and Ba are depleted with respect to UCC while Sr remained mostly enriched. The concentration of high field strength elements (HFSEs) such as Y, Zr, Hf, Th, and U are depleted with reference to UCC values. The studied sediments have average Cr and Ni values (Cr = 74.06 ppm; Ni = 41.61 ppm; Cr/Ni = 1.99). The Cr contents of the meta-sediments are lesser than those of PAAS (110 ppm), NASC (125 ppm) and UCC (92 ppm). However, the Cr/Ni ratios of many of the studied samples demonstrate a wide range (1.06–5.02) and on an average they are lower than NASC (2.15) but concurrent with PAAS (2) and UCC (1.96). Additionally, the Sc and Co abundances in the MIT core are comparable with the UCC values.
In terms of the total REE content (∑REE), the MIT sediments range from 2.61 to 110.27 ppm (avg. ∑REE = 64 ppm). In the chondrite normalised REE patterns (Figure 9d) the studied sediments have fractionated REE patterns, with (La/Yb)N ratio varying from 2.69 to 142.13, (La/Sm)N from 2.47 to 4.55. The (Gd/Yb)N ratios (0.75–1.90) suggest a reasonably flat HREE pattern. The MIT samples show weak Eu anomalies (0.83–0.98) except for two samples which have positive Eu anomalies (1.29 and 1.48). The lack of prominent negative Eu anomaly is attributed to the partial weathering of plagioclase feldspar, suggesting the robustness of REE during weathering. In PAAS-normalised (Figure 9e) REE patterns with a majority of samples showing slight LREE depletion [(La/Sm)N = 0.56–1.04], HREE enrichment (Gd/Yb)N = 0.54–3.56; except for a three samples), (La/Yb)N ratios varying from 0.28–2.25 and positive Eu anomalies (Eu/Eu* = 1.27–2.25).
Discussion
The paleohydrological and palaeoclimatic variabilities are usually ascribed by changes in terrestrial contribution, sediment provenance, weathering intensities, redox conditions and paleoproductivity that has been prominently assessed by geochemical variations in the sediments and sedimentary rocks (Claes et al., 2019; Sun et al., 2012). The ionic potential (IP), redox potential (Eh), distribution coefficient and pH are the prime governing factors that controls the mobility of elements in the sediments and sedimentary environment (Banerji et al., 2022a; Bjørlykke, 2010). Further, the climatic parameters substantially supports in the transportation of the elements released during weathering which is either mediated by river or wind while the currents also plays a role in case for marine environment (Bertrand et al., 1996; Sageman and Lyons, 2003).
Sediment maturity and paleoweathering
Major oxides and their ratios have been used for the classification of sediments and to evaluate their maturity and sorting (Bhatia and Crook, 1986; Gupta et al., 2024). Sorting of minerals according to their size, density, and shape can strongly modify the relative abundance of different minerals, and consequently, the concentrations of chemical elements (trace elements) hosted in their lattice (Garzanti et al., 2010, 2011; McLennan et al., 1993a). Das and Haake (2003) suggested that the chemistry of the sediments is dependent on quartz and clay fraction which is related to grain size wherein there is an increase in Al2O3 in finer sediments (clays) and SiO2 in the sands. A negative correlation between SiO2 and Al2O3 suggests control of aluminous clay and quartz content (Figure 8). Major oxide plots against Al2O3 abundances make it possible to compare if the sediments are derived from significantly different source rocks, as Al2O3 is likely to be immobile during weathering and diagenesis. The behavioural pattern of most of the major oxides in MIT sediment of southern Saurashtra mudflat show linear positive trends with increasing Al2O3 (Figure 8), suggesting major influences of hydraulic fractionation (Roy and Roser, 2012) or might suggest their detrital origin and association with clay, feldspar and micaceous minerals (Babeesh et al., 2017).
The decrease in feldspar and clay content is suggestive of increased maturity and is reflected in the increased in SiO2/Al2O3 ratios (Bouchez et al., 2011). The SiO2 is generally more mobile compared to Al2O3 under tropical to subtropical environmental conditions, and the degree of weathering of soil profile can be differentiated by SiO2/Al2O3 ratio (Nathan and Bijilal, 2016) and the SiO2/Al2O3 has been used to decipher the textural maturity and sorting of sediments wherein the igneous and basic rocks have values between 3 and 5 while values greater than 5 suggest the maturity of sediments (Roser et al., 1996). Herron (1988) used Fe2O3/K2O values as indicators of mineral stability while Cox and Lowe (1995) used the ratio of K2O and Al2O3 wherein clays range from 0 to 0.3 while the range between 0.3 and 0.9 represents felspars thereby delineating the composition of original sediments. The studied sediments from the mudflats of southern Saurashtra have SiO2/Al2O3, Fe2O3/K2O and K2O/Al2O3 ratios ranging from 3.34 to 8.63 (avg. = 5.82), 4.84 to 8.21 (avg. = 6.79) and 0.09 to 0.20 (avg. = 0.13) respectively, suggesting low maturity of the studied sediments which contain a considerable amount of clay.
The chemical composition of sediments rests on several geological and geomorphological factors (Joshi et al., 2021a; McLennan et al., 1993b). Various proxies like mineralogy, bulk chemistry and heavy minerals combined with statistical approaches have been widely used to distinguish the source lithologies of sediments and sedimentary rocks as well as for understanding the paleoweathering conditions (Banerji et al., 2022b; Joshi et al., 2021). Additionally, elemental ratios like Mg/Al, Ti/Al and Ti/Ca have been often used as terrestrial and well as proxies for clastic fluxes (Clift et al., 2014; Warrier and Shankar, 2009; Wayne Nesbitt and Markovics, 1997).
The intensity of the chemical weathering in the region has been inferred based on various weathering parameters like CIA, PIA and ICV (Cox et al., 1995; Fedo et al., 1995, 1996; Nesbitt and Young, 1982). The studied samples have a CIA range from 47 to 61 (Figure 10a and b) with an average value of 55 (CIA for UCC = 51 and PAAS = 70–75; Figure 8; McLennan et al., 1993; Sharma et al., 2013) suggesting that the source rocks have undergone weak chemical weathering. The CIA values of MIT sediments are well within the range and in agreement with the Mahi River sediments that drain Deccan traps (37–59; Sharma et al., 2013) and from the previously studied mudflats of Diu (CIA = 48.6–53.3; Banerji et al., 2021). The low CIA values for Mahi River sediments have been attributed to tectonically active and water-starved semi-arid conditions (Sharma et al., 2013) while low values for Diu mudflats were suggested to be caused by the absence of perennial rivers and climate. As the site of the present study is in the vicinity to the previously studied Diu mudflats, the possible reason for low CIA values could be attributed to the lack of perennial rivers and the prevailing climatic conditions. A number of studies have evaluated the degree of weathering and post-depositional K-addition using Al2O3–(CaO + Na2O)-K2O (A-CN-K) ternary diagram (Joshi et al., 2021c). As Al-rich sediments are released during weathering, the sediments tend to move sub-parallel to the A-CN join from their un-weathered source rocks which in turn can be used to decipher the composition of the source lithologies. The weathering trends for the studied samples indicate that it plots close to the feldspar line when extrapolated backward suggesting an intermediate sediment source (Figure 10a). The weathering trend of the mudflat sediment samples was similar to the Mahi River sediments that drain basaltic terrain (Sharma et al., 2013) and the previously studied Diu mudflats (Banerji et al., 2021).

(a) A-CN-K ternary diagram for the MIT core wherein PAAS (Taylor and McLennan, 1995), Saurashtra Basalt (Melluso et al., 2006); Mahi River, average Tholeiite (Sharma et al., 2013); Marine Sediment off southern Saurashtra (Kurian et al., 2013) and sediments from Diu mudflats (Banerji et al., 2021) are plotted for comparison; (b) CIA versus ICV plot for MIT sediment suggest its textural immaturity and witnessed initial chemical weathering; (c) Th/U versus Th plot for the MIT core (McLennan, 1993); (d) Al2O3 versus TiO2 diagram for the Jaffrabad mudflats, southern Saurashtra; (e) A-CNK-FM ternary diagram for Jaffrabad mudflats, southern Saurashtra. PAAS (Taylor and McLennan, 1995), Saurashtra Basalt (Melluso et al., 2006); Mahi River, average Tholeiite (Sharma et al., 2013) and Marine Sediment off southern Saurashtra (Kurian et al., 2013) and sediments from Diu mudflats (Banerji et al., 2021) are plotted for comparison.
The moderate CIA and high ICV values (ICV: 2.66–7.97 with an average of 5, CIA: 47–63 with an average of 56) from the studied sediments also corroborate limited chemical weathering of the source which is also attested by a low textural maturity of the studied sediments. The CIA versus ICV plot for the MIT sediment samples falls in the zone of weak weathering which reaffirms poor weathering of source rocks and generation of immature sediments (Figure 10b). The Th/U ratios of sediments have also been utilised to understand the weathering and recycling (Khan and Sarma, 2023; Ramos-Vázquez and Armstrong-Altrin, 2019). Generally, U+4 oxidises to U+6 due to the loss of U as removable (UO)−2 during recycling and weathering. This process increases the Th/U ratio in sediments which is >4 in case of intense weathering and recycling (McLennan et al., 1995). The studied sediments have Th/U ratios ranging from 0.92 to 2.93 which is much less than UCC (3.89) and PAAS (4.84) suggesting weak weathering and recycling (Figure 10c).
Geochemical imprints of sediment source
Major and trace elements, their ratios, U-Pb isotopes and inclusions in heavy minerals etc. have been frequently used to decipher the source of sediments (Hifzurrahman et al., 2023; Slabunov et al., 2017, 2024). In the UCC normalised major and trace element plot, the marginal depletion of SiO2, K2O and enrichment of MgO and FeO as compared to the UCC reveals the characteristics of the source rock to be mafic to intermediate in composition. The depletion in K2O and Na2O may be associated with a relatively low concentration of feldspars in the sediments. The sediment being derived from a mafic to intermediate source is also revealed by comparable concentrations of Co and Ni to the UCC. The higher partition coefficient of Ti-oxides and Fe-bearing accessory phases can lead to enriched TiO2 and Fe2O3 (Greber and Dauphas, 2019; Joshi et al., 2014) while the enrichment in mobile elements such as MgO and CaO reveals a poor intensity of chemical weathering in the source. Cai et al. (2022) suggested that Al, Ti oxides, Zr and hydroxides have low solubility in low-temperature aqueous medium and the ratio of the sample between Al2O3 and TiO2 is very close to their source. The studied sample from the mudflat of southern Saurashtra is plotted on the Al2O3 versus TiO2 variation diagram (Figure 10d) and the ratio between Al2O3 and TiO2 ranges from 5 to 11 which suggests the basic to an intermediate/mix source (Roser and Korsch, 1988). Additionally, the Al2O3 versus TiO2 diagram (Figure 10d) illustrates that most samples fall in the basic to intermediate source field. Heterogeneity in the source sediments for the southern Saurashtra mudflats can also be suggested in the Fe2O3/K2O versus log ratios of SiO2/ Al2O3 diagram (Figure 9a) where a range from Fe-shale to Fe-sand is noted. The source of the sediments can be further confirmed from the Al2O3-(CaO*+ Na2O + K2O)-(FeO + MgO) (A-CNK-FM; Figure 10e) ternary diagram (Lamine Malick and Ishiga, 2016; Nesbitt and Young, 1989) wherein the studied sediments plot near the feldspar and biotite join near the field of basalts and average tholeiite which further reaffirms limited chemical weathering of the mafic to an intermediate source.
The enrichment and depletion of trace elements in the UCC normalised trace element diagram (Figure 9c) also suggest derivation from a mafic to the intermediate source. As Th has high partition coefficient for felsic lithologies while Sc and Co prefer mafic rocks, their ratios can provide insights into source lithologies. Mafic source lithologies have Th/Sc ratio of <0.8 and Th/Co ratios ranging between 0.04 and 1.40 while the Th/Sc ratio >1 and Th/Co ratio ranging from 0.67 to 19.4 are typical of the upper crust. The studied sediments have Th/Sc and Th/Co values ranging from 0.04 to 0.42 and 0.02 to 0.35 which suggests their derivation from mafic lithologies. This is also supported by the Th/Cr versus La/Cr, Th/Sc versus La/Sc, La/Th versus Hf and Co/Th versus La/Sc plots (Figure 11) wherein the studied sediments fall well below the UCC and felsic compositions. The samples are characterised by low Th/Cr, Th/Sc and high La/Th concentrations thereby suggesting derivation from a mafic to the intermediate source. Furthermore, the studied sediments plot well below the average values of PAAS and UCC in the field of depleted mantle sources suggesting to be derived from mafic sources (Figure 10c).

Trace element discriminant plots of (a) Th/Cr versus La/Cr; (b) Th/Sc versus La/Sc (Bhatia and Crook, 1986); (c) La/Th versus Hf (Floyd and Leveridge, 1987); and (d) La/Sc versus Co/Th plot (Gu et al., 2002) for the MIT core from the Jaffrabad mudflats, southern Saurashtra. Potential sources are plotted for comparison (Condie, 1993).
The relative immobility of REEs during weathering and diagenesis makes them impeccable source indicators. The REE patterns of mafic rocks have low total REEs, (La/Yb)N ratios and smaller or no Eu anomalies whereas those of more silicic rocks usually show higher (La/Yb)N ratios and prominent negative Eu anomalies (Cullers et al., 1997; Cullers and Graf, 1984; McLennan, 1989; Šoster et al., 2020). Low concentrations of total REE and lack of prominent Eu anomalies in the studied sediments indicate a contribution from a mafic source. The low UCC normalised LREE/HREE (<0.1) values further support that the sediments were sourced from mafic lithologies. Similar REE (Figure 10c) patterns were earlier observed from sediments derived from mafic sources in Mahi River sediments (Sharma et al., 2013) and Cauvery catchment in southern India (Sensarma et al., 2008). The above discussion supporting mafic rocks as a dominant source, thereby suggesting that the studied sediments have been majorly derived from the Deccan basalts of Saurashtra.
Holocene climate, periodicities and global linkage
The climate variability of the recent interglacial period, ‘the Holocene Epoch’, has been a prime focus not only due to the climate perturbations but also perceived the growth and development of human society. Recently, the Holocene epoch has been broadly ratified into the Greenlandian (11.8−8.2 ka), Northgrippian (8.2−4.2 ka) and Meghalayan (4.2 ka−Present) Stages by the International Commission on Stratigraphy (ICS; Cohen et al., 2013; Walker et al., 2018). The present section aims to demonstrate the monsoon reconstruction, the interplay of dominant periodicities responsible for the monsoon variability and its global linkages.
Paleomonsoon variability
Present day, the study area is solely influenced by ISM variability causing improved hydrological conditions leading to enhanced terrestrial contribution. A high terrestrial contribution during 10,650−5000 cal yr BP followed by a prominent decline after 5000 cal yr BP can be inferred from the variations in the Al2O3, TiO2, Na2O, K2O, Fe2O3 and trace elements (Mn, Cu, Co; Figure 5). The increased terrestrial contribution is supported by the TOC/TN values>12 indicating dominance of terrestrially derived TOC mediated by improved hydrological environment (Figure 4). Further, the region also witnessed improved chemical weathering accompanied by the deposition of clastic sediments from the hinterland during 10,650−5000 cal yr BP, as indicated by the higher values of CIA and PIA. The simultaneous prevalence of enhanced detrital flux and intense chemical weathering underscores warm and humid conditions during 10,650−5000 cal yr BP. The climatic period corroborated well with the globally recognised Holocene Climate Optimum (HCO) which co-occurred with the strengthening of ISM. The strong ISM strengthening has been demonstrated from Wadhwana Lake, mainland Gujarat during 7500−5560 cal yr BP based on a multiproxy approach (Prasad et al., 2014). Likewise, the palynological investigation of the Nakta Lake sediment, Chhattisgarh revealed warm and humid conditions during 11,700−8500 cal yr BP (Quamar and Kar, 2020). The persistent depletion in the δ18Osw during 12−3 ka from the western Bay of Bengal suggested intensified ISM (Govil and Divakar Naidu, 2011) which corroborated well with the observations made from the eastern Arabian Sea (Govil and Naidu, 2010). Various marine and terrestrial (Banerji et al., 2022c; Haridas et al., 2022; Misra et al., 2019; Quamar et al., 2024) records have underscored strong ISM strength during the globally recognised HCO.
The lowering of elemental concentration of Al2O3, TiO2, Na2O, K2O and other trace elements such as Mn, Cu and Co during 5500−2700 cal yr BP (Figure 5) corroborated well with the reduction in the clay content. While the deposition of coarser fraction during 5500−2700 cal yr BP can be inferred from increased TiO2/Al2O3, SiO2/Al2O3, Zr and Zr/Al2O3 (Figure 6) supported by higher silt content and increased ICV values. Additionally, weak monsoon and reduced humidity can be inferred from reduced chemical weathering demonstrated by low CIA and PIA. Further, during 4000−3300 cal yr BP, the region also witnessed increased in-situ productivity as indicated by enhanced TOC, along with an increase Cu/Al2O3 and Ni/Al2O3 which accounts for improved paleoproductivity (Figure 7). In addition to enhanced in-situ productivity, the peaking of ICV and Zr/Al2O3 (Figure 6) during 4000−3300 cal yr BP demonstrates the climate amelioration during 4000−3300 cal yr BP. The multiproxy study from the mudflats of Vasoj and Diu of the southern Saurashtra coast suggested a warm and humid climate during 4710−2825 cal yr BP (Banerji et al., 2015) and 4105−2640 cal yr BP (Banerji et al., 2017) respectively. Likewise, the study from Nal Sarovar, mainland Gujarat demonstrated the wettest climate during 4.8−3 ka which trended towards arid climate till 2 ka (Prasad et al., 1997) and the slack water deposits from the Mahi River basin elucidated enhanced episodes of monsoon during >5, 4.6 and 4.6−1.7 ka (Sridhar, 2007). The palynological investigation on the swamp, in the NW of Amarkantak, Madhya Pradesh suggested a warm and moderately humid climate during 3600−2761 cal yr BP associated with the strengthened ISM (Chauhan, 2015). An extensive mangrove forest during 7220−4770 years BP followed by its gradual decline since 3500 years BP has been inferred from the Konkan coast based on the palynological study (Limaye and Kumaran, 2012). An enhanced ISM during 10,000−5800 years BP followed and preceded by weak ISM during 11,500–10,500 years BP and 5000–2000 years BP has been demonstrated from the Sanai Tal lake, central Ganga Plain, Uttar Pradesh (Sharma et al., 2004). Thus, it can be suggested that the region witnessed reduced ISM during 5500−2700 cal yr BP with an intermittent period of climate amelioration during 4000−3300 cal yr BP resulted due to improved hydrological conditions in the study area led by the influence of monsoonal activity.
Wavelet analysis
The study of cyclicities in geological processes tend to seek periodicities in data series and explain them in terms of known natural phenomena (Preston and Henderson, 1964). Previous studies have shown cyclic patterns in marine organisms and atmospheric carbon dioxide over the past 500 Ma (Rohde and Muller, 2005). Prokoph and Puetz (2015) noted period-tripling patterns in geological and palaeobiological events on a timescale of ~30–1600 Ma. Several studies have suggested a cyclicity of ~700–900 Myr for the formation of supercontinents or continental growth events (Chen and Cheng, 2018; Joshi et al., 2022). The identification of the periodicities of the geological events in nature provides clues on past cyclic events and delineates future projections (Chen and Cheng, 2018). Generally, a Cross Wavelet Transform (XWT) is performed to reconnoitre the consistency in the phase relationship within the regions in time-frequency space with large common power. The XWT aids in revealing their common power and relative phase in time-frequency space. The continuous wavelet transform (CWT), cross wavelet transform (XWT) and wavelet coherence (WTC) analysis for CIA, ICV and TOC/TN are performed (Figures 12 and 13). The statistically significant common features in the XWT of CIA, ICV and TOC/TN time series (Figure 13) stand out at a 5% error.

The continuous wavelet transforms for (a) CIA, (b) ICV and left parenthesis (c) TOC/TN for the Jaffrabad mudflats. The thick black contour designates the 95% significance level against red noise and the cone of influence (COI) where edge effects might distort the picture is shown as a lighter shade.

(a and b) Cross wavelet transform for CIA-ICV and their Wavelet coherence, (c and d) Cross wavelet transform for CIA and TOC/TN and their Wavelet coherence for the MIT sediment core from southern Saurashtra mudflat. The 95% significance level against red noise is shown as a thick contour. The relative phase relationship is shown as arrows (with in-phase pointing right, and anti-phase pointing left).
In the CWTs, the CIA, ICV and TOC/TN show prominent periodicities of 2048 years and 1024 years from ~5500 to 7500 years BP and ~4000–9000 years BP, respectively (Figure 12). The intensity of the periodicities for CIA and TOC/TN is comparatively higher than ICV. A short pulse of ~256 years periodicity during 7000−8000 cal yr BP and a continuous periodicity of ~512 years during 5000−8000 cal yr BP has been indicated by the CWT of TOC/TN (Figure 12c). The periodicity of ~2050 years corroborated well with the ~2200 years periodicity observed in the sediment core from off Oman and was associated with the oceanic circulation and 14C variation in the atmosphere (Naidu and Malmgren, 1995). Further, the ~1024 years cycle is in line with the ~950 years cyclicity reported from speleothems of Oman (Neff et al., 2001) and lake sediments from Alaska (Hu et al., 2003). The northern hemispheric ∆14C data of tree ring represented both 2200 years and 950 years cyclicity (Lean, 2001) while the atmospheric variation of 14C is attributed to solar activity (Stuiver and Braziunas, 1993) thus, both the presence of periodicities underscores the role of solar forcing. The XWT of CIA and ICV (Figure 13a) demonstrated that they were in-phase only at ~2050 years cycle during 6000−7000 cal yr BP other than this, throughout it remained anti-phased with a cycle of ~1024 years. In view of the present study area, the probable reason for their being anti-phased could be because the CIA is associated with chemical weathering which generally intensifies during warm and humid conditions usually caused by strong ISM while the ICV could be due to activation of local rivers (which otherwise are seasonal) by monsoonal activity led by WD as well. The enhanced winter precipitation has been witnessed by Gujarat during 3600−3400 years BP (Prasad et al., 2007) and during 503−193 cal yr BP (Banerji et al., 2019) wherein the latter has been linked with the enhanced western disturbance. The in-phase behaviour of the CIA and ICV during 6000−7000 cal yr BP could be due to the interplay of both ISM and the WD triggered by solar forcing (~2050 years cyclicity) and reduced TSI (Steinhilber et al., 2009). Nevertheless, throughout the intensity of ISM remained dominant as depicted from the CWT of CIA (Figure 12a). Additionally, the WTC (Figure 13b) indicated a prominent anti-phased scenario of CIA and ICV with cyclicity of ~256 years during 3000−4000 cal yr BP. The time range of 3000−4000 cal yr BP corresponds well with the lesser-known Minoan Warm Period (MWP: ~3300 years BP; Easterbrook, 2016). The periodicity of ~256 years possibly corresponds to ~282 years cyclicity of North Atlantic cold event registered from the Chesapeake Bay, eastern North America which in turn has been linked with the solar variations (Willard et al., 2005). There is a possibility that the study area witnessed the influence of WD linked with the North Atlantic cold climate during 3000−4000 cal yr BP. The possibility of winter rainfall is further supported by the occurrence of ~256 years cyclicity within the 95% significance level against red noise observed in the case of ICV while CIA though represented the same periodicity but remained <95 % significance level. The plausibility of winter precipitation during 3000−4000 cal yr BP can be supported by the previous study on palynological investigation of the Kothiyakhad sedimentary sequence, mainland Gujarat which revealed enhanced winter rainfall during 3600−3400 years BP (Prasad et al., 2007). However, the ~256 years cyclicity could also be representing ~250 years cycle of orbital forcings such as eccentricity, obliquity and precession (Loutre et al., 1992) or ~205 years of sun spot cycle (Solanki et al., 2004). Thus, an interplay of both ISM and WD could have played a significant role during 3000−4000 cal yr BP reframing the monsoonal variability. However, a detailed investigation is needed to delineate the possible forcing mechanism responsible for the climate perturbation during 3000−4000 cal yr BP and its impacts on the Indian monsoon system.
The XWT of CIA and TOC/TN (Figure 13c) was in-phase during 11,000−6500 cal yr BP with a periodicity of ~1024 years while after 6500 cal yr BP it suddenly shifted to anti-phased relationship continuing with the similar cycle (~1024 years). Strengthened monsoonal conditions resulted in enhanced weathering and terrestrially derived TOC causing an in-phased behaviour of CIA and TOC/TN with ~1024 years cyclicity underscoring the control of solar forcing. While after 6500 cal yr BP, the in-situ TOC was dominant thereby leading to anti-phased behaviour of CIA and TOC/TN. Moreover, the WTC (Figure 13d) indicated episodic in-phase events for CIA and TOC/TN with cyclicity of ~128−256 years (4000−5000 cal yr BP), ~256 (7000−8000 cal yr BP), 64−128 years (8000−9000 cal yr BP and 9000−10,000 cal yr BP). The possibility of in-phase events of CIA and TOC/TN could be due to the ISM variability episodes as the periodicities such as ~64 years, ~128 years and ~256 years are in-line with the periodicities obtained for ISM such as ~54 years, ~95 years, ~125 years and ~250 years from the sedimentary records of NE Arabian Sea thereby representing typically of solar forcing (Agnihotri et al., 2002).
Even though the wavelet analysis in the present study represented various periodicities originating from solar, and orbital and also showed the possible occurrence of winter precipitation, however, solar forcing remained the prime driving mechanism resulting in major variations in the monsoon thereby impacting the regional as well as the global climate system.
Global teleconnection
Unlike glacial-interglacial transitions, the Holocene climate remained nearly consistent yet the convoluted global climate system and its association with natural forcings remains poorly understood which restricted the climate modelling and their future climate projections. To decrypt the present climate variability and future projections, the high-resolution climate reconstruction and its association with climate variables need to be addressed on a regional and global scale. In view of this, the present section aims to address the linkage of monsoon variability engrossed in the mudflats of western India, with the global climate system and natural forcing mechanism.
The climate system during the Holocene has been largely regulated by Orbital, Solar and Volcanic forcings while other anthropogenic factors such as vegetation cover, Green House Gases (GHGs) and stratospheric ozone has played a prominent role during the last few centuries leading to a complicated feedback mechanism (Wanner et al., 2008). The enhanced events of the volcanic eruption during 9000−7000 yrs BP ago are prominently seen in the GISP2 ice core of Greenland (Zielinski and Mershon, 1997) resulted due to the crustal re-adjustments caused by the unloading of the glaciers during the termination of the glacial period (Grove, 1976; Zielinski et al., 1996). While the non-existence of such frequent volcanic episodes in the Antarctic ice core (Castellano et al., 2004) further supports the role of crustal readjustments during the early Holocene rather than its role in controlling the climate mechanism. The Greenland ice core (Figure 14a) demonstrated nearly 85% of the volcanic events matched with the documented volcanic eruptions during the last two millennia while only 30% of such events corroborated with documented volcanism during 2−9 ka (Zielinski et al., 1994; Zielinski and Mershon, 1997). In the MIT core, the climate for the last 2700 cal yr BP remains unregistered, thus the implication of volcanic forcing on monsoonal episodes recorded by the MIT core remains unaddressed. However, the previous study on the mudflat sediment core from Rohisa, western Gujarat revealed the possible role of volcanic event occurred after 800 cal yr BP in triggering the onset of the Little Ice Age (Banerji et al., 2019) thereby underscoring the role of volcanic forcing on monsoon variability of western India.

Comparison of the global natural forcings such as Volcanic forcing from GISP2 ice core sulphate records (a; Castellano et al., 2005; Zielinski and Mershon, 1997), Orbital forcing (b; Berger and Loutre, 1994), Solar Insolation, Total Solar Insolation, and 10Be (b-c; Berger and Loutre, 1991; Finkel and Nishiizumi, 1997; Steinhilber et al., 2009), red colour intensity representing ENSO events from the Laguna Pallcacocha (d; Moy et al., 2002), Ice rafted debris (IRD) from North Atlantic (e; Bond et al., 2001), reconstructed the NAO index (e; Olsen et al., 2012) and AMOC strength (e; McManus et al., 2004), Siberian high and Icelandic low represented by K and Na content in the GISP2 ice core (f; Mayewski et al., 1997), ITCZ migration from the Cariaco Basin (g; Haug et al., 2001; Hughen et al., 1996) and Temperature reconstruction through GISP2 icecore (h; Alley, 2004) with the geochemical proxies of the present study depicting chemical weathering (i) and terrestrial contribution (j-k).
Generally, the orbital parameters (obliquity – At 9 ka: 24.230°; present day: 23.45°, perihelion – At 9 ka: 30 July; present day: 3 January and eccentricity – At 9 ka: 0.0193; present day: 0.0167) collectively (Figure 14b) led to differences for the solar radiation during July that was beyond 7% and 25–35 W/m2 over a broader region of latitudes (Berger, 1978; Kutzbach and Guetter, 1986). Overall, a high obliquity and enhanced solar insolation during the commencement of the Holocene possibly led to the warming of the northern hemisphere (Lawrence et al., 2006; Mantsis et al., 2011) with the simultaneous strengthening of the ISM, as the higher solar forcings have been frequently linked with the monsoon strengthening episodes, especially during the early-Holocene (Banerji et al., 2020, 2022c). The conspicuous decline in the cosmogenic isotope such as 10Be (Figure 14c) further supports the improvement in the solar forcing (Finkel and Nishiizumi, 1997; Lal et al., 2005). Thus, the monsoon intensification from 10,650 to 5000 cal yr BP could be associated with the enhanced solar forcing with an interplay of the orbital parameters. Further, there has been a strong association of the ISM with the El Niño Southern Oscillation (ENSO) and North Atlantic Oscillation (NAO) of the Pacific Ocean and Atlantic Ocean respectively (Banerji et al., 2020, 2022c; Srivastava et al., 2017). During the Holocene, a shift towards increased El Niño activity was observed during the Late-Holocene which has been demonstrated by a coupled ocean-atmospheric model (Clement et al., 2000) and the lithic flux in off-Peru marine sediment (Rein et al., 2004) while low activities (Figure 14d) were observed during early Holocene (Moy et al., 2002). In general, the weak ISM or drought-like conditions prevails during the El Niño years (Pant and Parthasarathy, 1981; Rasmusson and Carpenter, 1983; Sikka, 1980). Thus, the ISM weakening after 5000 cal yr BP indicated by terrestrial and weathering proxies observed in the MIT agrees with increased El Niño like conditions. The prominent shift towards weak ISM indicated by the terrestrial and weathering proxies during ~5500 cal yr BP is in-line with the increased ice-rafted debris (IRD; Figure 14e; Bond et al., 2001) and weakening of Atlantic Meridional Overturning (AMOC; Figure 14e and f; McManus et al., 2004) thereby providing a linkage of ISM with the North Atlantic climate. As the ISM variability is a manifestation of the lateral migration of the ITCZ, the study from the Cariaco Basin (Figure 14g) suggested that the ITCZ has marginally shifted during the early-mid-Holocene compared to its prominent shift towards south during the Late-Holocene (Haug et al., 2001). The southward shift of the ITCZ has been associated with the weak ISM and can lead to strong WD resulting in winter precipitation. The plausibility of winter monsoon has been indicated during 3660–3400 years BP based on palynological investigation on the Kothiyakhad sedimentary sequence, Mahi estuary, Gujarat (Prasad et al., 2007). Further, the northern hemisphere temperature as indicated by the temperature reconstruction based on the GISP2 ice core of Greenland (Figure 14h) corroborated well with the trends of the geochemical proxies of the MIT sediment core (Figure 14i–k). The low chemical weathering and reduced detrital contribution in MIT possibly support the ISM weakening thereby underscoring the fact that the mudflat of Jaffrabaad, southern Saurashtra, western India responded in tandem with the regional as well as global climate variability.
Possible signatures of high sea level during early-mid-Holocene
A rapid rise of ~20 m/1000 years in the sea level has been deciphered for the western coast of India during 10,000−7000 years BP while during last 7000 years BP the sea level more or less fluctuated along the present level (Hashimi et al., 1995). In the present study, the occurrence of broken shells and intact shells of bivalve and gastropods between 26 cm (5130 cal yr BP) and 65 cm (10,650 cal yr BP) along with the presence of silty clay from 26 to 60 cm (~5130−10,440 cal yr BP) suggests the plausibility of high sea level which led to flooding of the mudflat and transforming it into a lagoon system. The presence of silty mud during 10,440−10,650 cal yr BP possibly underscores the fact that the rapid rise in the sea level during the onset of the Holocene submerged the coastal mudflat and transformed it into lagoon system. The rapid rise in sea level led to the sudden submergence of coastal plants which resulted in poor water ventilation causing oxygen-deficient conditions. In the sediment core MIT, the oxygen-deficient conditions can be deciphered from the depleted Mn/Al2O3 along with enhanced TOC, V/Al2O3, Zn/Al2O3, Mo/Al2O3, Cr/Al2O3 and Cu/Al2O3 (Figure 7). The present-day elevation of the MIT core site is ~3 m whereas the dating of the oyster shell collected from the Jaffrabad region demonstrated a high sea level of ~3 m during 5.4 ± 0.05 ka (Juyal et al., 1995). A multiproxy approach on the relict mudflat from the Vasoj village, NE of Diu Island indicated a high sea level of ~2 m inclusive of ~1 m of tectonic component (Banerji et al., 2015). As the Saurashtra peninsula witnessed neo-tectonism (Banerji et al., 2015; Bhatt and Bhonde, 2006; Pant and Juyal, 1993), the high sea level of ~3 m must have been affected by tectonic component which remains beyond the scope of the present work. However, it can be prominently underscored that a high sea level prevailed in the region during the early to mid-Holocene period followed by a gradual regression.
Conclusion
The mudflats on the southern coast of Saurashtra get terrestrial contributions only when seasonal rivers are activated. This presents a unique opportunity to study the region’s paleomonsoon, paleo-sediment supply, and paleoweathering. With relatively small contributions from intermediary sources, the major and trace elements geochemistry of the MIT core indicates poor maturity and sorting of the sediments, which mostly generated from restricted chemical weathering of Deccan basalts from the hinterland.
Despite the region experiencing strong monsoonal events during the early Holocene, the chemical weathering remained less operative thereby causing the deposition of immature sediments in the mudflat. Additionally, the sediment’s geochemical values revealed an intermediate to mafic signature, which was likely caused by the basalt from the Deccan traps in the hinterland, which was a major source of terrestrial sediment. Other sources contributed additional material, resulting in signatures with an intermediate or mixed composition, which will require further investigation in future studies.
The hypothesis that the ISM fluctuates in combination with other global climatic events and natural forcing mechanisms is supported by the mudflat records from Jaffrabad in southern Saurashtra, Gujarat, which show a connectivity between the ISM and global climate dynamics. The geochemical analysis of the MIT sediment core revealed ISM intensification with a warm and humid climate during 10,650−5000 cal yr BP followed by a prominent decline in the monsoon. The Holocene climate system was evidently controlled by the solar and orbital forcings, as supported by the trends of geochemical proxies. Further, the monsoon weakening after 5500 cal yr BP is in-line with the increased El Niño events and southward migration of the ITCZ. On the other hand, the existence of both whole and broken shell fragments provided evidence for the high sea level that formed the coastal lagoon system that eventually gave rise to a mudflat around 5500 years ago. However, further geomorphological studies are necessary to elucidate the evolutionary phases of the mudflats of Jaffrabad, southern Saurashtra coast. Additionally, the wavelet analysis also revealed interplay between ISM and WD during 3000−4000 cal yr BP which corresponded to the MWP, but it also invokes the need to explore other evidence for the occurrence of WD during the MWP.
Supplemental Material
sj-xlsx-1-hol-10.1177_09596836241266398 – Supplemental material for Geochemical records of mudflat sediments from southern Saurashtra, Western India: Implications for Holocene climate and global teleconnection
Supplemental material, sj-xlsx-1-hol-10.1177_09596836241266398 for Geochemical records of mudflat sediments from southern Saurashtra, Western India: Implications for Holocene climate and global teleconnection by Upasana S Banerji, Ravi Bhushan, Kumar Batuk Joshi, Ankur Dabhi, AK Sudheer, Chandra Prakash Dubey, Rakesh Kumar Panda, Nayana V Haridas and Mahesh Gaddam in The Holocene
Footnotes
Acknowledgements
USB and RB are thankful to Dr Navin Juyal, Dr. Balaji D. and Mr. J.P. Bhavsar for their kind help and support during the fieldwork. USB and RB are grateful to Prof. J.S. Ray for his support and help during the trace element measurement by Q-ICP-MS. USB, KBJ, and NVH thank the Director of NCESS for his encouragement and support. USB and NVH are thankful to Dr D Padmalal and Dr K. Maya for their encouragement and support. KBJ is thankful to Mr Nishanth for the help in the XRF measurements. NVH acknowledges the research fellowship through the MoES Research Fellow Program (MRFP). Authors extend special thanks to the anonymous reviewers for their constructive comments and suggestions which significantly improved the manuscript.
Author contributions
Funding
The author(s) disclosed receipt of the following financial support for the research, authorship, and/or publication of this article: This research was funded by the Ministry of Earth Sciences, Govt. of India under GEOTRACES Project.
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.
