Skip to main content

Main menu

  • Home
  • Articles
    • Current
    • Special Feature Articles - Most Recent
    • Special Features
    • Colloquia
    • Collected Articles
    • PNAS Classics
    • List of Issues
  • Front Matter
    • Front Matter Portal
    • Journal Club
  • News
    • For the Press
    • This Week In PNAS
    • PNAS in the News
  • Podcasts
  • Authors
    • Information for Authors
    • Editorial and Journal Policies
    • Submission Procedures
    • Fees and Licenses
  • Submit
  • Submit
  • About
    • Editorial Board
    • PNAS Staff
    • FAQ
    • Accessibility Statement
    • Rights and Permissions
    • Site Map
  • Contact
  • Journal Club
  • Subscribe
    • Subscription Rates
    • Subscriptions FAQ
    • Open Access
    • Recommend PNAS to Your Librarian

User menu

  • Log in
  • My Cart

Search

  • Advanced search
Home
Home
  • Log in
  • My Cart

Advanced Search

  • Home
  • Articles
    • Current
    • Special Feature Articles - Most Recent
    • Special Features
    • Colloquia
    • Collected Articles
    • PNAS Classics
    • List of Issues
  • Front Matter
    • Front Matter Portal
    • Journal Club
  • News
    • For the Press
    • This Week In PNAS
    • PNAS in the News
  • Podcasts
  • Authors
    • Information for Authors
    • Editorial and Journal Policies
    • Submission Procedures
    • Fees and Licenses
  • Submit
Research Article

Multidecadally resolved polarity oscillations during a geomagnetic excursion

Yu-Min Chou, Xiuyang Jiang, Qingsong Liu, Hsun-Ming Hu, Chung-Che Wu, Jianxing Liu, Zhaoxia Jiang, Teh-Quei Lee, Chun-Chieh Wang, Yen-Fang Song, Cheng-Cheng Chiang, View ORCID ProfileLiangcheng Tan, View ORCID ProfileMahjoor A. Lone, Yongxin Pan, Rixiang Zhu, Yaoqi He, Yu-Chen Chou, An-Hung Tan, Andrew P. Roberts, Xiang Zhao, and Chuan-Chou Shen
  1. aHigh-Precision Mass Spectrometry and Environment Change Laboratory, Department of Geosciences, National Taiwan University, 10617 Taipei, Taiwan, Republic of China;
  2. bDepartment of Ocean Science and Engineering, Southern University of Science and Technology, 518055 Shenzhen, People’s Republic of China;
  3. cKey Laboratory of Humid Subtropical Eco-Geographical Processes, Ministry of Education, College of Geography Science, Fujian Normal University, 350117 Fuzhou, People’s Republic of China;
  4. dLaboratory for Marine Geology, Qingdao National Oceanography Laboratory for Science and Technology, 266237 Qingdao, People’s Republic of China;
  5. eKey Laboratory of Marine Sedimentology and Environmental Geology, First Institute of Oceanography, State Oceanic Administration, 266061 Qingdao, People’s Republic of China;
  6. fCollege of Marine Geosciences, Ocean University of China, 266100 Qingdao, People’s Republic of China;
  7. gInstitute of Earth Sciences, Academia Sinica, 11529 Taipei, Taiwan, Republic of China;
  8. hNational Synchrotron Radiation Research Center, 30076 Hsinchu, Taiwan, Republic of China;
  9. iState Key Laboratory of Loess and Quaternary Geology, Institute of Earth Environment, Chinese Academy of Sciences, 710061 Xi’an, People’s Republic of China;
  10. jInstitute of Global Environmental Change, Xi’an Jiaotong University, 710049 Xi’an, People’s Republic of China;
  11. kInstitute of Geology and Geophysics, Chinese Academy of Sciences, 100029 Beijing, People’s Republic of China;
  12. lCollege of Tourism and Air Service, Guizhou Minzu University, 550025 Guiyang, People’s Republic of China;
  13. mDepartment of Mechanical Engineering, Chien Hsin University of Science and Technology, 32097 Taoyuan, Taiwan, Republic of China;
  14. nResearch School of Earth Sciences, Australian National University, Canberra, ACT 2601, Australia;
  15. oResearch Center for Future Earth, National Taiwan University, 10617 Taipei, Taiwan, Republic of China

See allHide authors and affiliations

PNAS September 4, 2018 115 (36) 8913-8918; first published August 20, 2018; https://doi.org/10.1073/pnas.1720404115
Yu-Min Chou
aHigh-Precision Mass Spectrometry and Environment Change Laboratory, Department of Geosciences, National Taiwan University, 10617 Taipei, Taiwan, Republic of China;
bDepartment of Ocean Science and Engineering, Southern University of Science and Technology, 518055 Shenzhen, People’s Republic of China;
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • For correspondence: river@ntu.edu.tw chouym@sustc.edu.cn
Xiuyang Jiang
cKey Laboratory of Humid Subtropical Eco-Geographical Processes, Ministry of Education, College of Geography Science, Fujian Normal University, 350117 Fuzhou, People’s Republic of China;
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Qingsong Liu
bDepartment of Ocean Science and Engineering, Southern University of Science and Technology, 518055 Shenzhen, People’s Republic of China;
dLaboratory for Marine Geology, Qingdao National Oceanography Laboratory for Science and Technology, 266237 Qingdao, People’s Republic of China;
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Hsun-Ming Hu
aHigh-Precision Mass Spectrometry and Environment Change Laboratory, Department of Geosciences, National Taiwan University, 10617 Taipei, Taiwan, Republic of China;
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Chung-Che Wu
aHigh-Precision Mass Spectrometry and Environment Change Laboratory, Department of Geosciences, National Taiwan University, 10617 Taipei, Taiwan, Republic of China;
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Jianxing Liu
dLaboratory for Marine Geology, Qingdao National Oceanography Laboratory for Science and Technology, 266237 Qingdao, People’s Republic of China;
eKey Laboratory of Marine Sedimentology and Environmental Geology, First Institute of Oceanography, State Oceanic Administration, 266061 Qingdao, People’s Republic of China;
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Zhaoxia Jiang
fCollege of Marine Geosciences, Ocean University of China, 266100 Qingdao, People’s Republic of China;
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Teh-Quei Lee
gInstitute of Earth Sciences, Academia Sinica, 11529 Taipei, Taiwan, Republic of China;
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Chun-Chieh Wang
hNational Synchrotron Radiation Research Center, 30076 Hsinchu, Taiwan, Republic of China;
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Yen-Fang Song
hNational Synchrotron Radiation Research Center, 30076 Hsinchu, Taiwan, Republic of China;
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Cheng-Cheng Chiang
hNational Synchrotron Radiation Research Center, 30076 Hsinchu, Taiwan, Republic of China;
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Liangcheng Tan
iState Key Laboratory of Loess and Quaternary Geology, Institute of Earth Environment, Chinese Academy of Sciences, 710061 Xi’an, People’s Republic of China;
jInstitute of Global Environmental Change, Xi’an Jiaotong University, 710049 Xi’an, People’s Republic of China;
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Liangcheng Tan
Mahjoor A. Lone
aHigh-Precision Mass Spectrometry and Environment Change Laboratory, Department of Geosciences, National Taiwan University, 10617 Taipei, Taiwan, Republic of China;
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Mahjoor A. Lone
Yongxin Pan
kInstitute of Geology and Geophysics, Chinese Academy of Sciences, 100029 Beijing, People’s Republic of China;
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Rixiang Zhu
kInstitute of Geology and Geophysics, Chinese Academy of Sciences, 100029 Beijing, People’s Republic of China;
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Yaoqi He
lCollege of Tourism and Air Service, Guizhou Minzu University, 550025 Guiyang, People’s Republic of China;
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Yu-Chen Chou
aHigh-Precision Mass Spectrometry and Environment Change Laboratory, Department of Geosciences, National Taiwan University, 10617 Taipei, Taiwan, Republic of China;
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
An-Hung Tan
mDepartment of Mechanical Engineering, Chien Hsin University of Science and Technology, 32097 Taoyuan, Taiwan, Republic of China;
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Andrew P. Roberts
nResearch School of Earth Sciences, Australian National University, Canberra, ACT 2601, Australia;
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Xiang Zhao
nResearch School of Earth Sciences, Australian National University, Canberra, ACT 2601, Australia;
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Chuan-Chou Shen
aHigh-Precision Mass Spectrometry and Environment Change Laboratory, Department of Geosciences, National Taiwan University, 10617 Taipei, Taiwan, Republic of China;
oResearch Center for Future Earth, National Taiwan University, 10617 Taipei, Taiwan, Republic of China
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • For correspondence: river@ntu.edu.tw chouym@sustc.edu.cn
  1. Edited by Lisa Tauxe, University of California, San Diego, La Jolla, CA, and approved July 23, 2018 (received for review November 22, 2017)

  • Article
  • Figures & SI
  • Info & Metrics
  • PDF
