(Table S1) Age determination of ODP Site 165-1002

A series of 14C measurements in Ocean Drilling Program cores from the tropical Cariaco Basin, which have been correlated to the annual-layer counted chronology for the Greenland Ice Sheet Project 2 (GISP2) ice core, provides a high-resolution calibration of the radiocarbon time scale back to 50,000 years before the present. Independent radiometric dating of events correlated to GISP2 suggests that the calibration is accurate. Reconstructed 14C activities varied substantially during the last glacial period, including sharp peaks synchronous with the Laschamp and Mono Lake geomagnetic field intensity minimal and cosmogenic nuclide peaks in ice cores and marine sediments. Simulations with a geochemical box model suggest that much of the variability can be explained by geomagnetically modulated changes in 14C production rate together with plausible changes in deep-ocean ventilation and the global carbon cycle during glaciation.

Paleointensity of various DSDP sites (Table 3)

The existing database for paleointensity estimates of the ancient geomagnetic field contains more than 1500 data points collected through decades of effort. Despite the huge amount of work put into obtaining these data, there remains a strong bias in the age and global distribution of the data toward very young results from a few locations. Also, few of the data meet strict criteria for reliability and most are of unknown quality. In order to improve the age and spatial distribution of the paleointensity database, we have carried out paleointensity experiments on submarine basaltic glasses from a number of DSDP sites. Of particular interest are the sites that provide paleointensity data spanning the time period 0.3-5 Ma, a time of relatively few high quality published data points. Our new data are concordant with contemporaneous data from the published literature that meet minimum acceptance criteria, and the combined data set yields an average dipole moment of 5.49 +/- 2.36*10**22 Am**2. This average value is comparable to the average paleofield for the period 5-160 Ma (4.2 +/- 2.3*10**22 Am**2) (Juarez et al., 1998, doi:10.1038/29746) and is substantially less than the value of approximately 8*10**22 Am**2 often quoted for the last 5 Myr (e.g. McFadden and McElhinny (1982) J. Geomagn. Geoelectr. 34, 163-189; Goguitchaichvili et al., 1999, doi:10.1016/S0012-821X(99)00010-2).

Age determination of three sediment cores

A total of five sediment cores from three sites, the Arctic Ocean, the Fram Strait and the Greenland Sea, yielded evidence for geomagnetic reversal excursions and associated strong lows in relative palaeointensity during oxygen isotope stages 2 and 3. A general similarity of the obtained relative palaeointensity curves to reference data can be observed. However, in the very detail, results from this high-resolution study differ from published records in a way that the prominent Laschamp excursion is clearly characterized by a significant field recovery when reaching the steepest negative inclinations, whereas only the N-R and R-N transitions are associated with the lowest values. Two subsequent excursions also reach nearly reversed inclinations but without any field recovery at that state. A total of 41 accelerator mass spectrometry (AMS) 14C ages appeared to allow a better age determination of these three directional excursions and related relative palaeointensity variations. However, although the three sites yielded more or less consistent chronological as well as palaeomagnetic results a comparison to another site, PS2644 in the Iceland Sea, revealed significant divergences in the ages of the geomagnetic field excursions of up to 4 ka even on basis of uncalibrated AMS 14C ages. This shift to older 14C ages cannot be explained by a time-transgressive character of the excursions, because the distance between the sites is small when compared with the size of and the distance to the geodynamo in the Earth's outer core. The most likely explanation is a difference of reservoir ages and/or mixing with old 14C-depleted CO2 from glacier ice expelled from Greenland at site PS2644.

(Table 1) Depth and revised depth from ODP Sites 169-1033 and 169-1034

