(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).

Geomagnetic paleointensity dating of South China Sea sediments for the last 130 kyr

Relative paleointensity records from the northern South China Sea, northwest Pacific Ocean were studied in two gravity piston cores. Continuous mineral magnetic and paleomagnetic measurements were made using discrete sediment samples. Detailed rock magnetic parameters, such as thermomagnetic and high-field hysteresis data, indicate that pseudo-single domain magnetite in a narrow range of grain-size and concentration is the main contributor to the remanent magnetization. The uniform magnetic mineralogy meets the commonly accepted criteria for establishing relative paleointensity records. The relative paleointensity (RPI) curves were constructed by normalizing the natural remanent magnetization (NRM) with isothermal remanent magnetization (IRM), both in the 20-60 mT demagnetization state. Dating constraints have been provided by radiocarbon ages in the upper 400 cm of both cores. Furthermore, we have correlated our paleointensity records with NAPIS-75, S.Atlantic-1089, Sint-200 and NOPAPIS-250 to determine the chronological RPI framework for the South China Sea (SCS-PIS). Although some temporal offsets of paleointensity features between the different records have been recognized, their similar shape suggests that relative paleointensity on the 10(3)-10(4) year scale is globally coherent and can provide an age framework for sediments independent of delta O-18 ages.

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.

Long-term variations in the magnetic field of the Sun and possible implications for terrestrial climate

Recent studies of the variation of geomagnetic activity over the past 140 years have quantified the "coronal source" or "open" magnetic flux F-s that leaves the solar atmosphere and enters the heliosphere and have shown that it has risen, on average, by 34% since 1963 and by 140% since 1900. This variation is reflected in studies of the heliospheric field using isotopes deposited in ice sheets and meteorites by the action of galactic comic rays. The variation has also been reproduced using a model that demonstrates how the open flux accumulates and decays, depending on the rate of flux emergence in active regions and on the length of the solar cycle. The cosmic ray flux at energies > 3 GeV is found to have decayed by about 15% during the 20(th) century (and by about 4% at > 13 GeV). We show that the changes in the open flux do reflect changes in the photospheric and sub-surface field which offers an explanation of why open flux appears to be a good proxy for solar irradiance extrapolation. Correlations between F-s, solar cycle length, L, and 11-year smoothed sunspot number, R-11, explain why the various irradiance reconstructions for the last 150 years are similar in form. Possible implications of the inferred changes in cosmic ray flux and irradiance for global temperatures on Earth are discussed.

Solar causes of the long-term increase in geomagnetic activity

We analyze the causes of the century-long increase in geomagnetic activity, quantified by annual means of the aa index, using observations of interplanetary space, galactic cosmic rays, the ionosphere, and the auroral electrojet, made during the last three solar cycles. The effects of changes in ionospheric conductivity, the Earth's dipole tilt, and magnetic moment are shown to be small; only changes in near-Earth interplanetary space make a significant contribution to the long-term increase in activity. We study the effects of the interplanetary medium by applying dimensional analysis to generate the optimum solar wind-magnetosphere energy coupling function, having an unprecedentedly high correlation coefficient of 0.97. Analysis of the terms of the coupling function shows that the largest contributions to the drift in activity over solar cycles 20-22 originate from rises in the average interplanetary magnetic field (IMF) strength, solar wind concentration, and speed; average IMF orientation has grown somewhat less propitious for causing geomagnetic activity. The combination of these factors explains almost all of the 39% rise in aa observed over the last three solar cycles. Whereas the IMF strength varies approximately in phase with sunspot numbers, neither its orientation nor the solar wind density shows any coherent solar cycle variation. The solar wind speed peaks strongly in the declining phase of even-numbered cycles and can be identified as the chief cause of the phase shift between the sunspot numbers and the aa index. The rise in the IMF magnitude, the largest single contributor to the drift in geomagnetic activity, is shown to be caused by a rise in the solar coronal magnetic field, consistent with a rise in the coronal source field, modeled from photospheric observations, and an observed decay in cosmic ray fluxes.

Ionospheric and geomagnetic responses to changes in IMF B_Z: a superposed epoch study

