Showing posts with label Exoplanets. Show all posts
Showing posts with label Exoplanets. Show all posts

Monday, 16 June 2025

Direct imaging of 14 Herculis c by the James Webb Space Telescope.

The planet 14 Herculis c was discovered by the ELODIE Planet Search Survey in 2005, using the radial velocity method, which it was detected by the gravitational effect it has on its host star, 14 Herculis A, as these cause the star to wobble slightly on its axis. 14 Herculis is the outermost of two known planets in the 14 Herculis system, the other being 14 Herculis b. The star 14 Herculis A is a K-type orange dwarf star, slightly smaller than the Sun, 58.4 light years from our Solar System in the constellation of Hercules. The planet 14 Herculis c orbits this star at about 20 AU (i.e. about 20 times the average distance between the Earth and the Sun), with a orbital period of 152.8 years. It is estimated to be about the same size as Jupiter, but much more massive. This is due to the distance at which it orbits its star, which makes it much cooler than Jupiter, making the gasses from which it is made more dense.

Although several subsequent studies have confirmed the existence of 14 Herculis c, the planet had never been directly imaged. This changed in 2025, when the James Webb Space Telescope's Near Infrared Camera was trained upon the 14 Herculis system. The Webb study also enabled the direct measurement of the planet's temperature for the first time, showing that it's average surface temperature may be around -3°C.

14 Herculis c. The view is mostly black, with very faint red splotches in the central region of the image. At the center of the image, there is a black circle, and in the center of that, there is a star symbol representing a real star. This black circle blocks the light from the host star. To the lower right of the circle is a fuzzy bright orange circle, which is the exoplanet. NASA/ESA/CSA/STScl/William Balmer/Daniella Bardalez Gagliuffi.

14 Herculis c was imaged at 4.44 microns, an infrared wavelength equivalent to a temperature of -3°C. Although this was the predicted temperature for the exoplanet, it was found to be much dimmer than expected. This is theorized to indicate a much more active atmosphere than predicted, with warmer gas molecules from the planet's interior being brought to the surface rapidly due to internal churning.

See also...

Saturday, 14 December 2024

Direct imaging of a possible planet orbiting the fast moving star HIP 36277.

In the past three decades almost six thousand exoplanets (planets orbititing stars other than our Sun) have been detected, using a variety of methods. Direct imaging has proven to be a useful technique for detecting planets with masses greater-than-or-equal-to that of Jupiter at distances of more than 10 AU from their host stars (i.e. more than ten times as far from their host stars as the Earth is from the Sun). Notable planets discovered in this way include 51 Eridani b, which has a mass 2.6 times that of Jupiter and orbits a star 96 light years from Earth in the constellation of Eridanus at a distance of 11.1 AU, HIP 65426 b (formally named Najsakopajk), which has a mass 7.1 times that of Jupiter, and orbits a star 385 light years from Earth in the constellation of Centaurus at a distance of 87 AU, and PDS 70 b and PDS 70 c, which have masses of 3.2 and 7.5 times that of Jupiter, and orbit a star 370 light years from Earth in the constellation of Centaurus at distances of 20.8 and 34.3 AU, respectively.

However, large planets at large separations from their host stars are relatively rare, which means that a large number of stellar systems have to be surveyed in this way in order to detect a few planets. Most giant planets known orbit their host stars at distances of 1-3 AU, which would make direct imaging them with current technology impossible if they are more than about 50 parsecs (163.1 light years) away from us. 

The radial velocity method uses the movement of stars to detect to infer the presence of companions. This has proven very effective as a  way to detect very large planets close to stars, such as Beta Pictoris c, which has a mass about nine times that of Jupiter, and orbits a star 63 light years from Earth in the constellation of Pictor at a distance of about 2.7 AU, or HD206893 c, which has a mass about 12.7 times that of Jupiter, and orbits a star 125 light years from  Earth in the constellation of Capricornus,  at a distance of 3.53 AU.

The proper motion anomaly method can identify potential companions to stars by measuring their parallax (the amount they move in a year because we are observing them from different points on the Earth's orbit) over several years; if the star moves more than predicted (i.e. anomalously), then this is likely to be because of an unseen companion moving the star. This method has been used to identify several potential planets which have subsequently been directly imaged. These include HIP 99770 b, which has a mass of about 16 times that of Jupiter, and which orbits a star 133 light years from Earth in the constellation of Cygnus at a distance of 17 AU, AF Leporis b, which is 2-5 times the mass of Jupiter and orbits a star 87.5 light years from Earth in the constellationof Lepus, as well as the brown dwarf HD21152 B, which has a mass 22-36 times that of Jupiter, and which orbits a star 150 light years from Earth in the constellation of Taurus at a distance of about 18 AU.

In a paper published in the Monthly Notices of the Royal Astronomical Society on 9 December 2024, a team of astronomers led by Dino Mesa of the Osservatorio Astronomico di Padova present the results of a study which targeted three stars in the Hipparcos-Gaia PMa catalogue identified as having proper motion anomalies with the SHARK-NIR coronagraphic camera and LMIRCam camera and coronagraph of the Large Binocular Telescope in Arizona. 

Because they were interested in planetary-sized companion bodies, rather than secondary stars, Mesa et al. looked for stars which showed small proper motion anomalies, and because they wished to be able to image bodies within 10 AU of their host star, they restricted themselves to objects within 50 parsecs (163 light years) of the Earth.

The first star selected, HIP 11696 A (also known as HD 15407 A) is an F-type (yellow-white dwarf) star with a mass about 1.40 times that of our Sun, 49.3 parsecs (160.8 light years) from Earth in the constellation of Persius. HIP 11696 is a young star, which has been estimated to be about 80 million years old, although it is also thought likely to be a member of the AAB Doradus Moving Group, which would make it between 125 and 149 million years old. Mesa et al. use an intermediate age of 137 million years for their calculations in their study.

HIP 11696 A has a companion star, HIP 1696 B, which is a K-type (orange dwarf) star with about 80% of the mass of our Sun, separated by about 1000 AU - far enough to be excluded from the field of view of the SHARK-NIR instrument. HIP 11696 A appears to be producing an unusual amount of light in the mid-infrared range, which may be indicative of a recent collision between rocky planets of planetary embryos in the inner part of the system. A debris disk has been detected at a distance of 0.6-1.0 AU from the star, which makes it unlikely that there are any massive planets orbiting close to the star. Nevertheless, an anomaly in the motion of HIP 11696 A which could not be explained by the presence of HIP 11969 B was detected. It has been suggested that this might be caused by a planet with a mass about 6.39 times that of Jupiter orbiting at 3 AU from HIP 11696 A, or a planet with a mass about 16.6 times that of Jupiter orbiting at about 30 AU from the star.

HIP 11696 A was observed with the SHARK-NIR and LMIRCam instrument on the night of 28 October 2023. Mesa et al. detected a bright object to the southeast of the star at a distance of 1.5" (1.5 arc seconds; the sky can be imagined as a sphere surrounding the Earth, divided into 360 degrees (°), with each degree divided into 60 arc minutes (') and each arc minute divided into 60 arc seconds (")). However, this object was also imaged previously by the Keck II telescope in November 2009, and the Gemini North Telescope in August 2013, with no movement relative to HIP 11696 A between these images, leading Mesa et al. to conclude that this is in fact a background object rather than a planetary companion to the star. Based upon this inability to image a planet close to the star, Mesa et al. calculate that if a planet is responsible for the observed wobble in HIP 11696 A's orbit, then this is likely to be between 2.5 and 28 AU from the star, and have a mass 4-16 times that of Jupiter.

(Top) Final image obtained for HIP 11696 using SHARK-NIR data. This image was obtained by applying a PCA method subtracting 5 principal components. (Bottom) Final image obtained for HIP 11696 using LMIRCam data. In this case, a PCA method subtracting 10 principal components was applied. In both cases, a bright candidate companion is visible South-East from the star. Because the image is looking up, the positions of east and west are reversed. Mesa et al. (2024).

The second star identified, HIP 47110 A (also known as HD 82939 A) is a G-type (yellow dwarf) star with about 98% of the mass of our Sun, 38.7 parsecs (126.2 light years) from Earth in the constellation of Leo Minor. HIP 47110 A has been identified as a possible member of the Pleiades Moving group, with an age of approximately 112 million years.

HIP 47110 A has a companion star, HIP 47110 B, which is a M-type (red dwarf) star with a separation of larger than 162" (interpreted to be more than 6280 AU), enabling it to be excluded from the field of vision. Again, there is an anomaly in the motion of HIP 47110 A which cannot be explained by the presence of HIP 447110 B, and which has been hypothesized to be caused by a planet. It has been suggested that this might be caused by a planet with a mass about 2.5 times that of Jupiter orbiting at between 5 and 10 AU from HIP 47110 A, or a planet with a mass about 11.35 times that of Jupiter orbiting at about 30 AU from the star.

HIP 47110 A was observed on the night of 20 February 2024, but no potential companion was observed. Based upon this, Mesa et al. exclude the possibility of a planet close to star, calculating that the observed orbitary wobble must be caused by a planet between 3 and 30 AU from the star with a mass of between 2 and 10 times that of Jupiter.

The third star in the study, HIP 36277, is a K-type (orange dwarf) star with a mass 0.67 times that of our Sun, located 46.3 parsecs (151 light years) from the Earth in the constellation of Dorado. HIP 36277 was identified as a young runaway star (star which has been ejected from the star cluster which birthed it, and which is therefore travelling at a high speed in a distance at odds with galactic rotation) with an age of about 41.2 million years. However, spectrographic analysis of the star has suggested a much older age, most probably more than a billion years old and possibly more than 10 billion years. 

An anomaly on the motion of HIP 36277, which has been interpreted as potentially due to a planet with a mass 2.3 times that of Jupiter at a distance of 5 AU from the star, 2.64 times that of Jupiter at 10 AU from the star, or 15.18 times the mass of Jupiter at 30 AU from the star.

HIP 36277 was observed on the night of 21 February 2024, with a bright object observed to the southeast of the star in both SHARK-NIR and LMIRCam images. This body could also be identified in images from the Gaia space telescope, with similar parallax and proper motion values, which demonstrates physical association with the star. The object is separated from the star by 1.9". The precise size of this object is difficult to calculate, given the uncertainty of the age of the star, but Mesa et al. calculate that if the star is 41 million years old, then it would have a mass of between 16.2 and 73.9 times that of Jupiter (with a median value of 37.8 Jupiter masses), making it most likely a brown dwarf companion to the star (brown dwarfs are objects intermediate to stars and planets in size; they are not large enough to fuse ordinary hydrogen in their cores, but are large enough to fuse the heavier isotope deuterium). However, if the star is about five billion years old, then the body is likely to have a mass about 0.1 times that of our Sun, making it a small M-type (red dwarf) star. 

A second object was also visible in the SHARK-NIR data, but not observed by LMIRCam. This object is to the south of the star, and separated by about 0.0625", which Mesa et al. calculate to be equivalent to about 28.9 AU. Again, the mass of such an object would be dependent on its age, with a 41-million-year-old object having a mass about 7.6 times that of Jupiter, making it a large planet, while at 5 billion years old it would have a mass between 65.9 and 72.1 times that of Jupiter, again indicative of a brown dwarf.

