## Abstract

It is common practice to estimate a mean diameter for spherical or sub-spherical particles or vesicles in a rock by multiplying the average diameter of the approximately circular cross-sections visible in thin section by a factor of 1.273. This number-weighted average may be dominated by the hard-to-measure fine tail of the size distribution, and is unlikely to be representative of the average particle diameter of greatest interest for a wide range of geological problems or processes. Average particle size can be quantified in a variety of ways, based on the mass or surface area of the particles, and here we provide exact relations of these different average measures to straightforward measurements possible in thin section, including an analysis of how many particles to measure to achieve a desired level of uncertainty. The use of average particle diameter is illustrated firstly with a consideration of the accumulation of olivine phenocrysts on the floor of the 135 m thick picrodolerite/crinanite unit of the Shiant Isles Main Sill. We show that the 45 m thick crystal pile on the sill floor could have formed by crystal settling within about a year. The second geological example is provided by an analysis of the sizes of exsolved Fe-rich droplets during unmixing of a basaltic melt in a suite of experimental charges. We show that the size distribution cannot be explained by sudden nucleation, followed by either Ostwald ripening or Brownian coalescence. We deduce that a continuous process of droplet nucleation during cooling is likely to have occurred.

## Introduction

Useful answers to many geological problems can be obtained from relatively simple calculations that provide time- or length-scales, correct within an order of magnitude, which can then be used to place constraints on the processes likely to have been involved in the problem in question. Good examples of this kind of approach are based on determinations of grain size, with the quantification of particle size in rocks (either grains or bubbles) providing an opportunity to make progress on many problems of petrological interest, such as magma solidification time-scales (Cashman and Marsh, 1988; Cashman, 1993; Higgins, 1996), crystallization during magma ascent (Cashman, 1992; Hammer *et al*., 1999), rates of production and coalescence of volatile-filled bubbles from magma (Herd and Pinkerton, 1997), buoyancy-driven particle migration or other fluid dynamical processes (Robertson and Barnes, 2015), rates of Ostwald ripening (Cabane *et al*., 2001, 2005) and pattern formation in metamorphic rocks (Holness, 1997).

Much recent work using grain size to quantify geological processes is based on a sophisticated treatment involving the characterization and interpretation of the particle size distribution (as introduced by Marsh, 1988). The accuracy of such an approach is enhanced by determination of the true 3D distribution of grain sizes by disaggregation, dissolution of the matrix (e.g. Holness, 1997) or tomographic analysis (e.g. Carlson and Denison, 1992; Denison and Carlson, 1997). However, given the limitations of the materials we work with, most studies using grain size are based on observations of thin sections, in which case stereological corrections are required to convert the range of grain intersection size to a 3D grain-size distribution (e.g. Cashman and Marsh, 1988; Johnson, 1994; Higgins, 2000).

For spherical particles, converting the distribution of circular cross-sections observed in thin section into an estimate for the true 3D distribution of particle diameters is mathematically well defined (Wicksell, 1925). For non-spherical particles, such as parallelepipeds, available numerical methods are based on the assumption of invariant particle shape regardless of size (e.g. Higgins, 1994; Sahagian and Proussevitch, 1998; Morgan and Jerram, 2006), which is not likely to be true for natural samples (e.g. Mock and Jerram, 2005; Duchêne *et al*., 2008). However, for many applications a mean particle diameter is often sufficient to provide order of magnitude estimates that can be used to constrain timescales of geologically interesting problems. The question then arises as to how one might obtain a mean particle diameter from thin-section observations.

For the particular case of a monodisperse population of spheres (i.e. one with a uniform particle size), the average diameter of circular cross-sections obtained by random cross-sections through the population is π/4 times the sphere diameter. The simplicity of this relationship has led to its common application to estimate an average 3D particle diameter for polydisperse (i.e. a population with a range of 3D particle sizes) as well as monodisperse particle distributions (Hughes, 1978; Cashman and Marsh, 1988; Kong *et al*., 1995; Herd and Pinkerton, 1997), although model-based maximum likelihood approaches are also used (Kong *et al*., 1995).

