Weather stations can represent local weather variability and extremes more reliably than gridded products and are therefore better suited for local climate impact applications like calculation of the Fire Weather Index (FWI), a multivariate index for wildfire danger assessment. However, the prediction at multiple sites poses the challenge of preserving spatial consistency across locations, requiring a suitable multi-site approach. This study evaluates the potential of Convolutional Neural Networks (CNNs) for statistical downscaling (SD) of FWI predictions across the Iberian Peninsula. We compare our CNN-Multi-Gaussian (CNN-MG) model against Generalized Linear Models (GLMs) and a benchmark single-site CNN approach. Our evaluation focuses on predictive accuracy, distributional congruence, spatial coherence and extreme events reproducibility using daily FWI data from 29 locations in Spain. The CNN-MG model, which integrates the covariance structure of the predictands, outperformed other methods in representing FWI distributions across both single and multisite scales. Moreover, our model provides greater physical interpretability via eXplainable Artificial Intelligence (XAI) techniques, while also emphasizing simplicity and ease of training.

Neural Weather Models (NWM) are novel data-driven weather forecasting tools based on neural networks that have recently achieved comparable deterministic forecast skill to current operational approaches using significantly less real-time computational resources. The short inference times required by NWMs allow the generation of a large ensemble potentially providing benefits in quantifying the forecast uncertainty, particularly for extreme events, which is of critical importance for various socio-economic sectors. Here we propose a novel ensemble design for NWMs spanning two main sources of uncertainty: epistemic —or model uncertainty,— and aleatoric —or initial condition uncertainty. For the epistemic uncertainty, we propose an effective strategy for creating a diverse ensemble of NWMs that captures uncertainty in key model parameters. For the aleatoric, we explore the “breeding of growing modes” for the first time on NWMs, a technique traditionally used for operational numerical weather predictions as an estimate of the initial condition uncertainty. The combination of these two types of uncertainty produces an ensemble of NWM-based forecasts that is shown to improve upon benchmark probabilistic NWM and is competitive with the 51-member ensemble of the European Centre for Medium-Range Weather Forecasts based on the Integrated Forecasting System (IFS) —a gold standard in weather forecasting,— in terms of both error and calibration. In addition, we report better probabilistic skill than the IFS over land for two key variables: surface wind and air surface temperature.