Beryllium concentrations and deposition rates in ODP holes from the north Atlantic

Increases in the production rate of cosmogenic radionuclides associated with geomagnetic excursions have been used as global tie-points for correlation between records of past climate from marine and terrestrial archives. We have investigated the relative timing of variations in 10Be production rate and the corresponding palaeomagnetic signal during one of the largest Pleistocene excursions, the Iceland Basin (IB) event (ca. 190 kyr), as recorded in two marine sediment cores (ODP Sites 1063 and 983) with high sedimentation rates. Variations in 10Be production rate during the excursion were estimated by use of 230Thxs normalized 10Be deposition rates and authigenic 10Be/9Be. Resulting 10Be production rates are compared with high-resolution records of geomagnetic field behaviour acquired from the same discrete samples. We find no evidence for a significant lock-in depth of the palaeomagnetic signal in these high sedimentation-rate cores. Apparent lock-in depths in other cores may sometimes be the result of lower sample resolution. Our results also indicate that the period of increased 10Be production during the IB excursion lasted longer and, most likely, started earlier than the corresponding palaeomagnetic anomaly, in accordance with previous observations that polarity transitions occur after periods of reduced geomagnetic field intensity prior to the transition. The lack of evidence in this study for a significant palaeomagnetic lock-in depth suggests that there is no systematic offset between the 10Be signal and palaeomagnetic anomalies associated with excursions and reversals, with significance for the global correlation of climate records from different archives.

Radiocarbon ages, grain sizes, paleomagentic properties and XRF data of Lake Nam Co sediments, Tibetan Plateau

In September 2008 several cores (68 cm-115 cm length) (water depth: 93 m) were retrieved from Lake Nam Co (southern-central Tibetan Plateau; 4718 m a.s.l.). This study focuses on the interpretation of high-resolution (partly 0.2 cm) data from three gravity cores and the upper part of a 10.4 m long piston core, i.e., the past 4000 cal BP in terms of lake level changes, hydrological variations in the catchment area and consequently variations in monsoon strength. A wide spectrum of sedimentological, geochemical and mineralogical investigations was carried out. Results are presented for XRF core-scans, grain size distribution, XRD-measurements and SEM-image analyses. These data are complemented by an age-depth model using 210Pb and 137Cs analyses as well as eleven AMS-14C-ages. This model is supported by excellent agreement between secular variations determined on one of the gravity cores to geomagnetic field models. This is a significant improvement of the chronology as most catchments of lacustrine systems on the Tibetan Plateau contain carbonates resulting in an unknown reservoir effect for radiocarbon dates. The good correlation of our record to the geomagnetic field models confirms our age-depth model and indicates only insignificant changes in the reservoir effect throughout the last 4 ka. High (summer-) monsoonal activity, i.e. moist environmental conditions, was detected in our record between approximately 4000 and 1950 cal BP as well as between 1480 and 1200 cal BP. Accordingly, lower monsoon activity prevails in periods between the two intervals and thereafter. This pattern shows a good correlation to the variability of the Indian Ocean Summer Monsoon (IOSM) as recorded in a peat bog ~1000 km in NE direction from Lake Nam Co. This is the first time that such a supra regional homogenous monsoon activity is shown on the Tibetan Plateau and beyond. Finally our data show a significant lake level rise after the Little Ice Age (LIA) in Lake Nam Co which is suggested to be linked to glacier melting in consequence of rising temperatures occurring on the whole Tibetan Plateau during this time.

(Table 2) Radiocarbon dates of sediment cores obtained during Bio Hespérides cruise SVAIS and OGS Explora cruise EGLACOM between 2007-2008