In this contribution we concentrate on systems containing spherical particles, such as bubbles, droplets in an emulsion, or equant mineral grains such as olivine or spinel, and argue that such a simple approach to determining the average particle diameter has three problems, which can be remedied easily. Firstly, the average value obtained using this method may be affected strongly by the smallest particles in the population, which is precisely the part of the size distribution that is most likely to be overlooked or not properly resolved. Secondly, this approach provides no estimate of the uncertainties in the result. Thirdly, and perhaps most importantly, it doesn't address a question of great significance for polydisperse particle populations which is absent for the monodisperse case, namely: which of the various ways of calculating the average diameter is most appropriate for the problem we are interested in?

Firstly we discuss the various merits of different ways of calculating the average for sphere diameters, present some simple exact results linking them to circular cross-sections, and provide simulated data to show how many grains need to be measured to achieve any required degree of accuracy. We then explore how sensitive these statistics are to ignoring the smallest cross-sections in a sample, and whether the proposed method can be applied to non-spherical, but equant grains (specifically, we look at cubes). Lastly we illustrate the usefulness of various measures of the average particle size to constrain timescales of settling of olivine grains on the floor of a basaltic sill, and the mechanisms of coarsening of an unmixed immiscible basaltic melt.

## Calculating the average

_{sph}spherical particles with diameters , where . One way to characterize this distribution is to find the ‘number-weighted’ mean diameter which we will term

*D*

_{1,0}. This is simply the sum of all the sphere diameters divided by the number of spheres: However, this simple average may be problematic, particularly if the grain size distribution is broad, as an abundance of very small grains will dominate this average so that

*D*

_{1,0}may be much smaller than we expect. Additionally, it may be difficult to identify all the tiniest particles (for example, within the limits of the X-ray tomographic analysis of any particular sample), so that even if we were happy to use

*D*

_{1,0}the calculated value might be prone to large errors.

*D*

_{4,3}and the surface-weighted mean diameter

*D*

_{3,2}: The volume-weighted mean diameter can be defined as follows: choose random points in space, and if the point lies within one of the particles of interest, record the diameter of that particle.

*D*

_{4,3}is the mean of all the diameters sampled this way. Because large spheres are much more likely to be sampled than small spheres, this average is dominated by the larger spheres in the distribution, and much less sensitive to under-representation of the tail of fine dust. The area-weighted mean diameter

*D*

_{3,2}can be thought of as taking a mean of diameters where spheres in the sample are chosen with a probability proportional to their surface area (rather than volume as in

*D*

_{4,3}). It also is less sensitive to the tail of fine particles than the simple number-weighted average.

The choice of average diameter depends on which captures the relevant properties of the system under investigation. *D*_{4,3} represents the size class around which most of the mass of the particles lies, and for that reason may be taken as a good measure of particle diameter from a compositional point of view. In contrast, if the problem under consideration involves processes controlled by interfacial area (for example the aggregation of crystals to make sintered clusters, or the adsorption of water onto the surface of soil particles) then would be the best measure of average particle dimension. This is because the amount of surface area *S* per unit volume of sample is simply , where ϕ is the volume fraction of spherical grains. If we are interested in the permeability of rocks then we note that the Kozeny-Carman relation (Carman, 1937) gives an approximate expression for the hydrodynamic permeability in terms of *S*, so that is once more a key quantity. Hydrodynamic permeability may also be relevant to sedimentation of concentrated suspensions, as the rate of sedimentation is likely to be determined by D'Arcy flow through the bed as a whole, rather than by particles settling individually. In Ostwald ripening, crystal growth is driven by interfacial energy, but in the scaling regime of LSW theory (Lifshitz and Soyolov, 1961; Wagner, 1961) *all* the mean diameters have the same cube root dependence on time, so there is no obvious preferred choice in this case. Comparison of several mean diameters can, however, be used to shed light on whether Ostwald ripening is the dominant growth mechanism. For settling of dilute suspensions, we show below that is the more relevant parameter to calculate. As a less geologically relevant aside, statistical studies suggest that when observers look at cross-sections of different sphere distributions, they tend to rank them by size according to (Alderliesten, 2008); a correlation which has yet to be given a rigorous (physiological) explanation.

*D*and is , where (Farr, 2013) and , then the logarithmic width σ (a simple measure of the range of particle sizes present in the rock) can be calculated from

## Sections through sphere distributions

*d*'s refer to two-dimensional sections and upper case