The Holocene section in Saanich Inlet, Vancouver Island, British Columbia, is 50-70 m thick. Cores from Saanich Inlet obtained during Leg 169S of the Ocean Drilling Program afford an excellent opportunity to obtain an ultrahigh-resolution paleomagnetic and environmental magnetic record for the Holocene and Late Pleistocene of western Canada. We have used an automated, long-core cryogenic magnetometer to study over 380 m of continuous u-channel samples from ODP Sites 1033 and 1034, the two sites that constitute Leg 169S. Holocene records of paleomagnetic inclination and intensity show excellent intra-site correlation and can be used to fine-tune the lithologic correlation among cores from each site. The Late Pleistocene magnetic records provide a means of intra-site correlation of the otherwise featureless marine clay. Near the Holocene/Late Pleistocene boundary, both sites contain a magnetic intensity feature that is interpreted as a Missoula-type flood event on the Fraser River. The composite Holocene inclination records from the two sites are quite similar and provide a means of comparing current age-models that are based on radiocarbon dating of material from each site. This comparison shows only minor differences in the available age-models. It also provides strong evidence that the sediments of Saanich Inlet represent a reliable record of geomagnetic field behavior.

Component NRM directions for sediment core HLY0503-06JPC

The Arctic oceans have been fertile ground for the recording of apparent excursions of the geomagnetic field, implying that the high latitude field had unusual characteristics at least over the last 1-2 Myrs. Alternating field demagnetization of the natural remanent magnetization (NRM) of Core HLY0503-6JPC from the Mendeleev Ridge (Arctic Ocean) implies the presence of primary magnetizations with negative inclination apparently recording excursions in sediments deposited during the Brunhes Chron. Thermal demagnetization, on the other hand, indicates the presence of multiple (often anti-parallel) magnetization components with negative inclination components having blocking temperatures predominantly, but not entirely, below ~ 350 °C. Thermo-magnetic tests, X-ray diffraction (XRD) and scanning electron microscopy (SEM) indicate that the negative inclination components are carried by titanomaghemite, presumably formed by seafloor oxidation of titanomagnetite. The titanomaghemite apparently carries a chemical remanent magnetization (CRM) that is partially self-reversed relative to the detrital remanent magnetization (DRM) carried by the host titanomagnetite. The partial self-reversal could have been accomplished by ionic ordering during oxidation, thereby changing the balance of the magnetic moments in the ferrimagnetic sublattices.

Atlantic 231Pa/230Th and biogenic opal compilation of the Holocene and the LGM

The strength and geometry of the Atlantic meridional overturning circulation is tightly coupled to climate on glacial-interglacial and millennial timescales, but has proved difficult to reconstruct, particularly for the Last Glacial Maximum. Today, the return flow from the northern North Atlantic to lower latitudes associated with the Atlantic meridional overturning circulation reaches down to approximately 4,000 m. In contrast, during the Last Glacial Maximum this return flow is thought to have occurred primarily at shallower depths. Measurements of sedimentary 231Pa/230Th have been used to reconstruct the strength of circulation in the North Atlantic Ocean, but the effects of biogenic silica on 231Pa/230Th-based estimates remain controversial. Here we use measurements of 231Pa/230Th ratios and biogenic silica in Holocene-aged Atlantic sediments and simulations with a two-dimensional scavenging model to demonstrate that the geometry and strength of the Atlantic meridional overturning circulation are the primary controls of 231Pa/230Th ratios in modern Atlantic sediments. For the glacial maximum, a simulation of Atlantic overturning with a shallow, but vigorous circulation and bulk water transport at around 2,000 m depth best matched observed glacial Atlantic 231Pa/230Th values. We estimate that the transport of intermediate water during the Last Glacial Maximum was at least as strong as deep water transport today.

Analysis of magnetoencephaoographic data as a nonlinear dynamical system

