Error propagation

All the methods and equations presented thus far have assumed that all
parameters are either known or measured with infinite precision. In
reality, however, the analytical equipment used to measure isotopic
compositions, elemental concentrations and radioactive half-lives is
not perfect. It is crucially important that we quantify the resulting
analytical uncertainty before we can reliably interpret the resulting
ages.

For example, suppose that the extinction of the dinosaurs has been dated at 65 Ma in one field location, and a meteorite impact has been dated at 64 Ma elsewhere. These two numbers are effectively meaningless in the absence of an estimate of precision. Taken at face value, the dates imply that the meteorite impact took place 1 million years after the mass extinction, which rules out a causal relationship between the two events. If, however, the analytical uncertainty is significantly greater than 1 Myr (e.g. 64 ± 2 Ma and 65 ± 2 Ma), then such of a causal relationship remains very plausible.

Suppose that our geochronological age (t) is calculated as a function (f) of some measurements (X and Y ):

| (10.1) |

Suppose that we have performed a large number (n) of replicate measurements of X and Y :

| (10.2) |

It is useful to define the following summary statistics:

- The mean:
(10.3) is a useful definition for the ‘most representative’ value of X and Y , which can be plugged into Equation 10.1 to calculate the ‘most representative’ age.

- The variance:
(10.4) with σ

_{X}and σ_{Y }the ‘standard deviations’, is used to quantify the amount of dispersion around the mean. - The covariance:
(10.5) quantifies the degree of correlation between variables X and Y .

X, Y , σ_{X}^{2}, σ_{Y }^{2} and cov(X,Y ) can all be estimated from the
input data (X,Y ). These values can then be used to infer σ_{t}^{2}, the
variance of the calculated age t, a process that is known as ‘error
propagation’. To this end, recall the definition of the variance (Equation
10.4):

| (10.6) |

We can estimate (t_{i} -t) by differentiating Equation 10.1:

| (10.7) |

Plugging Equation 10.7 into 10.6, we obtain:

σ_{t}^{2} | = ∑
_{i=1}^{n}^{2} | (10.8) |

= σ_{X}^{2}^{2} + σ_{
Y }^{2}^{2} + 2 cov(X,Y ) | (10.9) |

This is the general equation for the propagation of uncertainty with two variables, which is most easily extended to more than two variables by reformulating Equation 10.9 into a matrix form:

| (10.10) |

where the innermost matrix is known as the variance-covariance matrix and the outermost matrix (and its transpose) as the Jacobian matrix. Let us now apply this equation to some simple functions.

Let X and Y indicate measured quantities associated with analytical uncertainty. And let a and b be some error free parameters.

- addition:
t = aX + bY ⇒ = a, = b ⇒ σ _{t}^{2}= a^{2}σ_{ X}^{2}+ b^{2}σ_{ Y }^{2}+ 2ab cov(X,Y )(10.11) - subtraction:
(10.12) - multiplication:
t = aXY ⇒ = aY, = aX ⇒ σ _{t}^{2}= (aY )^{2}σ_{ X}^{2}+ (aX)^{2}σ_{ Y }^{2}+ 2a^{2}XY cov(X,Y )⇒ ^{2}=^{2}+^{2}+ 2(10.13) - division:
(10.14) - exponentiation:
(10.15) - logarithms:
(10.16) - power:
(10.17)

Recall the definition of the arithmetic mean (Equation 10.3):

Applying the equation for the error propagation of a sum (Equation 10.11):

| (10.18) |

where we assume that all n measurements were done independently, so
that cov(X_{i},X_{j}) = 0∀i,j. The standard deviation of the mean is known as
the standard error:

| (10.19) |

This means that the standard error of the mean monotonically decreases
with the square root of sample size. In other words, we can arbitrarily
increase the precision of our analytical data by acquiring more data.
However, it is important to note that the same is generally not the case for
the accuracy of those data. The difference between precision and accuracy is
best explained by a darts board analogy:

Whereas the analytical precision can be computed from the data using the error propagation formulas introduced above, the only way to get a grip on the accuracy is by analysing another sample of independently determined age. Such test samples are also known as ‘secondary standards’.