Persistent or repeated surface habitability on Mars
during the Late Hesperian - Amazonian
EdwinS.Kite1,*,JonathanSneed1,DavidP.Mayer1,SharonA.Wilson2
1.UniversityofChicago(*[email protected])
2.SmithsonianInstitution.
Abstract.
LargealluvialfandepositsonMarsrecordrelativelyrecenthabitablesurfaceconditions
(≲3.5Ga,LateHesperian–Amazonian).Wefindnetsedimentationrate<(4-8)μm/yrinthe
alluvial-fandeposits,usingthefrequencyofcratersthatareinterbeddedwithalluvial-fan
depositsasafluvial-processchronometer.Consideringonlytheobservedinterbedded
craterssetsalowerboundof>20Myronthetotaltimeintervalspannedbyalluvial-fan
aggradation,>103-foldlongerthanpreviouslowerlimits.Amorerealisticapproachthat
correctsforcratersfullyentombedinthefandepositsraisesthelowerboundto>(100-300)
Myr.Severalfactorsnotincludedinourcalculationswouldfurtherincreasethelowerbound.
Thelowerboundrulesoutfan-formationbyabriefclimateanomaly.Therefore,duringthe
LateHesperian–AmazonianonMars,persistentorrepeatedprocessespermittedhabitable
surfaceconditions.
1.Introduction.
LargealluvialfansonMarsrecordoneormoreriver-supportingclimateson≲3.5GaMars
(e.g.Moore&Howard2005,Grant&Wilson2012,Kiteetal.2015,Williamsetal.2013)(Fig.
1).Thisclimatepermittedprecipitation-sourcedrunoffproductionof>0.1mm/hrthatfed
riverswithdischargeupto102m3/s(Morganetal.2014).Suchlargewaterdischarges
probablyrequireawateractivitythatwouldpermitlife.Thelarge(>10km2)alluvialfans
with≲2°slopescorrespondtoarelativelyrecent(Grant&Wilson2011)epochofMars
surfacehabitability(Williamsetal.2013).(Althoughsmall(<10km2)alluvialfanswith>10°
slopesformed<5MyaonMars(e.g.Williams&Malin2008),inthispaperwefocusonthe
largealluvialfans).Didtheselargefansresultfromasingleanomalousburstofwet
conditions,suchasmightresultfromanimpactorvolcaniceruption?Or,dothefansrecord
persistentorrepeatedwetconditions,forexampleastheresultofasustainedwarmer
climateregime?Betterconstraintsonthetimespanofalluvialfanformationwould
constrainmodelsofLateHesperian–Amazonianclimate.
Previousworkonthedurationoftheintervaloffanbuild-upusedsedimentologytoestimate
thetimeoverwhichsedimenttransportoccurred(e.g.,Armitage2011,Williamsetal.2011,
Palucisetal.2014).Thesesedimentologicmethodsrequireassumptionsaboutflow
intermittencyorsediment:waterratio,which(foralmostallMarsalluvialfans)arepoorly
constrained.Therefore,thelowerlimitsobtainedfromsedimentologicalmethodsareshort–
forexample,~3600years(Morganetal.2014).Moreover,brief(1-100yr)aggradation
intervalshavebeenproposedfordeltasonMarsthataresimilartothealluvialfansinage
1
andvolume(Kleinhansetal.2010,Hauberetal.2013).Anotherapproachtoestimatingthe
intervalofalluvial-fanbuild-upistomeasurethedensityofcraterssuperimposedon
differentfans,andusethespreadofcrater-retentionagesforthefansurfacesasaproxyfor
therangeoffan-formationages.Thismethodisnotreliable,becausecratercountsforareas
<103km2cannotdistinguishbetweenagesof(forexample)2.5Gaand3.0Ga(Warneretal.
2015).Therefore,thetimespanofhabitableclimatesintheLateHesperian-Amazonian
remainsanopenquestion.
Figure1.LargealluvialfandepositonMars(28°S333°E).The≲3.5Gaageofthefansis
shownbytheirrelativelylowcraterdensity.Aeolianerosionexposeslayersandchannels
withinthedeposit.
2
Figure2.Idealizedcross-sectionthroughanalluvialfandeposit.Craterdensitywithin
alluvialfansmaybeestimatedusingthefrequencyofvisibleinterbeddedcratersatthe
exposedsurface.CaseA:Iffluvialsedimentaccumulationwasmodest,somecratersthat
formedlateinthehistoryofalluvial-fanaggradationareonlypartiallyburiedandarevisible
today(thoughoutnumberedbypostfluvialcraters;notshown).CaseB:Iffluvial
accumulationwasrapid,thetemporalwindowforpartialcraterburialwascorrespondingly
brief,andfewinterbeddedcratersarevisibleattheexposedsurface.CaseC:Ifpostfluvial
erosionwassevere,thearealdensityofexposedembeddedcratersofagivendiameteris
stillproportionaltothevolumetricdensityofthosecraters.Ineithererosionalorwellpreservedfanscenarios,assumingsteadyaggradation,acountofinterbeddedcraters
constrainsthefanaggradationrate.
Togetamoreaccurateestimateoftheintervaloverwhichalluvialfansformed,weusedthe
embedded-cratermethod(Hartmann1974,Kiteetal.2013a).Forpost-Noachiandepositson
Mars,thismethodworksasfollows.Craterdensityonastableplanetarysurface(asurface
thatisneithererodingnoraggrading)isproportionaltoexposureduration.Thismethod
maybeextendedtothreedimensions:thetotalnumberofcratersinterbeddedwithina
3
sedimentarydepositisproportionaltothetimespannedbyactivesedimentation(including
anyhiatuses)(Fig.2).Thegreaterthevolumetricdensityofcraters,thelongerthetimespan
ofdeposition.Ofcourse,manyormostoftheinterbeddedcratersmaybecompletelyburied
(Fig.2);thereforeacompletecountofinterbeddedcratersisusuallyimpossible.However,if
animpactoccursneartheendofactivesedimentation,thenthecorrespondingcratermaybe
onlypartiallyburied.Smallercratersaremorereadilyburied,andlargercratersrequire
moresedimenttobecompletelyobscured.Thetimeneededtoaccumulatethepopulationof
visiblyembedded(synfluvial)cratersis
τraw,D=ND/(fDa)
[1]
whereτraw,Distheminimumtimerequiredtobuilduptheobservedpopulationofembedded
craters(withminimumdiameterD),NDisthenumberofobservedembeddedcraters,fDis
thepastcraterflux(#/km2/yr),andaisthecountarea(km2).Prefluvialcraters(whichare
overlainbyfandeposits,butthatformedbeforethestartoffluvialdeposition;Irwinetal.
2015)areexcluded.Thisproceduregivesastrictlowerlimitontheintervaloffanformation;
itdoesnotaccountforcraterswithinthedepositthatarefullyentombed.Tocorrectforfully
entombedcraters,wecanassumesteadyaggradationanddividefanthicknessZbybest-fit
aggradationrateWtogetdurationoffanformationτsteady:
WD≈1.33Dφ/τD
[2]
τsteady,D=Z/WD
[3]
Thenumeratorin[2]correspondstotherequiredburialdepthforobliteration.Theamount
ofburialthatisrequiredisconstrainedbythegeometryofsmallimpactcraters(Melosh
1989,Wattersetal.2015).φistheobliterationdepthfractionforagivencrater,expressed
relativetocraterdiameter.(Inthispaper,wedefineacrateras“obliterated”ifitcanno
longerbeidentifiedinahigh-resolutionopticalimage;Kite&Mayer2017).Thefactorof
1.33correctsforthefactthat,foranyminimum-diameterD,themediandiameterinacount
willexceedD–thus,Dφisanunderestimateoftherequiredburialdepth.Thiscorrection
dependsonthecraterproductionfunctionused.Thecorrectionisrelativelysmall(1.3×–
1.5×)inoursizerangeofinterest,becausecraterfrequencyfallsoffsteeplywithincreasing
diameter,andwerepresentitbyafixedfactor.Post-depositionalerosionoffandeposits
doesnotaffectthevalidityof[2-3].Arandomlyorientedcutthroughthecrater-containing
volumeintersectseachcraterwithaprobabilityproportionaltothatcrater’ssize;the
sampleofcraterspartiallyexhumedatanerosionalsurfaceisbiasedtowardslarger
impacts.Justaswithpartialburial,thenumberofcratersthatareexposedisproportionalto
thevolumetriccraterdensity(andinverselyproportionaltoaggradationrate).Thus,the
volumetricdensityofinterbeddedcratersmaybeestimatedfromsurfacecountsbothon
pristinefansurfacesandforfansthataredeeplyeroded(Kiteetal.2013a).Therefore,
embedded-cratercountscanbeusedasaMarsfluvialprocessspeedometer.
2.Methodsandresults.
Inordertosetalowerboundonthetimespanofalluvial-fanformation,wesearched
1.7×104km2ofpreviously-cataloguedfans(correspondingtomostofthesurfaceareaof
4
largealluvialfansonMars;Wilsonetal.2012)using6m-per-pixelCTXimagestoscoutfor
candidateembeddedcraters.Candidatecratersshowpossibleevidenceofinterbeddingwith
paleochannelsorotherfandeposits.Eachcandidatefeaturewasreviewedbythreeofthe
authors(E.S.K,D.P.M.,andJ.S.)forfinalclassification.Whereavailable,25cm-per-pixel
HiRISEimagesandanaglyphs,plusCTXDigitalTerrainModels(DTMs)wereusedtoinspect
candidatesflaggedintheinitialCTXsurvey.Eachcandidatefeaturewascategorizedas
qualitylevel1,2,3(representingdecreasingconfidencethatthecraterwasembedded),orit
wasdiscarded.Atotalof25embeddedcraterswerefoundat<30°latitude(D=0.08-5.0km;
Fig.3,Figs.S1-S2).Thelatitudecutwasselectedtoavoidglaciatedcraters.Thoseembedded
craterswerethenclassified(usuallywiththeaidofCTXDTMs)as“synfluvial,”“uncertain”,
or“prefluvial.”Thesecratersconstituteasmallfractionofthetotalnumberofcratersonthe
surfacesofthefans,mostofwhichappeartobepostfluvial(Grant&Wilson2012).
Figure3.EmbeddedsynfluvialcraterswithinalluvialfansonMarsfromourcatalog(Fig.
S2).
AsupplementaryHiRISE-onlysurvey(570km2)wascarriedouttocheckforresolution
effects.13embeddedcraterswerefound(D=0.05–0.22km;FigureS3).For0.08km<D<
0.16km,HiRISEembedded-craterdensitiesND/aarethesame(withinPoissonerror)asCTX
embedded-craterdensities.Therefore,theHiRISEcheckprovidesnoevidencethatour
conclusionswouldbechangedbyaHiRISEre-surveyofthefullareacoveredbyourCTX
survey.
Diametermeasurementerrorwasestimatedbyblindlyremeasuringcratersandfoundtobe
negligible.Craterdiameterchangeduringdegradationisignored.
Thecontributionoffalsepositivestoourcatalogislikelysmall.Althoughpolygonalfaulting
inEarthmarinesedimentscanproducecrater-likeconcentriclayering(Tewkseburyetal.
2014),thisisunlikelyforMarsalluvial-fandeposits.Forexample,embeddedcratersare
isolated(notspace-filling),andfrequentlyshowpreservedrims.However,thereare
certainlyfalsenegativesinoursurveyarea:re-surveyofacraterofinterestfoundseveral
additionalcandidateswithscores≤3onpanelinspection.Itispossiblethatmanyembedded
cratersareeasilyidentifiableasimpactcratersinCTXimagery,buthaveanexpressionthat
isindistinguishablefrompostfluvialcratersatCTXscale.Therefore,ourembeddedcrater
densitiesarelowerlimits.
5
WeestimatedfanthicknessesbydifferencingCTXDTMprofilesacrossfansandanalogous
profilesacrosspartsofthesamefan-hostingcraters(n=17)thatlackedfans(Palucisetal.
2014)(Table1).Wefoundmaximumfanthicknessof1.1km,withthicknesses~1km
common.
Crater
Name
46W0N OrsonWelles
45W24S
Degana
37W18S
Luba
34W26S
Holden
28W26S
Ostrov
Estimatedmaximum
fanthickness
900m
200m
>200m
1.1km
1km
27W28S “SEofOstrov”
15W25S
-
67E,22S
Harris
73E22S
Saheki
74E23S “SEofSaheki”
76E26S
Runanga
84E29S
-
84E33S
Majuro
134E1N
-
142E12N
Eddie
170E19N
Kotka
180E5S
-
>300m
150m
300m
>1.1km
500m
250m
900m
300m–400m
>500m
250m
300m
200m
Comments
Flatfloorapproximation.
Conservativeestimate
Flatfloorapproximation.
Flatfloorapproximation.
Conservativeestimate
Flatfloorapproximation
Flatfloorapproximation.
Alternativecontrolgivesmax.thickness450m.
TakenfromMangoldetal.2012
Table1.Alluvialfanthicknessestimates(crater-by-cratermaximumestimatedthickness).
Thicknesses~1kmarehighlighted.
3.Analysis.
Theusualprocedureforestimatingcrater-countingerroristousePoissonstatistics(Michael
etal.2016).Theresultsofthisprocedureareshownbythebluelinesandblueerrorbarsin
Fig.4.Togeneratetheseresults,weassumedafixedcraterflux(Michaeletal.2013),no
changeinatmosphericscreeningfromtoday’sMars,astrong-rocktargetstrength,anda
fixedobliterationdepthfractionφ=0.1.
Thetrueerrorislargerthanthisbecauseofuncertaintyin:-truecraterflux(Johnsonetal.
2016);targetstrength;filteringbyapotentiallythickerpastatmosphere;thetimeof
formationofthefans;andtheamountofburialorerosion(expressedasafractionofthe
crater’sdiameter)thatisneededtopreventthecraterfrombeingdetectedatCTXresolution.
Therefore,weadoptedconservativepriorprobabilitiesontheseparametersinaMonte
CarlosimulationofourlowerboundthatalsoincludesPoissonerror(detailsaregiveninthe
SupplementaryInformation).Specifically,weassumed(1)afactor-of-4uncertaintyincrater
flux(log-uniformuncertaintybetween0.5×and2×theMichaeletal.2013fluxes);(2)loguniformuncertaintyintargetstrengthbetweenlimitsof65kPaand10MPa(Dundasetal.
2010);(3)log-uniformuncertaintyinpaleo-atmosphericpressurebetweenlimitsof6mbar
and1000mbar;(4)auniformuncertaintybetweenfanformation2.0Ga(low-endfancrater
6
retentionage)and3.6Ga(agesoflargecratersthathostfans);and(5)alog-uniformprior
forobliterationdepthfraction(expressedasafractionofD)from0.05(rimburial;Melosh
1989)and0.2(originalcraterdepth;Wattersetal.2015).ForeachMonteCarlotrial,the
effectofPoissonerroriscalculatedanalytically.Giventheobservationsandtherandomlysampledparameters,eachMonteCarlotrialyieldsananalyticprobabilityforeachcandidate
age(oreachcandidateaggradationrate)ineachsizebin.Theseprobabilitiesaresummed
over103MonteCarlotrialsandnormalized.Resultsareshownbythegraybandsandblack
starsinFig.4.
Thebest-fitlowerlimitonthetime-spanoffanaggradation(Fig.4a)increaseswith
increasingdiameter,asexpected.Bins≥1.4kmcontainonly1crater,andthePoisson
uncertaintyconstitutesmostofthetotaluncertaintyinthelowerlimit.Forthesmaller
diameters,thecounting-statisticserrorissmallcomparedtothetotaluncertaintyinthe
lowerlimit,butsystematicundercountingofembeddedcratersismostlikelyforcratersthat
aresmaller(andthusmoreeasilyburiedandmodified).Forexample,for10mdiameter
craters,thenominaltimespantobuilduptheobservedcratersis~1Myr,buttheamountof
burialneededtoobliteratethecraterisonly~1m,andgivenfanthicknesses~1km(Table
1),itisveryunlikelythattherearezerocratersentombedwithinthelower99.9%ofthe
fans.Thelargest>1km-diameterbincontains2embeddedcratersandyieldsa95%lower
limitof>17Myr,whichweroundto>20Myr.Usingthesingle~5kmcraterfoundinour
surveygivesa>54Myrlowerbound.
(a)
7
(b)
Figure4.(a).Minimumtime-spanofsedimentationbasedonobservedsynfluvialcraters,
binnedbyminimumdiameter(cumulativeplot).Blueerrorbarsbracketthe90%confidence
intervalsonlowerlimit(byPoissonestimation).FullMonteCarlofitcorrespondstothegray
band.Blackzoneisexcludedwith>95%confidence.Whitezoneisexcludedin<5%oftrials.
BlackasteriskscorrespondtothemedianoutcomeoftheMonteCarloprocedure.Red
dashedlinesassumeaconstantaggradationrateandthatburialby10%ofacrater’s
diameterissufficienttoobscurethecraterfromorbiter-imagesurveys.(b)Corresponding
aggradationrateestimatesforalluvialfans,binnedbyminimumdiameter(cumulativeplot).
Errorbarsarethesameasin(a).
Turningtotherateplot(Fig.4b),constantaggradationratesof<(4-8)μm/yearmatchour
data.Differentcraterdiametersprobedifferentdepthrangesandthusaggradationrateover
differenttimescales,butourdatadonotrequireachangeinaggradationratewithtimescale
(Jerolmack&Sadler2007).Dividingtypicalfanthicknessesbythisrategives125-250Myr
(100-300Myrto1s.f.).Forbothplots,theMonteCarloproceduregiveslimitsthataremore
uncertain,andslightlymorepermissive(Fig.4).
Thetimespanofliquidwaterestimatedfromsteadyaggradation(>(100-300)Myr;Fig.4b)
isamorerealisticestimateofthetruefan-formingintervalthanthetimespanestimated
fromobservedsynfluvialcraters(>20Myr;Fig.4a).Forexample,iftheobservedembedded
cratersrepresentthetotalembedded-craterpopulation,thenfanaggradationmusthave
startedveryfastanddecreasedsharplyneartheendoffanbuild-up.Suchanaccumulation
historywouldfavorsmall-craterpreservationrelativetolarge-craterpreservation–
oppositewhatisobserved.Furthermore,afterfluvialdepositionstopped,manyofthefan
depositsunderwentaeolianerosion(e.g.Fig.1).Becausepostfluvialerosionwoulddestroy
8
someembeddedcraters,theobservedembeddedcratersareveryunlikelytorepresentthe
totalpopulationofcratersthatformedembeddedwithinthefans.
Resultsareplottedexcludingcratersforwhichsynfluvialversusprefluvialstatuscouldnot
bedetermined,butincludingcratersofqualityscore3.Thesechoiceshavelittleeffectonour
conclusions.Thatisbecausethequality-3embeddedcraters,andthecraterswhose
synfluvialstatusisuncertain,correspondtosmalldiameterbinsforwhichourcountsare
probablyincomplete.
4.Discussion.
4.1.Factorsnottakenintoaccountwouldincreaseourlowerlimit.
Byfittingasingletime-span(andseparatelyfittingasingleerosionrate)tothefandeposits,
weimplicitlyassumethattheprobabilityoffindinganembeddedcraterisspatiallyuniform.
However,thisassumptionisinconsistentwiththeclumpinessoftheobservedspatial
distributionofembeddedcraters.Forexample,Crater“W”(Kraaletal.2008)has20%ofthe
embeddedcraters(5outof25)eventhoughthefanareaatthatsiteisonly500km2(3%of
total).Thisclusterisunlikelytobeduetochance(assumingcrateringisaPoissonprocess):
iffanareaisdividedinto(17000km2)/(500km2)=34equal-areasites,theexpectednumber
ofcraterspersiteisλ=25/34=0.73.Theprobabilityoffinding≥5embeddedcratersin≥1
sitesisthenonly1–(1–Σx=5∞f(x|0.73))34=3%,wherefisthePoissonprobabilitydistribution
function.Tobeconsistentwithdataatthe50%level,wemustincreaseλby200%.This
mightcorrespondtouniformfanagewithspatiallynon-uniformdetectabilityofembedded
craters.AlternativelyCrater“W”mightrecordanomalouslyslowaggradation.Ineithercase,
ourbestestimateofthetimespanoffan-formingclimatesis200%longerthaninourlower
limit.
Spatialstaggeringoffanaggradationwouldincreaseourlowerlimit.Time-varyingorbital
forcingwouldfavorsnowmelt(e.g.Kiteetal.2013b)indifferentplacesatdifferenttimes.
Localizedprecipitationwouldnotbegloballycorrelated(Kiteetal.2011).
Ifthefandepositsunderwenterosionduringtheperiodofnetaggradation,thenthiswould
destroysomeembeddedcraters.Therefore,iferosionoccurredduringfanaggradation,our
upperlimitwouldincreasefurther.
MarscraterfluxesareextrapolatedfromLunarradiogenicagesandcorrespondingcrater
densities.Thosecraterdensitieshavebeenarguedtobeincorrect(Robbins2014).Adopting
Robbins’chronologywouldreducefluxuncertaintyfrom(1-12)×modern(Neukum
chronologyfunction)to(1.5–2.2)×modern.BecausehighaggradationratesinourMonte
Carlorunsalwayscorrespondtohighcraterfluxes,includingRobbins’chronologywould
raiseourlowerlimit.
9
4.2.Implicationsforpaleohydrologyandclimate.
D<100membeddedcratersplaceanupperlimitonpaleoatmosphericpressure(Kiteetal.
2014).SinceD<50mimpactcratersareextremelyrareonEarth(atmosphericcolumn
density104kg/m2),theexistenceofD<100membeddedcratersintheMarsfanssuggests
atmosphericcolumndensity<2×104kg/m2,i.e.P<1bararoundthetimeoffanaggradation.
Previousanalysesoftheintervaloverwhichalluvialfansformedhavedividedfanvolumeby
theinferredfluvialsedimenttransportflux(e.g.Jerolmacketal.2004).Thisdurationisa
lowerboundontheintervaloverwhichalluvialfansformed,becausenotallyearsneed
producerunoff.Ourlowerboundexceedssedimentologicallowerboundsby>1000×.Many
alluvialfansare~1kmthick.Supposeafan:alcovearearatioof0.5.Typicalwater:sediment
ratiosonEarthare103:1.2000kmofwaterat0.5msnowmelt/yeargives4Myr.These
calculationsarehighlyuncertain.Forexample,iftheamountofsnowmeltislimitedbysnow
supplyto10cm/yr,thenthefanformationtimeis20Myr,equaltoourstrictlowerlimiton
thetotaltimespanofalluvialfanformation.However,givenquasi-periodicorbital
variability,fine-tuningofMars’hydrologicalcycleisrequiredtoproducesmallamountsof
runoffeveryyear,especiallyforourpreferredlowerlimitof>(100-300)Myr.Slownet
aggradationratesinareasofsteeprelief(Fig.1)suggestintermittency.
Intermittencyinalluvial-fan-formingclimateisfurthersuggestedbycombiningourdata
withotherconstraints.Thepaucityofmineralogicevidenceforin-situalterationoffan
deposits(McKeownetal.2013),thepresenceofhydratedsilica(possiblyopal;Carteretal.
2012),andthepersistenceofolivine(Stoparetal.2006),whencombinedwiththe>20Myr
spanofsurfaceliquidwaterrequiredbyourdata,suggestthatclimateconditionswerecold
andthatintermittencyfurtherreducedliquidwaterinteractionwithsoil.Coldconditionsare
alsosuggestedbysedimentary-depositmineralogyatGale(McLennanetal.2014,Siebachet
al.2017).IntermittencyisalsosuggestedbymultiplepulsesoffanformationatHolden
(Irwinetal.2008),Gale,andMelasChasma(Williams&Weitz2014).
Insummary,thedataexcludeanyexplanationthatproducesasingleburstofhabitabilityof
<20Myrduration.Forexample,thedataexclude triggeringbythethermalpulsecausedby
theimpactsthatformedthelargecraterswhichhostthealluvialfans.Thedatadisfavor
fluvialsedimenttransporteveryyearfor>20Myr. Amongotherpossibilities,thedatapermit
a long-livedhabitableenvironment(snowmeltorrainfall);achaostrigger(Bakeretal.
1991);orobliquity-pacedfluvialintermittency(Kiteetal.2013b).
10
Figure5.GeomorphichistoryofMars(Howard2007,Fassett&Head2011).FSV=Fresh
ShallowValleys(Wilsonetal.2016).RSL=RecurringSlopeLineae.Pre-valleynetwork
fluvialsedimenttransportisfromIrwinetal.(2013).
4.3.Implicationsforhabitability.
Thetimespanofalluvial-fanaggradation(§3)isaproxyforthetimespanofspatially
associatedpaleolakes.Thosepaleolakesincludecandidateplayadepositsatthetoesoffans
(e.g.Morganetal.2014),a>100mdeeppaleolakeinCrater“P”(Kraaletal.2008)suggested
byacommonfan-frontal-scarpelevationof-2700m(inourCTXDTMs),andtheEberswalde
paleolake(whichsharesadrainagedividewithfansatHolden;Irwinetal.2015).In
addition,riversandlakesoccurredduringtheLateHesperian-AmazonianinVallesMarineris
(e.g.Mangoldetal.2004)andArabiaTerra(Wilsonetal.2016).Lakedepositshavegood
biosignaturerecoverypotential(Summonsetal.2011),andbiosignaturerecoveryfrom
Proterozoiclakedepositsisroutine(Petersetal.2005).However,biosignatureswouldhave
beendestroyediflakewaterswereoxidizing.
5.Conclusions.
ExplainingyoungalluvialfansonMarsisachallengetoclimatemodels.Todeterminethe
timespanofalluvialfanformingconditions,wecountedembeddedcraters.Wefoundahigh
densityofembeddedcraters,whichrequiresthattheriver-permittingclimate(s)spanned
>20Myr.Ifaggradationwassteadyat<(4-8)μm/yr,whichisconsistentwithourdata,then
fanbuild-uprequired>(100-300)Myr(Table1). Thedatamakethechallengeofexplaining
thealluvialfansmoresevere,becausetheyexcludeasingleshort-livedanomalyasthecause
ofthealluvialfans.
11
Acknowledgements.
WethankAlanHoward,MarisaPalucis,BeckyWilliams,RossIrwin,Jean-PierreWilliams,Bill
Dietrich,andBrianHynek.Wethanktwoanonymousreviewersfortheiruseful
recommendations,andAndrewDombardforeditorialhandling.Alldatamaybeobtainedby
contactingtheleadauthor([email protected]).Wearegratefulforfinancialsupportfrom
theUStaxpayerviaNASAgrantNNX16AG55G.
References.
Armitage,JohnJ.,Warner,NicholasH.,Goddard,Kate,Gupta,Sanjeev,2011.Timescalesofalluvialfan
developmentbyprecipitationonMars.Geophys.Res.Lett.38(17).CiteIDL17203.
Baker,V.R.,Strom,R.G.,Gulick,V.C.,Kargel,J.S.,Komatsu,G.1991.Ancientoceans,icesheetsandthe
hydrologicalcycleonMars.Nature352,589-594.
Carter,J.,etal.,2012.CompositionofDeltasandAlluvialFansonMars.In:43rdLunarandPlanetary
ScienceConference.March19–23,2012,Woodlands,TX.LPIContributionNo.1659,id.1978.
Dundas,C.M.,Keszthelyi,L.P.,Bray,V.J.&McEwen,A.S.Roleofmaterialpropertiesinthecratering
recordofyoungplaty-ridgedlavaonMars.Geophys.Res.Lett.37,L12203(2010).
Fassett,C.I.,andJ.W.Head,2011,SequenceandtimingofconditionsonearlyMars.Icarus,211,
1204-1214.
Grant,J.A.,Wilson,S.A.,2011.LatealluvialfanformationinsouthernMargaritiferTerra,Mars.
Geophys.Res.Lett.38(8).CiteIDL08201.
Grant,J.A.,Wilson,S.A.,2012,Apossiblesynopticsourceofwaterforalluvialfanformationin
southernMargaritiferTerra,Mars,Planetary&SpaceScience
Grindrod,P.M.;Heap,M.J.;Fortes,A.D.;Meredith,P.G.;Wood,I.G.;Trippetta,F.;Sammonds,P.R.2010,
Experimentalinvestigationofthemechanicalpropertiesofsyntheticmagnesiumsulfatehydrates:
ImplicationsforthestrengthofhydrateddepositsonMars,J.Geophys.Res.115,CiteIDE06012.
Hartmann,W.K.,1974,GeologicalobservationsofMartianarroyos.,J.Geophys.Res.79,3951–3957.
Hartmann,W.K.,2005.Martiancratering8:isochronrefinementandthechronologyofMars.Icarus
174,294–320.
Hauber,E.;Platz,T.;Reiss,D.;LeDeit,L.;Kleinhans,M.G.;Marra,W.A.;Haas,T.;Carbonneau,P.,
2013,AsynchronousformationofHesperianandAmazonian-ageddeltasonMarsandimplications
forclimate,J.Geophys.Res.:Planets,118,pp.1529-1544.
Howard,A.D.,2007.SimulatingthedevelopmentofMartianhighlandlandscapesthroughthe
interactionofimpactcratering,fluvialerosion,andvariablehydrologicforcing.Geomorphology91,
332–363.
Irwin,R.P.;Grant,J.A.;Howard,A.D.,2008,TheAlluvialFanComplexinHoldenCrater:Implications
fortheEnvironmentofEarlyMars,39thLunarandPlanetaryScienceConference,abstractid.1869
Irwin,R.P.,Lewis,K.W.,Howard,A.D.,Grant,J.A.,2015.PaleohydrologyofEberswaldecrater,Mars.
Geomorphology.http://dx.doi.org/10.1016/j.geomorph.2014.10.012.
12
Jerolmack,DouglasJ.;Mohrig,David;Zuber,MariaT.;Byrne,Shane,2004,Aminimumtimeforthe
formationofHoldenNortheastfan,Mars,GeophysicalResearchLetters,31,CiteIDL21701.
Jerolmack,D.J.,Sadler,P.,2007.Transienceandpersistenceinthedepositionalrecordofcontinental
margins.J.Geophys.Res.(EarthSurface)112,F03S13.
Johnson,B.C.;Collins,G.S.;Minton,D.A.;Bowling,T.J.;Simonson,B.M.;etal.2016,Spherulelayers,
craterscalinglaws,andthepopulationofancientterrestrialimpactors,Icarus,271,350-359.
Kite,E.S.,Michaels,T.I.,Rafkin,S.C.R.,Manga,M.,&W.E.Dietrich,2011.Localizedprecipitationand
runoffonMars,J.Geophys.Res.Planets,116,E07002.
Kite,E.S.,Lucas,A.,Fassett,C.I.,2013a,PacingearlyMarsriveractivity:embeddedcratersinthe
AeolisDorsaregionimplyriveractivityspanned>∼(1-20)Myr.Icarus225,850–855.
Kite,E.S.,I.Halevy,M.A.Kahre,M.Manga,andM.Wolff,2013b,Seasonalmeltingandtheformation
ofsedimentaryrocksonMars,Icarus,223,181–210.
Kite,E.S.,Williams,J.-P.,Lucas,A.,Aharonson,O.,2014.Lowpalaeopressureofthemartian
atmosphereestimatedfromthesizedistributionofancientcraters.Nat.Geosci.7,335–339.
Kite,E.S.,Howard,A.,Lucas,A.S.,Armstrong,J.C.,Aharonson,O.,Lamb,M.P.,2015.Stratigraphyof
AeolisDorsa,Mars:stratigaphiccontextofthegreatriverdeposits.Icarus253,223-242.
Kite,E.S.;Mayer,D.P.,2017,Marssedimentaryrockerosionratesconstrainedusingcratercounts,
withapplicationstoorganicmatterpreservationandtotheglobaldustcycle,Icarus,286,212-222.
Kleinhans,M.G.,vandeKasteele,H.E.,Hauber,E.,2010.Palaeoflowreconstructionfromfandelta
morphologyonMars.EarthPlanet.Sci.Lett.294,378–392.
Kraal,E.R.,etal.,2008,Catalogueoflargealluvialfansinmartianimpactcraters,Icarus,194,101110.
Mangold,N.,etal.,2004.EvidenceforprecipitationonMarsfromdendriticvalleysintheValles
MarinerisArea.Science305,78–81.
McLennan,S.M.,etal.,2014.ElementalgeochemistryofsedimentaryrocksatYellowknifeBay,Gale
Crater,Mars.Science343(6169).http://dx.doi.org/10.1126/science.1244734
Melosh,H.J.,1989.Impactcratering:ageologicprocess.In:ResearchsupportedbyNASA,253.Oxford
UniversityPress,NewYork,p.11.(OxfordMonographsonGeologyandGeophysics,No.11),1989.
Michael,G.G.;Kneissl,T.;Neesemann,A.,2016,Planetarysurfacedatingfromcratersize-frequency
distributionmeasurements:Poissontiminganalysis,Icarus,277,279-285.
Moore,J.M.,&Howard,A.D.,2005,LargealluvialfansonMars,J.Geophys.Res.,110,E4,E04005.
Michael,G.G.,2013.Planetarysurfacedatingfromcratersize-frequencydistributionmeasurements:
multipleresurfacingepisodesanddifferentialisochronfitting.Icarus,226,885–890.
Morgan,A.M.,etal.,2014.Sedimentologyandclimaticenvironmentofalluvialfansinthemartian
SahekicraterandacomparisonwithterrestrialfansintheAtacamaDesert.Icarus,229,131–156.
Neukum,G.;Ivanov,B.A.;Hartmann,W.K.,2001,CrateringRecordsintheInnerSolarSystemin
RelationtotheLunarReferenceSystem,SpaceScienceReviews,v.96,Issue1/4,p.55-86.
13
Okubo,C.H.,2007,Strengthanddeformabilityoflight-tonedlayereddepositsobservedbyMER
Opportunity:EagletoErebuscraters,Mars.GeophysicalResearchLetters34,CiteIDL20205.
Palucis,M.,etal.2014,TheoriginandevolutionofthePeaceVallisfansystemthatdrainstothe
Curiositylandingarea,GaleCrater,Mars,J.Geophys.Res.Planets.,119,Issue4,pp.705-728
Peters,K.E.,etal.(eds.),2005,Thebiomarkerguide:volume2,CambridgeUniversityPress.
Pilorget,C.,&F.Forget,2016,FormationofgulliesonMarsbydebrisflowstriggeredbyCO2
sublimation,NatureGeoscience9,65-69.
Robbins,StuartJ.,2014,Newcratercalibrationsforthelunarcrater-agechronology,Earthand
PlanetaryScienceLetters,Volume403,p.188-198.
Siebach,K.L.,M.B.Baker,J.P.Grotzinger,S.M.McLennan,R.Gellert,L.M.Thompson,J.A.Hurowitz,
2017,SortingoutCompositionalTrendsinSedimentaryRocksoftheBradburyGroup(AeolusPalus),
GaleCrater,Mars,J.Geophys.Res.Planets,DOI:10.1002/2016JE005195
Stopar,J.D.,G.J.Taylor,V.E.Hamilton,L.Browning(2006)Kineticmodelofolivinedissolutionand
extentofaqueousalterationonMars,GCA,70,p.6136-6152
Summons,R.E.,etal.,2011.PreservationofMartianorganicandenvironmentalrecords:finalreport
oftheMarsBiosignatureWorkingGroup.Astrobiology11,157–181
Tewksbury,B.J.;Hogan,J.P.;Kattenhorn,S.A.;Mehrtens,C.J.;Tarabees,E.A.,2014,Polygonalfaults
inchalk:InsightsfromextensiveexposuresoftheKhomanFormation,WesternDesert,Egypt
Geology,vol.42,issue6,pp.479-482.
Warner,N.H.;Gupta,S.;Calef,F.;Grindrod,P.;Boll,N.;Goddard,K.,2015,Minimumeffectiveareafor
highresolutioncratercountingofmartianterrains,Icarus,Volume245,p.198-240.
Watters,W.A.,Geiger,L.M.,Fendrock,M.,Gibson,R.,2015.Morphometryofsmallrecentimpact
cratersonMars:sizeandterraindependence,short-termmodification.J.Geophys.Res.(Planets)120,
226–254.
Williams,RebeccaM.E.;Malin,MichaelC.,2008,Sub-kilometerfansinMojaveCrater,Mars,Icarus,
198,365-383.
Williams,R.M.E.,DeanneRogers,A.,Chojnacki,M.,Boyce,J.,Seelos,K.D.,Hardgrove,C.,Chuang,F.,
2011.EvidenceforepisodicalluvialfanformationinfarwesternTerraTyrrhena,Mars,Icarus211
(1),222–237.
Williams,R.M.E.etal.,2013,MartianFluvialConglomeratesatGaleCrater,Science,Volume340,
Issue6136,pp.1068-1072.
Williams,R.M.E.,Weitz,C.M.,2014.Reconstructingtheaqueoushistorywithinthesouthwestern
MelasBasin,Mars:cluesfromstratigraphicandmorphometricanalysesoffans.Icarus242,19–37.
Wilson,S.A.;Grant,J.A.;Howard,A.D.,2012,DistributionofIntracraterAlluvialFansandDeltaic
DepositsintheSouthernHighlandsofMars,43rdLunar&Planet.Sci.Conf.,abstractid.2462.
Wilson,S.A.;Howard,A.D.;Moore,J.M.;Grant,J.A.,2016,Acold-wetmiddle-latitudeenvironmenton
MarsduringtheHesperian-Amazoniantransition:EvidencefromnorthernArabiavalleysand
paleolakes,J.Geophys.Res.:Planets,121,1667-1694.
14