Source author record

Yuk L. Yung

Yuk L. Yung appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

16works
6topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

16 published item(s)

preprint2022arXiv

A bimodal distribution of haze in Pluto's atmosphere

Pluto, Titan, and Triton make up a unique class of solar system bodies, with icy surfaces and chemically reducing atmospheres rich in organic photochemistry and haze formation. Hazes play important roles in these atmospheres, with physical and chemical processes highly dependent on particle sizes, but the haze size distribution in reducing atmospheres is currently poorly understood. Here we report observational evidence that Pluto's haze particles are bimodally distributed, which successfully reproduces the full phase scattering observations from New Horizons. Combined with previous simulations of Titan's haze, this result suggests that haze particles in reducing atmospheres undergo rapid shape change near pressure levels ~0.5Pa and favors a photochemical rather than a dynamical origin for the formation of Titan's detached haze. It also demonstrates that both oxidizing and reducing atmospheres can produce multi-modal hazes, and encourages reanalysis of observations of hazes on Titan and Triton.

preprint2022arXiv

Seasonal variations of chemical species and haze in Titan's upper atmosphere

Seasonal variation is significant in Titan's atmosphere due to the large change of solar insolation resulting from Titan's 26.7° axial tilt relative to the plane of Saturn's orbit. Here we present an investigation of hydrocarbon and nitrile species in Titan's upper atmosphere at 400-1200 km, which includes the mesosphere and the lower thermosphere, over more than one fourth of Titan's year (2006-2014, LS=318°-60°), using eighteen stellar occultation observations obtained by Cassini/UVIS. Vertical profiles of eight chemical species (CH4, C2H2, C2H4, C2H6, C4H2, C6H6, HCN, HC3N) and haze particles are retrieved from these observations using an instrument forward model, which considers the technical issue of pointing motion. The Markov-chain Monte Carlo (MCMC) algorithm is used to obtain the posterior probability distributions of parameters in the retrieval, which inherently tests the extent to which species profiles can be constrained. The results show that no change of the species profiles is noticeable before the equinox, while the decrease of atmospheric temperature and significant upwelling in the summer hemisphere are found five terrestrial years afterwards. Altitude of the detached haze layer decreases towards the vernal equinox then it disappears, and no reappearance is identified within the time range of our data, which is consistent with observations from Cassini/ISS. This study provides observational constraints on the seasonal change of Titan's upper atmosphere, and suggests further investigations of the atmospheric chemistry and dynamics therein.

preprint2022arXiv

Vertical distribution of cyclopropenylidene and propadiene in the atmosphere of Titan

Titan's atmosphere is a natural laboratory for exploring the photochemical synthesis of organic molecules. Significant recent advances in the study of the atmosphere of Titan include: (a) detection of C$_3$ molecules: C$_3$H$_6$, CH$_2$CCH$_2$, c-C$_3$H$_2$, and (b) retrieval of C$_6$H$_6$, which is formed primarily via C$_3$ chemistry, from Cassini-UVIS data. The detection of $c$-C$_3$H$_2$ is of particular significance since ring molecules are of great astrobiological importance. Using the Caltech/JPL KINETICS code, along with the best available photochemical rate coefficients and parameterized vertical transport, we are able to account for the recent observations. It is significant that ion chemistry, reminiscent of that in the interstellar medium, plays a major role in the production of c-C$_3$H$_2$ above 1000 km.

preprint2021arXiv

Global Climate Model Occultation Lightcurves Tested by August 2018 Ground-Based Stellar Occultation

Pluto's atmospheric profiles (temperature and pressure) have been studied for decades from stellar occultation lightcurves. In this paper, we look at recent Pluto Global Climate Model (GCM) results (3D temperature, pressure, and density fields) from Bertrand et al. (2020) and use the results to generate model observer's plane intensity fields (OPIF) and lightcurves by using a Fourier optics scheme to model light passing through Pluto's atmosphere (Young, 2012). This approach can accommodate arbitrary atmospheric structures and 3D distributions of haze. We compared the GCM model lightcurves with the lightcurves observed during the 15-AUG-2018 Pluto stellar occultation. We find that the climate scenario which best reproduces the observed data includes an N2 ice mid-latitude band in the southern hemisphere. We have also studied different haze and P/T ratio profiles: the haze effectively reduces the central flash strength, and a lower P/T ratio both reduces the central flash strength and incurs anomalies in the shoulders of the central flash.

