ORIGINAL RESEARCH article

Front. Astron. Space Sci., 01 November 2019

Sec. Earth and Planetary Materials

Volume 6 - 2019 | https://doi.org/10.3389/fspas.2019.00062

On the Role of the Urey Reaction in Extracting Carbon From the Earth's Atmosphere and Adding It to the Continental Crust

  • Department of Earth and Planetary Sciences, University of California, Davis, Davis, CA, United States

Abstract

Urey (,, ) introduced the Urey reaction to explain the origin of carbonates in the continental crust. The reaction extracts CO2 from the atmosphere in acid rain that reacts with calcium silicates, the products are transported to the oceans where organic and inorganic processes results in the deposition of calcium carbonates. The origin of the CO2 can be (1) the early atmosphere and (2) the mantle. The reaction removes almost all the atmospheric CO2 that survived the formation of the Earth. If surface volcanism introduces more CO2 than is lost by subduction, the CO2 entering the atmosphere will be lost to the continental crust by the Urey reaction. Studies of carbon fluxes between surface reservoirs have demonstrated the importance of the Urey reaction in controlling the concentration of CO2 in the atmosphere. At the present time the only surface reservoir with a large mass of carbon is the continental crust. In order to quantify the rate at which the Urey reaction removes CO2 from the atmosphere we utilize data from the Paleocene-Eocene thermal maximum (PETM). This was a period of elevated global temperatures (ΔT = 4 − 5°C) at 56 Ma attributed to a pulse of volcanism and associated CO2 in the north Atlantic. We utilize the decay time of this thermal anomaly to quantify the rate at which CO2 is extracted from the atmosphere by the Urey reaction.

1. Introduction

We first give a brief introduction to the distribution of carbon in the Earth. In doing this we will utilize the concept of carbon reservoirs. In some cases the definition of a carbon reservoir is well-defined. The Earth's atmosphere is an example. Carbon in the atmosphere is primarily in CO2 and its concentration is relatively uniform. In other cases the definition of a carbon reservoir is arbitrary. The Earth's mantle is an example. The relatively uniform composition of mid-ocean ridge basalts (including CO2) is taken as evidence of a near-uniform upper mantle reservoir. However, the variability of ocean-island and hot-spot basalts is evidence for a complex heterogeneous lower-mantle.

Houghton () considered the present distribution of carbon in near-surface reservoirs including the atmosphere, oceans, ocean sediments, vegetation and soil, and fossil fuels. However, the near-surface reservoirs contain only a small fraction of the carbon in the Earth. The major reservoirs for carbon at the present are the core, the mantle, and the continental crust. It is generally accepted that a large fraction of the carbon resides in the core reservoir. Both during early core formation and accretion, and during core formation after the moon-forming giant impact, carbon in the solid Earth was largely partitioned into the core. Dasgupta () estimates that the mass of carbon in the core is Gt (1 Gt = 1012 kg) but there is certainly a large error bar. We suggest that a reasonable estimate for the present mass of carbon in the core is Gt.

A number of authors have studied the global carbon flux at mid-ocean ridges. These values can be used to estimate the composition of carbon in the mantle source region. Michael and Graham () give a mass fraction and Rosenthal et al. () give a value . We suggest that a reasonable estimate for the composition of carbon in this source region is . Assuming the entire mantle to have this concentration, the mass of carbon in the mantle is Gt.

The third large reservoir for carbon at the present time is the continental crust. About 1% of the carbon in the Earth is in the continental crust; a large fraction of this carbon is in carbonates (limestones, marble, etc.) Wedepohl () has given a comprehensive review of the composition of the continental crust with an emphasis on carbon. He gives a total carbon mass in the continental crust Gt. Of this total he attributes about 80% to sediments, almost entirely in carbonates. Hayes and Waldbauer () have reviewed the literature on carbon in the crust and have given a range of masses and gave a preferred total carbon mass Gt. In this paper we will accept the value given by Wedepohl () as a representative value for the present mass of carbon in the continental crust Gt, but suggest that a reasonable range of acceptable values is – 2 × 108 Gt.

