We develop a new retrieval scheme for obtaining two-dimensional surface maps of exoplanets from scattered light curves. In our scheme, the combination of the L1-norm and Total Squared Variation, which is one of the techniques used in sparse modeling, is adopted to find the optimal map. We apply the new method to simulated scattered light curves of the Earth, and find that the new method provides a better spatial resolution of the reconstructed map than those using Tikhonov regularization. We also apply the new method to observed scattered light curves of the Earth obtained during the two-year DSCOVR/EPIC observations presented by Fan et al. (2019). The method with Tikhonov regularization enables us to resolve North America, Africa, Eurasia, and Antarctica. In addition to that, the sparse modeling identifies South America and Australia, although it fails to find the Antarctica maybe due to low observational weights on the poles. Besides, the proposed method is capable of retrieving maps from noise injected light curves of a hypothetical Earth-like exoplanet at 5 pc with noise level expected from coronagraphic images from a 8-m telescope. We find that the sparse modeling resolves Australia, Afro-Eurasia, North America, and South America using 2-year observation with a time interval of one month. Our study shows that the combination of sparse modeling and multi-epoch observation with 1 day or 5 days per month can be used to identify main features of an Earth analog in future direct imaging missions such as the Large UV/Optical/IR Surveyor (LUVOIR).

Solar filaments are dark structures on the solar disk, with an elongated spine and several barbs extending out from the spine. When appearing above the solar limb, a filament is called a prominence, with several feet extending down to the solar surface. It was generally thought that filament barbs are simply the prominence feet veering away from the spine and down to the solar surface. However, it was recently noticed that there might be another dynamic type of barbs, which were proposed to be due to filament thread longitudinal oscillation. If this is the case, the dynamic barbs would not extend down to the solar surface. With the quadrature observations of a filament barb on 2011 June 5 from the {\it Solar Dynamics Observatory} and the {\it STEREO} satellites, we confirm that the filament barb is due to filament thread longitudinal oscillations. Viewed from the side, the filament barb looks like an appendix along the spine of the prominence, and does not extend down to the solar surface as a foot.

Resolving spatially-varying exoplanet features from single-point light curves is essential for determining whether Earth-like worlds harbor geological features and/or climate systems that influence habitability. To evaluate the feasibility and requirements of this spatial feature resolving problem, we present an analysis of multi-wavelength single-point light curves of Earth, where it plays the role of a proxy exoplanet. Here, ~10,000 DSCOVR/EPIC frames collected over a two-year period were integrated over the Earth's disk to yield a spectrally-dependent point source and analyzed using singular value decomposition. We found that, between the two dominant principal components (PCs), the second PC contains surface-related features of the planet, while the first PC mainly includes cloud information. We present the first two-dimensional (2D) surface map of Earth reconstructed from light curve observations without any assumptions of its spectral properties. This study serves as a baseline for reconstructing the surface features of Earth-like exoplanets in the future.

Cassini/UVIS FUV observations of stellar occultations at Titan are well suited for probing its atmospheric composition and structure. However, due to instrument pointing motion, only five out of tens of observations have been analyzed. We present an innovative retrieval method that corrects for the effect of pointing motion by forward modeling the Cassini/UVIS instrument response function with the pointing motion value obtained from the SPICE C-kernel along the spectral dimension. To illustrate the methodology, an occultation observation made during flyby T52 is analyzed, when the Cassini spacecraft had insufficient attitude control. A high-resolution stellar model and an instrument response simulator that includes the position of the point source on the detector are used for the analysis of the pointing motion. The Markov Chain Monte-Carlo method is used to retrieve the line-of-sight abundance profiles of eleven species (CH4, C2H2, C2H4, C2H6, C4H2, C6H6, HCN, C2N2, HC3N, C6N2 and haze particles) in the spectral vector fitting process. We obtain tight constraints on all of the species aside from C2H6, C2N2 and C6N2, for which we only retrieved upper limits. This is the first time that the T52 occultation was used to derive abundances of major hydrocarbon and nitrile species in Titan's upper and middle atmosphere, as pointing motion prohibited prior analysis. With this new method, nearly all of the occultations obtained over the entire Cassini mission could yield reliable profiles of atmospheric composition, allowing exploration of Titan's upper atmosphere over seasons, latitudes, and longitudes.

Over the course of last decade, the Nice model has dramatically changed our view of the solar system's formation and early evolution. Within the context of this model, a transient period of planet-planet scattering is triggered by gravitational interactions between the giant planets and a massive primordial planetesimal disk, leading to a successful reproduction of the solar system's present-day architecture. In typical realizations of the Nice model, self-gravity of the planetesimal disk is routinely neglected, as it poses a computational bottleneck to the calculations. Recent analyses have shown, however, that a self-gravitating disk can exhibit behavior that is dynamically distinct, and this disparity may have significant implications for the solar system's evolutionary path. In this work, we explore this discrepancy utilizing a large suite of Nice odel simulations with and without a self-gravitating planetesimal disk, taking advantage of the inherently parallel nature of graphic processing units. Our simulations demonstrate that self-consistent modeling of particle interactions does not lead to significantly different final planetary orbits from those obtained within conventional simulations. Moreover, self-gravitating calculations show similar planetesimal evolution to non-self-gravitating numerical experiments after dynamical instability is triggered, suggesting that the orbital clustering observed in the distant Kuiper belt is unlikely to have a self-gravitational origin.

On 2015 November 13, the small artificial object designated WT1190F entered the Earth atmosphere above the Indian Ocean offshore Sri Lanka after being discovered as a possible new asteroid only a few weeks earlier. At ESA's SSA-NEO Coordination Centre we took advantage of this opportunity to organize a ground-based observational campaign, using WT1190F as a test case for a possible similar future event involving a natural asteroidal body.

The New Horizons flyby of Pluto confirmed the existence of hazes in its atmosphere. Observations of a large high- to low- phase brightness ratio, combined with the blue color of the haze, suggest that the haze particles are fractal aggregates, analogous to the photochemical hazes on Titan. Therefore, studying the Pluto hazes can shed light on the similarities and differences between the Pluto and Titan atmospheres. We model the haze distribution using the Community Aerosol and Radiation Model for Atmospheres assuming that the distribution is shaped by sedimentation and coagulation of particles originating from photochemistry. Hazes composed of both purely spherical and purely fractal aggregate particles are considered. Agreement between model results and occultation observations is obtained with aggregate particles when the downward flux of photochemical products is equal to the column-integrated methane destruction rate ~1.2 $\times$ 10$^{-14}$ g cm$^{-2}$ s$^{-1}$, while for spherical particles the mass flux must be 2-3 times greater. This flux is nearly identical to the haze production flux of Titan previously obtained by comparing microphysical model results to Cassini observations. The aggregate particle radius is sensitive to particle charging, and a particle charge to radius ratio of 30 e-/{\mu}m is necessary to produce ~0.1-0.2 {\mu}m aggregates near Pluto's surface, in accordance with forward scattering measurements. Such a particle charge to radius ratio is 2-4 times higher than those previously obtained for Titan. Hazes composed of spheres with the same particle charge to radius ratio have particles that are 4 times smaller. These results further suggest that the haze particles are fractal aggregates. We also consider the effect of condensation of HCN, and C$_{2}$-hydrocarbons on the haze particles, which may play an important role in shaping their distributions.