Ice-Core Layer Counting

Greenland ice is counted layer by layer to 60,000 years and no further; Antarctic ice is not counted at all. Inside the counted range, up to five independent seasonal chemical signals agree with each other and with argon, uranium and tree-ring clocks to within a few per cent.

Draft:  August 2026

"It is relatively easy to count annual layers downward from the surface through considerable depths in the Greenland ice sheet." — Larry Vardiman, Institute for Creation Research, 19921

"Fundamentally, in counting any annual marker, we must ask whether it is absolutely unequivocal, or whether nonannual events could mimic or obscure a year." — Richard Alley and colleagues, GISP2 team, 19972

"In the upper halves of the Greenland ice cores, creationist and uniformitarian age assignments agree to within a factor of two or three of one another." — Jake Hebert, Institute for Creation Research, 20213

Summary

Snow that falls on Greenland in July is chemically different from snow that falls in January.  It carries different dust, different sea salt, different acidity, different oxygen isotopes.  Bury it, and if nothing smears the difference away, the ice preserves an alternating signal that can be counted like tree rings.

That is the whole method.  Everything else is a question about whether the signal survives, whether it is really annual, and how far down you can still see it.

Eight conclusions, stated first.

  1. The counting stops at 60,000 years. The Greenland Ice Core Chronology 2005 (GICC05) is continuously counted layer by layer to 60,000 years before AD 2000, and below that the chronology comes from an ice-flow model, not from counting.4 The famous "800,000 years of ice" from Antarctica is not counted at all — East Antarctic accumulation is far too low for annual layers to survive, and those ages come entirely from flow models tuned to orbital templates.3 Anyone citing 800,000 counted layers has misread the literature, and the young-earth critique is right about that.
  2. Inside the counted range, up to five physically independent signals have to agree. Sodium peaks in winter (marine air masses), calcium in spring (continental dust), ammonium in late spring (North American biomass burning), nitrate in summer, and the water isotopes trace a sine wave with summer maxima.5 These come from four different source regions and four different transport pathways. A single storm can fake one of them. Faking all five, in the right order, is a different proposition.
  3. The counters publish their own error bars, and they are not small. Ambiguous layers in GICC05 are counted as half a year plus or minus half a year, and the ambiguity is summed down the core. The result — the maximum counting error — is "on average 1% in the Holocene and 5% in the Glacial," reaching 950 years over the 41,800–60,000-year section.4 GISP2's team give roughly 1–10% in the top 2,500 m, about 20% below that, and declare the core unusable below 2,800 m.6
  4. The count has been checked against three different decay systems and lands within a few per cent of all of them. Argon-argon on the Laschamp geomagnetic excursion: 40,400 ± 1,000 years, against 41,250 ± 814 counted — 2.1%.4 Argon-argon on the North Atlantic Ash Zone II: 54,500 ± 1,000 against 55,380 ± 1,184 counted — 1.6%. Uranium-thorium on Chinese cave stalagmites, across the whole 60,000 years: agreement "within 800 years."4
  5. The chronology was caught being wrong, by seven years, and fixed. The 774–775 CE burst of solar particles is fixed to a calendar year by tree rings on five continents. Its beryllium-10 signature in Greenland and Antarctic ice sat at 768 CE on the chronologies then in use.7 One named cause was a sulfate layer wrongly attributed to Vesuvius in AD 79 and used as a fixed tie point; another, in the authors' own words, was "use of manual-layer interpretation techniques that may favor interpretations consistent with a priori knowledge or existing chronologies (WDC)."7 That is a real admission of the exact bias the young-earth critique alleges, and it should be quoted. It should also be quoted with what followed: the corrected chronologies are good to better than ±5 years over 2,500 years, and 24 of 32 independently recorded historical dust-veil years now fall within ±3 years of an ice sulfate spike, where on the old chronology there was no statistically significant association at all.7
  6. The young-earth case is that the deep "annual" cycles are storm cycles. Michael Oard's model: in a 700-year post-Flood ice age with enormous snowfall, one real year deep in the core "may be interpreted by uniformitarians as 100 or even 1,000 'annual cycles'."8 The mechanism is genuine — sub-annual oscillations do exist, and the GISP2 team say so themselves. What has never been published is a quantitative demonstration that they produce a factor of fourteen, consistently, in five signals at once.
  7. The Lost Squadron argument should be retired. P-38 fighters abandoned on Greenland in 1942 were recovered in 1992 under about 268 feet of snow and firn. The figure is real; the inference is not, because the aircraft sit near the southeast coast, the wettest part of the ice sheet, while the deep cores are drilled at the dry interior divide. Vardiman, who introduced the argument, immediately warned that "life is never as simple as implied above."1 Jake Hebert's 2018 treatment concedes the apples-to-oranges charge. Creation Ministries International's own article has carried an editorial correction for well over a decade.
  8. The strongest unanswered problem for the young-earth reading is in the upper core, not the lower. Hebert concedes that the top halves of the Greenland cores agree with creationist estimates "to within a factor of two or three."3 But the GISP2 depth–age table puts the Holocene boundary at 1,678 m and 11,651 counted years. On a post-Flood model, essentially all of that ice fell in the roughly 4,000 years after the ice age ended, under something like modern climate, where the high-snowfall over-counting mechanisms do not apply. That is an over-count of about 2.9 in the part of the core where the layers are thickest and easiest to read — and about 5.0 in the ice below 500 m if the model is granted a perfectly counted top 500 m.9 No mechanism has been offered for either.

The thing that is actually surprising

Larry Vardiman's 1992 article is the founding young-earth treatment of ice cores, and it concedes more than most later ones do.  It states that Greenland layer counting is "relatively easy… through considerable depths," that historic volcanic eruptions cross-check the count "with a high degree of accuracy," and that Hammer and colleagues had already identified about a dozen historical eruptions in a central Greenland core.1  His argument is not that the top of the core is uncountable.  It is that the bottom is, and that the deep chronology depends on models rather than counting — which, for Antarctica, is correct.


Part 1: What is actually being counted

The five signals

Snow arriving at the Greenland summit carries a chemical load that changes with the season, because its sources change with the season.  The GICC21 team list them:

  • Sodium (Na⁺) peaks in winter, "because of increased advection of marine air masses." Winter storms bring sea spray inland.
  • Calcium (Ca²⁺) peaks in spring, "because of enhanced transport from terrestrial reservoirs" — the continental dust season.
  • Ammonium (NH₄⁺) peaks in late spring or summer, "because of enhanced biogenic activity in the North American continent," including boreal wildfire.
  • Nitrate (NO₃⁻), also biogenic, peaks in summer.
  • Water isotopes (δ¹⁸O, δD) "show a sinusoidal pattern with winter valleys and summer peaks, mainly representing temperature variations at the drill site."5

To these the GISP2 programme adds two physical signals:  visual stratigraphy, in which summer depth-hoar layers scatter light differently from winter wind slab, and laser-light scattering from dust, which resolves the dust season at millimetre scale.6

The important point is not that there are five signals but that they are driven by different physics in different places.  Sea salt comes from the North Atlantic.  Dust comes from central Asia.  Ammonium comes from North American fires.  Isotopes come from the local temperature.  A blizzard that puts a spurious peak into the sodium record has no reason to put a matching spurious peak into the ammonium record from a continent away.

How a year is identified