Superposed epoch studies have been carried out in order to determine the ionospheric response at mid-latitudes to southward turnings of the interplanetary magnetic field (IMF). This is compared with the geomagnetic response, as seen in the indices K p, AE and Dst. The solar wind, IMF and geomagnetic data used were hourly averages from the years 1967–1989 and thus cover a full 22-year cycle in the solar magnetic field. These data were divided into subsets, determined by the magnitudes of the southward turnings and the concomitant increase in solar wind pressure. The superposed epoch studies were carried out using the time of the southward turning as time zero. The response of the mid-latitude ionosphere is studied by looking at the F-layer critical frequencies, f o F2, from hourly soundings by the Slough ionosonde and their deviation from the monthly median values, δf o F2. For the southward turnings with a change in B z of δB z > 11.5 nT accompanied by a solar wind dynamic pressure P exceeding 5 nPa, the F region critical frequency, f o F2, shows a marked decrease, reaching a minimum value about 20 h after the southward turning. This recovers to pre-event values over the subsequent 24 h, on average. The Dst index shows the classic storm-time decrease to about −60 nT. Four days later, the index has still to fully recover and is at about −25 nT. Both the K p and AE indices show rises before the southward turnings, when the IMF is strongly northward but the solar wind dynamic pressure is enhanced. The average AE index does register a clear isolated pulse (averaging 650 nT for 2 h, compared with a background peak level of near 450 nT at these times) showing enhanced energy deposition at high latitudes in substorms but, like K p, remains somewhat enhanced for several days, even after the average IMF has returned to zero after 1 day. This AE background decays away over several days as the Dst index recovers, indicating that there is some contamination of the currents observed at the AE stations by the continuing enhanced equatorial ring current. For data averaged over all seasons, the critical frequencies are depressed at Slough by 1.3 MHz, which is close to the lower decile of the overall distribution of δf o Fl values. Taking 30-day periods around summer and winter solstice, the largest depression is 1.6 and 1.2 MHz, respectively. This seasonal dependence is confirmed by a similar study for a Southern Hemisphere station, Argentine Island, giving peak depressions of 1.8 MHz and 0.5 MHz for summer and winter. For the subset of turnings where δB z > 11.5 nT and P ≤ 5 nPa, the response of the geomagnetic indices is similar but smaller, while the change in δf o F2 has all but disappeared. This confirms that the energy deposited at high latitudes, which leads to the geomagnetic and ionospheric disturbances following a southward turning of the IMF, increases with the energy density (dynamic pressure) of the solar wind flow. The magnitude of all responses are shown to depend on δB z . At Slough, the peak depression always occurs when Slough rotates into the noon sector. The largest ionospheric response is for southward turnings seen between 15–21 UT.

Field-aligned plasma flow in the quiet, mid-latitude ionosphere deduced from topside soundings

A method for quantifying diffusive flows of O+ ions in the topside ionosphere from satellite soundings is described. A departure from diffusive equilibrium alters the shape of the plasma scale-height profile near the F2-peak where ion-neutral frictional drag is large. The effect enables the evaluation of , the field-aligned flux of O+ ions relative to the neutral oxygen atom gas, using MSIS model values for the neutral thermospheric densities and temperature. Upward flow values are accurate to within about 10%, the largest sources of error being the MSIS prediction for the concentration of oxygen atoms and the plasma temperature gradient deduced from the sounding. Downward flux values are only determined to within 20%. From 60,000 topside soundings, taken at the minimum and rising phase of the solar cycle, a total of 1098 mean scale-height profiles are identified for which no storm sudden commencement had occurred in the previous 12 days and for which Kp was less than 2o, each mean profile being an average of about six soundings. A statistical study ofdeduced from these profiles shows the diurnal cycle of O+ flow in the quiet, topside ionosphere at mid-latitudes and its seasonal variations. The differences betweenand ion flux observations from incoherent scatter radars are considered using the meridional thermospheric winds predicted by a global, three-dimensional model. The mean interhemispheric flow from summer to winter is compared with predictions by a numerical model of the protonospheric coupling of conjugate ionospheres for up to 6 days following a geomagnetic storm. The observed mean (of order 3 × 1016 ions day−1 along a flux tube of area 1 m2 at 1000 km) is larger than predicted for day 6 and the suggested explanation is a decrease in upward flows from the winter, daytime ionosphere between the sixth and twelfth days.

Extracting planetary waves from geomagnetic time series using Empirical Mode Decomposition

Empirical Mode Decomposition is presented as an alternative to traditional analysis methods to decompose geomagnetic time series into spectral components. Important comments on the algorithm and its variations will be given. Using this technique, planetary wave modes of 5-, 10-, and 16-day mean periods can be extracted from magnetic field components of three different stations in Germany. In a second step, the amplitude modulation functions of these wave modes can be shown to contain significant contribution from solar cycle variation through correlation with smoothed sunspot numbers. Additionally, the data indicate connections with geomagnetic jerk occurrences, supported by a second set of data providing reconstructed near-Earth magnetic field for 150 years. Usually attributed to internal dynamo processes within the Earth's outer core, the question of who is impacting whom will be briefly discussed here.

Oceanic long-distance navigation : do experienced migrants use the earth's magnetic field?

Albatrosses and sea turtles are known to perform extremely long-distance journeys between disparate feeding areas and breeding sites located on small, isolated, oceanic islands or at specific coastal sites. These oceanic journeys, performed mainly over or through apparently featureless mediums, indicate impressive navigational abilities, and the sensory mechanisms used are still largely unknown. This research used three different approaches to investigate whether bi-coordinate navigation based on magnetic field gradients is likely to explain the navigational performance of wandering albatrosses in the South Atlantic and Indian Oceans and of green turtles breeding on Ascension Island in the South Atlantic Ocean. The possibility that magnetic field parameters can potentially be used in a bi-coordinate magnetic map by wandering albatrosses in their foraging area was investigated by analysing satellite telemetry data published in the literature. The possibilities for using bi-coordinate magnetic navigation varied widely between different areas of the Southern Oceans, indicating that a common mechanism, based on a bi-coordinate geomagnetic map alone, was unlikely for navigation in these areas. In the second approach, satellite telemetry was used to investigate whether Ascension Island green turtles use magnetic information for navigation during migration from their breeding island to foraging areas in Brazilian coastal waters. Disturbing magnets were applied to the heads and carapaces of the turtles, but these appeared to have little effect on their ability to navigate. The only possible effect observed was that some of the turtles with magnets attached were heading for foraging areas slightly south of the control turtles along the Brazilian coast. In the third approach, breeding female green turtles were deliberately displaced in the waters around Ascension Island to investigate which cues these turtles might use to locate and return to the island; the results suggested that cues transported by wind might be involved in the final stages of navigation.