Your new experience awaits. Try the new design now and help us make it even better

REVIEW article

Front. Remote Sens., 08 August 2025

Sec. Multi- and Hyper-Spectral Imaging

Volume 6 - 2025 | https://doi.org/10.3389/frsen.2025.1635015

This article is part of the Research TopicEarth Observations from the Deep Space: 10 Years of the DSCOVR MissionView all 4 articles

The single beacon: progresses in understanding Earth as an exoplanet using DSCOVR/EPIC observations

Xun Jian,Xun Jian1,2Lixiang GuLixiang Gu1Siteng Fan
Siteng Fan1*Stuart J. BartlettStuart J. Bartlett2Jiani YangJiani Yang2Jonathan H. JiangJonathan H. Jiang3Yangcheng LuoYangcheng Luo4Yuk L. Yung,Yuk L. Yung2,3
  • 1Department of Earth and Space Sciences, Southern University of Science and Technology, Shenzhen, Guangdong, China
  • 2Division of Geological and Planetary Sciences, California Institute of Technology, Pasadena, CA, United States
  • 3Jet Propulsion Laboratory, California Institute of Technology, Pasadena, CA, United States
  • 4Laboratoire de Météorologie Dynamique (LMD)/Institut Pierre-Simon Laplace (IPSL), École Normale Supérieure (ENS), Université Paris Sciences et Lettres (PSL), Centre National de la Recherche Scientifique (CNRS), Institut Polytechnique de Paris, Sorbonne Université, Paris, France

Almost 6000 exoplanets have thus far been confirmed, revolutionizing our understanding of planetary habitability. Yet, despite the identification of Earth-like exoplanets, definitive evidence of extraterrestrial life remains elusive. Studying Earth, the only confirmed habitable and inhabited planet, as a proxy exoplanet provides critical insights for interpreting forthcoming exoplanet direct-imaging data. Observations from the Deep Space Climate Observatory/Earth Polychromatic Imaging Camera (DSCOVR/EPIC), located at the first Sun-Earth Lagrangian point (L1), offer a unique opportunity to analyze Earth’s full-disk, single-point multi-spectrum light curves. Here, we review progress that treat EPIC data as if Earth were an unresolved, distant world. These studies reveal information about planetary rotation, cloud patterns, and surface types. Autocorrelation of the time series recovers the 24 h rotation period, while principal component analysis (PCA) highlights the land-ocean spectral contrast, enabling the reconstruction of a coarse two-dimensional surface map. Modeling studies further quantify the contributions of different planetary surfaces and clouds to Earth’s observable brightness, with low-level clouds playing a dominant role. Additionally, the effects of Earth’s atmosphere, particularly within strong oxygen bands, have been simulated and evaluated. The rich temporal–spectral “light-curve complexity” produced by its heterogeneous surface and dynamic atmosphere has emerged as a practical, observation-based metric of habitability. Comparisons with simulations and other solar system planets demonstrate that Earth’s light curves exhibit the highest complexity, underscoring its unique status as the only known habitable and inhabited exoplanet. These findings provide a valuable observational baseline for future exoplanet studies, refining our ability to recognize life-supporting worlds beyond the Solar System.

1 Introduction

Detecting and characterizing habitable planets across interstellar distances is a central goal of exoplanet science and astrobiology (Catling et al., 2018). Since the first exoplanet was discovered (Campbell et al., 1988), nearly 6000 exoplanets have been detected and confirmed (NASA Exoplanet Archive1). A few of these Earth-sized planets such as Proxima Centauri b (Anglada-Escudé et al., 2016), TRAPPIST-1 e, f, and g (Gillon et al., 2017), and LHS 1140b (Dittmann et al., 2017) are located within their habitable zones (HZs; Kasting et al., 1993; Kopparapu et al., 2013). These are postulated to be potentially habitable (Fujii et al., 2018), in the sense that they could support water-based, Earth-like life (Lammer et al., 2019). The search for putative life on exoplanets currently centers on the detection of biosignatures, i.e., gaseous, surface, and temporal indicators (Meadows, 2005; 2008), that may suggest biological activity. Current biosignature detection strategies incorporate multiple approaches, including the detection of specific chemical proxies (e.g., oxygen and methane; Meadows, 2017), atmospheric chemical disequilibrium sustained by biogenic gas fluxes (Wogan and Catling, 2020), and technosignatures such as chlorofluorocarbons that rarely form abiotically (Lustig-Yaeger et al., 2023). However, despite extensive efforts to detect biosignatures on exoplanets, no definitive evidence of extraterrestrial life has been found to date. Earth remains the only known example of a habitable world that sustains life (e.g., Meadows et al., 2018; Robinson and Reinhard, 2018), serving as the single reference (data point) for understanding life’s characteristics (Seager and Bains, 2015). Investigations of Earth’s biosignatures will provide essential guidance for NASA’s Habitable Worlds Observatory (HWO2), which is designed to characterize 25 exo-Earth candidates (EECs) while advancing astrophysical research covering ultraviolet, optical, and near-infrared (UV/O/NIR) wavelengths. Consequently, despite the ongoing challenge of detecting extraterrestrial life, Earth remains an indispensable reference for refining detection strategies and interpreting prospective biosignatures.

Ongoing observations from the James Webb Space Telescope (JWST) are investigating atmospheric retention on Earth-like exoplanets orbiting M-dwarf stars (e.g., Greene et al., 2023; Lustig-Yaeger et al., 2023b). However, providing a detailed atmospheric characterization of these planets with JWST remains challenging (e.g., Morley et al., 2017; Krissansen-Totton et al., 2018). Looking ahead, the next-generation of 40-m ground-based Extremely Large Telescopes (ELTs) will enable direct detection of potentially habitable planets via both thermal emission and reflected light (e.g., Quanz et al., 2015; Kasper et al., 2021). Achieving contrast levels of ∼10−8 at 15 milliarcseconds and ∼10−9 at 100 milliarcseconds with limiting magnitudes of I ∼ 27 (Kasper et al., 2021), instruments like ELT-PCS (Planetary Camera and Spectrograph) will be capable of imaging only a handful of nearby Earth-like planets around M-dwarfs, such as Proxima Centauri (Anglada-Escudé et al., 2016), Barnard’s Star (Ribas et al., 2018), Lalande 21185 (Díaz et al., 2019), and Teegarden’s Star (Zechmeister et al., 2019). In addition, high-eccentricity planets such as Gl 514 b (Damasso et al., 2022) and HD 20794 d (Nari et al., 2025), which may sustain temperate climates during parts of their orbits, could also fall within the contrast detection limits. The HWO concept (Feinberg et al., 2024), building upon the concept of the LUVOIR (Large UV/Optical/IR Surveyor; The LUVOIR Team, 2019) and HabEx (Habitable Exoplanet Observatory; Gaudi et al., 2020) missions, aims to obtain spectra of at least 25 EECs, as the Astro2020 Decadal Survey recommended, to search for biosignatures and enable transformative astrophysics across ultraviolet, optical, and near-infrared wavelengths. Besides, recent advancements in hybrid interferometric telescopes may offer alternative pathways for high-resolution direct imaging. Instruments such as the 74-m Colossus and the 20–25 m ExoLife Finder (ELF) might enable two-dimensional reconstructions of exoplanetary surface features (Kuhn et al., 2014; Moretto et al., 2014; Berdyugina et al., 2018; Berdyugina and Kuhn, 2019). However, no existing or planned observatories can achieve the spatial resolution required for direct exoplanet imaging, thus posing significant challenges for characterizing exoplanetary geophysical processes, geodynamic activity, and surface features (Jiang et al., 2018). Furthermore, characterizing planets’ temporal variations necessitates continuous monitoring through sustained direct imaging processes of the exoplanetary system. Therefore, obtaining time-resolved, single-point photometric is crucial for probing exoplanetary surfaces and dynamics.

In this paper, we present a review of advancements in understanding Earth as a proxy exoplanet using DSCOVR/EPIC observations and exploring their implications for future exoplanet research. The organization of this paper is as follows. In Section 2, we provide a brief review of previous studies that have examined Earth from an exoplanetary perspective. In Section 3, we briefly describe the observational characteristics and technical specifications of the DSCOVR/EPIC instrument, including its spectral coverage, imaging configuration, and data products. In Section 4, we introduce the analysis of disk-integrated multi-spectral EPIC observations, along with the subsequent temporal analysis and two-dimensional map reconstruction. In Section 5, we discuss extended approaches that incorporate physical modeling of Earth’s reflectance using spatially resolved surface and atmospheric data to interpret the principal components of the light curves. In Section 6, we review a framework based on information theory and complexity analysis to evaluate potential agnostic biosignatures. Finally, in Section 7, we summarize the key findings and discuss their implications for future exoplanet detection missions and biosignature interpretation strategies.

2 Earth within the context of exoplanets

Early foundational work (Sagan et al., 1993) established Earth as the prototype for exoplanet biosignature detection, using Galileo spacecraft data to demonstrate how atmospheric oxygen/methane, ocean glint, and the so-called vegetation red edge (VRE) could collectively indicate biological activity. This pioneering Earth-as-an-exoplanet approach was later advanced by the Deep Impact mission (EPOXI; Livengood et al., 2011), which developed rotational mapping techniques using diurnal brightness variations (Cowan et al., 2009; Fujii and Kawahara, 2012). EPOXI observations were also used to validate 3D spectral models (Robinson et al., 2010; 2011) that predicted ocean reflection signatures for Earth-like exoplanets. Nevertheless, these studies were limited by either short observation times or insufficient spatial resolution, therefore lacking the detection and interpretation potential of long-term variations in light curves. In this context, Earth observations from the Deep Space Climate Observatory (DSCOVR) satellite’s Earth Polychromatic Imaging Camera (EPIC3) provide a unique opportunity, offering both high spatial resolution and continuous multi-year temporal coverage at multi-spectral wavelengths for studying Earth’s time-varying spectral characteristics. Through disk-integrated observations of Earth and multi-year temporal sampling, EPIC has provided continuous, multi-band photometric time series data of Earth as a single-point light source, thus enabling the systematic spectral analysis presented in subsequent sections.

Until recently, there was no exoplanet optical imaging mission concept that could resolve Earth-like exoplanets, as doing so would require a telescope or array with an effective aperture of ∼2.3 km (Jiang et al., 2018), well beyond current technological capabilities. EPIC images from the DSCOVR mission provide a unique opportunity to monitor Earth in detail over extended periods. The resulting radiance images (Figure 1) show that the Earth appears brightest in the 443 and 551 nm bands due to the Sun’s peak irradiance as reported by the World Climate Research Program (WCRP). Various surfaces, with different physical and chemical state, as well as their roughness (Henderson-Sellers and Wilson, 1983), reflect the solar radiation differently: reflected light from clouds dominates all 10 wavelengths observed by EPIC; ice and snow generally exhibit relatively high reflectance across the 10 wavelengths (e.g., the ice cap of Antarctic; Shields et al., 2013); vegetation is more discernible in the red and near-infrared bands due to strong internal leaf scattering and minimal absorption in this spectral region (Knipling, 1970; Jacquemoud and Baret, 1990); soil reflectance also depends on wavelength, with higher reflectance in the near-infrared, influenced by factors such as moisture, composition, and surface roughness (Cierniewski and Verbrugghe, 1997); ocean reflectance appears relatively enhanced in the 388 nm and blue–green wavelengths primarily due to low water absorption and the contribution of surface Fresnel reflection (Pope and Fry, 1997). However, ocean does not exhibit distinct spectral peaks in reflectance (Shields et al., 2013), and the observed signal can vary with environmental conditions such as wind speed (Enomoto, 2007). The reflectance differences across wavelengths suggest the potential to identify Earth-like features on exoplanets through spectral signatures.