This thesis presents the results from an investigation into the merits of analysing Magnetoencephalographic (MEG) data in the context of dynamical systems theory. MEG is the study of both the methods for the measurement of minute magnetic flux variations at the scalp, resulting from neuro-electric activity in the neocortex, as well as the techniques required to process and extract useful information from these measurements. As a result of its unique mode of action - by directly measuring neuronal activity via the resulting magnetic field fluctuations - MEG possesses a number of useful qualities which could potentially make it a powerful addition to any brain researcher's arsenal. Unfortunately, MEG research has so far failed to fulfil its early promise, being hindered in its progress by a variety of factors. Conventionally, the analysis of MEG has been dominated by the search for activity in certain spectral bands - the so-called alpha, delta, beta, etc that are commonly referred to in both academic and lay publications. Other efforts have centred upon generating optimal fits of "equivalent current dipoles" that best explain the observed field distribution. Many of these approaches carry the implicit assumption that the dynamics which result in the observed time series are linear. This is despite a variety of reasons which suggest that nonlinearity might be present in MEG recordings. By using methods that allow for nonlinear dynamics, the research described in this thesis avoids these restrictive linearity assumptions. A crucial concept underpinning this project is the belief that MEG recordings are mere observations of the evolution of the true underlying state, which is unobservable and is assumed to reflect some abstract brain cognitive state. Further, we maintain that it is unreasonable to expect these processes to be adequately described in the traditional way: as a linear sum of a large number of frequency generators. One of the main objectives of this thesis will be to prove that much more effective and powerful analysis of MEG can be achieved if one were to assume the presence of both linear and nonlinear characteristics from the outset. Our position is that the combined action of a relatively small number of these generators, coupled with external and dynamic noise sources, is more than sufficient to account for the complexity observed in the MEG recordings. Another problem that has plagued MEG researchers is the extremely low signal to noise ratios that are obtained. As the magnetic flux variations resulting from actual cortical processes can be extremely minute, the measuring devices used in MEG are, necessarily, extremely sensitive. The unfortunate side-effect of this is that even commonplace phenomena such as the earth's geomagnetic field can easily swamp signals of interest. This problem is commonly addressed by averaging over a large number of recordings. However, this has a number of notable drawbacks. In particular, it is difficult to synchronise high frequency activity which might be of interest, and often these signals will be cancelled out by the averaging process. Other problems that have been encountered are high costs and low portability of state-of-the- art multichannel machines. The result of this is that the use of MEG has, hitherto, been restricted to large institutions which are able to afford the high costs associated with the procurement and maintenance of these machines. In this project, we seek to address these issues by working almost exclusively with single channel, unaveraged MEG data. We demonstrate the applicability of a variety of methods originating from the fields of signal processing, dynamical systems, information theory and neural networks, to the analysis of MEG data. It is noteworthy that while modern signal processing tools such as independent component analysis, topographic maps and latent variable modelling have enjoyed extensive success in a variety of research areas from financial time series modelling to the analysis of sun spot activity, their use in MEG analysis has thus far been extremely limited. It is hoped that this work will help to remedy this oversight.

The diagenesis and palaeomagnetism of permian and triassic sediments from central Spain

Detailed diagenetic and palaeomagnetic studies have been made of Permian and Triassic rocks from the Iberian Cordillera, Spain. Five stratigraphical units comprising the Autunian, Saxonian, Buntsandstein, Muschelkalk, and Keuper have been studied in a number of sections which have been well documented sedimentologically. Autuninan rocks have a characteristic remanence which is exclusively reversed and corresponds to the Kiaman Interval. The pole position identified is consistent with previous studies, which indicate the rotation of Iberia in post-Triassic times. The Saxonian facies show complex multicomponent magnetizations; no polarity zonation can be resolved. The Buntsandstein is remarkable in that much of it is remagnetised in a direction similar to the present-day local geomagnetic field direction. The secondary remanence is carried by fine-grained haematite which has been formed as a result of carbonate dissolution associated with structural inversion of the Iberian Cordillera. Dating of diagenetic events associated with this remagnetization is also possible. Fragments of primary remanence are preserved in some fine-grained lithologies of the Buntsandstein. These indicate that normal and reversed zones of magnetization were originally present. The magnetization of the Muschelkalk and Keuper carbonates is also complex; secondary components similar to those of the Buntsandstein are present but there is evidence that primary components were predominantly normal during Karnian times.

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.