1 Introduction
The inner planets of the solar system are considered to have formed from a swarm of planetesimals; meteorite parent bodies are witnesses of these early planetesimals. One of the major goals in studying meteorites is to define the accretion time scale and thermal history of their parent bodies and thus to elucidate the processes of planet formation.
Acapulcoites belong to the group of primitive achondrites and are believed to originate jointly with the closely related lodranites from a common parent body (McCoy et al., 1993). These meteorites experienced high temperatures and are highly metamorphosed and/or were partially molten. The temperature peak, ∼1150 °C, is significantly higher than that estimated for the most highly equilibrated chondrites but lower than that estimated for differentiated meteorites. Acapulco is the only fall of this class. Its bulk chemical composition is similar to that of H chondrites, but it lacks chondrules and exhibits a coarsely granoblastic texture (Christophe Michel-Levy et al., 1978, 1979; Palme et al., 1981). The meteorite is unshocked and does not show evidence for disturbance after crystallization (McCoy et al., 1996). Although Acapulco experienced partial melting, it exhibits high abundances of volatiles and trapped noble gases (Palme et al., 1981). Indications for non-equilibrium, e.g. different nitrogen isotope compositions between metal and silicates and LREE- and U- excesses in phosphates with respect to silicates have been described (Kim and Marti, 1992, 1993; Zipfel et al., 1995).
The thermal history of the Acapulco parent body is heavily debated in the literature. Its cooling rate has been studied with a variety of mineralogical tools: Ca zoning in olivine (Köhler and Brey, 1990; Zipfel et al., 1995), Fe/Mg ordering in orthopyroxene (Zema et al., 1996), metallographic observations and Ni content in taenite (McCoy et al., 1996). Based on these investigations, the cooling of Acapulco was rapid, at least 103 °C/Ma in the temperature range between 650–350 °C (Min et al., 2003) and this fast cooling has been interpreted as being the consequence of the break-up of the Acapulco parent body (Zema et al., 1996).
In contrast, no consensus exists concerning the cooling rate below 350 °C. Low cooling rates, ∼1–20 °C/Ma, have been deduced from 244Pu fission track density and metallographic studies (Marti et al., 1994; Pellas and Fieni, 1992; Pellas et al., 1997). Also long-lived chronometers, U/Pb, K/Ar, and Sm/Nd, have been applied in order to constrain the thermal evolution (Amelin, 2005; Göpel et al., 1992; McCoy et al., 1996; Pellas et al., 1997; Prinzhofer et al., 1992; Renne, 2000). The 39Ar/40Ar age, 4514 ± 16 Ma, calculated with the 40K decay constant recommended by Steiger and Jäger (1977), is ∼50 Ma younger than the 207Pb/206Pb age of the phosphate (Pellas et al., 1997). This age difference defines a slow cooling rate corresponding to ∼5 °C/Ma for the temperature interval 550–300 °C. However, Renne (2000) argued that the 39Ar/40Ar age becomes 4554 Ma, ∼50 Ma older when the 40K decay constant from Endt and van der Leun (1973) is used. This older age translates into a faster cooling rate and better agrees with the petrologic observations that support a rapid cooling. Trieloff et al. (2001) contested this approach, whereas Renne (2001) in response reiterated his arguments. Recently, new Hf/W data were published for two acapulcoites (Touboul et al., 2009). Touboul et al. (2009) estimate the Hf/W closure temperature at ∼975 °C for cooling rates ∼102 °C/Ma. These authors translate the relative Hf/W age into an absolute age, 4563.1 ± 0.8 Ma, using the U/Pb and Hf/W anchor of angrites. This age corresponds to the end of Hf/W reequilibration between metal and pyroxene.
The goal of this study is to evaluate the chemical history of the Acapulco meteorite with U/Th/Pb systematics and to better constrain the thermal history of its parent body at high temperatures, between ∼1000–350 °C. The 207Pb/206Pb dating method is of particular interest since it permits age determinations with an uncertainty of ∼1 Ma, allowing us to anchor this cooling history into an absolute time scale. This study is part of a collaborative work organized by K. Marti.
2 Sample preparation and analytical procedure
Two chips of Acapulco (1.67 and 0.89 g) were broken in a boron carbide mortar, metal was eliminated and the crushed fraction was sieved. This procedure was repeated until most of the sieved material was in the 100–200 and 200–400 mesh grain size ranges. The ground material was composed of a magnetic fraction (29.3 wt %), 100–200 mesh fraction (49.3 wt %), 200–400 mesh fraction (21.5 wt %) and a < 400 mesh fraction (0.3 wt %). Two aliquots from the 100–200 and one from the 200–400 mesh fraction were reserved for “bulk” U/Pb analysis. The rest of the material was processed on a Frantz magnetic separator and mineral separates were prepared from this aliquot.
The feldspar separate was washed prior to dissolution with acetone and leached with 0.5 M HBr. We used a mixed 205Pb/233U tracer for all samples, the U/Pb chemistry that was applied is described in Göpel et al. (1985) and Manhès et al. (1984). Total chemistry blanks correspond to 6.2 × 10−12 g Pb and 0.5 × 10−12 g U. Isotope measurements were performed on a Thomson CSF thermal ionization mass spectrometer (Birck and Allègre, 1973).
The Th content is not directly determined by isotopic dilution, it is derived from the measured Pb isotope composition, the measured U content and the 207Pb/206Pb age, assuming a Pb evolution in a closed system:
(1) |
λ and λ” are the decay constants of 238U and 232Th, respectively, “meas” corresponds to the measured Pb isotopic ratio, “ini” to the initial composition of Acapulco. This initial isotopic composition is defined in this study and will be discussed in more detail in Section 3.2.
3 Results
U and Pb concentrations, Pb isotopic compositions as well as 207Pb/206Pb model ages, 238U/204Pb and 232Th/238U ratios of all analyzed samples are reported in Table 1.
Systématique U–Pb des échantillons d’Acapulco.
Sample | U | Pb | 238U/204Pb | 232Th/238U | 206Pb/204Pb | 204Pb/206Pb | 207Pb/206Pb | 208Pb/206Pb | Pb/Pb age | ρ 4/6–7/6 | ρ 4/6–8/6 |
Weight | (ppb) + | (ppb) + | +,† | *,† | $ | + | + | + | (Ma) # | ||
feldspar | 0.05 | 2407.1 | 0.0009 | 9.346 | 0.107070 | 1.10839 | 3.1536 | 6030 | −0.378 | −0.521 | |
6.2 mg | ± 0.0007 | ± 0.006 | ± 0.000079 | ± 0.00051 | ± 0.0023 | ± 120 | |||||
fr 100 a | 82.3 | 440.4 | 12.761 | 1.03 | 24.759 | 0.040362 | 0.80431 | 1.3474 | 4562.0 | −0.506 | −0.263 |
15.1 mg | ± 0.218 | ± 0.019 | ± 0.000033 | ± 0.00029 | ± 0.0011 | ± 1.5 | |||||
fr 100 b | 63.8 | 466.4 | 9.19 | 1.55 | 23.125 | 0.043231 | 0.81767 | 1.5029 | 4563.8 | −0.535 | −0.613 |
11.5 mg | ± 0.07 | ± 0.018 | ± 0.000034 | ± 0.00030 | ± 0.0010 | ± 1.6 | |||||
fr 200 | 157.2 | 786.5 | 13.32 | 0.84 | 24.092 | 0.041496 | 0.80934 | 1.3507 | 4561.8 | −0.079 | −0.171 |
7.5 mg | ± 0.15 | ± 0.025 | ± 0.000043 | ± 0.00045 | ± 0.0012 | ± 1.9 | |||||
phosphate 1 | 2760.9 | 5106.2 | 137.11 | 0.56 | 148.02 | 0.006751 | 0.65106 | 0.32864 | 4556.9 | 0.200 | 0.068 |
6.1 mg | ± 1.13 | ± 0.18 | ± 0.000008 | ± 0.00067 | ± 0.00037 | ± 1.7 | |||||
phosphate 2 | 2317.2 | 4362.2 | 127.38 | 0.58 | 138.63 | 0.007208 | 0.6525 | 0.3457 | 4555.3 | 0.581 | 0.746 |
6.1 mg | ± 0.72 | ± 0.41 | ± 0.000022 | ± 0.0020 | ± 0.0010 | ± 5.0 | |||||
troilite | 3.7 | 108.9 | 1.72 | ||||||||
5.6 mg | |||||||||||
leachate feldspar | 0.03 | 2.9 | 0.44 | 1.35 | 13.779 | 0.0794 | 0.979 | 2.426 | 4556 | 0.991 | 0.994 |
6.2 mg | ± 0.70 | ± 0.051 | ± 0.0048 | ± 0.023 | ± 0.065 | ± 30 |
3.1 Budget of U, Th and Pb in Acapulco
The phosphate separate was dissolved in two steps. Therefore, the successive dissolution steps, “Phos 1” and “Phos 2”, each one representing roughly half of the total sample, should be combined. The abundances for the combined phosphate are: 9.5 ppm Pb, 2.8 ppm Th and 5.1 ppm U; the 238U/204Pb ratio and the 232Th/238U ratio correspond to 132.5 and 0.57, respectively. The Th/U ratio allows us to discriminate between the two phosphates that are present in meteorites. Apatite exhibits low ratios (< 1) whereas whitlockite is characterized by high (∼10) values (Crozaz, 1974; Göpel et al., 1994). Based on this record we deduce that the mineral separate studied here consists mainly of apatite, which is in agreement with Zipfel et al. (1995).
The feldspar is characterized by a high Pb and a low U abundance, translating into a low 238U/204Pb ratio.
Because of an abnormally high procedural blank, the Pb isotope composition of the troilite separate is not well determined and will not be considered in the following. U and Pb abundances determined for this mineral, however, allow us to use them as upper limits in the budget calculation of the bulk meteorite presented below.
The grain size fractions display variable Pb, U, Th concentrations but similar Pb isotope compositions; this can be explained by slight variations in their modal mineralogy. In order to minimize this mineralogical effect, we calculate a more representative bulk sample for Acapulco. It is calculated based on the distribution of the different grain size fractions taking their weight proportion into account (details in section “sample preparation”). We assumed that the metal that was eliminated during the grinding step, contains negligible amounts of U, Th and Pb. Then this reconstituted bulk Acapulco material exhibits the following concentrations: 0.39 ppm Pb, 70.1 ppb U and 72.8 ppb Th. The atomic 238U/204Pb and 232Th/238U ratios correspond to 12.1 and 1.07, respectively. These U and Th concentrations are in agreement with the “Acapulco average” defined by Palme et al. (1981) and Zipfel et al. (1995). Compared to CI carbonaceous chondrites (Tatsumoto et al., 1976), the reconstituted bulk material of Acapulco is enriched in U and Th by a factor of 8.5 and 2.6, respectively, whereas 204Pb is depleted by a factor of 10.5. The 232Th/238U ratio, ∼1, is unusual for undifferentiated meteorites, in that their Th/U ratios range between 3.5 and 4.
Understanding the distribution of U, Th and nonradiogenic Pb among the different mineral phases is essential in interpreting the chronometric information. Two groups published the modal mineralogy of this meteorite (Palme et al., 1981; Zipfel et al., 1995). For our budget calculation, the following U, Th, Pb values are applied: bulk sample: 4.9 ppb 204Pb, 70.1 ppb U, 72.8 ppb Th; phosphate: 33 ppb 204Pb, 5.1 ppm U, 2.8 ppm Th; feldspar: 47 ppb 204Pb, 0.05 ppb U; troilite: < 1.8 ppb 204Pb, 3.7 ppb U.
Based on the mineral composition published by Palme et al. (1981) phosphates account for ∼100% U, 54% Th and 9.3% 204Pb. Feldspar contains 84% of 204Pb and < 0.01% U, troilite < 1% 204Pb and < 0.2% U. Using the modal mineral composition defined by Zipfel et al. (1995), we find that phosphate contains 73% U, 39% Th and 6.6% 204Pb, feldspar accounts for 90% of 204Pb and < 0.01 % U, troilite holds less than 2.9% 204Pb and less than 0.4% U.
Both evaluations show that phosphate and feldspar jointly represent the major hosts for U, nonradiogenic Pb and for a significant part of Th in the Acapulco meteorite. Therefore, the interpretation of the U/Pb chronology can be based on the U/Pb systematics measured in the bulk material, the feldspar and phosphate.
3.2 U/Pb isotope systematics
The Pb isotope composition of the feldspar lies very close to the primordial Pb composition that was measured in the Canyon Diablo iron meteorite (Göpel et al., 1985; Tatsumoto et al., 1973) (Table 1). This is the first time that such an unradiogenic composition is directly measured in undifferentiated high temperature material. The Pb isotope evolution of the solar system is almost linear during the first hundred million years. The composition of the Acapulco feldspar falls close to this “instantaneous” Pb evolution
that is issued from the primordial Pb composition and may be the result from a short (∼10 Ma) evolution of the primordial Pb composition. The nominal value of the initial Acapulco Pb composition, representing the Pb composition between minerals when the internal U/Pb equilibration occurred, is defined as the intersection between the ∼4.57 Ga Pb evolution line starting from the primordial Pb isotope composition with the internal Acapulco isochron: 206Pb/204Pb = 9.3336, 207Pb/204Pb = 10.3483, 208Pb/204Pb = 29.4533.
The Pb isotope compositions of the phosphate samples are radiogenic with 206Pb/204Pb values > 138. The 207Pb/206Pb model ages, assuming a primordial Pb composition (Tatsumoto et al., 1973) for the initial Pb, are similar: 4556.9 ± 1.7 Ma for Phos 1, 4555.3 ± 5.0 Ma for Phos 2. The 207Pb/206Pb model ages are slightly higher than the corresponding single stage 238U/206Pb and 235U/207Pb model ages indicating a discordance of 1.4% (Phos 1) and 1.0% (Phos 2). Two phosphates out of six published by Amelin (2005) also show a discordant behavior, that lies in the same order of magnitude.
The bulk material shows slightly radiogenic Pb compositions. Its 207Pb/206Pb model age is calculated assuming a single stage evolution starting from the primordial Pb composition. The 207Pb/206Pb age of the reconstituted Acapulco bulk material corresponds to 4562.0 ± 1.6 Ma and dates the closure of the U/Pb system at the whole rock scale. The 238U/206Pb ages (4815–5917 Ma) and 235U/207Pb model ages (4638–4945 Ma) are higher than the corresponding 207Pb/206Pb model ages, indicating an excess of radiogenic Pb. It corresponds in the reconstituated bulk Acapulco to 18%.
3.3 Internal 207Pb/206Pb age of Acapulco
The Pb isotope compositions of all samples define a line in the 207Pb/206Pb versus 204Pb/206Pb space (Fig. 1 and Table 2). The bulk sieved fraction plot on the isochron that is defined by the phosphate and feldspar. This shows that all minerals in Acapulco, beside the feldspar and phosphate, exhibit an identical 207Pb/206Pb age and that the addition of radiogenic Pb, present in excess in the bulk material, must have occurred recently, most probably during the excavation of the meteorite from its parent body (Manhès et al., 1984).
Isochrones définies par régression linéaire des analyses d’Acapulco.
Samples used for regression calculation | Slope | Intercept | MSWD | Probability | Age [Ma] |
fr100 a, fr200, fr100 b, phos1, phos2, feldspar | 4.5593 | 0.62033 | 0.65 | 0.63 | 4556.4 |
± 0.0093 | ± 0.00046 | ± 1.1 | |||
fr100 a, fr200, fr100 b, phos1, phos2, feldspar, 6 phoshates* | 4.5622 | 0.62014 | 0.85 | 0.58 | 4555.9 |
± 0.0067 | ± 0.00026 | ± 0.6 | |||
phos1, phos2, feldspar | 4.5594 | 0.62022 | 0.39 | 0.53 | 4556.1 |
± 0.0094 | ± 0.00066 | ± 1.5 |
Different linear regressions, including or not including the grain size fractions into the calculation, are reported in Table 2. The similarity between the different ages confirms the isotopic consistency between bulk material and mineral separates. The 207Pb/206Pb age based on all samples analyzed in this study corresponds to 4556.4 ± 1.1 Ma; the age calculated only with phosphate data from Amelin (2005) corresponds to 4556.5 ± 1.3 Ma. When the phosphate data from Amelin (2005) and our analyses are used together in this regression, the 207Pb/206Pb age is defined with a smaller uncertainty, 4555.9 ± 0.6 Ma. This age records the end of the U/Pb re-equilibration between the minerals of Acapulco and more specifically it indicates the end of the Pb diffusion outside of the phosphate.
4 Discussion
The Pb isotope systematics of Acapulco bulk material, phosphate and feldspar are consistent: U/Pb closure at the whole rock scale occurred at ∼4562 Ma, whereas the internal U/Pb re-equilibration occurred or ended at 4556 Ma, 6 Ma later. In the following, we discuss the different events that are registered by the U/Pb system.
4.1 Enrichment of U and Th in bulk Acapulco
Zipfel et al. (1995) point out that the Acapulco bulk exhibits large excesses of LREEs and U and show that phosphates are the main carriers for these elements. However, these authors explicitly rule out a direct addition of phosphate. Such an addition would have concomitantly induced major gains in Ca and P, but the Ca/Mg and P/Mg ratios of bulk Acapulco are approximately chondritic. Therefore, they propose that a melt, rich in incompatible elements, infiltrated Acapulco at such low temperature (∼800 °C) that equilibrium among different minerals could not be reached. LREEs and U were extracted by phosphates from this melt while the remaining melt, that was rich in all other incompatible elements, subsequently drained away and left Acapulco (Zipfel et al., 1995). Using Zipfel's model, we estimate the amount of U and Th that were added to the bulk sample. Compatible refractory elements (Sc, Al, Y) are enriched in Acapulco by a factor of 1.43 with respect to CI material. We assume that Acapulco prior to the melt intrusion is equally enriched in Th and U and characterized by a Th/U of 3.5 similar to CI. The calculated abundances of the Acapulco precursor correspond to 11.7 ppb U and 40.9 ppb Th. The amount of U and Th that had been added by the melt is calculated by subtracting these abundances from the measured U (70.1 ppb) and calculated Th abundances (72.8 ppb) and corresponds to 58.4 ppb U and 31.9 ppb Th. Both elements, U and Th, have been added to the bulk rock with a low Th/U ratio (0.55); this ratio is similar to that found in the Acapulco phosphates. Therefore, we propose that the melt selectively mobilized U, Th and LREEs out of the phosphates from the Acapulco parent body, these elements migrated and were incorporated during the formation of apatite. No other incompatible elements have migrated into Acapulco.
4.2 207Pb/206Pb age of Acapulco
Two questions arise when one tries to elucidate the early events that led to the Pb depletion of the Acapulco bulk material:
- • Did the melt that infiltrated Acapulco add Pb together with U, Th and LREEs? No experimental data are available that allow us to denote the behavior of Pb.
- • Is the Pb depletion related to nebular processes and/or to Pb loss during the high temperature regime in the Acapulco parent body?
Based on the U/Pb record alone, we cannot answer this question. The age of 4562 Ma of bulk Acapulco may be related to nebular processes (partial condensation or partial evaporation of solar dust), or to planetary processes (thermal closure during the high temperature regime of Acapulco, incorporation of U and eventually Pb by melt percolation). If the melt did not incorporate Pb at all, the 207Pb/206Pb age of Acapulco bulk material dates the melt percolation, because the majority of U was added by this melt. On the other hand, if the melt mobilized Pb as well as U and if Acapulco incorporated Pb and U in similar proportions, this 207Pb/206Pb age does not register the melt percolation but the age is related to a U/Pb fractionation that occurred previously.
The comparison between U/Pb and 182Hf/182W data obtained on two acapulcoites partially resolves the ambiguity concerning the meaning of 207Pb/206Pb age of the Acapulco bulk (Touboul et al., 2009). The “absolute” Hf/W age of these samples corresponds to 4563.1 ± 0.8 Ma, this age indicates the end of the Hf/W equilibration between metal and silicate (Touboul et al., 2009). Consequently, the younger 207Pb/206Pb age of the Acapulco bulk is not related to the end of equilibration with the nebular gas but it dates the U/Pb closure at the whole rock scale or it dates the melt intrusion during the cooling of the Acapulco parent body.
4.3 Cooling history of the Acapulco parent body
Comparing the Hf/W and U/Pb thermochronometries with mineralogical cooling rates allows us to better constrain the thermal history of the Acapulco parent for the temperature range ∼1000–350 °C. This is shown in Fig. 2.
The internal 207Pb/206Pb age of Acapulco corresponds to the end of Pb diffusion out of the apatite. This closure temperature is deduced from Pb diffusion experiments on terrestrial apatites (Cherniak et al., 1991) taking the crystal size of the phosphates (∼200 μm) into account. For cooling rates of 1, 10, 102 and 103 °C/Ma the closure temperature corresponds to ∼445, 490, 540 and 585 °C, respectively. The closure temperature of the U/Pb system lies between 650 °C and 350 °C, where high cooling rates, at least 103 °C/Ma, are recorded by mineralogical investigations. Therefore, the internal 207Pb/206Pb age precisely dates when Acapulco cooled rapidly. This is not the case for the other chronologies (Hf/W, K/Ar and U/Th–He) whose closure temperatures lie outside the temperature range 650–350 °C where mineralogical tools record the fast cooling. The mean cooling rate between the closure of the Hf/W system at 4563.1 ± 0.8 Ma and the U/Pb system at 4555.9 ± 0.6 Ma can be calculated. This time interval (7.2 ± 1.0 Ma) and the closure temperatures of the Hf/W system (985 ± 50 °C) and U/Pb system (525 ± 50 °C) define a mean cooling rate of ∼ 65 ± 15 °C/Ma. This cooling rate is more than one order of magnitude slower than the minimum cooling rate of 103 °C/Ma suggested by Min et al. (2003). Such distinct cooling rates in an overlapping temperature range (650–520 °C) are incompatible with a model of a continuously cooling and undisrupted parent body. Combining the information derived from the U/Pb and Hf/W system with mineralogically determined cooling rates requires the break-up of the Acapulco parent body, as proposed by Min et al. (2003) and Renne (2000). Therefore, the cooling rate of Acapulco before the disrupture of the parent body is not constrained and may vary from a few up to ∼102 °C/Ma. The size of the Acapulco parent body and the position of Acapulco inside this parent body (i.e. near the surface or deeply embedded) is less constrained than suggested by Touboul et al. (2009) whose model is only based on radiochronometric records.
The time when the parent body broke up can be precisely established from the U/Pb closure age (4555.9 ± 0.6 Ma) and the cooling rate (103 °C/Ma). This event may have occurred 0.4 Ma before the closure of the U/Pb (if Acapulco's cooling rate corresponded to ∼10 °C/Ma) or at the time of the U–Pb closure. Therefore the break-up of the Acapulco parent body is dated at 4555.9 +0.8/−0.6 Ma.
Acapulco may have continued to cool rapidly until 120 °C inside a small fragment as supported by Min et al. (2003) and Renne (2000). On the other hand, when the ages obtained from 244Pu fission track studies and Ar ages calculated with the decay constant from Steiger and Jäger (1977) are taken into account, then Acapulco must have been buried in a re-agglomarated large body in order to explain the low cooling rate, < 2 °C/Ma, between 350 °C and the retention of 244Pu fission tracks at ∼120 °C.
5 Conclusion
The U–Th–Pb systematics of separated minerals and in bulk fractions of Acapulco define major events of the chemical and thermal history of its parent body.
For an undifferentiated meteorite, Acapulco exhibits high U and Th abundances with an unusual Th/U ratio of ∼1. It was proposed that the main part of these elements was added early to the Acapulco precursor by a low temperature melt. We support this idea and suggest that this melt selectively extracted U, Th and LREEs from phosphates. Subsequently, these elements migrated into Acapulco and were incorporated by phosphates when they crystallized. The low Th/U ratio determined for the Acapulco bulk, is the result of this melting event. It occurred between 4562 Ma and 4556 Ma.
The feldspar contains the frozen Pb composition of Acapulco at 4556 Ma, it lies close to the primordial Pb composition (Göpel et al., 1985; Tatsumoto et al., 1973). The small difference with this primordial composition can be explained by the decay of U and Th during ∼10 Ma, before 4556 Ma, the end of U/Pb re-equilibration in Acapulco.
The comparison between Hf/W ages determined in acapulcoites and the U/Pb age in Acapulco as well as cooling rates determined with mineralogical tools for Acapulco support the idea that the break-up of the parent body occurred during the cooling period. This event is precisely dated with the U/Pb system at 4556 ± 1 Ma. The initial size of the parent body and/or the position of Acapulco inside this parent body are poorly constrained.
Acknowledgements
The authors would like to thank two anonymous reviewers for their comments that helped to improve the final version of the manuscript. This work was financially supported by ANR-08Blan-0260. This is IPGP contribution number 2560.