ALBERT

Description: Monthly tropical sea surface temperature (SST) data are used as predictors to make statistical forecasts of cold season (November–March) precipitation and temperature for the contiguous United States. Through the use of the combined-lead sea surface temperature (CLSST) model, predictive information is discovered not just in recent SSTs but also from SSTs up to 18 months prior. We find that CLSST cold season forecast anomaly correlation skill is higher than that of the North American Multimodel Ensemble (NMME) and the SEAS5 model from the European Centre for Medium-Range Weather Forecasts (ECMWF) when averaged over the United States for both precipitation and 2-m air temperature. The precipitation forecast skill obtained by CLSST in parts of the Intermountain West is of particular interest because of its implications for water resources. In those regions, CLSST dramatically improves the skill over that of the dynamical model ensembles, which can be attributed to a robust statistical response of precipitation in this region to SST anomalies from the previous year in the tropical Pacific.

Description: Forecast skill of numerical weather prediction (NWP) models for precipitation accumulations over California is rather limited at subseasonal time scales, and the low signal-to-noise ratio makes it challenging to extract information that provides reliable probabilistic forecasts. A statistical postprocessing framework is proposed that uses an artificial neural network (ANN) to establish relationships between NWP ensemble forecast and gridded observed 7-day precipitation accumulations, and to model the increase or decrease of the probabilities for different precipitation categories relative to their climatological frequencies. Adding predictors with geographic information and location-specific normalization of forecast information permits the use of a single ANN for the entire forecast domain and thus reduces the risk of overfitting. In addition, a convolutional neural network (CNN) framework is proposed that extends the basic ANN and takes images of large-scale predictors as inputs that inform local increase or decrease of precipitation probabilities relative to climatology. Both methods are demonstrated with ECMWF ensemble reforecasts over California for lead times up to 4 weeks. They compare favorably with a state-of-the-art postprocessing technique developed for medium-range ensemble precipitation forecasts, and their forecast skill relative to climatology is positive everywhere within the domain. The magnitude of skill, however, is low for week-3 and week-4, and suggests that additional sources of predictability need to be explored.

Description: Downscaling precipitation fields is a necessary step in a number of applications, especially in hydrological modeling where the meteorological forcings are frequently available at too coarse resolution. In this article, we review the Gibbs sampling disaggregation model (GSDM), a stochastic downscaling technique originally proposed by Gagnon et al. The method is capable of introducing realistic, weather-dependent, and possibly anisotropic fine-scale details, while preserving the mean rain rate over the coarse-scale pixels. The main developments compared to the former version are (i) an adapted Gibbs sampling algorithm that enforces the downscaled fields to have a similar texture to that of the analysis fields, (ii) an extensive test of various meteorological predictors for controlling specific aspects of the texture such as the anisotropy and the spatial variability, and (iii) a review of the regression equations used in the model for defining the conditional distributions. A perfect-model experiment is conducted over a domain in the southeastern United States. The metrics used for verification are based on the concept of gridded, stratified variogram, which is introduced as an effective way of reproducing the abilities of human eyes for detecting differences in the field texture. Results indicate that the best overall performances are obtained with the most sophisticated, predictor-based GSDM variant. The 600-hPa wind is found to be the best year-round predictor for controlling the anisotropy. For the spatial variability, kinematic predictors such as wind shear are found to be best during the convective periods, while instability indices are more informative elsewhere.

Description: Characteristics of the European Centre for Medium-Range Weather Forecast’s (ECMWF’s) 00 UTC diagnosed 2-meter temperatures (T2m) from 4D-Var and global ensemble forecasts initial conditions were examined in 2018 over the contiguous US at 1/2-degree grid spacing. These were compared against independently generated, upscaled high-resolution T2m analyses that were created with a somewhat novel data assimilation methodology, an extension of classical optimal interpolation (OI) to surface data analysis. The analysis used a high-resolution, spatially detailed climatological background and was statistically unbiased. Differences of the ECMWF 4D-Var T2m initial states from the upscaled OI reference were decomposed into a systematic component and a residual component. The systematic component was determined by applying a temporal smoothing to the time series of differences between the ECMWF T2m analyses and the OI analyses. Systematic errors at 00 UTC were commonly 1 K or more and larger in the mountainous western US, with the ECMWF analyses cooler than the reference. The residual error is regarded as random in character and should be statistically consistent with the spread of the ensemble of initial conditions after inclusion of OI analysis uncertainty. This analysis uncertainty was large in the western US, complicating interpretation. There were some areas suggestive of an over-spread initial ensemble, others under-spread. Assimilation of more observations in the reference OI analysis would reduce analysis uncertainty, facilitating more conclusive determination of initial-condition ensemble spread characteristics.

