Introduction

The Late Devonian Frasnian/Famennian (F/F) boundary interval records one of the most important biotic crisis, interpreted to be one of the big five mass extinction events in Earth’s history (McLaren 1982; Raup and Sepkoski 1982; McGhee 1996; Stanley 2016). Numerous publications have shown that many Late Devonian fossil groups were impacted including trilobites, brachiopods, ostracods, ammonoids and tentaculitoids and even specific fossil communities such as coral reef ecosystems (Stearn 1987; Copper 2002; Ma et al. 2002; Olempska 2002; Racki 2005; Bond 2006; Wei et al. 2012; Crônier and François 2014; Ma et al. 2016; Mottequin and Poty 2016). The F/F interval is generally characterised by two dark, organic-rich lithologies that were dependent on the facies setting (Carmichael et al. 2019). These are known as the Lower and Upper Kellwasser events (LKE, UKE) that represent a stepwise ocean anoxia and extinction sequence (Walliser 1996; Gerecke and Schindler 2012). The section described herein lacks the diagnostic Kellwasser type lithologies. Here, the LKE took place during the Palmatolepis rhenana conodont zone of the late Frasnian, while the UKE is associated with the F/F boundary (Palmatolepis linguiformis conodont zone; Ziegler and Sandberg 1990). Even though numerous studies on more than 150 Kellwasser localities worldwide have been published, causes and mechanisms for the Late Devonian mass extinctions are still being debated. These events have been linked to an extraterrestrial impact (McGhee 1996; Wang et al. 1996), glaciation (Streel et al. 2000; Song et al. 2017), a cooling event that was linked with abrupt warming and sea-level rise (Thompson and Newton 1988; Algeo and Scheckler 1998; Riquier et al. 2005, 2006; Algeo and Scheckler 2010; Huang et al. 2018b), multiple arc volcanic events (Racki et al. 2018), and eutrophication and the development of suboxic to anoxic conditions present during both events as a result of high biological productivity (Algeo and Scheckler 2010; Whalen 2015; Carmichael and Kido 2014, Carmichael et al. 2016). Possible causes and triggers for anoxia and extinction of Late Devonian events were recently reviewed by Carmichael et al. (2019). These authors documented that the current knowledge concerning the F/F crisis is limited by a significant sampling bias, due to the fact that most previous studies sampled epicontinental seaways or passive continental shelves, primarily from localities across Europe and North America. Recently, more data on the F/F interval were published from southeast Asia (e.g. Huang and Gong 2016; Königshof et al. 2017; Huang et al. 2018a, b; Komatsu et al. 2018; Racki et al. 2019). Here, we report new biostratigraphical, sedimentological and magnetic susceptibility records across the F/F boundary interval exposed in the upper part of the Xom Nha Formation, Central Vietnam.

Geological setting

The F/F boundary is known from two sections in Vietnam. One section is exposed in the Si Phai Pass, Dong Van area (Ha Giang Province) in the northeast, close to the Chinese border (Tạ Hoa 2002, 2007; Königshof et al. 2017; Komatsu et al. 2018) and structurally is associated with the South China Block, within the Early Palaeozoic East Bac Bo Intracontinental Orogenic Belt. The other, the Xom Nha section in Central Vietnam (Figs. 1 and 2; Minh Hoa District, Quang Binh Province), belongs to the Indochina Block (Truong Son Belt, Viet Lao basin group; Thanh et al. 2006; Usuki et al. 2009; Trị et al. 2011). The Xom Nha section is exposed in a steep cliff at an entrance of a cave close to the small village of Xom Nha (N 17°48′26,9″, E 105°56′41,1″; Figs. 2 and 3).

Fig. 1
figure 1

Simplified map of structural units in Vietnam. The studied section belongs to the Truong Son Belt, Central Vietnam (small rectangle, more details are provided in the figure above)

Fig. 2
figure 2

Location and geological context of the studied section in Central Vietnam (Minh Hoa District, Quang Binh Province). Asterisk marks the position of the Xom Nha section (N 17°48′26,9″, E 105°56′41,1″; ± 3m)

Fig. 3
figure 3

