Statistics on Venus: Craters and Catastrophes (?) Steven A. Hauck, II

Transcription

1 Statistics on Venus: Craters and Catastrophes (?) Steven A. Hauck, II Department of Terrestrial Magnetism Carnegie Institution of Washington

2 Acknowledgements Roger Phillips Washington University Maribeth Price South Dakota School of Mines and Technology Sean Solomon Carnegie Institution of Washington

3 The big question What does it mean for the evolution of a planet if the spatial distribution of impact craters on its surface cannot be distinguished from a completely spatially random distribution?

4 Outline Why Venus? Why impact craters? Dating with craters. Geology in brief. Monte Carlo models and statistical tests. Implications for Venus.

5 The Basics 2 nd planet from Sun Mean radius = 6052 km ( = 6371 km) Mean density = 5243 kg/m 3 ( = 5515 kg/m 3 ) 1 Venus year = 225 days 1 Venus day = 243 days (retrograde) Surface pressure = 91 atmospheres Surface temperature = 740 K

6

7

8 Magellan SAR Mosaic of Venus

9

10 Why study impact craters?

11

12 50 km

13 Motivation Can we learn something about the history of Venus from the distribution of impact craters on the surface?

14 Relevance Surface history places a constraint on the evolution of the whole planet. Ultimately provides a contrast to the Earth which is comparable in size and presumably composition.

15 Venusian Impact Craters

16 Martian Impact Craters Craters > 4km from Barlow database over Mars shaded relief from MOLA

17 Terrestrial Impact Craters Space Imagery Center:

18 Craters Surface Ages 1) Assume the rate of impact crater formation is approximately constant (only to first-order) The rate has an impact size-dependence 2) Assume that cratering process is spatially and temporally random 3) Divide the surface into units based upon geologic criteria (e.g., morphology, superposition relationships) 4) Calculate area density of points (craters) within units 5) Relative differences give relative ages Convert to absolute age if an estimate of mean surface age is available

19 Absolute Ages Calibration points: Earth and moon Other bodies? Assumption that Mars, Venus, and Mercury have some multiple of the lunar impactor population Comparison of present day minor planets with (asteroids) known oribital elements with planetary orbits Uncertainty abounds Venus has the additional problem of its thick atmosphere

20 Crater Ages Production Age: Number of craters superposed on a geologic unit reflect the time since the unit was emplaced. Retention Age: Number of craters within a geologic unit reflect a competition between crater emplacement and removal.

21 More Background ~1000 impact craters on the surface Early analysis showed that the spatial distribution of impact craters cannot be distinguished from one that is completely spatially random [CSR] Most craters appear pristine. Dense atmosphere has a profound filtering effect Surface mean crater production age ~750 Myr Refs: Phillips et al., 1992; Schaber et al., 1992; Herrick and Phillips, 1994; McKinnon et al., 1997

22 Venusian Impact Craters

23 The big question What does it mean for the evolution of a planet if the spatial distribution of impact craters on its surface cannot be distinguished from a completely spatially random distribution?

24 Early Models Based on the notion that Venus impact craters are randomly distributed, two end-member models were proposed : The equilibrium resurfacing model (ERM) [Phillips et al., 1992] The catastrophic resurfacing model (CRM CRM) [Phillips et al., 1992; Schaber et al., 1992; Bullock et al., 1993; Strom et al., 1994]

25 Large-scale Geology Distinct morphologic units can be defined at the 1:8,000,000 scale (C1-MIDR). [Price and Suppe, 1994, 1995; Tanaka et al., 1997] The volcanic plains are the areally most extensive unit covering ~65% of the planetary surface. Plains can be divided into sub-units based upon dominant flow morphology and radar brightness. [Price, 1995; Tanaka et al., 1997]

26 SAR Images of Type PL1 and PS PL1 PS 100 km 200 km

27 Venusian Plains Units N 0 0 PL1 PL2 PL3 PS Impact Crater 45 S Plains units after Price [1995]

