Dynamical effects of multiple impacts: Large impacts on a Mars-like planet

Thomas Ruedas; Doris Breuer

doi:10.1016/j.pepi.2019.01.003

Dynamical effects of multiple impacts: Large impacts on a Mars-like planet

Ruedas, Thomas;Breuer, Doris 2019-01-27 00:00:00 1,2,3 Thomas Ruedas Doris Breuer Dynamical eects of multiple impacts: Large impacts on a Mars-like planet nal version 9 January 2019 published in: Phys. Earth Planet. Inter. 287, pp. 76{92 (2019) Institute of Planetary Research, German Aerospace Center (DLR), Berlin, Germany Museum fur Naturkunde Berlin, Germany Institute of Planetology, Westf alische Wilhelms-Universit at, Munster, Germany The version of record is available at http://dx.doi.org/10.1016/j.pepi.2019.01.003. This author post-print version is shared under the Creative Commons Attribution Non-Commercial No Derivatives License (CC BY-NC-ND 4.0). arXiv:1901.09463v1 [astro-ph.EP] 27 Jan 2019 Dynamical eects of multiple impacts: Large impacts on a Mars-like planet Thomas Ruedas Institute of Planetary Research, German Aerospace Center (DLR), Berlin, Germany Museum fur Naturkunde Berlin, Germany Institute of Planetology, Westfalische Wilhelms-Universitat, Munster, Germany Doris Breuer Institute of Planetary Research, German Aerospace Center (DLR), Berlin, Germany Highlights • Combined eects are not simply linear superpositions of individual impacts. • Many details of impact history are not preserved in the long-term evolution. • Global signatures of dierent impact histories converge in the long term. • Lithospheric delamination may prolong melt generation, aect core. Abstract The earliest stage of the evolution of a fully assembled planet is profoundly aected by a number of basin-forming impacts large enough to change the dynamics of its deeper interior. These impacts are in some cases quite closely spaced and follow one another in short time intervals, so that their eects interact and result in behavior that may dier from a simple sum of the eects of two individual and isolated impacts. We use two-dimensional models of mantle convection in a Mars-like planet and a simple parameterized representation of the principal eects of impacts to study some of the dynamical eects and interactions of multiple large impacts. In models of only two impacts, we con rm that the dynamical eects of the impacts reinforce each other the closer they are in space and time but that the eects do not always correspond to straightforward superpositions of those of single, isolated impacts. In models with multiple (4{8) impacts with variable sizes, distances, and frequencies, the global response of the mantle is as variable as the impact sequences in the short term, but in the long term the dierent evolutionary paths converge for several indicator variables such as the mean ow velocity, temperature, or heat ow. Nonetheless, beyond a certain impact frequency and energy, lithospheric instabilities triggered by large impacts occur on a global scale, reinvigorate mantle dynamics for long time spans, and entail a late stage of melt production in addition to the initial melting stage that is not observed in one- or two-impact models. After one or several very large impacts, some lithospheric material may founder and sink to the core{mantle boundary, and if enough of it accumulates there, it enhances the heat ux out of the core for several hundred millions of years, with possible eects on dynamo activity. Mars; mantle convection; impact 1 Introduction In recent years, the interactions between very large meteorite impacts and planetary interiors have received increasing attention, and several authors have studied a variety of aspects of impact-induced mantle dynamics especially on Mars (e.g., Reese et al., 2002, 2004; Watters et al., 2009; Roberts et al., 2009; Arkani-Hamed and Ghods, 2011; Golabek et al., 2011; Roberts and Arkani-Hamed, 2012; Ruedas and Breuer, 2017, 2018), but also on Mercury (Roberts and Barnouin, 2012; Padovan et al., 2017), Venus (e.g., Gillmann et al., 2016), Earth (e.g., O'Neill et al., 2017), and the Moon (Rolf et al., 2017). Most of these studies concentrate on the eects of a single large impact in order to isolate the physical eects more clearly. However, in the earliest history of the terrestrial planets, basin-forming impacts have followed each other in close succession, and at times also in close proximity, which raises the question whether and how the eects of subsequent large impacts on the deep interior of the target planet may in uence each other. Roberts et al. (2009); Roberts and Arkani-Hamed (2012) and Rolf et al. (2017) have already explored this possibility to some extent for multiple impacts on Mars and the Moon, respectively, and observed that a succession of large impacts leaves a cumulative signature on physical variables such as the heat ow. Corresponding author: T. Ruedas, Institute of Planetary Research, German Aerospace Center (DLR), Berlin, Germany (Thomas.Ruedas@dlr.de) 2 In this study we consider dierent successions of impacts, both idealized ones and examples derived from observed martian giant impacts. Our principal goal is to elucidate some more details of how impacts in close temporal or spatial proximity may in uence each other and whether the cumulative eects of a succession of several impacts are additive or lead to some sort of saturation, e.g., because the aected mantle is so depleted in fusible mineral components that it ceases to produce melt readily. 2 Method and model set-up The evolution of the martian interior including concomitant processes such as melting is represented as a two- dimensional fully dynamical model using the convection code STAGYY (Tackley, 1996, 2008) with the additions and modi cations described in detail in Ruedas and Breuer (2017, 2018) on a two-dimensional spherical annulus grid (Hernlund and Tackley, 2008) with 512×128 points. The convection algorithm solves the equations of mass, mo- mentum, and energy conservation in the anelastic, compressible approximation, making use of a detailed petrological model that includes, among others, melting-related density changes and mineral solid-phase transitions. Following the preferred interior models by Konopliv et al. (2011) and Rivoldini et al. (2011), we de ne our models to consist essentially of two layers, i.e., an upper mantle with an olivine-dominated mineralogy and a lower mantle with an assemblage mostly of a high-pressure olivine polymorph and majorite garnet (Bertka and Fei, 1997); there is no basal (pv+fp) layer. The chemical model assumed for the mantle is the one by W anke and Dreibus (1994), and for the core, we use a sulfur content derived from Rivoldini et al. (2011). The viscosity of the mantle is described by an extended Arrhenius-type law for diusion creep of olivine that includes dependencies on temperature T , pressure p, and chemical composition, in particular water content (Hirth and Kohlstedt, 2003; Zhao et al., 2009); its value at surface pressure and the initial potential temperature of the mantle for the bulk water and iron contents of unmolten mantle material is also the scaling viscosity of the model. Internal heating is provided by the radioactive decay of 40 232 235 238 K, Th, U, and U, which are tracked with tracer particles and are redistributed by melting, by which they become concentrated in the melt and in the crust forming from it; the movement and partitioning of water is treated in the same manner. The surface and core{mantle boundary (CMB) are implemented as isothermal free-slip bound- aries, whereby the temperature at the CMB is allowed to evolve with time according to a simple energy balance and evolution model of the core; the core model largely follows the approach by Nimmo et al. (2004) and Williams and Nimmo (2004) and is described in some detail in Ruedas et al. (2013, App. B). The model is initialized on the basis of a radially symmetric adiabatic state prevailing at 4.5 Ga with a core superheated by 150 K and a surface at 218 K; the surface temperature is kept constant in all models over the entire evolution, as it is not thought to have changed strongly enough to have had a feedback eect on the interior. The actual initial state is derived from this con guration by setting the supersolidus region to the solidus, extracting the corresponding amount of melt to form an initial crust of 50 km, and redistributing the depletion of the mantle residue randomly throughout the mantle; the latter lets the mantle appear well-mixed and compositionally homogeneous on large scales. The thickness of this primordial crust is not well constrained, but with the extraction threshold and the mantle solidus xed, it is mostly dependent on temperature. Subsequent regular melting in the course of the evolution is controlled by the extent to which the temperature exceeds the local solidus; if the amount of newly formed melt passes a certain threshold value and has a supersolidus vertical path to the base of the lithosphere, it is assumed to be extractable, and the excess amount is instantaneously moved to the top to build the crust. The extraction threshold has been estimated to lie between 0.1% (e.g., McKenzie, 1985) and 1{2% (e.g., Faul, 2001), and so we chose a value of 0.7%. The solidus is the same as in Ruedas and Breuer (2017, App. A). The most important model parameters are summarized in Table 1. Further technical detail may be found in Ruedas and Breuer (2017, 2018) and their supplementary material. The impact process is not included via coupling or interfacing the convection algorithm with a fully dynamical hydrocode simulation, but is described using a combination of simpler analytical or semi-empirical approximations. As the impacts we consider are modeled after real impact basins on Mars, we start from the known diameter of the nal crater, D (cf. Table 2), and use scaling laws for complex craters to work back to determine the diameter of the impactor: 0:427 0:28 g % 1:1346 0:1475 D = 0:6918D D (1) imp f sc 0:56 v % imp z;imp (e.g., Melosh, 2011; Ruedas and Breuer, 2018). Here, v = v sin is the vertical component of the impactor z;imp imp velocity, and D is the diameter of the transition between simple and complex craters; in the absence of more speci c sc information, the angle is always set to 45°, which is the most probable impact angle for an isotropic impactor ux (e.g., Gault and Wedekind, 1978). The density and the absolute value of the impactor velocity are set to those of a homogeneous rocky asteroidal object hitting Mars (cf. Table 1; Carry, 2012; Ivanov, 2001). With regard to the energy input from the impact, we follow the approach outlined by Watters et al. (2009) in many respects, i.e., we base the estimate of the shock pressure on a one-dimensional impedance-match calculation and derive the waste heat from that value as a function of distance from the impact center using the \inverse-r" pressure decay formula from Ruedas (2017). The temperature in the shock-heated volume is forced to a value just above the solidus after the production and extraction of the impact melt, for which an upper limit is imposed. The melt extracted is instantaneously added to the top of the model. 3 Table 1: Main model parameters Mars Planetary radius, R 3389.5 km Archinal et al. (2011) Mass 6:41712 10 kg Konopliv et al. (2011); Jacobson (2010) Surface gravity, g 3:717 m/s Surface temperature 218 K Catling (2015) Mantle thickness, z 1659.5 km after Rivoldini et al. (2011); Konopliv et al. (2011) Initial potential temperature, T 1700 K pot Initial core superheating 150 K Surface porosity, ' 0.2 Cliord (1993) surf Melt extraction threshold 0.007 Simple/complex crater transition, D 5.6 km Robbins and Hynek (2012, Tab. 3) sc Bulk silicate Mars Mg# 0.75 W anke and Dreibus (1994); Taylor and McLennan (2009) Initial bulk silicate water content 36 ppm W anke and Dreibus (1994) Present-day K content 305 ppm W anke and Dreibus (1994) Present-day Th content 56 ppb W anke and Dreibus (1994) Present-day U content 16 ppb W anke and Dreibus (1994) Core S content 16 wt.% Rivoldini et al. (2011) Impactor Density, % 2720 kg/m Carry (2012) imp Velocity, v 9.6 km/s Ivanov (2001) imp In addition to a reference model without any impacts, this study comprises two model subsets with a total of 32 models. The rst subset consists of models with two impacts of identical size, of which the rst strikes at 4 Ga (i.e., at t = 500 Myr model time) and whose spatial and temporal separation is varied systematically in order to study possible eects. The impacts of the two series of this subset each have the magnitude of either an Isidis-forming (I) or of a Utopia-forming (U) event, respectively. The spatial separation x is measured in multiples of the isobaric core diameter D and is 1, 2, 5, or 10 in the models of either series; the isobaric core diameter is de ned based on ic the position of the in ection point of the \inverse-r" pressure decay curve from Ruedas (2017) (cf. Eq. 17). The temporal separation t is given in multiples of an estimated \decay time" t needed for certain dynamical eects caused by a single impact of a given magnitude to decay to levels close to the pre-impact state. We chose the mean ow velocity v as the main indicator for the dynamical state and consider the decay to be essentially complete when rms the dierence of v between the disturbed model and the impact-free reference drops below 5% of the maximum rms disturbance at the time of impact; the exact choice is somewhat arbitrary, and we rounded the time to full million years for convenience. In the model series, the second impact follows the rst at 0.5, 1, and 2 decay times. x and t are thus functions of the size of the impact, and as Isidis and Utopia have quite dierent magnitudes, so have D ic and t : for Isidis, we have D = 135:4 km and t = 6 Myr, for Utopia, we have D = 385:6 km and t = 12 Myr. In d ic d ic d addition to these two-impact models, each series also includes an additional model with a single impact identical to the rst impact as a further aid in identifying the particular eects of the second impact. The second subset of models includes six models with a succession of four to eight large impacts on Mars taken from the list of 20 large basins by Frey (2008) (cf. Roberts et al., 2009), whereby the impact sequences are chosen such that all events in a given model lie approximately on a great circle (GC) (Fig. 1). The spacing of the impacts corresponds to the approximate position of the basin centers on the great circle, and their timing follows the ages given by Frey (2008); the method of mapping them is outlined in A. This model set thus combines a certain element of randomness with the constraints from observations on Mars. 3 Results All models in this study are derived from an impact-free reference model whose setup is identical with the 36 ppm H O reference models of Ruedas and Breuer (2017, 2018) except that it starts at 4.5 Ga instead of 4.4 Ga. The reference model starts with a regular pattern of upwellings and downwellings that transitions into a more chaotic mode after a few tens of millions of years. After 500{600 million years, the ow reorganizes itself into a more regular large-scale pattern again over the course of several hundreds of millions of years, ending up with a long-lived nal con guration with ve approximately regularly spaced large plumes while cooling steadily and forming a thick thermal boundary layer at the surface. The initially superheated core cools quite rapidly and loses its excess heat in the rst few hundreds of millions of years, but it never begins to crystallize. The divergence of the evolutionary paths of the models thus begins with the rst impact striking in a given model. For the two-impact models, this occurs early in the transition from the chaotic to the nal quasi-static con guration, for the great-circle models it starts at some point during the chaotic stage. Hence, all impacts, and especially the earliest ones, strike a planet with a still thin lithosphere and hot, globally vigorously convecting mantle whose uppermost part is undergoing some regular global melting, albeit at quite low intensity. 4 5 Table 2: Impact basins considered in this study in order of descending age, after Frey (2008) and Roberts et al. (2009). The impactor diameter D and the diameter D imp ic and depth z of the isobaric core are estimates from scaling laws derived from the observed nal basin diameter D . The ages and positions in the list apply only to impacts ic f in the GC models; the position given in the \Position" column are the geographical coordinates, the single degree entries in the GC columns are the positions on the great circle of the GC model determined as described in A, measured clockwise from the bottom of the circle in the cross-section images. Basin/Code Age D D z D Position Model sets imp ic ic f (Myr) (km) (km) (km) (km) GC1 GC2 GC3 GC4 GC5 GC6 U/I-* Amenthes (Ame) 4218 188.3 103.8 51.8 1070 110.6°E/0.9°S 73.0° 271.5° Zephyria (Z) 4210 213.2 117.4 58.6 1193 164.3°E/12.4°S 255.9° 126.9° 257.6° Daedalia (Dae) 4199 526.9 290.4 144.0 2639 228.3°E/26.5°S 188.1° 243.4° Sirenum (Sir) 4196 188.0 103.6 51.8 1069 205.3°E/67.4°S 200.0° SW Daedalia (SWD) 4176 231.9 127.8 63.3 1278 213.9°E/29.4°S 175.3° 235.6° Ares (Are) 4167 680.5 375.0 185.4 3300 343.9°E/4.0°N 84.0° 303.8° 83.5° Amazonis (Ama) 4154 580.9 320.2 158.5 2873 187.9°E/27.1°N 299.5° 350.7° 302.1° In Amazonis (IA) 4152 207.2 114.2 56.5 1156 192.5°E/29.3°N 302.9° 354.8° 306.6° Solis (Sol) 4148 311.1 171.4 85.4 1663 275.3°E/23.8°S 117.5° 230.4° N Tharsis (NTh) 4143 244.8 134.8 67.3 1347 243.6°E/17.6°N 66.0° 290.0° 42.6° Chryse (C) 4140 324.3 178.6 89.0 1725 318.0°E/25.0°N 51.0° Hematite (Hem) 4138 189.1 104.2 51.5 1065 357.8°E/3.2°N 87.6° 93.5° Scopolus (Sc) 4133 439.5 242.2 120.2 2250 81.8°E/6.9°N 276.4° 42.9° Acidalia (Aci) 4132 629.6 347.1 172.0 3087 342.7°E/59.8°N 29.7° 9.6° North Polar (NP) 4124 297.7 164.0 81.8 1600 195.2°E/80.0°N 350.6° Utopia (U) 4111 699.5 385.4 190.5 3380 115.5°E/45.0°N 324.3° (X) SE Elysium (SEE) 4107 256.2 141.1 70.5 1403 170.3°E/3.7°N 272.9° 273.5° Hellas (Hel) 4065 399.5 220.1 109.4 2070 66.4°E/42.3°S 226.6° 133.7° 170.6° Argyre (Arg) 4043 238.2 131.2 65.5 1315 317.5°E/49.0°S 158.6° Isidis (I) 3810 247.5 136.4 67.5 1352 87.8°E/13.4°N 284.9° 46.8° (X) North Polar 60˚ 60˚ Acidalia Utopia In Amazonis Chryse N Tharsis Isidis Amazonis Hematite SE Elysium Scopolus 0˚ 0˚ Ares Amenthes Zephyria Solis SW Daedalia Daedalia Hellas 2 4 Argyre Sirenum −60˚ −60˚ Figure 1: Schematic map in Hammer projection of the locations of the large impact basins on Mars after Frey (2008) listed in Table 2 and the great circles (colored curves) by which they are assigned to 2D models. 3.1 Two-impact models The quasi-instantaneous input of energy into the planetary interior by an impact produces a steep change in the temporal evolution of several variables that serve as indicators of the dynamical state of the mantle, e.g., the global means of the root-mean-square velocity, v , the mantle temperature, T , and the heat ows through the surface rms mean and the CMB, q , and q . This step-like change is followed by an initially rapid rebound towards the pre-impact t CMB state over the next few millions of years (Fig. 2). The signal of the second impact is added to the decaying signature of the rst in dierent ways for dierent variables. To compare the second impact with the rst, we determine the maximum change of these variables due to each impact and calculate the ratio of the second signal to the rst; e.g., for T the second-to- rst ratio of impact eect amplitudes would be mean T (t ) T (t ) mean 2 mean;sngl 2 = ; (2) T (t ) T (t ) mean 1 mean;ref 1 where \sngl" and \ref " indicate the corresponding models with a single impact and the impact-free reference model, respectively, and t and t indicate the times of the rst and the second impact. These ratios are expressed as functions 1 2 0 0 of the normalized distance, x = x=D , and time lapse, t = t=t , respectively. If this ratio is smaller than 1, ic d i.e., if the change of the indicator variable is smaller in the second impact than in the rst, that would suggest that the eect is not simply additive but may tend towards a state of saturation. If the ratio were equal to 1, impacts with the same parameters would have identical eects, implying that the target has no memory of the eects of preceding impacts. is therefore a measure of the superposition eciency of impacts with respect to the indicator variable X . The models show that the superposition eciency for v (t) tends to increase with decreasing x and also, rms although much less clearly, with t (Fig. 3), because the lingering thermal anomaly from the rst impact boosts the upwelling triggered by the second. The physical cause for this boost may be that in closely spaced anomalies whose temperatures decrease away from the central axis, the superposition leads to particularly high peak temperatures near the axis, which have a non-linear positive eect on the viscosity and thus on ow velocities. A simpli ed analytical model for pipe-like ow with temperature-dependent viscosity indicates that a strong focusing of the ow and a speedup by more than an order of magnitude relative to the isoviscous case is possible for peak temperature anomalies 0 0 of a few hundred kelvins (cf. B). (x ; t ) is not linear in both variables, and the data suggest that there may be a slight increase beyond a certain x in the Utopia-size models, which may be due to the increase of the aected total volume. Except for closely spaced impact pairs, < 1. Later peaks in v (t) in Fig. 2 are not directly related v rms to the impacts and re ect other mantle dynamical processes such as rising plumes or delamination of crust. They do not appear in the impact-free reference model and are an indirect, longer-term consequence of local lithospheric instabilities induced by the impacts. To assess the thermal eect, we consider rst a simpli ed setup consisting of two identical spherical anomalies of radius R whose centers are separated by x; for x < 2R, these anomalies would overlap to some extent. Especially for impacts in short succession, the thermal anomaly from the rst strike is still so strong that the joint anomaly after the second one can be treated as a homogeneous volume as an approximation. The magnitude of the thermal anomaly formed by overlapping spheres is proportional to its volume V (x), and the functional form of V points to a sublinear increase with x up to x 2R (cf. C); for larger osets, it should be constant, as both anomalies are independent 6 10 Δt=0.5t Δt=2t d d Reference Single Isidis Double Isidis, Δx=D ic Double Isidis, Δx=10D ic Single Utopia Double Utopia, Δx=D ic Double Utopia, Δx=10D ic 100 500 550 600 650 500 550 600 650 t t (Myr) (Myr) Figure 2: Root-mean-square velocity, mean temperature, and surface heat ow for selected models with two impacts of either Isidis-like or Utopia-like magnitude. The impact-free and the single-impact models are also shown for reference. The Utopia models have been smoothed at t > 530{550 My by averaging over 2 My-intervals for better readability. The ticmarks on the top abscissas mark the times at which impacts occurred in any of the models. A collection of plots for all models can be found in the Supplementary Material. T (K) 2 mean q ) (mW/m v (cm/yr) rms 1.2 Δt'=0.5, Δx'=1, I I 1.15 Δt'=0.5, Δx'=1, U U Δt'=1, Δx'=2, I I Δt'=1, Δx'=2, U U 1.1 Δt'=2, Δx'=5, I I Δt'=2, Δx'=5, U U Δx'=10, I 1.05 Δx'=10, U 0.95 0.9 1.1 0.9 0.8 0.7 0.6 0.5 1.1 0.9 0.8 0.7 0.6 0.5 0.4 0 2 4 6 8 10 0 0.5 1 1.5 2 Δx' Δt' Figure 3: Second-to- rst ratio of impact eect amplitudes (superposition eciencies) of v , T , and q as functions rms mean t 0 0 of x (left column) and t (right column). α α q,t T from each other. Estimating the cooling time of a spherical anomaly from the analytical solution for diusive cooling 2 2 of a solid sphere, one nds that the temporal decay is governed by a term of the form exp( t=R ), for a diusivity (e.g., Yeom, 2015). The level from which the mantle is reheated by the second impact depends on how much the rst anomaly has cooled by the time it strikes, and so we expect a decrease of the thermal superposition eciency that approximately reproduces that trend, whereby the eective diusivity is enhanced by convection within the anomaly and entrainment of cooler material from outside. The trends in Fig. 3 con rm that grows approximately linearly with spatial separation and decreases approximately linearly with temporal separation; the apparently linear form of the latter is due to the limited range of t considered, whose lowest value of 0.5 lies already after the stage of strongest cooling. The increase with x is clearer for the Utopia-size impacts, presumably because they penetrate more deeply into the sublithospheric mantle and are hence less aected by the cold lid. As < 1 in almost all of the two-impact models, the mantles move towards a state of saturation in which further impacts ideally would not entail more heating but transform their energy otherwise, e.g., by melt production. The signature of an impact on the global surface heat ow is proportional to the surface area aected, whereby the strongest eects are expected within the nal crater. Under the assumptions made in our models, we expect that the signatures of multiple impacts are only additive if the distance between two impacts exceeds a certain value, because otherwise the aected areas will overlap to some extent. Impact modi cation of a surface region that has already been modi ed by a previous impact is reduced, because processes like pore lling of regolith by volcanism can only occur once. In general, the superposition eciency decreases approximately linearly with increasing spatial and q;t temporal separation (see Fig. 3 and D); however, for the large Utopia impacts, there is a clear minimum at x = 5 at all t , which indicates that other processes also leave their mark on q (t). The CMB heat ow, by contrast, shows only minute variations as a direct consequence of the impact in our models, and no meaningful dependence of the 0 0 peak amplitudes on x or t is found (cf. Fig. 6, left column). The shock-heated volumes are also depleted in fusible components and are thus compositionally distinct and less dense than pristine mantle material, which contributes to their buoyancy and the reinforcement of convection caused by impacts. On the other hand, a mantle volume that was strongly depleted by extensive melting in the rst impact cannot experience much further depletion by melting from a second nearby impact. Thus, the buoyancy of the ow after the latter should be expected to be more thermally dominated in models with small x , but the potentially diminished compositional contribution to the total buoyancy is apparently not signi cant enough to compensate the other reinforcing eects on v discussed above. As the impact-generated anomalies ascend and spread out beneath rms the lid, they may not only in uence each other dynamically, as seen in the v (t) curves, but will also collide and rms merge, especially if they are closely spaced. In models with widely spaced impacts in close succession, both anomalies develop somewhat independently for a certain time until they have spread far enough below the lid to run into each other. In such cases a piece of normal mantle can get caught between them and induce a downwelling due to its relatively higher density, especially for smaller impacts with less vigorous dynamics. An example for this can be seen 0 0 in the double-Utopia model with x = 10 and t = 0:5 in the lower right block of Fig. 4. Heat and enhanced convection can also result in increased production of melt and crust, but the extreme depletion of the shallower mantle in the impact-aected region counteracts this eect. Global crustal thickness evolution maps h(x; t) and the corresponding cumulative distributions H (h; t) (Figs. 5, S3{6) show that impacts broaden the range of thickness by creating thinner and thicker regions, especially in the large (Utopia-sized) impact model series. In the following, we will normalize h with the mean crustal thickness h at the time t of the rst impact, and we express the 0 0 time elapsed since t in multiples of the decay time t in order to achieve a scale-free representation that might reveal 0 d common trends for impacts of dierent magnitudes. The cumulative distribution function of the crustal thickness H (h; t) quanti es how large the fraction of crust with a thickness no larger than h is at the time t; for example, in the scale-free plots of some selected cases in Fig. 5, H (1:1; t + 0:5t ) = 0:87 means that at half a decay time after 0 d the rst impact, 87% of the entire crust has a thickness of 1:1h or less. The curves in the left column of the gure show H (h) at that time after the second impact at which a part of the crust reaches near-maximum thickness and has a diverse thickness distribution. As the curve for the reference model shows, the thickness distribution in general is essentially unimodal and shows little scatter around the mean value h=h = 1, i.e., the crust has a rather uniform thickness. The unimodality remains largely valid in the models with impacts as well, but the steep increase at 1 is not quite as high, whereas the anks of the curve become longer and higher, re ecting the broadening of the thickness distribution due to production of anomalous crust by impacts. The color-coded maps in the middle and right column of Fig. 5 show how this cumulative distribution evolves through time in a time window around (mostly after) the impact time; the curves in the left column are thus sections of these maps in the vertical direction and are marked as vertical lines in the map plots. The spatial pattern and evolution of post-impact crustal thickness depends on spacing, timing, and size of the impacts. Comparison of the impact-free and the single-impact model shows clearly that the additional melt extraction after the impact results in a marked and lasting increase in h in some regions, especially near the site of a Utopia-sized impact (cf. Fig. S3). In the Utopia series, the second impact then increases the fraction of thickened crust, possibly including some more growth of already thickened regions, but at least in models in which the impacts are further apart in space and time, the most extreme eect is not as long-lived. On the other hand, some parts of the crust become thinner a certain time after the impacts (usually 5t to 15t after t ), which seems to d d 0 be related to delamination of very thick crustal roots, as the increase in the fraction of thinner-than-average crust is related to the gradual disappearance of some of the thickest crust in the H diagram. Note that the rather discrete, step-like thickness increments are partly an artifact of the grid discretization, which limits the accuracy with which 9 double-Isidis double-Utopia 2 2 2 2 2 1 1 1 1 1 2 2 2 2 1 1 1 1 Figure 4: Snapshots of temperature (T , left column) and depletion (f , right column) for the models with the closest and widest impact spacing and the shortest time interval. The impact locations are marked in the upper rows of each block as 1 and 2. The rst impact occurs at 500 My in all models, the second impact follows at 503 My in the double-Isidis models and at 506 My in the double-Utopia models, respectively (upper rows of each block). The evolution immediately following the second impact is shown in the snapshots in the lower rows of each block. 10 the crustal thickness can be determined. In the Isidis series, these features are less clear, and the thinning is much less signi cant as less delamination takes place. The occurrence of extra-thick crust due to the second impact, however, is also visible but has a much more transient character and fairly uniform lifetime of about 10t in the entire parameter range covered. As a consequence of the interaction of the two anomalies, the location of the thickest post-impact crust does not necessarily coincide exactly with the impact site in both series. In all cases investigated, the dierences between the models diminish with time, and the impact signature in the temporal evolution of the system's dynamical variables fades and has disappeared long before the present. In particular, thermal anomalies fade with time and leave no signal in the present-day heat ux. Compositional anomalies, however, are preserved, as pointed out by Ruedas and Breuer (2017), but it would be dicult to draw a sharp boundary between the region of in uence of one impact or the other if the impacts are close enough to overlap. Variations in post-impact crust formation are also re ected in the crustal thickness and may thus preserve a long-term record of impact-induced mantle dynamics. In addition to the eects of impacts at the surface or in the shallow mantle, it has been proposed that impacts also aect a core dynamo that may be active (e.g., Roberts and Arkani-Hamed, 2014). We did not include any direct impact eects such as shock heating of the core in our treatment as did those authors and therefore neglect the potential eects of an impact-generated hot layer at the top of the core whose appearance was proposed by them. However, given that even the largest impacts in our models are somewhat smaller than the one modeled by Roberts and Arkani-Hamed (2014) and that shock heating decreases to nearly zero towards the CMB in our models, we do not expect a signi cant direct eect of the shock on the core. Impact-related variations of the heat ow through the CMB would therefore mostly be induced by post-impact mantle processes. The left panel of Fig. 6 shows the total entropy production of the core, which is a good proxy for q (t) in our models, for the rst billion years in the double-Utopia CMB 0 0 models with the longest time interval between impacts (t = 2); the curves for smaller t look very similar, whereas in the double-Isidis models the deviations from the reference are negligible. At the times of the impacts, there are tiny downward excursions that indicate a short period of reduced entropy production and heat ow out of the core during the presence of the thermal anomaly in the mantle; this eect would be compensated if shock heating also aected the core. The more signi cant eect, however, sets in only several decay times after the second impact and consists in an increased entropy production and heat ow out of the core. It is caused by the arrival at the CMB of cooler lithospheric material into which small fragments of eclogitized crust are embedded and that had become unstable in the aftermath of the impacts, delaminated, and sunk to the bottom of the mantle where it cools the core from above. However, in the lower part of the martian mantle, eclogite has a signi cantly smaller density excess than in the shallow mantle (e.g., Aoki and Takahashi, 2004), and as the crustal fragments form a close union with the lighter harzburgitic residue that went down along with them, the old lithosphere does not accumulate permanently at the CMB. Instead, it is heated enough by its internal sources in the eclogite and by conduction from the surrounding warm mantle to become unstable and be mixed back into the mantle again by the general mantle ow. Hence, while spacing and time interval between impacts are unimportant for post-impact eects on the CMB, the size and, to some degree, the number of the impacts are decisive, because they control the extent of lithospheric instability. 3.2 Great-circle models As the models with multiple impacts on a great circle show, a succession of various dierent large impacts produces a strongly variable depletion pattern in the uppermost mantle (Fig. 7). This pattern still re ects the diversity of the impacts that produce it, but the vigorous post-impact dynamics and merging of the individual anomalies precludes the distinction of clear boundaries between the traces of the discrete events. As the individual impacts in these models lie further apart from each other and span a longer time interval than those in the previous sets, the more indirect eects as a whole, especially those on the lithosphere, also last longer and extend to the entire planet: speci cally, crustal delamination as a long-term consequence of the impacts becomes a planet-wide process and lasts much longer than in the previous sets. Figure 8 shows the root-mean-square convection ow velocity, the mantle mean temperature, and the surface heat ow for models GC3, GC4, and GC5, which represent the models with the highest, an intermediate, and the lowest number of impacts on a great circle in the model set. The basic shape of the curves resembles those from the idealized double-impact models in the previous section. The impacts are again marked by sharp peaks with subsequent gradual decays of these variables; later narrow peaks in v correspond to the quick motion of cold lithospheric instabilities rms formed in the aftermath of the impacts, mostly related to the formation of thick crust. Although we do not normalize times here with the decay time t because there is no single t applicable to the entire duration of a given model, it d d is worth noting that estimates for the t for individual impacts seem to follow the aforementioned relation between t d d 2 2 and the size of the thermal anomaly, exp( t=R ). The center panel of Figure 8 shows that the T of models mean with particularly intense impact activity even drops below the T of the impact-free reference model some time mean after the impacts. The reason for this maybe unexpected behavior is that the heat from the impact shocks has largely diused away after some tens of millions of years, but the copious amounts of cool crust that were produced in their aftermath are becoming unstable and partly sink into the mantle, where they form ecient heat sinks; the crust is so cold, because the melt from which it is formed has been set to the surface temperature upon eruption. At some point between 1.5 Gy and 2.5 Gy, i.e., outside the plotted range, all of the models return to and remain at T values above the reference as the delaminated crust has been warmed up by its own internal radioactive heat mean 11 1 2.2 2.2 No No impact impact reference reference 2 2 0.8 1.8 1.8 0.6 1.6 1.6 1.4 1.4 0.4 1.2 1.2 0.2 1 1 Reference model @ t 0 0.8 0.8 0.8 1 1.2 1.4 1.6 1.8 2 2.2 0 5 10 15 20 25 30 0 5 10 15 20 25 30 h/h (t-t )/t (t-t )/t 0 d 0 d 1 2.2 2.2 Isidis, Utopia, single single 2 2 0.8 1.8 1.8 0.6 1.6 1.6 0.4 1.4 1.4 1.2 1.2 0.2 Single +6t Isidis @ t 0 d 1 1 Single Utopia @ t 0.8 1 1.2 1.4 1.6 1.8 2 2.2 0.8 0.8 0 5 10 15 20 25 30 0 5 10 15 20 25 30 h/h (t-t )/t (t-t )/t 0 d 0 d 2.2 2.2 Double Double Isidis, Utopia, Δx'=1, Δx'=1, Δt'=0.5 Δt'=0.5 2 2 0.8 1.8 1.8 0.6 1.6 1.6 0.4 1.4 1.4 1.2 1.2 0.2 Double +1.5t Isidis, Δx'=1, Δt'=0.5 @ t 0 d 1 1 Double +0.5t Utopia, Δx'=1, Δt'=0.5 @ t 0 d 0.8 0.8 0.8 1 1.2 1.4 1.6 1.8 2 2.2 0 5 10 15 20 25 30 0 5 10 15 20 25 30 h/h (t-t )/t (t-t )/t 0 d 0 d 1 2.2 2.2 Double Double Isidis, Utopia, Δx'=1, Δx'=1, Δt'=2 Δt'=2 2 2 0.8 1.8 1.8 0.6 1.6 1.6 1.4 1.4 0.4 1.2 1.2 0.2 Double +2t Isidis, 1 1 Δx'=1, Δt'=2 @ t 0 d Double +2.5t Utopia, Δx'=1, Δt'=2 @ t 0 d 0.8 0.8 0.8 1 1.2 1.4 1.6 1.8 2 2.2 0 5 10 15 20 25 30 0 5 10 15 20 25 30 h/h (t-t )/t (t-t )/t 0 d 0 d 2.2 2.2 Double Double Isidis, Utopia, Δx'=10, Δx'=10, Δt'=0.5 Δt'=0.5 2 2 0.8 1.8 1.8 0.6 1.6 1.6 0.4 1.4 1.4 1.2 1.2 0.2 Double +6t Isidis, Δx'=10, Δt'=0.5 @ t 0 d 1 1 Double +0.5t Utopia, Δx'=10, Δt'=0.5 @ t 0 d 0.8 0.8 0.8 1 1.2 1.4 1.6 1.8 2 2.2 0 5 10 15 20 25 30 0 5 10 15 20 25 30 h/h (t-t )/t (t-t )/t 0 d 0 d 2.2 2.2 Double Double Isidis, Utopia, Δx'=10, Δx'=10, Δt'=2 Δt'=2 2 2 0.8 1.8 1.8 0.6 1.6 1.6 1.4 1.4 0.4 1.2 1.2 0.2 1 1 Double +6t Isidis, Δx'=10, Δt'=2 @ t 0 d Double +2t Utopia, Δx'=10, Δt'=2 @ t 0 d 0.8 0.8 0.8 1 1.2 1.4 1.6 1.8 2 2.2 0 5 10 15 20 25 30 0 5 10 15 20 25 30 h/h 0 (t-t )/t (t-t )/t 0 d 0 d 0 0 0 0 0 0 0 0 0 0 0 0 0.2 0.2 0.2 0.2 0.2 0.2 0.2 0.2 0.2 0.2 0.2 0.2 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.6 0.6 0.6 0.6 0.6 0.6 0.6 0.6 0.6 0.6 0.6 0.6 0.8 0.8 0.8 0.8 0.8 0.8 0.8 0.8 0.8 0.8 0.8 0.8 1 1 1 1 1 1 1 1 1 1 1 1 H H H H H H H H H H H H Figure 5: Cumulative distribution function H of the crustal thickness of the two-impact model series from Fig. 2 as a function of crustal thickness h and time t at and after the impacts. The left column shows H (h) at that time after the second impact at which the crust reaches its maximum thickness. These curves are thus sections of the H (h; t) maps in the other two columns parallel to the ordinate; they are marked there as vertical lines in the corresponding color. The crustal thickness is normalized with the mean crustal thickness h immediately before the rst impact, i.e., at t = 500 Myr, which can be calculated from the impact-free reference model. The abscissa in the middle and right column gives the time since the rst impact in multiples of the decay timescale t , which is 6 Myr for the Isidis-sized impacts (left column) and 12 Myr for the Utopia-sized impacts (right column). The impact-free and the single-impact models are also shown for comparison. The full set of gures is provided in the Supplementary Material. H H H H H H h/h h/h h/h h/h h/h h/h 0 0 0 0 0 0 h/h h/h h/h h/h h/h h/h 0 0 0 0 0 0 35 Reference Reference Single GC1 double Δx=D Great-circle GC2 Utopia, models Δt=2t d ic Δx=2D GC3 ic Δx=5D GC4 ic Δx=10D GC5 ic GC6 -5 0 200 400 600 800 1000 0 200 400 600 800 1000 t t (Myr) (Myr) Figure 6: Total entropy production of the core for the double-Utopia models with the longest time interval between impacts (left) and for the great-circle models (right); a dynamo can exist if S > 0. The corresponding heat ows at the CMB can be found in the Supplementary Material. The ticmarks on the top abscissa mark the times at which impacts occurred in any of the models. SWD Dae Sir Hel Z Sol SWD Dae GC3 GC4 GC5 Ame Ame NTh Are Sc NTh Ama Aci IA Figure 7: Snapshots of the present-day depletion patterns for the great-circle models 3, 4, and 5. The northernmost points of the respective great circles lie at the six o'clock position. Snapshots of all models are presented in the Supplementary Material. S (MW/K) c production. In all curves, the largest impacts dominate strongly the total deviation from the impact-free reference. In the Introduction, we had speculated that at some point, a planetary mantle could reach a \saturated" state in the sense that additional energy from another impact does not reinvigorate the dynamics, for instance because the shallow mantle is so hot and/or depleted throughout that some sort of stable layering is established that suppresses further whole-mantle convection. However, in none of the models do we nd any clear indication of such a \saturation", which suggests that Mars was always far enough away from it, if such a state is possible at all. Nonetheless, and even though the larger spatial and temporal separation of the individual events diminishes the direct interactions between impacts, the dierent histories of the GC models induce dierent dynamical states of the mantle that lead to dierent signatures between models for a given impact. This becomes most readily visible in the two late impacts, Hellas (GC1, GC4, GC6) and Isidis (GC1, GC3). Both events leave signals of dierent magnitude in the dierent models in which they occur as a consequence of the dierent evolutionary paths (Fig. S7). In all six models, however, the global characteristics return to a state close to that of the reference on timescales that are substantially shorter than the age of the planet. Given the relatively short lifetime of those temperature-related anomalies, it is of interest to consider how longer- lived anomalies of other characteristics compare between great-circle models. Hence, we look again at the statistics of crustal thickness as expressed by H (h; t) (Fig. 9). As in the two-impact models, we observe that thick crust is formed immediately after the impacts but is partly destroyed by delamination on timescales of a few ten to a few hundred million years. However, the more perturbed dynamics of models with several impacts result in more extensive lithospheric instability in the later evolution, which in turn stirs the mantle up and enables prolonged melt and crust production in some places by transporting relatively fresh or re-enriched material into the melting zone; the refertilization of depleted mantle would be a consequence of the entrainment of delaminated crust (e.g., Rosenthal et al., 2014). As a consequence, there is a second era of (limited and localized) crustal production that begins considerable time after the end of the era of the large impacts and generates some sites of stable thick crust that persist to the present, which can be seen in the h(x; t) plots in the left part of Fig. 9. The resurgent crust production is also visible in some of the H (h; t) plots as a new increase in the maximum thicknesses. For the multi-impact sequences of the great-circle models, it is also worthwhile to take a look at q (t), because CMB in most of them we observe substantial lithospheric instability and subsequent transient accumulation of relatively cool material at the CMB. The right panel of Fig. 6 shows that indeed in most of these models core entropy production is increased substantially for several hundred million years after the impacts relative to the impact-free reference model, and q doubles in this time interval. The only exception, in which the increase is much less pronounced, is model CMB GC5, which had only four impacts and was also the least productive one in terms of post-impact crust formation. 4 Discussion Although the very setup of our models is arti cial by design in the two-impact models or due to the restriction in impact count and geometry imposed by the choice of the great circles and the two-dimensional grid, some lessons can be learned with regard to the geological implications of the martian impact record and its uncertainties. There are dierent compilations of large impact basins in the literature, and they show a substantial variance with regard to the timing, the size, and even the very number of impacts, because many of the old basins that predate Utopia are eroded or otherwise obliterated and therefore dicult to identify and characterize. At the one extreme, the survey by Frey and Mannoia (2013) identi ed some 30 large basins falling in the time interval between the formation of the dichotomy and the creation of Utopia, some of them of a magnitude comparable to the Utopia event. At the other extreme lies a recent analysis by Bottke and Andrews-Hanna (2017) that strongly indicates that such a large ux of giant impactors would almost certainly have hit some part of the dichotomy boundary, but as clear evidence for this seems to be lacking, it concludes that the actual number of such events should be less than 12. We decided to use the list by Frey (2008), because of its intermediate number of basin-forming impacts with potential eects on mantle dynamics and for easier comparison with the study by Roberts et al. (2009). Moreover, Werner (2008) identi ed partly dierent basins, assigned them ages that dier by tens or even more than hundred million years from those given in other compilations, and also determined dierent diameters. On the other hand, the global characteristics we chose to assess the dynamical state tend to converge towards the impact-free reference in the long run; for instance, the nal T in the great-circle models is only about 11{17 K higher and q less than 1.3 mW/m lower than in the impact-free mean t reference model. This suggests that the uncertainties in the dierent compilations are not a very strong concern with regard to the overall dynamical and thermal evolution of Mars. In the short-term aftermath of the impacts, however, there are substantial dierences, which would leave traces in the geological record that would still be visible today. With regard to the long-term preservation of the impact record, the triggering of lithospheric instability is of particular importance. A crucial factor in this process is the formation of crust that grows thick enough for its base to cross the basalt{eclogite transition and become gravitationally unstable. With our parameters, this transition would generally take place somewhere around 140 km depth in Mars, which is a bit less than twice the global mean thickness reached after a few hundred million years in the impact-free reference model. That global mean value is rather on the high side of other, independent crustal thickness estimates for Mars and would be lower if a lower initial mantle temperature or a higher melt extraction threshold had been chosen. A lower crustal thickness in our models would imply that lithospheric instability becomes less likely and less extensive or does not occur at all. Likewise, the cooling 14 7 Reference GC3 GC4 GC5 1750 0 300 350 400 450 500 550 600 650 700 t (Myr) Figure 8: Root-mean-square velocity, mean temperature, and surface heat ow for the great-circle models 3, 4, and 5, which include 8, 6, and 4 impacts of varying magnitudes, respectively. The impact-free model is also shown for reference. A collection of plots for all models can be found in the Supplementary Material. T (K) mean q ) (mW/m v (cm/yr) rms Ame Dae Sir SWD Are IA Ama Sol NTh Sc Aci Hel I No impact reference 360 2.2 1.8 1.6 1.4 1.2 0 0.8 0 500 1000 1500 2000 2500 3000 0 500 1000 1500 2000 2500 3000 t-t t-t (Myr) (Myr) GC3 360 2.2 Are 1.8 Sol 1.6 Dae SWD 1.4 Z 120 1.2 Ame 60 Sc I 0 0.8 0 500 1000 1500 2000 2500 3000 0 500 1000 1500 2000 2500 3000 t-t t-t (Myr) (Myr) 0 0 GC4 360 2.2 NTh 1.8 Dae SWD 1.6 Sir 1.4 Hel 120 1.2 Aci 0 0.8 0 500 1000 1500 2000 2500 3000 0 500 1000 1500 2000 2500 3000 t-t t-t (Myr) (Myr) 0 0 GC5 360 2.2 IA Ama Ame 1.8 1.6 1.4 1.2 NTh 0 0.8 0 500 1000 1500 2000 2500 3000 0 500 1000 1500 2000 2500 3000 t-t t-t (Myr) (Myr) 0 0 0.8 0.8 0.8 0.8 1 1 1 1 1.2 1.2 1.2 1.2 1.4 1.4 1.4 1.4 1.6 1.6 1.6 1.6 1.8 1.8 1.8 1.8 2 2 2 2 2.2 2.2 2.2 2.2 0 0 0 0 0.2 0.2 0.2 0.2 0.4 0.4 0.4 0.4 0.6 0.6 0.6 0.6 0.8 0.8 0.8 0.8 1 1 1 1 h/h h/h h/h h/h H H H H Figure 9: Normalized crustal thickness h(x; t)=h (left column) and cumulative distribution function H (h; t) of the crustal thickness (right column) of the impact-free reference model and the great-circle models 3, 4, and 5 (cf. Fig. 8 at and after the impacts. The crustal thickness is normalized with the mean crustal thickness h immediately before the rst impact in Table 2, i.e., Amenthes at 280 Myr, which can be calculated from the impact-free reference model. The abscissa gives the time since the rst impact; contrary to Fig. 5, no normalization has been applied, because the impacts have dierent sizes. The full set of gures is provided in the Supplementary Material. x x x x (°) (°) (°) (°) h/h h/h h/h h/h 0 0 0 0 Ame Ame Z Sir Dae SWD Dae IA Ama SWD Are NTh Sol NTh Aci Sc Hel I eect of the delaminated crust would be absent, and models with impacts would always be hotter at all times than the impact-free reference. Thus, post-impact dynamics also depends crucially on pre-impact crust formation and especially on the thickness of the primordial crust, which is poorly known. This highlights once more the importance of obtaining tight constraints on the thickness of the crust, e.g., from the seismic experiment on the InSight mission. Moreover, at the impact sites themselves the actual extent of extraction of melt generated directly by the impacts determines to a large degree the thickness of the newly formed post-impact crust. Our assumption of very ecient extraction is likely an extreme case and de nes an upper limit more conducive to crustal delamination. By contrast, Padovan et al. (2017) assumed no extraction of impact melt and formed their post-impact crust exclusively from melt produced in the aftermath of the impact by the dynamics of the evolving thermal anomaly, which results in a substantially thinner post-impact crust. The true extent of melt extraction in impacts is currently unknown and requires further study. A straightforward criterion for the occurrence of a second melting era, e.g., in terms of a threshold number of impacts or a minimum size does not seem to exist, even though we nd that the two-impact models and the great- circle model with the fewest impacts (GC5 with only four events) have no such era whereas the other GC models do. The absence of marked thermal anomalies at the time the second crust-forming era begins suggests that its cause is not primarily thermal; indeed, it is only hundreds of millions of years later that the average mantle temperature of the GC models is again higher than that of the impact-free reference model. The critical ingredient is therefore presumably the availability of easily fusible material such as delaminated crust, which however would rst have to be produced from very thick crust by extensive melting. This can happen more easily if the impacts shock-heat an already hot target, which requires them to penetrate more or less the entire lithosphere. Fig. 10 shows the depth reach of the impacts in some of the great-circle models as a function of time and in relation to the depth where 90% of the solidus temperature of fertile martian peridotite is reached; this value serves as a proxy for the shallowest depth at which sublithospheric conditions are reached or at which the mantle is so hot that it will melt without a lot of further energy input. The gure reveals that of the scanty series of impacts in GC5, only one (Amazonis) reaches rmly into the sublithospheric mantle, whereas the other series aect the sublithospheric mantle more strongly and over a longer timespan. On the other hand, the strong eect of the double-Utopia models with two deep-reaching impacts in rather short succession shows that the magnitude alone is not sucient. Instead, a longer gap between the individual events is more likely to ensure that the convective activity is maintained over a longer period and can stir and mix delaminated crustal material back eciently. Like all crust-related features, however, this mechanism and hence the very occurrence of the second crust-forming era depends on the eciency of melt extraction and crust formation during and directly after the impact, as discussed in the previous paragraph. The enhancement of q after several large impacts induced by foundered lithospheric material raises the possi- CMB bility that a core dynamo is reinforced or reactivated by this process, because that cool material causes locally a steep temperature gradient at the CMB that can potentially drive thermal convection in the core. Such an inhomogeneous temperature pattern at the CMB, albeit on a larger scale, had been envisioned by Stanley et al. (2008) as a possible cause for the dichotomy of the magnetic eld of Mars. If the martian crustal dichotomy has been caused by a giant impact (e.g., Golabek et al., 2011), however, it is conceivable that similar crust production and eclogite formation would occur on a correspondingly larger scale than in our models, and the eect on the CMB would also be more substantial and global in nature. The models by Golabek et al. (2011) did not display this eect, but as discussed above, this may be due to the choice of model parameters such as depth to the basalt{eclogite transition or the thickness of the initial crust, which seem to make delamination less likely in their models; as those authors point out, other parameters, e.g., the extraction eciency or the initial mantle temperature, are also highly uncertain and thus add to the uncertainty in the crustal thickness. Apart from this longer-term, indirect eect of an impact on q it is worth mentioning that some workers also CMB observed a slight enhancement of the CMB heat ux over several dozen millions of years directly after the impact as a consequence of the merging of the cores of the impactor and the target planet (Monteux and Arkani-Hamed, 2014). In principle, modeling this eect requires a more fully dynamical numerical model of core merging than the simple parameterized calculation of heat eects on the mantle used in our models and was therefore omitted; impactors in models based on scaling laws are homogeneous by design. However, the example in gure 16 of Monteux and Arkani-Hamed (2014) suggests that even for impactors that are comparable in size or even larger than the largest ones we considered, the eect would not exceed 5 mW/m and is thus of a magnitude comparable to that induced by foundering lithospheric material. In impactors of the size and with the physical properties we considered, a core is not expected to exceed a few dozen kilometers in size, and its mass would therefore be at least four orders of magnitude lower than that of the martian core. An order-of-magnitude estimate of the heat from the release of gravitational energy and from the shock that would be transported to the martian core by an intact impactor core suggests that the temperature rise of the core would not exceed a few kelvins even under very favorable conditions. In the more likely scenario of a dispersed impactor core (Kendall and Melosh, 2016), the amount of material that eventually ends up in the core will rather be much smaller. Therefore, even in basin-forming impacts core merging is not expected to be signi cant in terms of the CMB heat ow. Williams and Nimmo (2004) used a simple entropy-based criterion for dynamo existence in a non-crystallizing core that considers a cooling rate-dependent entropy term S related to the speci c heat of the core, the entropy production from radioactive heating S , and the entropy change due to conduction along the core adiabat S . Their criterion R k states that a dynamo can exist if S + S > S . The dynamo in our impact-free reference model becomes extinct s R k rather quickly, sooner than commonly thought if the cratering-based timescales for Mars are correct. If some process 17 0 Isidis @ 3810 Myr Ame SWD 150 Sol Sc Dae GC3 Are 350 Sir SWD NTh 150 Hel Dae GC4 Aci IA Ame NTh 150 z +R ic ic GC5 Ama ic z(0.9T ) sol 4250 4200 4150 4100 4050 4000 Age (Myr) Figure 10: Depth range aected by impacts of the great-circle models 3, 4, and 5 in relation to the depth to 0:9T , sol the depth at which 90% of the solidus temperature of fertile martian peridotite is reached. The red spike indicates the depth z of the center of the isobaric core, the blue spike the depth to the deepest point inside the isobaric core, ic z + R , as a measure of the maximum depth extent of shock-heating. The Isidis event, which occurred much later ic ic than all others, is indicated in GC3 with arrows, because it falls outside the plotted time interval. The full set of models can be found in the Supplementary Material. z z z (km) (km) (km) such as massive crust production and subsequent delamination result in the accumulation of eclogite at the CMB, as is the case as a consequence of impacts in our models, the evolution shown in Fig. 6 implies that its in uence on the CMB heat ux and hence on core entropy production could extend the lifetime of a dynamo or revive an already defunct one. The resulting timing would also be in better agreement with accepted dynamo durations for Mars, i.e., until 4:1 Ga (e.g., Lillis et al., 2013; Vervelidou et al., 2017). In the context of discussing global average characteristics of model variables, a few remarks about the characteristics considered here (and in various other studies) will be useful. The issue is the calculation of globally averaged physical properties or variables such as v , T , q , or h in two-dimensional models and how they relate to the rms mean t;CMB ave actual three-dimensional geometry of real planets. Ruedas and Breuer (2017, App. B) already pointed out that the calculation of variables averaged over the surface such as h or q results in an overestimate of the signatures of ave t local anomalies such as a crater, because the linear anomaly in a two-dimensional model occupies a larger fraction of the (linear) surface of the model than the two-dimensional crater area does in the two-dimensional surface of a three- dimensional model. The conversion factor for an anomaly spanning an arc from the 2D model to three dimensions is = 1 cos ; (3) i.e., it is not constant but depends on size and is greater for anomalies up to about three quarters of the circumference in size. Likewise, anomalous volumes in the mantle aect the average disproportionately in two-dimensional models in 2 2 2 comparison to three-dimensional reality. A spherical anomaly with radius R occupies a fraction f = R =(R R ) a 2 a P c 3 3 3 in a 2D model but a fraction f = R =(R R ) in a 3D model. The conversion factor for body variables or anomalies a c is thus 2 2 R R P c = R : (4) 23 a R R For Mars with the structure assumed here, = 2:5163 10 R (with R in kilometers), i.e., for anomalies of the size 23 a a of interest here, the amplitudes should be reduced by 0.5 to 1.5 orders of magnitude; correction factors for the other bodies in the inner Solar System are within a factor of about 2. On the other hand, in three dimensions the larger number of impacts would compensate to some extent the downward correction for the geometry eect by as far as the sum eect is concerned, because a 3D model with all 20 impacts listed in Table 2 would have between 2.5 and 5 times as many impacts as the 2D models carried out here. This should be kept in mind when comparing absolute values with estimates based on 3D models. 5 Conclusions The thermal and compositional evolution of terrestrial planets is in uenced in various ways by large, basin-forming impacts. These impacts are not isolated events but can in uence each other via their dynamical eects, if they are close enough to each other in space and time. Very closely spaced impacts occurring shortly after one another can almost appear like a single larger impact, whereas the interaction is less direct and more complex as the distance and/or the time interval between events grows. Under the assumption of complete extraction of the melt generated by the impact, very large impacts can also trigger lithospheric instabilities that modify the convection ow eld further, and beyond some threshold these instabilities can become so extensive that they stir the mantle on a global scale and potentially reinforce crustal production, including a late stage of long-lasting but low-volume and localized volcanism. Impact-triggered volcanism also enhances the variability of crustal thickness. The remixing of cool crustal material can lead to an intermediate cooling of the mantle that begins some time after the impacts and lasts several hundreds of millions of years before internal heating of this material partly restores previous temperatures. At this scale, the accumulation of lithospheric material at the CMB may also have a reinforcing eect on core dynamics, in particular on a core dynamo. The most realistic but also most expensive way of modeling mantle dynamics is by a three-dimensional spherical model. Among the consequences of models of a two-dimensional cross-section of a sphere are geometric eects that lead to an overestimate of global averages of certain variables or properties in models which feature strong local anomalies. Two-dimensional models can thus accurately simulate physical processes locally and globally, but in order to derive global average characteristics for a 3D planet such as the mean global heat ow from them, the signatures of strong local processes have to be corrected for the geometry eect. With this caveat, they can still provide useful insights into the evolution of three-dimensional systems. Acknowledgments Constructive reviews by two anonymous referees helped us to clarify the paper and consider some additional topics of interest. We appreciate helpful comments on size and age estimates for large basins by Stephanie Werner and on some issues related to impact dynamics by Kai Wunnemann. TR was partly supported by DFG (Deutsche Forschungsgemeinschaft) grant Ru 1839/1-1 and supplementary funding from the DFG programme SFB-TRR 170 and DFG grant Ru 1839/2-1. DB was supported by SFB-TRR 170. This is TRR 170 publication no. 56. The numerical 19 calculations were carried out on the computational resource ForHLR II at the Steinbuch Centre for Computing, Karlsruhe Institute of Technology, funded by the Ministry of Science, Research and the Arts Baden-Wurttem berg and DFG. Fig. 1 was made with the Generic Mapping Tools (Wessel et al., 2013). A Projection of basin positions onto great circles Real crater locations are given by a longitude and a latitude on a spherical surface, but in two-dimensional models, the impacts are implicitly assumed to occur on a great circle. Here we describe how a suitably chosen subset of impact sites can be tted with a great circle as closely as possible to determine impact site locations for use in 2D models. We assume that the selected impacts are close enough to the great circle for them to be mapped onto it in a useful manner. For a great circle with a known pole (in the northern hemisphere by convention), this can be achieved by rst determining the coordinates of the circle's northernmost point: = ; ' = ' ; (5) N P N P where and ' are longitude and latitude, and the subscripts P and N indicate the pole of the great circle and its northernmost point, respectively. The equations of the great circle and its derivative with respect to are then ' = arctan[tan ' cos( )] (6a) N N d' tan ' sin( ) tan ' sin( ) N N N N = = : (6b) d 1 + [tan ' cos( )] 1 + tan ' N N The distance from the actual impact location ( ; ' ) to a point (; ') on the great circle is given by the orthodrome c c arc = arccos[sin ' sin ' + cos ' cos ' cos( )]; (7) c c c and should be minimal for the projection point on the great circle, i.e., its derivative with respect to the longitude should be zero: h i d' d' sin ' cos ' cos ' sin ' cos( ) + cos ' sin( ) c c c c d d d = p 1 [sin ' sin ' + cos ' cos ' cos( )] c c c h i d' d' sin ' cos ' cos ' sin ' cos( ) + cos ' sin( ) c c c c d d = = 0; (8) sin d' where is taken from Eq. 6b. The two existing solutions correspond to the maximum and the minimum distance. The desired minimum solution is found numerically, whereby the maximum distance can easily be excluded by applying the root- nding algorithm in a suciently narrow interval around , assuming that the great circle is a good enough 0 0 t. The coordinate for use in the model can then be given as the arc distance of the projection point ( ; ' ) from c c the northernmost point of the great circle as 0 0 0 s = arccos[sin ' sin ' + cos ' cos ' cos( )]: (9) N N N c c c The placement of the impacts on a great circle in our GC model series uses this relation. Hence, a given impact will have dierent positions in dierent models. B Poiseuille ow with temperature-dependent rheology In order to construct a model simple enough for an analytical solution that captures some features of the upwelling in a hot plume generated by the thermal anomaly of an impact, we assume: 1. that the ow is steady and occurs in a vertical cylindrical pipe with rigid walls at R as in the classical pipe- ow (Poiseuille ow) problem; 2. that the temperature anomaly T is only a function of distance r from the central axis of the pipe; 3. that the decisive in uence on viscosity is the temperature and that it can be described by a Frank-Kamenetskii law, i.e., (T (r)) = exp(bT (r)) with constant parameters and b; 4. that the timescale of the ow process is much shorter than the thermal diusion timescale, such that we can assume a time-independent T (r) that is zero for all r R and T at r = 0. For our max purposes, T (r) = T [1 (r=R) ] is a useful choice that ful lls these criteria. max We start with the same assumptions as in the derivation for the classical Poiseuille ow problem as given by Turcotte and Schubert (1982), i.e., by stating that the force resulting from the pressure dierence p that drives the ow and the shear force on the walls of the cylindrical control volume of length l and radius r cancel out: r dp r p = 2rl , = : (10) 2 dz 20 In our setup, is directly related to the ow velocity v by the ow law, i.e., = dv=dr. The driving pressure results from the thermal buoyancy of the anomaly, which causes a contrast in the vertical lithostatic pressure gradient of dp=dz = g%(r). The density variation is controlled by the thermal expansion due to the temperature anomaly, i.e., % = % T (r); % and are the background density and thermal expansivity, respectively, and are constant. 0 0 Setting Eq. 10 equal to the ow law and inserting the relations for the viscosity, the density and temperature anomaly, and the pressure gradient, we arrive at a dierential equation that can be solved by separation of variables: 2 2 r dv g% T r r 0 max exp bT 1 = 1 (11) 0 max R dr 2 R 2 2 g% R r r ) v(r) = 1 bT 1 exp bT 1 1 : (12) max max 4 b T R R 0 max This ow pro le resembles a Gaussian bell curve and is more focused at the axis for a given R than the classical Poiseuille ow pro le, which is parabolic. The focusing increases with the magnitude of the temperature anomaly and with increasing b. C Simpli ed anomaly volumes In order to estimate the dependence of the impact-generated thermal anomalies on their spatial and temporal sepa- ration, we need the dependence of the anomaly volumes on the separation of the two impacts, x. As we are only interested in a general order-of-magnitude estimate, we approximate the anomalies by two overlapping spheres of radius R whose centers are x apart from each other, i.e., the anomalies are centered at x = R and x = R + x and intersect at x = R + x=2; R can be any reasonable length scale of the anomaly, likely a small multiple of the isobaric core radius. Our anomaly outline is thus de ned by the outline 2 2 2 R (x R) = 2xR x ; 0 x R + x=2 a(x) = (13) 2 2 R (x R x) ; R + x=2 < x 2R + x=2 and zero elsewhere. With respect to the comparison of two- and three-dimensional models, we derive the solutions for both two- and three-dimensional anomalies. As both anomalies are assumed to be of equal size, all integrals are only calculated on the interval [0; R + x=2], and the result is doubled. We only need to consider separations up to 2R, because beyond that distance, there are simply two separate spherical anomalies. The \volume" of the anomaly in two-dimensional models such as our numerical models is the area between a(x) and a(x): R+ " # x x x V = 4 2xR x dx = R 1 + 2 arcsin + R 4R 2R " # x x 2 x = A 1 + 1 + arcsin ; (14) circ R 4R 2R where A is the area of a circle with radius R. circ In a three-dimensional real planet, the idealized anomaly is given by rotating a(x) around the x axis. The rotational body has the volume x x R+ R+ 2 2 Z Z 2 2 V = 2 a (x) dx = 2 (2xR x ) dx 0 0 " # 4 1 3 x 1 x 3 2 3 = R + R x x = V 1 + ; (15) sph 3 12 4 R 16 R where V is the volume of a sphere with radius R. Both functions are strictly monotonic on the interval of interest, sph 0 x=R 2. In their nal forms, they are cast as products of the volume of a simple sphere (in two or three dimensions) and a form factor. The form factor functions for volumes dier by at most 5% between the two- and three-dimensional cases, which suggests that eects that depend on volume follow comparable trends in their dependence on impact spacing in two and three dimensions. In this discussion, we have also assumed that the spheres are fully immersed in the target planet, i.e., they are not truncated by the surface. In the models, we have set the depth of the isobaric core according to the estimate 0:5 imp 0:361 z = 0:1524D v (16) ic imp imp (Ruedas and Breuer, 2018, eq. 6), which is an extension of an empirical relation by Pierazzo et al. (1997). Comparison of the isobaric core diameters and depths in Table 2 shows that this simpli cation is well justi ed for the impacts considered here if R is taken to be of the order of R . ic 21 D Eect of crater overlap on heat ow We consider only complex craters and neglect the curvature of the target body for simplicity. Taking the diameter of the nal crater as the measure of the area in which the heat ow is modi ed, two (identical) impacts will be independent if the distance of their centers is greater than D . For a given impact, our measure of the isobaric core is given by n 1 D = 2r = D b ; (17) ic in imp n + 1 where b and n are empirical parameters that depend on v (Ruedas, 2017). Solving for D , equating with Eq. 1, imp imp and rearranging gives 0:8814 0:377 0:4972 % n + 1 imp imp 0:8814 D = 1:3836 D ; (18) ic 0:2486 0:13 % g D b(n 1) sc 0:8814 0:8814 i.e., D D ; for our setup b = 2:463 and n = 1:424, and hence D 17:9D . Thus, all of our double-impact f f ic ic models except the Utopia-sized ones with x = 10 are expected to lie in the overlap region. It should be kept in mind, however, that D D does not imply that the entire thermal anomaly is shielded by the modi ed crust; the f ic shock-heated region is several times larger than the isobaric core, in which the highest shock pressures are reached. Neglecting curvature, we can use Eq. 14 with R being replaced by D =2 to estimate the eective total impact- modi ed area for overlapping impacts on a sphere: " # D 2x x x A = + 1 + 2 arcsin : (19) 4 D D D f f For our two-dimensional models, the modi ed \surface" is an arc of length x D , and so the eective modi ed \area" increases linearly with x=D up to 2D . f f References Aoki, I., Takahashi, E., 2004. Density of MORB eclogite in the upper mantle. Phys. Earth Planet. Inter. 143{144, 129{143. doi:10.1016/j.pepi.2003.10.007. Archinal, B.A., A'Hearn, M.F.A., Bowell, E., Conrad, A., Consolmagno, G.J., Courtin, R., Fukushima, T., Hestroer, D., Hilton, J.L., Krasinsky, G.A., Neumann, G., Oberst, J., Seidelmann, P.K., Stooke, P., Tholen, D.J., Thomas, P.C., Williams, I.P., 2011. Report of the IAU Working Group on Cartographic Coordinates and Rotational Elements: 2009. Celest. Mech. Dyn. Astr. 109, 101{135. doi:10.1007/s10569-010-9320-4. erratum in 110, 401{403. Arkani-Hamed, J., Ghods, A., 2011. Could giant impacts cripple core dynamos of small terrestrial planets? Icarus 212, 920{934. doi:10.1016/j.icarus.2011.01.020. Bertka, C.M., Fei, Y., 1997. Mineralogy of the Martian interior up to core-mantle boundary pressures. J. Geophys. Res. 102, 5251{5264. doi:10.1029/96JB03270. Bottke, W.F., Andrews-Hanna, J.C., 2017. A post-accretionary lull in large impacts on early Mars. Nature Geosci. 10, 344{348. doi:10.1038/ngeo2937. Carry, B., 2012. Density of asteroids. Planet. Space Sci. 73, 98{118. doi:10.1016/j.pss.2012.03.009. Catling, D.C., 2015. Planetary atmospheres, in: Spohn, T. (Ed.), Physics of Terrestrial Planets and Moons. 2nd ed.. Elsevier. volume 10 of Treatise on Geophysics. chapter 10.13, pp. 429{472. doi:10.1016/B978-0-444-53802-4. 00185-8. Cliord, S.M., 1993. A model for the hydrologic and climatic behavior of water on Mars. J. Geophys. Res. 98, 10973{11016. doi:10.1029/93JE00225. Faul, U.H., 2001. Melt retention and segregation beneath mid-ocean ridges. Nature 410, 920{923. doi:10.1038/ Frey, H., 2008. Ages of very large impact basins on Mars: Implications for the late heavy bombardment in the inner solar system. Geophys. Res. Lett. 35. doi:10.1029/2008GL033515. Frey, H.V., Mannoia, L.M., 2013. A revised, rated and dated inventory of very large candidate impact basins on Mars. Lunar Planet. Sci. 44. Gault, D.E., Wedekind, J.A., 1978. Experimental studies of oblique impact. Proc. Lunar Planet. Sci. Conf. 9, 3843{ 22 Gillmann, C., Golabek, G.J., Tackley, P.J., 2016. Eect of a single large impact on the coupled atmosphere-interior evolution of Venus. Icarus 268, 295{312. doi:10.1016/j.icarus.2015.12.024. Golabek, G.J., Keller, T., Gerya, T.V., Zhu, G., Tackley, P.J., Connolly, J.A.D., 2011. Origin of the martian dichotomy and Tharsis from a giant impact causing massive magmatism. Icarus 215, 346{357. doi:10.1016/j.icarus.2011. 06.012. Hernlund, J.W., Tackley, P.J., 2008. Modeling mantle convection in the spherical annulus. Phys. Earth Planet. Inter. 171, 48{54. doi:10.1016/j.pepi.2008.07.037. Hirth, G., Kohlstedt, D.L., 2003. Rheology of the upper mantle and the mantle wedge: A view from the experimen- talists, in: Eiler, J. (Ed.), Inside the Subduction Factory. American Geophysical Union, Washington, D.C.. volume 138 of AGU Geophysical Monograph, pp. 83{105. Ivanov, B.A., 2001. Mars/Moon cratering rate ratio estimates. Space Sci. Rev. 96, 87{104. doi:10.1023/A: Jacobson, R.A., 2010. The orbits and masses of the martian satellites and the libration of Phobos. Astron. J. 139, 668{679. doi:10.1088/0004-6256/139/2/668. Kendall, J.D., Melosh, H.J., 2016. Dierentiated planetesimal impacts into a terrestrial magma ocean: Fate of the iron core. Earth Planet. Sci. Lett. 448, 24{33. doi:10.1016/j.epsl.2016.05.012. Konopliv, A.S., Asmar, S.W., Folkner, W.M., Karatekin, O., Nunes, D.C., Smrekar, S.E., Yoder, C.F., Zuber, M.T., 2011. Mars high resolution gravity elds from MRO, Mars seasonal gravity, and other dynamical parameters. Icarus 211, 401{428. doi:10.1016/j.icarus.2010.10.004. Lillis, R.J., Robbins, S., Manga, M., Halekas, J.S., Frey, H.V., 2013. Time history of the martian dynamo from crater magnetic eld analysis. J. Geophys. Res. 118, 1488{1511. doi:10.1002/jgre.20105. McKenzie, D., 1985. The extraction of magma from the crust and mantle. Earth Planet. Sci. Lett. 74, 81{91. doi:10.1016/0012-821X(85)90168-2. Melosh, H.J., 2011. Planetary Surface Processes. Number 13 in Cambridge Planetary Science, Cambridge University Press. Monteux, J., Arkani-Hamed, J., 2014. Consequences of giant impacts in early Mars: Core merging and martian dynamo evolution. J. Geophys. Res. 119, 480{505. doi:10.1002/2013JE004587. Nimmo, F., Price, G.D., Brodholt, J., Gubbins, D., 2004. The in uence of potassium on core and geodynamo evolution. Geophys. J. Int. 156, 363{376. doi:10.1111/j.1365-246X.2003.02157.x. O'Neill, C., Marchi, S., Zhang, S., Bottke, W., 2017. Impact-driven subduction on the Hadean Earth. Nature Geosci. 10, 793{797. doi:10.1038/ngeo3029. Padovan, S., Tosi, N., Plesa, A.C., Ruedas, T., 2017. Impact-induced changes in source depth and volume of magmatism on Mercury and their observational signatures. Nat. Comm. 8. doi:10.1038/s41467-017-01692-0. Pierazzo, E., Vickery, A.M., Melosh, H.J., 1997. A reevaluation of impact melt production. Icarus 127, 408{423. doi:10.1006/icar.1997.5713. Reese, C.C., Solomatov, V.S., Baumgardner, J.R., 2002. Survival of impact-induced thermal anomalies in the martian mantle. J. Geophys. Res. 107. doi:10.1029/2000JE001474. Reese, C.C., Solomatov, V.S., Baumgardner, J.R., Stegman, D.R., Vezolainen, A.V., 2004. Magmatic evolution of impact-induced Martian mantle plumes and the origin of Tharsis. J. Geophys. Res. 109. doi:10.1029/2003JE002222. Rivoldini, A., Van Hoolst, T., Verhoeven, O., Mocquet, A., Dehant, V., 2011. Geodesy constraints on the interior structure and composition of Mars. Icarus 213, 451{472. doi:10.1016/j.icarus.2011.03.024. Robbins, S.J., Hynek, B.M., 2012. A new global database of Mars impact craters 1 km: 2. Global crater properties and regional variations of the simple-to-complex transition diameter. J. Geophys. Res. 117. doi:10.1029/2011JE003967. Roberts, J.H., Arkani-Hamed, J., 2012. Impact-induced mantle dynamics on Mars. Icarus 218, 278{289. doi:10.1016/ j.icarus.2011.11.038. Roberts, J.H., Arkani-Hamed, J., 2014. Impact heating and coupled core cooling and mantle dynamics on Mars. J. Geophys. Res. 119, 729{744. doi:10.1002/2013JE004603. Roberts, J.H., Barnouin, O.S., 2012. The eect of the Caloris impact on the mantle dynamics and volcanism of Mercury. J. Geophys. Res. 117. doi:10.1029/2011JE003876. 23 Roberts, J.H., Lillis, R.J., Manga, M., 2009. Giant impacts on early Mars and the cessation of the Martian dynamo. J. Geophys. Res. 114. doi:10.1029/2008JE003287. Rolf, T., Zhu, M.H., Wunnemann, K., Werner, S.C., 2017. The role of impact bombardment history in lunar evolution. Icarus 286, 138{152. doi:10.1016/j.icarus.2016.10.007. Rosenthal, A., Yaxley, G.M., Green, D.H., Hermann, J., Kov acs, I., Spandler, C., 2014. Continuous eclogite melting and variable refertilisation in upwelling heterogeneous mantle. Sci. Rep. 4. doi:10.1038/srep06099. Ruedas, T., 2017. Globally smooth approximations for shock pressure decay in impacts. Icarus 289, 22{33. doi:10. 1016/j.icarus.2017.02.008. Ruedas, T., Breuer, D., 2017. On the relative importance of thermal and chemical buoyancy in regular and impact- induced melting in a Mars-like planet. J. Geophys. Res. 122, 1554{1579. doi:10.1002/2016JE005221. Ruedas, T., Breuer, D., 2018. \Isocrater" impacts: Conditions and mantle dynamical responses for dierent impactor types. Icarus 306, 94{115. doi:10.1016/j.icarus.2018.02.005. Ruedas, T., Tackley, P.J., Solomon, S.C., 2013. Thermal and compositional evolution of the martian mantle: Eects of phase transitions and melting. Phys. Earth Planet. Inter. 216, 32{58. doi:10.1016/j.pepi.2012.12.002. Stanley, S., Elkins-Tanton, L., Zuber, M.T., Parmentier, E.M., 2008. Mars' paleomagnetic eld as the result of a single-hemisphere dynamo. Science 321, 1822{1825. doi:10.1126/science.1161119. Tackley, P.J., 1996. Eects of strongly variable viscosity on three-dimensional compressible convection in planetary mantles. J. Geophys. Res. 101, 3311{3332. doi:10.1029/95JB03211. Tackley, P.J., 2008. Modelling compressible mantle convection with large viscosity contrasts in a three-dimensional spherical shell using the yin-yang grid. Phys. Earth Planet. Inter. 171, 7{18. doi:10.1016/j.pepi.2008.08.005. Taylor, S.R., McLennan, S.M., 2009. Planetary Crusts. Cambridge University Press. Turcotte, D.L., Schubert, G., 1982. Geodynamics | Applications of Continuum Physics to Geological Problems. 2nd ed., John Wiley & Sons. Vervelidou, F., Lesur, V., Grott, M., Morschhauser, A., Lillis, R.J., 2017. Constraining the date of the martian dynamo shutdown by means of crater magnetization signatures. J. Geophys. Res. 122, 2294{2311. doi:10.1002/ 2017JE005410. W anke, H., Dreibus, G., 1994. Chemistry and accretion history of Mars. Phil. Trans. R. Soc. Lond. A 349, 285{293. doi:10.1098/rsta.1994.0132. Watters, W.A., Zuber, M.T., Hager, B.H., 2009. Thermal perturbations caused by large impacts and consequences for mantle convection. J. Geophys. Res. 114. doi:10.1029/2007JE002964. Werner, S.C., 2008. The early martian evolution { Constraints from basin formation ages. Icarus 195, 45{60. doi:10. 1016/j.icarus.2007.12.008. Wessel, P., Smith, W.H.F., Scharroo, R., Luis, J., Wobbe, F., 2013. Generic Mapping Tools: Improved version released. Eos Trans. AGU 94, 409{410. doi:10.1002/2013EO450001. Williams, J.P., Nimmo, F., 2004. Thermal evolution of the Martian core: Implications for an early dynamo. Geology 32, 97{100. doi:10.1130/G19975.1. Yeom, J., 2015. Rs01b1t0 solid sphere with suddenly applied temperature at outer surface. URL: http://exact. unl.edu/exact/contents/contents.php. Zhao, Y.H., Zimmerman, M.E., Kohlstedt, D.L., 2009. Eect of iron content on the creep behavior of olivine: 1. Anhydrous conditions. Earth Planet. Sci. Lett. 287, 229{240. doi:10.1016/j.epsl.2009.08.006. http://www.deepdyve.com/assets/images/DeepDyve-Logo-lg.png Astrophysics arXiv (Cornell University) http://www.deepdyve.com/lp/arxiv-cornell-university/dynamical-effects-of-multiple-impacts-large-impacts-on-a-mars-like-Ws7kFb4ExN

Loading next page...

References (62)

S. Robbins, B. Hynek (2012)
A new global database of Mars impact craters ≥1 km: 2. Global crater properties and regional variations of the simple‐to‐complex transition diameter
Journal of Geophysical Research, 117
F. Nimmo, G. Price, J. Brodholt, D. Gubbins (2003)
The influence of potassium on core and geodynamo evolution
Geophysical Journal International, 156
J. Roberts, J. Arkani‐Hamed (2013)
Impact heating and coupled core cooling and mantle dynamics on Mars
Journal of Geophysical Research: Planets, 119
J. Monteux, J. Arkani‐Hamed (2014)
Consequences of giant impacts in early Mars: Core merging and Martian dynamo evolution
Journal of Geophysical Research: Planets, 119
T. Rolf, Meng‐Hua Zhu, K. Wünnemann, S. Werner (2016)
The role of impact bombardment history in lunar evolution
Icarus, 286
H. Wänke, G. Dreibus (1994)
Chemistry and accretion history of Mars
Philosophical Transactions of the Royal Society of London. Series A: Physical and Engineering Sciences, 349
(2015)
10.13 – Planetary Atmospheres
, 10
C. Reese, V. Solomatov, J. Baumgardner (2002)
Survival of impact-induced thermal anomalies in the Martian mantle
Journal of Geophysical Research, 107
G. Hirth, D. Kohlstedt (2013)
Rheology of the Upper Mantle and the Mantle Wedge: A View from the Experimentalists
Geophysical monograph, 138
J. Roberts, R. Lillis, M. Manga (2009)
Giant impacts on early Mars and the cessation of the Martian dynamo
Journal of Geophysical Research, 114
P. Tackley (2008)
Modelling compressible mantle convection with large viscosity contrasts in a three-dimensional spherical shell using the yin-yang grid
Physics of the Earth and Planetary Interiors, 171
J. Kendall, H. Melosh (2016)
Differentiated planetesimal impacts into a terrestrial magma ocean: Fate of the iron core
Earth and Planetary Science Letters, 448
A. Rosenthal, G. Yaxley, D. Green, J. Hermann, I. Kovács, C. Spandler (2014)
Continuous eclogite melting and variable refertilisation in upwelling heterogeneous mantle
Scientific Reports, 4
P. Tackley (1996)
Effects of strongly variable viscosity on three‐dimensional compressible convection in planetary mantles
Journal of Geophysical Research, 101
(2004)
2004), and as the crustal fragments form a close union with the lighter harzburgitic
(2013)
A revised, rated and dated inventory of very large candidate impact basins on Mars
Jean‐Pierre Williams, F. Nimmo (2004)
Thermal evolution of the Martian core: Implications for an early dynamo
Geology, 32
(2015)
Rs01b1t0 solid sphere with suddenly applied temperature at outer surface
A. Rivoldini, T. Hoolst, O. Verhoeven, A. Mocquet, V. Dehant (2011)
Geodesy constraints on the interior structure and composition of Mars
Icarus, 213
R. Jacobson (2010)
THE ORBITS AND MASSES OF THE MARTIAN SATELLITES AND THE LIBRATION OF PHOBOS
The Astronomical Journal, 139
Philip Bland, S. Schwenzer (2018)
Planetary Surface Processes
P. Wessel, W. Smith, R. Scharroo, J. Luis, F. Wobbe (2013)
Generic Mapping Tools: Improved Version Released
Eos, Transactions American Geophysical Union, 94
U. Faul (2001)
Melt retention and segregation beneath mid-ocean ridges
Nature, 410
D. Turcotte, G. Schubert (1982)
Geodynamics : applications of continuum physics to geological problems
H. Frey (2008)
Ages of very large impact basins on Mars: Implications for the late heavy bombardment in the inner solar system
Geophysical Research Letters, 35
C. Bertka, Y. Fei (1997)
Mineralogy of the Martian interior up to core‐mantle boundary pressures
Journal of Geophysical Research, 102
S. Werner (2008)
The early martian evolution—Constraints from basin formation ages
Icarus, 195
Ivanov (2001)
Mars/Moon cratering rate ratio estimates
Space Sci. Rev., 96
W. Bottke, J. Andrews‐Hanna (2017)
A post-accretionary lull in large impacts on early Mars
Nature Geoscience, 10
C. O'Neill, S. Marchi, Siqi Zhang, W. Bottke (2017)
Impact-driven subduction on the Hadean Earth
Nature Geoscience, 10
C. Reese, V. Solomatov, J. Baumgardner, D. Stegman, A. Vezolainen (2004)
Magmatic evolution of impact-induced Martian mantle plumes and the origin of Tharsis
Journal of Geophysical Research, 109
Y. Zhao, M. Zimmerman, D. Kohlstedt (2009)
Effect of iron content on the creep behavior of olivine: 1. Anhydrous conditions
Earth and Planetary Science Letters, 287
T. Ruedas (2017)
Globally smooth approximations for shock pressure decay in impacts
Icarus, 289
T. Ruedas, P. Tackley, S. Solomon (2012)
Thermal and compositional evolution of the martian mantle: Effects of phase transitions and melting
Physics of the Earth and Planetary Interiors, 216
S. Padovan, N. Tosi, A. Plesa, T. Ruedas (2017)
Impact-induced changes in source depth and volume of magmatism on Mercury and their observational signatures
Nature Communications, 8
J. Roberts, O. Barnouin (2011)
The effect of the Caloris impact on the mantle dynamics and volcanism of Mercury
Journal of Geophysical Research, 117
G. Golabek, G. Golabek, T. Keller, T. Gerya, G. Zhu, P. Tackley, J. Connolly (2011)
Origin of the martian dichotomy and Tharsis from a giant impact causing massive magmatism
Icarus, 215
J. Roberts, J. Arkani‐Hamed (2011)
Impact-induced mantle dynamics on Mars
Icarus, 218
Planetary Crusts (2007)
Planetary Crusts
(2009)
dynamics and for easier comparison with the study
D. McKenzie (1985)
The extraction of magma from the crust and mantle
Earth and Planetary Science Letters, 74
R. Lillis, S. Robbins, M. Manga, J. Halekas, H. Frey (2013)
Time history of the Martian dynamo from crater magnetic field analysis
Journal of Geophysical Research: Planets, 118
A. Konopliv, S. Asmar, W. Folkner, Ö. Karatekin, D. Nunes, S. Smrekar, C. Yoder, M. Zuber (2011)
Mars high resolution gravity fields from MRO, Mars seasonal gravity, and other dynamical parameters
Icarus, 211
D. Gault, J. Wedekind (1978)
Experimental studies of oblique impact.
, 3
C. Gillmann, G. Golabek, P. Tackley (2016)
Effect of a single large impact on the coupled atmosphere-interior evolution of Venus
Icarus, 268
J. Arkani‐Hamed, J. Arkani‐Hamed, A. Ghods (2011)
Could giant impacts cripple core dynamos of small terrestrial planets
Icarus, 212
Y. Zhao, M. Zimmerman, D. Kohlstedt (2017)
Effect of iron content on the creep behavior of Olivine: 2. Hydrous conditions
Physics of the Earth and Planetary Interiors, 278
W. Watters, M. Zuber, B. Hager (2003)
Supplementary materials for Paper 2007JE002964 Thermal perturbations caused by large impacts and consequences for mantle convection
B. Archinal, M. A’Hearn, E. Bowell, A. Conrad, G. Consolmagno, R. Courtin, T. Fukushima, D. Hestroffer, J. Hilton, G. Krasinsky, G. Neumann, J. Oberst, P. Seidelmann, P. Stooke, D. Tholen, P. Thomas, I. Williams (2011)
Report of the IAU Working Group on Cartographic Coordinates and Rotational Elements: 2009
Celestial Mechanics and Dynamical Astronomy, 109
F. Vervelidou, V. Lesur, M. Grott, A. Morschhauser, R. Lillis (2016)
Constraining the Date of the Martian Dynamo Shutdown by Means of Crater Magnetization Signatures
Journal of Geophysical Research: Planets, 122
B. Carry (2012)
Density of asteroids
Planetary and Space Science, 73
I. Aoki, E. Takahashi (2004)
Density of MORB eclogite in the upper mantle
Physics of the Earth and Planetary Interiors, 143
J. Hernlund, P. Tackley (2008)
Modeling mantle convection in the spherical annulus
Physics of the Earth and Planetary Interiors, 171
R. Redman (1950)
Planetary Atmospheres
Nature, 165
S. Clifford (1993)
A model for the hydrologic and climatic behavior of water on Mars
Journal of Geophysical Research, 98
T. Ruedas, D. Breuer (2018)
"Isocrater" impacts: Conditions and mantle dynamical responses for different impactor types
Icarus, 306
S. Stanley, L. Elkins‐Tanton, M. Zuber, E. Parmentier (2008)
Mars' Paleomagnetic Field as the Result of a Single-Hemisphere Dynamo
Science, 321
B. Ivanov (2001)
Mars/Moon Cratering Rate Ratio Estimates
Space Science Reviews, 96
E. Pierazzo, A. Vickery, H. Melosh (1997)
A Reevaluation of Impact Melt Production
Icarus, 127
M. Davies, V. Abalakin, C. Cross, R. Duncombe, H. Masursky, B. Morando, T. Owen, P. Seidelmann, A. Sinclair, G. Wilkins, Y. Tjuflin (1980)
Report of the IAU Working Group on Cartographic Coordinates and Rotational Elements of the Planets and Satellites
Celestial mechanics, 22
T. Ruedas, D. Breuer (2017)
On the relative importance of thermal and chemical buoyancy in regular and impact‐induced melting in a Mars‐like planet
Journal of Geophysical Research: Planets, 122
(2004)
surface” is an arc of length x ≥ Df , and so the effective modified “area” increases linearly with ∆x/Df up to 2Df

