The rise and decline of many complex, pre-European maize-farming cultures in the American Southwest coincides with the warm, climatically quiescent Medieval Climate Anomaly (MCA, ca. 850-1350 CE) and transition to the cool, hydrologically variable Little Ice Age (LIA, ca. 1350-1850 CE). The effects of drought on early subsistence agriculture in the Southwest is well studied, but the impact of temperature stress and variability on the growth of maize crops and which areas were most resilient to such stress remain open questions. We statistically downscaled outputs from a paleo-climate model experiment (CESM1 LME) to map changes to cumulative growing degree days for maize (GDD, 30/10°C) and precipitation over Utah between 850 and 1449 CE, and downscaled GDD changes to local Fremont Culture archaeological site occupations from radiocarbon-dated contexts mapped as spatially discrete kernel density estimates of summed probability distributions (SPDs). We then analyzed correspondences between Fremont SPDs and GDD/precipitation between 850 and 1449 CE. In general, we found (1) high Fremont occupation intensity coincident with GDD that is less volatile than the long term average, and low occupation intensity coincident with, or following, periods of volatile GDD; (2) intensified occupation of high-elevation sites during the MCA, followed by a retreat to lower elevation sites coincident with a sudden rise in annual temperature volatility and increasing drought conditions; and (3) these occupation changes occurred in spite of the greater temperatures and variability in GDD at low-elevation sites. We found evidence that increased inter-annual variability of growing seasons prior to the onset of the LIA, was likely a determinant of Fremont subsistence strategy decision making, and high-elevation site occupation. The most resilient Fremont occupations in the face of these challenges were sited where growing season lengths were least variable.

The American Southwest has experienced a series of severe droughts interspersed with strong wet episodes over the past decades, prompting questions about future climate patterns and potential intensification of weather disruptions under warming conditions. Here we show that interannual hydroclimatic variability in this region has displayed a significant level of non-stationarity over the past millennium. Our tree ring-based analysis of past drought indicates that the Little Ice Age (LIA) experienced high interannual hydroclimatic variability, similar to projections for the 21st century. This is contrary to the Medieval Climate Anomaly (MCA), which had reduced variability and therefore may be misleading as an analog for 21st century warming, notwithstanding its warm (and arid) conditions. Given past non-stationarity, and particularly erratic LIA, a 'warm LIA' climate scenario for the coming century that combines high precipitation variability (similar to LIA conditions) with warm and dry conditions (similar to MCA conditions) represents a plausible situation that is supported by recent climate simulations. Our comparison of tree ring-based drought analysis and records from the tropical Pacific Ocean suggests that changing variability in El Niño Southern Oscillation (ENSO) explains much of the contrasting variances between the MCA and LIA conditions across the American Southwest. Greater ENSO variability for the 21st century could be induced by a decrease in meridional sea surface temperature gradient caused by increased greenhouse gas concentration, as shown by several recent climate modeling experiments. Overall, these results coupled with the paleo-record suggests that using the erratic LIA conditions as benchmarks for past hydroclimatic variability can be useful for developing future water-resource management and drought and flood hazard mitigation strategies in the Southwest.