ARTICLE IN PRESS
JQSR2725_proof ■ 20 April 2010 ■ 1/13
Quaternary Science Reviews xxx (2010) 1e13
Contents lists available at ScienceDirect
Quaternary Science Reviews
journal homepage: www.elsevier.com/locate/quascirev
OF
Holocene climate trend, variability, and shift documented by lacustrine stable
isotope record in the northeastern United States
Cheng Zhao a, *, Zicheng Yu a, Emi Ito b, Yan Zhao c
a
Department of Earth and Environmental Sciences, Lehigh University, Bethlehem, PA, 18015, USA
Limnological Research Center, Department of Geology and Geophysics, University of Minnesota, 310 Pillsbury Drive SE, Minneapolis, MN 55455, USA
c
MOE Key Laboratory of Western China’s Environmental System, College of Earth and Environmental Sciences, Lanzhou University, Lanzhou 730000, China
RO
b
a b s t r a c t
Article history:
Received 6 April 2009
Received in revised form
1 March 2010
Accepted 31 March 2010
Earlier studies indicated that the general pattern of the Holocene climate in the northeastern United
States changed from cool and dry (11.6e8.2 ka; 1 ka ¼ 1000 cal yr BP) to warm and wet (8.2e5.4 ka) to
warm and dry (5.4e3 ka) to cool and wet (after 3 ka). A new w35-year resolution stable isotope record of
endogenic calcite from a sediment core for Lake Grinnell in northern New Jersey provided a chance to
examine the Holocene climate variations of the region in a finer detail. After the Younger Dryas cold
climate reversal, the d18O fluctuated around a constant value of 7.4& until 5.8 ka, thereafter shifted to
a steadily decreasing trend to the most recent value of 8.2&. Responding to this shift, the widely
observed hemlock decline in the northeastern USA occurred about w350e500 (143.5) years later.
Detrended d18O and d13C records show a clear covariance at 910-year periodicity. The amplitudes of
centennial-scale d18O variations became much smaller after 4.7 ka. At the same time, the dominant
frequency of these variations changed from 330 to 500 years. We suggest that a non-linear response of
atmospheric circulation to the gradual decrease in insolation is responsible for the shift in the climate
trend at 5.8 ka as indicated by the deceasing d18O values. A dominant frequency shift in solar forcing and
the decreased seasonal contrast of insolation might have caused the change in climate variability at
4.7 ka through modulating ocean and atmosphere circulations.
Ó 2009 Elsevier Ltd. All rights reserved.
1. Introduction
CO
RR
EC
TE
DP
a r t i c l e i n f o
Stable isotopes from lacustrine carbonate can be a useful tool for
reconstructing long-term climate changes and investigating
possible causes and mechanisms of these climate changes. In the
northeastern United States, many efforts have been made to
reconstruct the Holocene climate history from lake sediments,
including lithologically derived lake-level reconstructions (e.g.
Dwyer et al., 1996; Mullins, 1998; Newby et al., 2000; Shuman et al.,
2004; Dieffenbacher-Krall and Nurse, 2005), pollen-inferred vegetation changes (e.g. Davis,1969; Watts, 1979; Davis et al.,1980; Suter,
1985; Maenza-Gmelch, 1997; Shuman et al., 2001), as well as stable
isotope records (e.g. Anderson et al., 1997; Huang et al., 2002; Kirby
et al., 2002; Ellis et al., 2004; Hou et al., 2006, 2007). This Holocene
climate history includes a cool and dry period at 11.6e8.2 ka
(1 ka ¼ 1000 cal year BP), a warm and wet climate at 8.2e5.4 ka,
a warm and dry interval at 5.4e3 ka, and a cool and wet condition
during the past 3 ka (Shuman et al., 2004). Superimposed on this
general trend, a series of millennial-, centennial- and multidecadal-
UN
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
* Corresponding author.
E-mail address: [email protected] (C. Zhao).
scale climate oscillations have been reported (e.g. Kirby et al., 2002;
Cronin et al., 2005; Viau et al., 2006; Li et al., 2007; Springer et al.,
2008). In addition, the mid-Holocene was a period of large climate
change. For examples, many lakes in eastern North America experienced a low lake-level from 5 to 3 ka (e.g. Dwyer et al., 1996; Yu
et al., 1997; Newby et al., 2000; Lavoie and Richard, 2000;
Almquist et al., 2001; Mullins and Halfman, 2001; Shuman et al.,
2001, 2004; Dieffenbacher-Krall and Nurse, 2005), which has been
attributed to broad-scale changes in atmosphere circulation (e.g. Yu
et al., 1997; Kirby et al., 2002). In response to the dry climate, an
abrupt decline of the hemlock (Tsuga canadensis) around 5.5 ka has
been documented from pollen data across eastern North America (e.
g. Allison et al., 1986; Foster and Zebryk, 1993; Fuller, 1998; Bennett
and Fuller, 2002; Foster et al., 2006) and the decline has recently
been suggested to be caused by increased aridity (e.g. Yu et al., 1997;
Shuman et al., 2004; Foster et al., 2006). However, there are only
a few stable isotope records with temporal resolution sufficient for
reconstructing the multidecadal- to centennial-scale Holocene
climate history (e.g. Kirby et al., 2002) and examining the possible
mechanisms of climate changes.
In this paper, we present an w35-year-resolution calcite stable
isotope record of a marl sediment core collected from a small,
0277-3791/$ e see front matter Ó 2009 Elsevier Ltd. All rights reserved.
doi:10.1016/j.quascirev.2010.03.018
Please cite this article in press as: Zhao, C., et al., Holocene climate trend, variability, and shift documented by lacustrine stable isotope...,
Quaternary Science Reviews (2010), doi:10.1016/j.quascirev.2010.03.018
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
ARTICLE IN PRESS
2
C. Zhao et al. / Quaternary Science Reviews xxx (2010) 1e13
hydrologically open hardwater lake in northern New Jersey, USA.
The objectives of this study are to reconstruct a detailed climate
history covering the last w12,000 years with a special focus on the
nature of the proposed mid-Holocene climate shift and to identify
and understand the possible forcing factors and how the climate
system responded to those forcings.
OO
PR
CO
RR
E
CT
ED
Lake Grinnell (41060 N, 74 380 W, 170 m above sea level) is
located in Sussex County in northern New Jersey (Fig. 1A). The lake
sits in a limestone terrain and is situated in the Kittatinny Valley of
the Ridge and Valley physiographic province, which is underlain
by northeast-trending folded and faulted blocks of dolostone,
limestone, slate and siltstone (Witte and Monteverde, 2006). The
Laurentide Ice Sheet retreated from this area around 14,000 years
ago (Witte, 2001). The regional climate has been strongly affected
by the moisture from the Gulf of Mexico and tropical Atlantic Ocean
(Lawrence et al., 1982; Peixoto and Oort, 1983). Based on the
instrumental data (1895e2006) from New Jersey Climate Division 1
(northern New Jersey) (http://www1.ncdc.noaa.gov/pub/data/cirs),
the regional mean annual temperature is 10.4 C with an average
F
2. Study region and site
summer (JuneeAugust) temperature of 21.5 C and an average
winter (DecembereFebruary) temperature of 1.1 C (Fig. 2A). The
mean annual precipitation is 1177 mm with almost uniform
monthly distribution (Fig. 2A). The average of monthly weighted
mean precipitation d18O values (collected from 5 stations in the
Global Network of Isotopes in Precipitation (GNIP), including Ste.
Agathe, Quebec; Ottawa, Ontario; Simcoe, Ontario; Coshocton,
Ohio; and Chicago, Illinois (http://isohis.iaea.org/userupdate/
GNIPMonthly.xls)) in eastern North America (see locations in
Fig. 1A) shows a high correlation with the mean monthly temperature in northern New Jersey (Fig. 2A)
Lake Grinnell is a small hardwater lake of w0.2 km2 in surface
area (Fig. 1B) and w7 km2 in catchment area. It has two basins
separated by a ridge where water is <3 m deep. The lake has a short
mean residence time of w80 days and is mostly recharged by
shallow groundwater and perhaps precipitation on the catchment
basin (Laura Nicholson, personal communication). The water depth
in the basins was raised about 2.4 m artificially by a dam to form the
present impoundment area, with a maximum water depth of about
10.4 m. The date of the first dam construction is not well known
except that it occurred before 1913 (Valgenti et al., 1950). Lake
Grinnell is the headwater lake of a branch of Wallkill River, which
UN
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
JQSR2725_proof ■ 20 April 2010 ■ 2/13
Fig. 1. Location and settings. (A) Map showing the location of our study site Lake Grinnell (red square) in northern New Jersey. Other locations (black dots): 1, Queechy Lake, New
York; 2, Green Lake, New York; 3, Inglesby Lake, Ontario; 4, Crawford Lake, Ontario; and 5, Pretty Lake, Indiana. Also shown the locations of five stations in the Global Network of
Isotopes in Precipitation (GNIP) in eastern North America (blue triangles): a, Ste. Agathe, Quebec; b, Ottawa, Ontario; c, Simcoe, Ontario; d, Coshocton, Ohio; and e, Chicago, Illinois.
(B) Air photo of the Lake Grinnell area. White dot represents the coring location of core GL05-1. White squares indicate the locations of water samples collected: a, Lake center;
b, Boat deck; c, Travis’ and Naby’s wells; d, Sincaglia’s well; e, Lake outlet (see Table 1 for detail). The water depth contours are in meters. (C) Core photo of core GL05-1. The first
drive is from 660 to 760 cm below the lake surface, and the last drive (Drive 9) is from 1460 to 1560 cm.
Please cite this article in press as: Zhao, C., et al., Holocene climate trend, variability, and shift documented by lacustrine stable isotope...,
Quaternary Science Reviews (2010), doi:10.1016/j.quascirev.2010.03.018
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
233
234
235
236
237
238
239
240
ARTICLE IN PRESS
JQSR2725_proof ■ 20 April 2010 ■ 3/13
C. Zhao et al. / Quaternary Science Reviews xxx (2010) 1e13
monitoring around White Lake from fall 1992 through spring 1994
showed that water table responded quickly to precipitation events,
suggesting local recharge of shallow aquifer located in the limestone (Nicholson, 1995). Given the high permeability of limestone
aquifers and hardwater composition of the lake, the hydrologic
budget of the lake likely includes significant groundwater component as well as surface inflow.
3. Methods
PR
OO
F
Water samples at 50 cm water depth from different locations in
Lake Grinnell (Fig. 1B) were collected mostly from 2006 to 2008
(Table 1). Several domestic wells around the lake were sampled in
2006 and 2008 by first pumping out about 1.5 L of water before
collecting in 500 ml bottles. One shallow well (Travis’ well) was
sampled in both years. Samples for hydrogen and oxygen isotopes
were sealed in 20 ml high density polyethylene plastic bottles with
scintillation caps right after collection. Samples for dissolved inorganic carbon (DIC) were filtered through a 0.45-micron filter paper
and ‘poisoned’ by adding two drops of CuSO4 solution to kill
organisms and sealed in Amber DIC glass bottles. All water and DIC
samples were sent to the University of Arizona Isotope Lab for
dDwater and d18Owater and d13CDIC analyses. Isotope values are
reported in conventional d notation relative to VSMOW (Vienna
Standard Mean Ocean Water) for water dD and d18O, and VPDB
(Vienna Pee Dee Belemnite) for DIC d13C. The analytical precisions
are 0.09& for d18O, 0.9& for dD, and 0.3& for d13C.
The 900-cm long sediment core (GL05-1; Fig. 1C) was taken on
30 January 2005 from near the center of the lake (see Fig. 1B for
coring location) at 660 cm water depth with a Livingstone-Wright
piston corer of 5 cm inside diameter. Sub-samples were taken
directly from the core at every 2 cm, and then homogenized for
loss-on-ignition (LOI), stable isotope and pollen analyses. LOI
analysis was carried out to estimate organic matter content based
on weight loss from combustion at 550 C about 1 h and carbonate
content at 1000 C about 1 h (Dean, 1974). Terrestrial plant
macrofossils and charcoals were picked from seven intervals and
dated using an accelerator mass spectrometry (AMS) at the
University of California, Irvine (Table 2). All the 14C dates were
calibrated to calendar ages using CALIB 5.01 program based on the
IntCal04 dataset (Reimer et al., 2004). The age model is derived
from a fourth-polynomial curve with a fixed surface age (55 cal
year BP ¼2005 AD; Fig. 3A).
Samples for calcite stable isotope analyses were air-dried at
room temperature. Plant macrofossils, mollusk and ostracode
shells, and other macroscopic fragments were removed and discarded under a stereo-microscope. Each of the 334 calcite samples
CT
ED
Fig. 2. Regional climate and water isotope data. (A) Mean monthly temperature and
precipitation between 1895 and 2006 from New Jersey Climate Division 1 (northern
New Jersey), showing a large seasonal temperature change but almost uniform
monthly precipitation. The monthly averages of weighted precipitation d18O values
from the five GNIP stations (average of 5 stations) in eastern North America. (B) Stable
isotope values of lake water and nearby wells at Lake Grinnell from different seasons,
plotting close to the global meteoric water line (GMWL; dD¼8d18O þ 10). The annual
weighted mean precipitation isotope data from 5 GNIP stations in eastern North
America were also included.
CO
RR
E
flows northeastward to Hudson River. The water chemistry
measurements made in 2002 show a pH of 8.9, an electrical
conductivity of 503.6 mS/cm, and total dissolved solids of 0.4 g/L.
The alkalinity of the lake water was 0.15 g/L in 1950 (Valgenti et al.,
1950). The lake bottom is almost completely covered with musk
grass Chara, an aquatic multicelluar algae that is common in
hardwater lakes (Valgenti et al., 1950).
Due to its location within the Kittatinny Valley, the narrow low
ridges on either side of Lake Grinnell are about 30e40 m higher
than the lake surface in elevation. White Lake (Sussex County,
Fig. 1B), about 500 m to the southwest of Lake Grinnell, is about 7 m
higher in its lake surface elevation. There is no channel connecting
White Lake to Lake Grinnell today but surface or groundwater
connection between the two lakes is possible. Groundwater table
Table 1
Water isotope measurements from and around Lake Grinnell, New Jersey.a
UN
241
242
243
244
245
246
247
248
249
250
251
252
253
254
255
256
257
258
259
260
261
262
263
264
265
266
267
268
269
270
271
272
273
274
275
276
277
278
279
280
281
282
283
284
285
286
287
288
289
290
291
292
293
294
295
296
297
298
299
300
301
302
303
304
305
3
Sample ID
Collection date (month/year)
dD (&, VSMOW)
d18O (&, VSMOW)
020623C
060211D
060618a
060618b
061010a
12302006
080517A
080517B
080908C
080908D
060618c
060618d
080517C
080908A
06/2002
02/2006
06/2006
06/2006
10/2006
12/2006
05/2008
05/2008
09/2008
09/2008
06/2006
06/2006
05/2008
09/2008
31
51
48
44
43
50
45
47
43
41
53
51
53
53
5.2
7.9
6.8
6.9
7.3
7.7
7.0
7.0
5.9
5.9
8.1
8.2
8.2
8.2
a
b
d13C of DIC (&, VPDB)
8.5
7.3
7.1
12.0
Sampling locationb
Boat deck
Boat deck
Boat deck
Lake center
Lake outlet
Lake outlet
Lake outlet
Boat deck
Boat deck
Lake outlet
Travis’ well (w1.5 m)
Naby’s well (w1.5 m)
Sincaglia’s well (w30 m)
Travis’ well (w1.5 m)
The analytical precisions are 0.09& for d18O, 0.9& for dD, and 0.3& for d13C as reported from the University of Arizona Isotope Lab.
See Fig. 1B for sampling location.
Please cite this article in press as: Zhao, C., et al., Holocene climate trend, variability, and shift documented by lacustrine stable isotope...,
Quaternary Science Reviews (2010), doi:10.1016/j.quascirev.2010.03.018
306
307
308
309
310
311
312
313
314
315
316
317
318
319
320
321
322
323
324
325
326
327
328
329
330
331
332
333
334
335
336
337
338
339
340
341
342
343
344
345
346
347
348
349
350
351
352
353
354
355
356
357
358
359
360
361
362
363
364
365
366
367
368
369
370
ARTICLE IN PRESS
4
C. Zhao et al. / Quaternary Science Reviews xxx (2010) 1e13
Table 2
Radiocarbon (14C) dates from Lake Grinnell, New Jersey.
Depth (cm, below
lake surface)
AMS lab IDa
14
C date error
(yr BP)
Median age
(cal BP 2s range)b
Materials dated
720e724
762e764
898e902
1022e1024
1150e1154
1242e1244
19482
16918
19483
16915
16914
16919
1330 30
1430 15
3495 15
5340 20
7625 30
9475 20
1265.5 36.5
1324 24
3769 63
6082.5 80.5
8417 41
10715.5 59.5
1278e1282
16921
9870 45
11,283 80
2 Larix needles, 1 Cyperaceae seed, plant macrofossils and charcoals
7 woody scales from Abies or Picea, plant macrofossils
plant macrofossils and charcoals
1 Scirpus seed, 1 woody scale, plant macrofossils and charcoals
1 woody scale, plant macrofossils and charcoals
1 Pinus seed wing, 4 Pinus needles, 5 pieces Pinus seed fragments
and plant macrofossils
8 Pinus needles, plant macrofossils
F
using the program AnalySeries (version 2.0.4), and calculated the
rates of changes (expressed as 100 D(d18O)/Dt) between each
adjacent pair of raw measurements. The raw d18O and d13C data
were fitted and detrended with a 3rd order-polynomial curve using
the program AutoSignal (version 1.0) to remove the long-term
trend. In order to evaluate the potential solar forcing on the
centennial-scale climate variability and to highlight these centennial-scale variations recorded in d18O time series, we also
performed spectral analysis using the program REDFIT (Schulz and
Mudelsee, 2002) using a Hanning window and band-pass filters in
the detrended d18O and raw residual D14C time series (Reimer et al.,
2004). The program AnalySeries was also used to perform crossspectral analysis of the detrended d18O and d13C time series in order
to evaluate their coherence, and to apply a band-pass filter (at
910 80 year periodicity) to show their corresponding changes at
millennial-scale.
CO
RR
E
CT
ED
was analyzed for oxygen and carbon isotopes in the Stable Isotope
Laboratory at the University of Minnesota using a Finnigan MAT
252 isotope ratio mass spectrometer coupled to a Kiel P carbonate
preparation device. All results are reported in standard d notation
relative to VPDB. The analytical precision is 0.06& for both d18O
and d13C. To evaluate the regional vegetation response to climate
changes, pollen analysis was done on 0.7 cm3 sub-samples at every
8-cm interval (Zhao et al., in press). The samples were prepared
using a modified standard acetolysis procedure (see details in Zhao
et al., in press).
The computer program RAMPFIT (Mudelsee, 2000) was used to
analyze raw d18O data to objectively detect the shift in the longterm trend. In order to visualize and evaluate the magnitude of
centennial-scale variability of the d18O time series, we obtained the
deviations of raw d18O values from the RAMPFIT trend, filtered the
raw d18O time series after applying a 200e500 year band-pass filter
OO
b
Sample lab ID for Keck-Carbon Cycle AMS facility at the University of California, Irvine.
Calibrated ages for radiocarbon dates were obtained based on the IntCal04 dataset using CALIB 5.01 (Reimer et al., 2004).
PR
a
UN
371
372
373
374
375
376
377
378
379
380
381
382
383
384
385
386
387
388
389
390
391
392
393
394
395
396
397
398
399
400
401
402
403
404
405
406
407
408
409
410
411
412
413
414
415
416
417
418
419
420
421
422
423
424
425
426
427
428
429
430
431
432
433
434
435
JQSR2725_proof ■ 20 April 2010 ■ 4/13
Fig. 3. Results from core GL05-1 at Lake Grinnell, New Jersey. (A) Age model based on seven calibrated 14C dates using a fourth-polynomial curve. (B) Organic matter content (%,
dotted curve) and carbonate content (%, solid curve) from LOI analysis. (C) d18O values from calcite. (D) d13C values from calcite.
Please cite this article in press as: Zhao, C., et al., Holocene climate trend, variability, and shift documented by lacustrine stable isotope...,
Quaternary Science Reviews (2010), doi:10.1016/j.quascirev.2010.03.018
436
437
438
439
440
441
442
443
444
445
446
447
448
449
450
451
452
453
454
455
456
457
458
459
460
461
462
463
464
465
466
467
468
469
470
471
472
473
474
475
476
477
478
479
480
481
482
483
484
485
486
487
488
489
490
491
492
493
494
495
496
497
498
499
500
ARTICLE IN PRESS
JQSR2725_proof ■ 20 April 2010 ■ 5/13
C. Zhao et al. / Quaternary Science Reviews xxx (2010) 1e13
OO
PR
5. Discussion
5.1. Interpretation of calcite stable isotopes at Lake Grinnell
ED
The water isotope values of the lake and surrounding wells are
within the range of amount-weighted isotope values in annual
precipitation from the five GNIP stations in eastern North America
(Fig. 1A) (http://isohis.iaea.org/userupdate/GNIPMonthly.xls). Also,
these water samples scatter around the global meteoric water line
(GMWL, Fig. 2B). The lake water DIC d13C values are higher (average
7.6&) than the groundwater d13C value of 12& (Table 1).
The 900-cm long core has about 230 cm clay sediment at the
bottom (1560e1330 cm below the lake surface). The rest of the core
is mostly marl (1330e670 cm), except for 10 cm of organic-rich
gyttja at the very top of the core (670e660 cm; Fig. 1C). This study
focuses only on the marl section. The marl contains >90% carbonate
and w6% organic matter (Fig. 3B). The seven calibrated 14C dates
from the core GL05-1 are all in order. The chronology indicates
the marl section covers the last 12500 calibrated years (Fig. 3A). The
accumulation rate of the marl is almost constant throughout the
core at w0.54 mm/year based on the age model.
The d18O values range from 8.9& to 6.7&, with centennialscale fluctuations (Fig. 3C). A major negative excursion at
12.5e11.3 ka reached a minimum of 8.9&. After 11.3 ka, the d18O
values fluctuate around a constant value of 7.40.4& until 5.8 ka,
thereafter steadily decreasing to most recent value of 8.2&
(Fig. 4A), which suggests an equilibrium precipitation of calcite
from present-day lake water. The high d18O values of about 6.7 to
6.8& occur at 11.1, 10.2, 9.5, 6e5.7 and 4.8 ka. The d18O deviations
from RAMPFIT data show w0.7& variations before and <0.5&
variations after 4.7 ka (Fig. 4B). The smoothed time series after
a 200e500 year band-pass filter also show w0.2 to 0.3& variations
before and w0.1& variations after 4.7 ka (Fig. 4C). The D(d18O)/Dt
results indicate highest rates of changes around 4.7 ka (Fig. 4D).
In addition, spectral analysis of the detrended d18O data shows
a 500-year periodicity above 95% confidence level (CL) after 4.7 ka
(Fig. 5C) and a 330-year periodicity above 95% CL before 4.7 ka
(Fig. 5D).
The d13C values range from 9.0& to 4.2& (Fig. 3D). A major
positive excursion between 12.5 and 11.3 ka reached 6& at
w11.7 ka. The d13C values were steady around 7.4& between 11.3
and 9 ka and thereafter steadily increased to 4.5& at 8.5 ka
(Fig. 3D). From 8.2 to 5.5 ka, d13C values averaged 6.0 0.4&.
From 5.5 to 4.7 ka, d13C increases from 6 to 5& with little
fluctuation, followed by a decreasing trend to 7& at 3.5 ka. From
3 ka to the present, d13C values remained nearly constant at
5.0 1&. The cross-spectral analysis of the detrended d18O and
d13C values (Fig. 6A, B) shows highest coherency at w0.8 (Fig. 6C)
and in-phase (w0 ) relationship (Fig. 6D) at the 910-year periodicity (although both curves have a very weak correlation, with an R
value of w0.05), which is also supported by the covariance of the
smoothed detrended d18O and d13C time series after applying
a 910 80 year band-pass filter (Fig. 6E).
F
4. Results
CO
RR
E
CT
Changes in calcite d18O values obtained at Lake Grinnell likely
reflect changes in climate conditions in northern New Jersey, rather
UN
501
502
503
504
505
506
507
508
509
510
511
512
513
514
515
516
517
518
519
520
521
522
523
524
525
526
527
528
529
530
531
532
533
534
535
536
537
538
539
540
541
542
543
544
545
546
547
548
549
550
551
552
553
554
555
556
557
558
559
560
561
562
563
564
565
5
Fig. 4. Time series analyses of calcite d18O time series from Lake Grinnell, New Jersey. (A) Raw d18O values with RAMPFIT line to show change in the long-term trend. (B) d18O
deviations from long-term fitted RAMPFIT line (detrended record). (C) Smoothed curve of raw d18O time series after a 200e500 year band-pass filter to show low-frequency
components. (D) The rate of change (100 D(d18O)/Dt) based on raw d18O time series.
Please cite this article in press as: Zhao, C., et al., Holocene climate trend, variability, and shift documented by lacustrine stable isotope...,
Quaternary Science Reviews (2010), doi:10.1016/j.quascirev.2010.03.018
566
567
568
569
570
571
572
573
574
575
576
577
578
579
580
581
582
583
584
585
586
587
588
589
590
591
592
593
594
595
596
597
598
599
600
601
602
603
604
605
606
607
608
609
610
611
612
613
614
615
616
617
618
619
620
621
622
623
624
625
626
627
628
629
630
ARTICLE IN PRESS
OO
PR
ED
CT
Fig. 5. Spectral analyses of the detrended d18O time series (by a third-polynomial curve) at Lake Grinnell and raw residual D14C time series (data from Reimer et al., 2004). (A)
Spectral power of residual D14C after 4.7 ka. (B) Spectral power of residual D14C before 4.7 ka. (C) Spectral power of detrended d18O after 4.7 ka. (D) Spectral power of detrended d18O
before 4.7 ka. The legends of spectral power, spectral background and different confidence levels (CL) were also indicated.
CO
RR
E
than a local hydrogeological behavior for the following reasons: (1)
the lake has a very short-residence time of about 80 days. Thus, the
lake water would be quickly replaced by inflowing ground and
surface water. (2) Groundwater table monitoring around White
Lake (Sussex County), 500 m south of Lake Grinnell within the
same valley, from fall 1992 to spring 1994 showed that the water
table responded quickly to precipitation events (Nicholson, 1995).
(3) Chara encrustations (see details below), the calcite materials
that we used for d18O analyses, mostly precipitate from late spring
to summer. Thus they only capture the summer water and climatic
signal. (4) Our measurements of summer lake water show higher
d18O values when comparing to well water and winter lake water
(Fig. 2B, Table 1), suggesting that precipitation contributes significantly to the lake water during summers. Therefore, we argue that
calcite d18O values at Lake Grinnell likely reflect the summer
conditions in the region.
The fine-grained calcite that we analyzed is comprised of mostly
disintegrated Chara encrustations as well as some endogenic
calcite. Chara uses dissolved CO2 from the lake water for photosynthesis, which is facilitated significantly by Chara calcification
(McConnaughey, 1991; McConnaughey and Falk, 1991), inducing
calcite precipitation as a result. The d18O values of Chara encrustations have been documented to be close to the expected equilibrium values (Coletta et al., 2001; Ito, 2002). Our stable isotope
measurements of Lake Grinnell water are consistent with prior
results. At our study region, the average warm season (Maye
September) temperature between 2006 and 2008 is about 21 C,
with a maximum mean monthly temperature of w27 C and
a minimum temperature of w15 C. The temperature of Chara
photosynthetic activity is probably close to the average warm
UN
631
632
633
634
635
636
637
638
639
640
641
642
643
644
645
646
647
648
649
650
651
652
653
654
655
656
657
658
659
660
661
662
663
664
665
666
667
668
669
670
671
672
673
674
675
676
677
678
679
680
681
682
683
684
685
686
687
688
689
690
691
692
693
694
695
C. Zhao et al. / Quaternary Science Reviews xxx (2010) 1e13
F
6
JQSR2725_proof ■ 20 April 2010 ■ 6/13
season temperature, considering the high heat capacity of water
assuming that the groundwater inflow mixes quickly with the lake
water. For an average temperature of 21 C and average warm
season d18Owater values of 6.6& (VSMOW) (see Table 1), the
calculated d18Ocalcite at equilibrium is 8.1& (VPDB) using the
fractionation equation published in Leng and Marshall (2004) as
modified from Kim and O’Neil (1997). The uppermost d18Ocalcite
value from the core is 8.2& (VPDB), which is very close to the
calculated value of 8.1&.
The d18O values of calcite in isotopic equilibrium reflect isotopic
composition and temperature of the lake water during calcite
precipitation. The water in Lake Grinnell is derived from groundwater and precipitation, with groundwater table responding
quickly to precipitation events suggesting a relatively shallow
aquifer and local recharge (see details in Section 2). The d18O values
of meteoric precipitation are controlled by moisture sources and air
temperature. Evaporative enrichment can also increase the lake
water d18O values. The lake/groundwater isotope values and
regional precipitation isotopic compositions fall on the GMWL
(Fig. 2B), suggesting that this short-residence time lake has
a minimal evaporation (Leng and Marshall, 2004 and references
therein). Therefore, it is unlikely that the changes in precipitation/
evaporation ratio control the d18Ocalcite signal. Changes in the
moisture source of precipitation could also influence the d18O
values of lake water. For example, Burnett et al. (2004) investigated
the moisture sources of winter precipitation in New York state and
concluded that the lake-effect precipitation had a much lower d18O
values (17.9&) than the moisture originating in the Gulf of Mexico
and in the nearby Atlantic Ocean (8.2&). However, the averaged
d18O values of annual precipitation in northern New Jersey is about
Please cite this article in press as: Zhao, C., et al., Holocene climate trend, variability, and shift documented by lacustrine stable isotope...,
Quaternary Science Reviews (2010), doi:10.1016/j.quascirev.2010.03.018
696
697
698
699
700
701
702
703
704
705
706
707
708
709
710
711
712
713
714
715
716
717
718
719
720
721
722
723
724
725
726
727
728
729
730
731
732
733
734
735
736
737
738
739
740
741
742
743
744
745
746
747
748
749
750
751
752
753
754
755
756
757
758
759
760
ARTICLE IN PRESS
JQSR2725_proof ■ 20 April 2010 ■ 7/13
OO
PR
ED
CT
Fig. 6. Covariance of the detrended d18O and d13C time series from Lake Grinnell. (A) Detrended d18O time series by its third-polynomial fitting curve. (B) Detrended d13C time series
by its third-polynomial fitting curve. (C) The spectral coherence between the detrended d18O and d13C time series. (D) The phase relationship between the detrended d18O and d13C
time series. (E) Smoothed curves of detrended d18O (solid) and d13C time series (dashed) after a 910 80 year band-pass filter.
CO
RR
E
8& (Bowen and Wilkinson, 2002). The shortest distance from the
Laurentide Great Lakes (Lake Erie) to Lake Grinnell is w300 km and
they are separated by the Allegheny Mountains which could act as
a sufficient orographic barrier against any lake-effect precipitation
from reaching northern New Jersey. Thus the major moisture source
in northern New Jersey is more likely from the Gulf of Mexico and
tropical Atlantic Ocean (Lawrence et al., 1982; Peixoto and Oort,
1983). Changes in lake water temperature of this small, shortresident-time lake likely follow the air temperature. Hence changes
in air temperature are the most likely cause of changes in d18O
values of Chara encrustations at Lake Grinnell. As the average of
monthly weighted mean precipitation d18O values in eastern North
America show a close correlation with the monthly air temperature
in northern New Jersey (Fig. 2A). The annual weighted mean
precipitation d18O values and air temperatures (http://isohis.iaea.
org/userupdate/GNIP2001yearly.xls) from the five GNIP stations
(Fig. 1A) are also correlated (d18Oprecipitation ¼ 0.7 Tair 14.9;
R2 ¼ 0.70) for the period 1962e2001. Considering both the positive
relationship between meteoric water d18O and air temperature at
0.7 & per C in eastern North America, and the negative relationship with water temperature during calcite precipitation at 0.24
& per C (Craig, 1965), we can use a simple carbonate d18O-air
temperature relation of 0.46 & per C as a first approximation to
estimate air temperature changes.
The d13C values of Chara encrustation depend mainly on local
factors, particularly the d13C values of DIC in lake waters (Ito, 2002).
The d13C of lake water DIC could be influenced by the atmospherelake exchange of CO2, the DIC isotope composition of groundwater
inflow, the decomposed or respired organic matter input into the
lake water, and especially the photosynthetic activity by Chara
(McKenzie, 1985; Cole et al., 1994; Ito, 2002; Leng and Marshall,
UN
761
762
763
764
765
766
767
768
769
770
771
772
773
774
775
776
777
778
779
780
781
782
783
784
785
786
787
788
789
790
791
792
793
794
795
796
797
798
799
800
801
802
803
804
805
806
807
808
809
810
811
812
813
814
815
816
817
818
819
820
821
822
823
824
825
7
F
C. Zhao et al. / Quaternary Science Reviews xxx (2010) 1e13
2004). The d13C values were unlikely influenced by interaction
with atmospheric CO2 because of the lake’s small size, short water
residence time and the relatively constant atmospheric CO2
concentration during the Holocene (with a range of about 25 ppmv;
Indermühle et al., 1999). The groundwater DIC d13C value of 12&
(Table 1) reflects about equal contributions from both limestone
bedrock (w0&) and C3 plants in this region (from about 32 to
20&, Leng and Marshall, 2004 and references therein). The
decomposed and respired organic matter could decrease the DIC
d13C value of lake water by releasing 12C from organic matter. Chara
will increase the d13C of the DIC by preferentially incorporating 12C
in photosynthesis (McConnaughey et al., 1997). The lake water DIC
d13C values range from 7 to 8.5& (Table 1), which are higher
than that of groundwater DIC, likely to be the result of extensive
photosynthetic activity by Chara that covers the lake floor.
5.2. Holocene climate history inferred from calcite isotope records
The negative excursion of w1.5 & in d18O at 12.511.3 ka
represents a 3e4 C cooling in air temperature, corresponding to
the Younger Dryas (YD) cold event (e.g. Cwynar and Levesque,
1995; Shuman et al., 2002; Yu, 2007). The 3e4 C temperature
decrease during the YD inferred at Lake Grinnell is comparable to
other temperature records in the northeastern United States,
including a similar cooling inferred from pollen records in New
Jersey and Connecticut (Peteet et al., 1993), the w5 C cooling
derived from d18O values in carbonate at White Lake (Warren
County), New Jersey (Yu, 2007), the 5.6 C cooling at Blood Pond in
Massachusetts (Hou et al., 2007), and the more than 5 C cooling
at Crooked Pond in Massachusetts (Huang et al., 2002) estimated
from dD values in palmitic acid. This cold period was locally
Please cite this article in press as: Zhao, C., et al., Holocene climate trend, variability, and shift documented by lacustrine stable isotope...,
Quaternary Science Reviews (2010), doi:10.1016/j.quascirev.2010.03.018
826
827
828
829
830
831
832
833
834
835
836
837
838
839
840
841
842
843
844
845
846
847
848
849
850
851
852
853
854
855
856
857
858
859
860
861
862
863
864
865
866
867
868
869
870
871
872
873
874
875
876
877
878
879
880
881
882
883
884
885
886
887
888
889
890
ARTICLE IN PRESS
haven’t been well documented in other paleotemperature records
in the northeastern United States. Some of the IRD events have been
documented by the moisture reconstructions for the past 7000
years from this region (Li et al., 2007; Springer et al., 2008). Thus,
we suggest that the IRD events in North Atlantic Ocean only have
influences on the moisture condition rather than the temperature
in the northeastern United States.
5.3. Climate shift in the mid-Holocene
OO
F
The d18O value shift at mid-Holocene seems to be a regional
signal and has been revealed by other d18O records from lake
sediment carbonate in eastern North America (Fig. 7). In New York,
the d18O values from Queechy Lake (Stuiver, 1970) and Green Lake
(Kirby et al., 2002) show a negative shift around 6e5 ka. In
southern Ontario, the d18O values from Crawford Lake (Yu et al.,
1997) and Inglesby Lake (Edwards and Fritz, 1988) also show
a negative shift around 5 ka. Similar negative shift was recorded in
the d18O values from Pretty Lake in Indiana (Stuiver, 1968).
Although the rate and specific timing of the shift are different from
site to site, all these shifts seem to imply a consistent climate
transition in eastern North America during the mid-Holocene. In
addition, a mid-Holocene decrease in dD values was found from
Crooked Pond in Massachusetts (Huang et al., 2002). We suggest
that the transition is likely caused by the same climate driver and
occurred at same time considering the poor dating and coarse
resolution for some of these published records.
The d18O deviations to the RAMPFIT output (Fig. 4B), the
smoothed time series by the band-pass filter (Fig. 4C), and spectral
analysis results of d18O data (Fig. 5C, D) indicate that the isotopic
record of Lake Grinnell had higher frequency and larger amplitude
variations during the early to mid-Holocene (before 4.7 ka) than
the late Holocene (after 4.7 ka; Figs. 4, 5). The change in the
amplitude and/or frequency of climate variability during the
Holocene has been reported in paleoclimate archives in North
America. For example, temperature reconstructions from pollen
data in eastern North America showed that the amplitude of
millennial-scale oscillations has been decreasing during the Holocene (Viau et al., 2006). Similarly, a clear change in variability (both
amplitude and frequency) at 5 ka was shown in the sea surface
temperature (SST) data obtained from the subpolar North Atlantic
(Berner et al., 2008). Decreased amplitude of millennial-scale
oscillations since the mid-Holocene in the western Canada was also
suggested by a model of climate fluctuations extended back to 15 ka
based on a 4000 year paleoclimate record in conjunction with the
Milankovitch cycle (Campbell et al., 1998). In addition, the rate of
d18O changes (Fig. 4D) indicated that the most rapid climate variations happened in the mid-Holocene around 4.7 ka.
The mid-Holocene hemlock decline was thought to be a rapid
vegetation response to climate change (Yu et al., 1997; Shuman
et al., 2004; Foster et al., 2006). Close examination of the difference in the timing of the shift in d18O and the decline in hemlock
pollen (Zhao et al., in press) at Lake Grinnell (both pollen and
isotope data from the same samples) (Fig. 8), indicates there is
a delay between the hemlock decline and the of climate change.
The climate shift occurred at 5.8 ka (at 1014 cm below the lake
surface), while the hemlock decline started between w5.45 and
5.3 ka (at 988e996 cm below the lake surface). We have radiocarbon ages from above and below the hemlock decline. The dating
error for both ages yields a maximum error of 143.5 year (also see
Table 2). The age difference between the d18O shift and hemlock
decline, based on our age-depth model, indicates the hemlock
decline lagged w350e500 (143.5) years behind the climate shift.
Similar age difference is obtained if we use linear interpolation to
construct the age model. If the mid-Holocene hemlock decline was
CO
RR
E
CT
ED
characterized by increased abundance of boreal taxa such as Pinus,
Picea, Betula, and Alnus (Zhao et al., in press). The positive excursion
of w3& in d13C at YD may reflect the decreased decomposition and
respiration rate of organic matter in the cold climate. The timing of
the YD event recorded at Lake Grinnell is 400e500 years younger
than other regional records (e.g., Shuman et al., 2002), probably due
to the age extrapolation in our age model from the oldest dating
control in the early Holocene.
From 11.3 to 5.8 ka, d18O values remain relatively high around
7.4&. For the period before 9 ka, other climate reconstructions
from eastern North America show a cool but gradually increasing
temperature (Suter, 1985; Anderson et al., 1997; Huang et al., 2002;
Viau et al., 2006) and a dry climate (e.g. Dwyer et al., 1996; Newby
et al., 2000; Lavoie and Richard, 2000; Almquist et al., 2001;
Shuman et al., 2001, 2004; Dieffenbacher-Krall and Nurse, 2005).
The cool temperatures would have caused lower precipitation d18O
values to Lake Grinnell region, but the dry climate might cause
higher lake water d18O values due to longer water residence time
leading to stronger evaporation effect than present. This interval
was also characterized by a vegetation transition from a pine forest
to an oak forest, and a decline of hemlock (Zhao et al., in press),
consistent with reduced moisture availability. The low d13C values
at 11.3e9 ka could be caused by the less limestone contributions to
the groundwater DIC induced by reduced groundwater recharge
under a dry climate.
A warm early to mid-Holocene climate indicated by the stable
d18O values from 9 to 5.8 ka is in agreement with many studies in
the northeastern United States and around the North Atlantic
region (e.g. Stuiver, 1970; Davis et al., 1980; Suter, 1985; Huang
et al., 2002; Ellis et al., 2004; McFadden et al., 2004; Hou et al.,
2006). At Crooked Pond in Massachusetts, dD values indicated an
annual temperature rise by 3e7 C between 8.2 and 5.4 ka (Huang
et al., 2002). Similar magnitude temperature rise was also reconstructed from the pollen-inferred January mean temperature at
Winneconnet Pond in Massachusetts (Suter, 1985). At Lake Grinnell,
this period was also characterized by the stabilization of oak forest,
with a high hemlock population (Zhao et al., in press). The relatively
high and stable d13C values could reflect the high photosynthesis
activity by Chara.
Starting at 5.8 ka, d18O values shifted from a stable state to
a gradual and steady decreasing trend reflecting a w2 C cooling
since the mid-Holocene if this shift were controlled by temperature
alone. A dry climate between 5 and 3 ka was reported from most
lake sediment records in eastern North America (e.g. Dwyer et al.,
1996; Mullins, 1998; Yu et al., 1997; Newby et al., 2000; Lavoie
and Richard, 2000; Almquist et al., 2001; Mullins and Halfman,
2001; Shuman et al., 2001, 2004; Dieffenbacher-Krall and Nurse,
2005). This dry interval corresponds with a period of low d13C
values centered at 3.5 ka from Lake Grinnell, again indicating less
limestone contributions to the groundwater DIC under the dry
climate similar to the situation at 11.3e9 ka.
The carbon-oxygen isotopic covariance at 910-year periodicity
(Fig. 6E) likely reflects the positive controls of air temperature on
both d13C and d18O values in temperate lakes, as suggested by
Drummond et al. (1995). Higher mean annual temperature will
cause higher precipitation d18O values, which in turn leads to
higher d18O of groundwater discharging into Lake Grinnell.
Assuming plenty of sunlight is available, the rate of photosynthesis
increases with temperature so that residual DIC becomes enriched
in 13C. The lack of such a covariance at shorter timescales might be
caused by the more complex and localized factors controlling the
carbonate d13C values.
Our data don’t show distinct climate change around 8.2 ka or
other millennial-scale ice-rafted debris (IRD) events as documented
in the North Atlantic (Bond et al., 2001). Also those IRD events
UN
891
892
893
894
895
896
897
898
899
900
901
902
903
904
905
906
907
908
909
910
911
912
913
914
915
916
917
918
919
920
921
922
923
924
925
926
927
928
929
930
931
932
933
934
935
936
937
938
939
940
941
942
943
944
945
946
947
948
949
950
951
952
953
954
955
C. Zhao et al. / Quaternary Science Reviews xxx (2010) 1e13
PR
8
JQSR2725_proof ■ 20 April 2010 ■ 8/13
Please cite this article in press as: Zhao, C., et al., Holocene climate trend, variability, and shift documented by lacustrine stable isotope...,
Quaternary Science Reviews (2010), doi:10.1016/j.quascirev.2010.03.018
956
957
958
959
960
961
962
963
964
965
966
967
968
969
970
971
972
973
974
975
976
977
978
979
980
981
982
983
984
985
986
987
988
989
990
991
992
993
994
995
996
997
998
999
1000
1001
1002
1003
1004
1005
1006
1007
1008
1009
1010
1011
1012
1013
1014
1015
1016
1017
1018
1019
1020
ARTICLE IN PRESS
JQSR2725_proof ■ 20 April 2010 ■ 9/13
OO
PR
ED
Fig. 7. Regional d18O records from lake sediments in eastern North America: (A) Lake Grinnell, New Jersey, with RAMPFIT line; (B) Queechy Lake, New York (Stuiver, 1970); (C) Green
Lake, New York (Kirby et al., 2002); (D) Crawford Lake, Ontario (Yu et al., 1997); (E) Inglesby Lake, Ontario (Edwards and Fritz, 1988); and (F) Pretty Lake, Indiana (Stuiver, 1968).
CT
really caused by climate change, our data suggest that the decline
took w350e500 years to respond. The w350e500 year lag is
similar to the suggestion from other sites that the hemlock decline
occurred >300 years after climate change (Shuman et al., 2009).
RR
E
5.4. Possible causes of Holocene climate trend, variability and shift
CO
Changes in insolation have been previously considered as
driving long-term climate trends in the northeastern United States
(Kirby et al., 2002) (Fig. 9A, B, and C). The presence of the remaining
Laurentide Ice Sheet (LIS; Fig. 9A) also affects the climate system.
For example, the cool and dry early Holocene was probably caused
by the effect of cold glacial anticyclone circulation from the
remaining LIS (COHMAP, 1988). The glacial anticyclone creates an
anomalous north-northeasterly flow on the south-east of the ice
sheet (Kutzbach et al., 1998; Webb et al., 1998). Thus, cloudiness
and precipitation south of the ice sheet were suppressed (Felzer
et al., 1996; Webb et al., 1998). This mechanism would be responsible for the cool and dry early Holocene condition in the northeastern United States (Shuman and Donnelly, 2006). An analysis of
spatial and temporal expression of early Holocene climate around
the North Atlantic region shows that locations close to LIS experienced maximum Holocene temperatures that lagged peak summer
insolation about 1000e3000 years (Kaplan and Wolfe, 2006). In
addition, reconstructed SSTs along the northwestern Atlantic slope
waters show a cool early Holocene close to the LIS but a warm early
Holocene farther south (Sachs, 2007). After the retreat of LIS,
insolation became the dominant control (Fig. 9A).
The mid-Holocene climate shift at 5.8 ka recorded at Lake
Grinnell (Fig. 9D) likely reflect a non-linear response to the gradual
decrease of insolation within the climate system. General circulation models (GCM) have shown that subtle changes in insolation,
aided with the positive feedbacks between atmosphere, ocean, sea
ice, and vegetation, could produce significant climate responses
UN
1021
1022
1023
1024
1025
1026
1027
1028
1029
1030
1031
1032
1033
1034
1035
1036
1037
1038
1039
1040
1041
1042
1043
1044
1045
1046
1047
1048
1049
1050
1051
1052
1053
1054
1055
1056
1057
1058
1059
1060
1061
1062
1063
1064
1065
1066
1067
1068
1069
1070
1071
1072
1073
1074
1075
1076
1077
1078
1079
1080
1081
1082
1083
1084
1085
9
F
C. Zhao et al. / Quaternary Science Reviews xxx (2010) 1e13
(e.g. Foley et al., 1994; Ganopolski et al., 1998; Claussen et al., 1999;
Kerwin et al., 1999). As summarized by Steig (1999) and Sandweiss
et al. (1999), the mid-Holocene climate shift has been documented
in various paleoclimate records and in cultural changes at global
scale. In addition, the fastest rate of annual insolation decline
(Fig. 9C) could promote a more unstable climate during the midHolocene through rapidly changing thermal gradient at different
latitudes. The insolation decrease could cause a significant atmosphere circulation change in eastern North America. Several
mechanisms of the atmosphere circulation change have been
proposed for the decreased precipitation d18O values observed in
eastern North America (Fig. 7), including the more frequent intrusions of dry Pacific air mass with low d18O values (Yu et al., 1997),
the increased cold season precipitation (Denniston et al., 1999;
Shuman and Donnelly, 2006), and the higher altitude of the
moisture condensation elevated by the cooler and drier western air
(Smith and Hollander, 1999). The contribution of Pacific air mass is
unlikely the cause of the d18O shift at Lake Grinnell because it can
bring very little moisture even into the continental interior region
like Iowa (Simpkins, 1995). In the case of the latter two mechanisms, the average air temperature would be lower at Lake Grinnell
region and thus decrease the d18O values.
The frequency of centennial-scale climate oscillations recorded
in Lake Grinnell sediment is possibly caused by a change in dominant frequency of solar variations. The results of spectral analysis
show that the 330-year periodicity was significant before 4.7 ka, but
the 500-year periodicity was significant after 4.7 ka in both d18O and
D14C time series (Fig. 5). In order to better show the relationship of
these two curves at the dominant frequencies, we smoothed both
time series after applying a 330 (30)-year band-pass filter before
4.7 ka and a 500 (50)-year filter after 4.7 ka on both the detrended
d18O and raw residual D14C time series. Both filtered curves seem to
show corresponding variations, especially after 6 ka (Fig. 9E). The
poor correlation in some intervals during the early Holocene may be
Please cite this article in press as: Zhao, C., et al., Holocene climate trend, variability, and shift documented by lacustrine stable isotope...,
Quaternary Science Reviews (2010), doi:10.1016/j.quascirev.2010.03.018
1086
1087
1088
1089
1090
1091
1092
1093
1094
1095
1096
1097
1098
1099
1100
1101
1102
1103
1104
1105
1106
1107
1108
1109
1110
1111
1112
1113
1114
1115
1116
1117
1118
1119
1120
1121
1122
1123
1124
1125
1126
1127
1128
1129
1130
1131
1132
1133
1134
1135
1136
1137
1138
1139
1140
1141
1142
1143
1144
1145
1146
1147
1148
1149
1150
ARTICLE IN PRESS
OO
PR
ED
CO
RR
E
CT
Fig. 8. Timing of mid-Holocene climate shift and hemlock decline. (A) Age model with two dating points adjacent to the hemlock decline. Dating error is indicated by the error bars,
with age models based on the fourth-polynomial curve used in this paper (solid) and on a linear interpolation (dashed). (B) The d18O values from Lake Grinnell with RAMPFIT line.
The gray band shows the climate shift point. (C) Hemlock pollen percentage from the same core as isotope data at Lake Grinnell (Zhao et al., in press). The gray band indicates the
hemlock decline between two pollen samples.
UN
1151
1152
1153
1154
1155
1156
1157
1158
1159
1160
1161
1162
1163
1164
1165
1166
1167
1168
1169
1170
1171
1172
1173
1174
1175
1176
1177
1178
1179
1180
1181
1182
1183
1184
1185
1186
1187
1188
1189
1190
1191
1192
1193
1194
1195
1196
1197
1198
1199
1200
1201
1202
1203
1204
1205
1206
1207
1208
1209
1210
1211
1212
1213
1214
1215
C. Zhao et al. / Quaternary Science Reviews xxx (2010) 1e13
F
10
JQSR2725_proof ■ 20 April 2010 ■ 10/13
Fig. 9. Regional climate records and climate controls. (A) Annual insolation at 40 N (Berger and Loutre, 1991) and change in Laurentide Ice Sheet (LIS) areas compared to full glacial
period (gray area). Half-way point of insolation changes is indicated by the dashed line. (B) Summer (solid) and winter (dashed) insolation at 40 N (Berger and Loutre, 1991). Halfway point of summer insolation changes is indicated by the dashed line. (C) The rate of change in annual insolation per 1000 year (kyr) at 40 N. The gray bar indicated the most
rapid decrease in annual insolation. (D) The d18O values from Lake Grinnell with RAMPFIT line. The seven calibrated 14C dates and errors are indicated on the left. (E) Smoothed
curves of detrended d18O (solid curve) and raw residual D14C time series (dashed curve; Reimer et al., 2004) after a 330 30 year band-pass filter before 4.7 ka and a 500 50 year
band-pass filter after 4.7 ka.
Please cite this article in press as: Zhao, C., et al., Holocene climate trend, variability, and shift documented by lacustrine stable isotope...,
Quaternary Science Reviews (2010), doi:10.1016/j.quascirev.2010.03.018
1216
1217
1218
1219
1220
1221
1222
1223
1224
1225
1226
1227
1228
1229
1230
1231
1232
1233
1234
1235
1236
1237
1238
1239
1240
1241
1242
1243
1244
1245
1246
1247
1248
1249
1250
1251
1252
1253
1254
1255
1256
1257
1258
1259
1260
1261
1262
1263
1264
1265
1266
1267
1268
1269
1270
1271
1272
1273
1274
1275
1276
1277
1278
1279
1280
ARTICLE IN PRESS
JQSR2725_proof ■ 20 April 2010 ■ 11/13
C. Zhao et al. / Quaternary Science Reviews xxx (2010) 1e13
1346
1347
1348
1349
1350
1351
1352
1353
1354
1355
1356
1357
1358
Acknowledgements
1359
1360
We thank Robert Moeller, Yongxiang Li, Long Li and EES 357
1361
Spring 2005 class at Lehigh University for field coring assistance,
1362
John Naby and his neighbors around Lake Grinnell for water
1363
sampling assistance, and two anonymous reviewers for helpful
1364
comments that improved the manuscript. This project was
1365
supported by American Chemical Society e Petroleum Research
1366
Fund (#42418-AC2) and US NSF (EAR 0518774) to Z.C. Yu. This is
1367
Contribution 09-22 of the Limnological Research Center of the
1368
University of Minnesota.
1369
1370
References
Q1 1371
1372
Allison, T.D., Moeller, R.E., Davis, M.B., 1986. Pollen in laminated sediments provides
1373
evidence for a mid-Holocene forest pathogen outbreak. Ecology 67, 1101e1105.
1374
Almquist, H., Dieffenbacher-Krall, A.C., Flanagan-Brown, R., Sanger, D., 2001. The
Holocene record of lake levels of Mansell Pond, central Maine, USA. The
1375
Holocene 11, 189e201.
1376
Anderson, W.T., Mullins, H.T., Ito, E., 1997. Stable isotope record from Seneca Lake,
1377
New York: evidence for a cold paleoclimate following the Younger Dryas.
1378
Geology 25, 135e138.
Asmerom, Y., Polyak, V., Burns, S., Rassmussen, J., 2007. Solar forcing of Holocene
1379
climate: new insights from a speleothem record, southwestern United States.
1380
Geology 35, 1e4.
1381
Bennett, K.D., Fuller, J.L., 2002. Determining the age of the mid-Holocene Tsuga
canadensis (hemlock) decline, eastern North America. The Holocene 12,
1382
421e429.
1383
Berger, A., Loutre, M.F., 1991. Insolation values for the climate of the last 10 million
1384
of years. Quaternary Science Reviews 10, 297e317.
Berner, K.S., Koc, N., Divine, D., Godtliebsen, F., Moros, M., 2008. A decadal-scale
1385
Holocene sea surface temperature record from the subpolar North Atlantic
1386
constructed using diatoms and statistics and its relation to other climate
1387
parameters. Paleoceanography 23. doi:10.1029/2006PA001339.
Bond, G., Kromer, B., Beer, J., Muscheler, R., Evans, M.N., Showers, W., Hoffmann, S.,
1388
Lotti-Bond, R., Hajdas, I., Bonani, G., 2001. Persistent solar influence on North
1389
Atlantic climate during the Holocene. Science 294, 2130e2136.
1390
Bowen, G.J., Wilkinson, B., 2002. Spatial distribution of d18O in meteoric precipitation. Geology 30, 315e318.
1391
Burnett, A.W., Mullins, H.T., Patterson, W.P., 2004. Relationship between atmo1392
18
spheric circulation and winter precipitation d O in central New York State.
1393
Geophysical Research Letters 31, L22209. doi:10.1029/2004GL021089.
1394
Campbell, I., Campbell, C., Apps, M.J., Rutter, N.W., Bush, A.B.G., 1998. Late Holocene
w1500 yr climatic periodicities and their implications. Geology 26, 471e473.
1395
Chapman, M.R., Shackleton, N.J., 2000. Evidence of 550-year and 1000-year
1396
cyclicities in North Atlantic circulation patterns during the Holocene. The
1397
Holocene 10, 287e291.
Claussen, M., Kubatzki, C., Brovkin, V., Ganpolski, A., 1999. Simulation of an abrupt
1398
change in Saharan vegetation in the mid-Holocene. Geophysical Research
1399
Letters 26, 2037e2040.
1400
COHMAP Members, 1988. Climatic changes of the last 18,000 years: observations
and model simulations. Science 241, 1043e1052.
1401
Cole, J.J., Caraco, N.F., Kling, G.W., Kratz, T.K., 1994. Carbon dioxide supersaturation
1402
in surface waters of lakes. Science 265, 1568e1570.
1403
Coletta, P., Pentecost, A., Spiro, B., 2001. Stable isotopes in charophyte incrustations:
relationships with climate and water chemistry. Palaeogeography, Palae1404
oclimatology, Palaeoecology 173, 9e19.
1405
Craig, H., 1965. The measurement of oxygen isotope paleotemperatures. In:
1406
Tongiorgi, E. (Ed.), Stable Isotopes in Oceanographic Studies and Paleotemperatures. Consiglio Nazionale delle Ricerche Laboratorio di Gelogia Nuclear, Pisa,
1407
pp. 161e182.
1408
Cronin, T.M., Thunell, R., Dwyer, G.S., Saenger, C., Mann, M.E., Vann, C., Seal II, R.R.,
1409
2005. Multiproxy evidence of Holocene climate variability from estuarine
1410
sediments, eastern North America. Paleoceanography 20, 4006e4027.
PR
OO
F
weakening ocean and atmosphere circulations. The fast decline
in annual insolation would also facilitate the amplitude of
centennial-scale climate variability during the mid-Holocene.
4. Our record from Lake Grinnell supports the notion that
a gradual change in climate forcing could trigger an shift in the
climate system. Also, changes in mean climate states are likely
associated with changes in climate variability, in terms of both
magnitude and frequency. This implies that a continued global
warming may induce changes in climate variability (including
climate extremes), especially if the climate system crosses
a threshold into a new mean state.
CT
ED
caused by the relatively larger dating errors at Lake Grinnell. Small
solar output variations would be amplified in the climate system
through several mechanisms, including the changes of ozone
production and stratospheric temperature structure caused by
changes in ultraviolet part of solar spectrum (Haigh, 1994; Shindell
et al., 1999), the changes in cosmic ray intensity on cloud formation
and precipitation (van Geel et al., 1999), the changes in global-scale
atmosphere circulation (Rind and Overpeck, 1994), and the changes
of sea ice and the North Atlantic deep water formation (Chapman
and Shackleton, 2000; Bond et al., 2001). Solar-forced centennialscale climate oscillations have also been widely hypothesized in
North America, for example, in the northern Great Plains (Yu and Ito,
1999), in the midcontinent (Tian et al., 2006), in the southwestern
United States (Asmerom et al., 2007), and in the sub-arctic Alaska
(Hu et al., 2003).
The solar forcing alone may not fully explain the mid-Holocene
shift in centennial-scale climate oscillations as the amplitude of
solar variations does not match the amplitude of the climate
variability (Fig. 9E). The observed amplitude shift of the centennialscale climate oscillations was possibly caused by seasonality
changes in insolation. For example, the contrast between summer
and winter insolation was larger during the early to mid-Holocene
than during the late Holocene (Fig. 9B). The strong seasonality
could cause intense internal oscillations in the oceanic and possibly
atmospheric circulations during the early to mid-Holocene, as
shown in the North Atlantic region (Schulz and Paul, 2002; Berner
et al., 2008) and in western Canada (Campbell et al., 1998). In
addition, the most rapid decrease in annual insolation during the
mid-Holocene might further facilitate climate oscillations (Fig. 9C).
Consequently, the early to mid-Holocene centennial-scale climate
variability was not only at higher frequency but also at higher
amplitude than the late Holocene, especially with the highest
amplitude during the mid-Holocene.
6. Conclusions and implication
CO
RR
E
1. The isotopic record at Lake Grinnell provides a detailed climate
history and landscape changes from the late-glacial to Holocene
for the northeastern United States. The negative excursion of
w1.5 & in d18O values at approximate 12.5e11.3 ka represents
a regional expression of the Younger Dryas cold event. The d18O
data show a high and stable mean value of about 7.4& during
the early to mid-Holocene from 11.3 to 5.8 ka and a negative
shift to a steady decreasing trend from 7.4& to 8.2& after
5.8 ka. The carbon-oxygen isotope covariance at 910-year periodicity likely reflects the controls of air temperature on both
d18O and d13C values in temperate lakes at millennial timescale.
2. The general Holocene climate trend recorded at Lake Grinnell
was dominantly controlled by the insolation. The remaining
Laurentide Ice Sheet would be responsible for the cool and dry
climate from 11.3 to 9 ka in eastern North America, as also
shown in other regional records. The 5.8 ka climate shift in the
long-term trend recorded at Lake Grinnell likely indicates a nonlinear response of atmosphere circulation to a gradual decline of
insolation during the mid-Holocene, as also documented in
other regional isotope records. Responding to this climate shift,
the hemlock declined in the northeastern United States, but
lagged w350e500 (143.5) years behind the climate shift.
3. The mid-Holocene shift in climate variability at centennialscale at 4.7 ka was detected not only in the periodicity shift
from 330 years to 500 years, but also in the decrease in the
amplitude. The increased periodicity at 4.7 ka was likely
controlled by the dominant frequency shift in variations of solar
output. The reduced amplitude at 4.7 ka was possibly caused by
the decrease in seasonal insolation contrast and accompanying
UN
1281
1282
1283
1284
1285
1286
1287
1288
1289
1290
1291
1292
1293
1294
1295
1296
1297
1298
1299
1300
1301
1302
1303
1304
1305
1306
1307
1308
1309
1310
1311
1312
1313
1314
1315
1316
1317
1318
1319
1320
1321
1322
1323
1324
1325
1326
1327
1328
1329
1330
1331
1332
1333
1334
1335
1336
1337
1338
1339
1340
1341
1342
1343
1344
1345
11
Please cite this article in press as: Zhao, C., et al., Holocene climate trend, variability, and shift documented by lacustrine stable isotope...,
Quaternary Science Reviews (2010), doi:10.1016/j.quascirev.2010.03.018
ARTICLE IN PRESS
OO
F
Lavoie, M., Richard, P.J.H., 2000. Postglacial water-level changes of a small lake in
southern Quebec, Canada. The Holocene 12, 621e634.
Lawrence, J.R., Gedzelman, S.D., White, J.W.C., Smiley, D., Lazov, P., 1982. Storm
trajectories in eastern US D/H isotopic composition of precipitation. Nature 296,
638e640.
Leng, M.J., Marshall, J.D., 2004. Palaeoclimate interpretation of stable isotope data
from lake sediment archives. Quaternary Science Reviews 23, 811e831.
Li, Y.X., Yu, Z.C., Kodama, K.P., 2007. Sensitive moisture response to Holocene
millennial-scale climate variations in the Mid-Atlantic region, USA. The Holocene 17, 3e8.
Maenza-Gmelch, T.E., 1997. Holocene vegetation, climate, and fire history of the
Hudson Highlands, southeastern New York, USA. The Holocene 7, 25e37.
McConnaughey, T.A., 1991. Calcification in Chara coralline: CO2 hydroxylation
generates protons for bicarbonate assimilation. Limnology and Oceanography
36, 619e628.
McConnaughey, T.A., Burdett, J., Whelan, J.F., Paull, C.K., 1997. Carbon isotopes in
biological carbonates: respiration and photosynthesis. Geochimica et Cosmochimica Acta 61, 611e622.
McConnaughey, T.A., Falk, R.H., 1991. Calcium-proton exchange during algal calcification. Biological Bulletin 180, 185e195.
McFadden, M.A., Mullins, H.T., Patterson, W.P., Anderson, W.T., 2004. Paleoproductivity of eastern Lake Ontario over the past 10,000 years. Limnology and
Oceanography 49, 1570e1581.
McKenzie, J.A., 1985. Carbon isotopes and productivity in the lacustrine and marine
environment. In: Stumm, W. (Ed.), Chemical Processes in Lakes. John Willey &
Sons, New York, pp. 99e118.
Mudelsee, M., 2000. Ramp function regression: a tool for quantifying climate
transitions. Computers & Geosciences 26, 293e307.
Mullins, H.T., 1998. Holocene lake level and climate change inferred from marl
stratigraphy of the Cayuga Lake Basin, New York. Journal of Sedimentary
Research 68, 569e578.
Mullins, H.T., Halfman, J.D., 2001. High-resolution seismic reflection evidence for
middle Holocene environmental change, Owasco Lake, New York. Quaternary
Research 55, 322e333.
Newby, P.E., Killoran, P., Waldorf, M.R., Shuman, B.N., Webb, R.S., Webb III, T., 2000.
14,000 years of sediment, vegetation, and water-level changes at the Makepeace
Ceder Swamp, southeastern Massachusetts. Quaternary Research 53, 352e368.
Nicholson, L.J., 1995. Ground-water flow in a buried valley and effects of quarry
dewatering on Howells Pond, Sussex County, New Jersey. Master Thesis, Drexel
University, USA.
Peixoto, J.P., Oort, A.H., 1983. The atmospheric branch of the hydrological cycle and
climate. In: Street-Perrott, A. (Ed.), Variations in the Global Water Budget.
Reidel, pp. 5e65.
Peteet, D.M., Daniels, R.A., Heusser, L.E., Vogel, J.S., Southon, J.R., Nelson, D.E., 1993.
Late-glacial pollen, macrofossils, and fish remains in northeastern U.S.A. e the
younger Dryas oscillation. Quaternary Science Reviews 12, 597e612.
Reimer, P.J., Baillie, M.G.L., Bard, E., Bayliss, A., Beck, J.W., Bertrand, C.J.H.,
Blackwell, P.G., Buck, C.E., Burr, G.S., Cutler, K.B., Damon, P.E., Edwards, R.L.,
Fairbanks, R.G., Friedrich, M., Guilderson, T.P., Hogg, A.G., Hughen, K.A.,
Kromer, B., McCormac, G., Manning, S., Ramsey, C.B., Reimer, R.W., Remmele, S.,
Southon, J.R., Stuiver, M., Talamo, S., Taylor, F.W., van der Plicht, J.,
Weyhenmeyer, C.E., 2004. INTCAL04 terrestrial radiocarbon age calibration,
0e26 cal kyr BP. Radiocarbon 46, 1029e1058.
Rind, D., Overpeck, J., 1994. Hypothesized causes of decade-to-century-scale climate
variability: climate model results. Quaternary Science Reviews 12, 357e374.
Sachs, J.P., 2007. Cooling of northwest Atlantic slope waters during the Holocene.
Geophysical Research Letters 34. doi:10.1029/2006GL028495.
Sandweiss, D.H., Maasch, K.A., Anderson, D.G., 1999. Transitions in the mid-Holocene. Science 283, 499e500.
Schulz, M., Mudelsee, M., 2002. REDFIT: estimating red-noise spectra directly from
unevenly spaced paleoclimatic time series. Computers & Geosciences 28,
421e426.
Schulz, M., Paul, A., 2002. Holocene climate variability on centennial-to-millennial
time scales: 1. Climate records from the North-Atlantic realm. In: Wever, G.E.A.
(Ed.), Climate Development and History of the North Atlantic Realm. Springer,
New York, pp. 41e54.
Shindell, D., Rind, D., Balachandran, N., Lean, J., Lonergan, P., 1999. Solar cycle
variability, ozone, and climate. Science 284, 305e308.
Shuman, B., Bravo, J., Kaye, J., Lynch, J.A., Newby, P., Webb III, T., 2001. Late
Quaternary water-level variations and vegetation history at Crooked Pond,
southeastern Massachusetts. Quaternary Research 56, 401e410.
Shuman, B., Donnelly, J.P., 2006. The influence of seasonal precipitation and
temperature regimes on lake levels in the northeastern United States during the
Holocene. Quaternary Research 65, 44e56.
Shuman, B., Newby, P., Huang, Y., Webb III, T., 2004. Evidence for the close climate
control of New England vegetation history. Ecology 85, 1297e1310.
Shuman, B., Webb III, T., Bartlein, P., Williams, J.W., 2002. The anatomy of a climatic
oscillation: vegetation change in eastern North America during the Younger
Dryas chronozone. Quaternary Science Reviews 21, 1777e1791.
Shuman, B.N., Newby, P., Donnelly, J.P., 2009. Abrupt climate change as an important
agent of ecological change in the Northeast U.S. throughout the past 15,000
years. Quaternary Science Reviews 28, 1693e1709.
Simpkins, W.W., 1995. Isotopic composition of precipitation in central Iowa. Journal
of Hydrology 172, 185e207.
CO
RR
E
CT
ED
Cwynar, L.C., Levesque, A.J., 1995. Chironomid evidence for late-Glacial climatic
reversals in Maine. Quaternary Research 43, 403e415.
Davis, M.B., 1969. Climatic changes in southern Connecticut recorded by pollen
deposition at Rogers Lake. Ecology 50, 409e422.
Davis, M.B., Spear, R.W., Shane, L.C.K., 1980. Holocene climate of New England.
Quaternary Science Reviews 14, 240e250.
Dean, W.E.J., 1974. Determination of carbonate and organic matter in calcareous
sediments and sedimentary rocks by loss-on-ignition: comparison with other
methods. Journal of Sedimentary Research 44, 242e248.
Denniston, R.F., Gonzalez, L.A., Asmerom, Y., Baker, R.G., Reagan, M.K., Bettis III, E.A.,
1999. Evidence for increased cool season moisture during the middle Holocene.
Geology 27, 815e818.
Dieffenbacher-Krall, A.C., Nurse, A.M., 2005. Late-glacial and Holocene record of
lake levels of Mathews Pond and Whitehead Lake, northern Marine, USA.
Journal of Paleolimnology 34, 283e310.
Drummond, C.N., Patterson, W.P., Walker, J.C.G., 1995. Climatic forcing of carbone
oxygen isotopic covariance in temperate-region marl lakes. Geology 23,
1031e1034.
Dwyer, T.R., Mullins, H.T., Good, S.C., 1996. Paleoclimatic implications of Holocene
lake-level fluctuations, Owasco Lake, New York. Geology 24, 519e522.
Edwards, T.W.D., Fritz, P., 1988. Stable isotope paleoclimate records from southern
Ontario: comparison of results from marl and wood. Canadian Journal of Earth
Sciences 25, 1379e1406.
Ellis, K.G., Mullins, H.T., Patterson, W.P., 2004. Deglacial to middle Holocene (16,600
to 6000 calendar years BP) change in the northeastern United States inferred
from multi-proxy stable isotope data, Seneca Lake, New York. Journal of Paleolimnology 31, 343e361.
Felzer, B., Oglesby, R.J., Webb III, T., Hyman, D.E., 1996. Sensitivity of a general
circulation model to changes in northeastern hemisphere ice sheets. Journal of
Geophysical Research 101, 19,077e19,092.
Foley, J.A., Kutzbach, J.E., Coe, M.T., Levis, S., 1994. Feedbacks between climate and
boreal forests during the Holocene epoch. Nature 371, 52e54.
Foster, D.R., Oswald, W.W., Faison, E.K., Doughty, E.D., Hansen, B.C.S., 2006. A
climatic driver for abrupt mid-Holocene vegetation dynamics and the hemlock
decline in New England. Ecology 87, 2959e2966.
Foster, D.R., Zebryk, T.M., 1993. Long-term vegetation dynamics and disturbance
history of a Tsuga-dominated forest in New England. Ecology 74, 982e988.
Fuller, J.L., 1998. Ecological impact of the mid-Holocene hemlock decline in
southern Ontario, Canada. Ecology 79, 2337e2351.
Ganopolski, A., Kubatzki, C., Claussen, M., Brovkin, V., Petoukhov, V., 1998. The
influence of vegetation-atmosphereeocean interaction on climate during the
mid-Holocene. Science 280, 1916e1919.
van Geel, B., Raspopov, O.M., Renssen, H., van der Plicht, J., Dergchev, V.A., Meijer, H.
A.J., 1999. The role of solar forcing upon climate change. Quaternary Science
Reviews 18, 331e338.
Haigh, J.D., 1994. The role of stratospheric ozone in modulating the solar radiative
forcing of climate. Nature 370, 544e546.
Hou, J., Huang, Y., Oswald, W.W., Foster, D.R., Shuman, B., 2007. Centennial-scale
compound-specific hydrogen isotope record of Pleistocene-Holocene climate
transition from southern New England. Geophysical Research Letters 34,
L19706. doi:10.1029/2007GL030303.
Hou, J., Huang, Y., Wang, Y., Shuman, B., Oswald, W.W., Faison, E., Foster, D.R., 2006.
Postglacial climate reconstruction based on compound-specific D/H ratios of
fatty acids from Blood Pond, New England. Geochemistry Geophysics
Geosystems 7, Q03008. doi:10.1029/2005GC001076.
Hu, F.S., Kaufman, D., Yoneji, S., Nelson, D., Shemesh, A., Huang, Y., Tian, J., Bond, G.,
Clegg, B., Brown, T., 2003. Cyclic variation and solar forcing of Holocene climate
in the Alaskan subarctic. Science 301, 1890e1893.
Huang, Y., Shuman, B., Wang, Y., Web III, T., 2002. Hydrogen isotope ratios of
palmitic acid in lacustrine sediments record late Quaternary climate variations.
Geology 30, 1103e1106.
Indermühle, A., Stocker, T.F., Joos, F., Fischer, H., Smith, H.J., Wahlen, M., Deck, B.,
Mastroianni, D., Tschumi, J., Blunier, T., Meyer, R., Stauffer, B., 1999. Holocene
carbon-cycle dynamics based on CO2 trapped in ice at Taylor Dome, Antarctica.
Nature 398, 121e126.
Ito, E., 2002. Application of stable isotope techniques to inorganic and biogenic
carbonates. In: Last, W.M., Smol, J.P. (Eds.), Developments in Paleoenvironmental Research Series. Tracking Environmental Change Using Lake Sediments, Physical and Geochemical Methods, vol. 2. Kluwer Academic Publishers,
The Netherlands, pp. 351e371.
Kaplan, M.R., Wolfe, A.P., 2006. Spatial and temporal variability of Holocene
temperature in the North Atlantic region. Quaternary Research 65, 223e231.
Kerwin, M.W., Overpeck, J.T., Webb, R.S., DeVernal, A., Rind, D.H., Healy, R.J., 1999.
The role of oceanic forcing in mid-Holocene Northern Hemispheric climate
change. Paleoceanography 14, 200e210.
Kim, S.T., O’Neil, J.R., 1997. Equilibrium and nonequilibrium oxygen isotope effects in
synthetic carbonates. Geochimica et Cosmochimica Acta 61, 3461e3475.
Kirby, M.E., Mullins, H.T., Patterson, W.P., Burnett, A.W., 2002. Late glacialeHolocene
atmospheric circulation and precipitation in the northeast United States
inferred from modern calibrated stable oxygen and carbon isotopes. Geological
Society of America Bulletin 114, 1326e1340.
Kutzbach, J.E., Gallimore, R., Harrison, S., Behling, P., Selin, R., Laarif, F., 1998. Climate
and biome simulations for the past 21,000 years. Quaternary Science Reviews
17, 473e506.
UN
1411
1412
1413
1414
1415
1416
1417
1418
1419
1420
1421
1422
1423
1424
1425
1426
1427
1428
1429
1430
1431
1432
1433
1434
1435
1436
1437
1438
1439
1440
1441
1442
1443
1444
1445
1446
1447
1448
1449
1450
1451
1452
1453
1454
1455
1456
1457
1458
1459
1460
1461
1462
1463
1464
1465
1466
1467
1468
1469
1470
1471
1472
1473
1474
1475
C. Zhao et al. / Quaternary Science Reviews xxx (2010) 1e13
PR
12
JQSR2725_proof ■ 20 April 2010 ■ 12/13
Please cite this article in press as: Zhao, C., et al., Holocene climate trend, variability, and shift documented by lacustrine stable isotope...,
Quaternary Science Reviews (2010), doi:10.1016/j.quascirev.2010.03.018
1476
1477
1478
1479
1480
1481
1482
1483
1484
1485
1486
1487
1488
1489
1490
1491
1492
1493
1494
1495
1496
1497
1498
1499
1500
1501
1502
1503
1504
1505
1506
1507
1508
1509
1510
1511
1512
1513
1514
1515
1516
1517
1518
1519
1520
1521
1522
1523
1524
1525
1526
1527
1528
1529
1530
1531
1532
1533
1534
1535
1536
1537
1538
1539
1540
ARTICLE IN PRESS
JQSR2725_proof ■ 20 April 2010 ■ 13/13
C. Zhao et al. / Quaternary Science Reviews xxx (2010) 1e13
Watts, W.A., 1979. Late Quaternary vegetation of central Appalachia and the New
Jersey coastal plain. Ecological Monographs 49, 427e469.
Webb III, T., Anderson, K.H., Webb, R.S., Bartlein, P.J., 1998. Late Quaternary climate
changes in eastern North America: a comparison of pollen-derived estimates
with climate model results. Quaternary Science Reviews 17, 587e606.
Witte, R.W., 2001. Late Wisconsian end moraines in northwestern New Jersey:
observations on their distributions, morphology, and composition. In: 66th
Field Conference of Pennsylvania Geologists, pp. 81e98.
Witte, R.W., Monteverde, D.H., 2006. Quaternary Geology and Geologic Material
Resources of the Newton East Quadrangle, Sussex County, New Jersey. New
Jersey Geological Survey, Bureau of Geology and Topography, Trenton, New
Jersey.
Yu, Z.C., 2007. Rapid response of forested vegetation to multiple climatic oscillations
during the last deglaciation in the northeastern United States. Quaternary
Research 67, 297e303.
Yu, Z.C., Ito, E., 1999. Possible solar forcing of century-scale drought frequency in the
northern Great Plains. Geology 27, 263e266.
Yu, Z.C., McAndrews, J.J., Eicher, U., 1997. Middle Holocene dry climate caused by
change in atmospheric circulation pattern: evidence from lake levels and stable
isotopes. Geology 25, 251e254.
Zhao, Y., Yu, Z.C., Zhao, C., Hemlock (Tsuga Canadensis) declines at 9800 and 5300
cal yr BP caused by Holocene climatic shifts in northeastern North America. The
Holocene, in press.
CO
RR
E
CT
ED
PR
OO
F
Smith, M.A., Hollander, D.J., 1999. Historical linkage between atmospheric circulation patterns and the oxygen isotopic record of sedimentary carbonates from
Lake Mendota, Wisconsin, USA. Geology 27, 589e592.
Springer, G.S., Rowe, H.D., Hardt, B., Edwards, R.L., Cheng, H., 2008. Solar forcing of
Holocene droughts in a stalagmite record from West Virginia in east-central North
America. Geophysical Research Letters 35, L17703. doi:10.1029/2008GL034971.
Steig, E., 1999. Mid-Holocene climate change. Science 286, 1485e1487.
Stuiver, M., 1968. Oxygen 18 content of atmospheric precipitation during the last
11,000 years in the Great Lakes region. Science 162, 994e997.
Stuiver, M., 1970. Oxygen and carbon isotope ratios of fresh-water carbonates as
climatic indicators. Journal of Geophysical Research 75, 5247e5257.
Suter, S.M., 1985. Late-glacial and Holocene vegetation history in southeastern
Massachusetts: a 14,000 year pollen record. Current Research in the Pleistocene
2, 87e88.
Tian, J., Nelson, D.M., Hu, F.S., 2006. Possible linkages of late-Holocene drought in
the North American midcontinent to Pacific Decadal oscillation and solar
activity. Geophysical Research Letters 33, L23702. doi:10.1029/2006GL028169.
Valgenti, F., Jr., Salvato, J.C., Johnson, H.S., Campbell, C.A., Long, W.R.M., Hand, J.W.,
Jeffers JR, H.W., MacDonald, R.D., String, A.W., Hart, D.H., Cubbage, C.M., 1950.
New Jersey fisheries survey, report number one, lakes and ponds.
Viau, A.E., Gajewski, K., Sawada, M.C., Fines, P., 2006. Millennial-scale temperature
variations in North America during the Holocene. Journal of Geophysical
Research 111. doi:10.1029/2005JD006031.
UN
1541
1542
1543
1544
1545
1546
1547
1548
1549
1550
1551
1552
1553
1554
1555
1556
1557
1558
1559
13
Please cite this article in press as: Zhao, C., et al., Holocene climate trend, variability, and shift documented by lacustrine stable isotope...,
Quaternary Science Reviews (2010), doi:10.1016/j.quascirev.2010.03.018
1560
1561
1562
1563
1564
1565
1566
1567
1568
1569
1570
1571
1572
1573
1574
1575
1576
1577
1578