A list of the estimated present masses of carbon in the reservoirs we consider is given in Table 1. The estimated total mass of carbon in all of the reservoirs is cM = 4.68 × 109 Gt. With the mass of the Earth, 5.97 × 1024 kg, the concentration of carbon in the bulk Earth by mass is about 0.1%. These values are consistent with previous estimates (Dasgupta, ; DePaolo, ). We do not give error bars because the uncertainties of some values are difficult to specify.

Table 1

ReservoirMass of carbon GtReferences
Core4 × 109Dasgupta,
Mantle2 × 108This paper
Continental crust4.2 × 107Wedepohl,
Oceans3.8 × 104Houghton,
Atmosphere8.5 × 102NOAA,
Total4.24 × 109

Masses of carbon in the Earth's carbon reservoirs considered in this paper (1 Gt = 1012 kg).

A large fraction of the carbon in the surface reservoirs is in the continental crust, ~99.99%. Much of this carbon is in the carbonates (limestones, marbles). Some 20–30% of the sediments in the continental crust are carbonates, independent of age (Mackenzie and Morse, ). Urey (,, ) proposed that a reaction between atmospheric CO2 and calcium silicate was a major source of carbon in the continental crust. This reaction is now known as the Urey reaction. The principal purpose of this paper is to evaluate the role of the Urey reaction in terms of our present knowledge of the distribution of carbon in the Earth.

Two sources of the carbon in the continental crust are available. The first of these is the atmosphere. Today the atmosphere has very little carbon in it, however the atmosphere of Venus has a large mass of carbon as CO2. It is quite possible that the early Earth also had a very high concentration of CO2, the Urey reaction could have removed a significant fraction of this carbon to the continental crust. The second source of the carbon in the continental crust is the mantle. If the flux of carbon out of the mantle in surface volcanism exceeds the loss of carbon in subduction, the mantle could have lost carbon to the continental crust. The Urey reaction could have played an important role in this process.

2. The Urey Reaction

Urey (,, ) introduced the Urey reaction to explain the origin of calcium carbonates (limestones, marbles) in the continental crust. The basic reaction combines atmospheric CO2 with a calcium silicate to generate a calcium carbonate plus silica. In order to include the role of acid rain in the process we write the Urey reaction in the form

The aqueous products are removed by calcium carbonate precipitation. Carbon dioxide (carbonic acid) in acid rain dissolves calcium silicate (wollastonite) sediments to give calcium, bicarbonate, and silica. The resulting calcium and bicarbonate ions flow in rivers to the oceans where either organic or inorganic precipitation produces the calcium carbonate. Blättler and Higgins () have considered Urey's carbonate-silicate cycle in some detail utilizing measured calcium isotope ratios.

The role of the Urey reaction in maintaining a low concentration of CO2 in the Earth's atmosphere was explicitly stated by Urey (). “As carbon dioxide was formed it reacted with silicates to form limestone. Of course the silicates may have been a variety of minerals but the pressure of CO2 was always kept at a low level by this reaction or similar reactions just as it is now.”

Although the mass of carbon in the atmosphere today is small (~850 Gt), the mass may have been much higher in the past. One of the major differences between Venus and the Earth is atmospheric composition. The atmospheric pressure on Venus is about a factor of 100 greater than the atmospheric pressure on Earth and is 96% carbon dioxide. The mass of carbon in the Venus atmosphere is Gt. Scaling the atmospheric carbon masses to the overall masses of Venus and the Earth gives an estimate of the mass of carbon in the early atmosphere of the Earth. The estimated value is Gt (Kasting and Ackerman, ).

Another major difference between the Earth and Venus is the giant impact that is hypothesized to have been responsible for the origin of the moon. During the final stages of the Earth's accretion, a collision between the Earth and an astronomical body about the size of Mars occurred 20–100 Myrs after the origin of the solar system. A detailed review of the influence of this collision on the first several hundred million years of Earth evolution has been given by Zahnle et al. (). The intensity of the collision melted and vaporized significant fractions of the Earth and the incoming body. Some fraction of the incoming body (and some fraction of the Earth) subsequently accreted and cooled to form Earth's moon. It is estimated that the hot silicate atmosphere cooled and condensed over about a thousand years. A global magma ocean cooled and solidified much more slowly, possibly over tens of millions of years. Sleep and Zahnle () and Sleep et al. (, ) have proposed that a substantial fraction of the CO2 in the atmosphere was absorbed by the solid surface crust of the global magma ocean. Episodic foundering subduction of this dense crust transported the carbon into the Earth's mantle, removing CO2 from the atmosphere.

