3. Results
Analysis of exoplanet habitability within circumstellar habitable zones reveals several critical insights, visualized through a series of figures that illustrate various aspects of the data. By linking these figures together, a narrative is constructed that enhances the understanding of exoplanetary habitability and aligns it with existing literature. Systematic examination of 5,595 confirmed exoplanets from the NASA Exoplanet Archive, applying equation (1) to calculate the average surface temperatures () and categorize their habitable zone status, was performed.
Figure 1 depicts the cumulative count history of all confirmed exoplanet discoveries, highlighting their habitable zone status. The exponential growth in exoplanet discoveries, particularly since the launch of the Kepler Space Telescope, underscores significant advancements in detection technologies and methodologies (Borucki et al., 2010). This trend illustrates our increasing ability to identify potentially habitable exoplanets, reflecting ongoing refinements in search strategies and expanding observational capabilities. The incomplete data for 2024 is indicative of ongoing discoveries, suggesting that the counts of habitable and non-habitable zone exoplanets continued growth, driven by new technologies enabling more advanced data analysis techniques and observation missions.
To understand the spatial distribution of exoplanet discoveries,
Figure 2 illustrates the distance-wise distribution of all confirmed exoplanets from the Solar System, categorized by their habitable zone status. The majority of exoplanets are concentrated within 1,000 light-years, highlighting observational biases where nearer stars are more frequently surveyed due to limitations in current detection technologies. This bias underscores the challenge of detecting distant exoplanets and emphasizes the need for next-generation telescopes capable of probing deeper into the galaxy to identify more distant, potentially habitable exoplanets (Gaudi et al., 2020). The clustering of discoveries within 1,000 light-years also reflects the limitations in the sensitivity and resolution of current instruments, as well as the prioritization of closer stars for detailed observation (Winn & Fabrycky, 2015).
Figure 3 presents the distribution of confirmed single-hosted exoplanets based on their habitable zone status. Among these exoplanets, 77.75% are categorized as "Too Hot," 8.04% as "Too Cold," 4.48% are within the HZ, and 9.73% have indeterminate status (N/A). The significant proportion of "Too Hot" exoplanets suggests an observational bias, as closer-in planets with shorter orbital periods are easier to detect using methods like the Transit method. This finding aligns with previous studies that highlight the detection bias towards short-period exoplanets, often leading to an underrepresentation of planets within the habitable zone (Seager & Mallén-Ornelas, 2003). The skew towards "Too Hot" exoplanets also reflects the challenges in detecting cooler, potentially habitable planets that lie farther from their host stars, where longer orbital periods and lower transit probabilities complicate their detection (Howard et al., 2012).
Delving deeper into the detection methods,
Figure 4 (a) shows the habitable zone status of single-hosted exoplanets discovered via the Transit and Transit Timing Variations methods. Specifically, 89.75% are "Too Hot," 0.90% are "Too Cold," 3.10% are within the HZ, and 6.25% are N/A. The overwhelming majority of "Too Hot" exoplanets discovered through these methods underscores the inherent observational bias towards detecting planets with shorter orbital periods, which are more likely to transit their host star frequently. This bias reflects the limitations of the Transit method in identifying habitable zone exoplanets, as it is more sensitive to planets that orbit closer to their stars (Petigura et al., 2013). The underrepresentation of "Too Cold" and HZ planets emphasizes the need for more sensitive instruments and extended observation periods to detect planets in wider orbits.
Figure 4 (b) presents the habitable zone status of single-hosted exoplanets discovered using the Radial Velocity method. The distribution is more balanced compared to the Transit method, with 49.45% "Too Hot," 35.32% "Too Cold," 11.70% within the HZ, and 3.53% N/A. The Radial Velocity method's sensitivity to planets at various distances from their host stars provides a broader view of exoplanetary systems, although it still shows a detection bias towards larger planets. This method's ability to detect planets in a wider range of orbits highlights its complementary role in identifying habitable zone exoplanets, addressing some of the limitations inherent in the Transit method (Mayor & Queloz, 1995). The broader distribution of HZ exoplanets detected by Radial Velocity indicates its potential in revealing more distant, possibly habitable planets that are missed by transit surveys. Note: Appendix-
Figure 1 extends the analysis presented in
Figure 4 (a) and 4 (b) to the Imaging method, revealing this technique’s detection capability for HZ exoplanets to still be quite limited in comparison.
Moving to the stellar classifications of host stars,
Figure 5 shows the proportions of host stars across all single-hosted exo-systems containing at least one HZ exoplanet. The breakdown is as follows: 26.32% are M-type stars, 29.82% are K-type stars, 35.53% are G-type stars, and 7.89% are F-type stars. This distribution indicates a higher prevalence of HZ exoplanets around G-type and K-type stars, aligning with the fact that these stars are prime targets for habitability studies due to their stable lifetimes and favorable conditions for liquid water (Kasting et al., 1993). The relatively lower proportion of HZ exoplanets around M-type stars, despite their abundance in the galaxy, reflects the challenges in detecting potentially habitable planets around these dimmer, cooler stars (Shields et al., 2016). M-type stars, while abundant, have smaller HZs closer to the star, making planets within these zones more susceptible to stellar activity and tidal locking, which could hinder habitability (Dressing & Charbonneau, 2015).
Figure 6 compares the relative abundance of M, K, G, F, A, and B stellar classes in the Milky Way galaxy to the stellar classes of confirmed exoplanet host stars in single-host systems. For exoplanet host stars, the proportions are 7.8% M-type, 24.7% K-type, 47.4% G-type, 19.4% F-type, 0.6% A-type, and 0.2% B-type. In contrast, the overall abundance in the Milky Way is 76.5% M-type, 12.1% K-type, 7.6% G-type, 3.0% F-type, 0.6% A-type, and 0.13% B-type. This discrepancy highlights a selection bias towards G-type stars, which are similar to our Sun and are often prioritized in exoplanet searches due to their potential for habitability (Reid et al., 2002). The underrepresentation of M-type star hosts further underscores the observational challenges and biases in current exoplanet surveys. The preference for G-type stars also reflects historical biases and the assumption that solar analogs are more likely to host habitable planets (Brown, 2015).
Examining exoplanet discovery methods further,
Figure 7 shows the proportion of single-hosted exoplanets categorized by these various methodologies. The breakdown is as follows: Transit 76.20%, Radial Velocity 17.81%, Pulsation Time Variation 0.04%, Pulsar Timing 0.12%, Orbital Brightness Modulation 0.18%, Microlensing 3.95%, Imaging 1.12%, Astrometry 0.04%, Disk Kinematics 0.02%, and Transit Timing Variations 0.53%.
The dominance of the Transit method reflects its efficiency in detecting exoplanets, particularly those close to their host stars. This method's predominance in exoplanet discovery highlights its strengths in surveying large areas of the sky and detecting numerous exoplanets, although it also emphasizes the need for complementary methods to provide a more complete picture of exoplanetary systems (Winn & Fabrycky, 2015). The significant presence of exoplanets discovered via microlensing and imaging methods showcases their role in identifying planets at greater distances and in different stellar environments (Bennett et al., 2014). Note: Appendix-
Figure 2 focuses the analysis of
Figure 7 specifically to HZ exoplanets, illustrating the co-dominance and near parity in terms of overall counts between the Transit and Radial Velocity methods.
Figure 8 shows the habitable zone width as a function of host star mass, including only those hosts with mass known to ≤10% uncertainty. The Sun is included for reference as a large yellow dot. The figure indicates that the width of the habitable zone increases with the mass of the host star and can be approximated as a power function. This relationship aligns with theoretical models where more massive stars have broader habitable zones due to their higher luminosities, which affect the range of distances at which liquid water could exist on a planet's surface (Kasting et al., 1993). This finding suggests that more massive stars may offer wider potential zones for habitability, though they also present challenges such as shorter lifespans and higher levels of stellar activity (Selsis et al., 2007). The broader HZ around massive stars implies that planets can orbit at greater distances while still maintaining surface conditions conducive to liquid water, but the increased stellar radiation and shorter stellar lifetimes may limit long-term habitability (Kopparapu et al., 2013).
Figure 9 illustrates the variation of habitable zone boundary distances from the host star as a function of host mass, overlaid with exoplanet semi-major axes corresponding to host mass. The inner (red dots) and outer (blue dots) boundaries of the HZ are shown alongside the semi-major axes of exoplanets. The Solar System’s planets, as well as the Sun’s similarly calculated HZ inner and outer HZ boundaries, are depicted for reference. This visualization demonstrates how the HZ boundaries expand outward with increasing host star mass, while the distribution of exoplanet semi-major axes suggests a tendency for planets to reside closer to their stars in lower-mass systems and further away in higher-mass systems. This pattern reflects the influence of stellar mass on planetary formation and orbital dynamics, providing insights into the distribution of potentially habitable exoplanets across different stellar environments (Kopparapu et al., 2013). The data suggest that in lower-mass systems, planets are more likely to be found closer to the star, within the narrower HZ, while in higher-mass systems, planets can be situated further out within the broader HZ, potentially offering a more stable environment for life (Lammer et al., 2009).
Figure 10 illustrates the relationship between the effective temperatures of single host stars and the surface temperatures of their corresponding exoplanets, with exoplanets grouped by type: Gas Giants, Neptunian Planets, Super-Earths, and Terrestrial Planets. The Earth's position is labeled for reference, along with those for the Solar System’s other planets, with the habitable zone range indicated within the vertical green dashed lines. The size of each point represents the relative size of the exoplanet. An exoplanetary system of particular note is TRAPPIST-1, consisting of seven confirmed planets and the host itself, an ultra-cool (effective surface temperature of 2,566 K) red dwarf much smaller in radius than the Sun. Notwithstanding, the TRAPPIST-1 system contains two HZ planets and three of terrestrial size, one of which is among the HZ pair as calculated using equation (1). Given its intriguing array of exoplanets, this system has been one of the more closely studied since its discovery in 2016.
This figure provides a comprehensive visualization of the relationship between stellar temperatures and exoplanet surface temperatures. It shows that few exoplanets fall within the HZ range, indicating the rarity of conditions suitable for surface liquid water. Higher stellar temperatures generally correspond to higher planetary surface temperatures, as evident from the upward trend of data points. Gas Giants and Neptunian Planets predominantly lie outside the HZ, while Super-Earths and Terrestrial Planets show a greater propensity to occupy or approach the HZ. This visualization underscores the need for advanced observational technologies and methodologies to discover and characterize exoplanets within the HZ. Future missions equipped with direct imaging capabilities and improved sensitivity are essential to identifying and studying Earth-like planets in habitable zones of a broader range of stellar types.
This figure also brings attention to assumption limitations introduced earlier in this paper. For the planets in our solar system, surface temperatures calculated with Equation (1) generally align with known mean temperatures (
https://science.nasa.gov/resource/solar-system-temperatures/) with a 6-37% difference. However, Venus is an exception. Although its surface is too hot for life as we know it, Equation (1) flags the planet as within the habitable zone. This discrepancy arises from the assumption of a standardized bulk temperature factor (k=1.13k) based on Earth's values when accounting for the atmospheric greenhouse effect. In reality, Venus has a very thick atmosphere composed primarily of CO
2, trapping heat and resulting in a much higher bulk temperature factor (k=3.17). This limitation is discussed earlier in section 2.6. Bracketing the inner Solar System-based atmospheric greenhouse assumption, this on the cooler end, is Mars. While the atmosphere of Mars is also predominantly composed of CO
2, it is far less dense and accordingly much less capable of trapping solar radiation. The particular exception of Venus indicates that variations in the atmospheric greenhouse effect will need to be further considered to better determine exoplanet surface temperatures.
Overall, the analysis presented above highlights several observational biases inherent in current exoplanet detection methods. The Transit method, which dominates exoplanet discoveries, is more likely to detect planets with shorter orbital periods, leading to an overrepresentation of "Too Hot" exoplanets. This bias is evident in the high percentage of "Too Hot" planets discovered via Transit methods (
Figure 3 and
Figure 4a). Conversely, the Radial Velocity method, which can detect exoplanets at various distances from their host stars, presents a more balanced distribution of habitable zone statuses (
Figure 4b).
The stellar classification of host stars (
Figure 5 and
Figure 6) shows a preference for G-type stars in exoplanet searches, despite M-type stars being by far the most common in the Milky Way. This selection bias may be due to the more stable and longer-lived nature of G-type stars, which are conducive to sustaining life-supporting environments over extended periods.
Figure 8 and
Figure 9 provide insights into the relationship between host star mass and habitable zone characteristics. The widening of the habitable zone with increasing host star mass suggests that more massive stars offer a larger range of distances where conditions might support surface liquid water. However, the semi-major axis distribution of exoplanets indicates a concentration of planets closer to lower-mass stars, potentially due to the higher likelihood of detecting such planets through current observational techniques – again echoing the dominance of the Transit method.
Figure 10 adds additional context by illustrating the relationship between the effective temperatures of host stars and the surface temperatures of their corresponding exoplanets. This figure highlights the challenges in finding planets within the HZ and underscores the importance of advanced observational technologies to overcome current limitations.
Overall, these results emphasize the need for continued development and deployment of diverse detection methods to achieve a more comprehensive understanding of exoplanetary systems. Future missions should aim to mitigate observational biases by targeting a broader range of stellar types and distances, thereby enhancing our ability to identify potentially habitable exoplanets and build a more complete understanding of planetary systems in general.