A high-resolution paleomagnetic and rock magnetic study has been carried out on sediment cores collected in glaciomarine silty-clay sequences from the continental shelf and slope of the southern Storfjorden trough-mouth fan, on the northwestern Barents Sea continental margin. The Storfjorden sedimentary system was investigated during the SVAIS and EGLACOM cruises, when 10 gravity cores, with a variable length from 1.03 m to 6.41 m, were retrieved. Accelerator mass spectrometry (AMS) 14C analyses on 24 samples indicate that the cores span a time interval that includes the Holocene, the last deglaciation phase and in some cores the last glacial maximum. The sediments carry a well-defined characteristic remanent magnetization and have a valuable potential to reconstruct the paleosecular variation (PSV) of the geomagnetic field, including relative paleointensity (RPI) variations. The paleomagnetic data allow reconstruction of past dynamics and amplitude of the geomagnetic field variations at high northern latitudes (75°-76° N). At the same time, the rock magnetic and paleomagnetic data allow a high-resolution correlation of the sedimentary sequences and a refinement of their preliminary age models. The Holocene PSV and RPI records appear particularly sound, since they are consistent between cores and they can be correlated to the closest regional stacking curves (UK PSV, FENNOSTACK and FENNORPIS) and global geomagnetic model for the last 7 ka (CALS7k.2). The computed amplitude of secular variation is lower than that outlined by some geomagnetic field models, suggesting that it has been almost independent from latitude during the Holocene.

Paleomagnetic and beryllium isotope record of West Equatorial Pacific sediments

The reliability of paleomagnetic records as proxies of the geomagnetic field intensity is still a matter of controversy since volcanic materials hardly provide continuous records, and marine sediments are suspected to carry a remanence biased by post-depositional realignments and/or by overprints. Such long standing debate emphasizes the need for the development of methods independent from paleomagnetism to decipher geomagnetic intensity variations. High resolution measurements of authigenic 10Be/9Be along with a detailed sedimentary record of directional and relative paleointensity variations evidence, over the 0.6-1.3 Ma time interval, frequent and recurrent excursions or short events in the late Matuyama and the early Brunhes epochs, among which two Brunhes-Matuyama reversal precursors and an intra-Jaramillo excursion. The results of this study confirm the idea of a highly unstable geomagnetic field as suggested by paleomagnetic evidences.

Energy estimation of cosmic rays with the Engineering Radio Array of the Pierre Auger Observatory

The Auger Engineering Radio Array (AERA) is part of the Pierre Auger Observatory and is used to detect the radio emission of cosmic-ray air showers. These observations are compared to the data of the surface detector stations of the Observatory, which provide well-calibrated information on the cosmic-ray energies and arrival directions. The response of the radio stations in the 30-80 MHz regime has been thoroughly calibrated to enable the reconstruction of the incoming electric field. For the latter, the energy deposit per area is determined from the radio pulses at each observer position and is interpolated using a two-dimensional function that takes into account signal asymmetries due to interference between the geomagnetic and charge-excess emission components. The spatial integral over the signal distribution gives a direct measurement of the energy transferred from the primary cosmic ray into radio emission in the AERA frequency range. We measure 15.8 MeV of radiation energy for a 1 EeV air shower arriving perpendicularly to the geomagnetic field. This radiation energy-corrected for geometrical effects-is used as a cosmic-ray energy estimator. Performing an absolute energy calibration against the surface-detector information, we observe that this radio-energy estimator scales quadratically with the cosmic-ray energy as expected for coherent emission. We find an energy resolution of the radio reconstruction of 22% for the data set and 17% for a high-quality subset containing only events with at least five radio stations with signal.

Measurement of the radiation energy in the radio signal of extensive air showers as a universal estimator of cosmic-ray energy.

We measure the energy emitted by extensive air showers in the form of radio emission in the frequency range from 30 to 80 MHz. Exploiting the accurate energy scale of the Pierre Auger Observatory, we obtain a radiation energy of 15.8 +/- 0.7 (stat) +/- 6.7 (syst) MeV for cosmic rays with an energy of 1 EeV arriving perpendicularly to a geomagnetic field of 0.24 G, scaling quadratically with the cosmic-ray energy. A comparison with predictions from state-of-the-art first-principles calculations shows agreement with our measurement. The radiation energy provides direct access to the calorimetric energy in the electromagnetic cascade of extensive air showers. Comparison with our result thus allows the direct calibration of any cosmic-ray radio detector against the well-established energy scale of the Pierre Auger Observatory.

(Table 1) Types of giant magnetofossils in ODP holes

