Study of the Moon and comparative planetology: Problems and perspectives

On the base of the comparative planetologic study of the Moon and terrestrial planets two fundamental features of their history and structure have been established.Firstly, shell-like structure of the terrestrial planets could be understood only in the terms of the heterogeneous accretion theory. At the final stages of major terrestrial planet formation the leading role belonged to the planetosimals of carbonaceous chondritic composition. Secondly, there are two types of the crust on the planetary surface. Their formation are considered to be independent and differing in the geological time. The primary planetary crust of predominantly feldspathic composition is considered to form during the pregeologic period at the final stage of planetary formation due to the impact-explosive processes. The hydrosphere and atmosphere is thought to appear contemporaneously. The basaltic planetary crust is forming later due to the radioactive decay and superimposed on the primary feldspathic crust.     

3D radar wavefield tomography of comet interiors

Answering fundamental questions about the origin and evolution of small planetary bodies hinges on our ability to image their surface and interior structure in detail and at high resolution. The interior structure is not easily accessible without systematic imaging using, e.g., radar transmission and reflection data from multiple viewpoints, as in medical tomography. Radar tomography can be performed using methodology adapted from terrestrial exploration seismology. Our feasibility study primarily focuses on full wavefield methods that facilitate high quality imaging of small body interiors. We consider the case of a monostatic system (co-located transmitters and receivers) operated in various frequency bands between 5 and 15?MHz, from a spacecraft in slow polar orbit around a spinning comet nucleus. Using realistic numerical experiments, we demonstrate that wavefield techniques can generate high resolution tomograms of comets nuclei with arbitrary shape and complex interior properties.     

Co-registration of Chang’E-1 stereo images and laser altimeter data with crossover adjustment and image sensor model refinement

Topographic information is of fundamental importance for various scientific investigations in lunar and planetary exploration. To provide high-precision, seamless mapping capability, it is critical to co-register image and altimetry data, the two major data sources for topographic modeling. This paper presents a method for co-registration of Chang’E-1 (CE-1) stereo images and laser altimeter (LAM) data with crossover adjustment and refinement of the image sensor model. The crossover adjustment is tested in a larger area (0–60°N, 50–0°W); the image refinement and co-registration with LAM data are tested in an area (46.2–50.0°N, 31.8–28.8°W) within the larger area. Experimental results demonstrate that this co-registration reduces the mean differences of inconsistency from more than 200 m to just 3.21 m in the Z direction of object space. In image space, the mean errors of homologous points both in the column and row directions are reduced to below 0.1 pixel. This indicates that the proposed crossover adjustment of LAM data and refinement of the CE-1 stereo image model can effectively improve co-registration of the two data sets.     

Antiparticle content in the magnetosphere

In the present work we assess the stable and transient antiparticle content of planetary magnetospheres, and subsequently we consider their capture and application to high delta-v space propulsion. We estimate the total antiparticle mass contained within the Earth’s magnetosphere to assess the expediency of such usage. Using Earth’s magnetic field region as an example, we have considered the various source mechanisms that are applicable to a planetary magnetosphere, the confinement duration versus transport processes, and the antiparticle loss mechanisms. We have estimated the content of the trapped population of antiparticles magnetically confined following production in the exosphere due to nuclear interactions between high energy cosmic rays (CR) and constituents of the residual planetary upper atmosphere.The galactic antiprotons that directly penetrate into the Earth’s magnetosphere are themselves secondary by its nature, i.e. produced in nuclear reactions of the cosmic rays passing through the interstellar matter. These antiproton fluxes are modified, dependent on energy, when penetrating into the heliosphere and subsequently into planetary magnetospheres. During its lifetime in the Galaxy, CR pass through the small grammage of the interstellar matter where they produce secondary antiprotons. In contrast to this, antiprotons generated by the same CR in magnetosphere are locally produced at a path length of several tens g/cm² of matter in the ambient planetary upper atmosphere. Due to the latter process, the resulting magnetically confined fluxes significantly exceed the fluxes of the galactic antiprotons in the Earth’s vicinity by up to two orders of magnitude at some energies.The radiation belt antiparticles can possibly be extracted with an electromagnetic-based “scoop” device. The antiparticles could be concentrated by and then stored within the superimposed magnetic field structure of such a device. In future developments, it is anticipated that the energy of the captured antiparticles (both rest energy and kinetic energy) can be adapted for use as a fuel for propelling spacecraft to high velocities for remote solar system missions.     

