The paper investigates the applicability of machine learning (ML) to weather prediction by building a reservoir computing-based, low-resolution, global prediction model. The model is designed to take advantage of the massively parallel architecture of a modern supercomputer. The forecast performance of the model is assessed by comparing it to that of daily climatology, persistence, and a numerical (physics-based) model of identical prognostic state variables and resolution. Hourly resolution 20-day forecasts with the model predict realistic values of the atmospheric state variables at all forecast times for the entire globe. The ML model outperforms both climatology and persistence for the first three forecast days in the midlatitudes, but not in the tropics. Compared to the numerical model, the ML model performs best for the state variables most affected by parameterized processes in the numerical model. a computationally efficient hybrid modeling approach. The present paper implements the parallel ML technique of Pathak, Wikner, et al. (2018) to build a model that predicts the weather in the same format as a global numerical model. We train and verify the model on hourly ERA5 reanalysis data from the European Centre for Medium-Range Weather Forecasts (Hersbach et al., 2019).The work presented here can also be considered an attempt to develop a ML model that can predict the evolution of the three-dimensional, multivariate, global atmospheric state. To the best of our knowledge, the only similar prior attempts were those by Scher (2018) and Scher and Messori (2019), but they trained their three-dimensional multivariate ML model on data that were produced by low-resolution numerical model simulations. In addition, Dueben and Bauer (2018) and Weyn et al., 2019Weyn et al., (2020 designed ML models to predict two-dimensional, horizontal fields of select atmospheric state variables. Similar to our verification strategy, they also verified the ML forecasts against reanalysis data. Compared to all of the aforementioned studies, an important new aspect of our work is that we employ reservoir computing (RC) (Jaeger, 2001;Lukoševičius & Jaeger, 2009;Lukoševičius, 2012;Maass et al., 2002) rather than deep learning (e.g., Goodfellow et al., 2016), which is primarily motivated by the significantly lower computer wall clock time required to train an RC-based model. This difference in training efficiency would allow for a larger number of experiments to tune the ML model at higher resolutions.The structure of the paper is as follows. Section 2 describes the ML model, while section 3 presents the results of the forecast experiments, using as benchmarks persistence of the atmospheric state, daily climatology, and numerical forecasts from a physics-based model of identical prognostic state variables and resolution. Section 4 summarizes our conclusions.