Combined Magnetostratigraphy From Three Localities of the Rainstorm Member of the Johnnie Formation in California and Nevada, United States Calibrated by Cyclostratigraphy: A 13 R/Ma Reversal Frequency for the Ediacaran

A combined magnetostratigraphy for the Rainstorm Member of the Ediacaran Johnnie Formation was constructed using the sediment accumulation rates determined by rock magnetic cyclostratigraphy for three localities of the Rainstorm Member to provide a high resolution, time-calibrated record of geomagnetic field reversal frequency at a critical time period in Earth history. Two previously reported magnetostratigraphy records from Death Valley, California, the Nopah Range and Winters Pass Hills (Minguez et al., 2015), were combined with new paleomagnetic and cyclostratigraphic results from the Desert Range locality of the Rainstorm Member in south central Nevada, United States . The Johnnie oolite marker bed is at the base of each of the three sections and allows their regional correlation. The Nopah Range and Desert Range localities have similar sediment accumulation rates of ∼5 cm/ka, so their stratigraphic sections can be combined directly. The Winters Pass Hills locality has a higher sediment accumulation rate of 8.4 cm/ka, therefore its stratigraphic positions are multiplied by 0.6 to combine with the Desert Range and Nopah Range magnetostratigraphy. The thermal demagnetization results from the Desert Range locality isolates characteristic remanent magnetizations that indicate two nearly antipodal east-west and shallow directions and a mean paleopole (11.7˚N, 348.4˚E) that is consistent with “shallow” Ediacaran directions. The Desert Range also yields a magnetic susceptibility rock magnetic cyclostratigraphy that records short eccentricity, obliquity, and precession astronomically-forced climate cycles in the Ediacaran. The high-resolution combined magnetostratigraphy with nearly meter-scale stratigraphic spacing (nominally 23 ka, based on the Desert Range sediment accumulation rate), indicates 11 polarity intervals in a cyclostratigraphy-calibrated duration of 849 ka, indicating a reversal frequency of 13 R/Ma. The Rainstorm Member records the Shuram carbon isotope excursion, hence its age is ∼574 Ma. Given the recent cyclostratigraphy-calibrated reversal frequency of 20 R/Ma from the Zigan Formation (Levashova et al., 2021) at 547 Ma, our results show that reversal frequency was high but fluctuated during the Ediacaran.


INTRODUCTION
Evidence has been accumulating that the reversal rate of the Earth's magnetic field was extremely high in the latest Precambrian/earliest Paleozoic. Magnetostratigraphy studies from Ediacaran rocks (e.g., Pavlov and Gallet, 2010;Bazhenov et al., 2016;Meert et al., 2016) or middle Cambrian rocks (Pavlov and Gallet, 2001;Gallet et al., 2019) show large numbers of polarity reversals in sedimentary stratigraphic sections. Although all these studies showed a high reversal frequency based on estimated sediment accumulation rates, only the most recent study of the Ediacaran Zigan Formation from the Urals (Levashova et al., 2021) presented a calibrated sediment accumulation rate using rock magnetic cyclostratigraphy (Kodama and Hinnov, 2014). The Zigan Formation is dated to be 547 Ma in the latest Ediacaran, just before the Cambrian explosion of life at about 540 Ma. The reversal frequency obtained from Levashova et al. (2021) is 20 R/Ma (reversals per million years). For comparison, the highest reversal frequency, over the past 170 Ma, based on the sea floor magnetic anomaly record of geomagnetic polarity, is 10-12 R/Ma for the latest Jurassic/early Cretaceous. The polarity stratigraphy record from the astronomically-calibrated Newark Basin rocks (Kent et al., 2017), provides a well-calibrated record of reversals earlier than the Jurassic/Cretaceous during the Triassic and indicates 66 polarity intervals in the 33.2 million year period from 199.5 Ma to 232.7 Ma or an average reversal frequency of about 2 R/Ma. Hüsing et al. (2014) provide an astronomically-calibrated magnetostratigraphy in Lower Jurassic marine sediments from the United Kingdom and also indicate a low (∼2-3 R/Ma) reversal frequency at about 200 Ma. Li et al. (2017) see polarity reversal patterns in the astronomically-calibrated magnetostratigraphy of Upper Triassic Xujiahe Formation of South China similar to those in the coeval Newark supergroup rocks.

