Reddish coloration and color diversity among icy bodies in the outer solar system are one of the significant clues for understanding the status and history of the solar system. However, the origin of color distribution remains debatable. Here, we demonstrate reddish coloration that is stable only at cryogenic temperatures in a laboratory experiment. The reddish coloration was produced on methanol- and water-containing ice irradiated with nitrogen-containing cryoplasma at 85 K. The reddish color visually faded and disappeared at 120-150 K as the ice was heated, unlike well-known refractory organic tholins that are stable even when heated to room temperature. This temperature dependence of reddish coloration under cryogenic conditions could provide a new possible explanation for the absence of ultra-red coloration closer to the Sun in the outer solar system. Our result implies that a reddish material specific to cryogenic environments is useful for the investigation of color diversity and formation mechanism of the outer solar system.

The so-called Lidov-Kozai oscillation is very well known and applied to various problems in solar system dynamics. This mechanism makes the orbital inclination and eccentricity of the perturbed body in the circular restricted three-body system oscillate with a large amplitude under certain conditions. It is widely accepted that the theoretical framework of this phenomenon was established independently in the early 1960s by a Soviet Union dynamicist (Michail L'vovich Lidov) and by a Japanese celestial mechanist (Yoshihide Kozai). A large variety of studies has stemmed from the original works by Lidov and Kozai, now having the prefix of "Lidov-Kozai" or "Kozai- Lidov." However, from a survey of past literature published in late nineteenth to early twentieth century, we have confirmed that there already existed a pioneering work using a similar analysis of this subject established in that period. This was accomplished by a Swedish astronomer, Edvard Hugo von Zeipel. In this monograph, we first outline the basic framework of the circular restricted three-body problem including typical examples where the Lidov-Kozai oscillation occurs. Then, we introduce what was discussed and learned along this line of studies from the early to mid-twentieth century by summarizing the major works of Lidov, Kozai, and relevant authors. Finally, we make a summary of von Zeipel's work, and show that his achievements in the early twentieth century already comprehended most of the fundamental and necessary formulations that the Lidov-Kozai oscillation requires. By comparing the works of Lidov, Kozai, and von Zeipel, we assert that the prefix "von Zeipel-Lidov-Kozai" should be used for designating this theoretical framework, and not just Lidov-Kozai or Kozai-Lidov.

To reproduce the orbits and masses of the terrestrial planets (analogs) of the solar system, most studies scrutinize simulations for success as a batch. However, there is insufficient discussion in the literature on the likelihood of forming planet analogs simultaneously in the same system (analog system). To address this issue, we performed 540 N-body simulations of protoplanetary disks representative of typical models in the literature. We identified a total of 194 analog systems containing at least three analogs, but only 17 systems simultaneously contained analogs of the four terrestrial planets. From an analysis of our analog systems, we found that, compared to the real planets, truncated disks based on typical outcomes of the Grand Tack model produced analogs of Mercury and Mars that were too dynamically cold and located too close to the Venus and Earth analogs. Additionally, all the Mercury analogs were too massive, while most of the Mars analogs were more massive than Mars. Furthermore, the timing of the Moon-forming impact was too early in these systems, and the amount of additional mass accreted after the event was too great. Therefore, such truncated disks cannot explain the formation of the terrestrial planets. Our results suggest that forming the four terrestrial planets requires disks with the following properties: 1) Mass concentrated in narrow core regions between ~0.7-0.9 and ~1.0-1.2 au; 2) an inner region component starting at ~0.3-0.4 au; 3) a less massive component beginning at ~1.0-1.2 au; 4) embryos rather than planetesimals carrying most of the disk mass; and 5) Jupiter and Saturn placed on eccentric orbits.

Context: Distant trans-Neptunian objects are subject to planetary perturbations and galactic tides. The former decrease with the distance, while the latter increase. In the intermediate regime where they have the same order of magnitude (the 'inert Oort cloud'), both are weak, resulting in very long evolution timescales. To date, three observed objects can be considered to belong to this category. Aims: We aim to provide a clear understanding of where this transition occurs, and to characterise the long-term dynamics of small bodies in the intermediate regime: relevant resonances, chaotic zones (if any), and timescales at play. Results: There exists a tilted equilibrium plane (Laplace plane) about which orbits precess. The dynamics is integrable in the low and high semi-major axis regimes, but mostly chaotic in between. From 800 to 1100 au, the chaos covers almost all the eccentricity range. The diffusion timescales are large, but not to the point of being indiscernible in a 4.5 Gyrs duration: the perihelion distance can actually vary from tens to hundreds of au. Orbital variations are favoured in specific ranges of inclination corresponding to well-defined resonances. Starting from uniform distributions, the orbital angles cluster after 4.5 Gyrs for semi-major axes larger than 500 au, because of a very slow differential precession. Conclusions: Even if it is characterised by very long timescales, the inert Oort cloud is much less inert than it appears. Orbits can be considered inert over 4.5 Gyrs only in small portions of the space of orbital elements, which include (90377) Sedna and 2012VP113. Effects of the galactic tides are discernible down to semi-major axes of about 500 au. We advocate including the galactic tides in simulations of distant trans-Neptunian objects, especially when studying the formation of detached bodies or the clustering of orbital elements.