Description: The spread of an ensemble of weather predictions initialized from an ensemble Kalman filter may grow slowly relative to other methods for initializing ensemble predictions, degrading its skill. Several possible causes of the slow spread growth were evaluated in perfect- and imperfect-model experiments with a two-layer primitive equation spectral model of the atmosphere. The causes examined were the covariance localization, the additive noise used to stabilize the assimilation method and parameterize the system error, and the model error itself. In these experiments, the flow-independent additive noise was the biggest factor in constraining spread growth. Preevolving additive noise perturbations were tested as a way to make the additive noise more flow dependent. This modestly improved the data assimilation and ensemble predictions, both in the two-layer model results and in a brief test of the assimilation of real observations into a global multilevel spectral primitive equation model. More generally, these results suggest that methods for treating model error in ensemble Kalman filters that greatly reduce the flow dependency of the background-error covariances may increase the filter analysis error and decrease the rate of forecast spread growth.

Description: Probabilistic quantitative precipitation forecasts (PQPFs) were generated from The Observing System Research and Predictability Experiment (THORPEX) Interactive Grand Global Ensemble (TIGGE) database from July to October 2010 using data from Europe (ECMWF), the United Kingdom [Met Office (UKMO)], the United States (NCEP), and Canada [Canadian Meteorological Centre (CMC)]. Forecasts of 24-h accumulated precipitation were evaluated at 1° grid spacing within the contiguous United States against analysis data based on gauges and bias-corrected radar data. PQPFs from ECMWF’s ensembles generally had the highest skill of the raw ensemble forecasts, followed by CMC. Those of UKMO and NCEP were less skillful. PQPFs from CMC forecasts were the most reliable but the least sharp, and PQPFs from NCEP and UKMO ensembles were the least reliable but sharper. Multimodel PQPFs were more reliable and skillful than individual ensemble prediction system forecasts. The improvement was larger for heavier precipitation events [e.g., 〉10 mm (24 h)−1] than for smaller events [e.g., 〉1 mm (24 h)−1]. ECMWF ensembles were statistically postprocessed using extended logistic regression and the five-member weekly reforecasts for the June–November period of 2002–09, the period where precipitation analyses were also available. Multimodel ensembles were also postprocessed using logistic regression and the last 30 days of prior forecasts and analyses. The reforecast-calibrated ECMWF PQPFs were much more skillful and reliable for the heavier precipitation events than ECMWF raw forecasts but much less sharp. Raw multimodel PQPFs were generally more skillful than reforecast-calibrated ECMWF PQPFs for the light precipitation events but had about the same skill for the higher-precipitation events; also, they were sharper but somewhat less reliable than ECMWF reforecast-based PQPFs. Postprocessed multimodel PQPFs did not provide as much improvement to the raw multimodel PQPF as the reforecast-based processing did to the ECMWF forecast. The evidence presented here suggests that all operational centers, even ECMWF, would benefit from the open, real-time sharing of precipitation forecast data and the use of reforecasts.

Description: This manuscript presents a detailed multiscale analysis—using observations, model analyses, and ensemble forecasts—of the extreme heat wave over Russia and historic floods over Pakistan during late July 2010, with an emphasis on the floods over northern Pakistan. The results show that recirculation of air and dynamically driven subsidence occurring with the intensification of the blocking anticyclone in late July 2010 were key factors for producing the exceptionally warm temperatures over western Russia. Downstream energy dispersion from the blocking region led to trough deepening northwest of Pakistan and ridge building over the Tibetan Plateau, thereby providing the linkage between the Russian heat wave and Pakistan flood events on the large scale, in agreement with previous studies. The extratropical downstream energy dispersion and enhanced convective outflow on the large scale associated with the active phase of the Madden–Julian oscillation facilitated the formation of an intense upper-level jet northwest of Pakistan. During this same period an intense southeasterly, low-level, barrier jet–like feature formed over northern Pakistan in conjunction with a westward-moving monsoon depression. This low-level jet and deep easterly flow on the equatorward flank of an anomalous anticyclone over the Tibetan Plateau further enhanced the transport of deep tropical moisture into Pakistan and produced a sustained upslope flow and an extended period of active convection, thereby providing an important contribution leading to the exceptional rainfall amounts. The deep easterly flow and intense low-level jet were features that were absent during previous convective episodes over northern Pakistan in 2010, and hence, were likely key factors in the increased severity of the late July event.

