We investigate the scaling properties of the sources of crackling noise in a fully dynamic numerical model of sedimentary rocks subject to uniaxial compression. The model is initiated by filling a cylindrical container with randomly sized spherical particles that are then connected by breakable beams. Loading at a constant strain rate the cohesive elements fail, and the resulting stress transfer produces sudden bursts of correlated failures, directly analogous to the sources of acoustic emissions in real experiments. The source size, energy, and duration can all be quantified for an individual event, and the population can be analyzed for its scaling properties, including the distribution of waiting times between consecutive events. Despite the nonstationary loading, the results are all characterized by power-law distributions over a broad range of scales in agreement with experiments. As failure is approached, temporal correlation of events emerges accompanied by spatial clustering. DOI: 10.1103/PhysRevLett.112.065501 PACS numbers: 61.43.Gt, 46.50.+a, 89.75.Da, 91.60.Ba Understanding the processes that lead to catastrophic failure of porous granular media is an important problem in a wide variety of applications, notably in Earth science and engineering [1][2][3][4][5][6][7][8]. Such failure is often preceded by detectable changes in mechanical properties (stress and strain) and in geophysical signals (elastic wave velocity, electrical conductivity, and acoustic emissions) measured remotely at the sample boundary [9]. In particular, acoustic emissions result from sources of internal damage due to sudden local dislocations in the form of tensile or shear microcracks whose origin time, location, orientation, duration, and magnitude can all be inferred from the radiated wave train [10]. Typically, only a very small proportion of the microcracks revealed by destructive thin sectioning after the test result in detectable acoustic emissions [11]. As a consequence, experimental data provide only limited insight into the complexity of the microscopic processes at work prior to failure, notably the probability distributions of the relevant parameters, their scaling properties, and their population dynamics.Theoretical approaches to the dynamics and statistics of rupture cascades have typically been based on stochastic fracture models comprising lattices of springs [12], beams [13,14], fuses [15,16], or fibers [17][18][19]. However, such lattice models involve a strong simplification of the material microstructure and the inhomogeneous stress field. For example, macroscopic laws of damage for cohesive elements are often implemented at the mesoscopic scale on a regular two-dimensional grid, avoiding the truly threedimensional microstructure of real porous media, and often using power-law rheology as an input. Here, we adopt a discrete element modeling (DEM) approach that relaxes all of these restrictions and allows a realistic investigation of the emergent properties of the dynamics, including the temporal and spatial...