AMHERST COLLEGE Department of Geology GEOLOGY 41 - Environmental and Solid Earth Geophysics Lab 2: Geochronology Solution U/Pb data aliquot # U (ppm) Pb (ppm) 206 Pb/ 204 Pb 207 Pb/ 204 Pb 208 Pb/ 204 Pb z1 228.4 23.21 1341.89 93.260 417.990 z2 436.2 57.03 544.604 46.980 269.640 z3 349.1 40.81 468.614 42.634 118.620 z4 542.7 58.82 573.22 47.697 275.262 1) The concentrations of U and Pb are given in mass units (ppm) and need to be converted to the number of atoms of U and Pb in each aliquot. To so this, divide the concentration of U and Pb of each aliquot by the atomic weights of U and Pb. (This method gives you the number of atoms in an arbitrary mass of zircon, but it works because the age equations deal with atom ratios). For this first aliquot this is: U = 228.4/238.0291 = 0.9595 Pb = 23.21/206.18 = 0.1126 for all of our aliquots we get: U (atoms) Pb (atoms) z1 0.9595 0.1126 z2 1.8326 0.2766 z3 1.4666 0.1979 z4 2.2800 0.2853 2) For each aliquot determine the number of atoms of each of the U ( 235 U, 238 U) and P ( 208 Pb, 207 Pb, 206 Pb, 204 Pb) isotopes. For U, this is relatively easy as the ratio of U isotopes in all U is constant. For Pb, this is a bit more complicated because 3 of the Pb isotopes are radiogenic and the abundance of each isotope must be calculated from the analytical data each aliquot. (Do you know how to do this?-if not see Peter) U = 238 U + 235 U U/ 235 U = 238 U/ 235 U + 1 = 1 + 1/0.00725670 =138.804 235 U/U = 1/138.804 = 0.00720442 238 U/U = 235 U/U x 238 U/ 235 U = 0.00720442 / 0.00725670 = 0.99279 235 U = U (atoms) x 235 U/U = U (atoms) x 0.0072044 238 U = U (atoms) x 238 U/U = U (atoms) x 0.99279
We can get the Pb ratios the same way, but we need to do this for each aliquot. Pb = 204 Pb + 206 Pb + 207 Pb + 207 Pb Pb/ 204 Pb =1 + 206 Pb/ 204 Pb + 207 Pb/ 204 Pb + 207 Pb/ 204 Pb 204 Pb/Pb =1/ (1 + 206 Pb/ 204 Pb + 207 Pb/ 204 Pb + 207 Pb/ 204 Pb) and like U 204 Pb = U (atoms) x 204 Pb/Pb this gives us 204 Pb, we can get the others by: 206 Pb/Pb = 204 Pb/Pb x 206 Pb/ 204 Pb 207 Pb/Pb = 204 Pb/Pb x 207 Pb/ 204 Pb 208 Pb/Pb = 204 Pb/Pb x 208 Pb/ 204 Pb 204 Pb for all of our aliquots this is: 206 Pb 207 Pb 208 Pb z1 0.0000607 0.0815 0.00566 0.0254 0.00691 0.0952 z2 0.0003208 0.175 0.0151 0.0865 0.0132 1.819 z3 0.0003137 0.147 0.0134 0.0372 0.0106 1.456 z4 0.0003180 0.182 0.0152 0.0875 0.0164 2.263 3) Some of the Pb in your zircon was present when it crystallized. You know that this is the case because there is some non-radiogenic 204 Pb in your samples. Use the common Pb isotopic ratios given below to determine the amount of radiogenic 207 Pb & 206 Pb in each aliquot (i.e., determine the amount of 207 Pb* & 206 Pb* in each aliquot). 235 U..OK this is getting fussy, but is really just algebra. 206 Pb* = 206 Pb - 206 Pb common = 206 Pb - 204 Pb x ( 206 Pb/ 204 Pb) common 206 Pb* = 206 Pb - 204 Pb x 17.915 238 U for all of our aliquots this is: 206 Pb* 207 Pb* z1 0.0804 0.00472 z2 0.0169 0.0101 z3 0.141 0.00849 z4 0.0177 0.0102 4) For each aliquot determine the 207Pb*/235U & 206Pb*/238U ratios and using the decay equations for 235U & 238U determine a 207Pb*/235U and 206Pb*/238U age. Are these ages concordant? Are the ages for each aliquot the same? This is just plug and chug using the U-Pb age equations. t 238 = 1 ln ( 206 Pb " + 1 )! 238 238 U or for aliquot z1:
t 238 = 1 ln ( 80.383!! 1 0 1.55125 x 10 952.634 + 1 )! 522. 2 x 10 6 years for all of our aliquots this is: 206 Pb* /238 U 207 Pb* /235 U 206 Pb /238 U 207 Pb /235 U age (Ma) age (Ma) z1 0.0844 0.682 522.2 528.1 z2 0.0929 0.763 572.5 575.7 z3 0.0971 0.803 597.5 598.7 z4 0.0780 0.622 484.3 490.9 The ages are not concordant and vary by 100 Ma from aliquot to aliquot
5) Plot the 207 Pb*/ 235 U & 206 Pb*/ 238 U ratios for each aliquot on a concordia diagram. Can you now make sense out of your zircon ages? 0.100 600 0.095 2 0 6 P b / 2 3 8 U 0.090 0.085 550 0.080 500 0.075 0.60 0.65 0.70 0.75 0.80 0.85 0.90 2 0 7 P b / 2 3 5 U Aliquots fall along a discordia array with an ~600 Ma upper intercept indicating a 600 Ma age for the igneous? Protolith of the gneiss. 6) What observations can you make about the tectonic evolution of the Pelham Hills based on your U/Pb zircon and 40 Ar/ ages? T closure (zircon) > 600 C 600 Ma T closure (hornblende) = 550 C 360 Ma T closure (muscovite) = 400 C 335 Ma Plutonism at ca. 600 Ma, followed by cooling to low? temperature. Metamorphism (and tectonism) at or before ca. 360 Ma followed by slow cooling (6 deg.ma) through 335 Ma (and younger times?). Isostatically driven unroofing? 40 Ar/ data T( C) of heating % / 40 Ar released 36 Ar/ 40 Ar
step this step 56.7 Ma biotite standard fusion 100 0.4122 0.00001 hornblende muscovite 600 12 0.0408 0.00041 700 17 0.0514 0.00046 800 28 0.0557 0.00027 900 18 0.0586 0.00007 1000 15 0.0573 0.00009 fusion 10 0.0531 0.00033 500 9 0.064 0.00002 600 23 0.0651 0.00001 700 14 0.0643 0.00003 800 28 0.0644 0.00003 fusion 26 0.0629 0.00017 1) In order to determine the 40 Ar/ ages of your samples, you must know the ratio 40 Ar*/, the ratio of radiogenic 40 Ar to. To do this, you must correct your measured ratios of 40 Ar/ for any Ar the was introduced during sample preparation. You know that some Ar must have been introduced into each aliquot that was analyzed because there is some 36 Ar in each aliquot. Assume that the source of 36 Ar is modern atmosphere ( 40 Ar/ 36 Ar = 295.5) and correct each of the measured values of 40 Ar/ for 40 Ar that was introduced during sample preparation and determine 40 Ar*/ for each aliquot. To correct for atmospheric 40 Ar we need to subtract it from the total 40 Ar. (a) 40 Ar! = 40 Ar " 40 Ar atm = 40 Ar " 295. 5 x 36 Ar We only have Ar isotope ratios, but that is OK, because we need 40 Ar*/, so lets divide (a) by to get: 40 Ar! = 40 Ar " 295. 5 x 36 Ar = 40 Ar " 295. 5 x 40 Ar 36 Ar 40 Ar
...for this biotite standard this is: 40 Ar! = 1 0. 4122 " 295. 5 x. 00001 x 1 0. 4122 = 2. 419 For all of our aliquots we get: standard fusion 2.419 hornblende muscovite 600 21.540 700 16.811 800 16.521 900 16.712 1000 16.988 fusion 16.996 total gas 17.324 500 15.533 600 15.316 700 15.414 800 15.390 fusion 15.100 total gas 15.314 2) Determine the neutron flux parameter J for your analyses using the Ar isotopic ratios collected from the 56.7 Ma biotite standard. This comes directly from the equation for J, the age of the standard, and the 40Ar*/39Ar ratio calculated above. 40 Ar! = 1 J ( e " t # 1 ) J = 1 / 40 Ar! or using the standard data we get J = 1 2. 941 ( e 5. 453 x 57. 6 x 10! 4! 1 ) = 0. 013199 ( e " t # 1 ) 3) Determine an 40 Ar/ age for each of the heating steps for the two samples (muscovite and hornblende). Also determine a total gas 40 Ar/ age for both samples. We can solve our equation for t and calculate ages for each aliqot. 40 Ar! = 1 J ( e " t 1 # 1 ) t =! ln ( J x 40 Ar " + 1 ) or for the 600 deg hornblende aliquot this is:
t = 1 5. 545 x 10! 1 0 ln ( 0. 013199 x 21. 45 + 1 ) = 451. 4 x10 6 years
For all of our aliquots we get: standard fusion 56.70 hornblende 600 451.44 700 361.54 800 355.89 900 359.62 1000 364.99 fusion 365.15 total gas 371.5 muscovite 500 336.46 600 332.17 700 334.12 800 333.65 fusion 327.88 total gas 332.13 4) Using the ages determined above, plot an Ar age spectrum (apparent age vs. cumulative percent of Ar) for each of the 2 samples.
400.00 390.00 380.00 370.00 Age (Ma) 360.00 350.00 340.00 hbl musc 330.00 320.00 310.00 300.00 100 80 60 40 20 0 cumulative Ar
5) Assuming a ±3 Ma error in the age of each release step, do your ages represent a plateau? YES