(Top) Final image obtained for HIP 36277 using SHARK-NIR data. This image was obtained by applying a PCA method subtracting 5 principal components. (Bottom) Final image obtained for HIP 36277 using LMIRCam data. In this case, a PCA method subtracting 10 principal components was applied. In both cases, a bright candidate companion is visible southeast from the star. Furthermore, in the SHARK-NIR image, a possible fainter object is visible just south of the star. Mesa et al. (2024).

See also...

Wednesday, 17 January 2024

A new image of the Beta Pictoris system.

The European Space Agency has released a new image of the Beta Pictoris, as seen by the James Webb Space Telescope, in a press release issued on 10 January 2024. The image shows the two main debris disks of the system, which are offset from one another by about 5°, as well as a new and previously unseen feature, which has been named the 'Cat's Tail'.

The Beta Pictoris System. European Space Agency.

Debris disks are rings of dust, rock and icy material left surrounding stars after planet formation has occurred (unlike protoplanetary disks, which are present around very young stars only, and which are thought to be largely consumed by planetary formation). Our Solar System has two such debris disks, the Asteroid Belt and the Kuiper Belt, and in recent years improved telescope technology has allowed astronomers to detect debris disks around about 80 of other stars. Most of these disks have been detected only at infrared wavelengths; imaging disks at visual wavelengths enables astronomers to study the morphology and composition of such disks in far greater detail, enabling them to predict the presence of planets that cannot be directly detected and better understand planetary formation around other stars.

Beta Pictoris is a young star 62.4 light years from Earth in the constellation Pictor. Although it is thought to be only 20-25 million years old, it is 1.75 times as massive as the Sun, and therefore 8.7 times as luminous. 

The main debris disk of Beta Pictoris was discovered in 1983 by the Infrared Astronomical Satellite, one of the first extra-Solar debris disks ever discovered. This disk is seen edge-on from the Earth, and extends about 1835 from the star (i.e 1835 times as far from the star as the Earth is from the Sun), with the northeast arm rotating towards us and the southwest arm rotating away.

In 2006 the Hubble Space Telescope discovered as much fainter secondary debris disk around Beta Pictoris, offset by about 5° from the main disk. This secondary disk extends about 130 AU from the star, and is thought to be caused by a massive planet deflecting material from the main disk onto an orbit aligned with its own. 

The precise nature of the Cat's Tail is unclear, but in 2014 the Atacama Large Millimeter/submillimeter Array detected a clump of carbon monoxide at a similar location. The high levels of radiation emitted by the star should break down carbon monoxide quite quickly, suggesting that the gas was a result of a recent event. The Cat's Tail is thought to be comprised of dusty material, with an equivalent mass similar to that of a large body in the Main Asteroid Belt of our Solar System, spread our over a distance of about 16 billion km (107 AU). This is probably only offset from the main debris disk by about 5°, but appears curved because of the rotation of the system. It is likely that the dust of the Cat's Tail and the carbon monoxide clump detected by the Atacama Large Millimeter/submillimeter Array are the result of a single, catastrophic event.

See also...

Friday, 4 December 2020

Simulating the orbital dynamics of the Kepler-1649 System and the climate of Kepler-1649c.

Over the past several decades, thousands of exoplanets have been discovered enabling statistical studies of exoplanet occurrence rate. These studies include the investigation of planets that lie within the Venus Zone and Habitable Zone. A primary purpose of such studies is to place our solar system within the context of planetary system architectures. One aspect of this context is determining the frequency of systems that have terrestrial planets in both the Venus Zone and Habitable Zone. Such systems may provide potential analogs for Venus and Earth, that in turn may provide key insights into the divergence of the Venus/Earth atmospheric evolutions. Examples of such systems include TRAPPIST-1, LHS 1140, and Kepler-186. Amongst the factors that influence the climate evolution of terrestrial exoplanets are the eccentricity variations, which may have contributed to the water loss of Venus. Thus, a pathway to understanding the major factors behind the Venus/Earth divergence is exploring the dynamical impacts of eccentricity variations that may have occurred (and continue to occur) in exoplanetary systems.

A particularly strong contributor to the currently known inventory of terrestrial exoplanets was the Kepler mission, from which a catalog of Habitable Zone planets was constructed. The Kepler-1649 system was discovered in 2017, which discovery identified the observed planet as a potential Venus analog due to the planetary size and similarity of insolation flux to that of Venus. The prospect of the planet as a Venus analog was strengthened by a set of detailed climate simulations that showed the rapid rise in surface temperature from a variety of initial conditions. The system became a strong comparative planetology target when a further terrestrial planet was discovered in the Habitable Zone earlier this year. Thus, Kepler-1649 has joined a select group of systems that are known to harbour terrestrial planets in both the Venus Zone and Habitable Zone. This class of system also provide opportunities to study the role of complex dynamical evolution on the orbital variation of planets, including those within the Habitable Zone.

In a paper published on the arXiv Database at Cornell University on 18 November 2020, Stephen Kane and Zhexing Li of the Department of Earth and Planetary Sciences at the University of California, Riverside, Eric Wolf of the Laboratory for Atmospheric and Space Physics at the University of Colorado, Boulder, and Colby Ostberg, and Michelle Hill, also of the Department of Earth and Planetary Sciences at the University of California, Riverside, present a dynamical and climate study of the Kepler-1649 system, describing the potential eccentricity variations that result from the interactions between the planets, and the possible impact on the climate of the Habitable Zone planet.

