The traces left by earthquakes in the unlithified sediments, recorded as soft-sediment deformation structures (SSDS), are well reconstructed as palaeo-seismic signals, while the origin of SSDS, seismic vs. Aseismic, is challenging. The present study discusses the origin of SSDS and its implications on palaeoceanography and sediment architecture. In the Middle Jurassic succession of Spiti Himalayan region in India, the topmost part of the Ferruginous Oolitic Formation (FOF) consists of four layers of SSDS and is underlain by the lower member of the Spiti Formation (SF). The sedimentary facies analysis documents the palaeogeographic shift from the middle shelf (carbonate-shale repository: FOF) to the outer shelf (black shale: lower member of SF). The SSDS layers, exhibiting load casts, ball and pillow structures, indicate gravitational instability, while syn-sedimentary faults and insitu breccia are the results of brittle deformation. The dominance of storms in depositional sites often argues for a possible triggering agent for SSDS. Therefore, it was necessary to distinguish between seismic vs. aseismic triggering agents. The lateral continuity, vertical repetition, confinement of SSDS at the top part of FOF and sharp change of facies assemblage indicate seismicity-induced syn-sediment deformation, i.e. seismite. The transition from middle shelf to outer shelf at the onset of seismite indicates that seismic impact possibly caused the rapid subsidence, resulting in the palaeogeographic shift. The rapid transgression is recorded as carbonate-shale repository to anoxic black shale. This study highlights the importance of sedimentological analysis to distinguish the seismite and its implications on palaeogeographic evolution and sedimentary architecture.