Figure 1
Ten disk images show Earth's radiance at wavelengths from 317 nanometers to 780 nanometers. Each image uses a color scale from blue to red indicating radiance levels. The scale ranges from zero to 0.30 watts per square meter.

Figure 1. DSCOVR EPIC radiance imagery of the Earth taken at 14:47 UTC, 2016 August 7.

Despite the valuable insights gained from studying Earth as a proxy exoplanet, such approaches are inherently limited by our Earth-centric definition of life (Schwieterman et al., 2018). Consequently, researchers often focus on identifying unambiguous molecular or chemical signals that presuppose a fundamental degree of similarity between alien biochemistry and Terran biochemistry (Bartlett et al., 2022), with further contextual information (e.g., Kiang et al., 2018; Meadows et al., 2018). However, alien life may manifest in molecular or chemical forms distinct from those found on Earth (Bartlett and Wong, 2020). Even though, some fundamental physicochemical constraints may still apply universally. For instance, Vladilo and Hassanali (2018) argue that any biochemistry capable of sustaining life must involve extensive hydrogen bonding, which imposes thermal stability requirements broadly consistent with those observed for terrestrial life. These constraints suggest that while the molecular architecture of alien life may differ, its existence may still be confined to specific environmental conditions. In response, some researchers have proposed alternative and agnostic biosignatures (e.g., Dorn et al., 2011; Johnson et al., 2018; Marshall et al., 2021; Walker et al., 2018). In this sense, some agnostic approaches, which particularly emphasize 1) the complexity of molecular components (Dorn et al., 2011; Marshall et al., 2021) 2) the complexity of chemical reaction networks (Walker et al., 2018), and 3) planetary thermodynamic disequilibria (Krissansen-Totton et al., 2016; 2018; 2022), are motivated by the observation that the primary distinction between living and non-living systems lies in the capacity of life to utilize information to perform complex chemical and physical processes (Baluška and Levin, 2016; Davies and Walker, 2016; Farnsworth et al., 2013; Walker et al., 2016). Therefore, given that the complexity of the living world has generally increased over time, sometimes through striking major transitions, the greater the complexity and potentially the further from equilibrium a molecule or reaction network is observed to be, the higher the probability that it originates from a biotic process (Adami et al., 2000; Lineweaver et al., 2013; Smith and Szathmary, 1997).

Based on the ideas above, new perspectives on habitability propose that future works should instead focus on identifying systems that integrate sources of free energy, a variety of material constituents, and informational complexity (Wong et al., 2022), or alternatively, that habitability is fundamentally a manifestation of the existence of life (Chopra and Lineweaver, 2016; Lenardic and Seales, 2021). Additional approaches rooted in specific entropy metrics have also shown promise for assessing the potential habitability of exoplanets (e.g., Bartlett et al., 2022; Segal et al., 2024).

3 DSCOVR/EPIC

DSCOVR (formerly known as Triana4) was launched on 11 February 2015 and has been located at the first Sun-Earth Lagrange point (L1), approximately 1.5 million km from Earth since 8 June 2015. The collaborative mission between NASA (National Aeronautics and Space Administration), NOAA (National Oceanic and Atmospheric Administration), and the USAF (U.S. Air Force) primarily monitors space weather (e.g., solar winds), which is essential for ensuring the accuracy and lead time of NOAA’s space weather alerts and forecasts. It is also equipped with two NASA Earth-viewing instruments: the National Institute of Standards and Technology Advanced Radiometer (NISTAR5) and EPIC, which peer back at Earth and the entire planet to detect changes in the planet’s albedo, ozone absorption, and clouds (Herman et al., 2018). NISTAR is a cavity radiometer designed to measure the absolute and spectrally integrated irradiance that is reflected and emitted from the entire sunlit face of the Earth in three broad wavelengths as a single pixel. EPIC is an imager that provides global spectral images of the entire sunlit face of Earth with high spatial resolution (about 18 × 18 km2, 10 wavelength bands).

The EPIC instrument consists of a 2048 × 2048 hafnium-coated charge-coupled device (CCD) camera with 12-bit readout electronics (Herman et al., 2018). The images are captured using 10 narrowband filters covering the ultraviolet (UV, 317.5, 325, 340, and 388 nm), visible (443, 551, 680, and 688 nm), and near-infrared (NIR, 764 and 779.5 nm) spectral ranges. Specifically, two of these channels (688 and 764 nm) are within the strongly absorbing oxygen B and A bands (Geogdzhayev and Marshak, 2018), respectively. Since 13 June 2015, EPIC has been capturing one full set of 10 wavelengths every 68–110 min, with exposure times ranging from 22 ms (551 nm) to 654 ms (317.5 nm). The camera has a 0.61° field of view (FOV) and a 1″07 angular sampling resolution and produces 2048 × 2048 images that are downsized to 1024 × 1024 pixels (except at 443 nm) with a nadir spatial resolution of about 18 km pixel-1.The spectral characteristics of the 10 filters, including their bandwidths, peak transmissions, and quantum efficiencies, are summarized in Table 1 (Geogdzhayev and Marshak, 2018; Herman et al., 2018; Jiang et al., 2018).

Table 1
www.frontiersin.org

Table 1. Instrument parameters, calibration factors, and solar irradiance of EPIC wavelength band (adapted from Jiang et al. (2018)).

The EPIC Level 1B (L1B6) data are provided in engineering units of counts per second (counts s-1), which are obtained by reading the CCD electronic signals through a 12-bit analog-to-digital (A/D) converter with a gain of 42 electrons per count and then dividing by the exposure time (Herman et al., 2018). The EPIC L1B images with the initial unit of counts per second (counts s-1) were converted into reflectance (R) and radiance (I) through a calibration process (Jiang et al., 2018), and thus the EPIC radiance imagery of Earth was acquired (Figure 1).

While DSCOVR’s primary mission did not include exoplanet studies, its high-resolution Earth observations still provide a unique opportunity to construct disk-integrated rotational light curves of our planet. These time-resolved photometric measurements serve as an empirical reference for modeling Earth-like exoplanets that would appear as unresolved point sources in direct imaging observations (Jiang et al., 2018). Such datasets are particularly valuable for upcoming next-generation telescopes such as the ELT-PCS (Kasper et al., 2021), which will soon be capable of directly detecting and characterizing exoplanets within the habitable zones.

4 Observation of single-point Earth

4.1 Spectrum and frequency

Disk-integrated brightness of the sunlit portion of Earth (and occasionally the Moon) observed by EPIC serves as an ideal proxy for point-source measurements, analogous to observations of distant exoplanets (Jiang et al., 2018). The 388 nm radiance (Figure 1) is predominantly influenced by highly reflective clouds, with minimal contribution from the darker surface background, thus serving as an effective proxy for cloud distribution. In contrast, the 780 nm channel captures significant signals from both clouds and surface features. Fourier analysis using Lomb–Scargle periodograms (Lomb, 1976; Scargle, 1982; Hans et al., 1999) reveals distinct periodic signals in these channels, capturing the temporal variability driven by Earth’s rotation, orbital configuration, and the spatial distribution of clouds and surface features.

The qualitative results (Jiang et al., 2018) show that the variability of Earth’s light curve is largely influenced by the distribution of landmasses, oceans, and clouds. In the UV and visible bands, the reflected radiance is dominated by clouds and enhanced whenever clouds are present in the EPIC FOV. Surface contributions (i.e., positive landmass and negative ocean signals) are mixed with cloud signals in the visible wavelengths and become more prominent in the near-infrared channels. Earth appears brighter during the Southern Hemisphere (SH) summer due to its orbital eccentricity (∼0.017), which brings it closer to the Sun during that season. Additionally, the prominent feature in the Fourier spectra is the 24 h (1 day) peak (Figure 2), reflecting Earth’s 24 h rotation period. A few notable features with periods shorter than 24 h appear across all wavelengths at 12, 8, and 6 h, with the 12 h peak especially strong in the UV (Figures 2b,d). These signals likely result from periodic patterns of the Earth’s surface and cloud systems, which repeat daily when the Earth rotates. The period of 12 h could be a combined signature from clouds, oceans, and landmass patterns. For example, Asia and North America enter the FOV roughly 12 h apart in local time, and the Western Pacific faces EPIC around 02:00–03:00 UTC, followed by the Atlantic Ocean around 14:00–15:00 UTC. The 8-h and 6-h components may arise from the spatial variation of the cloud system and the continents, such as the ∼8 h difference between western Australia and northeastern Africa, and the ∼6 h interval between South America and southern Africa as they sequentially enter the EPIC FOV.

Figure 2
Time series and power spectrum plots. Panels (a) and (b) show the \(I_{388}\) time series and its power spectrum; (c) and (d) show \(I_{790}\) time series and power spectrum. Panels (e) and (f) illustrate PC1 time series and power spectrum; (g) and (h) depict PC2 time series and power spectrum. Each power spectrum features frequencies in days\(^{-1}\) on a logarithmic scale with power density, showing several peaks.

Figure 2. (a) Disk-averaged radiance spectra of I388 from DSCOVR EPIC, with the envelopes of daily maxima and minima indicated by black lines. (b) Fourier power spectrum of the I388 time series, with black dashed lines marking the frequencies corresponding to annual, semi-annual, diurnal, and semi-diurnal cycles. (cf), as well as (g,h), follow the same format as (a,b), respectively, but correspond to I780, PC1, and PC2. Panels (ad) are modified from Fan et al. (2019) (eh) are adapted from Fan et al. (2019). All panels are based on the same dataset as used in Figure 3 of Fan et al. (2019) and Figure 4 of Jiang et al. (2018).

4.2 Principal component and reflector