*D*'s to three-dimensional spheres), we can use these data to define the various mean circle diameters, in the same way as for the different mean sphere diameters. For example and would be Provided the number of measured grains in the sample is very large, simple integration (see Appendix) shows that the mean sphere diameters can be obtained from different mean circle diameters thus: and the number,

*N*, of spheres per unit volume can be deduced from the number,

*n*, of circles per unit area through where is the harmonic mean of the circle diameters (the reciprocal of the arithmetic mean of the reciprocals of the diameters).

Although these relations are independent of the size distribution of spheres, they are only exact in the limit of an infinite number of individual measurements. For a finite number of measurements of individual particles, there will be some scatter in results if different particle populations are measured (e.g. different parts of the same thin section, or different thin sections of the same sample), and there may also be some systematic error in the mean taken over many realizations of the experiment. As an extreme example of systematic bias, if only one grain is measured (i.e. if ), all the different averages , etc. would be identical, leading to the absurd prediction , which can never be true for any sphere size distribution.

## Confidence intervals: how many circular sections to measure?

In this section we address the question of how many particles should be measured to achieve a given level of confidence in the estimates for average grain diameter. The results of our analysis depend on the distribution of sizes and, for simplicity, we assume this to be lognormal. The different lognormal distributions we consider are shown in Fig. 1.

Suppose we measure the diameters of particles in a thin section, and then calculate estimates for , and using equations 11 and 12. If this procedure were performed several times, each time measuring a different population of particles in the thin section, the results would have some scatter, due to statistical fluctuations, which (together with any systematic bias) will give an estimate of the uncertainty in the result.

Figure 2 shows results where we have used a computer to generate monodisperse spheres in random locations in space and generated circular sections from a plane drawn through this distribution. No account is taken of sphere overlaps, so the simulations strictly represent the dilute limit. The plots show the 68% confidence intervals for the predicted quantities compared to their true values (which are known in this case). The confidence intervals mean that 68% of the results lie in the interval; 68% being chosen because for a normal distribution this would represent plus or minus one standard deviation.

If we define the fractional error of the method to be the difference between the least accurate point in the confidence interval and the true value, divided by the true value, then we can plot this fractional error also for different lognormal sphere size distributions. This is done in Fig. 3, where we see that in general the fractional error is proportional to , as would be expected from the central limit theorem, but overall it is harder to accurately measure parameters for the wider sphere size distributions.

## Sensitivity to ignoring the smallest circular cross sections

The smallest cross-sections in a sample may fall below the limit of resolution in an image, so it is important to know how sensitive are the statistics we propose to the omission of the small tail of the circle distribution. Figure 4 shows similar data to Fig. 3, but where some of the small circles have been omitted from the statistics, specifically, all those circles with diameters less than some fraction α of . We see that estimates of itself are not affected materially even when α is as large as 0.2, while estimates for are a little more sensitive, and require for there to be no measurable effect. This lack of sensitivity, of both the volume-weighted and area-weighted averages, to an under-representation of the smallest particles is unsurprising. Therefore for the statistics we propose here (, and their ratio), a sensible rule of thumb would be to check *post hoc* that all circles larger than one fifth of the calculated value of have been included in the averages, and preferably all those larger than one tenth of .

## Does the procedure work for cubic crystals?

For non-spherical particles, we define the equivalent diameter of a particle as the diameter of a sphere that has the same volume. Similarly for the cross-sections, we define the equivalent circle diameter as the diameter of a circle with the same area as the cross-section of the particle (which for a cubic particle will be a polygon with 3 to 6 sides; see Higgins (1994) and Morgan and Jerram (2006) for examples).

Figure 5 shows the errors incurred when the proposed procedure in this paper (derived for spheres) is applied to a random distribution of randomly oriented cubes. While the errors do not die out at , and never fall below ∼3%, this procedure can indeed be used to obtain reasonable estimates of the equivalent mean diameters for the cube population.

We note however that although values for equivalent , and their ratio are well predicted, it would not be appropriate to apply the sphere result of to estimate the specific surface area in the system. The correct expression for cubes results in a surface area about 24% higher (the surface area of a cube being 24% higher than a sphere of the same volume). Instead, the standard stereological method (Russ, 1986) that the specific surface area is times the perimeter per unit area in the cross section would be the appropriate method for this statistic.

## Application to gravitational settling of a polydisperse grain population