Soft-sediment deformation structures (SSDS), typically form during or shortly after sedimentation, but prior to lithification, are observed in various sedimentary environments from mountain to deep sea (Sims, 1973, 1975; Quia et al.1994; Sarkar et al.1995; Owen, 1996; Owen et al.2011; Owen and Moretti, 2011; He et al.2014; Sarkar et al.2014; Laborde-Casadaban et al.2021). The primary processes for SSDS formation, liquefaction and fluidization, governs by various factors, including deformation mechanisms, driving forces, and triggering agents (Owen et al.2011). Liquefaction and fluidization processes are primarily induced by reverse density or gravitational instability (Owen, 1987; Quia et al.1994; Moretti et al. (2016). Several significant research in deep-time sediment sequence relates the SSDS with ancient earthquake (seismite; Seilacher, 1969). The correlation of SSDS with earthquake magnitude and its relation with maximum distance of liquefaction from earthquake epicentre provide a better understanding of tectonics activity (Silva et al.1997; Vanneste et al.1999; Rodríguez-Pascua et al.2000; Greb and Archer, 2007; Salomon et al.2018; Morsilli et al.2020). Recent investigations have revealed that the formation of SSDS can also be triggered by aseismic forces such as storm waves, tidal shearing, rapid sedimentation and overloading (Quai et al.1994; Moretti and Sabato, 2007; Owen et al.2011; Van Loon and Dechan, 2013; Van Loon and Pisarska-Jamroży, 2014; Jamil et al.2021). Studies exploring the relationship between seismic and aseismic triggers of SSDS suggest that neither the deformation structures nor their characteristics can uniquely identify the triggering agents (Owen et al.2011). The researchers engaged in palaeo-seismicity study propose that seismic deformed beds may exhibit lateral persistence, vertical repetition and change of SSDS morphology in lateral continuity (Seth et al.1990; Montenat et al.2007; VanLoon, 2009; VanLoon and Pisarska-Jamroży, 2014; Owen et al.2011; Quai et al.2013; Zhong et al.2022). However, a number of researchers records these characteristics of SSDS beds in storm- and tide-dominated environments (Alfaro et al.2002; Owen et al.2011; Shanmugam, 2016 and references therein). Thus, discerning the seismic repercussions within rock sequences, spanning significant geological periods, always presents a formidable task.

On this controversy, recent studies have put forth the proposition that the seismic activity can lead to changes in facies and facies associations (Owen et al.2011; Moretti and Ronchi, 2011; Sarkar et al.2014; Guo et al. 2023). Highlighting the importance of comprehending the deformation mechanisms and driving forces involved, a number of studies have documented the change in facies and facies association at onset of seismite, such as fluvial to lacustrine (Moretti and Ronchi, 2011), shallow to deep marine (Sarkar et al.2014) and marine to lacustrine (Chen et al.2020; Laborde-Casadaban et al.2021; Yu et al.2022). Despite these findings, there has been limited exploration on the change in facies and facies association for correlating palaeo-seismites. Therefore, the impact of seismic activity on basin palaeoceanography and sedimentation style is poorly constrained.

Here, we focus on the Middle Jurassic succession of Spiti Himalaya, in and around Langza, India, from Ferruginous Oolitic Formation (FOF) to the lower member of Spiti Formation (SF), which was deposited in the southern hemisphere of Tethys along the north-eastern boundary of the Indian plate (Fig. 1). The depositional sequence represents a transition from limestone to anoxic black shale (Bhargava and Singh, 2019; Fursich et al.2021; Albert et al.2021). Several studies focus on the change in palaeoclimate and ocean chemistry (Bhargava and Singh, 2019; Fursich et al.2021; Singh et al.2021), while the cause of facies transition is yet to be explored. The present study records a number of beds with soft-sediment deformations on top of the FOF. The detailed sedimentological study provides insights into the sediment architecture, depositional environment and sea-level fluctuations. The objectives of this study are (a) origin of SSDS, (b) role of seismic vs. aseismic processes on the SSDS formation and (c) impact of seismicity on palaeogeographic evolution and sedimentary architecture. The integrated study of facies analysis, characteristics of SSDS and their origin documents the role of seismicity on basin palaeogeographic evolution and sedimentation architecture.

The Himalayan sedimentary succession deposited along the north-eastern boundary of India in the southern hemisphere of Tethys holds crucial information on Eocambrian to Eocene (Fig. 2a; Bhargava and Singh, 2019; Fursich et al.2021; Albert et al.2021; Sun et al.2021). The Jurassic succession is mainly preserved at the inner part of the higher Himalayas, often representing the scenario during Gondwana’s disintegration. While the Triassic-Jurassic boundary has not been drawn yet, the overall depositional scheme of the Jurassic succession consists of Tangling Formation (TF), FOF and SF, respectively (Bhargava, 2008; Bhargava and Singh, 2019; Fig. 2a). The lithostratigraphic distribution of the Jurassic succession and their respective age is in Table 1. The TF, resting on the Para Formation, consists of the greyish-blue carbonate deposited during Early Jurassic directly transit to FOF with a disconformity in-between (Fig. 2a; Bhargava, 2008; Bhargava and Singh, 2019; Fig. 2a). The FOF consists of alternating beds of oolitic limestone and shale, deposited from Bajocian to Early Callovian (Krishna et al.1982; Bhargava and Singh, 2019 and reference therein; Fig. 2a). The SF is composed predominantly of shale, directly transit from FOF and gradationally passes into the Giumal Formation with sand-shale alternation (Bhargava and Singh, 2019). The lithological variation divided the entire sequence into three major divisions: lower, middle and upper members (Fig. 2a). The fossils constrain Middle Callovian to Late Oxfordian, Kimmeridgian to Early Tithonian, and Early Tithonian to Berriasian age for lower, middle and upper member of SF, respectively (Bhargava and Bassi, 1998; Pandey et al.2013; Singh et al.2021; Fig. 2a; Table 1).

Facies and facies association

The studied sedimentary sequence is dominantly composed of limestone and shale, distinguished as two major parts with laterally persistence-continuous soft-sediment deformation in-between FOF and lower member of SF (Fig. 2b). The lower part of the interval zone is composed of alternate beds of limestone and shale, whereas the upper part consists of thick black shale with occasional calc-arenite beds (Fig. 2b). The lower and upper parts of the SSDS interval identified as facies association I and II, respectively. Although most of the primary sedimentary structures are erased, the preservation of SSDS within limestone suggests that these beds as part of F.A-I. The detailed facies in each facies association is as follows.

Facies Association I: Facies distribution in FOF

The facies association I (F.A-I) represents the FOF of the Kioto Group (Fig. 2a). It is characterized by shallow-water limestone grainstone beds alternating with variable thickness of mud (Fig. 2b). The shale beds mainly consist of planar lamination, while primary sedimentary structures in limestone include massive, hummocky cross-stratified, trough cross-stratification, ripple laminated and planar laminated beds (Fig. 3). The limestone beds consist of abundant broken fragments of fossil invertebrates, with sharp erosional base gradationally transiting to massive, hummocky cross-stratification, trough cross-stratification and ripple laminae (Fig. 3a-d). The hummocky cross-stratification is characterized by series of hummocks and swells, showing erosional contact with planar laminated shale beds (Fig. 3d). The wavelength of hummocks varies from 80 cm to 35 cm. The cross-stratifications are identified by very low-angle cross-bedding (Fig. 3e). The cross-stratification dips 10°–25° towards NW. The tops of the beds are ripple laminated. The planar laminated limestone at the top of the sequence is dominated by micrite (Fig. 3f). The shale beds are mainly grey in colour, preserving planar laminae (Fig. 3g).

Facies Association II: Facies distribution in lower member of SF

The facies association II (F.A-II) represents the lower member of SF, consisting of a thick shale sequence with occasional calc-arenite beds (Fig. 4a-b). The shale is mainly dark black with few inter-bedding of grey calc-arenite beds (Fig. 4a). This planar laminated shale is intervened by calc-arenite beds with sharp and erosional scour base (Fig. 4c). The massive beds show haphazard orientation of fossils like belemnite and bivalves (Fig. 4c). The calc-arenite beds are massive in places, with a thickness varying from 10 to 25 cm. These beds are composed of micrite with abundance belemnite fossils (Fig. 4c-d). The massive beds are graded to ripple lamination in places (Fig. 4d). The planar laminated dark shale shows bed-parallel orientation of belemnites (Fig. 4d).

Characteristics of SSDS

The studied sediment sequence preserves the SSDS at the topmost part of FOF (Fig. 2b). The SSDS beds are composed of fine- to medium-grained carbonate, which are overlain and underlain by shale. The SSDS beds have appeared four times in the stratigraphic sequence, which were documented in two sections (Fig. 2b; Fig. 5a). The SSDS displays a wide range of scale and geometry. The SSDS are mainly syn-sedimentary faults, insitu breccia, load casts and ball-pillow structures (Fig. 5b-f; 6). Detailed description of SSDS are as follows.

Syn-sedimentary faults: Small-scale faults, which vertically transit to normal depositional laminae, are prevalent features within sediment sequences. They predominantly manifest as high-angle normal faults (Fig. 5b), occasionally exhibiting reverse faulting. The fault offsets typically range from a few cm to around 15 cm. These faults are present in a series of sub-parallel flexure associated with slide planes (Fig. 5c). The pole plots of these faults show two prominent orientations: one trending towards the northwest and the other towards the southeast (Fig. 5c). This layer is overlain by undeformed planar laminated shale (Fig. 5b, c). Furthermore, these faults commonly coincide with faulted breccia.

Load casts: Load casts are recognized in heterolithic succession at the interface between carbonate and shale bed (Fig. 5d). They display flat tops and irregular concave up base (Fig. 5d, e). Load structures show slight to deep penetration into the underlying shale layer and range in depth from few cm to 15 cm. Few structures show single loads, while many of them display multiple load collages (Fig. 5e). The multi-loads are asymmetric in fashion. The single and multi-load structures are associated with flame structures (Fig. 5e). The overall thickness of the bed is 40 cm, which shows a massive appearance.

Ball and Pillow structures: The ball and pillow structures occur within limestone-shale inter-bedding (Fig. 5f). The thickness of this deformed bed is about 45 cm. The deformed bed exhibits isolated ball- and pillow-like structures. Entire bed shows detached ball and pillow features and trace of original laminae is rarely preserved. The maximum diameter of ball/pillow is around 30 cm.

Faulted breccia: The breccia consists of irregular, sub-angular to angular carbonate clasts suspended within a matrix of mudstone (Fig. 6a-b). The clasts exhibit a predominantly random orientation, varying in length from ∼2 cm to very large (up to 1 m). The larger breccia clasts display a distinctive rectangular shape, retaining remnants of the original bedding and lamination, providing clues to their source (Fig. 6a-b). The jagged boundaries of clasts show interlocking relationship, revealing the original bedding patterns, which often become disrupted along faulted boundaries (Fig. 6a-b). The average thickness of breccia bed is ca. 1.5 m. This bed occurs on the undeformed limestone and is covered by thin (∼5 cm) shale bed.

Depositional environment and Palaeobathymetry

The studied sections of Jurassic Spiti Himalaya consist of two facies association with soft-sediment deformation beds at their respective interval. The deformed beds rest on F.A-I and are composed of alternative carbonate and shale beds. The presence of massive beds, hummocky cross-stratification and cross-stratification within the limestone indicate high-energy flow. The transition of massive bed to ripple lamination through hummocky cross-stratification within a single bed infers a waning flow. The sharp erosional base of the limestone also supports the high-energy storm. The ripple laminated limestone may have been deposited in the waning flow of storm or wave influence. The deposition of shale and carbonate in the same environment infers depositional condition was fluctuating abruptly. Considering the underlying thick carbonate sequence of Tagling Formation, we infer the depositional environment within the storm-dominating middle shelf setting.

The dominance of black shale with occasional calc-arenite beds in F.A-II suggests greater depth than the F.A-I, where the planar laminated shale deposited from suspension falls out in a calm and quiet environment. The earlier studies on palaeoredox suggested that the depositional site was dysoxic-anoxic (Fursich et al.2021; Albert et al.2021). The absence of coarse-grained sediment in an anoxic zone infers depositional site below the wave and storm base. In the absence of wave and storm, the sharp scour base of massive calc-arenite beds suggests that these are the deposits of slope failure or turbidite (Bouma, 1962; Shanmugam, 2021). These massive beds may also be related to hybrid beds like turbidites-debrite couplets (Jamil et al.2021). The development of hybrid beds is highly influenced by frontal and lateral lobe settings. The frontal lobes initially deposit coarse sandstones, which is followed by the gradual evolution of flow conditions (Jamil et al.2021). The haphazard orientation of fossils in the massive bed also supports the context of rapid sedimentation. The thin ripple lamination on massive beds may be related to the spreading of sediments in lateral lobes starting from an off-axis area that provides a small avenue for the transition of flow and the deposition of hybrid beds (Davis et al.2009; Fonnesu et al.2016, 2018; Jamil et al.2021). In general, progradation of hydride bed movement in deep-sea progress on the erosional surface (Hofstra et al.2016; Kane et al.2017; Jamil et al.2021). Considering the thick black shale deposits, we infer that F.A-II is deposition of deep marine, outer shelf.

Origin of SSDS: Seismic vs. Aseismic triggering processes

The studied sequence consists of shale-limestone repository (F.A-I) and thick black shale (F.A-II), with laterally extended four repetitive SSDS beds on top of the F.A-I (Fig. 2b, 5a). The SSDS include syn-sedimentary faults, load casts, ball-pillow structures and insitu breccia. Each of these deformed beds are overlain and underlain by the undeformed beds. The syn-sedimentary faults, load casts, and insitu breccia occur in both the studied sections, while ball and pillow structure is observed in one of the sections.

SSDS form in unlithified sediments during or soon after deposition, exhibit plastic and/or brittle deformations (Owen, 1987; Moretti and Sabato, 2007). The deformation mechanism enables to overcome the yield strength of solid materials, which temporally behave like a liquid manner and deform as a viscous fluid (Owen, 1987), ultimately resulting the liquefaction and fluidization. Liquefaction and fluidization processes obscure the primary lamination and fabric of the sediment deposits (Obermeier, 1996). The driving forces responsible for the SSDS include gravitational body forces, unstable density gradients, shear stresses, and unevenly distributed loads, reflecting the SSDS morphology (Owen, 1987).

Deformation mechanism and driving force

The deformation mechanism classifies the distribution of SSDS in the studied succession into two groups: liquefied and fluidized structures (load casts, ball-pillow structures) and plastic or brittle deformation (syn-sedimentary faults and insitu breccia). The deposition of high-density limestone bed on low-density muds created an inverted density gradient (Anketell et al.1970; Owen, 1987; Quai et al.1994; Owen et al.2011). The liquefaction induced by the unstable density at the limestone-shale interface formed the load structures (Owen, 1987). In a similar fashion, in the joint action of density gradient and gravity instability, liquefaction and fluidization of underlying beds leads to the formation of the ball-pillow by the separation of the loads (Moretti et al.2001; Quai et al.2013). The alignment of flames along a specific direction suggests strong shear force direction. Ball and pillow structures commonly form by the liquefaction-fluidization processes, triggered by seismic shocks (Owen, 1996; Molina et al.1998; Moretti et al.1999; Tian Hongshui et al.2016a, b; Zhong et al.2019).

Syn-sedimentary faulting in limestone beds demonstrates brittle behaviour, as indicated by the discontinuities caused by faulting in alternating beds of mud and sand (Owen, 1987). The inclined micro-fault blocks exhibit a parallel stair-step pattern, which reflects the influence of gravitational forces (specifically lateral stress) and/ or slumping on the deformation process (Salomon et al.2018). During the slump initiation, layer-parallel normal and reverse faults form at the head and toe of the slump, respectively (Salomon et al.2018). The fluidization pressure during deformation damaged the harder limestone beds and formed breccia (Quai et al.1994; Yu et al.2022). The semi-plastic deformation mechanism of the shale and limestone layers forms breccia.

Argument of seismic vs. aseismic origin

Furthermore, although both brittle-ductile deformation and liquefaction/fluidization contribute to the studied SSDS, analysis of the triggering agents in marine settings indicates that these structures can form due to sediment overloading (Allen, 1982; Moretti et al.2001; Moretti and Sabato, 2007), tidal shear (Tessier and Terwindt, 1994; Greb and Archer, 2007), tsunamis (Takashimizu and Masuda, 2000; Nanayama et al.2000; Schnyder et al.2005; Le Roux et al.2008), wave-storm pounding (Molina et al.1998) and seismic shocks (Sim, 1975; Owen et al.2011; Moretti and Ronchi, 2011; Sarkar et al.2014; He and Qiao, 2015). The preservation of undeformed beds on top of SSDS beds suggests that deformation must have been triggered at the sediment-water interface; otherwise, liquefaction and fluidization would have affected the overlying beds (Owen, 1996). This discounts the possibility of overloading as the main trigger. The tidal shearing also can form a similar type of SSDS via storm pounding-induced liquefaction and fluidization, but the absence of any tidal signal in the depositional site also discards the tidal shearing effect. The dominance of wave and storm signatures in the depositional environment argues as a possible triggering agent. However, the confinement of the SSDS layer at the topmost part of the FOF infers that seismicity must have triggered deformation. We interpret the soft sediment deformation structures in the Jurassic succession of Spiti Himalaya as seismite. The seismicity favours capturing lateral continuity of the deformation structures over several hundred metres, vertical repetition, and change of facies assemblages from F.A-I to II (Owen et al.2011; Moretti and Ronchi, 2011; Sarkar et al.2014; He et al.2014; Chen et al.2020; Guo et al. 2023).

Relationship of the seismicity and SSDS: Magnitude of Earthquake

While the differential liquefaction and fluidization is the main agent to form the soft-sediment deformation, the final morphology of SSDS is controlled by the driving force. SSDS not only identify the palaeo-seismite but also measure the earthquake magnitude and faulting activity (Quia et al.2013; He et al.2014, 2015; Zhong et al.2022). The SSDS in the studied sediment sequence provide the intensity of palaeo-earthquake as well (Allen, 1986; Rodríguez-Pascua et al.2000). Some authors stated that intensity of an earthquake greater than 5 (M≥5) can produce liquefaction and fluidization (Atkinson, 1984; Audemard and De Santis, 1991; Galli, 2000; Wheeler, 2002; Zhong, 2017; Zhong et al.2022; Guo et al. 2023). The computer modelling suggests M≈5 produce liquefaction feature within 4 km radius of the earthquake epicentre (Atkinson, 1984), but a few studied debates on M≈4.5, M≈5.5, and M≈7 for liquefaction and fluidization structures depending on the grain size (Marco and Agnon, 1995; Ambraseys, 1988; Obermeier, 1996). The preservation of load casts, ball-pillow structures indicates that the studied sediment was affected by at least an earthquake magnitude of 5.0–7.0 (Rodríguez-Pascua et al.2000). SSDS in four repetitive beds suggests the sediment succession may have affected four seismic pulses.

Response of basin on seismic impact: Carbonate shelf to black shale transition

The vertical transition from F.A-I to F.A-II provides clear evidence of a distinct shift from a mixed siliciclastic (shale)-carbonate repository to black shale. The overall increase of mud and absence of carbonate in the F.A-II sequence define the demise of carbonate platform (Fig. 7a). Increasing mud and the disappearance of carbonate suggest that the depositional site may have entered below the photic zone via transgression (Fig. 7a). The onset of sea-level rise is marked by the change from mixed carbonate repository to black shale. The occurrence of the seismite at the junction between two contrasting depositional setups suggests the role of seismicity on the environmental change. The facies immediately overlying on the seismites represent the deepest among all facies associations, implying a rapid rise in relative sea level to its maximum extent. This rapid transgression is defined as an event of basin subsidence, including the formation of syn-sedimentary faults and breccia. This subsidence resulted in a sudden shift from the middle to the outer shelf, evident through carbonate setback and thick anoxic black shale (Fig. 7b-c). The regional distribution of the SSDS, coupled with the evidence of a sudden and sharp increase in sea level, can only be satisfactorily explained by intense tectonic activity, including earthquakes that deformed the sedimentary layers at the surface, leading to rapid basin subsidence (Fig. 7b-c; He et al.2014; Sarkar et al.2014). It is important to note that sediment supply remained low during this period, resulting in the dominance of fine-grained suspension fall-out and the absence of coarse sediment in the depositional site. The abrupt palaeogeographic changes, observed during the onset of seismites, are further supported by turbidity flows. The occurrence of turbidites in F.A-II suggests that the tectonic subsidence may have created some slopes, along which calc-arenite beds were transported from shallow-water environments and preserved in the deep marine setting (Fig. 7c). The significant positive shift observed across the seismite interval provides additional force for the water-level rise and develops anoxic conditions, leading to sedimentation (Fig. 7c).

Implication in region study and global perspective

The records of seismite in the studied Jurassic succession indicate the occurrence of palaeo-seismic belt close to the depositional setting (Quai et al.2013). The seismic impact subsided the basin and caused the transition of the depositional environment from mixed siliciclastic-carbonate repository to anoxic black shale. The rapid subsidence increased the sea level without affecting the eustatic sea level (Moretti and Roonchi, 2011; Sarkar et al.2014; He et al.2014). On the deformed Himalayan belt, the basin-wide correlation of sedimentary succession is always challenging. However, the seismite at the boundary between FOF and the lower member of SF is very useful for the stratigraphic correlation. The seismic pulses often directly correlate with the convergence and divergence of plates (Chen et al.2020; Laborde-Casadaban et al.2021; Mueller et al.2023). Considering the Middle Jurassic Tethys sea-floor deposition, we infer that the Middle Jurassic Gondwana breakup and synchronous active rifting (Lovecchio et al.2020) may be linked to the studied seismic pulses.

The seismite in the present study includes the impact of seismicity across the basin and is associated with a change in the palaeogeography, which is reflected in facies distribution and sedimentation architecture. Therefore, the sedimentological study should be the key for the identification of seismite for basin-wide correlation.

Considering the SSDS distribution and sediment dynamics pattern at onset of the SSDS beds, we draw the following conclusion.

  • Laterally persistent four repetitive deformed beds at the topmost part of FOF result from the seismic impact. The liquefaction and fluidization were induced by reverse density/ gravity instability (for load casts, and ball-pillow structures) and plastic or brittle deformation (syn-sedimentary faults and insitu breccia) provoked by seismic shocks.

  • The dominance of strong waves and storm in the depositional site argue against the seismic impact, while the sharp change of facies distribution interprets seismicity is the only possibility. Recurrence of SSDS beds at four intervals suggests that the depositional site was affected by earthquake at least four times. The association of SSDS features indicates moderate to strong magnitude earthquakes.

  • The seismic impact rapidly subsided the basin, which increased the sea level. Basin subsidence caused the shift of the depositional environment from middle shelf to outer shelf. Further, the deepening of basin with climatic calamity enhanced anoxicity, causing transition from carbonate-shale repository to black shale deposition.

  • The seismite beds may be useful for the stratigraphic correlation, demarcating the boundary between FOF and SF along the palaeo-Tethyan coast of India. The Gondwana breakup during Middle Jurassic and related rifting may be the reasons for earthquakes. The study highlights the importance of sedimentological context of the stratigraphic record for distinguishing seismic vs. aseismic impacts.

The authors acknowledge the infrastructure support provided by their host institutes. SM, AS, PRU and RSN acknowledge Director, Birbal Sahni Institute of Palaeosciences, for providing fieldwork and chemical analysis funds. SM, AS and RSN acknowledge the in-house project-2 of Birbal Sahni Institute of Palaeosciences, Lucknow, India. SB is thankful to the Indian Institute of Technology Bombay. The BSIP manuscript no. is BSIP/RDCC/Publication no. 10/2023-24.

The authors declare no conflict of interest.

This is an Open Access article, distributed under the terms of the Creative Commons Attribution licence (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted re-use, distribution and reproduction, provided the original article is properly cited.