tmpDE6B TMP
tmpDE6B TMP
tmpDE6B TMP
DOI: 10.1038/ncomms10018
OPEN
The Intertropical Convergence Zone (ITCZ) encompasses the heaviest rain belt on the Earth.
Few direct long-term records, especially in the Pacic, limit our understanding of long-term
natural variability for predicting future ITCZ migration. Here we present a tropical
precipitation record from the Southern Hemisphere covering the past 282,000 years, inferred
from a marine sedimentary sequence collected off the eastern coast of Papua New Guinea.
Unlike the precession paradigm expressed in its East Asian counterpart, our record shows
that the western Pacic ITCZ migration was inuenced by combined precession and obliquity
changes. The obliquity forcing could be primarily delivered by a cross-hemispherical thermal/
pressure contrast, resulting from the asymmetric continental conguration between Asia and
Australia in a coupled East AsianAustralian circulation system. Our nding suggests that the
obliquity forcing may play a more important role in global hydroclimate cycles than previously
thought.
1 CAS
Key Laboratory of Crust-Mantle Material and Environment, School of Earth and Space Science, University of Science and Technology of China,
Hefei 230026, China. 2 High-Precision Mass Spectrometry and Environment Change Laboratory (HISPEC), Department of Geosciences, National Taiwan
University, No. 1, Sec. 4, Roosevelt Road, Taipei 10617, Taiwan ROC. 3 State Key Laboratory of Loess and Quaternary Geology, Institute of Earth Environment,
Chinese Academy of Sciences, Xian 710075, China. 4 CAS Center for Excellence in Tibetan Plateau Earth Sciences, Beijing 100101, China. 5 Department of
Earth Sciences, National Taiwan Normal University, Taipei 11677, Taiwan ROC. 6 Department of Physics, University of Arizona, Tucson, Arizona 85721, USA.
7 Research Center for Environmental Changes, Academia Sinica, Taipei 11529, Taiwan ROC. 8 Department of Geography and Anthropology, Louisiana State
University, Baton Rouge, Louisiana 70803, USA. 9 Department of Earth Sciences, University of Cambridge, Cambridge CB2 3EQ, UK. * These authors
contributed equally to this work. w Present address: Department of Earth Sciences, University of Cambridge, Cambridge CB2 3EQ, UK (L.L.); Japan Agency for
Marine-Earth Science and Technology (JAMSTEC), Yokosuka 2370061, Japan (H.A.). Correspondence and requests for materials should be addressed to C.C.S. (email: river@ntu.edu.tw).
NATURE COMMUNICATIONS | 6:10018 | DOI: 10.1038/ncomms10018 | www.nature.com/naturecommunications
ARTICLE
60N
30N
30S
60S
60E
90E
0
120E
30
150E
60
90
180
120
150W
150
180
120W
210
90W
240
60W
270
300
30W
30E
330
Figure 1 | Precipitation map and study site. Map of mean annual precipitation (cm per year; 19882004; data source: http://jisao.washington.edu/data/
gpcp/). Solid and dashed white lines represent the mean positions of the ITCZ in July and January, respectively. Symbols denote locations of Chinese
caves6,3032 (green circles), marine sediment cores MD05-2925 in this study (blue circle) and 54MC of ref. 23 (white circle), GBR coral24 (white triangle),
PNG coastal coral29 (green triangle), PNG coastal seawater27 (white diamond), surface seawater of the Coral Sea25 (black diamond) and Gregory Lakes38
(orange circle). Simulated precipitation results in sectors of PNG (512S and 130160E, blue solid lines) and North Australia (1220S and 120160E,
blue dashed lines) are given in Fig. 3 and Supplementary Fig. 8.
2
ARTICLE
1.0
1.0
log10 sample/PAAS
1.5
1.5
2.0
PNG G. ruber
PNG G. ruber
2.5
East Pacific
G. sacculifer
2.0
PNG river
2.5
6.0
7.0
3.0
PNG coastal
seawater
PNG coral
GBR coral
3.5
Souhwest pacific
seawater
8.0
4.0
La
Pr Sm Gd Dy Er Yb
Ce Nd Eu Tb Ho Tm Lu
La Pr Sm Gd Dy Er Yb
Ce Nd Eu Tb Ho Tm Lu
La Pr Sm Gd Dy Er Yb
Ce Nd Eu Tb Ho Tm Lu
Figure 2 | REE patterns of waters and marine carbonates. (a) Averaged MD05-2925 G. ruber data with temporal variability range over the past 282 kyr in
this study (grey area). (b) PNG river17. (c) PNG coastal seawater (depth of 40 m at station EUC-Fe 27 from ref. 27, hollow diamonds) and open-ocean
surface seawater of the Coral Sea in the southwest Pacic (depth of 0200 m at station SA-7 from ref. 25, black diamonds). Comparison of REE patterns
from (d) MD05-2925 G. ruber in this study (blue circles), (e) PNG coastal coral29 (green triangles), (f) East Pacic core-top planktonic foraminifera
G. sacculifer (site 54MC of ref. 23; hollow circles), and (g) GBR coral24 (hollow triangles). Site locations are plotted in Fig. 1 and Supplementary Fig. 1.
Arrows depict the trend of the REE patterns. The REEs are shale normalized26.
ARTICLE
560
520
12
440
10
400
8
ASM
6
Modelled PNG
precipitation (mm per day )
PNG precipitation
1.0
0.5
10
0.0
8
520
6
480
440
25
Obliquity ()
G. ruber Nd/Ca
(mol mol1)
1.5
400
24
Chinese stalagmite
18O ()
480
23
22
0
50
100
150
200
250
Discussion
The obliquity effect on SH tropical Pacic precipitation is most
likely associated with its control on the meridional thermalpressure contrast. Modelling results by FOAM suggest that high
obliquity is responsible for the establishment of a strong Siberian
high cell (Supplementary Fig. 8a) and East Asian winter monsoon
system22,33,34. In an experiment using the Geophysical Fluid
Dynamics Laboratory modelling, climate feedbacks and seasonal
response may outcompete the local radiative forcing of obliquity
and induce complicated response of northern high-latitude
climate35. However, FOAM-inferred atmospheric response
(Supplementary Fig. 8) is supported by other simulation results
using Community Climate System Model version 3 (ref. 36;
Supplementary Fig. 10) and Community Earth System version 1
(Supplementary Fig. 11).
Similar with the Siberian high, the Australian low, the
counterpart of the meridional circulation loop, is also affected
by obliquity (Supplementary Fig. 8c). Although, precessional
forcing dominates local land-ocean thermal contrasts and
inuences the Australian low, the simulated Australian low does
not rigidly follow precession. Extreme low pressures are always
induced by high obliquity during the past 282 kyr and
4
ARTICLE
and was recovered in June 2005 during the IMAGES XIII-PECTEN (Past Equatorial Climate: Tracking El Nino) cruise on board the R.V. Marion Dufresne of the
French Polar Institute (IPEV). The core site is located at the southern margin of the
Western Pacic Warm Pool, 110 km to Fergusson Island, 50 km off southeastern
tip of PNG (Fig. 1 and Supplementary Fig. 1).
The core sediment is composed of a mixture of biogenic carbonate and silty
clay40. The chlorophyll level of 0.2 mg m 3 (ref. 41) for surrounding surface water
in eastern PNG suggests low regional productivity. The dissolved-oxygen
concentrations are high (43 ml l 1) through the whole water column including
bottom waters of eastern PNG42. The local benthic oxygen ux, reecting organic
matter remineralization, is only 0.1 mol m 2 per year (ref. 43). It is lower than the
values of 0.8 mol m 2 per year for the reducing margins (notably in the eastern
boundary upwelling systems and North Indian Ocean)43. These data indicate an
oxidative condition at this study site. The upper 1,510 cm was used in this study.
Age model. The age model was established based on accelerator mass spectrometry (AMS) radiocarbon (14C) dates (Supplementary Table 1) and oxygen isotope
stratigraphy (Supplementary Fig. 12). A series of planktonic foraminiferal AMS 14C
dates at 19 different depths, including 200 individuals of Globigerinoides sacculifer
(4500 mm) each, from the upper 292 cm of the core were measured. Dates were
calibrated to calendar ages (before 1950 AD) using CALIB 6.0.1 software44 with a
reservoir age difference (DR) estimated from the Marine Reservoir Correction
Database (http://calib.qub.ac.uk/marine/). The calculated weighted mean DR value
is 6423 years for the selected four sites around the Solomon Sea45. The
chronology was based on linear interpolation between calibrated 14C dates
(Supplementary Table 1).
For the depths 4292 cm, the age model was developed by correlating the
composite benthic foraminiferal oxygen isotopic data of core MD05-2925 to the
LR04 stack record46 (Supplementary Fig. 12). Composite benthic foraminiferal
oxygen isotope data are established with benthic foraminifera (4250 mm, 24
individuals each depth), including the Uvigerina spp. (201 samples), Cibicidoides
wuellerstor (11 samples) and Bulimina spp. (7 samples) at core depths of
1571,897 cm (Supplementary Fig. 13). Measurement of d18O data, relative to
Vienna Pee Dee Belemnite carbonate standard, was performed on a Micromass
IsoPrime isotope ratio mass spectrometer with 1s reproducibility of 0.05%
(ref. 47). d18O offsets of C. wuellerstor ( 0.64%)48 and Bulimina spp.
( 0.11%)49 from Uvigerina spp. were corrected. This age model is supported
by the last occurrence of G. ruber (pink) occurred at depths of
830835 cm, corresponding to 129.8 kyr BP (Supplementary Fig. 12), consistent
with the observation in the southern South China Sea50.
Screening for diagenesis. Scanning electron microscopy images of 30 uncleaned
individuals of planktonic foraminifera G. ruber (white, s.s. 250300 mm) at six
depths of 477 (50.1 kyr BP) and 617 cm (81.6 kyr BP) with low REE content, 527
(56.8 kyr BP) and 577 cm (73.2 kry BP) with high REE content and 877 (135.0 kyr
BP) and 917 cm (146.1 kyr BP) with moderate REE content (Supplementary
Fig. 14) were carefully screened. Thirty more uncleaned individuals picked from six
depths (87, 267, 787, 1,087, 1,317 and 1,477 cm), respectively, at marine isotope
stages 1, 2, 5, 6, 7 and 8 were also checked with scanning electron microscopy. No
nodules of Mn-Fe oxides were noticeable and all shell walls were intact and
primitive (Supplementary Fig. 14). Additional careful inspection under microscope
did not observe Mn-Fe oxides for 1,200 tests from the selected 12 depths. For
conservative consideration, we still applied a full cleaning procedure on all samples.
Measurement of foraminiferal trace elements. REE contents of down-core
planktonic foraminifera G. ruber (white, s.s. 250300 mm) were measured
(Supplementary Fig. 2). Although no Mn-Fe nodules were noticeable
(Supplementary Fig. 14), G. ruber tests were cleaned with a full cleaning procedure
for foraminiferal trace metal analysis, modied from refs 51,52. About 20
foraminiferal individuals were gently crushed, placed in a Teon vial and washed
sequentially with the following reagents (all at pH 8.59.0): (i) ethanol H2O,
ARTICLE
(ii) 1% H2O2, (iii) 0.56 M NH4Cl and (iv) 0.43 M NH2OH. Cleaned tests, polished
with 10 3 M HNO3 to dissolve a possible thin post-depositional magnesium-rich
surface layer51,52, were rinsed with ultrapure water three times to wash off the
residues of chemicals and then dissolved in 5% HNO3 for instrumental analyses.
All chemical procedures were performed on a class-100 laminar-ow bench in a
class-10,000 clean room in the High-precision Mass Spectrometry and
Environment Change Laboratory (HISPEC), Department of Geosciences, National
Taiwan University.
REE/Ca ratios were calculated using the ion beams of 46Ca, 139La, 140Ce, 141Pr,
146Nd, 147Sm, 153Eu, 160Gd, 159Tb, 163Dy, 165Ho, 166Er, 172Yb and 175Lu, detected
on an ICP-SF-MS, Thermo Fisher ELEMENT II, equipped with a dry introduction
Cetac ARIDUS20 system. Two-month 2s reproducibility is 1.96.5%. Mg/Ca,
Mn/Ca and Fe/Ca ratios with respective 2s errors of 0.23%, 0.68% and
2.7% were determined on the same ICP-SF-MS, equipped with a quartz Scotttype double-pass spray chamber53.
An insignicant correlation between Mg/Ca and Fe/Ca data (Supplementary
Fig. 3a) indicates the effectiveness of the cleaning techniques. Moreover, the
measured REE/Ca patterns (Fig. 2) are different from shale-like patterns for
uncleaned foraminifera with greater light REE (LREE) contents enrichment and
unclear Ce anomalies54. We also tested our cleaning procedure/analytical technique
by an interlaboratory comparison for analysing REE/Ca ratios of benthic
foraminifera C. wuellerstor sample from core GGC-15 (ref. 20). The results
showed that our REE data replicate measurements using a REE cleaning method at
the University of Cambridge (Fig. 5 of ref. 20). Detailed instrumentation and
delity of our methodology for foraminiferal test REE/Ca determination are
described in ref. 20.
Nd isotopic measurement. Planktonic foraminifera G. ruber and sediment
(o63 mm) samples were collected from two depth intervals of 472477 cm
(49.550.1 kyr BP, 580 individuals, 4250 mm) and 537542 cm (58.860.6 kyr BP,
250 individuals, 4250 mm) of core MD05-2925 (Supplementary Fig. 2). The picked
planktonic foraminifera samples were cleaned with the same protocol for REE/Ca
ratio analysis and then dissolved in 2 M HNO3. The sediment samples were rst
cleaned with 10% CH3COOH to remove carbonate, and subsequently cleaned with
a reductive reagent (1 M NH2OH HCl in 25% CH3COOH) to remove possible
Fe-Mn phases on the sample surface55. The cleaned sediment samples were
decomposed in a mixed solution of HF, HClO4 and HNO3, and then dissolved
in 2 M HNO3.
Neodymium in the 2 M HNO3 dissolved samples was extracted by a two-stage
column separation56. The REE fraction in the solution was puried from the
remaining major and trace elements using Eichrom RE resin. Neodymium was
subsequently separated from the other REE with Eichrom Ln resin.
Neodymium isotopic compositions were measured by a multi-collector
ICP-MS, Thermo Fisher Neptune, in the HISPEC. The measured 143Nd/144Nd
ratios were normalized to 146Nd/144Nd 0.7219 using an exponential law.
La Jolla standard was measured at 0.5118110.000014 (or 0.27 e; 2s, n 13).
All 143Nd/144Nd ratios were calibrated to the reported value relative to the
La Jolla standard value of 0.511858 (ref. 57). Sample 143Nd/144Nd ratios
[(143Nd/144Nd)sample] are expressed as e notation dened by an equation
of eNd [(143Nd/144Nd)sample/(143Nd/144Nd)CHUR 1] 104, where the
143Nd/144Nd ratio of CHUR standard for Chondritic Uniform Reservoir
[(143Nd/144Nd)CHUR] is 0.512638 (ref. 58).
Modelling simulation. The simulated precipitation and climatological records
used in this study are from an orbital-accelerated transient run using FOAM
conducted by Kutzbach et al.21 and re-analysed by Shi et al.22. FOAM, a fully
coupled, mixed-resolution, and high-throughput general circulation model,
provides a good simulation of mean condition and variability59. With a factor of
100, FOAM was integrated for 2,820 years under orbital forcing only to obtain
climate evolution over the past 282 kyr. Changes in global ice volume/sea level and
greenhouse gases were not considered. The spatial resolution is set to 4 7.5 for
atmosphere and 1.4 2.8 for ocean. Because of the limitation of orbital
acceleration, it is difcult for the deep ocean to reach equilibrium so that the full
potential of the deep ocean feedback cannot be achieved. However, in previous
studies21,22,60, the responses of monsoon precipitation, mostly considered as a
response to the changes in the atmosphere-surface ocean system, to the orbital
insolation can be successfully retrieved in the annual variability. A detailed
description on the transient experiment is available in ref. 21.
References
1. Waliser, D. E. & Gautier, C. Satellite-derived climatology of the ITCZ. J. Clim 6,
21622174 (1993).
2. Cai, W. et al. More extreme swings of the South Pacic convergence zone due
to greenhouse warming. Nature 488, 365369 (2012).
3. Haug, G. H. et al. Climate and the collapse of Maya Civilization. Science 299,
17311735 (2003).
4. Yancheva, G. et al. Inuence of the intertropical convergence zone on the East
Asian monsoon. Nature 445, 7477 (2007).
6
ARTICLE
Acknowledgements
This study used samples offered by the Marine Core Repository and Laboratory, Taiwan
Ocean Research Institute (TORI), Taiwan ROC. The location of MD05-2925 site was
suggested by M.-T. Chen and M.-Y. Lee, and collected during the IMAGES PECTEN
Cruise, led by L. Beaufort and M.-T. Chen. We thank the scientic party of WEPAMA,
MD122/IMAGES VII for coring. A. Chiang helped pick foraminiferal specimens.
Community Climate System Model version 3 data were downloaded from East Asian
Paleoenvironmental Science Database (http://paleo-data.ieecas.cn/en/). Constructive and
comprehensive reviews by ve anonymous reviewers over the past two years signicantly
improved this paper. We are thankful for nancial support provided by Taiwan ROC
MOST (952611-M-002-019, 962611-M-002-019, 992611-M-002-005, 1002116-M002-009, 1012611-M-002-003, 1022611-M-001-006, 1032119-M-002-022, 1032611M-001-004, and 1042119-M-002-003), National Taiwan University (101R7625),
National Key Basic Research Program of China (2013CB956102 and 2013CB955904),
State Key Laboratory of Loess and Quaternary Geology, Institute of Earth Environment,
CAS (SKLLQG1411), National Science Foundation of China (41572148 and 41290255).
Data of the reported planktonic foraminifera G. ruber records of trace metals and oxygen
stable isotopes and benthic foraminiferal d18O over the past 282 kyr are available in the
Supplementary Materials.
Author contributions
C.-C.S. directed the project; C.-C.S., Y.L., K.-Y.W., and Z.P. designed the research.
Y.L., L.L., C.-J.C., C.-C.W. and H.A. performed foraminiferal test REE/Ca and Nd
isotope analyses. H.-S.M., Y.-C.C., L.L. and C.-K.C. did d18O analyses. Z.S. and
S.-Y.L. was responsible for model simulations and G.S.B. for AMS 14C dating.
Y.L., C.-C.S. and L.L. prepared the draft and all authors contributed to the completion of
manuscript.
Additional information
Supplementary Information accompanies this paper at http://www.nature.com/
naturecommunications
Competing nancial interests: The authors declare no competing nancial interests.
Reprints and permission information is available online at http://npg.nature.com/
reprintsandpermissions/
How to cite this article: Liu, Y. et al. Obliquity pacing of the western Pacic Intertropical
Convergence Zone over the past 282,000 years. Nat. Commun. 6:10018
doi: 10.1038/ncomms10018 (2015).
This work is licensed under a Creative Commons Attribution 4.0
International License. The images or other third party material in this
article are included in the articles Creative Commons license, unless indicated otherwise
in the credit line; if the material is not included under the Creative Commons license,
users will need to obtain permission from the license holder to reproduce the material.
To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/