Loading

Significance

A stalagmite-based paleomagnetic record of the post-Blake excursion reveals details of repeated centennial–millennial interhemispheric polarity drifts and saw-tooth inclination oscillations during periods of low geomagnetic field intensity at 100 thousand years before present. One surprisingly abrupt centennial reversal transition occurred in 144 ± 58 years (2σ) and provides unprecedented evidence that raises fundamental questions about the speed of geomagnetic field shifts. Such rapid polarity changes could severely affect satellites and human society in the future if the current geomagnetic field intensity continues to decrease.

Abstract

Polarity reversals of the geomagnetic field have occurred through billions of years of Earth history and were first revealed in the early 20th century. Almost a century later, details of transitional field behavior during geomagnetic reversals and excursions remain poorly known. Here, we present a multidecadally resolved geomagnetic excursion record from a radioisotopically dated Chinese stalagmite at 107–91 thousand years before present with age precision of several decades. The duration of geomagnetic directional oscillations ranged from several centuries at 106–103 thousand years before present to millennia at 98–92 thousand years before present, with one abrupt reversal transition occurring in one to two centuries when the field was weakest. These features indicate prolonged geodynamo instability. Repeated asymmetrical interhemispheric polarity drifts associated with weak dipole fields likely originated in Earth’s deep interior. If such rapid polarity changes occurred in future, they could severely affect satellites and human society.

  • multidecadally resolved polarity oscillations
  • geomagnetic excursion
  • abrupt reversal transition
  • asymmetrical interhemispheric polarity drifts
  • stalagmite U-Th dating

The geomagnetic field, which is generated by convection in Earth’s liquid outer core (1), has existed for at least 3.45 billion y (2). It shields Earth from direct impact of solar wind and cosmic radiation, and its intensity and direction vary on timescales from milliseconds to tens of millions of years (3). Geomagnetic polarity reversals have occurred irregularly throughout Earth history (4), with the most recent polarity reversal defining the Matuyama–Brunhes boundary (MBB) that is astronomically dated at 773 ± 1 thousand years ago (ka) (5). At least 13 short-lived geomagnetic excursions have occurred since the MBB (e.g., the Laschamp excursion at 41 ka) (6).

Polarity reversals and excursions occur during periods with low geomagnetic field intensities (7, 8) (SI Appendix, Fig. S1), which could severely disorder the ionosphere, atmosphere, and biosphere (9, 10). Infrastructure and facilities for communication, transportation, and livelihoods can be threatened by such conditions (11). The contemporary geomagnetic field intensity is decreasing, and it has been suggested that a field reversal is impending (11, 12). Understanding details of transitional field behavior and obtaining precise estimates of transition duration are important for the Earth sciences, and they are also potentially significant for society and the biosphere.

Since the discovery of geomagnetic reversals in igneous rocks in the 1920s (13), diverse geological archives, such as lava flows, sedimentary rocks, and marine/lacustrine sediments, have been studied over recent decades (14, 15). Igneous rocks can provide ages of reversals and excursions via radioisotopic dating techniques (16, 17). However, intermittent eruption of igneous rocks hinders development of continuous records for studying transitional field configurations (18). Sedimentary records offer continuity (19) but are limited intrinsically by difficulties in obtaining radioisotopic chronologies and by smoothing of sedimentary paleomagnetic signals (20). Constraining details of transitional field structure during polarity reversals and excursions need accurate age control.

Speleothems have proven to be highly suitable for paleomagnetic reconstructions (21⇓–23) because of their continuous formation, short magnetization lock-in time, and minor postdepositional complications. Speleothem records have been obtained for the Laschamp excursion in North America (23) at 41.10 ± 0.35 ka and the Blake excursion (BE) in northern Spain (22) at 119–112 ka. Ongoing refinement of high-sensitivity superconducting magnetometry (24) and radioisotopic U-Th chronology (25, 26) offers the possibility of reconstructing geomagnetic excursions at centennial scale or better. We here report a multidecadally resolved Chinese stalagmite paleomagnetic record for the 107- to 91-ka interval with closely spaced U-Th ages. The record provides important insights into ancient geomagnetic field behavior of excursional events that reveal repeated centennial–millennial interhemispheric polarity drifts and saw-tooth inclination oscillations with one abrupt centennial reversal transition.

U-Th Dating and Paleomagnetic Analyses

