  An important problem in modern applied science is characterizing the behavior
of systems with complex internal dynamics subjected to external forcings from
their environment. While a great variety of techniques has been developed to
analyze such non-autonomous systems, many approaches rely on the availability
of ensembles of experiments or simulations in order to generate sufficient
information to encapsulate the external forcings. This makes them unsuitable to
study important classes of natural systems such as climate dynamics where only
a single realization is observed. Here, we show that operator-theoretic
techniques previously developed to identify slowly decaying observables of
autonomous dynamical systems provide a powerful means for identifying trends
and persistent cycles of non-autonomous systems using data from a \emph{single}
trajectory of the system. Using systematic mathematical analysis and prototype
examples, we demonstrate that eigenfunctions of Koopman and transfer operators
provide coordinates that simultaneously capture nonlinear trends and coherent
modes of internal variability. In addition, we apply our framework to two
real-world examples from present and past climate dynamics: Variability of sea
surface temperature (SST) over the industrial era and the mid-Pleistocene
transition (MPT) of Quaternary glaciation cycles. Our results provide a
nonparametric representation of SST and surface air temperature (SAT) trends
over the industrial era, while also capturing the response of the seasonal
precipitation cycle to these trends. In addition, our paleo-climate analysis
reveals the dominant glaciation cycles over the past 3 million years and the
MPT with a high level of granularity.