The Planetary Laboratory for Image Analysis (PLIA)

The observation, characterization, and understanding of planetary atmospheres are key components of Solar System exploration. Their study requires, among other observations, the use of image data acquired from spacecrafts. The basic tasks that are generally used to work with planetary image data are: image navigation, projection, image processing, limb-darkening correction, composition of multiwavelength images and mosaicing. We illustrate the mathematical techniques required for the study of planetary images obtained from spacecrafts using as an example our own software package (PLIA). We present examples of analysis of Galileo SSI, Cassini ISS and Venus Express VIRTIS observations covering a large range of image peculiarities. The paper is organized as a guide for anyone willing to replicate one or several of the capabilities incorporated in the PLIA software. An open source version of PLIA for astronomical image processing is available.     

Crater detection via genetic search methods to reduce image features

Recent approaches to crater detection have been inspired by face detection’s use of gray-scale texture features. Using gray-scale texture features for supervised machine learning crater detection algorithms provides better classification of craters in planetary images than previous methods. When using Haar features it is typical to generate thousands of numerical values from each candidate crater image. This magnitude of image features to extract and consider can spell disaster when the application is an entire planetary surface. One solution is to reduce the number of features extracted and considered in order to increase accuracy as well as speed. Feature subset selection provides the operational classifiers with a concise and denoised set of features by reducing irrelevant and redundant features. Feature subset selection is known to be NP-hard. To provide an efficient suboptimal solution, four genetic algorithms are proposed to use greedy selection, weighted random selection, and simulated annealing to distinguish discriminate features from indiscriminate features. Inspired by analysis regarding the relationship between subset size and accuracy, a squeezing algorithm is presented to shrink the genetic algorithm’s chromosome cardinality during the genetic iterations. A significant increase in the classification performance of a Bayesian classifier in crater detection using image texture features is observed.     

The first new application of the Mathematical Theory of Stochastic Processes to lunar and planetary science: Topography-Profile Diagrams of Mars

The Mathematical Statistics Theory (MST) and the Mathematical Theory of Stochastic Processes (MTSP) are different branches of the more general Mathematical Probability Theory (MPT) that can be used to investigate physical processes through mathematics. Each model of a stochastic process, according to MTSP, can provide one or more interpretations in the MST domain. A large body of work on impact crater statistics according to MST exists, showing cumulative crater frequency (N km⁻²) as a function of age (years) for some particular crater diameter. However, this is only one possible representation in the MST domain of the bombardment of the planetary surface modeled as a stochastic process according to MTSP. The idea that other representations are possible in the MST domain of the same stochastic process from MTSP has been recently presented. The importance of the approach is that each such mathematical-based interpretation can provide a large amount of new information. Coupled with MOLA data, Topography-Profile Diagrams (TPD) are one of the many examples that can provide a large amount of new information regarding the history of Mars. TPD consists of: (1) Topography-Profile Curve (TPC), which is a representation of the planet’s topography, (2) Density-of-Craters Curve (DCC), which represents density of craters, (3) Filtered-DCC (FDCC), which represents DCC filtered by a low-pass filter, included with the purpose of reducing the noise, and (4) Level-of-Substance-Over-Time Curve (LSOTC), which represents interpretation of the influence on the distribution of craters shown by FDCC. TPC uniquely corresponds to the computation of TPD, whereas DCC depends on algorithms for computing the elevation of each crater according to the topography, center coordinates, and radius of impact crater, and FDCC relies on the architecture of the custom designed low-pass filter for filtering DCC. However, all variations of DCC and FDCC, which includes the various impact crater data sets, showed a correlation among the density of craters and elevation over 70–80% of the planet surface. Additionally, if we assume that the ocean primarily caused the noted correlation, LSOTC offers a mathematical approach for estimating topographic change of the ocean’s extent over time. Accordingly, TPD is the first new practical application of MTSP to lunar and planetary sciences, showing correlation of topography to a physical process.     

