Annual resolution sediment layers, known as varves, can provide continuous and high-resolution chronologies of sedimentary sequences. In addition, varve counting is not burdened with the high laboratory costs of geochronological analyses. Despite a more than 100-year history of use, many existing varve counting techniques are time consuming and difficult to reproduce. We present countMYvarves, a varve counting toolbox which uses sliding-window autocorrelation to count the number of repeated patterns in core scans or outcrop photos. The toolbox is used to build an annually-resolved record of sedimentation rates, which are depth-integrated to provide ages. We validate the model with repeated manual counts of a high sedimentation rate lake with biogenic varves (Herd Lake, USA) and a low sedimentation rate glacial lake (Lago Argentino, Argentina). In both cases, countMYvarves is consistent with manual counts and provides additional sedimentation rate data. The toolbox performs multiple simultaneous varve counts, enabling uncertainty to be quantified and propagated into the resulting age-depth model. The toolbox also includes modules to automatically exclude non-varved portions of sediment and interpolate over missing or disrupted sediment. CountMYvarves is open source, runs through a graphical user interface, and is available online for download for use on Windows, macOS or Linux at https://doi.org/10.5281/zenodo.4031811.

Deposits of at least three glaciations are present in New Jersey and the New York City area. The oldest deposits are magnetically reversed. Pollen and stratigraphic relations suggest that they are from the earliest Laurentide advance at ~2.4 Ma. Deposits of a second advance are overlain by peat dated to 41 ka and so are pre-Marine Isotope Stage (pre-MIS) 2. Their relation to marine deposits indicates that they predate MIS 5 but postdate MIS 11 and may postdate MIS 7 or 9, suggesting an MIS 6 age. The most recent deposits are of MIS 2 (last glacial maximum [LGM]) age. Radiocarbon dates and varve counts tied to glacial-lake events indicate that LGM ice arrived at its terminus at 25 ka, stood at the terminus until ~24 ka, retreated at a rate of 80 m/yr until 23.5 ka, and then retreated at a rate of 12 m/yr to 18 ka. At 18 ka the retreat record connects to the base of the North American Varve Chronology at Newburgh, New York. The 25–24 ka age for the LGM is slightly younger than, but within the uncertainty of, cosmogenic ages; it is significantly older than the oldest dated macrofossils in postglacial deposits in the region.

Age-depth models provide essential temporal frameworks in paleoenvironmental science. We use a sample of 80 recently-published age-depth models to comment on current practices in building and reporting radiocarbon-based age-depth models. We address options for model building, sampling strategies, dating densities, and best practices for reporting age-depth models and associated data. Our review reveals incomplete reporting of 14C ages, model-building methods, age-depth models and associated meta-data in many recent studies. All information needed to evaluate, reproduce and update an age-depth model should accompany every published model. We also present a case study of building age-depth models for a lake sediment core that has both 14C ages and an independent varve chronology. The case study illustrates that choosing the ‘best model’ is not a simple task, and that model accuracy is ultimately controlled by differences between 14C ages and true age that likely occur in many late Quaternary records.

Central Asia is highly vulnerable to large earthquakes, yet existing records of past seismic activity in this area are still insufficient to reliably assess regional earthquake hazard on longer timescales. Within this study, the sediments of Sary Chelek, a mountain lake in the western Kyrgyz Tian Shan, were investigated to explore its potential as a natural paleoseismic archive. The lacustrine deposits are characterized by a succession of annually laminated (varved) sediments overlying event deposits that consist of large-scale turbidites and distorted lake sediments, similar to earthquake-related deposits described from other lake sediment records. Microscopic sediment analysis furthermore revealed distorted varves in the laminated sequence that closely resemble earthquake-related soft-sediment deformation structures. Varve counting and radiometric dating determine the formation of the distorted varves and the emplacement of the large-scale event deposits to the early 1990s and mid-1940s, respectively. This is in good temporal agreement with the occurrence of two large earthquakes that struck western Kyrgyzstan in AD 1992 and AD 1946. These results and particularly the precise age control of the Sary Chelek sediment record highlight its potential for establishing a long and precisely dated record of regional earthquake activity.

Formation of annually laminated sediments in Summit Lake, White Pass, British Columbia is controlled by runoff generated by snowpack and glacier melt and major rainfall events. The 700-yr varve record is divided into two subannual series (early and late) based on sedimentological criteria and sedimentary structures within each varve. A comparison of recent subannual laminae with nearby meteorological records supports the interpretation they are formed by river discharge events generated by major snow and glacier melt events and large late-summer rainfall events. A significant correlation exists between the late subannual thickness series and the size of the largest rainfall events in late summer. The long record indicates there was an abrupt increase in the thickness and frequency of major rainfall-induced sedimentary events at the end of the seventeenth century. In addition, the frequency of laminae generated by early runoff events also increased. However, early subannual varve thickness component remains statistically the same as the thickness prior to the end of the seventeenth century. This suggests the change in varve thickness at this time is due to increases in major late-summer rainfall frequency rather than increased sediment availability caused by regional Little Ice Age glacier advances.