The Kepler-1649 system is particularly fascinating from a comparative planetology aspect, since it contains both potential analogs to Venus and Earth within the same system. The stellar properties of Kepler-1649A include its mass (0.1977 times that of the Sun), radius (0.2317 times that of the Sun), luminosity (0.00516 times that of the Sun), and effective temperature (3240 K, compared to 5778 K for our Sun). The radii of the known planets in the system are 1.017 times that of the Earth for Kepler-1649b and 1.06 times that of the Earth for Kepler-1649c. To calculate the semi-major axes of the planets (average distance at which the planets orbit), Kane et al. used information provided in the paper which described the discovery of Kepler-1649c, resulting in semi-major axes of 0.0842 AU (i.e. 8.42% of the distance at which the Earth orbit's the Sun) for Kepler-1649b, and 0.0827 AU for Kepler-1649c. The Conservative Habitable Zone boundaries are defined by the runaway greenhouse limit at the inner edge and the maximum greenhouse limit at the outer edge. The Optimistic Habitable Zone boundaries are determined empirically based on the assumption that Venus may have had surface liquid water as recently as 1 billion years ago, and that Mars may have surface liquid water 3.8 billion years ago. Based on these system properties, Kane et al. calculated the extent of the Conservative Habitable Zone and Optimistic Habitable Zone boundaries for Kepler-1649 as 0.075–0.147 AU and 0.059–0.155 AU respectively. The outer edge of the Venus Zone corresponds to the inner edge of the Conservative Habitable Zone (the runaway greenhouse limit) at 0.075 AU.

 
A top-down view of the Kepler-1649 system, where the star is location at the center of the crosshairs and the orbits of the known planets are indicated by solid circles. The extent of the Conservative Habitable Zone is shown as the light green region and the Optimistic Habitable Zone extension to the Conservative Habitable Zone is shown as the dark green region. Kane et al. (2020).

Previous descriptions of the Kepler-1649 system have not constrained the eccentricities in any significant way, though smaller planets tend to have small eccentricities. It has further been suggested that there could be an additional planet in an orbit that lies between the two known planets. This suggestion was partially motivated within the context of mean motion resonances that may exist within the system. The known b and c planets are close to the 9:4 mean motion resonance, which is a relatively weak resonance in terms of planetary perturbations. However, an additional planet between the known planets at a location of 0.063 AU would create a 3:2 resonant chain with high stability potential. Furthermore, an analysis of Kepler systems has found that planetary architectures may preferentially contain similar size planets, with the vast majority of Kepler planets having a spacing greater than 10 mutual Hill radii (the Hill radius of a body is the region in which it dominates the attraction of satellites). Assuming an Earth mass for all three planets, a planet located at the 3:2 resonant chain location between planets b and c would be 12.5 mutual Hill radii away from either planet, consistent with the observed spacings in other Kepler systems.

To investigate the prospect of an additional planet as a function of eccentricity, Kane et al. conducted a dynamical simulation using the REBOUND orbital dynamics package. Such an investigation requires an estimate of the planetary masses. There are numerous planetary mass-radius relationships that have been developed from which such an estimate may be derived For their dynamical simulations, Kane et al. adopted the mass estimates provided by the Forecaster tool, which produces masses estimates of 1.06 and 1.21 times that of the Earth for planets b and c respectively.

The dynamical simulation was carried out using the MEGNO (Mean Exponential Growth of Nearby Orbits) package within REBOUND with the symplectic integrator WHFast. MEGNO is a chaos indicator that is useful in dynamical simulations to distinguish quasi-periodic or chaotic orbital time evolution of planetary systems. The final numerical value returned by a MEGNO simulation indicates whether the system would end up in a chaotic state or not, where a chaotic state is less likely to maintain long-term stability. For their simulation, Kane et al. explored the possibility of an additional planet in between the orbits of planet b and c by testing how varying eccentricity and semi-major axis values for an injected planet of one Earth mass affects the stochasticity of the system. Kane et al. conducted four sets of simulations that varied eccentricity, each providing a grid of orbital parameters for the injected planet’s eccentricity (0.0–0.3) and semi-major axis (0.0482–0.0827 AU). Since the eccentricities of the two known planets are unconstrained, Kane et al. tested planet b and c eccentricities of 0.0, 0.1, 0.2, and 0.3. Each simulation duration was 10 million orbits for the outer most planet with a time step of 0.001 years (0.365 days). The integration was set to stop and return a large MEGNO value if the planet was ejected beyond 100 AU. 

There is a large parameter space between approximately 0.052 AU and 0.075 AU with low eccentricity values where the injected Earth-mass planet could be located with long-term stability. The largest range of eccentricities for the injected planet that allow a stable configuration occur at 0.063 AU which, as stated above, is the location of the 3:2 mean motion resonance chain. The additional spikes visible in the stability distribution, most obvious in the top-left panel, also correspond to mean motion resonance locations, such as the 5:4 resonance location at 0.071 AU. However, as the eccentricities for both planet b and c increase to higher values, the stable locations for the injected planet indicated gradually diminishes and the locations shift to higher eccentricity values. When the eccentricities for both planet b and c exceed 0.3, the stable locations for the injected planet vanish. Even for an eccentricity of 0.1, the system is most stable when the pericenters of the planets are approximately aligned. It is worth noting that the described simulations assumed system coplanarity, and the introduction of mutual inclinations may present additional regions of stability, particularly for scenarios that mitigate potential overlap of orbits.

A non-zero eccentricity for the known planets can place significant restrictions on the presence of an additional terrestrial planet between them. Kane et al. next investigate the time-dependent dynamical effects of non-zero eccentricities. Such investigations have been carried out in numerous other works in a more general context. Kane et al.'s N-body integrations were performed with the Mercury Integrator Package. The simulations used a hybrid symplectic/Bulirsch-Stoer integrator with a Jacobi coordinate system to provide increased accuracy for multi-planet systems. The time resolution of the simulations were set to 0.2 days in order to comply with the recommendation that resolutions be 1/20 of the shortest orbital period (8.69 days for planet b). 

Kane et al. performed a series of dynamical simulations using the above criteria that explore a range of eccentricities. Specifically, they started each simulation with a circular orbit for planet b and eccentricities for planet c that ranged from 0.0 to 0.5 in steps of 0.05. Each simulation was allowed to run for 1 million years, or 18 million orbits of planet c, with orbital parameters for both planets out-put every 10 simulation years. The suite of simulations that explored the full range of starting eccentricities for planet c were conducted using Earth masses for both planets b and c, and also using the Forecaster masses.

The full set of dynamical simulations demonstrated that the system is able to remain stable for 1 million years up to a starting eccentricity for planet c of 0.325, beyond which the system rapidly loses dynamical integrity. The 0.325 eccentricity upper limit was true for both the Earth mass and Forecaster mass scenarios. In particular, Kane et al. found that the eccentricity of planet c induces a near-equivalent eccentricity in planet b, whereby each  planet oscillates in eccentricity in order to conserve the system angular momentum.

To investigate this further, Kane et al. conducted a fourier analysis of each of the dynamical simulation outputs described above. The periods of the eccentricity oscillations, for both Forecaster and Earth masses, as a function of starting planet c eccentricity. There are several aspects to note about this. Firstly, the period of the oscillations decreases with increasing eccentricity. Secondly, the period of the oscillations is systematically smaller for the larger Forecaster masses. Thirdly, the dynamical instability that results beyond eccentricities of 0.3 is clearly visible and is likely related to the diminishing oscillation period with increasing eccentricity. These simulation results agree well with similar calculations that may be performed using the analytical expressions derived by Sam Hadden and Yoram Lithwick. It is worth noting that the period of the eccentricity oscillations and the eccentricity at which instability occurs is likely influenced by the proximity of the planets to the 9:4 mean motion resonance. Additional factors include the mutual inclinations of the planets combined with non-zero stellar obliquity. Regardless, the high frequency of eccentricity oscillations may result in significant climate consequences for the planets. 

Kane et al. discuss the implications of the variable eccentricity for the possible climate scenarios of planet c. Given the system parameters produced the insolation flux received by planet c at its semi-major axis is about 75% of the mean Earth insolation flux. The eccentricity oscillations described will result in a correspondingly variable insolation flux for the planet, whose amplitude and period depend on the eccentricity. For example, consider the case of a starting eccentricity of 0.2 for planet c. Since, the eccentricity of planet c never reaches zero, the maximum insolation flux is always larger than that for a circular orbit. The maximum insolation in this case oscillates approximately around an Earth equivalent flux with a variation of ±20%.

The Resolving Orbital and Climate Keys of Earth and Extraterrestrial Environments with Dynamics (ROCKE-3D), is a three-dimensional climate system model developed at the NASA Goddard Institute for Space Studies specifically for studying planetary and exoplanetary atmospheres. Kane et al. used ROCKE-3D to constrain the climatic effects of oscillations in the orbital-rotational properties of Kepler-1649c. ROCKE-3D has previously been used for studying the climate and habitability of terrestrial extrasolar planets, including studies of Proxima Centauri b, moist and runaway greenhouse states, the effects of planet orbital-rotational properties, and the effects of ocean circulation. While the atmospheric composition and surface properties of Kepler-1649c remain unconstrained, here we have conducted a set of simulations using ROCKE-3D to illuminate the possible role of eccentricity evolution on climate. For a starting eccentricity of 0.3, Kepler-1649c could experience eccentricity oscillations on 1200 year timescales. To constrain the bounds of the potential climate states experienced by Kepler-1649c during its dynamical evolution, Kane et al. simulated scenarios with 0 eccentricity and synchronous rotation, and with an eccentricity of 0.3 and a 3:2 resonant rotation state. Note, for these cases Kane et al.

were not simulating the continuous evolution of Kepler-1649c over many thousands of years, but rather were simulating the equilibrium climate states for potential end-member orbital-rotational configurations. However, if the spin dynamics become chaotic, these end-member rotation states may not be fully realised. Thus, Kane et al. also conducted simulations that assume a constant, non-resonant asynchronous spin rate in order to approximate chaotic spin-dynamics. For these cases, they assumed a 16 Earth-day rotation period, which is an intermediate value between the synchronous rotation period of 19.535 days, and the 3:2 resonant rotation period of 13.0235 days. Each of Kane et al.'s simulations achieve radiative equilibrium after 100–200 years, well within the timescales of possible eccentricity variations described. This is also comparable to the timescales provided by the continuously variable eccentricity approach.

Kepler-1649c lies within the Conservative Habitable Zone, close to the inner edge. In all climate simulations conducted, Kane et al. assume a globally ocean covered planet with no emergent continents. While this represents only one of many potential surface configurations for Kepler-1649c, the study of ocean covered worlds provides a first order approximation for the study of climate evolution on potentially habitable worlds.

Kane et al. tested two different assumptions for the ocean layer. 'Slab' ocean simulations assume a 100 m deep ocean layer with no horizontal transports, but with the thermodynamic exchange of heat and moisture between the ocean and atmosphere. 'Dynamic' ocean simulations utilise a 900 m deep fully coupled ocean circulation model capable of predicting horizontal and vertical transports, in addition to thermodynamic exchanges with the atmosphere. Both ocean configurations allow for sea-ice formation. Salinity is held equal to that of the modern Earth’s ocean. Atmospheric compositions also remain unconstrained from observations and therefor constitute a large uncertainty on the climate state. Kane et al. chose two nominal atmospheric compositions which are 'Earth-like' (1 bar total with 400 ppm carbon dioxide and an nitrogen background), and a 1 bar carbon dioxide atmosphere. For all climate simulations, water vapor and clouds are variable constituents in the atmosphere, controlled by temperature, evaporation, convection, advection, and condensation. For planets residing in the middle of the habitable zone, as does Kepler-1649c, an Earth-like atmospheric composition will result in a cold climate state, dominated by ice and where ocean transport processes are of particular importance, while a 1 bar CO2 atmosphere will result in a warm climate state, dominated by the hydrological cycle, clouds, and atmospheric transport. Thus the selection of these two atmospheric compositions allow Kane et al. to capture some diversity of habitable zone climates. A wide ranging parameter study of the infinite variety of possible surface and atmospheric compositions for Kepler-1649c is beyond the scope of this work. In total, we have conducted 16 different climate simulations of Kepler-1649c

Kane et al. produced surface temperature maps from our climate simulations with ROCKE-3D. The results used for these are from an average of the last 100 orbits of Kepler-1649c (i.e. 1950 Earth days, 5.3 Earth years). The bold white line indicates where the surface temperature is 273 K, and thus serves as a close proxy for the sea ice margin. The columns show climate results from different rotation states, including synchronous, resonant, and chaotic. Synchronous rotation assumes zero eccentricity and a rotation period equaling the orbital period. Resonant rotation assumes a rotation period in a 3:2 resonance with the orbital period and an eccentricity of 0.3. Chaotic rotation, as mentioned previously, is approximated by a 16 day asynchronous rotation period. The global mean surface temperature and surface temperature maps are largely unchanged for chaotic rotators with eccentricities of zero or 0.3.

 
Surface temperatures maps from Kane et al.'s climate simulations of Kepler-1649c. The left column assumes synchronous rotation with zero eccentricity, the center column assumed 3:2 resonant rotation with an eccentricity of 0.3, and the right column assumes chaotic rotation approximated by a 16 day asynchronous orbit. For the chaotic states, eccentricity does not have a significant effect on the climate and thus zero and 0.3 eccentricity results are averaged together. Rows (labeled) illustrate different surface an atmospheric configurations. All simulations performed with ROCKE-3D. Kane et al. (2020).

Global mean surface temperatures range between 244.7 K and 304.8 K, and surface liquid water is maintained over at least some of the planet surface in all cases simulated. Thus, all climates simulated are at least partially habitable, with many states being globally habitable. If the orbital-rotational state of Kepler-1649c can fully oscillate between synchronous and resonant modes, significant temporal changes to the global mean surface temperatures and the total habitable area of the planet, which here Kane et al. use the open ocean fraction as a proxy for the habitable surface area, would occur. The shift from synchronous rotation to resonant rotation affects climate both by changing the longitudinal patterns of stellar flux incident on the planet, and also by changing the Coriolis parameter which feeds back on atmospheric and ocean circulations. The breaking of synchronicity and the increase in rotation rate tends to result in increased zonal transports resulting in longitudinally banded patterns, as opposed to the classic eyeball climate states found for slow and synchronous rotators. The inclusion of non-zero eccentricity with a 3:2 orbital resonance produces a characteristic double hot-spot evident in temporal mean climate statistics, due to the fact that two opposite faces of the planet are periodically irradiated at perihelion. The double hot spot pattern is clearly evident in slab ocean cases. The orbital period of Kepler-1649c is short (19.5 days) compared to thermal timescales of the ocean, thus opposing sides of the planet remain permanently warm. For dynamic ocean cases, a faint signal of the double hot spot phenomenon persists, but zonal heat transports by the ocean results in nearly zonally uniform patterns of surface temperature.

For corresponding end-member scenarios of synchronous and resonant modes, global mean surface temperatures could change by as little as 8.9 K for an Earth-like atmosphere with dynamic ocean circulation, or by as much as 28.4 K for a carbon dioxide dominated planet with a slab ocean. For both atmospheric compositions, the climate sensitivity between end-member orbital-rotational states is reduced when considering a full dynamic ocean compared to the slab ocean cases. Perhaps surprisingly, the optically thicker 1 bar carbon dioxide simulations are more sensitive to the orbital-rotational state than are the Earth-like atmosphere. While further investigation is warranted, we surmise that this is due to more vigorous water vapor and cloud feedbacks that occur for our 1 bar carbon dioxide atmospheres, due to their generally warmer climate states compared to the Earth-like atmospheres which remain cold and ice dominated.

However, as noted above, these end-member scenarios may never be fully realised. The relatively rapid changes in orbital-rotational states may instead drive Kepler-1649c into chaotic rotation, which would result in asynchronous planet rotation and not into any favoured resonant state. In such a scenario, the climate patterns of Kepler-1649c become zonally uniform for all combinations of atmospheric composition and ocean assumptions, as star light would irradiate all longitudes of the planet equally with no preferred faces receiving more insolation on time-average than any other.

In all cases simulated, the change in surface temperature caused by the dynamical evolution is relatively small on the substellar hemisphere, but is much larger for the antistellar hemisphere. Kane et al. define substellar and antistellar hemispheres based on the synchronous rotating orbital-rotational reference frame. Naturally, for non-synchronous cases (resonant and chaotic), all sides of the planet receive periodic stellar flux, meaning that the opposing hemisphere is no longer permanently cloaked in darkness and enjoys its moment in the starlight, resulting in considerable warming. The substellar-antistellar dichotomy is most strongly seen for the slab ocean cases. In both synchronous and resonant rotating instances, the inclusion of dynamic ocean circulations effectively transports heat from the day-side to the night-side of the planet, muting the substellar-antistellar dichotomy.

Kane et al.'s 3-D climate simulations reveal that the modulation of orbital-rotational states does not trigger true climate catastrophes for Kepler-1649c (i.e. runaway glaciation or runaway greenhouse states) given the atmospheric compositions explored. If Keperly-1649c oscillates between synchronous and resonant states, significant millennial scale climate changes would occur that are far in excess of anything that the Earth has experienced over such a short timescale. Given the timescales of variation, life as we know it may struggle to adapt to such wide temperature variations indicated by Kane et al.'s simulations. Extant lifeforms may remain confined to the substellar hemisphere due its relative climatic stability compared to the significantly more intense variations that occur on the antistellar hemisphere. However, if instead Kepler-1649c settles into chaotic rotation, millennial scale climate changes would be negligible.

Other factors not explored here may also affect the sensitivity of Kepler-1649c’s climate to oscillations in its orbital-rotational state. Kane et al. have modeled only selected states in equilibrium. If, the orbital-rotational state of Kepler-1649c instead evolves monotonically between our illustrated end-member cases, the climates states may instead gradually smear between synchronous and resonant. Chaotic rotation states may also experience time-evolution of the rotation period. However, differences in non-resonant asynchronous rotation periods are not likely to drive significant changes in climate, compared to toggling into, and out of, synchronous and resonant states.

The atmospheric composition also presents a significant unconstrained parameter. While the atmospheric compositions studied here all result in temperate climates, if the atmosphere of Kepler-1649c were considerably different, the climate sensitivity between states may deviate considerably. For instance, if Kepler-1649c possesses an overwhelmingly dense, cloud and haze filled atmosphere, the surface may receive little direct starlight and become decoupled from changes in the patterns of stellar flux, and surface temperatures may trend towards a globally uniform state regardless of rotation state. Alternatively, if Kepler-1649c instead has a dry, clear, optically transparent atmosphere, the change in rotation state may also cause little difference to the surface temperature as the planet’s energy balance and thus climate would be governed simply by the albedo and emissivity properties of the planet’s solid surface, along with the instantaneous stellar flux input. Finally, our results indicate that the ocean treatment makes a considerable difference in the modulation of surface temperatures between synchronous and resonant states, but makes little difference for the chaotic rotation state. Still, changes in the ocean depth, and the presence and location of continents would significantly influence the results shown here through the modulation of ocean circulations. Jérémy Leconte argued that tidally locked planets with emergent continents will tend to have their continents migrate to the substellar or antistellar points. If continents cluster near the substellar point, recent work with ROCKE-3D suggests that this would mute the effects of ocean circulations. The varied location of continents remains a significant uncertainty in modern terrestrial exoplanet climate modeling, due to their tight coupling to ocean circulations, ocean heat transport, and ultimately climate.

A major focus of exoplanetary studies lies in the exploration of the full diversity of planetary system architectures. However, some of these architectures are similar enough to our own solar system such that they may shed insight into the orbital and atmospheric evolution of our terrestrial planets. Thus, systems that contain terrestrial planets within both the Venus Zone and the Habitable Zone are of particular interest for the investigation of potential Venus/Earth analogs that exist within similar yet different contexts of age, composition, host star type, and orbital dynamics.

The Kepler-1649 system is one of the few known systems to harbour such a potential Venus/Earth analog, where the similar age and, presumably, composition of the planets allows the removal of several free parameters from which other comparative studies may be conducted. Kane et al. have presented the results of dynamical analysis in which they have demonstrated that additional terrestrial planets could exist in orbits that lie between the two known planets. This is not an unlikely scenario given the dynamical packing observed for compact planetary systems and the necessity of coplanarity for transits to be observed. Additionally, the locations of significance resonance, such as the potential 3:2 >mean motion resonance chain, lend credence to a possible planet located between the known b and c planets. Kane et al.'s dynamical simulations further demonstrate that varying eccentricity can be used as a tool to exclude either eccentricities of the known planets or the presence of the hypothetical additional planet. 

Kane et al. also explored the effect of eccentricity on the dynamical evolution of the two known planets, the results of which show that the system is stable up to eccentricities of 0.3 and that the frequency of the eccentricity oscillations can be exceptionally high (1000–3000 years). Such high frequency eccentricity oscillations are a fundamental component of Milankovitch cycles for exoplanets and may play an important role in climate evolution.

Kane et al.'s ROCKE-3D simulations for the climate of Kepler-1649c under various eccentricity and ocean circulation assumptions demonstrate that in fact the surface temperature variation remains relatively small, even for eccentricity at the high end. The region of highest climate stability consistently lies in the substellar hemisphere for all of our simulations, indicating a potential longitude dependence of metabolic processes for any life that may exist on such a planet. Thus, the overall combination of both dynamical and climate evolution studies for this planetary system show that even though their orbits may change far more rapidly than that observed within the Solar System, atmospheric characterisation may yet reveal robust environments that maintain temperate surface conditions. Such dynamically active systems should therefore not be ruled out in the continuing search for potentially habitable worlds.

See also...















Follow Sciency Thoughts on Facebook.

Follow Sciency Thoughts on Twitter.

 

Tuesday, 21 April 2020

Understanding the Earth's Archean atmosphere.

The environment of the Archean Eon from 4 to 2.5 billion years ago has to be understood to appreciate biological, geological, and atmospheric evolution on our planet and Earth-like exoplanets. Its most distinguishing characteristic was negligible oxygen, unlike today’s air, which contains, by dry volume, 21% oxygen, 78% nitrogen, 0.9% argon, and 0.1% other gases. With its radically different atmosphere and lack of macroscopic, multicellular life, the Archean world was alien. However, at that time, the beginnings of modern Earth emerged. For example, Cyanobacteria probably evolved during this period, and these oxygenic photoautotrophs eventually oxygenated the air, setting the stage for later, complex life, including us.

In a review paper published in the journal Science Advances on 26 February 2020, David Catling of the Department of Earth and Space Sciences and cross-campus Astrobiology Program at the University of Washington, and Kevin Zahnle of the Space Sciences Division at the NASA's Ames Research Center, discuss our current understanding of the Archean atmosphere and climate.

The earliest well-preserved sedimentary and volcanic rocks are Archean and provide insights into atmospheric composition, climate, and life. These perspectives are unavailable for the Hadean eon, from about 4.6 to 4 billion years ago, which generally lacks these rocks. For context, the Archean precedes the Proterozoic eon of 2.5 billion to 541 million years ago, and Archean eras provide a timeline for Catling and Zahnle's discussion: the Eoarchean (4 to 3.6 billion years ago), Palaeoarchean (3.6 to 3.2 billion years ago), Mesoarchean (3.2 to 2.8 billion years ago), and Neoarchean (2.8 to 2.5 billion years ago).

The Archean was originally conceived to span the time from after the originof life to the advent of free oxygen. While the origin of life dates back to before 3.5 to 3.8 billion years ago or earlier, newer information puts atmospheric oxygenation after 2.4 billion year ago, inside the Proterozoic. Therefore, considering the Archean in the older sense, the origin of life falls outside this Eon, while the story of oxygen’s rise falls within.

Data about the Archean atmosphere come from how individual gases, or the air as a whole, affected chemical and physical phenomena (e.g., the composition of aerosols, chemical reactions in soils, raindrop terminal velocity, isotopic fractionations, etc.) that were recorded in rocks. So, after a brief discussion of the Hadean, Catling and Zahnle review what the Archean atmosphere was made of. However, because of limited proxy data, major uncertainties remain about the exact levels of atmospheric gases over time.

Catling and Zahnle's discussion also considers evidence for early life and its possible global influence. They assume that metabolically useful gases would have been consumed, while waste gases would have been excreted, as they are today.

Oxygenic photosynthesis produced the most impactful waste gas. The oxygen from Cyanobacterial ancestors flooded the atmosphere rapidly at a time between 2.4 and 2.3 billion years ago, with the transition marked in the rocks by the sudden disappearance of mass-independent fractionation of sulphur isotopes. The Great Oxidation Event thus began, which ended about 2.1 to 2.0 billion years ago. Although this switch to an oxygenated atmosphere and shallow ocean occurred in the Palaeoproterozoic era (2.5 to 1.6 billion years ago), the weakly reducing atmosphere that was eliminated typified the Archean. In this context, 'weakly reducing' means minor levels of reducing gases, such as carbon monoxide, hydrogen, and methane, in an anoxic atmosphere of bulk oxidized gases, carbon dioxide and nitrogen. A major outstanding question concerns how trends of biological and geological evolution relate to the Great Oxidation Event.

Atmospheric composition, in turn, affected Archean climate. At 4 billion years ago, solar luminosity was 25 to 30% lower than today, but Archean Earth was not persistently frozen because abundant evidence shows an active hydrological cycle. Liquid water under a fainter Sun likely implies more abundant greenhouse gases than today. 

In addition to a gradual increase in solar luminosity, slow changes in the solid Earth over time provided boundary conditions for atmospheric evolution. On geological time scales, volcanic and metamorphic gases replenish atmospheric volatiles that escape to space or are chemically sequestered into solid materials.

The composition of the Hadean atmosphere is obscured by a lack of well-preserved rocks, but analysis of zircons, crystals of zirconium silicate, suggests that continents, oceans, and perhaps life all originated in the Hadean. Zircons are tiny (less thab 0.5 mm) durable pieces of continental crust. Elevated oxygen¹⁸/oxygen¹⁶ in 4.3-billion-year-old zircons was possibly inherited from oxygen¹⁸-enriched, weathered surface rocks that were later buried and melted, which implies the presence of surficial liquid water and even land. In addition, graphite inside a 4.1-billion-year-old zircon has a biogenic-like carbon¹³/carbon¹² ratio of −24 per mil compared to a standard reference material, although lack of context means that an abiotic origin cannot be eliminated.

Impact bombardment would have affected Hadean and subsequent Archean environments. The lunar record implies a decay of terrestrial impact bombardment extending into the Archean. The estimated median age of the last impact big enough to vaporize the entire ocean is about 4.3 billion years ago, which provides a crude upper age limit on the origin of life. An origin of life during the period 4.3 to 4.0 billion years ago is consistent with phylogenetic inferences. Later, the Late Heavy Bombardment is a hypothesized interval of enhanced bombardment superposed on the general decline, which occurred between 4.2 and 4.0 to about 3.5 billion years ago, based on the ages of lunar rocks and meteorite shocks, althoughmany dispute that the Late Heavy Bombardment was a discrete event. Regardless, an Late Heavy Bombardment would likely not sterilise Earth; any microbial life would have rebounded.

The mantle contains excess highly siderophile elements (elements that will dissolve readily in iron either as solid solutions or in the molten state, which tend to sink into the core) relative to concentrations expected after Earth’s iron core formed, which removed highly siderophile elements. Similarities of isotopes and relative proportions of these highly siderophile elements to those in enstatite chondrite and achondrite meteorites suggest that this highly reducing meteoritic material was delivered late in Earth’s accretion. Thus, carbon and nitrogen were supplied in graphite and nitrides. Therefore, the Hadean atmosphere and mantle were probably initially highly reducing, before subsequent oxidation either by hydrogen escape or disproportionation of mantle iron oxide accompanied by iron loss to the core. In any case, iron-cored impactors would reduce seawater to hydrogen and create transient, highly reducing atmospheres that may have been important for the origin of life. 

Because a liquid ocean likely existed by about 4.4 billion years ago, feedbacks in the geologic carbon cycle probably stabilised the long-term climate. However, consumption of carbon dioxide in the weathering of impact ejecta by carbonic acid suggests a cool early Hadean surface near 0°C under the faint Sun.

Estimates of nitrogen bound in today’s solid Earth range a few to about 40 bar equivalent, which allows for considerable nitrogen in the Hadean atmosphere unless nitrogen was incorporated into a reducing, deep layer of magma, a magma ocean, formed after the Moon-forming impact about 4.5 billion years ago. So, although nitrogen was one of the bulk atmospheric gases, its Hadean level, higher than today or lower, remains unclear. 

Precambrian events and atmospheric change. Catling & Zahnle (2020).

Proxies constrain Archean atmospheric composition. Gases reacted with the seafloor or land, leaving chemical traces in seafloor minerals or in soils that became palaeosols. In addition, atmospheric particles carried isotopic signatures into sediments that were diagnostic of atmospheric composition. Occasionally, fluid inclusions in rocks trapped seawater with dissolved air or even microbial gases.

Sometimes the physical environment affected rocks and minerals. Their preservation allows estimates of environmental temperature and barometric pressure. 

The strongest constraint on Archean atmospheric composition is that the ground-level mixing ratio of oxygen was less than one millionth of present atmospheric level or less thab 0.2–parts per million by volume oxygen for air of 1 bar, indicated by the presence of sulphur isotope mass-independent fractionation in Archean sedimentary minerals. An oft-quoted limit of one hundred thousandth  present atmospheric oxygen derives from an earlier photochemical model that could not address oxygen levels lower than this. Usually, isotope fractionation is proportional to the mass difference between isotopes; e.g., in diffusive separation of sulfur-containing gases, sulphur³⁴ becomes about half as abundant, relative to sulphur³², as sulphur³³. However, some particular photochemical reactions produce mass-independent fractionation that, by definition, deviates from proportionality to mass.

Archean sulphur isotope mass-independent fractionation is tied to the production of elemental sulphur from photochemistry in anoxic air; in contrast, in an oxic atmosphere, sulphur isotope mass-independent fractionation nearly disappears. Photochemistry imparts sulphur isotope mass-independent fractionation and produces elemental sulphur, starting with reactions that photolyse volcanic sulphur dioxide. This photolysis occurs when short-wavelength ultraviolet penetrates Earth’s troposphere in the absence of a stratospheric ozone layer. Scarce ozone implies negligible molecular oxygen, from which ozone derives.

In anoxic air, sulfur ends up in insoluble octasulphur aerosols and watersoluble sulphate and sulphur dioxide, unlike today’s atmospheric sulphur, which almost entirely oxidises to sulphate. In the Archean, as well as anoxia, gases such as methane or hydrogen produced sufficiently reducing conditions that sulphur gases persisted and polymerized into octasulphur. Reactions that dominantly imparted sulphur isotope mass-independent fractionation are debated: They include polysulphur formation, sulphur dioxide photolysis, and other reactions. In any case, when octasulphur particles fell to Earth’s surface, their sulphur isotope mass-independent fractionation isotopic composition complemented that of sulphate particles, allowing preservation of different sulphur isotope mass-independent fractionation signs in these phases as sedimentary pyrite (iron sulphide) and barite (barium sulphate). Sulphate, of course, could later be microbially transformed to pyrite.

Numerous redox-sensitive tracers corroborate negligible Archean oxygen. Anoxia proxies that have long been recognised include the lack of pre-Great Oxidation Event continental sediments stained by red ferric oxides (redbeds), detrital grains from well-aerated rivers of siderite (iron carbonate), uraninite (uranium dioxide), or pyrite that would oxidise and dissolve or rust at high oxygen pressures, and palaeosols with iron washed out by anoxic rainwater. Furthermore, Archean marine sediments have low concentrations of elements that enter rivers during oxidative continental weathering. Conversely, glacial sediments contain continental materials lacking oxidative weathering loss of molybdenum.

Iron formations, which are marine chemical sedimentary rocks rich in iron and silica (15 to 40 percent by weight iron and 40 to 60 percent by weight silicon dioxide), indicate that the deep Archean ocean contained aqueous (disoloved) iron and so was anoxic. In 'Superior-type iron formations' that formed near-shore, rare earth elements show that dissolved iron was partly sourced from seafloor hydrothermal vents and upwelled onto continental shelves where the iron precipitated. In today’s deeply oxygenated oceans, oxidised iron instead precipitates locally around vents. Archean shallow-water iron formations constrain atmospheric oxygen at less than 0.0024 bar.

The absence or presence ofmass-dependent fractionation of various isotopes can also indicate anoxic versus oxic conditions. In the case of sulphur, today’s oxidative weathering of continental sulphides produces soluble sulphate, which rivers carry to the ocean. When Bacteria reduce sulphate to pyrite in seafloor sediments, they impart mass-dependent sulphur isotope fractionation if sulphate is present at sufficient concentrations, as in the modern oceans, but usually not in Archean seawater. The absence of this isotope fractionation indicates little Archean seawater sulphate and implies anoxic air. Using similar arguments of oxygen-sensitive weathering, transport, and fractionation, isotopes of copper, chromium, iron, uranium, molybdenum, and selenium, indicates an anoxic Archean atmosphere.

Atmospheric oxidation of 2.7-billion-year-old iron-nickel micrometeorites has been used to argue for oxygen near-modern levels above about 75 km altitude. However, given copious evidence for an anoxic Archean atmosphere, an alternative explanation is that high carbon dioxide levels (perhaps over 70%) oxidised the micrometeorites.

Nonetheless, even under a globally anoxic atmosphere, lakes and shallow seawater inhabited by oxygenic photosynthesisers could have become 'oxygen oases', local or regional areas with elevated Oxygen, although estimates for disolved oxygen levels in Archean oxygen oases only range from 0.4 to 7% of present disolved oxygen levels in surface seawater.

When exactly oxygenic photosynthesis began and dominated over anoxygenic photosynthesis is debated, but signs of biological carbon fixation appear early. Graphite in a roughly 3.7-billion-year-old outcrop of sedimentary rock is carbon¹²-enriched in amounts typical of photosynthetic microbes. Then, at 3.52 billion years ago, the proportion of carbon¹³ in kerogen and associated marine carbonate of found in Australia is similar to biological isotope fractionation in modern oceans.

Fossil evidence for Archean Cyanobacteria has been reported. Light organic carbon isotopes and structures like those made by filamentous Cyanobacteria found within Stromatolites or other microbially induced sedimentary structures are consistent with Cyanobacteria by 3.2 to 2.7 billion years ago. Cyanobacteria could be corroborated by biomarkers, which are remnant organic molecules from particular organisms. However, putative Archean biomarkers have been plagued by younger contamination.

Instead, Cyanobacterial interpretations are strengthened by geochemical data suggesting oxygen oases at 3.2 to 3.0 billion years ago. Fractionated iron and molybdenum isotopes and levels of redox-sensitive metals suggest marine photic zone oxygen (chromium isotope data have been used to argue for the existence of approximately 3-billion-year-old terrestrial oxygen, but they were probably caused by modern oxidative weathering).

Then, by 2.8 to 2.6 billion years ago, increasing concentrations and isotopic fractionation of molybdenum and sulphur in marine shales suggest that oxygen proximal to Cyanobacterial mats and Stromatolites on land oxidised sulphides and boosted sulphate and molybdenum riverine fluxes to the oceans. These trends are consistent with isotopic evidence for Neoarchean methanotrophy and oxidative nitrogen cycling.

Concentration spikes at 2.5 to 2.66 billion years ago in molybdenum, selenium, and rhenium and isotopic excursions of molybdenum, selenium, uranium, and nitrogen have been interpreted as arising from oxygen transients or 'whiffs' of oxygen. Critics argue that the data derive from post-Great Oxidation Event alteration. Alternatively, oxygen oases of Cyanobacteria within soils or lakebeds may have mobilised these elements into rivers and then the sea.

Phylogenetic analyses mostly suggest that Cyanobacteria originated by 2.8 billion years ago. Molecular clocks must be calibrated by physical evidence, and phylogenetic methods are themselves debated. While some argue that oxygenic photosynthesis evolved only 50 to 100 million years before the Great Oxidation Event, most studies suggest an earlier Palaeoarchean or Mesoarchean age.

Amajor unresolved issue is how the Great Oxidation Event was related to underlying geological or biological trends. Life on its own cannot change the net redox state of the global environment because each biological oxidant is complemented by a reductant. In oxygenic photosynthesis, a mole of organic carbon accompanies every mole of molecular oxygen.

Consequently, for long-term molecular oxygen accumulation, some organic carbon must be segregated from the oxygen and buried. Alternatively, oxygen molecules are liberated if microbes use the organic carbon to make other reductants, such as sulphides from sulphates, that are buried. However, Archean seawater sulphate concentrations were small, so organic carbon burial is the oxyen flux that matters for the Great Oxidation Event. Atmospheric oxidation also occurs when hydrogen escapes to space after the photochemical breakdown of gases such as molecular hydrogen and methane, which ultimately derive from water and are relatively abundant in anoxic air.

Schematic history of atmospheric oxygen. Colored arrows represent known oxygen constraints, but the black line is speculative. An Archean upper bound of under 0.2-mbar oxygen (blue) is for photochemistry that generates octasulphur aerosols, preserving observed mass-independent isotope fractionation in sulphur compounds. The size and shape of an oxygen overshoot during the Great Oxidation Event are highly uncertain; a lower bound (red arrow) comes from iodine incorporation into carbonates. In the Proterozoic, a lower bound (light green) of 0.006 bar is required for an oxygen-rich atmosphere to be photochemically stable. However, oxygen levels likely remained low for most of the Proterozoic. Neoproterozoic oxygenation began around 800 million years ago. From about 600 million years ago, a lower bound of over 0.02-bar oxygen (dark green) is from plausible oxygen demands of macroscopic Ediacaran and Cambrian biota. Charcoal since 400 million years ago ago implies a lower bound of over 0.15 bar (purple). The post-Devonian black line for oxygen evolution approximately represents curves from calculations of carbon and sulphur isotopic mass balance. Catling & Zahnle (2020).

Atmospheric oxygen is determined by net redox fluxes into and out of the atmosphere. This simple truth is deceptive because these redox fluxes are themselves controlled by less easily constrained oxidative weathering (both seafloor and continental) and volcanic and metamorphic degassing, as well as hydrogen escape to space. Permanent atmospheric oxygenation requires an oxygen level of 0.03 x present atmospheric level or more to prevent a destabilising positive feedback of photochemical destruction of tropospheric oxygen that otherwise occurs when an incipient ozone column is still transparent to far ultraviolet light.

These oxygen levels would have only been attained when the oxygen flux from the burial of organic carbon exceeded the kinetically efficient sink from oxygen-consuming gases (carbon monoxide, hydrogen, hydrogen sulphide, and sulphur dioxide) from volcanism and metamorphism plus fluxes of reducing cations such as iron²⁺ from seafloor vents. A minority assume that such a flux imbalance applied as soon as oxygenic photosynthesis evolved, mandating a rapid rise of atmospheric oxygen. In the consensus assessment that oxygenic photosynthesis evolved long before the Great Oxidation Event, efficient consumption of oxygen initially suppressed oxygen levels.

Ideas favoring increased oxygen fluxes from organic burial appeal to more continental shelf area available for burial, more phosphorus to stimulate photosynthesis, or subduction of organic carbon relative to ferric iron. Because organic carbon burial extracts carbon¹² and leaves inorganic carbonates carbon¹² depleted, it is difficult to reconcile these hypotheses with the remarkable constancy of the carbon isotope record, which indicates little change in average organic burial rates between the Archean and Proterozoic. However, an increase in organic burial might have occurred if negligible oxidative weathering of carbon¹²-rich organics on land or a sink of seafloor carbonate meant that the operation of the carbon cycle or the isotopic composition of carbon input into the surface environment differed from today.

Hypotheses for a slowly decreasing oxygen sink to the Great Oxidation Event tipping point rely on a decline in the ratio of reduced-to-oxidised species from volcanic, metamorphic, and hydrothermal sources. Some emphasize the role of hydrogen escape to space in oxidising solid Earth, lowering Earth’s capacity to release oxygen-consuming reductants. New evidence from xenon isotopes supports rapid Archean hydrogen escape.

Molecular nitrogen dominates today’s atmosphere, and three lines of evidence have begun to constrain Archean nitrogen levels. First, the largest size of 2.7-billion-year-old fossil raindrop imprints provides a conservative limit of paleopressure of under 2.1 bar and a probable limit of under 0.52 to 1.2 bar. Second, the molecular nitrogen/argon³⁶ ratio in fluid inclusions indicates a nirogen presure of less than 1.0 bar at 3.3 billion years ago and less than 1.1 bar at 3.5 to 3.0 billion years ago. Third, vesicle volumes in 2.7-billion-year-old basaltic lava flows erupted at sea level imply a 0.23 bar palaeopressure.

The inferred history of nitrogen presure depends on how the geologic nitrogen cycle has changed over time. Nitrogen can only greatly accumulate as atmospheric nitrogen or in rocks as ammonium, amide, nitride, or organic nitrogen.Under typical mantle temperatures and redox conditions, volcanic gases contain molecular nitrogen, not ammonia, and because molecular nitrogen is unreactive, it enters the air. Today, molecular nitrogen is also produced when oxidative weathering of organic matter on the continents makes nitrate that undergoes rapid biological denitrification into molecular nitrogen. Within large uncertainties, volcanic and oxidative weathering inputs of atmospheric nitrogen are comparable. Using sedimentary carbon/nitrogen data, it has been argued that the sum of these nitrogen sources was balanced over the Phanerozoic primarily by nitrogen burial in organic matter, so that the Phanerozoic partial pressure of nitrogen varied little.

Another theory states that low Archean palaeopressure arose because today’s long-term nitrogen atmospheric input from oxidative weathering and denitrification was absent. If so, nitrogen presure would have risen at the Great Oxidation Event, and nitrogen in today’s air must have been in solid phases previously. Certainly on modern Earth, nitrification, by humankind’s addition of nitrate to land and sea, has enhanced denitrification, indicated by increased atmospheric nitrous oxide. On the other hand, a model can be constructed where nitrogen presure declines after the Great Oxidation Event if burial of organic nitrogen increased. 

In the Hadean, nitrogen presure either started high and then diminished or was initially low if nitrogen partitioned into a very reducing magma ocean. However, low nitrogen/arbon in today’s mid-ocean ridge source basalts suggests considerable nitrogen degassing once the upper mantle became oxidised because then nitrogen became insoluble in magmas and upper mantle fluids. Even today, some of the upper mantle lies within the stability field for ammonium, so that increased oxidation of the early mantle and mantle wedge could have caused more subducted nitrogen to outgas as molecular nitrogen.

Marine phyllosilicates (sheet silicate minerals) at 3.8 billion years ago are ammonium enriched, which probably came from porewater ammonium derived from degraded organics, and these data have been used to argue that Archean nitrogen was sequestered into solid phases after an early advent of biological nitrogen fixation. In the early ocean, ammonium would have been the stable form of dissolved nitrogen unlike today’s nitrate. Consequently, a postulated drawdown of Archean nitrogen involves biological fixation, organic burial, and subduction of ammonium in refractory minerals. The rate of organic burial must have been relatively high for a time-integrated loss to affect nitrogen presure significantly, which is not necessarily inconsistent with carbon isotopic constraints because early high degassing of carbon required more carbon to be buried.

Nitrogen isotopes appear to confirm that biologically fixed nitrogen entered the Archean mantle. Sedimentary organics have 7 parts per thousand more nitrogen¹⁵ as a proportion of total nitrogen; compared to modern atmospheric nitrogen, whereas mantle nitrogan has 2 parts per thousand less nitrogen¹⁵ than modern atmospheric nitrogen. Fractionation mostly arises when denitrification preferentially converts nitrate or nitrite nitrogen¹⁴ into gaseous molecular nitrogen. Thus, heavy nitrogen¹⁵ in 3.1 to 3.5-billion-year-old mantle-derived diamonds may be a sedimentary component.

Ammonium substitutes for potassium, and breakdown of previously subducted ammonium-containing minerals in magmas at oceanic islands releases moleciular nitrogen and radiogenic Argon⁴⁰ derived from Potasium⁴⁰. The ratio Argon⁴⁰/moleciular nitrogen in plume-related lavas scatters by a factor of about 4 to 5, and higher values (older from more Argon⁴⁰) correlate with smaller, Archean-like values of nitrogen¹⁵ ratios, consistent with a history of ammonium subduction. Because moleciular nitrogen is uncorrelated with non-radiogenic Argon³⁸ or Argon³⁶, nitrogen in the current mantle is not primordial but recycled. 

Some models have used nitrogen ingassing versus outgassing fluxes to infer past atmospheric nitrogen presure, however estimates of modern nitrogen outgassing vary a great deal, making such models dubious.

Unlike molecular nitrogen, other nitrogen-bearing Archean gases would have been present only in trace quantities. With only tiny atmospheric fluxes of nitrate or nitrite, nitrous oxide from denitrification would have been negligible, perhaps restricted to lakes. Lightning production of nitric oxide, followed by hydrogen addition and nitroxyl dissolution and decomposition, might maintain ground-level nitrous oxide to a few parts per billion by volume, compared to a pre-industrial levels of about 270 parts per billion nitrous oxide. The other nitrogen oxides (nitric oxide and nitrogen dioxide) would have been at trace levels because their lightning production in carbon dioxide-nitrogen air is inefficient.

Ammonia levels of 10 to 100 parts per billion would provide a greenhouse effect to counteract the faint young Sun, but these levels of ammonia cannot be sustained against ultra violet photolysis. Possibly, a stratospheric organic haze, such as that on Saturn’s moon, Titan, shielded tropospheric ammonia from ultra violet light. However, whether this shielding actually occurred depends on whether hazes actually existed and on the size and radiative properties of haze particles, which remain uncertain.

Another nitrogen-bearing gas, hydrogen cyanide is more stable photochemically than ammonia and made in reducing atmospheres by lightning, impacts, or ultra violet-driven photochemistry. In particular, nitrogen atoms from molecular nitrogen photolysis in the upper atmosphere can mix to lower levels and react with methane photolysis products to make hydrogen cyanide. At Archean biogenic methane levels of about 3000 parts per million, hydrogen cyanide concentrations reach 100 parts per million.

Since the early Hadean, carbon dioxide has probably always been Earth’s most important non-condensable greenhouse gas. Carbon dioxide also affects seawater pH and influences the carbon cycle through the formation of carbonates and organic matter. However, direct evidence for Archean carb dioxide on levels remains scanty. 

Palaeosols provide some estimates. Acid leaching in Archean soils arose from carbon dioxide dissolved in rainwater. Mass-balance calculations give 10 to 50 times present atmospheric levels of carbon dioxide at 2.7 billion years ago and 8 to 69 times present atmospheric levels of carbon dioxide at billion years ago. However, these analyses assume that all the carbon dioxide that entered the soils caused dissolution, so carbon dioxide pressure could have been higher if only a fraction of the carbon dioxide had been used. Another study used an analysis with an estimate of the composition of temperature-dependent aqueous solutions during weathering and, because of a weaker dependence of weathering on carbon dioxide pressure, obtained higher carbon dioxide pressure of 85 to 510 times present atmospheric levels at 2.77 billion years ago, 78 to 2500 times present atmospheric levels at 2.75 billion years ago, and 160 to 490 times present atmospheric levels at 2.46 billion years ago.

High Archean carbon dioxide pressure does not have to induce acidic seawater and dissolve marine carbonates. Instead, an increase in calcium ion concentrations could maintain an ocean saturated in calcium carbonate at alkaline pH; alternatively, seawater pH could be slightly lower than today, but calcium carbonate would remain saturated. In fact, sedimentary marine carbonates appear from 3.52 billion years ago onward.

Calcium isotopes might provide insight into coupled Archean seawater pH and carbon dioxide pressure. These two variables and carbonate alkalinity define a system where any two variables imply the third. In evaporating seawater, calcium isotopes could undergo Rayleigh distillation (isotope fractionation) under alkaline conditions, but limestones from 2.6-billion-year-old the Campbellrand marine evaporites of South Africa show no spread in calcium isotope values, which might imply a pH of 6.4 to 7.4 for a Neoarchean sea.

Siderite (iron carbonate) in Archean iron formations has been proposed as a carbon dioxide pressure proxy if it precipitated in equilibrium with the atmosphere. However, data suggest that this siderite was diagenetic (formed within the rocks at some time after they were deposited) so Catling and Zahnle omit this proxy.

Potentially, two negative feedbacks control long-term carbon dioxide pressure. First, the net consumption of carbon dioxide in acid weathering of continents or the seafloor ends up making carbonates. Seafloor weathering occurs when water in the permeable abyssal plains dissolves basaltic minerals, releasing calcium ions and precipitating calcium carbonate in veins and pores. Second, 'reverse weathering' reactions have been proposed. This reaction consumes aqueous silicon dioxide and uses cations to make aluminosilicate clays instead of carbonates that consume carbon, so that carbon dioxide stays in the atmosphere. High dissolved seawater silica and pH are hypothesised to promote reverse weathering and create clay minerals.

If atmospheric carbon dioxide pressure is high, pH is low, and temperature is warm, the formation of carbonates is favored, and carbon dioxide pressure falls, in negative feedback on carbon dioxide pressure and climate. However, if carbon dioxide pressure is low and pH is high, reverse weathering may be an alternative to carbonate formation and a negative feedback on low carbon dioxide pressure.

Estimates of today’s reverse weathering flux vary by an order of magnitude, the rate coefficient (the rate at which a reaction occurs) for reverse weathering reactions varies by many orders of magnitude, and the solubility of authigenic phyllosilicates (needed to calculate reverse weathering) varies by over an order of magnitude. Consequently, a self-consistent, coupled carbon-silica cycle model since 4 billion years ago shows that reverse weathering can be important or unimportant for the Proterozoic climate depending on parameter choice, while reverse weathering in the Archean is muted because of the probably lower land fraction and sedimentation rate. Considering these factors, estimates of Archean carbon dioxide pressure and seawater pH that Catling and Zahnle present are based on carbonate-silicate cycle models without reverse weathering.

In the Archean, with greater seafloor production than today, seafloor weathering could have been comparable to continental weathering, and negative feedback likely maintained average Archean surface temperatures between 0° and 40°C with seawater pH 6.4 to 7.4. The corresponding carbon dioxide pressure would have been 0.006 to 0.6 bar 4 billion years ago, assuming that about 10 000 parts per million by volume methane also contributed to the greenhouse effect.

Anoxic Archean air could hold 1000 s of parts per million by volume of methane if a microbial flux of methane was comparable to today’s. Phylogenetically, methanogens date back to more than 3.5 billion years ago. In contrast, the modern oxygenated atmosphere destroys reducing gases rapidly, limiting tropospheric methane and hydrogen abundances to 1.8 and 0.55 parts per million by volume, respectively.

Evidence points to high levels of Archean methane. First, signs of methanogens and methanotrophs from light carbon isotopes in Archean organics imply methane’s presence. Second, Archean sulphur isotope mass-independent fractionation requires more than 20 parts per million by volume methane  to generate particulate sulphur. Third, the deuterium-to-hydrogen ratio of 3.7-billion-year-old seawater estimated from serpentine minerals is 2.5% lighter than today, which could be explained by rapid Archean escape of hydrogen and isotopic fractionation. This hydrogen was likely derived from ultra violet photolysis of methane in the upper atmosphere.  Fourth, globally extensive glaciations during the Great Oxidation Event provide circumstantial evidence for high Archean methane. At the tipping point, air flips from anoxic to oxic in only about 10 000 years, causing a roughly 10°C temperature drop by oxidizing methane. This chemical transition is far faster than the 100 000-to-1 000 000-year response of the carbonate-silicate thermostat.

Another carbon-containing gas, carbon monoxide, was probably not abundant in the presence of an Archean microbial biosphere. The Last Universal Common Ancestor was likely capable of anaerobic carbon monoxide consumption, which involves water as a substrate and catalysts such as iron sulphide that were probably widespread.With this reaction, microbes would draw carbon monoxide down to 1 to 100 parts per million by volume. However, episodic carbon monoxide levels at much higher levels may have occurred when large impacts delivered cometary carbon monoxide ice or organic matter that was oxidized.

Because of relatively abundant methane, a high-altitude Archean organic haze might have formed. This idea was first suggested in the 1980s. Empirically, if the methane-carbon dioxide ratio exceeds about 10% in an ultra violet-irradiated carbon dioxide-nitrogen-methane mixture, radicals from methane photolysis polymerise into organic particles. 

When and whether an organic haze formed are uncertain. A haze could have affected tropospheric sulfur gases by blocking ultra violet photons. Consequently, the structure of sulphur isotope mass-independent fractionation variations in Archean sedimentary minerals and their correlation with light, organic carbon¹³ have been attributed to episodic hazes driven by variable atmospheric methane-carbon dioxide ratios. However, given only a hazy understanding of which species and reactions are important for sulphur isotope mass-independent fractionation, interpretations of episodic hazes are speculative rather than definitive.

The Archean lower atmosphere is unlikely to have been hydrogen-rich given the antiquity of methanogens, ome of which convert hydrogen into methane. Anoxygenic photosynthesis also consumes hydrogen. In models with methanogens and hydrogen-based photosynthesisers, atmospheric hydrogen mixing ratios depend on assumed hydrogen outgassing and biological productivity but generally are no more than 0.0001. These levels preclude hydrogen as an important Archean greenhouse gas. Detrital magnetite carried in rivers 3.0 to 2.7 billion years ago would have dissolved at high hydrogen pressures via microbial reduction of iron to its soluble form using hydrogen, providing an upper limit of atmospheric hydrogen pressure of  0.01 bar.

Changes in atmospheric xenon isotopes through the Archean stop after the Great Oxidation Event (analogous to sulphur isotope mass-independent fractionation) and potentially tell us about ancient oxygen, methane, and hydrogen levels, including in the otherwise hidden Hadean. The trend also relates to how much hydrogen escaped from Earth and hence Earth’s total oxidation over time.

Xenon in fluid inclusions in Archean rocks becomes isotopically heavier through the Archean relative to an initial solar composition until the fractionation reaches that of modern air around 2.1 to 2.0 billion years ago. Xenon dragged out into space by escaping hydrogen during the Archean and Hadean best explains the progressive mass fractionation. An alternative explanation from trapping of positively charged xenon ions in organic hazes has the problem that weathering or microbial processing of buried organics would release xenon and modulate the xenon isotopes in post-Archean air, which is not observed. Today, the nine atmospheric xenon isotopes (124 to 136) have a huge fractionation of 4.2% per atomic mass unit relative to solar or chondritic sources, whereas the six, lighter krypton isotopes (76 to 86) are barely fractionated. Earth’s xenon/krypton ratio is also 4 to 20 times less than meteorites, implying selective xenon loss.

Unlike lighter krypton, xenon is easily ionized by solar ultra violet or charge exchange with positive hydrogen ions, so positive xenon ions can be dragged out to space by escaping positive hydrogen ions. Whereas positive xenon ions are unreactive with atomic hydrogen, molecular hydrogen, or carbon dioxide, any positive krypton ions are neutralised by reactions with molecular hydrogen, explaining the lack of krypton isotope fractionation. Ions are tethered to Earth’s magnetic field lines, but a 'polar wind' of hydrogen ions escapes along open field lines at the poles, accounting for about 15% of all hydrogen that escapes today. Positive xenon ions could be dragged by a vigorous ancient polar wind. That requires copious hydrogen to be derived from relatively abundant methane and/or molecular hydrogen in the lower atmosphere. Ultra violet decomposes methane in an anoxic upper atmosphere, releasing hydrogen.

The hypothesized xenon escape works only in anoxic air, so, like sulphur isotope mass-independent fractionation, xenon isotopes record the Great Oxidation Event. Oxic air destroys molecular hydrogen and methane, making their abundances too low to supply enough hydrogen to drag along xenon. In addition, molecular would remove positive xenon ions in a resonant charge exchange reaction; the oxygen molecule can donate an electron to the xenon ion, creating a positively charged oxygen molecule and a xenon atom (electrons have a negative charge).

Catling and Zahnle's preliminary model finds that xenon can escape when the total hydrogen mixing ratio exceeds about 1% for the solar extreme ultra violet flux expected around 3.5 billion years ago. If nearly all Archean hydrogen was biologically converted into methane, Catling and Zahnle we interpret the xenon escape constraint of more than 1% of the hydrogen mixing ratio to be more than 0.5% methane, which is a rare empirical constraint on Archean methane.

The inferred escape of a strong reductant, hydrogen, would oxidise the entire Earth. The total oxidation during the Archean is equivalent to oxygen from a tenth of more of an ocean.

Sulphur isotope mass-independent fractionation proves that Archean sulphur-containing gases existed but tells us little about their concentrations. All sulphur gases apart from octasulphur, which condenses, are susceptible to ultra violet photolysis and are short lived and interconverting. Consequently, atmospheric sulphur is sequestered into stable sulphuric acid aerosols by oxidation or octasulphur aerosols by reduction. For a wide range of Archean atmospheric redox conditions, both types of aerosol should form and precipitate.

Although sulfur gases absorb ultra violet, they probably did not shield Earth’s surface. Surface temperatures greater than 50°C would be required to produce enough atmospheric octasulphur vapour to shield surface life from ultra violet radiation.

Related to atmospheric composition, is another basic problem: How the Archean climate remained clement under a fainter Sun. In the Sun’s core, nuclear reactions fuse four protons into helium nuclei, increasing the mean molecularmass and decreasing the number of particles per unit volume. In response, the weight of the overlying column presses inward, and the core temperature rises. A greater radial temperature gradient drives more outward radiation flux, so the Sun brightens. Solar evolution models show that the Sun was 25 to 30% fainter at 4 billion years ago.

One solution to the Faint Young Sun problem is that the young Sun was not faint but more massive and therefore as bright as today. The Sun would then need to lose enough mass over time so that a declining weight of the column above the core matched the pressure loss from particles fused in the core. However, observations are unsupportive: Fast mass loss only happens in the first 200 million years after Sun-like stars form.

As well as solar luminosity, Earth’s mean global temperature depends on the Bond albedo (the reflectivity over all wavelengths) and greenhouse effect. Some Faint Young Sun hypotheses invoke clouds: that the Archean Bond albedo was low because of substantially different cloud cover or that abundant high-altitude cirrus clouds warmed early Earth. However, no plausible combination of few low clouds (which tend to reflect sunlight) and many high, icy clouds (which tend to warm the surface with downwelling infrared) solves the Faint Young Sun problem. Three-dimensional climate models with a lower solar flux produce weaker evaporation and fewer low clouds, but a strong greenhouse effect is still required because the albedo decrease is small.

The most plausible solution to the Faint Young Sun is a big greenhouse effect. Water vapor is an important greenhouse gas but, on its own, cannot solve the Faint Young Sun problem. Water vapor condenses into rain and snow, so its abundance is limited by the saturation vapour pressure that depends only on temperature, which is set by non-condensable greenhouse gases, such as carbon dioxide and methane. 

Consequently, substantial carbon dioxide is the most obvious driver of an enhanced Archean greenhouse effect. Catling and Zahnle expect Archean carbon dioxide pressure to have been high because of Earth’s carbonate-silicate cycle thermostat. If global temperatures and carbon dioxide are low, carbon dioxide removal via rainfall is slow, as are rates of continental and seafloor silicate weathering. Then, geological emissions of carbon dioxide raise atmospheric carbon dioxide levels, increasing global temperatures. Conversely, if the climate warms too much, greater rainfall and silicate weathering consume carbon dioxide and cool the Earth.

These negative feedbacks would have moderated the Archean climate to amean global temperature of 0° to 40°C (or 0° to 50°C without land and only seafloor weathering), carbon dioxide alone could have solved the Faint Young Sun problemwith 0.004 to 0.03 bar carbon dioxide at 2.5 billion years ago and 0.024 to 1 bar carbon dioxide at 4 billion years ago, where the spread comes from uncertain carbonate-silicate model parameters.

However, without oxygen in the atmosphere, a methane greenhouse bears consideration. Methane levels of 103 to 104 parts per million by volume produce roughly 10° to 15°C of Archean greenhouse warming, which lowers the carbon dioxicde pressure needed to warm the Archean Earth. Some of this warming comes from a few parts per million by volume of ethane, which is derived from methane.

Both substantial methane and carbon dioxide may be necessary if lower bounds on nitrogen pressure are valid. Although nitrogen itself is not an effective greenhouse gas, it pressure broadens line and continuum infrared absorption, enhancing the greenhouse effect. If Archean nitrogen pressure was about half of today, a few degrees would be lost from the greenhouse effect. 

However, too much methane relative to carbon dioxide would create the high altitude organic haze, which may cool Earth by up about 20°C by absorbing incoming solar radiation and radiating energy back to space in a so-called 'anti-greenhouse effect'. A cooling limit exists because as a haze becomes more ultra violet absorbing, it shields methane molecules from the photolysis needed for further haze formation. Thus, attenuation of sunlight saturates.

An organic haze might protect tropospheric ammonia from ultra violet, but whether ammonia was a viable greenhouse gas is unclear. The composition of organic aerosols in the Archean where sulphur was present and the carbon/oxygen ratio was less than 1 was different from particles on Titan where sulphur is absent and the carbon/oxygen ratio is more than 1. Consequently, the ultra violet attenuation properties of early Earth haze particles remain uncertain.

Another proposed early greenhouse gas is molecular hydrogen, if it attained percentage abundances. Molecular collisions produce temporary charge separation, allowing hydrogen to undergo collision-induced absorption of infrared photons. Archean hydrogen is unlikely to have been abundant. Instead, hydrogen may be relevant for the Hadean greenhouse, if percentage hydrogen levels arose from hydrogen outgassing or impacts.

Sulphur gases were likely unimportant for warming. Sulphur dioxide is a strong greenhouse gas, but tropospheric sulphur dioxide dissolves in rainwater, and Archean stratospheric sulphur dioxide photochemically disproportionates into ultra violet-stable polysulphur and sulphuric acid aerosols. The latter raises the albedo and more than offsets any sulphur dioxide greenhouse warming. Hydrogen sulphide is a weak greenhouse gas that overlaps in the absorption with of methane and is ultra violet fragile. Carbonyl sulphide has been suggested as an Archean greenhouse gas because ultra violet photolysis of carbonyl sulphide was posited to explain observed relationships in sulphur isotope mass-independent fractionation. Carbonyl sulphide persists today because it is not oxidised by hydroxide radicals that cleanse the troposphere of other pollutants. However, in the absence of an ozone shield, carbonyl sulphide is susceptible to photolysis and cannot build up significantly.

Last, an alternative to the carbonate-silicate cycle stabilisation of Earth’s climate is the 'Gaia hypothesis', which proposes that life interacts with inorganic Earth as a self-stabilising system that maintains habitable conditions. The Archean biosphere was surely important for atmospheric composition, which intertwines Earth’s inhabitance and habitability. However, debate continues about if and how adaptive evolution of the biosphere can produce a more stable climate than an abiotic Earth.

Catling and Zahnle next consider proxy data for the Archean climate. Glacial rocks are one line of evidence. These include dropstones in sediments from melted icebergs, unsorted clasts mixed with silt and clay deposited at the base of glaciers that are preserved as diamictites, and striations made by rocks embedded in moving glaciers that scoured underlying surfaces. In the Archean, any indication of polar environments from glacial rocks would suggest a relatively cool world given that since about 34 million years ago, a climate with poles covered in ice has required average global temperatures below about 20°C.

Glacial rocks are reported from 3.5, 2.9, and 2.7 billion years ago. The oldest, from the Barberton of South Africa, includes clasts in finely laminated sediments, interpreted as dropstones, found below diamictites. Later, in the rounghly 2.9-billion-year-old Pongola in South Africa, which was 43° to 48° palaeolatitude, diamictites contain striated clasts, while associated silty laminates have dropstones with splash-up of substrata. Then, at 2.7 billion years ago, diamictites with dropstones occur in India, and Montana.

Another line of evidence about palaeotemperature are ratios of oxygen¹⁸/oxygen¹⁶ in marine cherts and carbonates. These ratios decline with increasing age, which, at face value, suggests Archean ocean temperatures of 50° to 85°C. Cherts and carbonates precipitated from seawater acquire less oxygen¹⁸ relative to the seawater as temperature increases because, at equilibrium, warmth allows stronger oxygen¹⁸ bonds to be broken and replaced with oxygen¹⁶. Not all Archean isotopic studies infer hot temperatures, however.Combined oxygen and hydrogen isotope ratios have suggested surface temperatures below 40°C, while oxygen isotopes in Archean phosphates have produced 26° to 35°C upper limits.

Most researchers have doubted the isotopic inferences of persistently high Archean surface temperatures because the climate was sometimes cold enough for the aforementioned glaciations, and chemical weathering of quartz from hot climates is lacking. Two alternative interpretations are as follows: (i) Archean surface temperatures were similar to today’s because the oxygen isotope composition of seawater, increased by about 15 parts per thousand oxygen¹⁸ since 3.5 billion years ago. (ii) Older cherts and carbonates have nonprimary oxygen¹⁸ proportions from high alteration temperatures or hydrothermal fluids.

The first explanation relies on changes in the relative rates of high and low temperature water-rock interactions. Hydrothermal seafloor interactions increase the proportion of oxygen¹⁸ in seawater, whereas low-temperature seafloor and continental weathering decrease the proportion of oxygen¹⁸ in seawater. Archean oceans with shallow seafloor hydrothermal circulation and an increase in Phanerozoic pelagic sediments that lowers seafloor weathering might cause secular increase in the proportion of oxygen¹⁸ in seawater. However, measurements find the ancient proportion of oxygen¹⁸ in seawater matches Vienna standard mean ocean water in 2-billion-year-old ophiolites, 3.8-billion-year-old serpentine minerals, Archean pillow basalts, and kerogens in Archean cherts. In addition, clumped isotopes show little change in Phanerozoic proportions of oxygen¹⁸ in seawater. Oxygen isotopes in iron oxides since 2 billion yeasr ago suggest a secular increase in the proportion of oxygen¹⁸ in seawater. 

The second possibility is that proportion of oxygen¹⁸ in cherts and carbonates is more altered with age. Microanalysis shows how chert replaced sedimentary carbonates and did not precipitate from seawater. In addition, stringent geochemical and petrographic selection criteria for palaeothermometry eliminate Archean cherts. The same issue applies for silicon isotopes in Archean cherts, which have been used to infer hot surface temperatures. Last, the triple oxygen isotope composition of Archean cherts cannot be reconciled with equilibrium precipitation from seawater and requires alteration.

Phylogenetic inferences of early thermophilic microbes have also been used to argue for hot Archean oceans However, other more convincing reconstructions and biochemical considerations suggest a mesophilic (living at below 50°C) common ancestor.

Temperature limits on mineral formation may also favor a temperate Palaeoarchean. Pseudomorphs of barite after gypsum are reported in 3.5-billion-year-old evaporites that are now partially silicified. In halitesaturated brines, gypsum forms at under 18°C; otherwise, anhydrite precipitates. It has been claimed from x-ray tomography that hydrothermal barite, not gypsum, was actually primary. However, three-dimensional universal stage petrography found interfacial angles diagnostic of gypsum at the original crystal faces, which are partially replaced by quartz and draped by quartz, and so cannot show up in x-ray density contrasts.

Additional factorsmay have affected the Archean climate. One external parameter was a faster rotation of early Earth. Today, the Moon recedes about 4 cm a year from Earth due to tides, so Earth is despinning over time to conserve angularmomentum in the Earth-Moon system.Going back in time, depositional cycles in the 2.45-billion-year-old banded iron Weeli Wolli Formation, Australia indicate a roughly 17 hour day.

In general, a faster Archean rotation rate reduces heat transport from equator to pole, particularly if the atmosphere had less mass. Thus,warm tropics and glaciated poleswould arise if the Archean
had a thin atmosphere.

External parameters that modulate the carbonate-silicate cycle would also affect the long-term climate. These include land fraction, carbon dioxide outgassing, and biological modification of weathering. In general, how carbon dioxide outgassing changedwith time dominates the uncertainty. 

The general view that emerges is an Archean atmosphere devoid of oxygen and enriched in carbon dioxide and methane greenhouse gases that countered the Faint Young Sun. Nitrogen was a bulk gas, but proxy data, although debated, suggest that nitrogen levels were similar to today or possibly a factor of a few lower, which, if correct, requires understanding of how the nitrogen cycle operated differently in an anoxic world.

An overview of post-Archean atmospheric evolution in the context of biological evolution and constraints on mean global temperature in the Archean in the context of the glacial record. (A) Uncertainties on gas concentrations are a factor of a few or more. Atmospheric nitrogen may have tracked oxygen levels due to an oxidative weathering and denitrification source of nitrogen, but nitrogen pressure changes are debated. Methane was oxidized as oxygen rose but could have been protected subsequently under an ozone layer, depending on post-Archean methane source fluxes. The secular decline of carbon dioxide is a feedback effect in the geological carbon cycle induced by decreasing solar luminosity. (B) Constraints on Archean mean global temperature. Vertical blue bars denote that glacial rocks exist, noting that the durations of glaciations in the early Proterozoic and earlier are poorly known. A proposed Mesoproterozoic glaciation is not plotted because its age is disputed and possibly Sturtian. Cainozoic glaciations only occur at a global mean temperature below about 20°C. Red arrows on the Archean glaciations are a more conservative 25°C upper limit, taking into account of the possible effects of different land configurations and lack of vegetation. Low carbon dioxide during the Phanerozoic (A) correlates with glaciations (B), such as the Carboniferous-Permian ones, 335 to 256 million years ago. Precambrian greenhouse gasesmust also have fluctuated, but the amount is unknown and so not reflected in (A). Catling & Zahnle (2020).

The Archean climate was probably mostly moderate. Although some argue for a hot Archean climate, glacial rocks at 2.7, 2.9, and 3.5 billion years ago and recent isotopic analyses call that idea into question; indeed, our current understanding of feedbacks in the geologic carbon cycle suggests surface temperatures within 0° to 40°C.

Reconstructing the history of oxygen is informed by proxies, but they are often indirect. Once cyanobacteria evolved, local or regional oxygen oases in lakes or shallow seawater are possible, and various redox-sensitive proxies suggest that these oases actually existed by 3.2 to 2.8 billion years ago under globally anoxic air. The newest relevant data are mass fractionations of xenon isotopes, which gradually increase through the Archean relative to initial solar values. These data are best explained by xenon escaping to space as an ion dragged out by hydrogen originating photolytically from an anoxic atmosphere enriched in methane until the Great Oxidation Event.

Despite improved knowledge of theArchean, a relative dearth of uncontested data that constrain basic environmental variables, such as Archean seawater pH, climatic temperature, barometric pressure, and changing levels of the greenhouse gases carbon dioxide and methane through time, means that more proxy development and measurements are critical. It is also implausible that the Archean, one-third of the history of Earth, had a constant climate. However, our knowledge of Archean climatic variability is meager.

The biosphere was surely a major influence on Archean atmospheric composition, as it is today. Consequently, resolving when key biological innovations evolved, such as nitrogen fixation, methanogenesis, anoxygenic photosynthesis, and oxygenic photosynthesis, and understanding their influence are essential for improving models of atmospheric evolution.

This understanding may help us interpret future exoplanet data because atmospheres on other rocky Earth-sized worlds would initially be anoxic. On Earth, oxygenic photosynthesis evolved only once, perhaps because of its biochemical complexity. Consequently, if life exists elsewhere, inhabited planets with Archean-like atmospheres may be the most common type. So, determining what the Archean atmosphere was made of, and the influence of life, could help us distinguish biogenic gases in exoplanet atmospheres and so find life elsewhere.

Finally, the connection between atmospheric evolution and debated trends in solid Earth evolution was outside the scope of Catling and Zahnle's review. However, quantifying temporal changes in land area, surface rock composition, weathering, the size of volcanic and metamorphic fluxes of gases, and the redox evolution of outgassing all need further development to understand their implications for Earth’s atmospheric evolution or to extrapolate to exoplanets.

See also...

https://sciencythoughts.blogspot.com/2020/04/identifying-worlds-oldest-impact.htmlhttps://sciencythoughts.blogspot.com/2020/04/microfossils-from-palaeoproterozoic.html
https://sciencythoughts.blogspot.com/2020/01/understanding-influence-of-large-bolide.htmlhttps://sciencythoughts.blogspot.com/2019/07/microfossils-from-archean-of-south.html
https://sciencythoughts.blogspot.com/2019/06/evaluating-possibility-that-iron-oxides.htmlhttps://sciencythoughts.blogspot.com/2015/06/microbial-sediments-from-early-to.html
Follow Sciency Thoughts on Facebook.