DOCSLIB.ORG

Arxiv:1410.0082V2 [Astro-Ph.EP] 21 Oct 2014 Gration, As Opposed to Once Per Time Step

Total Page：16

File Type：pdf, Size：1020Kb

Download full-text PDF Read full-text
Arxiv:1410.0082V2 [Astro-Ph.EP] 21 Oct 2014 Gration, As Opposed to Once Per Time Step Draft version October 22, 2014 Preprint typeset using LATEX style emulateapj v. 5/2/11 NEAR-EARTH ASTEROID SATELLITE SPINS UNDER SPIN-ORBIT COUPLING Shantanu P. Naidu1, Jean-Luc Margot1,2 (Dated:) Draft version October 22, 2014 ABSTRACT We develop a fourth-order numerical integrator to simulate the coupled spin and orbital motions of two rigid bodies having arbitrary mass distributions under the influence of their mutual gravitational potential. We simulate the dynamics of components in well-characterized binary and triple near- Earth asteroid systems and use surface of section plots to map the possible spin configurations of the satellites. For asynchronous satellites, the analysis reveals large regions of phase space where the spin state of the satellite is chaotic. For synchronous satellites, we show that libration amplitudes can reach detectable values even for moderately elongated shapes. The presence of chaotic regions in the phase space has important consequences for the evolution of binary asteroids. It may substantially increase spin synchronization timescales, delay BYORP-type evolution, and extend the lifetime of binaries. The variations in spin rate due to large librations also affect the analysis and interpretation of lightcurve and radar observations. Subject headings: minor planets, asteroids: general { minor planets, asteroids: individual (2000 DP107, 1999 KW4, 2002 CE26, 2004 DC, 2003 YT1, Didymos, 1991 VH, 2001 SN263, 1994 CC), planets and satellites: dynamical evolution and stability 1. INTRODUCTION semi-major axis, eccentricity, and component sizes. In Binary near-Earth asteroids (NEAs) are numerous in practice, only systems observed with radar fall in this the asteroid population. Both radar and lightcurve data class. Over 35 binary NEAs have been observed with have shown that ∼ 16% of NEAs larger than ∼200 m di- radar, but only about ten have sufficient data to yield ameter have satellites (Pravec et al. 1999; Margot et al. mutual orbits and component size estimates. We apply 2002; Pravec et al. 2006). It is now widely accepted our technique to these systems. that binary NEAs form by a spin-up process (Margot Sections 2 and 3 describe the implementation of our et al. 2002) and that the specific spin-up mechanism is coupled spin-orbit integrator and cover energy and an- the YORP torque (Rubincam 2000). Binary NEA sys- gular momentum conservation properties. Section 4 ex- tems exhibit interesting post-fission and spin-orbit dy- plains different kinds of satellite spin librations and sets namics (e.g. Ostro et al. 2006; Scheeres et al. 2006; Fahne- up the notation used in subsequent sections. In section 5, stock and Scheeres 2008; McMahon and Scheeres 2013) we examine the spin-orbit coupling effect and compare that profoundly affect their evolution (e.g. Jacobson and numerical and analytical estimates of libration ampli- Scheeres 2011a; Fang and Margot 2012a; Jacobson et al. tudes. Section 6 introduces surface of section plots which 2014), but the range of dynamical regimes has not been are used to identify resonant, chaotic, and non-resonant fully explored. quasi-periodic trajectories. We examine the trajectories In this paper, we develop a method for simulating the and spin configurations of satellites in well-characterized coupled spin and orbital motions of two rigid bodies with binary and triple near-Earth asteroid systems in section 7 arbitrary mass distributions. This technique is signifi- and show that large chaotic zones exist in the phase space cantly faster than a similar implementation by Fahne- of known asynchronous satellites. We also compute libra- stock and Scheeres (2006), because in our implementa- tion amplitudes for synchronous satellites. We discuss tion the computationally expensive volume integrals over implications of the results in section 8. the two bodies are computed only once before the inte- arXiv:1410.0082v2 [astro-ph.EP] 21 Oct 2014 gration, as opposed to once per time step. We use our 2. NUMERICAL INTEGRATION technique to perform a survey of the dynamics of all well- We numerically investigate the coupled spin and or- characterized binary NEA systems and map the range of bital dynamics of two extended rigid objects under their dynamical behaviors, including the spin configurations mutual gravitational influence. We neglect the transla- of asteroid satellites. These results provide important tional motion of the system barycenter and use the 6 insights for modeling the physical properties of binaries first-order differential equations of motion (EOMs) de- and for understanding the long term evolution of the bi- rived by Maciejewski (1995). Here we express these nary systems. EOMs in the body-fixed frame of the primary: The sample of well-characterized binaries includes all NEA systems with known estimates of system mass, @V P P_ = P × Ω − ; R_ = R × Ω + ; 1 @R 1 m 1 Department of Earth, Planetary, and Space Sciences, Uni- _ _ versity of California, Los Angeles, CA 90095, USA Γ2 = Γ2 × Ω1 + µ2; Γ1 = Γ1 × Ω1 + µ1; 2 Department of Physics and Astronomy, University of Cali- fornia, Los Angeles, CA 90095, USA S_ = SΩ^ 2 − Ω^ 1S; S_1 = S1Ω^ 1: (1) 2 Here R and P are the relative position and linear mo- inertia integrals are listed below: mentum vectors of the secondary with respect to the pri- Z mary, respectively, m = mpms=(mp + ms) is the reduced 2 (−A + B + C) Ix2 = x dm = ; mass of the system, where mp and ms are the masses of 2 the primary and secondary, respectively, V is the mutual Z 2 (A − B + C) gravitational potential, µ's are the torque vectors acting Iy2 = y dm = ; on the two components, Ω's are their angular velocity 2 Z vectors, and Γ's are their angular momentum vectors. 2 (A + B − C) I 2 = z dm = : (4) Subscripts 1 and 2 denote quantities that refer to the z 2 primary and the secondary, respectively. Further, S and S1 are attitude rotation matrices: the former mapping Here, A ≤ B ≤ C are the principal moments of inertia from the secondary frame to the primary frame, and the of the object about the x, y, and z axes, respectively. latter mapping from the primary frame to the inertial We use the Cash-Karp method (Cash and Karp 1990) frame. A hat (^) symbol above a vector specifies an op- to integrate equations (1). It is a fifth-order Runge-Kutta erator that maps a 3-vector (e.g., v = [vx; vx; vz]) to an integrator with adaptive stepsize control which uses an antisymmetric 3x3 matrix, as follows: embedded fourth-order Runge-Kutta formula to compute errors. We use the implementation provided by Press " 0 −vz vy # et al. (1992) and set the fractional error tolerance to −15 v^ = vz 0 −vx : (2) 10 . −vy vz 0 In all simulations, we assume a planar system, i.e., both bodies are in principal axis rotation about their The term @V=@R is the gradient of the mutual gravita- z (shortest) axes and their equatorial planes are aligned tional potential, which is the gravitational force (vector) with the mutual orbit at all times. We start all simu- between the two components. All vectors in equations lations at the pericenter of the osculating mutual orbit (1) are expressed in the body-fixed frame of the primary. and with the longest axis of each body pointing towards However, when computing Ω^ 2, one must express Ω2 in each other. The system parameters and initial osculating the body-fixed frame of the secondary. mutual orbital parameters for all simulations are given The gravitational force and torques are computed at in Table 1. each time step according to the method detailed in Ashenberg (2007). They are functions of R, S, and the 3. ENERGY AND ANGULAR MOMENTUM inertia integrals of the two bodies. The inertia integrals CONSERVATION encode the mass distribution information of a body and In this section, we describe results of tests designed to are of the form: evaluate the energy and angular momentum conservation Z properties of the integrator. Figure 1 illustrates a repre- p q r sentative test run with the parameters given in the first Ixpyq zr = x y z dm; (3) B line of Table 1. For this test case, we used triaxial ellip- soids with principal axis half-lengths of a = 600, b = 500, where dm is a mass element of body B at body-fixed co- and c = 400 m for the primary, and a = 252, b = 229, ordinates (x, y, z) and the integral is a volume integral c = 190 m for the secondary. The initial spin periods of over the entire body. The body-fixed coordinate system the primary and the secondary are 2.775 h and 32.59 h, is aligned with the principal axes and its origin is at the respectively. center of mass of B. The exponents p, q, and r are ei- The total energy is conserved at a level of 10−2 Joules ther 0 or positive integers, such that p + q + r > 0. We per year, which is about 10−11 times the mean orbital use inertia integrals up to fourth order in the integra- energy and less than 10−8 times the magnitude of the tions, where the order of an inertia integral is given by energy exchanged between the binary components and the sum of exponents, i.e., p + q + r. The inertia inte- the mutual orbit. grals depend only on the mass distribution of the object Angular momentum is conserved at a level of and remain constant throughout the integration, so we 220 kg m2 s−1 per year, which is less than 10−11 compute them only once before the integration. At each times the total angular momentum of the system time step, current values of R and S from the integrator (∼ 1014 kg m2 s−1) and less than 10−8 times the angular are passed as arguments to the modules that compute momentum exchanged between the component spins and force and torques.
Load more

Recommended publications
  • Supplementary Information For
    Supplementary Information for Vertical Angular Momentum Constraint on Lunar Formation and Orbital History ZhenLiang Tian and Jack Wisdom Z. Tian, J. Wisdom E-mail: [email protected], [email protected] This PDF file includes: Supplementary text Figs. S1 to S4 Table S1 References for SI reference citations www.pnas.org/cgi/doi/10.1073/pnas.2003496117 ZhenLiang Tian and Jack Wisdom 1 of 11 Supporting Information Text Model Comparison. Our model, which is the same as in (1), differs from the model of (2) in the following aspects: 1. Probably the most important difference is the tidal model. We use the conventional Darwin-Kaula constant Q model (3). In this model the tidal potential is expanded in a Fourier series and the response to each term is given a phase delay to model dissipation, by analogy to the damped harmonic oscillator. The explicit form builds in considerable constraints from orbital mechanics. There is little observational evidence to constrain the phase delays, but results from the analysis of lunar laser ranging (LLR) data are consistent with the Darwin-Kaula choice to set all the phase delays to be equal (4). Ćuk et al. (2) stated that there is a problem with the constant Q model: that it predicts a decay of orbital eccentricity that, to leading order, is independent of eccentricity. However, this is not true. Using the expressions in (3), one can show that the constant Q model gives a decay of eccentricity that is, to leading order, proportional to eccentricity. Nevertheless, given this misunderstanding of the constant Q model, Ćuk et al.
    [Show full text]
  • The Chaotic Rotation of Hyperion*
    ICARUS 58, 137-152 (1984) The Chaotic Rotation of Hyperion* JACK WISDOM AND STANTON J. PEALE Department of Physics, University of California, Santa Barbara, California 93106 AND FRANGOIS MIGNARD Centre d'Etudes et de Recherehes G~odynamique et Astronomique, Avenue Copernic, 06130 Grasse, France Received August 22, 1983; revised November 4, 1983 A plot of spin rate versus orientation when Hyperion is at the pericenter of its orbit (surface of section) reveals a large chaotic zone surrounding the synchronous spin-orbit state of Hyperion, if the satellite is assumed to be rotating about a principal axis which is normal to its orbit plane. This means that Hyperion's rotation in this zone exhibits large, essentially random variations on a short time scale. The chaotic zone is so large that it surrounds the 1/2 and 2 states, and libration in the 3/2 state is not possible. Stability analysis shows that for libration in the synchronous and 1/2 states, the orientation of the spin axis normal to the orbit plane is unstable, whereas rotation in the 2 state is attitude stable. Rotation in the chaotic zone is also attitude unstable. A small deviation of the principal axis from the orbit normal leads to motion through all angles in both the chaotic zone and the attitude unstable libration regions. Measures of the exponential rate of separation of nearby trajectories in phase space (Lyapunov characteristic exponents) for these three-dimensional mo- tions indicate the the tumbling is chaotic and not just a regular motion through large angles. As tidal dissipation drives Hyperion's spin toward a nearly synchronous value, Hyperion necessarily enters the large chaotic zone.
    [Show full text]
  • Turbulence, Entropy and Dynamics
    TURBULENCE, ENTROPY AND DYNAMICS Lecture Notes, UPC 2014 Jose M. Redondo Contents 1 Turbulence 1 1.1 Features ................................................ 2 1.2 Examples of turbulence ........................................ 3 1.3 Heat and momentum transfer ..................................... 4 1.4 Kolmogorov’s theory of 1941 ..................................... 4 1.5 See also ................................................ 6 1.6 References and notes ......................................... 6 1.7 Further reading ............................................ 7 1.7.1 General ............................................ 7 1.7.2 Original scientific research papers and classic monographs .................. 7 1.8 External links ............................................. 7 2 Turbulence modeling 8 2.1 Closure problem ............................................ 8 2.2 Eddy viscosity ............................................. 8 2.3 Prandtl’s mixing-length concept .................................... 8 2.4 Smagorinsky model for the sub-grid scale eddy viscosity ....................... 8 2.5 Spalart–Allmaras, k–ε and k–ω models ................................ 9 2.6 Common models ........................................... 9 2.7 References ............................................... 9 2.7.1 Notes ............................................. 9 2.7.2 Other ............................................. 9 3 Reynolds stress equation model 10 3.1 Production term ............................................ 10 3.2 Pressure-strain interactions
    [Show full text]
  • Long-Term Evolution of Orbits About a Precessing Oblate Planet: 1. The
    Published in: “Celestial Mechanics and Dynamical Astronomy” (2005) 91 : 75 108 − Long-term evolution of orbits about a precessing∗ oblate planet: 1. The case of uniform precession. Michael Efroimsky US Naval Observatory, Washington DC 20392 USA e-mail: me @ usno.navy.mil August 2, 2018 Abstract It was believed until very recently that a near-equatorial satellite would always keep up with the planet’s equator (with oscillations in inclination, but without a secular drift). As explained in Efroimsky and Goldreich (2004), this misconception originated from a wrong interpretation of a (mathematically correct) result obtained in terms of non-osculating orbital elements. A similar analysis carried out in the language of osculating elements will endow the planetary equations with some extra terms caused arXiv:astro-ph/0408168v4 22 May 2010 by the planet’s obliquity change. Some of these terms will be nontrivial, in that they will not be amendments to the disturbing function. Due to the extra terms, the variations of a planet’s obliquity may cause a secular drift of its satellite orbit inclination. In this article we set out the analytical formalism for our study of this drift. We demonstrate that, in the case of uniform precession, the drift will be extremely slow, because the first-order terms responsible for the drift will be short-period and, thus, will have vanishing orbital averages (as anticipated 40 years ago by Peter Goldreich), while the secular terms will be of the second order only. However, it turns out that variations of the planetary precession make the first-order terms secular.
    [Show full text]
  • Early Evolution of the Earth-Moon System with a Fast-Spinning
    Icarus 256 (2015) 138–146 Contents lists available at ScienceDirect Icarus journal homepage: www.elsevier.com/locate/icarus Early evolution of the Earth–Moon system with a fast-spinning Earth Jack Wisdom ⇑, ZhenLiang Tian Massachusett Institute of Technology, Cambridge, MA 02139, USA article info abstract Article history: The isotopic similarity of the Earth and Moon has motivated a recent investigation of the formation of the Received 14 October 2014 Moon with a fast-spinning Earth (Cuk, M., Stewart, S.T., [2012]. Science, doi:10.1126/science.1225542). Revised 21 February 2015 Angular momentum was found to be drained from the system through a resonance between the Moon Accepted 22 February 2015 and Sun. They found a narrow range of parameters that gave results consistent with the current angular Available online 7 March 2015 momentum of the Earth–Moon system. However, a tidal model was used that was described as approx- imating a constant Q tidal model, but it was not a constant Q model. Here we use a conventional constant Keywords: Q tidal model to explore the process. We find that there is still a narrow range of parameters in which Earth angular momentum is withdrawn from the system that corresponds roughly to the range found earlier, Moon Dynamics but the final angular momentum is too low to be consistent with the Earth–Moon system. Exploring a broader range of parameters we find a new phenomenon, not found in the earlier work, that extracts angular momentum from the Earth–Moon system over a broader range of parameters. The final angular momentum is more consistent with the actual angular momentum of the Earth–Moon system.
    [Show full text]
  • 1. INTRODUCTION Imations, We Thought It Worthwhile to Verify the Evolution the Origin of the Moon Is Still Unclear
    THE ASTRONOMICAL JOURNAL, 115:1653È1663, 1998 April ( 1998. The American Astronomical Society. All rights reserved. Printed in the U.S.A. RESONANCES IN THE EARLY EVOLUTION OF THE EARTH-MOON SYSTEM JIHAD TOUMA McDonald Observatory, University of Texas, Austin, TX 78712; touma=harlan.as.utexas.edu AND JACK WISDOM Department of Earth, Atmospheric, and Planetary Sciences, Massachusetts Institute of Technology, Cambridge, MA 02139; wisdom=mit.edu Received 1997 November 25; revised 1997 December 26 ABSTRACT Most scenarios for the formation of the Moon place the Moon near Earth in low-eccentricity orbit in the equatorial plane of Earth. We examine the dynamical evolution of the Earth-Moon system from such initial conÐgurations. We Ðnd that during the early evolution of the system, strong orbital resonances are encountered. Passage through these resonances can excite large lunar orbital eccentricity and modify the inclination of the Moon to the equator. Scenarios that resolve the mutual inclination problem are presented. A period of large lunar eccentricity would result in substantial tidal heating in the early Moon, providing a heat source for the lunar magma ocean. The resonance may also play a role in the formation of the Moon. Key words: celestial mechanics, stellar dynamics È Earth È Moon regression of the MoonÏs node. With so many approx- 1. INTRODUCTION imations, we thought it worthwhile to verify the evolution The origin of the Moon is still unclear. The classic sce- of the system with numerical simulation. Our hope was that narios for the formation of the Moon include coaccretion, through direct numerical simulation we would uncover Ðssion, and capture.
    [Show full text]
  • Evolution of a Terrestrial Multiple Moon System
    THE ASTRONOMICAL JOURNAL, 117:603È620, 1999 January ( 1999. The American Astronomical Society. All rights reserved. Printed in U.S.A. EVOLUTION OF A TERRESTRIAL MULTIPLE-MOON SYSTEM ROBIN M. CANUP AND HAROLD F. LEVISON Southwest Research Institute, 1050 Walnut Street, Suite 426, Boulder, CO 80302 AND GLEN R. STEWART Laboratory for Atmospheric and Space Physics, University of Colorado, Campus Box 392, Boulder, CO 80309-0392 Received 1998 March 30; accepted 1998 September 29 ABSTRACT The currently favored theory of lunar origin is the giant-impact hypothesis. Recent work that has modeled accretional growth in impact-generated disks has found that systems with one or two large moons and external debris are common outcomes. In this paper we investigate the evolution of terres- trial multiple-moon systems as they evolve due to mutual interactions (including mean motion resonances) and tidal interaction with Earth, using both analytical techniques and numerical integra- tions. We Ðnd that multiple-moon conÐgurations that form from impact-generated disks are typically unstable: these systems will likely evolve into a single-moon state as the moons mutually collide or as the inner moonlet crashes into Earth. Key words: Moon È planets and satellites: general È solar system: formation INTRODUCTION 1. 1000 orbits). This result was relatively independent of initial The ““ giant-impact ÏÏ scenario proposes that the impact of disk conditions and collisional parameterizations. Pertur- a Mars-sized body with early Earth ejects enough material bations by the largest moonlet(s) were very e†ective at clear- into EarthÏs orbit to form the Moon (Hartmann & Davis ing out inner disk materialÈin all of the ICS97 simulations, 1975; Cameron & Ward 1976).
    [Show full text]
  • Interview with Peter Goldreich
    PETER GOLDREICH (b. 1939) INTERVIEWED BY SHIRLEY K. COHEN March, April and November 1998 ARCHIVES CALIFORNIA INSTITUTE OF TECHNOLOGY Pasadena, California Preface to the LIGO Series Interviews The interview of Peter Goldreich (1998) was originally done as part of a series of 15 oral histories conducted by the Caltech Archives between 1996 and 2000 on the beginnings of the Laser Interferometer Gravitational-Wave Observatory (LIGO). Many of those interviews have already been made available in print form with the designation “The LIGO Interviews: Series I.” A second series of interviews was planned to begin after LIGO became operational (August 2002); however, current plans are to undertake Series II after the observatory’s improved version, known as Advanced LIGO, begins operations, which is expected in 2014. Some of the LIGO Series I interviews (with the “Series I” designation dropped) have now been placed online within Caltech’s digital repository, CODA. All Caltech interviews that cover LIGO, either exclusively or in part, will be indexed and keyworded for LIGO to enable online discovery. The original LIGO partnership was formed between Caltech and MIT. It was from the start the largest and most costly scientific project ever undertaken by Caltech. Today it has expanded into an international endeavor with partners in Europe, Japan, India, and Australia. As of this writing, 760 scientists from 11 countries are participating in the LSC—the LIGO Scientific Collaboration. http://resolver.caltech.edu/CaltechOH:OH_Goldreich_P Subject area Physics, geology, planetary science, astronomy, astrophysics, LIGO Abstract Interview in five sessions in March, April, and November 1998 with Peter Goldreich, Lee A.
    [Show full text]
  • Planète 1 Planète
    Planète 1 Planète Une planète est un corps céleste orbitant autour du Soleil ou d'une autre étoile de l'Univers et possédant une masse suffisante pour que sa gravité la maintienne en équilibre hydrostatique, c'est-à-dire sous une forme presque sphérique. La définition scientifique rigoureuse, d'une part veut qu'on réserve le nom de « planète » aux objets du système solaire (dans les autres cas, on parle d'exoplanètes), d'autre part exige que ce corps céleste ait éliminé tout corps rival se déplaçant sur une orbite proche (cela peut signifier soit en faire un de ses satellites, soit provoquer sa destruction par collision). Ce dernier critère ne s'applique pas aux exoplanètes[1] . On estime que le nombre de planètes dans notre seule galaxie est de 1000 milliards[2] . Selon cette définition (qui fut approuvée le 24 août 2006 en clôture de la 26e Assemblée Générale de l'Union astronomique internationale (UAI) huit planètes ont été recensées dans le système solaire : Mercure, Vénus, la Terre, Mars, Jupiter, Saturne, Uranus et Neptune (les quatre premières étant des planètes telluriques alors que les quatre suivantes sont des géantes gazeuses). En même temps que la définition d'une planète était clarifiée, l'UAI définissait comme étant une planète naine un objet céleste répondant à tous les critères, sauf l'élimination des corps sur une orbite proche. Contrairement à ce que suggère l'usage habituel d'un adjectif, une planète naine n'est pas une planète. On compte actuellement cinq planètes naines dans le système solaire : Cérès, Pluton, Makemake, Haumea et Éris[3] .
    [Show full text]
  • Long Term Evolution of the Solar System Jack Wisdom
    LONG TERM EVOLUTION OF THE SOLAR SYSTEM JACK WISDOM Massachusetts Institute of Technology Abstract. The mapping method of Wisdom (1982) has been generalized to encompass all n-body problems with a dominant central mass (Wisdom and Holman, 1991). The new mapping method is presented as well as a number of initial applications. These include billion year integrations of the outer planets, a number of 100 million year integrations of the whole solar system, and a systematic survey of test particle stability in the outer solar system. 1· Introduction A new mapping method for the numerical integration of planetary orbits is pre- sented in Wisdom and Holman (1991). This mapping method is typically about an order of magnitude faster than other methods of integration which have been used in studies of the long-term evolution of planetary orbits. Here, the basic idea of the new η-body mappings is first presented. This is followed by a presentation of new billion year integrations of the outer planets carried out with the new mapping. The new integrations are compared to the 845 million year integration of the outer planets performed with a conventional integration technique on the Digital Orrery. Preliminary results of a number of 100 million year integrations of the whole solar system are then presented, followed by a quick summary of a systematic survey of test particle stability in the outer solar system. 2. Mapping Method A complete presentation of the new mapping method can be found in Wisdom and Holman (1991). The new mapping method is a generalization of the mapping method of Wisdom (1982).
    [Show full text]
  • Linda T. Elkins-Tanton Publications
    Linda T. Elkins-Tanton Publications Submitted 69. Scheinberg, A., L.T. Elkins-Tanton, S. Zhong, Timescale and morphology of Martian mantle overturn immediately following magma ocean solidification, submitted. 68. Erkaev, N. V., H. Lammer, L. Elkins-Tanton, P. Odert, K. G. Kislyakova, Yu. N. Kulikov, M. Leitzinger, M. Gu#del, Escape of the martian protoatmosphere, submitted. 67. Fu, Roger R. and L. T. Elkins-Tanton, The fate of magmas in planetesimals and the retention 1 of primitive chondritic crusts, submitted. 66. Black, B.A, J.-F. Lamarque, C. Shields, L. Elkins-Tanton, J. Kiehl, Acid rain and ozone depletion from pulsed Siberian Traps magmatism, submitted. 65. Black, B.A., E.H. Hauri, L.T. Elkins-Tanton, Sulfur isotopic evidence for sources of volatiles in Siberian Traps magmas, submitted. 2013 64. Elkins-Tanton, L.T., Evolutionary dichotomy for rocky planets, Nature (News and Views), 2013. 63. Lammer H., M. Blanc, W. Benz, M. Fridlund, V. Coudé du Foresto, M. Güdel, H. Rauer, S. Udry, R.-M. Bonnet, M. Falanga, D. Charbonneau, R. Helled, W. Kley, J. Linsky, L. T. Elkins-Tanton, Y. Alibert, E. Chassefière, T. Encrenaz, A. P. Hatzes, D. Lin, R. Liseau, W. Lorenzen, S. N. Raymond, The Science of Exoplanets and their Systems. Accepted at Astrobiology. 62. Elkins-Tanton, L.T., What makes a habitable planet? Eos, Transactions of the American Geophysical Union 94, 149-150, 2013. 61. Mandler, B.E. and L. T. Elkins-Tanton, The origin of eucrites, diogenites and olivine diogenites: magma ocean crystallization and shallow magma chamber processes on Vesta, Meteoritics & Planetary Science, 1-17, doi:10.1111/maps.12135, 2013.
    [Show full text]
  • I. Seasonal Changes in Titan's Cloud Activity II. Volatile
    I. Seasonal Changes in Titan’s Cloud Activity II. Volatile Ices on Outer Solar System Objects Thesis by Emily L. Schaller In Partial Fulfillment of the Requirements for the Degree of Doctor of Philosophy California Institute of Technology Pasadena, California 2008 (Defended April 28, 2008) ii c 2008 Emily L. Schaller All Rights Reserved iii Acknowledgements I would first like to thank my parents Meg and Dave Schaller. Starting with the picture book of the solar system when I was five years old and continuing throughout my PhD they have always been supportive of my love of science. This thesis is dedicated to them. I next want to thank Dave Lamb for his love and support throughout this process. I’ve actually really enjoyed being a grad student and I think that has been in large part due to the people in GPS. I thank my advisor Mike Brown, for his support and willingness to let me dive in and explore different areas of planetary science that were interesting to me — even if they weren’t exactly what I was supposed to be doing. I really value all that I have learned from him and the tools he has given me over these past few years. Thanks to all the people who have passed through Mike’s group including Antonin Bouchez, Chad Trujillo, Henry Roe, Lindsey Malcolm, Sarah Horst, Kris Barkume, Darin Ragozzine, and Meg Schwamb: I have learned something from each of them. Thanks to Kris Barkume for her friendship and for those late night bottles of wine. I thank my roommates (Rowena Lohman, Colette Salyk, Amy Hofmann, and, of course, Harry the cat) for putting up with me.
    [Show full text]