An almost universal process occurring during the solidification of basaltic magma is the relative movement of crystals and residual liquid under the influence of gravity. It is this process which is the fundamental driver for fractionation. Here we discuss a simple treatment of settling under gravity of an initially dilute suspension of crystals entrained in a basaltic magma. If we assume that emplacement of the crystal-bearing magma is essentially instantaneous, and that the crystal-bearing magma has a Newtonian rheology, the initial buildup rate of the thickness of the layer of crystals on the magma chamber floor for a dilute suspension of polydisperse spheres in a non-convecting Newtonian liquid can be calculated as follows:

*N*spheres per unit volume in suspension and let the fraction of the number of these spheres that have diameters between

*D*and be , so that is a normalized probability density: . These spheres are initially randomly distributed in a liquid column above the floor, and we assume that their growth during settling is negligible. If all the spheres were of diameter

*D*then due to Stokes' law (the suspension is dilute) they would all settle at a speed where η is the fluid viscosity, g the acceleration due to gravity and the difference in density between spheres and liquid. Thus, in unit time, the volume of spheres whose centres pass through a horizontal plane in the fluid per unit area of that plane will be which represents the instantaneous buildup of sphere volume on the floor. For a range of sphere sizes, this rate will be The number

*N*of spheres per unit volume is related to the volume fraction in suspension through , while the rate of buildup of the layer thickness is related to through , where is the volume fraction of spheres in the sediment at the time they are deposited. Combining these with equation 20 and 2 we obtain . However, the size distribution applies to the spheres in suspension, while we only have access to the size distribution in the sediment. Suppose there are spheres per unit volume in the suspension with sizes close to . The number per unit volume in the sediment will be increased in proportion to the speed of sedimentation. Therefore in the suspension (at this moment in time) will be equal to in the sediment, and so we obtain where is now obtained from the distribution of sphere sizes in the sediment (at the time of deposition). Note that this is a local measure at this point in the accumulating sediment, and a more global analysis capturing variations with sediment depth could be performed, based on the conservation laws of sedimentation in partial differential form (Simakin

*et al*., 1998).

As settling proceeds, the suspension becomes less dilute near the floor and the particles become closer together. During the final stages of sedimentation therefore, Stokes' Law no longer holds and sedimentation rates become controlled more by the permeability of the particle accumulation.

## Olivine settling in the Shiant Isles Main Sill

To illustrate our treatment of settling timescales using average particle size we focus on the Shiant Isles Main Sill, which is the largest of the four separate Tertiary alkaline basalt sills exposed on the Shiant Isles (Outer Hebrides, Scotland), and intruded into Jurassic sediments (Gibb and Henderson, 1984). The Shiant Isles Main Sill is 165 m thick (Gibb and Henderson, 1984) and is a composite body (Drever and Johnston, 1959; Gibb and Henderson, 1989, 1996). The bulk of the sill (135 m of stratigraphy) formed from a single pulse of olivine-phyric magma that contained 10 vol.% olivine phenocrysts together with 1–2 wt.% Cr-spinel and a small amount of plagioclase (Gibb and Henderson, 2006). The olivine phenocrysts settled to the (contemporary) floor of the intrusion to form a picrodolerite ∼45 m thick, leaving an essentially aphyric magma that crystallized to form the remainder (a crinanite, dominated by plagioclase and augite, with interstitial olivine, Fe-Ti oxides and analcime). Here we concentrate on the olivine accumulation that forms the picrodolerite. The olivine grains are generally equant and rounded, commonly forming clusters and loose chains in which the individual grains are joined by small areas of grain boundary.

*ImageJ*software, and the diameters of circles with the same areas were then calculated. A total of 517 olivine crystals were identified. The statistics of the population are shown in Table 1. From these, using equations 14 and 15 we obtain and . Equations 16 to 18 then give the following estimates: As , which is much smaller than one fifth, we can be confident that these estimates are not affected by any smaller grains which were not resolved. We also note from equations A12 and A10 (see Appendix) that so, assuming has the same relative error as and , If we assume that the only mineral accumulating on the sill floor was olivine, we can calculate the original volume fraction of settled grains if we also assume that the interstitial liquid crystallized as a mixture of interstitial plagioclase, augite and overgrowths on the existing olivine grains. This means that the observed grains are all somewhat larger than the original settled population, and so a correction must be made for this post-accumulation overgrowth.