A quantitative analysis of the individual contributions from these components is essential for a comprehensive understanding of the light curve. While the spectral properties of known surface types can be estimated under different stellar spectral energy distributions (Shields et al., 2013), the specific surface composition of exoplanets remains unknown. Therefore, to maintain generality in the absence of such information, all types of land surfaces are treated as equivalent and represented with an “average” land surface (Fan et al., 2019). Acknowledging that future exoplanet observations may span different wavelength ranges, the reflectance time series is scaled to have zero mean and unit variance prior to further analysis. This normalization ensures equal weighting across spectral channels, which is particularly important for multivariate techniques such as singular value decomposition (SVD) and Gradient Boosted Regression Trees (GBRT; Friedman, 2001), which are sensitive to differences in data scale. The input dataset consisted of 10 reflectance channels along with 2 fractional coverage labels (land and cloud fraction, Figure 3) at each time step. In contrast, unnormalized reflectance values can be retained to preserve physical brightness information from different wavelengths. Additionally, EPIC Level 2 products, which are derived based on well-characterized properties of Earth (e.g., the distribution of surface types, the composition of the surface and atmosphere, and the spectral signatures of various spatial features), provide an opportunity to construct physically grounded models of planetary reflectance. The spatially resolved EPIC-view land cover (NASA/LARC/SD/ASDC 2017) and cloud data (NASA/LARC/SD/ASDC 2018), classified as the more specific fractional contributions of six distinct spatial features (i.e., ocean, desert, snow/ice, vegetation, and low and high clouds) were incorporated as additional indicators, resulting in a time series of 10 reflectance channels and 6 fractional coverage labels at each time step (Gu et al., 2021).

Figure 3
Line graph showing the relative importance of different elements across ten principal components (PC1 to PC10). Elements include low clouds, ocean, high clouds, vegetation, desert, snow/ice, and references from Fan et al., 2019.

Figure 3. Relative importance of different spatial features for the first ten principal components. The line plot shows importance values (ranging from 0 to 1) of six spatial features: Low Clouds (orange), Ocean (blue), High Clouds (green), Vegetation (purple), Desert (golden yellow), and Snow/Ice (light blue). These values were derived from the Gradient Boosted Regression Trees (GBRT) model (modified from Gu et al. (2021)). The black and gray lines correspond to the Land/Ocean and Cloud model results reported by Fan et al. (2019), respectively; methodological details of the PCA and GBRT methods can be found in Section 3 of that work.

The GBRT results (Figure 3) as well as the eigenvalues of PCs (Figure 4a) suggest that PC1 and PC2 are the most important components with distinct contributions to land and cloud features, while the contributions from other PCs are relatively negligible (Figure 3; Fan et al., 2019). For land fraction, PC2 is dominant with a weight of 0.88, whereas PC1 contributes only 0.03. The variation in land/ocean fraction is thus primarily captured by PC2, showing a strong correlation of 0.91. Both PCs contribute substantially to cloud fraction, with weights of 0.28 for PC1 and 0.49 for PC2. The GBRT results also show that PC1 is strongly associated with low cloud, showing a flat spectral shape from UV to NIR wavelengths, consistent with cloud spectra, and a prominent peak corresponding to low clouds in observed data (∼70% relative importance) (Gu et al., 2021). The time series of PC1 is highly correlated with low cloud fraction (R = −0.8). PC2 captures information about land and ocean surfaces, with a spectrum that increases in reflectance with wavelength, matching their spectral behavior (Figure 4). GBRT shows peak importance for land/ocean labels (>84%), except for snow/ice, which may be due to its relatively low coverage on Earth’s surface (∼10%), confinement to high latitudes (resulting in smaller weight for disk integration), and/or slow temporal change. PC2 is strongly correlated with land/ocean contrast (R = −0.96). PC4 contains high cloud information, showing a flat spectrum with oxygen absorption features (688 and 764 nm). While GBRT indicates high cloud dominance (>71%), with a strong correlation (R = −0.82) between high cloud fraction and PC4, PC3 remains unclear but may relate to atmospheric Rayleigh scattering or ozone absorption due to its steep UV slope. PC5 shows peaks for vegetation/desert labels, suggesting it reflects the vegetation-desert contrast. Its spectrum (Figure 4b) is consistent with the VRE, and in experiments with only vegetation or desert, PC5 decreases significantly, providing weak evidence of its link to vegetation.

Figure 4
Graphical analysis of eigenvalues and principal component (PC) spectra. Panel (a) displays eigenvalue trends versus PC index with multiple color-coded lines representing different data sources: Fan et al. (2019), Gu et al. (2021, 2022) for observations and simulations. Panel (b) shows various PC spectra across wavelengths from 300 to 800 nm, with multiple line graphs indicating different PCs. Horizontal axes represent wavelength or PC index, while vertical axes indicate eigenvalue or spectral values. Various line styles and markers differentiate the data sources and PCs.

Figure 4. (a) Singular values (logarithmic scale) of principal components from multiple datasets: Fan et al. (2019) (orange diamonds), Gu et al. (2021) observations (red triangles) and synthetic data (blue inverted triangles), Gu et al. (2022) simulations (green right-pointing triangles), entire observations (purple circles), and selected observations (orange squares). Subsampling and normalization procedures follow the corresponding dataset specifications. (b) Eigenvectors of the first six principal components from the Gu et al. (2021) observations and simulations, and Gu et al. (2022) observations and simulations. To facilitate comparison, PC3 and PC4 are swapped between observations and simulations, following Gu et al. (2021), Gu et al. (2022). The legend in (b) corresponds exactly to that used in (a).

4.3 Surface mapping

Given the strong correlation between PC2 and surface features (Figure 3), and the comparable importance of PC1 and PC2 for clouds, it is inferred that clouds consist of two components: one independent of surface type, and the other associated with surface features. This interpretation aligns with previous studies that incorporate surface type into cloud albedo calculations (Thompson and Barron, 1981; Biasiotti et al., 2022). Specifically, PC2 exhibits smoother and more consistent daily changes (Figures 2g,h), while PC1 shows greater variability across days (Figures 2e,f). PC2 also tends to peak during summer, when the sunlit side of the Northern Hemisphere, with extensive landmass, is visible to both the Sun and the spacecraft. Its consistent diurnal variation likely reflects the combined effect of Earth’s rotation and the asymmetry in surface distribution, supporting the conclusion that the cloud variability is temporally correlated with surface features.

PC2 primarily captures the surface information about Earth with a strong linear correlation, thus making it possible to reconstruct a two-dimensional (2D) surface map. To facilitate generalization for future observations of Earth-like exoplanets, a set of simplified assumptions was adopted (Fan et al., 2019; Fan and Yung, 2020): the spectral properties of reflective surfaces were treated as unknown, the incoming solar flux was assumed to be uniform and known, and the entire planetary surface was modeled as a Lambertian reflector. Although Earth’s oceans are known to be strongly non-Lambertian, this assumption allows the method to be applied more broadly when surface characteristics are uncertain. The viewing geometry was assumed to be known and derived from DSCOVR navigation data over the 2-year observation period. Under these assumptions, the surface reconstruction becomes a linear regression problem (Fan et al., 2019). The reconstructed map (Figure 5a) shows PC2 values are positively correlated with land fraction. Coastlines are defined using the median PC2 value, aligning with the assumption that the total land coverage is unknown. Compared with the global land/ocean map of Earth (Figure 5b), the reconstructed map successfully captures all major continental features.

Figure 5
(a) Map showing 2D reconstruction with smooth color transitions ranging from red to blue, representing varying data values. (b) Global land/ocean map with distinct red for land and dark blue for ocean, illustrating geographical distribution. Both maps include color scales.

Figure 5. (a) Two-dimensional (2D) surface map of Earth, reconstructed using the time series of the PC2. The black contour line indicates the median value of PC2, which is used to represent coastlines. The regularization parameter of λ = 10−3 was used to construct the map. (b) Global land/ocean map of Earth for comparison (adapted figure from Fan et al. (2019)).

Building upon the framework (Fan et al., 2019; Fan and Yung, 2020), further studies have proposed enhancements to improve mapping fidelity and flexibility. Since global mapping is inherently ill-posed, Aizawa et al. (2020) applied sparse modeling techniques to improve the spatial resolution of inferred maps. Besides, Kawahara (2020) formulated a single inverse problem that unifies spin-orbit tomography and NMF (nonnegative matrix factorization)-based spectral unmixing, enabling the simultaneous retrieval of both spectral and geographical information from disk-integrated exoplanet light curves. Moreover, to account for temporal variability, Kawahara and Masuda (2020) developed a dynamic mapping approach by extending the Bayesian framework of Farr et al. (2018), utilizing analytical solutions to the Bayesian inverse problem with (multivariate) Gaussian priors and their isomorphic representations, thereby enabling the modeling of time-varying geography.

5 Numerical simulation

5.1 Reflecting surfaces

Since exoplanets will remain spatially unresolved point sources for the foreseeable future, quantifying the contribution of each spatial component to the integrated light curve is crucial for characterizing their habitability. Although PC2 is primarily interpreted (Fan et al., 2019) in terms of land and cloud fractions, Earth, as a well-observed and well-understood planet, provides a unique opportunity for more detailed modeling beyond such general correlations. Leveraging the availability of high-resolution data on Earth’s surface composition and atmospheric constituents, it becomes feasible to simulate its reflective signatures through physical models that account for diverse surface types, cloud structures, and atmospheric scattering effects.

As discussed above, the EPIC Level 2 products provide the opportunity to clarify 6 spatial features. By using a two-dimensional moving average applied to the pixel-level reflectance data within the phase space defined by the solar and viewing zenith angles, the reflected intensity for each spatial feature at a given wavelength is calculated. This enables the construction of a temporally varying Earth image model, which can then be used to reconstruct DSCOVR observations, analyze their relationship with the spatial features, and interpret the physical meaning of the principal components (Gu et al., 2021).

Assumed that the entire sunlit disk is uniformly covered by each of the six spatial features individually, the corresponding disk-averaged reflection spectra for a synthetic Earth (Figures 6a,b) provides characteristic spectral profiles associated with reflecting surfaces (Gu et al., 2021). For each time step, the reflectance contribution of a given spatial feature was calculated by dividing its total reflectance by the disk-integrated reflectance and a daily average of these contributions was then applied (Figure 6c). The slopes in the first three UV channels result from Rayleigh scattering and ozone absorption. Absorption features from the oxygen A and B bands at 764 nm and 688 nm influence the spectral shape, leading to minimal spectral differences among spatial features. Oxygen absorption is weaker in high clouds compared to low clouds (or the surface) due to shorter optical path lengths, suggesting the possibility of distinguishing high and low clouds in the disk-integrated spectrum. The spectral shape of ice/snow reflectance closely matches that of low clouds, making it difficult to identify their contribution to individual light curves due to smaller coverage. Vegetation is the only feature showing a monotonic increase in reflectance at the longest wavelengths, consistent with the VRE (i.e., 680, 688, and 764 nm). In the first four UV channels (Figure 6c), ocean and cloud contributions dominate (around 40%), with atmospheric scattering and ozone absorption influencing the reflectance (Jiang et al., 2018), making it hard to distinguish surface features. At longer wavelengths, low clouds contribute up to ∼60%, while high clouds are more significant at 688 and 764 nm due to stronger oxygen absorption in the lower atmosphere. Vegetation shows higher contributions in the NIR channels, consistent with the VRE. The seasonal cycle, with more vegetation and less ocean in the northern summer, is also apparent in the time series.

Figure 6
Charts of reflectance data across different wavelengths from 300 to 800 nm. (a) Shows disk-integrated reflectance spectra with varying patterns for vegetation, desert, snow/ice, ocean, high clouds, low clouds, and Earth mean. (b) Displays spectral deviation from Earth mean. (c) Depicts temporal contributions for different wavelengths, highlighting proportions of ocean, desert, vegetation, snow/ice, high clouds, and low clouds over months.

