Are there connections between the Earth's magnetic field and climate?

Understanding climate change is an active topic of research. Much of the observed increase in global surface temperature over the past 150 years occurred prior to the 1940s and after the 1980s. The main causes invoked are solar variability, changes in atmospheric greenhouse gas content or sulfur due to natural or anthropogenic action, or internal variability of the coupled ocean–atmosphere system. Magnetism has seldom been invoked, and evidence for connections between climate and magnetic field variations have received little attention. We review evidence for correlations which could suggest such (causal or non-causal) connections at various time scales (recent secular variation ∼ 10–100 yr, historical and archeomagnetic change ∼ 100–5000 yr, and excursions and reversals ∼ 10³–10⁶ yr), and attempt to suggest mechanisms. Evidence for correlations, which invoke Milankovic forcing in the core, either directly or through changes in ice distribution and moments of inertia of the Earth, is still tenuous. Correlation between decadal changes in amplitude of geomagnetic variations of external origin, solar irradiance and global temperature is stronger. It suggests that solar irradiance could have been a major forcing function of climate until the mid-1980s, when “anomalous” warming becomes apparent. The most intriguing feature may be the recently proposed archeomagnetic jerks, i.e. fairly abrupt (∼ 100 yr long) geomagnetic field variations found at irregular intervals over the past few millennia, using the archeological record from Europe to the Middle East. These seem to correlate with significant climatic events in the eastern North Atlantic region. A proposed mechanism involves variations in the geometry of the geomagnetic field (f.i. tilt of the dipole to lower latitudes), resulting in enhanced cosmic-ray induced nucleation of clouds. No forcing factor, be it changes in CO₂ concentration in the atmosphere or changes in cosmic ray flux modulated by solar activity and geomagnetism, or possibly other factors, can at present be neglected or shown to be the overwhelming single driver of climate change in past centuries. Intensive data acquisition is required to further probe indications that the Earth's and Sun's magnetic fields may have significant bearing on climate change at certain time scales.     

Trends in the F2 ionospheric layer due to long-term variations in the Earth's magnetic field

The Earth's magnetic field presents long-term variations with changes in strength and orientation. Particularly, changes in the dip angle (I) and, consequently, in the sin(I)cos(I) factor, affect the thermospheric neutral winds that move the conducting plasma of the ionosphere. In this way, a lowering or lifting of the F2-peak (hmF2) is induced together with changes in foF2, depending on season, time and location. A simple theoretical approximation, developed in a previous work, is extended to a worldwide latitude–longitude grid to assess hmF2 and foF2 trends due to Earth's magnetic field secular variations. Compared to the greenhouse gases effects over the ionosphere, the Earth's magnetic field may be able to produce stronger trends which vary with season, time and location. However, to elucidate the origin of F2-region trends, long-term variations in the three possible known mechanisms should be considered altogether—greenhouse gases, geomagnetic activity and Earth's magnetic field.     

The palaeogeomagnetic field strength,variations in reversal frequency,and geomagnetic dynamo models

Abstract

The geomagnetic field and its frequent polarity reversals are generally attributed to magnetohydrodynamic (MHD) processes in the Earth's metallic and fluid core. But it is difficult to identify convincingly any MHD timescales with that over which the reversals occur. Moreover, the geological record indicates that the intervals between the consecutive reversals have varied widely. In addition, there have been superchrons when the reversals have been frequent, and at least two, and perhaps three, 35-70 Myr long superchrons when they were almost totally absent. The evaluation of these long-term variations in the palaeogeophysical record can provide crucial constraints on theories of geomagnetism, but it has generally been limited to only the directional or polarity data. It is shown here that the correlation of the palaeogeomagnetic field strength with the field's protracted stability during a fixed polarity superchron provides such a constraint. In terms of a strong field dynamo model it leads to the speculation that the magnetic Reynolds number, and the toroidal field, increase substantially during a superchron of frequent reversals.     

Agreement between Earth's rotation and mass displacement as detected by GRACE