preprint2020arXiv

A Statistical Estimation of the Occurrence of Extraterrestrial Intelligence in the Milky Way Galaxy

In the field of Astrobiology, the precise location, prevalence and age of potential extraterrestrial intelligence (ETI) have not been explicitly explored. Here, we address these inquiries using an empirical galactic simulation model to analyze the spatial-temporal variations and the prevalence of potential ETI within the Galaxy. This model estimates the occurrence of ETI, providing guidance on where to look for intelligent life in the Search for ETI (SETI) with a set of criteria, including well-established astrophysical properties of the Milky Way. Further, typically overlooked factors such as the process of abiogenesis, different evolutionary timescales and potential self-annihilation are incorporated to explore the growth propensity of ETI. We examine three major parameters: 1) the likelihood rate of abiogenesis (λA); 2) evolutionary timescales (Tevo); and 3) probability of self-annihilation of complex life (Pann). We found Pann to be the most influential parameter determining the quantity and age of galactic intelligent life. Our model simulation also identified a peak location for ETI at an annular region approximately 4 kpc from the Galactic center around 8 billion years (Gyrs), with complex life decreasing temporally and spatially from the peak point, asserting a high likelihood of intelligent life in the galactic inner disk. The simulated age distributions also suggest that most of the intelligent life in our galaxy are young, thus making observation or detection difficult.

preprint2020arXiv

Earth as a Proxy Exoplanet: Deconstructing and Reconstructing Spectrophotometric Light Curves

Point source spectrophotometric ("single-point") light curves of Earth-like planets contain a surprising amount of information about the spatial features of those worlds. Spatially resolving these light curves is important for assessing time-varying surface features and the existence of an atmosphere, which in turn is critical to life on Earth and significant for determining habitability on exoplanets. Given that Earth is the only celestial body confirmed to harbor life, treating it as a proxy exoplanet by analyzing time-resolved spectral images provides a benchmark in the search for habitable exoplanets. The Earth Polychromatic Imaging Camera (EPIC) on the Deep Space Climate Observatory (DSCOVR) provides such an opportunity, with observations of ~5000 full-disk sunlit Earth images each year at ten wavelengths with high temporal frequency. We disk-integrate these spectral images to create single-point light curves and decompose them into principal components (PCs). Using machine learning techniques to relate the PCs to six preselected spatial features, we find that the first and fourth PCs of the single-point light curves, contributing ~83.23% of the light curve variability, contain information about low and high clouds, respectively. Surface information relevant to the contrast between land and ocean reflectance is contained in the second PC, while individual land sub-types are not easily distinguishable (<0.1% total light curve variation). We build an Earth model by systematically altering the spatial features to derive causal relationships to the PCs. This model can serve as a baseline for analyzing Earth-like exoplanets and guide wavelength selection and sampling strategies for future observations.

preprint2016arXiv

Constraints on the Microphysics of Pluto's Photochemical Haze from New Horizons Observations

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-/μm is necessary to produce ~0.1-0.2 μ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.

preprint2016arXiv

Hypotheses for near-surface exchange of methane on Mars