Figure 6. (a) Disk-integrated, full-phase reflection spectra of six spatial features (colored lines) compared with that of the globally averaged Earth (black line). (b) Spectral differences between each spatial feature and the global average, based on (a). (c) Daily averaged time series of fractional reflectance contributions from the six spatial features (adapted from Gu et al. (2021)).

The eigenvalues (Figure 4a) of PC1 and PC2 indicate that they together explain over ∼97% of the total variance, with values nearly an order of magnitude larger than that of PC3 (Gu et al., 2021). PC3–6 form a secondary group, potentially capturing weaker yet meaningful spatial signals. Components beyond PC6, each contributing less than ∼0.02% of the variance, are treated as noise and excluded from further analysis. Despite slightly lower eigenvalues, the synthetic single-point light curves successfully reproduce all characteristic spectral features, except for minor switching between the observed and reconstructed eigenvectors for PC3 and PC4 (Figure 4b). The result is promising, especially considering that the reconstruction was achieved using only six spatially distinct yet temporally invariant surface and atmospheric components.

5.2 Atmosphere

Although a non–radiative transfer (non-RT) Earth model was employed (Gu et al., 2021) and successfully demonstrated that information related to low clouds, surface type distribution, and high clouds is encoded in different PCs of Earth’s light curves, the absence of radiative transfer processes in such models still limited their ability to assess the influence of atmospheric composition and scattering. Previous studies (e.g., Tinetti et al., 2006a; Tinetti et al., 2006b; Robinson et al., 2011; Feng et al., 2018; Batalha et al., 2019) have attempted to model Earth’s light curves using RT models, yet these efforts lack validation against observations extending beyond a few days. Consequently, the development of a comprehensive RT model of Earth, rigorously tested against long-term and diverse observations, remains both critical and pressing.

To investigate the influence of atmospheric radiative transfer on Earth’s disk-integrated reflectance, the Earth Spectrum Simulator (ESS) is developed (Gu et al., 2022). The ESS model incorporates the two-stream-exact-single-scattering (2S-ESS) line-by-line RT model that combines precise single-scattering calculations with a two-stream approximation for multiple scattering and offers a balance between computational efficiency and accuracy (Spurr and Natraj, 2011) and has been widely applied in Earth remote sensing of greenhouse gases and aerosols (Xi et al., 2015; Zhang et al., 2015; 2016; Zeng et al., 2017; 2018; Huang et al., 2020; Natraj et al., 2021), thus accounting for surface reflection, cloud and aerosol scattering, and gaseous absorption.

In the 2S-ESS model, surface reflectance was simulated through BRDFs (Bidirectional Reflectance Distribution Functions, Figure 7a), and then land surfaces were categorized into seven types (NASA/LARC/SD/ASDC, 2017). For snow- and ice-covered pixels, a Lambertian surface albedo was assumed. Additionally, the Rayleigh scattering and molecular absorption were simulated considering only the influence of O2, H2O, and O3, due to their dominant influence in the narrow spectral bands of EPIC. The vertical water vapor and temperature were assumed spatially invariant while surface pressure varied spatially and temporally. The surface pressure (Figure 7c) data was derived from MERRA-2 (Modern-Era Retrospective analysis for Research and Application, Version 2). The 3D ozone distributions were derived (Figure 7d), and all atmospheric fields were sampled on the 5th, 15th, and 25th days of each month to reduce computational demands. Aerosol optical properties (Figure 7e) were calculated using Mie scattering theory, except for dust, where nonsphericity is considered (Meng et al., 2010; Crisp et al., 2021). Clouds (Figure 7f) were further classified as either liquid or ice, with one type per pixel and fixed vertical distributions and their optical properties were based on Mie theory (Hess et al., 1998) for liquid clouds and a cirrus habit mixture (Baum et al., 2011) for ice clouds. For pixels with mixed cloud types at degraded resolution, multiple RT simulations were performed and averaged based on cloud fraction to obtain the final reflectance.

Figure 7
Six globe images demonstrate different environmental data: (a) Land cover with diverse colors indicating vegetation density. (b) Ocean productivity with gradients in color. (c) Surface pressure displaying intense reds. (d) Ozone levels with variations across regions. (e) Aerosol distribution in vivid patterns. (f) Cloud cover with intricate textures. Each image includes a color bar legend on the side.

Figure 7. (a) Land reflectance map calculated using the Ross–Thick–Li–Sparse model (Roujean et al., 1992). Model coefficients for the 443, 551, 680, and 780 nm channels were obtained from the DSCOVR MAIAC Version 02 dataset (NASA/LARC/SD/ASDC, 2018) and the viewing geometry information is provided by the EPIC Level 1B dataset (NASA/LARC/SD/ASDC, DSCOVR EPIC Level 1B Version 3). (b) Ocean reflectance map from the Cox–Munk model (Cox and Munk, 1954). (c) Surface pressure map from MERRA-2 reanalysis (Global Modeling and Assimilation Office, 2015b). (d) Column-integrated ozone map generated following the approach described in Crisp et al. (2021). (e) Total aerosol optical depth. The mixing ratios of dust, sulfate, sea salt, black carbon, and organic carbon are extracted from MERRA-2 (Global Modeling and Assimilation Office, 2015a). (f) Total cloud optical depth. Cloud optical depth, effective pressure, and effective radius (for ice clouds) are obtained from EPIC composite data (NASA/LARC/SD/ASDC, 2017) (adapted from Gu et al. (2022)).

Generally, the ESS can reproduce the major features observed accurately. In the shortest wavelength channels, strong Rayleigh scattering results in generally bright images, with the 340 nm channel being the brightest (Figure 8), which is consistent with the brightness variations caused by the enhanced ozone absorption between 340 nm and 317 nm. In the longer wavelength channels, surface information becomes more clearly visible, except for the two oxygen absorption channels (688 nm and 764 nm). In particular, due to the VRE, the reflectance of vegetated pixels increases significantly in the 780 nm channel, forming a sharp contrast with desert pixels, whose reflectance gradually increases across the 551 nm, 680 nm, and 780 nm channels. These spectral differences among different surface types benefit from the detailed modeling of land surface BRDF. In terms of atmospheric scattering features, the use of MAIAC data enables effective simulation of the backscattering peak observed under EPIC viewing geometry; meanwhile, the adoption of the Cox–Munk surface model allows the simulation to capture the feature of ocean glint (Figure 7b).

Figure 8
Ten global maps displaying various intensity levels at different nanometers (317nm to 780nm) using a color scale. The scale ranges from blue (low intensity) to red (high intensity), representing data over the Earth's surface, particularly over oceans and landmasses.

Figure 8. EPIC synthetic reflectance images in 10 channels at 14:47 UTC, 2016 August 7 (adapted from Gu et al. (2022)).

Specifically, the ESS model successfully reproduces the time-averaged spectra with errors generally below 5%, and within 7% overall, except at the two oxygen absorption bands (688 nm and 764 nm), where simulated reflectance exceeds observations by up to 19.53%. The model exhibits slightly larger temporal variability than observations. In the 317–443 nm range, the simulated reflectance tends to be slightly darker than observed, while simulations at longer wavelengths generally appear brighter. A pronounced seasonal cycle is captured in the Southern Hemisphere, with higher reflectance during summer, likely due to the closer Earth–Sun distance and the presence of the Antarctic ice sheet. A smaller enhancement is also observed in the Northern Hemisphere summer, which may be attributed to increased land surface reflectance and cloud cover driven by the South Asian monsoon (Jiang et al., 2018).

The comparison between the simulated and observed PCs (Figure 4) includes an additional subsampled dataset of approximately 300 time points, selected to align with the simulation output. This subset captures the main temporal variability of the light curves, as indicated by the relatively small amplitudes of higher-order singular values (Figure 4a; Gu et al., 2021). At the same time points, the singular values from the simulations are generally larger than those from the observations, suggesting stronger temporal variability in the simulated light curves. Notably, PC3 and PC4 appear to be interchanged between the simulation and observation (Figure 4b). The simulation successfully reconstructs PC1-5; in particular, the first two PCs show near-perfect agreement, demonstrating that the ESS model accurately reproduces the two most dominant signals, namely, low cloud and surface reflectance patterns, which is in line with the previous findings (Gu et al., 2021). The observed PC4 and PC3 also correspond well to the simulated PC3 and PC4, respectively. Since observed PC4 is primarily associated with high clouds (Gu et al., 2021), and ESS tends to simulate higher reflectance for high clouds, the corresponding eigenvalue increases due to the enhanced variance contribution from high cloud coverage.

Although the details of cloud features show some differences, mainly because the simulation uses the DSCOVR composite cloud dataset, which does not exactly match the instantaneous cloud distribution observed by EPIC, the ESS model still serves as a valuable tool for generating light curves of Earth-like exoplanets and provides an important platform for combining with global climate models (GCMs) to support exoplanet characterization studies with its high fidelity, robustness, and efficiency.

6 Habitability metric

6.1 Complexity

In the context of measuring exoplanet complexity, two key aspects must be considered: 1) how to formally quantify complexity, and 2) how to use remotely detectable signals to assess the complexity of a given exoplanet. To address the first, the formal epsilon machine (EM) framework has been in development for several decades (Crutchfield, 2012), and was applied to the challenge of exoplanet characterization (Bartlett et al., 2022). EMs are discrete mathematical objects (finite state machines) that can compactly represent both simple random and simple deterministic processes (as contrasted with minimal Turing machine programs, which ascribe high complexities to simple random processes). They also have the capacity to reveal the underlying phase-space structure of nonlinear systems (Brodu, 2011; Packard et al., 1980; Crutchfield et al., 1986). EMs are particularly well suited to the analysis of planetary light curves, as they are designed to fulfill several key criteria: 1) reproduce deterministic features of the input time series; 2) reproduce random stochastic features of the input time series; 3) serve as the optimal predictor of the generating process; and 4) represent the internal structure of the generating process in the most compact form (Brodu, 2011; Crutchfield, 2012). These kinds of algorithms take a time series as input and produce a minimal optimized model (the EM) capable of reproducing time series that are, subject to parameter tuning and data quality, etc., statistically equivalent to the original. The size of the resulting model, particularly the information content of its state-space distribution, defines a metric known as statistical complexity.

To address the second challenge, it is hypothesized that planets exhibiting the highest complexity are the most likely to host biospheres, with Earth serving as the initial reference. Even if lifeless planets with high complexity are identified, the approach remains valuable as a characterization tool, given the largely unknown space of exoplanet types, which could vary widely in complexity. For example, the Moon exhibits very low complexity, while the Sun and Jupiter display somewhat richer structures and variability. Earth, with its active atmosphere, hydrosphere, and biosphere, appears to be the most complex body in the Solar System (Kasting and Siefert, 2002; Dyke et al., 2011; Olejarz et al., 2021), although this has yet to be quantitatively confirmed. By using EMs derived from reflectance time series, exoplanet complexity levels can be estimated, which could serve as a potential biosignature.