Integrated method for crater detection from topography and optical images and the new PH9224GT catalogue of Phobos impact craters

In a large majority of lunar and planetary surface images, impact craters are the most abundant geological features. Therefore, it is not surprising that crater detection algorithms (CDAs) are one of the most studied subjects of image processing and analysis in lunar and planetary science. In this work we are proposing an Integrated CDA, consisting of: (1) utilization of DEM (digital elevation map)-based CDA; (2) utilization of an optical-based CDA; (3) re-projection of used datasets and crater coordinates from normal to rotated view and back; (4) correction of the brightness and contrast of a used optical image; and (5) tile generation for the optical-based CDA and an assembling of results with an elimination of multiple detections, in combination with a pyramid approach down to the resolution of the available DEM image; and (6) a final integration of the results of DEM-based and optical-based CDAs, including a removal of duplicates. The proposed CDA is applied to one specific asteroid-like body, the small Martian moon Phobos. The experimental evaluation of the proposed CDA is done by a manual verification of crater-candidates and a search for uncatalogued craters. The evaluation has shown that the proposed CDA was used successfully for cataloging Phobos craters. The major result of this paper is the PH9224GT – currently the most complete global catalogue of the 9224 Phobos craters. The possible applications of the new catalogue are: (1) age estimations for any selected location; and (2) comparison/evaluation of the different chronology and production functions for Phobos. This confirms the practical applicability of the new Integrated CDA – an additional result of this paper, which can be used in order to considerably extend the current crater catalogues.     

Thermal protection system development,testing, and qualification for atmospheric probes and sample return missions: Examples for Saturn,Titan and Stardust-type sample return

The science community has continued to be interested in planetary entry probes, aerocapture, and sample return missions to improve our understanding of the Solar System. As in the case of the Galileo entry probe, such missions are critical to the understanding not only of the individual planets, but also to further knowledge regarding the formation of the Solar System. It is believed that Saturn probes to depths corresponding to 10 bars will be sufficient to provide the desired data on its atmospheric composition. An aerocapture mission would enable delivery of a satellite to provide insight into how gravitational forces cause dynamic changes in Saturn’s ring structure that are akin to the evolution of protoplanetary accretion disks. Heating rates for the “shallow” Saturn probes, Saturn aerocapture, and sample Earth return missions with higher re-entry speeds (13–15 km/s) from Mars, Venus, comets, and asteroids are in the range of 1–6 KW/cm². New, mid-density thermal protection system (TPS) materials for such probes can be mission enabling for mass efficiency and also for use on smaller vehicles enabled by advancements in scientific instrumentation. Past consideration of new Jovian multiprobe missions has been considered problematic without the Giant Planet arcjet facility that was used to qualify carbon phenolic for the Galileo probe. This paper describes emerging TPS technologies and the proposed use of an affordable, small 5 MW arcjet that can be used for TPS development, in test gases appropriate for future planetary probe and aerocapture applications. Emerging TPS technologies of interest include new versions of the Apollo Avcoat material and a densified variant of Phenolic Impregnated Carbon Ablator (PICA). Application of these and other TPS materials and the use of other facilities for development and qualification of TPS for Saturn, Titan, and Sample Return missions of the Stardust class with entry speeds from 6.0 to 28.6 km/s are discussed.     

