Self-organized criticality 1 is one of the key concepts to describe the emergence of complexity in natural systems. The concept asserts that a system self-organizes into a critical state where system observables are distributed according to a power law. Prominent examples of self-organized critical dynamics include piling of granular media 2 , plate tectonics 3 and stick-slip motion 4 . Critical behaviour has been shown to bring about optimal computational capabilities 5 , optimal transmission 6 , storage of information 7 and sensitivity to sensory stimuli 8-10 . In neuronal systems, the existence of critical avalanches was predicted 11 and later observed experimentally 6,12,13 . However, whereas in the experiments generic critical avalanches were found, in the model of ref. 11 they only show up if the set of parameters is fine-tuned externally to a critical transition state. Here, we demonstrate analytically and numerically that by assuming (biologically more realistic) dynamical synapses 14 in a spiking neural network, the neuronal avalanches turn from an exceptional phenomenon into a typical and robust self-organized critical behaviour, if the total resources of neurotransmitter are sufficiently large.In multi-electrode recordings on slices of rat cortex and neuronal cultures 6,12 , neuronal avalanches were observed whose sizes were distributed according to a power law with an exponent of −3/2. The distribution was stable over a long period of time.Variations of the dynamical behaviour are induced by application or wash-out of neuromodulators. Qualitatively identical behaviour can be reached in models such as those in refs 11,15 by variations of a global connectivity parameter. In these models, criticality only shows up if the interactions are fixed precisely at a specific value or connectivity structure.Here, we study a model with activity-dependent depressive synapses and show that existence of several dynamical regimes can be reconciled with parameter-independent criticality. We find that synaptic depression causes the mean synaptic strengths to approach a critical value for a certain range of interaction parameters, whereas outside this range other dynamical behaviours are prevalent, see Fig. 1. We analytically derive an expression for the average coupling strengths among neurons and the average inter-spike intervals in a mean-field approach. The mean-field approximation is applicable here even in the critical state, because the quantities that are averaged do not exhibit power laws, but unimodal distributions. These mean values obey a self-consistency equation that allows us to identify the mechanism that drives the dynamics of the system towards the critical regime. Moreover, the critical regime induced by the synaptic dynamics is robust to parameter changes.Consider a network of N integrate-and-fire neurons. Each neuron is characterized by a membrane potential 0 < h i (t ) < θ. The neurons receive external inputs by a random process ξ τ ∈ {1,. . .,N} that selects a neuron ξ τ (t ) = i at a rate τ and a...