Reconstruction of annual solar irradiance over the last three millennia
Pith reviewed 2026-05-22 03:59 UTC · model grok-4.3
The pith
A physics-based model produces the first annual total solar irradiance record over the last three millennia.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
The authors produce a continuous physics-based total solar irradiance record at annual resolution covering the last three millennia by extending the SATIRE-T model with annually resolved sunspot numbers from cosmogenic isotope records, resulting in a maximum difference of 1.04 with uncertainties of -0.2 to +0.14 W m^{-2} between the maximum and minimum of the 50-year running mean values over the full interval.
What carries the argument
The SATIRE-T model, which calculates irradiance changes from the evolution of solar surface magnetic features, applied to sunspot numbers reconstructed from cosmogenic isotopes.
If this is right
- The record supplies a uniform annual TSI series that bridges the satellite era and earlier times for direct use in climate studies.
- It quantifies the amplitude of long-term solar variability as no more than about one watt per square meter in smoothed values.
- It demonstrates that a physics-based approach can be carried backward in time once suitable sunspot input series exist.
Where Pith is reading between the lines
- Climate simulations could now ingest this TSI series to test how much of observed past temperature changes can be attributed to solar forcing alone.
- If later work confirms the small amplitude, it would tighten upper limits on solar contributions to centennial-scale climate shifts.
- The method opens a path for similar annual reconstructions of other solar parameters such as ultraviolet flux once suitable input data become available.
Load-bearing premise
The model relating sunspot numbers to irradiance changes, which was tuned on recent observations, continues to hold when applied to sunspot numbers from thousands of years ago.
What would settle it
An independent TSI proxy or reconstruction from a different method that shows a range in 50-year averages substantially larger than 1.04 W m^{-2} over the same three-millennia span.
Figures
read the original abstract
Solar irradiance measurements are limited to the last few decades, requiring reconstructions to assess solar variability on longer timescales and its impact on Earth's climate. We present the first physics-based reconstruction of total solar irradiance (TSI) at annual resolution over the last three millennia. The reconstruction is obtained by extending the SATIRE-T model beyond the telescopic era using recently published, annually resolved sunspot number series derived from cosmogenic isotope records. This yields a continuous, physics-based TSI record extending from the satellite era back over the last three millennia, with annual resolution throughout the pre-telescopic period. Over the full three-millennia interval, the reconstructed TSI exhibits a maximum difference of $1.04_{-0.2}^{+0.14}\,\mathrm{W\,m^{-2}}$, defined as the difference between the maximum and minimum of the 50-yr running mean values.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper extends the SATIRE-T semi-empirical model, previously calibrated on satellite-era TSI and telescopic sunspot data, to produce an annual-resolution reconstruction of total solar irradiance (TSI) over the last three millennia. Input sunspot numbers are taken from an inversion of cosmogenic 14C and 10Be production rates. The central quantitative result is a maximum difference of 1.04_{-0.2}^{+0.14} W m^{-2} between the maximum and minimum of the 50-year running-mean TSI values over the full interval.
Significance. A validated annual TSI reconstruction spanning three millennia would be a useful input for climate-sensitivity studies and for testing solar-dynamo models on centennial timescales. The physics-based character of SATIRE-T and the use of annually resolved proxy sunspot numbers are strengths relative to purely empirical long-term reconstructions. However, the significance is limited by the absence of explicit validation that the modern-calibrated flux-to-irradiance mapping remains unchanged when applied to the statistical properties of the pre-telescopic SSN series.
major comments (2)
- [Abstract and §2] Abstract and §2 (model description): the central claim that the 1.04 W m^{-2} range is robust rests on the assumption that SATIRE-T parameters (spot/facular contrasts and area-filling factors) calibrated on post-1978 data remain valid for the amplitude, cycle length, and grand-minimum depth statistics of the isotope-derived SSN series. No sensitivity tests that vary these parameters within their modern uncertainties or that compare the resulting TSI against independent proxies (e.g., 10Be-based reconstructions or tree-ring 14C) are reported.
- [§3] §3 (reconstruction results): the reported asymmetric uncertainties on the 1.04 W m^{-2} difference appear to reflect only the uncertainty in the 50-yr running-mean extrema; they do not propagate the uncertainty arising from possible secular changes in the SSN-to-magnetic-flux conversion or in the contrast ratios. A quantitative assessment of this additional systematic uncertainty is required for the quoted range to be interpreted as a total error budget.
minor comments (2)
- [Abstract] The definition of the reported maximum difference (difference between max and min of the 50-yr running mean) should be stated explicitly in the abstract and again in the results section for clarity.
- [Figures] Figure captions should indicate the time intervals covered by the satellite calibration data versus the proxy-driven reconstruction to avoid visual conflation of the two regimes.
Simulated Author's Rebuttal
We thank the referee for their constructive comments, which help clarify the robustness of our reconstruction. We respond to each major comment below.
read point-by-point responses
-
Referee: [Abstract and §2] Abstract and §2 (model description): the central claim that the 1.04 W m^{-2} range is robust rests on the assumption that SATIRE-T parameters (spot/facular contrasts and area-filling factors) calibrated on post-1978 data remain valid for the amplitude, cycle length, and grand-minimum depth statistics of the isotope-derived SSN series. No sensitivity tests that vary these parameters within their modern uncertainties or that compare the resulting TSI against independent proxies (e.g., 10Be-based reconstructions or tree-ring 14C) are reported.
Authors: SATIRE-T is a physics-based model in which spot and facular contrasts are determined from the wavelength-dependent radiative properties of magnetic features, quantities that are intrinsic to the Sun and independent of the statistical properties of the input sunspot series. The area-filling factors are likewise tied to the observed magnetic flux. Because the underlying physics does not change on millennial timescales, the same parameter set calibrated on the satellite era is applied to the isotope-derived SSN series. Nevertheless, we agree that an explicit demonstration of robustness is valuable. In the revised manuscript we will add sensitivity tests that vary the contrast ratios and filling-factor scaling within the uncertainties reported in the original SATIRE-T calibration papers and will include a brief comparison of the resulting TSI envelope with independent 10Be-based reconstructions. revision: yes
-
Referee: [§3] §3 (reconstruction results): the reported asymmetric uncertainties on the 1.04 W m^{-2} difference appear to reflect only the uncertainty in the 50-yr running-mean extrema; they do not propagate the uncertainty arising from possible secular changes in the SSN-to-magnetic-flux conversion or in the contrast ratios. A quantitative assessment of this additional systematic uncertainty is required for the quoted range to be interpreted as a total error budget.
Authors: The asymmetric uncertainties quoted for the 1.04 W m^{-2} range are obtained solely from the extrema of the 50-year running mean of the central reconstruction. We acknowledge that this does not yet incorporate systematic contributions from possible secular changes in the SSN-to-flux conversion or in the adopted contrasts. In the revised version we will add a quantitative assessment of these systematics, for example by recomputing the TSI range for the upper and lower bounds of the conversion factor and contrast uncertainties, and will present the resulting total error budget alongside the existing statistical range. revision: yes
Circularity Check
No significant circularity in the reconstruction chain
full rationale
The paper computes TSI by applying the pre-existing SATIRE-T model (with parameters fixed from modern satellite and telescopic data) to an externally derived, annually resolved sunspot-number time series obtained by inverting cosmogenic isotope production rates. The reported 1.04 W m^{-2} maximum difference is obtained by taking the max-minus-min of the 50-year running mean of this forward-modelled TSI series. This is a standard forward application of a calibrated model to independent proxy input; the output series is not redefined as or fitted to the same quantities that entered the model. No load-bearing step reduces by construction to the paper's own fitted values or to a self-citation chain whose validity is assumed inside the present work. The derivation therefore remains self-contained against external benchmarks.
Axiom & Free-Parameter Ledger
free parameters (1)
- SATIRE-T model parameters
axioms (1)
- domain assumption The physical relationship between sunspot number and total solar irradiance encoded in SATIRE-T holds for pre-telescopic solar conditions.
Reference graph
Works this paper leans on
-
[1]
Cyclic sunspot activity during the first millennium. A&A , author =. 2026 , pages =. doi:10.1051/0004-6361/202659038 , language =
-
[2]
doi:https://doi.org/10.1038/s41586-024-07679-4 , year =
Extreme solar storms and the quest for exact dating with radiocarbon , journal =. doi:https://doi.org/10.1038/s41586-024-07679-4 , year =
-
[3]
Type. ApJ , author =. 2017 , pages =. doi:10.3847/1538-4357/aa6afe , language =
-
[4]
Monthly Notices of the Royal Astronomical Society , author =
Missing. Monthly Notices of the Royal Astronomical Society , author =. 2020 , anote =. doi:10.1093/mnras/staa059 , abstract =
-
[5]
Monthly Notices of the Royal Astronomical Society , author =
How much hydrogen is in. Monthly Notices of the Royal Astronomical Society , author =. 2022 , anote =. doi:10.1093/mnras/stac088 , abstract =
-
[6]
Dorn-Wallenstein, Trevor Z. and Levesque, Emily M. and Davenport, James R. A. and Neugent, Kathryn F. and Morris, Brett M. and Bostroem, K. Azalee , month = jun, year =. The
-
[7]
Monthly Notices of the Royal Astronomical Society , author =
Towards a. Monthly Notices of the Royal Astronomical Society , author =. 2020 , anote =. doi:10.1093/mnras/staa3029 , abstract =
-
[8]
Core-Collapse Supernovae from 9 to 120 Solar Masses Based on Neutrino-powered Explosions
Core-. ApJ , author =. 2016 , anote =. doi:10.3847/0004-637X/821/1/38 , abstract =
work page internal anchor Pith review doi:10.3847/0004-637x/821/1/38 2016
-
[9]
Annual Review of Nuclear and Particle Science , author =
Explosion. Annual Review of Nuclear and Particle Science , author =. 2012 , anote =. doi:10.1146/annurev-nucl-102711-094901 , abstract =
-
[10]
A. A&A , author =. 2017 , anote =. doi:10.1051/0004-6361/201730698 , abstract =
-
[11]
Monthly Notices of the Royal Astronomical Society , author =
A simple approach to the supernova progenitor–explosion connection , volume =. Monthly Notices of the Royal Astronomical Society , author =. 2016 , pages =. doi:10.1093/mnras/stw1083 , abstract =
-
[12]
Monthly Notices of the Royal Astronomical Society , author =
Progenitor. Monthly Notices of the Royal Astronomical Society , author =. 2021 , anote =. doi:10.1093/mnras/stab1800 , abstract =
-
[13]
Baraffe, I. and Clarke, J. and Morison, A. and Vlaykov, D. G. and Constantino, T. and Goffrey, T. and Guillet, T. and Saux, A. Le and Pratt, J. , month = jan, year =. A study of convective core overshooting as a function of stellar mass based on two-dimensional hydrodynamical simulations , url =
-
[14]
The cosmic carbon footprint of massive stars stripped in binary systems , volume =. ApJ , author =. 2021 , anote =. doi:10.3847/1538-4357/ac2f44 , abstract =
-
[15]
A statistical study of 233 pulsar proper motions , volume =. Monthly Notices of the Royal Astronomical Society , author =. 2005 , anote =. doi:10.1111/j.1365-2966.2005.09087.x , abstract =
-
[16]
ApJ , author =. 2009 , pages =. doi:10.1088/0004-637X/703/2/2205 , language =
-
[17]
Sykes, B. and Mueller, B. and Cordero-Carrión, I. and Cerdá-Durán, P. and Novak, J. , month = oct, year =. A
-
[18]
Nature Methods , author =. 2020 , anote =. doi:10.1038/s41592-019-0686-2 , abstract =
-
[19]
How much. Monthly Notices of the Royal Astronomical Society , author =. 2012 , anote =. doi:10.1111/j.1365-2966.2012.20464.x , abstract =
-
[20]
Progenitors of. ApJ , author =. 2019 , anote =. doi:10.3847/1538-4357/ab4ad7 , abstract =
-
[21]
How. ApJ , author =. 2003 , anote =. doi:10.1086/375341 , abstract =
-
[22]
A. ApJ , author =. 2008 , pages =. doi:10.1086/590053 , abstract =
-
[23]
Stellar evolution , url =
-
[24]
Zeitschrift Angewandte Mathematik und Mechanik , author =
7. Zeitschrift Angewandte Mathematik und Mechanik , author =. 1925 , anote =. doi:10.1002/zamm.19250050212 , urldate =
-
[25]
Monthly Notices of the Royal Astronomical Society , author =
The excess of cool supergiants from contemporary stellar evolution models defies the metallicity-independent. Monthly Notices of the Royal Astronomical Society , author =. 2021 , pages =. doi:10.1093/mnras/stab383 , abstract =
-
[26]
The joy of writing a paper , language =
Sterk, P J and Rabe, K F , pages =. The joy of writing a paper , language =
-
[27]
Monthly Notices of the Royal Astronomical Society , author =
The search for failed supernovae with the. Monthly Notices of the Royal Astronomical Society , author =. 2015 , anote =. doi:10.1093/mnras/stv776 , abstract =
- [28]
-
[29]
Monthly Notices of the Royal Astronomical Society , author =
Numerical experiments to help understand cause and effect in massive star evolution , volume =. Monthly Notices of the Royal Astronomical Society , author =. 2022 , anote =. doi:10.1093/mnras/stac538 , abstract =
-
[30]
Stellar evolution with rotation , abstract =
Meynet, G and Maeder, A , pages =. Stellar evolution with rotation , abstract =
-
[31]
Nucleosynthesis in massive stars and the. Physics Reports , author =. 1993 , anote =. doi:10.1016/0370-1573(93)90058-L , abstract =
-
[32]
The. ApJ , author =. 2020 , pages =. doi:10.3847/1538-4357/ab6739 , abstract =
-
[33]
Astronomy and Astrophysics, Vol
Mass-dependent mass loss rates of. Astronomy and Astrophysics, Vol. 220, p. 135-143 (1989) , author =. 1989 , pages =
work page 1989
-
[34]
The. ApJ , author =. 1993 , anote =. doi:10.1086/172886 , abstract =
-
[35]
The. ApJ , author =. 1995 , anote =. doi:10.1086/175963 , abstract =
-
[36]
Synthesis of the. Reviews of Modern Physics , author =. 1957 , anote =. doi:10.1103/RevModPhys.29.547 , urldate =
-
[37]
An. ApJ , author =. 1993 , anote =. doi:10.1086/173117 , abstract =
-
[38]
and Black, Diane , month = jul, year =
Chamba, Nushkia and Knapen, Johan H. and Black, Diane , month = jul, year =. How to plan your astronomy research paper in ten steps , url =
- [39]
-
[40]
Prialnik, Dina , month = oct, year =. An
-
[41]
Different to the core:. A&A , author =. 2021 , pages =. doi:10.1051/0004-6361/202140506 , abstract =
-
[42]
Kozyreva, Alexandra and Klencki, Jakub and Filippenko, Alexei V. and Baklanov, Petr and Mironov, Alexey and Justham, Stephen and Chiavassa, Andrea , month = jul, year =. The circumstellar material around the
-
[43]
, year = 1955, month = jan, volume =
The. ApJ , author =. 1955 , pages =. doi:10.1086/145971 , abstract =
-
[44]
Astronomy and Astrophysics Review , author =
The yellow hypergiants , volume =. Astronomy and Astrophysics Review , author =. 1998 , pages =. doi:10.1007/s001590050009 , abstract =
- [45]
-
[46]
The. Annual Review of Astronomy and Astrophysics , author =. 2009 , pages =. doi:10.1146/annurev.astro.46.060407.145222 , abstract =
-
[47]
2019, ApJS, 243, 10, doi: 10.3847/1538-4365/ab2241
Modules for. ApJ Supplement Series , author =. 2019 , anote =. doi:10.3847/1538-4365/ab2241 , abstract =
-
[48]
Modules for. ApJ Supplement Series , author =. 2018 , anote =. doi:10.3847/1538-4365/aaa5a8 , abstract =
work page internal anchor Pith review doi:10.3847/1538-4365/aaa5a8 2018
-
[49]
Modules for Experiments in Stellar Astrophysics (MESA): Binaries, Pulsations, and Explosions
Modules for. ApJ Supplement Series , author =. 2015 , anote =. doi:10.1088/0067-0049/220/1/15 , abstract =
work page internal anchor Pith review doi:10.1088/0067-0049/220/1/15 2015
-
[50]
, archivePrefix = "arXiv", eprint =
Modules for. ApJ Supplement Series , author =. 2011 , anote =. doi:10.1088/0067-0049/192/1/3 , abstract =
-
[51]
Low-. ApJ , author =. 2005 , anote =. doi:10.1086/428642 , language =
-
[52]
The. ApJ , author =. 2020 , anote =. doi:10.3847/1538-4357/ab8cc1 , abstract =
-
[53]
Frontiers in Astronomy and Space Sciences , author =
Full-disc. Frontiers in Astronomy and Space Sciences , author =. 2022 , pages =. doi:10.3389/fspas.2022.1038949 , abstract =
-
[54]
Evolution of. ApJ , author =. 2010 , anote =. doi:10.1088/2041-8205/717/1/L62 , abstract =
-
[55]
The. ApJ , author =. 1996 , anote =. doi:10.1086/176778 , abstract =
-
[56]
The. ApJ , author =. 2013 , anote =. doi:10.1088/0004-637X/767/1/3 , abstract =
-
[57]
Monthly Notices of the Royal Astronomical Society , author =
Comparing compact object distributions from mass- and presupernova core structure-based prescriptions , volume =. Monthly Notices of the Royal Astronomical Society , author =. 2022 , anote =. doi:10.1093/mnras/stab3797 , abstract =
-
[58]
Monthly Notices of the Royal Astronomical Society , author =
Whatever happened to the progenitors of supernovae 2008cn, 2009kr and 2009md? , volume =. Monthly Notices of the Royal Astronomical Society , author =. 2015 , anote =. doi:10.1093/mnras/stu2658 , abstract =
-
[59]
Monthly Notices of the Royal Astronomical Society , author =
Dependence of convective boundary mixing on boundary properties and turbulence strength , volume =. Monthly Notices of the Royal Astronomical Society , author =. 2019 , pages =. doi:10.1093/mnras/stz312 , abstract =
-
[60]
Turbulent. ApJ , author =. 2007 , pages =. doi:10.1086/520318 , abstract =
-
[61]
A solar coronal loop in a box:. A&A , author =. 2022 , pages =. doi:10.1051/0004-6361/202141451 , abstract =
-
[62]
Impact of spatially correlated fluctuations in sunspots on metrics related to magnetic twist , volume =. A&A , author =. 2022 , pages =. doi:10.1051/0004-6361/202243357 , abstract =
-
[63]
Monthly Notices of the Royal Astronomical Society , author =
Modelling of σ. Monthly Notices of the Royal Astronomical Society , author =. 2014 , anote =. doi:10.1093/mnras/stu885 , abstract =
-
[64]
The. ApJ , author =. 1995 , anote =. doi:10.1086/175779 , abstract =
-
[65]
Strotjohann, Nora L. and Ofek, Eran O. and Gal-Yam, Avishay and Sollerman, Jesper and Chen, Ping and Yaron, Ofer and Zackay, Barak and Rehemtulla, Nabeel and Gris, Phillipe and Masci, Frank J. and Rusholme, Ben and Purdum, Josiah , month = feb, year =. Direct detection of supernova progenitor stars with
-
[66]
Observational. ApJ Letters , author =. 2023 , pages =. doi:10.3847/2041-8213/acb704 , abstract =
-
[67]
The. ApJ , author =. 2014 , anote =. doi:10.1088/0004-637X/783/1/10 , abstract =
-
[68]
Spin-down by dynamo action in simulated radiative stellar layers , volume =. Science , author =. 2023 , anote =. doi:10.1126/science.abk2169 , abstract =
-
[69]
Rotating massive main-sequence stars. A&A , author =. 2011 , pages =. doi:10.1051/0004-6361/201016113 , abstract =
-
[70]
Presupernova. ApJ , author =. 2000 , anote =. doi:10.1086/308158 , abstract =
-
[71]
Fallback and. ApJ , author =. 2008 , anote =. doi:10.1086/526404 , abstract =
-
[72]
Compact. ApJ , author =. 2012 , anote =. doi:10.1088/0004-637X/749/1/91 , abstract =
-
[73]
The. ApJ , author =. 2015 , anote =. doi:10.1088/0004-637X/801/2/90 , abstract =
-
[74]
A. ApJ , author =. 2018 , anote =. doi:10.3847/1538-4357/aac2da , abstract =
-
[75]
The. ApJ , author =. 2019 , anote =. doi:10.3847/1538-4357/ab1589 , abstract =
-
[76]
Royal Society Open Science , author =
Chemical element transport in stellar evolution models , volume =. Royal Society Open Science , author =. 2017 , anote =. doi:10.1098/rsos.170192 , abstract =
-
[77]
Binary. Science , author =. 2012 , pages =. doi:10.1126/science.1223344 , abstract =
-
[78]
The Astronomical Journal , author =. 2009 , pages =. doi:10.1088/0004-6256/137/2/3358 , language =
-
[79]
A. ApJ , author =. 2007 , pages =. doi:10.1086/522073 , abstract =
-
[80]
Progenitor-. ApJ , author =. 2012 , anote =. doi:10.1088/0004-637X/757/1/69 , abstract =
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.