Magnetotactic bacteria biomineralize magnetic minerals with precisely controlled size, morphology, and stoichiometry. These cosmopolitan bacteria are widely observed in aquatic environments. If preserved after burial, the inorganic remains of magnetotactic bacteria act as magnetofossils that record ancient geomagnetic field variations. They also have potential to provide paleoenvironmental information. In contrast to conventional magnetofossils, giant magnetofossils (most likely produced by eukaryotic organisms) have only been reported once before from Paleocene-Eocene Thermal Maximum (PETM; 55.8 Ma) sediments on the New Jersey coastal plain. Here, using transmission electron microscopic observations, we present evidence for abundant giant magnetofossils, including previously reported elongated prisms and spindles, and new giant bullet-shaped magnetite crystals, in the Southern Ocean near Antarctica, not only during the PETM, but also shortly before and after the PETM. Moreover, we have discovered giant bullet-shaped magnetite crystals from the equatorial Indian Ocean during the Mid-Eocene Climatic Optimum (~40 Ma). Our results indicate a more widespread geographic, environmental, and temporal distribution of giant magnetofossils in the geological record with a link to "hyperthermal" events. Enhanced global weathering during hyperthermals, and expanded suboxic diagenetic environments, probably provided more bioavailable iron that enabled biomineralization of giant magnetofossils. Our micromagnetic modelling indicates the presence of magnetic multi-domain (i.e., not ideal for navigation) and single domain (i.e., ideal for navigation) structures in the giant magnetite particles depending on their size, morphology and spatial arrangement. Different giant magnetite crystal morphologies appear to have had different biological functions, including magnetotaxis and other non-navigational purposes. Our observations suggest that hyperthermals provided ideal conditions for giant magnetofossils, and that these organisms were globally distributed. Much more work is needed to understand the interplay between magnetofossil morphology, climate, nutrient availability, and environmental variability.

(Table 1) Age markers and errors for the last 1215 ka in IODP Exp302

Despite its importance in the global climate system, age-calibrated marine geologic records reflecting the evolution of glacial cycles through the Pleistocene are largely absent from the central Arctic Ocean. This is especially true for sediments older than 200 ka. Three sites cored during the Integrated Ocean Drilling Program's Expedition 302, the Arctic Coring Expedition (ACEX), provide a 27 m continuous sedimentary section from the Lomonosov Ridge in the central Arctic Ocean. Two key biostratigraphic datums and constraints from the magnetic inclination data are used to anchor the chronology of these sediments back to the base of the Cobb Mountain subchron (1215 ka). Beyond 1215 ka, two best fitting geomagnetic models are used to investigate the nature of cyclostratigraphic change. Within this chronology we show that bulk and mineral magnetic properties of the sediments vary on predicted Milankovitch frequencies. These cyclic variations record ''glacial'' and ''interglacial'' modes of sediment deposition on the Lomonosov Ridge as evident in studies of ice-rafted debris and stable isotopic and faunal assemblages for the last two glacial cycles and were used to tune the age model. Potential errors, which largely arise from uncertainties in the nature of downhole paleomagnetic variability, and the choice of a tuning target are handled by defining an error envelope that is based on the best fitting cyclostratigraphic and geomagnetic solutions.

(Table 2) U-channel paleomagnetic record of ODP Site 178-1098

Thick Holocene sedimentary sections (>45 m) cored in the Palmer Deep by the United States Antarctic Program (USAP) and during Ocean Drilling Program (ODP) Leg 178 provide the first opportunity to examine past geomagnetic field behavior at high southern latitudes. After removal of a low-coercivity drilling overprint the sediments display a stable, single-component remanent magnetization. Two short cores that recovered the uppermost 2.6 m of sediment have inclinations that fluctuate about the present day inclination (-57°) measured at Faraday Station, and several features with wavelengths of 10 to 20 cm appear to be correlative. However, shipboard measurements of inclination fluctuations on split-core samples from three holes drilled at ODP Site 1098 do not correlate well with each other, even though the intensity and susceptibility data correlate very well and the overall mean inclination for cores from each hole is consistent with the expected geocentric axial dipole (GAD) inclination. The correlation is improved dramatically by using inclinations measured on u-channels taken from the pristine center of a split core. Consequently, the anomalous directions and the resulting poor between-hole correlation of inclinations obtained from shipboard data can be attributed to coring-induced deformation, which is common on the outer edge of ODP piston cores, and/or measurement artifacts in the split-core data. Our preferred inclination record is thus derived from u-channel results. The upper ~25 m represents continuous sedimentation over the past 9000 yr, with an average sedimentation rate exceeding 250 cm/kyr (0.25 cm/yr). Given that remanence measurements on u-channels average over an interval <7 cm long, we obtained independent measurements of the paleo-geomagnetic field that average over only ~30 yr. This high-resolution record is characterized by an inclination that fluctuates within +/-15° of the current GAD inclination.