Stalagmite SX11 was collected in a chamber 800 m from the entrance of Sanxing Cave (107°11′ E, 27°22′ N, 2 km in length, 720-m altitude), Guizhou Province, southwestern China (Fig. 1 A and B), in 2011. The stalagmite is 1 m in length and 8 cm in diameter, and it has a candle-like shape and yellow to dark brown color (Fig. 1C). At the micrometer scale, irregularly shaped detrital particles are observed within laminae (SI Appendix, Fig. S2). Subsamples with 5-mm thickness were used for U-Th dating to provide 2σ precision of ±10s of years (Methods Summary). The stalagmite, which was deposited from 106.66 to 91.23 ka, contains a hiatus from 100.81 to 98.77 ka (Fig. 1C and SI Appendix, Figs. S3 and S4).

Fig. 1.
  • Download figure
  • Open in new tab
  • Download powerpoint
Fig. 1.

Sampling location and photograph of Sanxing stalagmite SX11. (A) Location map of Sanxing Cave (red circle). (B) The sampling site of SX11 in the cave chamber (red star). (C) Stalagmite photograph. The white line is a tape measure, and the dashed brown curved line denotes a hiatus at a depth of 45.5 cm.

One hundred ninety-four subsamples were cut along horizontal growth laminae to avoid possible distorted paleomagnetic directions at stalagmite slope edges (Methods Summary). The natural remanent magnetization (NRM) ranges from 10−3 to 10−5 A/m and averages 10−4 A/m. Alternating field (AF) demagnetization results for most (180 of 194) subsamples indicate the presence of a weak viscous remanent magnetization (VRM) that is removed after AF treatment at 10–25 mT. At high demagnetization fields, the NRM is dominated by one stable paleomagnetic component that generally decays linearly to the origin of demagnetization diagrams and that defines the characteristic remanent magnetization (ChRM) with AF < 50 mT (Methods Summary, Fig. 2 A and C, SI Appendix, Fig. S5, and Dataset S2). Thermal demagnetization of a three-axis isothermal remanent magnetization (IRM) (27) indicates that the soft component (Bcr < 50 mT) is fully thermally demagnetized at 575 °C, which indicates that the ChRM carrier is dominantly magnetite (Fig. 2 B and D and SI Appendix, Fig. S6). The ratio of anhysteretic remanent magnetization (ARM) to saturation isothermal remanent magnetization applied in a 1-T pulse field (SIRM1T) suggests that the magnetite particles are in the single domain to vortex state size range (SI Appendix, Fig. S7). These properties are typical of the pedogenic and detrital magnetite particles that are likely to be delivered to the stalagmite by drip waters. The stability of the recorded paleomagnetic directions throughout the stalagmite (SI Appendix, Fig. S5) confirms that the magnetite particles must be in this magnetically stable particle size range. A one-component ChRM with a signal dominated by magnetite indicates that the magnetic particles acquired a depositional remanent magnetization without chemical alteration/precipitation of magnetic carriers during stalagmite growth (28).

Fig. 2.
  • Download figure
  • Open in new tab
  • Download powerpoint
Fig. 2.

Orthogonal projections of AF demagnetization results for the NRM and thermal demagnetization results for a three-axis IRM. Typical orthogonal vector component plots for depth intervals of (A) 565–570 and (C) 970–975 mm with vector endpoints projected onto the vertical (white squares) and horizontal (black squares) planes. Letters of E, W, S, N, U, and D denote east, west, south, north, up, and down, respectively. ChRM line fits (red lines) indicate the presence of one stable component with NRM values of 6.10 × 10−3 A/m at 970–975 mm and as low as 5.30 × 10−5 A/m at 565–570 mm. Thermal demagnetization curves for a three-axis IRM at depth intervals of (B) 560–580 and (D) 960–980 mm imparted with fields of 2.5 (red circles), 0.5 (yellow circles), and 0.05 T (blue triangles). Unblocking up to 575 °C indicates that magnetite is consistently the dominant magnetic mineral.

Geomagnetic excursions are often defined as a deviation of the virtual geomagnetic pole (VGP) from the normal paleosecular variation range (NPSVR) (29). Stalagmite SX11 was not oriented azimuthally, and VGPs could not be calculated from paleomagnetic directions. NPSVR for inclination at the site was calculated with a paleosecular variation model for VGP scatter (30). The gradient color bands in Fig. 3B indicate inclination distributions from modeled paleosecular variation for normal (yellow in Fig. 3B) and reversed (blue in Fig. 3B) polarities at 27°22′ N, where the most intense color expresses the highest probability; geocentric axial dipole field inclinations for both polarities ±45.8° are also plotted with dashed lines in Fig. 3B (Methods Summary and SI Appendix, Fig. S8). Two excursions, post-Blake 1 (PB1) and PB2, which are defined based on inclination, are recorded in the stalagmite inclination record at 105.8–103.5 and 98.5–92.0 ka, respectively (Fig. 3). These excursions coincide with the PB excursion (PBE) events at 105 and 95 ka, respectively, which were discovered in Lac du Bouchet, France (31). 40Ar/39Ar dating of lava flows from Lipari Island and New Mexico indicates that the beginning of the PBE was at 104–105 ka (Lipari: 105 ± 1 ka; New Mexico: 104 ± 12 ka) (17). The PBE is also recorded in Icelandic lava flows at 94 ± 8 ka (32).

Fig. 3.
  • Download figure
  • Open in new tab
  • Download powerpoint
Fig. 3.

Paleomagnetic record for Sanxing stalagmite SX11 at 107–91 ka. (A) Two excursions, PB1 and PB2, are indicated in the magnetostratigraphic column. (B) Inclination record with 95% confidence limit of ChRM (calculated following ref. 51 and shaded in gray). 230Th ages and 2σ errors are plotted. Numbers with uncertainties are 230Th ages per data point (unit: years before AD 1950). Ages in italics were calculated using the linear age model to interpolate between two adjacent determined 230Th ages with propagated errors. Horizontal dashed lines are inclinations for a geocentric axial dipole field at the site. Yellow and blue color gradient bands denote the probability distributions of modeled paleosecular variation at 27°22′ N for normal and reversed polarity, respectively. The letters a–f and numbers 1–4 represent different intervals of geomagnetic field behavior described in the text. (C) Relative declination record with 95% confidence limit of ChRM (shaded in gray). (D) Virtual axial dipole moment (VADM) records of the SINT-800 (blue line), PISO-1500 (red line), and PADM2M (black line) RPI stacks with 95% confidence limits (39⇓–41). The vertical light gray bar indicates the location of the hiatus.