The progress in GRACE data processing should improve the estimation of low degree spherical harmonics which are expected to agree better with Earth's rotation observations. The polar motion and length-of-day excitations determined from the spherical harmonics of the GRGS latest release (RL02) are explored and compared to the previous release (RL01). The RL02 gives best fit of the observed annual variations than RL01 and geophysical models do. However, the observed residual signal obtained after removing annual and semiannual oscillations is still better explained by the geophysical models even if RL02 estimates are improved at these frequencies scales. Linear trends are also estimated over study period (2003–2008). The linear trends of χ₁ based on GRACE RL01/02 and EOP are similar but they are very different for χ₂. Further studies with longer time series of GRACE and future gravimetric missions could help better interpret the long term variations and the effects of ice sheet mass loss or post glacial rebound. Concerning LOD variations, GRACE/LAGEOS mass displacement information brings better agreement with EOP observations, compared to the pressure term estimated by models, however the RL02 has not shown significant improvement.     

Scaling properties of the runoff variations in the arid and semi-arid regions of China: a case study of the Yellow River basin

We analyzed long daily runoff series at six hydrological stations located along the mainstem Yellow River basin by using power spectra analysis and multifractal detrended fluctuation analysis (MF-DFA) technique with aim to deeply understand the scaling properties of the hydrological series in the Yellow River basin. Research results indicate that: (1) the runoff fluctuations of the Yellow River basin exhibit self-affine fractal behavior and different memory properties at different time scales. Different crossover frequency (1/f) indicates that lower crossover frequency usually corresponds to larger basin area, and vice versa, showing the influences of river size on higher frequency of runoff variations. This may be due to considerable regulations of river channel on the runoff variations in river basin of larger basin size; (2) the runoff fluctuations in the Yellow River basin exhibit short-term memory properties at smaller time scales. Crossover analysis by MF-DFA indicates unchanged annual cycle within the runoff variations, implying dominant influences of climatic changes on changes of runoff amount at longer time scales, e.g. 1 year. Human activities, such as human withdrawal of freshwater and construction of water reservoirs, in different reaches of the Yellow River basin may be responsible for different scaling properties of runoff variations in the Yellow River basin. The results of this study will be helpful for hydrological modeling in different time scales and also for water resource management in the arid and semi-arid regions of China.     

A Tikhonov regularization method to estimate Earth's oblateness variations from global GPS observations

Earth's oblateness is varying due to the redistribution of Earth's fluid mass and the interaction of various components in the Earth system. Nowadays, continuous Global Positioning System (GPS) observations can estimate Earth's oblateness (J₂) variations with the least squares method, but are subject to ill-conditioned equations with limited GPS observations and aliasing errors from truncated degrees. In this paper, a Tikhonov regularization method is used to estimate J₂ variations from global continuous GPS observations. Results show that the J₂ has been better estimated from GPS observations based on a Tikhonov regularization method than the usual least squares method when compared to SLR solutions. Furthermore, the amplitudes and phases of the annual and semi-annual J₂ variations are closer to the SLR results with truncated degrees from 2 to 5. Higher truncated degrees will degrade the J₂ estimate. Annual J₂ variations are best estimated from GPS observations with truncated degree 4 and semi-annual J₂ variations are best estimated with truncated degree 2.     

Low frequency 1/f-like fluctuations of the AE-index as a possible manifestation of self-organized criticality in the magnetosphere

Low frequency stochastic variations of the geomagnetic AE-index characterized by 1/f^b-like power spectrum (where f is a frequency) are studied. Based on the analysis of experimental data we show that the B_z-component of IMF, velocity of solar wind plasma, and the coupling function of Akasofu are insufficient factors to explain these behaviors of the AE-index together with the 1/f^b fluctuations of geomagnetic intensity. The effect of self-organized criticality (SOC) is proposed as an internal mechanism to generate 1/f^b fluctuations in the magnetosphere. It is suggested that localized spatially current instabilities, developing in the magnetospheric tail at the initial substorm phase can be considered as SOC avalanches or dynamic clusters, superposition of which leads to the 1/f^b fluctuations of macroscopic characteristics in the system. Using the sandpile model of SOC, we undertake numerical modeling of space-localized and global disturbances of magnetospheric current layer. Qualitative conformity between the disturbed dynamics of self-organized critical state of the model and the main phases of real magnetospheric substorm development is demonstrated. It is also shown that power spectrum of sandpile model fluctuations controlled by real solar wind parameters reproduces all distinctive spectral features of the AE fluctuations.     

Impacts of riparian zone plant water use on temporal scaling of groundwater systems

