Application of CCM SOCOL-AERv2-BE to cosmogenic beryllium isotopes: description and validation for polar regions

The short-living cosmogenic isotope 7Be, which is produced by cosmic rays in the atmosphere, is often used as a tracer for atmospheric dynamics, with precise and high-resolution measurements covering the recent decades. The long-living isotope 10Be, as measured in polar ice cores with an annual resolution, is a proxy for long-term cosmic-ray variability, whose signal can, however, be distorted by atmospheric transport and deposition that need to be properly modeled to be accounted for. While transport of 7Be can be modeled with high accuracy using the known meteorological fields, atmospheric transport of 10Be was typically modeled using case-study-specific simulations or simplified box models based on parameterizations. Thus, there is a need for a realistic model able to simulate atmospheric transport and deposition of beryllium with a focus on polar regions and (inter)annual timescales that is potentially able to operate in a self-consistent mode without the prescribed meteorology. Since measurements of 10Be are extremely laborious and hence scarce, it is difficult to compare model results directly with measurement data. On the other hand, the two beryllium isotopes are believed to have similar transport and deposition properties, being different only in production and lifetime, and thus the results of 7Be transport can be generally applied to 10Be. Here we present a new model, called CCM SOCOL-AERv2-BE, to trace isotopes of 7Be and 10Be in the atmosphere based on the chemistry–climate model (CCM) SOCOL (SOlar Climate Ozone Links), which has been improved by including modules for the production, deposition, and transport of 7Be and 10Be. Production of the isotopes was modeled for both galactic and solar cosmic rays by applying the CRAC (Cosmic Ray Atmospheric Cascade) model. Transport of 7Be was modeled without additional gravitational settling due to the submicron size of the background aerosol particles. An interactive deposition scheme was applied including both wet and dry deposition. Modeling was performed using a full nudging to the meteorological fields for the period of 2002–2008 with a spin-up period of 1996–2001. The modeled concentrations of 7Be in near-ground air were compared with the measured ones at a weekly time resolution in four nearly antipodal high-latitude locations: two in the Northern (Finland and Canada) and two in the Southern (Chile and the Kerguelen Islands) Hemisphere. The model results agree with the measurements in the absolute level within error bars, implying that the production, decay, and lateral deposition are correctly reproduced. The model also correctly reproduces the temporal variability of 7Be concentrations on annual and sub-annual scales, including the presence and absence of the annual cycle in the Northern and Southern Hemisphere, respectively. We also modeled the production and transport of 7Be for a major solar energetic particle event (SPE) on 20 January 2005, which appears insufficient to produce a measurable signal but may serve as a reference event for historically known extreme SPEs. Thus, a new full 3D time-dependent model, based on CCM SOCOL, of 7Be and 10Be atmospheric production, transport, and deposition has been developed. Comparison with real data on the 7Be concentration in the near-ground air validates the model and its accuracy.

Neogene continental denudation and the beryllium conundrum

Reconstructing Cenozoic history of continental silicate weathering is crucial for understanding Earth’s carbon cycle and greenhouse history. The question of whether continental silicate weathering increased during the late Cenozoic, setting the stage for glacial cycles, has remained controversial for decades. Whereas numerous independent proxies of weathering in ocean sediments (e.g., Li, Sr, and Os isotopes) have been interpreted to indicate that the continental silicate weathering rate increased in the late Cenozoic, beryllium isotopes in seawater have stood out as an important exception. Beryllium isotopes have been interpreted to indicate stable continental weathering and/or denudation rates over the last 12 Myr. Here we present a Be cycle model whose results show that variations in the 9Be weathering flux are counterbalanced by near-coastal scavenging while the cosmogenic 10Be flux from the upper atmosphere stays constant. As a result, predicted seawater 10Be/9Be ratios remain nearly constant even when global denudation and Be weathering rates increase by three orders of magnitude. Moreover, 10Be/9Be records allow for up to an 11-fold increase in Be weathering and denudation rates over the late Cenozoic, consistent with estimates from other proxies. The large increase in continental weathering indicated by multiple proxies further suggests that the increased CO2 consumption by continental weathering, driven by mountain-building events, was counterbalanced by other geological processes to prevent a runaway icehouse condition during the late Cenozoic. These processes could include enhanced carbonate dissolution via pyrite weathering, accelerated oxidation of fossil organic carbon, and/or reduced basalt weathering as the climate cooled.