Characteristics and evolution of inertinite abundance and atmospheric pO2 during China’s coal-forming periods

Dong-Dong Wang,Lu-Sheng Yin,Long-Yi Shao,Da-Wei Lyu,Hai-Yan Liu,Shuai Wang and Guo-Qi Dong

Abstract Coal, especially the inertinite in it, is highly sensitive to climate changes, showing an obvious response to paleoclimate conditions, in particular, to paleo-oxygen concentration (pO2). In this study, the inertinite abundance data of typical coal-forming periods in China were systematically collected and analyzed. Its characteristics and control factors were studied, and its evolution was established. Based on inertinite abundance data, pO2 evolution curves of various coal-forming periods in China were established, which fluctuated between 15% and 30% during the entire Phanerozoic. The inertinite abundance in coal deposits during Paleozoic in China was basically consistent with that of other areas of the world, while it was quite different globally from the Mesozoic to the Cenozoic. The results show that the inertinite abundance in coal deposits is controlled by pO2 and other factors including climatic zones, plant differentiation, sedimentary environments, and tectonic activities. The inertinite abundance in coal deposits in China during the Jurassic was high, suggesting dry paleoclimate of inland China.

Keywords: Inertinite abundance, pO2,Paleoclimate,Sedimentary environment, Coal-forming period, Comparative

1 Introduction

Coal is considered to be highly sensitive to climate changes. The maceral types and composition reflect paleoclimatic conditions during coal forming periods(Stach 1990), particularly the inertinite abundance (Diessel 2010). Inertinite is formed by fusainization of coalforming materials during peatification stage, and generally reflects relatively dry conditions (Stach 1990; Teichmüller 1989). At present, there are many disagreements on the genesis of inertinite, related to oxidation (Corrêada-Silva and Wolf 1980; Falcon 1989; Hunt 1989; Taylor et al. 1989; Silva and Kalkreuth 2005; Silva et al. 2008;Van Niekerk et al. 2010; Hower et al. 2011; Richardson et al. 2012; Wang et al. 2016) and fire (Bustin 1997;Scott 2000; Scott and Glasspool 2007; Diessel 2010;Hudspith et al. 2012; Shao et al. 2017). Diessel (2010)examined the overall distributions of global stratigraphy,and found that incomplete combustion was the main source of inertinite in coal. The genesis of inertinite requires a high oxygen content in atmosphere during coalforming periods. Therefore, relatively dry oxidation conditions would be favorable for inertinite forming.

In accordance with the findings presented by Scott(2000, 2002), Scott et al. (2000), and Scott and Glasspool(2006, 2007), a correlation exists between the inertinite content in coal and the atmospheric oxygen levels at the time of formation. The forming of inertinite is affected by atmospheric paleo-oxygen concentrations (pO2) during the coal-forming periods (Robinson 1989; Scott and Glasspool 2007; Diessel 2010; Sen et al. 2016). At a global scale, the pO2is considered to be the main factor controlling the inertinite abundance in coal. And local areas, paleoclimate conditions may also lead to different inertinite abundance (Diessel 2010). Therefore, based on the correlations between inertinite and pO2, researchers have speculated that the inertinite abundance in coal since the Paleozoic can be used to quantify the pO2during coal-forming periods. Consequently, a quantitative calculation formula was put forward (Robinson 1991;Wildman et al. 2004; Scott and Glasspool 2007; Glasspool and Scott 2010; Glasspool et al. 2015).

At present, researches involved in characteristics and evolution of inertinite abundance and pO2during coalforming periods in geological history are relatively less,and no unified understanding has been acquired globally. China, with a large area, experienced many coalforming periods over a long span of time, and therefore,is favorable for case studies of coal formation. In this study, the inertinite abundance data in coal during each coal-forming period in China were systematically collected and analyzed. Characteristics of the inertinite were analyzed and evolution pathways were successfully determined, and the evolution curves of pO2were established as well. Comparative studies have been carried out around the world to further enrich basic theories.We hope that our findings would help supplement and better understand the previous research results and in the global range to a certain extent.

2 Data sources and pO2 reconstruction method

2.1 Data sources

In this study, the inertinite abundance data during each coal-forming period were obtained from a wide range of sources and samples as follows: Atlas for Coal Petrography of China (Yang 1996); Coal Petrology of China

(Han 1996); published data in available academic resources; national coal resources potential evaluation report (Sun, 2013); parts of the provincial coal resources potential evaluation reports (Sun, 2013); and parts of the internal research reports, as well as portions of our own laboratory test results.

In order to avoid any error which could be potentially caused by different test conditions and descriptions of test data from different sources, the abundance data were first uniformly sorted and classified. For example, the inertinite abundance of coal was determined as the percentage of inertinite in total organic matter of coal. It is assumed that the coal is mineral-matter free. In order to ensure wide coverage and high representativeness of the data points, within the scope of small errors, the similar data points from the same region and during the same period were combined, and their mean values were selected. The data points with large changes in value during the same period and in the same area should be fully reflected as representatives.

In addition, statistical methods were used to systematically collate the inertinite abundance data from different sources during coal-forming periods, which include:mean abundance of inertinite; sample numbers; maximum values; minimum values; and standard deviations,as well as literary sources and regional distribution of data (Tables 1, 2, 3, 4, 5 and 6). The stratigraphic divisions and chronology framework used in this study referenced the 2014 edition of the China Stratigraphic Chart proposed in the attached table of the All China Commission of Stratigraphy (2015).

2.2 Reconstruction method of pO2 during coal-forming periods

There are only a few studies about calculating pO2from inertinite abundance in coal. In this study, the mosthighly-recognized pO2reconstruction model and an exponential function formula (Glasspool and Scott 2010)were applied, based on the following known conditions:

1) When pO2≤15%, charcoal is basically not produced (which means that the inertinite abundance is close to 0%) (Belcher and McElwain 2008). Under such conditions, in a typical coal rock,only one charcoal particle out of ≥500 counts, was found, and therefore, the inertinite abundance was considered to be 0.2% (Taylor et al. 1998).

2) The oxygen concentration in the modern atmosphere is approximately 21%. Based on a series of peat tests produced under different ecological,climatic, and geographical conditions, the inertinite abundance has been determined to be approximately 4.3% (Glasspool and Scott 2010).