The smallest and largest equivalent diameters in the sample are and .

*et al*. (2000)). The proportion of Sample SC459 that is olivine is 48 vol.%, which suggests that if the volume fraction of olivine in an olivine-only sediment at the time of deposition were then , so and the change in volume fraction to remove the overgrowth is thus . From this and equations 28 and 26 we estimate that removing the overgrowth reduces the value of , but only by the insignificant amount of 3 µm, so we can continue to use the original value of .

The corrected olivine mode in SC459 is lower than the 56–54 vol.% expected for a random loose packing of cohesionless monodisperse spheres (Onoda and Liniger, 1990; Ciamarra and Coniglio, 2008; Zamponi, 2008; Farrell *et al*., 2010), and lower still than random loose packings achieved for polydisperse particles (Epstein and Young, 1962; Jerram *et al*., 2003), suggesting that olivine was not the only phase settling from the incoming magma. However, the efficiency of random loose packing is reduced for strongly cohesive particles, for which a stable distribution can be achieved at lower volume fractions (Dong *et al.*, 2006; Yang *et al.*, 2007). The presence of highly non-spherical, loose clustered chains of olivine will reduce this still further (Campbell *et al*., 1978; Jerram *et al*., 2003), with packings as low as 37 vol.% observed for settled accumulations of loose chains and clusters of olivine and magnetite (Campbell *et al*., 1978),. If we assume that the crystal pile was not densified by compaction, shear or shaking, it is therefore plausible that the accumulated olivine grains in SC459 preserve a randomly loose packed, mechanically stable framework of crystals and loose crystal clusters and chains with an overall solid fraction of ∼45 vol.%.

_{80}(Gibb and Henderson, 1996), which has a density of 3.29 g/cm

^{3}. Gibb and Henderson (2006) provide two suggested compositions for the parental magma that formed the 135 m picrodolerite/crinanite unit. Assuming an intrusion temperature of 1200°C (Gibb and Henderson, 2006), a confining pressure of 0.5 kbar and a H

_{2}O content of 1 wt.%, the two parental magmas have a density of 2.6 g/cm

^{3}, calculated according to the method of Bottinga and Weill (1970), using volume and thermal expansion data from Lange and Carmichael (1987); Kress and Carmichael (1991). Therefore the density difference is . The largest uncertainty relates to the viscosity of the liquid phase through which sedimentation occurs. From Giordano

*et al*. (2008) the two compositions lead to predicted viscosities of 3.9 and 8.5 Pa. Using the above properties in equation 22, we obtain two values for the rate of buildup of the sediment in the magma chamber: The layer of sedimented olivine is 45 m thick (Gibb and Henderson, 1989, 1996) so, provided convection did not significantly delay the process, we predict a time of around 280 or 600 days for the layer to accumulate. This is much faster than the cooling time of the sill: a simple upper bound on the time for solidification can be obtained by assuming the entire thickness of the 165 m sill is intruded into cold country rock, which results in an estimate of 90 years from emplacement to complete crystallization (using the simple calculation developed by Holness

*et al*., 2012).

Note that although the volume fraction of the settled olivine is ∼45 vol.% in the lower part of the picrodolerite, the olivine mode decreases upwards towards the crinanite (Gibb and Henderson, 1996). It is probable that this reduction is matched by an increase in the accumulation of another phase, such as plagioclase. The implications of this will be explored in a future contribution.

## Mechanisms of droplet growth in phase-separating magmas

Silicate liquid immiscibility in basaltic systems, first recognized by Roedder and Weiblen (1971) (with significant further observations by Philpotts (1979, 1982)), is recognized increasingly as an important factor controlling fractionation and the compositions of erupted magmas (Veksler *et al*., 2007; Charlier *et al*., 2011). The potential for immiscibility to affect the liquid line of descent on the scale of a magma chamber depends on the ease with which the two conjugate liquids can separate under the influence of gravity (e.g. Holness *et al*., 2011; VanTongeren and Mathez, 2012) and this, in turn, is affected strongly by the coarseness of the emulsion (e.g. Chung and Mungall, 2009). At present almost nothing is known about the kinetics of emulsion coarsening (e.g. Martin and Kushiro, 1991; Veksler *et al*., 2010). However the size distribution of the droplets potentially carries information about the mechanisms of their formation and subsequent growth.