The Curiosity rover recently detected a background of 0.7 ppb and spikes of 7 ppb of methane on Mars. This in situ measurement reorients our understanding of the Martian environment and its potential for life, as the current theories do not entail any geological source or sink of methane that varies sub-annually. In particular, the 10-fold elevation during the southern winter indicates episodic sources of methane that are yet to be discovered. Here we suggest a near-surface reservoir could explain this variability. Using the temperature and humidity measurements from the rover, we find that perchlorate salts in the regolith deliquesce to form liquid solutions, and deliquescence progresses to deeper subsurface in the season of the methane spikes. We therefore formulate the following three testable hypotheses. The first scenario is that the regolith in Gale Crater adsorbs methane when dry and releases this methane to the atmosphere upon deliquescence. The adsorption energy needs to be 36 kJ/mol to explain the magnitude of the methane spikes, higher than existing laboratory measurements. The second scenario is that microorganisms convert organic matter in the soil to methane when they are in liquid solutions. This scenario does not require regolith adsorption, but entails extant life on Mars. The third scenario is that deep subsurface aquifers produce the bursts of methane. Continued in situ measurements of methane and water, as well as laboratory studies of adsorption and deliquescence, will test these hypotheses and inform the existence of the near-surface reservoir and its exchange with the atmosphere.

preprint2016arXiv

Reflected Light Curves, Spherical and Bond Albedos of Jupiter- and Saturn-like Exoplanets

Reflected light curves observed for exoplanets indicate bright clouds at some of them. We estimate how the light curve and total stellar heating of a planet depend on forward and backward scattering in the clouds based on Pioneer and Cassini spacecraft images of Jupiter and Saturn. We fit analytical functions to the local reflected brightnesses of Jupiter and Saturn depending on the planet's phase. These observations cover broad bands at 0.59-0.72 and 0.39-0.5 μm, and narrow bands at 0.938 (atmospheric window), 0.889 (CH4 absorption band), and 0.24-0.28 μm. We simulate the images of the planets with a ray-tracing model, and disk-integrate them to produce the full-orbit light curves. For Jupiter, we also fit the modeled light curves to the observed full-disk brightness. We derive spherical albedos for Jupiter, Saturn, and for planets with Lambertian and Rayleigh-scattering atmospheres. Jupiter-like atmospheres can produce light curves that are a factor of two fainter at half-phase than the Lambertian planet, given the same geometric albedo at transit. The spherical albedo is typically lower than for a Lambertian planet by up to a factor of 1.5. The Lambertian assumption will underestimate the absorption of the stellar light and the equilibrium temperature of the planetary atmosphere. We also compare our light curves with the light curves of solid bodies: the moons Enceladus and Callisto. Their strong backscattering peak within a few degrees of opposition (secondary eclipse) can lead to an even stronger underestimate of the stellar heating.

preprint2015arXiv

A Ground-Based Albedo Upper Limit for HD 189733b from Polarimetry

We present 50 nights of polarimetric observations of HD 189733 in $B$ band using the POLISH2 aperture-integrated polarimeter at the Lick Observatory Shane 3-m telescope. This instrument, commissioned in 2011, is designed to search for Rayleigh scattering from short-period exoplanets due to the polarized nature of scattered light. Since these planets are spatially unresolvable from their host stars, the relative contribution of the planet-to-total system polarization is expected to vary with an amplitude of order 10 parts per million (ppm) over the course of the orbit. Non-zero and also variable at the 10 ppm level, the inherent polarization of the Lick 3-m telescope limits the accuracy of our measurements and currently inhibits conclusive detection of scattered light from this exoplanet. However, the amplitude of observed variability conservatively sets a $3 σ$ upper limit to the planet-induced polarization of the system of 58 ppm in $B$ band, which is consistent with a previous upper limit from the POLISH instrument at the Palomar Observatory 5-m telescope (Wiktorowicz 2009). A physically-motivated Rayleigh scattering model, which includes the depolarizing effects of multiple scattering, is used to conservatively set a $3 σ$ upper limit to the geometric albedo of HD 189733b of $A_g < 0.37$. This value is consistent with the value $A_g = 0.226 \pm 0.091$ derived from occultation observations with HST STIS (Evans et al. 2013), but it is inconsistent with the large $A_g = 0.61 \pm 0.12$ albedo reported by (Berdyugina et al. 2011).

preprint2015arXiv

A Multiple Scattering Polarized Radiative Transfer Model: Application to HD 189733b