Rainstorm Member of the Johnnie Formation
Our study has recovered a high reversal frequency from the Rainstorm Member of the Ediacaran Johnnie Formation. The Rainstorm Member is the most widespread member of the Johnnie Formation outcropping throughout the southern Death Valley-Kingston Range area in California, the Spring Mountains, the Nevada Test Site and in the Desert Range of Nevada. It contains three units, the siltstone unit, the carbonate unit and the siltstone and quartzite unit. The basal siltstone unit is thin, typically less than 30 m thick. The overlying carbonate unit is comprised of pale-red to grey-red siltstone and limey silt. It ranges from 50-80 m in thickness (Stewart, 1970) and is the main unit sampled in this study. The base of our sampled section is marked by the Johnnie Oolite, a 2-3 m thick oolitic dolomite which occurs near the middle or in the upper half of the siltstone unit (Stewart, 1970). It serves as a marker bed to which the base of the stratigraphic sections at different localities are tied. This is important because it allows the stratigraphic positions of paleomagnetic data collected at various localities of the Rainstorm Member to be very accurately combined, critical for the magnetostratigraphy that we will present in this paper. The top portions of the sections are correlated by extending the sediment accumulation rates derived from the magnetic susceptibility cyclostratigraphy. In the section sampled in this study the Rainstorm Member strikes 326˚and dips 45˚to the northeast. Structural deformation in the area is related to displacements along Mesozoic and Tertiary thrust systems and strike-slip faults (Stewart, 1970).
We initially studied the Rainstorm Member of the Johnnie Formation because it carries a record of the Shuram carbon isotope excursion (LeGuerroue et al., 2006). The Shuram is the largest negative carbon isotope excursion in Earth history with a nadir of approximately −12 δ 13 C in inorganic carbon occurring globally just before the Cambrian explosion of multi-cellular life (Knoll and Carroll, 1999). Minguez et al. (2015) used magnetic susceptibility cyclostratigraphy and magnetostratigraphy to determine the duration of the Shuram excursion in the Johnnie Formation rocks. Minguez et al.'s (2015) rock magnetic, magnetic susceptibility, cyclostratigraphic study of the Rainstorm Member at the Nopah Range and Winters Pass Hills observed spectral peaks that were identified as short eccentricity, obliquity and precession Milankovitch cycles. We also studied the duration of the Shuram in the Wonoka Formation in Australia (Minguez and Kodama, 2017) and the Doushantuo Formation in south China (Gong et al., 2017). In all three studies we found a duration of 8-9 Ma for the Shuram carbon isotope excursion using rock magnetic, either magnetic susceptibility or anhysteretic remanent magnetization (ARM), cyclostratigraphy. We also found a normal polarity interval at the nadir of the carbon isotope excursion in Death Valley, California and in South Australia providing evidence that the excursion was synchronous globally (Minguez and Kodama, 2017). Investigating this relationship between field polarity and the nadir of the excursion was not possible in China because the Doushantuo Formation was found to be remagnetized in the Mesozoic (Gong et al., 2017;Gong et al., 2019).
We studied the Rainstorm Member of the Johnnie Formation at two localities in the Death Valley region, the Nopah Range and Winters Pass Hills (Minguez et al., 2015). (Van Alstine and Gillette 1979) reported a magnetostratigraphic study of the Rainstorm Member from the Desert Range of south-central Nevada about 130 km northeast of the Nopah Range and Winters Pass Hills localities. Our initial intent was to conduct a magnetostratigraphy and rock magnetic cyclostratigraphy study of the Rainstorm Member in the Desert Range to investigate the relationship of a normal polarity event with nadir of the Shuram carbon isotope excursion.
One critical aspect of our study is the age of the reversal frequency we observed in the Rainstorm Member. The Rainstorm Member records the onset of the Shuram carbon isotope excursion including its nadir. The end of the Shuram excursion has been dated at 550.55 ± 0.75 Ma (Condon et al., 2005) from an ash bed that was initially interpreted to be at the very top of the Doushantuo Formation. However, , work on the upper part of the Doushantuo Fm in the Three Gorges area suggests that the euxinic shale that records the end of the Shuram excursion (Member 4 of the Doushantuo Formation), is not directly below the ash layer dated by Condon et al. (2005), and hence the end of the Shuram may be much older than previously thought. Rooney et al. (2020) report Re-Os ages from different continents of the onset (574 ± 4.7 Ma) and end (567 ± 3.0 Ma) of the Shuram carbon isotope excursion. This would suggest the Rainstorm rocks we studied are approximately 574 Ma in age.
In this paper, we combine magnetostratigraphy results from the Rainstorm Member of the Johnnie Formation collected at three localities: the Nopah Range and Winters Pass Hills in Death Valley and the Desert Range in Nevada. The magnetostratigraphies are accurately combined using magnetic susceptibility cyclostratigraphy (Kodama and Hinnov, 2014) and the resulting reversal frequency is time calibrated using the rock magnetic cyclostratigraphy.

