pith. sign in

arxiv: 2605.30242 · v2 · pith:5UZE676Rnew · submitted 2026-05-28 · 📊 stat.AP

Multi-source land-use emissions reveal rising airborne fraction

J. Eduardo Vera-Valdes This is my paper

Pith reviewed 2026-06-28 23:44 UTC · model grok-4.3

classification 📊 stat.AP
keywords airborne fractionland-use emissionsmixed-effects modelcarbon cycleCO2 emissionstrend estimationGlobal Carbon Budget
0
0 comments X

The pith

A mixed-effects model applied to all land-use emission series finds the airborne fraction of CO2 rose from 0.40 to 0.47 between 1959 and 2024.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper applies a mixed-effects statistical model to every available land-use and land-cover change emission series to estimate trends in the airborne fraction of anthropogenic CO2. It concludes that this fraction increased from about 0.40 to about 0.47 over 1959-2024. A sympathetic reader would care because the airborne fraction determines how much emitted CO2 accumulates in the atmosphere versus being taken up by land and ocean sinks. The random-effects structure explains why earlier studies reached weak or inconclusive results by not fully handling variation across emission estimates. The result remains after checks that exclude the final year and propagate denominator uncertainty.

Core claim

Using a mixed-effects model with random intercepts and slopes by LULC series on all measurement series from the Global Carbon Budget 2025, the airborne fraction increased over 1959-2024 from about 0.40 to about 0.47. This conclusion is robust to excluding the final year and to alternative specifications that explicitly propagate denominator uncertainty.

What carries the argument

Mixed-effects model with random intercepts and slopes by LULC series, which estimates a shared trend across emission measurement series while allowing each series its own intercept and slope.

If this is right

  • A rising share of anthropogenic CO2 emissions remains in the atmosphere.
  • Land and ocean sinks absorb a declining proportion of total emissions.
  • Carbon-budget assessments must reflect the higher accumulation rate.
  • Near-term mitigation requirements become more stringent to meet temperature targets.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • If the trend continues, atmospheric CO2 concentrations would rise faster for any given emission path.
  • The same modeling approach could be used to examine trends in other carbon-cycle fluxes measured with high inter-series uncertainty.
  • Updated emission inventories released in future Global Carbon Budget reports provide a direct test of whether the increase persists.

Load-bearing premise

The random intercepts and slopes by LULC series fully capture systematic differences and uncertainties across measurement series without biasing the estimated common trend.

What would settle it

Incorporating new land-use emission series through 2025 or later into the same mixed-effects model yields a flat or declining airborne fraction trend.

read the original abstract

The airborne fraction is the share of anthropogenic carbon dioxide emissions that remains in the atmosphere and is a key indicator of carbon-cycle response and remaining carbon budgets under continued emissions. Whether this share is rising remains debated because inference is sensitive to uncertainty in land-use and land-cover change (LULC) emissions. Here we use all available LULC measurement series from Global Carbon Budget 2025 and estimate airborne-fraction trends with a mixed-effects model with random intercepts and slopes by LULC series. We find that the airborne fraction increased over 1959-2024, from about 0.40 to about 0.47, and that this conclusion is robust to excluding the final year and to alternative specifications that explicitly propagate denominator uncertainty. These results clarify why earlier studies reported weak or inconclusive trend evidence and strengthen support for the view that an increasing share of emitted carbon dioxide is accumulating in the atmosphere rather than being taken up by land and ocean sinks, with implications for carbon-budget assessment and near-term mitigation requirements.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

1 major / 2 minor

Summary. The manuscript claims that the airborne fraction (AF) of anthropogenic CO2 emissions increased from approximately 0.40 to 0.47 over 1959-2024. This conclusion is obtained by applying a mixed-effects model with random intercepts and slopes by LULC series to all available land-use and land-cover change emission series from the Global Carbon Budget 2025; the trend is reported as robust to exclusion of the final year and to alternative specifications that propagate denominator uncertainty.

Significance. If the central trend estimate is unbiased, the multi-source approach strengthens evidence for a rising airborne fraction relative to prior single-series analyses, with direct implications for carbon-budget accounting and near-term mitigation targets. The use of publicly available GCB datasets and explicit robustness checks on denominator uncertainty are positive features.