6.2 Shannon entropy

As an initial statistical assessment, the Shannon entropy was calculated for each proxy exoplanet reflectance time series (Bartlett et al., 2022). The reflectance values were discretized and normalized, and the entropies were computed using the standard formula based on the frequency of occurrence for each reflectance value. These entropies were then averaged over the ten EPIC wavelength channels. On the complexity-entropy plane (Figure 9a), groups are observed along a monotonic slope of increasing complexity and entropy despite inherent noise and data sampling issues. The simplest synthetic Earths (cloudless, mono-surface worlds) exhibit the smallest Shannon entropy, while the highest entropy is observed in the most qualitatively and quantitatively complex case (the original Earth). Within the cloudy worlds, there is a relatively small distinction between the mono-surface and multi-surface groups, but it is promising to observe any distinction despite the cloud cover’s masking effect.

Figure 9
Four scatter plots labeled (a), (b), (c), and (d) display variations in statistical complexity versus entropy metrics. Panel (a) compares Shannon entropy to statistical complexity across different environments and conditions, denoted by various shapes and colors. Panel (b) illustrates Shannon entropy against statistical complexity for Earth and Jupiter at different wavelengths. Panel (c) shows average differential entropy versus average statistical complexity, with distinctions in surface types and cloud presence. Panel (d) plots joint differential entropy of eigencolors versus average statistical complexity, with markers representing diverse environmental conditions and cloud statuses.

Figure 9. (a) Statistical complexity versus Shannon entropy averaged across all wavelength bands for the original and synthetic Earth data, as analyzed in Bartlett et al. (2022) (b) Complexity–entropy comparison between Earth and Jupiter time series after processing. The Jupiter datasets were obtained from the Cassini mission, which collected data over approximately 44 Earth days (∼107 Jupiter days) across nine wavelength channels. Three of these channels were selected for comparison with Earth data: (1) 443 nm (Earth) and 450.9 nm (Jupiter), (2) 551 nm (Earth) and 568.2 nm (Jupiter), and (3) 764 nm (Earth) and 750.5 nm (Jupiter). (c) Statistical complexity versus differential entropy. Consistent with the findings of Bartlett et al. (2022), cloudy worlds exhibit the highest values. A relatively small distinction is observed between the mono-surface and multi-surface groups within the cloudy category. (d) Statistical complexity as a function of joint differential entropy of eigencolors. (a,b) are modified from Bartlett et al. (2022), and (c,d) are modified from Segal et al. (2024).

To further explore the potential for distinguishing living from non-living planets, the gas giant, Jupiter, which is generally assumed to be sterile and known for its complex atmosphere that exhibits both periodic and stochastic features, was compared with Earth (Bartlett et al., 2022; Figure 9b). For each wavelength considered, Earth exhibited higher complexity and entropy than Jupiter, with an average difference of approximately 1.25 bits for complexity and 1 bit for entropy, although these values vary depending on the wavelength. Even though the lower image resolution of the Cassini data may have affected the information content of the Jupiter data, these results still suggest that Earth’s planetary complexity exceeds that of Jupiter, which is consistent with the hypothesis that planetary complexity correlates with the presence of life. However, further investigation is required to fully understand the causal influences leading to this observed correlation.

6.3 Joint differential entropy

When analyzing an object using discrete data, it’s possible that important details in the original measurements are lost during the discretization (binning) step. In contrast, methods applicable to continuous time series data may reveal important, perhaps subtle features that could aid in planetary characterization. This was the motivation for Segal et al. (2024) to employ an approximation for the entropy of a continuous variable to DSCOVR EPIC reflectance time series. Specifically, they used the Vasicek estimate of the differential entropy (Vasicek, 1976). PCA applied to the spatially unresolved multiwavelength reflectance data for each synthetic Earth used by Bartlett et al. (2022) and the average differential entropy (Figure 9c) was derived. This was found to be approximately proportional to the joint differential entropy (JDE). Figure 9d shows that the average JDE of eigencolors, obtained via PCA, reveals a new level of distinctions between the synthetic Earth time series data. Unlike the original wavelength channels, eigenvectors are linearly independent which means that cross-wavelength mutual information, such as that introduced by cloud cover, is only accounted for once (Segal et al., 2024). This is supported by the average Pearson correlation coefficient among the 10 eigencolors, which is on the order of 10–17 for both Earth and the reconstructed worlds. In contrast, when the correlation is computed across the original EPIC narrow band wavelengths, the average value is approximately R = 0.57 for the original Earth data and decreases to R = 0.08 for the cloud-removed version. These results are consistent with previous findings that reflectance variability due to clouds is across all 10 EPIC bands (Jiang et al., 2018).

7 Summary

This paper presents recent milestones in studying Earth as a proxy for habitable exoplanets, utilizing single-point, multi-spectral data from the DSCOVR/EPIC instrument. By integrating observations across ultraviolet, visible, and near-infrared bands, these measurements provide essential insights into the temporal and spectral variability of Earth’s disk-integrated reflectance. Such information is crucial for guiding the interpretation of remote sensing data from spatially unresolved observations of terrestrial exoplanets.

Earth’s photometric signature, when treated as that of a distant point source, captures significant diurnal and seasonal variations. These variations are shaped by the dynamics of atmospheric and surface features. Initial analyses demonstrated that averaging EPIC’s full-disk images could replicate the observational constraints expected in exoplanet missions. Temporal analysis techniques revealed periodic signals tied to planetary rotation, cloud cover, and surface distribution. To extract meaningful structure from these light curves, principal component analysis was employed to identify dominant modes of variability. The leading principal components captured distinct physical processes, including cloud dynamics and land-ocean contrast. These results enabled the reconstruction of a simplified surface map of Earth without predefined surface information, offering a method adaptable to diverse planetary systems.

Further investigations introduced physically based models that quantified the individual roles of atmospheric layers and surface types. The influence of low clouds proved especially strong in the ultraviolet and visible bands, while contributions from vegetation, high-altitude clouds, and oceans became more prominent at longer wavelengths. These studies clarified how different components contribute to the overall photometric signal. To simulate Earth’s reflectance with greater fidelity, a radiative transfer model was developed through the Earth Spectrum Simulator. This tool incorporated surface reflectance properties, atmospheric scattering, and absorption processes and its simulations closely matched observations, providing strong validation of the model’s accuracy and its capacity to replicate key spectral and temporal patterns.

Beyond chemical heuristics for biosignature detection, new methodologies based on information theory and complexity science were described, which offer new means to assess exoplanets and agnostically search for life. These approaches revealed that Earth’s reflectance time series displays significantly higher complexity than our gas giant companion (Jupiter). This complexity is proposed as a potential agnostic biosignature, reflecting dynamic processes that could be causally linked to biological activity.

The combined use of observational analysis, physical modeling, and complexity metrics forms a comprehensive framework for characterizing terrestrial exoplanets, integrating various scientific disciplines to enhance our understanding of diverse planetary systems. By leveraging high-resolution multi-spectral observations, physical models, and advanced techniques for assessing complexity, this research not only refines our ability to study exoplanetary environments but also contributes to the development of diagnostic tools for evaluating habitability and signs of life.

Author contributions

XJ: Writing – review and editing, Writing – original draft. LG: Writing – review and editing, Methodology. SF: Project administration, Writing – review and editing, Methodology. SJB: Methodology, Writing – review and editing. JY: Writing – review and editing. JHJ: Writing – review and editing, Conceptualization. YL: Validation, Writing – review and editing. YY: Supervision, Writing – review and editing.

Funding

The author(s) declare that financial support was received for the research and/or publication of this article. SF acknowledges funding from the Stable Support Plan Program for the Higher Education Institutions of the Shenzhen Science and Technology Innovation Commission through grant No. 20231115103030002. JHJ acknowledged the support from the Jet Propulsion Laboratory, California Institute of Technology, sponsored by NASA. JHJ and SJB acknowledge support from the NASA XRP (Exoplanets Research Program) grant no. 22-XRP22_2-0044.

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Generative AI statement

The author(s) declare that no Generative AI was used in the creation of this manuscript.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Footnotes

1https://exoplanetarchive.ipac.caltech.edu/

2https://science.nasa.gov/astrophysics/programs/habitable-worlds-observatory/

3http://epic.gsfc.nasa.gov

4http://www.nesdis.noaa.gov/DSCOVR/

5http://www.nesdis.noaa.gov/DSCOVR/spacecraft.html

6https://eosweb.larc.nasa.gov/project/dscovr/dscovr_table/

References

Adami, C., Ofria, C., and Collier, T. C. (2000). Evolution of biological complexity. Proc. Natl. Acad. Sci. 97 (9), 4463–4468. doi:10.1073/pnas.97.9.4463

PubMed Abstract | CrossRef Full Text | Google Scholar

Aizawa, M., Kawahara, H., and Fan, S. (2020). Global mapping of an exo-Earth using sparse modeling. Astrophysical J. 896 (1), 22. doi:10.3847/1538-4357/ab8d30

CrossRef Full Text | Google Scholar

Anglada-Escudé, G., Amado, P. J., Barnes, J., Berdiñas, Z. M., Butler, R. P., Coleman, G. A., et al. (2016). A terrestrial planet candidate in a temperate orbit around Proxima Centauri. nature 536 (7617), 437–440. doi:10.1038/nature19106

PubMed Abstract | CrossRef Full Text | Google Scholar

Baluška, F., and Levin, M. (2016). On having no head: cognition throughout biological systems. Front. Psychol. 7, 902. doi:10.3389/fpsyg.2016.00902

PubMed Abstract | CrossRef Full Text | Google Scholar

Bartlett, S., and Wong, M. L. (2020). Defining lyfe in the universe: from three privileged functions to four pillars. Life 10 (4), 42. doi:10.3390/life10040042

PubMed Abstract | CrossRef Full Text | Google Scholar

Bartlett, S., Li, J., Gu, L., Sinapayen, L., Fan, S., Natraj, V., et al. (2022). Assessing planetary complexity and potential agnostic biosignatures using epsilon machines. Nat. Astron. 6 (3), 387–392. doi:10.1038/s41550-021-01559-x

CrossRef Full Text | Google Scholar

Batalha, N. E., Marley, M. S., Lewis, N. K., and Fortney, J. J. (2019). Exoplanet reflected-light spectroscopy with PICASO. Astrophysical J. 878 (1), 70. doi:10.3847/1538-4357/ab1b51

CrossRef Full Text | Google Scholar

Baum, B. A., Yang, P., Heymsfield, A. J., Schmitt, C. G., Xie, Y., Bansemer, A., et al. (2011). Improvements in shortwave bulk scattering and absorption models for the remote sensing of ice clouds. J. Appl. Meteorology Climatol. 50 (5), 1037–1056. doi:10.1175/2010jamc2608.1

CrossRef Full Text | Google Scholar