An alternative to the mechanism proposed by Sleep and Zahnle () is that carbon was directly transferred to the mantle during accretion (Li et al., ) or that carbon was transferred by tectonic processes during the Hadean. Shaw () has provided a comprehensive review of CO2 evolution in the Hadean to the early Proterzoic.

An important question concerning carbon is the origin of the carbon in the continental crust. One hypothesis is that they were extracted from the atmosphere by the Urey reaction relatively early in the Earth's history. This hypothesis was put forward by Kramers () and by Lowe and Tice (). It should be noted that the potential mass of carbon in the early atmosphere, 1.57 × 108 Gt, is substantially larger than the estimated present mass of carbon in the continental crust Gt. Thus, a substantial fraction of the carbon in the continental crust could have been transferred from the atmosphere.

In order for the Urey reaction to extract CO2 from the atmosphere the early Earth must have had continental crust in order to generate surface deposits of calcium silicates. In addition the Earth must of had oceans in order for the acid rain to catalyze the Urey reaction between atmospheric CO2 and the surface deposits of calcium silicates. Little data is available for timing the initiation of the extraction of CO2 from the atmosphere. We will assume that the process began at a time t0 after the solidification of the magma ocean at about 4.5 Ga. We further assume that the Urey reaction extracted carbon from the atmosphere at a constant rate until the concentration of CO2 in the atmosphere was reduced to a very low level. The rate of loss of carbon from the atmosphere is given by

We require at t = t0 and obtain

This transport is completed at t = τ when where

Taking Gt we give the dependence of the carbon flux on the flux time τ − t0 in Figure 1. Combining Equations (3) and (4) we find the time dependence of the mass of carbon in the atmosphere

for the period t0 < t < t0 + τ. For t0 + τ < t < t0 + tp, we have added to the continental crust.

Figure 1

Again taking Gt we give the mass of carbon in the continental crust as a function of time t in Figure 2 for t0 = 1 Gyr and τ = 1 and 3 Gyr. Carbon is lost from the atmosphere to the continental crust after plate tectonics provides the calcium silicate required for the Urey reaction to proceed.

Figure 2

In discussing the role of the Urey reaction we focus our attention on its role in the generation of carbonate sediments in the continental crust. Above we discussed how it could be directly responsible for the extraction of primordial CO2 from the atmosphere. We now turn our attention to its indirect role in the addition of carbon from the mantle into the continental crust. Carbon is extracted from the mantle by volcanism and is returned to the mantle by subduction. If the rate of addition exceeds the rate of loss then the difference may be responsible for the origin of a large fraction of the carbon in the continental crust. Order of magnitude estimates of this process given by Hayes and Waldbauer () demonstrated the plausibility of this mechanism.

A comprehensive model for the variability of atmosphere CO2 over Phanerozoic times has been given by Berner and Kothavala (). This model, GEOCARB III, is complex and involves both organic and inorganic processes. Transport of carbon between the atmosphere, oceans, and continental crust is quantified on the million year time scale. The balance is dominated by the exchange of carbon between carbonates in the continental crust and carbon in the surficial reservoirs (oceans and atmosphere) and organic carbon (Berner and Caldeira, ). This exchange occurs through the Urey reaction given in Equation (1). When erosion is high, the Urey reaction extracts CO2 from the atmosphere adding carbonates to the continental crust. High erosion rates are associated with low sea level and large continental areas. When erosion is low, the Urey reaction operates in the opposite direction (from right to left in Equation 1) with carbonates decomposing to give CO2. An example of this metamorphic process is the subduction of carbonate sediments and the generation and return of CO2 to the atmosphere in subduction zone volcanics (Frezzotti et al., ).

