Suggested Reading
There is a lively field within rock magnetism that attempts to exploit the dependence of rock magnetic parameters on concentration, grain size and mineralogy for the purpose of gleaning information about past (and present) environments. Applications in applied rock magnetism (environmental magnetism) run from detection of industrial pollution to characterizing climatic change across major climatic events to constraining rainfall variations across Asia during the Quaternary. In this lecture we will review the basic tool-kit used by environmental magnetists and illustrate various applications with examples.
There are four basic methodologies involved in most applied rock magnetism studies: imaging of magnetic separates, hysteresis parameter estimation, thermomagnetic measurements (including Curie Temperature determination and low temperature measurements) and anhysteretic remanence (ARM) measurements. Imaging is done using optical, scanning electron and transmission electron microscopes (see e.g., Figure 8.1a) on magnetic separates, or thin sections. Hysteresis measurements (including magnetic susceptibility) are made on vibrating sample magnetometers (VSMs), alternating gradient force magnetometers (AGFMs, see Lecture 7), and susceptibility meters (Figure 8.1b) of various sorts. These measurements can be done as a function of frequency or temperature. Thermomagnetic measurements are made on a “Curie Balance” (Figure 8.1c) which measures saturation magnetization as a function of temperature. ARMs are measured using an instrument that applies a large, alternating field (an AF demagntetizer) in the presence of a small DC bias field (see Lecture 5).
|
|
Images of magnetic phases are used to constrain the origin of the magnetic phases. Igneous (Figure 8.2a), detrital or aeolian (Figure 8.2b), authigenic (Figure 8.2c), biogenic (Lecture 6), anthropogenic (Figure 8.2d) and cosmic (Figure 8.2e) sources all have distinctive ear-marks, so actually looking at the particles in question can provide invaluable information.
Hysteresis behavior is strongly controlled by mineralogy and grain size; hence hysteresis loops have the potential to help constrain the makeup of a given rock specimen. The hysteresis loop of a given sample will be the sum of all the curves generated by the individual grains. Each population of grains with a consistent coercivity spectrum will leave its imprint on the resulting loop.
We have already encountered hysteresis loops in Lecture 7) and many of the associated parameters that characterize them. There are a few more, however, that are useful in environmental magnetism (see Figure 8.3).
|
|
The slope relating magnetization and applied (low) fields is called the initial magnetic susceptibility (
i) (see Lecture 1) and Lecture 3). This is a reversible measurement if the field is low enough because the magnetization will return to its initial state when the field is turned off. As the applied magnetic field increases, individual particles will reach their flipping fields, or undergo some other irreversible reorganization of spin states (rearranging domain walls, collapsing vortex structures, etc.).
Because the response to an external field is greatly enhance if a particle is superparamagnetic, SP grains are hugely more susceptible than an equivalently sized SD grain. The definition of whether a given grain is SP or not depends on the time scale of observation, so a grain can behave in a superparamagnetic manner over a long period (and come into equilibrium with the applied field) but be SD over a shorter time scale (and have only a sluggish response to a small applied field). Therefore
i is strongly frequency dependent (as well as being strongly temperature dependent). Some instruments allow the measurement of
i at various frequencies allowing the definition of the so-called frequency dependent susceptibility or
fd. This is often used to estimate the contribution of SP particles to the total susceptibility.
Saturation magnetization (Ms) is the magnetization measured in the presence of a saturating field (Bs). This measurement must often be “corrected” for the contribution of paramagnetic minerals whose high field susceptibility
hf must be subtracted. Fortunately, paramagnetic behavior is linear up to several tesla so can usually be estimated and removed. If we subtract the high field susceptibility (which is only the paramagnetic contribution)
from the initial susceptibility, we can estimate the contribution of the ferrimagnetic (sensu lato) particles or
ferri.
Susceptibility can also be measured as a function of the orientation of the specimen with respect to the applied magnetic field. If the susceptibility is independent of orientation, it is said to be isotropic. Anisotropic orientations of magnetic minerals can lead to an anisotropic magnetic susceptibility response which in turn can be interpreted in terms of preferred orientation of magnetic phases. This topic will be addressed in later lectures.
The portion of the hysteresis loop that is recorded while the field is ramping up is called the ascending loop and the return portion recorded as the field is ramping down is the descending loop. Once the field is high enough, irreversible changes in the magnetization of the sample take place and the magnetization will no longer return to its initial state after the field is switched off; it displays hysteresis (see Lecture 7). The magnetization thus acquired is an IRM (see Lecture 5). The remanence remaining after application of a saturating field was termed saturation remanence Mr in Lecture 7) (also known as Mrs or Msr in the literature). It is also synonymous with the saturation IRM (sIRM).
As mentioned in Lecture 5, the coercive field (Bc) is that field required to reduce the net magnetization to zero and the bulk coercivity of remanence (Bcr) is the field necessary to flip half the magnetic moments (so when the field is turned off, the remanence is reduced to zero). Two ways of estimating Bcr were described in Lecture 5(Bcr and B'cr). A third way is the intersection
method described in Lecture 7) (Bcr*). A fourth way is the
M method illustrated in Figure 8.3b whereby the difference between the ascending and descending loops (
M) from Figure 8.3a is plotted versus applied field. The field at which the value of
M
is 50% of the maximum is here called Bcr**.
Robertson and France (1994) suggested that if populations of magnetic materials have generally log-normally distributed coercivity spectra and if the IRM is the linear sum of all the contributing grains, then an IRM acquisition curve could be “unmixed” into the contributing components. The basic idea is illustrated in Figure 8.4 whereby two components each with log normally distributed coercivity spectra (see dashed and dashed-dotted lines in the inset) create the IRM acquisition curve shown. Thus by obtaining a very well determined IRM acquisition plot (the “linear acquisition plot” or LAP in Figure 8.4 using the terminology of Kruiver et al., 2001), one could first differentiate it to get the “gradient acquisition plot” or GAP (heavy solid line in the inset to Figure 8.4 ). This then can be “unmixed” to get the parameters of the contributing components such as the mean and standard deviation of the log-normal distribution (called B1/2 and DP respectively by Robertson and France, 1994). Note that B1/2 is synonymous with Bcr if there is only one population of coercivities. Also, other forms of magnetic remanence (e.g., ARM), demagnetization as well as acquisition, and other distributions are also possible as are fancier methods of inversion (see e.g., Egli, 2003).
|
|
Another very useful technique for characterizing the magnetic mineralogy in a sample is the “3D IRM” technique of Lowrie (1990). Some important magnetic phases in geological materials (Table 6.1; Lecture 6) are magnetite (maximum blocking temperature of ~580oC, maximum coercivity of about 0.3 T), hematite (maximum blocking temperature of ~ 675oC and maximum coercivity larger than several tesla), goethite (maximum blocking temperature of ~ 125oC and maximum coercivity of much larger than 5 T), and various sulfides. The relative importance of these minerals in bulk samples can be constrained by a simple trick that exploits both differences in coercivity and unblocking temperature (Lowrie, 1990).
The “3D IRM” technique proceeds as follows:
• Apply an IRM along three orthogonal directions in three different fields. The first field, applied along X1, should be sufficient to saturate all the minerals within the sample and is usually the largest field achievable in the laboratory (say 2 T). The second field, applied along X2, should be sufficient to saturate magnetite, but not to realign high coercivity phases, such as goethite or fine-grained hematite (say 0.4 T). The third IRM, applied along X3, should target low coercivity minerals and the field chosen is typically something like 0.12 T.
• The composite magnetization can be characterized by determining the blocking temperature spectra for each component. This is done by heating the sample in zero field to successively larger temperatures, cooling then measuring the remaining magnetization. The magnitude of the three cartesian components (x1,x2,x3) of the remaining remanence is then plotted versus demagnetizing temperature.
An example of 3D IRM data are shown in Figure 8.5. The curve is dominated by a mineral with a maximum blocking temperature of between 550o and 600oC and has a coercivity less than 0.4 T, but greater than 0.12 T. These properties are typical of magnetite (Table 6.1; Lecture 6). There is a small fraction of a high coercivity (>0.4 T) mineral with a maximum unblocking temperature > 650oC, which is consistent with the presence of hematite (Table 6.1; Lecture 6).
|
|
|
|
|
|
IRM and ARM acquisition and demagnetization curves could be a rich source of information about the magnetic phases in rocks. However, these are extremely time consuming to measure taking hours for each one. Hysteresis loops on the other hand are quick, taking about 10 minutes to measure the outer loop. In principle, the same information could be had from hysteresis loops as in the IRM acquisition curves.
Hysteresis loops, like IRM acquisition curves are the sum of all the contributing particles in the sample. There are several basic types of loops which are recognized the “building blocks” of the hysteresis loops we measure on geological materials. We illustrate some of the building blocks of possible hysteresis loops in Figure 8.6. Figure 8.6a shows the negative slope typical of diamagnetic material such as carbonate or quartz, while Figure 8.6b shows a paramagnetic slope. Such slopes are common when the sample has little ferromagnetic material and is rich in iron-bearing phases such as biotite or clay minerals.
When grain sizes are very small, a sample can display superparamagnetic “hysteresis” behavior (Figure 8.6c). The SP curve follows a Langevin function L(
) (see Lecture 7) where
is MsvB/kT, but integrates over the distribution of v in the sample.
Above some critical volume, grains will have relaxation times that are sufficient to retain a stable remanence (Lecture 5). As discussed in Lecture 7), populations of randomly oriented stable grains can produce hysteresis loops with a variety of shapes, depending on the origin of magnetic anisotropy and domain state. We show loops from samples that illustrate representative styles of hysteresis behavior in Figure 8.6d-f. Figure 8.6d shows a loop characteristic of samples whose remanence stems from SD magnetite with uniaxial anisotropy. In Figure 8.6e, we show data from specular hematite whose anisotropy is magnetocrystalline in origin (hexagonal within the basal plane). Note the very high Mr/Ms ratio of nearly one. Finally, we show a loop that has lower Mr/Ms ratios than single domain, yet some stability. Loops of this type have been characterized as pseudo-single domain PSD (Figure 8.6f). We now know that PSD behavior is typical of vortex remanence state particles.
|
|
In the messy reality of geological materials, we often encounter mixtures of several magnetic phases and/or domain states. Such mixtures can lead to distorted loops, such as those shown in Figure 8.7. In Figure 8.7a, we show a mixture of hematite plus SD-magnetite. The loop is distorted in a manner that we refer to as goose-necked. Another commonly observed mixture is SD plus SP magnetite which can result in loops that are either wasp-waisted (see Figure 8.7b) or pot-bellied (see Figure 8.7c).
Considering the loops shown in Figure 8.7, we immediately notice that there are two distinct causes of loop distortion: mixing two phases with different coercivities and mixing SD and SP domain states. We differentiate the two types of distortion as “goose-necked” and “wasp-waisted” (see Figure 8.7) because they look different and they mean different things.
Jackson et al. (1990) suggested that the
M curve (see Figure 8.3) could be differentiated to reveal different coercivity spectra contained in the hysteresis loop. The
M curve and its derivitive (d
M/dH) are sensitive only to the remanence carrying phases, and not, for example, to the SP fraction, we can use these curves to distinguish the two sources of distortion.
In Figure 8.8, we show several representative loops, along with the
M and d
M/dH curves. Distortion resulting from two phases with different coercivities (e.g., hematite plus magnetite or two distinct grain sizes of the same mineral) results in a “two humped” d
M/dH curve, whereas wasp-waisting which results from mixtures of SD + SP populations have only one “hump”.
Jackson et al. (1990) also suggested a way to deal with noisy data using Fourier smoothing. This treatment is described in the appendix.
One holy grail of applied rock magnetism is a diagnostic set of measurements that will yield unambiguous grain size information. To this end, large amounts of hysteresis data have been collected on a variety of minerals that have been graded according to size and mode of formation. The most complete set of data are available for magnetite, as this is the most abundant magnetic phase in the world. There are three sources for magnetite typically used in these experiments: natural crystals that have been crushed and seived into grain size populatins, crystals that were grown by a glass ceramic technique and crystals grown from hydrothermal solution. In Figure 8.9a-c we show a compilation of grain size dependence of coercive force, remanence ratio, coercivity of remanence respectively. There is a profound dependence not only on grain size, but on mode of formation as well. Crushed particles tend to have much higher coercivities and remanence ratios than grown crystals, presumably because of the increased dislocation density which stabilizes domain walls in much the same manner as do voids. These abnormally high values disappear to a large extent when the particles are annealed at high temperature - a procedure which allows the dislocations to “relax” away.
The behavior of initial magnetic susceptibility is shown in Figure 8.9d. There is no strong trend with grain size over the entire range of grain sizes from single domain to multi-domain magnetite. Susceptibility is predicted to be sensitive to the SD/SP domain state transition, however, because in SP particles, the magnetization is unconstrained by magnetocrystalline or shape anisotropy energies, hence has a larger response to an applied field by a factor of ln(C
). Taking C to be 1010s-1 and
to be order 100 s, we find a factor of ~ 28 enhancement of magnetic susceptibility for SP grains over an SD grain of the same volume. This is the basis for using frequency dependence to detect the contribution of SP grains to a population. Because SP behavior depends on the time scale of observation, particles may behave SP at lower frequencies and not at higher frequencies. However, because
is exponentially sensitive to temperature
(T) often
yields much more information than
(f).
Grain size trends in ARM are shown in Figure 8.9e. This trend is very poorly constrained because ARM is also a strong function of concentration and the method by which the particles were prepared. Some (e.g. Banerjee et al., 1981) have suggested that the ratio of ARM (normalized by the DC field applied to get the so-called ARM susceptibility or
ARM to
could be used to determine grain size in magnetite, but there are substantial practicle difficulties with this method, unless a great deal is known about concentration,
origin of the minerals and magnetic mineralogy.
|
|
There is a bewildering array of parameter ratios that are in popular use in the applied rock and mineral magnetism literature: Mr/Ms,Bcr/Bc,ARM/
,ARM/Mr,Mr/
,IRM(x)/Mr.
The ratios Mr/Ms and Bcr/Bc are sensitive to remanence state (SP, SD, flower, vortex, MD) and the source of magetic anisotropy (cubic, uniaxial, defects), hence reveal something about grain size and shape. For this reason Day et al., (1977) began plotting these ratios on a diagram known as the “Day” plot (see e.g., Figure 8.10).
Day plots are divided into regions of nominally SD, “PSD” and MD behavior using some theoretical bounds as guides (see Lecture 7). The designation PSD stands for pseudo-single domain and has Mr/Ms ratios in between those characteristic of SD behavior (0.5 or higher) and MD (.05 or lower). In practice nearly all geological materials plot in the PSD box so the usefulness of the Day plot is limited.
|
|
|
|
Of course Mr/Ms could also be plotted against Bc (e.g, Neél, 1955). This type of plot has been characterized to some extent using micromagnetic modeling techniques (e.g., Tauxe et al., 2002) to aide in the interpretation of hysteresis data in terms of domain state and origin of magnetic anisotropy energy.
Ratios involving ARM and
can be complicated because both of these parameters have ratio complicated behaviors themselves. ARM is a strong function of concentration and not monotonic with grain size.
is a “garbage can” parameter that reflects everything in the sample to some extent. Under certain uncomplicated conditions, both of these parameters can be quite useful, but care should be exercised in their interpretation.
Ratios of a lower back-field IRM (IRM(-x) to saturation IRM (Mr or sIRM), the so-called S-ratio, can be used to quantify the ratio of hard (magnetized at saturation) to soft (remagnetized in the back field direction) minerals in a sample. IRM is less affected by particle interactions so behaves more linearly with concentration.
Rock magnetic parameters are relatively quick and easy to measure, compared to geochemical, sedimentological and paleontological data. When used judiciously, they can be enormously helpful in constraining a wide variety of climatic and environmental changes. There are two basic types of plots of the rock and mineral magnetic parameters discussed in this lecture: bi-plots and depth plots. Bi-plots, for example ARM versus
have been in use since Banerjee et al. (1981) (see e.g., Figure 8.12). Biplots can be useful for detecting changes in grain size, concentration, mineralogy, etc. If for example, the data in a plot of Mr
verus
plot on a line, it may be appropriate to interpret the dominant control on the rock magnetic parameters as changes in concentration alone.
Depth plots are useful for core correlation, variations in concentration, mineralogy and grain size as a function of depth. An example of a paleoceanographic application of rock magnetism from Hartl et al. (1995) is illustrated in Figures 8.13 and 8.14. Trends in isotopes and carbonate are shown in Figure 8.13. These indicate a major change in the environment at the end of the Eocene.
|
|
IRM and
were measured as well. The pronounced changes in carbon isotopes appeared to be mirrored in the rock magnetic variations. It was tempting to attribute these variations to be simply related to the complimentary changes in % CaCO3 because both IRM and
are approximately linear functions of concentration. However, a look at the bi-plot of IRM versus
(Figure 8.14) shows
a more complicated and interesting picture. If the variations in these parameters were only caused by changes in concentration, the bi-plot would show a straight line. Instead, the data plot along two lines with the Eocene data having a different slope and an overall larger constribution to
than the Oligocene group. As mentioned earlier in the lecture,
is a very strong function of the fraction of superparamagnetic grains to the population and Hartl et al. (1995) argue that there is a shift in grain size associated with the Eocene/Oligocene boundary in this case, with a greater fraction of SP grains in the Eocene than
in the Oligocene. A change in magnetic grain size can be the result of changes in the pore water chemistry resulting from changes in organic carbon delivery. This mechanism is consistent with the carbon isotopic variations shown in Figure 8.13.
|
|
In practice, hysteresis measurements may yield rather noisy data. Jackson et al. [1990] suggested that noisy hysteresis data could be filtered using a Fourier transform. The advantages of Fourier smoothing are that the calculated hysteresis parameters are less sensitive to noise and that the
M and d
M/dH curves are more readily interpreted.
The steps involved in Fourier smoothing of hysteresis loops are as follows (see Figure A1):
• First, the contribution of paramagnetic and diamagnetic phases must be removed. Figure A1a shows some typical data from carbonate rich sediments. These samples have a strong diamagnetic (negative high field slope) contribution. We remove the diamagnetic contribution by calculating a best-fit line using linear regression for the data at high fields (after the ferromagnetic phases have reached saturation) and removing its contribution by subtraction (see Figure A1b).
• In order to ensure uniformity of data treatment, Jackson et al. [1990] recommend truncating the data at some fixed percentage of Ms (after slope adjustment). We truncate the data at 99.9% of Ms in Figure A1b.
• A Fourier transform requires data with a single y value for every x value and hysteresis data, as normally plotted are not suitable. The loops can be mapped into a suitable form for Fourier analysis by transforming the field values into radians, as shown in Figure A1c. The unfolded loop starts at the point when the descending curve intersects the y axis (Mr). From H = 0
-Hmax,
H is mapped linearly to radians (H' = 0
/2). From H = -Hmax
0, H is mapped to H' =
/2
. From H = 0
+Hmax, we map H to H' =
3
/2, and finally, for H = +Hmax
0, H is converted to H' = 3
/2
2
.
• The “unfolded” data can then be subjected to a Fourier Transform as described by Jackson et al. [1990]. The data can be smoothed by retaining only a specified number of terms (see Figure A1d). Finally, hysteresis parameters can be calculated from the reconstituted loop and
M and d
M/dH curves can be plotted (see Figure A1e-f).
|
|
Banerjee, S. K., King, J. & Marvin, J. (1981), ‘A rapid method for magnetic granulometry with applications to environmental studies’, Geophys. Res. Lett. 8, 333-336.
Day, R., Fuller, M. D. & Schmidt, V. A. (1977), ‘Hysteresis properties of titanomagnetites: grain size and composition dependence’, Phys. Earth Planet. Inter. 13, 260-266.
Dunlop, D. (2002), ‘Theory and application of the Day plot (Mrs/Ms versus Hcr/Hc) 2. Application to data for rocks, sediments, and soils’, J. Geophys. Res 107, doi:10.1029/2001JB000487.
Dunlop, D. & Argyle, K. (1997), ‘Thermoremanence, anhysteretic remanence and susceptibility of submicron magnetites: Nonlinear field dependence and variation with grain size’, J. Geophys. Res 102, 20199-20210.
Egli, R. (2003), ‘Analysis of the field dependence of remanent magnetization curves’, Journal of Geophysical Research-Solid Earth 108(B2).
Evans, M. & Heller, F. (2003), Environmental Magnetism: Principles and Applications of Enviromagnetics, Academic Press.
Hartl, P., Tauxe, L. & Herbert, T. (1995), ‘Eartliest Oligocene increase in South Atlantic productivity as interpreted from ”rock magnetics” at DSDP Site 522’, Paleoceanography 10, 311-325.
Heider, F., Zitzelsberger, A. & Fabian, K. (1996), ‘Magnetic susceptibility and remanent coercive force in grown magnetite crystals from 0.1
m to 6mm’, Phys. Earth Planet. Inter. 93, 239-256.
Hunt, C. P., Moskowitz, B. M. & Banerjee, S. K. (1995), ‘Rock Physics and Phase Relations, A Handbook of Physical Constants’, pp. 189-204.
Jackson, M., Worm, H. U. & Banerjee, S. K. (1990), ‘Fourier analysis of digital hysteresis data: rock magnetic applications’, Phys. Earth Planet. Inter. 65, 78-87.
Kruiver, P. P., Dekkers, M. J. & Heslop, D. (2001), ‘Quantification of magnetic coercivity components by the analysis of acquisition curves of isothermal remanent magnetisation’, Earth and Planetary Science Letters 189(3-4), 269-276.
Lowrie, W. (1990), ‘Identification of ferromagnetic minerals in a rock by coercivity and unblocking temperature properties’, Geophys. Res. Lett. 17, 159-162.
Maher, B. A. & Thompson, R., eds (1999), Quaternary Climates, Environments and Magnetism, Cambridge University Press.
Néel, L. (1955), ‘Some Theoretical Aspects of Rock-Magnetism’, Adv. Phys 4, 191-243.
Robertson, D. J. & France, D. E. (1994), ‘Discrimination of remanence-carrying minerals in mixtures, using isothermal remanent magnetisation acquisition curves’, Physics of The Earth and Planetary Interiors 82(3-4), 223-234.
Tauxe, L., Bertram, H. & Seberino, C. (2002), ‘Physical interpretation of hysteresis loops: Micromagnetic modelling of fine particle magnetite’, Geochem., Geophys., Geosyst. 3, DOI 10.1029/ 2001GC000280.
Tauxe, L., Mullender, T. A. T. & Pick, T. (1996), ‘Potbellies, wasp-waists, and superparamagnetism in magnetic hysteresis’, Jour. Geophys. Res. 101, 571-583.
|
[ Home ]
[ Databases ]
[ Events ]
[ Tools ]
[ Publications ]
[ Links ]
[ Google Site Search ]
[ Metadata ] [ Who's Who ] [ Browsers and Plugins ] [ Database Indexes ] | |
| This page was last updated on 04-Apr-2008 Sponsored by NSF EAR 0000998 |
Supported by the San Diego Supercomputer Center and the Scripps Institution of Oceanography |