major comments (1)
  1. [Methods (mixed-effects model)] Methods section (mixed-effects model): The specification AF_{t,s} = β0 + β1 year_t + u_{0,s} + u_{1,s} year_t + ε assumes conditional independence across the s LULC series given the random effects. Because every series is constructed from the identical atmospheric CO2 growth and fossil-fuel emissions in the numerator of the AF ratio, the series are deterministically linked; any unmodeled common measurement error or shock in the shared components propagates identically across groups. With the small number of series typical of GCB releases, the random-effects covariance cannot absorb this cross-series dependence, which risks attenuation or inflation of the fixed slope β1. The reported robustness checks address only denominator uncertainty and year exclusion, not this structural non-independence.
minor comments (2)
  1. [Abstract] Abstract and Methods: State the exact number of LULC series retained after any exclusions and report the estimated random-effects covariance matrix (or its condition number) to allow assessment of whether the random-slope structure is well-identified.
  2. [Results] Figure 1 or equivalent results figure: Add pointwise confidence bands that incorporate the full estimated covariance structure rather than only the fixed-effect standard error.

Simulated Author's Rebuttal

1 responses · 0 unresolved

We thank the referee for their constructive comments on the manuscript. We respond point-by-point to the major comment below.

read point-by-point responses

  1. Referee: Methods section (mixed-effects model): The specification AF_{t,s} = β0 + β1 year_t + u_{0,s} + u_{1,s} year_t + ε assumes conditional independence across the s LULC series given the random effects. Because every series is constructed from the identical atmospheric CO2 growth and fossil-fuel emissions in the numerator of the AF ratio, the series are deterministically linked; any unmodeled common measurement error or shock in the shared components propagates identically across groups. With the small number of series typical of GCB releases, the random-effects covariance cannot absorb this cross-series dependence, which risks attenuation or inflation of the fixed slope β1. The reported robustness checks address only denominator uncertainty and year exclusion, not this structural non-independence.

    Authors: We agree that the AF series share identical atmospheric growth and fossil-fuel emissions, inducing cross-series dependence through the common numerator. However, this shared component affects the AF level uniformly across all series at each time t and does not differentially influence the series-specific trends. The mixed-effects specification is designed to recover the population-average trend β1 while permitting series-specific random intercepts and slopes to capture heterogeneity arising from the LULC denominators. Because the common shock is absorbed into the fixed effects (or would appear as a common residual that does not correlate with the year trend in a manner that biases the slope), the point estimate of β1 remains consistent. The random-effects structure already accounts for between-series variation; with the small number of GCB series the model remains parsimonious and our existing robustness checks (denominator uncertainty propagation and year exclusion) test sensitivity to the series-specific components. We therefore maintain that the reported trend is not biased by this dependence structure. revision: no

Circularity Check

0 steps flagged

No circularity: trend estimated from external data via standard model

full rationale

The paper computes AF ratios from public Global Carbon Budget datasets (atmospheric CO2, fossil emissions, and multiple LULC series) and fits the mixed-effects model AF_t,s = β0 + β1*year_t + u0,s + u1,s*year_t + ε to recover the fixed slope β1. This produces the reported increase (0.40 to 0.47) as a model output rather than by definition or self-citation. No load-bearing step reduces to a fitted parameter renamed as prediction, no uniqueness theorem imported from the same authors, and no ansatz smuggled via citation. The derivation chain is self-contained against external benchmarks.

Axiom & Free-Parameter Ledger

0 free parameters · 0 axioms · 0 invented entities

Abstract-only review; no explicit free parameters, axioms, or invented entities are described. The work relies on standard mixed-model assumptions and external GCB data.

pith-pipeline@v0.9.1-grok · 5696 in / 1010 out tokens · 22186 ms · 2026-06-28T23:44:14.523995+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Forward citations

Cited by 1 Pith paper

Reviewed papers in the Pith corpus that reference this work. Sorted by Pith novelty score.

  1. Global Warming Has Been Accelerating Since At Least 1990

    stat.ME 2026-06 unverdicted novelty 5.0

    Statistical analysis of global temperature records finds evidence of supralinear acceleration since at least 1990 that strengthens with more recent data.

Reference graph

Works this paper leans on