We will next provide estimates for the rate at which surface volcanism transfers carbon from the mantle to the surface reservoirs. Volumetrically, surface volcanism is dominated by mid-ocean ridge volcanism. A number of authors have studied the global carbon flux from the mantle at mid-ocean ridges. Some studies involve direct measurements of CO2 concentrations and others utilize correlations with other elements, such as barium. Hayes and Waldbauer () summarized the early literature and gave a mid-ocean ridge global carbon flux Mtyr−1. Dasgupta and Hirschmann () summarized the available data and gave a global carbon flux Mtyr−1. Michael and Graham () obtained an estimate of the carbon flux from mid-ocean ridges using measurements of CO2/Ba ratios with published compilations of trace elements. They suggest that the global flux of carbon out of the mantle at mid-ocean ridges is Mtyr−1. Rosenthal et al. () used CO2/Ba and CO2/Nb systematics and conclude that the global flux of carbon is Mtyr−1. Chavrit et al. () carried out CO2 vesicularity measurements on some 400 mid-ocean ridge glasses. These results were interpreted to give a global carbon flux Mtyr−1. Since not all CO2 would be expected to be in cavities, the actual carbon flux would be expected to be somewhat higher. Based on these previous studies we suggest that a reasonable choice for the carbon flux from the mantle at mid-ocean ridges is Mtyr−1.

Flux estimates of carbon in ocean island basalts and hot spot volcanics have been reviewed by Hayes and Waldbauer (). Volumetrically this volcanism is about 10% of the mid-ocean ridges volcanism but some of this volcanism may be tapping a primitive lower mantle reservoir with higher concentrations of carbon. Considering the uncertainties we will neglect this source of carbon to the atmosphere in our modeling.

Just as carbon is lost from the mantle at mid-ocean ridges, carbon is returned to the mantle at subduction zones. A detailed study of carbon fluxes at subduction zones has been carried out by Kelemen and Manning (). These authors suggest that the downward flux of carbon at global subduction zones is 53 ± 13 Mtyr−1. However, a substantial fraction of this carbon never makes it to the mantle. Carbon is returned to the oceans and atmosphere in subduction zone volcanism. Kelemen and Manning () estimate the loss of carbon due to subduction zone volcanics is in the range 30 ± 12 Mtyr−1. Taking the difference of these values, the present subduction flux of carbon to the mantle is Mtyr−1. Dasgupta and Hirschmann () have also provided a comprehensive study of carbon subduction and conclude that Mtyr−1 reach the mantle. An earlier review of carbon subduction was given by Hayes and Waldbauer (). These authors concluded that the input flux of carbon is 36 ± 12 Mtyr−1 and that the flux to the mantle is Mtyr−1.

We obtain an estimate for the net flux of carbon from the mantle required to produce a large mass of carbon in the continental crust assuming that this flux has been a constant during the period tpt0 to tp (present time). We use the relation

Taking the carbon mass transported from the mantle to the continental crust Gt we give the dependence of the carbon flux on the flux time tpt0 in Figure 3. The results are essentially identical to those given in Figure 1 for the transfer from the atmosphere. This is not surprising since Equation (6) is basically the same as Equation (3).

Figure 3

The time dependence of the mass transferred from the mantle to the continental crust is given by

Again taking Gt we give the dependence of the mass of carbon added to the continental crust as a function of time t in Figure 4 for t0 = 1.4 Gyr.

Figure 4

When excess CO2 in the atmosphere has been depleted by the Urey reaction an approximate steady state balance is established between the volcanic input of CO2 and extraction by the Urey reaction. It has been argued by many authors that this balance is controlled by the Urey reaction. We approximate this balance by the relation

where is the rate of volcanic input of carbon into the atmosphere. We assume that the extraction rate is proportional to the mass of carbon in the atmosphere . The rate at which the Urey reaction rate processes carbon τacc takes account of the rate at which acid rain can interact with calcium silicate sediments, clearly τacc can be a function of time.

3. PETM