Itokawa particles, which are the recovered samples from the S-complex asteroid 25143 Itokawa by the Hayabusa spacecraft, demonstrate that S-complex asteroids are parent bodies of ordinary chondrite meteorites. Furthermore, they clarify that the space weathering age of the Itokawa surface is of the order of several thousand years. Traditionally, Q-type asteroids have been considered fresh-surfaced. However, as the space weathering timescale is approximately three orders of magnitude lesser than the conventionally considered age, the previously proposed formation mechanisms of Q-type asteroids cannot sufficiently explain the surface refreshening. In this study, we propose a new hypothesis on the surface state of Q-type asteroids: Q-type asteroids have a non-fresh weathered surface with a paucity of fine particles. For verifying this hypothesis, laboratory experiments on the space weathering of ordinary chondrites are performed. Based on the results of these experiments, we found that large (more than 100 {\mu}m) ordinary chondritic particles with space weathering exhibit spectra consistent with Q-type asteroids.

The near-Earth asteroid 3200 Phaethon (1983 TB) is an attractive object not only from a scientific viewpoint but also because of JAXA's DESTINY+ target. The rotational lightcurve and spin properties were investigated based on the data obtained in the ground-based observation campaign of Phaethon. We aim to refine the lightcurves and shape model of Phaethon using all available lightcurve datasets obtained via optical observation, as well as our time-series observation data from the 2017 apparition. Using eight 1-2-m telescopes and an optical imager, we acquired the optical lightcurves and derived the spin parameters of Phaethon. We applied the lightcurve inversion method and SAGE algorithm to deduce the convex and non-convex shape model and pole orientations. We analysed the optical lightcurve of Phaethon and derived a synodic and a sidereal rotational period of 3.6039 h, with an axis ratio of a/b = 1.07. The ecliptic longitude (lambda) and latitude (beta) of the pole orientation were determined as (308, -52) and (322, -40) via two independent methods. A non-convex model from the SAGE method, which exhibits a concavity feature, is also presented.

We present a visible multi-band photometry of trans-Neptunian objects (TNOs) observed by the Subaru Telescope in the framework of Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP) from March in 2014 to September in 2016. We measured the five broad-band (g, r, i, z, and Y) colors over the wavelength range from 0.4 um to 1.0 um for 30 known TNOs using the HSC-SSP survey data covering ~500 deg2 of sky within +/-30 deg of ecliptic latitude. This dataset allows us to characterize the dynamical classes based on visible reflectance spectra as well as to examine the relationship between colors and the other parameters such as orbital elements. Our results show that the hot classical and scattered populations share similar color distributions, while the cold classical population has a reflective decrease toward shorter wavelength below the i band. Based on the obtained color properties, we found that the TNO sample examined in the present work can be separated into two groups by inclination (I), the low-I population consisting of cold classical objects and high-I population consisting of hot classical and scattered objects. The whole sample exhibits an anti-correlation between colors and inclination, but no significant correlation between colors and semi-major axis, perihelion distance, eccentricity, or absolute magnitude. The color-inclination correlation does not seem to be continuous over the entire inclination range. Rather, it is seen only in the high-I population. We found that the low- and high-I populations are distinguishable in the g-i vs. eccentricity plot, but four high-I objects show g-i colors similar to those of the low-I population. If we exclude these four objects, the high-I objects show a positive correlation between g-i and eccentricity and a negative correlation between g-i and inclination with high significance levels.

