Nanoscale and mesoscopic disorder and associated local hysteretic responses underpin the unique properties of spin and cluster glasses, phase-separated oxides, polycrystalline ferroelectrics, and ferromagnets alike. Despite the rich history of the field, the relationship between the statistical descriptors of hysteresis behavior such as Preisach density, and micro and nanostructure has remained elusive. By using polycrystalline ferroelectric capacitors as a model system, we now report quantitative nonlinearity measurements in 0.025-1 μm 3 volumes, approximately 10 6 times smaller than previously possible. We discover that the onset of nonlinear behavior with thickness proceeds through formation and increase of areal density of micron-scale regions with large nonlinear response embedded in a more weakly nonlinear matrix. This observation indicates that large-scale collective domain wall dynamics, as opposed to motion of noninteracting walls, underpins Rayleigh behavior in disordered ferroelectrics. The measurements provide evidence for the existence and extent of the domain avalanches in ferroelectric materials, forcing us to rethink 100-year old paradigms.piezoelectric | ferroelectric | piezoelectric force microscopy | lead zirconate titanate | thin films H ysteretic response to external stimuli is a ubiquitous feature of disordered systems ranging from friction, capillary condensation, metal-insulator transitions, to magnetic, ferroelectric, and ferroelastic materials. Hysteresis in ferroic systems is utilized in many applications, including magnetic and ferroelectric memory technologies (1) and shape-memory and high-toughness materials (2). The traditional objects of study in the field of ferroic materials driven by information technology applications are strong-field hysteresis loops corresponding to complete switching of material between the limiting states. However, equally important are hysteretic responses in weak (subcoercive) fields, which directly underpin the enhanced properties often observed in disordered and multidomain materials (3 and 4). For ferroelectrics, enhanced electromechanical coupling and dielectric responses enabled by domain wall motion strongly influence electromechanical energy conversion in medical transducers, actuators for precise positioning, and piezoelectric and composite magnetostrictive sensors (5). At the same time, local hysteresis is the primary mechanism of energy losses and dissipation. Design of materials with high coupling coefficients and acceptable losses necessitates deciphering fundamental mechanisms behind lowfield hysteresis behavior.In disordered magnetic systems, the weak-field hysteresis has been a field of active study since the seminal 1870 work by Rayleigh (6), who described the universal parabolic law for ferromagnetic minor loops,where MðHÞ is the field-dependent magnetization, χ init is the initial susceptibility, H max is the maximum applied magnetic field, and α is the irreversible Rayleigh constant. The same constant defines the field-dependent suscep...