Berdyugina, S. V., and Kuhn, J. R. (2019). Surface imaging of proxima b and other exoplanets: albedo maps, biosignatures, and technosignatures. Astronomical J. 158 (6), 246. doi:10.3847/1538-3881/ab2df3

CrossRef Full Text | Google Scholar

Berdyugina, S. V., Kuhn, J. R., Langlois, M., Moretto, G., Krissansen-Totton, J., and Catling, D. (2018). “The Exo-Life Finder (ELF) telescope: new strategies for direct detection of exoplanet biosignatures and technosignatures,”Ground-based Airborne Telesc. VII, 10700, 1453–1466. doi:10.1117/12.2313781

CrossRef Full Text | Google Scholar

Biasiotti, L., Simonetti, P., Vladilo, G., Silva, L., Murante, G., Ivanovski, S., et al. (2022). EOS-ESTM: a flexible climate model for habitable exoplanets. Mon. Notices R. Astronomical Soc. 514 (4), 5105–5125. doi:10.1093/mnras/stac1642

CrossRef Full Text | Google Scholar

Brodu, N. (2011). Reconstruction of epsilon-machines in predictive frameworks and decisional states. Adv. Complex Syst. 14 (05), 761–794. doi:10.1142/s0219525911003347

CrossRef Full Text | Google Scholar

Campbell, B., Walker, G. A. H., and Yang, S. (1988). “A search for planetary mass companions to nearby stars,” in Bioastronomy—the next steps: proceedings of the 99th colloquium of the international astronomical union held in balaton (Netherlands: Springer), 83–90. Hungary, June 22-27, 1987.

Google Scholar

Catling, D. C., Krissansen-Totton, J., Kiang, N. Y., Crisp, D., Robinson, T. D., DasSarma, S., et al. (2018). Exoplanet biosignatures: a framework for their assessment. Astrobiology 18 (6), 709–738. doi:10.1089/ast.2017.1737

PubMed Abstract | CrossRef Full Text | Google Scholar

Chopra, A., and Lineweaver, C. H. (2016). The case for a Gaian bottleneck: the biology of habitability. Astrobiology 16 (1), 7–22. doi:10.1089/ast.2015.1387

PubMed Abstract | CrossRef Full Text | Google Scholar

Cierniewski, J., and Verbrugghe, M. (1997). Influence of soil surface roughness on soil bidirectional reflectance. Int. J. remote Sens. 18 (6), 1277–1288. doi:10.1080/014311697218412

CrossRef Full Text | Google Scholar

Cowan, N. B., Agol, E., Meadows, V. S., Robinson, T., Livengood, T. A., Deming, D., et al. (2009). Alien maps of an ocean-bearing world. Astrophysical J. 700 (2), 915–923. doi:10.1088/0004-637x/700/2/915

CrossRef Full Text | Google Scholar

Cox, C., and Munk, W. (1954). Measurement of the roughness of the sea surface from photographs of the sun’s glitter. J. Opt. Soc. Am. 44 (11), 838–850. doi:10.1364/josa.44.000838

CrossRef Full Text | Google Scholar

Crisp, D., O’Dell, C., Eldering, A., Fisher, B., Oyafuso, F., and Payne, V. (2021). Orbiting carbon observatory (OCO-2) level 2 full physics algorithm theoretical basis document, Version 3.0–Rev, 1. Available online at: https://docserver.gesdisc.eosdis.nasa.gov/public/project/OCO/OCO_L2_ATBD.pdf.

Google Scholar

Crutchfield, J. P. (2012). Between order and chaos. Nat. Phys. 8 (1), 17–24. doi:10.1038/nphys2190

CrossRef Full Text | Google Scholar

Crutchfield, J. P., Farmer, J. D., Packard, N. H., and Shaw, R. S. (1986). Chaos. Sci. Am., 255(1), 46–57. doi:10.1038/scientificamerican1286-46

CrossRef Full Text | Google Scholar

Damasso, M., Perger, M., Almenara, J. M., Nardiello, D., Pérez-Torres, M., Sozzetti, A., et al. (2022). A quarter century of spectroscopic monitoring of the nearby M dwarf Gl 514-A super-Earth on an eccentric orbit moving in and out of the habitable zone. Astronomy and Astrophysics 666, A187. doi:10.1051/0004-6361/202243522

CrossRef Full Text | Google Scholar

Davies, P. C., and Walker, S. I. (2016). The hidden simplicity of biology. Rep. Prog. Phys. 79 (10), 102601. doi:10.1088/0034-4885/79/10/102601

PubMed Abstract | CrossRef Full Text | Google Scholar

Díaz, R. F., Delfosse, X., Hobson, M. J., Boisse, I., Astudillo-Defru, N., and Bonfils, X. (2019). The SOPHIE search for northern extrasolar planets-XIV. A temperate (Teq∼ 300 K) super-earth around the nearby star Gliese 411. Astronomy and Astrophysics 625, A17. doi:10.1051/0004-6361/201935019

CrossRef Full Text | Google Scholar

Dittmann, J. A., Irwin, J. M., Charbonneau, D., Bonfils, X., Astudillo-Defru, N., Haywood, R. D., et al. (2017). A temperate rocky super-Earth transiting a nearby cool star. Nature 544 (7650), 333–336. doi:10.1038/nature22055

PubMed Abstract | CrossRef Full Text | Google Scholar

Dorn, E. D., Nealson, K. H., and Adami, C. (2011). Monomer abundance distribution patterns as a universal biosignature: examples from terrestrial and digital life. J. Mol. Evol. 72, 283–295. doi:10.1007/s00239-011-9429-4

PubMed Abstract | CrossRef Full Text | Google Scholar

Dyke, J. G., Gans, F., and Kleidon, A. (2011). Towards understanding how surface life can affect interior geological processes: a non-equilibrium thermodynamics approach. Earth Syst. Dyn. 2 (1), 139–160. doi:10.5194/esd-2-139-2011

CrossRef Full Text | Google Scholar

Enomoto, T. (2007). Ocean surface albedo in AFES. JAMSTEC Rep. Res. Dev., 6, 21–30. doi:10.5918/jamstecr.6.21

CrossRef Full Text | Google Scholar

Fan, S., and Yung, Y. L. (2020). Surface mapping of earth-like exoplanets using single point light curves. J. Vis. Exp. 159, e60951. doi:10.3791/60951

PubMed Abstract | CrossRef Full Text | Google Scholar

Fan, S., Li, C., Li, J. Z., Bartlett, S., Jiang, J. H., and Natraj, V. (2019). Earth as an exoplanet: a two-dimensional alien map. Exoplanets in our backyard: solar system and exoplanet synergies on planetary formation. Evol. Habitability 2195, 3013. doi:10.3847/2041-8213/ab3a49

CrossRef Full Text | Google Scholar

Farnsworth, K. D., Nelson, J., and Gershenson, C. (2013). Living is information processing: from molecules to global systems. Acta biotheor. 61, 203–222. doi:10.1007/s10441-013-9179-3

PubMed Abstract | CrossRef Full Text | Google Scholar

Farr, B., Farr, W. M., Cowan, N. B., Haggard, H. M., and Robinson, T. (2018). Exocartographer: a Bayesian framework for mapping exoplanets in reflected light. Astronomical J. 156 (4), 146. doi:10.3847/1538-3881/aad775

CrossRef Full Text | Google Scholar

Feinberg, L., Ziemer, J., Ansdell, M., Crooke, J., Dressing, C., and Mennesson, B. (2024). “The Habitable Worlds Observatory engineering view: status, plans, and opportunities,”Space Telesc. Instrum. 2024 Opt. Infrared, Millim. Wave, 13092. 511–522. doi:10.1117/12.3018328

CrossRef Full Text | Google Scholar

Feng, Y. K., Robinson, T. D., Fortney, J. J., Lupu, R. E., Marley, M. S., Lewis, N. K., et al. (2018). Characterizing earth analogs in reflected light: atmospheric retrieval studies for future space telescopes. Astronomical J. 155 (5), 200. doi:10.3847/1538-3881/aab95c

CrossRef Full Text | Google Scholar

Friedman, J. H. (2001). Greedy function approximation: a gradient boosting machine. Ann. statistics 29, 1189–1232. doi:10.1214/aos/1013203451

CrossRef Full Text | Google Scholar

Fujii, Y., and Kawahara, H. (2012). Mapping earth analogs from photometric variability: spin–orbit tomography for planets in inclined orbits. Astrophysical J. 755 (2), 101. doi:10.1088/0004-637x/755/2/101

CrossRef Full Text | Google Scholar

Fujii, Y., Angerhausen, D., Deitrick, R., Domagal-Goldman, S., Grenfell, J. L., Hori, Y., et al. (2018). Exoplanet biosignatures: observational prospects. Astrobiology 18 (6), 739–778. doi:10.1089/ast.2017.1733

PubMed Abstract | CrossRef Full Text | Google Scholar

Gaudi, B. S., Seager, S., Mennesson, B., Kiessling, A., Warfield, K., and Cahoy, K. (2020). The Habitable Exoplanet Observatory (HabEx) mission concept study final report. arXiv preprint arXiv:2001.06683.

Google Scholar

Geogdzhayev, I. V., and Marshak, A. (2018). Calibration of the DSCOVR EPIC visible and NIR channels using MODIS Terra and Aqua data and EPIC lunar observations. Atmos. Meas. Tech. 11 (1), 359–368. doi:10.5194/amt-11-359-2018

PubMed Abstract | CrossRef Full Text | Google Scholar

Gillon, M., Triaud, A. H., Demory, B. O., Jehin, E., Agol, E., Deck, K. M., et al. (2017). Seven temperate terrestrial planets around the nearby ultracool dwarf star TRAPPIST-1. Nature 542 (7642), 456–460. doi:10.1038/nature21360

PubMed Abstract | CrossRef Full Text | Google Scholar

Greene, T. P., Bell, T. J., Ducrot, E., Dyrek, A., Lagage, P. O., and Fortney, J. J. (2023). Thermal emission from the Earth-sized exoplanet TRAPPIST-1 b using JWST. Nature 618 (7963), 39–42. doi:10.1038/s41586-023-05951-7

PubMed Abstract | CrossRef Full Text | Google Scholar

Gu, L., Fan, S., Li, J., Bartlett, S. J., Natraj, V., Jiang, J. H., et al. (2021). Earth as a proxy exoplanet: deconstructing and reconstructing spectrophotometric light curves. Astronomical J. 161 (3), 122. doi:10.3847/1538-3881/abd54a

CrossRef Full Text | Google Scholar

Gu, L., Zeng, Z. C., Fan, S., Natraj, V., Jiang, J. H., Crisp, D., et al. (2022). Earth as a proxy exoplanet: simulating DSCOVR/EPIC observations using the Earth spectrum simulator. Astronomical J. 163 (6), 285. doi:10.3847/1538-3881/ac5e2e

CrossRef Full Text | Google Scholar

Hans, P., Van, D., Olofsen, E., VanHartevelt, J. H., and Kruyt, E. W. (1999). Searching for biological rhythms: peak detection in the periodogram of unequally spaced data. J. Biol. rhythms 14 (6), 617–620. doi:10.1177/074873099129000984

