(2.2)
where (41K/39K)reference is the average of the two bracketing analyses of the K reference solution prior to and after the analysis of the sample. The sample solutions were introduced to the IsoProbe using using an Aridus II (Cetac Technologies) desolvation system with a Teflon nebulizer with a measured uptake of 60 microlitres/minute. Approximately 150 ng of K was used per isotopic analysis run. The typical uncertainty in δ 41K is conservatively given as ± 0.3‰ and is based on the combined uncertainties in the measured 41K/39K of the bracketing reference and samples. Further details on the mass spectrometry technique can be found in Richter et al. (2011).
2.3. RESULTS
2.3.1. Major Element Diffusion Profiles
Major‐element diffusion profiles from two experiments that differ in duration (2.5 versus 6 hours) are shown in Fig. 2.3. Note that the x‐coordinates of both profiles have been shifted such that the phonolite‐rhyolite interfaces overlap and the profiles are framed within each panel. Additionally, in the 6‐hour experiment the capsule was loaded with more phonolite powder (0.07 g) than rhyolite powder (0.04 g), and consequently, the rhyolite side of the couple is short, occupying only about one‐third of the capsule length.
The major element profiles exhibit features that are interesting, although we cannot fully explain some of them. We observe uphill diffusion (i.e., diffusion of a component against its own concentration gradient) in all the major‐element oxides, including SiO2. The SiO2 and Al2O3 profiles change relatively little from the starting materials, which is expected because of the low diffusivities of these components and the coupling between Si and Al in the silicate framework components of the liquids. One thing that is unexpected is that the absolute concentrations of SiO2 and Al2O3 on the rhyolite end of the profiles are low relative to the starting values. This should be regarded as an analytical inconsistency as it cannot be explained by a dilution effect. The absolute concentrations of Na2O are even more problematic, being too high relative to the starting materials. The simplest explanation is that the starting materials actually had Na2O concentrations that were about 0.5% higher than reported. This is unexpected, but possible, since it is well known that Na can be volatilized during EMP analysis, leading to incorrectly low measured concentrations. These issues do not significantly affect the main conclusions that we draw from the Ca and K isotopic results.
A notable feature of the data is the dramatic rearrangement of the Na, Ca, and K components in the 2.5 hour run. The K profile is strongly affected by rapid diffusion of K from the phonolite to the rhyolite, producing a steep K concentration profile in the rhyolite. The large increase in K in the rhyolite, especially close to the contact with phonolite, is accommodated mostly by rapid diffusion of Ca and some Na from the rhyolite to the phonolite. Apparently, as the Ca gradient becomes steep within the rhyolite, Na diffuses both into the phonolite, and toward the rhyolite end of the capsule to maintain a near‐constant Na/Ca within the rhyolite. We have not attempted to capture this coupling in our modeling of Ca (see Section 2.5), but since Na diffuses rapidly, its movement probably does not greatly affect the diffusion of Ca. The Na/Ca ratio (or Ca/(Ca+Na) ratio) is nearly constant within the rhyolite, and within the phonolite, which is a clear indication that both Na and Ca diffuse quickly to what must be thermodynamically favored Na/Ca ratios that are presumably determined largely by the SiO2 (± Al2O3) of the two liquids.
2.3.2. Ca and K Isotopes
The phonolite‐rhyolite diffusion couples produced large diffusive isotope effects in both Ca (Fig. 2.4a) and K (Fig. 2.4b). For Ca isotopes, the felsic side is enriched in the heavy isotope whereas the mafic side is enriched in the light isotope, consistent with a net flux of CaO from the felsic to the mafic liquid despite a negligible initial compositional contrast. In the 2.5‐hour run, the δ 44Ca is equal to the initial value (within error) on the felsic side. It is unclear if the δ 44Ca has changed from the initial value on the mafic side because the endmost wafer was not successfully extracted during sample preparation. The total range in δ 44Ca is about 2 in the 2.5‐hour run. In the 6‐hour run, the shape of the profile is similar but more time‐evolved. The Ca isotopic composition has changed significantly from the initial value at both ends of the capsule with a total range in δ 44Ca of about 1.5‰.