1. Introduction
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, Reference Sims1973, Reference Sims1975; Quia et al. Reference Quia, Song, Gao, Peng, Li, Gao, Song and Zhang1994; Sarkar et al. Reference Sarkar, Banerjee and Chakraborty1995; Owen, Reference Owen1996; Owen et al. Reference Owen, Moretti and Alfaro2011; Owen and Moretti, Reference Owen, Moretti and Alfaro2011; He et al. Reference He, Quia, Jiao, Xu, Cai, Guo and Zhang2014; Sarkar et al. Reference Sarkar, Choudhuri, Banerjee, van Loon and Bose2014; Laborde-Casadaban et al. Reference Laborde-Casadaban, Homberg, Schnyder, Borderie and Raine2021). The primary processes for SSDS formation, liquefaction and fluidization, governs by various factors, including deformation mechanisms, driving forces, and triggering agents (Owen et al. Reference Owen, Moretti and Alfaro2011). Liquefaction and fluidization processes are primarily induced by reverse density or gravitational instability (Owen, Reference Owen, Jones and Preston1987; Quia et al. Reference Quia, Song, Gao, Peng, Li, Gao, Song and Zhang1994; Moretti et al. (Reference Moretti, Alfaro and Owen2016). Several significant research in deep-time sediment sequence relates the SSDS with ancient earthquake (seismite; Seilacher, Reference Seilacher1969). 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. Reference Silva, Cañaveras, Sánchez-Moral, Lario and Sanz1997; Vanneste et al. Reference Vanneste, Meghraoui and Camelbeeck1999; Rodríguez-Pascua et al. Reference Rodríguez-Pascua, Calvo, De Vicente and GómezGras2000; Greb and Archer, Reference Greb and Archer2007; Salomon et al. Reference Salomon, Grasemann, Plan, Gier and Schöpfer2018; Morsilli et al. Reference Morsilli, Giona Bucci, Gliozzi, Lisco and Moretti2020). 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. Reference Quia, Song, Gao, Peng, Li, Gao, Song and Zhang1994; Moretti and Sabato, Reference Moretti and Sabato2007; Owen et al. Reference Owen, Moretti and Alfaro2011; Van Loon and Dechan, Reference Van Loon, Han and Han2013; Van Loon and Pisarska-Jamroży, Reference Van Loon and Pisarska-Jamroży2014; Jamil et al. Reference Jamil, Siddiqui, Umar, Usman, Ahmed, Rahman and Zaidi2021). 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. Reference Owen, Moretti and Alfaro2011). 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. Reference Seth, Sarkar and Bose1990; Montenat et al. Reference Montenat, Barrier, Ott d’Estevou and Hibsch2007; VanLoon, Reference Van Loon2009; VanLoon and Pisarska-Jamroży, Reference Van Loon2014; Owen et al. Reference Owen, Moretti and Alfaro2011; Quai et al. Reference Quia and Guo2013; Zhong et al. Reference Zhong, Jiang, Li, Su, Xu, Liang and Fan2022). However, a number of researchers records these characteristics of SSDS beds in storm- and tide-dominated environments (Alfaro et al. Reference Alfaro, Delgado, Estévez, Molina, Moretti and Soria2002; Owen et al. Reference Owen, Moretti and Alfaro2011; Shanmugam, Reference Shanmugam2016 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. Reference Owen, Moretti and Alfaro2011; Moretti and Ronchi, Reference Moretti and Ronchi2011; Sarkar et al. Reference Sarkar, Choudhuri, Banerjee, van Loon and Bose2014; Guo et al. Reference Guo, He and Li2023). 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, Reference Moretti and Ronchi2011), shallow to deep marine (Sarkar et al. Reference Sarkar, Choudhuri, Banerjee, van Loon and Bose2014) and marine to lacustrine (Chen et al. Reference Chen, Zhong, Ogg, van Loon, Chen, Yang, Liu and Xu2020; Laborde-Casadaban et al. Reference Laborde-Casadaban, Homberg, Schnyder, Borderie and Raine2021; Yu et al. Reference Yu, Liang, Tian, Han, Wang and Zhao2022). 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, Reference Bhargava and Singh2019; Fursich et al. Reference Fürsich, Alberti, Pandey, Chasker and Bhosale2021; Albert et al. Reference Alberti, Fürsich, Pandey, Andersen, Garbe-Schönberg, Bhosale, Chaskar and Habermann2021). Several studies focus on the change in palaeoclimate and ocean chemistry (Bhargava and Singh, Reference Bhargava and Singh2019; Fursich et al. Reference Fürsich, Alberti, Pandey, Chasker and Bhosale2021; Singh et al. Reference Singh, Deori, Pandey, Shekhawat, Verma, Banerjee and Sarkar2021), 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.
Figure 1. Global paleogeography during Early Jurassic (a, modified from Damborenea, Reference Damborenea2002, note the northern boundary of Indian continental margin along Tethys marked by blue rectangle), geological sketch map of Himalaya and position within India (b, map modified after Bhargava and Bassi, Reference Bhargava and Bassi1998). Note the enlarged view of the Spiti region and road map showing location of studied sections, marked by the red stars.
2. Geological setting
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, Reference Bhargava and Singh2019; Fursich et al. Reference Fürsich, Alberti, Pandey, Chasker and Bhosale2021; Albert et al. Reference Alberti, Fürsich, Pandey, Andersen, Garbe-Schönberg, Bhosale, Chaskar and Habermann2021; Sun et al. Reference Sun, Richoz, Krystyn, Grasby, Chen, Banerjee and Joachimski2021). 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, Reference Bhargava2008; Bhargava and Singh, Reference Bhargava and Singh2019; 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, Reference Bhargava2008; Bhargava and Singh, Reference Bhargava and Singh2019; Fig. 2a). The FOF consists of alternating beds of oolitic limestone and shale, deposited from Bajocian to Early Callovian (Krishna et al. Reference Krishna, Kumar and Singh1982; Bhargava and Singh, Reference Bhargava and Singh2019 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, Reference Bhargava and Singh2019). 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, Reference Bhargava and Bassi1998; Pandey et al. Reference Pandey, Pathak and Krishna2013; Singh et al. Reference Singh, Deori, Pandey, Shekhawat, Verma, Banerjee and Sarkar2021; Fig. 2a; Table 1).
Figure 2. Lithostratigraphy and biostratigraphy distribution in parts of Spiti Himalaya sequence (a, note the age of studied succession marked by pink color), lithological log in two best locations from FOF to lower member of SF (b, note the SSDS zone at top part of FOF, distance from one section to another section ∼1 km).
Table 1. Lithostratigraphy and age of Jurassic sediment succession (Modified after Fursich et al. Reference Fürsich, Alberti, Pandey, Chasker and Bhosale2021)
3. Result
3.a. 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.
3.a.1. 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).
Figure 3. Facies distribution in FOF (Facies association-I): a single limestone bed preserved massive character, hummocky cross-stratification and planar lamination in ascending order with sharp, scour base on shale (a, hammer length 35 cm, note scour base marked by red arrows), sketch of Fig. 6a (b), abundant fossils in massive bed with sharp scour base (c, camera cover diameter 5 cm), hummocks and swales are marked by blue and red arrows respectively (d, pen length 15 cm), cross-stratified limestone with sharp, erosional base (e, hammer length 35 cm), planar laminated limestone facies (f, camera cover diameter 5 cm), alternating shale and limestone beds (g, marker length 15 cm).
3.a.2. 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).
Figure 4. Facies distribution in lower member of SF (Facies association-II): thick shale sequence with occasional calc-arenite beds (a, hammer length 35 cm, note the sharp, erosional base of calcareous beds are marked by arrows), log of Fig. 4a and primary sedimentary structures distribution in association of calcareous bed (b, note count of fossils), massive bed with sharp, scour base (c, lens cover diameter 5 cm), ripple laminated and planar laminated shale (d, lens cover diameter 5 cm).
3.b. 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.
Figure 5. Soft-sediment deformation structures zone in studied section(a, bar length 5m; note the SSDS zone marked by red dotted line), syn-sedimentary faults (b, bar length 30 cm); abundance slide planes at the fault plane (c, bar length 30 cm, note the pole of slide inset), laterally continuous load cast bed (d, bar length 1m; note the loads are protruded in underlying shale layer; see asymmetric load in inset); enlarge view of load casts marked by yellow rectangle in Fig. 5d (e, note flame structures are marked by yellow arrow), ball-pillow structures (f, bar length 15cm; ball structures are marked by white dotted lines).
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.
Figure 6. Soft-sediment deformation structure: insitu breccia (a, bar length 1 m; enlarge view in inset); sketch of insitu breccia bed (b).
4. Discussion
4.a. 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. Reference Fürsich, Alberti, Pandey, Chasker and Bhosale2021; Albert et al. Reference Alberti, Fürsich, Pandey, Andersen, Garbe-Schönberg, Bhosale, Chaskar and Habermann2021). 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, Reference Bouma1962; Shanmugam, Reference Shanmugam2021). These massive beds may also be related to hybrid beds like turbidites-debrite couplets (Jamil et al. Reference Jamil, Siddiqui, Umar, Usman, Ahmed, Rahman and Zaidi2021). 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. Reference Jamil, Siddiqui, Umar, Usman, Ahmed, Rahman and Zaidi2021). 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. Reference Davis, Haughton, McCaffrey, Scott, Hogg and Kitching2009; Fonnesu et al. Reference Fonnesu, Patacci, Haughton, Felletti and McCaffrey2016, Reference Fonnesu, Felletti, Haughton, Patacci and McCaffrey2018; Jamil et al. Reference Jamil, Siddiqui, Umar, Usman, Ahmed, Rahman and Zaidi2021). In general, progradation of hydride bed movement in deep-sea progress on the erosional surface (Hofstra et al. Reference Hofstra, Pontén, Peakall, Flint, Nair and Hodgson2016; Kane et al. Reference Kane, Pontén, Vangdal, Eggenhuisen, Hodgson and Spychala2017; Jamil et al. Reference Jamil, Siddiqui, Umar, Usman, Ahmed, Rahman and Zaidi2021). Considering the thick black shale deposits, we infer that F.A-II is deposition of deep marine, outer shelf.
4.b. 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, Reference Owen, Jones and Preston1987; Moretti and Sabato, Reference Moretti and Sabato2007). 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, Reference Owen, Jones and Preston1987), ultimately resulting the liquefaction and fluidization. Liquefaction and fluidization processes obscure the primary lamination and fabric of the sediment deposits (Obermeier, Reference Obermeier1996). 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, Reference Owen, Jones and Preston1987).
4.b.1. 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. Reference Anketell, Cegla and Dzulynski1970; Owen, Reference Owen, Jones and Preston1987; Quai et al. Reference Quia, Song, Gao, Peng, Li, Gao, Song and Zhang1994; Owen et al. Reference Owen, Moretti and Alfaro2011). The liquefaction induced by the unstable density at the limestone-shale interface formed the load structures (Owen, Reference Owen, Jones and Preston1987). 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. Reference Moretti, Soria, Alfaro and Walsh2001; Quai et al. Reference Quia and Guo2013). 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, Reference Owen1996; Molina et al. Reference Molina, Alfaro and Moretti1998; Moretti et al. Reference Moretti, Alfaro, Caselles and Canas1999; Tian Hongshui et al. Reference Hongshui, Van Loon, Hailin, Shenhe and Jiewang2016a, b; Zhong et al. Reference Zhong, Jiang, Li, Xu, Shi, Zhang and Wei2019).
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, Reference Owen, Jones and Preston1987). 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. Reference Salomon, Grasemann, Plan, Gier and Schöpfer2018). During the slump initiation, layer-parallel normal and reverse faults form at the head and toe of the slump, respectively (Salomon et al. Reference Salomon, Grasemann, Plan, Gier and Schöpfer2018). The fluidization pressure during deformation damaged the harder limestone beds and formed breccia (Quai et al. Reference Quia, Song, Gao, Peng, Li, Gao, Song and Zhang1994; Yu et al. Reference Yu, Liang, Tian, Han, Wang and Zhao2022). The semi-plastic deformation mechanism of the shale and limestone layers forms breccia.
4.b.2. 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, Reference Allen1982; Moretti et al. Reference Moretti, Soria, Alfaro and Walsh2001; Moretti and Sabato, Reference Moretti and Sabato2007), tidal shear (Tessier and Terwindt, Reference Tessier and Terwindt1994; Greb and Archer, Reference Greb and Archer2007), tsunamis (Takashimizu and Masuda, Reference Takashimizu and Masuda2000; Nanayama et al. Reference Nanayama, Shigeno, Satake, Shimokawa, Koitabashi, Miyasaka and Ishii2000; Schnyder et al. Reference Schnyder, Baudin and Deconinck2005; Le Roux et al. Reference Le Roux, Nielsen, Kemnitz and Henriquez2008), wave-storm pounding (Molina et al. Reference Molina, Alfaro and Moretti1998) and seismic shocks (Sim, Reference Sims1975; Owen et al. Reference Owen, Moretti and Alfaro2011; Moretti and Ronchi, Reference Moretti and Ronchi2011; Sarkar et al. Reference Sarkar, Choudhuri, Banerjee, van Loon and Bose2014; He and Qiao, Reference He and Qiao2015). 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, Reference Owen1996). 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. Reference Owen, Moretti and Alfaro2011; Moretti and Ronchi, Reference Moretti and Ronchi2011; Sarkar et al. Reference Sarkar, Choudhuri, Banerjee, van Loon and Bose2014; He et al. Reference He, Quia, Jiao, Xu, Cai, Guo and Zhang2014; Chen et al. Reference Chen, Zhong, Ogg, van Loon, Chen, Yang, Liu and Xu2020; Guo et al. Reference Guo, He and Li2023).
4.c. 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. Reference Quia and Guo2013; He et al. Reference He, Quia, Jiao, Xu, Cai, Guo and Zhang2014, Reference He, Quia, Zhang, Tian, Cai, Chen and Zhang2015; Zhong et al. Reference Zhong, Jiang, Li, Su, Xu, Liang and Fan2022). The SSDS in the studied sediment sequence provide the intensity of palaeo-earthquake as well (Allen, Reference Allen1986; Rodríguez-Pascua et al. Reference Rodríguez-Pascua, Calvo, De Vicente and GómezGras2000). Some authors stated that intensity of an earthquake greater than 5 (M≥5) can produce liquefaction and fluidization (Atkinson, Reference Atkinson1984; Audemard and De Santis, Reference Audemard and de Santis1991; Galli, Reference Galli2000; Wheeler, Reference Wheeler, Ettensohn, Rast and Brett2002; Zhong, Reference Zhong2017; Zhong et al. Reference Zhong, Jiang, Li, Su, Xu, Liang and Fan2022; Guo et al. Reference Guo, He and Li2023). The computer modelling suggests M≈5 produce liquefaction feature within 4 km radius of the earthquake epicentre (Atkinson, Reference Atkinson1984), 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, Reference Marco and Agnon1995; Ambraseys, Reference Ambraseys1988; Obermeier, Reference Obermeier1996). 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. Reference Rodríguez-Pascua, Calvo, De Vicente and GómezGras2000). SSDS in four repetitive beds suggests the sediment succession may have affected four seismic pulses.
4.d. 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. Reference He, Quia, Jiao, Xu, Cai, Guo and Zhang2014; Sarkar et al. Reference Sarkar, Choudhuri, Banerjee, van Loon and Bose2014). 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).
Figure 7. Conceptual model: vertical section showing F.A-I passes upward to F.A-II with a seismite layer at the top of F.A-I (a, note fining up sequence and transgression-regression cycle), depositional scenarios for F.A-I (b, note the sea level-1), depositional scenarios for F.A-II (c, note sea level-2 provoked by seismicity induced subsidence).
4.e. 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. Reference Quia and Guo2013). 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, Reference Moretti and Ronchi2011; Sarkar et al. Reference Sarkar, Choudhuri, Banerjee, van Loon and Bose2014; He et al. Reference He, Quia, Jiao, Xu, Cai, Guo and Zhang2014). 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. Reference Chen, Zhong, Ogg, van Loon, Chen, Yang, Liu and Xu2020; Laborde-Casadaban et al. Reference Laborde-Casadaban, Homberg, Schnyder, Borderie and Raine2021; Mueller et al. Reference Mueller, Tamburelli, Menegoni, Perozzo, Amadori, Crispini, Federico, Seno and Maino2023). Considering the Middle Jurassic Tethys sea-floor deposition, we infer that the Middle Jurassic Gondwana breakup and synchronous active rifting (Lovecchio et al. Reference Lovecchioa, Rohaisb, Josephb, Bolattia and Ramos2020) 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.
5. Conclusions
Considering the SSDS distribution and sediment dynamics pattern at onset of the SSDS beds, we draw the following conclusion.
-
(a) 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.
-
(b) 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.
-
(c) 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.
-
(d) 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.
Acknowledgements
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.
Competing interests
The authors declare no conflict of interest.
Open access
