Many soil processes and functions, including carbon storage, water and element cycling, and crop productivity,
are directly or indirectly affected by soil aggregates. As a key property, aggregate tensile strength is usually
quantified by the Weibull distribution. This distribution was found to depend on soil texture, organic matter
content, and aggregate shape; however, effects of cementing agents (i.e., clay mineral and organic compounds)
on aggregate tensile strength cannot be empirically studied independent of shape. The objective was to apply
Discrete Element Modeling (DEM) to simulate the rupture of soil aggregates of different shapes but the same
assumed distributions of cementing agents. The hypothesis was that the Weibull modulus of aggregate tensile
strength could be used to describe the heterogeneity of the particle binding forces within the aggregate. Multiple
DEM simulations of aggregate rupture were carried out for sets of randomly generated binding force distributions
between particles of the same diameter. Aggregate properties and shapes were obtained from those observed in
arable Gleyic Cambisol topsoil from plots with and without vegetation of a soil compaction trial; additional
aggregate shapes were simulated to increase the predictive range. The simulation results showed that (i) the
sensitivity of Weibull modulus of aggregate tensile strength to the heterogeneity of particle binding force distribution
increased with aggregate sphericity, (ii) the value of the Weibull modulus decreased with increasing
heterogeneity of this distribution, and (iii) internal heterogeneity also reduced the crack surface area produced
by aggregate rupture. Results corresponded to friability theory implying that small values of Weibull modulus
indicated more complex soil structural organization. The DEM simulations could help explaining complex interactions
between aggregate shape cementing agents and stability. Our DEM approach also uncovered the
importance of Weibull modulus to assess structural features of soil aggregates and to describe the heterogeneity
of the particle binding forces within the aggregate.