Centennial changes in the heliospheric magnetic field and open solar flux: the consensus view from geomagnetic data and cosmogenic isotopes and its implications

Svalgaard and Cliver (2010) recently reported a consensus between the various reconstructions of the heliospheric field over recent centuries. This is a significant development because, individually, each has uncertainties introduced by instrument calibration drifts, limited numbers of observatories, and the strength of the correlations employed. However, taken collectively, a consistent picture is emerging. We here show that this consensus extends to more data sets and methods than reported by Svalgaard and Cliver, including that used by Lockwood et al. (1999), when their algorithm is used to predict the heliospheric field rather than the open solar flux. One area where there is still some debate relates to the existence and meaning of a floor value to the heliospheric field. From cosmogenic isotope abundances, Steinhilber et al. (2010) have recently deduced that the near-Earth IMF at the end of the Maunder minimum was 1.80 ± 0.59 nT which is considerably lower than the revised floor of 4nT proposed by Svalgaard and Cliver. We here combine cosmogenic and geomagnetic reconstructions and modern observations (with allowance for the effect of solar wind speed and structure on the near-Earth data) to derive an estimate for the open solar flux of (0.48 ± 0.29) × 1014 Wb at the end of the Maunder minimum. By way of comparison, the largest and smallest annual means recorded by instruments in space between 1965 and 2010 are 5.75 × 1014 Wb and 1.37 × 1014 Wb, respectively, set in 1982 and 2009, and the maximum of the 11 year running means was 4.38 × 1014 Wb in 1986. Hence the average open solar flux during the Maunder minimum is found to have been 11% of its peak value during the recent grand solar maximum.

Reconstruction of geomagnetic activity and near-Earth interplanetary conditions over the past 167 yr – Part 2: A new reconstruction of the interplanetary magnetic field

We present a new reconstruction of the interplanetary magnetic field (IMF, B) for 1846–2012 with a full analysis of errors, based on the homogeneously constructed IDV(1d)composite of geomagnetic activity presented in Part 1 (Lockwood et al., 2013a). Analysis of the dependence of the commonly used geomagnetic indices on solar wind parameters is presented which helps explain why annual means of interdiurnal range data, such as the new composite, depend only on the IMF with only a very weak influence of the solar wind flow speed. The best results are obtained using a polynomial (rather than a linear) fit of the form B = χ · (IDV(1d) − β)α with best-fit coefficients χ = 3.469, β = 1.393 nT, and α = 0.420. The results are contrasted with the reconstruction of the IMF since 1835 by Svalgaard and Cliver (2010).

The heliospheric magnetic field

The heliospheric magnetic field (HMF) is the extension of the coronal magnetic field carried out into the solar system by the solar wind. It is the means by which the Sun interacts with planetary magnetospheres and channels charged particles propagating through the heliosphere. As the HMF remains rooted at the solar photosphere as the Sun rotates, the large-scale HMF traces out an Archimedean spiral. This pattern is distorted by the interaction of fast and slow solar wind streams, as well as the interplanetary manifestations of transient solar eruptions called coronal mass ejections. On the smaller scale, the HMF exhibits an array of waves, discontinuities, and turbulence, which give hints to the solar wind formation process. This review aims to summarise observations and theory of the small- and large-scale structure of the HMF. Solar-cycle and cycle-to-cycle evolution of the HMF is discussed in terms of recent spacecraft observations and pre-spaceage proxies for the HMF in geomagnetic and galactic cosmic ray records.

Reconstruction of geomagnetic activity and near-Earth interplanetary conditions over the past 167 yr - Part 4: Near-Earth solar wind speed, IMF, and open solar flux