Manually:  a scientist looks at the parallel signals plotted against depth and marks where the sequence winter–spring–summer completes.  Increasingly also automatically:  the StratiCounter algorithm treats the multi-proxy record as a hidden semi-Markov model and "learns to recognize the specific annual pattern," returning both a layer count and a probability distribution over layer boundaries.5

Both are used, and GICC21's authors are explicit that neither alone suffices:  "A layer counting bias was found in all ice cores because of site-specific signal disturbances;  therefore the manual comparison of all ice cores was deemed necessary to increase timescale accuracy."5

That sentence is worth pausing on.  It says that automatic counting is biased at every site, and that the fix is human judgement across multiple cores.  A critic is entitled to note that human judgement is exactly what the Sigl 2015 correction blamed for a seven-year error.  The answer is that the residual error is measured, and the measurement is the next section.

How the uncertainty is computed

GICC05's convention is deliberately conservative.  Where the counters cannot decide whether a feature is a year boundary, they record it as half a year, plus or minus half a year, and sum the halves down the core.  The accumulated ambiguity is the maximum counting error (MCE), treated as a 2σ bound.4

The resulting numbers, from the paper's own text and Table 1:

horizon GICC05 age (yr b2k) 1σ maximum counting error
Younger Dryas / Preboreal transition 11,703 ±50 ~99 yr (0.85%)
Onset of GI-1 (Bølling) 14,692 ±93 ~186 yr
Onset of GI-2 23,340 ±298 ~596 yr (2.6%)
Onset of GI-8 38,220 ±725 ~1,450 yr
Laschamp excursion 41,250 ±814 ~1,630 yr (3.9%)
Onset of GI-17 59,440 ±1,287 ~2,574 yr (4.3%)

Summarised in the paper's own words:  "The maximum counting error of the time scale is on average 1% in the Holocene and 5% in the Glacial."4

The authors also flag an error they cannot quantify, and it should be reported:  "There is potentially an additional uncertainty in the time scale arising from a systematic bias due to the assumed shape of annual layers…  Such a bias is, however, difficult to quantify and it is, therefore, not included in the uncertainty estimates given."4

For the last 3,800 years GICC21 does better:  its uncertainty rate is "1 year per century," and it revises GICC05 younger by about 13 years at 3,835 years ago.5


Part 2: Where the counting stops

This is the section most likely to surprise a reader who has met the subject through popular writing on either side.

Greenland, GICC05:  continuously counted to 60,000 years b2k.  Holocene from DYE-3, GRIP and NorthGRIP;  glacial from NorthGRIP alone.  Below 60,000 years, ages come from an ice-flow model.4

Greenland, GISP2:  counted to roughly 110,000 years at 2,800 m, with a stated error of about 20% between 2,500 and 2,800 m, and with the ice below 2,800 m deformed and excluded.6

Antarctica, WAIS Divide (WD2014):  annual-layer counted through 2,850 m of core, in a high-accumulation coastal-influenced setting.7

Antarctica, EPICA Dome C and Vostok:  not counted.  Accumulation at Dome C is a few centimetres of water equivalent a year;  annual layers are destroyed by diffusion long before they reach depth.  Ages come from flow models and orbital tuning.  Hebert states this correctly and it is worth quoting him:  "if accumulation rates are too low, as is the case for the deep cores in East Antarctica, direct counting of layers is not possible, so glaciologists use theoretical age-date models to assign ages to different depths within the ice."3

His further charge — that those models "implicitly assume 'millions of years' by taking the height of the ice sheet to be constant or nearly constant" — is a fair description of the standard assumption.  Whether it is circular depends on whether the model is being used to establish an old age or to interpolate within one established otherwise.  In practice the deep Antarctic chronologies are tuned to orbital and gas-record targets, which is a separate argument treated in Orbital Cycles in Rock.

None of the deep-ice ages are used as evidence in this series.  The counting evidence stops at 60,000 years, and that is sufficient for the question at hand.


Part 3: The cross-checks

A layer count could, in principle, be systematically wrong.  It could not be systematically wrong in a way that happens to track four unrelated clocks.

Against argon-argon

Two horizons in the Greenland glacial section have independent radiometric ages from lavas elsewhere:

  • The Laschamp geomagnetic excursion, a near-collapse of the dipole field recorded worldwide, dated by ⁴⁰Ar/³⁹Ar on the Laschamp and Olby lava flows at 40,400 ± 1,000 years. Counted position in NorthGRIP: 41,250 ± 814. Difference 850 years, 2.1%, inside the combined uncertainty.4
  • North Atlantic Ash Zone II, an Icelandic tephra, ⁴⁰Ar/³⁹Ar 54,500 ± 1,000; counted 55,380 ± 1,184. Difference 880 years, 1.6%.

Potassium-argon is an electron-capture decay.  Layer counting involves no decay at all.  These two have no physical reason to agree unless both are approximately right.

Against uranium-thorium

The Hulu Cave stalagmites in eastern China are dated by ²³⁰Th/²³⁴U, an alpha-decay chain with entirely different chemistry and different half-lives from potassium-argon.  Their oxygen-isotope record shows the same abrupt warmings (Dansgaard–Oeschger events) as Greenland ice.  Line the two up and:

"GICC05 compares well to the time scale of that record with absolute age differences of less than 800 years throughout the 60 ka period."4

and, on the interstadial onsets specifically:

"the two time scales agree within 500 yrs throughout the 60 ka period.  Such a good agreement would be unlikely if one of the records had a significant dating error."4

Two other caves agree — Kleegruben in Austria to within 300 years, Villars in France placing the onset of GI-12 at 46.6 ± 0.5 ka against the counted 46,860 ± 956.  One cave does not:  Moomi Cave on Socotra matches GI-12 to within 350 years but deviates by more than 1,000 years on adjacent interstadials, and the authors say so.4

Against cosmogenic isotopes and tree rings

This is the sharpest check, and it is also where the ice chronology was caught out.  It is described in Part 4.

Against history

Big eruptions leave sulfuric acid in the ice, and some of them are in the written record.  GICC21's Table 3 lines the two up:5

eruption historical date GICC21 counted date
Katmai 1912 1912 ± 2
Laki 1783 1783 ± 2
Öræfajökull 1362 1362 ± 2
Samalas 1257 1259 ± 2
Tianchi (Changbaishan) 946 946 ± 3.7
Eldgjá 939 939 ± 3.8
Okmok 43 BCE 43 ± 5.6 BCE

And in the other direction:  32 independent historical accounts of a dimmed or discoloured Sun — the oldest a Babylonian note from 255 BCE that "the disk of the sun looked like that of the moon" — were compared against volcanic years in the NEEM core.  24 of 32 fall within ±3 years, with a Monte Carlo p-value below 0.001.  On the old chronology, the same test returned no significant association at all.7


Part 4: The seven-year error, told properly

This episode is the single most useful thing in the ice-core literature for a reader trying to judge whether the field checks itself, and it should be told without spin in either direction.

The setup.  In 2012 Fusa Miyake and colleagues found that Japanese cedar rings recorded a sudden 12‰ jump in carbon-14 between 774 and 775 CE, "about 20 times larger than the change attributed to ordinary solar modulation."10  A spike that abrupt is atmospheric and therefore global.  It also produces beryllium-10, which falls out onto ice sheets.

