In all panels, blue points and bars denote experimental data analyzed and aggregated at the family level (Caporaso et al., 2011). Red lines and shading denote best fit model predictions as the mean and standard deviation, respectively, across 20 random instances of the best fit ensemble level parameters, . (A–D) Illustrations of various time series statistics in (E–L). (A) The distribution of richness , the number of consumers present at a sampling time, and its mean are well fit by the model. (B) The variance and mean over time of each family’s abundance (abu) scale as a power law with exponent . Here, in experimental data and in simulations. (C) The distribution of log10(abundance change) across all families is well fit by an exponential with standard deviation . The gray line denotes the best fit exponential distribution, and is largely overlapping with the model prediction in red. (D) The distribution of restoring slopes , defined based on the linear regression between the abundance change and the relative abundance for a species across time, is tightly distributed around a mean that reflects the environmental restoring force. Best fit values of model parameters were determined by minimizing errors in , , , and (E–H, respectively). Using these values, our model also reproduced the distribution of prevalences (fraction of sampling times in which a consumer is present, I), the relationship between prevalence and mean abundance (J), the distributions of residence and return times (durations of sustained presence or absence, respectively, as illustrated in D) (K), and the rank distribution of abundances (L).