In the concluding paper of this tetralogy, we here use the different geomagnetic activity indices to reconstruct the near-Earth interplanetary magnetic field (IMF) and solar wind flow speed, as well as the open solar flux (OSF) from 1845 to the present day. The differences in how the various indices vary with near-Earth interplanetary parameters, which are here exploited to separate the effects of the IMF and solar wind speed, are shown to be statistically significant at the 93% level or above. Reconstructions are made using four combinations of different indices, compiled using different data and different algorithms, and the results are almost identical for all parameters. The correction to the aa index required is discussed by comparison with the Ap index from a more extensive network of mid-latitude stations. Data from the Helsinki magnetometer station is used to extend the aa index back to 1845 and the results confirmed by comparison with the nearby St Petersburg observatory. The optimum variations, using all available long-term geomagnetic indices, of the near-Earth IMF and solar wind speed, and of the open solar flux, are presented; all with ±2sigma� uncertainties computed using the Monte Carlo technique outlined in the earlier papers. The open solar flux variation derived is shown to be very similar indeed to that obtained using the method of Lockwood et al. (1999).

Reconstruction of geomagnetic activity and near-Earth interplanetary conditions over the past 167 yr - Part 3: improved representation of solar cycle 11

Svalgaard (2014) has recently pointed out that the calibration of the Helsinki magnetic observatory’s H component variometer was probably in error in published data for the years 1866–1874.5 and that this makes the interdiurnal variation index based on daily means, IDV(1d), (Lockwood et al., 2013a), and the interplanetary magnetic field strength derived from it (Lockwood et al., 2013b), too low around the peak of solar cycle 11. We use data from the modern Nurmijarvi station, relatively close to the site of the original Helsinki Observatory, to confirm a 30% underestimation in this interval and hence our results are fully consistent with the correction derived by Svalgaard. We show that the best method for recalibration uses the Helsinki Ak(H) and aa indices and is accurate to ±10 %. This makes it preferable to recalibration using either the sunspot number or the diurnal range of geomagnetic activity which we find to be accurate to ±20 %. In the case of Helsinki data during cycle 11, the two recalibration methods produce very similar corrections which are here confirmed using newly digitised data from the nearby St Petersburg observatory and also using declination data from Helsinki. However, we show that the IDV index is, compared to later years, too similar to sunspot number before 1872, revealing independence of the two data series has been lost; either because the geomagnetic data used to compile IDV has been corrected using sunspot numbers, or vice versa, or both. We present corrected data sequences for both the IDV(1d) index and the reconstructed IMF (interplanetary magnetic field).We also analyse the relationship between the derived near-Earth IMF and the sunspot number and point out the relevance of the prior history of solar activity, in addition to the contemporaneous value, to estimating any “floor” value of the near-Earth interplanetary field.

Open solar flux estimates from near-Earth measurements of the interplanetary magnetic field: comparison of the first two perihelion passes of the Ulysses spacecraft

Results from all phases of the orbits of the Ulysses spacecraft have shown that the magnitude of the radial component of the heliospheric field is approximately independent of heliographic latitude. This result allows the use of near- Earth observations to compute the total open flux of the Sun. For example, using satellite observations of the interplanetary magnetic field, the average open solar flux was shown to have risen by 29% between 1963 and 1987 and using the aa geomagnetic index it was found to have doubled during the 20th century. It is therefore important to assess fully the accuracy of the result and to check that it applies to all phases of the solar cycle. The first perihelion pass of the Ulysses spacecraft was close to sunspot minimum, and recent data from the second perihelion pass show that the result also holds at solar maximum. The high level of correlation between the open flux derived from the various methods strongly supports the Ulysses discovery that the radial field component is independent of latitude. We show here that the errors introduced into open solar flux estimates by assuming that the heliospheric field’s radial component is independent of latitude are similar for the two passes and are of order 25% for daily values, falling to 5% for averaging timescales of 27 days or greater. We compare here the results of four methods for estimating the open solar flux with results from the first and second perehelion passes by Ulysses. We find that the errors are lowest (1–5% for averages over the entire perehelion passes lasting near 320 days), for near-Earth methods, based on either interplanetary magnetic field observations or the aa geomagnetic activity index. The corresponding errors for the Solanki et al. (2000) model are of the order of 9–15% and for the PFSS method, based on solar magnetograms, are of the order of 13–47%. The model of Solanki et al. is based on the continuity equation of open flux, and uses the sunspot number to quantify the rate of open flux emergence. It predicts that the average open solar flux has been decreasing since 1987, as Correspondence to: M. Lockwood (m.lockwood@rl.ac.uk) is observed in the variation of all the estimates of the open flux. This decline combines with the solar cycle variation to produce an open flux during the second (sunspot maximum) perihelion pass of Ulysses which is only slightly larger than that during the first (sunspot minimum) perihelion pass.