Early Cretaceous Terrestrial Milankovitch Cycles in the Luanping Basin, North China and Time Constraints on Early Stage Jehol Biota Evolution

This research analyzes the cyclostratigraphy of the lacustrine Dabeigou Formation (DBG) of early Jehol Biota age ( ∼ 130–135 Ma) in the Luanping Basin, northern China. A high-resolution (2 cm interval), 117.82-m-long magnetic susceptibility (MS) stratigraphic series was measured along the Yushuxia section. MS is positively correlated with thorium, potassium and uranium concentrations associated with gamma ray intensity, and represents a proxy for detrital inﬂux to the Luanping Basin. Power spectral analysis identiﬁes a hierarchy of sedimentary cycles with wavelengths of 16.38 m, 5.85–3.28 m, 1.88–1.33 m, and 0.98–0.7 m, which are interpreted to represent Earth’s orbital eccentricity, obliquity and precession index cycles. Objective testing of the MS series supports the interpretation of Milankovitch cycles, indicating an average sedimentation rate of 4.642–4.723 cm/kyr. A ﬂoating astronomical time scale with a duration of 2478 kyr is established from interpreted 405 kyr long orbital eccentricity cycles along the MS series. The 405-kyr tuned DBG MS time series closely matches the predicted orbital eccentricity of the La2004 astronomical solution from 130.787 to 133.265 Ma, providing independent temporal constraints on early stage Jehol Biota evolution. Finally, this estimated time interval for the DBG MS time series indicates that it occurred entirely within the Weissert Event.


