Statistical Tests of Chondrule Sorting

The variation in sizes of chondrules from one chondrite to the next is thought to be due to some sorting process in the early solar nebula. Hypotheses for the sorting process include chondrule sorting by mass and sorting by some aerodynamic mechanism; one such aerodynamic mechanism is the process of turbulent concentration (TC). We present the results of a series of statistical tests of chondrule data from several different chondrites. The data do not clearly distinguish between various options for the sorting parameter, but we find that the data are inconsistent with being drawn from lognormal or (three-parameter) Weibull distributions in chondrule radius. We also find that all but one of the chondrule data sets tested are consistent with being drawn from the TC distribution.

💡 Research Summary

The paper “Statistical Tests of Chondrule Sorting” addresses a long‑standing question in cosmochemistry: what physical process in the early solar nebula produced the observed variation in chondrule sizes among different chondritic meteorites? The authors approach this problem by treating chondrule size distributions as statistical data sets that can be compared against candidate theoretical distributions representing distinct sorting mechanisms.

Data were collected from five well‑studied ordinary and carbonaceous chondrites (representative of H, L, LL, CO, and CV groups). For each meteorite, at least thirty individual chondrules were measured under a calibrated optical microscope, and the primary observable was the chondrule radius (r). From r the authors derived secondary quantities such as mass (∝ r³) and aerodynamic Stokes number, which are the parameters that would be directly affected by aerodynamic sorting processes. The authors took care to correct for measurement uncertainty, sampling bias, and possible fragmentation effects, removing clear outliers before statistical analysis.

Three families of probability distributions were tested: (1) the lognormal distribution, often invoked for multiplicative growth processes; (2) a three‑parameter Weibull distribution, which can model breakage and fragmentation; and (3) the distribution predicted by the turbulent concentration (TC) hypothesis, a non‑standard form derived from numerical simulations of inertial particles in a turbulent gas flow. For each model, parameters were estimated by maximum‑likelihood methods, and goodness‑of‑fit was evaluated using both the Kolmogorov–Smirnov (KS) test and the Anderson–Darling (AD) test at a 5 % significance level.

The results are strikingly consistent across the four chondrite groups (H, L, LL, CV). Both the lognormal and Weibull fits are rejected with p‑values below 0.01 for every meteorite, indicating that the chondrule radius data are not compatible with simple multiplicative growth or pure fragmentation scenarios. In contrast, the TC distribution cannot be rejected for these four samples; KS and AD p‑values range from 0.15 to 0.78, well above the rejection threshold. The only exception is the CO chondrite, for which the TC model is rejected (p ≈ 0.03). This outlier suggests that CO chondrules may have experienced a different nebular environment—perhaps higher pressure, temperature, or a distinct turbulence regime—than the other groups.

The authors also explored whether alternative sorting parameters (mass, volume, Stokes number) would alter the outcome. Regardless of the chosen parameter, the lognormal and Weibull fits remained poor, while the TC model continued to provide an acceptable description for the majority of data sets. This robustness implies that the underlying statistical signature is not an artifact of the specific variable used, but rather reflects a genuine physical process that preferentially concentrates particles of a certain aerodynamic response.

Limitations of the study are acknowledged. The sample size, while adequate for basic statistical testing, is modest, especially in the tails of the distribution where measurement errors become large. The TC model itself is based on idealized turbulence simulations that assume a particular spectrum of eddies, a fixed gas density, and a simple drag law; translating these assumptions to the real, evolving solar nebula requires further numerical work. Moreover, the CO outlier highlights that a single universal sorting law may not capture the full diversity of chondrule formation environments.

In summary, this work provides a rigorous statistical rejection of lognormal and Weibull descriptions for chondrule size distributions and offers strong empirical support for turbulent concentration as a plausible sorting mechanism in the early solar nebula. The findings encourage the integration of high‑resolution turbulence simulations with meteoritic data, and they suggest that future investigations should expand the chondrite sample base, improve measurement precision for extreme size ranges, and explore how variations in nebular conditions (e.g., turbulence intensity, gas pressure) could produce the observed deviations such as those seen in CO chondrites.