As a further test of our approach to the use of the Urey reaction to explain variations in atmospheric CO2 concentrations we will consider the Paleocene-Eocene thermal maximum (PETM). This was a period of elevated global temperatures that occurred at about 56 Ma. The elevated temperatures are attributed to increased concentrations of CO2 and possibly methane in the atmosphere. The maximum increase in global surface temperatures is in the range 4−5°C and the period of excess temperatures lasted some 106 yrs. The occurrence of this period of high atmospheric CO2 and temperatures was first documented by Kennett and Stott () from an observation of large carbon and oxygen isotope ratios in an ocean drilling core.

A comprehensive review of the subsequent studies of this event was given by McInerney and Wing (). By the time this review paper was published some 400 research papers on a wide variety of aspects of PETM were included. The onset of the period of high atmosphere CO2 and temperature was dated at 56.3 Ma, the onset lasted <10 ka, and the subsequent decay lasted some 120–220 ka. Observational data used to constrain the structure of the thermal pulse involves a variety of isotope measurements in drilling cores. One example is the temporal dependence of the observed δ13C anomaly in a drilling core (Cui et al., ). A second example is the use of boron isotope data as a proxy for seawater (Gutjahr et al., ). Jones et al. () have carried out a detailed study of the warming associated with the PETM. They conclude that the global mean surface anomaly is in the range 4–5°C.

Several sources of the carbon associated with the PETM have been proposed. These include methane clathrates and melting of permafrost. Storey et al. () have made a very strong case for associating the PETM with the flood basalt volcanism resulting from the opening of the north Atlantic. This volcanism is documented in east Greenland and the Faroe Islands. The earliest of these flood basalts have an age of 56.1±0.5 Ma. We support the hypothesis that a substantial fraction of the carbon responsible for the PETM was CO2 entering the atmosphere from these flood basalts.

Studies of isotopic anomalies associated with the PETM systematically show a rapid increase at about 55.8 Ma with a subsequent decay of the anomaly over about 105 yrs. There is considerable scatter between drill sites. As an example of an isotopic anomaly we give data from the Deep Sea Drilling Project (DSDP) site 401 in the northeast Atlantic Ocean (Gutjahr et al., ). A strong negative carbon isotope anomaly was found in foraminifera Morozovella subbotinae. The observed δ13C anomaly is given as a function of time tPETM relative to the onset of the anomaly in Figure 5.

Figure 5

). Also shown is the least squares best fit of our relaxation model from Equation (12) with relaxation time τacc = 105 kyr and .

Gutjahr et al. () have carried out an extensive modeling study using the available isotope data to provide estimates for the carbon generated during the PETM. We will utilize these values and take the total carbon added to the Earth's surface to be 12,000 Gt. Before and after the PETM we take the background carbon mass in the atmosphere to be Gt and the peak mass of carbon in the atmosphere during the PETM to be Gt.

We now carry out an approximate analysis of the decay of the PETM due to the loss of CO2 from the atmosphere by the Urey reaction. We extend the balance given in Equation (8) to include the transient removal of carbon from the atmosphere and write

From Equation (8) the background mass of carbon in the atmosphere is given by

We prescribe an initial mass of carbon in the atmosphere at , with and solve Equation (10) taking τacc to be constant with the result

The excess mass of carbon in the atmosphere decays exponentially as the Urey reaction extracts carbon from the atmosphere.

In order to estimate the value of the characteristic decay time τacc we assume that the negative carbon isotope anomaly is proportional to the mass of carbon in the atmosphere and write Equation (11) as

where is the background isotope ratio before and after PETM and is the minimum isotope ratio prescribed at t = 0. Taking the dependence of δ13C on tPETM from Equation (12) is compared with the isotope anomaly values given in Figure 5. The least squares fit of the exponential decay to the data points in given taking τacc = 105 kyr and .

We next use this value of τacc to obtain the dependence of atmosphere carbon mass on time during PETM based on the model dependence given in Equation (11). Taking the values Gt and Gt as given by Gutjahr et al. () with τacc = 105 kyr the model results are given in Figure 6.

Figure 6

