Abstract. The South American Summer Monsoon (SASM) is the main driver of regional hydroclimate variability across tropical and subtropical South America. It is best recorded on paleoclimatic timescales by stable oxygen isotope proxies, which are more spatially representative of regional hydroclimate than proxies for local precipitation alone. Network studies of proxies that can isolate regional influences lend particular insight into various environmental characteristics that modulate hydroclimate, such as atmospheric circulation variability and changes in the regional energy budget as well as understanding the climate system sensitivity to external forcings. We extract the coherent modes of variability of the SASM over the last millennium (LM) using a Monte Carlo empirical orthogonal function (MCEOF) decomposition of 14 δ18O proxy records and compare them with modes decomposed from isotope-enabled climate model data. The two leading modes reflect the isotopic variability associated with (1) thermodynamic changes driving the upper-tropospheric monsoon circulation (Bolivian High–Nordeste Low waveguide) and (2) the latitudinal displacement of the South Atlantic Convergence Zone (SACZ). The spatial characteristics of these modes appear to be robust features of the LM hydroclimate over South America and are reproduced both in the proxy data and in isotope-enabled climate models, regardless of the nature of the model-imposed external forcing. The proxy data document that the SASM was characterized by considerable temporal variability throughout the LM, with significant departures from the mean state during both the Medieval Climate Anomaly (MCA) and the Little Ice Age (LIA). Model analyses during these periods suggest that the local isotopic composition of precipitation is primarily a reflection of upstream rainout processes associated with monsoon convection. Model and proxy data both point to an intensification of the monsoon during the LIA over the central and western parts of tropical South America and indicate a displacement of the South Atlantic Convergence Zone (SACZ) to the southwest. These centennial-scale changes in monsoon intensity over the LM are underestimated in climate models, complicating the attribution of changes on these timescales to specific forcings and pointing toward areas of important model development.