Scaling relationships characterize complex systems at criticality. In the brain, these relationships are evident in scale-invariant activity cascades, so-called neuronal avalanches, quantified by power laws in avalanche size and duration. At the cellular level, neuronal avalanches are identified in spatially distributed groups of neurons that participate in cascades of coincident action potential firing. Such spatiotemporal synchronization is central to theories on brain function, yet scaling relationships in avalanche synchronization have been challenging to study when only a fraction of neurons is observed, underestimating avalanche properties. Here, we study these biases from fractional sampling in an all-to-all, balanced network of excitatory and inhibitory neurons with critical branching process dynamics. We focus on the growth of mean avalanche size with avalanche duration. For parabolic avalanches, this growth is quadratic, quantified by the scaling exponent, , which signifies rapid spatial expansion of coincident firing within a relatively short period of time. In contrast, for fractionally sampled networks. We show that temporal coarse-graining combined with a threshold for the minimally required coincident firing in the network recovers , even when sampling as few as 0.1% of the neurons. In contrast, a commonly proposed 'crackling noise' approach fails to recover under those conditions. Our approach robustly identifies for ongoing neuronal activity in frontal cortex of awake mice using cellular 2-photon imaging. Our findings demonstrate how to correct scaling bias from fractional sampling and identifies rapid, scale-invariant synchronization of cell assemblies in the brain.
Keywords: 2-photon imaging; E/I balanced neural network; Neuronal avalanches; criticality; frontal cortex; mouse; ongoing activity; scaling exponent; subsampling; thresholding.