The GSTools package provides a Python-based platform for geoostatistical applications. Salient features of GSTools are its random field generation, its kriging capabilities and its versatile covariance model. It is furthermore integrated with other Python packages, like PyKrige, ogs5py or scikit-gstat, and provides interfaces to meshio and PyVista. Four presented workflows showcase the abilities of GSTools.