Collection, analysis, and age-dating of sediment cores from Herring River wetlands and other nearby wetlands in Wellfleet, Massachusetts, 2015–17

A piston coring system was used to collect sediment cores from coastal wetlands in 2015, 2016, and 2017. The PVC core liner (diameter 11 cm, 1-2 m in length) was fitted with a gasketed piston that was placed on the sediment surface. The clear, sharpened core liner was pushed down into the marsh subsurface, while the piston was maintained at the marsh surface via tension on the piston. We visually observed the sediment surface to ensure that the soil column did not compact during collection. Once the core reached the desired depth, the core liner and piston were removed from the marsh with a pulley system. The process date reflects the last occurrence of the work.

At the end of each core collection date, the sediment cores were immediately returned to the USGS Woods Hole Coastal and Marine Science Center. They were placed in a refrigerator for 1 to 3 days, then split vertically, sectioned at 1 or 2 cm intervals, frozen, and then freeze dried for 7 days, until sediment weights did not change further. Dry bulk density was determined as the dry weight of a known volume of sample. Approximately 5 g of dried sediment sample was blended and homogenized prior to sealing in a jar for a minimum of three weeks and then placed on a planar-type gamma counter for 24 to 48 hours to measure 7Be, 137Cs, 210Pb, and 226Ra at 477, 662, 46.5 and 352 kiloelectronvolts (KeV) energies respectively (Canberra Inc., USA). Detector efficiency was determined from EPA standard pitchblende ore in the same geometry as the samples. Excess 210PB was calculated as the decay-corrected difference between total 210Pb and supported 210Pb (considered to be equal to 226Ra). Activities of 7Be, 137Cs, and excess 210Pb were decay corrected to time of collection, using their respective half-lives. Suppression of low energy peaks by self-absorption was corrected for according to Cutshall and others, 1983. Gamma spectroscopy detection limits were determined in APTEC software for each sample; refer to the attribute accuracy section of this metadata for further details. Values reported are above this limit, while values below are reported as 0. Core sections not analyzed are reported as blank cells. Sediment ages and accretion rates were calculated with the Constant Rate of Supply (CRS) excess 210Pb age model, a variant on the advection-decay equation (Appleby and Oldfield, 1978; Goldberg, 1963). This model assumes that 210Pb supply to the sediment surface is constant through time, but allows for changing sedimentation rates, in addition to decay, to control the down-core activity of excess 210Pb. The common form of the CRS model, as derived by Appleby and Oldfield (1978), solves for age based on the distribution of exces 210Pb in the sediment record. Prior to application of the age model, excess 210Pb profiles were evaluated to ensure they were sufficiently resolved to apply the CRS model without bias towards ages that are too old or accretion rates that are too low at depth (Binford, 1990). All gamma analyses were ongoing from 2015 and completed in 2018.
Appleby, P.G., and Oldfield, F., 1978, The calculation of lead-210 dates assuming a constant rate of supply of unsupported 210Pb to the sediment: Catena, v. 5, issue 1, p. 1–8, https://doi.org/10.1016/S0341-8162(78)80002-2.
Binford, M.W., 1990, Calculation and uncertainty analysis of 210 Pb dates for PIRLA project lake sediment cores: Journal of Paleolimnology, v. 3, issue 3, p. 253-267, https://doi.org/10.1007/BF00219461.
Cutshall, N.H., Larsen, I.L., and Olsen, C.R., 1983, Direct analysis of 210 Pb in sediment samples—Self-absorption corrections: Nuclear Instruments and Methods in Physics Research, v. 206, issues 1–2, p. 309–312, https://doi.org/10.1016/0167-5087(83)91273-5.
Goldberg, E.D., 1963, Geochronology with 210 Pb, in Miller, J.A., convener, Radioactive dating: International Atomic Energy Agency Symposium on Radioactive Dating, Athens, Greece, November 19-23, 1962, [Proceedings], p. 121-131.

A subsample of freeze dried sediment (0.5 g) was ball-milled to a fine powder and packaged for analysis at the USGS Woods Hole Coastal and Marine Science Center. Samples were analyzed for Total Carbon and Nitrogen at the University of California-Davis (UC Davis) Stable Isotope Facility (SIF). Weight percent carbon and nitrogen and the isotopic signature of carbon and nitrogen in organic matter was analyzed with an Elementar Vario EL Cube or Micro Cube elemental analyzer interfaced to a PDZ Europa 20-20 isotope ratio mass spectrometer. During analysis, samples are interspersed with several replicates of at least four different laboratory reference material; these reference materials have a long-term standard deviation of 0.2 per mil for 13C and 0.3 per mil for 15N. Analyses occurred concurrently with analyses in the previous process step; after radionuclide analysis was completed on a core, we would begin the total carbon and nitrogen analysis (as other cores were still in the process of being analyzed for radionuclides).

Land surface elevations were measured by Trimble Real-Time Kinematic (RTK) GPS across the collection area after core collection. Surveyors returned to the core collection site on the dates listed below, navigated to the core collection locations, and surveyed several points. HROC and HRHT: 01/11/2017; MYN, MYUR, and HRCH: 12/11/2017; HRP and HRHK: 07/11/2018. Average vertical accuracy of this method is +/- 5 centimeters. No formal positional accuracy tests were conducted.

Raw data were entered into an Excel spreadsheet where all calculations were completed. The Excel spreadsheet was exported as a CSV file from Excel for Mac version 15.33. The CSV file was processed in MATLAB to round calculated values to appropriate place values and again exported as a comma separated comma separated text file (*.csv). Note that all calculations were performed prior to rounding and truncating values reported in this data release.