One possible scenario for emulsion formation is that there is an initial interval when many nuclei form, followed by evolution of the structure without further nucleation. This could happen, for example, with spinodal decomposition, or if the system is between the binodal and spinodal lines and there are many potential sites for heterogeneous nucleation. Once droplets have formed in such a system, they can grow by various mechanisms. In a non-convecting liquid, growth can happen either by Ostwald ripening, or through coalescence as droplets diffuse and collide under Brownian forces. Both of these mechanisms lead to a state where all the mean diameters grow as the cube root of time (Crist and Nesarikar, 1995), but each has a characteristic size distribution, which can be probed by moment-based averages. From the derivations in Appendix B, we see that Ostwald ripening leads asymptotically to , whereas the size distribution resulting from coalescence under Brownian forces is characterized by .

We analysed the sizes of droplets in experimental charges described previously by Charlier and Grove (2012). The charges, containing material representative of compositions of tholeiitic basalts from the Sept Iles intrusion (charges SI-5, SI-8 and SI-13, Table 2) and Mull (charges M-5, M-6 and M-9, Table 2), were cooled at 1°C per hour from a starting temperature of 1100°C to a range of temperatures at which they were held for periods of up to 96 hours before quenching. At the end of each experiment the charges contained solid phases (detailed in Table 2 of Charlier and Grove (2012)) together with Si-rich glass containing exsolved quenched droplets of Fe-rich liquid (Fig. 7). Using back-scatter images we measured the diameters of isolated droplets in large regions of glass, avoiding those in direct contact with the mineral phases (thus avoiding droplets that may have nucleated heterogeneously on the mineral surfaces).

Sample | Final T (°C) | Equilibration time (hours) | D_{4,3}/µm | D_{3,2}/µm | D_{4,3}/D_{3,2} | |
---|---|---|---|---|---|---|

SI-5 | 282 | 1006 | 96 | 3.94 ± 1.54 | 2.22 ± 0.87 | 1.77 ± 0.69 |

SI-8 | 484 | 963 | 48 | 2.26 ± 0.30 | 1.59 ± 0.21 | 1.42 ± 0.19 |

SI-13 | 299 | 1020 | 64 | 1.36 ± 0.04 | 1.28 ± 0.04 | 1.06 ± 0.03 |

M-5 | 3205 | 1005 | 92 | 2.42 ± 0.10 | 1.80 ± 0.07 | 1.35 ± 0.05 |

M-6 | 621 | 963 | 48 | 1.85 ± 0.11 | 1.51 ± 0.09 | 1.23 ± 0.07 |

M-9 | 2017 | 1020 | 64 | 2.22 ± 0.09 | 1.75 ± 0.07 | 1.27 ± 0.05 |

Sample | Final T (°C) | Equilibration time (hours) | D_{4,3}/µm | D_{3,2}/µm | D_{4,3}/D_{3,2} | |
---|---|---|---|---|---|---|

SI-5 | 282 | 1006 | 96 | 3.94 ± 1.54 | 2.22 ± 0.87 | 1.77 ± 0.69 |

SI-8 | 484 | 963 | 48 | 2.26 ± 0.30 | 1.59 ± 0.21 | 1.42 ± 0.19 |

SI-13 | 299 | 1020 | 64 | 1.36 ± 0.04 | 1.28 ± 0.04 | 1.06 ± 0.03 |

M-5 | 3205 | 1005 | 92 | 2.42 ± 0.10 | 1.80 ± 0.07 | 1.35 ± 0.05 |

M-6 | 621 | 963 | 48 | 1.85 ± 0.11 | 1.51 ± 0.09 | 1.23 ± 0.07 |

M-9 | 2017 | 1020 | 64 | 2.22 ± 0.09 | 1.75 ± 0.07 | 1.27 ± 0.05 |

The number of droplets measured is . The equilibration time gives the time for which the charge was held at the final temperature after having been cooled from a starting temperature of 1100°C at 1°C/hr. Charges SI-5, SI-8 and SI-13 have a bulk composition identical to that of a dyke cutting the Sept Iles intrusion and charges M-5, M-6 and M-9 have a bulk composition typical of an intermediate basalt from the Mull Tertiary volcano (for further details see Charlier and Grove (2012)).