The test.  Michael Sigl and colleagues measured beryllium-10 in Greenland and Antarctic cores and looked for the spike.  They found it — in both hemispheres, "exceeding the natural background concentration (>150%; 6σ) for one-to-two consecutive years."  On the existing chronologies it sat at 768 CE, seven years off.  A second anomaly compatible with the 994 CE event sat at 987 CE.7

The diagnosis.  Five named causes, verbatim:

"1) low effective resolution of some ice core measurements (NGRIP, GRIP);  2) use of only single (or few) parameters for annual-layer interpretation (GRIP, Dye-3);  3) intra-annual variations in various ice-core parameters falsely interpreted as layer boundaries (e.g., caused by summer melt in Dye-3);  4) use of tephra believed to originate from the 79 CE Vesuvian eruption as a fixed reference horizon to constrain the Greenland ice-core dating;  5) use of manual-layer interpretation techniques that may favor interpretations consistent with a priori knowledge or existing chronologies (WDC)."7

Causes 2, 3 and 5 are precisely the failure modes the young-earth literature has been alleging for thirty years.  They are here conceded, in Nature, by the people who built the chronologies.

A correction that has to be made.  Vesuvius AD 79 is often cited, including in some old-earth writing, as a confirmation of ice-core dating.  It is the opposite.  Sigl and colleagues rejected the attribution:  "in view of the new annual-layer dating results, we concluded that this layer dates to 87/88 CE," and note that the documentary evidence has Vesuvian ash blowing toward the eastern Mediterranean, not Greenland.7  Do not use Vesuvius as an ice-core tie point.

The repair, and what it bought.  The revised chronologies (NS1-2011 for Greenland, WD2014 for Antarctica) are constrained by the 775 CE beryllium spike plus three historically dated dust-veil years and the 1257 Samalas eruption.  They differ from the old ones "by up to 11 years (GICC05, Greenland) and 14 years (WDC06A-7, Antarctica)," and are "accurate to within less than ±5 years during the past 2,500 years."7

The validations not used in building them.  Three:  the Changbaishan tephra found in Greenland ice at 946–947 CE against a documented winter 946/47 eruption and a radiocarbon wiggle-match of 946 ± 3 CE on a tree killed by it;  the 994 CE beryllium rise appearing at 993 in the ice, matching Japanese tree-ring carbon-14 which places it between the 993 and 994 growing seasons;  and the 32 dust-veil accounts above.7

The scale of the error, in the terms of this debate.  Seven years in 1,240 is 0.6%.  If the deep layers were storms at the rate a young-earth model requires — roughly fourteen per year — the 775 CE spike would not have been seven layers deep at 1,240 counted years.  It would have been about seventeen thousand layers deep.


Part 5: The young-earth case at full strength

Vardiman 1992: the founding argument

Larry Vardiman's Impact article set the shape of the creationist position, and it is more careful than its reputation.  Its structure:

  1. The Lost Squadron aircraft were buried under ~250 feet in ~50 years, so ice can accumulate fast (raised, then immediately hedged — see below).
  2. At Camp Century, annual layer thickness runs from about fourteen inches near the surface to less than two inches near the bottom; take the mean and a 4,000-foot ice sheet forms in under 6,000 years.1
  3. Greenland counting is easy near the top and cross-checks well against historic eruptions, but confidence falls with depth; Hammer and colleagues expected ambiguity to limit counting to "less than about 8,500" layers.
  4. Antarctic ages of 160,000 years come from oxygen-isotope-to-temperature-to-accumulation-rate modelling, not from counting, and depend on the assumption that accumulation has not varied — which a post-Flood model denies. "Perhaps as much as 95% of the ice near the poles could have accumulated in the first 500 years or so after the Flood."1

Points 3 and 4 are substantially correct as descriptions of the method.  Point 2 is arithmetic on a real depth–age relation.  Point 1 is the weak one, and Vardiman knew it.

Oard 2001: the deep core is storms

Michael Oard's TJ paper is the technical core of the position.  Its argument:

  • Creationists and glaciologists agree about the top of the core. "Annual layers in the very top section of the core are easily interpreted from the ratio of oxygen isotopes, and creationists agree with these interpretations."8
  • Below that the models diverge. In a post-Flood ice age with roughly 6 m water-equivalent accumulation per year in central Greenland, an annual layer at 2,000 m depth would be about 2,000 mm thick rather than the conventional 20 mm — a factor of a hundred.8
  • The number of measurements made down the core is chosen according to the expected layer thickness, so the resolution of the instruments is set by the model. "Significantly, their interpretation of annual layers from these variables has been determined by the thickness of the annual layers that they expect, based on their model."8
  • Sub-annual oscillations demonstrably exist: a storm has warm and cold sectors, snow dunes redistribute accumulation, and multi-week warm and cold spells are routine at high latitudes.
  • His authority for this is the GISP2 team itself:

"Fundamentally, in counting any annual marker, we must ask whether it is absolutely unequivocal, or whether nonannual events could mimic or obscure a year.  For the visible strata (and, we believe, for any other annual indicator at accumulation rates representative of central Greenland), it is almost certain that variability exists at the subseasonal or storm level, at the annual level, and for various longer periodicities (2-year, sunspot, etc.).  We certainly must entertain the possibility of misidentifying the deposit of a large storm or a snow dune as an entire year or missing a weak indication of a summer and thus picking a 2-year interval as 1 year."2

That quotation is genuine and accurately reproduced.  It is also incomplete in a way that matters:  the sentences that follow it in Alley et al. 1997 go on to quantify the resulting error, at roughly 1% in the Holocene, 5% in the Wisconsinan, and under 2% against an independent uranium-thorium marker at 30,500 years.2  The GISP2 team raise the storm problem in order to bound it, not to concede that it is unbounded.

Woodmorappe 2002: two cores 28 km apart that disagree

The sharpest specific objection in the young-earth ice-core literature is not about the deep core at all.  GRIP and GISP2 were drilled 28 km apart in central Greenland.  Over a few centuries around the thirteenth millennium before present, their counted layers do not match, and the reconciliation study says so in terms:

"then the problem is not missing core or other 'block' data loss.  Rather, the GRIP core lacks about half the annual layers throughout this interval, or the GISP2 ice contains many subannual structures which mimic annual bands, or the layers are in fact annual but one of the counts is erroneous."11

Woodmorappe's conclusion:  the correlation between cores is achieved by "ad hoc contractions and expansions (i.e. 'accordioned')," and "The lack of close correspondence between the two cores, especially in large segments, in spite of their geographic proximity, favours a catastrophic storm-dominated accumulation."11

The observation is correct and it is not disputed.  What has happened since is that the field agreed with it and rebuilt the chronology.  GICC05 does not rest on GISP2 or on single-parameter counting;  it is built from DYE-3, GRIP and NorthGRIP together, with up to five parallel chemical signals, and its ambiguous layers are carried explicitly as the maximum counting error rather than resolved by fiat.  GICC21's abstract states the problem in Woodmorappe's own terms — "A layer counting bias was found in all ice cores because of site-specific signal disturbances;  therefore the manual comparison of all ice cores was deemed necessary" — and then quantifies what is left at one year per century.5

So the right reading of this episode is neither that the cores are unreliable nor that the discrepancy was imaginary.  It is that a real inter-core discrepancy of a few hundred years, in a 14,500-year record, drove a methodological rebuild whose residual error is now measured against external markers at 0.6%.

Hebert 2014–2021: resolution and circularity

