AbstractThe present study applies the empirical orthogonal function (EOF) method to investigate the interannual covariations of East Asian–Australian land precipitation (EAALP) during boreal winter based on observational and reanalysis datasets. The first mode of EAALP variations is characterized by opposite-sign anomalies between East Asia (EA) and Australia (AUS). The second mode features an anomaly pattern over EA similar to the first mode, but with a southwest–northeast dipole structure over AUS. El Niño–Southern Oscillation (ENSO) is found to be a primary factor in modulating the interannual variations of land precipitation over EA and western AUS. By comparison, the Indian Ocean subtropical dipole mode (IOSD) plays an important role in the formation of precipitation anomalies over northeastern AUS, mainly through a zonal vertical circulation spanning from the southern Indian Ocean (SIO) to northern AUS. In addition, the ENSO-independent cold sea surface temperature (SST) anomalies in the western North Pacific (WNP) impact the formation of the second mode. Using the atmospheric general circulation model ECHAM5, three 40-yr numerical simulation experiments differing in specified SST forcings verify the impacts of the IOSD and WNP SST anomalies. Further composite analyses indicate that the dominant patterns of EAALP variability are largely determined by the out-of-phase and in-phase combinations of ENSO and IOSD. These results suggest that in addition to ENSO, IOSD should be considered as another crucial factor influencing the EAALP variability during the boreal winter, which has large implications for improved prediction of EAALP land precipitation on the interannual time scale.