The formation of the four terrestrial planets of the solar system is one of the most fundamental problems in the planetary sciences. However, the formation of Mercury remains poorly understood. We investigated terrestrial planet formation by performing 110 high-resolution N-body simulation runs using more than 100 embryos and 6000 disk planetesimals representing a primordial protoplanetary disk. To investigate the formation of Mercury, these simulations considered an inner region of the disk at 0.2-0.5 au (the Mercury region) and disks with and without mass enhancements beyond the ice line location, aIL, in the disk, where aIL = 1.5, 2.25, and 3.0 au were tested. Although Venus and Earth analogs (considering both orbits and masses) successfully formed in the majority of the runs, Mercury analogs were obtained in only nine runs. Mars analogs were also similarly scarce. Our Mercury analogs concentrated at orbits with a ~ 0.27-0.34 au, relatively small eccentricities/inclinations, and median mass m ~ 0.2 Earth masses. In addition, we found that our Mercury analogs acquired most of their final masses from embryos/planetesimals initially located between 0.2 and ~1-1.5 au within 10 Myr, while the remaining mass came from a wider region up to ~3 au at later times. Although the ice line was negligible in the formation of planets located in the Mercury region, it enriched all terrestrial planets with water. Indeed, Mercury analogs showed a wide range of water mass fractions at the end of terrestrial planet formation.

Apollo-type NEA (3200) Phaethon, classified at the B/F-type taxonomy, probably the main mass of the Phaethon-Geminid stream complex (PGC), can be the most metamorphic C-complex asteroid in our solar system, since it is heated up to ~1000 K by the solar radiation around its perihelion passages. Hence, its surface material may be easily decomposed in near-sun environment. Phaethon's spectrum exhibits extremely blue-slope in the VIS-NIR region (so-called Phaethon Blue). Another candidate large member of the PGC, (155140) 2005 UD, shows a B/F-type color, however with a C-type-like red color over its ~1/4 rotational part, which implies an exposition of less metamorphosed primordial internal structure of the PGC precursor by a splitting or breakup event long ago. If so, some rotational part of Phaethon should show the C-type color as well as 2005 UD. Hence, we carried out the time-series VIS-spectroscopic observations of Phaethon using 1-m telescope in order to detect such a signature. Also, R-band photometries were simultaneously performed in order to complement our spectroscopy. Consequently, we obtained a total of 68 VIS-spectrophotometric data, 78% of which show the B-type blue-color, as against the rest of 22% showing the C-type red-color. We successfully acquired rotationally time-resolved spectroscopic data, of which particular rotational phase shows a red-spectral slope as the C-type color, as 2005 UD does, suggesting longitudinal inhomogeneity on Phaethon's surface. We constrained this C-type red-colored area in the mid-latitude in Phaethon's southern hemisphere based on the rotationally time-resolved spectroscopy along with Phaethon's axial rotation state, of which size suggests the impact-induced origin of the PGC. We also surveyed the meteoritic analog of Phaethon's surface blue-color, and found thermally metamorphosed CI/CM chondrites as likely candidates.

The Karin family is a very young asteroid family created by an asteroid breakup 5.8 Myr ago. Since the members of this family probably have not experienced significant orbital or collisional evolution yet, it is possible that they still preserve properties of the original family-forming event in terms of their spin state. As we carried out a series of photometric observations of the Karin family asteroids, here we report an analysis result of lightcurves including the rotation period of eleven members. The mean rotation rate of the Karin family members turned out to be much lower than those of NEAs or smaller MBAs (diameter D<12 km), and even lower than that of larger MBAs (D>130 km). We investigated a correlation between the peak-to-peak variation magnitude reduced to zero solar phase angle and the rotation period of the eleven Karin family asteroids, and found a possible trend that elongated members have lower spin rate, and less elongated members have higher spin rate. However, this has to be confirmed by another series of future observations.

We review previously published and newly obtained crater size-frequency distributions in the inner solar system. These data indicate that the Moon and the terrestrial planets have been bombarded by two populations of objects. Population 1, dominating at early times, had nearly the same size distribution as the present-day asteroid belt, and produced the heavily cratered surfaces with a complex, multi-sloped crater size-frequency distribution. Population 2, dominating since about 3.8-3.7 Ga, has the same size distribution as near-Earth objects (NEOs), had a much lower impact flux, and produced a crater size distribution characterized by a differential -3 single-slope power law in the crater diameter range 0.02 km to 100 km. Taken together with the results from a large body of work on age-dating of lunar and meteorite samples and theoretical work in solar system dynamics, a plausible interpretation of these data is as follows. The NEO population is the source of Population 2 and it has been in near-steady state over the past ~3.7-3.8 gigayears; these objects are derived from the main asteroid belt by size-dependent non-gravitational effects that favor the ejection of smaller asteroids. However, Population 1 were main belt asteroids ejected from their source region in a size-independent manner, possibly by means of gravitational resonance sweeping during giant planet orbit migration; this caused the so-called Late Heavy Bombardment (LHB). The LHB began some time before ~3.9 Ga, peaked and declined rapidly over the next ~100 to 300 megayears, and possibly more slowly from about 3.8-3.7 Ga to ~2 Ga. A third crater population (Population S) consists of secondary impact craters that can dominate the cratering record at small diameters.