Jake Hebert's contribution is the most recent and the most methodological.  His strongest single point:

"uniformitarian glaciologists increased their assigned age for the bottom of the GISP2 core from 85,000 years to 111,000 years simply by increasing the resolution of their dust measurements in the bottom 500 m of the core!"3

The fact is correct.  The inference — that the increase was driven by a desire to match the orbitally tuned marine record — is not obviously wrong and deserves a straight answer.  Here is one:  annual layers at that depth are about 4 mm thick, and the earlier dust measurements were made at 8 mm resolution.  An instrument that samples every 8 mm cannot resolve a 4 mm layer;  it must under-count, and by roughly a factor of two, which is roughly the size of the revision.  Increasing resolution to 1 mm was not a fudge, it was the necessary fix.  But the concern was legitimate and deserved that answer rather than a wave-off.

Hebert also states plainly what the young-earth programme still lacks:

"I think creationists have refuted claims that ice cores demand an old Earth.  However, we have done very little actual modelling of the ice sheets."3

There is also a mechanistic proposal that deserves noting, because it is one of the few in this literature that tries to explain an ice-core observation rather than to doubt it.  Larry Vardiman's 1997 paper argues that the long decline in the oxygen-isotope ratio through the "ice age" section of the cores need not record a hundred thousand years of ocean temperature change:  it can be produced in thousands of years by isotopic fractionation and trajectory dispersion as an ice shelf grows outward, since snow falling far from the shelf edge is isotopically lighter.12  Whether the mechanism can carry the observed amplitude has not been tested by anyone;  it is a hypothesis with a physical basis and no follow-up.

The existing model is Vardiman's 1994 one-dimensional analytic treatment, which assumes high ice-age accumulation decaying exponentially to modern rates and treats thinning as proportional to current thickness rather than solving for stress.3  Hebert has since computed annual-layer thicknesses within it and has published two further steps toward a physics-based treatment — a shallow isothermal ice dome with a frozen base in 2022, and an ice ridge in 2023.13  A 2025 paper by Kady Kunisaki takes up the sediment and ice-core records together within the same framework.14  That is honest, ongoing work, and it should be watched.  One result from it is worth reporting because it is counter-intuitive and comes from their own model:  "in creation-based ice models the thinnest annual layers are not necessarily found at the very bottom of the ice, but can occur at intermediate depths."3

The Lost Squadron

In July 1942, six P-38 fighters and two B-17 bombers made forced landings on the Greenland ice cap.  In 1992 one P-38 was recovered from beneath roughly 268 feet of snow and ice.  The argument writes itself:  five feet a year, so a 10,000-foot ice sheet takes 2,000 years.

Why it does not work.  The aircraft came down near the southeast coast, in the wettest sector of the Greenland ice sheet, where accumulation exceeds 1,500 mm of water equivalent per year — a figure mapped in the Journal of Glaciology the year before the recovery.15  The deep cores are drilled at the interior divide, which receives roughly a sixth of that.  The 268 feet is mostly firn, not ice:  about 53 m water equivalent, which is an ordinary fifty years for that location.  At the SE-Dome site 220 km away, the 1963 nuclear-test tritium horizon lies at 81 m depth — the same burial in the same fifty-odd years, dated three independent ways.16

Who declines the argument is the tell.  Vardiman raised it in 1992 and warned in the next paragraph that "life is never as simple as implied above," listing compaction, horizontal thinning, and site-to-site precipitation differences.1  Oard's technical ice-core papers do not use it.  Hebert's 2018 treatment concedes the apples-to-oranges charge and demotes the aircraft to an illustration of post-Flood snowfall rather than a refutation of ice-core dating.  Creation Ministries International's own article has carried an editorial note for well over a decade conceding that published accumulation rates make the burial depth unsurprising.17

It survives only in popular presentations.  It should be dropped.


Part 6: The problem in the upper core

The debate usually concentrates on the deep, thin, ambiguous layers.  The harder problem for a young-earth reading is in the thick, easy, unambiguous ones.

Take the GISP2 depth–age table, which is public.9  The Holocene–Pleistocene boundary sits at 1,678 m, with 11,651 counted years above it.  In an Oard-style chronology the ice age ended about 4,000 years ago;  everything above that boundary therefore fell in the last ~4,000 years, under climate not very different from today's.

That is the regime where all sides agree the counting works.  The layers are thick.  The isotope cycle is clean.  Vardiman calls counting there "relatively easy."  Hebert allows agreement "to within a factor of two or three."  But a factor of two or three is not enough:  11,651 counted years in ~4,000 real years is a factor of about 2.9, and it has to be achieved in ice that fell after the high-snowfall, high-storm conditions that the model invokes to explain over-counting had ended.

And the requirement gets worse if the model is given the part it most obviously deserves.  The same table puts 500 m depth at 2,097 counted years.  Grant that the top 500 m is counted perfectly — recent ice, thick layers, historic volcanic markers all through it — and the remaining 9,554 counted years have to fit into 4,000 − 2,097 = 1,903 real years.  That is a required over-count of 5.0 in the ice between 500 m and 1,678 m.9

More snowfall does not help.  Thicker layers are easier to resolve, not harder.  The mechanism the model needs — many storm cycles miscounted as years — works best where layers are thin and ambiguous, which is exactly where this ice is not.

This is, in the assessment of this article, the strongest single unanswered objection to the young-earth reading of the ice cores, and it is computable from a table anyone can download.


Objections and responses

"Antarctic ice is 800,000 years old, so the earth is old." Not by counting it is not.  East Antarctic accumulation is too low for annual layers to survive to depth;  those ages come from flow models and orbital tuning.  The young-earth critique is correct on this point and it should be conceded plainly.  The counting evidence is Greenland, and it stops at 60,000 years.

"Layers can be counted wrong, and the team admits it." They do admit it, they quantify it, and the quantification is the point.  Ambiguous layers are counted at half weight;  the accumulated ambiguity is 1% in the Holocene and 5% in the glacial;  and the residual error against external clocks has been measured at 0.6% over the last 1,240 years and about 2% against argon-argon at 41,000.  A method whose error is measured is not thereby refuted.

"A single storm can produce a false layer." Yes, in one signal.  The counting requires the signals to agree, and they come from different continents and different physical processes.  A North Atlantic storm has no mechanism for putting a spurious ammonium peak into snow whose ammonium comes from North American wildfire.

"Snowfall was far higher after the Flood, so the deep layers are storms." This is the real argument, and it needs a number.  Sixty thousand counted layers in about 4,300 years requires roughly fourteen countable oscillations per year, sustained, in five parallel signals, for the entire depth of the core.  No published creationist model produces that figure, and Hebert states that ice-sheet modelling on the creationist side has barely begun.

"The count was raised from 85,000 to 111,000 years to match the marine record.  That is circular." The revision followed an eightfold increase in dust-measurement resolution, from 8 mm to 1 mm, in a section where layers are about 4 mm thick.  The original measurement physically could not resolve those layers.  The suspicion was reasonable;  the explanation is mechanical.

"The Vesuvius layer proves the chronology works." It does not, and this claim should be dropped by whoever is making it.  The Vesuvius attribution was withdrawn in 2015;  the layer dates to 87/88 CE, and using it as a tie point was one of the named causes of the seven-year error.