51 extracted references · 4 canonical work pages · cited by 1 Pith paper · 1 internal anchor

  1. [1]

    Canadell, J. G.et al.Contributions to accelerating atmospheric CO2 growth from economic activity, carbon intensity, and efficiency of natural sinks.Proceedings of the National Academy of Sciences104, 18866–18870 (2007)

    2007

  2. [2]

    R.et al.Global and regional drivers of accelerating CO2 emissions

    Raupach, M. R.et al.Global and regional drivers of accelerating CO2 emissions. Proceedings of the National Academy of Sciences104, 10288–10293 (2007)

    2007

  3. [3]

    Is the airborne fraction of anthropogenic CO2 emissions increasing? Geophysical Research Letters36, L21710 (2009)

    Knorr, W. Is the airborne fraction of anthropogenic CO2 emissions increasing? Geophysical Research Letters36, L21710 (2009)

    2009

  4. [4]

    P., Alden, C

    Ballantyne, A. P., Alden, C. B., Miller, J. B., Tans, P. P. & White, J. W. C. Increase in observed net carbon dioxide uptake by land and oceans during the past 50 years.Nature488, 70–72 (2012)

    2012

  5. [5]

    Le Quéré, C.et al.Trends in the sources and sinks of carbon dioxide.Nature Geoscience2, 831–836 (2009)

    2009

  6. [6]

    & Koopman, S

    Bennedsen, M., Hillebrand, E. & Koopman, S. J. On the evidence of a trend in the CO2 airborne fraction.Nature616, E1–E3 (2023)

    2023

  7. [7]

    & Koopman, S

    Bennedsen, M., Hillebrand, E. & Koopman, S. J. A regression-based approach to the CO2 airborne fraction.Nature Communications15, 8507 (2024)

    2024

  8. [8]

    F., Salawitch, R

    Bennett, B. F., Salawitch, R. J., McBride, L. A., Hope, A. P. & Tribett, W. R. Quantification of the Airborne Fraction of Atmospheric CO2 Reveals Stability in Global Carbon Sinks Over the Past Six Decades.Journal of Geophysical Research: Biogeosciences129, e2023JG007760 (2024)

    2024

  9. [9]

    Vera-Valdés, J. E. & Grivas, C. Robust estimation of carbon dioxide airborne fraction under measurement errors.Environmental Research Communications7, 031009 (2025)

    2025

  10. [10]

    Friedlingstein, P.et al.Global carbon budget 2025.Earth System Science Data 18, 3211–3288 (2026). 9

    2025

  11. [11]

    Hansis, E., Davis, S. J. & Pongratz, J. Relevance of methodological choices for accounting of land use change carbon fluxes.Global Biogeochemical Cycles29, 1230–1246 (2015)

    2015

  12. [12]

    Gasser, T.et al.Historical CO2 emissions from land use and land cover change and their uncertainty.Biogeosciences17, 4075–4101 (2020)

    2020

  13. [13]

    Qin, Z.et al.Global spatially explicit carbon emissions from land-use change over the past six decades (1961-2020).One Earth7, 835–847 (2024)

    1961

  14. [14]

    Henderson, C. R. Estimation of variance and covariance components.Biometrics 9, 226–252 (1953)

    1953

  15. [15]

    M.et al.Generalized linear mixed models: A practical guide for ecology and evolution.Trends in Ecology & Evolution24, 127–135 (2009)

    Bolker, B. M.et al.Generalized linear mixed models: A practical guide for ecology and evolution.Trends in Ecology & Evolution24, 127–135 (2009)

    2009

  16. [16]

    A.Measurement Error Models

    Fuller, W. A.Measurement Error Models. (Wiley, New York, 1987)

    1987

  17. [17]

    J., Ruppert, D., Stefanski, L

    Carroll, R. J., Ruppert, D., Stefanski, L. A. & Crainiceanu, C. M.Measurement Error in Nonlinear Models: A Modern Perspective. (Chapman & Hall/CRC, Boca Raton, 2006)

    2006

  18. [18]

    & Thoning, K

    Lan, X., Tans, P. & Thoning, K. W. Trends in globally-averaged CO2 de- termined from NOAA global monitoring laboratory measurements. (2026) doi:10.15138/9N0H-ZH07

    work page doi:10.15138/9n0h-zh07 2026

  19. [19]

    & Tubiello, F

    Conchedda, G. & Tubiello, F. N. Drainage of organic soils and GHG emissions: Validation with country data.Earth System Science Data12, 3113–3137 (2020)

    2020

  20. [20]

    & Joos, F

    Müller, J. & Joos, F. Committed and projected future changes in global peat- lands – continued transient model simulations since the last glacial maximum. Biogeosciences18, 3657–3687 (2021)

    2021

  21. [21]

    Science Advances7, eabf1332 (2021)

    Qiu, C.et al.Large historical carbon emissions from cultivated northern peatlands. Science Advances7, eabf1332 (2021)

    2021

  22. [22]

    Haverd, V.et al.A new version of the CABLE land surface model (subversion revision r4601) incorporating land use and land cover change, woody vegetation demography, and a novel optimisation-based approach to plant coordination of photosynthesis.Geoscientific Model Development11, 2995–3026 (2018)

    2018

  23. [23]

    Melton, J. R.et al.CLASSIC v1.0: The open-source community successor to the canadian land surface scheme (CLASS) and the canadian terrestrial ecosystem model (CTEM) – part 1: Model framework and site-level performance.Geoscien- tific Model Development13, 2825–2850 (2020)

    2020

  24. [24]

    Lawrence, D. M.et al.The community land model version 5: Description of new features, benchmarking, and impact of forcing uncertainty.Journal of Advances in Modeling Earth Systems11, 4245–4287 (2019)

    2019

  25. [25]

    A.et al.Taking off the training wheels: The properties of a dynamic vegetation model without climate envelopes, CLM4.5(ED).Geoscientific Model Development8, 3593–3619 (2015)

    Fisher, R. A.et al.Taking off the training wheels: The properties of a dynamic vegetation model without climate envelopes, CLM4.5(ED).Geoscientific Model Development8, 3593–3619 (2015). 10

    2015

  26. [26]

    Tian, H.et al.North american terrestrial CO2 uptake largely offset by CH4 and N2O emissions: Toward a full accounting of the greenhouse gas budget.Climatic Change129, 423–426 (2015)

    2015

  27. [27]

    Geoscientific Model Development15, 1971–1994 (2022)

    Ma, L.et al.Global evaluation of the ecosystem demography model (ED v3.0). Geoscientific Model Development15, 1971–1994 (2022)

    1971

  28. [28]

    & Hoffman, F

    Yang, X., Thornton, P., Ricciuto, D., Wang, Y. & Hoffman, F. Global evaluation of terrestrial biogeochemistry in the energy exascale earth system model (E3SM) and the role of the phosphorus cycle in the historical terrestrial carbon balance. Biogeosciences20, 2813–2836 (2023)

    2023

  29. [29]

    Needham, J.et al.Vertical canopy gradients of respiration drive plant carbon budgets and leaf area index.New Phytologist246, 144–157 (2025)

    2025

  30. [30]

    Felzer, B. S. & Jiang, M. Effect of land use and land cover change in context of growth enhancements in the united states since 1700: Net source or sink?Journal of Geophysical Research: Biogeosciences123, 3439–3457 (2018)

    2018

  31. [31]

    Z.et al.The carbon budget of china: 1980–2021.Science Bulletin69, 114–124 (2024)

    Xia, J. Z.et al.The carbon budget of china: 1980–2021.Science Bulletin69, 114–124 (2024)

    1980

  32. [32]

    Yue, X.et al.Development and evaluation of the interactive model for air pollution and land ecosystems (iMAPLE) version 1.0.Geoscientific Model Development17, 4621–4642 (2024)

    2024

  33. [33]

    K., Koven, C

    Shu, S., Jain, A. K., Koven, C. D. & Mishra, U. Estimation of permafrost SOC stock and turnover time using a land surface model with vertical heterogeneity of permafrost soils.Global Biogeochemical Cycles34, e2020GB006585 (2020)

    2020

  34. [34]

    H.et al.JSBACH 3 - the land component of the MPI earth system model: Documentation of version 3.2

    Reick, C. H.et al.JSBACH 3 - the land component of the MPI earth system model: Documentation of version 3.2. (2021) doi:10.17617/2.3279802

    work page doi:10.17617/2.3279802 2021

  35. [35]

    C., Hodson, E

    Poulter, B., Frank, D. C., Hodson, E. L. & Zimmermann, N. E. Impacts of land cover and climate data selection on understanding terrestrial carbon dynamics and the CO2 airborne fraction.Biogeosciences8, 2027–2036 (2011)

    2027

  36. [36]

    Smith, B.et al.Implications of incorporating n cycling and n limitations on pri- mary production in an individual-based dynamic vegetation model.Biogeosciences 11, 2027–2054 (2014)

    2027

  37. [37]

    Schaphoff, S.et al.LPJmL4 – a dynamic global vegetation model with managed land – part 1: Model description.Geoscientific Model Development11, 1343–1375 (2018)

    2018

  38. [38]

    & Joos, F

    Lienert, S. & Joos, F. A bayesian ensemble data assimilation to constrain model parameters and land-use carbon emissions.Biogeosciences15, 2909–2930 (2018)

    2018

  39. [39]

    Vuichard, N.et al.Accounting for carbon and nitrogen interactions in the global terrestrial ecosystem model ORCHIDEE (trunk version, rev 4999): Multi-scale evaluation of gross primary production.Geoscientific Model Development12, 4751–4779 (2019). 11

    2019

  40. [40]

    Walker, A. P.et al.The impact of alternative trait-scaling hypotheses for the maximum photosynthetic carboxylation rate (v-cmax) on global gross primary production.New Phytologist215, 1370–1386 (2017)

    2017

  41. [41]

    & Yamagata, Y

    Kato, E., Kinoshita, T., Ito, A., Kawamiya, M. & Yamagata, Y. Evaluation of spatially explicit emission scenario of land-use change and biomass burning using a process-based biogeochemical model.Journal of Land Use Science8, 104–122 (2013)

    2013

  42. [42]

    Disequilibrium of terrestrial ecosystem CO2 budget caused by disturbance- induced emissions and non-CO2 carbon export flows: A global model assessment

    Ito, A. Disequilibrium of terrestrial ecosystem CO2 budget caused by disturbance- induced emissions and non-CO2 carbon export flows: A global model assessment. Earth System Dynamics10, 685–709 (2019)

    2019

  43. [43]

    Gloor, M.et al.What can be learned about carbon cycle climate feedbacks from the CO2 airborne fraction?Atmospheric Chemistry and Physics10, 7739–7752 (2010)

    2010

  44. [44]

    Oehlert, G. W. A note on the delta method.The American Statistician46, 27–29 (1992)

    1992

  45. [45]

    Aitken, A. C. On least squares and linear combination of observations.Proceedings of the Royal Society of Edinburgh55, 42–48 (1935)

    1935

  46. [46]

    Huang, B.et al.Extended reconstructed sea surface temperature, version 5 (ERSSTv5): Upgrades, validations, and intercomparisons.Journal of Climate https://doi.org/10.1175/JCLI-D-16-0836.1 (2017) doi:10.1175/JCLI-D-16-0836.1

    work page doi:10.1175/jcli-d-16-0836.1 2017

  47. [47]

    M., Meehl, G

    Ammann, C. M., Meehl, G. A., Washington, W. M. & Zender, C. S. A monthly and latitudinally varying volcanic forcing dataset in simulations of 20th century climate.Geophysical Research Letters30, (2003). Methods Our primary estimator is a linear mixed-effects trend model with random intercepts and random slopes for each LULC measurement series14,15. As a r...

    2003

  48. [48]

    Linearizegaround the mean vector(µG,t,µC,t) =E[( ˆGt, ˆCt)]: g(Xt)≈g(µG,t,µC,t) +∇g(µG,t,µC,t)⊤(Xt−E[Xt])

  49. [49]

    Compute the gradient: ∇g(G,C) = [ ∂g/∂G ∂g/∂C ] = [ 1/C −G/C2 ]

  50. [50]

    Write the covariance matrix of(ˆGt, ˆCt): Σt = [ σ2 G,t σGC,t σGC,t σ2 C,t ] . 16

  51. [51]

    Apply the delta-method variance formula Var(g(Xt))≈∇g(µG,t,µC,t)⊤Σt∇g(µG,t,µC,t), and with plug-in evaluation at(Gt,Ct), this becomes Var(AFt)≈ ( 1 Ct )2 σ2 G,t + ( Gt C 2 t )2 σ2 C,t−2Gt C 3 t σGC,t. In practice, we use plug-in estimates (replace unknown moments by their empirical counterparts): σ2 G,t≈ˆVar( ˆGt), σ 2 C,t≈ˆVar( ˆCt), σ GC,t≈0, where we s...

    work page internal anchor Pith review Pith/arXiv arXiv 1959