METHODS
We received permission from Nellis Air Force Base and from the Desert National Wildlife Refuge to sample the Rainstorm Member for 2 days in May 2018. Both entities control the sampling locality. The target area was that initially sampled by Van Alstine and Gillett (1979). Figure 1 shows the two locations in this area where we were able to locate exposures of the Rainstorm Member. Location 1 (Sites JO-7; Figure 1) gave us access to the bottom 7 m of the Rainstorm Member, including the Johnnie Oolite. Oriented hand samples were collected every meter of section. 25 mm diameter cores were drilled from these hand samples for paleomagnetic measurement and demagnetization at Lehigh University. At least three cores were drilled from each hand sample. Location 2 (Sites 16.5-45; Figure 1) allowed us to sample from 16.5 to 45 m above the Johnnie Oolite, which was also exposed at this location and thus allowed an accurate stratigraphic position for our sites. We sampled this location with a gasoline-powered modified chain saw drill (Pomeroy). We drilled three independently oriented samples per horizon and sampled the section at nominally 2 m stratigraphic intervals.
Unoriented rock samples approximately 3 × 3 × 3 cm in size were collected every 25 cm of section for the rock magnetic cyclostratigraphy study at both locations 1 and 2. The section was sampled nearly continuously for 45 m, with the exception of the stratigraphic gap between 7 and 16.5 m, where no exposures could be located in the 2 days of access we were granted.
The oriented paleomagnetic core samples were stepwise thermally demagnetized in a single chamber ASC Model TD-48SC thermal demagnetizer in 10-20 steps between 100-200˚C and 680˚C. Remanence was measured in a 2G Enterprises 755 superconducting rock magnetometer. Both thermal demagnetization and remanence measurement were conducted in a magnetostatically-shielded room built by Lodestar Magnetics (ambient field intensity ∼350 nT) in the Paleomagnetism Laboratory of Lehigh University. The bulk susceptibility of the unoriented cyclostratigraphy samples was measured using an Agico KLY-3s Kappabridge at Lehigh University. Three measurements were made on each sample and averaged. The susceptibility measurements were normalized by mass. The Kappabridge was also used to measure the temperature dependence of susceptibility from liquid nitrogen temperature (−196°C) to 0°and from 20°C to 700°C in order to identify magnetic mineralogy and the relative proportion of ferromagnetic and paramagnetic mineral contributions to the bulk susceptibility.
Principal component analysis was carried out on vector endpoint diagrams (Kirschvink, 1980) using the PmagPy software (Tauxe et al., 2016). Spectral analysis of the magnetic susceptibility time series was conducted using the 2π multi-taper method (MTM, Thomson, 1982). Before spectral analysis the time series was resampled at a 0.36 m even spacing using linear interpolation and a Lowess detrending model with a 35% window was applied. Spectral peak significance was estimated by robust red noise and its confidence limits at 90, 95, 99% calculated from Mann and Lees (1996). These calculations were carried out using the Acycle (v 2.0) software .
Carbon isotope (δ 13 C) values were measured for 24 samples from the Rainstorm member in order to determine if the rocks had recorded the Shuram carbon isotope excursion. Initially the objective was to determine the stratigraphic position of the nadir of the excursion for comparison to the magnetostratigraphy, but