In this work, we study groundwater system temporal scaling in relation to plant water use and near‐river‐stage fluctuations in riparian zones where phreatophytes exist. Using detrended fluctuation analysis (DFA), we investigate the influence of regular diurnal fluctuations due to phreatophyte water use on temporal scaling properties of groundwater level variations. We found that groundwater use by phreatophytes, at the field site on the Colorado River, USA, results in distinctive crossovers (slope changes when the plots are fitted with straight lines) in the logarithm plots of root‐mean‐square fluctuations of the detrended water level time series versus time scales of groundwater level dynamics. For groundwater levels monitored at wells close to the river, we identified one crossover at ～1 day in the scaling characteristics of groundwater level variations. When time scale exceeds 1 day, the scaling properties decrease from persistent to close to 1/f noise, where f is the frequency. For groundwater levels recorded at wells further away from the river, the slope of the straight line fit (i.e. scaling exponent) is smallest when the time scale is between 1 and 3 days. When the time scale is < 1 day, groundwater variations become persistent. When the time scale is between 1 and 3 days, the variations are close to white noise, but return to persistent when the time scale is > 3 days. Copyright © 2011 John Wiley & Sons, Ltd.     

A Quaternary geomagnetic instability time scale

Reversals and excursions of Earth's geomagnetic field create marker horizons that are readily detected in sedimentary and volcanic rocks worldwide. An accurate and precise chronology of these geomagnetic field instabilities is fundamental to understanding several aspects of Quaternary climate, dynamo processes, and surface processes. For example, stratigraphic correlation between marine sediment and polar ice records of climate change across the cryospheres benefits from a highly resolved record of reversals and excursions. The temporal patterns of dynamo behavior may reflect physical interactions between the molten outer core and the solid inner core or lowermost mantle. These interactions may control reversal frequency and shape the weak magnetic fields that arise during successive dynamo instabilities. Moreover, weakening of the axial dipole during reversals and excursions enhances the production of cosmogenic isotopes that are used in sediment and ice core stratigraphy and surface exposure dating. The Geomagnetic Instability Time Scale (GITS) is based on the direct dating of transitional polarity states in lava flows using the ⁴⁰Ar/³⁹Ar method, in parallel with astrochronologic age models of marine sediments in which oxygen isotope and magnetic records have been obtained. A review of data from Quaternary lava flows and sediments gives rise to a GITS that comprises 10 polarity reversals and 27 excursions that occurred during the past 2.6 million years. Nine of the ten reversals bounding chrons and subchrons are associated with ⁴⁰Ar/³⁹Ar ages of transitionally-magnetized lava flows. The tenth, the Gauss-Matuyama chron boundary, is tightly bracketed by ⁴⁰Ar/³⁹Ar dated ash deposits. Of the 27 well-documented geomagnetic field instabilities manifest as short-lived excursions, 14 occurred during the Matuyama chron and 13 during the Brunhes chron. Nineteen excursions have been dated directly using the ⁴⁰Ar/³⁹Ar method on transitionally-magnetized volcanic rocks and these form the backbone of the GITS. Excursions are clearly not the rare phenomena once thought. Rather, during the Quaternary period, they occur nearly three times as often as full polarity reversals.     

Geomagnetic polarity reversals in a turbulent core

The behavior of the main magnetic field components during a polarity transition is investigated using the α²-dynamo model for magnetic field generation in a turbulent core. It is shown that rapid reversals of the dipole field occur when the helicity, a measure of correlation between turbulent velocity and vorticity, changes sign. Two classes of polarity transitions are possible. Within the first class, termed component reversals, the dipole field reverses but the toroidal field does not. Within the second class, termed full reversals, both dipole and toroidal fields reverse. Component reversals result from long term fluctuations in core helicity; full reversals result from short term fluctuations. A set of time-evolution equations are derived which govern the dipole field behavior during an idealized transition. Solutions to these equations exhibit transitions in which the dipole remains axial while its intensity decays rapidly toward zero, and is regenerated with reversed polarity. Assuming an electrical conductivity of 3 × 10⁵ mho m^?1 for the fluid core, the time interval required to complete the reversal process can be as short as 7500 years. This time scale is consistent with paleomagnetic observations of the duration of reversals. A possible explanation of the cause of reversals is proposed, in which the core's net helicity fluctuates in response to fluctuations in the level of turbulence produced by two competing energy sources—thermal convection and segregation of the inner core. Symmetry considerations indicate that, in each hemisphere, helicity generated by heat loss at the core-mantle boundary may have the opposite sign of helicity generated by energy release at the inner core boundary. Random variations in rates of energy release can cause the net helicity and the α-effect to change sign occasionally, provoking a field reversal. In this model, energy release by inner core formation tends to destabilize stationary dynamo action, causing polarity reversals.     