3) Several studies have determined the upper limit of Phanerozoic pO2and its corresponding inertinite.Combustion experiments by Belcher and McElwain(2008) revealed that the upper limits of fully-burnt oxygen range between 25% and 35%. Therefore,pO2with a median value of 30% was taken as the upper limit of fully-burnt oxygen (Glasspool et al.2015). Lenton and Watson (2000) determined that an inertinite abundance of 44.4% was related to a pO2of 30%.

Therefore, based on the above-mentioned information, this study estimated and calibrated the pO2based on inertinite abundance as follows: 15% pO2corresponds to 0.2% inertinite abundance; 21% pO2corresponds to 4.3% inertinite abundance; and 30%pO2corresponds to 44.4% inertinite abundance. Subsequently, an exponential function of pO2and inertinite abundance was established using data processing software, with a high correlation coefficient R2value(R2=0.99). The equation is as follows (Glasspool and Scott 2010):

Table 1 The inertinite abundance and pO2 characteristics in Middle Devonian coals in China. IAVG: Mean inertinite abundance; N:Samples number; R:Percentage range; S:Standard deviation

According to the above equation, the mean inertinite abundance during each coal-forming period could be converted into pO2. In this way, the pO2evolution curve during coal-forming periods of the entire Phanerozoic Eon was successfully reconstructed.This method supplements the research of global pO2evolution in geological history.

3 Inertinite abundance and pO2 characteristics during different coal-forming periods

Atmospheric pO2is closely related to the evolution of lignocellulosic plants (Mills et al. 2016), and plays an indispensable role in biochemical cycles of the Earth (Berner 1999). The emergence of terrestrial plants, along with the development of lignin at the end Silurian, facilitated the formation of large-scale peat and coal deposits(Wang 2002). Subsequently, inertinite-containing strata began to appear and continued to the present day (Premovic 2006). Therefore, the long-term changes in atmospheric pO2in geological history may potentially explain the distribution and abundance of inertinite in coal deposits(Diessel 2010).

Statistical tables of inertinite abundance and pO2during different coal-forming periods in China were completed in Tables 1, 2, 3, 4, 5 and 6, and the spatial distributions of related data points were plotted in Fig.1.

3.1 Devonian

During the Devonian, plants evolved systematically, and a large number of aquatic plants moved towards land.The diversity of plant systems, and their amazing abilities to reproduce and grow, laid the foundation for the formation of higher plants on land. It has been discovered that during the Late Silurian (Rimmer and Scott 2006) and Early Devonian (Glasspool et al. 2006), an inertinite in the form of charcoal appeared for the first time in the world.

The Devonian is considered to be one of the most important coal-forming periods in China, during which both southern and northern China were characterized by tropical climate conditions (Li and Jiang 2013). The Early Devonian plants were in the early stages of growth and were considered to be higher terrestrial plants. They grew in clusters around coastal wetland areas with limited coverage, and generally formed thin coal seams.During the Middle Devonian, sea level fluctuated frequently, and plants were further propagated in low-lying areas, such as coastal boughs and lagoons. Coal seams were also formed in those areas (Han et al. 1993). Tectonic movements and transgression in southern China had profound influence on the development of coal basins and the formation of coal deposits (Zhang 1995; Li et al. 2018). It has been determined that during the Devonian, coal-bearing strata in China were mainly distributed in the southwestern regions, which mainly developed during the Dongganglingian (Givetian) of Middle Devonian, such as the Haikou Formation in Panzhihua of Sichuan Province, and the Qujing Formation in Luquan of Yunnan Province.

Devonian coal seams are generally thin, ranging between 0.1 m and 0.3 m, characterized by the low abundance and limited distribution of inertinite. The majority inertinite abundance ranges between 0.1% and 3.0%(Table 1). Correspondingly, the overall pO2was not high. It was approximately 15% during Early Devonian,and increased significantly during Middle and Late Devonian, and reached approximately 20% by Late Devonian.

3.2 Carboniferous-Permian

During the Carboniferous-Permian period, the southern and northern China were under tropical climate conditions (Tabor and Poulsen 2008) with significantlyincreased diversity of terrestrial plants (mainly ferns and gymnosperms), leading to profound changes in inertinite type and abundance (Table 2).

3.2.1 Early carboniferous

During Early Carboniferous, coal-bearing strata were widely distributed in the southern China. Under humid and rainy tropical climates, plants as Calamites and Sublepidodendron formed dense forests (Yang 1996), which mainly distributed in coastal lowlands, and did not show obvious climatic differentiation and palaeogeographical isolation (Liu and Quan 1996). The preservation condition for lignin in the plants had been continuously improved, which was conducive to the formation and preservation of coal seams (Diessel 2010). In southern China, the coal seams were mainly developed in deltaic environments, followed by meandering river-lakeenvironments, alluvial fan-braided river environments,and tidal flat-coastal environments, respectively (Liu 1990; Zhang 1995; Zheng 2008).

Table 2 The inertinite abundance and pO2 characteristics in Carboniferous-Permian coals in China. IAVG: Mean inertinite abundance; N: Samples number; R: Percentage range; S:Standard deviation

Table 2 The inertinite abundance and pO2 characteristics in Carboniferous-Permian coals in China. IAVG: Mean inertinite abundance; N: Samples number; R:Percentage range;S:Standard deviation(Continued)

Table 2 The inertinite abundance and pO2 characteristics in Carboniferous-Permian coals in China. IAVG: Mean inertinite abundance; N: Samples number; R:Percentage range;S:Standard deviation(Continued)

Table 3 The inertinite abundance and pO2 characteristics in Late Triassic coal-bearing strata in China. IAVG: Mean inertinite abundance; N:Samples number; R: Percentage range; S: Standard deviation

The coal-bearing strata of the Early Carboniferous were mainly distributed in the Ceshui Formation of Dewuan (Serpukhovian) Stage in central Hunan Province; Simen Formation of Visean Stage in Guangxi Province; and Wanshoushan Formation of Visean Stage in Yunnan Province. In addition, it has also been determined that in northwestern China,such locations as Chouniugou Formation of Visean Stage in Jingyuan, Gansu were also conducive to coal deposits.

The mean inertinite abundance in coals of the southern China during Early Carboniferous ranged between 5.6% and 20.2%, with the largest variation (5.0%-22.3%)observed in Ceshui Formation in different areas of Hunan Province. In Datang area of Guangxi, the value ranged between 13.0% and 15.0%, with an average of 14.0%. The pO2during Early Carboniferous had increased compared to that of Devonian, and remained stable around 25%.

3.2.2 Late Carboniferous-Early Permian

