Electronic Supplementarymaterial for the Paleohydrology of Sluice Pond, NE Massachusetts

Electronic Supplementarymaterial for the Paleohydrology of Sluice Pond, NE Massachusetts

Electronic SupplementaryMaterial for “The paleohydrology of Sluice Pond, NE Massachusetts, and its regional significance”

J. Bradford Hubeny, Francine M.G. McCarthy, Jonathan Lewis, Matea Drljepan, Cameron Morissette, John W. King, Mark Cantwell, Nicole M. Hudson, Mary Lynne Crispo

Age Constraints

The age model for core SP09KC2 is reported in Fig. 3. Here we offer supplemental information on the chronology, including 14C calibration and correlation to other dated records (210Pb and paleomagnetics). A detailed list of chronologic control points used in the age model is presented in Table ESM1.

Five AMS 14C dates from core SP09KC2 and oneAMS 14C date from core SP07PC4 were measured on terrestrial macrofossils (Table ESM2). Each sample wastreated with an acid-base-acid (ABA) procedure (Olsson 1986). Analyses were conducted at the NSF-Arizona AMS Laboratory. Calibrations were calculated using the CLAM 2.2 program (Blaauw 2010) and Intcal13 calibration data set (Reimer et al. 2013).

Additional chronologic control of core SP09KC2 was provided by paleomagneticdirectional and intensity correlations to well-dated regional records. Paleomagnetic data were measured at 1-cm resolution on u-channels sub-sampled from the basin core, using a 2-G® cryogenic magnetometer. Progressive alternating field (AF) demagnetization was performed; a stable signal was noted after the 150 Oe step, and these data were used for correlation to regional dated paleomagnetic records. Correlations were based on similarities of major features in each record, and all interpretations were consistent within the framework of the radiocarbon constraints. Paleomagnetic secular variation records typically correspond with centennial-scale variability(King and Peck 2001), and therefore the error on correlations was assigned a value of ± 100 years. Paleomagnetic inclination was correlated to the dated NE inclination record of King and Peck (2001) (Fig.ESM1). Relative paleomagnetic intensity was calculated as Natural Remanent Magnetization normalized to Anhysteretic Remanent Magnetization (NRM/ARM). The record was correlated to the dated record of St-Onge et al. (2003) (Fig.ESM2).

The recent chronology was constrained with a 210Pb and 137Cs CIC age model (Crescenzi et al. 2010). This chronology was quantified on a short surface core, SP10KC1. Core SP09KC2 was correlated to SP10KC1, and subsequently to its age model, by matching the volume magnetic susceptibility records (Fig.ESM3), confirmed with visual lithologic descriptions. Long core SP09KC2 has a compressed upper section and is missing <10 cm of sediment from the top as a result of the coring process, which used heavy weight to achieve deep penetration. The presence of three distinct peaks in the record enabled stratigraphic correlation between the cores.

The final age model, reported in Fig. 3, was calculated using 18 age constraints (Table ESM1) as a fourth-order polynomial regression (Blaauw 2010). One thousand iterations were calculated, and the calculated goodness-of-fit was 34.21.

Cyclic organic matter deposition

The high-resolution (1-cm) record of organic matter from LOI analysis (Fig. 5) enabled us to examine cyclic components of the record through spectral analysis of the time series. Spectral analysis was calculated on the LOI-OC timeseries in kSpectra 3.4,using the multi-taper method with 3 tapers. Because this technique requires evenly spaced time series, data were resampled in AnalySeries 2.0 (Paillard et al. 1996)at a temporal resolution of 24.8 yr (average sampling interval) prior to analysis. The 99%, 95%, and 90% confidence intervals were calculated for the spectral analysis, assuming a “red noise” first-order autoregressive (AR(1)) background (Mann and Lees 1996). Results are reported in Fig.ESM4, and peaks extending above the AR(1) background curves represent statistically significant frequencies within the record.