Description: NOAA’s second-generation reforecasts are approximately consistent with the operational version of the 2012 NOAA Global Ensemble Forecast System (GEFS). The reforecasts allow verification to be performed across a multidecadal time period using a static model, in contrast to verifications performed using an ever-evolving operational modeling system. This contribution examines three commonly used verification metrics for reforecasts of precipitation over the southeastern United States: equitable threat score, bias, and ranked probability skill score. Analysis of the verification metrics highlights the variation in the ability of the GEFS to predict precipitation across amount, season, forecast lead time, and location. Beyond day 5.5, there is little useful skill in quantitative precipitation forecasts (QPFs) or probabilistic QPFs. For lighter precipitation thresholds [e.g., 5 and 10 mm (24 h)−1], use of the ensemble mean adds about 10% to the forecast skill over the deterministic control. QPFs have increased in accuracy from 1985 to 2013, likely due to improvements in observations. Results of this investigation are a first step toward using the reforecast database to distinguish weather regimes that the GEFS typically predicts well from those regimes that the GEFS typically predicts poorly.

Description: The ability of five operational ensemble forecast systems to accurately represent and predict atmospheric rivers (ARs) is evaluated as a function of lead time out to 10 days over the northeastern Pacific Ocean and west coast of North America. The study employs the recently developed Atmospheric River Detection Tool to compare the distinctive signature of ARs in integrated water vapor (IWV) fields from model forecasts and corresponding satellite-derived observations. The model forecast characteristics evaluated include the prediction of occurrence of ARs, the width of the IWV signature of ARs, their core strength as represented by the IWV content along the AR axis, and the occurrence and location of AR landfall. Analysis of three cool seasons shows that while the overall occurrence of ARs is well forecast out to a 10-day lead, forecasts of landfall occurrence are poorer, and skill degrades with increasing lead time. Average errors in the position of landfall are significant, increasing to over 800 km at 10-day lead time. Also, there is a 1°–2° southward position bias at 7-day lead time. The forecast IWV content along the AR axis possesses a slight moist bias averaged over the entire AR but little bias near landfall. The IWV biases are nearly independent of forecast lead time. Model spatial resolution is a factor in forecast skill and model differences are greatest for forecasts of AR width. This width error is greatest for coarser-resolution models that have positive width biases that increase with forecast lead time.

Description: High-quality, high-resolution, hourly unbiased surface (2 m) temperature analyses are needed for many applications, including training and validation of statistical postprocessing applications. These temperature analyses are often generated through data assimilation procedures, whereby a background short-range gridded forecast is adjusted to newly available observations. Even with frequent updates to newly available observations, surface-temperature analysis errors and biases can be comparatively large relative to errors and biases of midtropospheric variables, especially over land, despite more near-surface in situ observations. Larger near-surface errors may have several causes, including biased background forecasts and the spatial heterogeneity of surface temperatures that results from subgrid-scale surface, vegetation, land-use, and terrain variations. Are biased raw background forecasts the predominant cause of surface temperature analysis errors? Part I of this two-part series describes a simple benchmark for evaluating the error characteristics of short-term (1 h) raw model background surface temperature forecasts. For stations with a relatively complete time series of data, it is possible to generate an hourly, diurnally, and seasonally dependent observation climatology at a station. The deviation of the current hour’s temperature observation with respect to this hour’s and Julian day’s climatology is added to the climatology for the next hour. For contiguous U.S. stations in July 2015, the station benchmark was lower in error than interpolated 1-h high-resolution numerical predictions of surface temperature from NOAA’s High-Resolution Rapid Refresh (HRRR) system, although not including full postprocessing. For August 2018, 1-h HRRR forecasts were much improved when tested against the station benchmark.