From Late Carboniferous to Early Permian, coal-bearing strata mainly developed in the basins of North China,with humid climate conditions (Chang and Gao 1993;Tan 2017), and extended stable tectonic setting. Coalforming processes were active in the following formations: Taiyuan Formation of Xiaoyiaoan-Longlinian(Kasimovian-Artinskian) Stages; Shanxi Formation of Zisongian-Longlinian (Asselian-Artinskian) Stages;and Liangshan Formation of Luodianian-Xiangboan(Kungurian-early Roadian) Stages (which spanned the Lower and Middle Permian). (1) In Taiyuan Formation of Xiaoyiaoan Stage (Kasimovian-Gzhelian), the mean inertinite abundance in coals ranges between 19.1% and 58.8%, with relatively high values occurring in the Ordos Basin. In the area of Jiangsu-Shandong-Henan, this value ranges between 8.3% and 28.6%. (2)In Taiyuan Formation of Zisongian-Longlinian (Asselian-Artinskian) Stages, the mean inertinite abundance remains high with significant variations, with the highest value of 51.6% occurring in Junggar Basin,and the lowest of 8.5% occurring in northern Jiangsu Province. (3) In Shanxi Formation of Zisongian-Longlinian (Asselian-Artinskian) Stages, the inertinite abundance slightly decreased when compared to the values of previous periods, ranging between 15.9%and 45.0%, (4) The inertinite abundance decreased from west to east from Taiyuan Formation to Shanxi Formation (Han 1996). (5) In the LiangshanFormation of Luodianian Stage (early Kungurian), the inertinite abundance was relatively low, ranging between 6.4% and 10.0%.

Table 4The inertinite abundance andpO2 characteristics in Early-Middle Jurassic coal-bearing strata in China. IAVG: Mean inertinite abundance; N: Samples number; R:Percentagerange;S: Standard deviation

Table 4 The inertinite abundance and pO2 characteristics in Early-Middle Jurassic coal-bearing strata in China. IAVG: Mean inertinite abundance; N: Samples number; R:Percentage range;S: Standard deviation(Continued)

Table 5 The inertinite abundance and pO2 characteristics in Early Cretaceous coal-bearing strata in China. IAVG: Mean inertinite abundance; N: Samples number; R: Percentage range; S:Standard deviation

Table 6 The inertinite abundance and pO2 characteristics in Paleogene-Neogene coal-bearing strata in China. IAVG: Mean inertinite abundance; N: Samples number; R:Percentage range;S: Standard deviation

