Boolean delay equations (BDEs) are equations with discrete variables evolving in continuous time. They serve as exploratory tools in the study of nonlinear and complex systems. In this review paper, we outline their formulation and illustrate their properties by application to a solid-earth problem and a climate one. The first problem is the seismotectonic description and prediction of earthquakes and of their clustering. The second one is that of the coupled atmosphere-ocean phenomenon of the El Niño-Southern Oscillation (ENSO). Both involve irregular behavior that is hard to predict, although some form of cyclicity is present in both, especially in ENSO. The paper concludes with broad perspectives on the further use of BDEs in the geosciences and elsewhere.