Localized changes in geomagnetic total intensity values prior to the 1995 Hyogo-ken Nanbu (Kobe) earthquake

Changes in total geomagnetic field intensity, of 2–3 nT, were observed prior to the 1995 Hyogo-ken Nanbu (Kobe) earthquake at the Amagase (AMG) site, located approximately 70 km from the epicenter. We examined whether the observed variations are local signals arising from the Earth's crust, or global variations that are unlikely to originate from the crust. To remove global-scale variations in total geomagnetic intensity data, we employed a regional geomagnetic field model. Using data recorded at five reference sites in Japan, we estimated global-scale variations in total geomagnetic intensity, and removed them from the observed total geomagnetic intensity at the AMG site. The reminder still showed variations during the period prior to the Kobe earthquake. In addition, these pre-seismic variations include two of the largest shifts recorded during the entire observation period at the AMG site, raising the possibility that these variations were indeed related to the earthquake. These variations cannot be interpreted as signals arising from the area close to the seismic source because of the large distance between the epicenter and the site. Therefore, our results raise the possibility that the physical state of the Earth's crust shows marked changes over a wide region in the lead-up to a seismic event.     

Variations in planetary convection via the effect of climate on damage

A new model for the generation of plate tectonics suggests an important interaction between a planet's climate and its lithospheric damage behavior; and thus provides a simple explanation for the tectonic difference between Earth and Venus. We propose that high surface temperatures will lead to higher healing rates (e.g. grain growth) in the lithosphere that will act to suppress localization, plate boundary formation, and subduction. This leads to episodic or stagnant lid convection on Venus because of its hotter climate. In contrast, Earth's cooler climate promotes damage and plate boundary formation. The damage rheology presented in this paper attempts to describe the evolution of grain size by allowing for grain reduction via deformational work input and grain growth via surface tension-driven coarsening. We explore the interaction of damage and healing in two-dimensional numerical convection simulations. We also develop a simple “drip-instability” model to test the hypothesis that the competition between damage and healing controls convective and plate tectonic style by modulating episodicity at subduction zones. At small values of damage, f_A, (or large values of healing, k_A) the lithosphere remains strong enough to resist subduction on time scales of billions of years. At intermediate values of f_A and k_A the lithosphere may become mobilized and allow for short bursts of tectonic behavior followed by periods of quiescence. At large (small) values of f_A (k_A ) the fineness is increased so that the viscosity of the plate boundary is reduced to allow for continuous, unimpeded subduction of lithosphere and plate-like deformation. The results suggest the feasibility of our proposed hypothesis that the interplay of climate and damage control the mode of tectonics on a planet.     

On the validity of a model for the reversals of the Earth's magnetic field

Abstract

We have used techniques of nonlinear dynamics to compare a special model for the reversals of the Earth's magnetic field with the observational data. Although this model is rather simple, there is no essential difference to the data by means of well-known characteristics, such as correlation function and probability distribution. Applying methods of symbolic dynamics we have found that the considered model is not able to describe the dynamical properties of the observed process. These significant differences are expressed by algorithmic complexity and Renyi information.     

The geomagnetic field at the Proterozoic-Paleozoic boundary and lower-mantle plumes

Data on the amplitude of variations in the direction of the geomagnetic field and the frequency of reversals in the Vendian-Cambrian are presented. It has been established from these data that (a) distributions of variations in the direction of the geomagnetic field S _p are bimodal (modes 9° and 11°); (b) the maximum of the average amplitude S _p takes place by 5–10 Myr later than the Vendian-Cambrian boundary; (c) S _p tends to increase as plume epicenters are approached; and (d) the plume formation is more often confined to intervals with different frequencies of geomagnetic reversals than to the interval of a stable state of the geomagnetic field without reversals (Vendian hyperchron). The listed features of the geomagnetic field behavior are repeated near all boundaries of geological eras of the Phanerozoic.     

Temporal variations and scaling of streamflow and baseflow and their nitrate-nitrogen concentrations and loads