Fig. 1 Spatial distribution of the data points of inertinite abundance in coal deposits during different coal-forming periods in China. The map of China is modified after the Standard Map Service of the National Administration of Surveying, Mapping and Geoinformation of China (http://bzdt.ch.mnr.gov.cn/)(No. GS(2016)1569).E.: Early;M.: Middle;L.: Late;and the same follows

During Late Carboniferous-Early Permian, the pO2increased significantly compared with the Early Carboniferous data. In Late Carboniferous, it reached the maximum value of almost 30%. It was above 27%throughout the Early Permian with no significant decrease and began to decrease since the late Early Permian,dropping down to approximately 24%.

3.2.3 Middle Permian

During the Middle Permian, coal-bearing strata of China were mainly developed in large-scale delta environments(Li et al. 2018), including the southern regions of northern China, west of Henan Province, southern and northern Anhui Province, Xuzhou area of Jiangsu Province,with semi-humid climate (Chang and Gao 1993; Tan 2017). Coal was mainly distributed in the following formations of Middle Permian: (1) Lower Shihezi Formation of Xiangboan Stage (late Kungurian-early Roadian),with mean inertinite abundance ranging between 7.9%and 45.3%; (2) Upper Shihezi Formation of Lengwuan Stage (Capitanian), with mean inertinite abundance ranging between 23.0% and 42.1%, and the inertinite mainly comprised of macrinite; (3) Tongziyan Formation of Lengwuan Stage (Capitanian), which was mainly distributed in Guangdong and Fujian provinces. The mean inertinite abundance varied greatly in different regions, ranging from 6.5%to 28.5%;(4)Liangshan Formation of Luodianian-Xiangboan Stages (Kungurian), of which the coalbearing strata mainly developed during Early-Middle Permian.Coal seams were thickest in the Middle Permian,with the mean inertinite abundance ranging from 13.8%to 58.3%.

During the early Middle Permian, the pO2increased slightly, and remained at approximately 28% until the late Middle Permian.

3.2.4 Late Permian

During the Late Permian, coal-bearing strata in China were mainly distributed in the southern regions (Liu 1990; Zhang 1995), dominated by a semiarid climate at the early stage, and an arid climate at the late stage(Chang and Gao 1993; Tan 2017). It was also formed in an epicontinental marine environment with humid climate conditions in the early Late Permian (Golonka 2011). The distribution of coal-bearing strata was directly controlled by palaeogeographical environment,with the transgression and regression causing continuous migration of coal-rich belts (Li et al. 2018).

Coal-bearing strata were mainly developed in Longtan Formation of Wuchiapingian Stage, and in Changxing and Wangjiazhai formations of the Changhsingian Stage.(1) In Longtan Formation, the mean inertinite abundance varied among different regions: 20.3%-20.7% in eastern Jiangxi;13.0%-15.6%in southern Anhui,southern Jiangsu and Zhejiang; and 11.8%-21.1% in Yunnan, Guizhou,and Sichuan.(2)In Changxing Formation,the value was 2.4%-33.3% in Sichuan and Chongqing with large variations. (3) In Wangjiazhai Formation, the value was 23.7%-42.4%in Guizhou,with an average of 29.3%.

During the early Late Permian, pO2did not vary much compared to the Middle Permian data,and remained stable at approximately 26%.During the mid Late Permian, it decreased to approximately 20%, and then increased once again to approximately 28%in late Late Permian.

3.3 Triassic

During the Triassic, gymnosperms quickly emerged and became dominant. At a global scale, the Triassic was characterized by dry and arid climate, although Middle Triassic was under increased rainfall and occasional humid conditions (Preto et al. 2010). The mass extinction event at the end Permian resulted in a 10 Ma coalforming gap during the Triassic, when most regions were covered by widespread desert conditions (Kutzbach and Gallimore 1989). The coal was scarce during the Middle Triassic(Retallack et al. 1996).

Coal-bearing strata mainly developed during the Late Triassic in the humid, hot, and rainy southern regions in China. In particular, it occurred in the Xujiahe Formation of Peikucuo Stage (Norian-Rhaetian) in Sichuan,Yunnan, Hubei, Guangdong, and Guizhou (Zhang 1995;Shao et al. 2014). Among those, the strata in Sichuan had the best coal-bearing properties, under coal-forming environments including coastal plains, coastal lake-delta plains, and coastal delta plains (Lu et al. 2008). However,during the same period, northern China was mostly under arid or semi-arid climate. It was only at the end of Late Triassic that coal-bearing strata of river and lake facies were formed in Wayaobao Formation of the Ordos Basin (Tian et al. 2011).

The inertinite abundance varied between the northern and southern China. For example, in the Wayaobao Formation in the north-central Ordos Basin, it was about 24.8%, and was mainly semi-fusinite. However, in southern China,the coal-bearing strata almost covered the entire southwestern regions, with the mean inertinite abundance in coals ranging between 5.8% and 26.3%. In regard to the Maantang Formation of Yazhiliangian Stage (Carnian), the value was 7.9%-19.5% in Zixing(Hunan Province), and 23.2%-28.4% in Sichuan Province, respectively. It varied most in Anyuan Formation of Jiangxi Province, ranging from 1.1% to 21.2%(Table 3).

During the Triassic, the pO2recorded a declining trend, from 25% in the early Late Triassic, to 18% at the end Late Triassic, which was the lowest value of the entire Mesozoic Era.

3.4 Jurassic

During the Jurassic, gymnosperms like cycads, conifers,and ginkgo were extremely abundant. The coal-bearing strata mainly developed in the northern and northwestern regions in China during the Early and Middle Jurassic, under a subtropical-warm temperate humid climate(Huang and Hou 1988). The coal-bearing basins were generally dominated by large- and medium-sized inland lake basins; while the coal-forming environments were dominated by alluvial-lake delta systems, followed by lacustrine systems(Li et al. 2018).

The coal-bearing strata during the Early Jurassic were mainly distributed in Badaowan Formation of Yongfeng Stage (Hettangian-Sinemurian) in the Junggar Basin of Xinjiang; Xishanyao Formation of Liuhuanggou Stage(Pliensbachian-Toarcian) in the Tuha and Yili basins,Xinjiang; as well as Xiahuayuan Formation of Shihezi Stage (Pliensbachian-Toarcian) in Hebei and Shandong regions (Shao et al. 2009; Shi et al. 2011). During that period, inertinite abundance in coal deposits varied significantly (Table 4), with mean values of 3.3%-14.8% in Badaowan Formation of Junggar Basin, and 4.7%-49.0%in Xishanyao Formation of Yili Basin.

During the Middle Jurassic, coal was well accumulated in entire northern China (Wu et al. 2008; Qin et al.2009). Most of the coal-bearing strata were found in Yan’an Formation of Shihezi Stage (Pliensbachian-Toarcian) in the Ordos Basin. In Inner Mongolia, the mean inertinite abundance in coals ranged between 31.9% and 84.0%, with maximum values of up to 86.9% in Longdong Coalfield of Gansu Province. In Xishanyao Formation of Shihezi Stage (Pliensbachian-Toarcian)(Xinjiang), the mean inertinite abundance varied greatly,ranging between 12.2% and 68.8%. In Yaojie Formation of Shihezi Stage (Pliensbachian-Toarcian) in Tianzhu county (Gansu), the value was between 1.1% and 3.4%,which was the lowest values during that period. In Toutunhe Formation of Manasi Stage (Bathonian-Callovian)in Minhe Basin (Qinghai), the value was between 24.0%and 66.5%.

Fig. 2 Summary of inertinite abundance relative to the geological age, atmospheric paleo-oxygen concentration, paleotemperature, and carbon dioxide concentration fluctuations. a Geochronological distribution of global inertinite abundance (Diessel 2010). Dashed lines refer to the three repetitions of inertinite abundance from high to low; mf: Mineral-free; b Geochronological distribution of inertinite abundance in China. Four major cycles are identified (marked by red lines); c Curve of atmospheric paleo-oxygen concentration in previous studies (Berner 1999, 2006;Berner and Kothavala 2001; Ward et al. 2006); d Simulation curve of atmospheric paleo-oxygen concentration in this study. Dashed lines in c and d refer to the pO2 value of the modern atmosphere (approximately 21%); e Variations of the carbonate-hosted oxygen isotopes, and the calculated temperatures (Shaviv and Veizer 2003). Dashed line refers to the base line of oxygen isotope; f Changes in atmospheric carbon dioxide concentration. R: Ratio of atmospheric CO2 concentration in geological history compared to the modern atmospheric CO2 concentration) (Royer et al.2004)

pO2changed drastically during the Jurassic compared to earlier periods. Since the Early Jurassic, it began to rise and reached approximately 28% by the end Early Jurassic. During the Middle Jurassic, it rapidly increased to 30% or more, which was the highest value ever recorded and was maintained throughout the Late Jurassic.

Fig.3 a Variation curves of inertinite abundance; b Prediction curve of the Phanerozoic atmospheric paleo-oxygen concentration

3.5 Cretaceous

The Cretaceous was an important period of plant evolution. Angiosperms appeared and flourished, which provided important materials for coal production. During the Early Cretaceous, a series of continental fault basins were formed due to rifting and faulting activities in northeastern China, which provided sites for coal formation (Li et al. 1987). The coal-bearing strata were mainly formed in alluvial fans, fan deltas, lakeside deltas, and lacustrine environments. The coal-rich belts were mainly located on the sides of main faults near basin margins,with their distribution directions consistent with basin trends (Cai et al. 2011; Shao et al. 2013).

During the Early Cretaceous, coal-bearing strata were mainly distributed in northeastern China, e.g., on the western side of Greater Khingan Mountains and the southern margin of Songliao Basin. They occurred in the Muling Formation, Yixian Formation, and Jiufotang Formation of Rehei Stage (lower Aptian) in Heilongjiang;Fuxin Formation and Shahai Formation of Liaoxi Stage(upper Aptian) in Hebei; and the Huolinhe Formation of Jibei Stage (Berriasian-Hauterivian) and the Yimin Formation of Rehei Stage (Barremian) in Inner Mongolia.The mean inertinite abundance in coal is relatively higher in the west than in the east of northeastern China(Table 5). It varied greatly in some coalfields in western Inner Mongolia, e.g., ranging between 26.1% and 55.8%in Yimin Formation of Jalainur Coalfield, with 42.4% on average. However, this value was relatively low in coalfields in eastern Inner Mongolia. For example, it ranges between 2.2%-16.5% in Fuxin and Tiefa Coalfields of Liaoning Province, and 9.4%-25.0% in Jixi and Hegang Coalfields of Heilongjiang Province. During the Early Cretaceous, pO2was slightly lower than that during the Jurassic, which, however, still kept at a relatively high level of approximately 25%.

