Igneous rock area and age in continental crust

Rock quantity and age are fundamental features of Earth’s crust that pertain to many problems in geoscience. Here we combine new estimates of igneous rock area in continental crust from the Macrostrat database (https://macrostrat.org/) with a compilation of detrital zircon ages in order to investigate rock cycling and crustal growth. We find that there is little or no decrease in igneous rock area with increasing rock age. Instead, igneous rock area in North America exhibits four distinct Precambrian peaks, remains low through the Neoproterozoic, and then increases only modestly toward the recent. Peaks in Precambrian detrital zircon age frequency distributions align broadly with peaks in igneous rock area, regardless of grain depositional age. However, detrital zircon ages do underrepresent a Neoarchean peak in igneous rock area; young grains and ca. 1.1 Ga grains are also overrepresented relative to igneous area. Together, these results suggest that detrital zircon age distributions contain signatures of continental denudation and sedimentary cycling that are decoupled from the cycling of igneous source rocks. Models of continental crustal evolution that incorporate significant early increase in volume and increased sedimentation in the Phanerozoic are well supported by these data.


INTRODUCTION
Quantitative constraints on the age-varying properties of rocks in Earth's crust are critical for generating and testing hypotheses about the long-term evolution of Earth systems. A priori expectations for the quantity-age distribution of some rock types can be formulated with assumptions about how geological processes operate. For example, a fundamental prediction of the sedimentary cycle is that surviving sediment quantity should decrease exponentially with increasing age (e.g., Mackenzie and Pigott, 1981). The same principles of rock cycling apply to igneous rocks in continental crust, but models are less firmly grounded in a steadystate world view. This is because it is accepted that while today continents occupy ∼30% of Earth's surface, at some early point in Earth's history there cannot have been any continental crust. Between these two constraints, nearly all possible models have been proposed, each with different preferences for the relative importance of cycling versus time-varying production (Armstrong, 1981;Roberts and Spencer, 2015;Puetz et al., 2017;Condie et al., 2018;Dhuime et al., 2018;Condie and Aster, 2010). Resolving these models and calibrating rock cycling has implications for how we interpret deep-time records and for generating and testing hypotheses for drivers of long-term changes in Earth systems (e.g., Hayes and Waldbauer, 2006;Husson and Peters, 2018).
Several attempts have been made to estimate continent-or global-scale rock quantity so as to constrain rock cycling and crustal growth models with minimum estimates of original volume. Some are based on geological maps (e.g., Blatt and Jones, 1975;Goodwin, 1996;Wilkinson et al., 2009), the most widely produced models for the lithology and age of rocks in Earth's crust, albeit only explicitly for a surface. Studies that integrate both surface and subsurface data provide a more complete description of crustal age and composition, but most have emphasized sediments (Ronov et al., 1980;Husson and Peters, 2017). Recent advances in high-throughput zircon U-Pb geochronology and geochemistry provide a proxy for crustal growth and recycling (e.g., Cawood et al., 2013;Payne et al., 2016;Korenaga, 2018;Rosas and Korenaga, 2018;Puetz and Condie, 2020), but these methods rely on several key assumptions, including that the frequency of crystallization ages among compilations of detrital zircon (DZ) is proportional to the quantity of igneous rocks that sourced the sediment.
We leverage advances in geoinformatics in order to provide new constraints on the areaage relationship of igneous rocks in continental crust. Our study is focused on North America, where surface and subsurface data are available, but we consider this record in the context of global map data.

DATA AND METHODS
Geologic maps in the Macrostrat database (https://macrostrat.org/; Peters et al., 2018) are grouped into four scales that combine sources into coherent two-dimensional representations. Here, we use the two scales that are globally complete for continents: "tiny" (∼1:20,000,000 scale) and "small" (∼1:5,000,000 scale; Fig. 1). The tiny-scale map derives from Chorlton (2007), and the small-scale map was composited from this and other sources (see Table S1 in the Supplemental Material 1 ). All bedrock maps in Macrostrat consist of polygons for unit boundaries, each of which minimally has chronostratigraphic age(s) and lithology descriptions linked to vocabularies (see https://macrostrat. org/api/defs/).
Two map scales are considered here to illustrate the effects of different temporal binning 1 Supplemental Material. Supplemental Figures S1 and S2, and Table S1 (data used in this analysis). Please visit https://doi .org/10.1130/GEOL.S.14772795 to access the supplemental material, and contact editing@geosociety.org with any questions. schemes and rock unit definitions. To make the scales as comparable as possible, polygons were clipped to the outline of land today, and major oceanic islands were removed. Map polygons containing igneous or metaigneous rock are shown in Figure 1A and 1B, along with outlines for other rock types. For an online interactive version, see https://macrostrat.org/map/.
We also include geologic column data that summarize the lithologies and ages of rocks in the surface and subsurface regionally. Macrostrat columns are not yet global in coverage, and here we focus on 949 columns in North America (Fig. 1C). Column rock units acquire an age model that incorporates correlations to chronostratigraphic bins and relative age constraints between units within bins . Thus, the ages of rock units in columns are typically more finely resolved than in maps. Columns can also include igneous rocks of different ages and lithologies through a thickness of crust that is covered by sediment, a more volumetrically relevant representation of igneous rock quantity than that provided by maps (Fig. 1).
Area versus age was calculated for 1 m.y. increments by summing the Cartesian area in square kilometers (World Geodetic System 1984 [WGS84] spheroid) of all polygons ( Fig. 1) containing igneous and/or metaigneous rocks with an intersecting age estimate. We also include concordant U-Pb "best ages" (Spencer et al., 2016) for 69,453 DZ from 746 samples matched to 392 Phanerozoic sedimentary units in North American columns. DZ measurements derive from multiple sources, most aggregated by Puetz (2018) and all of which are accessible via Macrostrat's application programming interface (API) and included in the Supplemental Material. DZ data were not used to construct Macrostrat age models.

RESULTS
Area-age results for igneous and metaigneous rocks in Macrostrat's global geological maps and North American columns ( Fig. 1) share many similarities over 3.5 b.y. (Fig. 2). First, the absolute values of the area estimates are similar. This coincidence occurs because igneous rocks compose 19.8% and 21.9% of the total global area of the tiny-and small-scale maps, respectively, whereas columns occupy 17.1% of the total global map area. Thus, North America stripped of sediments to reveal all igneous rocks in the surface and subsurface has approximately the same total area as surface-exposed igneous rocks do globally; North America does have proportionally more igneous rock at the surface compared to globally (22.8% and 29.6% of the tiny-and small-scale map area in North America is igneous).
The more salient similarities between igneous rock area in global maps and North American columns involve temporal patterns (Fig. 2), including shared late Archean and late Paleoproterozoic peaks followed by a decrease into the Mesoproterozoic and then a smaller mid-to late Phanerozoic rise. The better temporal resolution of the small-scale map makes patterns more apparent, but results are consistent between map scales. Importantly, neither the global geological map nor North American column data exhibit a sustained increase in igneous area toward the present (Fig. 2A). To further assess this long-term trend and to address potential overweighting of poorly time-resolved rock units, we normalized the area of each igneous polygon by its estimated duration in millions of years. Such normalization does impact the temporal trajectory (Fig. 2B). Notably, the Archean and Paleoproterozoic peaks are lower and there is a large increase in area toward the recent during the Phanerozoic. Normalization by duration may, however, introduce bias by increasing the area per million years of units from extant igneous systems that will range into the future while decreasing that of igneous units that formed over an area for a long duration. Despite such distortions, and regardless of which estimate is used, there is little or no long-term decrease in igneous rock area with increasing age for most of the past 3.5 b.y. (Fig. 2).
Detrital zircon age frequencies in North America share many similarities with igneous rock area estimated from columns in the same region (Fig. 3A). This is true for combined DZ ages and for Precambrian-aged grains when they are subdivided by the depositional age of their host sediments. Notably, both igneous area and DZ ages exhibit Neoarchean and late Paleoproterozoic peaks, but their relative magnitudes are different. DZ age frequency also aligns broadly with a peak in column area at ca. 1.4 Ga, and there is a shared late Mesoproterozoic peak, albeit one that is larger and somewhat younger in DZ. Both igneous area and DZ age frequencies are low through the Neoproterozoic and Paleozoic and then increase in the Mesozoic-Cenozoic, with DZ becoming richer in grains relative to igneous rock area toward the recent (Fig. 3B).

DISCUSSION
Data on igneous-metaigneous rock area in the continental surface and subsurface combined with DZ age distributions have several implications. The simplest is that North America is sufficiently large and tectonically diverse to capture a signal of igneous rock quantity with parallels seen globally (Fig. 2). This result is consistent with the finding that the sedimentary record of North America contains a global signal (Ronov et al., 1980;Peters and Husson, 2017) and probably reflects the fact that "global tectonics" is quantitatively expressed in all such syntheses of rocks from large samples of continental crust. Nevertheless, there are clear differences between the global and regional data. For example, the small-scale global map shows a Neoproterozoic increase in igneous area whereas igneous area declines to a minimum in North America (Fig. 2), a likely signal of the Gondwanan Pan-African orogeny.
Another useful result is that both igneous rock area in North America and DZ age frequencies in the same region have similar temporal variation (Fig. 3A). Thus, both records are likely detecting the same quantity-age property of the crust. More interesting, however, are the differences in these records. Notably, Archean igneous rocks are more abundant in direct measures of quantity (Fig. 2) than suggested by DZ age frequencies. Indeed, there is an overall decrease in grain frequency relative to igneous area with increasing age (Fig. 3B). One hypothesis for this discrepancy is that zircon fertility is lower in older igneous rocks because they are more mafic on average (Moecher and Samson, 2006;Lee et al., 2016). However, the mean felsic-tomafic ratio of igneous rocks in Macrostrat columns is not markedly different in the Archean (Fig. S1 in the Supplemental Material), and major elemental data from sediments also indicate that Archean crust had a compositional diversity similar to that of modern continents (Lipp et al., 2021). Correcting DZ abundance for changes in the zircon saturation of magmas (Keller et al., 2017) does significantly increase DZ density estimates in the Archean (Fig. S2), but not enough to account for the discrepancy (Fig. 3B).
In the absence of a compositional shift in igneous source rocks, another explanation for the observed change in DZ abundance relative to igneous rock area is that DZ grains undergo attrition in a way that is much faster than that of their igneous source rocks. This is expected if older grains are more likely to have undergone metamictization and Pb loss as well as physical destruction during transport (e.g., Markwitz et al., 2017;Andersen et al., 2019). The DZ record may also be inherently overprinted by its sedimentary origin (e.g., Andersen et al., 2016), with igneous rocks yielding fewer grains per unit area as they age due to increasing isolation from active orogens and lowering of their mean elevation (Spencer et al., 2018).
On shorter time scales, a notable difference between DZ age frequency and igneous area occurs at ca. 1.1 Ga, where the well-known peak in DZ ages associated with the assembly of Rodinia (Condie and Aster, 2010) is larger and younger than the peak in igneous area. It is possible that this offset reflects error in the column age model or overweighting by columns of small but widespread intrusive rocks, such as the Mackenzie dike swarm (western Canada) spike at ca. 1260 Ma (Fig. 2). Indeed, the ca. 1.2 Ga peak in igneous area is driven primarily by mafic rocks (Fig. S1), which may account for some of the difference with DZ ages. It is also possible that the offset reflects the unroofing of zircon-rich igneous rocks from a narrow, active stretch of the eastern North American margin (Park et al., 2010) or input of grains from outside of North America.
Importantly, Precambrian-aged DZ age densities are similar relative to each other and to igneous area, no matter the depositional age of their host sediment (Fig. 3A). Post-Cambrian Paleozoic sediments do contain more Ediacaran grains, a signature that reflects the docking of Ediacaran-aged arcs and terranes along the eastern margin of Laurentia (Park et al., 2010), and Cambrian sediments are richer in Archean grains, but the overall differences are small.
The similarities in Precambrian DZ age frequencies among depositional cohorts, combined with clear evidence of early Paleozoic reburial of most if not all of the igneous source rocks in North America (Fig. 3; Peters and Gaines, 2012;Keller et al., 2019), raises the possibility that much of the Precambrian DZ grain population was introduced en masse to the surface environment during late Precambrian continental denudation. Under this model, Precambrianaged DZ grains in post-Cambrian Phanerozoic sedimentary rocks did not come predominantly from exposed igneous rocks, but rather from sediments that were recycled from the margins of a largely denuded continent. This pool of DZ was then spread back over the continent during progressive Phanerozoic reburial of the Great Unconformity, beginning with a thin blanket of predominately marine sediments that covered essentially all of Laurentia by the end-Ordovician. The possibility that a similar episode of continental denudation followed by reburial occurred during the Paleoproterozoic Keller et al., 2019) provides an intriguing, if speculative, hypothesis for the low abundance of Archean DZ. In this scenario, Archean grains in post-Cambrian sediments would have undergone two such major cycles of continental exhumation and reburial.
Finally, regardless of sampling approach or geographic scale (Fig. 2), igneous rock quantity does not decrease exponentially with increasing age, as predicted by most basic models of rock cycling. Scaling igneous area by unit duration (Fig. 2B) does produce a decrease, but only over the recent to mid-Paleozoic. This is the same time scale over which there is a large decrease in nonmarine sediment , a higher elevation, faster-cycling component of the sedimentary system that does include igneous rocks. Thus, the rock age distributions reported here provide minimum bounds on net continental crustal quantities that account for rocks in the subsurface and in more rapidly cycling geomorphic systems. Covariance in old igneous and sedimentary rock quantity and the shifting abundance of sediment in continental crust (Fig. 3C) reinforces crustal evolution as a critical driver of long-term changes in Earth  = 11,196). KD for all concordant grains is also shown. KD scaling is identical between groups but arbitrary relative to rock area. KD lines encompass 95% confidence limits from 100 bootstrap replicate samples of grains. (B) Ratio of DZ KD and saturation-corrected KD (DZ density) to North American igneous rock area; envelopes incorporate bootstrap confidence limits for KD and igneous area (from A). Inset map shows location of all DZ samples (black dots), columns (Fig. 1C), and land area today. (C) Sedimentary rock area from columns (as in A). Geologic time abbreviations are as in Figure 2. systems. Critical tests of this hypothesis and many new insights will be gained by expanding the geographic footprint and temporal acuity of surface-subsurface column data globally.

ACKNOWLEDGMENTS
Macrostrat development was supported by U.S. National Science Foundation grant EAR-1150082 and EarthCube grant ICER-1440312. C. Walton acknowledges the UK Natural Environment Research Council (NERC) and UK Research and Innovation (UKRI) for support through a NERC Doctoral training partnership (DTP) studentship, grant NE/L002507/1. O. Shorttle acknowledges support from NERC grants NE/T012633/1 and NE/T00696X/1. We thank Alan Carroll and the Adam Maloof Lab at Princeton University (USA) for discussion, and Christopher Spencer, Justin Payne, and an anonymous reviewer for insightful reviews.