The patterns of temporal variations of precipitation (P), streamflow (SF) and baseflow (BF) as well as their nitrate-nitrogen (nitrate) concentrations (C) and loads (L) from a long-term record (28 years) in the Raccoon River, Iowa, were analyzed using variogram and spectral analyses. The daily P is random but scaling may exist in the daily SF and BF with a possible break point in the scaling at about 18 days and 45 days, respectively. The nitrate concentrations and loads are shown to have a half-year cycle while daily P, SF, and BF have a one-year cycle. Furthermore, there may be a low-frequency cycle of 6–8 years in C. The power spectra of C and L in both SF and BF exhibit fractal 1/f scaling with two characteristic frequencies of half-year and one-year, and are fitted well with the spectrum of the gamma distribution. The nitrate input to SF and BF at the Raccoon watershed seems likely to be a white noise process superimposed on another process with a half-year and one-year cycle.     

Tracking the lithium isotopic evolution of the mantle using carbonatites

Carbonatites are mantle-derived, intraplate magmas that provide a means of documenting isotopic variations of the Earth's mantle through time. To investigate the secular Li isotopic evolution of the mantle and to test whether Li isotopes document systematic recycling of material processed at or near the Earth's surface into the mantle, we analyzed the Li isotopic compositions of carbonatites and spatially associated mafic silicate rocks. The Li isotopic compositions of Archean (2.7 Ga) to Recent carbonatites (δ⁷Li = 4.1 ± 1.3 (n = 23, 1σ)) overlap the range typical for modern mantle-derived rocks, and do not change with time, despite ongoing crustal recycling. Thus, the average Li isotopic composition of recycled crustal components has not deviated greatly from the mantle value (~ + 4) and/or Li diffusion is sufficiently fast to attenuate significant heterogeneities over timescales of 10⁸ years. Modeling of Li diffusion at mantle temperatures suggests that limited δ⁷Li variation in the mantle through time reflects the more effective homogenization of Li in the mantle compared to radiogenic isotope systems. The real (but limited) variations in δ⁷Li that exist in modern mantle-derived magmas as well as carbonatites studied here may reflect isotopic fractionation associated with shallow-level processes, such as crustal assimilation and diffusive isotopic fractionation in magmatic systems, with some of the scatter possibly related to low-temperature alteration.     

Relativistic electron enhancements,magnetic storms,and substorm activity

The relativistic electron fluxes of the Earth's outer radiation belt are subjected to strong temporal variations. The most prominent changes are initiated by fast solar wind streams impinging upon the magnetosphere, which often also cause enhanced substorm activity and magnetic storms. Using 4 years of data from the particle detector REM aboard the UK satellite Strv-1b in a GTO, we investigated the relation between these different appearances of geomagnetic activity. A typical sequence is that there is a drop in the relativistic electron intensity during the main phase of the magnetic storm and a successive enhancement during the recovery phase which sometimes leads to much higher than pre-storm fluxes. Whereas the flux drop is well correlated with the magnetic storm intensity and is mainly due to the deceleration and loss of particles caused by the ring-current-induced magnetic field changes, there is only a bad correlation between the post-storm electron flux and Dst. As we show here, it is much more the level of substorm activity during the whole event which determines the size of the flux enhancements.     

Long-range dependence in the Cenozoic reversal record

The frequency distribution of intervals between Cenozoic geomagnetic reversals approximates a power law, while their occurrence over time shows temporal clustering of short and long intervals. Application of the Aggregate Variance and Absolute Value methods suggest long-range dependence in this time series, a possible indication that the geodynamo operates as a self-organised complex system. This hypothesis may allow the Cretaceous superchron to be considered as an integral part of the ordinary reversal regime attested in the Cenozoic record.     

Intensity of the geomagnetic field from recent Italian lavas using a new paleointensity method

A new method has been tested on Etna historic lavas for determining the geomagnetic field intensity (F) using the thermoremanent magnetization of volcanic rocks. The procedure involves a number of very short duration heatings above the Curie point, to produce successive laboratory TRM. Thus, it is possible to check the variations in the TRM-acquiring capacity of the samples with the time of heating (t). The curve J = f(t) is then extrapolated towards t = 0, leading to a virtual value of TRM without any laboratory heating, i.e., without the changes that currently occur when the lavas are heated to produce the TRM. Using such virtual values of TRM, satisfactory results of F had been derived from the majority of the samples studied. These results are consistent with Thellier's archeomagnetic data: they show a nearly constant intensity of the geomagnetic field during the last three centuries and, further back into the past, a significant increase of this.