3.6 Paleogene to Neogene

Since the Cenozoic, evergreen and deciduous broadleafed plants of modern angiosperms have gradually flourished, which indicated the obvious seasonal climate changes, as well as the rapid development of modern plants which are more adapted to severe climatic conditions (Zhao et al. 1995).

From Paleogene to Neogene, coal-bearing strata mainly developed in small continental basins in the northeastern and southwestern China. The coal-forming processes mainly occurred in the lake and delta marsh environments (Zhang 1995;Li et al. 2018).

During the Paleogene, the coal-bearing basins were mainly distributed to the east of Greater Khingan Mountains and north of Qinling Mountains and southwest of Guangxi. The coal-forming periods include Eocene and Oligocene. (1) The Eocene coalbearing strata mainly included Hunchun and Shulan formations of Ashantou Stage (Ypresian-Lutetian) in northeastern China; Nadu Formation of Yuanquan Stage in Baise and Nanning basins in Guangxi; and Lijiaya Formation of Yuanquan (Bartonian) in Liangjia Coal Mine in Shandong. (2) The Oligocene coalbearing strata included the Pinghu Formation of late Caijiachongian Stage (Priabonian) within Xihu Depression of East China Sea, and the Yacheng Formation of Tabenbuluckian Stage (Chattian) in Qiongdongnan Basin. The inertinite abundance in coal was generally low (Table 6), with the mean values ranging between 1.1% and 11.4%. (3) During the Oligocene, a few coal seams developed in Huerjing Formation of Wulanbulagean Stage (Rupelian) in Erlianhot, Inner Mongolia, within which the inertinite abundance was very high, ranging from 23.8% to 38.7%, with a mean value of 30.4%, which was significantly higher than that in other regions during the same period.

During the Neogene, coal-bearing strata were mainly distributed in the Miocene and Pliocene coal basins in the eastern coastal and the southwestern areas of southern China (Zhang 1995), mainly in continental lacustrine basins. The Miocene coal-bearing strata included the Xiaolongtan Formation of Tonggurian Stage (Langhian-Serravallian) in Yunnan Province and the Nanzhuang Formation of Shanwangian Stage (Burdigalian) in the northwestern basins of Taiwan. The Pliocene coalbearing strata mainly included the Shagou Formation of Mazegouan Stage (Piacenzian) in Shaotong-Qujing Basin of Yunnan. During the Neogene, the total inertinite abundance was relatively low, with the mean value ranging between 0.2% and 10.2%.

During the Quaternary, coal seams were generally developed in Yuanma Formation of Nihewanian Stage(Gelasian-Calabrian) in Tengchong, Yunnan. The inertinite was mainly composed of filamentous bodies and fungi, and was of low abundance (mean value of 4.0%)(Jin and Qin 1989).

Throughout the Cenozoic, pO2fluctuated around 21%,which was roughly the same as that in the current atmosphere.

3.7 Inertinite abundance and pO2 evolution

During the entire Phanerozoic, the abundance and distribution of inertinite in coal varied greatly. In terms of time,there were four major cycles in which the inertinite abundance first increased and then decreased, i.e., during Early Devonian-Late Permian; Late Triassic-Early Jurassic; Middle Jurassic-Late Cretaceous; and Paleogene-Neogene, respectively (Figs. 2, 3a). During these periods, the evolution process of pO2also changed a lot (Figs. 2,3b).

3.7.1 Early Devonian-Late Permian

During the Paleozoic, the inertinite abundance in coal was at a very low level in the Early Devonian. During Middle Devonian, this value was as low as 0.1% in Qujing Formation of Dongganglingian Stage in Yunnan. From Early Carboniferous to Late Permian, accompanied by the unprecedented diversity occurring in the plant kingdom, this value increased rapidly and reached its peak of 58.8% in Taiyuan Formation of Xiaoyiaoan Stage in Junggar Coalfield, Inner Mongolia. It remained high throughout the Permian,e.g., in Changxing Formation of Changxingian Stage in Guizhou Province. Until Late Permian, it slowly decreased to 30.7%. From Devonian to Permian, the pO2kept increasing and reached a maximum value of approximately 29% during the Middle Permian.Throughout Permian-Triassic, the pO2values fluctuated several times. However, prior to the Early Triassic, the pO2values had rapidly decreased to 23%.

3.7.2 Late Triassic-Early Jurassic

During the Mesozoic, the inertinite abundance in coalbearing strata of China was quite different from that of the Paleozoic. During the early Late Triassic, it rapidly decreased to 1.1%, the lowest value of the entire Phanerozoic, in Anyuan Formation of Peikucuo Stage in Leping (Pingxiang, Jiangxi). Immediately after the global coal gap of the Early and Middle Triassic, the inertinite abundance of the Late Triassic increased rapidly until reaching the same level as before the Early Triassic, with the peak value of 28.4% in Sichuan. From Late Triassic to Early Jurassic, the value decreased rapidly, only 4.7%in Xishanyao Formation of Liuhuanggou Stage in Tuha Basin (Xinjiang). Correspondingly, during the early Late Triassic, pO2was at the lowest level of only 18%, which,however, increased rapidly since the Early Jurassic,reaching up to 26%by the end of Early Jurassic.

3.7.3 Middle Jurassic-Late Cretaceous

During the entire Jurassic, the inertinite abundance displayed an increasing trend. It was not high at the Early Jurassic,but increased rapidly since the early Middle Jurassic, and reached 84%, i.e., the maximum value of the entire geological history at Yan’an Formation of Shihezi Stage in Gansu Province. However, the value significantly decreased since the Early Cretaceous. By the Late Cretaceous, it reduced to 28.3% in Fuxin Formation of Liaoxi Stage in Wanquan Coalfield,Hebei Province.Correspondingly, from Middle Jurassic to Late Cretaceous,pO2also changed greatly, which first increased to over 30% during the Middle Jurassic, and then slowly decreased but still remained as over 25%.

3.7.4 Paleogene-Neogene