Field photographs of the Xom Nha Formation bracketing the F/F interval: a limestones of the Len Don mountain and entrance to the cave, b taking samples with a rock saw at the steep cliff, c monotonous marl and limestone above the F/F boundary, d details of the F/F boundary within layer 18 (boundary marked in red): sampled beds are shown with yellow numbers; white numbers are from earlier field studies that were correlated with the new sampling in 2017. The F/F boundary is placed within bed 18 (18b contains the first conodont Palmatolepis triangularis)

The calcareous rocks exposed here belong to the upper part of the Xom Nha Formation and we have sampled a 2.6 m thick section which covers the F/F stage boundary interval. The Xom Nha section is mainly composed of fossiliferous, brecciated and laminated limestones, while shale beds occur in a subordinate number. The lithological characteristics of the F/F interval, such as black limestones and black shales known from many sections worldwide, are lacking as a result of facies, as described from other sections in Asia such as from Thailand (Königshof et al. 2012; Racki et al. 2019). The overall fauna is rather scarce in many layers except for conodonts, which occur in large numbers, particularly in the Frasnian part of the section. Conodonts from the sampled section range from the Frasnian Palmatolepis hassi Zone (Ziegler and Sandberg 1990; Ji and Ziegler 1993) to the Famennian Palmatolepis crepida Zone (Spalletta et al. 2017; Fig. 4).

Fig. 4
figure 4