Catalysis of peptide formation by inorganic oxides: high efficiency of alumina under mild conditions on the Earth-like planets.

Amino acid condensation catalyzed by inorganic oxides is a widely recognized way for prebiotic peptide formation. Silica and alumina are widely distributed in the Earth-like planets' crust as minerals of different complexity, and thus are attractive model catalysts for the studies of abiotic peptide synthesis. Experiments performed in other laboratories have shown that this process can be efficient at > 80 degrees C, which is not easy to find on the planetary surface in combination with sufficient concentrations of amino acids and necessary catalysts. In the present work we tested catalytic activity of three forms of alumina (which proved to be an efficient catalyst for this process) in the intermolecular condensation of L-alanine. We expanded the temperature interval down to 55 degrees C and used the simplest permanent heating procedure, without employing fluctuating drying/wetting conditions. The most important finding is that even under the lowest temperature considered (i.e. 55 degrees C), short peptide formation can be detected already after 10-30 days of heating. This fact implies that the abiotic peptide formation might occur in a wide variety of planetary environments, without need for high temperatures, given the presence of amino acid building blocks and alumina-containing minerals.     

Detection of sub-kilometer craters in high resolution planetary images using shape and texture features

Counting craters is a paramount tool of planetary analysis because it provides relative dating of planetary surfaces. Dating surfaces with high spatial resolution requires counting a very large number of small, sub-kilometer size craters. Exhaustive manual surveys of such craters over extensive regions are impractical, sparking interest in designing crater detection algorithms (CDAs). As a part of our effort to design a CDA, which is robust and practical for planetary research analysis, we propose a crater detection approach that utilizes both shape and texture features to identify efficiently sub-kilometer craters in high resolution panchromatic images. First, a mathematical morphology-based shape analysis is used to identify regions in an image that may contain craters; only those regions - crater candidates - are the subject of further processing. Second, image texture features in combination with the boosting ensemble supervised learning algorithm are used to accurately classify previously identified candidates into craters and non-craters. The design of the proposed CDA is described and its performance is evaluated using a high resolution image of Mars for which sub-kilometer craters have been manually identified. The overall detection rate of the proposed CDA is 81%, the branching factor is 0.14, and the overall quality factor is 72%. This performance is a significant improvement over the previous CDA based exclusively on the shape features. The combination of performance level and computational efficiency offered by this CDA makes it attractive for practical application.     

Planetary landings with terrain sensing and hazard avoidance: A review

The site selection and certification processes for planetary landers equipped with hazard avoidance capability are reviewed. The prior (large-scale) ‘landability’ of the target areas determined from orbital remote sensing can be only a few tens of per cent (～40 % for Perseverance, ～80–90 % for Apollo 11 and 60–70 % for subsequent Apollo landers) because on-board sensing is able to find safe areas at smaller scales (meters to tens of meters) than the delivery ellipse which may be several to several tens of kilometers across. This contrasts with the ‘blind’ landings of unguided missions, where safe terrain occupying 95–99 % of the landing ellipse are typically sought. The particulars of Apollo 11 and Perseverance/Ingenuity are discussed, together with the similar Chang-E-3 and Tianwen-1 Moon and Mars landers, and the Hayabusa-2 and OSIRIS-REx asteroid contacts, since these missions all used on-board terrain-relative navigation to steer relative to hazards either mapped previously or detected in real-time. These missions set the context for the application of these techniques to the Dragonfly mission to Titan, which has a more austere remote sensing basis on which to select a landing site, but whose rotor propulsion allows substantial divert capability.     

Planetary protection policy overview and application to future missions.

The 1967 treaty on the peaceful uses of outer space reflected both concerns associated with the unknown nature of the space environment and the desire of the world scientific community to preserve the pristine nature of celestial objects until such times as they could be studied in an effective manner. Since 1967, NASA has issued policy directives that have adopted the guidelines of COSPAR for protecting the planets from contamination by Earth organisms and for protecting the Earth from the unknown. This paper presents the current status of planetary protection (quarantine) policy within NASA, and a prospectus on how planetary protection and back contamination issues might be addressed in relation to future missions envisioned for development by NASA either independently, or in cooperation with the space agencies of other nations.     

