A Steeply-inclined Trajectory For The Chicxulub Impact
PENA steeply-inclined trajectory for the Chicxulubimpact1234567890():,;G.S. Collins 1 , N. Patel 1, T.M. DavisonExpedition 364 Science Party*1,A.S.P. Rae1,2, J.V. Morgan1 & S.P.S. Gulick3,IODP-ICDPThe environmental severity of large impacts on Earth is influenced by their impact trajectory.Impact direction and angle to the target plane affect the volume and depth of origin ofvaporized target, as well as the trajectories of ejected material. The asteroid impact thatformed the 66 Ma Chicxulub crater had a profound and catastrophic effect on Earth’senvironment, but the impact trajectory is debated. Here we show that impact angle anddirection can be diagnosed by asymmetries in the subsurface structure of the Chicxulubcrater. Comparison of 3D numerical simulations of Chicxulub-scale impacts with geophysicalobservations suggests that the Chicxulub crater was formed by a steeply-inclined (45–60 tohorizontal) impact from the northeast; several lines of evidence rule out a low angle ( 30 )impact. A steeply-inclined impact produces a nearly symmetric distribution of ejected rockand releases more climate-changing gases per impactor mass than either a very shallow ornear-vertical impact.1 Department Earth Science and Engineering, Imperial College London, London SW7 2AZ, UK. 2 Institute of Geology, University of Freiburg, Freiburg 79104,Germany. 3 Institute for Geophysics and Department of Geological Sciences, University of Texas at Austin, Austin, TX 78758, USA. *A list of members andtheir affiliations are listed at the end of the paper. email: g.collins@imperial.ac.ukNATURE COMMUNICATIONS (2020)11:1480 https://doi.org/10.1038/s41467-020-15269-x www.nature.com/naturecommunications1
ARTICLENATURE COMMUNICATIONS st8,13. However, these large-scale, peripheral features areall likely to be pre-existing features, unrelated to the impact12.As models of the subsurface based on potential field data haveinherently poor resolution and suffer from non-uniqueness, themost robust evidence of asymmetry comes from seismic reflectionand refraction data14. High-resolution seismic reflection imagesalong a concentric arc outside the crater rim clearly show that thenortheastern gravity low in the offshore half of the crater occursin an area where Cretaceous and Cenozoic sedimentary rocks areparticularly thick, and the northwestern gravity high occurswhere this sedimentary sequence is thinnest and basement rocksare closer to surface (ref. 3, Fig. 3). Given the observed correlationbetween the gravity signature and depth to basement outside thecrater, the thickness of the sedimentary sequence is the mostlikely control on the offshore gravity anomaly. There is thereforeno evidence that the azimuthal asymmetry in the outer gravitysignature is impact related.A nominal geographical centre of the crater (21.29 N, 89.53 W)is defined by the geometric centre of the crater rim demarcated byboth a circular high in horizontal gravity gradient and the prominent cenote ring11,13 (Fig. 1). Relative to this point, the centre(21.24 N, 89.58 W) of both the central gravity high, attributed tothe uplift of dense lower-crustal rocks11, and the surroundingannular gravity low, which underlies the inner edge of the peakring, are shifted several km to the west-southwest (Fig. 1 andSupplementary Fig. 1). In contrast, three-dimensional (3D) seismiche 66 Ma asteroid impact event that formed the Chicxulubcrater, Mexico, marks the end of the Mesozoic Era of Earthhistory and has been attributed as the cause of the contemporaneous mass extinction1. Numerical impact simulationscombined with geophysical investigation of subsurface structurehave constrained the kinetic energy of the impact under thesimplifying assumption of a vertical trajectory2–4. The trajectoryangle and direction of the Chicxulub impact are not known, but anear-vertical impact is unlikely. Only one quarter of impactsoccur at angles between 60 and the vertical and only 1 in 15impacts is steeper than 75 .For constant impactor size and speed, a shallower impact angleproduces a smaller crater, a more asymmetric dispersal of ejectedmaterial5 and partitions more impact energy at shallower depths6.As a result, impact direction and trajectory angle to the targetplane are important impact parameters that determine, amongother things, the direction of most severe environmentalconsequences and the volume and depth of origin of vaporisedtarget7,8, as well as ejecta5 and crater asymmetries9.Since the discovery of the Chicxulub impact structure based ondiagnostic evidence of shock metamorphism and geophysicalanomalies10, several asymmetries in the geophysical character ofthe crater have been noted8,11,12, which may result from obliqueimpact8,11 and/or impact in a heterogeneous target3,12. Amongthe most obvious of these (Fig. 1) are radially oriented gravitylows to the south and northeast, and a radial gravity high to theGravity (mGal)–16–12–8–404812162024283236404448N22 Mantle upliftcenterInner rimftoten ngltifauExM0077APeak ringcenterCratercenterMérida21 35 15 –100 –90 –80 –89 Fig. 1 Asymmetries of the geophysical signature of the Chicxulub crater. Background colourmap shows Bouguer gravity anomaly map in the vicinity ofthe crater (gravity data courtesy of Hildebrand and Pilkington). The red circle marks the nominal position of the crater centre; the green circle marks thecentre of maximum mantle uplift; the blue circle marks the centre of the peak ring (as defined by the annular gravity low surrounding the central high); thewhite triangle marks the location of the Expedition 364 drill site through the peak ring (Hole M0077A). The coastline is displayed with a thin white line;cenotes and sinkholes with white dots, and the city of Mérida with a white square. The dotted lines offshore mark the approximate location of the innercrater rim and the extent of faulting as imaged by seismic data14. Inset depicts the regional setting, with red rectangle outlining the region shown in thegravity map. Adapted from ref. 14.2NATURE COMMUNICATIONS (2020)11:1480 https://doi.org/10.1038/s41467-020-15269-x www.nature.com/naturecommunications
NATURE COMMUNICATIONS https://doi.org/10.1038/s41467-020-15269-xvelocity data indicate that the maximum uplift of the mantlebeneath the crater occurs at 21.38 N, 89.52 W, ca. 10 km to thenorth-northeast of the crater centre4 (Fig. 1, Supplementary Fig. 1and ref. 4, Fig. 4d).The south-westerly offset of the central gravity high relative tothe crater centre was previously interpreted as indicating impactfrom the southwest, on the premise that central uplift motionwould be directed uprange11. An alternative interpretation, of atrajectory from the southeast, was proposed on the basis of anorthwest–southeast elongation of the central gravity high andmagnetic anomaly, and the northwest truncation of the annulargravity low8. However, seismic reflection and refraction datareveal that the zone of structural uplift is not elongated towardsthe northwest14 (see also Supplementary Fig. 1) and that thetruncation of the gravity low in the northwest is a pre-impactfeature of the regional anomaly caused by the shallow depth tobasement in this direction. The short-wavelength component ofthe magnetic anomaly shows a slight (10%) elongation in thenorthwest–southeast direction15,16, but is also offset to thesouthwest of the crater centre (Supplementary Fig. 1). The shortwavelength and steep gradients of this anomaly both suggest ashallow source, probably related to the melt sheet and impactbreccia and not to structural crater asymmetry15,16. On the otherhand, the long-wavelength component of the magnetic anomalyis a magnetic high elongated and offset along a direction southwest of the crater centre (Supplementary Fig. 1), consistent with azone of uplifted basement rocks southwest of the crater centre15.Here we use 3D numerical modelling to examine the relationship between impact angle and structural crater asymmetries in aChicxulub-scale peak-ring crater in a flat-layered target withoutlateral pre-impact asymmetry. We show that the observed asymmetry in the positions of the central uplift, peak-ring centre andmaximum mantle uplift, relative to the crater centre, can beattributed to the angle and azimuth of the impact trajectory.Comparison of our simulation results with geophysically constrained models of the Chicxulub crater structure is used to inferthe likely trajectory and angle of the impact. The recent jointInternational Ocean Discovery Program (IODP) and InternationalContinental Scientific Drilling Program (ICDP) Expedition 364recovered 600 m of peak-ring rocks from the Chicxulub crater17that provide additional constraints to discriminate between impactscenarios. Our simulations also reveal azimuthal variation in peakring material properties, which provide context for IODP-ICDPExpedition 364 core analysis.Results and discussionNumerical simulation results. We performed a series of 3Dsimulations of impacts that produce a Chicxulub-scale crater,using the iSALE3D shock physics code18,19. The simulationsassumed a flat, two-layer target comprising crust and mantleand considered four different impact angles (90 (vertical), 60 ,45 and 30 ) and two impact speeds (12 and 20 km/s). Furthermodel details are described in the “Methods” section. Oursimulations provide insight into crater asymmetries diagnostic ofimpact angle and trajectory in the absence of any target asymmetry (Figs. 2 and 3, Supplementary Figs. 3–5 and SupplementaryMovies 1–4).In our vertical impact simulation (Supplementary Fig. 3),crater formation is axially symmetric and consistent withprevious two-dimensional (2D) numerical simulations thatemployed an axially symmetric geometry2,3,17. Collision of theasteroid with the target surface generates a detached shockwavethat propagates symmetrically from the impact site. In the firstminute after impact, an excavation flow initiated by theshockwave produces a deep, bowl-shaped cavity, often termedARTICLEthe transient crater. The material flow depresses the crustmantle boundary beneath the transient crater, uplifts the crustin the transient crater wall and expels the unvaporized portionof the 3-km-thick sedimentary rock sequence from thetransient crater as part of the ejecta curtain (SupplementaryFig. 3 and Supplementary Movie 1).The transient crater is unstable and collapses dramatically toproduce a much flatter, broader final crater. In the vertical impactsimulation, collapse manifests as uplift of the crater floor anddownward and inward collapse of the transient crater rim and asurrounding collar of sedimentary rocks. Floor uplift beginsdirectly beneath the transient crater centre and proceeds verticallyupward, overshooting the pre-impact surface to form a largecentral uplift. At the same time, rim collapse occurs symmetrically at all azimuths, converging towards, and helping to driveup, the central uplift. Finally, the overheightened central uplift ofcrustal rocks collapses downward and outward, overthrusting thecollapsed transient crater rim to form an uplifted ring ofcrystalline basement, overlying inwardly slumped sedimentaryrocks from outside the transient crater. Although the spatialresolution of the numerical simulations is insufficient to resolvethe characteristic sharp-peaked topography of the inner ringobserved in extraterrestrial peak-ring craters, we are able toidentify the position and structure of the material that forms thepeak ring in the numerical simulations as a 10-km-wide collararound the central uplift (Supplementary Fig. 2). This modelof peak-ring crater formation is supported by geophysicaldata20,21 and recent geological drilling17 at Chicxulub, as wellas remote-sensing data from the Schrödinger peak-ring crater onthe Moon22.Impacts at progressively shallower angles to the horizontalresult in an increasingly asymmetric development of the crater,internally (Figs. 2 and 3, Supplementary Fig. 4; SupplementaryMovies 2–4), while the planform of the final impact basin remainsapproximately circular (Supplementary Table 1; SupplementaryFig. 9). As impact angle decreases, the downrange offset of thecrater centre from the impact point increases; less uplift of thetransient crater rim occurs in the uprange direction; and moreuplift occurs in the downrange direction (Figs. 2 and 3,Supplementary Fig. 4). Relative to the final crater centre, thedeepest part of the transient crater (and depressed mantle) alsoshifts with decreasing impact angle, first to the uprange direction(Fig. 2), then back towards the centre (Supplementary Fig. 4 andFig. 3). The collapse phase of crater formation is also modified byimpact angle. Uplift of the crater floor during crater collapsebegins uprange of the crater centre, but has a downrangecomponent such that the central uplift axis is tilted downrangeand the centre of the uplift prior to its collapse is downrangeof the crater centre (Figs. 2 and 3, Supplementary Fig. 4;Supplementary Movies 2–4). Conversely, downward and outwardcollapse of the central uplift occurs preferentially in the uprangedirection, resulting in enhanced overthrusting of the central uplifton top of transient crater rim in the uprange direction. The netresult of the downrange-directed rise and uprange-directed fall ofthe central uplift is a simulated peak ring with a geometric centreonly modestly offset in the downrange direction (Figs. 4 and 5).The impact simulations shown in Figs. 2 and 3 employ animpact speed of 12 km/s, only slightly larger than the minimumpossible speed—Earth’s escape velocity of 11.2 km/s. While theseresults are likely to be representative of the 25% of all impactsthat occur at speeds below 15 km/s, we also conducted anothersuite of simulations with a more probable impact speed of 20 km/s (close to Earth’s mean and median asteroid impact speed23) toexamine the sensitivity of our results to impactor speed. Thehigher-speed impacts produced similar offsets in mantle-upliftcentre and simulated peak-ring centre (compare Figs. 2 and 3NATURE COMMUNICATIONS (2020)11:1480 https://doi.org/10.1038/s41467-020-15269-x www.nature.com/naturecommunications3
ARTICLENATURE COMMUNICATIONS https://doi.org/10.1038/s41467-020-15269-xa20T 0s010102030405060Tracer peak pressure [GPa]z 0602080T 20 s10z 0206080T 180 s10z 0206080T 300 s10z [km]0–10–20–30–40–80–60–40–200x [km]20406080Fig. 2 Development of the Chicxulub crater for a 60 impact. The impact scenario depicted is for a 17-km diameter impactor with a density of2630 kg m3 and a speed of 12 km/s. Evolution of the crater up to 5 min after impact is depicted. Shown are cross-sections through the numericalsimulation along the plane of trajectory, with x ¼ 0 defined at the crater centre (measured at the pre-impact level; z ¼ 0); the direction of impact isfrom right to left. The upper 3 km of the pre-impact target, corresponding to the average thickness of sedimentary rocks at Chicxulub, is tracked bytracer particles (sandy brown). Deformation in the crust (mid-grey) and upper mantle (dark grey) is depicted by a grid of tracer particles (black).Tracer particles within the peak-ring material are highlighted based on the peak shock pressure recorded (white–blue colour scale); melted targetmaterial ( 60 GPa) is highlighted in red.4NATURE COMMUNICATIONS (2020)11:1480 https://doi.org/10.1038/s41467-020-15269-x www.nature.com/naturecommunications
ARTICLENATURE COMMUNICATIONS https://doi.org/10.1038/s41467-020-15269-xa20T 0s010102030405060Tracer peak pressure [GPa]z 0602080T 20 s10z 0206080T 180 s10z 0206080T 300 s10z [km]0–10–20–30–40–80–60–40–200x [km]20406080Fig. 3 Development of the Chicxulub crater for a 30 impact. The impact scenario depicted is for a 21-km diameter impactor with a density of 2630 kg m3and a speed of 12 km/s. Evolution of the crater up to 5 min after impact is depicted. Shown are cross-sections through the numerical simulation along theplane of trajectory, with x ¼ 0 defined at the crater centre (measured at the pre-impact level); the direction of impact is from right to left. Colours andshading of material and tracer particles are the same as Fig. 2.with Supplementary Figs. 6 and 7), and the same trends in offsetswith impact angle (Fig. 5).An important consequence of higher impact speed isenhanced melt production caused by higher shock pressuresclose to the impact site (e.g., compare Fig. 2a and SupplementaryFig. 6a). The larger melt volume complicates the interpretationof peak-ring structure in the 20 km/s simulations as thedynamics of the melt are not expected to be well captured,given the 500-m spatial resolution of the 3D simulations, andwould likely continue long after the simulation end time.NATURE COMMUNICATIONS (2020)11:1480 https://doi.org/10.1038/s41467-020-15269-x www.nature.com/naturecommunications5
ARTICLENATURE COMMUNICATIONS 05060Tracer peak pressure [GPa]a2010z –60–40–2002040608020406080b2010z [km]0–10–20–30–400x [km]Fig. 4 Final simulated Chicxulub crater for a 30 and 60 impact angle. Shown are cross-sections, along the plane of trajectory, through the finalsimulated craters formed by a 30 (a) and 60 (b) impact, with x ¼ 0 defined at the crater centre (measured at the pre-impact level); the direction ofimpact is from right to left. Sandy-brown tracers indicate the final position of the upper 3 km of the pre-impact target (sedimentary rock); red tracersindicate the position of melt; tracers with blue–white shading indicate shock pressures of simulated peak-ring materials. The geometric centre of the craterrim defines the coordinate origin (x ¼ 0); negative x-values are downrange.Offset from crater center/Dc0.100.05Obs. mantle uplift0.00Obs. peak ring–0.05Peak ring center (12 km/s)Mantle uplift center (12 km/s)Peak ring center (20 km/s)–0.10Mantle uplift center (20 km/s)3040507060Impact angle [ ]8090Fig. 5 Offset of structural crater features relative to crater centre. Theoffsets of the centre of mantle uplift (squares) and the centre of thesimulated peak ring (circles), relative to the crater centre, are shown as afunction of impact angle to the horizontal. Grey bands denote theapproximate relative offsets of the peak-ring and mantle-uplift centres atChicxulub, taking into account the uncertainty in crater diameter andlocations of the different features (see Supplementary Fig. 1).6Nevertheless, the lateral distribution of the melt material relativeto the peak-ring material at the end of the simulations(Supplementary Figs. 8 and 10) suggests that below an impactangle of 45 there is a high concentration (thick sheet) ofsurficial melt in the downrange quadrant of the crater, which islikely to hinder or prevent formation of a topographic peak ringat these azimuths. Our results therefore support the idea thathorse-shoe shaped peak-ring planforms are indicative ofshallow-angle impacts, with the gap in the peak ring diagnosticof the downrange direction24.Comparison with observations. Asymmetry in crater development produces differences in central crater structure in theuprange and downrange directions. While the centre of thesimulated peak ring appears to be consistently offset downrange ofthe crater centre by 5% of the crater diameter in the three obliqueimpacts, the centre of the mantle uplift is offset uprange of thecrater centre in the 60 impact and, to a lesser extent, the 45 impact; and is offset downrange in the 30 impact (Fig. 5). Thispattern of mantle-uplift offset relative to the final crater rim is aconsequence of the corresponding change in the offset of thedeepest part of the transient crater relative to the centre of the finalcrater. Geophysical observations at Chicxulub suggest the peakring and mantle-uplift centres are offset in different, approximately opposite directions from the crater centre (Fig. 1).Uncertainty in the precise locations of the centres of the crater,NATURE COMMUNICATIONS (2020)11:1480 https://doi.org/10.1038/s41467-020-15269-x www.nature.com/naturecommunications
NATURE COMMUNICATIONS https://doi.org/10.1038/s41467-020-15269-xpeak ring and mantle uplift (Supplementary Fig. 1), as well asuncertainty in the crater diameter25, contribute to an approximateuncertainty of 26% and 48% for the relative offset of the peak ringand mantle uplift, respectively (grey bands in Fig. 5). Comparisonof these observations with our simulation results suggests that theobservedconfigurationismostsimilartothe60 impact simulations (or possibly the 45 impact simulation at20 km/s; Fig. 5).Tracer particles that track the history of material in thesimulation afford analysis of the provenance of peak-ringmaterials and their variation with azimuth. The mean depth oforigin of peak-ring materials is 10–12 km for the 45 , 60 and 90 impacts, only dropping significantly, to 8 km, in the 30 impact.In the 30 scenario a significant fraction of the simulated peakring originates from the sedimentary sequence in the uprangedirection (Fig. 4); the presence of significant amounts ofsedimentary material in the simulated peak ring is not consistentwith geophysical interpretations or results from Expedition36417,26,27.We also observe a systematic change in the up/downrangedifference in subsurface structure of simulated peak rings withimpact angle (Fig. 4). Similar to the situation in a vertical impact,at 60 the simulated peak ring is formed of overthrusted graniticcrustal rocks from the central uplift above down-slumpedsedimentary rocks from the transient crater wall, in all directions.However, the sedimentary rocks are deeper and extend fartherbeneath the simulated peak ring in the uprange directioncompared with the downrange direction (Fig. 4). At 45 and30 this difference is more pronounced: on the downrange side ofthe crater, the inwardly slumped sedimentary rocks do not extendunder the simulated peak ring (Fig. 4) owing to enhancedtransient crater rim uplift in this direction. This downrangeconfiguration is inconsistent with geophysical interpretations atChicxulub, which suggest sedimentary slump blocks lie beneaththe outer portion of the peak ring at all azimuths offshore12,25.However, pre-impact asymmetries in sediment thickness, waterdepth, particularly in the northeast part of the crater (andpotentially in the crust), may also affect structure beneath thepeak ring3.A proposed indicator of shallow-angle impact is thetruncation of the peak ring in the downrange direction24. Ournumerical simulations at typical terrestrial impact speeds(20 km/s) are consistent with the production of a gap in thepeak ring in the downrange direction for impact anglesshallower than 45 (Supplementary Figs. 8 and 10). However,a prominent gap in the Chicxulub peak ring that might indicatea shallow-angle impact is not supported by the geophysical data.The topographic expression of the peak ring is clearly resolvedin all radial seismic reflection lines through the offshore portionof the crater28 and is particularly prominent in the northwestseismic reflection line Chicx-B28, the downrange directionaccording to shallow-angle impact hypothesis proposed bySchultz and d’Hondt8. While the onshore portion of the craterhas not been seismically imaged, the annular negative gravityanomaly that has been shown to correlate with peak-ringposition offshore is well-pronounced and continuous in thisregion, with no break that might indicate an abundance of meltor change in the character of the peak ring. The continuity of thegeophysical signature of the peak ring therefore also supports amore steeply inclined impact trajectory.In summary, our numerical simulations of oblique Chicxulubscale impacts appear to be most consistent with the internalstructure of the Chicxulub crater for a steeply inclined impactangle of 45–60 to the horizontal. If the observed asymmetries inthe Moho uplift, central uplift and peak ring of the Chicxulubimpact structure are attributable to impact trajectory, the impliedARTICLEdirection of impact is northeast-to-southwest. This is the oppositedirection to that proposed by Hildebrand et al.11 based on theoffset of the central uplift relative to the crater centre. Our resultsindicate that uplift of the crater floor occurs in a downrangerather than uprange direction, consistent with numericalsimulations of complex crater formation19 and geologicalinterpretation of eroded structural uplifts at terrestrial complexcraters9,29.Analyses of Venusian craters have not shown a clear linkbetween asymmetries in central crater features and direction ofimpact. A slight tendency for the peak-ring centre to be offset inthe downrange direction was observed, but the results wereinconclusive, in part owing to the relatively small numberof craters used in the study30. The magnitude of the offset(0.03–0.07 D) is, however, consistent with our numericalsimulation results. In contrast, there is no correlation betweenimpact trajectory direction and the offset from the crater centreof central peaks in small complex craters31. While we did notsimulate central peak formation in this work, our resultsprovide a possible explanation for the absence of correlation. Atsteep angles, uplift of the crater floor initiates uprange of thecrater centre, while at shallow angles uplift initiates downrange.If central peaks represent frozen central uplifts, offsets in eitheruprange or downrange direction might therefore be expected atmoderately oblique angles 30–60 .Implications of a steeply inclined Chicxulub impact. Impactsthat occur at a steep angle of incidence are more efficient atexcavating material and driving open a large cavity in the crustthan shallow incidence impacts5,19. Our preferred impact angle ofca. 60 is close to the most efficient, vertical scenario, whichsuggests that previous estimates of impactor kinetic energy basedon high-resolution 2D vertical impact simulations2,17 do not needto be revised dramatically based on impact angle.Steeply inclined impacts favour a more symmetric distributionof material ejected from the crater among both proximal anddistal ejecta5. Asymmetry in the distribution of ejecta wasoriginally used by Schultz and d’Hondt8 as an argument for ashallow impact angle towards the northwest. This was based onthe observation that both the particle size and layer thicknesswere relatively large in North American K–Pg sites. Subsequentwork has shown that number and size of shocked quartz grainspresent in the global ejecta layer decreases with distance fromChicxulub, and is independent of azimuth32–34. In addition, the1–3-cm-thick double layer in North America is also observed tothe south and southeast of Chicxulub in Colombia35 and theDemerara Rise36 at equivalent paleodistances from Chicxulub.The global K–Pg boundary layer therefore has a more-or-lesssymmetric ejecta distribution, consistent with our preferred steepimpact angle.Impact angle has an important influence on the mass ofsedimentary target rocks vaporised by the Chicxulub impact37.Recent complementary numerical simulations of impact vapourproduction in oblique impacts using the SOVA shock physicscode showed that a trajectory angle of 30–60 constitutes theworst-case scenario for the high-speed ejection of CO2 andsulfur by the Chicxulub impact37. At this range of impactangles, the ejected mass of CO2 is a factor of two-to-three timesgreater than in a vertical impact and approximately an order ofmagnitude greater than a very shallow-angle (15 ) scenario37.An absence of evaporites in the IODP-ICDP Expedition 364drill core is consistent with highly efficient vaporisation ofsedimentary rocks at Chicxulub27. Our simulations thereforesuggest that the Chicxulub impact produced a near-symmetricdistribution of ejecta and was among the worst-case scenariosNATURE COMMUNICATIONS (2020)11:1480 https://doi.org/10.1038/s41467-020-15269-x www.nature.com/naturecommunications7
ARTICLENATURE COMMUNICATIONS https://doi.org/10.1038/s41467-020-15269-xfor the lethality of the impact by the production of climatechanging gases.MethodsNumerical simulations. The Chicxulub impact was simulated using the iSALE3Dshock physics code18,19. Tabular equations of state generated using ANEOS38 withinput parameters for dunite39 and granite40 were used to describe the thermodynamic response of the mantle and crust, respectively. The impactor was alsomodelled as a granite sphere, with a density of 2650 kg/m3, because of a currentlimitation of iSALE3D that does not allow for more than one boundary betweenmaterials per grid cell. The actual Chicxulub impactor density is not known.Although a carbonaceous chondrite composition has been proposed41,42, the bulkporosity of the Chicxulub asteroid prior to impact is undetermined. The Murchison(CM2) carbonaceous chondrite meteorite has a bulk density only 10% less than ourimpactor density, suggesting that our assumed impactor density is reasonable.Moreover, for a given impactor mass our simulation results are not expected to besensitive to the assumed impactor density (or other impactor material properties).Material strength was modelled using an approach appropriate for geologicalmaterials43. The choice of model parameters was based on previous vertical impactsimulations using iSALE2D3,4,17 and the similar SALEB code44,45 and obliqueimpact simulations of the early stages of the Chicxulub impact2,46.A flat, two-layer target was employed, with a crustal thickness of 33 km.Material number limitations precluded inclusion of a rheologically distinctsedimentary layer in the target; however, Lagrangian tracer particles allowedmaterial at this strat
Impact direction and angle to the target plane affect the volume and depth of origin of vaporized target, as well as the trajectories of ejected material. The asteroid impact that formed the 66 Ma Chicxulub crater had a profound and catastrophic effect on Earth’s environment, but the i
Bruksanvisning för bilstereo . Bruksanvisning for bilstereo . Instrukcja obsługi samochodowego odtwarzacza stereo . Operating Instructions for Car Stereo . 610-104 . SV . Bruksanvisning i original
(like hiking and dining) or different transportation modes, such as walking and driving. We show examples of trajectory classification in Section 7. Trajectory Outlier Detection: Different from trajectory patterns that frequently occur in trajectory data, trajectory ou
10 tips och tricks för att lyckas med ert sap-projekt 20 SAPSANYTT 2/2015 De flesta projektledare känner säkert till Cobb’s paradox. Martin Cobb verkade som CIO för sekretariatet för Treasury Board of Canada 1995 då han ställde frågan
service i Norge och Finland drivs inom ramen för ett enskilt företag (NRK. 1 och Yleisradio), fin ns det i Sverige tre: Ett för tv (Sveriges Television , SVT ), ett för radio (Sveriges Radio , SR ) och ett för utbildnings program (Sveriges Utbildningsradio, UR, vilket till följd av sin begränsade storlek inte återfinns bland de 25 största
Hotell För hotell anges de tre klasserna A/B, C och D. Det betyder att den "normala" standarden C är acceptabel men att motiven för en högre standard är starka. Ljudklass C motsvarar de tidigare normkraven för hotell, ljudklass A/B motsvarar kraven för moderna hotell med hög standard och ljudklass D kan användas vid
LÄS NOGGRANT FÖLJANDE VILLKOR FÖR APPLE DEVELOPER PROGRAM LICENCE . Apple Developer Program License Agreement Syfte Du vill använda Apple-mjukvara (enligt definitionen nedan) för att utveckla en eller flera Applikationer (enligt definitionen nedan) för Apple-märkta produkter. . Applikationer som utvecklas för iOS-produkter, Apple .
March 10, 2014 The task of determining the net force acting upon an . in your Hilroy scribbler. P12 Inclined Planes.notebook 20 March 10, 2014. P12 Inclined Planes.notebook 21 March 10, 2014 Friction and Inclined Planes The force of friction must e
Simple machines have few or no moving parts and use energy to make work easier. There are the six types of simple machines: Wedge Inclined plane Lever Screw Pulley Wheel & axle In this builder the focus will be on the wedge, inclined plane, and lever. Three Simple Machines Inclined plane An inclined plane is a flat .