We examined the hydroclimatic signal in a record of annual lamina (varve) thickness from High Arctic Lake A, Ellesmere Island (83°00.00′N, 75°30.00′W). In this unglacierized catchment, nival melt is the dominant source for meltwater and transport of sediment to the lake, and autumn snowfall is highly influential on varve thickness through the amount of snow available for melt in the following year. For the period during which climatic data are available, varve thickness in Lake A was significantly correlated (r = 0.50, p < 0.01) with the cumulative snowfall from August to October (ASO) during the previous year and to a lesser extent, ASO mean daily temperature (r = 0.39, p < 0.01) at Alert, Nunavut (175 km east). The varve thickness record, interpreted as a proxy record of ASO snowfall and by extension, ASO temperature, indicated above-mean conditions during five periods of the past millennium, including most of the 20th century. These results corresponded well to other available high-resolution proxy climate records from the region, with some discrepancies prior to AD 1500 and during the period AD 1700–1900.

We present a varve thickness chronology from glacier-dammed Iceberg Lake in the southern Alaska icefields. Radiogenic evidence confirms that laminations are annual and record continuous sediment deposition from A.D. 442 to A.D. 1998. Varve thickness is positively correlated with Northern Hemisphere temperature trends, and more strongly with a local, ∼600 yr long tree ring width chronology. Varve thickness increases in warm summers because of higher melt, runoff, and sediment transport (as expected), but also because shrinkage of the glacier dam allows shoreline regression that concentrates sediment in the smaller lake. Varve thickness provides a sensitive record of relative changes in warm season temperatures. Relative to the entire record, temperatures implied by this chronology were lowest around A.D. 600, warm between A.D. 1000 and A.D. 1300, cooler between A.D. 1500 and A.D. 1850, and have increased dramatically since then. Combined with stratigraphic evidence that contemporary jökulhlaups (which began in 1999) are unprecedented since at least A.D. 442, this record suggests that 20th century warming is more intense, and accompanied by more extensive glacier retreat, than the Medieval Warm Period or any other time in the last 1500 yr.

An integrated seismo- and chronostratigraphic investigation at Lac Dasserat, northwestern Quebec, identified 74 separate failures within eight event horizons. Horizons E and B, and H and G have strong or moderately-strong multi-landslide signatures, respectively, composed of 11-23 failures, while horizons F, D, C, and A have minor landslide signatures consisting of a single or pair of deposit(s). Cores collected at six sites recovered glacial Lake Ojibway varve deposits that are interbedded with the event horizons. The correlation of the varves to the regional Timiskaming varve series allowed varve ages or ranges of varve ages to be determined for the event horizons. Horizons H, G, E, and B are interpreted to be evidence of paleoearthquakes with differing levels of interpretative confidence, based on the relative strength of the multi-landslide signatures, the correlation to other disturbed deposits of similar age in the region, and the lack or possibility of alternative aseismic mechanisms. The four interpreted paleoearthquakes occurred between 9770 ± 200 and 8470 ± 200 cal yr BP, when glacial Lake Ojibway was impounded behind the Laurentide Ice Sheet during deglaciation. They probably represent an elevated period of seismicity at deglaciation that was driven by crustal unloading.

Areal mapping of the middle Wisconsin varved clay site along Sixmile Creek near Ithaca, New York, has revealed a much more widespread and varied array of sediments than previously thought. Lacustrine clays, some varved, are interbedded with sands and gravels interpreted as sub-aqueous fan deposits, and both are overlain by a deformation till. Nine radiocarbon dates indicate a 34–37 14C ka BP age for the lacustrine sediment, with the deformation till less than a few thousand years younger. Beneath this sequence is a deposit dated at ± 42 14C ka BP. Both strata represent a tundra climate with a mean July temperature of about 10°C. The Sixmile Creek deformation till must correlate with the 35 14C ka BP till along the Genesee River, 125 km to the NW, and defines a Cherrytree stade glacial advance into the Appalachian Plateau, much further south than what has generally been accepted. Such an advance would require drainage from a proglacial lake in the western Ontario basin to flow westward instead of northeastward. The Sixmile strata suggest a colder than accepted middle Wisconsin stage. Recent data indicate that this stage is one of progressive cooling, with large climatic fluctuations.

The varve record from High Arctic, proglacial Bear Lake reveals a regionally coherent hydroclimatic signal as well as complexities due to changing hydroclimatic and limnologic conditions. Varve formation is strongly dependent on underflows that exhibit variability in strength during the past 750 yr. Periods with reduced underflow sedimentation and accumulation rates fail to produce varves in the distal part of the lake. Isolated coarse silt and sand grains occur in 80% of the varves and are interpreted to be niveo-aeolian in origin. Coarse (>500 μm) sand grains deposited on the lake ice by strong winter winds are notably less common since A.D. 1850, likely due to reduced storminess. Regression of the varve thickness record with meteorological records indicates high correlations with autumn (September and October) temperatures and total monthly snowfall. These correlations are best at times when underflow activity is sufficiently strong to produce varves throughout the lake. The close association with warmer temperatures and snow-bearing synoptic systems moving north in Baffin Bay suggests that the primary climate signal in the varves is varying autumn snow pack that controls nival discharge in the following year. The similarity between the other records of melt season temperature and sea-ice cover and the Bear Lake record suggests that summer and autumn conditions were generally similar across the Baffin Bay region through much of the last millennium.