Modelling of cometary nuclei: Planetary missions preparation

The planning of planetary missions requires the knowledge as much as possible accurate of the target. This knowledge, which cannot be obtained only by ground based observations, can be supported by theoretical modeling. This is particularly true in the case of Rosetta, a cornerstone ESA mission that will be launched in January 2003, and of its target, the comet 46P/Wirtanen. In this paper we show how, using a nucleus thermal evolution model, it is possible to foresee the temperature of the nucleus surface and the activity level along the orbit and define the range of possible values for these properties. Activity level close to the aphelion depends on the presence of CO. Surface temperatures depend on the physical properties attributed to the dust: temperatures reached by a crust mainly composed by silicatic grains differ from the temperatures reached by a crust mainly composed by organic grains.     

Planetary penetrators: Their origins,history and future

Penetrators, which emplace scientific instrumentation by high-speed impact into a planetary surface, have been advocated as an alternative to soft-landers for some four decades. However, such vehicles have yet to fly successfully. This paper reviews in detail, the origins of penetrators in the military arena, and the various planetary penetrator mission concepts that have been proposed, built and flown. From the very limited data available, penetrator developments alone (without delivery to the planet) have required ∼$30M: extensive analytical instrumentation may easily double this. Because the success of emplacement and operation depends inevitably on uncontrollable aspects of the target environment, unattractive failure probabilities for individual vehicles must be tolerated that are higher than the typical ‘3-sigma’ (99.5%) values typical for spacecraft. The two pathways to programmatic success, neither of which are likely in an austere financial environment, are a lucky flight as a ‘piggyback’ mission or technology demonstration, or with a substantial and unprecedented investment to launch a scientific (e.g. seismic) network mission with a large number of vehicles such that a number of terrain-induced failures can be tolerated.     

An introduction to the lunar and planetary science activities in Korea

Korea is planning a series of lunar space programs in 2020 starting with a lunar orbiter and a lander with a rover. Compared to other countries, Korea has a relatively brief history in space and planetary sciences. With the expected Korean missions on the near-term horizon and the relatively few Korean planetary scientists, Korea Institute of Geoscience and Mineral Resources (KIGAM) has established a new planetary research group focusing on development of prospective lunar instruments, analysis of the publicly available planetary data of the Moon, organizing nationwide planetary workshops, and initiating planetary educational programs with academic institutions. Korea has also initiated its own rocket development program, which could acquire a rocket-launch capability toward the Korean lunar mission. For the prospective Korea’s lunar science program, feasibility studies for some candidate science payloads have been started since 2010 for an orbiter and a lander. The concept design of each candidate instrument has been accomplished in 2012. It is expected that the development of science payloads may start by 2014 as Phase A. Not only developing hardware required for the lunar mission but also educational activities for young students are high priorities for Korea. The new plan of the Korean lunar mission can be successfully accomplished with international cooperative outreach programs in conjunction with internationally accessible planetary data system (PDS). This paper introduces the KIGAM’s international cooperative planetary research and educational programs and also summarizes other nationwide new developments for Korean lunar research projects at Kyung Hee University and Hanyang University.     

Environmental degradation analysis using NOAA/AVHRR data