CrossRef Full Text | Google Scholar

Henderson-Sellers, A., and Wilson, M. F. (1983). Albedo observations of the earth’s surface for climate research. Philosophical Trans. R. Soc. Lond. Ser. A, Math. Phys. Sci. 309 (1508), 285–294. doi:10.1098/rsta.1983.0042

CrossRef Full Text | Google Scholar

Herman, J., Huang, L., McPeters, R., Ziemke, J., Cede, A., and Blank, K. (2018). Synoptic ozone, cloud reflectivity, and erythemal irradiance from sunrise to sunset for the whole earth as viewed by the DSCOVR spacecraft from the earth–sun Lagrange 1 orbit. Atmos. Meas. Tech. 11 (1), 177–194. doi:10.5194/amt-11-177-2018

CrossRef Full Text | Google Scholar

Hess, M., Koepke, P., and Schult, I. (1998). Optical properties of aerosols and clouds: the software package OPAC. Bull. Am. meteorological Soc. 79 (5), 831–844. doi:10.1175/1520-0477(1998)079<0831:opoaac>2.0.co;2

CrossRef Full Text | Google Scholar

Huang, Y., Natraj, V., Zeng, Z. C., Kopparla, P., and Yung, Y. L. (2020). Quantifying the impact of aerosol scattering on the retrieval of methane from airborne remote sensing measurements. Atmos. Meas. Tech. 13 (12), 6755–6769. doi:10.5194/amt-13-6755-2020

CrossRef Full Text | Google Scholar

Jacquemoud, S., and Baret, F. (1990). PROSPECT: a model of leaf optical properties spectra. Remote Sens. Environ. 34 (2), 75–91. doi:10.1016/0034-4257(90)90100-z

CrossRef Full Text | Google Scholar

Jiang, J. H., Zhai, A. J., Herman, J., Zhai, C., Hu, R., Su, H., et al. (2018). Using deep space climate observatory measurements to study the Earth as an exoplanet. Astronomical J. 156 (1), 26. doi:10.3847/1538-3881/aac6e2

CrossRef Full Text | Google Scholar

Johnson, S. S., Anslyn, E. V., Graham, H. V., Mahaffy, P. R., and Ellington, A. D. (2018). Fingerprinting non-terran biosignatures. Astrobiology 18 (7), 915–922. doi:10.1089/ast.2017.1712

PubMed Abstract | CrossRef Full Text | Google Scholar

Kasper, M., Urra, N. C., Pathak, P., Bonse, M., Nousiainen, J., Engler, B., et al. (2021). PCS--A roadmap for exoearth imaging with the ELT. arXiv Prepr. arXiv:2103.11196. doi:10.48550/arXiv.2103.11196

CrossRef Full Text | Google Scholar

Kasting, J. F., and Siefert, J. L. (2002). Life and the evolution of Earth's atmosphere. Science 296 (5570), 1066–1068. doi:10.1126/science.1071184

PubMed Abstract | CrossRef Full Text | Google Scholar

Kasting, J. F., Whitmire, D. P., and Reynolds, R. T. (1993). Habitable zones around main sequence stars. Icarus 101 (1), 108–128. doi:10.1006/icar.1993.1010

PubMed Abstract | CrossRef Full Text | Google Scholar

Kawahara, H. (2020). Global mapping of the surface composition on an exo-Earth using color variability. Astrophysical J. 894 (1), 58. doi:10.3847/1538-4357/ab87a1

CrossRef Full Text | Google Scholar

Kawahara, H., and Masuda, K. (2020). Bayesian dynamic mapping of an exo-Earth from photometric variability. Astrophysical J. 900 (1), 48. doi:10.3847/1538-4357/aba95e

CrossRef Full Text | Google Scholar

Kiang, N. Y., Domagal-Goldman, S., Parenteau, M. N., Catling, D. C., Fujii, Y., Meadows, V. S., et al. (2018). Exoplanet biosignatures: at the dawn of a new era of planetary observations. Astrobiology 18 (6), 619–629. doi:10.1089/ast.2018.1862

PubMed Abstract | CrossRef Full Text | Google Scholar

Knipling, E. B. (1970). Physical and physiological basis for the reflectance of visible and near-infrared radiation from vegetation. Remote Sens. Environ. 1 (3), 155–159. doi:10.1016/s0034-4257(70)80021-9

CrossRef Full Text | Google Scholar

Kopparapu, R. K., Ramirez, R., Kasting, J. F., Eymet, V., Robinson, T. D., Mahadevan, S., et al. (2013). Habitable zones around main-sequence stars: new estimates. Astrophysical J. 765 (2), 131. doi:10.1088/0004-637x/765/2/131

CrossRef Full Text | Google Scholar

Krissansen-Totton, J., Bergsman, D. S., and Catling, D. C. (2016). On detecting biospheres from chemical thermodynamic disequilibrium in planetary atmospheres. Astrobiology 16 (1), 39–67. doi:10.1089/ast.2015.1327

PubMed Abstract | CrossRef Full Text | Google Scholar

Krissansen-Totton, J., Garland, R., Irwin, P., and Catling, D. C. (2018). Detectability of biosignatures in anoxic atmospheres with the James Webb Space Telescope: A TRAPPIST-1e case study. Astronomical J. 156 (3), 114. doi:10.3847/1538-3881/aad564

CrossRef Full Text | Google Scholar

Krissansen-Totton, J., Thompson, M., Galloway, M. L., and Fortney, J. J. (2022). Understanding planetary context to enable life detection on exoplanets and test the Copernican principle. Nat. Astron. 6 (2), 189–198. doi:10.1038/s41550-021-01579-7

CrossRef Full Text | Google Scholar

Kuhn, J. R., Berdyugina, S. V., Langlois, M., Moretto, G., Thiébaut, E., Harlingten, C., et al. (2014). “Looking beyond 30m-class telescopes: the Colossus project,”, 9145. SPIE, 533–540. Ground-Based Airborne Telesc. V. doi:10.1117/12.2056594

CrossRef Full Text | Google Scholar

Lammer, H., Sproß, L., Grenfell, J. L., Scherf, M., Fossati, L., Lendl, M., et al. (2019). The role of N2 as a geo-biosignature for the detection and characterization of Earth-like habitats. Astrobiology 19 (7), 927–950. doi:10.1089/ast.2018.1914

PubMed Abstract | CrossRef Full Text | Google Scholar

Lenardic, A., and Seales, J. (2021). Habitability: a process versus a state variable framework with observational tests and theoretical implications. Int. J. Astrobiol. 20 (2), 125–132. doi:10.1017/s1473550420000415

CrossRef Full Text | Google Scholar

C. H. Lineweaver, P. C. Davies, and M. Ruse (2013). Complexity and the arrow of time (Cambridge University Press).

Google Scholar

Livengood, T. A., Deming, L. D., A'hearn, M. F., Charbonneau, D., Hewagama, T., Lisse, C. M., et al. (2011). Properties of an Earth-like planet orbiting a Sun-like star: earth observed by the EPOXI mission. Astrobiology 11 (9), 907–930. doi:10.1089/ast.2011.0614

PubMed Abstract | CrossRef Full Text | Google Scholar

Lomb, N. R. (1976). Least-squares frequency analysis of unequally spaced data. Astrophysics space Sci. 39, 447–462. doi:10.1007/bf00648343

CrossRef Full Text | Google Scholar

Lustig-Yaeger, J., Fu, G., May, E. M., Ceballos, K. N. O., Moran, S. E., Peacock, S., et al. (2023). A JWST transmission spectrum of the nearby Earth-sized exoplanet LHS 475 b. Nat. Astron. 7 (11), 1317–1328. doi:10.1038/s41550-023-02064-z

CrossRef Full Text | Google Scholar

Lustig-Yaeger, J., Meadows, V. S., Crisp, D., Line, M. R., and Robinson, T. D. (2023). Earth as a transiting exoplanet: a validation of transmission spectroscopy and atmospheric retrieval methodologies for terrestrial exoplanets. Planet. Sci. J. 4 (9), 170. doi:10.3847/psj/acf3e5

CrossRef Full Text | Google Scholar

Marshall, S. M., Mathis, C., Carrick, E., Keenan, G., Cooper, G. J., Graham, H., et al. (2021). Identifying molecules as biosignatures with assembly theory and mass spectrometry. Nat. Commun. 12 (1), 3033. doi:10.1038/s41467-021-23258-x

PubMed Abstract | CrossRef Full Text | Google Scholar

Meadows, V. S. (2005). Modelling the diversity of extrasolar terrestrial planets. Proc. Int. Astronomical Union 1 (C200), 25–34. doi:10.1017/s1743921306009033

CrossRef Full Text | Google Scholar

Meadows, V. S. (2008). Planetary environmental signatures for habitability and life. Exopl. Detect. Form. Prop. habitability, 259–284. doi:10.1007/978-3-540-74008-7_10

CrossRef Full Text | Google Scholar

Meadows, V. S. (2017). Reflections on O2 as a biosignature in exoplanetary atmospheres. Astrobiology 17 (10), 1022–1052. doi:10.1089/ast.2016.1578

PubMed Abstract | CrossRef Full Text | Google Scholar

Meadows, V. S., Reinhard, C. T., Arney, G. N., Parenteau, M. N., Schwieterman, E. W., Domagal-Goldman, S. D., et al. (2018). Exoplanet biosignatures: understanding oxygen as a biosignature in the context of its environment. Astrobiology 18 (6), 630–662. doi:10.1089/ast.2017.1727

PubMed Abstract | CrossRef Full Text | Google Scholar

Meng, Z., Yang, P., Kattawar, G. W., Bi, L., Liou, K. N., and Laszlo, I. (2010). Single-scattering properties of tri-axial ellipsoidal mineral dust aerosols: a database for application to radiative transfer calculations. J. Aerosol Sci. 41 (5), 501–512. doi:10.1016/j.jaerosci.2010.02.008

CrossRef Full Text | Google Scholar

Moretto, G., Kuhn, J. R., Thiébaut, E., Langlois, M., Berdyugina, S. V., Harlingten, C., et al. (2014). “New strategies for an extremely large telescope dedicated to extremely high contrast: the Colossus project,”, 9145. SPIE, 590–598. Ground-based Airborne Telesc. V. doi:10.1117/12.2055797

CrossRef Full Text | Google Scholar

Morley, C. V., Kreidberg, L., Rustamkulov, Z., Robinson, T., and Fortney, J. J. (2017). Observing the atmospheres of known temperate Earth-sized planets with JWST. Astrophysical J. 850 (2), 121. doi:10.3847/1538-4357/aa927b

CrossRef Full Text | Google Scholar

Nari, N., Dumusque, X., Hara, N. C., Mascareño, A. S., Cretignier, M., Hernández, J. G., et al. (2025). Revisiting the multi-planetary system of the nearby star HD 20794-Confirmation of a low-mass planet in the habitable zone of a nearby G-dwarf. Astronomy and Astrophysics 693, A297. doi:10.1051/0004-6361/202451769