A 2200-yr long, high-resolution (∼5 yr) record of drought variability in northwest Montana is inferred from diatoms and δ18O values of bio-induced carbonate preserved in a varved lacustrine core from Foy Lake. A previously developed model of the diatom response to lake-level fluctuations is used to constrain estimates of paleolake levels derived from the diatom data. High-frequency (decadal) fluctuations in the de-trended δ18O record mirror variations in wet/dry cycles inferred from Banff tree-rings, demonstrating the sensitivity of the oxygen-isotope values to changes in regional moisture balance. Low frequency (multi-centennial) isotopic changes may be associated with shifts in the seasonal distribution of precipitation. From 200 B.C. to A.D. 800, both diatom and isotope records indicate that climate was dry and lake level low, with poor diatom preservation and high organic carbon: nitrogen ratios. Subsequently, lake level rose slightly, although the climate was drier and more stable than modern conditions. At A.D. 1200, lake level increased to approximately 6 m below present elevation, after which the lake fluctuated between this elevation and full stage, with particularly cool and/or wetter conditions after 1700. The hydrologic balance of the lake shifted abruptly at 1894 because of the establishment of a lumber mill at the lake's outlet. Spectral analysis of the δ18O data indicates that severe droughts occurred with multi-decadal (50 to 70 yr) frequency.

Sedimentological analyses of 289 years (AD 1718–2006) of varved sediment from Shadow Bay, southwest Alaska, were used to investigate hydroclimate variability during and prior to the instrumental period. Varve thicknesses relate most strongly to total annual discharge (r2 = 0.75, n = 43, p < 0.0001). Maximum annual grain size depends most strongly on maximum spring daily discharge (r2 = 0.63, n = 43, p < 0.0001) and maximum annual daily discharge (r2 = 0.61, n = 43, p < 0.0001), while varve thickness is poorly correlated with maximum annual grain size (r2 = 0.004, n = 287, p = 0.33). Relations between varve thickness and annual climate variables (temperature, precipitation, North Pacific (NP) and Pacific Decadal Oscillation (PDO) indices) are insignificant. On multidecadal timescales, however, regime shifts in varve thickness and total annual discharge coincide with shifts in NP and PDO indices. Periods with increased varve thickness and total annual discharge were associated with warm PDO phases and a strengthened Aleutian Low. The varve-inferred record of PDO suggests that any periodicity in the PDO varied over time, and that the early 19th century marked a transition to a more frequent or detectable shifts.

During its maximum extension, the Saalian ice cap reached into the central Netherlands, where glacier tongues excavated over 100 m deep basins in the unconsolidated Middle and Early Pleistocene sediments. The basins are filled by relatively thick successions of Late Saalian, Eemian and Weichselian sediments. The fill of the Amsterdam glacial basin is among the best known and studied in the Netherlands. The Late Saalian sediments consist mainly of warves and ill-bedded clays and silts with, along its southern margin, influxes of sands from the surrounding ice-pushed ridges. During deposition of these sediments, the Amsterdam basin formed part of a large lake extending into the present North Sea. Draining of this lake at the end of the Late Saalian left small, shallow pools at the site of the glacial basins.

Late Saalian and Eemian sediments are probably separated by a short break, although sedimentation may have been continuous in the deepest part of the basin. The Eemian deposits consist in main lines of a thin, diatom-rich sapropel at the base, overlain by an up to 30 m thick clay-rich sequence covered by a wedge of sand that measures more than 20 m in the northern part of the basin and that peters out southwards. As appears from the fauna, most of the clays were deposited in a lagoonal setting shielded behind a threshold and/or barrier. The rate of sediment supply was low so that lagoonal conditions were maintained over a long timespan. Sands derived from the surrounding ice-pushed ridges and transported by longshore drift and tidal currents formed a spit at the northern margin of the basin, which moved southward after eustatic sea-level rise stabilized and the lagoon was filled by clay. Loading of this clay-rich sequence by the spit and its washover fans induced subsidence, however, because of compaction, so that marine conditions were maintained until after the Eemian highstand. Fluvial and eolian sediments of Weichselian age, locally reaching a thickness of almost 10 m, eventually levelled the Amsterdam glacial basin.

This update on radiocarbon calibration results from the 19th International Radiocarbon Conference at Oxford in April 2006, and is essential reading for all archaeologists. The way radiocarbon dates and absolute dates relate to each other differs in three periods: back to 12400 cal BP, radiocarbon dates can be calibrated with tree rings, and the calibration curve in this form should soon extend back to 18000 cal BP. Between 12400 and 26000 cal BP, the calibration curves are based on marine records, and thus are only a best estimate of atmospheric concentrations. Beyond 26000 cal BP, dates have to be based on comparison (rather than calibration) with a variety of records. Radical variations are thus possible in this period, a highly significant caveat for the dating of middle and lower Paleolithic art, artefacts and animal and human remains.