To accurately predict weather and climate numerical weather prediction (NWP) models are being used (Bauer et al., 2015). Recently, machine learning methods have received more attention as an alternative approach for weather and climate prediction. For example, in Bihlo (2021), Bauer (2018), andWeyn et al. (2021) neural networks are trained on reanalysis data to produce purely data-driven weather forecasts. While still in its infancy, if successful these data-driven approaches would enable issuing weather forecasts at several orders of magnitude faster than conventional NWP models (Pathak et al., 2022).Owing to the inherent uncertainty in the atmospheric dynamical system (Lorenz, 1963) it has been recognized early on in the development of NWP that a measure of uncertainty of a numerical weather forecast can substantially enhance the value of these forecasts. This gave rise to the field of ensemble weather prediction (Leutbecher & Palmer, 2008), which aims to quantify the various sources of uncertainty in NWP models, chief of which are uncertainty in the initial conditions and errors in the numerical model formulation. To overcome these uncertainties, in addition to the single deterministic weather forecast, an ensemble of perturbed forecasts is generated whose overall divergence, or spread, ideally will provide a measure of the uncertainty in the deterministic prediction. In some applications one is interested in the potential worst cases scenarios . Here, the ensemble members themselves are needed as the spread only would not provide the spatial correlation between the different points. The main limiting factor in generating the ensemble is still computational in nature as each ensemble run takes up computational resources thereby limiting the total number of such ensembles, typically less then 100, that can be computed on an operational basis.