Equatorial and mid-latitude records of the last ... - ScienceDirect

3 downloads 0 Views 866KB Size Report
Equatorial and mid-latitude records of the last geomagnetic reversal from the Atlantic Ocean. Jean-Pierre Valet 1, Lisa Tauxe 2 and Bradford Clement 3 i Centre ...
Earth and Planetary Science Letters, 94 (1989) 371-384 Elsevier Science Publishers B.V., A m s t e r d a m - Printed in The Netherlands

371

[51

Equatorial and mid-latitude records of the last geomagnetic reversal from the Atlantic Ocean J e a n - P i e r r e Valet 1, Lisa T a u x e 2 and B r a d f o r d C l e m e n t 3 i Centre des Faibles Radioactivit~s, Laboratoire mixte CNRS-CEA, Avenue de la Terrasse, 91198 Gifsur Yvette Cedex (France) 2 Scripps Institution of Oceanography, GRD, A-020, La Jolla, CA 92093 (U.S.A.) 3 Department of Geology, Florida International University, Miami, FL 33199 (U.S.A.) Received January 25, 1989; revised version accepted June 8, 1989 Three records of the Matuyama-Brunhes reversal have been obtained from O D P deep-sea cores distributed along the West African margin in the Atlantic Ocean. These studies in addition to the record from DSDP Site 609B provide a latitudinal transect extending from the equator to 50 ° N. Simulations of natural smoothing by post-depositional processes show that characteristic features of a transitional field geometry dominated by axisymmetry should be preserved in the records, especially at the equator. The results do not seem to favour the hypothesis that axisymmetrical terms would dominate during this transition.

1. Introduction Many models have been proposed in attempts to account for the non dipolar structure of the geomagnetic field during reversals. The first compilations of paleomagnetic records across reversal boundaries appeared to support models in which the transitional field was dominated by low-order zonal terms (e.g. [1]). Subsequent results [2-5], however, suggested that in most cases the transitional field cannot be purely axisymmetric and non-zonal components were thus incorporated in the models [6-8]. The most recent reversal is certainly the most appropriate for studies of the geometry of the reversing field because it is easily identified in sedimentary sequences which have undergone the least chemical alteration. The first attempts to analyze the harmonic content of the field from the existing records of MatuyamaBrunhes were mostly consistent with the presence of zonal axisymmetrical components [7,9,10]. In particular, the approach followed by Williams and Fuller [11] focussed on the relationship between the site latitude the duration and the signature of the inclination during the transition. However, most of the records used were poorly documented and the geographical distribution of the sites was not very adequate. Studies of Matuyama-Brunhes involving latitudinal transects as well as records from the equator and the southern hemisphere are 0012-821X/89/$03.50

© 1989 Elsevier Science Publishers B.V.

crucial for testing whether the transitional field is indeed dominated by axisymmetrical terms. In this paper we present the results of three records of Matuyama-Brunhes obtained from Holes 659C (18 o 04' N, 21° 01' W), 664D (0 o 06' N, 23°16'W °) and 665B drilled during ODP Leg 108. The complete results of the magnetostratigraphy [12] as well as detailed paleomagnetic studies of the first two records are reported elsewhere [13] and we include here data from Hole 665B located also at the equator (23°16'W). Thus, together with Hole 609B ( 5 0 ° N , 2 5 ° W ) drilled during DSDP Leg 94 [14] these records provide a latitudinal transect extending from the equator to 50 ° N along a mean longitude of 20 ° W.

2. Sampling and laboratory techniques The sediments from Leg 108 are composed of light grey to dark grey (659C, 664D) or brown (665B) nannofossil ooze. No apparent change in the lithology was observed through the transition zones which were identified from the results of shipboard measurements at 23.2 meters below sea floor (mbsf) (659C), 28 mbsf (664D) and 13.70 mbsf (665B). Sampling was accomplished with " U " channels [15] 50-100 cm long. The " U " channels were then sliced into 0.5-1 cm thick wafers giving two or three specimen at the same level for Holes 659C and 664D; only one speci-

372

664D-4-S-100A

664-4-5-20A

664D-4-5- 50A

UP -W

UP

W

UP

~I

- -

50mT

'

' -"-"~N

, I5, ' " T , ,3qnq./m-~'5~. , s O' ~w : ~ - I : : - : 1 ~ "

DOWN • Dec • Inc

NRN 37,S3 :0E-B A/M

659C-3-4-25

659C-3-3-130

5

w uP

W UP

~,,

665-4

NRM [ 1 3 7

10E-6

'~OmT

s

35

N

63mT .

t

l

oDec • Inc

DOWN

NRM 5 0 8

u~/~

'tOE-6

A/M

665B-20

665B-19

A/M

up/w665B-23 7

N/N

~/'N

2o _

_

N/N

Fig. 1. Demagnetization diagrams relative to samples (most of them with intermediate directions) from Holes 659C, 664D and 665B. The samples 665B-19 and 20 are stratigraphically consecutive and adjacent to the transition.

373

men per level was taken in 665B as the magnetization intensity of the sediment was lower. Coercivities of the natural remanent magnetization (NRM) lower than 60 roT, blocking temperatures of the N R M of under 580 o C and saturation of isothermal remanent magnetization (IRM) attained below 0.2 T suggest that the magnetic mineralogy is dominated by magnetite in every section. Attempts to determine the range of grain size using the method of the low-field susceptibility versus anhysteretic susceptibility (ARM) [16] are consistent with fine magnetite grains. The downcore variations in susceptibility and ARM do not show changes by more than a factor of 2 (see Fig. 3 and Valet et al. [13]). These results establish the uniformity of the magnetic properties along the transitional intervals. The measurements were done in the shielded rooms of the CFR and Scripps paleomagnetic laboratories using LETI and CTF three-axis cryogenic magnetometers. A total of 273 samples including all those located near or within the transitions were stepwise demagnetized by alternating fields (AF) and a procedure of double demagnetization used as a routine when anomalous directions were observed between two successive steps; thermal demagnetization was also performed on selected samples. The remaining 231 specimens were demagnetized at three successive steps. Typical demagnetization diagrams relative to the three sites are shown in Fig. 1. Samples from 659C and 664D have a soft component of magnetization removed by fields between 5 mT (664D) and 15 mT (659C) and then a characteristic component decreasing linearly to the origin. Samples from 665B show also a secondary viscous component with a direction more or less close to the present field entirely removed by AF demagnetization at 15 roT. A total of 46 samples which yielded ambiguous directions (three of which were within polarity reversal) were rejected; no sample from 665B was rejected.

3. Paleomagnetic results The directional changes obtained from 256 individual samples at site 659C are plotted as vector means for each stratigraphic level in Fig. 2. Be-

cause the core was not oriented the declinations have been corrected by adjusting the mean normal value to 0 °. Using the same correction for the entire section, the mean reversed declination is 180 o. The declination changes associated with the reversal identified at 23.15 mbsf are characterized by an abrupt variation from 180 ° to 0 ° and a relatively gradual change in inclination. The mean normal (33 +_ 9.8 °) and reverse (29.7 + 7 °) inclinations are in good agreement with the geocentered dipole field inclination of 33 ° at the site latitude. The absence of data points within 23.37-23.44 mbsf is due to the lack of stable pseudo-single domain grains. In Fig. 1 we show the demagnetization diagrams associated with the uppermost reversed (665B-19) and the lowermost normal (665B-20) consecutive samples from Hole 665B. It is clear from these diagrams that no intermediate direction was obtained. The complete record of the section shown in Fig. 2 is characterized by a sudden change in declination at 13.7 mbsf associated with a large decrease of the magnetization intensity while the inclinations remain close to 0 o. The declinations associated with the reverse and normal polarities are perfectly antipodal. The most detailed record was obtained from Hole 664D and is also plotted in Fig. 2 as vector means for each stratigraphic level. The inclinations are on average close to 0 o in agreement with the equatorial position of the site. The transition identified at 27.8 mbsf is characterized by easterly transitional declinations while the inclination values remain close to 0°. The upper part of the U channel (80 cm above the transition record) show anomalously high inclinations with a general increase in scatter of the directions. These levels are also characterized by a change in sediment colour and they appear to be physically disturbed. The magnetization intensities from Holes 664D and 665B have been normalized with respect to variations in the content of magnetic minerals by the anhysteretic remanent magnetization (0.1 mT direct field, 100 mT AF field) demagnetized at 20 roT. The N R M / A R M variations (Fig. 3) obtained at Hole 664D start decreasing before the directional changes and reach a minimum of 20% during the reversal. The N R M intensity (after 20 mT) and the N R M / A R M variations plotted in Fig. 3 show identical variations. The N R M / A R M de-

374

9di INCLINATION

INCLINATION

80 (50

665B

664D

40

~,

~J .I

20

oo

- ; =

H

0

"~mJ~mlHr m

--20 --40

DEPTH in METERS

13..70

13.60

--80

13.80

DEPT t

273

I DECLI NATION

28.1

27.7

28.5

DECLINATION

180I

ol METERS

13,60

13.65

13.70

13.75

(~

1380

80

.,

27.3

27,7

28.1

,

METERS

28.5

659C

60 ,40

N ~

20

o --40 --60 --80 t

23`2

25O

,

a

23.6

DE PTH

j

24

200 150

100 50

0

~5o23



i

23.2

23.4

i

23,6 DEPTH

i

23.8

METERS

j

24

Fig. 2. Records of the declination and inclination variations across the M a t u y a m a - B r u n h e s reversal obtained from O D P Holes 659C, 664D and 665B after stepwise demagnetization.

375

creases symmetrically around the polarity change represented by the latitude of the virtual geomagnetic pole (VGP) and exhibits a minimum synchronously with the directional jump.

4. Resolution and reliability of the records Before discussing these results in terms of transitional field behavior we must first investigate

NRM/ARM

664D

i

i

J I

I I I

27.3

27.6

28

28.4

METERS'

665B

.

.

13.60

.

.

.

V,"

13.70

,

,

,

DEPTH

13.80

Fig. 3. Variations of the N R M (20 m T ) / A R M ( 2 0 mT) relative to the records from Holes 664D and 665B. The directional changes (represented as the V(3P latitudes) and the N R M changes relative to Hole 665B are shown also for comparison.

376 carefully the quality of these sediments as paleomagnetic recorders. The possible influence of post-depositional remanent magnetization (PDRM) on recordings of reversals have been discussed by several authors [17-19] and there is a common agreement among the studies that P D R M may introduce a smoothing of the transitional directions which could explain in particular the overall longitudinal confinement of the trajectories. It is thus important to attempt to distinguish specific features which may be compatible with post-depositional processes from those inherent to the field behaviour. A first approach to these questions is given by the record from Hole 665B characterized by an absence of transitional directions while the decrease and the recovery of the magnetization intensity are symmetric around the polarity change. Such characteristics might be the consequence of a sedimentary hiatus. However, assuming that the N R M / A R M variations would then be representative of the relative variations of the paleofield intensity, the relationship between the directional changes and the N R M / A R M plot shown in Fig. 3 would imply that the hiatus occurred exactly when the field direction was reversing. Although such an hypothesis cannot definitely be ruled out, the absence of any visible discontinuity in the sedimentary column and the constancy in the mineralogy and chemical composition of the sediment suggest that no major break occurred. Another hypothesis is that this record is the result of mixing of grains with opposite polarities within the individual samples due to post-depositional reorientations. When the number of grains oriented along the normal polarity becomes predominant the declination jumps abruptly from reverse to normal polarity. The downcore changes in the magnetization intensity reflect the variations in the amount of grains with antipodal polarities within the samples. It must be emphasized that this interpretation implies a very rapid transition at the scale of the deposition rate quite compatible with the low value (1.5 c m / k y r ) obtained at Hole 665B. It is also required that both polarities be perfectly antipodal as seen effectively in Fig. 2. Therefore, the changes in the intensity of magnetization are symmetric with respect to the directional changes (Fig. 3) and the depth over which intensity variations are recorded corre-

sponds to the depth over which P D R M was active, i.e. 22 cm at Hole 665B. Interestingly, this result agrees quite well with values reported in the literature based on laboratory experiments [20,21]. If we consider this hypothesis as valid in the case of Hole 665B we must also examine the influence of time averaging due to post depositional reorientations on the other records. Based on the value of 22 cm obtained above for the depth of the lock-in zone we have thus examined possible effects of P D R M on a record similar to 664D characterized by a deposition rate between 3.8 and 4.4 c m / k y r . As the main question underlying this paper deals with the predominance of axisymmetrical components we used the field intensity and directional variations predicted by the model of Fuller et al. [1] for the BrunhesMatuyama transition as seen from the equator by assuming that 20, 30 and 50% of the total energy of the dipole term are respectively redistributed to the gO, g03 and g4° Gauss coefficients. Fuller et al. point out that the form of the synthetic record depends on how much energy goes to each harmonic and their signs. Any other model in which only zonals are included could also have been choosen (as well as the most recent model developed by this group [7] combining low-order zonals with a drifting non-dipole field). A striking and common feature is that (with the exception of a pure octupole) they all are characterized by inclinations passing through + 90 ° or - 90 °. We assumed that the process of remanence smoothing could be simulated by a weighted running average with a Gaussian window (other shapes of window led to identical observations) performed on the vector data X, Y, Z individually. The results and conclusions given below must thus be considered within the limits of this assumption. Variations of the inclination, declination and intensity obtained by using successive windows with lengths of 7, 15 and 22 cm for a transitional interval 7 cm long are shown in Fig. 4a. The main feature is that the amplitude of the inclination peak of the original record remains totally unchanged while the intensity variations are spread over a larger zone and the transition is displaced downwards (wider windows lead to same characteristics). These results are the consequence of competing effects produced by averaging directions restricted within the N - S plane. As the

377 (b)

(a)

l

-56

cm -40

-24

- 10

360 .4

27o DECLINATION

DEPTH -56

-40

-24

-lO

180

-40

-24

DEPTHincm ~10

Fig. 4. (a) Simulation of PDRM smoothingby a weighted running average with Gaussian windows of 7, 15 and 22 cm on a transition entirely dominated by axisymmetrical components. The inclination peak is persistent whatever the width of the window. (b) Same calculations as in figure a for a transitional field with zonal and non zonal components.The declination variations are shortened and the inclination changes still characterized by a significant peak. declination can only be 180 ° or 0 ° the horizontal components are cancelled by averaging and consequently the vertical components are enhanced. The introduction of non-zonal components as (schematically) illustrated in Fig. 4b by a gradual westward rotation of the declination from 180 ° to 360 °, of course, limits these effects and smooths the inclination more effectively. The main feature, however, is that the smoothing does not become stronger for windows wider than 7 cm so that a significant peak of more than 40 ° in amplitude persist in the records. The declination variations appear shorter than on the original record. It is also interesting that the duration of the transition is not significantly affected by P D R M especially if seen from the intermediate directions deviating by more than 30 ° from full polarity. On the contrary, the duration of the intensity changes is considerably enhanced and cannot be reconciled with the original variations if the window of smoothing is larger than the transition length. This result precludes any interpretation of the intensity variations from sediments when recorded with low deposition rates.

The main conclusion emanating from this simulation is that the pattern of directions characteristic for axisymmetric terms dominating during the transition could not be eliminated by the P D R M (if the smoothing in Nature occurs by a mechanism compatible with our calculations). In particular, there is no way of observing significant deviations from the N - S axis (i.e. significant deviations of the declination from 0 ° or 180 °) if the transitional field remained confined close to this axis and the recording of vertical inclinations could not be smoothed out. Because the full polarity states are strongly dominated by zonal terms, these results suggest that zonal harmonics (for records of reversals obtained with low deposition rates) cannot be eliminated but rather advantaged by PDRM. Since these calculations have been performed for a record similar to 664D we can now distinguish whether zonal harmonics were dominating the transition recorded at this site. We are well aware however that the absence of these components would not mean that the record was not affected by post-depositional processes.

378

5. Discussion

would induce variations of the inclination at the equator. It must be noted that no evidence of a perfect axisymmetrical field during a reversal has been reported so far. Recently Merrill and McFadden [22] have proposed that the condition of occurrence of reversals is that the coupling between the dipole family (characterized by the sum of degree n and order m of the spherical harmonic being odd) and the quadrupole families (n + m = even) be critically enhanced at a time when the relative magnitude of the quadrupole terms is high with respect to the dipole family contribution. The present results do not seem to support the presence of zonal harmonics issued from the quadrupole family.

5.1. Inclination variations The equatorial position of Hole 664D is ideal for testing the hypothesis of axisymmetry. Because the field must pass through the vertical, an equatorial site would have the largest variation in inclination (from 0 o to either + 90 o or - 90 ° ) if quadrupolar terms are dominating. Octupolar terms or in a more general sense fields which are symmetric about the equator would have no change in inclination at the equator. Recently, Williams et al. [7] proposed a model combining low-order zonal harmonics with a drifting non-dipole field of intensity and drift similar to the present non-dipole field. The authors show that the directional changes resulting from this configuration are controlled by the phase relationship of the variation of strength of the zonal harmonics and the drifting non-dipole field. However, high inclinations should still be observed at some time during the transition in equatorial records, either related to the zonal terms or to the passage of features of the non-dipole field over the site. The absence of marked changes in inclination (which remained close to 0 °) during the transition in the record from 664D (Fig. 2) is thus inconsistent with the hypothesis that low-order zonal harmonics dominated. The alternative hypothesis, that an octupolar field family dominated the reversal may be difficult to account for since it requires the absence of any other component (in particular the declination should j u m p from 180 ° to 0 ° ) . Indeed, the presence of other components even of weak intensity combined with an octupolar field

5.2. Duration of the transition Another interesting aspect of our data set is based on a comparison of the duration of the reversal at various latitudes. Assuming for the moment that the transition of field is dominated by octupolar terms, durations of the reversal at equatorial latitudes should be significantly shorter than at mid-northern latitudes. Such a latitudinal dependency has been reported for the Matuyam a - B r u n h e s reversal by Clement and Kent [10] from deep-sea cores of the Pacific Ocean. Their study dealt, however, with records with low sedimentation rates (1 c m / k y r and less) which were thus very sensitive to the criteria used in the determinations of the length of the transitional intervals as well as to variations in the deposition rate. Before calculating durations, therefore, it is useful first to examine the consistency between the lengths of the transitional intervals and the sedi-

TABLE 1 Durations of the t r a n s i t i o n - - c a l c u l a t i o n s of the transition duration from the lengths of the transitional intervals and estimates of the sedimentation rate associated with each record Site

Location

Transition length (cm)

Deposition rate ( c m / k y r )

Duration

min.

max.

min.

max.

(kyr) 0.6-2.2 kyr d = 1.4 kyr 2.0-4.4 kyr d = 3.2 kyr 1.8-2.9 kyr d = 2.3 kyr 2.0-2.6 kyr d = 2.3 kyr

Hole 659C

18°0'N, 21°W

2

6

2.7

3.1

Hole 664D

0 ° 0 6 ' N , 23 ° W

9

17

3.8

4.4

Hole 609B

49 ° 8 ' N , 25 ° W

15

17

5.8

8.2

Lake Tecopa

35 ° 5 ' N, 116 o W

10

12

4.5

5.0

379 mentation rate. The boundaries of the transition zones may be defined by using such criteria as the level at which the directions exceed the maximal amplitude of the pre- and post-reversal variations or the occurrence of VGPs latitudes lower than 45 o or 60 o. In order to account for the uncertainties related to the use of different criteria we considered the minimum and maximum thickness obtained for each record using the various methods (Table 1). We assume that differences between the values are representative of the error. Values of the deposition rates using the magnetostratigraphic results obtained for the Brunhes chron and the M a t u y a m a - U p p e r Jaramillo interval were both taken into account although they are similar. We calculated durations for the data from Holes 659C, 664D, DSDP Hole 609B (all located along the same longitude between 0 ° and 50 ° N) and Lake Tecopa [13]. It can be seen in Fig. 5 that there is a good correlation between the lengths of the transitional zones and the deposition rates. Because no transitional direction was recorded at Site 665B, no reliable determination of the transition length can be obtained although the results are consistent with a very short duration of the reversal. Within the limits of the errors described above, the best-guess durations obtained for the O D P Holes and DSDP Hole 609B (Fig. 6) do not show any dependency with site latitude and thus does not support the concentration that octupolar terms be dominant. The mean duration of the reversal is about 2.3 kyr and this value agrees also with the record of Matuyama-Brunhes from Lake Tecopa in California [13] obtained from a very different kind of sediment. This result indicates that the duration of Brunhes-Matuyama was very likely short. Our value differs by almost a factor of two with values usually reported in the literature for other reversals [3,4,9,11]. We have shown above that the interpretation of the intensity changes is difficult. According to the results in Fig. 3 the dipole intensity (Hole 664D) would have been reduced for about 4-5 kyr (about twice as long as the transition defined by directional changes) and this duration may be exaggerated due to post depositional processes. Larger secular variation can be expected during periods of low intensity before and after the transition. Both declination and inclination records are effectively characterized by distinctive oscillations but

0 -10" -20 UJ

-30"

~E t-

-40 -50" -60

i

i

200

400

i



600 AGE in Kyr

z

800

1000

a

DEPOSITION RATE in cmlkyr

l 659C

.1

I

665B

TRANSITION 5

10

15

LENGTH in cm

Fig. 5. (a) Age versus depth of the transition boundaries for every core. The sedimentation rates do not show major changes. (b) Relationship between the lengths of the transitional zones and the deposition rates. There is a good correlation between the data sets. The deposition rate corresponding with Lake Tecopa was deduced from the thickness of sediment between the Matuyama-Brunhes transition and the lower and upper limits obtained on the position of the upper Jaramillo from preliminary magnetostratigraphy (no estimate of the thickness relative to the Brunhes period can be obtained from the section).

the time resolution of the record as well as the number of directions associated with full polarity states are too short to draw a firm conclusion.

5.3. Rotated directions D', I ' and VGP paths We calculated the rotated directions D ' , I ' [9] and the VGP paths associated with each record from the declination, inclination data presented in

380 5i DURAT ION in Kyr

3

TECOPA

664D

:

2

609B

±

659C[

li

~



2~

,LATITUDE 5o~

Fig, 6. Plot of the durations calculated for each record of Matuyama-Brunhes with error bars versus site latitude. The duration associated with the record from Lake Tecopa is shown in dashed lines.

Fig. 2. The few intermediate directions at Hole 659C result in a poorly defined trajectory with three intermediate VGPs (latitudes > 30 ° from the rotational axis) in the northern hemisphere and no VGPs in the southern hemisphere (Fig. 7a). The equal area plot of the intermediate " r o t a t e d " directions (Fig. 7b) show three intermediate directions restricted to the lower hemisphere with declinations between 125 ° and 180 °. The V G P path recorded at Hole 664D (Fig. 8a) is described by 17 intermediate points as defined above distributed along a mean longitude 70 °E, i.e. 90 ° to the east of the site longitude. The D ' , I ' plots shown in Fig. 8b are characterized by 22 "intermediate" rotated directions with declinations close to the east axis.

5. 4. Geometry of the transitional field Any attempt to test possible field geometries during the M a t u y a m a - B r u n h e s reversal requires a large number of detailed paleomagnetic records. The last reversal can easily be identified (in principle) in sedimentary sequences and consequently m a n y distinct records have been published. In a review paper on reversals, a classification of records in terms of quafity A, B and C was published by Fuller et al. [1]. Owing to the subsequent developments in Paleomagnetism additional requirements (including stepwise demagnetization of all the transitional samples and analysis of the magnetic mineralogy) have been proposed [13]. We are aware that such criteria may be considered subjective and their use is a matter of debate. It is frustrating, however, that after selecting records consistent with these requirements,only the two

records of M a t u y a m a - B r u n h e s from Lake Tecopa in California and that from O D P Hole 664D are retained. In fact most of the other studies are excluded either because the samples have not been stepwise demagnetized, the V G P paths are too poorly defined or the presence of a primary magnetization is not clearly established. The results obtained from Lake Tecopa and Hole 664D, considered separately, lead to different interpretations. The lack of inclination variations as well as the pattern of the directional changes at Hole 664D indicate that axisymmetrical components were not recorded at the equator. In contrast the results from Lake Tecopa taken alone could be interpreted in terms of axisymmetry. If we deal strictly with these two records we are thus somewhat tantalized by ambiguity. The record obtained from Hole 609B by Clement and Kent [14] stands out from the rest of the studies because the transitional samples have been demagnetized at at least three steps. Another reason why we consider this record is that its location is at the same longitude as Hole 664D but much further north. Indeed, similar V G P paths are predicted from sites located on the same longitude, if the transitional field is dominated by zonal terms. We observe in Fig. 9 that the two trajectories not only do not follow a c o m m o n path but are (surprisingly) constrained along almost antipodal longitudes 90 ° west and 90 ° east from the sites meridian respectively. Moreover, a n o r m a l - t r a n s i tional-normal rebound is apparent on the 609B record which is not observed at the equator. Similar observations are drawn from the equal area plots of the rotated directions D ' , I ' (Fig. 8) which clearly indicate that the discrepancy results from the declination variations rotating eastwards in the first case (664D) and westwards in the other case (609B). We are thus faced with puzzling observations which could leave us somewhat suspicious as to the reliability of the present records. However, no artefact due to incomplete cleaning of secondary components of magnetization or problems linked to the magnetic mineralogy can be strongly suspected. We have also discussed above possible effects of P D R M on these records. The global appearance of the VGP paths (as, e.g., their longitudinal confinement) may be the consequence of smoothing effects and the influence of this factor

381

(cl)

N

S

(b)

659C

Fig. 7. VGP paths and rotated directions D', 1' obtained from ODP Hole 659C.

certainly cannot be underestimated but it is difficult to imagine that antipodal transitional directions would have resulted from smoothing if the field variations were identical at both Sites 664D and 609B. It is also difficult to explain why the structure of the directional changes would be so different. Obviously, high-resolution records are required to observe in detail the transitional field variations. Unfortunately, none of these records could be duplicated and this condition seems to be necessary in order to firmly establish the reliability of paleomagnetic studies of reversals. Such a requirement is however extremely difficult to be

met especially when dealing with marine cores from ODP protected by a severe (but maybe justified) sampling policy.

6. Conclusion The present results give insights on field behavior during the last reversal by suggesting that zonal terms did not dominate during the transitional period. It seems that the characteristics of the records would be consistent with a complicated structure of the transitional field involving harmonics of at least degree three but the

382

S

(b)

(c)

664D

609B

Fig. 8. VGP paths and rotated directions D ' , I ' from O D P Hole 664D. For comparison we show also the D ' , f from DSDP Hole 609B.

record obtained

383

60°N

o ,

o

60

• 664D

120

180

240

TECOPA

300

0

'

' 6o60°S

609B

Fig. 9. VGP paths relative to the records of Matuyama-Brunhes from ODP Hole 664D, Lake Tecopa and DSDP Hole 609B. The location of the sites is shown on the planisphere.

resolution and the number of records are too poor to draw firm conclusions. In a paper related to a statistical model for geomagnetic reversals, Constable [23] suggests that during historical and archeological times the geodynamo exhibited complex field variations and there is no apparent reason for the field to be simpler during transitions. The model proposed is based on a statistical description of the geomagnetic secular variation over the last 5 Myr that has been extended to include reversals by introducing a time dependence so that any Gauss coefficient at any instant in time depends on its value at present time plus a random component. Reversals are simulated by allowing the axial dipole to decay and grow back with opposite sign while letting the rest of the secular variation continue as usual. Many of the features of reversal records are predicted, such as looping and loitering in the VGP paths and a fairly rapid transition between northern and southern latitudes. There is no requirement that the VGP path be related to the site longitude but appears to be located within a region determined by the field configuration at the onset of the reversal and for any given reversal may vary from site to site. A model of interacting dipoles undergoing spontaneous reversals proposed recently by Mazaud and Laj [24] predict VGP trajectories not related to the site but mostly constrained in longitude. Such features could be compatible with the present results. Further tests of these models require a large number of detailed

records with a special attention to periods immediately preceeding the transition.

Acknowledgements The paper benefited greatly from the reviewers comments (R. Merrill and two anonymous reviewers). We are pleased to thank C. Laj for very valuable suggestions. We enjoyed discussions with R. Weeks and C. Constable. A. Mennel and D. Clark are acknowledged for help in measurements and technical assistance. The authors enjoyed discussions with P. Tucholka, R. Coe, and K. Hoffman. Support was provided for this study by NATO, the CNRS and the CEA for J.-P.V. and by NSF grants EAR 8515743 and Texas A & M 777348 to L.T. This is CFR Contribution No. 1024.

References 1 M. Fuller, I. Williams and K.A. Hoffman, Paleomagnetic records of geomagnetic field reversals and the morphology of the transitional fields, Rev. Geophys. Space Phys. 17, 179-203, 1979. 2 M. Prrvot, E.A. Mankinen, C.S. Gromm6 and R.S. Coe, The Steens Mountain polarity transition: field intensity variations and discussion of reversal models, J. Geophys. Res. 90, 10417-10448, 1985. 3 J.-P. Valet and C. Laj, Invariant and changing transitional field configuration in a sequence of geomagnetic reversals, Nature 311, 552-555, 1984. 4 E. Herrero-Bervera and F. Theyer, Non axisymmetric behaviour of Olduvai and Jararnillo polarity transitions r e -

384

5

6

7

8

9

10

11

12

13

14

corded in north-central Pacific deep-sea sediments, Nature 322, 159-162, 1986. C. Laj, S. Guitton, C. Kissel and A. Mazaud, Complex behaviour of the geomagnetic field during three successive polarity reversals, 11-12 m.y.B.P., J. Geophys. Res. 93, 11655-11665, 1988. S. Bogue and R.S. Coe, Transitional paleointensities from Kauai, Hawai, and geomagnetic reversals models, J. Geophys. Res. 89, 10341-10354, 1984. I. Williams, R. Weeks and M. Fuller, A model for transition fields during geomagnetic reversals, Nature 332, 719-720, 1988. F. Theyer, E. Herrero-Bervera and V. Hsu, The zonal harmonic models of polarity transitions: a test using successive reversals, J. Geophys. Res. 90, 1963-1982, 1985. K.A. Hoffman, Geomagnetic reversals and excursions: their paleomagnetic record and implications for the geodynamo, Rev. Geophys. 21,614-620, 1983. B.M. Clement and D.V. Kent, A detailed record of the lower Jaramillo polarity transition from a southern hemisphere deep-sea core, Nature 310, 488-491, 1984. I. Williams and M. Fuller, Zonal harmonics models of reversal transition fields, J. Geophys. Res. 86, 11657-11665, 1981. L. Tauxe, J.P. Valet and 1. Bloemendal, Magnetostratigraphy of Leg 108 sediments, in: Scientific Results ODP Leg 108, B. Ruddiman, M. Sarntheim and J. Baldauf, eds., in press, 1989. J.P. Valet, L. Tauxe and .D.C. Clark, The MatuyamaBrunhes transition recorded from lake Tecopa sediments, Earth Planet. Sci. Lett. 87, 463-472, 1988. B.M. Clement and D.V. Kent, Geomagnetic polarity transition records from five hydraulic piston core sites in the north Atlantic, in: B. Ruddiman, Init. Rep. DSDP 94, 1987.

15 L. Tauxe, J.L. LaBrecque, R. Dodson and M. Fuller, A new technique for paleomagnetic analysis of hydraulic piston cores (abstract), EOS 64, 219, 1983. 16 J.W. King, S.K. Banerjee, J. Marvin and O. Ozdemir, A comparison for determining the relative grain size of magnetite in natural materials: some results for lake sediments, Earth Planet. Sci. Lett. 59, 404-419, 1982. 17 J.C. SaUoway, The structure of the geomagnetic field during polarity transitions, 296 pp., Ph.D. Thesis, University of Edinburgh, 1983. 18 J.-P. Valet, Inversions g~omagn&iques du Miocene suprrieur en Cr&e. Modalitrs du renversement et caractrristiques du champ de transition, 195 pp., Thrse d'rtat, Universit6 d'Orsay, 1985. 19 K.A. Hoffman and B.S. Slade, Polarity transition records and the acquisition of remanence: a cautionary note, Geophys. Res. Lett. 13, 483-486, 1986. 20 R. Lovlie, The intensity pattern of post-depositional remanence acquired in some marine sediments deposited during a reversal of the external magnetic field, Earth Planet. Sci. Lett. 30, 209-214, 1976. 21 I. Katsura and S. Sasajima, Note on the experimental acquisition of remanence in sediment during a field reversal, Mem. Fac. Sci. Kyoto Univ., Ser. Geol. Mineral., pp. 67-74, 1983. 22 R.T. Merrill and P.L. McFadden, Secular variation and the origin of geomagnetic field reversals, J. Geophys. Res., 93, 11589-11597, 1988. 23 C. Constable, A simple statistical model for geomagnetic reversals (abstract), EOS 69, 1988 and J. Geophys. Res, in preparation. 24 A. Mazaud and C. Laj, Simulation of geomagnetic reversals by a model of interacting dipole sources, Earth Planet. Sci. Lett. 92, 299-306, 1989.