We present lightcurve observations and their multiband photometry for near-Earth object (NEO) 2011 XA3. The lightcurve has shown a periodicity of 0.0304 +- 0.0003 day (= 43.8 +- 0.4 min). The fast rotation shows that 2011 XA3 is in a state of tension (i.e. a monolithic asteroid) and cannot be held together by self-gravitation. Moreover, the multiband photometric analysis indicates that the taxonomic class of 2011 XA3 is S-complex, or V-type. Its estimated effective diameter is 225 +- 97 m (S-complex) and 166 +- 63 m (V-type), respectively. Therefore, 2011 XA3 is a candidate for the second-largest fast-rotating monolithic asteroid. Moreover, the orbital parameters of 2011 XA3 are apparently similar to those of NEO (3200) Phaethon but F/B-type. We computed the orbital evolutions of 2011 XA3 and Phaethon. However, the results of the computation and distinct taxonomy indicate that neither of the asteroids is of common origin.

The dependence on an applied electric field of the ionization current produced by an energetic electron stopped in liquid helium can be used to determine the spatial distribution of secondary electrons with respect to their geminate partners. An analytic expression relating the current and distribution is derived. The distribution is found to be non-Gaussian with a long tail at larger distances.

The scintillation properties of liquid helium upon the recoil of a low energy helium atom are discussed in the context of the possible use of this medium as a detector of dark matter. It is found that the prompt scintillation yield in the range of recoil energies from a few keV to 100 keV is somewhat higher than that obtained by a linear extrapolation from the measured yield for an 5 MeV alpha particle. A comparison is made of both the scintillation yield and the charge separation by an electric field for nuclear recoils and for electrons stopped in helium.