Paleomagnetic Analysis
Principal component analysis of the thermal demagnetization vector endpoint diagrams (Zijderveld, 1967) isolated a characteristic remanent magnetization (ChRM) typically between the temperatures of 580°C and 680°C ( Figure 2). Maximum angular deviations (MADs) of the fit ranged from 1.3˚to, in one case, as high as 17.4˚, but averaged 6.1˚for 60 samples from 22 horizons. An equal area net plot of the ChRMs (Figure 3) clearly shows two nearly antipodal groups of directions indicating that in most cases a polarity assignment could be made. One direction is nearly due west and flat (D 268.7˚, I 11.7˚, N 31, α 95 17.0˚) which is interpreted to be normal polarity, based on North America's paleogeography in the late Precambrian (Cocks and Torsvik, 2011). The other direction is 155.6˚away pointing northeast and flat (D 65.0˚, I -4.8˚, N 29, α 95 9.6˚). Even given the uncertainties, the directions do not pass a reversals test, but polarity assignments can be made based on the great circle distance between the virtual geomagnetic poles (VGPs) calculated for these directions and the mean VGP (paleomagnetic pole) determined from these data. In addition to the two nearly antipodal groups of ChRMs, there are other groupings of directions ( In about 60% of the samples, the low temperature magnetization component is observed and is clearly separate from the high temperature ChRM. The low temperature component is typically isolated between 200°C and 500°-600°C (Figures 2A,B) and has a stratigraphic mean direction (D 9°, I 43°) similar to a low temperature component observed by Minguez et al. (2015) in the Rainstorm member of the Johnnie Formation in the Nopah Range, Death Valley, California. Minguez et al. (2015) attributed the low temperature Nopah Range component to Miocene age overprinting, based on its similarity to the North American Miocene paleopole (Torsvik et al., 2012). We suggest the same origin for the Desert Range low temperature component.
The paleomagnetic pole for the Rainstorm Member of the Johnnie Formation collected at the Desert Range is 11.7˚N,    (Bono et al., 2014) to be the primary magnetization of the Sept-Île intrusive rocks. A plot of the VGP paleolatitude (great circle distance from the paleopole-90°) ( Figure 4) as a function of stratigraphic position can be used to determine the magnetostratigraphy of the Rainstorm Member at the Desert Range locality. Multiple samples were collected at each horizon (nominally three samples) to check the robustness of the polarity interpretation for each horizon. VGP latitudes between 90˚and 45°are interpreted to be normal polarity, latitudes between −45°and −90°are interpreted to be reversed polarity. A 45°cutoff for polarity determination is  Stewart and Barnes (1966). (bottom) MTM spectrum of the magnetic susceptibility data series calculated by Acycle . Black line is robust red noise model (Mann and Lees, 1996) with 90% (thin red line), 95% (dashed red line), and 99% (dashed blue line) confidence intervals. Significant spectral peaks are outlined in boxes consistent with (Berger et al., 1992) theoretical peaks for short eccentricity (125/95 ka), obliquity (29 ka) and precession (18.7/16 ka) at 500 Ma. . The stratigraphic periods of the highly significant spectral peak assigned to short eccentricity are labeled for each plot. The Desert Range and the Nopah Range both have short eccentricity peaks close to a 5 m period, while the Winters Pass Hills has a short eccentricity peak with a longer stratigraphic period indicating a higher SAR. These relationships in SAR will be used to produce a combined stratigraphy from these three localities of the Rainstorm Member.
Frontiers in Earth Science | www.frontiersin.org November 2021 | Volume 9 | Article 764714 based on (Deenen et al., 's 2011) analysis of paleosecular variation. If at least one sample in a horizon falls within these bounds a polarity interpretation is made. If none of the samples in a horizon can be assigned a polarity under these conditions, the horizon's polarity is designated indeterminate. At the bottom of the section, the Johnnie Oolite's polarity is indeterminate. From 1 to 7 m the section is normal polarity. From 7 m to 16.5 m in the section is the 9.5 m stratigraphic gap in outcrop. From 18 m to 26 m the section is reversed polarity, followed by a normal polarity interval from 28 m to 30 m. The section is reversed polarity from 33 m to 42 m and returns to normal polarity at the top of the 45 m section. The magnetostratigraphy is fairly straightforward, but has a large 9.5 m gap and a nominal 2 m spacing between horizons with one completely indeterminate horizon directly above the 9.5 m stratigraphic gap.

Rock Magnetic Cyclostratigraphy
The mass normalized magnetic susceptibility (χ)data series measured for the unoriented samples has values typically of 4 × 10 −8 m 3 /kg but high values can reach 1.5 × 10 −7 m 3 /kg. (Figure 5). The data series shows both high frequency sub-meter/meter scale variability and lower frequency several meter scale variability. Multitaper method (MTM) analysis of the time series using Acycle  with robust red noise modeling (Mann and Lees, 1996) to identify significant spectral peaks is used to calculate the power spectrum. Spectral peaks that rise above the 95% confidence level of the robust red noise (dashed red line in Figure 5) are observed at 0.2 cycles/meter (5 m stratigraphic period), 0.65 cycles/m (1.5 m stratigraphic period) and multiple peaks between 0.9 and 1.15 cycles/m (1.1-0.9 m stratigraphic periods). These peaks can be matched quite easily to the Milankovitch peaks for short eccentricity (125-95 ka), obliquity (∼29 ka) and precession (18.7 FIGURE 7 | Magneti susceptibility vs temperature plots for Rainstorm Member samples. The top plots show the magnetic susceptibility behavior between room temperature (nominally 20°C) and 700°C. Red data is the heating curve, blue is the cooling curve. Note magnetic susceptibility shows evidence of a decrease during heating at about 680°C suggesting hematite as the main magnetic mineral carrying the susceptibility. The cooling curves both show a decrease at about 580°C, suggesting magnetite may have been destroyed during the heating. The bottom two plots show the magnetic susceptibility behavior between liquid nitrogen temperatures (−196°C) and 0°C. The poor fit of the 1/T paramagnetic hyperbola (green) to the data suggests that most of the susceptibility is carried by ferromagnetic, rather than paramagnetic, minerals. The data shown in the top plots would suggest that ferromagnetic mineral is hematite.
Frontiers in Earth Science | www.frontiersin.org November 2021 | Volume 9 | Article 764714 and 16 ka) calculated for 500 Ma by Berger et al. (1992). Using COCO (Correlation Coefficient) analysis from Acycle, a significant fit of these spectral peaks to the Berger et al. (1992) theoretical cycles yields a sediment accumulation rate (SAR) for the Rainstorm member of 5.3 cm/ka indicating that the 45 m of section was deposited in 849 ka. χ data series from Minguez et al. (2015) indicates strong, significant short eccentricity peaks at 4.8 m stratigraphic period for the Nopah Range and 8.4 m stratigraphic period for the Winters Pass Hills ( Figure 6). The Desert Range Rainstorm Member eccentricity peak by comparison is 5.2 m when the center of the peak is chosen. These results show that the SARs are nearly the same for the Nopah Range and Desert Range localities of the Rainstorm Member and approximately 1.67 times faster for the Winters Pass Hills locality. We will use these results to combine the magnetostratigraphic and isotopic data from all three localities.

Rock Magnetic Analysis
In order to identify the magnetic minerals in the Rainstorm Member and determine the magnetic minerals carrying the predominance of the magnetic susceptibility, magnetic susceptibility versus temperature experiments were conducted on a subset of Rainstorm Member cyclostratigraphy samples. In the high temperature experiments, samples heated from room temperature up to 700°C, show the signature of a magnetic mineral that loses its magnetism at about 680°C, which is most likely hematite (Figure 7). This observation is in agreement with the red color of the rocks in the Desert Range. However, the magnetic mineralogy of the Rainstorm Member can vary, as well as its color. In Winters Pass Hills the main ferromagnetic mineral is hematite, but at the Nopah Range the ChRMs are carried mainly by magnetite (Minguez et al., 2015). On cooling from 700°C a small decrease compared to the heating curve just below 580°C suggests that magnetite may have been destroyed by the heating. In the low temperature experiments, in which samples warmed from liquid nitrogen temperatures up to 0°C, when the magnetic susceptibility curve is fit with a 1/T hyperbola (Hrouda, 1994;Hrouda et al., 1997), it is clear that almost all the susceptibility is carried by the ferromagnetic minerals in the rock, i.e., hematite, rather than the paramagnetic clay minerals (Figure 7). Of the six samples measured the lowest ferromagnetic component was ∼58%, but most were above 90% ferromagnetic ( Table 1).

Carbon Isotope Analysis
The

DISCUSSION
The paleomagnetic and the rock magnetic results indicate that hematite carries the ChRM and the magnetic susceptibility (χ). The χ cyclostratigraphy results indicate that the hematite records global climate cycles. The paleomagnetic pole calculated from the ChRMs agrees with the "shallow" Ediacaran directions that (Bono et al., 2014) show is the primary Ediacaran direction. These observations and the record of a magnetostratigraphy in the Rainstorm Member supports a primary magnetization for the Rainstorm Member. Specifically, the ChRMs are isolated between 580°C and 680°C, the magnetic susceptibility versus temperature experiments show that the magnetic susceptibility persists until 680°C and that most of the susceptibility is carried by ferromagnetic minerals. The spectral analysis of the magnetic susceptibility data series shows significant spectral peaks that can be matched to short eccentricity, obliquity, and precession Milankovitch cycles in the Ediacaran. Correlation coefficient analysis (COCO) indicates a sediment accumulation rate (SAR) of 5.3 cm/ka for the Rainstorm Member in the Desert Range. These results strongly support an early origin, likely depositional, for the ferromagnetic mineral, hematite, carrying the magnetic susceptibility. The cyclostratigraphy at the Desert Range, the Nopah Range, and Winters Pass Hills localities can be used to construct a combined magnetostratigraphy from these three exposures of the Rainstorm Member. The short eccentricity (95-125 ka) peaks at the Nopah Range and the Desert Range both have a stratigraphic thickness period of ∼5 m, and hence these sections have similar SARs (∼5m/100 ka). Therefore, their independent magnetic polarity stratigraphies can be combined without any adjustments to the stratigraphic positions of their sampling horizons. The Winters Pass Hills locality has a short eccentricity peak with an 8.4 m stratigraphic thickness period indicating a faster SAR (∼8.4m/100 ka). Therefore, the stratigraphic positions of horizons in the Winters Pass Hills section should be multiplied by 5/8.4 0.6 to place its sampling horizons in an equivalent time scale before its magnetic polarity stratigraphy can be combined with the Desert Range and the Nopah Range magnetostratigraphies. The carbon isotope results from all three localities can be combined in this manner to check this approach for combining the three locality sections stratigraphically using the cyclostratigraphy. The combined carbon-isotope stratigraphy shows a good match of carbon-isotope values between the three sections ( Figure 8). Also, it is apparent that the 9.5 m stratigraphic gap in the Desert Range section is filled by data from the Nopah Range and Winters Pass Hills.
When the magnetic polarity stratigraphy from the three Rainstorm Member localities are combined in this manner, using the cyclostratigraphy, a high-resolution magnetostratigraphy results (Figure 9). The sampling interval is decreased to slightly more than 1 m. Thirty-eight horizons are spaced over 45 m of section, with an average spacing of 1.2 m. Using the SAR of 5.3 cm/ka indicates that this stratigraphic spacing is a temporal spacing of about 23 ka. The same rules used to determine the magnetostratigraphy for the Desert Range section (Figure 4) can be applied to the combined magnetostratigraphy. In the resulting combined magnetostratigraphy, there are some horizons in which all the samples are indeterminate. These horizons typically occur at polarity interval boundaries. They could have resulted from unresolved magnetization overprinting, or they may have deposited during polarity transitions. In some instances, samples with opposite polarities occur at nearly the same stratigraphic level, but this behavior occurs at polarity interval boundaries. It can be attributed to the stratigraphic position inaccuracies of combining the stratigraphic sections. FIGURE 9 | Combined magnetostratigraphy from the Desert Range, Nopah Range and Winters Pass Hills localities. Same color scheme for the data points from the localities as in Figure 8. VGP paleolatitude (GCD-90°) is plotted as a function of stratigraphic position combined in the same manner is in Figure 8 (carbon isotope plot). The rules for determining the polarity of each horizon are given in the text. Note the nearly meter stratigraphic spacing for the combined stratigraphy yielding a high resolution magnetostratigraphy that records 11 polarity intervals in the 45 m section. The magnetic susceptibility data series for the Desert Range is shown on the right (in red), with the data filtered (Gaussian filter, calculated by Acycle) at short eccentricity plotted on top (in blue). Assuming an SAR of 5.3 cm/ka from the Desert Range cyclostratigraphy, the duration of the 45 m section is 849 ka. This yields a reversal rate of 13 R/Ma. The resulting combined magnetostratigraphy is based on the Desert Range locality with an SAR of 5.3 cm/ka determined by COCO from the cyclostratigraphy. This rate can be used to estimate the duration of the 45 m section. The section was deposited over 849 ka. Therefore, 11 polarity intervals observed with the high-resolution, combined magnetostratigraphy in 849 ka indicates a reversal frequency of 13 R/Ma. This is slightly higher than the highest frequency observed in the past 170 Ma of 10-12 R/ Ma which occurred during the seafloor magnetic anomalycalibrated Jurassic/early Cretaceous time period.
Rooney et al.'s age estimate for the Shuram carbon isotope excursion between 574 Ma and 567 Ma suggests that the Rainstorm member, which records the early stages of the Shuram excursion is about 574 Ma in age. This age estimate makes the reversal frequency of 13 R/Ma about 27 Ma older than the cyclostratigraphy-calibrated reversal frequency of 20 R/Ma that Levashova et al. (2021) estimated for the Zigan Formation. It shows that reversal frequencies were high, but also highly variable during the Ediacaran. It is interesting to note that our reversal frequency estimate is close in age to the 565 Ma very low field intensity reported by Bono et al. (2019) for the Sept-Île intrusive rocks.

CONCLUSION
Rock magnetic cyclostratigraphy and paleomagnetic polarity stratigraphy at three localities in the southwestern United States of the Rainstorm Member of the Ediacaran Johnnie Formation produced a high-resolution, combined magnetostratigraphy for the age of the Shuram carbon isotope excursion at 574 Ma. This reversal stratigraphy is precisely time-calibrated by rock magnetic cyclostratigraphy and yields a reversal frequency of 13 R/Ma for the time period just before the very low field intensity observed in the Sept-Île intrusive rocks at 565 Ma (Bono et al., 2019). Given the recent cyclostratigraphy-calibrated reversal frequency of 20 R/Ma from the Zigan Formation (Levashova et al., 2021) at 547 Ma, our results show that reversal frequency was high but fluctuated during the Ediacaran.

DATA AVAILABILITY STATEMENT
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found below: doi.org/10.7288/V4/MAGIC/19237.

AUTHOR CONTRIBUTIONS
KK planned the study, organized and participated in the collection of the samples, measured the paleomagnetic samples, analyzed the data and wrote the paper.

FUNDING
This work was funded in part of United States National Science Foundation grant EAR-1322002 to the author.