"They tuned the ice chronology to the tree rings, so of course they agree." The direction of the correction refutes this.  The chronology was not tuned to the tree rings before 2015 — that is why it was seven years wrong.  The tree-ring marker was applied from outside, found a discrepancy, and forced a revision.  A tuned system does not do that.

"The Lost Squadron planes were under 268 feet after fifty years." They were, at a site receiving more than 1,500 mm water equivalent a year.  The cores are drilled where accumulation is about a sixth of that.  The argument has been withdrawn or hedged by every creationist author who has examined it closely.


Issues and unanswered questions

  • The "shape of annual layers" bias is unquantified, by the counters' own statement, and is not included in the stated uncertainties.4
  • Around 30,000 years the Greenland scale carries about 1,000 uncertain years, and around 40,000 there is an unresolved conflict of about 1,000 years between a floating New Zealand kauri sequence placed by beryllium-10 against the ice and the same sequence placed by carbon-14 against Hulu Cave. The specialists conclude the ice placement is the wrong one.18
  • GICC05 systematically over-counts the Holocene by about 65 years — 0.55%, in the direction of too many layers, not too few.18 Two independently counted chronologies in opposite hemispheres differ by 0.63% in accumulated years, and IntCal20 applies a blunt ×1.0063 correction to reconcile them.
  • Moomi Cave, Socotra, does not fit. It agrees on one interstadial and deviates by more than a thousand years on the adjacent ones. No explanation has been published.
  • Alley et al. 1997 has not been read in full for this article. The quotation Oard uses and the following sentences that quantify the error are verified; the rest of the paper has not been obtained directly.

Verdict

Ice-core layer counting is a good clock over a limited range and a poor one outside it, and both halves of that sentence get lost in popular argument.

The young-earth critique gets three things right.  Antarctic deep ice is not counted.  Greenland counting is genuinely difficult at depth and its error was genuinely underestimated.  And the field has, on the record, tuned chronologies to markers that turned out to be misidentified.

It gets the conclusion wrong for a reason that can be stated in one line:  the residual error, after all of that, has been measured against three unrelated decay systems and one absolutely dated cosmic-ray event, and it comes out between 0.6% and 5%.  A young-earth chronology needs it to be about 1,300%.

And the sharpest problem is not in the disputed deep ice at all.  It is in the top 1,678 metres, where everyone agrees the counting works, and where a post-Flood model still has to explain 11,651 layers in about 4,000 years.


Further reading