The inertinite abundance in coal-bearing strata decreased during the Cenozoic when compared with prior periods, but was relatively stable. However, this value increased since the beginning of the Paleocene.During Oligocene, it reached 30.4%, the maximum value of the entire period, in the Huerjing Formation of Wulanbulagean Stage in Erenhot, Inner Mongolia,and then decreased again. During Neogene, it decreased to 1.8% in Shagou Formation of Mazegouan Stage. During the entire Cenozoic, pO2did not change significantly, remaining between 20% and 22%,which was equivalent to the current oxygen levels in the atmosphere.

4 Genetic analysis and comparison of inertinite abundance and pO2 characteristics during different coal-forming periods

China has experienced many coal-forming periods over a long span of time, The data generated from the coal deposits are important for the global study on inertinite abundance and pO2characteristics in coal-bearing strata.In this study, by examining the characteristics of inertinite abundance in China’s coal-bearing strata during each coal-forming period, and comparing the results with worldwide researches conducted by Diessel (2010),it was found that the basic trends were the same (Fig. 2).For example, during the Early Paleozoic, coal-forming processes were weak, and the inertinite abundance in coal was relatively low. Then, during the Late Paleozoic,this abundance showed regular increases or decreases on a global scale during different coal-forming periods.These trends were manifested as the high inertinite abundance from Late Carboniferous to Early Permian.The Middle Permian values were lower, and there was a rapid increase again during the Late Permian. It is also determined that from the Mesozoic to the Cenozoic,and from the Early Triassic to the Middle Triassic, the inertinite abundance changed greatly. There were generally three repetitions of inertinite abundance form high to low, which covered the time span of: (1) from Late Triassic to Early Jurassic; (2) from Middle Jurassic to Late Cretaceous; and (3) from Paleogene to Neogene, respectively. However, differing from the global research results of Diessel (2010), the inertinite abundance of the Middle Jurassic coal in northern and northwestern China was very high, reaching the peaks of each of the coal-forming periods. The atmospheric pO2in each of the coal-forming periods, which was simulated by the inertinite abundance in coal, also displayed corresponding characteristics (Fig. 2).

4.1 Early Devonian-Late Permian

During the Paleozoic in China, the inertinite abundance in coal-bearing strata displayed a trend of alternating increase and decrease. This was particularly obvious in regard to the high inertinite abundance from Late Carboniferous to Early Permian, showing low values during Middle Permian and rapid increase again during Late Permian. The variations were consistent with those observed for the global coal deposits, indicating that the formation of inertinite in coal during that period was significantly influenced by global factors (for example,atmospheric paleo-oxygen concentrations), which went beyond the influence of regional environments and their specific flora (Diessel 2010).

During the Early Carboniferous, the existing plants mainly grew in coastal lowland areas, and the climatic differentiation and palaeogeographical isolation were not obvious. In Late Carboniferous, latitudinal zonation occurred in the plant kingdom in order to adapt to climatic differentiation (Liu and Quan 1996). In the Carboniferous, coal deposits were mainly distributed in the southern continental Laurasia-Russia and the western Baltic Plate under tropical environments. Coal deposits were also located in the northwestern regions of continental Gondwana in Africa; and the northern parts of China and the Siberian Plate of the northern temperate zones. The first few regions were tropical year-round moist plant zones, while the Siberian Plate belonged to the temperate plant zones (Li and Jiang 2013). During the Carboniferous, terrestrial plants were diverse and vascular systems rapidly developed. The lignin tissues of plants were prone to smoldering, which was conducive to forming inertinite (Diessel 2010). The climatic zoning and geographical division of the plants were still similar to those of the Carboniferous (Li and Jiang 2013). During the Late Permian, the Earth’s vegetation was rich in lignin and was more prone to smoldering than full combustion (Robinson 1989). The tendency of lignin to form coke (for example, inert components) in the pyrolysis process may have contributed greatly to the increases in the content of inertinite components at the end Permian(Diessel 2010). Furthermore, during the Late Permian,local volcanic eruptions in northern China triggered massive wildfires (Zhao et al. 2010), which also contributed to the increases of inertinite abundance in China’s coal deposits.

Fig. 4 Reconstruction map of global palaeogeographic environment during Middle-Late Jurassic (175-146 Ma) (Li and Jiang 2013). Note: The time node of plate reconstruction is 160 Ma. Orogenic belt (Golonka 2011); Shallow sea (Golonka 2011); Ocean currents (Scotese 2002);Paleoclimate(Sellwood and Valdes 2008)

There is a close relationship between the oxygen content in atmosphere and the metabolism of plants and animals (Robinson 1989). The distribution of oxygen within the atmosphere of an area is particularly related to the variations in abundance of inertinite (Berner 2006; Ward et al. 2006). The growth of forest vegetation,as well as plant assimilation and transpiration led to the removal of a large amount of carbon dioxide from the atmosphere (Diessel 2010). The carbon extracted from atmospheric carbon dioxide was photosynthesized into plants, and much of that became buried and stored in the form of fossil fuel. Marshall et al. (2020) found a large number of mutant plant sporopollen based on mass extinction events during the Devonian-Carboniferous. The variation of plant debris reflectivity distribution indicates a higher charcoal content, indicating that atmospheric oxygen does not decrease. During the Carboniferous-Permian coal-forming periods, the pO2changed greatly. The diversity of terrestrial plants increased significantly during the Early Carboniferous,resulting in a continuous increase in paleo-oxygen concentrations throughout the Carboniferous. Until the Mesozoic, the pO2remained at or above 26%, which was consistent with the predictions based on plant carbon isotopic fractionation (Lenton and Watson 2000). During the Middle Permian, the pO2first decreased, and then fluctuated. However, it continued to remain at high levels. The aforementioned fluctuation pattern of the pO2was very similar to the bimodal distribution model of pO2established in previous study (Bergman et al.2004). Diessel (2010) found a relationship between the high inertinite content, the high atmospheric oxygen and low atmospheric carbon dioxide levels in Late Paleozoic. However, the inertinite abundance was found to be more complex than the relatively undifferentiated oxygen curves and closer to the δO18paleotemperature model.

In previous comprehensive analyses, the global coalforming plant types, sedimentary environments, and the climatic zones of the Paleozoic were found to be relatively stable. The inertinite abundance in global coal deposits were remarkably consistent, and were mainly controlled by the changes of global atmospheric pO2.The differences in the abundance and dispersion degrees of inertinite in coal deposits located in different regions also reflected the fact that they had been influenced by such secondary factors as climate, deposition, and structure.

4.2 Late Triassic-Early Jurassic

