In this post I’d like to outline some of the more important aspects of Ar/Ar thermochronology. I am going to start with a brief outline of the K-Ar technique, but spend most of the post on the more popular Ar/Ar variant. This post is a long time coming; as far as published reliable radiometric dates, I’d bet that Ar/Ar thermochronology is the most used chronometer in earth science (although I did recently hear Mark Harrison claim he had “age characterized” 100,000 zircons for U-Pb geochronology, so I could be way way off). There are a few primary reasons for the popularity of the method, some I’ll discuss later. But first the basics.
Potassium (K) is the 7th most abundant element in the earth’s crust, and is a major component of many common rock forming minerals. Very common minerals like biotite and muscovite mica and potassium feldspar contain weight percent levels of K. Naturally occurring K is made of 3 isotopes, 39K (93.258%), 40K (0.012%) and 41K (6.730%). 39K and 41K are stable, 40K is radioactive. Natural Ar is fairly abundant in the earth’s atmosphere (~1% of the air we breathe) and is composed of three isotops, 40Ar (99.600%), 38Ar (0.063%), and 36Ar (0.337%). 40K decay is an interesting process, unlike the U and Th systems I’ve blogged about earlier which undergo a chain decay, 40K has a branched or dual decay scheme. That means that 89.52% of the time, 40K decays by beta emission to 40Ca. 10.48% of time time, 40K decays by extranuclear electron capture to 40Ar. The 40K to 40Ar process is the one we will concern ourselves with in this post. 40K – 40Ca is rarely used as a geochronometer because Ca is so abundant in nature. It is hard to differentiate in minerals what Ca formed in the mineral from decay of K and what Ca was included in the original mineral structure. Argon, a noble gas, is rarely present in significant amounts within crystals, so we do not have to worry as much about the original 40Ar component (although you will see later that the 40Ar/39Ar technique lets us deal even with “inherited” 40Ar). 40K has a half-life of 1.25 billion years (Ga), which makes it perfect for addressing a wide range of geologic problems. The graph below shows the evolution of 40K, 40Ca, and 40Ar abundance versus time.
Like in all graphs that show parents and daughters, the curves intersect at the half-life. Focusing on the 40Ar then, we can see that the daughter to parent ratio, or 40Ar/40K changes over time, like this
We can measure the 40Ar and 40K in a sample, and from that ratio, calculate an age. The age equation (shown below)
then relates the amount of parent (40K), daughter (40Ar) and time. Like with all geochronologic systems, actually calculating the age is fairly straightforward. The actual K-Ar age equation is a little more involved because of the branched decay I mentioned earlier, so the final form looks like thiswhere
is the 40K decay constant (total)
is the decay constant for 40Ca production from 40K. The term in the age equation just before the daughter to parent ratio therefore describes the fraction of 40K decays that result in the formation of a 40Ar atom.
The K-Ar geochronometer came into widespread use in the 1940’s, once the basics of the decay scheme were worked out. There are two main issues with the technique, however, that are kind of a pain, and are the reasons that more often that not the 40Ar/39Ar technique is much more popular today. First, you cannot measure K and Ar on the same machine. Typically, samples would first be heated (often using an induction furnace) to very high temperatures (~1400°C). At these temperatures the Ar trapped in the sample is liberated and can then be measured. The 40Ar from the sample would be spiked with 38Ar, and then measured on a mass spectrometer. The degassed sample would then be retrieved from the furnace, and analyzed separately for 40K. With the amount of 40Ar and 40K, you culd then calculated an age.
The second problem with the K/Ar technique is determining where the Ar came from. As I mentioned earlier, Ar is fairly abundant in the atmosphere, making up about 1% of what we breathe. It is also fairly common in groundwater and geo-fluids (fluids associated with igneous intrusions and metamorphism). Some of this Ar could become trapped in the crystals, therefore rendering the calculated age kind of meaningless. We know the atmospheric 40Ar/36Ar ratio very well, so by measuring the 36Ar we can correct for this “trapped Ar”, but we have to assume that the composition is atmospheric. In some cases this is a perfectly reasonable assumption, but with the advent of the 40Ar/39Ar technique, we are actually able to evaluate that assumption.
The basis of the 40Ar/39Ar technique is the fact that when K is bombarded with neutrons (in a nuclear reactor), some important reactions take place. Namely, 39K is converted into 39Ar. We know the natural 39K/40K ratio, so if we can measure the efficiency of this transformation in the reactor, we can use 39Ar as a proxy for 40K. This would allow us to measure an age by simply measuring the 40Ar/39Ar ratio of a sample. What is even better is that 39Ar does not exist naturally, so we do not have to worry about any inheritance or atmospheric contamination.
39Ar production can be described by this equation
is the neutron flux at energy Eis the neutron capture cross section at energy E for the reaction 39K(n,p)39Ar, and
is the duration of the irradiation. In practice, instead of measuring these values directly we figure them out using well defined age standards, producing this equation
where 40Ar* is the 40Ar derived from radioactive decay, 39Ar K is the 39Ar produced in the reactor from 39K, and J is a term that incorporates all of the irradiation variables but can be calculated from standards of known age (determined with other methods), whereThe advent and refinement of the 40Ar/39Ar method solved some of the larger problems with the K/Ar technique. First, sample analysis became much more straightforward, measuring different isotopes of the same element is much easier that measuring different elements. Measuring a ratio of two isotopes is also easier that measuring abundances of two elements. The 40Ar/39Ar frees us from uncertainty in the spikes used during analysis.
Second, it allows the geochronologist to evaluate the assumption that the “extraneous Ar” in the sample is atmospheric in composition. This is done by progressively degassing the sample at higher and higher temperatures (called step-heating experiments). The first temperature steps are at low temperatures, and tend only to liberate very loosely bound Ar from the crystal, perhaps Ar adsorbed to the surface of the crystal, or held near the surface. With increasing temperature we access and analyze gas that is more tightly bound up in the crystal structure. In the end we then have a series of steps, each with separate isotopic values. If the gas came only from the atmosphere, and was not related to any K, then 39Ar/40Ar would equal 0. If the gas was entirely radiogenic, and had no atmospheric input, then 36Ar/40Ar would equal 0. In truth, we typically have some mixture of inherited (or atmospheric) Ar and radiogenic Ar. If we plot up these values on an inverse isochron diagram we can determine what the composition of the extraneous Ar reservoir is, and therefore correct our values to determine the real age.
as shown above. If your sample is a mixture of inherited Ar from a single reservoir and radiogenic Ar of a single age, then your individual temperature steps will define a mixing line between those two compositions. We can use the inherited Ar composition to correct our 40Ar values, and get the real age. It is also common to see the data displayed as an age spectra, shown below. Here, each box represents an indivual temperature step from a single sample. The y-axis is the age (here in Ma), corrected for extraneous Ar. The height of the box is proportional to the uncertainty of the measurements. The x-axis is the cumulative percent of Ar released. The steps are in order, so the first low temperature steps are on the left, and progressively higher temperatures towards the right. The width of the box is proportional to the amount of radiogenic Ar (40Ar*) released during that step. The large flat area in the middle of the graph is called a plateau and indicates that ages from different temperature steps agree. There is no real physical basis for the definitions of a plateau, and variations abound. Realistically, if our measurement precision was high enough, we would never see plateaus, because natural volume diffusion should round out the profile and make the early temperature steps slightly younger. But, age spectra are very common in literature, so it is important to know how to read them.
The inverse isochron diagram is much more powerful, and I am suspect of Ar papers that do not include inverse isochron diagrams. If the sample is very radiogenic, and there is little 36Ar, then the corrections aren’t significant, but still that information should be included. Plateaus and age spectra diagrams are interesting on their own, almost from a Psychology of Science perspective. Even though they are not the most powerful Ar/Ar tool, they are visually pleasing, and telling someone you got a nice plateau usually convinces them of your data quality faster than discussing inverse isochron diagrams. The definitions of a proper plateau usually involve some minimum number of consecutive steps representing some fraction of the total gas that are within error of each other for age and have no statistically discernible slope. You can imagine that if your measurements were kind of crappy, and the age uncertainty on each step high, then getting consecutive steps to agree isn’t very difficult. As labs drive for higher and higher precision, and individual step ages have uncertainties less than 1% and official plateau get tougher and tougher to create, even when the age and sample are perfectly acceptable.
40Ar/39Ar thermochronology is a very common and useful tool for geologists, and is particulary interesting because it allows the direct examination of a number of the assumptions involved with radiogenic isotope geochronology. Later posts will examine the power of the technique in more detail.
Note – the inverse isochron diagram I used is one I made up just as a teaching tool, the age spectra is real data I collected a few years ago.