Large-scale modes of variability are present in the climate system. These modes are known to have influences on each other but are usually viewed as linear influences. The nonlinear connections among a set of key climate indices are explored here using tools from information theory, which allow us to characterize the causality between indices. It was found that quadratic nonlinear dependencies between climate indices are present at low frequencies, reflecting the complex nature of their dynamics.