The general trend of OC in Sluice Pond mirrors that of regional temperature variability in the Holocene (Huang et al. 2002) (Fig. 5). Higher frequency components of the OC record, as determined through spectral analysis, assist interpretation of climatic controls for Sluice Pond (Fig. ESM4). A range of significant multidecadal and centennial periodicities were observed, demonstrating the dynamic nature of the Sluice Pond record. Similar periodic components have been found in numerous other Holocene climate proxy records, including reconstructions of the circumpolar vortex(Kirby et al. 2002), solar variability(Solanski et al. 2004) and ocean-atmospheric teleconnection patterns such as the Atlantic Multidecadal Oscillation(Delworth and Mann 2000; Hubeny et al. 2006; Knight et al. 2005). Although this study was unable to unravel the ultimate controlling factorson Sluice Pond cyclic variability, the significant cycles observed here corroborate and expand upon other research that focuses on such variability.

Table ESM1Details and list of 18 age constraints from core SP09KC2 used for the generation of CLAM 2.2 (Blaauw 2010)age model reported in Fig. 3

Core Depth (cm) / Date Type / Calibratedage
IntCal13 / Error / 14C 2probability / Reference Information
0 / 210Pb Correlation / -35 / 5 / Fig. ESM3
12 / 210Pb Correlation / 80 / 5 / Fig. ESM3
86 / AMS 14C / 1839 / 12 / 7.7 / Table ESM2
1924 / 64 / 87.1
101 / Inclination Correlation / 2350 / 100 / Fig. ESM1
102 / Intensity Correlation / 2150 / 100 / Fig. ESM3
130 / Intensity Correlation / 3300 / 100 / Fig. ESM3
166 / AMS 14C / 3619 / 1 / 0.2 / Table ESM2
3859 / 227 / 94.8
166 / Intensity Correlation / 4100 / 100 / Fig. ESM3
181 / Intensity Correlation / 4600 / 100 / Fig. ESM3
224 / Intensity Correlation / 5850 / 100 / Fig. ESM3
285 / Intensity Correlation / 7050 / 100 / Fig. ESM3
286 / AMS 14C / 7364 / 41 / 26.5 / Table ESM2
7456 / 50 / 68.4
290 / Inclination Correlation / 6830 / 100 / Fig. ESM1
305 / Inclination Correlation / 7450 / 100 / Fig. ESM1
339 / Inclination Correlation / 8620 / 100 / Fig. ESM1
366 / AMS 14C / 9146 / 140 / 95 / Table ESM2
452 / Inclination Correlation / 11240 / 100 / Fig. ESM1
467 / AMS 14C / 11409 / 210 / 94.9 / Table ESM2
11685 / 1 / 0.1

Table ESM2Sluice Pond AMS radiocarbon data. All dates from terrestrial macrofossils and analyses conducted at the NSF – Arizona AMS Laboratory. Calibrations were calculated using the CLAM 2.2 program (Blaauw 2010) and Intcal13 calibration data set (Reimer et al. 2013)

Sample Name / Composite Depth (cm) / Material / Lab ID / 13C (‰) / 14C age
(yr BP) / 2 calibrated age ranges
(cal BP) / Probability
SP09KC2S2, 10cm / 86 / Leaf / AA87772 / -28.3 / 1957 ± 36 / 1827 – 1851
1860 - 1988 / 7.7
87.1
SP09KC2S3 0cm / 166 / Leaf / AA99261 / -23.3 / 3551 ± 89 / 3618 – 3620
3632 – 4086 / 0.2
94.8
SP09KC2S5, 10cm / 286 / Leaf / AA87773 / -29.3 / 6515 ± 44 / 7323 – 7404
7406 – 7505 / 26.5
68.4
SP09KC2S6 0cm / 366 / Leaf / AA99262 / -29.1 / 8172 ± 58 / 9006 – 9286 / 95.0
SP09KC2CC / 467 / Seed casings / AA87774 / -27.1 / 9905 ± 73 / 11199 – 11619
11684 – 11685 / 94.9
0.1
SP07PC4Sec2, 88.5cm / 189 / Wood / AA92031 / -27.7 / 6982 ± 99 / 7625 – 7637
7653 - 7979 / 0.8
94.1
/ Fig. ESM1Paleomagnetic inclination from sediment core SP09KC2 and its correlation to the dated NE inclination record of King and Peck (2001). Correlations are consistent within the framework of all other chronologic constraints. Specific dates used in the age model are reported in Table ESM1
/ Fig. ESM2 Relative paleomagnetic intensity (NRM/ARM) from sediment core SP09KC2 and its correlation to the dated record of St-Onge et al. (2003) Correlations are consistent within the framework of all other chronologic constraints. Specific dates used in the age model are reported in Table ESM1
/ Fig. ESM3 Correlation of volume magnetic susceptibility of cores SP09KC2 and SP10KC1. SP10KC1 was dated with a 210Pb and 137Cs CIC age model (Crescenzi et al. 2010)and affords chronologic constraint to the correlated portion of SP09KC2. Specific dates used in the age model are reported in Table ESM1
Fig. ESM4 Spectral analysis of the LOI organic carbon record. Gray lines illustrate 99%, 95%, and 90% confidence intervals calculated based on AR(1) red noise assumption(Mann and Lees 1996). All significant periodicities are labeled in bold