For the PB1 excursion, at least three submillennial to millennial directional drift cycles (Fig. 3, a and b, 650 y; c and d, 1,100 y; e and f, 500 y) with inclinations that change from the NPSVR to shallow excursional angles are identified (Fig. 3). The first cycle (PB1-a and PB1-b) involves a rapid inclination deviation from 27.3° to 0.8° and a westward declination drift of 60° in 100 y (a in Fig. 3B; inclination enlarged in SI Appendix, Fig. S9) followed by a 150-y northward inclination shift back to NPSVR values and an eastward declination drift of 60°. Directions then remain in the NPSVR for 400 y (b in Fig. 3B; inclination enlarged in SI Appendix, Fig. S9). The second cycle (PB1-c and PB1-d in Fig. 3B; inclination enlarged in SI Appendix, Fig. S9) has a similar pattern, where the inclination deviated from 49.7° to −2.9° in 70 y (PB1-c), rebounded back to NPSVR values in 70 y, and then remained in the NPSVR for 900 y. The counterpart declination change involved a rapid 120° eastward movement followed by westward drift back to NPSVR values (PB1-d). The structure of the third cycle (PB1-e and PB1-f) differs from that of the first two cycles. The inclination deviated slowly from 27.3° to 8.8° in 350 y (PB1-e) and then returned in 150 y (PB1-f). The declination changed within a 60° range. The PB1-e and PB1-f record has large uncertainties with low stalagmite growth rates, which prevent clear visualization of field behavior. After 103 ka, inclinations remained within the NPSVR, and relative declinations were roughly constant for at least two millennia.

The PB2 excursion has two large inclination drift cycles (Fig. 3, 1 and 2, 2,300 y; 3 and 4, 3,900 y) from 98 to 92 ka (Fig. 3). The cycle at 98–96 ka involved an abrupt drift to full opposite polarity with inclination changing from 30.0° to −64.6° (one data point at 98.22 ka) and a westward declination shift of 120° from 98,360 ± 39 to 98,216 ± 43 y (1 in Fig. 3B; inclination enlarged in SI Appendix, Fig. S9) followed by a 1,400-y period with excursional inclinations from −19.5° to −0.1° and an inclination return to NPSVR over 800 y, with declination drifts of >60° (PB2-2). The second cycle at 96–92 ka records a reversal process with a gradual 40° inclination shift with four 250-y fluctuations, each recorded by multiple data points, from 62.3° to −20.9°. During this 1,050-y interval of PB2-3, the declination rotated eastward by >60° over the first 950 y and rotated westward by >60° during the final 100 y. Through the 2,950-y PB2-4 interval, excursional inclinations between −9.0° and 3.5° occurred over a 2,350-y period, with declination drifts of >30°. Directions then returned to NPSVR values slowly, with declination rotating >60° eastward in 300 y. The geomagnetic structure of these multimillennial cycles is similar to that of the first two submillennial cycles in the PB1 excursion (a–b and c–d), although with different timescales.

Asymmetrical Interhemispheric Polarity Oscillations

The multicycle variations during the PBE illustrated in Fig. 3 are also observed for the BE over the 119- to 112-ka interval in a stalagmite from northern Spain (22) that lacks precise age control. Such a detailed multidecadally resolved drift structure is only precisely dated in the Sanxing stalagmite. The multicentennial to millennial inclination drift cycles (Fig. 3B) are characterized by an abrupt southward shift and slow relaxation back to the NPSVR (except PB1-e and PB1-f), which are expressed as an asymmetrical saw-tooth oscillation. For the PB1 excursion, the three observed inclination events only drift to ∼0°, whereas one fully reversed polarity direction is observed in the PB2 excursion. The longest saw-tooth inclination oscillation is 4,000 y for the cycle from PB2-3 to PB2-4 (Fig. 3B), and the shortest is only 650 y for the period from PB1-a to PB1-b. Even with large differences in oscillation durations, a similar geomagnetic structure is depicted for the multicentennial to millennial cycles with an abrupt southward shift and slower rebound (Fig. 3B). The periodicity and amplitude of directional deviations of saw-tooth oscillations are amplified in PB2 with respect to PB1. The stalagmite BE record from northern Spain has a low inclination departure a few thousand years before the full reversed polarity event (22). The PBE record from SX11 also has a low inclination PB1 period before PB2. The low inclination event could be a precursor before the large amplitude geomagnetic excursion.

The most noticeable features of our record are the occurrence of two dynamical changes toward negative inclination at PB2. The one-millennium PB2-3 drift period from 96.0 to 94.9 ka is characterized by four 250-y fluctuations (Fig. 3B). Another abrupt movement occurred in 144 ± 58 y during PB2-1 (Fig. 3B and SI Appendix, Fig. S9). These polarity transition periods are much shorter than previously estimated millennial-scale reversal durations from volcanic records (33). The most surprising observation is the abrupt centennial polarity shift at PB2-1, which possibly occurred within one century. Such an extremely rapid polarity drift has not been shown before. Previous hotly debated rates from lava (34) and lacustrine sedimentary (35) records suggested extreme directional changes of 1° per week (34) or less than a century for the transition between polarity states (35). However, the rapid reversal indicated by lava flows was rescinded in later studies (36), and the lacustrine records (35) have been suggested (37) and then shown to be affected by chemical remagnetization (38).

Mechanisms of Asymmetrical Polarity Drifts

Transitional field morphologies and underlying physical mechanisms have never been established fully due to the lack of detailed geomagnetic excursion records with precise dating. Our Sanxing stalagmite record offers a radioisotopically dated dynamic view of geomagnetic field behavior on multidecadal to millennial timescales associated with reversals or excursions (Fig. 3 and SI Appendix, Fig. S9). The saw-tooth inclination oscillations represent asymmetrical interhemispheric polarity drifts at 105.8–103.5 and 98.5–92.0 ka, respectively (Fig. 3). A rapid reversal event within a century at 98 ka (PB2-1) (SI Appendix, Fig. S9) coincides with a global relative paleointensity (RPI) minimum (39⇓–41) (Fig. 3D), which indicates a strong link between the amplitude of the directional deviation and RPI.

A hypothesis for geomagnetic excursions has suggested that rapid field changes (Fig. 3) are associated with an intrinsically oscillating dynamo during periods with low dipole intensity (42), although exact timescales are unknown. Our stalagmite record provides a precisely radioisotopically dated excursion time series, which reveals that directional geomagnetic stability correlated inversely with field intensity (Fig. 3 B and D) and that the geodynamo was intrinsically less stable when dipole field intensities were low. Model simulations (43) indicate that, under a low-friction scenario with small Ekman number, disturbances can propagate in Earth’s liquid outer core. The enhanced outer core convective flow can generate a geomagnetic field with a particular polarity. This gradually increasing field may become unstable and then eventually collapse as depicted by inclination dynamics during PB2-3 and PB2-1. Ongoing convection continues to produce dynamo action and generates magnetic fields that could have the opposite polarity. Inclination variability shown in Fig. 3 suggests that such a nonlinear process probably repeats several times on centennial to millennial timescales in response to chaotic geomagnetic variations. Our results support the argument that the geodynamo was unstable and that it undergoes continuous change to reestablish a full polarity state between occasional collapses (43).