28 Age of the Plains Unit Area Craters Relative Age Estimated Age (Ma) PL ± 0.36 T 589 ± 270 PL ± 0.14 T 649 ± 106 PL ± 0.16 T 857 ± 116 PS ± 0.16 T 983 ± 120 PL1+PL ± 0.13 T 641 ± 99 PL3+PS ± 0.11 T 923 ± 84 SAP ± 0.09 T 829 ± 65 A unit of area is 10 6 km 2. Errors listed are 2σ. Note that both PL2 and PL1+PL2 have relative ages that do not overlap within 2σ of the single-age plains (SAP) model, suggesting that the younger plains have distinct ages that are statistically significant. The mean surface production age, used to calculate the last column, is estimated as T = 750 Ma [McKinnon et al., 1997].

29 Modeling > 200 Monte Carlo simulations Density of craters within a unit prescribed Modeling done with ArcView GIS Results post-processed to measure distances to all neighbors Mean distances to nearest neighbors compared to Venus observations using Mth nearest neighbor analysis.

30 Resurfacing Models Nominal MB1 MB2 SAP DAP DAP2 TAP Each unit has the observed age PL3-2σ, PS+2σ PL2-2σ, PS+1.5σ Single age for all plains units Combine young and old units as PL1+PL2 and PL3+PS DAP young + 2σ, DAP old - 2σ Divide units as PL1, PL2, and PL3+PS

31 Model

32 QQ and PP Plots QQ Plot of Venus Resurfacing Model with UnitsofDistinctCraterProductionAges PP Plot of Venus Resurfacing Model with Units of Distinct Crater Production Ages Observed [Deg] ExpectedDistance [Deg] Observed Distance [Deg/Deg] Expected Distance [Deg/Deg]

33 Tests Distance based Nearest Neighbor Analysis (and M th Nearest Neighbor ) compare mean distance from each crater to the 1 st,2 nd,,m th nearest neighbor to the expected distance. Density based Binomial probability probability of finding the number of craters that are observed in each unit if the hypothesis that distribution of craters in the plains is controlled only by a single random process is true. Chi-squared goodness-of-fit test compare the observed number of craters in each plains unit to the number expected by a particular model.

34 Two-sided p values of Testing the Hypothesis that Plains Resurfacing Models Represent Venus 1.0 p value Nominal MB 1 MB 2 SAP DAP DAP 2 TAP CSR CSR vs. Random M th Nearest Neighbor

35 Statistical Results Binomial Probability Unit P PL1 1.7 x 10-1 PL2 3.0 x 10-2 PL3 1.3 x 10-2 PS 4.5 x 10-7 PL1+PL2 1.4 x 10-2 PL3+PS 6.5 x SAP 4.2 x 10-6 Chi-squared goodness-of-fit Model P SAP 2.0 x 10-4 DAP 4.9 x 10-1 TAP 5.2 x 10-1 CRM 4.0 x 10-6

36 Results M th Nearest Neighbor Analysis None of the models presented (including a CSR population) can be distinguished from Venus crater distribution. Binomial probability The hypothesis that variations in the crater distribution are due to a single random process for the planet can be rejected for all units except PL1. Chi-squared goodness-of-fit test It is extremely unlikely that a SAP or CRM could result in the observed number of craters in each plains unit. Dual- or tri-age plains models cannot be rejected.

37 Conclusions CSR cannot be used as a constraint on resurfacing or geodynamic models because it is a non-unique interpretation of the crater distribution. None of the resurfacing models can be rejected as being representative of Venus based upon M th nearest neighbor analysis. Chi-squared test on crater populations within the plains units suggests that both the single-age plains and single-age planet (CSR) models can be rejected as being representative of Venus.

38 Conclusions II Binomial probability tests on plains crater populations suggest that the sub-unit ages are significant. The spread in plains ages on the order of one-half the mean production age of the surface is significant and suggests that Venus has been geologically active more recently than believed in the past. Hypotheses such as CRM and episodic resurfacing [Turcotte, 1993;1995] are unnecessary to explain the crater distribution of Venus.

39 50 km