Summary Statistical functionals are called elicitable if there exists a loss or scoring function under which the functional is the optimal point forecast in expectation. While the mean and quantiles are elicitable, it has been shown in Heinrich (2014) that the mode cannot be elicited if the true distribution can follow any Lebesgue density. We strengthen this result substantially, showing that the mode cannot be elicited if the true distribution can be any strongly unimodal distribution with continuous Lebesgue density, i.e., a continuous density with only one local maximum. Likewise, the mode fails to be identifiable relative to this class.

The Greater Horn of Africa (GHA) is highly vulnerable to climate and weather hazards such as drought, heat waves, and floods. There is a need for accurate seasonal forecasts to prepare for risks (such as crop failure and reduced grazing opportunities) and take advantage of favorable conditions (rains arrive on time and where they are needed) when they arise. As such, information at finer spatial scales than current state‐of‐the‐art global prediction models can provide is needed. Dynamical downscaling is one method employed to obtain information at finer scales. However, providers of seasonal forecasts over the GHA are hampered by limited computational resources and time constraints that restrict the number of global model ensemble members that can be downscaled. Some ensemble subselection criteria must be employed. Currently, providers take an uninformed (or random) approach. Specifically, forecasters simply take the first ensemble member of the global model seasonal forecast ensemble. However, recent work, focused on decadal prediction, has shown that subselecting global model ensemble members in an informed way, that is, according to their ability to reproduce key features of the climate system, results in improved predictions. This emerges from the fact that the climate system is likely more predictable than our models would have us believe. Seeing an opportunity for improvement, we apply the same thinking to the seasonal context and assess several procedures for subselecting ensemble members from seasonal predictions with exchangeable members. Such informed subselections have the potential to take advantage of information in an ensemble of global simulations that might be missed by random selection. Three subselection methods are investigated, with a focus on seasonal predictions for rainfall over GHA. We demonstrate that informed subselection leads to systematically higher skill than random subselection. We find that (1) for small subsample sizes, such as would be chosen for dynamical downscaling and/or downstream impact modeling, informed subselection nearly always outperforms random subselection, (2) subselecting based on well‐known teleconnections benefits those seasons in which such pathways are active, such as OND and JJAS, and (3) ‐means subselection outperforms random selection for small ensemble sizes throughout all seasons, including the notoriously difficult to predict MAM season. These techniques require only input that is available at the time of the forecast release and are easy to apply operationally.

Agricultural food production and natural ecological systems depend on a range of seasonal climate indicators that describe seasonal patterns in climatological conditions. This paper proposes a probabilistic forecasting framework for predicting the end of the freeze-free season, or the time to a mean daily near-surface air temperature below 0 • C (here referred to as hard freeze). The forecasting framework is based on the multi-model seasonal forecast ensemble provided by the Copernicus Climate Data Store and uses techniques from survival analysis for time-to-event data. The original mean daily temperature forecasts are statistically post-processed with a mean and variance correction of each model system before the time-to-event forecast is constructed. In a case study for a region in Fennoscandia covering Norway for the period 1993-2020, the proposed forecasts are found to outperform a climatology forecast from an observation-based data product at locations where the average predicted time to hard freeze is less than 40 days after the initialization date of the forecast on October 1.

We introduce a class of proper scoring rules for evaluating spatial point process forecasts based on summary statistics. These scoring rules rely on Monte-Carlo approximations of expectations and can therefore easily be evaluated for any point process model that can be simulated. In this regard, they are more flexible than the commonly used logarithmic score and other existing proper scores for point process predictions. The scoring rules allow for evaluating the calibration of a model to specific aspects of a point process, such as its spatial distribution or tendency towards clustering. Using simulations we analyze the sensitivity of our scoring rules to different aspects of the forecasts and compare it to the logarithmic score. Applications to earthquake occurrences in northern California, USA and the spatial distribution of Pacific silver firs in Findley Lake Reserve in Washington, USA highlight the usefulness of our scores for scientific model selection.

<p>In the agricultural sector there is a high interest for forecasts that predict relevant agroclimatic indicators related to heat accumulation and frost characteristics. The forecasts can simplify agricultural decisions related to planting and harvest timing. Motivated by this, we propose a probabilistic forecasting framework for predicting the end of the freeze-free season, or the time to a mean daily near-surface air temperature below 0 <sup>°</sup>C (here referred to as hard freeze). The forecasts are constructed based on a multi-model seasonal temperature forecast ensemble provided by the Copernicus Climate Data Store. The raw temperature forecast is statistically post-processed through a mean and variance correction. The resulting ensemble is next used as input to a survival analysis model. Survival analysis is a broad statistical field that is commonly used in the field of biostatistics, but rarely used in meteorology.</p><p>The forecasting framework is evaluated by predicting the time to hard freeze from October 1 for 1993-2020 for a region in Fennoscandia that covers Norway and parts of Sweden, Finland and Russia. We find that the proposed forecast outperforms a climatology forecast from an observation-based data product at locations where the average predicted time to hard freeze is less than 40 days after the initialization date.</p><p>Our work also forms an entry point showing how survival models can be used in general to construct seasonal forecasts for other meteorological events, e.g. the onset of the rainy season or the time to the next drought.</p>