Geographic bias effects on interpretations of secular trends of Hf isotope times series in zircons

The analysis of εHf(t) time series data from zircon reveals notable discrepancies based on sample type (igneous versus detrital), statistical weighting methodology, and geographic sampling bias. These differences warrant caution when interpreting data in the context of tectonic settings and the history of supercontinents. In terms of tectonic setting, accretionary orogens dominate in both sedimentary basins that host detrital zircons and in igneous zircon sources. Because of differences between the various time series, emphasis in this study is on peaks, valleys and secular trends and not on absolute εHf(t) values. Specifically, the igneous time series for εHf(t) in zircons shows more peaks and valleys than corresponding detrital time series for both weighted and unweighted data (weighting corrects for disproportionate geographic sampling). Also, sample-based (each sample considered separately) and site-based (samples grouped by geographic location and age) results align closely to the igneous time series, whereas the site-based detrital series displays more negative εHf(t) values. Regardless of the type of time series, the failure to compensate for disproportionate geographic sampling increases the prospects of producing an unrepresentative time-series. Nine zircon age peaks (both detrital and igneous) have corresponding εHf(t) peaks (3200, 2700, 2500, 2150, 1500, 1100, 750 Ma) and two have corresponding age valleys (1800–2000, 550 Ma). With exception of a geographically widespread 1500 Ma peak, most of the εHf(t) peaks and valleys are controlled by specific geographic regions and are likely not be global in extent.

Two distinct periods (200–0 and 1800–1600 Ma) display εHf(t) signatures that rise steadily for 100–200 Myr, coinciding with the final stages of supercontinent assembly and the transition to the retreat of exterior orogens. An εHf(t) peak at 750 Ma and a high at 1400–1100 Ma partly overlap with supercontinent breakup and valleys at 550 Ma and 900 Ma with supercontinent assembly. A large εHf(t) valley at 2000–1800 Ma corresponds with the onset of craton collisions that led to the final assembly of Columbia at 1800–1600 Ma. The steep rise in εHf(t) in the last 200 Myr in both igneous and detrital zircons is controlled by sites in Circum-Pacific orogens in North and South America and Southwest Asia, and it parallels the breakup of Pangea. The general increase in zircon εHf(t) in the last 500 Myr in both detrital and igneous data reflects an increase in the proportion of isotopically juvenile components in accretionary orogens.