Arxiv:1802.10422V1 [Astro-Ph.EP] 28 Feb 2018 Can Preclude Melting by Decreasing the Strength of Green- Uid Water for Sediment Induration at Low Latitudes
Total Page:16
File Type:pdf, Size:1020Kb
Effect of Mars Atmospheric Loss on Snow Melt Potential in a 3.5-Gyr Mars Climate Evolution Model Megan Mansfield∗ and Edwin S. Kite University of Chicago Department of Geophysical Sciences Michael A. Mischna Jet Propulsion Laboratory, California Institute of Technology (Dated: March 1, 2018) Post-Noachian Martian paleochannels indicate the existence of liquid water on the surface of Mars after about 3.5 Gya [27, 53]. In order to explore the effects of variations in CO2 partial pressure and obliquity on the possibility of surface water, we created a zero-dimensional surface energy balance model. We combine this model with physically consistent orbital histories to track conditions over the last 3.5 Gyr of Martian history. We find that melting is allowed for atmospheric pressures corresponding to exponential loss rates of dP=dt / t−3:73 or faster, but this rate is within 0:5σ of the rate calculated from initial measurements made by the Mars Atmosphere and Volatile EvolutioN (MAVEN) mission, if we assume all the escaping oxygen measured by MAVEN comes from atmospheric CO2 [44, 62]. Melting at this loss rate matches selected key geologic constraints on the formation of Hesperian river networks, assuming optimal melt conditions during the warmest part of each Mars year [27, 36, 37, 60]. The atmospheric pressure has a larger effect on the surface energy than changes in Mars's mean obliquity. These results show that initial measurements of atmosphere loss by MAVEN are consistent with atmospheric loss being the dominant process that switched Mars from a melt-permitting to a melt-absent climate [29], but non-CO2 warming will be required if < 2 Gya paleochannels are confirmed, or if most of the escaping oxygen measured by MAVEN comes from H2O. I. INTRODUCTION Hecht [23] used a similar model to determine the nec- essary surface properties, such as albedo and conductiv- Large (> 10 km2), late-stage Martian alluvial fans and ity, to allow melting on Mars. Additionally, Hecht em- river deltas provide evidence for surface liquid water on phasized the importance of evaporative cooling at low post-Noachian Mars after about 3.5 Gya [20, 27, 53, 67]. atmospheric pressures, as first pointed out by Ingersoll However, by about 3.5 Gya most of the conditions fa- [26]. Models which have not included evaporative cool- vorable to the existence of surface liquid water no longer ing produce much more melting on the surface of Mars, existed: much of Mars's atmosphere was lost before the which indicates its importance in setting melting condi- end of the Noachian and the Martian dynamo shut down tions [9, 65]. The high rate of evaporative cooling at low around the mid-Noachian [43, 54]. Additionally, the atmospheric pressure precludes melting of pure water ice amount of greenhouse gases released by volcanism on on present-day Mars, which has an average atmospheric Mars is less than on Earth because the mantle is more pressure of 600 Pa. reducing [59], and by 3.5 Gya the rate of volcanic de- Although atmospheric pressure can have a strong ef- gassing had slowed down significantly [32]. Because or- fect on the surface energy balance, changes in Mars's biters detect only minor post-Noachian carbonate, which orbital parameters, especially large changes in its mean may be due to SO2 or acidity [5, 21], it is difficult to jus- obliquity, can also affect melt conditions. Kite et al. [33] tify post-Noachian carbonate sequestration of more CO2 determined the conditions for snowmelt at a variety of at- than exists in the present, thin atmosphere [15]. mospheric pressures, obliquities, and eccentricities. They Atmospheric pressure has a strong effect on permitting found that combining these factors could produce inter- surface liquid water because lower atmospheric pressure mittent melting that could provide small amounts of liq- arXiv:1802.10422v1 [astro-ph.EP] 28 Feb 2018 can preclude melting by decreasing the strength of green- uid water for sediment induration at low latitudes. house warming and increasing the amount of evaporative While models that study moments in time can pro- cooling. Clow [8] showed these effects in a model which vide some insight into the past Martian climate, models determined the minimum atmospheric pressure necessary showing evolution over time give a better understanding to melt dusty snow on the surface of Mars. They found of the timing of melting in Mars's past. Manning et al. that melting could occur at relatively low pressures if [47] addressed the time evolution of the Martian climate the snowpack was assumed to be thick enough that melt by examining the relative sizes of various CO2 reservoirs accumulated at the base of the thick snow layer while over the history of Mars. Their model determined that the upper layers of snow were at colder temperatures. different climate states are stable for different obliqui- ties, but they only considered a subset of the full range of obliquities that Mars has experienced. ∗ meganmansfi[email protected] A model that simultaneously studies the effects of 2 time evolution, orbital variations, and changing energy obliquity jumps at different times in their timelines, but balance with changing atmospheric pressure is needed which all end with obliquities within 12◦ of the modern- to more fully understand the timing and intermittency day value of 25:19◦. Several obliquity tracks were con- of melting on Mars. To address this, we created a sidered because the chaotic nature of the Solar System zero-dimensional energy balance model that spans over means that the exact value of Mars's obliquity is not 3.5 Gyr of Mars history and examines how changes in or- known beyond a few hundred million years ago. There- bital parameters, solar luminosity, and atmospheric pres- fore, we studied a distribution of possible orbital histo- sure impact the melt conditions of snow on the surface ries to statistically constrain the influence of the orbital of Mars. We consider snowmelt as a source of water parameters on melting. The eight orbital histories all as opposed to catastrophic flows after impacts because produced similar energy histories, so only two will be ex- recent work suggests lakes lasted for at least a few thou- amined in detail in this paper. Obliquity histories for sand years, and this relatively long lifetime implies a cli- tracks 1 and 2, which will be examined in more detail mate favorable to melting and not just a sudden catas- later, are shown in Figure2. trophic event that allowed melting for a brief period [20, 27, 53, 63, 67]. We also ignore the effects of im- pacts and eruptions because changes due to atmospheric loss, solar brightening, and orbital variations are more well understood, and because the melt potential of im- pacts themselves depends on the atmospheric heat capac- ity and therefore the atmospheric pressure. Our model is zero-dimensional because this allows fast computation of a full 3.5-Gyr energy balance, whereas previous 3D mod- FIG. 1. A flowchart showing how key aspects of the model els have been unable to include time evolution because of were incorporated. Boxes outlined in red indicate new contri- computational limits [69]. This model is timely because butions from our model. Atmospheric loss, solar brightening, the Mars Atmosphere and Volatile EvolutioN (MAVEN) and obliquity histories were calculated as described in Section II.1 and combined in an energy balance described in section mission is currently in orbit around Mars collecting data II.3. on rates of atmospheric loss, and is presumably providing better constraints on the atmospheric pressure over time on Mars which we can use in our model to constrain the conditions that explain melting on post-Noachian Mars II.2. Loss of Atmospheric Pressure [44]. We describe our model in SectionII. In Section III we present our results, which indicate that atmospheric pressure has a more dominant effect on surface melt con- The chaotic orbital histories were combined with a ditions than obliquity variations. In SectionIV we dis- model of solar brightening over time and loss of atmo- cuss assumptions of our model and possible future exten- spheric pressure [2]. Although there are several possible sions, and we summarize our findings in SectionV. loss mechanisms for CO2 in past Mars, such as chemical fixation of CO2 in deep aquifers [7] and basal melting of a CO2 ice cap [39], we chose to focus on escape to space and assumed all atmospheric loss was due to escape to II. MODEL DESCRIPTION space. Additionally, we neglected carbonate formation because it is thought to be relatively ineffective in the II.1. Chaotic Orbital Histories post-Noachian, although this is debated [15, 25]. The current atmospheric pressure on Mars is approxi- To investigate the influence of orbital variability on mately 600 Pa, but at the temperatures that would allow post-Noachian surface liquid water, we created a model liquid water on the surface of Mars, the approximately that combined the effects of orbital variations, atmo- 600 Pa of additional CO2 currently trapped in Martian spheric loss, and solar brightening, as illustrated in Fig- polar caps would be released into the atmosphere, and so ure1. Physically consistent orbital histories were con- we assume a modern-day atmospheric pressure of 1200 Pa structed using the mercury6 N-body integrator and an [4]. Atmospheric loss to space was parameterized as a obliquity code [1,6]. The orbital histories recorded the power law orbital parameters of the simulated Mars at 200-year in- tervals over the last 3.5 Gyr and included quasi-periodic dP t α = −k 0 (1) variations in orbital parameters [35]. Most importantly, dt t they included changes in obliquity, which for Mars has probably jumped from low mean values around 20◦ to where P is the atmospheric pressure in Pa, t is time ◦ high mean values around 40 between two and nine times in Gyr since the Sun formed, t0 = 4:56 Gyr, and k = −1 in its history on ' 200-Myr timescales [41].