We present a multiple scattering vector radiative transfer model which produces disk integrated, full phase polarized light curves for reflected light from an exoplanetary atmosphere. We validate our model against results from published analytical and computational models and discuss a small number of cases relevant to the existing and possible near-future observations of the exoplanet HD 189733b. HD 189733b is arguably the most well observed exoplanet to date and the only exoplanet to be observed in polarized light, yet it is debated if the planet's atmosphere is cloudy or clear. We model reflected light from clear atmospheres with Rayleigh scattering, and cloudy or hazy atmospheres with Mie and fractal aggregate particles. We show that clear and cloudy atmospheres have large differences in polarized light as compared to simple flux measurements, though existing observations are insufficient to make this distinction. Futhermore, we show that atmospheres that are spatially inhomogeneous, such as being partially covered by clouds or hazes, exhibit larger contrasts in polarized light when compared to clear atmospheres. This effect can potentially be used to identify patchy clouds in exoplanets. Given a set of full phase polarimetric measurements, this model can constrain the geometric albedo, properties of scattering particles in the atmosphere and the longitude of the ascending node of the orbit. The model is used to interpret new polarimetric observations of HD 189733b in a companion paper.

preprint2015arXiv

Helium Atmospheres on Warm Neptune- and Sub-Neptune-Sized Exoplanets and Applications to GJ 436 b

Warm Neptune- and sub-Neptune-sized exoplanets in orbits smaller than Mercury's are thought to have experienced extensive atmospheric evolution. Here we propose that a potential outcome of this atmospheric evolution is the formation of helium-dominated atmospheres. The hydrodynamic escape rates of Neptune- and sub-Neptune-sized exoplanets are comparable to the diffusion-limited escape rate of hydrogen, and therefore the escape is heavily affected by diffusive separation between hydrogen and helium. A helium atmosphere can thus be formed -- from a primordial hydrogen-helium atmosphere -- via atmospheric hydrodynamic escape from the planet. The helium atmosphere has very different abundances of major carbon and oxygen species from those of a hydrogen atmosphere, leading to distinctive transmission and thermal emission spectral features. In particular, the hypothesis of a helium-dominated atmosphere can explain the thermal emission spectrum of GJ 436 b, a warm Neptune-sized exoplanet, while also consistent with the transmission spectrum. This model atmosphere contains trace amounts of hydrogen, carbon, and oxygen, with the predominance of CO over CH4 as the main form of carbon. With our atmospheric evolution model, we find that if the mass of the initial atmosphere envelope is 1E-3 planetary mass, hydrodynamic escape can reduce the hydrogen abundance in the atmosphere by several orders of magnitude in ~10 billion years. Observations of exoplanet transits may thus detect signatures of helium atmospheres and probe the evolutionary history of small exoplanets.

preprint2015arXiv

Stability of CO2 Atmospheres on Desiccated M Dwarf Exoplanets

We investigate the chemical stability of CO2-dominated atmospheres of desiccated M dwarf terrestrial exoplanets using a 1-dimensional photochemical model. Around Sun-like stars, CO2 photolysis by Far-UV (FUV) radiation is balanced by recombination reactions that depend on water abundance. Planets orbiting M dwarf stars experience more FUV radiation, and could be depleted in water due to M dwarfs' prolonged, high-luminosity pre-main sequences (Luger & Barnes 2015). We show that, for water-depleted M dwarf terrestrial planets, a catalytic cycle relying on H2O2 photolysis can maintain a CO2 atmosphere. However, this cycle breaks down for atmospheric hydrogen mixing ratios <1 ppm, resulting in ~40% of the atmospheric CO2 being converted to CO and O2 on a time scale of 1 Myr. The increased O2 abundance leads to high O3 concentrations, the photolysis of which forms another CO2-regenerating catalytic cycle. For atmospheres with <0.1 ppm hydrogen, CO2 is produced directly from the recombination of CO and O. These catalytic cycles place an upper limit of ~50% on the amount of CO2 that can be destroyed via photolysis, which is enough to generate Earth-like abundances of (abiotic) O2 and O3. The conditions that lead to such high oxygen levels could be widespread on planets in the habitable zones of M dwarfs. Discrimination between biological and abiotic O2 and O3 in this case can perhaps be accomplished by noting the lack of water features in the reflectance and emission spectra of these planets, which necessitates observations at wavelengths longer than 0.95 microns.