We find values of between 1.06 and 1.8, implying a size distribution in all cases except one (and including the estimated uncertainty) significantly more broad than either of these two mechanisms would predict. We interpret this discrepancy as evidence for continuous nucleation while the existing droplets are ripening, consistent with the design of the experiments in which the temperature was reduced at a steady rate into the binodal.

## Conclusions

Moment-based methods for particle size characterization provide a simple way to describe a population of (sub-) spherical particles (crystals, sedimentary clasts, emulsion droplets or bubbles), and have the advantage that exact results allow the different averages of the three dimensional population to be deduced (with estimated error bars) from two dimensional sections. Which of the mean diameters to use depends on the phenomena of interest, but we suggest that and give a good first characterization of the population (including the spread of diameters), while can be useful for accumulations of sedimenting grains or rising bubbles. Their ratio is a measure of the width of the size distribution, and this can carry information about growth mechanisms of inclusions.

## Acknowledgements

We thank Fergus Gibb and Michael Henderson for their help to access material from the Shiant Isles Main Sill. We acknowledge the loan of Shiant samples from the British Geological Survey and thank Michael Togher for his helpfulness and efficiency at facilitating the loans. We are grateful to Bernard Charlier for providing images of his experimental charges, and for the loan of the charges themselves to enable us to create further images. Comments from two anonymous reviewers greatly improved the manuscript. V.C.H. is supported by a Natural Environment Research Council studentship. M.B.H. acknowledges support from the Natural Environment Research Council [grant number NE/J021520/1].

## References

_{2}O

_{3}and the effect of composition, temperature, oxygen fugacity and pressure on their redox states

_{2}O-K

_{2}O-MgO-FeO-Fe

_{2}O

_{3}-Al

_{2}O

_{3}-TiO

_{2}-SiO

_{2}liquids: New measurements and derived partial molar properties

_{2}O-FeO-Fe

_{2}O

_{3}-Al

_{2}O

_{3}-SiO

_{2}and implications for the kinetics of silicate melt unmixing

### Appendix A: Solution for the moments

In what follows, we shall use upper case symbols to refer to three-dimensional quantities, and lower case to refer to two-dimensional quantities. Suppose we have a random distribution of spheres in space, where there are *N* spheres per unit volume. Further, suppose that the fraction of the number of spheres that have diameters in the range *D* to is , so that the probability density function is normalized: . We then imagine passing a plane section through this distribution, which produces an infinite collection of circular cross-sections. Let there be *n* such cross-sections per unit area of the plane, and let the fraction of the number of circles which have diameters between *d* and be , so that is also normalized: .

A sphere will only intersect the plane if it lies in a volume close to the plane, in particular if the perpendicular distance *y* of the sphere centre from the plane is such that . Thus the number of circles per unit area will be .

*y*is distributed uniformly, the probability density of circle diameters will be proportional to . Taking into account the normalization and , we obtain the distribution of the number of circles per unit area with different diameters for the monodisperse spheres case: Equation A6 then generalizes immediately to the case where the spheres have a distribution of sizes, as the numbers of circles per unit area add up in a linear manner: We can now simply write down the moments of the two dimensional distribution according to equation A2 with the only ‘trick’ required being to reverse the order of integration (see Fig. 8) between lines 2 and 3 below: The relations below follow immediately: and from these the relations between the different sphere and circle average diameters [equations 9 to 12] follow by means of equations A3 and A4.

### Appendix B: Droplet size distributions for different growth mechanisms

*u*is the ratio of the droplet radius to this critical radius, the normalized number density of droplets is given by where and zero otherwise (Meinders

*et al.*, 2001). This distribution can be integrated numerically to obtain moments of the distribution, and specifically we find Alternatively, droplets can grow through coalescence as they diffuse under Brownian motion. If we accept Smoluchowski's approximation to Brownian aggregation (Swift and Friedlander, 1964) then, starting with monodisperse droplets of number density at time , there will be a distribution of sizes at a later time

*t*. Specifically, let be the number of droplets per unit volume that are

*k*times the volume of the initial droplets, then Smoluchowski's approximation gives where and is Boltzmann's constant,

*T*the temperature and η the viscosity of the continuous phase. In the limit of large time and large

*k*, this reduces to Now, if the initial droplets have a volume then the diameter of a droplet formed from

*k*of these is . Thus the number of droplets per unit volume with diameters in the range

*D*to

*D*+ δ

*D*is where From this we find