A U-Pb zircon age constraint on the oldest-recorded air-breathing land animal

The oldest-known air-breathing land animal is the millipede Pneumodesmus newmani, found in the Cowie Harbour Fish Bed at Stonehaven, Scotland. Here we report the youngest, most concordant 238U-206Pb zircon age from ash below the fish bed of 413.7±4.4 Ma (±2σ), whereas the youngest age from a tuffaceous sandstone above the fish bed is statistically indistinguishable at 414.3±7.1 Ma. The Cowie Harbour Fish Bed thus appears to be lowermost Devonian (Lochkovian), contrary to the previously accepted mid-Silurian age based on palynomorphs from adjacent exposures. This has implications for the evolutionary timetable of land colonization, as the Cowie ages overlap late Lochkovian zircon ages reported elsewhere for andesite below the nearby (~50 mi) Rhynie Chert, which has more advanced terrestrial biota. The results postdate the possible late Silurian Ludford Lane locality in Shropshire, England. Pneumodesmus newmani is thus not the earliest air-breathing land animal, unless the Ludford Lane locality is younger than presently assigned.


Introduction
Evaluating the rate and character of land colonization requires precise dating of early land biotas-but this is often difficult because of the problem of correlating and dating non-marine deposits. The supposed oldest air-breathing land animal found to date, the millipede Pneumodesmus newmani, is from the Cowie Harbour Fish Bed at Cowie, Aberdeenshire, Scotland [1,2] (Fig 1). It is currently assigned a late Wenlock to early Ludlow age (mid-Silurian,~427 Ma) based on plant spores from supposedly correlative, but tectonically isolated, and undescribed small exposures of sediment at Carron wood and Carron water inland (Fig 1) [1,3,4]. Fish and arthropod fossils from the Cowie Harbour Fish Bed itself suggest a Pridoli (latest Silurian) age of 423-419 Ma [5]. No spores are recorded from the Cowie Harbour Fish Bed, which is unexpected if correlative with the inland exposures due to its close proximity.
No radiometric dates exist from the succession, thus we applied U-Pb Laser Ablation-Inductively Coupled Plasma Mass Spectrometry (LA-ICP-MS) to zircon grains extracted from adjacent beds above and below the Cowie Harbour Fish Bed (Fig 2). Zircon ages from these units have important implications for the development and rate of evolution of early land a1111111111 a1111111111 a1111111111 a1111111111 a1111111111 the area is permitted. This is also stated specifically in [8]. The field studies did not involve endangered or protected species.
Sample 16COW1 is from a fine-grained green tuff at the top of the Cowie Harbour Conglomerate member, whereas sample 16COW2 is from a coarse bed within the Red Sandstone member directly below the fish bed. Sample 16COW3 is from a fine-grained tuffaceous sandstone located within the Cowie Harbour Siltstone member 20m above the fish bed. Using standard heavy-mineral separation methods, sample 16COW2 surprisingly contained no zircons, but the others supplied >100 grains that are typically 50-100μm in length (Fig 3). One explanation for the absence of zircon in 16COW2 is that small, fine sand-sized zircons remained in suspension in the higher energy environment indicated by the coarse grain-size of the red sandstone [9]. The zircon absence in 16COW2 and the abundant small grain sizes seen in 16COW1 and 16COW3, also suggest that only fine-grained zircons may have been available for transport, and thus Cowie may have been a long distance from the eruption source.
Forty zircon grains each from sample 16COW1 and 16COW3 were selected and dated using an Element2 High Resolution (HR)-ICP-MS with an Excimer (192 nm) laser ablation system instrumentation in the Geo-Thermochronometry lab at the University of Texas at Austin. Each zircon was mounted on double-sided tape, and was examined optically during the mounting process to select euhedral grains and eliminate analysis of cracked or metamict zircons (Fig 3). In this dating method, an ablated dry aerosol is introduced into the High Resolution ICP-MS using ultra-high purity He carrier gas for 238 U, 235 U, 232 Th and 206 Pb-208 Pb isotopic measurements using ion-counting. Each analysis consisted of a 2-pulse cleaning ablation, a background measurement taken with the laser off, a 30 sec measurement with the laser firing, and a 30 sec cleaning cycle. The laser beam was 30μm diameter. Common Pb was corrected using the measured 204 Pb (Hg-corrected) and assuming initial composition reported by [10]. Elemental and isotopic fractionation of Pb/U and Pb isotopes, respectively, is corrected by interspersed analysis of primary and secondary zircon standards with a known age (GJ1 and Plesovice) [11,12]. The common unknown to standard measurement ratio was generally 3:1 or 4:1. Uncertainty resulting from calibration correction is 1-2% for both 206 Pb/ 207 Pb and 206 Pb/ 238 U. We attempted to date the youngest grains in each sample a second time, but the results were discordant. The youngest grains in Fig 3 have two laser ablation spots.
Data is summarized in Tables 1 and 2, with details provided in S1 Table. We report the 238 U-206 Pb age if the zircon was less than 850 Ma, and the 207 Pb-206 Pb age if the grain was older than 850 Ma. Ages were filtered for concordance, and 238 U-206 Pb ages that are >10% discordant are not reported (N = 9 out of 83 total). In addition to the age data, U (ppm) and U/Th for each analysis were obtained. Note that the whole zircon grains, as opposed to exposed zircon cross sections, were analyzed due to their small size. The time-resolved laser ablation technique allows us to evaluate data that is obtained as the laser drills continuously into the grain, thus identifying zones of different age and limiting the possibility of mixed ages. This process allows us to distinguish between zircon cores and rims, as different age populations as easily identified. All ages are reported with ±2σ uncertainty.

