Potassium feldspar present in global mineral aerosol (<5%) plays a disproportionate role in modulating the microphysics of mixed-phase cloud. Via exceptional ice nucleation properties, it is capable of changing cloud properties and behaviour. Here we identify times of substantial and abrupt change in the global availability of potassium feldspar since 600 Ma. Normally, weathering and vegetation cover contribute to low availability, with clay dominating mineral aerosol. Periods of maximum availability are reasoned to follow the emplacement and remobilization of ejecta blankets from major meteorite impact events, before returning to background after some hundreds to thousands of years. We review the 44 largest confirmed craters and evaluate the potassium feldspar content of their target rocks, which range from c. 0 to >30%. By combining crater size and tectonic reconstructions, we are able to provide a quantitative and self-consistent assessment of changes to global potassium feldspar availability. Considerable differences in potassium feldspar availability following meteorite impact events are revealed. Different impact events generated dust containing different amounts of potassium feldspar. Differing levels of influence upon climate are hypothesized, and should now be tested by looking at stratigraphic records of these events to reveal the sensitivity of climate to different dust mineralogy.
Supplementary material: Figures showing palaeolatitude determinations using various methods and reconstruction models, and estimates of alkali feldspar abundance and distribution across contemporary continental landmasses for meteorite impact events are available at: https://doi.org/10.6084/m9.figshare.c.4253312
Cloud dominates albedo and controls the distribution of water and sunshine across the Earth. Even subtle changes to cloud influence the wider Earth system. For example, sub-per cent per year changes in cloud cover over Greenland control the mass of its ice sheet by mediating the primary addition rate (Hofer et al. 2017). Variations in optical properties and behaviour of cloud are integral to weather, the hydrological cycle and climate, yet they remain difficult to parameterize and predict (Tan et al. 2016; Storelvmo 2017; Xie et al. 2017).
Ice nucleation is a microphysical process that dramatically changes optical properties (brightness, affecting radiative transfer) and behaviour (shortens lifetime by promoting rain) of cloud (Lohmann & Feichter 2005; Hoose & Möhler 2012). Whereas liquid water can supercool to −35°C in the atmosphere, ice formation is triggered at much warmer temperatures on the surfaces of aerosols acting as ice nucleating particles, INPs (Murray et al. 2012; Nenes et al. 2014; Herbert et al. 2015; Vali et al. 2015; Koop & Murray 2016). The amount of INPs in the atmosphere is critical to cloud behaviour and climate (Vergara-Temprado et al. 2018), as is the efficiency of those INPs.
Mineral dust accounts for around half the total aerosols in the modern atmosphere (Textor et al. 2006). This aerosol is composed mainly of clay (Atkinson et al. 2013) emitted from arid areas of the Earth (Gieré & Querol 2010). Mineral aerosol has been found to account for a large proportion of all ice nuclei in Earth's modern atmosphere, yet clay is poor at nucleating ice (Atkinson et al. 2013). In the past 5 years, research first identified feldspar as having exceptional ice nucleating properties (Atkinson et al. 2013). Further investigation revealed that only the alkali feldspars, not plagioclase, exhibit these properties (Zolles et al. 2015; Harrison et al. 2016; Peckhaus et al. 2016).
Of the alkali feldspars (Afs; mineral abbreviations from Siivola & Schmid (2007) and Whitney & Evans (2010) are used throughout), more potassic samples appeared, overall, to be better at nucleating ice than sodic (Harrison et al. 2016). Recently, an explanation as to why has been reported by Whale et al. (2017): the most efficient Afs are those that exhibit perthitic texture whose nanoscale topography, which includes etch pits and tunnels as a result of solution–feldspar reactions, probably acts to promote ice nucleation by providing suitable settings at the scale of a critical ice nucleus (the minimum cluster size of water molecules required for the spontaneous growth of ice) (Whale et al. 2017). Perthite is potassium feldspar (Kfs) with exsolved albitic (Ab) lamellae, and is prevalent in granitic rock types (see Parsons et al. (2015), to whom we direct the reader for a recent review of alkali feldspar microtextures). Consequently, Afs, and particularly perthitic Kfs, is considered to play a role in the global atmospheric and hydrological cycle that is highly disproportionate to its abundance. It is therefore Kfs that we focus our attention upon here.
Although laboratory-derived knowledge of specific mineral ice nucleation is growing (see Kiselev et al. 2017; Whale et al. 2017), determining real effects is challenging (Storelvmo 2017). This is partly because cloud properties are difficult to measure directly in real time (Tan et al. 2016). Furthermore, cloud properties are impossible to measure, at least directly, in the deep past.
Knowledge of the abundance and type of aerosols in the past is therefore important to improving understanding of modern and future atmosphere dynamics (see Buseck & Pósfái 1999; Andreae 2007), because it is the indirect effects of aerosols (i.e. via changes to cloud properties and behaviour; Komurcu et al. 2014) that contribute the largest uncertainties in global climate models (Stocker 2014). Currently, baseline knowledge concerning pre-industrial aerosols is poor (Carslaw et al. 2017).
Is a proxy record of ice nucleation regime possible?
Pankhurst (2017) provided an Earth-history scale overview of the availability of potassium feldspar (Kfs), to form the first step towards such a proxy record for an ice nucleation regime. No attempt was made to review mineral texture, or chemistry, as more work is required to determine links between ice nucleation efficiency, chemistry and texture of the feldspar group. It was simply reasoned that formation of granite and even more potassic intrusive and metamorphic equivalents (syenite, alkali feldspar granite, gneiss) fundamentally controls Kfs crustal distribution on a global scale, and that this lends the widest possible context for future study of links between mineral dust and cloud in the past.
Key points from Pankhurst (2017) include the following: (1) the accumulation of granites during construction and evolution of continental crust has increased Kfs primary availability, yet processes such as soil formation (feldspar to clay) and changing vegetative or ice cover probably play a more important role in modulating how much Kfs is able to be aerosolized on shorter timescales; (2) it is demonstrated that the availability of Kfs and its contribution to aerosol has not been constant across the averaged Earth's surface in time; (3) dust leaves a physical record (e.g. loess and other aeolian deposits) in which mineralogy can be measured; (4) a Kfs factor (KFF), which is simply the volumetric per cent abundance (i.e. pure Kfs has a KFF of 100), can be assigned to atmospheric dust, current land surfaces, rocks and aeolian deposits, which provides a shorthand nomenclature for linking multidisciplinary study.
Climate forcers such as pCO2 (Booth et al. 2017) and sulphate aerosol (Langmann 2014) have records composed of both historical measurements and often much longer term proxies (Booth et al. 2017; Wolfe et al. 2017). When linked to effects via laboratory-derived understanding and natural evidence, such proxy records provide the basis for model validation (Schmidt et al. 2011; Booth et al. 2017). Because the origins of dust in today's atmosphere can be fingerprinted to specific regions (Zhao et al. 2018), and a KFF can be derived from each (Pankhurst 2017), it becomes possible to make predictions about the effects of specific mineral dust events.
However, if the importance of Kfs to cloud behaviour is to be addressed with the aid of a pre-industrial and potentially much longer proxy record, what is needed first is to identify the best possible time periods in the rock record for multidisciplinary study. Owing to the low concentration of aerosol particles under pre-industrial conditions, global mean cloud albedo during much of Earth history could have been twice as sensitive to changes in aerosol emission as it is today (Carslaw et al. 2017). As such, any effects of changing mineral aerosol volume or variety may be more pronounced, and therefore evident.
Targeting the largest dust-refreshing events through time
The mineralogy of dust sources varies considerably across the planet (e.g. Scheuvens et al. 2013; Smith et al. 2014), yet in normal circumstances clay minerals dominate the global mineral aerosol budget (Atkinson et al. 2013) owing to the continual process of chemical weathering (White et al. 2001). Weathering (i.e. from feldspars to clays) forms a thin barrier between ‘fresh rock’ and ‘fresh air’. This regime has been in place since at least the Neoproterozoic (Eriksson et al. 2013), and is reasoned to both suppress and buffer mixed-phase cloud ice-nucleation efficiency in the Earth system (Pankhurst 2017). Changing this regime requires fresh minerals to cover the clay barrier, and therefore be available to become aerosolized. To affect global aerosol mineralogy, this natural barrier must be overcome rapidly and dramatically.
Meteorite impact events do just this. The rock and mineral fragments excavated from deep in the crust during major crater formation are distributed as ejecta (Collins et al. 2005; Osinski et al. 2011, 2012), blanketing the pre-existing strata. ‘Impact winter’ describes a scenario in which the sheer volume of dust in the atmosphere blocks sunlight and causes drastic cooling of land surfaces (Covey et al. 1994). This means that among other, longer-term effects (see Walkden & Parker 2008) we note that the potential to change the mineralogy available to the atmosphere exists after the dust has settled, until weathering and removal returns the surface veneer to the normal clay-rich state.
On Earth, Kfs is spatially heterogeneous owing to a range of tectonic, magmatic and land surface processes that control the distribution of typical upper crustal rock types (Pankhurst 2017; see Fig. 1). Underneath Earth's surface veneer, Kfs and the textures it forms are heterogeneously distributed (as are all minerals). Some rock types such as basalt and limestone contain ≪1% Kfs, whereas some granitoids can exceed 60%. Via the excavation of local-scale rocks from the Earth's crust and redistribution across a global-scale surface, major meteorite impact events hold the potential to change the Kfs content of dust available to the atmosphere in a geological instant. Has this ever occurred, and are there strikes that produced similar amounts of dust, but without Kfs?
Constructing a database of dust-refreshing events
Our aim is to collate the available data concerning scale and Kfs content of ejecta from meteorite impact events in a manner that allows first-order comparisons to be made. We begin with the most complete resource available, the Earth Impact Database (PASSC 2017), and later discuss inherent limitations to this approach. The largest crater diameter according to the Earth Impact Database is 160 km (Vredefort). All the confirmed bolide impacts with collapsed crater diameter >16 km were collated and for the purposes of this work are considered ‘major’ as they represent the top order of magnitude of crater diameters (D). From these major strikes a subset from 600 to 0 Ma was selected (n = 43), to which we added the Alamo impact event (D = 65 km; Pinto & Warme 2007). We restricted our time window to exclude times >600 Ma because the post-600 Ma geological record is more complete, tectonic reconstructions are more accurate, it is the period of Earth history in which soil formation and vegetative cover predominate, and this interval afforded us a dataset of useful size. These 44 impacts were assessed for (1) their target rocks’ KFF, and (2) the volume of the ejecta blanket resulting from the impactor. For craters with age determinations better than ±20 myr (n = 35), we were able to determine (3) their palaeogeographical location at the time of impact, which, when combined with (4) the lateral spread, thickness and particle sizes of the ejecta blanket, allowed us to estimate (5) the post-strike KFF of the global average land surface. We chose a cut-off of ±20 myr to limit the uncertainty produced by tectonic reconstructions and still retain c. 75% of the database.
Assessment of Kfs abundance in meteorite impact event target rocks and average global surface veneer
The level of detail available concerning mineralogy of target rocks in the literature is highly variable. Kfs abundance, if not directly measured, was determined from target rock descriptions. A robust assessment of Kfs abundance can even be made from field descriptions. This is because Kfs is (1) absent or negligible in many major rock types (e.g. basalt, limestone), and (2) a diagnostic rock-forming mineral used in quantitative terms within the nomenclature of granitoid suites (e.g. granitic or granodioritic; high and low Kfs respectively; see Streckeisen 1974). Where available, whole-rock geochemical data from granitoids (or metamorphic equivalents) were used to generate CIPW norms.
Terrigenous sediments are composed of between 5 and 15% feldspar, of which Kfs has been regarded as the most common owing to its higher availability and stability than plagioclases (Pl) under normal surface conditions (Folk 1980). The largest single continent-scale dataset of soil mineralogy gathered to date (Smith et al. 2014), however, returns an average Kfs:Pl of 1:1.7. Authigenic Kfs is common in marine sediments, but only in rare cases do abundances exceed 2% (Buyce & Friedman 1975). We assumed a low total Kfs content in sedimentary target rocks where mineralogical data were not available, as high Kfs content is reported only from immature sediments derived from granites or felsic large igneous provinces (e.g. Eriksson et al. 1994). Furthermore, the sandstone packages reviewed here are often intercalated with clay- and carbonate-dominated rock types, which dilute Kfs content in the overall sequence.
The most complete overview of meteorite impact event target rocks to date (Osinski et al. 2007) demonstrates that target geology ranges from highly homogeneous to heterogeneous. Volumes of individual rock types, and by extension the Kfs content, can be estimated with accuracy if outcrop maps (e.g. Lambert 2010) and/or subsurface observations, such as from geophysical surveys and drillcores (e.g. Tsikalas et al. 1998; Bron & Gostin 2012) are available. The corollary is that poor exposure and limited subsurface information leads to a broader range of plausible Kfs content.
A summarized description of target rocks and minerals from each crater is accompanied by an assessment of KFF in Table 1. Because the level of detailed lithological description available in the literature is variable and cannot be presented with equal confidence, the abundances are reported as a range of plausible values rather than a value with an uncertainty. Descriptions characterized by complex, heterogeneous lithologies, but without well-understood proportions, resulted in a larger range of plausible KFF. Lithology of target rocks can be matched to their ejecta blankets where they are sufficiently preserved (Gostin et al. 1986; Horton & Izett 2005; Thackrey et al. 2009), which validates this approach.
A global average KFF of five is derived from the most exhaustive database of soil mineralogy published to date (Smith et al. 2014). These soil data are used as the best available mineralogical surface veneer dataset for an entire, continental-sized landmass with a long and mixed geological history; as such they are the most representative available. A global atmospheric aerosol KFF of three (derived from Atkinson et al. 2013) is consistent with this value, given Kfs is present mainly in the silt, rather than clay size fraction (Claquin et al. 1999), and is more likely to be under-represented in aerosolized soil.
Ejecta blanket calculations
The scale of ejecta blankets is fundamentally linked to the volume of the transient cavity formed during the first stages of impact cratering, which is difficult to calculate with certainty. All the impacts considered here formed complex craters (including collapse of the transient cavity inwards), and have undergone variable modification and erosion post-impact. Owing to the considerable sources of uncertainty in any single case (see section below) and to maintain equal treatment, in the first instance we used the measured crater value of D to approximate the diameter of the transient cavity (Dtc) by a fixed ratio of 1:1.75. This value is based on detailed analysis of evidence from Chicxulub (Morgan et al. 2002), which, as its target rocks were sedimentary cover over crystalline basement, is both the most studied and most typical geological scenario in our database.
Equations within the Earth Impacts Effects program http://impact.ese.ic.ac.uk/ImpactEffects/ (Collins et al. 2005) were used to calculate the diameter and depth of the transient cavity (following terminology of French 1998), and ejecta thickness with distance for the subset of craters where appropriate (see palaeogeographical reconstructions and area calculations below) . To check if our use of a fixed ratio caused significant differences in model outputs for each impact, parameters were fixed with geological evidence where available (e.g. target rock type and thus density and presence of water). In all instances where the following parameters are not available, an impact trajectory of 45° (statistically most likely), bolide density of 3000 kg m−3 (dense rock) and pre-atmospheric entry velocity of 17 km s−1 were used. For a fair comparison between these events, we did not attempt to reconcile strewn fields or measured ejecta blanket thickness with our calculations, for these are subject to post-deposition modification, and such measurements are limited to a minority of these events.
The Dtc value given by consistent application of the 1:1.75 relationship taken from the largest case (Chicxulub, where D = 150 km) is normally an underestimation of the Dtc given by the program. In the smallest cases (i.e. Ames, D = 16 km) Dtc is underestimated by <20%. For the purposes of comparing Kfs availability post-strike between events, larger errors at the smaller end of the range of crater sizes are more acceptable than the opposite. Transient cavity depth is especially important for craters exhibiting pre-strike strata with differing Kfs abundance with depth. For consistency, we maintain use of the Dtc from the ratio of 1:1.75 described above. Around one-third to a half of the fragmented material either does not escape or falls back into the crater (Collins et al. 2005; Osinski et al. 2012); up to two-thirds could escape. As such, ejecta are dominated by material from the excavated zone as opposed to the (deeper) displaced zone (Osinski et al. 2011), yet rocks across the cavity depth can be found in ejecta (Osinski et al. 2008, 2012). It remains poorly known how much material is actually ejected and this may vary considerably on a case-by-case basis (Melosh 1989; Osinski et al. 2011). In our calculations of ejecta volume we use one-half of the transient cavity volume, to again maintain an underestimated final result.
In examples where thick, low-KFF sedimentary cover overlies high-KFF basement, we assume material from only the top half of the transient cavity contributed to the ejecta, in the absence of mass-balance calculations from the impacted rocks themselves (e.g. Hörz 1982). This step prevents the gross overestimation of Kfs, which owing to the normal low background would be more erroneous than underestimating by a little. This step produces differences in Kfs availability only for craters with diameter at the smaller end of the range, and applies to <20% of the dataset. Furthermore, because smaller craters produce less ejecta and so inherently have a smaller capacity to change the KFF of the global surface veneer, no first-order bias is anticipated in the later comparisons.
Melting has the potential to reduce the amount of Kfs released to the ejecta blanket. The proportion of melt volume to ejecta volume increases with energy transfer, which is a function of impactor size and velocity (Cintala & Grieve 1998). Even for the largest strike considered here, however, this melt volume is <10% of the transient cavity volume, and does not change substantially if comet densities (much lower) and velocities (potentially much higher) are used. It is worth noting that the processes of melt generation and emplacement are more complicated than first thought (Osinski et al. 2011). For instance, the El'gygytgyn impact crater provides evidence that more extensive melting can occur (Pittarello & Koeberl 2013a) and so illustrates the potential for uncertainty. Because average ejecta thickness with distance is more meaningful to global Kfs availability (see below) than the comparatively small difference that target rock melting of this scale would make to the ejecta blanket KFF, we do not include impact melt volume in our assessment.
Understanding of impact cratering is incomplete, and all individual model outputs should be considered highly uncertain. Nevertheless, the calculations of ejecta volume (see Table 2) are most sensitive to the transient cavity volume. Applying the same set of scaling laws (e.g. Croft 1984; Holsapple 1993, 2017; Collins et al. 2005) consistently across the largest 44 known impacts present in the rock record results in large relative differences in KFF and volume.
Dust volume calculations
Particles in the <100–200 µm range are mobilized by desert winds on ballistic trajectories close to the surface, and break apart or eject smaller particles on striking the soil (Mahowald et al. 2014). These smaller particles (<c. 50 µm) are able to be entrained into the upper boundary layers, after which they can be transported long distances and enter the atmosphere (Mahowald et al. 2014). Hypervelocity impacts produce proportionally more fine material than surface or subterranean explosions (O'Keefe & Ahrens 1985). According to a regression based on scaled laboratory experiments, c. 10% of the fragmented ejecta would be <200 µm (O'Keefe & Ahrens 1985). If the volume of the ejecta blanket is 50% of the transient crater volume (as discussed above), it follows that the available volume of primary dust meaningful to the climate system immediately post-impact can be estimated as 5% of the transient cavity volume. This estimate of total dust production is presented for each impact event relative to modern dust emission as the lower horizontal axis in Figure 1c. This value of 5% is considered highly uncertain in any individual case, but does serve as the most reasonable estimate at present. The volume and flux of dust mobilized from the ejecta layer as part of normal aerosolization mechanisms is another flux that is difficult to estimate, which is why we consider the fraction of continental area to be useful (see below).
Palaeogeographical reconstruction and area calculations
Palaeogeographical context of a meteorite impact event has the potential to be as important as target rock KFF. For instance, if an ejecta blanket of size X composed of 20% Kfs spanned a large continental area (above sea level), the available Kfs across the surface of the Earth would be comparable with that caused by an ejecta blanket of size X with 40% Kfs that straddled a coastline. Approximately half the material (and Kfs) in this latter case would be deposited directly into the ocean and take no part in later aerosolization processes.
Palaeolatitude was estimated for the majority of meteorite impact events using the calculator of van Hinsbergen et al. (2015), and compared with both palaeolatitude and palaeolongitude estimated using geological maps and plate reconstruction models (Matthews et al. 2016, strikes <300 Ma; Scotese 2016, >300 Ma; see Fig. 2). This allowed both an estimate of ejecta blanket reach around the globe, and palaeogeographical context (ejecta coverage of continental v. ocean area) using GPlates (https://www.gplates.org/; Gurnis et al. 2012). Small circles were drawn in GPlates with radii corresponding to average ejecta thickness, derived from the same model as used to calculate the transient cavity volume (Collins et al. 2005).
To calculate the proportion of the Earth's continental area covered by ejecta blanket, 8-bit binary images (255: continents; 0: oceans) using an equal-area projection were exported from GPlates at high resolution (10 000 pixels width) and analysed for area using ImageJ (https://imagej.nih.gov/ij/; Rasband 2015). Images at the same pixel resolution were exported with coloured lines (c. 3 pixels wide, value between 0 and 255), to which a threshold was applied and deleted. This removal of the fine lines isolated regions of continent inside and outside the relevant circles (ejecta minimum thickness of 10 cm to 10 µm). The area of each region was measured, and used to calculate the area within each isopach. This value was then expressed as a proportion of the total continental area. Using this approach, a self-consistent prediction of global Kfs abundance immediately after a meteorite impact event was derived.
Current limitations and future considerations
The Earth Impact Database, at present, is not itself peer reviewed, but rather collates information from the ‘hard’ literature and seeks to present it in a consistent and accessible format. Debate concerning accuracy of some of the information listed is currently continuing. For instance, recent refinement to a few of the crater diameter data (Osinski & Ferrière 2016) suggests that improvements to our results can be expected in future studies. Nevertheless, for our purposes here (see section ‘Constructing a database of dust-refreshing events’) it is necessary to consider, and treat, a large number of impact events consistently.
The primary amount and distribution of dust available to be remobilized after emplacement of an ejecta blanket calculated here are gross oversimplifications in each case. We supply the outputs attached with the following caveats.
The volumes calculated are internally consistent and are most sensitive to a parameter (transient cavity diameter) that has also been calculated and presented consistently (Earth Impact Database, PASSC 2017), yet these values include their own individual uncertainties that are not considered here.
The calculations of ejecta thickness with distance from a crater (Collins et al. 2005) were applied consistently and are helpful for visualizing and comparing the scale of events. However, the isopachs are not considered to be spatially precise, but rather approximate an average value of thickness at that distance (see Collins et al. 2005). These maps provide for an internally consistent method to assess continental v. oceanic distribution, which allows us to attach proportionately more relevance to regions closer to a crater than further away from it.
This relevance is calculated by using plate reconstruction models, which also include uncertainty we do not explicitly factor in. However, the comparison between independently derived palaeolatitudes (Fig. S1) lends confidence that these models have sufficiently small uncertainties for our purposes. Even when combined with the caveat attached to isopachs above, at the very least we can safely assume our approach introduces far more accuracy to our final global Kfs availability calculation than it removes.
We note that a portion of the continents as defined by GPlates polygons was subaqueous at each time considered. This portion is not estimated because of large differences in uncertainties across the 600 myr considered. The proportion of subaerial continental area covered by ejecta into continental interiors is likely to be underestimated slightly as a result, as shallow sea and shelf areas are typical of continental margins. This error is probably only a few per cent from small strikes and less for larger strikes. The possibility of simultaneous meteorite impact events (e.g. Spray et al. 1998), which would complicate these calculations, is low in comparison with the size of the dataset, and is not considered here.
The KFF of lithologies hosting the impact crater is considered a good guide to the KFF in the associated ejecta blanket, an assumption supported by studies that successfully link craters with their blankets using shared lithology and mineralogy (Hörz 1982; Gostin & Zbik 1999; Thackrey et al. 2009). However, determining an accurate abundance is difficult without detailed lithological maps, mineralogical notes and an assessment of melt percentage in distal ejecta. Presenting a range of plausible KFF values rather than a single value with an uncertainty is considered to better reflect this limitation.
Aerosolization of the ejecta blanket and return to background
Aerosolization potential (see Shao et al. 2011) after primary deposition is governed by landscape processes that are not within the scope of this work. What kind of primary spatial distribution is actually meaningful, and estimating the timescales of return to background Kfs availability, could be highly sensitive to climatic conditions at the time of individual meteorite impact events. Therefore how an instantaneous, voluminous dust source is redistributed after its primary deposition compares with the present dust cycle is unknown, and difficult to estimate owing to numerous feedbacks (Shao et al. 2011).
The rock record
The meteorite impact event record is incomplete, although all craters still present (not subducted) on Earth >85 km in diameter are likely to have been recognized by now (Kenkmann et al. 2015). We note that there are examples of low, average and high KFF across the range of crater size. Because much smaller craters than 85 km are calculated to effect considerable changes to global KFF, this missing record must be taken into account by future efforts to compare this parameter against other data in the rock record. Finally, the record is more complete with decreasing age; any real relationship observed between surface veneer KFF and other data should become stronger through time to the present.
We include only meteorite impact events that are reported with an age determination of ±20 myr in our assessments. This threshold provides reasonable accuracy when linking time of event to tectonic context. Better age-dating of impact events will lead to both (1) more events being able to be considered for Kfs availability and (2) increased confidence in any interpretations arising.
The pre-strike distribution of Kfs across the Earth is also not considered in this study. The possibility remains that, for instance, high-Kfs dust generated by a meteorite impact event covered land surfaces containing zero Kfs or very high-Kfs dust, meaning the upsurge calculated is under- or overestimated for each scenario respectively. This potential source of uncertainty becomes important for small strikes generating extremely high or low Kfs dust.
Results: target rocks, scale of ejecta blankets and primary global distribution
KFFs typical of major lithologies (Fig. 1a), and those estimated in ejecta from the most relevant subset of Earth's cratering record (Fig. 1c; n = 44) are presented alongside that of modern soil (Fig. 1b), which we use as the best available guide to the normal regime (see also Table 1). There is no systematic relationship between the scale of meteorite impact event and KFF of target lithology(ies) (Fig. 1c), which is unsurprising because of the randomized location of strikes and heterogeneous distribution of Kfs. Sixteen meteorite impact events mobilized low-KFF rocks. Half of the strikes mobilized Kfs to at least greater than the mean global surface value. Around a quarter contain Kfs that probably exceeded average upper continental crust (Taylor & McLennan 1995). Targets with mixed lithologies generally resulted in estimated KFF at or slightly higher than average surface availability.
During desert wind events, particles up to 100–200 µm are broken up into smaller fragments (<c. 50 µm) by saltation, before being entrained and carried long distances in the atmosphere (Mahowald et al. 2014). Because the proportion of ejecta fragments generated by a hypervelocity impact <200 µm can be estimated (c. 10 vol%; O'Keefe & Ahrens 1985), simple comparisons of Kfs availability to the atmosphere can be made between normal conditions and those after individual meteorite impact events.
Selected palaeogeographical reconstructions of ejecta distribution, a guide to dust generation, are presented in Figure 2 (all are provided in supplementary material Fig. S2). The assessment of ejecta distribution across continent v. ocean (Fig. 2 and Fig. S2) allows a quantitative estimate for the globally averaged surface veneer KFF to be calculated, and is plotted in Figure 3 as a timeline since 600 Ma and using KFF = 5 as a baseline (Pankhurst 2017). Each post-strike value is simply the weighted average of the KFF of continental surface area inside the 10 µm isopach (the KFF of the target rocks), and outside it (KFF = 5). We assume the area within the 10 µm isopach represents a meaningful, and at least internally consistent, approximation of the surficial extent of particles that can be readily aerosolized shortly after primary deposition.
This is because, although smaller particles (ejected over a larger area) are more easily mobilized, the ejecta blanket's capacity to act as a continuing dust source increases towards the strike location (thicker blanket, yet smaller area). Our calculations, which are conservative at each step, show that many meteorite impact events generated a volume of particles <200 µm between two and three orders of magnitude greater than that of modern annual dust emission by mass (Table 2; also shown in Fig. 1), although in many instances a significant portion would probably be deposited directly in the oceans (Fig. 2 and supplementary material Fig. S2). Owing to the global scale of ejecta blankets, we gauge them to represent highly important instantaneous dust sources, whose persistence is discussed below.
Toward understanding alkali feldspars' role in the climate system
There are many areas for improvement in both accuracy and precision of our calculations of Kfs availability to the atmosphere. However, the review of target rocks shows that KFF ranges from zero to >60. The volumetric estimations of the dust portion of ejecta blankets often exceed that of modern annual dust emission (1–3 Pg; Astitha et al. 2012) by two to three orders of magnitude (Table 2). Improved assessment of crater diameter (Osinski & Ferrière 2016) leading to changes of c. 10–20 KFF do influence our results, although considerably so only at the higher end of the range.
For instance, using Osinski & Ferrière's (2016) data, ejecta volume from Manicouagan would be c. 50% larger, although it would have a minimal influence on globally averaged surface veneer KFF. Chicxulub ejecta volume would be closer to double, yet its influence already reaches around the globe and in this case the most important factor becomes the KFF of the target rocks. Therefore, and despite the caveats detailed above, our attempt at quantifying Kfs availability demonstrates that a considerable range in both absolute abundance and KFF existed in ejecta blankets.
Context for future studies
Kfs is common in the Earth's upper crust. It accounts for c. 20.1 wt% (Taylor & McLennan 1995). Kfs is the key constituent of many granitoid types (in particular, syenite, alkali feldspar granite and granite sensu stricto) and some extrusive equivalents, which are heterogeneously distributed through tectonic and magmatic processes (Pankhurst et al. 2011a, 2011b, 2013). The distribution of different rock types across the continents results in major differences in the KFF of dust sources at local scales (Scheuvens et al. 2013; Smith et al. 2014).
Arid areas contribute a disproportionate amount of both total dust (Marticorena 2014) and Kfs owing to the low soil surface stability and lower rates of chemical weathering respectively. KFF in dust has been measured between 0.18 and 25, which reflects local upwind dust sources (Kandler et al. 2009, 2011; Díaz-Hernández et al. 2011). On average, however, Kfs contributes c. 3 wt% to the global mineral aerosol budget, and becomes comparatively well mixed in the atmosphere (Atkinson et al. 2013). This means that the mineralogy of atmospheric dust reflects an averaged global surface veneer mineralogy.
Background Kfs availability to the atmosphere, at a global scale, is buffered at low levels owing to chemical weathering and heterogeneous distribution of both high- and low-Kfs rock types (Pankhurst 2017). Instances of zero levels of atmospheric Kfs, such as prior to the first true granite formation and instance(s) of Snowball Earth, follow a logic that can be reasoned without any direct measurement (Pankhurst 2017). Placing upper limits on Kfs availability is more challenging to attempt with any useful degree of accuracy. However, we show here that global-scale dust-producing events can be quantified, and time-stamped, which implies that a way forward is possible.
For example, the emplacement of felsic large igneous provinces (Pankhurst et al. 2011a), particularly those typified by extensive, shallow-emplaced and highly evolved magmas, has the potential to cause at least regional spikes in Kfs availability (Pankhurst 2017). Glacial periods are linked to greater atmospheric dustiness (Fuhrer et al. 1999; Lambert et al. 2008; Vallelonga & Svensson 2014). Rapid ice-sheet retreat has the potential to increase the abundance and proportion of Kfs available to the atmosphere by exposing fresh, finely milled rock powder (Pankhurst 2017). Compared with the timescales of these terrestrial processes, however, meteorite impact events are instantaneous, and have the potential to effect more dramatic change in surface KFF by distributing local rocks over global scales.
Multidisciplinary study is needed
It is reasonable to hypothesize that global cloud behaviour, and by extension climate, is sensitive to Kfs availability. Attributing measurable changes to the scale of possible impacts on weather and climate this may have, however, is non-trivial (Storelvmo et al. 2011; Boucher et al. 2013; Yang et al. 2015; Tan et al. 2016; Storelvmo 2017). This is partly because globally averaged land surface KFF today is neither at a high nor a low with respect to the rest of Earth history (Pankhurst 2017). The rock record preserves an abundance of information that detailed multidisciplinary study can extract, and use to reconstruct palaeoclimate dynamics (e.g. Shields et al. 2003; Trela et al. 2018). It is possible that Kfs availability is important in palaeoclimate dynamics, yet recognition has not occurred, as the importance of Kfs to cloud nucleation was discovered only recently (Atkinson et al. 2013; Harrison et al. 2016; Whale et al. 2017). Comparatively well-dated examples of extreme availability of Kfs should be useful in testing hypotheses concerning any relationship between changes to atmospheric Kfs and climate.
Using bolide impacts to calculate time-stamped Kfs availability
Owing to atmospheric mixing, natural mineralogical variation in the atmosphere that allows effects to be linked and quantified is better captured by regional patterns within similar climatic conditions (Glaccum & Prospero 1980; Prospero et al. 1981; Jeong 2008; Atkinson et al. 2013). Yet this global mixing also raises an opportunity. Any sudden temporal changes to global surface mineralogy should be reflected in global aerosol mineralogy and represent a chance to gauge potential effects. The absolute age of a meteorite impact event need not be known precisely for these studies to be conducted. Strata above and below ejecta layers contain relative age information inherently linked to proxies of environmental conditions, which is more important than absolute age.
The scale of meteorite impact events is of equal importance to the KFF of target rocks when calculating global Kfs availability. For example, Manicouagan, Tookoonooka, Alamo and Kara are in the top 10 craters for size, yet cannot be said to have changed global Kfs availability considerably. Palaeo-continental context is also important (Fig. 2 and supplementary material Fig. S2). Despite Kara-Kul (D = 52) mobilizing 10× the ejecta volume of Rochechouart (D = 23 km), both are calculated to have raised surface veneer KFF to 8–9. This is due to c. 40% of Kara-Kul's ejecta blanket being deposited across ocean, and Rochechourt containing c. 2–3× the KFF in its target rocks (see Fig. S2).
Granites appear to be the key
Understanding why Kfs ice nucleation is so efficient under mixed-phase cloud conditions has been advanced by identification of active ice nucleation sites (Zolles et al. 2015), ab-initio modelling (Pedevilla et al. 2016) and direct observation of ice formation on Kfs surfaces (Kiselev et al. 2017). Crystal topography related to phase separation has recently been linked to exceptional ice nucleating properties (Whale et al. 2017). This helps explain why samples of sanidine (high-temperature Kfs polymorph) in early experiments did not nucleate ice well, and potentially why some samples of albite (anti-perthite; the Ab is the host of Kfs lamellae) are excellent ice nucleators (Whale et al. 2017).
Sanidine is common in volcanic rocks with similar bulk chemistry to granites, often observed in explosive super-eruption deposits where they are quenched from magmatic temperatures (e.g. Christiansen 2001). Based on present evidence it appears that volcanic eruptions, although overcoming the clay barrier, are less likely to be a source of alkali feldspar with highly efficient ice nucleation properties because rapid cooling rates make them less likely to form perthitic varieties of Kfs. In contrast, the Kfs reported from the rock record and discussed here is overwhelmingly composed of low(er) temperature, and slower forming polymorphs such as orthoclase and microcline, which are constituents of granite (sensu stricto) (see Table 1).
Timescales of potential impact
The actions of weathering, wash-out and burial must have returned Kfs to background levels after an upsurge in availability. Timescales depend upon regional land surface processes, which are difficult to estimate with enough accuracy to warrant comparison between events identified here. However, because modern annual dust recycling is estimated at c. 75% (Shao et al. 2011), meaning c. <1 Pg a−1 is transported to the oceans, a return to background levels of Kfs after being changed by major impact ejecta blankets may take place over thousands of years. These timescales are able to be resolved from some stratigraphic records (Galeotti et al. 2004).
What might be the effect of a ‘Kfs-rich’ atmospheric dust regime?
Increased ice nucleation efficiency would have different influences depending upon whether brighter clouds on average or shorter-lived clouds on average have more influence on the Earth's climate at the time of the event. Higher radiative transfer could lead to cooling, should the amount of cloud remain similar, or at least retain its role as the dominant albedo modulator. Shorter-lived clouds might decrease the average flux of water from oceans to the interiors of continents, and thereby also decrease the total amount of average cloud. Environments during times of super-continent tectonic regime would be expected to be more sensitive to this particular scenario of global hydrological change.
These two possible effects of increased ice nucleation efficiency on cloud properties and behaviour raise an intriguing set of tests that can be carried out when comparing similar-sized meteorite impact events in similar tectonic regimes, but different Kfs availability post-impact. Our study allows us to hypothesize that major environmental disruption via the hydrosphere is expected to have occurred in the years to thousands of years after a change to a high-Kfs regime, in addition to short-lived ‘impact winter’ effects expected via dust in the atmosphere. In contrast, similar-sized dust-refreshing events without Kfs are expected to have influences only by those mechanisms attached to impact winter. Hydrological effects are tractable and quantifiable using observations from sedimentary strata, and magnitude and type of disruption may possibly be determined using the palaeontological record.
Major meteorite impact events in the past 600 myr have mobilized instantaneous dust volumes exceeding modern annual global dust emission by 1–3 orders of magnitude. Emissions from the resulting ejecta blankets probably dominated atmospheric mineralogy for hundreds to thousands of years after the impact.
A subset of these ejecta blankets were high in Kfs; they hold the greatest potential to change the global atmospheric ice nucleation regime so far recognized in geological history. Other ejecta blankets of similar size cannot be said to have changed the ice nucleation regime from background. This distinction provides the basis for assessing the impacts of direct v. indirect effects of mineral aerosol. Comparison with other data from stratigraphic records presents an opportunity to further understanding of links within the Earth system via cloud behaviour and properties.
G. Osinski, L. Pittarello and H. J. Melosh are thanked for providing thoughtful reviews, which improved the paper. I. Fitzsimons is thanked for his editorial role
Scientific editing by Ian Fitzsimons