This paper proposes a particular approach to assess information about soil degradation, based on a methodology to calculate soil color from NOAA/AVHRR data. As erosive processes change physical and chemical properties of the soil, altering, consequently, the superficial color, monitoring the change in color over time can help to identify and analyze those processes. A relationship among the soil color (described in the Munsell Color System, i.e., in terms of Hue, Value and Chroma), vegetation indices, surface temperature and emissivity has been established, which is based on the theoretical model. The methodology has three main phases: determination of the regression models among soil color and vegetation indices, emissivity and surface temperature; generation of digital soil color models; and statistical evaluation of the estimated color. The tests showed that the methodology is efficient in determining soil color using the various vegetation indices (i.e., Normalized vegetation index NDVI, Modified soil adjusted vegetation index MSAVI). One vegetation index, i.e., Purified adjusted vegetation index (PAVI) is proposed to subsidies the effect of vegetation over the soil. Best results were obtained for the Hue color component. To further test the methodology, the estimated digital color models were compared with the characteristic color of soil classes in the test area. The results of this application confirmed the methodology’s capacity to determine the soil color from NOAA/AVHRR data. This type of study is quite helpful to know the erosion of soil as well as some abrupt change in soil due to natural hazards by space borne or air-borne sensors.     

我国中低纬地区上空低电离层中的行星波

本文通过5年的电离层吸收观测资料与平流层增温事件对比及吸收资料的谱分析,得出以下几点初步结论:1)极区平流层增温事件的影响可能通过子午环流和行星波传播,经过5—9天后到达中低纬地区,从而引起那里的电离层吸收变化;2)冬季行星波沿子午方向的平均速度大约在10m/s到15m/s之间变化;3)全年均有周期为32天、18天、10天、8天和2天的行星波出现,它对大气湍流系数有明显影响。计算得出行星波扰动引起中层的NO浓度偏离未扰值可高达40%。      

On the relationship between clear-sky planetary and surface albedos: a parameterization for simple energy balance climate models

In simple energy balance climate models all physical and dynamical processes are parameterized in terms of the single unknown variable: the surface temperature. To simulate the ice-albedo feedback, the surface albedo is usually assumed to be a function of surface temperature. But to compute the absorbed solar radiation in such models one requires the top-of-the atmosphere albedo: the planetary albedo. In the present study, a simple linear relationship is derived between planetary albedo and surface albedo for the case of clear skies. The relationship is based upon a regression equation derived from simulations and has a standard error of estimate of 0.028. The estimation of planetary albedo from surface albedo is checked by comparing zonally averaged clear-sky planetary albedos, estimated from zonally averaged surface albedos, to satellite determinations of zonally averaged minimum albedos for monthly mean conditions. The minimum albedos are assumed to be representative of the clear-sky planetary albedos. The results show root-mean square differences of 0.05 between the estimated clear-sky planetary albedos and the minimum albedos.More accurate relationships can be obtained if one uses an additional parameter - the solar zenith angle. In this case, the standard errors of estimate are reduced to 0.017 for a zenith angle of 0°, 0.018 for a zenith angle of 60° and 0.021 for a zenith angle of 85°.     

Bacterial cellulose may provide the microbial-life biosignature in the rock records

Bacterial cellulose (BC) is a matrix for a biofilm formation, which is critical for survival and persistence of microbes in harsh environments. BC could play a significant role in the formation of microbial mats in pristine ecosystems on Earth. The prime objective of this study was to measure to what extent spectral and other characteristics of BC were changed under the performance of BC interaction with the earthly rock – anorthosite – via microorganisms. The spectral analyses (Fourier Transform Infrared FT-IR, spectroscopy, and atomic absorption spectroscopy) showed unprecedented accumulation of chemical elements in the BC-based biofilm. The absorption capacity of IR by BC was shielded a little by mineral crust formed by microorganisms on the BC-based biofilm surface, especially clearly seen in the range of 1200–900 cm⁻¹ in FT-IR spectra. Confocal scanning laser microscopy analysis revealed that elements bioleached from anorthosite created surface coats on the BC nanofibril web. At the same time, the vibrational spectra bands showed the presence of the characteristic region of anomeric carbons (960–730 cm⁻¹), wherein a band at 897 cm⁻¹ confirmed the presence of β-1, 4-linkages, which may serve as the cellulose fingerprint region. Results show that BC may be a biosignature for search signs of living organisms in rock records.