Full geomagnetic polarity reversals are relatively rare (3). The intensity of Earth’s magnetic field varies continuously, and it has decreased by ∼10% over the past century, which has led to suggestions of an impending reversal (11, 12), although it is controversial (44). Our Sanxing stalagmite record indicates that short-lived but rapid field changes occur on centennial timescales when the geomagnetic field intensity decreases into a reversal or excursion. Such rapid geomagnetic changes can disturb animal navigation (45), damage satellites, and disturb radio communications (46). A much-improved view of past geomagnetic fluctuations without the smoothing or temporally aliased spot readings provided by sedimentary or volcanic archives is, thus, important for obtaining additional insights into geomagnetic transition processes.

Methods Summary

Subsampling.

The SX11 stalagmite was cut into halves along its growth axis. One-half of the stalagmite was cut into 8-cm3 cubic subsamples along its central axis for thermal demagnetization of a three-axis IRM. One thin section, 20–30 μm in thickness, was prepared at a depth interval of 53–55 cm along the central stalagmite axis for screening micrometer-scale minerals. Another halved stalagmite was cut into a column with a 4-cm2 cross-section along its central axis perpendicular to horizontal growth laminae. The column was sliced into 194 blocks, with dimensions of 2 (length) × 2 (width) × 0.5 (thickness) cm, containing 5–10% curved laminae for paleomagnetic analyses. After the analyses, 67 block subsamples with 5-mm thickness were selected for U-Th dating.

U-Th Dating and Age Model.

U-Th chemistry was conducted in a class-10,000 metal-free clean room with class-100 benches at the High-Precision Mass Spectrometry and Environment Change Laboratory, Department of Geosciences, National Taiwan University (47). Determinations of all isotopic compositions and concentrations were made on a Thermo-Finnigan NEPTUNE multicollector inductively coupled plasma mass spectrometer with Faraday-cup protocols (26). U-Th uncertainties were calculated at the 2σ level and include corrections for procedure blanks, abundance sensitivity, mass discrimination, and the occurrence of isotopes of interest in spike solution. Half-lives of U-Th nuclides used are given in ref. 26. Details of 230Th age calculation are described in ref. 48. U-Th isotopic compositions and 230Th dates are given in Dataset S1. Agreement of duplicated 230Th ages between coeval subsamples at 10 depths (SI Appendix, Fig. S3) indicates the fidelity of our dating approach. The chronology was built by linear interpolation between 230Th ages (before AD 1950). Ages shown in the figures are the determined 230Th ages with 2σ level (Dataset S1) or were calculated using the linear age model and two adjacent determined ages with propagated errors.

Synchrotron X-Ray Microscope, Rock Magnetism, and Paleomagnetism.

On the prepared thin section, micrometer-scale minerals within calcium carbonate matrix were inspected at high resolution (up to 30 nm) with synchrotron transmission X-ray microscopic techniques (49) using beamlines BL01A1 and BL01B1 at the Taiwan Light Source, National Synchrotron Radiation Research Center (Hsinchu, Taiwan, Republic of China). For thermal demagnetization of a three-axis IRM (27), fields of 2.5, 0.5, and 0.05 T were applied along three orthogonal sample axes on an MC-1 Pulse Magnetizer to identify magnetic minerals using coercivity and unblocking temperatures (Fig. 2 B and D and SI Appendix, Fig. S6).

The stalagmite was not oriented azimuthally; relative declinations were measured with respect to the orientation of the cut half surface. NRM data for the cut block subsamples were obtained using automated stepwise AF demagnetization with a high-resolution cryogenic magnetometer with automatic sample changer (background 2 × 10−7 A/m2) in a magnetically shielded paleomagnetic laboratory at the Institute of Geology and Geophysics, Chinese Academy of Sciences (Beijing, People’s Republic of China). Subsamples were demagnetized with AFs of 2.5, 5.0, 7.5, 10.0, 12.5, 15.0, 17.5, 20, 25, 30, 35, 40, 45, 50, 55, and 60 mT; test subsamples were also treated at additional steps of 70, 80, 90, and 100 mT. Orthogonal vector component plots were used to determine the ChRM for each subsample. A VRM was removed between 10 and 25 mT for most samples. The principal component analysis method of refs. 50 and 51 was used to estimate ChRM directions with 95% confidence limits for each sample. Demagnetization plots for 180 subsamples decay linearly and generally (155 of 180) to the origin. Unanchored fits were used following ref. 50. Fourteen subsamples with erratic demagnetization behavior were excluded. An ARM was applied with a 100-μT direct current bias field and a 100 mT AF. An SIRM1T was imparted in a 1-T field with an MMPM-10 impulse magnetizer. ARM and SIRM1T were only measured for the uppermost 80 subsamples from stalagmite SX11.

The NPSVR for the sampling site was calculated using the TK03 model (30) and the software of ref. 52. A histogram of inclinations from 20,000 simulations in TK03 (SI Appendix, Fig. S8) at the sample site (latitude 27°22′ N) indicates inclination distributions from modeled paleosecular variation for normal (yellow in SI Appendix, Fig. S8) and reversed (blue in SI Appendix, Fig. S8) polarities (where intense colors indicate high probability and white indicates zero probability). Possible VGP paths during the PBE were calculated by rotating relative declination values counterclockwise by 150° to give a mean value of 0° for the nontransitional interval at 103–101 ka (SI Appendix, Fig. S10).

Acknowledgments

We thank D. Heslop of the Australian National University for help with characteristic remanent magnetization calculations with uncertainties. Constructive reviews by Dr. Ron Shaar and two anonymous reviewers significantly improved this study. This research was supported mainly by Science Vanguard Research Program of the Ministry of Science and Technology (MOST) Grant 106-2628-M-002-013 (to C.-C.S.); National Taiwan University Grant 105R7625 (to C.-C.S.); and Higher Education Sprout Project of the Ministry of Education, Taiwan, Republic of China Grant 107L901001 (to C.-C.S.). It was also partially supported by Southern University of Science and Technology Grant Y01316111 (to Y.-M.C.). The diamond wire for sample preparation was supported by MOST Grant 104-2745-8-231-001 (to A.-H.T.). Stalagmite collection and field trips were supported by National Natural Science Foundation of China (NSFC) Grants 41372189 (to X.J.) and 41672170 (to X.J.). Paleomagnetic and rock magnetic analyses were supported by NSFC Grants 41374073 (to Q.L.) and 41430962 (to Q.L.). X-ray image acquisition and analysis were supported by National Synchrotron Radiation Research Center, Taiwan, Republic of China Grant NSRRC-2015-2-059-2 (to Y.-M.C.).

