A reconstruction of radiocarbon production and total solar irradiance from the Holocene <sup>14</sup>C and CO<sub>2</sub> records: implications of data and model uncertainties
Abstract. Past atmospheric CO2 concentrations reconstructed from polar ice cores combined with its Δ14C signature as conserved in tree-rings provide important information both on the cycling of carbon as well as the production of radiocarbon (Q) in the atmosphere. The latter is modulated by changes in the strength of the magnetic field enclosed in the solar wind and is a proxy for past changes in solar activity. We perform transient carbon-cycle simulations spanning the past 21 kyr using Bern3D-LPX, a fully featured Earth System Model of Intermediate Complexity (EMIC) with a 3-D ocean, sediment and a dynamic vegetation model. Using the latest atmospheric IntCal09/SHCal04 radiocarbon records, we reconstruct the Holocene radiocarbon fluxes and the total production rate. Our carbon-cycle based modern estimate of Q ≈ 1.7 atoms cm−2 s−1 is lower than previously reported by Masarik and Beer (2009) and more in line with Kovaltsov et al. (2012). Q is then translated into the solar modulation potential (Φ) using the latest geomagnetic field reconstruction and linked to a recent reanalysis of early instrumental data. In contrast to earlier reconstructions, our record suggests that periods of high solar activity (>600 MeV) were quite common not only in recent millennia but throughout the Holocene. Solar activity in our decadally-smoothed record is during 28% of the time higher than the modern average of 650 MeV during the past 9 ka. But due to considerable uncertainties in the normalization of Φ to instrumental data, the absolute value of Φ remains weakly constrained. Further, our simulations with a spatially resolved model (taking the interhemispheric Δ14C gradient into account) show that reconstructions that rely on the Northern Hemisphere 14C record only are biased towards low values during the Holocene. Notable deviations on decadal-to-centennial time scales are also found in comparison with earlier reconstructions. In a last step, past total solar irradiance (TSI) is quantified using a recently published Φ-TSI relationship yielding small changes in Holocene TSI of order 1 W m−2 with a Maunder Minimum irradiance reduction of 0.85 ± 0.17 W m−2. Future extension of TSI using autoregressive modeling suggest a declining solar activity in the next decades towards average Holocene conditions. Past TSI changes are finally translated into changes in surfaces atmospheric temperature (SAT) by forcing the Bern3D-LPX model with our new TSI record, yielding SAT anomalies of less than 0.1 K.