Accelerating mountain glacier recession in a warming climate threatens the sustainability of mountain water resources. The extent to which groundwater will provide resilience to these water resources is unknown, in part due to a lack of data and poorly understood interactions between groundwater and surface water. Here we address this knowledge gap by linking climate, glaciers, surface water, and groundwater into an integrated model of the Shullcas Watershed, Peru, in the tropical Andes, the region experiencing the most rapid mountain‐glacier retreat on Earth. For a range of climate scenarios, our model projects that glaciers will disappear by 2100. The loss of glacial meltwater will be buffered by relatively consistent groundwater discharge, which only receives minor recharge (~2%) from glacier melt. However, increasing temperature and associated evapotranspiration, alongside potential decreases in precipitation, will decrease groundwater recharge and streamflow, particularly for the RCP 8.5 emission scenario.

The relationship between monthly mean sea surface temperature (SST) anomalies in the commonly used El Niño regions and precipitation for 44 stations in Perú is documented for 1950–2002. Linear lag correlation analysis is employed to establish the potential for statistical precipitation forecasts from SSTs. Useful monthly mean precipitation anomaly forecasts are possible for several locations and calendar months if SST anomalies in El Niño 1+2, Niño 3.4, and Niño 4 regions are available. Prediction of SST anomalies in El Niño regions is routinely available from Climate Prediction Center, NOAA, with reasonable skill in the El Niño 3.4 region, but the prediction in El Niño 1+2 region is less reliable. The feasibility of using predicted SST anomalies in the El Niño 3.4 region to predict SST anomalies in El Niño 1+2 region is discussed.

Extreme precipitation events in the Peruvian Andes have significant socioeconomic impacts, yet their atmospheric dynamics are poorly understood. Here austral summer (December–March) wet and dry spells and their continental- and large-scale teleconnections are analyzed using reanalysis, gridded, and in situ precipitation data. Dry and wet spells in the Peruvian Andes show a pervasive dipole pattern with precipitation anomalies of the opposite sign over northeastern Brazil. Composite anomalies of various atmospheric fields during extreme precipitation events indicate that this dipole is related to large-scale adjustments in the upper-tropospheric Bolivian high–Nordeste low system, which in turn are modulated by northward-propagating extratropical Rossby wave trains. At upper- and midtropospheric levels, westerly wind anomalies over the Peruvian Andes suppress moisture flux from the Amazon during dry events, while wet events are characterized by opposite conditions. Yet, while easterly wind anomalies appear to be a prerequisite for heavy precipitation events in the region, they are not a sufficient forcing, as dry days can still occur during such periods. Dry spells in the Peruvian Andes appear to be linked to weakened convective activity over the western tropical Pacific, consistent with the previously documented El Niño influence over the region. Extreme dry and wet spells in northeastern Brazil only show a weak link to precipitation anomalies of the opposite sign over Peru but are strongly coupled with changes in the position and strength of the Nordeste low and the South Atlantic convergence zone.