preprint2015arXiv

Tracing the fate of carbon and the atmospheric evolution of Mars

The climate of Mars likely evolved from a warmer, wetter early state to the cold, arid current state. However, no solutions for this evolution have previously been found to satisfy the observed geological features and isotopic measurements of the atmosphere. Here we show that a family of solutions exist, invoking no missing reservoirs or loss processes. Escape of carbon via CO photodissociation and sputtering enriches heavy carbon (13C) in the Martian atmosphere, partially compensated by moderate carbonate precipitation. The current atmospheric 13C/12C and rock and soil carbonate measurements indicate an early atmosphere with a surface pressure <1 bar. Only scenarios with large amounts of carbonate formation in open lakes permit higher values up to 1.8 bar. The evolutionary scenarios are fully testable with data from the MAVEN mission and further studies of the isotopic composition of carbonate in the Martian rock record through time.

preprint2013arXiv

Bimodal Distribution of Sulfuric Acid Aerosols in the Upper Haze of Venus

The upper haze (UH) of Venus is variable on the order of days and it is populated by two particle modes. We use a 1D microphysics and vertical transport model based on the Community Aerosol and Radiation Model for Atmospheres to evaluate whether interaction of upwelled cloud particles and sulfuric acid particles nucleated in situ on meteoric dust are able to generate the two size modes and whether their observed variability are due to cloud top vertical transient winds. Nucleation of photochemically produced sulfuric acid onto polysulfur condensation nuclei generates mode 1 cloud droplets that then diffuse upwards into the UH. Droplets generated in the UH from nucleation of sulfuric acid onto meteoric dust coagulate with the upwelled cloud particles and cannot reproduce the observed bimodal size distribution. The mass transport enabled by cloud top transient winds are able to generate a bimodal size distribution in a time scale consistent with observations. Sedimentation and convection in the middle and lower clouds causes the formation of large mode 2 and mode 3 particles. Evaporation of these particles below the clouds creates a local sulfuric acid vapor maximum that causes upwelling of sulfuric acid back into the clouds. If the polysulfur condensation nuclei are small and their production rate is high, coagulation of small droplets onto larger droplets in the middle cloud may result in sulfuric acid "rain" below the clouds once every few Earth months. Reduction of the polysulfur condensation nuclei production rate destroys this oscillation and reduces the mode 1 particle abundance in the middle cloud by two orders of magnitude, though it better reproduces the sulfur-to-sulfuric-acid mass ratio in the cloud and haze droplets. In general we find satisfactory agreement between our results and observations, though improvements could be made by incorporating sulfur microphysics.

preprint2010arXiv

High Temperature Photochemistry in the Atmosphere of HD189733b

Recent infrared spectroscopy of hot exoplanets is beginning to reveal their atmospheric composition. Deep with in the planetary atmosphere, the composition is controlled by thermochemical equilibrium. Photochemistry becomes important higher in the atmosphere, at levels above ~1 bar. These two chemistries compete between ~1-10 bars in hot Jupiter-like atmospheres, depending on the strength of the eddy mixing and temperature. HD189733b provides an excellent laboratory in which to study the consequences of chemistry of hot atmospheres. The recent spectra of HD189733b and HD209458b contain signatures of CH4, CO2, CO and H2O. Here we identify the primary chemical pathways that govern the abundances of CH4, CO2, CO and H2O in the cases of thermochemical equilibrium chemistry, photochemistry, and their combination. Our results suggest that the abundance of these species can be photochemically enhanced above or below the thermochemical equilibrium value, so some caution must be taken when assuming that an atmosphere is in strict thermochemical equilibrium.