it has long been recognized that low stellar UV fluxes promote
N
2
O accumulation (e.g., Segura et al. 2003,2005; Grenfell
et al. 2014; Rugheimer & Kaltenegger 2018)and that higher
N
2
O surface fluxes would result in a more detectable
biosignature (Kaltenegger 2017; Schwieterman et al. 2018).
It has been hypothesized that biological fluxes of N
2
O during
the Proterozoic Eon (∼2500–540 million years ago, Ma)may
have been dramatically larger than at present, due to the limited
availability of copper catalysts in euxinic (anoxic and sulfur-
rich)oceans, which would have effectively short-circuited the
last metabolic step in the denitrification cycle (Buick 2007;
Figure 1). Higher atmospheric mixing ratios of N
2
O could
therefore also have contributed to greenhouse warming at the
time (Roberson et al. 2011), although lower O
2
concentrations
would have somewhat muted the impact of higher fluxes,
due to the reduced shielding of photolyzing UV radiation
(Roberson et al. 2011; Stanton et al. 2018). Under certain
conditions, N
2
O production on the earlier Earth could also have
been augmented by chemodenitrification—that is, the process
of N
2
O production via abiotic reduction of nitric oxide (NO)by
ferrous iron (Samarkin et al. 2010; Stanton et al. 2018).
Previous studies have not comprehensively examined the full
range of plausible N
2
Ofluxes over a range of pO
2
values and
stellar types, including the end-member scenario, in which the
nitrous oxide reductase enzyme that facilitates the last step in
the denitrification process simply does not evolve (Pauleta et al.
2013). In such a scenario, we will show that N
2
O can
accumulate to high concentrations—even for planets orbiting
FGK stars—with implications for the detectability of this
biosignature gas with current and upcoming observatories.
Here, we conduct a systematic photochemical and spectral
investigation of N
2
O as an exoplanet biosignature and place
upper limits on the N
2
O abundances and detectability from a
productive biosphere. In Section 2, we use the Earth system
(biogeochemical)model “cGENIE”to calculate denitrification
fluxes for an Earthlike marine biosphere as a function of
atmospheric oxygenation levels (pO
2
)and concentrations of
bioavailable phosphorous (as PO
4
3−
). We evaluate our results
against literature values for the Earth and generate realistic
bounds for plausible intermediate and maximum N
2
Ofluxes. In
Section 3, we calculate the photochemical stability and steady-
state mixing ratios of N
2
O given a large range of fluxes,
inclusive of those calculated in Section 2, with a variety of
oxygenation states and for stellar hosts that span the main
sequence (F4V to M8V). In Section 4, we generate emission
and transmission spectra for a subset of the scenarios
investigated in Section 3. Finally, we calculate the number of
transits of TRAPPIST-1e that will be required to detect N
2
O
with JWST for three N
2
Oflux scenarios, and find that detecting
N
2
O with NIRSpec is plausible for production fluxes near the
biospheric maxima. We discuss the implications and potential
false positives in Section 5. We conclude in Section 6.
2. Circumscribing Plausible Global N
2
O Fluxes with a
Biogeochemical Model
To map out how the total oceanic denitrification (and hence
the potential maximum N
2
O production)rate varies as a
function of pO
2
and phosphate availability (PO
4
3−
), we use the
cGENIE Earth system model of intermediate complexity.
cGENIE consists of a 3D ocean circulation model plus a 2D
energy balance and moisture model and a 2D sea-ice model.
The 2D grid is split into 36 ×36 equal-area cells, while we
adopt 16 depth layers in the ocean, following Cao et al. (2009).
cGENIE simulates a 3D marine biosphere, including phos-
phorous and nitrogen-limited primary production, and a set of
metabolisms, including aerobic respiration, anaerobic respira-
tion, methanogenesis, and aerobic methanotrophy (Ridgwell
et al. 2007; Olson et al. 2016). A simple 2D (not vertically
resolved)gridded calculation of basic atmospheric chemical
reactions is also included (Reinhard et al. 2020). cGENIE has
been leveraged to explore the coupled evolution of Earth’s
biosphere, atmosphere, and climate system over the entire
geologic timescale, including the Archean (e.g., Olson et al.
2013), Proterozoic (e.g., Olson et al. 2016; Reinhard et al.
2016,2020), and Phanerozoic (e.g., Kirtland Turner &
Ridgwell 2016). cGENIE has most recently been used to
explore the relationship between planetary obliquity, nutrient
cycling, and the consequent potential for atmospheric oxygena-
tion on exoplanets (Barnett & Olson 2022).
The biological N cycle in cGENIE includes diazotrophy (the
biological reduction of N
2
to NH
4
+
, which can then be
incorporated into biomass), nitrification (the oxidation of NH
4
+
to NO
3
−
), and denitrification (the biological reduction of NO
3
−
to N
2
). These processes are highlighted in Figure 1.
Diazotrophy occurs only when N is scarce relative to phosphate
(PO
4
3−
)and N:P <16 (the “Redfield Ratio”)within the photic
zone. Consequently, excess PO
4
3−
availability relative to N will
drive greater diazotrophy, until the global rates of N fixation
balance N loss. The rates of nitrification and denitrification are
both sensitive to atmospheric pO
2
, which directly influences
surface and benthic oxygen concentrations, but their relation-
ships to oxygen differ dramatically. Nitrification requires O
2
,
whereas denitrification occurs in the absence of O
2
. Denitri-
fication additionally requires reduced organic material and is a
multistep process, with several intermediate N species between
NO
3
−
and N
2
, such as N
2
O. The configuration of cGENIE
employed here neglects this complexity and, by default,
assumes the complete reduction of NO
3
−
to N
2
when organic
matter and NO
3
−
are in sufficient abundance and local dissolved
O
2
is low, following Naafs et al. (2019).
We estimate an upper bound on the possible N
2
Oflux arising
from incomplete denitrification for a given atmospheric pO
2
and ocean nutrient inventory by assuming that the entire
denitrification flux results in the evolution of N
2
O(rather than
going directly to N
2
). This could occur, for instance, if the
nitrous oxide reductase, the enzyme that facilitates the last step
of the denitrification process (Pauleta et al. 2013), has not
evolved or if dissolved copper, which is key to the functioning
of this catalyst, was severely rate-limiting in abundance
(Buick 2007). Nitrous oxide reductase can also be substantially
inhibited in community settings by other biological products,
including C
2
H
2
, CO, NO, N
3
−
, and CN
−
(Kristjansson &
Hollocher 1980; Koutný & Kučera 1999). Our goal here is not
to be overly prescriptive of the specific scenarios in which
denitrification is incomplete, but instead to examine plausible
maxima in N
2
O production by Earthlike biospheres.
Figure 2shows the total denitrification flux from Earth’s
marine biosphere as a function of pO
2
(relative to the present
atmospheric level, or PAL)and phosphate availability (PO
4
3−
,
relative to the present ocean level, or POL). We simulate pO
2
levels of 0%–100% and phosphate availability between one
and 2 times POL. The upper range of phosphate availability
represents a planet with higher nutrient availability from
enhanced continental weathering or a larger crustal abundance
3
The Astrophysical Journal, 937:109 (22pp), 2022 October 1 Schwieterman et al.