CrossRef Full Text | Google Scholar

Natraj, V., Luo, M., Blavier, J. F., Payne, V. H., Posselt, D. J., Sander, S. P., et al. (2021). Simulated multispectral temperature and atmospheric composition retrievals for the JPL GEO-IR.

Google Scholar

Olejarz, J., Iwasa, Y., Knoll, A. H., and Nowak, M. A. (2021). The Great Oxygenation Event as a consequence of ecological dynamics modulated by planetary change. Nat. Commun. 12 (1), 3985. doi:10.1038/s41467-021-23286-7

PubMed Abstract | CrossRef Full Text | Google Scholar

Packard, N. H., Crutchfield, J. P., Farmer, J. D., and Shaw, R. S. (1980). Geometry from a time series. Phys. Rev. Lett. 45 (9), 712–716. doi:10.1103/physrevlett.45.712

CrossRef Full Text | Google Scholar

Pope, R. M., and Fry, E. S. (1997). Absorption spectrum (380–700 nm) of pure water. II. Integrating cavity measurements. Appl. Opt. 36 (33), 8710–8723. doi:10.1364/ao.36.008710

PubMed Abstract | CrossRef Full Text | Google Scholar

Quanz, S. P., Crossfield, I., Meyer, M. R., Schmalzl, E., and Held, J. (2015). Direct detection of exoplanets in the 3–10 μm range with E-ELT/METIS. Int. J. Astrobiol. 14 (2), 279–289. doi:10.1017/s1473550414000135

CrossRef Full Text | Google Scholar

Ribas, I., Tuomi, M., Reiners, A., Butler, R. P., Morales, J. C., Perger, M., et al. (2018). A candidate super-Earth planet orbiting near the snow line of Barnard’s star. Nature 563 (7731), 365–368. doi:10.1038/s41586-018-0677-y

PubMed Abstract | CrossRef Full Text | Google Scholar

Robinson, T. D., and Reinhard, C. T. (2018). Earth as an exoplanet. Planet. Astrobiol. 379. doi:10.48550/arXiv.1804.04138

CrossRef Full Text | Google Scholar

Robinson, T. D., Meadows, V. S., and Crisp, D. (2010). Detecting oceans on extrasolar planets using the glint effect. Astrophysical J. Lett. 721 (1), L67–L71. doi:10.1088/2041-8205/721/1/l67

CrossRef Full Text | Google Scholar

Robinson, T. D., Meadows, V. S., Crisp, D., Deming, D., A’hearn, M. F., Charbonneau, D., et al. (2011). Earth as an ex trasolar planet: earth model validation using EPOXI Earth observations. Astrobiology 11, 393–408. doi:10.1089/ast.2011.0642

PubMed Abstract | CrossRef Full Text | Google Scholar

Roujean, J. L., Leroy, M., and Deschamps, P. Y. (1992). A bidirectional reflectance model of the Earth's surface for the correction of remote sensing data. J. Geophys. Res. Atmos. 97 (D18), 20455–20468. doi:10.1029/92jd01411

CrossRef Full Text | Google Scholar

Sagan, C., Thompson, W., Carlson, R., Gurnett, D., and Hord, C. (1993). A search for life on Earth from the Galileo spacecraft. Nature 365, 715–721. doi:10.1038/365715a0

PubMed Abstract | CrossRef Full Text | Google Scholar

Scargle, J. D. (1982). Studies in astronomical time series analysis. II-Statistical aspects of spectral analysis of unevenly spaced data. Astrophysical J. Part 1, 835–853. doi:10.1086/160554

CrossRef Full Text | Google Scholar

Schwieterman, E. W., Kiang, N. Y., Parenteau, M. N., Harman, C. E., DasSarma, S., Fisher, T. M., et al. (2018). Exoplanet biosignatures: a review of remotely detectable signs of life. Astrobiology 18 (6), 663–708. doi:10.1089/ast.2017.1729

PubMed Abstract | CrossRef Full Text | Google Scholar

Seager, S., and Bains, W. (2015). The search for signs of life on exoplanets at the interface of chemistry and planetary science. Sci. Adv. 1 (2), e1500047. doi:10.1126/sciadv.1500047

PubMed Abstract | CrossRef Full Text | Google Scholar

Segal, G., Parkinson, D., and Bartlett, S. (2024). Planetary complexity revealed by the joint differential entropy of eigencolors. Astronomical J. 167 (3), 114. doi:10.3847/1538-3881/ad20cf

CrossRef Full Text | Google Scholar

Shields, A. L., Meadows, V. S., Bitz, C. M., Pierrehumbert, R. T., Joshi, M. M., and Robinson, T. D. (2013). The effect of host star spectral energy distribution and ice-albedo feedback on the climate of extrasolar planets. Astrobiology 13 (8), 715–739. doi:10.1089/ast.2012.0961

PubMed Abstract | CrossRef Full Text | Google Scholar

Smith, J. M., and Szathmary, E. (1997). The major transitions in evolution. Oxford: OUP.

Google Scholar

Spurr, R., and Natraj, V. (2011). A linearized two-stream radiative transfer code for fast approximation of multiple-scatter fields. J. Quantitative Spectrosc. Radiat. Transf. 112 (16), 2630–2637. doi:10.1016/j.jqsrt.2011.06.014

CrossRef Full Text | Google Scholar

Thompson, S. L., and Barron, E. J. (1981). Comparison of Cretaceous and present earth albedos: implications for the causes of paleoclimates. J. Geol. 89 (2), 143–167. doi:10.1086/628577

CrossRef Full Text | Google Scholar

Tinetti, G., Meadows, V. S., Crisp, D., Fong, W., Fishbein, E., Turnbull, M., et al. (2006a). Detectability of planetary characteristics in disk-averaged spectra. I: the Earth model. Astrobiology 6 (1), 34–47. doi:10.1089/ast.2006.6.34

PubMed Abstract | CrossRef Full Text | Google Scholar

Tinetti, G., Meadows, V. S., Crisp, D., Kiang, N. Y., Kahn, B. H., Bosc, E., et al. (2006b). Detectability of planetary characteristics in disk-averaged spectra II: synthetic spectra and light-curves of earth. Astrobiology 6 (6), 881–900. doi:10.1089/ast.2006.6.881

PubMed Abstract | CrossRef Full Text | Google Scholar

Vasicek, O. (1976). A test for normality based on sample entropy. J. R. Stat. Soc. Ser. B Stat. Methodol. 38 (1), 54–59. doi:10.1111/j.2517-6161.1976.tb01566.x

CrossRef Full Text | Google Scholar

Vladilo, G., and Hassanali, A. (2018). Hydrogen bonds and life in the universe. Life 8 (1), 1. doi:10.3390/life8010001

PubMed Abstract | CrossRef Full Text | Google Scholar

Walker, S. I., Kim, H., and Davies, P. C. (2016). The informational architecture of the cell. Philosophical Trans. R. Soc. A Math. Phys. Eng. Sci. 374 (2063), 20150057. doi:10.1098/rsta.2015.0057

PubMed Abstract | CrossRef Full Text | Google Scholar

Walker, S. I., Bains, W., Cronin, L., DasSarma, S., Danielache, S., Domagal-Goldman, S., et al. (2018). Exoplanet biosignatures: future directions. Astrobiology 18 (6), 779–824. doi:10.1089/ast.2017.1738

PubMed Abstract | CrossRef Full Text | Google Scholar

Wogan, N. F., and Catling, D. C. (2020). When is chemical disequilibrium in Earth-like planetary atmospheres a biosignature versus an anti-biosignature? Disequilibria from dead to living worlds. Astrophysical J. 892 (2), 127. doi:10.3847/1538-4357/ab7b81

CrossRef Full Text | Google Scholar

Wong, M. L., Bartlett, S., Chen, S., and Tierney, L. (2022). Searching for life, mindful of lyfe’s possibilities. Life 12 (6), 783. doi:10.3390/life12060783

PubMed Abstract | CrossRef Full Text | Google Scholar

Xi, X., Natraj, V., Shia, R. L., Luo, M., Zhang, Q., Newman, S., et al. (2015). Simulated retrievals for the remote sensing of CO 2, CH 4, CO, and H 2 O from geostationary orbit. Atmos. Meas. Tech. 8 (11), 4817–4830. doi:10.5194/amt-8-4817-2015

CrossRef Full Text | Google Scholar

Zechmeister, M., Dreizler, S., Ribas, I., Reiners, A., Caballero, J. A., and Bauer, F. F. (2019). The CARMENES search for exoplanets around M dwarfs-Two temperate Earth-mass planet candidates around Teegarden’s Star. Astronomy and astrophysics 627, A49. doi:10.1051/0004-6361/201935460

CrossRef Full Text | Google Scholar

Zeng, Z. C., Zhang, Q., Natraj, V., Margolis, J. S., Shia, R. L., Newman, S., et al. (2017). Aerosol scattering effects on water vapor retrievals over the Los Angeles Basin. Atmos. Chem. Phys. 17 (4), 2495–2508. doi:10.5194/acp-17-2495-2017

CrossRef Full Text | Google Scholar

Zeng, Z. C., Natraj, V., Xu, F., Pongetti, T. J., Shia, R. L., Kort, E. A., et al. (2018). Constraining aerosol vertical profile in the boundary layer using hyperspectral measurements of oxygen absorption. Geophys. Res. Lett. 45 (19), 10–772. doi:10.1029/2018gl079286

CrossRef Full Text | Google Scholar

Zhang, Q., Natraj, V., Li, K. F., Shia, R. L., Fu, D., Pongetti, T. J., et al. (2015). Accounting for aerosol scattering in the CLARS retrieval of column averaged CO2 mixing ratios. J. Geophys. Res. Atmos. 120 (14), 7205–7218. doi:10.1002/2015jd023499

CrossRef Full Text | Google Scholar

Zhang, Q., Shia, R. L., Sander, S. P., and Yung, Y. L. (2016). XCO2 retrieval error over deserts near critical surface albedo. Earth Space Sci. 3 (2), 36–45. doi:10.1002/2015ea000143

CrossRef Full Text | Google Scholar

Keywords: exoplanet, habitability, remote sensing, DSCOVR, light curves, atmospheric variability, surface composition

Citation: Jian X, Gu L, Fan S, Bartlett SJ, Yang J, Jiang JH, Luo Y and Yung YL (2025) The single beacon: progresses in understanding Earth as an exoplanet using DSCOVR/EPIC observations. Front. Remote Sens. 6:1635015. doi: 10.3389/frsen.2025.1635015

Received: 25 May 2025; Accepted: 19 July 2025;
Published: 08 August 2025.

Edited by:

Alexander Marshak, National Aeronautics and Space Administration, United States

Reviewed by:

Alfonso Delgado Bonal, Universities Space Research Association (USRA), United States
Lorenzo Biasiotti, Trieste Astronomical Observatory, Italy

Copyright © 2025 Jian, Gu, Fan, Bartlett, Yang, Jiang, Luo and Yung. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Siteng Fan, ZmFuc3RAc3VzdGVjaC5lZHUuY24=

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.