The end-Permian mass extinction (Virgili 2008) halted coal formation worldwide for approximately 10 Ma. This is widely believed to result from the chain reaction caused by eruption of Siberian igneous provinces (Diessel 2010; Zhao et al. 2020). The coal seams appeared again in the Late Triassic with high inertinite abundance.After that period, the inertinite abundance displayed a trend of gradually decreasing until it reached its lowest point in the middle to late parts of the Early Jurassic.

The geographical distributions of global supercontinents determined the climatic and environmental characteristics of various regions during the Triassic. For example, from the Late Permian, the glaciers gradually began to disappear, and as the sea level rose, global warming occurred. The majority of Pangaea was located in arid and tropical zones, while the remainder was located in temperate zones. The Late Triassic coalforming regions were mainly distributed in the arid and tropical climatic zones located around the Northern Hemisphere’s Paleotethys Ocean (Li and Jiang 2013).From the Late Triassic to Early Jurassic, extensive volcanic activities in the central Atlantic magma province led to the release of large amounts of carbon dioxide and/or methane. This had resulted in extreme greenhouse effects, with wetter climatic conditions and increased storm and lightning activities (Petersen and Lindstr?m 2012). These storm and lightning events frequently led to a large number of fires (Belcher et al.2010), which produced high concentrations of charcoal(inertinite), resulting in a relatively high abundance of inertinite in the coal deposits.

Compared to the Permian, the records of Triassic charcoal are actually very sparse, reflecting the relatively low oxygen content in the atmosphere at that time. According to the simulations of inertinite abundance during the Late Triassic, the atmospheric pO2during that period had decreased from approximately 25% to 18% in China, globally even dropped below 13% (Berner 2002,2005). Paleoredox and phosphorus speciation data from the ocean depth profile of Svalbard were for this time provided by Schobben et al. (2020). Large areas of the ocean experienced hypoxia prior to the Permian-Triassic boundary and caused some species to disappear before the major extinction at the Permian-Triassic boundary (Huey and Ward 2005). A large amount of CO2was released into atmosphere by the explosions in Siberian large igneous provinces, which subsequently led to decrease in the pO2(Zhao et al. 2020). However, the terrestrial plants caused rise in pO2in the atmosphere significantly (Rimmer et al. 2015). The low pO2is also reflected by the slow recovery of land vegetation following the Permian-Triassic mass extinction event.

4.3 Middle Jurassic-Late Cretaceous

From the late Early Jurassic to the Late Cretaceous, the abundance of inertinite in China’s coal had displayed a trend of gradual decrease. Previous studies on a global scale suggested that the inertinite abundance in coal decreased from the Late Triassic to Late Jurassic (Diessel 2010). The coal forming process during the middle and late Late Jurassic in China was not favourable. During the Late Cretaceous, both the coal-forming intensity and the abundance of the inertinite in coal decreased gradually.

During the Jurassic, coal was mainly developed in the north and south latitudes 40° polar region. The climate at that time was dominated by temperate type, whereas frigid climate zone was distributed in the northern of Siberia. The development and distribution of coal seams had symbiotic relationships with the survival of plants at that time. The plants surviving near the equator were the most singular, with conifers and cycads mainly developing into scattered forests. In the middle latitudes,conifers, cycads, and ferns flourished with large numbers of species and genera. However, compared with the equatorial regions, the higher latitudes mainly included broad-leafed plants such as Ginkgo biloba. From a global perspective, only the coal-forming areas in northern and northwestern China were considered to be located in inland arid climatic zones. Meanwhile, the coal-forming areas in other regions were mainly located in the coastal areas of temperate zones, and a few in the coastal areas of cold zones (Li and Jiang 2013) (Fig. 4). For example,the Middle Jurassic Yan’an Formation in Ordos Basin was formed under relatively dry climatic conditions(Huang et al. 2010). The surfaces of peat bogs during coal-forming periods were under oxidative conditions for long periods of time (Zhang and Wu 1996; Chen et al. 2013). The simultaneous occurrences of charcoal fragments and high concentrations of pyrolytic paths in the Yan’an Formation coal seams of the Ordos Basin provided evidence of peatland wildfires. The cycles of peatland wildfires were caused by seasonal changes driven by the changes in rainfall. The low reflectivity of inertinites in all of the coal seams may have resulted from low-temperature surface peatland wildfires with relatively high atmospheric paleo-oxygen concentrations of approximately 29% (Zhang et al. 2020). The pO2in the Jurassic global atmosphere should have been consistent and did not vary from place to place. The differences in climatic zones, plant types, and sedimentary environments had significant impacts on the inertinite abundance in coal deposits. This may also have been an important reason why the inertinite abundance in the Jurassic coal in China was much higher than that in other parts of the world.

During the late Early Jurassic, the pO2increased rapidly (～28%); reached its maximum value (＞30%); and then gradually decreased (～25%). Previous studies suggested that the pO2in the atmosphere had decreased to approximately 12% during the Early and Middle Jurassic(Berner 1999, 2006; Berner and Kothavala 2001; Ward et al. 2006). In fact, some researchers have also questioned the effectiveness of such lower atmospheric pO2.Therefore, based on combustion experiments under the fully-controlled conditions of various materials, it is now believed that at least a 15% atmospheric pO2was required for inertinite to produce the Mesozoic wildfires(Belcher and McElwain 2008).

During the Cretaceous, the abundance of inertinite in coal deposits fluctuated, which was closely related to the climatic conditions of the same period. Also, during the Cretaceous, a high-temperature environment existed throughout the world, with the temperate zones in the Northern Hemisphere reaching 60° north latitude and in the Southern Hemisphere reaching 70° south latitude.Therefore, with intense volcanic activities around the world, the atmospheric carbon dioxide content reached the highest level in history, which exacerbated the greenhouse effects and played a key role in the prosperity of terrestrial animals and plants. The coal-forming regions of the world during Cretaceous were mainly located in the tropical and subtropical climatic zones between approximately 30° and 60° north latitudes in the Northern Hemisphere, and some were located within the cold zones where the coal mainly formed in coastal environments and subordinately in the inland areas (Li and Jiang 2013). The climate of the Cretaceous changed frequently, during which there were significant drought or rainfall events (Spicer 2003; Belcher and McElwain 2008;Belcher et al. 2010; Brown et al. 2012)which affected the abundance of inertinite in coal. Early Cretaceous anoxic event (Huang et al.2008)had profound effects on the atmosphere-ocean system, as well as the accompanying global cooling and the enhancement of the ocean oxidative capacity. Boudinot and Sepúlveda (2020) believed that the Mid-Cretaceous oceanic anoxic event would lead to a decrease in pCO2and an increase in pO2in the atmospheric and oceanic system, and trigger mountain fires, indicating that the oxygen content in the Mid-Cretaceous was not low. Atmospheric pO2was calculated by simulating the abundance of inertinite in coal during this period. The atmospheric pO2in the Early Cretaceous was high and displayed a gradually decreasing trend, which was significantly higher than the previous studies (Berner 1999, 2006; Berner and Kothavala 2001; Ward et al. 2006).

