The detection and characterization of potentially habitable exoplanets is one of the chief goals of astrophysics for the coming decades. Imaging in reflected light is well suited for characterizing Earth-like planets, as much can be learned about these planets in this wavelength range (i.e., ~0.3-2 μm). Several studies have been conducted to determine the abilities and limitations of reflectance spectroscopy, but most previous studies assumed a homogeneous atmospheric and surface composition. Here we investigate how heterogeneities in the atmosphere and surface of an Earth-like planet impact retrieval results. We extend the ExoReL retrieval framework to include a step function for retrieving wavelength varying surface albedo. We then use it to retrieve on visible-to-near-infrared spectra of realistic 3D Earth models with different surface features in view and varying cloud types/distributions synthesized with the Planetary Spectrum Generator. Including the ability to fit for wavelength dependent albedo mitigates degeneracies that arise when using 1D models to analyze 3D planets, and we recover an Earth-like planet in all cases. We detect surface albedo steps at ~0.7 and ~1.1 μm des
Characterizing the surface and atmosphere of Earth-like planets in reflected light is a key goal for upcoming direct imaging surveys. NASA's next flagship-class astrophysics mission concept, the Habitable Worlds Observatory (HWO), is a space-based Ultraviolet/Optical/Near-Infrared observatory with a mission design requirement to reach the $10^{-10}$ contrast necessary to characterize Earth-like planets around Sun-like stars. While reflected light from planetary surfaces provides a unique opportunity to constrain the coverage of surface materials and biopigments, detailed predictions of HWO's ability to retrieve surface fractions are necessary but have not been conducted. Here, we model photon-counting noise from astrophysical, instrumental, and post-processing sources for the HWO Exploratory Analytic Case 5 design equipped with a charge-6 vector-vortex coronagraph. By combining our photon-counting noise with five distinct modern Earth models at quadrature, we simulate single-visit HWO observations and perform spectral retrievals using the open-source code $\texttt{POSEIDON}$ to assess our ability to constrain both the surface and atmospheric composition. We find that degeneracies b
The search for potentially habitable exoplanets centers on detecting biosignature molecules in Earth-like atmospheres, which makes it essential to understand their detectability under biologically and geologically influenced conditions. In this study, we model the reflection and thermal emission spectra of such atmospheres across the UV/VIS/NIR and mid-IR regions and simulate their detectability with future mission concepts such as the Habitable Worlds Observatory (HWO) and the Large Interferometer for Exoplanets (LIFE). We employ Numerical Weather Prediction (NWP) model data, based on Earth's atmosphere, to derive temperature pressure profiles and couple them with a 1D photochemical model to assess the detectability of these molecules in Earth analogs located 10 parsecs away. We investigate the dominant reaction pathways and their contributions to the atmospheric composition of an Earth analog, with a focus on how they shape the resulting molecular signatures. We also examine the role of surface boundary conditions, which indirectly trace the effects of biological and geological processes, on the detectability of these molecules using HWO- and LIFE-type mission concepts. Our findi
The geomagnetic field has undergone hundreds of polarity reversals over Earth's history, at a variable pace. In numerical models of Earth's core dynamics, reversals occur with increasing frequency when the convective forcing is increased past a critical level. This transition has previously been related to the influence of inertia in the force balance. Because this force is subdominant in Earth's core, concerns have been raised regarding the geophysical applicability of this paradigm. Reproducing the reversal rate of the past million years also requires forcing conditions that do not guarantee that the rest of the geomagnetic variation spectrum is reproduced. These issues motivate the search for alternative reversal mechanisms. Using a suite of numerical models where buoyancy is provided at the bottom of the core by inner-core freezing, we show that the magnetic dipole amplitude is controlled by the relative strength of subsurface upwellings and horizontal circulation at the core surface. A relative weakening of upwellings brings the system from a stable to a reversing dipole state. This mechanism is purely kinematic because it operates irrespectively of the interior force balance.
Sulfur is a redox active element that may have helped mediate an electron flow that kickstarted life and which presently is an essential element for all life on Earth. Despite current uncertainties in global sulfur fluxes, modeling sulfur's abiotic cycling through Earth's deep history is important for understanding the impact of a planet wide biosphere on sulfur geochemical cycling and availability and vice versa. We present here an open-source, dynamical box model for estimating global sulfur fluxes and concentrations among surface and deep Earth reservoirs over Earth history, allowing tracking and estimation of the sulfur distribution in planetary reservoirs over deep time in the absence of life. While the main model presented here does not take into account the abrupt evolution of redox-shunting biosynthetic pathways such as oxygenic photosynthesis, we also modeled the abiotic sulfur cycle before and after a Great Oxidation Event-like transition on Earth-like planets. Our results suggest a considerably distinct chemical makeup of sulfur content in marine sediments in the absence of life on an Earth-like planet, leading to a marine sediment sulfate content two orders of magnitude
Future telescopes will characterize rocky exoplanets in reflected light, revealing their albedo, which depends on surface, cloud, and atmospheric properties. Identifying these features is crucial for assessing habitability. We present reference spectra and phase curves for an unresolved Earth-like exoplanet in reflected and polarized light, showing how phase- and wavelength-dependent reflectance reveals key planetary properties. Using the 3D radiative transfer code MYSTIC, we enhance surface and cloud modeling with validated, wavelength-dependent albedo maps of Earth's seasonal and spectral features, alongside a novel treatment of sub-grid cloud variability using ERA5 reanalysis data. Our models, incorporating high-resolution 3D cloud structures, show that sub-grid cloud variability reduces total reflectance and increases phase curve variability, especially at large phase angles where ocean glint dominates. We also find that neglecting wavelength-dependent albedo maps overestimates the vegetation red edge in spectra. Comparing an Ocean planet to an Earth-like planet with seasonal cloud variability, we show that polarization is more sensitive than intensity alone in distinguishing b
Barchans are eolian dunes of crescent shape found on Earth, Mars and other celestial bodies. Among the different types of barchan-barchan interaction, there is one, known as chasing, in which the dunes remain close but without touching each other. In this paper, we investigate the origins of this barchan-barchan dune repulsion by carrying out grain-scale numerical computations in which a pair of granular heaps is deformed by the fluid flow into barchan dunes that interact with each other. In our simulations, data such as position, velocity and resultant force are computed for each individual particle at each time step, allowing us to measure details of both the fluid and grains that explain the repulsion. We show the trajectories of grains, time-average resultant forces, and mass balances for each dune, and that the downstream barchan shrinks faster than the upstream one, keeping, thus, a relatively high velocity although in the wake of the upstream barchan. In its turn, this fast shrinkage is caused by the flow disturbance, which induces higher erosion on the downstream barchan and its circumvention by grains leaving the upstream dune. Our results help explaining the mechanisms be
Spectral characterization of near-Earth asteroids (NEAs) has revealed a continuum of space-weathered states for the surfaces of S-complex NEAs, with Q-class NEAs, an S-complex subclass, most closely matching the un-weathered surfaces of ordinary chondrite meteorites. Dynamical calculations of the orbital evolution of S-complex NEAs revealed that Q-class NEAs tend to have close encounters with terrestrial planets, suggesting that planetary tides may play a role in refreshing NEA surfaces. However, the exact physical mechanism(s) that drive resurfacing through tidal encounters and the encounter distance at which these mechanisms are effective, has remained unclear. Through the lens of the upcoming (99942) Apophis encounter with Earth in 2029, we investigate the potential for surface mobilization through tidally-driven seismic shaking over short-timescales during encounter and subsequent surface slope evolution over longer-timescales driven by tumbling. We perform multi-scale numerical modeling and find that the 2029 encounter will induce short-term tidally-driven discrete seismic events that lead to high-frequency (>0.1 Hz) surface accelerations that reach magnitudes similar to Ap
Earth's surface is deficient in available forms of many elements considered limiting for prebiotic chemistry. In contrast, many extraterrestrial rocky objects are rich in these same elements. Limiting prebiotic ingredients may, therefore, have been delivered by exogenous material; however, the mechanisms by which exogeneous material may be reliably and non-destructively supplied to a planetary surface remains unclear. Today, the flux of extraterrestrial matter to Earth is dominated by fine-grained cosmic dust. Although this material is rarely discussed in a prebiotic context due to its delivery over a large surface area, concentrated cosmic dust deposits are known to form on Earth today due to the action of sedimentary processes. Here we combine empirical constraints on dust sedimentation with dynamical simulations of dust formation and planetary accretion to show that localized sedimentary deposits of cosmic dust could have accumulated in arid environments on early Earth, in particular glacial settings that today produce cryoconite sediments. Our results challenge the widely held assumption that cosmic dust is incapable of fertilizing prebiotic chemistry. Cosmic dust deposits may
The Hadean, once thought to be uninhabitable and tumultuous, has more recently been recontextualized as a clement time in which oceans, land, and life likely appeared on Earth. This non-exhaustive chapter follows multiple threads from planet formation to the origin of life. We place significant emphasis on the solar system context for the Earth, the timing and nature of crustal formation and the evolution of the surface and atmosphere. Several scenarios for prebiotic chemistry are also discussed including atmospheric photochemistry, wet-dry and freeze-thaw cycles, and hydrothermal vent systems. We attempt to draw connections between the large-scale, planetary processes and various origin of life pathways to illustrate possible overlaps and correlations. In detail, we conclude with and discuss the "impact of impacts" to show how asteroid and comet impacts during the Hadean may have affected many of these processes and scenarios, from generating land to altering the chemical composition and oxidation state of the early Earth's atmosphere and surface.
The instability with respect to global glaciation is a fundamental property of the climate system caused by the positive ice-albedo feedback. The atmospheric carbon dioxide concentration at which this Snowball bifurcation occurs changes through Earth's history because of the slowly increasing solar luminosity. Quantifying this critical CO$_2$ level is not only interesting from a climate dynamics perspective, but also a prerequisite for understanding past Snowball Earth events as well as the conditions for habitability on Earth and other planets. Earlier studies are limited to investigations with simple climate models for Earth's entire history, or studies of individual time slices carried out with a variety of more complex models and for different boundary conditions, making comparisons and the identification of secular changes difficult. Here we use a coupled climate model of intermediate complexity to trace the Snowball bifurcation of an aquaplanet through Earth's history in one consistent model framework. We find that the critical CO$_2$ concentration decreases more or less logarithmically with increasing solar luminosity until about 1 billion years ago, but drops faster in more
Asteroids smaller than about 100 meters are observed to rotate very fast, with periods often much shorter than the critical limit of 2.2 h. Some of these super-fast rotators can also achieve a very large semi-major axis drift induced by the Yarkovsky effect, that in turn, is determined by internal and surface physical properties. We consider the small super-fast rotating near-Earth asteroid 2016 GE1. This object rotates in just 34 seconds, and a large Yarkovsky effect has been determined from astrometry. Here we aim to constrain the thermal inertia of the surface of this extreme object. We used a recently developed statistical method to determine the thermal properties of near-Earth asteroids. The method is based on the comparison between the observed and the modelled Yarkovsky effect, and the thermal conductivity (inertia) is determined by a Monte Carlo approach. Parameters of the Yarkovsky effect model are either fixed if their uncertainty is negligible, modelled with a Gaussian distribution of the errors if they are measured, or deduced from general properties of the population of near-Earth asteroids when they are unknown. Using a well-established orbit determination procedure,
In this paper, we compare the chemistry and the emission spectra of nitrogen-dominated cool, warm, and hot ultra-short-period (USP) super-Earth atmospheres in and out of chemical equilibrium at various surface pressure scenarios ranging from 0.1 to 10 bar. We link the one-dimensional VULCAN chemical kinetic code, in which thermochemical kinetic and vertical transport and photochemistry are taken into account, to the one-dimensional radiative transfer model, PETITRADTRANS, to predict the emission spectra of these planets. The radiative-convective temperature-pressure profiles were computed with the HELIOS code. Then, using PANDEXO noise simulator, we explore the observability of the differences produced by disequilibrium processes with the JWST. Our grids show how different surface pressures can significantly affect the temperature profiles, the atmospheric abundances, and consequently the emission spectra of these planets. We find that the divergences due to disequilibrium processes would be possible to observe in cooler planets by targeting HCN, C2H4, and CO, and in warmer planets by targeting CH4 with HCN, using the NIRSpec and MIRI LRS JWST instruments. These species are also fo
In the pursuit of lunar exploration and the investigation of water presence on the lunar surface, a comprehensive understanding of plasma-surface interactions is crucial since the regolith's space weathering can create H$_2$O. However, the Moon is in the Earth's magnetotail for nearly 20\% of its orbit, which could affect this water creation on the side facing the Earth if this condition shields it from the solar wind. The objective of this study is to understand how the passage of the Moon in the Earth's magnetotail affects the plasma delivery near the lunar surface. The Particle-In-Cell Electromagnetic (EM) Relativistic Global Model, known as IAPIC, is employed to kinetically simulate the Solar Wind-Magnetosphere-Ionosphere-Moon Coupling. The Earth's magnetotail does not prevent the influx of solar wind ions and ionospheric ions into the solar environment; therefore the space weathering of the regolith is not stopped in these conditions. In addition, the charge separation of solar wind ions and electrons happens is modeled, leading to electric fields and charging of the lunar surface that can be validated by observations. The study of the Sun-Earth-Moon system provides insight in
The ice shelves buttressing the Antarctic ice sheet determine the rate of ice-discharge into the surrounding oceans. The geometry of ice shelves, and hence their buttressing strength, is determined by ice flow as well as by the local surface accumulation and basal melt rates, governed by atmospheric and oceanic conditions. Contemporary methods resolve one of these rates, but typically not both. Moreover, there is little information of how they changed in time. We present a new method to simultaneously infer the surface accumulation and basal melt rates averaged over decadal and centennial timescales. We infer the spatial dependence of these rates along flow line transects using internal stratigraphy observed by radars, using a kinematic forward model of internal stratigraphy. We solve the inverse problem using simulation-based inference (SBI). SBI performs Bayesian inference by training neural networks on simulations of the forward model to approximate the posterior distribution, allowing us to also quantify uncertainties over the inferred parameters. We demonstrate the validity of our method on a synthetic example, and apply it to Ekström Ice Shelf, Antarctica, for which newly acq
Mercury is the closest planet to the Sun, possesses a weak intrinsic magnetic field and has only a very tenuous atmosphere (exosphere). These three conditions result in a direct coupling between the plasma emitted from the Sun (namely the solar wind) and Mercury's surface. The planet's magnetic field leads to a non-trivial pattern of plasma precipitation onto the surface, that is expected to contribute to the alteration of the regolith over geological time scales. The goal of this work is to study the solar wind plasma precipitation onto the surface of Mercury from a geographical perspective, as opposed to the local-time-of-day approach of previous precipitation modeling studies. We employ solar wind precipitation maps for protons and electrons from two fully-kinetic numerical simulations of Mercury's plasma environment. These maps are then integrated over two full Mercury orbits (176 Earth days). We found that the plasma precipitation pattern at the surface is most strongly affected by the upstream solar wind conditions, particularly by the interplanetary magnetic field direction, and less by Mercury's 3:2 spin-orbit resonance. We also found that Mercury's magnetic field is able t
Understanding processes that determine the global circulation of the atmosphere is necessary for long-term weather forecasting and climate studies which are critical for ensuring energy security. Processes in the atmosphere depend on many factors including the rotation of the Earth. The vertical component of the emerging centripetal acceleration is taken into account by introducing a geopotential height, and the influence of the horizontal component is usually not taken into account, although these components are close in magnitude. Here, I describe a mechanism of circulation formation due to a non-hydrostatic pressure gradient resulting from centripetal acceleration, which is created by the rotation of the Earth. The effect is similar to the curvature of a water surface in a rotating vessel. This causes the unevenly heated air to move in a meridional direction. The circulation consists of three oppositely directed flows separated vertically by the polar and tropical tropopause.
A purposely built deep learning algorithm for the Verification of Earth-System ParametERisation (VESPER) is used to assess recent upgrades of the global physiographic datasets underpinning the quality of the Integrated Forecasting System (IFS) of the European Centre for Medium-Range Weather Forecasts (ECMWF), which is used both in numerical weather prediction and climate reanalyses. A neural network regression model is trained to learn the mapping between the surface physiographic dataset plus the meteorology from ERA5, and the MODIS satellite skin temperature observations. Once trained, this tool is applied to rapidly assess the quality of upgrades of the land-surface scheme. Upgrades which improve the prediction accuracy of the machine learning tool indicate a reduction of the errors in the surface fields used as input to the surface parametrisation schemes. Conversely, incorrect specifications of the surface fields decrease the accuracy with which VESPER can make predictions. We apply VESPER to assess the accuracy of recent upgrades of the permanent lake and glaciers covers as well as planned upgrades to represent seasonally varying water bodies (i.e. ephemeral lakes). We show t
We argue that the reflection of relic neutrinos from the surface of the Earth results in a significant local $ν-\barν$ asymmetry, far exceeding the expected primordial lepton asymmetry. The net fractional electron neutrino number $\frac{n_{ν_e}-n_{\barν_e}}{n_{ν_e}}$ is up to $\mathcal{O}(10^5) \sqrt{\frac{m_ν}{0.1~\text{eV}}}$ larger than that implied by the baryon asymmetry. This enhancement is due to the weak 4-Fermi repulsion of the $ν_e$ from ordinary matter which slows down the $ν_e$ near the Earth's surface, and to the resulting evanescent neutrino wave that penetrates below the surface. This repulsion thus creates a net $ν_e$ overdensity in a shell $\sim 7~\text{meters} \sqrt{\frac{0.1~\text{eV}}{m_ν}}$ thick around the Earth's surface. Similarly the repulsion between $\barν_μ$ or $\barν_τ$ and ordinary matter creates an overdensity of $\barν_{μ, τ}$ of similar size. These local enhancements increase the size of $\mathcal{O}(G_F)$ torques of the $CνB$ on spin-polarized matter by a factor of order $10^5$. In addition, they create a gradient of the net neutrino density which naturally provides a way out of the forty-year-old ``no-go'' theorems on the vanishing of $\mathcal{O}
A critical early stage for the origin of life on Earth may have involved the production of hydrogen cyanide (HCN) in a reducing, predominantly H$_2$ atmosphere. HCN is crucial for the origin of life as it is a possible precursor to several biomolecules that make up RNA and proteins including nucleobases, nucleotides, amino acids, and ribose. In this work, we perform an in depth experimental and theoretical investigation of HCN production in reducing atmospheric conditions (89-95% H$_2$) possibly representing the earliest stages of the Hadean eon, ~4.5-4.3 billion years ago. We make use of cold plasma discharges - a laboratory analog to shortwave UV radiation - to simulate HCN production in the upper layers of the atmosphere for CH$_4$ abundances ranging from 0.1-6.5%. We then combine experimental mass spectrum measurements with our theoretical plasma models to estimate the HCN concentrations produced in our experiments. We find that upper atmospheric HCN production scales linearly with CH$_4$ abundance with the relation [HCN] = 0.13 $\pm$ 0.01[CH$_4$]. Concentrations of HCN near the surface of the Hadean Earth are expected to be about 2-3 orders of magnitude lower. The addition of