Implications of two Holocene time-dependent geomagnetic models for comogenic nuclide production rate scaling

The geomagnetic field is a major influence on in situ cosmogenic nuclide production rates at a given location (in addition to atmospheric pressure and, to a lesser extent, solar modulation effects). A better understanding of how past fluctuations in these influences affected production rates should allow more accurate application of cosmogenic nuclides. As such, this work explores the cosmogenic nuclide production rate scaling implications of two recent time-dependent spherical harmonic geomagnetic models spanning the Holocene. Korte and Constable (2011, Phys. Earth Planet. Inter. 188, 247-259) and Korte et al. (2011, Earth Planet. Sci. Lett. 312, 497-505) recently updated earlier spherical harmonic paleomagnetic models with new paleomagnetic data from sediment cores in addition to new archeomagnetic and volcanic data. These updated models offer improved resolution and accuracy over the previous versions, in part due to increased temporal and spatial data coverage. In addition, Pavon-Carrasco et al. (2014, Earth Planet. ScL Lett. 388, 98-109) developed another time-dependent spherical harmonic model of the Holocene geomagnetic field, based solely on archeomagnetic and volcanic paleomagnetic data from the same underlying paleomagnetic database as the Korte et al. models, but extending to 14 ka. With the new models as input, trajectory-traced estimates of effective vertical cutoff rigidity (R-C - the standard method for ordering cosmic ray data) yield significantly different time-integrated scaling predictions when compared to each other and to results using the earlier models. In addition, predictions of each new model using R-C are tested empirically using recently published production rate calibration data for both Be-10 and He-3, and compared to predictions using corresponding time-varying geocentric dipolar R-C formulations and a static geocentric axial dipole (GAD) model. Results for the few calibration sites from geomagnetically sensitive regions suggest that the Pavon-Carrasco et al. (2014) time-varying dipolar model tends to predict sea level, high latitude production rates more in line with those from calibration sites not affected by geomagnetic variations. This suggests that uncertainties arising from hemispheric and temporal sampling biases in the Holocene spherical harmonic models considered here, combined with the currently limited spatial and temporal distribution of production rate calibration sites as empirical tests, limit the robustness of the non-dipole aspects of these models for production rate scaling. These analyses should be revisited as such models improve and additional calibration sites become available. (C) 2015 Elsevier B.V. All rights reserved. The geomagnetic field is a major influence on in situ cosmogenic nuclide production rates at a given location (in addition to atmospheric pressure and, to a lesser extent, solar modulation effects). A better understanding of how past fluctuations in these influences affected production rates should allow more accurate application of cosmogenic nuclides. As such, this work explores the cosmogenic nuclide production rate scaling implications of two recent time-dependent spherical harmonic geomagnetic models spanning the Holocene. Korte and Constable (2011, Phys. Earth Planet. Inter. 188, 247-259) and Korte et al. (2011, Earth Planet. Sci. Lett. 312, 497-505) recently updated earlier spherical harmonic paleomagnetic models with new paleomagnetic data from sediment cores in addition to new archeomagnetic and volcanic data. These updated models offer improved resolution and accuracy over the previous versions, in part due to increased temporal and spatial data coverage. In addition, Pavon-Carrasco et al. (2014, Earth Planet. ScL Lett. 388, 98-109) developed another time-dependent spherical harmonic model of the Holocene geomagnetic field, based solely on archeomagnetic and volcanic paleomagnetic data from the same underlying paleomagnetic database as the Korte et al. models, but extending to 14 ka. With the new models as input, trajectory-traced estimates of effective vertical cutoff rigidity (R-C - the standard method for ordering cosmic ray data) yield significantly different time-integrated scaling predictions when compared to each other and to results using the earlier models. In addition, predictions of each new model using R-C are tested empirically using recently published production rate calibration data for both Be-10 and He-3, and compared to predictions using corresponding time-varying geocentric dipolar R-C formulations and a static geocentric axial dipole (GAD) model. Results for the few calibration sites from geomagnetically sensitive regions suggest that the Pavon-Carrasco et al. (2014) time-varying dipolar model tends to predict sea level, high latitude production rates more in line with those from calibration sites not affected by geomagnetic variations. This suggests that uncertainties arising from hemispheric and temporal sampling biases in the Holocene spherical harmonic models considered here, combined with the currently limited spatial and temporal distribution of production rate calibration sites as empirical tests, limit the robustness of the non-dipole aspects of these models for production rate scaling. These analyses should be revisited as such models improve and additional calibration sites become available. (C) 2015 Elsevier B.V. All rights reserved.