Footnotes

  • ↵1Y.-M.C., X.J., and Q.L. contributed equally to this work.

  • ↵2To whom correspondence may be addressed. Email: river{at}ntu.edu.tw or chouym{at}sustc.edu.cn.
  • Author contributions: C.-C.S. directed this project; Y.-M.C., Q.L., T.-Q.L., and C.-C.S. designed research; Y.-M.C., X.J., H.-M.H., C.-C. Wu, J.L., Z.J., C.-C. Wang, Y.-F.S., C.-C.C., L.T., M.A.L., Y.P., R.Z., Y.H., Y.-C.C., A.-H.T., and C.-C.S. performed research; Y.-M.C., X.J., Q.L., H.-M.H., C.-C. Wu, J.L., Z.J., T.-Q.L., C.-C. Wang, Y.-F.S., C.-C.C., L.T., M.A.L., Y.P., R.Z., Y.H., Y.-C.C., A.-H.T., A.P.R., X.Z., and C.-C.S. analyzed data; and Y.-M.C., X.J., Q.L., H.-M.H., C.-C. Wu, T.-Q.L., A.P.R., X.Z., and C.-C.S. wrote the paper.

  • The authors declare no conflict of interest.

  • This article is a PNAS Direct Submission.

  • This article contains supporting information online at www.pnas.org/lookup/suppl/doi:10.1073/pnas.1720404115/-/DCSupplemental.

Published under the PNAS license.