INTRODUCTION
The Early Cretaceous terrestrial Jehol Biota in East Asia is well-known for its exceptionally well-preserved feathered dinosaurs, early birds, mammals, pterosaurs, amphibians, insects and early flowering plants, and provides a unique window for understanding the evolution of Early Cretaceous terrestrial ecosystems (Zhou et al., 2003;Zhou, 2014). Geographically, the Jehol Biota is found in the Liaoning Province and adjacent areas in northeastern China. Stratigraphically, it is preserved mainly in the Dabeigou/Huajiying, Dadianzi, Yixian, and Jiufotang formations in ascending order in the Liaoning and Hebei provinces (Figure 1; Xu et al., 2019), with an age range of ∼135-120 Ma (Swisher et al., 1999(Swisher et al., , 2002He et al., 2006;Yang et al., 2007;Chang et al., 2009).
It has been increasingly recognized that the biota preserved in the Dabeigou Formation (DBG) reflect the origin and early stage Jehol Biota evolution (e.g., Tian et al., 2004;Ji et al., 2006;Wang and Ji, 2009;Zhou et al., 2009;Huang et al., 2015;Niu et al., 2015;Wang et al., 2015). The discovery of the earliest occurrence of the insect Ephemeropsis trisetalis and ostracods in the DBG, together with other fossils, suggests that the earliest Jehol Biota might already have appeared at the time when the lowermost part of Member 2 of the DBG was deposited (Zhou, 2014;Xu et al., 2019). The geochronology of the DBG indicates a range from ∼135 to ∼130 Ma (Liu et al., 2003;Zhang et al., 2005;He et al., 2006;Gao et al., 2018), but has large uncertainties, which impedes further progress in understanding early stage Jehol Biota evolution.
Cyclostratigraphy has been instrumental in establishing a 405-kyr-scale astronomical time scale (ATS) for the Cretaceous Period from well-preserved marine and continental successions (Hinnov, 2018;Huang, 2018). This has been possible due to the stability of the 405-kyr orbital eccentricity cycle over hundreds of millions of years (Laskar et al., 2004(Laskar et al., , 2011aKent et al., 2018). Cyclostratigraphy of the Cretaceous terrestrial strata from the Songliao Basin in northeastern China was interpreted with Milankovitch cycles, and established a ∼28 Myr-long ATS for the upper Turonian-Maastrichtian stages (Wu et al., 2009(Wu et al., , 2013a(Wu et al., , 2014. Subsequently, Wu et al. (2013b) analyzed cyclostratigraphy of the Lower Cretaceous Yixian Formation in the Sihetun Basin, northeastern China, estimating an average sedimentation rate of ∼1.70 cm/kyr for the Jianshangou Beds that host "feathered" dinosaur/primitive bird fossils of the Jehol Biota. Most recently, Liu et al. (2020) found that the cyclic alluvial-fluvial synrift deposits of the Lower Cretaceous (Valanginian-Hauterivian) Shahezi Formation in Songliao Basin were controlled by astronomical forcing, and reconstructed an ∼11.14 Myr-long floating astrochronology based on interpreted 405-kyr cycles.
The Lower Cretaceous DBG in the Yushuxia section, Luanping Basin (LPB), Hebei Province, the focus of this study, consists of 213 m of dominantly cyclic continental strata in wellexposed outcrops with abundant fossils (Liu et al., 2002;Tian et al., 2004;Wu et al., 2004;Zhang et al., 2007;Zhou et al., 2009;Niu et al., 2010). It provides a unique opportunity to develop an astrochronology to enable characterization of terrestrial climate change and environmental evolution of the Early Cretaceous Period. In this study, we conducted detailed time series analysis and modeling of the high-resolution (2 cm interval) magnetic susceptibility (MS) series of the DBG. The objectives were to search for and identify astronomical signals in the DBG, to develop an astrochronology for the early stage Jehol Biota, and to test the reliability of the astronomical solutions for early Cretaceous time.

GEOLOGICAL SETTING
The Luanping Basin (LPB) is located in northeastern Hebei Province, northern China, near the northern margin of the Yanshan orogenic belt (Figure 1). The paleo-latitude of LPB was ∼40 • N during the Cretaceous Period, i.e., approximately the same as today (Wang, 2013;Wang C. S. et al., 2013). Here, the Volcanic Sedimentary Basin Group was formed by the active tectonics and frequent volcanic eruptions during the Late Jurassic to Early Cretaceous . The LPB is one of the small, terrestrial extensional basins in the Yanshan structural belt with a complete sedimentary sequence (Wu et al., 2004;Zhang et al., 2007). Its sedimentary evolution involved three successive stages: (1) "vigorous volcanic eruptions, " (2) "extensional subsidence, " and (3) "sedimentary infilling" (Zhang et al., 2007). The Jurassic and Cretaceous terrestrial strata in LPB are divided into six lithologic formations (in descending order): Jiufotang, Xiguayuan, Dadianzi, DBG, Zhangjiakou and Houcheng formations (Tian et al., 2008;Wang and Ji, 2009;Wang S. E. et al., 2013).
The DBG of the 213 meter-thick Yushuxia section (Figure 1) is well-known for its continuous non-marine Lower Cretaceous deposits and preservation of early stage Jehol Biota, including plants, pollen, spores, spinicaudatans, gastropods, bivalves, ostracods, insects, and fishes (Liu et al., 2002;Tian et al., 2004;Li et al., 2004;Zhou et al., 2009;Wang et al., 2015). The stratigraphy is characterized by alternating deposition of mudstones, shales, and siltstones, interbedded with sandstones and conglomerates of semi-deep lake facies in the lower part and deep lake facies in front of fan deltas in the upper part (Figure 2; Liu et al., 2001;Wang and Ji, 2009).
The DBG is divided into 3 members and 46 layers according to a detailed description of sedimentology and stratigraphy Qin et al., 2018). Member 1 (49.39 m) consists of tuffaceous, pebbly and coarse-grained sandstones, and medium to fine-grained sandstones intercalated with tuffaceous siltstone, sedimentary tuff and thin volcanic ashes. Member 2 (98.44 m) consists of siliceous mudstones and shales, interlayered with siltstones, fine to medium-grained sandstones; abundant fossils occur in this member. Member 3 (76.68 m) consists of mudstones, silty mudstones, calcareous mudstones and pelitic siltstones, and abundant fossils. The DBG represents shorelineshallow lake to fan delta front-semi-deep lake facies in the lower part, and deep lake to fan delta front facies in the upper part (Liu et al., 2001).
The age of the DBG has been disputed for a long time. The paleontological data, including ostracods, spinicaudatans, palynotaxa, and plant macrofossils indicate that the DBG ranges from Late Jurassic to Early Cretaceous (e.g., Wan et al., 2016;Qin et al., 2018;Xi et al., 2019). Radioisotopic dating provides further constraints: in the Yushuxia section, SHRIMP zircon U-Pb dating of the tuffs in Layers 34 and 21, and the volcanic rocks of the underlying Zhangjiakou Formation indicates ages of 133.9 ± 2.5 Ma, 129.9 ± 1.2 Ma, and 135.4 ± 1.6 Ma, respectively (Figure 2; Liu et al., 2003;Gao et al., 2018). Zhang et al. (2005) obtained a LA-ICP-MS zircon U-Pb age of 135.2 ± 2.3 Ma from the Zhangjiakou Formation in the same section (Figure 2). 39 Ar/ 40 Ar dating of a tuff layer in the upper DBG of the Jiecaigou section, near LPB, yielded an age of 130.7 ± 1.2 Ma (He et al., 2006). Thus, the DBG can be constrained roughly from ∼135 to ∼130 Ma. The most recent constraints for Valanginian-Hauterivian time

DATA
Magnetic susceptibility (MS) is a measure of the degree of magnetization of a material when subject to an external magnetic field. In sediments and sedimentary rocks, MS is affected by the concentration, grain size and shape of magnetic minerals (Kodama, 2012). Many studies have shown that MS is an effective indicator for Milankovitch-scale signals in stratigraphy, and it has been extensively used to characterize both marine and terrestrial sediments (e.g., Wu et al., 2013b;Zhong et al., 2018;Kodama, 2019).
MS was measured at 0.02 m intervals (on average) along the recently exposed Yushuxia section using a portable SM-30 MS detector. A total of 117.82 meters was measured, for a total of  Frontiers in Earth Science | www.frontiersin.org TABLE 1 | Orbital eccentricity, obliquity, and precession index periodicities from the periodogram of the La2004 astronomical solution (Laskar et al., 2004), between 129.55 and 134.55 Ma, i.e., a 5 Myr-long interval with a median age 132.05 Ma, calculated by Acycle for the TimeOpt astronomical target (orbital eccentricity and precession index) and for the ASM and COCO astronomical frequency targets (all three astronomical parameters).

Astronomical parameter
Periodicity ( Figure 2B). Gamma ray intensity was measured at 0.1 m intervals along the section using a portable RS-230 GR detector; in this study measurements are presented for Layer 19 only for comparison (Figure 3).

Pre-processing
The Prior to analysis, the data of 16 ash beds was removed, and MS series was interpolated to the average sample rate of 0.02 m. A low-pass filter was applied to the MS stratigraphic series using tanerfilter.m to retain astronomical frequencies while removing the very high frequencies. A cut-off frequency of 5 cycles/m and roll-off rate of 10 12 were defined to reject wavelengths shorter than 0.2 m. The low-pass filtered MS series was then smoothed with a 40-meter-long window with the Matlab function smooth.m with the "lowess" option to estimate irregular long-term trends that could interfere with the detection of low frequency orbital eccentricity cycles. This smoothed curve was then subtracted from the MS series ( Figure 2C).

Spectral Analysis
Multi-taper method (MTM) power spectral analysis (Thomson, 1982) with pmtm.m, and the evolutionary Fast Fourier Transform (FFT) spectrogram with evofft.m (Kodama and Hinnov, 2015) were used to characterize the frequency content of the processed MS series. The MTM F-ratio test was used to determine significant harmonic lines. The MTM harmonic F-test used in the average spectral misfit analysis (see below) was performed using the "eha" function of the Astrochron package in R (Meyers, 2014); otherwise F-test significance values were computed with Acycle . The application of frequency ratios, e.g., 20:5:2:1 for long orbital eccentricity (405 kyr), short orbital eccentricity (100 kyr), obliquity (36.6 kyr), and precession (20 kyr) served as a preliminary test for astronomical frequencies in the MS stratigraphic series (Huang et al., 1992;Mayer and Appel, 1999;Laskar et al., 2004).

Sedimentation Rate Modeling
Three objective methods were used to test jointly for the presence of astronomical frequencies and the most probable sedimentation rate for the DBG MS stratigraphic series: 1. Average spectral misfit (ASM) analysis (Meyers and Sageman, 2007) estimates an optimal sedimentation rate given a set of astronomical target frequencies that is compared with the set of statistically significant frequencies of the stratigraphic data series for a specified range of sedimentation rates. The sedimentation rate with the lowest misfit between data and target frequencies is the "optimal sedimentation rate" (Supplementary Material).

Time Optimization (TimeOpt) analysis evaluates orbital
eccentricity-like variations in the data, together with supporting evidence from amplitude modulations of the precession index band, across a test range of sedimentation rates (Meyers, 2015(Meyers, , 2019   estimates the correlation coefficient between power spectra of a stratigraphic proxy series and an astronomical solution in the time domain for a range of sedimentation rates. The evolutionary correlation coefficient (eCOCO) procedure was used to investigate changes in sedimentation rate along the stratigraphic series . COCO and eCOCO analysis with 5000 Monte Carlo simulations was performed on test sedimentation rates ranging from 1 to 10 cm/kyr (Supplementary Material).
The ASM procedure was carried out using the Astrochron package in R (Meyers, 2014); TimeOpt, TimeOptSim, COCO, and eCOCO procedures were performed using the software ACycle v 2.1 . The astronomical target frequencies used in these three procedures are based the La2004 astronomical solution periodogram for the 5-myr-long interval 129.55-34.55 Ma (median age of 132.05 Ma) ( Table 1; Laskar et al., 2004).

The La2004 Astronomical Solution
The accuracy of astronomical solutions declines rapidly prior to 50 Ma, and solutions other than the well-known La2004, e.g., La2010 and La2011, are thought to incorporate parameterizations and initial conditions that are superior to those used in La2004 (Laskar et al., 2011a,b), and more closely fit geological data around and prior to 50 Ma, e.g., La2010c and ZB18a (Zeebe and Lourens, 2019). The main issue centers on modeling the fundamental frequencies that involve the orbits of Earth (g 3 , s 3 ) and Mars (g 4 , s 4 ), and the expectation for chaotic interactions between these two planetary orbits through deep time. For example, La2004 from 92 to 90 Ma and 87 to 85 Ma, predicts chaotic transitions between the Earth and Mars orbits. Recently, it was discovered that these two transitions closely fit Late Cretaceous cyclostratigraphic patterns from the Western Interior Seaway (Ma et al., 2017(Ma et al., , 2019. Here we will test whether the La2004 solution also fits cyclostratigraphic patterns of the Early Cretaceous time represented by the DBG.

MS Variations
MS values range from −2.4 × 10 −8 SI to 4.2 × 10 −7 SI, with an average of 7.02 × 10 −8 SI ( Figure 2B). The MS stratigraphic series shows a sustained cyclic pattern that is consistent with the decimeter-to meter-scale lithological cycles (Figures 2, 3). MS is positively correlated with gamma ray intensity and elevated 232 Th, 40 K, and 238 U (Figure 3). Th contributes ∼60% of the GR signal, which identifies GR as a primarily clay proxy ("shale log") (Crain, 2019) that tracks the detrital influx from surrounding terrestrial source areas (Ohta et al., 2011). Recently, Pei et al. (2019) identified primary magnetic remanence and constructed a magnetostratigraphic framework for the upper DBG in LPB, further supporting the validity of the MS series as an important paleoenvironmental and paleoclimatic proxy. We propose that DBG MS is a proxy for the terrestrial hydrologic cycle, weathering and fluvial intensity, and if these processes were forced by insolation, represent Milankovitch cycles.
The DBG Stratigraphic Spectrum 2π MTM power spectral analysis of the MS stratigraphic series indicates spectral peaks at wavelengths of 16.38 m, 5.85-4.82 m, 3.28-1.42 m, and 0.98-0.7 m above the 90% F-test value (Figures 4A,B). The ratio of the major cycle bands of 20.5-16.4 m: 5.85-4.82 m: 3.28-1.42 m: 0.98-0.7 m is similar to the ratio of the Early Cretaceous astronomical periodicities of long orbital eccentricity (405 kyr) : short orbital eccentricity (132, 124, 99, and FIGURE 5 | ASM analysis of the DBG MS series, applied over a range of sedimentation rates from 1 to 10 cm/kyr by 0.0225 cm/kyr steps, using the La2004 astronomical target frequencies for 129.55-134.55 Ma ( Table 1). All harmonics in the MS series exceeding 85% significance level, a total of 42, are evaluated. The astronomical target used is based on the list in Table 1, but has reduced the number of short eccentricity terms from four (1/132.129, 1/124.1212, 1/99.9024, 1/94.1609) to two (1/132.129, 1/95.25), which was found to yield a significantly lower H 0 (see Supplementary Material for ASM analysis with all 10 terms). (A) The ASM statistic between data and target spectrum. (B) The results of null hypothesis testing: H 0 is the probability of no target astronomical frequencies in the data. The optimal sedimentation rate occurs at the first joint lowest value ASM statistic and statistically significant H 0 minimum in the tested sedimentation rate range. (C) The number of astronomical terms used in the statistic; a total of 8 terms (frequencies) are included in the target (see above). (D) Comparison of significant spectral lines (dotted black lines) evaluated at the optimal sedimentation rate identified in (C) and the astronomical target lines (solid red lines). Vertical dash-dot red lines indicate solutions for optimal sedimentation rates. 95 kyr) : obliquity (37 kyr) : precession (23.0, 22.5, 21.8, 18.7, and 18.5 kyr) (Laskar et al., 2004). These cycles may be correlated with Milankovitch cycles of short orbital eccentricity, obliquity and precession according to the cycle length ratios. The evolutionary FFT spectrogram of the MS stratigraphic series ( Figure 4A) indicates a sudden shift of power to higher frequencies in the 35-65 m interval, followed at ∼65 m by an equally sudden shift back to the lower frequencies that is then maintained to the base of the series. The question is whether these shifts indicate sudden changes in sedimentation rate, or evolution of the astronomical frequencies.

ASM Analysis
ASM analysis indicates a statistically significant (H0∼1%) solution at 4.723 cm/kyr (Figure 5). However, only five of the total of eight astronomical target frequencies were evaluated for this result, with the long orbital eccentricity at 409.6 kyr, and short precession index terms at 18.08 and 17.9256 kyr not successfully fitted to within the half-Rayleigh spacing required by the ASM statistic (Equation 1 in Meyers and Sageman, 2007), as indicated by three vertical red lines that do not coincide with any of the black vertical dashed lines in Figure 5D.

TimeOpt/TimeOptSim Analysis
TimeOpt analysis (Figure 6) indicates a r 2 opt maximum at a sedimentation rate of 4.642 cm/kyr, for which TimeOptSim analysis estimates a p-value of 0.119 ( Figure 6I). This is close to the 4.723 cm/kyr sedimentation rate solution obtained by ASM (see above). The r 2 envelope assessment ( Figure 6A) indicates a much faster 9.224 cm/kyr with a very low p-value of 0.05348 ( Figure 6G); r 2 power ( Figure 6A) indicates 6.339 cm/kyr, and a much higher p-value of 0.47128 ( Figure 6H). Adopting 4.642 cm/kyr as the optimal sedimentation rate, the MS series has a periodogram with a strong orbital eccentricity signature ( Figure 6F).

COCO/eCOCO Analysis
COCO analysis indicates three maxima at 3.32, 4.72, and 6.34 cm/kyr, the largest at 4.72 cm/kyr with a correlation coefficient value exceeding 0.4 ( Figure 7A). All three sedimentation rates have H 0 significance levels lower than 0.01 ( Figure 7B). The results of eCOCO analysis shows that the sedimentation rates vary from 4.72 to 6.3 cm/kyr ( Figure 7A).

Sedimentation Rate Variations in the DBG MS Series
The optimal sedimentation rates indicated by the ASM, TimeOpt and COCO methods are summarized in Table 3. Only one sedimentation rate is shared by all three methods: 4.642 cm/kyr (TimeOpt) to 4.72 cm/kyr (ASM, COCO). Other sedimentation rates arising at 3.4 and 6.34 cm/kyr can be explained as the result of sedimentation rate variations, as revealed by eCOCO (Figure 7). The TimeOpt precession index amplitude envelope model result indicating 9.224 cm/kyr (red circles, Figure 6A) with an extremely low P-value of 0.05348 (Figure 6G), is twice the joint optimal sedimentation rate of 4.642 cm/kyr that has a higher P-value of 0.11902, but shows an excellent fit of orbital eccentricity model frequencies in the periodogram (from 0 to 0.12 cycles/kyr) ( Figure 6F). It is worth further considering this periodogram, and rescaling the frequencies by a factor of 2 (for 9.224 cm/kyr). This would shift the data spectral peaks that are currently at ∼100 kyr (0.01 cycles/kyr) to ∼50 kyr, and ∼400 kyr (0.0025 cycles/kyr) to ∼200 kyr, which do not coincide with the orbital eccentricity terms. Indeed, r 2 power , which compares the orbital eccentricity model and sedimentation rate-calibrated data, records a low value at 9.224 cm/kyr ( Figure 6A). The joint modeling by r 2 opt steers the TimeOpt solution toward  Figure 2B) that were identified with elevated statistical significance by the MTM harmonic F-ratio test ( Figure 4A) and bandpass filters ( Figure 2D).

Astrochronology of the DBG MS Series
The astronomical frequencies identified in the DBG are supported by coeval evidence for astronomical forcing in the Valanginian-Hauterivian marine systems (e.g., Fiet et al., 2006;Martinez et al., 2015;Aguirre-Urreta et al., 2019). The optimization procedures reveal that the 20. 5-16.4 m, 5.85-4.82 m, 3.28-1.42 m, and 0.98-0.75 m wavelengths in the MS stratigraphic series are associated with the periodicities of the long and short orbital eccentricity, obliquity and precession index, respectively. Notably, assuming that the 16-20 m thick MS cycles along the section are 405 kyr cycles indicates sedimentation rates that are closely aligned with the objective testing (Figure 7: vertical dashed and solid black lines in the eCOCO color maps). The recognition of 405-kyr cycles is important for the ongoing initiative to define cyclostratigraphy in terms of the 405-kyr g 2 -g 5 orbital eccentricity metronome for the geologic time scale (Kent et al., 2018;Hinnov, 2018). Applying the 405kyr-based chronology model ( Table 2) to convert the DBG MS stratigraphic series to a time series results in a spectrum with power concentrated in the short orbital eccentricity (influenced in part by three 100-kyr scale time points in Table 2), obliquity [1/(40.9 kyr), 1/(36.9 kyr) and 1/(27.6 kyr)] and precession index [1/(22.1 kyr) and 1/(18.9 kyr)] bands ( Figure 4C). The spectrogram of the 405-kyr tuned MS series also indicates that the long and short orbital eccentricity trade positions in dominance at ∼400-500 kyr intervals along the section (Figure 4D).
The DBG MS cyclostratigraphy and the La2004 astronomical solution from 135 to 130 Ma are shown together in Figure 8. The yellow shading indicates a proposed correlation interval between data and model, from 133.265 to 130.787 Ma. This specific correlation was selected according to the singular coincidence at 131.538 Ma of an obliquity maximum and an orbital eccentricity minimum in both the La2004 astronomical solution and the 405kyr tuned DBG time series (green vertical line in Figure 8). The obliquity-filtered DBG series (Figure 8D) also shows a pronounced amplitude modulation cycle that repeats every 6-7 cycles, which could be evidence for the recently described 173kyr s 3 -s 6 metronome, that is also present in the La2004 obliquity ( Figure 8B; Boulila et al., 2018;Hinnov, 2018).

Paleoclimate Implications: The Weissert Event
The Weissert Event, characterized by a positive marine d 13 C excursion and global cooling, is the first major Early Cretaceous Earth system perturbation in the geologic record (Erba et al., 2004;Gröcke et al., 2005;Gréselle et al., 2011;Föllmi, 2012;Bajnai et al., 2017;Price et al., 2018). It is thought to span as much as 5.85  from the start of the Late Valanginian to the end of the Early Hauterivian (Martinez et al., 2015). By this measure, the DBG, from 133.265 to 130.787 Ma (Figure 8), falls entirely within the Weissert Event. Geochemical analysis of the DBG indicates that the formation was deposited under relatively low-intensity hinterland weathering related to the boreal or semiarid regions; the upper unit of the DBG records hinterland weathering that is comparable to that in temperate and humid mid-latitude regions (Ohta et al., 2011). Palynoflora data also indicate that the paleoclimate changed from semi-humid to relatively temperate and humid conditions during DBG time (Wang et al., 2016).
The change toward more humid conditions during the late Berriasian and intensification during the Valanginian Weissert Event evidently influenced terrestrial biota, and especially the evolution of herbivore vertebrates. A study of oxygen isotope composition of vertebrate apatite in Liaoning Province, China, indicates global cooling during Jehol Biota time, demonstrated by an estimated air temperature of 10 ± 4 • C (Amiot et al., 2011). Zhou (2014) also found temperate-to-cool-temperate fossil wood genus Xenoxylon in the Jehol Lagerstätte and an absence of crocodilians in northeastern China.

CONCLUSION
The DBG in the LPB, Hebei Province, northeastern China represents the prelude to the Jehol Biota, the spectacular Lagerstätte in the Lower Cretaceous. However, the DBG also occurs in a geologic interval that is chronostratigraphically very poorly constrained. This state of affairs motivated measuring the DBG for MS to search for astronomical signals during DBG time, with the larger goal of developing an astrochronology for the formation. The results are as followings: • A high-resolution (2 cm measuring interval) MS stratigraphic series was collected from the predominantly cyclic continental (

DATA AVAILABILITY STATEMENT
The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation, to any qualified researcher.