期刊名称:Proceedings of the National Academy of Sciences
印刷版ISSN:0027-8424
电子版ISSN:1091-6490
出版年度:2022
卷号:119
期号:35
DOI:10.1073/pnas.2202387119
语种:English
出版社:The National Academy of Sciences of the United States of America
摘要:Li et al. (
1) contend that the marine
10Be/
9Be system presents a conundrum, recording constant global weathering and denudation rates over the past 12 My (
2), whereas other marine sedimentary isotope records (e.g., Li, Sr, and Os) collectively seem to indicate an increase in late Cenozoic weathering and denudation. They resolve this conundrum with a model focused on demonstrating that oceanic
10Be/
9Be is insensitive to changes in weathering and denudation rates. We believe this conclusion is untenable for three reasons.
First, the model in ref.
1 contains several design errors. In particular, equations 8 and 11 are dimensionally incorrect. Further, α (erosion to denudation; figure S1 of ref.
1) obtained by a logarithmic regression is allowed to exceed 1, yielding negative weathering rates at high denudation. Also, equation 10 assumes the coastal reservoir to be proportional to coastal length, although coastal length is highly fractal in nature. Finally, both their methods to estimate Φ
del, the efficiency with which Be escapes the coastal trap, result in a larger mismatch between measured and modeled ocean
10Be/
9Be than those reported in our original study (
3), even with unchanged parameters and ocean basin
10Be/
9Be data, implying their mass balance may be invalid.
Second, to explore the dependence of
9Be delivery on sediment yield using dimensionally correct equations, we revert to the original equations 5 and A13 of ref.
3 and find that it is neither unique nor probable that ocean
10Be/
9Be is insensitive to denudation over three orders of magnitude (
Fig. 1). Ref.
1 assumes that riverine and coastal particle concentration is linearly dependent on sediment yield (equations 9 and 11). However, regressing the global river’s particle concentration against sediment yield results in a power law exponent of 0.76 (
4). This relationship results in
10Be/
9Be that is sensitive to denudation (
Fig. 1). Further, a particle concentration effect (
5) controlling the Be partition coefficient cannot be discounted—unlike in ref.
1—if colloids transport
9Be into the ocean, resulting in an even greater sensitivity (model e in
Fig. 1) to denudation.
Fig. 1.
Dependence of ocean
10Be/
9Be on sediment yield. The red curve (a) is from Li et al. (
1); their model incorrectly assumes that coastal particle concentration is linearly dependent on sediment yield and that a particle concentration effect (
5) is absent. The other curves are calculated using the equations and data from the original model (
3) and global ocean values without the South Pacific (table 1 of ref.
3). In curve b, Φ
del is constant at 0.063. In all other curves Φ
del depends on partition coefficient
K
d
coast
B
e
and on coastal particle concentration (equation A13 of ref.
3) that scales with sediment yield by a power law exponent of 0.76 derived by regressing global river data (
4). In curve c,
K
d
B
e
(
coast
)
is 10
4. In curve d,
K
d
B
e
(
coast
)
is 10
5. In curve (e),
K
d
B
e
(
coast
)
depends on the particle concentration effect (
5). The black point shows global ocean data as calibration to curves b and e without the South Pacific for dearth of data (
3), whereas the ref.
1 model includes the South Pacific. All models except the Li et al. model (
1) show sensitivity of
10Be/
9Be to sediment yield and hence to weathering flux.
Third, “boundary exchange”—the release of “reactive” terrigenous Be from particles into seawater during early marine diagenesis—presents another input trajectory of
9Be (
Fig. 2) (
3). Ref.
1 ignores this pathway, although it has been shown to operate for several radiogenic and stable isotope systems (
6), including beryllium (
7).
Fig. 2.
The main entry pathways of
9Be into the open ocean after ref.
3. (
A) The dissolved Be input pathway where Φ
del is an estimate of the efficiency by which dissolved Be escapes coastal trapping; this is addressed by the model in ref.
1. (
B) Input of “reactive” (sediment-bound)
9Be and redissolution during “boundary exchange” (
6). In this scenario the extent of the Be leakeage from sediment (blue arrows) sets the value of Φ
del. Although we still lack models that quantify Φ
del associated with pathway (
B), we know that in this case the input of terrigenous Be (low
10Be/
9Be) into the ocean (high
10Be/
9Be) does not depend on the extent of coastal trapping (
3).
Even if we do not yet fully understand the
10Be/
9Be system, we know that
87Sr/
86Sr and
178Os/
186Os do not allow one to disentangle changes in flux from changes in lithological source (
8) and that
7Li/
6Li does not track continental mass fluxes (
9). Only if one accepts the alternative model (
1) that Neogene weathering fluxes have increased does the resulting imbalance in the carbon cycle require poorly constrained compensatory fluxes like sulfide weathering, organic carbon oxidation, or reduced basalt weathering. In contrast, the constant weathering flux compatible with
10Be/
9Be (
2) ensures balanced carbon fluxes; in turn, global cooling resulted from an increasingly reactive land surface driven by a change in the global ratio between denudation and weathering (
10) as indicated by
7Li/
6Li (
9). Given the simplicity of this scenario one may ask why the other tracers cited do not present the conundrum, rather than
10Be/
9Be, which in this group of weathering tracers is the only mass flux proxy.