The newly formed giant planets may have migrated and crossed a number of mutual mean motion resonances (MMRs) when smaller objects (embryos) were accreting to form the terrestrial planets. We investigated the effects of the planetesimal-driven migration of Jupiter and Saturn, and the influence of their mutual 1:2 MMR crossing on terrestrial planet formation for the first time, by performing N-body simulations. These simulations considered distinct timescales of MMR crossing and planet migration. In total, 68 high-resolution simulation runs using 2000 disk planetesimals were performed, which was a significant improvement on previously published results. Even when the effects of the 1:2 MMR crossing and planet migration were included in the system, Venus and Earth analogs (considering both orbits and masses) successfully formed in several runs. In addition, we found that the orbits of planetesimals beyond a ~1.5-2 AU were dynamically depleted by the strengthened sweeping secular resonances associated with Jupiter's and Saturn's more eccentric orbits (relative to present-day) during planet migration. However, this depletion did not prevent the formation of massive Mars analogs (planets with more than 1.5 times Mars' mass). Although late MMR crossings (at t > 30 Myr) could remove such planets, Mars-like small mass planets survived on overly excited orbits (high e and/or i), or were completely lost in these systems. We conclude that the orbital migration and crossing of the mutual 1:2 MMR of Jupiter and Saturn are unlikely to provide suitable orbital conditions for the formation of solar system terrestrial planets. This suggests that to explain Mars' small mass and the absence of other planets between Mars and Jupiter, the outer asteroid belt must have suffered a severe depletion due to interactions with Jupiter/Saturn, or by an alternative mechanism (e.g., rogue super-Earths).

Tens of thousands of sparse-in-time lightcurves from astrometric projects are publicly available. We investigate these data and use them in the lightcurve inversion method to derive new asteroid models. By having a greater number of models with known physical properties, we can gain a better insight into the nature of individual objects and into the whole asteroid population. We use sparse photometry from selected observatories from the AstDyS database, either alone or in combination with dense lightcurves, to determine new asteroid models by the lightcurve inversion method. We present 80 new asteroid models derived from combined data sets where sparse photometry is taken from the AstDyS database and dense lightcurves are from the Uppsala Asteroid Photometric Catalogue (UAPC) and from several individual observers. For 18 asteroids, we present updated shape solutions based on new photometric data. The addition of new models increases the total number of models derived by the lightcurve inversion method to ~200. We also present a simple statistical analysis of physical properties of asteroids where we look for possible correlations between various physical parameters with an emphasis on the spin vector. We present the observed and de-biased distributions of ecliptic latitudes with respect to different size ranges of asteroids as well as a simple theoretical model of the latitude distribution and then compare its predictions with the observed distributions. From this analysis we find that the latitude distribution of small asteroids (D < 30 km) is clustered towards ecliptic poles and can be explained by the YORP thermal effect while the latitude distribution of larger asteroids (D > 60 km) exhibits an evident excess of prograde rotators, probably of primordial origin.

Recent lunar crater studies have revealed an asymmetric distribution of rayed craters on the lunar surface. The asymmetry is related to the synchronous rotation of the Moon: there is a higher density of rayed craters on the leading hemisphere compared with the trailing hemisphere. Rayed craters represent generally the youngest impacts. The purpose of this paper is to test the hypotheses that (i) the population of Near-Earth asteroids (NEAs) is the source of the impactors that have made the rayed craters, and (ii) that impacts by this projectile population account quantitatively for the observed asymmetry. We carried out numerical simulations of the orbital evolution of a large number of test particles representing NEAs in order to determine directly their impact flux on the Moon. The simulations were done in two stages. In the first stage we obtained encounter statistics of NEAs on the Earth's activity sphere. In the second stage we calculated the direct impact flux of the encountering particles on the surface of the Moon; the latter calculations were confined within the activity sphere of the Earth. A steady-state synthetic population of NEAs was generated from a debiased orbital distribution of the known NEAs. We find that the near-Earth asteroids do have an asymmetry in their impact flux on the Moon: apex-to-antapex ratio of 1.32 +/- 0.01. However, the observed rayed crater distribution's asymmetry is significantly more pronounced: apex-to-antapex ratio of 1.65 +/- 0.16. Our results suggest the existence of an undetected population of slower (low impact velocity) projectiles, such as a population of objects nearly coorbiting with Earth; more observational study of young lunar craters is needed to secure this conclusion.

The quasi-Hilda comets (QHCs), being in unstable 3:2 Jovian mean motion resonance, are considered a major cause of temporary satellite capture (TSC) by Jupiter. Though the QHCs may be escaped Hilda asteroids, their origin and nature have not yet been studied in sufficient detail. Of particular interest are long TSCs/orbiters. Orbiters -- in which at least one full revolution about the planet is completed -- are rare astronomical events; only four have been known to occur in the last several decades. Every case has been associated with a QHC: 82P/Gehrels 3; 111P/Helin-Roman-Crockett; P/1996 R2 (Lagerkvist); and the possibly QHC-derived D/1993 F2 (Shoemaker-Levy 9, SL9). We focus on long TSC/orbiter events involving QHCs and Jupiter. Thus we survey the known QHCs, searching for further long TSCs/orbiters over the past century. First, we confirmed the long TSC/orbiter events of 82P, 111P, and 1996 R2 in order to test our method against previous work, applying a general N-body Newtonian code. We then used the same procedure to survey the remaining known QHCs and search for long TSC/orbiter events. We newly identified another long TSC/orbiter: 147P/Kushida-Muramatsu from 1949 May 14 (+97days/-106 days)--1961 July 15. Our result is verified by integrations of 243 cloned orbits which take account of the present orbital uncertainty of this comet. This is the third long TSC and the fifth orbiter to be found, thus long TSC/orbiter events involving Jupiter have occurred once per decade. Two full revolutions about Jupiter were completed and the capture duration was 12.17 (+0.29/-0.27) years; both these numbers rank 147P as third among long TSC/orbiter events, behind SL9 and 111P. This study also confirms the importance of the QHC region as a dynamical route into and out of Jovian TSC, via the Hill's sphere.

Although it is more complicated to search for near-Earth object (NEO) families than main belt asteroid (MBA) families, since differential orbital evolution within a NEO family can cause current orbital elements to drastically differ from each other, we have found that Apollo asteroids (1566) Icarus and the newly discovered 2007 MK6 are almost certainly related. Specifically, their orbital evolutions show a similar profile, time shifted by only ~1000 yr, based on our time-lag theory. The dynamical relationship between Icarus and 2007 MK6 along with a possible dust band, the Taurid-Perseid meteor swarm, implies the first detection of an asteroidal NEO family, namely the "Icarus asteroid family".

A large disruption in the main asteroid belt can cause a large flux, an "asteroid shower", on the terrestrial planets. We quantitatively examine the hypothesis that such an event was the cause of the lunar late heavy bombardment (LHB). We performed numerical integrations of about 20000 test particles starting in the vicinity of the nu6 secular resonance in the main asteroid belt. The purpose of these integrations is to calculate, for each of the terrestrial planets, the collision probability of asteroids coming from an asteroid break-up event in the inner part of the main belt. Compared with previous studies, we simulate nearly two orders of magnitude larger number of particles, and we include the orbital effects of the eight planets, Mercury to Neptune. We also examined in detail the orbital evolution of asteroid fragments once they enter the Earth's activity sphere, including the effect of the Earth-Moon orbit. We obtained the collision probability, the distributions of impact velocities, impact positions, and impact angles of asteroid fragments on the Moon. The collision probability with the Moon (~0.1%) suggests that a fairly large parent body, 1000-1500 km in diameter, is required if the LHB event is to be ascribed to a single asteroid disruption. An even larger parent body is required for less favorable initial conditions than we investigated here. We conclude that an asteroid disruption event is not a viable explanation for the LHB.

New insights into the history of the inner solar system are derived from the impact cratering record of the Moon, Mars, Venus and Mercury, and from the size distributions of asteroid populations. Old craters from a unique period of heavy bombardment that ended $\sim$3.8 billion years ago were made by asteroids that were dynamically ejected from the main asteroid belt, possibly due to the orbital migration of the giant planets. The impactors of the past $\sim$3.8 billion years have a size distribution quite different from the main belt asteroids, but very similar to the population of near-Earth asteroids.

We have studied the physical and chemical properties of a quiescent starless core L1521E with various molecular lines. It is found that there exists a compact dense core traced by the H^13CO^+, HN^13C, CCS, and HC_3N lines; their distributions have a single peak at the same position. The core radius is as small as 0.031 pc, whereas the H_2 density at the peak position is as high as (1.3-5.6)times10^5 cm^-3. Although the density is high enough to excite the inversion transitions of NH_3, these lines are found to be very faint in L1521E. The distributions of NH_3 and CCS seem to be different from those of well-studied starless cores, L1498 and L1544, where the distribution of CCS shows a shell-like structure while that of NH_3 is concentrated at the center of the core. Abundances of carbon-chain molecules are higher in L1521E than the other dark cloud cores, and especially those of sulfur-bearing molecules C_nS are comparable to the cyanopolyyne peak of TMC-1. Our results suggest that L1521E would be in a very early stage of physical and chemical evolution.

We present observations of the 3P1-3P0 fine structure transition of atomic carbon [CI], the J=3-2 transition of CO, as well as of the J=1-0 transitions of 13CO and C18O toward DR15, an HII region associated with two mid-infrared dark clouds (IRDCs). The 13CO and C18O J=1-0 emissions closely follow the dark patches seen in optical wavelength, showing two self-gravitating molecular cores with masses of 2000 Msun and 900 Msun, respectively, at the positions of the catalogued IRDCs. Our data show a rough spatial correlation between [CI] and 13CO J=1-0. Bright [CI] emission occurs in relatively cold gas behind the molecular cores, neither in highly excited gas traced by CO J=3-2 emission nor in HII region/molecular cloud interface. These results are inconsistent with those predicted by standard photodissociation region (PDR) models, suggesting an origin for interstellar atomic carbon unrelated to photodissociation processes.

Large scale mapping observations of the 3P1-3P0 fine structure transition of atomic carbon (CI, 492 GHz) and the J=3-2 transition of CO (346 GHz) toward the Orion A molecular cloud have been carried out with the Mt. Fuji submillimeter-wave telescope. The observations cover 9 square degrees, and include the Orion nebula M42 and the L1641 dark cloud complex. The CI emission extends over almost the entire region of the Orion A cloud and is surprisingly similar to that of 13CO(J=1-0).The CO(J=3-2) emission shows a more featureless and extended distribution than CI.The CI/CO(J=3-2) integrated intensity ratio shows a spatial gradient running from the north (0.10) to the south (1.2) of the Orion A cloud, which we interpret as a consequence of the temperature gradient. On the other hand, the CI/13CO(J=1-0) intensity ratio shows no systematic gradient. We have found a good correlation between the CI and 13CO(J=1-0) intensities over the Orion A cloud. This result is discussed on the basis of photodissociation region models.