Abstract. While climate change mitigation targets necessarily concern maximum mean state changes, understanding impacts and developing adaptation strategies will be largely contingent on how climate variability responds to increasing anthropogenic perturbations. Thus far Earth system modeling efforts have primarily focused on projected mean state changes and the sensitivity of specific modes of climate variability, such as the El Niño–Southern Oscillation. However, our knowledge of forced changes in the overall spectrum of climate variability and higher-order statistics is relatively limited. Here we present a new 100-member large ensemble of climate change projections conducted with the Community Earth System Model version 2 over 1850–2100 to examine the sensitivity of internal climate fluctuations to greenhouse warming. Our unprecedented simulations reveal that changes in variability, considered broadly in terms of probability distribution, amplitude, frequency, phasing, and patterns, are ubiquitous and span a wide range of physical and ecosystem variables across many spatial and temporal scales. Greenhouse warming in the model alters variance spectra of Earth system variables that are characterized by non-Gaussian probability distributions, such as rainfall, primary production, or fire occurrence. Our modeling results have important implications for climate adaptation efforts, resource management, seasonal predictions, and assessing potential stressors for terrestrial and marine ecosystems.