How old is the Earth?
Although orthodox science is confident that the Earth originated 4.56 billion years ago, the planet has no birth certificate. In order to date a rock one has to have a dating system, and any such system depends on a dating theory. The theory, in turn, needs to be tested against the things being dated. A coroner may perform an autopsy and deduce that the corpse before him is of a woman in her 70s, but if her whole appearance suggests she died in her teens, we should suspect there is something wrong with his dating method. Common sense plays a role.
To date a rock we must identify something in it that has changed over time at a regular rate. This is what radioisotope dating does. Certain elements are radioactive: they emit particles and decay into elements whose atomic number is consequently lower. Though we cannot predict exactly when a particular atom will emit a particle, overall the decay takes place at a statistically constant rate, and if we can deduce what the parent-daughter ratio was at the time the rock crystallised, we can calculate how much time must have elapsed to account for the parent-daughter ratio now. Obviously enough, the method applies only to rocks that contain the products of radioactive decay.
It also applies only if the rate of radioactive decay has always been the same. That is a crucial question. Suppose the rate had gradually decreased over Earth history, tailing off to zero. Would we be able to distinguish between these possibilities? As the currently measured value of the decay rate (or half-life) of an element has no theoretical basis, this can only be answered by testing the method against the primary evidence. It needs a ‘reality check’, much as the conclusion of our imaginary coroner needed a reality check. Unfortunately, this is never attempted. If a geological sequence suggests that it built up more rapidly than present-day decay rates suggest, the dating system is not questioned. Instead, the sequence is assumed to include long intervals which left no geological record. The reasoning is essentially circular.
The age of the Earth is a scientific question. Thus neither religious dogma (”biblical genealogies require a short timespan”) nor atheistic dogma (”evolution requires a long timespan”) should be allowed to dictate an answer. It requires the same approach as when considering to what extent organisms have evolved. If the fossil record shows species descending from other species by innumerable fine gradations and all converging upon a common ancestor, that is what we should conclude. If it shows only limited evolution, we should not pretend it shows something else. Science is not about imposing on the evidence interpretations that make up for its not showing what we hoped for.
Milankovitch cycles
Geology comes closest to verifying radioisotope dating when its results are compared to the time taken for Earth to complete certain astronomical cycles. While the association is speculative, it is thought the cycles may generate long-term swings in climate strong enough to be expressed in sedimentary sequences. Thus if a sequence has a pattern consistent with the repetition of astronomical cycles, it may be possible to calculate its entire length. The approach is similar to tree-ring dating. Working back from the present, the aim is to link such sequences into one master sequence, and compare the astronomical timescale with the one based on radioisotope decay.
The shortest astronomical cycles are the day, month, and year. In the case of the year, the cycle is defined by the seasons. With the Earth’s axis being at an angle relative to its orbit round the Sun, a hemisphere experiences summer when it is tilted towards the Sun and winter when it is tilted away. The seasons constitute the primary astronomical cycle. Three others affect the intensity of the seasons.
The Earth’s tilt – known as its obliquity – is currently 23.5 degrees. The greater the tilt, the greater the difference between summer and winter. Over the time that accurate measurements have been possible the tilt has decreased at a rate of 1.19 metres per century. It is thought to oscillate between a minimum of around 22 degrees and a maximum of 24.5 degrees. A complete cycle is calculated to take 41,000 years, with the period very slowly lengthening over time.
The shape of the Earth’s orbit round the Sun also varies, from a near-perfect circle to an ellipse. This pattern of change – called the eccentricity cycle – extends over 98,000 years. When the orbit is markedly elliptical, summers which coincide with the shortest distance between the Sun and the Earth will be warmer than those which coincide with the longest distance.
Finally, there is the phenomenon known as precession, where the elliptical orbit of the Earth turns around itself, and the Earth’s axis wobbles around itself, like the axis of a spinning top. Precession of the Earth’s orbit increases the seasonal contrast in one hemisphere while decreasing it in the other. At this present stage in the cycle the Earth is closest to the Sun when the northern hemisphere is in winter, so northern winters are less cold and summers less warm. Precession of the Earth’s axis (thought to be caused by the pull of the Moon and the Sun on the equatorial bulge) causes a corresponding shift in the celestial poles. These days the North Star is in the constellation Ursa Minor, but 5,000 years ago the North Star was in the constellation Draco. The precession cycle varies between 19,000 and 23,000 years – again, slowly increasing as time passes.
The three cycles of precession, obliquity and eccentricity are known as Milankovitch cycles, after the scientist who proposed that these long-term variations have a noticeable effect on climate. Milutin Milankovitch published his theory in 1941, but it was not until the 1980s that geologists began applying it to sedimentary sequences. These seemed to bear him out. Where deposits were stratified in couplets and sedimentation seemed to have been fairly continuous, one couplet could often be interpreted as indicating a warm climate and the other a cool climate.
An example is chalk sequences. Chalk consists predominantly of detritus from calcareous warm-water algae known as coccolithophores. Although from a distance it can look like a monotonous mass of white or grey, on closer inspection it proves to be bedded. In some places, beds of relatively pure chalk alternate with beds of marl – chalk mixed with clay. In others, the alternation is between beds that are compact and densely burrowed and beds that that are not so compact or densely burrowed. Often the tops of beds are marked by nodules of flint.
The interpretation given to the chalk-marl couplets is that they reflect the precession cycle, with the chalk accumulating when the climate was warmer and the marl when the climate was colder and wetter. Clay, a product of erosion, eroded from the land more quickly in the wetter climate. Oscillations in oxygen isotope ratios (δ18O) appear to confirm the interpretation, since they follow oscillations in temperature and trace a similar pattern. Climate change thus provides an explanation for the occurrence of bedding in chalk.
If the Milankovitch theory is correct, such cyclicity can be used to test the radioisotope timescale. For example, if radioisotope dating attributes 3 million years to a sequence and the precession cycle during the period was 20,000 years, the sequence should comprise 150 such couplets. Substantially more or fewer than 150 would be evidence that the timescale was in error. In practice, the two chronological systems rarely give the same result. Sometimes even approximate agreement is achieved only on the premise that the couplets reflect obliquity rather than precession – or a mixture. Some studies claim to have detected cyclicity corresponding to all three Milankovitch cycles. In such instances the smallest cycle is visible as repeating couplets and the larger cycles as ‘bundles’ of couplets – a bundle of four or five couplets being seen as evidence of the eccentricity cycle. Thus, in relation to the Biancone Formation, in the Alps, Mayer and Appel (1999) report:
Cycle periods of 45 cm, 80 cm and 180 cm likely correspond to dominant precession, obliquity and eccentricity cycles. Owing to the inaccuracy of the Cretaceous time scale, periods cannot be matched exactly, but cycle ratios are extremely close to expected ratios so that Milankovitch climate cycles could be positively identified in this Early Cretaceous section.
Since the two chronological systems do not agree, and the timescale based on radioisotope dating is not as precise as might be supposed, it is the latter which is judged to be inaccurate. Milankovitch cycles offer higher resolution, down to 20,000 years instead of uncertainties of a million years or more, and geochronologists are now using them to ‘tune’ the radioisotope timescale. The latest timescale published (Gradstein et al. 2004) tunes radioisotope dates back to the Miocene. Discontinuous portions of the Triassic, Jurassic, Cretaceous and Palaeocene have also been redated by this method.
Astronomical tuning is a rash move, for if the only criterion for determining whether a sedi- mentary sequence reflects Milankovitch cycles is the degree of closeness to the radioisotope timescale, there is no independent control on either method. The Milankovitch interpretation is validated by producing a result approximately consistent with the radioisotope dates, and the radioisotope timescale is then adjusted to it. The reasoning remains circular.
If the danger is not clear enough from the Biancone Formation, it is well illustrated by the debate over a succession in the western Dolomites, where 470 metres encompass around 600 bedding cycles. After making allowance for variable rates of sedimentation, researchers thought they could detect signals relating to both precession and eccentricity. According to Preto and Hinnov (2003):
All of the principal periodicities related to the precession index and eccentricity, calculated for 220 Ma [middle Triassic], are present: P1 (21.9 ky); P2 (17.8 ky); E1 (400 ky), E2 (95 ky), and E3 (125 ky), along with a peak at a frequency double that of the precession, which is a predicted feature of orbitally forced insolation at the equator. Components possibly related to Earth’s obliquity at ca. 35 ky and ca. 46 ky are present as well.
(P1 and P2 designate the different precession cycles, E1-E3 the different eccentricity cycles, ‘Ma’ is millions of years ago and ‘ky’ is a thousand years.) The succession preserved ‘the oldest pristine Milankovitch signature yet observed in the geologic record’. Other workers rejected the interpretation, finding that radioisotope ages from interbedded volcanic-ash layers indicated a much shorter interval. The procedures involved in analysing cyclical bedding could generate Milankovitch-like periodicities (once calibrated to the radioisotope timescale) that were nothing of the sort: the appearance of a close fit with the astronomical cycles was illusory.
There may also be mathematical grounds for thinking so. According to Mensur Omerbashich (2006), the procedures which result in claims that bedding cycles reflect multiple kinds of astronomical cycle are highly questionable. In the case he analysed, the alleged 99% confidence level was ‘mostly meaningless in the context of total-information quality’, owing to prior editing of the data and failure to distinguish between signal and noise.
Discrepancies between the two timescales, some minor, some major, can be documented many times over (e.g. Hilgen & Langereis 1989, Karner & Muller 2000, cf Ehrlich 2007). In order to save the theory, it has even been suggested that obliquity can dominate in one locality and precession, over the same interval, in another (Prokoph et al. 2001, Krijgsman et al. 2004, Nederbragt et al. 2007). Like radioisotope dating, Milankovitch theory is an example of a scientific paradigm (Kuhn 1962), an all-embracing explanation which provides the framework within which ‘normal science’ can operate but which is not itself tested or questioned. Problems are identified and solved, so far as possible, within its terms. Anomalies with the potential to undermine the paradigm are put to one side, ascribed to errors in the data, or accommodated by modifying the paradigm, but faith in the paradigm itself remains strong. Thus, science can be both critically questioning in the domain where normal science takes place and unquestioning at the fundamental level where what is at stake is its world-view and the sense of belonging to a brotherhood of common values and beliefs.
An alternative to Milankovitch
That said, the phenomenon of cyclical bedding is real, and if it does not reflect Milankovitch cycles, it must reflect some other rhythm. To return to chalk as an example, beds are often very extensive, and in one case have have been traced from outcrop to outcrop over more than one and half thousand miles (Gale et al. 1999). Chalk cyclicity calls for climatic explanations that apply on a global scale.
What if the warm-cool cycle of each couplet was of summer and winter – the prime astronomical cycle – rather than precession, so that we were looking at the deposit of just one year? Figure 1 (from Gale et al.) illustrates a series of chalk-marl couplets from the Cenomanian part of the Upper Cretaceous. The grey-scale fluctuations track the proportion of clay in the chalk, from 5-15% in the light-coloured beds to 10-30% in the darker beds. The top of each cycle is an omission surface, where sedimentation ceased for a time and animals made burrows. Although burrowing (‘bioturbation’) is reported to be pervasive, in the chalk beds it appears to be intensive only downward from the bedding surface. Most of the marls are totally bioturbated, suggesting that they provided a more oxygenated or nutrient-rich environment for deposit-feeders. This is also evidenced by the lack of significant bioturbation from the marl into the chalk.
Dense bioturbation does not take thousands of years. Laboratory experiments have shown that a small population of deposit-feeders, such as those that made the crustacean burrows known as Thalassinoides (Fig. 1), could churn a 10 cm depth of sediment in 42 days (Gingras et al. 2008). This is the kind of timescale that should be considered when interpreting the moderate bioturbation of the upper few cm of the (compacted) chalk beds. The lower level of activity in the lower part of the beds suggests even shorter periods. The coincidence of extensive burrowing at the top with the virtual cessation of sedimentation suggests that bioturbation did not occur lower down in the bed because the rate of sedimentation was too high (and/or activity was seasonal, related to life span and life cycle). The detail of the fabric also points to short periods, as when the walls of later generations of burrows are more sharply defined than earlier ones, indicating an initially high water content followed by dewatering and compaction (Jablonski & Bottjer 1983).
Figure 2 shows the δ18O profile of three chalk-marl couplets in the ‘white cliffs’ between the English towns of Dover and Folkestone, also of Cenomanian age. The zig-zag pattern follows that of the greyness scale (Fig. 1), and since low oxygen isotope ratios in sea-floor oozes correlate with warm sea-surface temperatures, and vice versa, explaining this as the effect of summer and winter does not seem unreasonable. The same pattern may be seen in the shells of individual organisms. A study of the tropical bottom-dwelling foram Cyclobirculina, for example, which has a life span of one year, recorded a variation in δ18O ratio from -0.5 in the summer to +1.5 in the winter (Wefer & Berger 1980 and Fig. 3). Growth rate increased with age, peaking in the spring. Similar patterns have been detected in other kinds of carbonate-secreting organism, both from the past (Purton & Brasier, 1999) and from the present (e.g. Dunbar & Wellington 1981, Hickson et al. 1999).
Since coccolithophores are much smaller than forams and have a life span of only days or weeks, the ooze composed of their remains would reflect the annual oscillation collectively, not individually. They are more prevalent in the chalk than the marl and are thought to represent productivity (Leary et al 1989). The chalks are dominated by microplankton that lived close to the surface, the marls by microplankton that lived in deep water or on the sea bottom. Temperatures rose steeply in the course of marl deposition, peaked near the base of the chalk beds, then fell less steeply (Fig. 2). Fossil escape burrows have been reported from some bases (Leary et al. 1989), indicating that the rate of deposition in the lower part of the chalk beds was sometimes extremely high. This would also account for the relatively sharp transition from marl to chalk, and for the chalk being rather purer nearer the base than nearer the top (Fig. 1). The overall difference in composition between chalk and marl would be due to higher rates of planktonic productivity during early summer.
In modern settings the rate of coccolith flux to the seafloor has also been observed to vary seasonally. Around the Canary Islands sedimentation peaks in February/March following phytoplankton blooms in the winter (Sprengel et al. 2002). In the Arabian Sea the settling rate peaks in January and declines to very low levels in the months of August to November (Andruleit et al. 2004). Thus the hiatus in sedimentation at the top of the chalk beds is directly comparable to the virtual cessation of coccolith flux in these settings from late summer to late autumn.
Similar patterns characterise the chalks that formed later in the Cretaceous, during the Maastrichtian. As with outcrops onshore, cores extracted from the North Sea show metre-scale cycles, but here the alternation is between chalk that consisted of porous, mostly unburrowed sediment and chalk, higher up, that is compact and intensively burrowed (Scholle et al. 1998); in the lower part of the cycle the sediment retains its primary laminated structure. There is little bioturbation other than the vertical escape burrows of animals momentarily buried by dumps of sedimentary snow. Along with the high porosity this obviously indicates a rapid rate of sedimentation. Further up, the rate slowed, allowing burrowing organisms to destroy the laminae. Calcite cementation began to harden the seafloor, enhancing burrow preservation, but insufficient time elapsed prior to the resumption of heavy sedimentation for true hardgrounds to form.
In the Central Graben of the North Sea the chalk cycles are of non-bioturbated, mainly laminated beds and slightly thicker bioturbated beds (Damholt & Surlyk 2004). The laminated beds consist of alternating millimetre-thick, graded, high-porosity laminae and non-graded, low-porosity laminae. The graded laminae are thought to have been deposited from small-volume turbidity currents and suspension clouds, the non-graded laminae from a rain of pelleted coccoliths produced directly from the plankton. In either case, only the briefest of timescales is required – minutes in the case of the graded laminae, hours in the case of the ungraded laminae. The bioturbated beds imply longer timescales, perhaps a few months in total followed by negligible sedimentation and then the commencement of the next cycle.
Are sedimentation rates of one cycle per year credible?
According to radioisotope dating, the Cretaceous period lasted 80 million years – a stretch of time that, accustomed though we are to such figures, is unimaginable. Is it also imaginary? Since the figure is based on data extraneous to the properties of the rocks themselves, those properties can provide a test of the timescale, and as we have seen, in the case of chalk sequences, they do not agree with ages of this magnitude. To put it another way, if we were to consider the primary evidence on its own merits, we would have no grounds for inferring the miniscule rates of sedimentation the radioisotope timescale imposes.
A typical chalk-marl couplet has a thickness of around 50 cm, sometimes less. Before compaction the thickness might have been double that, around 100 cm. Averaged over 20,000 years, this gives a sedimentation rate of 0.05 mm, or half a hair’s breadth, per year. That too is difficult to imagine. By contrast, if the couplets were annual, the sedimentation rate would have averaged 2.7 mm per day – higher during productivity blooms (sufficient, as we have seen, to bury organisms alive), lower at other times.
Measurements of modern sedimentation rates tend to be expressed in milligrams per square centimetre. With solid chalk having a density of around 2,500 mg (2.5 g) per cc, the uncompacted density would be about 1,250 mg per cc. Thus, if chalk-marl couplets represented 20,000 years, the sedimentation rate would be 6.25 mg per year or 0.017 mg per day; if the couplets represented one year, the rate would be 125 g per year, or 340 mg per day. Modern rates vary according to the same factors that ancient rates vary, but are closer to 6.25 mg than to 125 g per year. For example, Broerse et al. (2000) reported a carbonate flux of 2.2 mg per cc per year, which is even lower than that calculated for the Cenomanian.
However, modern rates do not provide a meaningful comparison. The escape burrows and the generally low level of bioturbation suggest that the chalk formed rapidly, mostly during plankton blooms that were much more frequent than in modern times. In normal modern conditions the export of carbonate to the seafloor is limited by a negative feedback loop, where 
low productivity is accompanied by low settling velocity, which results in high rates of consumption by animal grazers and high rates of chemical dissolution. In bloom conditions, by contrast, not only does the rate of production increase by orders of magnitude but faecal-pellet clumping leads to greatly increased settling velocity, with consequently low rates of consumption by grazers and low rates of dissolution (Wal et al. 1995). Coccolithophores, it should be noted, are ‘the most productive calcifying organisms on earth’ and can reach a density during blooms of some 115 million per litre. The great chalk formations bear witness to the exceptional ‘hypercalcifying activity of the coccolithophorids’ (Aubry 2005).
The main factors governing the dissolution of calcium carbonate are temperature, acidity and pressure. Calcium carbonate dissolves more readily where the water is colder, more acidic and deeper. Acidity decreases as temperature increases because warm water holds less carbon dioxide and thus less carbonic acid. In the warm, shallow seas of the Cretaceous, where the chalks were laid down, the conditions for calcium carbonate precipitation may have been optimal. Not only were the seas warm because they were shallow (most of Britain, for example, was inundated at this time), but the global climate itself was exceptionally warm (Huber et al. 2002).
Another point about the Cretaceous is that the ratio of magnesium to calcium in seawater was much lower than in modern seas. High levels of Mg have been found to inhibit growth rates. In experiments involving three species of coccolithophore in low-Mg seawater all showed dramatic increases in population growth rates; two of the species reproduced at rates approximately triple those observed in modern seawater (Stanley et al. 2005). The species dominating Cenomanian chalks are unlikely to have been less responsive. The ubiquitous Watznaueria barnesae, for example, is thought to have been highly prolific (Lees et al. 2005). Although the data are less clear, the same seems to have been true of Biscutum constans (or ellipticum).
Thus several lines of evidence additional to the purely sedimentary evidence support the inference that the Chalk was deposited rapidly compared to modern rates of deposition:
- The Chalk consists almost entirely of the remains of coccolithophores, which in the right conditions can reproduce at a prodigious rate.
- The warmth and shallowness of the waters where the coccoliths accumulated meant that a much smaller percentage of platelets shed by the coccolithophores dissolved before they reached the bottom.
- The low Mg/Ca ratio of Cretaceous seawater favoured coccolithophore reproduction.
- Some of the coccolithophore species that dominated chalk assemblages are known to have been prolific examples of their kind.
The role of faster plate tectonics
Coccolithophore blooms can occur only if there is a sufficient supply of inorganic nutrients, calcium ions and dissolved carbonate. Once the supply is exhausted the blooms cease, and they will not have the potential to recur until the minerals are replenished. Thus in the context of sedimentation rates of around 1 metre per year the final question to consider is how the minerals could be replenished that quickly.
Whereas modern seas are relatively poor in calcium, brines from the Cretaceous have high Ca/Na ratios and low Mg/Ca ratios, indicating that seawater was relatively rich in calcium. Seawater chemistry is controlled primarily by the flux of hydrothermal brine from mid-ocean ridges, and Cretaceous seas were rich in some minerals because rates of ocean crust production were higher then (Hardie 1996).
On the assumption of unchanging rates of radioactive decay, the rate of ocean crust production in the Middle Cretaceous is estimated to have been almost twice the modern rate, declining progressively thereafter to present-day values. But the seafloor spreading rates remain tiny, of the order of centimetres per year. On the other hand, if radioactive decay was higher in the past, enhanced radioactivity would have had the effect of speeding up geological processes, including activity at the mid-ocean ridges. The flux of hydrothermal brine into the ocean would then also have been higher. Large quantities of mineral nutrients would have been released directly into global current systems (which were vigorous), and large amounts of carbon dioxide would have been transferred from the mid ocean crust to the shallow epicontinental seas indirectly via atmospheric circulation. As is well known, the oceans are a major sink of carbon dioxide, on a par with tropical forests. The minerals removed from the oceans by coccolithophore blooms would have been quickly replenished in these circumstances. Along with levels of atmospheric carbon dioxide ten times higher than now (Royer et al 2004), the whole phenomenon of chalk may therefore be seen as indirect evidence of elevated levels of plate-tectonic activity.
Conclusions
Absolute dates in geology are inferences based on radioactive decay, where the key assumption is that the rate of decay has been constant over time. These dates command universal acceptance because they accord with a world-view in which long ages are required. The suggestion that they ought to be subjected to critical, independent testing is so heretical that it is, one suspects, inwardly suppressed before it can be entertained, let alone discussed in the literature.
Although Milankovitch cycles have the potential to test radioisotope dating, they come with their own presuppositions concerning the age of the solar system, the imperturbability of the Earth-Moon system over that time, the linearity of secular change in astronomical cycles, and the likelihood that these should dominate sedimentary processes. They are part of a world-view which assumes the constancy of isotopic half-lives over time, and in practice they are used only to ‘tune’ the radioisotope timescale, not test it. Tuning is necessary, for ‘very few rhythmically bedded geological sections yield clear-cut power spectra when calculated against stratigraphic depth, whatever the age of the section’ (Nederbragt et al. 2007). Discrepancies between the two chronologies are minimised in a variety of ways, including (1) interpreting some cycles as relating to precession, others to obliquity, (2) postulating that some cycles are not geologically preserved, and (3) proposing adjustments to the radioisotope timescale that leave the basis of radioisotope dating unchallenged.
As well as being self-consistent, radioisotope dating needs to be consistent with the nature of what is being dated. Our contention is that it is fundamentally inconsistent. The Earth History website sets out several instances of conflict between it and direct indications of process timescales (for links, see the Key Concepts section). Here we have focused on the possibility that the warm-cool cycles that characterise relatively tranquil offshore settings in the Cretaceous were annual cycles, and therefore represent sedimentation rates orders of magnitude higher than those inferred from the radiometric timescale. Perhaps the most obvious evidence is the pattern of bioturbation through the chalk beds. Given sedimentation rates of 0.05 mm per year, mobile deposit feeders should have repeatedly churned up their entire depth; there should never have been a time when any significant depth was not bioturbated. The very limited amount of bioturbation suggests that the cyclicity was annual.
Assuming that cycles attributed to precession (the shortest Milankovitch period) were annual, a simplistic calculation would divide the 80 million radioisotope years of the Cretaceous period by 20,000 to arrive at a true duration of around 4,000 years. However, Milankovitch cycles do not necessarily provide a good fit for sections that appear to be well dated, and the appearance of chronological precision can be illusory – to give an example, estimates of the duration of Ocean Anoxic Event II (latest Cenomanian) vary from 320,000 years to fully three times that figure. Short intervals are chronologically more difficult to quantify than long intervals.
Based on the growth rates of rudist bivalves, carbonate reef successions suggest that the length of a warm-cool cycle, if it records the annual cycle, is several times less than 20,000 radioisotope years and the length of the whole Cretaceous correspondingly greater than 4,000 true years. If the disparity between radioisotope and true time decreases as one approaches the historical period, the rhythmicity that continues to be visible in sedimentary successions far into the Cenozoic should prove to be even less amenable to Milankovitch explanations.