We now return to Equation (10). This result relates the background atmospheric carbon mass to the background rate of volcanic input of CO2 carbon into the atmosphere and the Urey decay time τacc. During the PETM we have taken the background carbon mass Gt. Taking τacc = 105 kyr we find from Equation (12) that Mtyr−1. This is an independent determination of the volcanic flux of carbon into the atmosphere. In Figure 4 we gave volcanic the flux required to generate a specific mass of carbon in the continental crust

4. Discussion

Urey () proposed the Urey reaction, Equation (1) to explain the origin of carbonates in the continental crust. He argued that the reaction would essentially remove all CO2 from the atmosphere. It is now accepted that in analogy to Venus, there may have been a large mass of carbon in the Earth's early atmosphere, as much as 108 Gt. However, only a fraction of this may have survived the moon forming impact. We give a very simplified model for the extraction of carbon from the atmosphere to the continental crust, taking the extraction rate to be constant. There are basically no constraints on the variation of this rate with time. If a significant fraction of the carbon in the continental crust was extracted from the atmosphere, it is likely that it occurred early in Earth's history as illustrated in Figure 1.

The Urey reaction also controls the equilibrium mass of carbon in the atmosphere after the removal of any large initial concentration. The input of carbon to the atmosphere is from volcanism and we show that the equilibrium mass of carbon in the atmosphere is proportional to the rate of volcanic injection divided by a characteristic Urey time τu. We quantify the value of τu by studying the observed relaxation of the Paleocene-Eocene thermal maximum, which occurred at 56 Ma, and the relaxation time is about yrs.

We also give a simplified model for the extraction of carbon from the atmosphere to the continental crust. If the volcanic loss of carbon from the mantle by volcanism exceeds the return of carbon by subduction, the difference is added to the continental crust. If the volcanic carbon enters the oceans organic precipitation creates carbonates. If the volcanic carbon enters the atmosphere it enters the continental crust by the Urey reaction. Current estimates of carbon fluxes from and to the mantle are sufficient to have produced all the carbon in the continental crust. At the present time it is not possible to quantify the relative importance of carbon addition to the continental crust from the early atmosphere and the mantle.

We have addressed two major questions concerning carbon in the atmosphere in this paper. The first is the origin of the carbon in the continental crust. We conclude that it is possible the carbon could have been extracted either from the early atmosphere or from the mantle over a longer period of time. Studies of the concentration of carbon in the atmosphere and continental crust over geologic time are required and should receive a high priority.

The second question we have addressed is the relaxation of injections of carbon into the atmosphere back to equilibrium values. We quantify this by studying the Paleocene-Eocene thermal maximum (PETM). This has obvious implications for the recovery from the process of anthropogenic injection of carbon into the atmosphere. We find the relaxation time to be about 50,000 yrs.

Statements

Data availability statement

All datasets generated for this study are included in the manuscript/supplementary files.

Author contributions

All authors listed have made a substantial, direct and intellectual contribution to the work, and approved it for publication.

Acknowledgments

The authors would like to acknowledge financial support from the Alfred P. Sloan Foundation under Grant G-2818-10087 through the Deep Carbon Observatory.

Conflict of interest

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

References

Summary

Keywords

Urey reaction, carbon, continental crust, atmosphere, early Earth

Citation

Kellogg LH, Turcotte DL and Lokavarapu H (2019) On the Role of the Urey Reaction in Extracting Carbon From the Earth's Atmosphere and Adding It to the Continental Crust. Front. Astron. Space Sci. 6:62. doi: 10.3389/fspas.2019.00062

Received

30 April 2019

Accepted

13 September 2019

Published

01 November 2019

Volume

6 - 2019

Edited by

Sami Mikhail, University of St Andrews, United Kingdom

Reviewed by

Fang Huang, Rensselaer Polytechnic Institute, United States; Nikolai Bagdassarov, Goethe University Frankfurt, Germany

Updates

Copyright

*Correspondence: Donald L. Turcotte

This article was submitted to Earth and Planetary Materials, a section of the journal Frontiers in Astronomy and Space Sciences

Disclaimer

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

Outline

Figures

Cite article

Copy to clipboard


Export citation file


Share article

Article metrics