4.4 Paleogene-Neogene

From the beginning of the Paleogene to the end of the Neogene, the abundance of inertinite increased slightly after the Paleocene, but still remained at a low level. Although some researchers reported forest fires during that period, the heavy rainfall events during the early Oligocene may have mitigated the impacts.

5 Discussion

The inertinite abundance in global coal during the Paleozoic displayed a fluctuating trend, particularly during the Late Carboniferous to Early Permian. It was higher in the Middle Permian coal, lower in the Late Permian coal, and more concentrated in the global coal.These findings indicated that the formation of inertinite during the Paleozoic was significantly influenced by global factors (such as atmospheric paleooxygen concentrations) (Diessel 2010), which went beyond the influence of regional environmental conditions and specific flora(Teichmüller 1952; Snyman 1961; Falcon 1975).

From the Mesozoic to the Cenozoic, the abundance of inertinite in coal changed greatly following the coal discontinuity which occurred from the Early to Middle Triassic. The general trend was that the inertinite abundance in the coal deposits repeated three times respectively, from high to low. On a global scale, the three replications were the Late Triassic to Late Jurassic; Late Jurassic to Late Cretaceous; and Late Cretaceous to Holocene (Diessel 2010; Glasspool and Scott 2010). In China, the three repetitions comprised the Late Triassic to Early Jurassic; Middle Jurassic to Late Cretaceous; and Paleogene to Neogene. The difference lies in the fact that the abundance of inertinite during the Middle to Late Jurassic in China was much higher than that in other parts of the world (Wang et al. 2019).

Diessel (2010) conducted a systematic study regarding the stratigraphic characteristics of global inertinite. The research results generally supported the theory that the global atmospheric pO2was the main controlling factor for the abundance of inertinite in global coal. In addition, the influences of climate conditions, sedimentary environments, differential subsidence, and other factors were believed to have led to regional and local changes in the abundance of inertinite in coal. This study also supported those viewpoints as a whole, and suggests that the global atmospheric pO2was the main controlling factor of the abundance of inertinite in global coal. During Paleozoic, due to the relatively single and stable coal-forming climatic zones, plant types, and sedimentary environmental conditions, the inertinite abundance in coal deposits was very obvious in response to the atmospheric pO2(Scott 2000, 2002; Glasspool et al.2015; Yan et al. 2019). However, since the Mesozoic, due to the diversification and complexity of coal-forming climatic zones, plant differentiation, coal-forming sedimentary environmental conditions, tectonic activities,sea level changes, and so on (Vail et al. 1977; Hunt and Smyth 1989), the local and dramatic changes in the abundance of inertinite in coal may have sometimes exceeded the influence of global atmospheric pO2.

6 Conclusions

In this investigation, the coal during each coal-forming period in China was taken as the research objects. The abundance of inertinite in coal during each coal-forming period was systematically analyzed. The atmospheric pO2characteristics during the coal-forming periods were successfully simulated and calculated. The following conclusions were reached in this study:

1) The systematic analysis of a large volume of data reveals the inertinite abundance and its evolution during the Phanerozoic coal-forming periods in China. Four major cycles were identified, in which the inertinite abundance first increased and then decreased. These are: Early Devonian to Late Permian; Late Triassic to Early Jurassic; Middle Jurassic to Late Cretaceous; and Paleogene to Neogene,respectively.

2) Based on a large amount of statistical data regarding the inertinite abundance in coal during each of the coal-forming periods in China, the atmospheric pO2characteristics of each of the coalforming periods were simulated. The evolution curves of atmospheric pO2during each coalforming period were quantitatively constructed,which provided a certain reference for further research regarding atmospheric pO2in geological history.

3) Through comparative studies on a global scale, the characteristics and evolution of the inertinite abundance and atmospheric pO2in various coalforming periods were further supplemented and improved, particularly in regard to the Jurassic.

4) It was further clarified in this study that the global atmospheric pO2was the main controlling factor of the inertinite abundance in global coal deposits.The relationship between the inertinite abundance in Paleozoic coal to the atmospheric pO2is established. However, since the Mesozoic, besides the diversity and complexity of coal-forming climatic zones,plant differentiation,coal-forming sedimentary environmental conditions, tectonic activities, and wildfire events, regional factors might have influenced global atmospheric pO2.

Abbreviations

pO2: Paleo-oxygen concentration; IAVG, %:Mean inertinite abundance;N: Number of samples; R: Percentage range (min-max), %; S: Standard deviation

Acknowledgements

We thank the editors and reviewers for their constructive comments to improve the manuscript.

Authors’ contributions

DDW proposed the main academic ideas of the manuscript, and participated in guiding the writing of the whole manuscript. LSY is mainly responsible for collecting data of inertinite abundance in coal in various geological periods,screening and processing data, and calculating paleooxygen content. LYS mainly makes a comprehensive analysis of the geological origin and evolution of the high inertinite abundance in coal. DWL mainly analyzes the development and distribution characteristics of inertinite abundance in coal in various geological periods. HYL mainly analyzes the reasons for the abnormally high inertinite abundance in Jurassic coal. SW mainly analyzes the geological reasons of inertinite abundance in Cenozoic coal. GQD assists in collecting geological data and drawing maps. All authors read and approved the final manuscript.

Funding

This work was supported by the National Natural Science Foundation of China (41972170, 41772096, 41402086), the Natural Science Foundation of Shandong Province (ZR2019MD021), and the SDUST Research Fund(2018TDJH101).

Availability of data and materials

All data generated or analyzed during this study are included in the manuscript.

Declarations

Competing interests

The authors declare that they have no competing interests.

Author details

1College of Earth Science and Engineering, Shandong University of Science and Technology, Qingdao 266590, Shandong Province, China.2College of Geoscience and Surveying Engineering, China University of Mining and Technology (Beijing), Beijing 100083, China.

Received: 19 July 2020 Accepted: 7 April 2021