References used in the Electronic Supplementary Material

Blaauw M (2010) Methods and code for 'classical' age-modelling of radiocarbon sequences. Quat Geochronol 5: 512-518

Crescenzi E, Hubeny JB, Moran SB, Kelly RP (2010) Historic trace metal accumulation in Sluice Pond, MA: Linkages between industry and the environment. GSA Abstr Prog 42: 183

Delworth TL, Mann ME (2000) Observed and simulated multidecadal variability in the Northern Hemisphere. Clim Dyn 16: 661-676

Huang Y, Shuman B, Wang Y, Webb T, III (2002) Hyrdrogen isotope ratios of palmitic acid in lacustrine sediments record late Quaternary climate variations. Geology 30: 1103-1106

Hubeny JB, King JW, Santos A (2006) Subdecadal to multidecadal cycles of Late Holocene North Atlantic climate variability preserved by estuarine fossil pigments. Geology 34: 569-572

King JW, Peck J (2001) Use of paleomagnetism in studies of lake sediments. In: Last WM, Smol JP (eds) Tracking environmental change using lake sediment Volume 1: Basin analysis, coring, and chronological techniques. Kluwer Academic Publishers, Dordrecht, The Netherlands, pp 371-389

Kirby ME, Mullins HT, Patterson WP, Burnett AW (2002) Late glacial-Holocene atmospheric circulation and precipitation in the northeast United States inferred from modern calibrated stable oxygen and carbon isotopes. Geol Soc Am Bull 114: 1326-1340

Knight JR, Allan RJ, Folland CK, Vellinga M, Mann ME (2005) A signature of persistent natural thermohaline circulation cycles in observed climate. Geophys Res Lett 32, L20708:

Mann ME, Lees JM (1996) Robust estimation of background noise and signal detection in climatic time series. Clim Change 33: 409-445

Olsson I (1986) Radiometric Methods. In: Berglund B (ed) Handbook of Holocene palaeoecology and palaeohydrology. John Wiley & Sons, Chichester, pp 273-312

Paillard D, Labeyrie L, Yiou P (1996) Macintosh program performs time-series analysis. Eos 77: 379

Reimer PJ, Bard E, Bayliss A, Beck JW, Blackwell PG, Bronk Ramsey C, Buck CE, Edwards RL, Friedrich M, Grootes PM, Guilderson TP, Haflidason H, Hajdas I, Hatté C, Heaton TJ, Hoffmann DL, Hogg AG, Hughen KA, Kaiser KF, Kromer B, Manning SW, Niu M, Reimer RW, Richards DA, Scott EM, Southon JR, Turney CSM, van der Plicht J (2013) IntCal13 and Marine13 radiocarbon age calibration curves, 0-50,000 years cal BP. Radiocarbon 55: 1869-1887

Solanski SK, Usoskin IG, Kromer B, Schüssler M, Beer J (2004) Unusual activity of the Sun during recent decades compared to the previous 11,000 years. Nature 431: 1084-1087

St-Onge G, Stoner JS, Hillaire-Marcel C (2003) Holocene paleomagnetic records from the St. Lawrence Estuary, eastern Canada: Centennial- to millennial-scale geomagnetic modulation of cosmogenic isotopes. Earth Planet Sci Lett 209: 113-130