a detail of the Frasnian/Famennian (F/F) boundary: microbioclastic wackestone is overlain by laminated limestone (vertical cross-section; numbers on the slab are a result of an older field campaign), no black shales or black limestone occur; b parallel-laminated lime mudstone overlying a wackestone; hardgrounds are characterised by a stepped escarpment on top of the hardground, forming a palaeorelief which is sometimes covered by a dark crust (manganese oxides) of variable thickness, sample VN-XN-18a, linguiformis Zone (equivalent of the Upper Kellwasser Limestone); c wackestone within the Palmatolepis rhenana rhenana Zone; bioclasts are preferably ostracods (some pelagic entomoceaceans with characteristic spines (upper left side) and cephalopods; sample VN-XN-15a, Pa. rhenana rhenana Zone; d microbioclastic wackestone with shell hash (mainly ostracods) and rare echinoderm ossicles, sample VN-XN-12, base of Pa. rhenana rhenana Zone; e bioclastic wackestone with coated bio- and lithoclasts; f wackestone with ostracod, cephalopod, trilobite and rare echinoderm and coral fragments, in sample VN-XN-11, Palmatolepis rhenana nasuta Zone; g beginning of a monomict stylobreccia as a result of pressure solution; h tectonic fault breccia, sample VN-XN-17a

Methods

Conodont biostratigraphy and microfacies analysis

Forty-six samples for conodont biostratigraphy from the 2.6 m thick section were processed. Conodont samples weighting 1 to 1.5 kg were collected from all beds. Some samples were barren or contained a very limited number of conodonts with a lack of zonal index taxa, so that a detailed stratigraphy for the middle Famennian was not possible. Limestones and marls were cleaned and crushed to about 2 cm3 and placed in wide-mesh sieves hung in buckets filled with warm tap water. Formic acid was added until the dilution reached a concentration of about 5% HCOOH. The process was repeated until samples were dissolved. The fractions (63 μm, 125 μm, 250 μm and 500 μm) were dried in an oven (~ 40°C) and the most promising conodont-sized fractions were separated by using liquid sodium polytungstate. Later, the conodonts were handpicked under a microscope and the most indicative species were photographed using an SEM. Given that it is common practice in high-resolution stratigraphic conodont studies, only Pa elements were identified. This is because many Late Devonian to Mississippian multielement reconstructions are still doubtful and incomplete.

Facies and microfacies types (FTs) were mainly compared on the basis of the classical models used by Wilson (1975) and Flügel (2004). In order to get a thorough overview on the section, we used thin sections that are at least 8 cm × 12 cm in size. These samples were usually collected at 10 cm intervals, but if necessary shorter intervals were used. Thin sections and conodont samples are stored at Senckenberg–Research Institute and Natural History Museum Frankfurt, Germany, under repository numbers SEM 70851 to SEM 70886 and SEM-XN-LD-01 to SEM-XN-LD-74, respectively.

Magnetic susceptibility measurements

From the 2.6 m thick Frasnian/Famennian Xom Nha section, 132 samples were collected for MS measurements from 24 layers (~2 cm/1 sample). Samples for magnetic susceptibility (MS) were measured in the Geomagnetism Department of the Institute of Geophysics laboratory at the Vietnam Academy of Sciences and Technology, Hanoi. The average distance between samples was 2 cm. All materials are “susceptible” to becoming magnetized in the presence of an external magnetic field. Low-field bulk, MS is an indicator of the strength of this transient magnetism (Ellwood and Gose 2006; Ellwood et al. 1995, 1997). MS is very different from remanent magnetism (RM), the intrinsic magnetization that accounts for the magnetostratigraphic polarity variations of materials. MS in stratigraphic sequences is generally considered to be an indicator of detrital iron-containing paramagnetic and ferrimagnetic grains, mainly clay minerals, and can be quickly and easily measured on small friable samples.

Results

Facies description and depositional environment

The sampled section (Fig. 3) is characterised by greyish well-bedded limestones, which exhibit four different microfacies (MF) types: lime-mudstone to wackestone, (MF 1), grey, fossiliferous wackestone to packstone (MF 2), brecciated grey limestone (MF 3) and alternating grey laminated limestone with marl layers (MF 4).

The overall facies setting points to outer shelf environments below fair-weather wave base. Beside conodonts, distinct layers are fossiliferous and contain tentaculitids, cephalopods, ostracods, trilobites and very rare crinoids and corals. Interestingly, the Xom Nha section does not exhibit characteristic black shales or black limestones around the F/F boundary (Fig. 4a) known from many places around the world, but instead show a similar lithology in comparison to other sections in Southeast Asia, particularly in Thailand and Vietnam (Königshof et al. 2012, 2017; Komatsu et al. 2018; Racki et al. 2019), which point to similar facies settings.

MF 1: microbioclastic lime mudstone to wackestone

Microfacies 1 occurs in Frasnian limestones and is dominated by greyish, thin-bedded lime mudstone to wackestone at the base of the section (Palmatolepis hassi Zone) and within the linguiformis Zone. The micritic matrix contains rare fossils, including conodonts, tentaculitids, benthic ostracods and small shell hash. In distinct layers, pelagic entomozoaceans also occur. The size of the bioclasts is small, burial solution occurs and some bioclasts are separated by clay seams. Shells are frequently dissolved and filled with sparry calcite. The fossil content generally varies between 5 and 10%. A mottled/nodular fabric occurs, along with styliolitization, but this is less frequent. Rarely, hard ground surfaces occur close to the F/F boundary (Fig. 4b). Similar sedimentological features around the F/F boundary are described from the Mae Sariang section in northwestern Thailand (Königshof et al. 2012).

Interpretation:

Hardgrounds are related to a combination of non-deposition or very low sedimentation rates. The overall fine-grained sediment and fossil content suggests a hemipelagic facies setting on an outer shelf environment with low sedimentation rates. The occurrence of pelagic entomozoaceans is restricted to distinct layers. Variable fabric ranging from thin-bedded to nodular/mottled can be explained as a result of bioturbation in a low-energy environment.

MF 2: fossiliferous wackestone to packstone

This rock type is composed of bioclastic wackestone to packstone which frequently occurs in the Palmatolepis rhenana rhenana and Pa. rhenana nasuta zones. Faunal elements comprise mainly conodonts, benthic ostracods, trilobites, tentaculitids, cephalopods and thin-shelled bivalves, which occur in descending order (Fig. 4c, d). Crinoid ossicles and small corals were observed, but occur with less frequency (Fig. 4e). Benthic ostracods are mainly disarticulated, but also complete ones occur. Pelagic ostracods are less frequent. The fossil content is higher than observed in MF 1 and can reach up to 30%. Single layers are packstone and bioturbation occurs in distinct layers. Fe–Mn crusts and microstromatolites around bioclasts and lithoclasts are observed in a few thin sections (Fig. 4f); the matrix is composed of lime mud. As a result of transport, the limestones show bedding and gradation, along with small-sized shall hash of up to 3 mm.

Interpretation:

This microfacies suggests a hemipelagic setting on an outer shelf environment with low sedimentation rates. Coated clasts can be interpreted as reworked hardground fragments (Zatoń et al. 2012; Lazar et al. 2013). Shallow-water bioclasts, such as corals, occur in distinct layers, and are very scarce. They occur in samples VN-XN-11 and VN-XN-12 (upper Pa. rhenana nasuta Zone/Pa. rhenana rhenana Zone), which appear to be linked with a sea-level fall. The increasing number of Icriodus specimens in sample VN-XN-11 may support this assumption.

MF 3: brecciated grey limestone

This facies type occurs in distinct layers within the Xom Nha section, and fauna is rare or absent. These limestones are characterised as a monomict stylobreccia (Fig. 4g). Another type of carbonate breccias occurs in the linguiformis Zone (sample VN-XN-17a) and exhibits strong fracturing of individual clasts as well as the whole rock. Microfractures differ in width and are filled with sparry calcite (Fig. 4h).

Interpretation:

Whereas the monomict stylobreccia is a result of pressure solution of lime mudstones, the second rock type is interpreted as a tectonic fault breccia.

MF 4: alternating grey limestone with marl layers

This microfacies occurs mainly in the Famennian part of the Xom Nha section and shows an alternation of laminated light-grey limestone and intercalated greyish marls. The marls are less frequent. These limestones and marls yielded rare fauna following the Upper Kellwasser Event. Occasionally, thin shell concentrations occur in the limestone. Limestones have a thickness of several cm and show parallel bedding as well as graded bedding with a relatively sharp lower erosional surface. These limestones are covered by thin-bedded marls.

Interpretation:

The fossiliferous limestones are interpreted to be deposited, in part, by distal turbidity currents and are typical of sediments found in basinal facies settings, such as lower slope facies sediments (Stow et al. 1996; Flügel 2004; Mulder 2011).

Conodont biostratigraphy and biofacies

Conodont biostratigraphy (Table 1, Figs. 5, 6, 7 and 8)

The results from the conodont stratigraphy in this section indicate that the upper part of the Xom Nha Formation ranges from the Palmatolepis hassi Zone to the Palmatolepis crepida Zone (Table 1). The crepida Zone is not shown in Figure 5, but this result is based on the earlier sampling (Tạ Hoa 2002). Zonation scheme follows Ziegler and Sandberg (1990), and Ji and Ziegler (1993) for the Frasnian conodont stratigraphy, and Hartenfels (2011) and Spalletta et al. (2017) for the Famennian conodont stratigraphy. Although the Montagne Noire Zonation (MN; Klapper 1989; Girard et al. 2005) provides a more detailed stratigraphy in comparison to the Frasnian Standard Zonation (Ziegler and Sandberg 1990; Ji and Ziegler 1993), we used the latter for the Frasnian because the conodont collection did not yield all naming MN zonal species.

Table 1 Conodont stratigraphy and conodont occurrence of the Xom Nha section in each sample. Biostratigraphically important conodonts are shown as well as other conodonts with long stratigraphical ranges
Fig. 5
figure 5

Important conodonts of the Xom Nha section (scale bares = 100 μm): a Palmatolepis hassi (sample VN-XN-2); b Ancyrodella nodosa (sample VN-XN-2); c Ancyrodella nodosa (sample VN-XN-2); d Ancyroides cf. leonis (broken specimen; sample VN-XN-3); e Ancyrodella nodosa (sample VN-XN-2); f Ancyrognathus triangularis (broken specimen, sample VN-XN-3); g Ancyrognathus primus (broken specimen, sample VN-XN-5); h Ancyrognathus triangularis (sample VN-XN-6); i Ancyrognathus triangularis (lower view, sample VN-XN-6); j Palmatolepis proversa (sample VN-XN-6); k Ancyrognathus aff. asymmetricus (sample VN-XN-6); l Polygnathus decorosus (sample VN-XN-8); m Ancyrodella nodosa (sample VN-XN-8); n Ancyrodella ioides (sample VN-XN-8); o Palmatolepis gigas gigas (sample VN-XN-6); Ancyrodella nodosa (sample VN-XN-8)

Fig. 6
figure 6

Important conodonts of the Xom Nha section (scale bares = 100 μm): a Ancyrodella ioides (sample VN-XN-9); b Polygnathus decorosus (sample VN-XN-11); c Ancyrognathus seddoni (sample VN-XN-11); d Palmatolepis unicornis (sample VN-XN-11); e Palmatolepis cf. foliacea (sample VN-XN-11); f Palmatolepis cf. hassi (sample VN-XN-11); g Palmatolepis subrecta (sample VN-XN-13); h Ancyrognathus seddoni (sample VN-XN-12); i Palmatolepis rhenana rhenana (sample VN-XN-12); j Ancyrodella nodosa (sample VN-XN-13); k Palmatolepis subrecta (sample VN-XN-14); l Palmatolepis eureka (sample VN-XN-13); m Palmatolepis juntianensis (sample VN-XN-13); n Palmatolepis juntianensis (sample VN-XN-14); o Palmatolepis subrecta (sample VN-XN-15); p Palmatolepis ederi (sample VN-XN-15)

Fig. 7
figure 7

Important conodonts of the Xom Nha section (scale bares = 100 μm): a Polygnathus tenellus (sample VN-XN-15); b Ancyrognathus seddoni (sample VN-XN-16b); c Palmatolepis linguiformis (sample VN-XN-16b); d Palmatolepis linguiformis (sample VN-XN-17a); e Ancyrodella curvata (sample VN-XN-17a); f Polygnathus webbi (sample VN-XN-17a); g Palmatolepis linguiformis (sample VN-XN-18a); h Ancyrodella curvata (sample VN-XN-18a); i Icriodus alternatus alternatus (sample VN-XN-18a); j Palmatolepis triangularis (sample VN-XN-18b); k Icriodus alternatus alternatus (sample VN-XN-18b); l Polygnathus brevilaminus (sample VN-XN-18b); m Icriodus alternatus alternatus (sample VN-XN-19); n Icriodus alternatus helmsi (sample VN-XN-19); o Palmatolepis triangularis (sample VN-XN-19); p Ancyrognathus sinelaminus (sample VN-XN-21c)

Fig. 8
figure 8

Important conodonts of the Xom Nha section (scale bares = 100 μm): a Palmatolepis cf. delicatula delicatula (sample VN-XN-21c); b Palmatolepis protorhomboidea (sample VN-XN-21c); c Polygnathus brevilaminus (sample VN-XN-21c); d Icriodus alternatus alternatus (sample VN-XN-21c); e Ancyrognathus sinelaminus (sample VN-XN-21c); f Icriodus alternatus helmsi (sample VN-XN-21c); g Palmatolepis crepida (sample X19, reillustrated from Ta Hoa (2002), figure II.4.18-9); h Palmatolepis quadrantinodosalobata (sample X20 reillustrated from Ta Hoa (2002), figure II.4.18-10)

Conodont zones were indicated by either the First Appearance Datum (FAD) or the presence of distinct taxa. The distribution of important conodont taxa is shown in Table 1 and most zonal indicative species are provided in Figs 5, 6, 7 and 8. Overall, the conodonts are well preserved. The conodont alteration index (CAI, Epstein et al. 1977) has a range in the Xom Nha section between CAI 4.5 and CAI 5.

The Late Devonian Frasnian and Famennian conodonts from the Xom Nha section in Central Vietnam contain mainly cosmopolitan species. Most sampled beds yielded conodonts (ranging from <10 conodont elements per sample to >300 elements per sample). The conodont record in the Frasnian is generally very good whereas the record in the Famennian is scarce, which is not a result of different sample sizes because conodont samples weighting 1 to 1.5 kg were collected from all beds. The base of the Xom Nha section starts in the Palmatolepis hassi Zone. In sample VN-XN-1, the conodont collection is dominated by Palmatolepis hassi and Pa. punctata. Other important species are Ancyrodella lobata, and Ancyrodella buckeyensis, which have their first occurrence in the early part of the Palmatolepis hassi Zone (Ji and Ziegler 1993). Sample VN-XN-3 did not yield many conodonts, but yielded some important conodonts, including Ancyrognathus triangularis and Ancyrodella nodosa. The first appearance of Ancyrognathus triangularis defines the next younger Ancyrognathus triangularis Zone. The entry of Ancyrodella nodosa was used by Ji and Ziegler (1993) as an important species, which has its first appearance close to the beginning of the Late hassi Zone (= Ancyrognathus triangularis Zone). The base of the Ancyrognathus triangularis conodont zone can be correlated to the base of the MN 8 zonation defined by the FAD of Palmatolepis proversa (Klapper 1989). Palmatolepis proversa was not found in sample VN-XN-3, but did appear later in the section in sample VN-XN-6. Whether the later appearance of this species is linked to sampling bias or the range of this species is diachronous in respect to conodont records in other localities remains questionable. As shown by Dzik et al. (2018), who investigated the F/F boundary in northeastern Vietnam, it seems obvious that even very fast or subtle changes in environmental conditions find expressions in the fossil record of species. Overlying conodont samples (VN-XN-5, VN-XN-6) yielded species that represent diagnostic species to identify this zone, including Ancyrognathus triangularis, and Palmatolepis proversa. Other important associated conodonts include Palmatolepis hassi, Ancyrognathus primus, Ancyrodella nodosa, Ancyrodella curvata and Ancyrodella lobata (Table 1). The next younger conodont zone in sample VN-XN-7 is defined by the FAD of Palmatolepis jamieae and Ancyrodella ioides. The latter appears at or close to the start of the Palmatolepis jamieae Zone (Ziegler and Sandberg 1990). The overlying sample VN-XN-8 yielded conodonts which correspond to the Palmatolepis rhenana nasuta Zone, which is defined by the FAD of the name-giving species Palmatolepis rhenana nasuta. This conodont zone ranges from sample VN-XN-8 to VN-XN-11. Important associated conodonts which were found in these samples are Palmatolepis hassi, Pa. jamieae, Pa. foliacea, Pa. plana, Pa. unicornis, Ancyrodella nodosa, Anc. lobata, Anc. buckeyensis, Ancyrognathus triangularis, Polygnathus decorosus and Po. webbi among others. Palmatolepis feisti which defines the MN 11 Zone (Klapper 1989) was not found in the Xom Nha section. Icriodontids are scarce (Table 1).

The next younger Palmatolepis rhenana rhenana Zone ranges from sample VN-XN-12 to VN-XN 16a (Table 1). The base of this interval is well defined by the FAD of Palmatolepis rhenana rhenana. In addition, Palmatolepis subrecta occurred for the first time in that sample. Palmatolepis juntianensis occurred a little bit later in this zone in the Xom Nha section, and was found in sample VN-XN-13. Palmatolepis cf. bogartensis, which represents the MN 13a Zone (Klapper 1989), shows up in sample VN-XN-15b, below the Palmatolepis linguiformis Zone. The next younger Palmatolepis linguiformis Zone is clearly defined by the FAD of Palmatolepis linguiformis in sample VN-XN-16b (Table 1) where a large number of individuals of this species were found. It seems likely that in some areas in Asia this species is common. Savage (2013) also reported a relatively large number of Pa. linguiformis from the Mae Sariang section in northwestern Thailand. The base of this conodont zone (Ziegler and Sandberg 1990; Ji and Ziegler 1993) corresponds to the MN Zone 13b of Girard et al. (2005). A further subdivision, suggested by Girard et al. (2005) (MN 13c = Ancyrognathus ubiquitus), was not possible to define in our section due to a lack of this marker species.

The Frasnian/Famennian boundary of the Xom Nha section corresponds to the base of sample VN-XN-18b (Table 1) and is based on the FAD of Palmatolepis triangularis (Ziegler and Sandberg 1990). In our section, we could not define the Palmatolepis subperlobata Zone which corresponds to the lowest part of the former Lower triangularis Zone (Hartenfels 2011; Spalletta et al. 2017) defined by the FAD of Palmatolepis subperlobata. This species does not occur in our sample VN-XN-18b, probably as a result of the scarcity of conodonts at that level. The conodont fauna shows low diversity, and Icriodus alternatus helmsi and Polygnathus brevilaminus were conodonts associated with this sample. The triangularis Zone (Hartenfels 2011; Spalletta et al. 2017) ranges from sample VN-XN-18a to sample VN-XN-21a. The next younger section corresponds to the Middle triangularis Zone of Ziegler and Sandberg (1990) and can be correlated with the Palmatolepis delicatula platys Zone (Spalletta et al. 2017). The useful conodont species for the zonal identification in the absence of the marker species is Ancyrognathus sinelaminus which was found in sample VN-XN-21c (Figs. 7p and 8e). This species enters at the base of this zone (Spalletta et al. 2017). Palmatolepis delicatula platys was not found which may support the observation by Schülke (1995) who reported that this species is rather rare outside Europe. Savage (2013) reported a small number of this species from a F/F boundary section in Thailand but this species was also described from the Toc Tat Formation in Northern Vietnam by Komatsu et al. (2018). The lack of this species in our section might be a result of facies or insufficient sample size. Other associated conodonts in this sample are Palmatolepis protorhomboidea, Pa. triangularis, Icriodus alternatus helmsi, Icr. alternatus alternatus and Polygnathus webbi. Stratigraphically, the section reaches up into the Palmatolepis crepida Zone, which is not shown herein, but was reported by Ta Hoa (2007).

Conodont biofacies

Based on the relatively abundant conodont assemblage of the Xom Nha section, conodont biofacies (e.g. Sandberg 1976; Weddige and Ziegler 1979; Sandberg and Ziegler 1979; Sandberg and Dreesen 1984; Sandberg et al. 1992) have been analysed. The whole section is mainly composed of lime mudstones and wackestones, suggesting an outer shelf environment below fair-weather wave base. Basically, the Xom Nha section is dominated by species of Palmatolepis and Polygnathus, which contain over 80% of the total P1 elements. Therefore, Palmatolepis elements exceed Polygnathus elements. Minor changes occur in some conodont zones (Table 1). The Frasnian part is generally dominated by a Palmatolepid-Polygnathid biofacies. In the lower part of the Xom Nha section, comprising the Palmatolepis hassi and Ancyrognathus triangularis zones, Palmatolepis species are dominant and comprise over 90% of the P1 elements. Ancyrodella species occur in subordinate numbers and the specimens are less frequent. Icriodus and Polygnathus specimens are present, although never numerous. In the uppermost part of the Ancyrognathus triangularis Zone, the number of Polygnathus increases and continues through the Palmatolepis rhenana nasuta Zone that is associated with an increasing number of Ancyrodella specimens. At the base of the Palmatolepis rhenana rhenana Zone, conodont biofacies is dominated by Palmatolepis and Ancyrodella. Polygnathus species occur but are very scarce. A similar biofacies occurs at the top of this conodont zone (Fig. 9).

Fig. 9
figure 9

Conodont Stratigraphy and conodont biofacies of the Xom Nha section. Pa. = Palmatolepis, Po. = Polygnathus, Anc. = Ancyrognathus, Ad. = Ancyrodella, and Icr. = Icriodus species

Within the Palmatolepis linguiformis Zone, the number of Polygnathus conodonts increases, and Icriodus specimens occur in small numbers below the F/F boundary. Icriodus specimens are present in some layers in the Frasnian, but they occur in very small numbers. In the Famennian, there is a marked increase of Icriodus and Polygnathus specimens, whereas the number of Palmatolepis specimens is observed to be decreasing. However, it should be noted that the total abundance of conodonts in the Famennian is scarce, particularly above sample VN-XN-22b. Thus, conodont biofacies analysis above this sample is not possible based on the present conodont record.

Magnetic susceptibility

From the 2.6 m thick Frasnian/Famennian Xom Nha section, 132 samples were collected for MS measurements and the MS data are given in Figure 10. The following interpretation convention for MS measurements is used. If the MS cyclic trends increase or decrease by a factor of two or more, and if the change is represented by two or more data points, then this change is assumed to be significant. This method is best employed when high-resolution data sets are being analysed (large numbers of closely spaced samples) and helps resolve variations associated with anomalous samples. Such variations may be due to weathering effects, secondary alteration and metamorphism. Longer-term trends are due to factors such as plate-driven eustacy, as opposed to shorter-term climate cycles, or event sequences such as impacts (Ellwood et al. 2004). In addition, variations in detrital input between localities, or a change in detrital sediment source, can often be compared between different localities. Abrupt MS shifts, such as that given at 2.0 m in Figure 10, may be interpreted as representing an unconformity, while cyclic trends are often interpreted as due to climate cycles, warm representing higher MS values due to erosion, while lows are generally due to cooler times with less rainfall and erosion. The Multi-Taper Method (MTM) uses MS data for time-series analysis (Fig. 11). The approach used here is to (a) collect closely and uniformly spaced samples in the field, (b) report MTM after normal detrending, (c) apply confidence limits to the MTM data, (d) place the important Milankovitch climate bands on the MTM diagrams for consistency checks and (e) establish a uniform model which can then be graphically compared to the MS cyclicity presented in Figure 10.

Fig. 10
figure 10

Magnetic susceptibility data from the F/F boundary interval of the Xom Nha section, Vietnam (red line marks the age of the F/F boundary at 372.2 Ma—see www.stratigraphy.org.)

Fig. 11
figure 11

MTM results show 2 major Milankovitch cycles exhibited in the MS data set: E2 (405 kyrs) and O2 (42 kyrs)

The MTM results show Milankovitch cycles with significant peaks. Using the work of Berger et al. (1992), we calculated two important climate cycles in the data set reported here, an eccentricity E2 (405 kyrs) cycle and an obliquity O2 (42 kyrs) cycle. Timing for the Xom Nha section is based on using the time-series analysis presented here (Fig. 11). Using the timing presented in the International Chronostratigraphic Chart, 2018 version with 372.2 Ma given for the F/F boundary (www.stratigraphy.org), we calculated the ages for the E2 and O2 main cycles for this section. From these results, we then calculated the interval of time represented by the MS samples in Figure 10 as 1.215 Ma. Therefore, the sediment accumulation rate (SAR) = 0.217 cm/kyrs for the collected section (Fig. 10).

The spectral power for MS data sets presented here was obtained using the MTM method after the data were detrended and subjected to a Hanning window to reduce spectral leakage and increase the dynamic range in the data (Jenkins and Watts 1968; Thompson 1982). Incidences of high-confidence peaks at the 90, 95 and 99% confidence limits in the resulting spectra were determined by employing the MTM analysis (Ghil et al. 2002), as calculated using the SSA-MTM toolkit (Dettinger et al. 1995). A null hypothesis of red noise was assumed (low frequency and high power in the spectrum, sloping towards lower values at high frequencies), using a three-taper model. As this method is prone to producing false positives, the use of statistical significance is limited here to its role in supporting (or not) the positions of multiple Milankovitch bands within the MTM dataset. The positions of these bands are mathematically fixed relative to each other, and thus a climate-forcing mechanism is supported by spectral analysis when the Milankovitch frequencies are also frequencies of high spectral power, as is the case here.

Conclusions

The overall microfacies and faunal characteristics shown here point to relatively deep and quiet water in the aphotic zone below storm wave base or current activity in an hemipelagic palaeoenvironment. Interestingly, the Xom Nha section does not show the characteristic black shale/limestone lithology around the F/F boundary, which is known from other hemipelagic sections around the world (Carmichael 2019, and references therein). The occurrence of discontinuous surfaces is related to a combination of non-deposition or low sedimentation rates. The main sedimentary process is settling from suspended matter. Based on the data presented herein, a sediment accumulation rate of 0.217 cm/kyrs was calculated. The conodont assemblage through the F/F boundary shows mainly cosmopolitan species, which is similar to many sections elsewhere. However, the section measured here differs from most sections in Europe, because the Xom Nha section in Central Vietnam contains relatively high numbers of the species Palmatolepis linguiformis, which has also been described from sections in Thailand. Deposition through the F/F boundary interval reported here for the Xom Nha section appears to have occurred over a period of ~1.2 Ma, and exhibits well-defined climate cyclicity.