Results
Overall, 38 zircon ages are reported from sample 16COW1 and 36 ages from 16COW3. Ages from 16COW1 range from 2684±32 Ma (5.66% discordant) to 410.4±7.2 Ma (3.37% discordant), whereas those from 16COW3 range from 1688±27 (4.32% discordant) to 414.3±7.1 Ma (0.65% discordant). We found ten Devonian zircon ages (N = 7 from 16COW1 and N = 3 from 16COW3) and eleven that yield Silurian results (N = 10 from 16COW1 and N = 1 from 16COW3). The average Devonian ages from 16COW1 is 414.2±8.7 Ma (MSWD = 0.82), whereas the Devonian results from sample 16COW3 average 416.0±8.4 Ma (MSWD = 0.64). One Devonian age in sample 16COW3 (415.0±12.0 Ma) is found as a rim on a 651±18 Ma core. This is the only zircon that contains an older core and younger rim; the remainder of the grains did not show age gradients using our filters for discordance. The majority of the zircon grains analyzed are older than Silurian (55% of all results from 16COW1 and 89% of all results   in 16COW3). We find no correlation between zircon age and U (ppm) or U/Th. The U/Th data is consistent with igneous zircons (e.g., [13]).

Discussion
The weighted means of zircon ages from volcanic ashes, which are frequently used to define bed ages, are averages of the ages of zircons blown out by the explosive eruption which produced the ash, and must include some older material from the volcanic pile and country rock. Ultimately, they do not give the age of the eruption, but only an older age, depending on the amount of earlier material incorporated [14,15]. They are meaningless in terms of the age of eruption of the ash that contains them [15]. Thus, we use the youngest, most concordant zircon ages as the maximum age of each bed, since a bed cannot be older than the youngest thing in it, although it can be younger (Fig 4; Tables 1 and 2). The youngest, most concordant U-Pb zircon age from the lowermost bed is 413.7±4.4 Ma (grain 6, Table 2), whereas the uppermost bed is 414.3±7.1 Ma (grain 4, Table 2). These are statistically indistinguishable, indicating rapid sedimentation, consistent with the paleo-environment. The Cowie Harbour Fish Bed thus appears to be~414 Ma in age, almost 13 million years younger than the currently accepted age. The younger zircon ages have the potential to be due to Pb loss from original mid-Silurian grains. However, zircons are often used to time major events in Earth's history because they retain their daughter isotopes under a range of extreme conditions. If the~414 Ma zircons were originally~427 Ma, the zircons would have had to experience~4% radiogenic Pb loss. Although this is a remote possibility, we selected clear, prismatic grains that lacked metamict features or textures suggestive of Pb loss. In addition, Fig 4C and 4D shows that a significant number of zircons yield the younger results, which requires they be considered when evaluating the age of Pneumodesmus newmani.
The uncertainty in the ages obtained using LA-ICP-MS is larger than ages that would be obtained using other approaches, such as isotope dilution thermal ionization mass spectrometry (ID-TIMS) analyses. The ID-TIMS approach applied to the youngest zircons would yield more precise ages since the two-sigma precision of individual ID-TIMS analyses is 0.1-0.3% compared with about 1-3% of individual LA-ICP-MS [16]. However, the only zircons available in the age range are small and aliquots of multiple grains would likely include the higher amounts of older material. We attempted a second LA-ICP-MS analysis of the youngest grains in the samples, but these yielded only discordant results. We anticipate that the TIMS analysis is almost certain to incorporate material from more than one zone and it is impossible to mechanically separate these in such small grains.
The Silurian-Devonian boundary is dated at 419.2±3.2 Ma [17]. Therefore, the Cowie Harbour Fish Bed and its air-breathing millipede are lowermost Devonian (Lochkovian) in age, contrary to the previously accepted Silurian age-on the current time scale. This has important implications for the evolutionary timetable of land colonization, as the Cowie age range statistically overlaps a U-Pb zircon age of 411.5±1.3 Ma (late Lochkovian, average of N = 4 grains, TIMS) for an andesite just below the Rhynie Chert with its more advanced and better preserved terrestrial biota [18] and is possibly younger than the supposed land centipedes and arachnid recorded from the Ludford Lane bone bed [19,20], which is correlated with the basal Pridoli (423±2.3 Ma). However, the biostratigraphic inferred ages of such early freshwater and land biotas are based on correlation with marine sections whose radiometric dating requires revision.
The fish/arthropod fauna of the Stonehaven Group, within which the Cowie Fish Bed is located, differs dramatically from those at structurally shallower levels of the Arbuthnott Group in the Midland Valley (see [21] for cross-section). Fault systems and unconformities may exist throughout the section in Fig 2 that are currently not exposed, and the Stonehaven Group may be of a separate provenance and depositional succession and structurally juxtaposed into its current position [22]. This observation leaves the significance of the dates here unaffected, but should be considered in interpretations for the structural development of the region.
Although attention has been paid to ascertaining ages to successions in Fig 2 via other means, consensus on the age of the units has not been achieved (e.g., [23], see Discussion regarding dating exposures of the Ringerike Group (Old Red Sandstone) in Norway). Many recent calibrations of the geological time scale have relied heavily on statistical techniques on conjunction with best-fit line techniques for estimating numerical ages for chronostratigraphic boundaries. For example, the Silurian time scale is based on calibrating a constrained optimization (CONOP) composite graptolite zonation to selected radiometric ages (e.g., [24]). But this, while useful, is no substitute for precise radiometric dating of ashes in fossiliferous sections if available.
Supporting information S1 Table. Description of the excel file: Detailed U-Pb LA-ICP-MS data. The first tab (Explanation) is a detailed key for the raw data presented under the second tab (Data). (XLSX) S1 Fig. Geological sketch-