References

  1. ↵
    1. Buffett BA
    (2000) Earth’s core and the geodynamo. Science 288:2007–2012.
    OpenUrlAbstract/FREE Full Text
  2. ↵
    1. Tarduno JA, et al.
    (2010) Geodynamo, solar wind, and magnetopause 3.4 to 3.45 billion years ago. Science 327:1238–1240.
    OpenUrlAbstract/FREE Full Text
  3. ↵
    1. Channell J,
    2. Kent D,
    3. Lowrie W,
    4. Meert J
    1. Lowrie W,
    2. Kent DV
    (2013) Geomagnetic polarity timescales and reversal frequency regimes. Timescales of the Paleomagnetic Field, eds Channell J, Kent D, Lowrie W, Meert J (American Geophysical Union, Washington, DC), pp 117–129.
  4. ↵
    1. McElhinny MW,
    2. McFadden PL
    (1999) Paleomagnetism: Continents and Oceans (Academic, New York).
  5. ↵
    1. Channell JET,
    2. Hodell DA,
    3. Singer BS,
    4. Xuan C
    (2010) Reconciling astrochronological and 40Ar/39Ar ages for the Matuyama‐Brunhes boundary and late Matuyama Chron. Geochem Geophys Geosyst 11:Q0AA12.
    OpenUrl
  6. ↵
    1. Singer BS
    (2014) A Quaternary geomagnetic instability time scale. Quat Geochronol 21:29–52.
    OpenUrl
  7. ↵
    1. Schubert G
    1. Laj C,
    2. Channell JET
    (2015) Geomagnetic excursions. Treatise on Geophysics, ed Schubert G (Elsevier, Amsterdam), 2nd Ed, Vol 5, pp 343–383.
    OpenUrl
  8. ↵
    1. Roberts AP
    (2008) Geomagnetic excursions: Knowns and unknowns. Geophys Res Lett 35:L17307.
    OpenUrlCrossRef
  9. ↵
    1. Kitaba I,
    2. Hyodo M,
    3. Katoh S,
    4. Dettman DL,
    5. Sato H
    (2013) Midlatitude cooling caused by geomagnetic field minimum during polarity reversal. Proc Natl Acad Sci USA 110:1215–1220.
    OpenUrlAbstract/FREE Full Text
  10. ↵
    1. Maffei ME
    (2014) Magnetic field effects on plant growth, development, and evolution. Front Plant Sci 5:445.
    OpenUrl
  11. ↵
    1. Constable C,
    2. Korte M
    (2006) Is Earth’s magnetic field reversing? Earth Planet Sci Lett 246:1–16.
    OpenUrl
  12. ↵
    1. Laj C,
    2. Kissel C
    (2015) An impending geomagnetic transition? Hints from the past. Front Earth Sci 3:61.
    OpenUrl
  13. ↵
    1. Matuyama M
    (1929) On the direction of magnetisation of basalt in Japan, Tyosen and Manchuria. Proc Imp Acad Jpn 5:203–205.
    OpenUrl
  14. ↵
    1. Opdyke ND,
    2. Channell JET
    (1996) Magnetic Stratigraphy (Academic, New York).
  15. ↵
    1. Clement BM
    (2004) Dependence of the duration of geomagnetic polarity reversals on site latitude. Nature 428:637–640.
    OpenUrlCrossRefPubMed
  16. ↵
    1. Cande SC,
    2. Kent DV
    (1995) Revised calibration of the geomagnetic polarity timescale for the Late Cretaceous and Cenozoic. J Geophys Res 100:6093–6095.
    OpenUrl
  17. ↵
    1. Singer BS,
    2. Guillou H,
    3. Jicha BR,
    4. Zanella E,
    5. Camps P
    (2014) Refining the Quaternary Geomagnetic Instability Time Scale (GITS): Lava flow recordings of the Blake and Post-Blake excursions. Quat Geochronol 21:16–28.
    OpenUrlCrossRef
  18. ↵
    1. van Zijl JSV,
    2. Graham KWT,
    3. Hales AL
    (1962) The palaeomagnetism of the Stormberg Lavas. II. The behaviour of the magnetic field during a reversal. Geophys J R Astron Soc 7:169–182.
    OpenUrlCrossRef
  19. ↵
    1. Channell JET,
    2. Lehman B
    (1997) The last two geomagnetic polarity reversals recorded in high-deposition-rate sediment drifts. Nature 389:712–715.
    OpenUrlCrossRef
  20. ↵
    1. Valet J-P,
    2. Meynadier L,
    3. Simon Q,
    4. Thouveny N
    (2016) When and why sediments fail to record the geomagnetic field during polarity reversals. Earth Planet Sci Lett 453:96–107.
    OpenUrl
  21. ↵
    1. Latham AG,
    2. Schwarcz HP,
    3. Ford DC,
    4. Pearce GW
    (1979) Palaeomagnetism of stalagmite deposits. Nature 280:383–385.
    OpenUrlCrossRef
  22. ↵
    1. Osete M-L, et al.
    (2012) The Blake geomagnetic excursion recorded in a radiometrically dated speleothem. Earth Planet Sci Lett 353–354:173–181.
    OpenUrl
  23. ↵
    1. Lascu I,
    2. Feinberg JM,
    3. Dorale JA,
    4. Cheng H,
    5. Edwards RL
    (2016) Age of the Laschamp excursion determined by U-Th dating of a speleothem geomagnetic record from North America. Geology 44:139–142.
    OpenUrlAbstract/FREE Full Text
  24. ↵
    1. Kirschvink JL,
    2. Kopp RE,
    3. Raub TD,
    4. Baumgartner CT,
    5. Holt JH
    (2008) Rapid, precise, and high‐sensitivity acquisition of paleomagnetic and rock‐magnetic data: Development of a low‐noise automatic sample changing system for superconducting rock magnetometers. Geochem Geophys Geosyst 11:Q05Y01.
    OpenUrl
  25. ↵
    1. Shen C-C, et al.
    (2012) High-precision and high-resolution carbonate 230Th dating by MC-ICP-MS with SEM protocols. Geochim Cosmochim Acta 99:71–86.
    OpenUrlCrossRef
  26. ↵
    1. Cheng H, et al.
    (2013) Improvements in 230Th dating, 230Th and 234U half-life values, and U–Th isotopic measurements by multi-collector inductively coupled plasma mass spectrometry. Earth Planet Sci Lett 371–372:82–91.
    OpenUrlCrossRefPubMed
  27. ↵
    1. Lowrie W
    (1990) Identification of ferromagnetic minerals in a rock by coercivity and unblocking temperature properties. Geophys Res Lett 17:159–162.
    OpenUrlCrossRef
  28. ↵
    1. Strauss BE, et al.
    (2013) The origin of magnetic remanence in stalagmites: Observations from electron microscopy and rock magnetism. Geochem Geophys Geosyst 14:5006–5025.
    OpenUrl
  29. ↵
    1. Vandamme D
    (1994) A new method to determine paleosecular variation. Phys Earth Planet Inter 85:131–142.
    OpenUrlCrossRef
  30. ↵
    1. Channell J,
    2. Kent D,
    3. Lowrie W,
    4. Meert J
    1. Tauxe L,
    2. Kent DV
    (2004) A simplified statistical model for the geomagnetic field and the detection of shallow bias in paleomagnetic tnclinations: Was the ancient magnetic field dDipolar? Timescales of the Paleomagnetic Field, eds Channell J, Kent D, Lowrie W, Meert J (American Geophysical Union, Washington, DC), pp 101–115.
  31. ↵
    1. Thouveny N,
    2. Creer KM,
    3. Blunk I
    (1990) Extension of the Lac du Bouchet palaeomagnetic record over the last 120,000 years. Earth Planet Sci Lett 97:140–161.
    OpenUrl
  32. ↵
    1. Jicha BR, et al.
    (2011) New age for the Skálamælifell excursion and identification of a global geomagnetic event in the late Brunhes chron. Earth Planet Sci Lett 310:509–517.
    OpenUrl
  33. ↵
    1. Valet J-P,
    2. Fournier A,
    3. Courtillot V,
    4. Herrero-Bervera E
    (2012) Dynamical similarity of geomagnetic field reversals. Nature 490:89–93.
    OpenUrlCrossRefPubMed
  34. ↵
    1. Coe RS,
    2. Prevot M,
    3. Camps P
    (1995) New evidence for extraordinarily rapid change of the geomagnetic field during a reversal. Nature 374:687–692.
    OpenUrlCrossRef
  35. ↵
    1. Sagnotti L, et al.
    (2014) Extremely rapid directional change during Matuyama-Brunhes geomagnetic polarity reversal. Geophys J Int 199:1110–1124.
    OpenUrlCrossRef
  36. ↵
    1. Coe RS,
    2. Jarboe NA,
    3. Le Goff M,
    4. Petersen N
    (2014) Demise of the rapid-field-change hypothesis at Steens mountain: The crucial role of continuous thermal demagnetization. Earth Planet Sci Lett 400:302–312.
    OpenUrl
  37. ↵
    1. Valet J-P,
    2. Fournier A
    (2016) Deciphering records of geomagnetic reversals. Rev Geophys 54:410–446.
    OpenUrl
  38. ↵
    1. Evans ME,
    2. Muxworthy AR
    (2018) A re-appraisal of the proposed rapid Matuyama–Brunhes geomagnetic reversal in the Sulmona Basin, Italy. Geophys J Int 213:1744–1750.
    OpenUrl
  39. ↵
    1. Guyodo Y,
    2. Valet J-P
    (1999) Global changes in intensity of the Earth’s magnetic field during the past 800 kyr. Nature 399:249–252.
    OpenUrlCrossRef
  40. ↵
    1. Channell JET,
    2. Xuan C,
    3. Hodell DA
    (2009) Stacking paleointensity and oxygen isotope data for the last 1.5 Myr (PISO-1500). Earth Planet Sci Lett 283:14–23.
    OpenUrl
  41. ↵
    1. Ziegler LB,
    2. Constable CG,
    3. Johnson CL,
    4. Tauxe L
    (2011) PADM2M: A penalized maximum likelihood model of the 0–2 Ma palaeomagnetic axial dipole moment. Geophys J Int 184:1069–1089.
    OpenUrlCrossRef
  42. ↵
    1. Busse FH,
    2. Simitev RD
    (2008) Toroidal flux oscillation as possible cause of geomagnetic excursions and reversals. Phys Earth Planet Inter 168:237–243.
    OpenUrlCrossRef
  43. ↵
    1. Zhang K,
    2. Gubbins D
    (2000) Is the geodynamo process intrinsically unstable? Geophys J Int 140:F1–F4.
    OpenUrlCrossRef
  44. ↵
    1. Brown M,
    2. Korte M,
    3. Holme R,
    4. Wardinski I,
    5. Gunnarson S
    (2018) Earth’s magnetic field is probably not reversing. Proc Natl Acad Sci USA 115:5111–5116.
    OpenUrlAbstract/FREE Full Text
  45. ↵
    1. Tian L-X,
    2. Pan Y-X,
    3. Metzner W,
    4. Zhang J-S,
    5. Zhang B-F
    (2015) Bats respond to very weak magnetic fields. PLoS One 10:e0123205.
    OpenUrlCrossRef
  46. ↵
    1. Stevenson DJ
    (2003) Planetary magnetic fields. Earth Planet Sci Lett 208:1–11.
    OpenUrl
  47. ↵
    1. Shen C-C, et al.
    (2003) Measurement of attogram quantities of 231Pa in dissolved and particulate fractions of seawater by isotope dilution thermal ionization mass spectroscopy. Anal Chem 75:1075–1079.
    OpenUrlCrossRefPubMed
  48. ↵
    1. Shen C-C, et al.
    (2008) Variation of initial 230Th/232Th and limits of high precision U–Th dating of shallow-water corals. Geochim Cosmochim Acta 72:4201–4223.
    OpenUrlCrossRef
  49. ↵
    1. Wang C-C, et al.
    (2015) Evolution and function of dinosaur teeth at ultramicrostructural level revealed using synchrotron transmission X-ray microscopy. Sci Rep 5:15202.
    OpenUrlCrossRefPubMed
  50. ↵
    1. Heslop D,
    2. Roberts AP
    (2016) Analyzing paleomagnetic data: To anchor or not to anchor? J Geophys Res 121:7742–7753.
    OpenUrl
  51. ↵
    1. Heslop D,
    2. Roberts AP
    (2016) Estimation and propagation of uncertainties associated with paleomagnetic directions. J Geophys Res 121:2274–2289.
    OpenUrl
  52. ↵
    1. Tauxe L, et al.
    (2016) PmagPy: Software package for paleomagnetic data analysis and a bridge to the Magnetics Information Consortium (MagIC) Database. Geochem Geophys Geosyst 17:2450–2463.
    OpenUrl