Sources


  1. 1.

    Vardiman, Larry.  “Ice Cores and the Age of the Earth.”  Impact.  Institute for Creation Research.  Apr 1992.  Article No. 226 in ICR's IMPACT series, by Larry Vardiman, Ph.D.  It argues that the great ages given to the bottom layers of the ice cores depend on models of past climate, not on direct counting.  It concedes that Greenland layers can be counted relatively easily through considerable depths, and it proposes that snowfall right after the Flood could have built the ice sheets quickly.  The article uses it as the founding young-earth treatment of ice cores, one that concedes more than most later ones do.

    1. 1a.

      Section Greenland Ice Cores:  “Although occasional ambiguities occur, it is relatively easy to count annual layers downward from the surface through considerable depths in the Greenland ice sheet.  This is possible because of the large precipitation rates in Greenland and the preservation of the annual effects.  It is also possible with a high degree of accuracy to cross check the counting of annual layers with occasional peaks in acidity and particulates from the fallout of historic volcanic events.” ↩ ↩

    2. 1b.

      Section Greenland Ice Cores:  Said immediately after presenting the Lost Squadron arithmetic:  “However, life is never as simple as implied above.  In making our calculations, we did not take into account the compaction of the snow into ice as it is weighted down by the snow above.  Neither did we consider the thinning of ice layers as the tremendous weight above forces the ice at lower levels to squeeze out horizontally.  More importantly, we did not consider the average precipitation rate and actual depths of ice for different locations on the Greenland ice sheet.” ↩ ↩

    3. 1c.

      Section Greenland Ice Cores:  “the average annual thickness of ice at Camp Century located near the northern tip of Greenland is believed to vary from about fourteen inches near the surface to less than two inches near the bottom (Hammer, et al., 1978).  If, for simplicity, we assume the average annual thickness to be the mean between the annual thickness at the top and at the bottom (about eight inches), this still gives an age of less than 6000 years for the 4000-foot-thick ice sheet.” ↩

    4. 1d.

      Section Antarctic Ice Cores:  “some creationist models predict significant quantities of snow immediately after the Flood (Oard, 1990).  Perhaps as much as 95% of the ice near the poles could have accumulated in the first 500 years or so after the Flood.” ↩

  2. 2.

    Alley, R. B. et al.  “Visual-stratigraphic dating of the GISP2 ice core: basis, reproducibility, and application.”  Journal of Geophysical Research 102(C12):26367–26381.  1997.  A GISP2 team paper on visual-stratigraphic dating, which counts the years in the ice core from the layers seen by eye, and on how well that method can be reproduced and applied.  The article uses it for the team's own limits on counting annual layers, including storms that can be mistaken for a year.

    1. 2a.

      Page 26378:  “Fundamentally, in counting any annual marker, we must ask whether it is absolutely unequivocal, or whether nonannual events could mimic or obscure a year.  For the visible strata (and, we believe, for any other annual indicator at accumulation rates representative of central Greenland), it is almost certain that variability exists at the subseasonal or storm level, at the annual level, and for various longer periodicities (2-year, sunspot, etc.).  We certainly must entertain the possibility of misidentifying the deposit of a large storm or a snow dune as an entire year or missing a weak indication of a summer and thus picking a 2-year interval as 1 year.” ↩ ↩

    2. 2b.

      The sentences that follow the passage quantify the resulting error at approximately 1% in the Holocene, 5% in the Wisconsinan, and under 2% against an independent uranium-thorium marker at 30,500 years. ↩

  3. 3.

    Hebert, Jake.  “The current state of creationist ice core research.”  Journal of Creation 35(3):12–14.  Dec 2021.  A short review by Jake Hebert, a young-earth writer on ice cores, in the creationist Journal of Creation.  It argues that thick ice sheets can form in thousands of years, that the deep Antarctic ages rest on models, and that creationist and mainstream ages for the top halves of the Greenland cores agree within a factor of two or three.  It also admits that creationists have done very little actual modelling of ice sheets, and the article uses it for those concessions.

    1. 3a.

      “Interestingly, in the upper halves of the Greenland ice cores, creationist and uniformitarian age assignments agree to within a factor of two or three of one another.  Only in the bottom halves of the Greenland cores do uniformitarian and creationist age estimates dramatically diverge.” ↩ ↩

    2. 3b.

      “if accumulation rates are too low, as is the case for the deep cores in East Antarctica, direct counting of layers is not possible, so glaciologists use theoretical age-date models to assign ages to different depths within the ice.  These models attempt to take into account thinning of the ice…  These models implicitly assume 'millions of years' by taking the height of the ice sheet to be constant or nearly constant.” ↩ ↩

    3. 3c.

      Hebert, on the age given to the bottom of the GISP2 core:  “This earlier pessimism is more realistic, as demonstrated when uniformitarian glaciologists increased their assigned age for the bottom of the GISP2 core from 85,000 years to 111,000 years simply by increasing the resolution of their dust measurements in the bottom 500 m of the core!” ↩

    4. 3d.

      “… I think creationists have refuted claims that ice cores demand an old Earth.  However, we have done very little actual modelling of the ice sheets.” ↩

    5. 3e.

      “In 1993 Larry Vardiman of the Institute for Creation Research published a one-dimensional mathematical model for the rapid formation of an ice sheet in the time since the Genesis Flood.  This model assumed high Ice Age accumulation rates that exponentially decayed to today's 'slow and gradual' rates.  This model did not directly calculate the stresses on a parcel of ice but implicitly assumed that thinning of the ice at any time was proportional to the current thickness of the ice sheet.”  Hebert's own extension of this model is CRSQ 57(3):175–185 (2021);  he reports working toward the higher-order Blatter–Pattyn model used in the Community Ice Sheet Model. ↩

    6. 3f.

      “One interesting result is that in creation-based ice models the thinnest annual layers are not necessarily found at the very bottom of the ice, but can occur at intermediate depths.”  Hebert describing his own CRSQ 57(3):175–185 computation. ↩

  4. 4.

    Svensson, A. et al.  “A 60 000 year Greenland stratigraphic ice core chronology.”  Climate of the Past 4:47–57.  2008.  A dating study that builds the Greenland Ice Core Chronology 2005 (GICC05), a timescale made by counting annual layers (yearly bands) in ice cores from DYE-3, GRIP and NorthGRIP out to 60,000 years.  It checks the count against a Chinese cave record (Hulu Cave), which agrees to within 800 years over 60,000 years, and against argon-argon dates, a radiometric clock, for two reference layers.  The article uses it for the counting limit and for the counters' own error bars.

    1. 4a.

      The chronology is continuously counted to 60,000 years before AD 2000: Holocene from DYE-3, GRIP and NorthGRIP ice cores, glacial from NorthGRIP alone.  NorthGRIP layer thicknesses over 41.8–60 ka are.  Page 50:  “1–1.5 cm thick in cold glacial stadials and 1.5–2.5 cm thick in mild interstadials.”  The 41.8–60 ka section was counted independently by two authors, whose counts agree within the stated error.  Page 49:  “although locally the difference sometimes exceeds the final counting error.” ↩ ↩

    2. 4b.

      The article's table gives 1σ ages taken from the paper's Table 1, and its MCE column is twice the 1σ value.  Page 49:  “The accumulated error obtained by summing up the uncertain annual layers is called the Maximum Counting Error (MCE) and is regarded as a 2σ error of the time scale (MCE = 2σ).”  Page 50:  “The MCE for the 41.8–60 ka section is 950 years or 5% which is similar to that of the 14.8–41.8 ka section.”  Page 54:  “The maximum counting error of the time scale is on average 1% in the Holocene and 5% in the Glacial.” ↩ ↩ ↩

    3. 4c.

      Page 49:  “There is potentially an additional uncertainty in the time scale arising from a systematic bias due to the assumed shape of annual layers…  Such a bias is, however, difficult to quantify and it is, therefore, not included in the uncertainty estimates given.” ↩ ↩

    4. 4d.

      Table 2:  The Laschamp event, a brief change in Earth's magnetic field: GICC05 gives 41,250 ± 814 years before AD 2000, at NorthGRIP depth 2118.66 m, against 40,400 ± 1,000 years by argon-argon dating (Guillou et al. 2004).  North Atlantic Ash Zone II: GICC05 gives 55,380 ± 1,184 years, at depth 2359.45 m, against 54,500 ± 1,000 years (Southon 2004).  The ¹⁴C radiocarbon ages in the same table are uncalibrated and must not be compared directly with counted ages; the paper's Figure 7 plots them against the calibration curves instead. ↩ ↩

    5. 4e.

      Section Abstract:  “Assuming that the Greenland climatic events are synchronous with those seen in the Chinese Hulu Cave speleothem record, GICC05 compares well to the time scale of that record with absolute age differences of less than 800 years throughout the 60 ka period.”  Page 52:  On interstadial onsets, excluding GI-12:  “the two time scales agree within 500 yrs throughout the 60 ka period.  Such a good agreement would be unlikely if one of the records had a significant dating error.”  Concluding section:  “The general agreement between the new Greenland chronology and many other independent time scales and reference horizons excludes the possibility of a significant hiatus in the Greenland ice cores and it supports that GICC05 has no significant long-term bias.” ↩ ↩ ↩

    6. 4f.

      Kleegruben Cave, Austria: GI-14 and GI-15 match GICC05 to.  Page 52:  “less than 300 yrs difference.”  Villars Cave, France: onset of GI-12 constrained to 46.6 ± 0.5 ka against GICC05's 46,860 ± 956, described as.  “excellent agreement.”  Moomi Cave, Socotra (core M1-2): GI-12 agrees within 350 yr but adjacent interstadials.  “deviate more than 1 ka from GICC05.” ↩

  5. 5.

    Sinnl, Giulia et al.  “A multi-ice-core, annual-layer-counted Greenland ice-core chronology for the last 3800 years: GICC21.”  Climate of the Past 18:1125–1150.  2022.  A revision of the Greenland ice-core timescale GICC05 for the last 3,835 years, built by synchronizing six deep and three shallow ice cores from the central Greenland ice sheet and counting their yearly layers.  It combines automatic counting on several chemical signals with manual checking, and it finds the new timescale younger than GICC05 by about 13 years at 3,835 years ago, with an uncertainty of about 1 year per century.  The article uses it for the counting method and for its finding that every core had a counting bias.

    1. 5a.

      Section 1.1:  “sodium (Na⁺) has a strong winter peak because of increased advection of marine air masses…  Relative to sodium, calcium (Ca²⁺) peaks in the spring because of enhanced transport from terrestrial reservoirs…  ammonium (NH₄⁺) has a maximum in the late spring or the summer because of enhanced biogenic activity in the North American continent;  nitrate (NO₃⁻), which is also related to biogenic processes, peaks in the summer…  Water isotopes (δ¹⁸O and δD) show a sinusoidal pattern with winter valleys and summer peaks, mainly representing temperature variations at the drill site.” ↩ ↩

    2. 5b.

      Section 1.2 Annual-layer counting methods:  “StratiCounter (SC) is a [program that identifies] annual layers in an ice-core multi-proxy dataset…  Based on a hidden semi-Markov model, the algorithm learns to recognize the specific annual pattern.  SC provides a layer count and a probability distribution of the recognized layer boundaries.” ↩

    3. 5c.

      Section Abstract:  “A layer counting bias was found in all ice cores because of site-specific signal disturbances;  therefore the manual comparison of all ice cores was deemed necessary to increase timescale accuracy.” ↩ ↩

    4. 5d.

      Section Abstract:  “After examining sources of error and their correlation lengths, the uncertainty rate was quantified to be 1 year per century.  The new timescale is younger than GICC05 by about 13 years at 3835 years ago.  The most recent 800 years are largely unaffected by the revision.” ↩

    5. 5e.

      Table 3:  Historically dated volcanic eruptions and cosmogenic events against GICC21.  Also listed: the 994 CE solar-particle event → 992 ± 3.6; the 775 CE event → 774 ± 4.1; "UE 88" (the layer formerly attributed to Vesuvius) → 89 ± 5.4; the 660 BCE event → 663 ± 6.8 BCE; Aniakchak (formerly attributed to Thera) → 1629 ± 7.3 BCE. ↩

  6. 6.

    Meese, D. A. et al.  “The Greenland Ice Sheet Project 2 depth-age scale: methods and results.”  Journal of Geophysical Research 102(C12):26411–26423.  1997.  The GISP2 team's paper on the depth-age scale of the GISP2 Greenland ice core, which links each depth in the core to an age, and on how that scale was built and what it gave.  The article uses it for the team's own error statements and for the signals the team counted.

    1. 6a.

      The GISP2 team's own uncertainty statement: roughly 1–10% in the top 2,500 m, about 20% between 2,500 and 2,800 m, and ice below 2,800 m deformed and excluded from the chronology. ↩ ↩

    2. 6b.

      The GISP2 count uses visual stratigraphy (summer depth hoar against winter wind slab), laser-light scattering from dust, electrical conductivity, and major-ion chemistry, with oxygen isotopes usable only in the upper part of the core. ↩

  7. 7.

    Sigl, M. et al.  “Timing and climate forcing of volcanic eruptions for the past 2,500 years.”  Nature 523(7562):543–549.  2015.  A Nature paper that dates volcanic eruptions over the past 2,500 years from Greenland and Antarctic ice cores, then checks the new dates against historical records and independent age markers.  It measures volcanic aerosol (tiny particles that eruptions put into the air) in the ice to date each eruption, and finds that large tropical and high-latitude eruptions were the main drivers of Northern Hemisphere temperature swings, with cooling lasting up to ten years after the largest.  The article uses it for the seven-year correction to the Greenland chronology and for the validation tests.

    1. 7a.

      Lines 94–103:  “Both polar ice sheets contain ¹⁰Be concentrations exceeding the natural background concentration (>150%; 6σ) for one-to-two consecutive years, compatible with the 775 CE event observed in tree rings.  Using the original ice-core age models, the ages of the ¹⁰Be maxima in NEEM-2011-S1, NGRIP, and WDC are 768 CE, offset by 7 years from the tree-ring event.  A further ¹⁰Be anomaly measured in NEEM-2011-S1 at 987 CE, compatible with the 994 CE event in tree rings, suggests a chronological offset was present by the end of the first millennium CE.” ↩ ↩

    2. 7b.

      Lines 598–606:  “use of manual-layer interpretation techniques that may favor interpretations consistent with a priori knowledge or existing chronologies (WDC).” ↩ ↩

    3. 7c.

      “In NS1-2011 we did not attribute acid layers to the historical eruptions Vesuvius 79 and Hekla 1104, due to a lack of corroborative tephra at these depths in this and a previous study.  Possible Vesuvian tephra was reported from the Greenland Ice Sheet Project (GRIP) ice core at 429.3 m depth, but in view of the new annual-layer dating results, we concluded that this layer dates to 87/88 CE…  Documentary sources also suggest that the main vector of ash transport following the Vesuvius 79 CE eruption was toward the eastern Mediterranean.” ↩

    4. 7d.

      The WAIS Divide WD2014 chronology is annual-layer counted to 2,850 m.  “Compared to the previous chronologies, age models differ prior to 1250 CE by up to 11 years (GICC05, Greenland) and 14 years (WDC06A-7, Antarctica).”  “We conclude that the revised ice-core timescales are accurate to within less than ±5 years during the past 2,500 years.” ↩ ↩

    5. 7e.

      “Using the new timescales, we found large sulfate signals in Greenland (e.g. 682, 574, 540 CE) between 500 and 2000 CE that frequently occurred within ±1 year from comparable – and independently dated – signals in Antarctica…  Back to 400 BCE, other large sulfate peaks (e.g., 44 BCE) were synchronous to within ±3 years.”  Thirty-two independent historical dust-veil observations back to 255 BCE; the Babylonian record reading.  “the disk of the sun looked like that of the moon”.  24 of 32 fall within ±3 yr of a NEEM volcanic year, Monte Carlo p < 0.001 over 10⁶ iterations.  “When placing volcanic event years on the original GICC05 timescale, we did not observe any statistically significant association with our independent validation points.” ↩ ↩

    6. 7f.

      Validation markers not used in constructing the chronology: Changbaishan/Tianchi tephra (volcanic ash) in NEEM-2011-S1 at 946–947 CE, against documentary evidence of a winter 946/47 eruption and a high-precision ¹⁴C wiggle-match of 946 ± 3 CE on a tree killed by it; the 994 CE ¹⁰Be rise appearing at 993 CE in the ice, consistent with Japanese tree-ring Δ¹⁴C placing the rise between the 993 and 994 growing seasons. ↩

  8. 8.

    Oard, Michael J.  “Do Greenland ice cores show over one hundred thousand years of annual layers?”  TJ (now Journal of Creation) 15(3):39–42.  2001.  A 2001 article by Michael J. Oard, a creationist who writes on the post-Flood Ice Age, disputing the claim that glaciologists count 110,000 annual layers down the GISP2 core.  GISP2 was drilled in central Greenland between 1989 and 1993, to bedrock at about 3,000 m.  Oard argues that the deep 'annual cycles' are variations within a single year, and that the counts follow the layer thickness the model expects.  The article uses it to show that creationists and glaciologists agree about the top of the core and disagree below it.

    1. 8a.

      Page 39:  The paper's abstract says:  “Annual layers in the very top section of the core are easily interpreted from the ratio of oxygen isotopes, and creationists agree with these interpretations.  Below this top section, the annual layers interpreted from the two models diverge significantly.” ↩

    2. 8b.

      Page 39:  “one annual layer deep in the ice sheet may be interpreted by uniformitarians as 100 or even 1,000 'annual cycles'.  However, instead of annual cycles, the oscillations simply represent variations within a single year.”  At 2,000 m in GISP2:  “the assumed annual layer thickness from the uniformitarian model would be about 20 mm, whereas the creationist estimate would be around 2,000 mm, 100 times more.” ↩ ↩

    3. 8c.

      “The assumed thickness of the annual layers is important because it determines how many measurements of each variable are made down the length of the core…  Significantly, their interpretation of annual layers from these variables has been determined by the thickness of the annual layers that they expect, based on their model.” ↩

  9. 9.

    NOAA National Centers for Environmental Information.  GISP2 Ice Core 1997 Layer-Counted Timescale.  1997.  A NOAA depth-age table for the GISP2 Greenland ice core (study 17835), with 11,919 rows down to 3,037 m, plus the GISP2 visual-stratigraphy log.  The article reads two rows from it: a depth of 1678 m with 11651 counted years at the Holocene boundary, and 499.500 m with 2097 counted years at 500 m.  The article uses them to test whether a post-Flood chronology could fit the counted layers into the last few thousand years. ↩ ↩ ↩

  10. 10.

    Miyake, Fusa et al.  “A signature of cosmic-ray increase in ad 774–775 from tree rings in Japan.”  Nature 486:240–242.  2012.  A study of tree rings from Japan that reports a signature of a cosmic-ray increase (a rise in high-energy particles from space) in ad 774–775.  The article uses it as the first step in showing that the spike was global.

    1. 10a.

      “a rapid increase of about 12‰ in the ¹⁴C content from ad 774 to 775, which is about 20 times larger than the change attributed to ordinary solar modulation.” ↩

  11. 11.

    Woodmorappe, John.  “Greenland ice cores: implicit evidence for catastrophic deposition.”  Journal of Creation 16(3):14–16.  Dec 2002.  A young-earth article in the creationist Journal of Creation by John Woodmorappe, who argues that one Greenland core must be given ad hoc contractions and expansions to fit the other, which he says favours catastrophic storm-dominated accumulation.  It is the sharpest specific objection in the young-earth ice-core literature, and it concerns the top of the cores rather than the deep section, which is why the article uses it.

    1. 11a.

      Woodmorappe quoting the GRIP/GISP2 reconciliation study he cites as his reference 2, on a few centuries of disagreement around the thirteenth millennium BP:  “then the problem is not missing core or other 'block' data loss.  Rather, the GRIP core lacks about half the annual layers throughout this interval, or the GISP2 ice contains many subannual structures which mimic annual bands, or the layers are in fact annual but one of the counts is erroneous.”  The passage is quoted as Woodmorappe quotes it. ↩

    2. 11b.

      Woodmorappe's conclusion:  “The variance in one core must be subject to ad hoc contractions and expansions (i.e. 'accordioned') in order to force it to fit the variance of other core.  The lack of close correspondence between the two cores, especially in large segments (blocs), in spite of their geographic proximity, favours a catastrophic storm-dominated accumulation of the water material in the Greenland ice cap.”  He also notes the GISP2 visual-layer count carries.  “constant, fine-scale counting uncertainties in the 1–2% range”.  and short core-loss gaps of usually under ten years between 3,000 and 9,000 years BP whose missing layers are interpolated. ↩

  12. 12.

    Vardiman, Larry.  “Rapid Changes in Oxygen Isotope Content of Ice Cores Caused by Fractionation and Trajectory Dispersion near the Edge of an Ice Shelf.”  Creation Ex Nihilo Technical Journal 11(1):52–60.  1997.  A 1997 paper in the creationist Creation Ex Nihilo Technical Journal by Larry Vardiman.  It proposes that the fall in ice-core oxygen-isotope ratios during the Ice Age can come from an ice shelf growing outward, through isotope fractionation and the dispersion of snow, and that this can happen in thousands of years.  The article uses it as one of the few mechanistic proposals in this literature, one that tries to explain an ice-core observation rather than doubt it.

    1. 12a.

      Section Abstract:  “The growth of ice shelves during the 'Ice Age' is shown to cause a decreased isotopic ratio at long distances from the edge of an ice shelf because of the fractionation of isotopes as a function of the vertical temperature distribution in the atmosphere and the dispersion of snow by crystal type, fall velocity, and wind fields.  If ice shelves grew slowly during the 'Ice Age' and melted rapidly during deglaciation, the trend observed in ice cores can be explained in thousands of years, consistent with a short interpretation of earth history.” ↩

  13. 13.

    Hebert, Jake.  “Towards a More Realistic Young-Earth Ice Sheet Model: A Shallow, Isothermal Ice Dome with a Frozen Base.”  Creation Research Society Quarterly.  2022.  “…A Shallow, Isothermal Ice Ridge with a Frozen Base”.  Two papers by the young-earth writer Jake Hebert, both attempts at a more realistic young-earth model of an ice sheet.  The first, from 2022, treats a shallow ice dome that is one temperature throughout and frozen to its base.  The second, from 2023, applies the same approach to a shallow, isothermal ice ridge with a frozen base.  The article cites them as the next steps toward a physics-based treatment of how creationist ice sheets form, and as work in progress. ↩

  14. 14.

    Kunisaki, Kadyn.  “Post-Flood Paleoclimate: Reexamining the Sediment and Ice Core Records in Light of the Global Flood.”  Creation Research Society Quarterly 62(2):89–101.  2025.  A 2025 paper in the creationist Creation Research Society Quarterly that reads published deep-sea sediment and Antarctic ice-core records of the Cenozoic (the last 65 million years) on a Biblical timeline, placing the Flood about 4,375 years ago.  It reinterprets those records as post-Flood cooling, and its disclaimer says, citing Mike Oard, that creation scientists have not thoroughly examined the deep-sea sediments it uses.  The article cites it as the most recent young-earth treatment of the ice-core and marine-sediment records together. ↩

  15. 15.

    Ohmura, Atsumu, and Niels Reeh.  “New precipitation and accumulation maps for Greenland.”  Journal of Glaciology 37(125):140–148.  1991.  A glaciology paper with maps of snow accumulation across the Greenland ice sheet, measured in millimetres of water equivalent (the depth of water the snow would make if melted).  It maps accumulation of over 1,500 mm of water equivalent per year in the southeast, against roughly a sixth of that at the interior divide where the deep cores are drilled.  The article uses it for that gradient, and the maps were published the year before the Lost Squadron was recovered. ↩

  16. 16.

    Iizuka, Yoshinori et al.  “A 60 year record of atmospheric aerosol depositions preserved in a high-accumulation dome ice core, Southeast Greenland.”  Journal of Geophysical Research: Atmospheres 122.  2017.  A 2017 study of a 60-year record of airborne particles preserved in the SE-Dome ice core, a high-accumulation dome in southeast Greenland.  The core lies 220 km from the Lost Squadron site, and its 1963 nuclear-test tritium horizon lies at 81 m depth.  The article uses it for that comparison: the same order of burial in the same fifty-odd years, dated three independent ways. ↩

  17. 17.

    Wieland, Carl.  “The lost squadron.”  Creation 19(3):10–14.  Jun 1997.  The popular creationist presentation of the Lost Squadron argument, by Carl Wieland of Creation Ministries International, in the magazine Creation.  It says the planes were buried under about 250 feet of ice that built up in less than 50 years, which it takes as a sign that ice can build up quickly.  The article uses it because the page has carried an editorial note since at least 2009 conceding that published accumulation rates make the burial depth unsurprising, and that concession comes from the same publisher. ↩

  18. 18.

    Muscheler, Raimund et al.  “Testing and Improving the IntCal20 Calibration Curve with Independent Records.”  Radiocarbon 62(4):1079–1094.  2020.  A 2020 research paper that tests the new IntCal20 radiocarbon calibration curve against records dated in other ways, such as volcanic eruptions dated by ice-core layer counting or by argon-argon dating.  It finds that these records largely support the curve, and its authors include Paula Reimer and Timothy Heaton, who also write the companion papers on IntCal20 in the same issue.  The article uses it for three limits on the ice-core clock:  about 1,000 uncertain years on the Greenland scale near 30,000 years ago, an unresolved conflict of about 1,000 years near 40,000, and a Holocene over-count of about 65 years.

    1. 18a.

      Around 40,000 years a floating New Zealand kauri sequence placed by matching its ¹⁰Be record to the ice cores comes out about 1,000 years younger than its placement against Hulu Cave ¹⁴C.  The authors conclude:  Page 1088:  “the ice core placement of the floating kauri sequence is likely too young.”  The paper also states:  “The Greenland ice-core time scale (GICC05) contains about 1000 uncertain years around 30,000 calendar yr BP.” ↩

    2. 18b.

      Page 1089:  “It indicates that GICC05 has the tendency to accumulate too many years (about 65) over the Holocene period.”  Sixty-five years over ~11,700 is 0.55%, and the sign is over-counting.  IntCal20 plots NGRIP δ¹⁸O with the GICC05 timescale.  “multiplied by 1.0063 as a first order correction to the offset between GICC05 and the West Antarctic Ice sheet (WAIS) Divide WD2014 chronology”.  Two independently counted chronologies, one per hemisphere, differ by 0.63%. ↩