ISSN: 0031-9201
eISSN: ARCH-3330
DOI: 10.1016/j.pepi.2019.01.003
Publisher site: See Article on Publisher Site

Abstract

1,2,3 Thomas Ruedas Doris Breuer Dynamical eects of multiple impacts: Large impacts on a Mars-like planet nal version 9 January 2019 published in: Phys. Earth Planet. Inter. 287, pp. 76{92 (2019) Institute of Planetary Research, German Aerospace Center (DLR), Berlin, Germany Museum fur Naturkunde Berlin, Germany Institute of Planetology, Westf alische Wilhelms-Universit at, Munster, Germany The version of record is available at http://dx.doi.org/10.1016/j.pepi.2019.01.003. This author post-print version is shared under the Creative Commons Attribution Non-Commercial No Derivatives License (CC BY-NC-ND 4.0). arXiv:1901.09463v1 [astro-ph.EP] 27 Jan 2019 Dynamical eects of multiple impacts: Large impacts on a Mars-like planet Thomas Ruedas Institute of Planetary Research, German Aerospace Center (DLR), Berlin, Germany Museum fur Naturkunde Berlin, Germany Institute of Planetology, Westfalische Wilhelms-Universitat, Munster, Germany Doris Breuer Institute of Planetary Research, German Aerospace Center (DLR), Berlin, Germany Highlights • Combined eects are not simply linear superpositions of individual impacts. • Many details of impact history are not preserved in the long-term evolution. • Global signatures of dierent impact histories converge in the long term. • Lithospheric delamination may prolong melt generation, aect core. Abstract The earliest stage of the evolution of a fully assembled planet is profoundly aected by a number of basin-forming impacts large enough to change the dynamics of its deeper interior. These impacts are in some cases quite closely spaced and follow one another in short time intervals, so that their eects interact and result in behavior that may dier from a simple sum of the eects of two individual and isolated impacts. We use two-dimensional models of mantle convection in a Mars-like planet and a simple parameterized representation of the principal eects of impacts to study some of the dynamical eects and interactions of multiple large impacts. In models of only two impacts, we con rm that the dynamical eects of the impacts reinforce each other the closer they are in space and time but that the eects do not always correspond to straightforward superpositions of those of single, isolated impacts. In models with multiple (4{8) impacts with variable sizes, distances, and frequencies, the global response of the mantle is as variable as the impact sequences in the short term, but in the long term the dierent evolutionary paths converge for several indicator variables such as the mean ow velocity, temperature, or heat ow. Nonetheless, beyond a certain impact frequency and energy, lithospheric instabilities triggered by large impacts occur on a global scale, reinvigorate mantle dynamics for long time spans, and entail a late stage of melt production in addition to the initial melting stage that is not observed in one- or two-impact models. After one or several very large impacts, some lithospheric material may founder and sink to the core{mantle boundary, and if enough of it accumulates there, it enhances the heat ux out of the core for several hundred millions of years, with possible eects on dynamo activity. Mars; mantle convection; impact 1 Introduction In recent years, the interactions between very large meteorite impacts and planetary interiors have received increasing attention, and several authors have studied a variety of aspects of impact-induced mantle dynamics especially on Mars (e.g., Reese et al., 2002, 2004; Watters et al., 2009; Roberts et al., 2009; Arkani-Hamed and Ghods, 2011; Golabek et al., 2011; Roberts and Arkani-Hamed, 2012; Ruedas and Breuer, 2017, 2018), but also on Mercury (Roberts and Barnouin, 2012; Padovan et al., 2017), Venus (e.g., Gillmann et al., 2016), Earth (e.g., O'Neill et al., 2017), and the Moon (Rolf et al., 2017). Most of these studies concentrate on the eects of a single large impact in order to isolate the physical eects more clearly. However, in the earliest history of the terrestrial planets, basin-forming impacts have followed each other in close succession, and at times also in close proximity, which raises the question whether and how the eects of subsequent large impacts on the deep interior of the target planet may in uence each other. Roberts et al. (2009); Roberts and Arkani-Hamed (2012) and Rolf et al. (2017) have already explored this possibility to some extent for multiple impacts on Mars and the Moon, respectively, and observed that a succession of large impacts leaves a cumulative signature on physical variables such as the heat ow. Corresponding author: T. Ruedas, Institute of Planetary Research, German Aerospace Center (DLR), Berlin, Germany (Thomas.Ruedas@dlr.de) 2 In this study we consider dierent successions of impacts, both idealized ones and examples derived from observed martian giant impacts. Our principal goal is to elucidate some more details of how impacts in close temporal or spatial proximity may in uence each other and whether the cumulative eects of a succession of several impacts are additive or lead to some sort of saturation, e.g., because the aected mantle is so depleted in fusible mineral components that it ceases to produce melt readily. 2 Method and model set-up The evolution of the martian interior including concomitant processes such as melting is represented as a two- dimensional fully dynamical model using the convection code STAGYY (Tackley, 1996, 2008) with the additions and modi cations described in detail in Ruedas and Breuer (2017, 2018) on a two-dimensional spherical annulus grid (Hernlund and Tackley, 2008) with 512×128 points. The convection algorithm solves the equations of mass, mo- mentum, and energy conservation in the anelastic, compressible approximation, making use of a detailed petrological model that includes, among others, melting-related density changes and mineral solid-phase transitions. Following the preferred interior models by Konopliv et al. (2011) and Rivoldini et al. (2011), we de ne our models to consist essentially of two layers, i.e., an upper mantle with an olivine-dominated mineralogy and a lower mantle with an assemblage mostly of a high-pressure olivine polymorph and majorite garnet (Bertka and Fei, 1997); there is no basal (pv+fp) layer. The chemical model assumed for the mantle is the one by W anke and Dreibus (1994), and for the core, we use a sulfur content derived from Rivoldini et al. (2011). The viscosity of the mantle is described by an extended Arrhenius-type law for diusion creep of olivine that includes dependencies on temperature T , pressure p, and chemical composition, in particular water content (Hirth and Kohlstedt, 2003; Zhao et al., 2009); its value at surface pressure and the initial potential temperature of the mantle for the bulk water and iron contents of unmolten mantle material is also the scaling viscosity of the model. Internal heating is provided by the radioactive decay of 40 232 235 238 K, Th, U, and U, which are tracked with tracer particles and are redistributed by melting, by which they become concentrated in the melt and in the crust forming from it; the movement and partitioning of water is treated in the same manner. The surface and core{mantle boundary (CMB) are implemented as isothermal free-slip bound- aries, whereby the temperature at the CMB is allowed to evolve with time according to a simple energy balance and evolution model of the core; the core model largely follows the approach by Nimmo et al. (2004) and Williams and Nimmo (2004) and is described in some detail in Ruedas et al. (2013, App. B). The model is initialized on the basis of a radially symmetric adiabatic state prevailing at 4.5 Ga with a core superheated by 150 K and a surface at 218 K; the surface temperature is kept constant in all models over the entire evolution, as it is not thought to have changed strongly enough to have had a feedback eect on the interior. The actual initial state is derived from this con guration by setting the supersolidus region to the solidus, extracting the corresponding amount of melt to form an initial crust of 50 km, and redistributing the depletion of the mantle residue randomly throughout the mantle; the latter lets the mantle appear well-mixed and compositionally homogeneous on large scales. The thickness of this primordial crust is not well constrained, but with the extraction threshold and the mantle solidus xed, it is mostly dependent on temperature. Subsequent regular melting in the course of the evolution is controlled by the extent to which the temperature exceeds the local solidus; if the amount of newly formed melt passes a certain threshold value and has a supersolidus vertical path to the base of the lithosphere, it is assumed to be extractable, and the excess amount is instantaneously moved to the top to build the crust. The extraction threshold has been estimated to lie between 0.1% (e.g., McKenzie, 1985) and 1{2% (e.g., Faul, 2001), and so we chose a value of 0.7%. The solidus is the same as in Ruedas and Breuer (2017, App. A). The most important model parameters are summarized in Table 1. Further technical detail may be found in Ruedas and Breuer (2017, 2018) and their supplementary material. The impact process is not included via coupling or interfacing the convection algorithm with a fully dynamical hydrocode simulation, but is described using a combination of simpler analytical or semi-empirical approximations. As the impacts we consider are modeled after real impact basins on Mars, we start from the known diameter of the nal crater, D (cf. Table 2), and use scaling laws for complex craters to work back to determine the diameter of the impactor: 0:427 0:28 g % 1:1346 0:1475 D = 0:6918D D (1) imp f sc 0:56 v % imp z;imp (e.g., Melosh, 2011; Ruedas and Breuer, 2018). Here, v = v sin is the vertical component of the impactor z;imp imp velocity, and D is the diameter of the transition between simple and complex craters; in the absence of more speci c sc information, the angle is always set to 45°, which is the most probable impact angle for an isotropic impactor ux (e.g., Gault and Wedekind, 1978). The density and the absolute value of the impactor velocity are set to those of a homogeneous rocky asteroidal object hitting Mars (cf. Table 1; Carry, 2012; Ivanov, 2001). With regard to the energy input from the impact, we follow the approach outlined by Watters et al. (2009) in many respects, i.e., we base the estimate of the shock pressure on a one-dimensional impedance-match calculation and derive the waste heat from that value as a function of distance from the impact center using the \inverse-r" pressure decay formula from Ruedas (2017). The temperature in the shock-heated volume is forced to a value just above the solidus after the production and extraction of the impact melt, for which an upper limit is imposed. The melt extracted is instantaneously added to the top of the model. 3 Table 1: Main model parameters Mars Planetary radius, R 3389.5 km Archinal et al. (2011) Mass 6:41712 10 kg Konopliv et al. (2011); Jacobson (2010) Surface gravity, g 3:717 m/s Surface temperature 218 K Catling (2015) Mantle thickness, z 1659.5 km after Rivoldini et al. (2011); Konopliv et al. (2011) Initial potential temperature, T 1700 K pot Initial core superheating 150 K Surface porosity, ' 0.2 Cliord (1993) surf Melt extraction threshold 0.007 Simple/complex crater transition, D 5.6 km Robbins and Hynek (2012, Tab. 3) sc Bulk silicate Mars Mg# 0.75 W anke and Dreibus (1994); Taylor and McLennan (2009) Initial bulk silicate water content 36 ppm W anke and Dreibus (1994) Present-day K content 305 ppm W anke and Dreibus (1994) Present-day Th content 56 ppb W anke and Dreibus (1994) Present-day U content 16 ppb W anke and Dreibus (1994) Core S content 16 wt.% Rivoldini et al. (2011) Impactor Density, % 2720 kg/m Carry (2012) imp Velocity, v 9.6 km/s Ivanov (2001) imp In addition to a reference model without any impacts, this study comprises two model subsets with a total of 32 models. The rst subset consists of models with two impacts of identical size, of which the rst strikes at 4 Ga (i.e., at t = 500 Myr model time) and whose spatial and temporal separation is varied systematically in order to study possible eects. The impacts of the two series of this subset each have the magnitude of either an Isidis-forming (I) or of a Utopia-forming (U) event, respectively. The spatial separation x is measured in multiples of the isobaric core diameter D and is 1, 2, 5, or 10 in the models of either series; the isobaric core diameter is de ned based on ic the position of the in ection point of the \inverse-r" pressure decay curve from Ruedas (2017) (cf. Eq. 17). The temporal separation t is given in multiples of an estimated \decay time" t needed for certain dynamical eects caused by a single impact of a given magnitude to decay to levels close to the pre-impact state. We chose the mean ow velocity v as the main indicator for the dynamical state and consider the decay to be essentially complete when rms the dierence of v between the disturbed model and the impact-free reference drops below 5% of the maximum rms disturbance at the time of impact; the exact choice is somewhat arbitrary, and we rounded the time to full million years for convenience. In the model series, the second impact follows the rst at 0.5, 1, and 2 decay times. x and t are thus functions of the size of the impact, and as Isidis and Utopia have quite dierent magnitudes, so have D ic and t : for Isidis, we have D = 135:4 km and t = 6 Myr, for Utopia, we have D = 385:6 km and t = 12 Myr. In d ic d ic d addition to these two-impact models, each series also includes an additional model with a single impact identical to the rst impact as a further aid in identifying the particular eects of the second impact. The second subset of models includes six models with a succession of four to eight large impacts on Mars taken from the list of 20 large basins by Frey (2008) (cf. Roberts et al., 2009), whereby the impact sequences are chosen such that all events in a given model lie approximately on a great circle (GC) (Fig. 1). The spacing of the impacts corresponds to the approximate position of the basin centers on the great circle, and their timing follows the ages given by Frey (2008); the method of mapping them is outlined in A. This model set thus combines a certain element of randomness with the constraints from observations on Mars. 3 Results All models in this study are derived from an impact-free reference model whose setup is identical with the 36 ppm H O reference models of Ruedas and Breuer (2017, 2018) except that it starts at 4.5 Ga instead of 4.4 Ga. The reference model starts with a regular pattern of upwellings and downwellings that transitions into a more chaotic mode after a few tens of millions of years. After 500{600 million years, the ow reorganizes itself into a more regular large-scale pattern again over the course of several hundreds of millions of years, ending up with a long-lived nal con guration with ve approximately regularly spaced large plumes while cooling steadily and forming a thick thermal boundary layer at the surface. The initially superheated core cools quite rapidly and loses its excess heat in the rst few hundreds of millions of years, but it never begins to crystallize. The divergence of the evolutionary paths of the models thus begins with the rst impact striking in a given model. For the two-impact models, this occurs early in the transition from the chaotic to the nal quasi-static con guration, for the great-circle models it starts at some point during the chaotic stage. Hence, all impacts, and especially the earliest ones, strike a planet with a still thin lithosphere and hot, globally vigorously convecting mantle whose uppermost part is undergoing some regular global melting, albeit at quite low intensity. 4 5 Table 2: Impact basins considered in this study in order of descending age, after Frey (2008) and Roberts et al. (2009). The impactor diameter D and the diameter D imp ic and depth z of the isobaric core are estimates from scaling laws derived from the observed nal basin diameter D . The ages and positions in the list apply only to impacts ic f in the GC models; the position given in the \Position" column are the geographical coordinates, the single degree entries in the GC columns are the positions on the great circle of the GC model determined as described in A, measured clockwise from the bottom of the circle in the cross-section images. Basin/Code Age D D z D Position Model sets imp ic ic f (Myr) (km) (km) (km) (km) GC1 GC2 GC3 GC4 GC5 GC6 U/I-* Amenthes (Ame) 4218 188.3 103.8 51.8 1070 110.6°E/0.9°S 73.0° 271.5° Zephyria (Z) 4210 213.2 117.4 58.6 1193 164.3°E/12.4°S 255.9° 126.9° 257.6° Daedalia (Dae) 4199 526.9 290.4 144.0 2639 228.3°E/26.5°S 188.1° 243.4° Sirenum (Sir) 4196 188.0 103.6 51.8 1069 205.3°E/67.4°S 200.0° SW Daedalia (SWD) 4176 231.9 127.8 63.3 1278 213.9°E/29.4°S 175.3° 235.6° Ares (Are) 4167 680.5 375.0 185.4 3300 343.9°E/4.0°N 84.0° 303.8° 83.5° Amazonis (Ama) 4154 580.9 320.2 158.5 2873 187.9°E/27.1°N 299.5° 350.7° 302.1° In Amazonis (IA) 4152 207.2 114.2 56.5 1156 192.5°E/29.3°N 302.9° 354.8° 306.6° Solis (Sol) 4148 311.1 171.4 85.4 1663 275.3°E/23.8°S 117.5° 230.4° N Tharsis (NTh) 4143 244.8 134.8 67.3 1347 243.6°E/17.6°N 66.0° 290.0° 42.6° Chryse (C) 4140 324.3 178.6 89.0 1725 318.0°E/25.0°N 51.0° Hematite (Hem) 4138 189.1 104.2 51.5 1065 357.8°E/3.2°N 87.6° 93.5° Scopolus (Sc) 4133 439.5 242.2 120.2 2250 81.8°E/6.9°N 276.4° 42.9° Acidalia (Aci) 4132 629.6 347.1 172.0 3087 342.7°E/59.8°N 29.7° 9.6° North Polar (NP) 4124 297.7 164.0 81.8 1600 195.2°E/80.0°N 350.6° Utopia (U) 4111 699.5 385.4 190.5 3380 115.5°E/45.0°N 324.3° (X) SE Elysium (SEE) 4107 256.2 141.1 70.5 1403 170.3°E/3.7°N 272.9° 273.5° Hellas (Hel) 4065 399.5 220.1 109.4 2070 66.4°E/42.3°S 226.6° 133.7° 170.6° Argyre (Arg) 4043 238.2 131.2 65.5 1315 317.5°E/49.0°S 158.6° Isidis (I) 3810 247.5 136.4 67.5 1352 87.8°E/13.4°N 284.9° 46.8° (X) North Polar 60˚ 60˚ Acidalia Utopia In Amazonis Chryse N Tharsis Isidis Amazonis Hematite SE Elysium Scopolus 0˚ 0˚ Ares Amenthes Zephyria Solis SW Daedalia Daedalia Hellas 2 4 Argyre Sirenum −60˚ −60˚ Figure 1: Schematic map in Hammer projection of the locations of the large impact basins on Mars after Frey (2008) listed in Table 2 and the great circles (colored curves) by which they are assigned to 2D models. 3.1 Two-impact models The quasi-instantaneous input of energy into the planetary interior by an impact produces a steep change in the temporal evolution of several variables that serve as indicators of the dynamical state of the mantle, e.g., the global means of the root-mean-square velocity, v , the mantle temperature, T , and the heat ows through the surface rms mean and the CMB, q , and q . This step-like change is followed by an initially rapid rebound towards the pre-impact t CMB state over the next few millions of years (Fig. 2). The signal of the second impact is added to the decaying signature of the rst in dierent ways for dierent variables. To compare the second impact with the rst, we determine the maximum change of these variables due to each impact and calculate the ratio of the second signal to the rst; e.g., for T the second-to- rst ratio of impact eect amplitudes would be mean T (t ) T (t ) mean 2 mean;sngl 2 = ; (2) T (t ) T (t ) mean 1 mean;ref 1 where \sngl" and \ref " indicate the corresponding models with a single impact and the impact-free reference model, respectively, and t and t indicate the times of the rst and the second impact. These ratios are expressed as functions 1 2 0 0 of the normalized distance, x = x=D , and time lapse, t = t=t , respectively. If this ratio is smaller than 1, ic d i.e., if the change of the indicator variable is smaller in the second impact than in the rst, that would suggest that the eect is not simply additive but may tend towards a state of saturation. If the ratio were equal to 1, impacts with the same parameters would have identical eects, implying that the target has no memory of the eects of preceding impacts. is therefore a measure of the superposition eciency of impacts with respect to the indicator variable X . The models show that the superposition eciency for v (t) tends to increase with decreasing x and also, rms although much less clearly, with t (Fig. 3), because the lingering thermal anomaly from the rst impact boosts the upwelling triggered by the second. The physical cause for this boost may be that in closely spaced anomalies whose temperatures decrease away from the central axis, the superposition leads to particularly high peak temperatures near the axis, which have a non-linear positive eect on the viscosity and thus on ow velocities. A simpli ed analytical model for pipe-like ow with temperature-dependent viscosity indicates that a strong focusing of the ow and a speedup by more than an order of magnitude relative to the isoviscous case is possible for peak temperature anomalies 0 0 of a few hundred kelvins (cf. B). (x ; t ) is not linear in both variables, and the data suggest that there may be a slight increase beyond a certain x in the Utopia-size models, which may be due to the increase of the aected total volume. Except for closely spaced impact pairs, < 1. Later peaks in v (t) in Fig. 2 are not directly related v rms to the impacts and re ect other mantle dynamical processes such as rising plumes or delamination of crust. They do not appear in the impact-free reference model and are an indirect, longer-term consequence of local lithospheric instabilities induced by the impacts. To assess the thermal eect, we consider rst a simpli ed setup consisting of two identical spherical anomalies of radius R whose centers are separated by x; for x < 2R, these anomalies would overlap to some extent. Especially for impacts in short succession, the thermal anomaly from the rst strike is still so strong that the joint anomaly after the second one can be treated as a homogeneous volume as an approximation. The magnitude of the thermal anomaly formed by overlapping spheres is proportional to its volume V (x), and the functional form of V points to a sublinear increase with x up to x 2R (cf. C); for larger osets, it should be constant, as both anomalies are independent 6 10 Δt=0.5t Δt=2t d d Reference Single Isidis Double Isidis, Δx=D ic Double Isidis, Δx=10D ic Single Utopia Double Utopia, Δx=D ic Double Utopia, Δx=10D ic 100 500 550 600 650 500 550 600 650 t t (Myr) (Myr) Figure 2: Root-mean-square velocity, mean temperature, and surface heat ow for selected models with two impacts of either Isidis-like or Utopia-like magnitude. The impact-free and the single-impact models are also shown for reference. The Utopia models have been smoothed at t > 530{550 My by averaging over 2 My-intervals for better readability. The ticmarks on the top abscissas mark the times at which impacts occurred in any of the models. A collection of plots for all models can be found in the Supplementary Material. T (K) 2 mean q ) (mW/m v (cm/yr) rms 1.2 Δt'=0.5, Δx'=1, I I 1.15 Δt'=0.5, Δx'=1, U U Δt'=1, Δx'=2, I I Δt'=1, Δx'=2, U U 1.1 Δt'=2, Δx'=5, I I Δt'=2, Δx'=5, U U Δx'=10, I 1.05 Δx'=10, U 0.95 0.9 1.1 0.9 0.8 0.7 0.6 0.5 1.1 0.9 0.8 0.7 0.6 0.5 0.4 0 2 4 6 8 10 0 0.5 1 1.5 2 Δx' Δt' Figure 3: Second-to- rst ratio of impact eect amplitudes (superposition eciencies) of v , T , and q as functions rms mean t 0 0 of x (left column) and t (right column). α α q,t T from each other. Estimating the cooling time of a spherical anomaly from the analytical solution for diusive cooling 2 2 of a solid sphere, one nds that the temporal decay is governed by a term of the form exp( t=R ), for a diusivity (e.g., Yeom, 2015). The level from which the mantle is reheated by the second impact depends on how much the rst anomaly has cooled by the time it strikes, and so we expect a decrease of the thermal superposition eciency that approximately reproduces that trend, whereby the eective diusivity is enhanced by convection within the anomaly and entrainment of cooler material from outside. The trends in Fig. 3 con rm that grows approximately linearly with spatial separation and decreases approximately linearly with temporal separation; the apparently linear form of the latter is due to the limited range of t considered, whose lowest value of 0.5 lies already after the stage of strongest cooling. The increase with x is clearer for the Utopia-size impacts, presumably because they penetrate more deeply into the sublithospheric mantle and are hence less aected by the cold lid. As < 1 in almost all of the two-impact models, the mantles move towards a state of saturation in which further impacts ideally would not entail more heating but transform their energy otherwise, e.g., by melt production. The signature of an impact on the global surface heat ow is proportional to the surface area aected, whereby the strongest eects are expected within the nal crater. Under the assumptions made in our models, we expect that the signatures of multiple impacts are only additive if the distance between two impacts exceeds a certain value, because otherwise the aected areas will overlap to some extent. Impact modi cation of a surface region that has already been modi ed by a previous impact is reduced, because processes like pore lling of regolith by volcanism can only occur once. In general, the superposition eciency decreases approximately linearly with increasing spatial and q;t temporal separation (see Fig. 3 and D); however, for the large Utopia impacts, there is a clear minimum at x = 5 at all t , which indicates that other processes also leave their mark on q (t). The CMB heat ow, by contrast, shows only minute variations as a direct consequence of the impact in our models, and no meaningful dependence of the 0 0 peak amplitudes on x or t is found (cf. Fig. 6, left column). The shock-heated volumes are also depleted in fusible components and are thus compositionally distinct and less dense than pristine mantle material, which contributes to their buoyancy and the reinforcement of convection caused by impacts. On the other hand, a mantle volume that was strongly depleted by extensive melting in the rst impact cannot experience much further depletion by melting from a second nearby impact. Thus, the buoyancy of the ow after the latter should be expected to be more thermally dominated in models with small x , but the potentially diminished compositional contribution to the total buoyancy is apparently not signi cant enough to compensate the other reinforcing eects on v discussed above. As the impact-generated anomalies ascend and spread out beneath rms the lid, they may not only in uence each other dynamically, as seen in the v (t) curves, but will also collide and rms merge, especially if they are closely spaced. In models with widely spaced impacts in close succession, both anomalies develop somewhat independently for a certain time until they have spread far enough below the lid to run into each other. In such cases a piece of normal mantle can get caught between them and induce a downwelling due to its relatively higher density, especially for smaller impacts with less vigorous dynamics. An example for this can be seen 0 0 in the double-Utopia model with x = 10 and t = 0:5 in the lower right block of Fig. 4. Heat and enhanced convection can also result in increased production of melt and crust, but the extreme depletion of the shallower mantle in the impact-aected region counteracts this eect. Global crustal thickness evolution maps h(x; t) and the corresponding cumulative distributions H (h; t) (Figs. 5, S3{6) show that impacts broaden the range of thickness by creating thinner and thicker regions, especially in the large (Utopia-sized) impact model series. In the following, we will normalize h with the mean crustal thickness h at the time t of the rst impact, and we express the 0 0 time elapsed since t in multiples of the decay time t in order to achieve a scale-free representation that might reveal 0 d common trends for impacts of dierent magnitudes. The cumulative distribution function of the crustal thickness H (h; t) quanti es how large the fraction of crust with a thickness no larger than h is at the time t; for example, in the scale-free plots of some selected cases in Fig. 5, H (1:1; t + 0:5t ) = 0:87 means that at half a decay time after 0 d the rst impact, 87% of the entire crust has a thickness of 1:1h or less. The curves in the left column of the gure show H (h) at that time after the second impact at which a part of the crust reaches near-maximum thickness and has a diverse thickness distribution. As the curve for the reference model shows, the thickness distribution in general is essentially unimodal and shows little scatter around the mean value h=h = 1, i.e., the crust has a rather uniform thickness. The unimodality remains largely valid in the models with impacts as well, but the steep increase at 1 is not quite as high, whereas the anks of the curve become longer and higher, re ecting the broadening of the thickness distribution due to production of anomalous crust by impacts. The color-coded maps in the middle and right column of Fig. 5 show how this cumulative distribution evolves through time in a time window around (mostly after) the impact time; the curves in the left column are thus sections of these maps in the vertical direction and are marked as vertical lines in the map plots. The spatial pattern and evolution of post-impact crustal thickness depends on spacing, timing, and size of the impacts. Comparison of the impact-free and the single-impact model shows clearly that the additional melt extraction after the impact results in a marked and lasting increase in h in some regions, especially near the site of a Utopia-sized impact (cf. Fig. S3). In the Utopia series, the second impact then increases the fraction of thickened crust, possibly including some more growth of already thickened regions, but at least in models in which the impacts are further apart in space and time, the most extreme eect is not as long-lived. On the other hand, some parts of the crust become thinner a certain time after the impacts (usually 5t to 15t after t ), which seems to d d 0 be related to delamination of very thick crustal roots, as the increase in the fraction of thinner-than-average crust is related to the gradual disappearance of some of the thickest crust in the H diagram. Note that the rather discrete, step-like thickness increments are partly an artifact of the grid discretization, which limits the accuracy with which 9 double-Isidis double-Utopia 2 2 2 2 2 1 1 1 1 1 2 2 2 2 1 1 1 1 Figure 4: Snapshots of temperature (T , left column) and depletion (f , right column) for the models with the closest and widest impact spacing and the shortest time interval. The impact locations are marked in the upper rows of each block as 1 and 2. The rst impact occurs at 500 My in all models, the second impact follows at 503 My in the double-Isidis models and at 506 My in the double-Utopia models, respectively (upper rows of each block). The evolution immediately following the second impact is shown in the snapshots in the lower rows of each block. 10 the crustal thickness can be determined. In the Isidis series, these features are less clear, and the thinning is much less signi cant as less delamination takes place. The occurrence of extra-thick crust due to the second impact, however, is also visible but has a much more transient character and fairly uniform lifetime of about 10t in the entire parameter range covered. As a consequence of the interaction of the two anomalies, the location of the thickest post-impact crust does not necessarily coincide exactly with the impact site in both series. In all cases investigated, the dierences between the models diminish with time, and the impact signature in the temporal evolution of the system's dynamical variables fades and has disappeared long before the present. In particular, thermal anomalies fade with time and leave no signal in the present-day heat ux. Compositional anomalies, however, are preserved, as pointed out by Ruedas and Breuer (2017), but it would be dicult to draw a sharp boundary between the region of in uence of one impact or the other if the impacts are close enough to overlap. Variations in post-impact crust formation are also re ected in the crustal thickness and may thus preserve a long-term record of impact-induced mantle dynamics. In addition to the eects of impacts at the surface or in the shallow mantle, it has been proposed that impacts also aect a core dynamo that may be active (e.g., Roberts and Arkani-Hamed, 2014). We did not include any direct impact eects such as shock heating of the core in our treatment as did those authors and therefore neglect the potential eects of an impact-generated hot layer at the top of the core whose appearance was proposed by them. However, given that even the largest impacts in our models are somewhat smaller than the one modeled by Roberts and Arkani-Hamed (2014) and that shock heating decreases to nearly zero towards the CMB in our models, we do not expect a signi cant direct eect of the shock on the core. Impact-related variations of the heat ow through the CMB would therefore mostly be induced by post-impact mantle processes. The left panel of Fig. 6 shows the total entropy production of the core, which is a good proxy for q (t) in our models, for the rst billion years in the double-Utopia CMB 0 0 models with the longest time interval between impacts (t = 2); the curves for smaller t look very similar, whereas in the double-Isidis models the deviations from the reference are negligible. At the times of the impacts, there are tiny downward excursions that indicate a short period of reduced entropy production and heat ow out of the core during the presence of the thermal anomaly in the mantle; this eect would be compensated if shock heating also aected the core. The more signi cant eect, however, sets in only several decay times after the second impact and consists in an increased entropy production and heat ow out of the core. It is caused by the arrival at the CMB of cooler lithospheric material into which small fragments of eclogitized crust are embedded and that had become unstable in the aftermath of the impacts, delaminated, and sunk to the bottom of the mantle where it cools the core from above. However, in the lower part of the martian mantle, eclogite has a signi cantly smaller density excess than in the shallow mantle (e.g., Aoki and Takahashi, 2004), and as the crustal fragments form a close union with the lighter harzburgitic residue that went down along with them, the old lithosphere does not accumulate permanently at the CMB. Instead, it is heated enough by its internal sources in the eclogite and by conduction from the surrounding warm mantle to become unstable and be mixed back into the mantle again by the general mantle ow. Hence, while spacing and time interval between impacts are unimportant for post-impact eects on the CMB, the size and, to some degree, the number of the impacts are decisive, because they control the extent of lithospheric instability. 3.2 Great-circle models As the models with multiple impacts on a great circle show, a succession of various dierent large impacts produces a strongly variable depletion pattern in the uppermost mantle (Fig. 7). This pattern still re ects the diversity of the impacts that produce it, but the vigorous post-impact dynamics and merging of the individual anomalies precludes the distinction of clear boundaries between the traces of the discrete events. As the individual impacts in these models lie further apart from each other and span a longer time interval than those in the previous sets, the more indirect eects as a whole, especially those on the lithosphere, also last longer and extend to the entire planet: speci cally, crustal delamination as a long-term consequence of the impacts becomes a planet-wide process and lasts much longer than in the previous sets. Figure 8 shows the root-mean-square convection ow velocity, the mantle mean temperature, and the surface heat ow for models GC3, GC4, and GC5, which represent the models with the highest, an intermediate, and the lowest number of impacts on a great circle in the model set. The basic shape of the curves resembles those from the idealized double-impact models in the previous section. The impacts are again marked by sharp peaks with subsequent gradual decays of these variables; later narrow peaks in v correspond to the quick motion of cold lithospheric instabilities rms formed in the aftermath of the impacts, mostly related to the formation of thick crust. Although we do not normalize times here with the decay time t because there is no single t applicable to the entire duration of a given model, it d d is worth noting that estimates for the t for individual impacts seem to follow the aforementioned relation between t d d 2 2 and the size of the thermal anomaly, exp( t=R ). The center panel of Figure 8 shows that the T of models mean with particularly intense impact activity even drops below the T of the impact-free reference model some time mean after the impacts. The reason for this maybe unexpected behavior is that the heat from the impact shocks has largely diused away after some tens of millions of years, but the copious amounts of cool crust that were produced in their aftermath are becoming unstable and partly sink into the mantle, where they form ecient heat sinks; the crust is so cold, because the melt from which it is formed has been set to the surface temperature upon eruption. At some point between 1.5 Gy and 2.5 Gy, i.e., outside the plotted range, all of the models return to and remain at T values above the reference as the delaminated crust has been warmed up by its own internal radioactive heat mean 11 1 2.2 2.2 No No impact impact reference reference 2 2 0.8 1.8 1.8 0.6 1.6 1.6 1.4 1.4 0.4 1.2 1.2 0.2 1 1 Reference model @ t 0 0.8 0.8 0.8 1 1.2 1.4 1.6 1.8 2 2.2 0 5 10 15 20 25 30 0 5 10 15 20 25 30 h/h (t-t )/t (t-t )/t 0 d 0 d 1 2.2 2.2 Isidis, Utopia, single single 2 2 0.8 1.8 1.8 0.6 1.6 1.6 0.4 1.4 1.4 1.2 1.2 0.2 Single +6t Isidis @ t 0 d 1 1 Single Utopia @ t 0.8 1 1.2 1.4 1.6 1.8 2 2.2 0.8 0.8 0 5 10 15 20 25 30 0 5 10 15 20 25 30 h/h (t-t )/t (t-t )/t 0 d 0 d 2.2 2.2 Double Double Isidis, Utopia, Δx'=1, Δx'=1, Δt'=0.5 Δt'=0.5 2 2 0.8 1.8 1.8 0.6 1.6 1.6 0.4 1.4 1.4 1.2 1.2 0.2 Double +1.5t Isidis, Δx'=1, Δt'=0.5 @ t 0 d 1 1 Double +0.5t Utopia, Δx'=1, Δt'=0.5 @ t 0 d 0.8 0.8 0.8 1 1.2 1.4 1.6 1.8 2 2.2 0 5 10 15 20 25 30 0 5 10 15 20 25 30 h/h (t-t )/t (t-t )/t 0 d 0 d 1 2.2 2.2 Double Double Isidis, Utopia, Δx'=1, Δx'=1, Δt'=2 Δt'=2 2 2 0.8 1.8 1.8 0.6 1.6 1.6 1.4 1.4 0.4 1.2 1.2 0.2 Double +2t Isidis, 1 1 Δx'=1, Δt'=2 @ t 0 d Double +2.5t Utopia, Δx'=1, Δt'=2 @ t 0 d 0.8 0.8 0.8 1 1.2 1.4 1.6 1.8 2 2.2 0 5 10 15 20 25 30 0 5 10 15 20 25 30 h/h (t-t )/t (t-t )/t 0 d 0 d 2.2 2.2 Double Double Isidis, Utopia, Δx'=10, Δx'=10, Δt'=0.5 Δt'=0.5 2 2 0.8 1.8 1.8 0.6 1.6 1.6 0.4 1.4 1.4 1.2 1.2 0.2 Double +6t Isidis, Δx'=10, Δt'=0.5 @ t 0 d 1 1 Double +0.5t Utopia, Δx'=10, Δt'=0.5 @ t 0 d 0.8 0.8 0.8 1 1.2 1.4 1.6 1.8 2 2.2 0 5 10 15 20 25 30 0 5 10 15 20 25 30 h/h (t-t )/t (t-t )/t 0 d 0 d 2.2 2.2 Double Double Isidis, Utopia, Δx'=10, Δx'=10, Δt'=2 Δt'=2 2 2 0.8 1.8 1.8 0.6 1.6 1.6 1.4 1.4 0.4 1.2 1.2 0.2 1 1 Double +6t Isidis, Δx'=10, Δt'=2 @ t 0 d Double +2t Utopia, Δx'=10, Δt'=2 @ t 0 d 0.8 0.8 0.8 1 1.2 1.4 1.6 1.8 2 2.2 0 5 10 15 20 25 30 0 5 10 15 20 25 30 h/h 0 (t-t )/t (t-t )/t 0 d 0 d 0 0 0 0 0 0 0 0 0 0 0 0 0.2 0.2 0.2 0.2 0.2 0.2 0.2 0.2 0.2 0.2 0.2 0.2 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.6 0.6 0.6 0.6 0.6 0.6 0.6 0.6 0.6 0.6 0.6 0.6 0.8 0.8 0.8 0.8 0.8 0.8 0.8 0.8 0.8 0.8 0.8 0.8 1 1 1 1 1 1 1 1 1 1 1 1 H H H H H H H H H H H H Figure 5: Cumulative distribution function H of the crustal thickness of the two-impact model series from Fig. 2 as a function of crustal thickness h and time t at and after the impacts. The left column shows H (h) at that time after the second impact at which the crust reaches its maximum thickness. These curves are thus sections of the H (h; t) maps in the other two columns parallel to the ordinate; they are marked there as vertical lines in the corresponding color. The crustal thickness is normalized with the mean crustal thickness h immediately before the rst impact, i.e., at t = 500 Myr, which can be calculated from the impact-free reference model. The abscissa in the middle and right column gives the time since the rst impact in multiples of the decay timescale t , which is 6 Myr for the Isidis-sized impacts (left column) and 12 Myr for the Utopia-sized impacts (right column). The impact-free and the single-impact models are also shown for comparison. The full set of gures is provided in the Supplementary Material. H H H H H H h/h h/h h/h h/h h/h h/h 0 0 0 0 0 0 h/h h/h h/h h/h h/h h/h 0 0 0 0 0 0 35 Reference Reference Single GC1 double Δx=D Great-circle GC2 Utopia, models Δt=2t d ic Δx=2D GC3 ic Δx=5D GC4 ic Δx=10D GC5 ic GC6 -5 0 200 400 600 800 1000 0 200 400 600 800 1000 t t (Myr) (Myr) Figure 6: Total entropy production of the core for the double-Utopia models with the longest time interval between impacts (left) and for the great-circle models (right); a dynamo can exist if S > 0. The corresponding heat ows at the CMB can be found in the Supplementary Material. The ticmarks on the top abscissa mark the times at which impacts occurred in any of the models. SWD Dae Sir Hel Z Sol SWD Dae GC3 GC4 GC5 Ame Ame NTh Are Sc NTh Ama Aci IA Figure 7: Snapshots of the present-day depletion patterns for the great-circle models 3, 4, and 5. The northernmost points of the respective great circles lie at the six o'clock position. Snapshots of all models are presented in the Supplementary Material. S (MW/K) c production. In all curves, the largest impacts dominate strongly the total deviation from the impact-free reference. In the Introduction, we had speculated that at some point, a planetary mantle could reach a \saturated" state in the sense that additional energy from another impact does not reinvigorate the dynamics, for instance because the shallow mantle is so hot and/or depleted throughout that some sort of stable layering is established that suppresses further whole-mantle convection. However, in none of the models do we nd any clear indication of such a \saturation", which suggests that Mars was always far enough away from it, if such a state is possible at all. Nonetheless, and even though the larger spatial and temporal separation of the individual events diminishes the direct interactions between impacts, the dierent histories of the GC models induce dierent dynamical states of the mantle that lead to dierent signatures between models for a given impact. This becomes most readily visible in the two late impacts, Hellas (GC1, GC4, GC6) and Isidis (GC1, GC3). Both events leave signals of dierent magnitude in the dierent models in which they occur as a consequence of the dierent evolutionary paths (Fig. S7). In all six models, however, the global characteristics return to a state close to that of the reference on timescales that are substantially shorter than the age of the planet. Given the relatively short lifetime of those temperature-related anomalies, it is of interest to consider how longer- lived anomalies of other characteristics compare between great-circle models. Hence, we look again at the statistics of crustal thickness as expressed by H (h; t) (Fig. 9). As in the two-impact models, we observe that thick crust is formed immediately after the impacts but is partly destroyed by delamination on timescales of a few ten to a few hundred million years. However, the more perturbed dynamics of models with several impacts result in more extensive lithospheric instability in the later evolution, which in turn stirs the mantle up and enables prolonged melt and crust production in some places by transporting relatively fresh or re-enriched material into the melting zone; the refertilization of depleted mantle would be a consequence of the entrainment of delaminated crust (e.g., Rosenthal et al., 2014). As a consequence, there is a second era of (limited and localized) crustal production that begins considerable time after the end of the era of the large impacts and generates some sites of stable thick crust that persist to the present, which can be seen in the h(x; t) plots in the left part of Fig. 9. The resurgent crust production is also visible in some of the H (h; t) plots as a new increase in the maximum thicknesses. For the multi-impact sequences of the great-circle models, it is also worthwhile to take a look at q (t), because CMB in most of them we observe substantial lithospheric instability and subsequent transient accumulation of relatively cool material at the CMB. The right panel of Fig. 6 shows that indeed in most of these models core entropy production is increased substantially for several hundred million years after the impacts relative to the impact-free reference model, and q doubles in this time interval. The only exception, in which the increase is much less pronounced, is model CMB GC5, which had only four impacts and was also the least productive one in terms of post-impact crust formation. 4 Discussion Although the very setup of our models is arti cial by design in the two-impact models or due to the restriction in impact count and geometry imposed by the choice of the great circles and the two-dimensional grid, some lessons can be learned with regard to the geological implications of the martian impact record and its uncertainties. There are dierent compilations of large impact basins in the literature, and they show a substantial variance with regard to the timing, the size, and even the very number of impacts, because many of the old basins that predate Utopia are eroded or otherwise obliterated and therefore dicult to identify and characterize. At the one extreme, the survey by Frey and Mannoia (2013) identi ed some 30 large basins falling in the time interval between the formation of the dichotomy and the creation of Utopia, some of them of a magnitude comparable to the Utopia event. At the other extreme lies a recent analysis by Bottke and Andrews-Hanna (2017) that strongly indicates that such a large ux of giant impactors would almost certainly have hit some part of the dichotomy boundary, but as clear evidence for this seems to be lacking, it concludes that the actual number of such events should be less than 12. We decided to use the list by Frey (2008), because of its intermediate number of basin-forming impacts with potential eects on mantle dynamics and for easier comparison with the study by Roberts et al. (2009). Moreover, Werner (2008) identi ed partly dierent basins, assigned them ages that dier by tens or even more than hundred million years from those given in other compilations, and also determined dierent diameters. On the other hand, the global characteristics we chose to assess the dynamical state tend to converge towards the impact-free reference in the long run; for instance, the nal T in the great-circle models is only about 11{17 K higher and q less than 1.3 mW/m lower than in the impact-free mean t reference model. This suggests that the uncertainties in the dierent compilations are not a very strong concern with regard to the overall dynamical and thermal evolution of Mars. In the short-term aftermath of the impacts, however, there are substantial dierences, which would leave traces in the geological record that would still be visible today. With regard to the long-term preservation of the impact record, the triggering of lithospheric instability is of particular importance. A crucial factor in this process is the formation of crust that grows thick enough for its base to cross the basalt{eclogite transition and become gravitationally unstable. With our parameters, this transition would generally take place somewhere around 140 km depth in Mars, which is a bit less than twice the global mean thickness reached after a few hundred million years in the impact-free reference model. That global mean value is rather on the high side of other, independent crustal thickness estimates for Mars and would be lower if a lower initial mantle temperature or a higher melt extraction threshold had been chosen. A lower crustal thickness in our models would imply that lithospheric instability becomes less likely and less extensive or does not occur at all. Likewise, the cooling 14 7 Reference GC3 GC4 GC5 1750 0 300 350 400 450 500 550 600 650 700 t (Myr) Figure 8: Root-mean-square velocity, mean temperature, and surface heat ow for the great-circle models 3, 4, and 5, which include 8, 6, and 4 impacts of varying magnitudes, respectively. The impact-free model is also shown for reference. A collection of plots for all models can be found in the Supplementary Material. T (K) mean q ) (mW/m v (cm/yr) rms Ame Dae Sir SWD Are IA Ama Sol NTh Sc Aci Hel I No impact reference 360 2.2 1.8 1.6 1.4 1.2 0 0.8 0 500 1000 1500 2000 2500 3000 0 500 1000 1500 2000 2500 3000 t-t t-t (Myr) (Myr) GC3 360 2.2 Are 1.8 Sol 1.6 Dae SWD 1.4 Z 120 1.2 Ame 60 Sc I 0 0.8 0 500 1000 1500 2000 2500 3000 0 500 1000 1500 2000 2500 3000 t-t t-t (Myr) (Myr) 0 0 GC4 360 2.2 NTh 1.8 Dae SWD 1.6 Sir 1.4 Hel 120 1.2 Aci 0 0.8 0 500 1000 1500 2000 2500 3000 0 500 1000 1500 2000 2500 3000 t-t t-t (Myr) (Myr) 0 0 GC5 360 2.2 IA Ama Ame 1.8 1.6 1.4 1.2 NTh 0 0.8 0 500 1000 1500 2000 2500 3000 0 500 1000 1500 2000 2500 3000 t-t t-t (Myr) (Myr) 0 0 0.8 0.8 0.8 0.8 1 1 1 1 1.2 1.2 1.2 1.2 1.4 1.4 1.4 1.4 1.6 1.6 1.6 1.6 1.8 1.8 1.8 1.8 2 2 2 2 2.2 2.2 2.2 2.2 0 0 0 0 0.2 0.2 0.2 0.2 0.4 0.4 0.4 0.4 0.6 0.6 0.6 0.6 0.8 0.8 0.8 0.8 1 1 1 1 h/h h/h h/h h/h H H H H Figure 9: Normalized crustal thickness h(x; t)=h (left column) and cumulative distribution function H (h; t) of the crustal thickness (right column) of the impact-free reference model and the great-circle models 3, 4, and 5 (cf. Fig. 8 at and after the impacts. The crustal thickness is normalized with the mean crustal thickness h immediately before the rst impact in Table 2, i.e., Amenthes at 280 Myr, which can be calculated from the impact-free reference model. The abscissa gives the time since the rst impact; contrary to Fig. 5, no normalization has been applied, because the impacts have dierent sizes. The full set of gures is provided in the Supplementary Material. x x x x (°) (°) (°) (°) h/h h/h h/h h/h 0 0 0 0 Ame Ame Z Sir Dae SWD Dae IA Ama SWD Are NTh Sol NTh Aci Sc Hel I eect of the delaminated crust would be absent, and models with impacts would always be hotter at all times than the impact-free reference. Thus, post-impact dynamics also depends crucially on pre-impact crust formation and especially on the thickness of the primordial crust, which is poorly known. This highlights once more the importance of obtaining tight constraints on the thickness of the crust, e.g., from the seismic experiment on the InSight mission. Moreover, at the impact sites themselves the actual extent of extraction of melt generated directly by the impacts determines to a large degree the thickness of the newly formed post-impact crust. Our assumption of very ecient extraction is likely an extreme case and de nes an upper limit more conducive to crustal delamination. By contrast, Padovan et al. (2017) assumed no extraction of impact melt and formed their post-impact crust exclusively from melt produced in the aftermath of the impact by the dynamics of the evolving thermal anomaly, which results in a substantially thinner post-impact crust. The true extent of melt extraction in impacts is currently unknown and requires further study. A straightforward criterion for the occurrence of a second melting era, e.g., in terms of a threshold number of impacts or a minimum size does not seem to exist, even though we nd that the two-impact models and the great- circle model with the fewest impacts (GC5 with only four events) have no such era whereas the other GC models do. The absence of marked thermal anomalies at the time the second crust-forming era begins suggests that its cause is not primarily thermal; indeed, it is only hundreds of millions of years later that the average mantle temperature of the GC models is again higher than that of the impact-free reference model. The critical ingredient is therefore presumably the availability of easily fusible material such as delaminated crust, which however would rst have to be produced from very thick crust by extensive melting. This can happen more easily if the impacts shock-heat an already hot target, which requires them to penetrate more or less the entire lithosphere. Fig. 10 shows the depth reach of the impacts in some of the great-circle models as a function of time and in relation to the depth where 90% of the solidus temperature of fertile martian peridotite is reached; this value serves as a proxy for the shallowest depth at which sublithospheric conditions are reached or at which the mantle is so hot that it will melt without a lot of further energy input. The gure reveals that of the scanty series of impacts in GC5, only one (Amazonis) reaches rmly into the sublithospheric mantle, whereas the other series aect the sublithospheric mantle more strongly and over a longer timespan. On the other hand, the strong eect of the double-Utopia models with two deep-reaching impacts in rather short succession shows that the magnitude alone is not sucient. Instead, a longer gap between the individual events is more likely to ensure that the convective activity is maintained over a longer period and can stir and mix delaminated crustal material back eciently. Like all crust-related features, however, this mechanism and hence the very occurrence of the second crust-forming era depends on the eciency of melt extraction and crust formation during and directly after the impact, as discussed in the previous paragraph. The enhancement of q after several large impacts induced by foundered lithospheric material raises the possi- CMB bility that a core dynamo is reinforced or reactivated by this process, because that cool material causes locally a steep temperature gradient at the CMB that can potentially drive thermal convection in the core. Such an inhomogeneous temperature pattern at the CMB, albeit on a larger scale, had been envisioned by Stanley et al. (2008) as a possible cause for the dichotomy of the magnetic eld of Mars. If the martian crustal dichotomy has been caused by a giant impact (e.g., Golabek et al., 2011), however, it is conceivable that similar crust production and eclogite formation would occur on a correspondingly larger scale than in our models, and the eect on the CMB would also be more substantial and global in nature. The models by Golabek et al. (2011) did not display this eect, but as discussed above, this may be due to the choice of model parameters such as depth to the basalt{eclogite transition or the thickness of the initial crust, which seem to make delamination less likely in their models; as those authors point out, other parameters, e.g., the extraction eciency or the initial mantle temperature, are also highly uncertain and thus add to the uncertainty in the crustal thickness. Apart from this longer-term, indirect eect of an impact on q it is worth mentioning that some workers also CMB observed a slight enhancement of the CMB heat ux over several dozen millions of years directly after the impact as a consequence of the merging of the cores of the impactor and the target planet (Monteux and Arkani-Hamed, 2014). In principle, modeling this eect requires a more fully dynamical numerical model of core merging than the simple parameterized calculation of heat eects on the mantle used in our models and was therefore omitted; impactors in models based on scaling laws are homogeneous by design. However, the example in gure 16 of Monteux and Arkani-Hamed (2014) suggests that even for impactors that are comparable in size or even larger than the largest ones we considered, the eect would not exceed 5 mW/m and is thus of a magnitude comparable to that induced by foundering lithospheric material. In impactors of the size and with the physical properties we considered, a core is not expected to exceed a few dozen kilometers in size, and its mass would therefore be at least four orders of magnitude lower than that of the martian core. An order-of-magnitude estimate of the heat from the release of gravitational energy and from the shock that would be transported to the martian core by an intact impactor core suggests that the temperature rise of the core would not exceed a few kelvins even under very favorable conditions. In the more likely scenario of a dispersed impactor core (Kendall and Melosh, 2016), the amount of material that eventually ends up in the core will rather be much smaller. Therefore, even in basin-forming impacts core merging is not expected to be signi cant in terms of the CMB heat ow. Williams and Nimmo (2004) used a simple entropy-based criterion for dynamo existence in a non-crystallizing core that considers a cooling rate-dependent entropy term S related to the speci c heat of the core, the entropy production from radioactive heating S , and the entropy change due to conduction along the core adiabat S . Their criterion R k states that a dynamo can exist if S + S > S . The dynamo in our impact-free reference model becomes extinct s R k rather quickly, sooner than commonly thought if the cratering-based timescales for Mars are correct. If some process 17 0 Isidis @ 3810 Myr Ame SWD 150 Sol Sc Dae GC3 Are 350 Sir SWD NTh 150 Hel Dae GC4 Aci IA Ame NTh 150 z +R ic ic GC5 Ama ic z(0.9T ) sol 4250 4200 4150 4100 4050 4000 Age (Myr) Figure 10: Depth range aected by impacts of the great-circle models 3, 4, and 5 in relation to the depth to 0:9T , sol the depth at which 90% of the solidus temperature of fertile martian peridotite is reached. The red spike indicates the depth z of the center of the isobaric core, the blue spike the depth to the deepest point inside the isobaric core, ic z + R , as a measure of the maximum depth extent of shock-heating. The Isidis event, which occurred much later ic ic than all others, is indicated in GC3 with arrows, because it falls outside the plotted time interval. The full set of models can be found in the Supplementary Material. z z z (km) (km) (km) such as massive crust production and subsequent delamination result in the accumulation of eclogite at the CMB, as is the case as a consequence of impacts in our models, the evolution shown in Fig. 6 implies that its in uence on the CMB heat ux and hence on core entropy production could extend the lifetime of a dynamo or revive an already defunct one. The resulting timing would also be in better agreement with accepted dynamo durations for Mars, i.e., until 4:1 Ga (e.g., Lillis et al., 2013; Vervelidou et al., 2017). In the context of discussing global average characteristics of model variables, a few remarks about the characteristics considered here (and in various other studies) will be useful. The issue is the calculation of globally averaged physical properties or variables such as v , T , q , or h in two-dimensional models and how they relate to the rms mean t;CMB ave actual three-dimensional geometry of real planets. Ruedas and Breuer (2017, App. B) already pointed out that the calculation of variables averaged over the surface such as h or q results in an overestimate of the signatures of ave t local anomalies such as a crater, because the linear anomaly in a two-dimensional model occupies a larger fraction of the (linear) surface of the model than the two-dimensional crater area does in the two-dimensional surface of a three- dimensional model. The conversion factor for an anomaly spanning an arc from the 2D model to three dimensions is = 1 cos ; (3) i.e., it is not constant but depends on size and is greater for anomalies up to about three quarters of the circumference in size. Likewise, anomalous volumes in the mantle aect the average disproportionately in two-dimensional models in 2 2 2 comparison to three-dimensional reality. A spherical anomaly with radius R occupies a fraction f = R =(R R ) a 2 a P c 3 3 3 in a 2D model but a fraction f = R =(R R ) in a 3D model. The conversion factor for body variables or anomalies a c is thus 2 2 R R P c = R : (4) 23 a R R For Mars with the structure assumed here, = 2:5163 10 R (with R in kilometers), i.e., for anomalies of the size 23 a a of interest here, the amplitudes should be reduced by 0.5 to 1.5 orders of magnitude; correction factors for the other bodies in the inner Solar System are within a factor of about 2. On the other hand, in three dimensions the larger number of impacts would compensate to some extent the downward correction for the geometry eect by as far as the sum eect is concerned, because a 3D model with all 20 impacts listed in Table 2 would have between 2.5 and 5 times as many impacts as the 2D models carried out here. This should be kept in mind when comparing absolute values with estimates based on 3D models. 5 Conclusions The thermal and compositional evolution of terrestrial planets is in uenced in various ways by large, basin-forming impacts. These impacts are not isolated events but can in uence each other via their dynamical eects, if they are close enough to each other in space and time. Very closely spaced impacts occurring shortly after one another can almost appear like a single larger impact, whereas the interaction is less direct and more complex as the distance and/or the time interval between events grows. Under the assumption of complete extraction of the melt generated by the impact, very large impacts can also trigger lithospheric instabilities that modify the convection ow eld further, and beyond some threshold these instabilities can become so extensive that they stir the mantle on a global scale and potentially reinforce crustal production, including a late stage of long-lasting but low-volume and localized volcanism. Impact-triggered volcanism also enhances the variability of crustal thickness. The remixing of cool crustal material can lead to an intermediate cooling of the mantle that begins some time after the impacts and lasts several hundreds of millions of years before internal heating of this material partly restores previous temperatures. At this scale, the accumulation of lithospheric material at the CMB may also have a reinforcing eect on core dynamics, in particular on a core dynamo. The most realistic but also most expensive way of modeling mantle dynamics is by a three-dimensional spherical model. Among the consequences of models of a two-dimensional cross-section of a sphere are geometric eects that lead to an overestimate of global averages of certain variables or properties in models which feature strong local anomalies. Two-dimensional models can thus accurately simulate physical processes locally and globally, but in order to derive global average characteristics for a 3D planet such as the mean global heat ow from them, the signatures of strong local processes have to be corrected for the geometry eect. With this caveat, they can still provide useful insights into the evolution of three-dimensional systems. Acknowledgments Constructive reviews by two anonymous referees helped us to clarify the paper and consider some additional topics of interest. We appreciate helpful comments on size and age estimates for large basins by Stephanie Werner and on some issues related to impact dynamics by Kai Wunnemann. TR was partly supported by DFG (Deutsche Forschungsgemeinschaft) grant Ru 1839/1-1 and supplementary funding from the DFG programme SFB-TRR 170 and DFG grant Ru 1839/2-1. DB was supported by SFB-TRR 170. This is TRR 170 publication no. 56. The numerical 19 calculations were carried out on the computational resource ForHLR II at the Steinbuch Centre for Computing, Karlsruhe Institute of Technology, funded by the Ministry of Science, Research and the Arts Baden-Wurttem berg and DFG. Fig. 1 was made with the Generic Mapping Tools (Wessel et al., 2013). A Projection of basin positions onto great circles Real crater locations are given by a longitude and a latitude on a spherical surface, but in two-dimensional models, the impacts are implicitly assumed to occur on a great circle. Here we describe how a suitably chosen subset of impact sites can be tted with a great circle as closely as possible to determine impact site locations for use in 2D models. We assume that the selected impacts are close enough to the great circle for them to be mapped onto it in a useful manner. For a great circle with a known pole (in the northern hemisphere by convention), this can be achieved by rst determining the coordinates of the circle's northernmost point: = ; ' = ' ; (5) N P N P where and ' are longitude and latitude, and the subscripts P and N indicate the pole of the great circle and its northernmost point, respectively. The equations of the great circle and its derivative with respect to are then ' = arctan[tan ' cos( )] (6a) N N d' tan ' sin( ) tan ' sin( ) N N N N = = : (6b) d 1 + [tan ' cos( )] 1 + tan ' N N The distance from the actual impact location ( ; ' ) to a point (; ') on the great circle is given by the orthodrome c c arc = arccos[sin ' sin ' + cos ' cos ' cos( )]; (7) c c c and should be minimal for the projection point on the great circle, i.e., its derivative with respect to the longitude should be zero: h i d' d' sin ' cos ' cos ' sin ' cos( ) + cos ' sin( ) c c c c d d d = p 1 [sin ' sin ' + cos ' cos ' cos( )] c c c h i d' d' sin ' cos ' cos ' sin ' cos( ) + cos ' sin( ) c c c c d d = = 0; (8) sin d' where is taken from Eq. 6b. The two existing solutions correspond to the maximum and the minimum distance. The desired minimum solution is found numerically, whereby the maximum distance can easily be excluded by applying the root- nding algorithm in a suciently narrow interval around , assuming that the great circle is a good enough 0 0 t. The coordinate for use in the model can then be given as the arc distance of the projection point ( ; ' ) from c c the northernmost point of the great circle as 0 0 0 s = arccos[sin ' sin ' + cos ' cos ' cos( )]: (9) N N N c c c The placement of the impacts on a great circle in our GC model series uses this relation. Hence, a given impact will have dierent positions in dierent models. B Poiseuille ow with temperature-dependent rheology In order to construct a model simple enough for an analytical solution that captures some features of the upwelling in a hot plume generated by the thermal anomaly of an impact, we assume: 1. that the ow is steady and occurs in a vertical cylindrical pipe with rigid walls at R as in the classical pipe- ow (Poiseuille ow) problem; 2. that the temperature anomaly T is only a function of distance r from the central axis of the pipe; 3. that the decisive in uence on viscosity is the temperature and that it can be described by a Frank-Kamenetskii law, i.e., (T (r)) = exp(bT (r)) with constant parameters and b; 4. that the timescale of the ow process is much shorter than the thermal diusion timescale, such that we can assume a time-independent T (r) that is zero for all r R and T at r = 0. For our max purposes, T (r) = T [1 (r=R) ] is a useful choice that ful lls these criteria. max We start with the same assumptions as in the derivation for the classical Poiseuille ow problem as given by Turcotte and Schubert (1982), i.e., by stating that the force resulting from the pressure dierence p that drives the ow and the shear force on the walls of the cylindrical control volume of length l and radius r cancel out: r dp r p = 2rl , = : (10) 2 dz 20 In our setup, is directly related to the ow velocity v by the ow law, i.e., = dv=dr. The driving pressure results from the thermal buoyancy of the anomaly, which causes a contrast in the vertical lithostatic pressure gradient of dp=dz = g%(r). The density variation is controlled by the thermal expansion due to the temperature anomaly, i.e., % = % T (r); % and are the background density and thermal expansivity, respectively, and are constant. 0 0 Setting Eq. 10 equal to the ow law and inserting the relations for the viscosity, the density and temperature anomaly, and the pressure gradient, we arrive at a dierential equation that can be solved by separation of variables: 2 2 r dv g% T r r 0 max exp bT 1 = 1 (11) 0 max R dr 2 R 2 2 g% R r r ) v(r) = 1 bT 1 exp bT 1 1 : (12) max max 4 b T R R 0 max This ow pro le resembles a Gaussian bell curve and is more focused at the axis for a given R than the classical Poiseuille ow pro le, which is parabolic. The focusing increases with the magnitude of the temperature anomaly and with increasing b. C Simpli ed anomaly volumes In order to estimate the dependence of the impact-generated thermal anomalies on their spatial and temporal sepa- ration, we need the dependence of the anomaly volumes on the separation of the two impacts, x. As we are only interested in a general order-of-magnitude estimate, we approximate the anomalies by two overlapping spheres of radius R whose centers are x apart from each other, i.e., the anomalies are centered at x = R and x = R + x and intersect at x = R + x=2; R can be any reasonable length scale of the anomaly, likely a small multiple of the isobaric core radius. Our anomaly outline is thus de ned by the outline 2 2 2 R (x R) = 2xR x ; 0 x R + x=2 a(x) = (13) 2 2 R (x R x) ; R + x=2 < x 2R + x=2 and zero elsewhere. With respect to the comparison of two- and three-dimensional models, we derive the solutions for both two- and three-dimensional anomalies. As both anomalies are assumed to be of equal size, all integrals are only calculated on the interval [0; R + x=2], and the result is doubled. We only need to consider separations up to 2R, because beyond that distance, there are simply two separate spherical anomalies. The \volume" of the anomaly in two-dimensional models such as our numerical models is the area between a(x) and a(x): R+ " # x x x V = 4 2xR x dx = R 1 + 2 arcsin + R 4R 2R " # x x 2 x = A 1 + 1 + arcsin ; (14) circ R 4R 2R where A is the area of a circle with radius R. circ In a three-dimensional real planet, the idealized anomaly is given by rotating a(x) around the x axis. The rotational body has the volume x x R+ R+ 2 2 Z Z 2 2 V = 2 a (x) dx = 2 (2xR x ) dx 0 0 " # 4 1 3 x 1 x 3 2 3 = R + R x x = V 1 + ; (15) sph 3 12 4 R 16 R where V is the volume of a sphere with radius R. Both functions are strictly monotonic on the interval of interest, sph 0 x=R 2. In their nal forms, they are cast as products of the volume of a simple sphere (in two or three dimensions) and a form factor. The form factor functions for volumes dier by at most 5% between the two- and three-dimensional cases, which suggests that eects that depend on volume follow comparable trends in their dependence on impact spacing in two and three dimensions. In this discussion, we have also assumed that the spheres are fully immersed in the target planet, i.e., they are not truncated by the surface. In the models, we have set the depth of the isobaric core according to the estimate 0:5 imp 0:361 z = 0:1524D v (16) ic imp imp (Ruedas and Breuer, 2018, eq. 6), which is an extension of an empirical relation by Pierazzo et al. (1997). Comparison of the isobaric core diameters and depths in Table 2 shows that this simpli cation is well justi ed for the impacts considered here if R is taken to be of the order of R . ic 21 D Eect of crater overlap on heat ow We consider only complex craters and neglect the curvature of the target body for simplicity. Taking the diameter of the nal crater as the measure of the area in which the heat ow is modi ed, two (identical) impacts will be independent if the distance of their centers is greater than D . For a given impact, our measure of the isobaric core is given by n 1 D = 2r = D b ; (17) ic in imp n + 1 where b and n are empirical parameters that depend on v (Ruedas, 2017). Solving for D , equating with Eq. 1, imp imp and rearranging gives 0:8814 0:377 0:4972 % n + 1 imp imp 0:8814 D = 1:3836 D ; (18) ic 0:2486 0:13 % g D b(n 1) sc 0:8814 0:8814 i.e., D D ; for our setup b = 2:463 and n = 1:424, and hence D 17:9D . Thus, all of our double-impact f f ic ic models except the Utopia-sized ones with x = 10 are expected to lie in the overlap region. It should be kept in mind, however, that D D does not imply that the entire thermal anomaly is shielded by the modi ed crust; the f ic shock-heated region is several times larger than the isobaric core, in which the highest shock pressures are reached. Neglecting curvature, we can use Eq. 14 with R being replaced by D =2 to estimate the eective total impact- modi ed area for overlapping impacts on a sphere: " # D 2x x x A = + 1 + 2 arcsin : (19) 4 D D D f f For our two-dimensional models, the modi ed \surface" is an arc of length x D , and so the eective modi ed \area" increases linearly with x=D up to 2D . f f References Aoki, I., Takahashi, E., 2004. Density of MORB eclogite in the upper mantle. Phys. Earth Planet. Inter. 143{144, 129{143. doi:10.1016/j.pepi.2003.10.007. Archinal, B.A., A'Hearn, M.F.A., Bowell, E., Conrad, A., Consolmagno, G.J., Courtin, R., Fukushima, T., Hestroer, D., Hilton, J.L., Krasinsky, G.A., Neumann, G., Oberst, J., Seidelmann, P.K., Stooke, P., Tholen, D.J., Thomas, P.C., Williams, I.P., 2011. Report of the IAU Working Group on Cartographic Coordinates and Rotational Elements: 2009. Celest. Mech. Dyn. Astr. 109, 101{135. doi:10.1007/s10569-010-9320-4. erratum in 110, 401{403. Arkani-Hamed, J., Ghods, A., 2011. Could giant impacts cripple core dynamos of small terrestrial planets? Icarus 212, 920{934. doi:10.1016/j.icarus.2011.01.020. Bertka, C.M., Fei, Y., 1997. Mineralogy of the Martian interior up to core-mantle boundary pressures. J. Geophys. Res. 102, 5251{5264. doi:10.1029/96JB03270. Bottke, W.F., Andrews-Hanna, J.C., 2017. A post-accretionary lull in large impacts on early Mars. Nature Geosci. 10, 344{348. doi:10.1038/ngeo2937. Carry, B., 2012. Density of asteroids. Planet. Space Sci. 73, 98{118. doi:10.1016/j.pss.2012.03.009. Catling, D.C., 2015. Planetary atmospheres, in: Spohn, T. (Ed.), Physics of Terrestrial Planets and Moons. 2nd ed.. Elsevier. volume 10 of Treatise on Geophysics. chapter 10.13, pp. 429{472. doi:10.1016/B978-0-444-53802-4. 00185-8. Cliord, S.M., 1993. A model for the hydrologic and climatic behavior of water on Mars. J. Geophys. Res. 98, 10973{11016. doi:10.1029/93JE00225. Faul, U.H., 2001. Melt retention and segregation beneath mid-ocean ridges. Nature 410, 920{923. doi:10.1038/ Frey, H., 2008. Ages of very large impact basins on Mars: Implications for the late heavy bombardment in the inner solar system. Geophys. Res. Lett. 35. doi:10.1029/2008GL033515. Frey, H.V., Mannoia, L.M., 2013. A revised, rated and dated inventory of very large candidate impact basins on Mars. Lunar Planet. Sci. 44. Gault, D.E., Wedekind, J.A., 1978. Experimental studies of oblique impact. Proc. Lunar Planet. Sci. Conf. 9, 3843{ 22 Gillmann, C., Golabek, G.J., Tackley, P.J., 2016. Eect of a single large impact on the coupled atmosphere-interior evolution of Venus. Icarus 268, 295{312. doi:10.1016/j.icarus.2015.12.024. Golabek, G.J., Keller, T., Gerya, T.V., Zhu, G., Tackley, P.J., Connolly, J.A.D., 2011. Origin of the martian dichotomy and Tharsis from a giant impact causing massive magmatism. Icarus 215, 346{357. doi:10.1016/j.icarus.2011. 06.012. Hernlund, J.W., Tackley, P.J., 2008. Modeling mantle convection in the spherical annulus. Phys. Earth Planet. Inter. 171, 48{54. doi:10.1016/j.pepi.2008.07.037. Hirth, G., Kohlstedt, D.L., 2003. Rheology of the upper mantle and the mantle wedge: A view from the experimen- talists, in: Eiler, J. (Ed.), Inside the Subduction Factory. American Geophysical Union, Washington, D.C.. volume 138 of AGU Geophysical Monograph, pp. 83{105. Ivanov, B.A., 2001. Mars/Moon cratering rate ratio estimates. Space Sci. Rev. 96, 87{104. doi:10.1023/A: Jacobson, R.A., 2010. The orbits and masses of the martian satellites and the libration of Phobos. Astron. J. 139, 668{679. doi:10.1088/0004-6256/139/2/668. Kendall, J.D., Melosh, H.J., 2016. Dierentiated planetesimal impacts into a terrestrial magma ocean: Fate of the iron core. Earth Planet. Sci. Lett. 448, 24{33. doi:10.1016/j.epsl.2016.05.012. Konopliv, A.S., Asmar, S.W., Folkner, W.M., Karatekin, O., Nunes, D.C., Smrekar, S.E., Yoder, C.F., Zuber, M.T., 2011. Mars high resolution gravity elds from MRO, Mars seasonal gravity, and other dynamical parameters. Icarus 211, 401{428. doi:10.1016/j.icarus.2010.10.004. Lillis, R.J., Robbins, S., Manga, M., Halekas, J.S., Frey, H.V., 2013. Time history of the martian dynamo from crater magnetic eld analysis. J. Geophys. Res. 118, 1488{1511. doi:10.1002/jgre.20105. McKenzie, D., 1985. The extraction of magma from the crust and mantle. Earth Planet. Sci. Lett. 74, 81{91. doi:10.1016/0012-821X(85)90168-2. Melosh, H.J., 2011. Planetary Surface Processes. Number 13 in Cambridge Planetary Science, Cambridge University Press. Monteux, J., Arkani-Hamed, J., 2014. Consequences of giant impacts in early Mars: Core merging and martian dynamo evolution. J. Geophys. Res. 119, 480{505. doi:10.1002/2013JE004587. Nimmo, F., Price, G.D., Brodholt, J., Gubbins, D., 2004. The in uence of potassium on core and geodynamo evolution. Geophys. J. Int. 156, 363{376. doi:10.1111/j.1365-246X.2003.02157.x. O'Neill, C., Marchi, S., Zhang, S., Bottke, W., 2017. Impact-driven subduction on the Hadean Earth. Nature Geosci. 10, 793{797. doi:10.1038/ngeo3029. Padovan, S., Tosi, N., Plesa, A.C., Ruedas, T., 2017. Impact-induced changes in source depth and volume of magmatism on Mercury and their observational signatures. Nat. Comm. 8. doi:10.1038/s41467-017-01692-0. Pierazzo, E., Vickery, A.M., Melosh, H.J., 1997. A reevaluation of impact melt production. Icarus 127, 408{423. doi:10.1006/icar.1997.5713. Reese, C.C., Solomatov, V.S., Baumgardner, J.R., 2002. Survival of impact-induced thermal anomalies in the martian mantle. J. Geophys. Res. 107. doi:10.1029/2000JE001474. Reese, C.C., Solomatov, V.S., Baumgardner, J.R., Stegman, D.R., Vezolainen, A.V., 2004. Magmatic evolution of impact-induced Martian mantle plumes and the origin of Tharsis. J. Geophys. Res. 109. doi:10.1029/2003JE002222. Rivoldini, A., Van Hoolst, T., Verhoeven, O., Mocquet, A., Dehant, V., 2011. Geodesy constraints on the interior structure and composition of Mars. Icarus 213, 451{472. doi:10.1016/j.icarus.2011.03.024. Robbins, S.J., Hynek, B.M., 2012. A new global database of Mars impact craters 1 km: 2. Global crater properties and regional variations of the simple-to-complex transition diameter. J. Geophys. Res. 117. doi:10.1029/2011JE003967. Roberts, J.H., Arkani-Hamed, J., 2012. Impact-induced mantle dynamics on Mars. Icarus 218, 278{289. doi:10.1016/ j.icarus.2011.11.038. Roberts, J.H., Arkani-Hamed, J., 2014. Impact heating and coupled core cooling and mantle dynamics on Mars. J. Geophys. Res. 119, 729{744. doi:10.1002/2013JE004603. Roberts, J.H., Barnouin, O.S., 2012. The eect of the Caloris impact on the mantle dynamics and volcanism of Mercury. J. Geophys. Res. 117. doi:10.1029/2011JE003876. 23 Roberts, J.H., Lillis, R.J., Manga, M., 2009. Giant impacts on early Mars and the cessation of the Martian dynamo. J. Geophys. Res. 114. doi:10.1029/2008JE003287. Rolf, T., Zhu, M.H., Wunnemann, K., Werner, S.C., 2017. The role of impact bombardment history in lunar evolution. Icarus 286, 138{152. doi:10.1016/j.icarus.2016.10.007. Rosenthal, A., Yaxley, G.M., Green, D.H., Hermann, J., Kov acs, I., Spandler, C., 2014. Continuous eclogite melting and variable refertilisation in upwelling heterogeneous mantle. Sci. Rep. 4. doi:10.1038/srep06099. Ruedas, T., 2017. Globally smooth approximations for shock pressure decay in impacts. Icarus 289, 22{33. doi:10. 1016/j.icarus.2017.02.008. Ruedas, T., Breuer, D., 2017. On the relative importance of thermal and chemical buoyancy in regular and impact- induced melting in a Mars-like planet. J. Geophys. Res. 122, 1554{1579. doi:10.1002/2016JE005221. Ruedas, T., Breuer, D., 2018. \Isocrater" impacts: Conditions and mantle dynamical responses for dierent impactor types. Icarus 306, 94{115. doi:10.1016/j.icarus.2018.02.005. Ruedas, T., Tackley, P.J., Solomon, S.C., 2013. Thermal and compositional evolution of the martian mantle: Eects of phase transitions and melting. Phys. Earth Planet. Inter. 216, 32{58. doi:10.1016/j.pepi.2012.12.002. Stanley, S., Elkins-Tanton, L., Zuber, M.T., Parmentier, E.M., 2008. Mars' paleomagnetic eld as the result of a single-hemisphere dynamo. Science 321, 1822{1825. doi:10.1126/science.1161119. Tackley, P.J., 1996. Eects of strongly variable viscosity on three-dimensional compressible convection in planetary mantles. J. Geophys. Res. 101, 3311{3332. doi:10.1029/95JB03211. Tackley, P.J., 2008. Modelling compressible mantle convection with large viscosity contrasts in a three-dimensional spherical shell using the yin-yang grid. Phys. Earth Planet. Inter. 171, 7{18. doi:10.1016/j.pepi.2008.08.005. Taylor, S.R., McLennan, S.M., 2009. Planetary Crusts. Cambridge University Press. Turcotte, D.L., Schubert, G., 1982. Geodynamics | Applications of Continuum Physics to Geological Problems. 2nd ed., John Wiley & Sons. Vervelidou, F., Lesur, V., Grott, M., Morschhauser, A., Lillis, R.J., 2017. Constraining the date of the martian dynamo shutdown by means of crater magnetization signatures. J. Geophys. Res. 122, 2294{2311. doi:10.1002/ 2017JE005410. W anke, H., Dreibus, G., 1994. Chemistry and accretion history of Mars. Phil. Trans. R. Soc. Lond. A 349, 285{293. doi:10.1098/rsta.1994.0132. Watters, W.A., Zuber, M.T., Hager, B.H., 2009. Thermal perturbations caused by large impacts and consequences for mantle convection. J. Geophys. Res. 114. doi:10.1029/2007JE002964. Werner, S.C., 2008. The early martian evolution { Constraints from basin formation ages. Icarus 195, 45{60. doi:10. 1016/j.icarus.2007.12.008. Wessel, P., Smith, W.H.F., Scharroo, R., Luis, J., Wobbe, F., 2013. Generic Mapping Tools: Improved version released. Eos Trans. AGU 94, 409{410. doi:10.1002/2013EO450001. Williams, J.P., Nimmo, F., 2004. Thermal evolution of the Martian core: Implications for an early dynamo. Geology 32, 97{100. doi:10.1130/G19975.1. Yeom, J., 2015. Rs01b1t0 solid sphere with suddenly applied temperature at outer surface. URL: http://exact. unl.edu/exact/contents/contents.php. Zhao, Y.H., Zimmerman, M.E., Kohlstedt, D.L., 2009. Eect of iron content on the creep behavior of olivine: 1. Anhydrous conditions. Earth Planet. Sci. Lett. 287, 229{240. doi:10.1016/j.epsl.2009.08.006.

Journal

Astrophysics – arXiv (Cornell University)

Published: Jan 27, 2019

Get 20M+ Full-Text Papers For Less Than $1.50/day. Start a 14-Day Trial for You or Your Team.

Learn More →

Dynamical effects of multiple impacts: Large impacts on a Mars-like planet

Dynamical effects of multiple impacts: Large impacts on a Mars-like planet

Get 20M+ Full-Text Papers For Less Than $1.50/day. Start a 14-Day Trial for You or Your Team.

Learn More →

Dynamical effects of multiple impacts: Large impacts on a Mars-like planet

Dynamical effects of multiple impacts: Large impacts on a Mars-like planet

References (62)

Abstract

Journal

Recommended Articles

There are no references for this article.

Our policy towards the use of cookies