PreviousNext
Back to top
Article Alerts
Email Article

Thank you for your interest in spreading the word on PNAS.

NOTE: We only request your email address so that the person you are recommending the page to knows that you wanted them to see it, and that it is not junk mail. We do not capture any email address.

Enter multiple addresses on separate lines or separate them with commas.
Multidecadally resolved polarity oscillations during a geomagnetic excursion
(Your Name) has sent you a message from PNAS
(Your Name) thought you would like to see the PNAS web site.
CAPTCHA
This question is for testing whether or not you are a human visitor and to prevent automated spam submissions.
Citation Tools
Multidecadally resolved polarity oscillations during a geomagnetic excursion
Yu-Min Chou, Xiuyang Jiang, Qingsong Liu, Hsun-Ming Hu, Chung-Che Wu, Jianxing Liu, Zhaoxia Jiang, Teh-Quei Lee, Chun-Chieh Wang, Yen-Fang Song, Cheng-Cheng Chiang, Liangcheng Tan, Mahjoor A. Lone, Yongxin Pan, Rixiang Zhu, Yaoqi He, Yu-Chen Chou, An-Hung Tan, Andrew P. Roberts, Xiang Zhao, Chuan-Chou Shen
Proceedings of the National Academy of Sciences Sep 2018, 115 (36) 8913-8918; DOI: 10.1073/pnas.1720404115

Citation Manager Formats

  • BibTeX
  • Bookends
  • EasyBib
  • EndNote (tagged)
  • EndNote 8 (xml)
  • Medlars
  • Mendeley
  • Papers
  • RefWorks Tagged
  • Ref Manager
  • RIS
  • Zotero
Request Permissions
Share
Multidecadally resolved polarity oscillations during a geomagnetic excursion
Yu-Min Chou, Xiuyang Jiang, Qingsong Liu, Hsun-Ming Hu, Chung-Che Wu, Jianxing Liu, Zhaoxia Jiang, Teh-Quei Lee, Chun-Chieh Wang, Yen-Fang Song, Cheng-Cheng Chiang, Liangcheng Tan, Mahjoor A. Lone, Yongxin Pan, Rixiang Zhu, Yaoqi He, Yu-Chen Chou, An-Hung Tan, Andrew P. Roberts, Xiang Zhao, Chuan-Chou Shen
Proceedings of the National Academy of Sciences Sep 2018, 115 (36) 8913-8918; DOI: 10.1073/pnas.1720404115
del.icio.us logo Digg logo Reddit logo Twitter logo CiteULike logo Facebook logo Google logo Mendeley logo
  • Tweet Widget
  • Facebook Like
  • Mendeley logo Mendeley

Article Classifications

  • Physical Sciences
  • Earth, Atmospheric, and Planetary Sciences
Proceedings of the National Academy of Sciences: 115 (36)
Table of Contents

Submit

Sign up for Article Alerts

Jump to section

  • Article
    • Abstract
    • U-Th Dating and Paleomagnetic Analyses
    • Asymmetrical Interhemispheric Polarity Oscillations
    • Mechanisms of Asymmetrical Polarity Drifts
    • Methods Summary
    • Acknowledgments
    • Footnotes
    • References
  • Figures & SI
  • Info & Metrics
  • PDF

You May Also be Interested in

Setting sun over a sun-baked dirt landscape
Core Concept: Popular integrated assessment climate policy models have key caveats
Better explicating the strengths and shortcomings of these models will help refine projections and improve transparency in the years ahead.
Image credit: Witsawat.S.
Model of the Amazon forest
News Feature: A sea in the Amazon
Did the Caribbean sweep into the western Amazon millions of years ago, shaping the region’s rich biodiversity?
Image credit: Tacio Cordeiro Bicudo (University of São Paulo, São Paulo, Brazil), Victor Sacek (University of São Paulo, São Paulo, Brazil), and Lucy Reading-Ikkanda (artist).
Syrian archaeological site
Journal Club: In Mesopotamia, early cities may have faltered before climate-driven collapse
Settlements 4,200 years ago may have suffered from overpopulation before drought and lower temperatures ultimately made them unsustainable.
Image credit: Andrea Ricci.
Steamboat Geyser eruption.
Eruption of Steamboat Geyser
Mara Reed and Michael Manga explore why Yellowstone's Steamboat Geyser resumed erupting in 2018.
Listen
Past PodcastsSubscribe
Birds nestling on tree branches
Parent–offspring conflict in songbird fledging
Some songbird parents might improve their own fitness by manipulating their offspring into leaving the nest early, at the cost of fledgling survival, a study finds.
Image credit: Gil Eckrich (photographer).

Similar Articles

Site Logo
Powered by HighWire
  • Submit Manuscript
  • Twitter
  • Facebook
  • RSS Feeds
  • Email Alerts

Articles

  • Current Issue
  • Special Feature Articles – Most Recent
  • List of Issues

PNAS Portals

  • Anthropology
  • Chemistry
  • Classics
  • Front Matter
  • Physics
  • Sustainability Science
  • Teaching Resources

Information

  • Authors
  • Editorial Board
  • Reviewers
  • Subscribers
  • Librarians
  • Press
  • Site Map
  • PNAS Updates
  • FAQs
  • Accessibility Statement
  • Rights & Permissions
  • About
  • Contact

Feedback    Privacy/Legal

Copyright © 2021 National Academy of Sciences. Online ISSN 1091-6490