Search

Reliable weather forecasts are valuable in a number of applications, such as agriculture, hydropower, and weather-related disease outbreaks. Global weather forecasts are widely used, but detailed evaluation over specific regions is paramount for users and operational centers to enhance the usability of forecasts and improve their accuracy. This study presents evaluation of the Global Forecast System (GFS) medium-range (1–15 day) precipitation forecasts in the nine subbasins of the Nile basin using NASA’s Integrated Multisatellite Retrievals (IMERG) Final Run satellite–gauge merged rainfall observations. The GFS products are available at a temporal resolution of 3–6 h and a spatial resolution of 0.25°, and the version-15 products are available since 12 June 2019. GFS forecasts are evaluated at a temporal scale of 1–15 days, a spatial scale from 0.25° to all the way to the subbasin scale, and for a period of one year (15 June 2019–15 June 2020). The results show that performance of the 1-day lead daily basin-averaged GFS forecast performance, as measured through the modified Kling–Gupta efficiency (KGE), is poor (0 < KGE < 0.5) for most of the subbasins. The factors contributing to the low performance are 1) large overestimation bias in watersheds located in wet climate regimes in the northern hemispheres (Millennium watershed, Upper Atbara and Setit watershed, and Khashm El Gibra watershed), and 2) lower ability in capturing the temporal dynamics of watershed-averaged rainfall that have smaller watershed areas (Roseires at 14 110 km2 and Sennar at 13 895 km2). GFS has better bias for watersheds located in the dry parts of the Northern Hemisphere or wet parts of the Southern Hemisphere, and better ability in capturing the temporal dynamics of watershed-average rainfall for large watershed areas. IMERG Early has better bias than GFS forecast for the Millennium watershed but still comparable and worse bias for the Upper Atbara and Setit and Khashm El Gibra watersheds. The variation in the performance of the IMERG Early could be partly explained by the number of rain gauges used in the reference IMERG Final product, as 16 rain gauges were used for the Millennium watershed but only one rain gauge over each Upper Atbara and Setit and Khashm El Gibra watershed. A simple climatological bias correction of IMERG Early reduces in the bias in IMERG Early over most watersheds, but not all watersheds. We recommend exploring methods to increase the performance of GFS forecasts, including postprocessing techniques through the use of both near-real-time and research-version satellite rainfall products.

A warm bias with a maximum value of over 4 K in the tropical tropopause layer (TTL) is detected in day-5 operational forecasts of the Global/Regional Assimilation and Prediction System (GRAPES) for global medium-range numerical weather prediction (GRAPES_GFS). In this study, the predicted temperature changes caused by different processes are examined, and the predicted cloud fractions are compared with the European Centre for Medium-Range Weather Forecasts ERA5 reanalysis data. It is found that the overprediction of the TTL cirrus fraction contributes to the warm bias due to cloud-radiative heating. The interactions among the ice nucleation, deposition/sublimation, and the large-scale condensation together determine the results of the TTL ice crystal content prediction. Moreover, a range of sensitivity experiments show that the TTL ice crystal content prediction is sensitive to the threshold relative humidity over ice (RHi) in the ice nucleation process. Then the uncertainties of the formulas for saturation vapor pressure over ice at very low temperatures are discussed. The RHi calculated based on the Magnus–Tetens formula is up to 10% higher than that based on the Goff–Gratch formula. As the Goff–Gratch formula is applicable over a broader range of 184–273 K, it is more suitable for the cold TTL. When the Goff–Gratch formula rather than the Magnus–Tetens formula is used in the microphysics scheme, the TTL cirrus forecasts are improved greatly, and the warm bias disappears completely. After investigating the interplay of the dynamical, microphysical, and radiative processes, we find a positive feedback mechanism that exacerbates the TTL cirrus prediction error.

Operational forecasting of tropical cyclone (TC) genesis has improved in recent years but still can be a challenge. Output from global numerical models continues to serve as a primary source of forecast guidance. Bulk verification statistics (e.g., critical success index) of TC genesis forecasts indicate that, overall, global models are increasingly able to predict TC genesis. However, as global model configurations are updated, TC genesis verification statistics will change. This study compares operational and retrospective forecasts from three configurations of NCEP’s Global Forecast System (GFS) to quantify the impact of model upgrades on TC genesis forecasts. First, bulk verification statistics from a homogeneous sample of model initialization cycles during the period 2013–14 are compared. Then, composites of select output fields are analyzed in an attempt to identify any key differences between hit and false alarm events. Bulk statistics indicate that TC genesis forecast performance decreased with the implementation of the 2015 version of the GFS, but then modestly recovered with the 2016 version of the model. In addition, the composite analysis suggests that false alarm forecasts in the 2015 version of the GFS may have been the result of inaccurately forecasting the location and/or strength of upper-level troughs poleward of the TC. There is also evidence of convective feedbacks occurring, such as ridging above the low-level circulation and upper-level convective outflow that were too strong, in this same set of false alarm forecasts. Overall, analyzing retrospective forecasts can assist forecasters in determining the strengths and weaknesses associated with a new configuration of a global model with respect to TC genesis.

Proactive quality control (PQC) is a fully flow dependent QC based on ensemble forecast sensitivity to observations (EFSO). Past studies showed in several independent cases that GFS forecasts can be improved by rejecting observations identified as detrimental by EFSO. However, the impact of cycling PQC in sequential data assimilation has, so far, only been examined using the simple Lorenz ’96 model. Using a low-resolution spectral GFS model that assimilates PrepBUFR (no radiance) observations with the local ensemble transform Kalman filter (LETKF), this study aims to become a bridge between a simple model and the implementation into complex operational systems. We demonstrate the major benefit of cycling PQC in a sequential data assimilation framework through the accumulation of improvements from previous PQC updates. Such accumulated PQC improvement is much larger than the “current” PQC improvement that would be obtained at each analysis cycle using “future” observations. As a result, it is unnecessary to use future information, and hence this allows the operational implementation of cycling PQC. The results show that the analyses and forecasts are improved the most by rejecting all the observations identified as detrimental by EFSO, but that major improvements also come from rejecting just the most detrimental 10% observations. The forecast improvements brought by PQC are observed throughout the 10 days of integration and provide more than a 12-h forecast lead-time gain. An important finding is that PQC not only reduces substantially the root-mean-squared forecast errors but also the forecast biases. We also show a case of “skill dropout,” where the control forecast misses a developing baroclinic instability, whereas the accumulated PQC corrections result in a good prediction.

The FV3GFS is the current operational Global Forecast System (GFS) at the National Centers for Environmental Prediction (NCEP), which combines a finite-volume cubed sphere dynamical core (FV3) and GFS physics. In this study, FV3GFS is used to gain understanding of rapid intensification (RI) of tropical cyclones (TCs) in shear. The analysis demonstrates the importance of TC structure in a complex system like Hurricane Michael, which intensified to a category 5 hurricane over the Gulf of Mexico despite over 20 kt (10 m s−1) of vertical wind shear. Michael’s RI is examined using a global-nest FV3GFS ensemble with the nest at 3-km resolution. The ensemble shows a range of peak intensities from 77 to 159 kt (40–82 m s−1). Precipitation symmetry, vortex tilt, moisture, and other aspects of Michael’s evolution are compared through composites of stronger and weaker members. The 850–200-hPa vertical shear is 22 kt (11 m s−1) in the mean of both strong and weak members during the early stage. Tilt and moisture are two distinguishing factors between strong and weak members. The relationship between vortex tilt and humidification is complex, and other studies have shown both are important for sheared intensification. Here, it is shown that tilt reduction leads to upshear humidification and is thus a driving factor for intensification. A stronger initial vortex and early evolution of the vortex also appear to be the key to members that are able to resist the sheared environment.

The Weather Research and Forecasting (WRF) Model and the Nonhydrostatic Icosahedral Atmospheric Model (NICAM) are forced with the Global Forecast System (GFS) data and run over the United Arab Emirates (UAE) for two 4-day periods: one in the cold season (16–18 December 2017) and another in the warm season (13–15 April 2018). The models’ performance is evaluated against four observational datasets: weather station observations, eddy-covariance flux measurements at Al Ain, microwave radiometer–derived temperature profile, and twice-daily radiosonde measurements at Abu Dhabi. An overestimation of the daily mean air temperature by 1°–3°C is noticed for both models and periods. This warm bias is attributed to the reduced cloud cover and resulting increased surface downward shortwave radiation flux. A comparison with the eddy-covariance data suggested that both models also underestimate the observed albedo. However, when the models predict heavier amounts of precipitation, they tend to be colder than observations, typically by 2°–3°C. NICAM and WRF overpredict the strength of the near-surface wind speed at all weather stations by roughly 1–3 m s−1, which has been attributed to a poor representation of its subgrid-scale fluctuations and surface drag parameterization. WRF tends to be wetter and NICAM drier than the station observations, possibly because of differences in the cloud microphysics schemes. While the performance of both models for the near-surface fields is comparable, NICAM outperforms WRF in the simulation of vertical profiles of temperature, relative humidity, and wind speed, being able to partially correct some of the biases in the GFS data.

In August 2018 and June 2019, NCEP upgraded the operational versions of the High-Resolution Rapid Refresh (HRRR) and Global Forecast System (GFS), respectively. To inform forecasters and model developers about changes in the capabilities and biases of these modeling systems over the western conterminous United States (CONUS), we validate and compare precipitation forecasts produced by the experimental, preoperational HRRRv3 and GFSv15.0 with the then operational HRRRv2 and GFSv14 during the 2017/18 October–March cool season. We also compare the GFSv14 and GFSv15.0 with the operational, high-resolution configuration of the ECMWF Integrated Forecasting System (HRES). We validate using observations from Automated Surface and Weather Observing System (ASOS/AWOS) stations, which are located primarily in the lowlands, and observations from Snowpack Telemetry (SNOTEL) stations, which are located primarily in the uplands. Changes in bias and skill from HRRRv2 to HRRRv3 are small, with HRRRv3 exhibiting slightly higher (but statistically indistinguishable at a 95% confidence level) equitable threat scores. The GFSv14, GFSv15.0, and HRES all exhibit a wet bias at lower elevations and neutral or dry bias at upper elevations, reflecting insufficient terrain representation. GFSv15.0 performance is comparable to GFSv14 at day 1 and superior at day 3, but lags HRES. These results establish a baseline for current operational HRRR and GFS precipitation capabilities over the western CONUS and are consistent with steady or improving NCEP model performance.

Atmospheric rivers, often associated with impactful weather along the west coast of North America, can be a challenge to forecast even on short time scales. This is attributed, at least in part, to the scarcity of eastern Pacific in situ observations. We examine the impact of assimilating dropsonde observations collected during the Atmospheric River (AR) Reconnaissance 2018 field program on the Navy Global Environmental Model (NAVGEM) analyses and forecasts. We compare NAVGEM’s representation of the ARs to the observations, and examine whether the observation–background difference statistics are similar to the observation error variance specified in the data assimilation system. Forecast sensitivity observation impact is determined for each dropsonde variable, and compared to the impacts of the North American radiosonde network. We find that the reconnaissance soundings have significant beneficial impact, with per observation impact more than double that of the North American radiosonde network. Temperature and wind observations have larger total and per observation impact than moisture observations. In our experiment, the 24-h global forecast error reduction from the reconnaissance soundings can be comparable to the reduction from the North American radiosonde network for the field program dates that include at least two flights.

There is a longstanding challenge in numerical weather and climate prediction to accurately model tropical wave variability, including convectively coupled equatorial waves (CCEWs) and the Madden–Julian oscillation. For subseasonal prediction, the European Centre for Medium-Range Weather Forecasts (ECMWF) Integrated Forecasting System (IFS) has been shown to be superior to the NOAA Global Forecast System (GFS) in simulating tropical variability, suggesting that the ECMWF model is better at simulating the interaction between cumulus convection and the large-scale tropical circulation. In this study, we experiment with the cumulus convection scheme of the ECMWF IFS in a research version of the GFS to understand which aspects of the IFS cumulus convection scheme outperform those of the GFS convection scheme in the tropics. We show that the IFS cumulus convection scheme produces significantly different tropical moisture and temperature tendency profiles from those simulated by the GFS convection scheme when it is coupled with other physics schemes in the GFS physics package. We show that a consistent treatment of the interaction between parameterized convective plumes in the GFS planetary boundary layer (PBL) and the IFS convection scheme is required for the GFS to replicate the tropical temperature and moisture profiles simulated by the IFS model. The GFS model with the IFS convection scheme, and the consistent treatment between the convection and PBL schemes, produces much more organized convection in the tropics, and generates tropical waves that propagate more coherently than the GFS in its default configuration due to better simulated interaction between low-level convergence and precipitation.

Since the implementation of all-sky radiance assimilation of the Advanced Microwave Sounding Unit-A (AMSU-A) in the operational hybrid 4D ensemble–variational Global Forecast System at NCEP in 2016, the all-sky approach has been tested to expand to the radiances of Advanced Technology Microwave Sounder (ATMS) in the Gridpoint Statistical Interpolation analysis system (GSI). Following the all-sky framework implemented for the AMSU-A radiances, ATMS radiance assimilation adopts similar procedures in quality control, bias correction, and model of observation error. Efforts have been focused on special considerations that are necessary because of the unique features of the ATMS radiances and water vapor channels, including surface properties based on fields of view size and shape, and taking care of large departures from the first guess (OmF) along coastlines and radiances affected by strong scattering. More importantly, it is shown that this work makes microwave radiance OmFs become more consistent among different sensors, and provides indications of the deficiencies in quality control procedures of the original ATMS and Microwave Humidity Sounder (MHS) clear-sky radiance assimilation. While the generalized tracer effect is noticed, the overall impact on the forecast skill is neutral. This work is included in the upcoming operational implementation in 2019.

Given the network of satellite and aircraft observations around the globe, do additional in situ observations impact analyses within a global forecast system? Despite the dense observational network at many levels in the tropical troposphere, assimilating additional sounding observations taken in the eastern tropical Pacific Ocean during the 2016 El Niño Rapid Response (ENRR) locally improves wind, temperature, and humidity 6-h forecasts using a modern assimilation system. Fields from a 50-km reanalysis that assimilates all available observations, including those taken during the ENRR, are compared with those from an otherwise-identical reanalysis that denies all ENRR observations. These observations reveal a bias in the 200-hPa divergence of the assimilating model during a strong El Niño. While the existing observational network partially corrects this bias, the ENRR observations provide a stronger mean correction in the analysis. Significant improvements in the mean-square fit of the first-guess fields to the assimilated ENRR observations demonstrate that they are valuable within the existing network. The effects of the ENRR observations are pronounced in levels of the troposphere that are sparsely observed, particularly 500–800 hPa. Assimilating ENRR observations has mixed effects on the mean-square difference with nearby non-ENRR observations. Using a similar system but with a higher-resolution forecast model yields comparable results to the lower-resolution system. These findings imply a limited improvement in large-scale forecast variability from additional in situ observations, but significant improvements in local 6-h forecasts.

An important issue in developing a forecast system is its sensitivity to additional observations for improving initial conditions, to the data assimilation (DA) method used, and to improvements in the forecast model. These sensitivities are investigated here for the Global Forecast System (GFS) of the National Centers for Environmental Prediction (NCEP). Four parallel sets of 7-day ensemble forecasts were generated for 100 forecast cases in mid-January to mid-March 2016. The sets differed in their 1) inclusion or exclusion of additional observations collected over the eastern Pacific during the El Niño Rapid Response (ENRR) field campaign, 2) use of a hybrid 4D–EnVar versus a pure EnKF DA method to prepare the initial conditions, and 3) inclusion or exclusion of stochastic parameterizations in the forecast model. The Control forecast set used the ENRR observations, hybrid DA, and stochastic parameterizations. Errors of the ensemble-mean forecasts in this Control set were compared with those in the other sets, with emphasis on the upper-tropospheric geopotential heights and vorticity, midtropospheric vertical velocity, column-integrated precipitable water, near-surface air temperature, and surface precipitation. In general, the forecast errors were found to be only slightly sensitive to the additional ENRR observations, more sensitive to the DA methods, and most sensitive to the inclusion of stochastic parameterizations in the model, which reduced errors globally in all the variables considered except geopotential heights in the tropical upper troposphere. The reduction in precipitation errors, determined with respect to two independent observational datasets, was particularly striking.

The FORMOSAT-7/COSMIC-2 Global Navigation Satellite System (GNSS) Radio Occultation (RO) satellite constellation was launched on June 2019 as a successor of the FORMOSAT-3/COSMIC mission. The Central Weather Bureau (CWB) of Taiwan has received FORMOSAT-7/COSMIC-2 GNSS RO data in real time from Taiwan Analysis Center for COSMIC. With the global numerical prediction system at CWB, a parallel semi-operational experiment assimilating the FORMOSAT-7/COSMIC-2 bending angle data with all other operational observation data has been conducted to evaluate the impact of the FORMOSAT-7/COSMIC-2 data. The first seven-month results show that the quality of the early FORMOSAT-7/COSMIC-2 data has been satisfactory for assimilation. Consistent and significant positive impacts on global forecast skills have been observed since the start of the parallel experiment, with the most significant impact found in the tropical region, reflecting the low-inclination orbital design of the satellites. The impact of the FORMOSAT-7/COSMIC-2 RO data is also estimated using the Ensemble Forecast Sensitivity to Observation Impact (EFSOI) method, showing an average positive impact per observation similar to other existing GNSS RO datasets, while the total impact is impressive by virtue of its large amount. Sensitivity experiments suggest that the quality control processes built in the Gridpoint Statistical Interpolation (GSI) system for RO data work well to achieve a positive impact by the low-level FORMOSAT-7/COSMIC-2 RO data, while more effort on observation error tuning should be focused to obtain an optimal assimilation performance. This study demonstrates the usefulness of the FORMOSAT-7/COSMIC-2 RO data in global numerical weather prediction during the calibration/validation period and leads to the operational use of the data at CWB.

Despite decades of research on the role of moist convective processes in large-scale tropical dynamics, tropical forecast skill in operational models is still deficient when compared to the extratropics, even at short lead times. Here we compare tropical and Northern Hemisphere (NH) forecast skill for quantitative precipitation forecasts (QPFs) in the NCEP Global Forecast System (GFS) and ECMWF Integrated Forecast System (IFS) during January 2015–March 2016. Results reveal that, in general, initial conditions are reasonably well estimated in both forecast systems, as indicated by relatively good skill scores for the 6–24-h forecasts. However, overall, tropical QPF forecasts in both systems are not considered useful by typical metrics much beyond 4 days. To quantify the relationship between QPF and dynamical skill, space–time spectra and coherence of rainfall and divergence fields are calculated. It is shown that while tropical variability is too weak in both models, the IFS is more skillful in propagating tropical waves for longer lead times. In agreement with past studies demonstrating that extratropical skill is partially drawn from the tropics, a comparison of daily skill in the tropics versus NH suggests that in both models NH forecast skill at lead times beyond day 3 is enhanced by tropical skill in the first couple of days. As shown in previous work, this study indicates that the differences in physics used in each system, in particular, how moist convective processes are coupled to the large-scale flow through these parameterizations, appear as a major source of tropical forecast errors.

The National Aeronautics and Space Administration (NASA) RapidScat scatterometer on board the International Space Station (ISS) provides observations of surface winds that can be assimilated into numerical weather prediction (NWP) forecast models. In this study, the authors assess the data quality of the RapidScat Level 2B surface wind vector retrievals and the impact of those observations on the National Oceanic and Atmospheric Administration (NOAA) Global Forecast System (GFS). The RapidScat is found to provide quality measurements of surface wind speed and direction in nonprecipitating conditions and to provide observations that add both information and robustness to the global satellite observing system used in NWP models. The authors find that with an assumed uncertainty in wind speed of around 2 m s−1, the RapidScat has neutral impact on the short-range forecast of surface wind vectors in the tropics but improves both the analysis and background field of surface wind vectors. However, the deployment of RapidScat on the ISS presents some challenges for use of these wind vector observations in operational NWP, including frequent maneuvers of the spacecraft that could alter instrument performance.

Wind energy requires accurate forecasts for adequate integration into the electric grid system. In addition, global atmospheric models are not able to simulate local winds in complex terrain, where wind farms are sometimes placed. For this reason, the use of mesoscale models is vital for estimating wind speed at wind turbine hub height. In this regard, the Weather Research and Forecasting (WRF) Model allows a user to apply different initial and boundary conditions as well as physical parameterizations. In this research, a sensitivity analysis of several physical schemes and initial and boundary conditions was performed for the Alaiz mountain range in the northern Iberian Peninsula, where several wind farms are located. Model performance was evaluated under various atmospheric stabilities and wind speeds. For validation purposes, a mast with anemometers installed at 40, 78, 90, and 118 m above ground level was used. The results indicate that performance of the Global Forecast System analysis and European Centre for Medium-Range Weather Forecasts interim reanalysis (ERA-Interim) as initial and boundary conditions was similar, although each performed better under certain meteorological conditions. With regard to physical schemes, there is no single combination of parameterizations that performs best during all weather conditions. Nevertheless, some combinations have been identified as inefficient, and therefore their use is discouraged. As a result, the validation of an ensemble prediction system composed of the best 12 deterministic simulations shows the most accurate results, obtaining relative errors in wind speed forecasts that are

This study discusses the development of the Hurricane Forecast Improvement Program (HFIP) Corrected Consensus Approach (HCCA) for tropical cyclone track and intensity forecasts. The HCCA technique relies on the forecasts of separate input models for both track and intensity and assigns unequal weighting coefficients based on a set of training forecasts. The HCCA track and intensity forecasts for 2015 were competitive with some of the best-performing operational guidance at the National Hurricane Center (NHC); HCCA was the most skillful model for Atlantic track forecasts through 48 h. Average track input model coefficients for the 2015 forecasts in both the Atlantic and eastern North Pacific basins were largest for the European Centre for Medium-Range Weather Forecasts (ECMWF) deterministic model and the National Centers for Environmental Prediction (NCEP) Global Forecast System (GFS) ensemble mean, but the relative magnitudes of the intensity coefficients were more varied. Input model sensitivity experiments conducted using retrospective HCCA forecasts from 2011 to 2015 indicate that the ECMWF deterministic model had the largest positive impact on the skill of the HCCA track forecasts in both basins. The most important input models for HCCA intensity forecasts are the Hurricane Weather Research and Forecasting (HWRF) Model and the Coupled Ocean–Atmosphere Mesoscale Prediction System-Tropical Cyclone (COAMPS-TC) model initialized from the GFS. Several updates were incorporated into the HCCA formulation prior to the 2016 season. Verification results indicate HCCA continued to be a skillful model, especially for short-range (12–48 h) track forecasts in both basins.

North Atlantic tropical cyclone (TC) forecasts from four ensemble prediction systems (EPSs) are verified using the National Hurricane Center’s (NHC) best tracks for the 2008–15 seasons. The 1–5-day forecasts are evaluated for the 21-member National Centers for Environmental Prediction (NCEP) Global Ensemble Forecast System (GEFS), the 23-member UKMO ensemble (UKMET), and the 51-member European Centre for Medium-Range Weather Forecasts (ECMWF) ensemble, as well as a combination of these ensembles [Multimodel Global (MMG)]. Several deterministic models are also evaluated, such as the Global Forecast System (GFSdet), Hurricane Weather Research and Forecasting Model (HWRF), the deterministic ECMWF model (ECdet), and the Geophysical Fluid Dynamical Laboratory model (GFDL).The ECdet track errors are the smallest on average at all lead times, but are not significantly different from the GEFS and ECMWF ensemble means. All models have a slow bias (90–240 km) in the along-track direction by 120 h, while there is little bias in the cross-track direction. Much of this slow bias is attributed to TCs undergoing extratropical transition (ET). All EPSs are underdispersed in the along-track direction, while the ECMWF is slightly overdispersed in the cross-track direction. The MMG and ECMWF track forecasts have more probabilistic skill than the ECdet and comparable skill to the NHC climatology-based cone forecast. TC intensity errors for the HWRF and GFDL are lower than the coarser models within the first 24 h, but are comparable to the ECdet at longer lead times. The ECMWF and MMG have comparable or better probabilistic intensity forecasts than the ECdet, while the GEFS’s weak bias limits its skill.

This study examines the performance of the NCEP Global Forecast System (GFS) surface layer parameterization scheme for strongly stable conditions over land in which turbulence is weak or even disappears because of high near-surface atmospheric stability. Cases of both deep snowpack and snow-free conditions are investigated. The results show that decoupling and excessive near-surface cooling may appear in the late afternoon and nighttime, manifesting as a severe cold bias of the 2-m surface air temperature that persists for several hours or more. Concurrently, because of negligible downward heat transport from the atmosphere to the land, a warm temperature bias develops at the first model level. The authors test changes to the stable surface layer scheme that include introduction of a stability parameter constraint that prevents the land–atmosphere system from fully decoupling and modification to the roughness-length formulation. GFS sensitivity runs with these two changes demonstrate the ability of the proposed surface layer changes to reduce the excessive near-surface cooling in forecasts of 2-m surface air temperature. The proposed changes prevent both the collapse of turbulence in the stable surface layer over land and the possibility of numerical instability resulting from thermal decoupling between the atmosphere and the surface. The authors also execute and evaluate daily GFS 7-day test forecasts with the proposed changes spanning a one-month period in winter. The assessment reveals that the systematic deficiencies and substantial errors in GFS near-surface 2-m air temperature forecasts are considerably reduced, along with a notable reduction of temperature errors throughout the lower atmosphere and improvement of forecast skill scores for light and medium precipitation amounts.

A new version of the Global Ensemble Forecast System (GEFS, v11) is tested and compared with the operational version (v10) in a 2-yr parallel run. The breeding-based scheme with ensemble transformation and rescaling (ETR) used in the operational GEFS is replaced by the ensemble Kalman filter (EnKF) to generate initial ensemble perturbations. The global medium-range forecast model and the Global Forecast System (GFS) analysis used as the initial conditions are upgraded to the GFS 2015 implementation version. The horizontal resolution of GEFS increases from Eulerian T254 (~52 km) for the first 8 days of the forecast and T190 (~70 km) for the second 8 days to semi-Lagrangian T574 (~34 km) and T382 (~52 km), respectively. The sigma pressure hybrid vertical layers increase from 42 to 64 levels. The verification of geopotential height, temperature, and wind fields at selected levels shows that the new GEFS significantly outperforms the operational GEFS up to days 8–10 except for an increased warm bias over land in the extratropics. It is also found that the parallel system has better reliability in the short-range probability forecasts of precipitation during warm seasons, but no clear improvement in cold seasons. There is a significant degradation of TC track forecasts at days 6–7 during the 2012–14 TC seasons over the Atlantic and eastern Pacific. This degradation is most likely a sampling issue from a low number of TCs during these three TC seasons. The results for an extended verification period (2011–14) and the recent two hurricane seasons (2015 and 2016) are generally positive. The new GEFS became operational at NCEP on 2 December 2015.

The current operational NCEP Global Forecast System (GFS) cumulus convection schemes are updated with a scale-aware parameterization where the cloud mass flux decreases with increasing grid resolution. The ratio of advective time to convective turnover time is also taken into account for the scale-aware parameterization. In addition, the present deep cumulus convection closure using the quasi-equilibrium assumption is no longer used for grid sizes smaller than a threshold value. For the shallow cumulus convection scheme, the cloud-base mass flux is modified to be given by a function of mean updraft velocity. A simple aerosol-aware parameterization where rain conversion in the convective updraft is modified by aerosol number concentration is also included in the update. Along with the scale- and aerosol-aware parameterizations, more changes are made to the schemes. The cloud-base mass-flux computation in the deep convection scheme is modified to use convective turnover time as the convective adjustment time scale. The rain conversion rate is modified to decrease with decreasing air temperature above the freezing level. Convective inhibition in the subcloud layer is used as an additional trigger condition. Convective cloudiness is enhanced by considering suspended cloud condensate in the updraft. The lateral entrainment in the deep convection scheme is also enhanced to more strongly suppress convection in a drier environment. The updated NCEP GFS cumulus convection schemes display significant improvements especially in the summertime continental U.S. precipitation forecasts.

Lake-effect snow (LES) is a cold-season mesoscale convective phenomenon that can lead to significant snowfall rates and accumulations in the Great Lakes region of the United States. While limited-area numerical weather prediction models have shown skill in prediction of warm-season convective storms, forecasting the sharp nature of LES precipitation timing, intensity, and location is difficult because of model error and initial and boundary condition uncertainties. Ensemble forecasting can incorporate and quantify some sources of forecast error, but ensemble design must be considered. This study examines the relative contributions of forecast uncertainties to LES forecast error using a regional convection-allowing data assimilation and ensemble prediction system. Ensembles are developed using various methods of perturbations to simulate a long-lived and high-precipitation LES event in December 2013, and forecast performance is evaluated using observations including those from the Ontario Winter Lake-Effect Systems (OWLeS) campaign. Model lateral boundary conditions corresponding to weather conditions beyond the Great Lakes region play an influential role in LES precipitation forecasts and their uncertainty, as evidenced by ensemble spread, particularly at lead times beyond one day. A strong forecast dependence on regional initial conditions was shown using data assimilation. This sensitivity impacts the timing and intensity of predicted precipitation, as well as band location and orientation assessed with an object-based verification approach, giving insight into the time scales of practical predictability of LES. Overall, an assimilation-cycling convection-allowing ensemble prediction system could improve future lake-effect snow precipitation forecasts and analyses and can help quantify and understand sources of forecast uncertainty.

The National Oceanic and Atmospheric Administration (NOAA) Global Data Assimilation System/Global Forecast System (GDAS/GFS) was extended to assimilate brightness temperatures from the Sondeur Atmosphérique du Profil d’Humidité Intertropicale par Radiométrie (SAPHIR) passive microwave water vapor sounder on board the Megha-Tropiques satellite. Quality control procedures were developed to assess the SAPHIR data quality for assimilating clear-sky observations over ocean surfaces, and to characterize observation biases and errors. A 6-week impact experiment was performed using the GDAS/GFS data assimilation system. The addition of SAPHIR observations on top of the current global observing system improved analysis and forecast humidity root-mean-square error (RMSE) results at the upper levels of the troposphere by about 6%, mostly at 100 hPa, when verified against European Centre for Medium-Range Weather Forecasts (ECMWF) analysis, though some degradation to the forecast humidity was seen at 150–200 hPa. The forecast impacts were predominant at earlier lead times between 24 and 96 h. Verification using global radiosonde observations also showed a reduction of the humidity RMSE from 4% to 6% between 500 hPa and the surface when assimilating SAPHIR, while temperature and wind speed RMSEs were reduced by up to 9% and 7% near the tropical tropopause, respectively. Other conventional forecast skill parameters including the 500-hPa geopotential height anomaly correlation showed neutral impact when assimilating SAPHIR.

Automated aircraft observations of wind and temperature have demonstrated positive impact on numerical weather prediction since the mid-1980s. With the advent of the Water Vapor Sensing System (WVSS-II) humidity sensor, the expanding fleet of commercial aircraft with onboard automated sensors is also capable of delivering high quality moisture observations, providing vertical profiles of moisture as aircraft ascend out of and descend into airports across the continental United States. Observations from the WVSS-II have to date only been monitored within the Global Data Assimilation System (GDAS) without being assimilated. In this study, aircraft moisture observations from the WVSS-II are assimilated into the GDAS, and their impact is assessed in the Global Forecast System (GFS). A two-season study is performed, demonstrating a statistically significant positive impact on both the moisture forecast and the precipitation forecast at short range (12–36 h) during the warm season. No statistically significant impact is observed during the cold season.

Given increasing computing power, an important question is whether additional computational resources would be better spent reducing the horizontal grid spacing of a convection-allowing model (CAM) or adding members to form CAM ensembles. The present study investigates this question as it applies to CAM-derived next-day probabilistic severe weather forecasts created by using forecast updraft helicity as a severe weather proxy for 63 days of the 2010 and 2011 NOAA Hazardous Weather Testbed Spring Forecasting Experiments. Forecasts derived from three sets of Weather Research and Forecasting Model configurations are tested: a 1-km deterministic model, a 4-km deterministic model, and an 11-member, 4-km ensemble. Forecast quality is evaluated using relative operating characteristic (ROC) curves, attributes diagrams, and performance diagrams, and forecasts from five representative cases are analyzed to investigate their relative quality and value in a variety of situations. While no statistically significant differences exist between the 4- and 1-km deterministic forecasts in terms of area under ROC curves, the 4-km ensemble forecasts offer weakly significant improvements over the 4-km deterministic forecasts over the entire 63-day dataset. Further, the 4-km ensemble forecasts generally provide greater forecast quality relative to either of the deterministic forecasts on an individual day. Collectively, these results suggest that, for purposes of improving next-day CAM-derived probabilistic severe weather forecasts, additional computing resources may be better spent on adding members to form CAM ensembles than on reducing the horizontal grid spacing of a deterministic model below 4 km.

The quality of ensemble streamflow forecasts in the U.S. mid-Atlantic region (MAR) is investigated for short- to medium-range forecast lead times (6–168 h). To this end, a regional hydrological ensemble prediction system (RHEPS) is assembled and implemented. The RHEPS in this case comprises the ensemble meteorological forcing, a distributed hydrological model, and a statistical postprocessor. As the meteorological forcing, precipitation, and near-surface temperature outputs from the National Oceanic and Atmospheric Administration (NOAA)/National Centers for Environmental Prediction (NCEP) 11-member Global Ensemble Forecast System Reforecast, version 2 (GEFSRv2), are used. The Hydrology Laboratory Research Distributed Hydrologic Model (HL-RDHM) is used as the distributed hydrological model, and a statistical autoregressive model with an exogenous variable is used as the postprocessor. To verify streamflow forecasts from the RHEPS, eight river basins in the MAR are selected, ranging in drainage area from ~262 to 29 965 km2 and covering some of the major rivers in the MAR. The verification results for the RHEPS show that, at the initial lead times (1–3 days), the hydrological uncertainties have more impact on forecast skill than the meteorological ones. The former become less pronounced, and the meteorological uncertainties dominate, across longer lead times (>3 days). Nonetheless, the ensemble streamflow forecasts remain skillful for lead times of up to 7 days. Additionally, postprocessing increases forecast skills across lead times and spatial scales, particularly for the high-flow conditions. Overall, the proposed RHEPS is able to improve streamflow forecasting in the MAR relative to the deterministic (unperturbed GEFSRv2 member) forecasting case.

The peak intensity occurrence frequency over the life cycles of parent cold-surge vortices (CSVs) for heavy rainfall/flood (HRF) events is classified into two types depending on their life cycles having two or three peak intensities, denoted as HRF2 or HRF3, respectively. The formation of an HRF2 event from its parent CSV(HRF2) formation is ≤5 days, while the formation of an HRF3 event is ≥6 days. The latter group contributes ~57% of the total number of HRF events. As a result of some model constraints, the formation and development of HRF3 events are not well forecasted by the Global Forecast System (GFS) and regional forecast models. The life cycle and second peak intensity for CSV(HRF3) allow for the introduction of a forecast advisory for HRF3 events. Identification of CSVs and two sufficient requirements for the formation and occurrence of HRF events were developed by previous studies. Nevertheless, two new necessary steps are now included in the proposed forecast advisory. The population ratio for CSV(HRF3) and the regular CSV is only about 15%. The occurrence optimum time to for the CSV(HRF3) second peak intensity from this vortex formation is about 3 days 6 h. The GFS forecast over to is utilized to identify CSV(HRF3). Then, the relay of the GFS forecast from the occurrence time of the CSV(HRF3) second peak is used to predict the formation/occurrence of HRF3 events. Six HRF3 events during cold seasons for 2013–16 are used to test the feasibility of this forecast advisory. Results clearly demonstrate this advisory is a success for the forecast of HRF3 events over the entire life cycles of their parent CSV(HRF3)s.

Since October 2013 a convective-scale weather prediction model has been used operationally to provide short-term forecasts covering large parts of the Nordic region. The model is now operated by a bilateral cooperative effort [Meteorological Cooperation on Operational Numerical Weather Prediction (MetCoOp)] between the Norwegian Meteorological Institute and the Swedish Meteorological and Hydrological Institute. The core of the model is based on the convection-permitting Applications of Research to Operations at Mesoscale (AROME) model developed by Météo-France. In this paper the specific modifications and updates that have been made to suit advanced high-resolution weather forecasts over the Nordic regions are described. This includes modifications in the surface drag description, microphysics, snow assimilation, as well as an update of the ecosystem and surface parameter description. Novel observation types are introduced in the operational runs, including ground-based Global Navigation Satellite System (GNSS) observations and radar reflectivity data from the Norwegian and Swedish radar networks. After almost two years’ worth of experience with the AROME-MetCoOp model, the model’s sensitivities to the use of specific parameterization settings are characterized and the forecast skills demonstrating the benefit as compared with the global European Centre for Medium-Range Weather Forecasts’ Integrated Forecasting System (ECMWF-IFS) are evaluated. Furthermore, case studies are provided to demonstrate the ability of the model to capture extreme precipitation and wind events.

This study evaluates the sensitivity of wind turbine hub-height wind speed forecasts to the planetary boundary layer (PBL) scheme, grid length, and initial condition selection in the Weather Research and Forecasting (WRF) Model over complex terrain. Eight PBL schemes available for the WRF-ARW dynamical core were tested with initial conditions sources from the North American Mesoscale (NAM) model and Global Forecast System (GFS) to produce short-term wind speed forecasts. The largest improvements in forecast accuracy primarily depended on the grid length or PBL scheme choice, although the most important factor varied by location, season, time of day, and bias-correction application. Aggregated over all locations, the Asymmetric Convective Model, version 2 (ACM2) PBL scheme provided the best forecast accuracy, particularly for the 12-km grid length. Other PBL schemes and grid lengths, however, did perform better than the ACM2 scheme for individual seasons or locations.

This study investigates the efficiency of the major operational global ensemble forecast systems of the world in capturing the spatiotemporal evolution of the forecast uncertainty. Using data from 2015, it updates the results of an earlier study based on data from 2012. It also tests, for the first time on operational ensemble data, two quantitative relationships to aid in the interpretation of the raw ensemble forecasts. One of these relationships provides a flow-dependent prediction of the reliability of the ensemble in capturing the uncertain forecast features, while the other predicts the 95th percentile value of the magnitude of the forecast error. It is found that, except for the system of the Met Office, the main characteristics of the ensemble forecast systems have changed little between 2012 and 2015. The performance of the UKMO ensemble improved in predicting the overall magnitude of the uncertainty, but its ability to predict the dominant uncertain forecast features was degraded. A common serious limitation of the ensemble systems remains that they all have major difficulties with predicting the large-scale atmospheric flow in the long (longer than 10 days) forecast range. These difficulties are due to the inability of the ensemble members to maintain large-scale waves in the forecasts, which presents a stumbling block in the way of extending the skill of numerical weather forecasts to the subseasonal range. The two tested predictive relationships were found to provide highly accurate predictions of the flow-dependent reliability of the ensemble predictions and the 95th percentile value of the magnitude of the forecast error for the operational ensemble forecast systems.

The quality of ensemble precipitation forecasts across the eastern United States is investigated, specifically, version 2 of the National Centers for Environmental Prediction (NCEP) Global Ensemble Forecast System Reforecast (GEFSRv2) and Short Range Ensemble Forecast (SREF) system, as well as NCEP’s Weather Prediction Center probabilistic quantitative precipitation forecast (WPC-PQPF) guidance. The forecasts are verified using multisensor precipitation estimates and various metrics conditioned upon seasonality, precipitation threshold, lead time, and spatial aggregation scale. The forecasts are verified, over the geographic domain of each of the four eastern River Forecasts Centers (RFCs) in the United States, by considering first 1) the three systems or guidance, using a common period of analysis (2012–13) for lead times from 1 to 3 days, and then 2) GEFSRv2 alone, using a longer period (2004–13) and lead times from 1 to 16 days. The verification results indicate that, across the eastern United States, precipitation forecast bias decreases and the skill and reliability improve as the spatial aggregation scale increases; however, all the forecasts exhibit some underforecasting bias. The skill of the forecasts is appreciably better in the cool season than in the warm one. The WPC-PQPFs tend to be superior, in terms of the correlation coefficient, relative mean error, reliability, and forecast skill scores, than both GEFSRv2 and SREF, but the performance varies with the RFC and lead time. Based on GEFSRv2, medium-range precipitation forecasts tend to have skill up to approximately day 7 relative to sampled climatology.

At the National Centers for Environmental Prediction, the global ensemble forecasts from the ensemble Kalman filter scheme in the Global Forecast System are applied in a regional three-dimensional (3D) and a four dimensional (4D) ensemble–variational (EnVar) data assimilation system. The application is a one-way variational method using hybrid static and ensemble error covariances. To enhance impact, three new features have been added to the existing EnVar system in the Gridpoint Statistical Interpolation (GSI). First, the constant coefficients that assign relative weight between the ensemble and static background error are now allowed to vary in the vertical. Second, a new formulation is introduced for the ensemble contribution to the analysis surface pressure. Finally, in order to make use of the information in the ensemble mean that is disregarded in the existing EnVar in GSI, the trajectory correction, a novel approach, is introduced. Relative to the application of a 3D variational data assimilation algorithm, a clear positive impact on 1–3-day forecasts is realized when applying 3DEnVar analyses in the North American Mesoscale Forecast System (NAM). The 3DEnVar DA system was operationally implemented in the NAM Data Assimilation System in August 2014. Application of a 4DEnVar algorithm is shown to further improve forecast accuracy relative to the 3DEnVar. The approach described in this paper effectively combines contributions from both the regional and the global forecast systems to produce the initial conditions for the regional NAM system.

A cycling run, which began 36 h before the model forecast, was employed to assimilate special Terrain-influenced Monsoon Rainfall Experiment (TiMREX) soundings, Global Telecommunications System (GTS) data, and Constellation Observing System for Meteorology, Ionosphere and Climate (COSMIC) global positioning system (GPS) radio occultation (RO) refractivity profiles to improve the model initial conditions provided by the National Centers for Environmental Prediction (NCEP) Global Forecast System (GFS) to study a coastal, heavy rainfall event over southwestern Taiwan during 15–16 June 2008. The 36-h cycling run with data assimilation (DA_ALL_DATA run) has a positive impact on the depiction of subsynoptic flow in the model initial conditions at 1200 UTC 15 June, including the warm moist tongue and southwesterly monsoon flow over the open ocean. Furthermore, the cold pool caused by the evaporative cooling of antecedent rains and orographic blocking over southwestern Taiwan are better resolved in the nested high-resolution domain in the DA_ALL_DATA run as compared to the initial conditions provided by the NCEP GFS. As a result, the heavy rainfall along the southwestern coast and afternoon localized heavy rainfall over northern Taiwan are better predicted in the DA_ALL_DATA run. Model sensitivity tests are also performed to diagnose the effects of terrain and rain-evaporative cooling on the intensity and depth of the cold pool and degree of orographic blocking on the southwesterly flow over southwestern Taiwan. It is apparent that including rain-evaporative cooling from antecedent rains and orographic effects in the model initial conditions are important to account for the predicted rainfall distribution of this coastal rainfall event.

The National Hurricane Center (NHC) has stated that guidance on tropical cyclone (TC) genesis is an operational forecast improvement need, particularly since numerical weather prediction models produce TC-like features and operationally required forecast lead times recently have increased. Using previously defined criteria for TC genesis in global models, this study bias corrects TC genesis forecasts from global models using multiple logistic regression. The derived regression equations provide 48- and 120-h probabilistic genesis forecasts for each TC genesis event that occurs in the Environment Canada Global Environmental Multiscale Model (CMC), the NCEP Global Forecast System (GFS), and the Met Office's global model (UKMET). Results show select global model output variables are good discriminators between successful and unsuccessful TC genesis forecasts. Independent verification of the regression-based probabilistic genesis forecasts during 2014 and 2015 are presented. Brier scores and reliability diagrams indicate that the forecasts generally are well calibrated and can be used as guidance for NHC’s Tropical Weather Outlook product. The regression-based TC genesis forecasts are available in real time online.

The capability of all-sky microwave radiance assimilation in the Gridpoint Statistical Interpolation (GSI) analysis system has been developed at the National Centers for Environmental Prediction (NCEP). This development effort required the adaptation of quality control, observation error assignment, bias correction, and background error covariance to all-sky conditions within the ensemble–variational (EnVar) framework. The assimilation of cloudy radiances from the Advanced Microwave Sounding Unit-A (AMSU-A) microwave radiometer for ocean fields of view (FOVs) is the primary emphasis of this study. In the original operational hybrid 3D EnVar Global Forecast System (GFS), the clear-sky approach for radiance data assimilation is applied. Changes to data thinning and quality control have allowed all-sky satellite radiances to be assimilated in the GSI. Along with the symmetric observation error assignment, additional situation-dependent observation error inflation is employed for all-sky conditions. Moreover, in addition to the current radiance bias correction, a new bias correction strategy has been applied to all-sky radiances. In this work, the static background error variance and the ensemble spread of cloud water are examined, and the levels of cloud variability from the ensemble forecast in single- and dual-resolution configurations are discussed. Overall, the all-sky approach provides more realistic simulated brightness temperatures and cloud water analysis increments, and improves analysis off the west coasts of the continents by reducing a known bias in stratus. An approximate 10% increase in the use of AMSU-A channels 1–5 and a 12% increase for channel 15 are also observed. The all-sky AMSU-A radiance assimilation became operational in the 4D EnVar GFS system upgrade of 12 May 2016.

Since November 2014, the Meteorological Services of Canada (MSC) has been running a real-time numerical weather prediction system that provides deterministic forecasts on a regional domain with a 2.5-km horizontal grid spacing covering a large portion of Canada using the Global Environmental Multiscale (GEM) forecast model. This system, referred to as the High Resolution Deterministic Prediction System (HRDPS), is currently downscaled from MSC’s operational 10-km GEM-based regional system but uses initial surface fields from a high-resolution (2.5 km) land data assimilation system coupled to the HRDPS and initial hydrometeor fields from the forecast of a 2.5-km cycle, which reduces the spinup time for clouds and precipitation. Forecast runs of 48 h are provided four times daily. The HRDPS was tested and compared to the operational 10-km system. Model runs from the two systems were evaluated against surface observations for common weather elements (temperature, humidity, winds, and precipitation), fractional cloud cover, and also against upper-air soundings, all using standard metrics. Although the predictions of some fields were degraded in some specific regions, the HRDPS generally outperformed the operational system for a majority of the scores. The evaluation illustrates the added value of the 2.5-km model and the potential for improved numerical guidance for the prediction of high-impact weather.

An evaluation of medium-range forecasts of tropical cyclones (TCs) is performed, covering the eastern North Pacific basin during the period 1 August–3 November 2014. Real-time forecasts from the Model for Prediction Across Scales (MPAS) and operational forecasts from the National Centers for Environmental Prediction (NCEP) Global Forecast System (GFS) are evaluated. A new TC-verification method is introduced that treats TC tracks as objects. The method identifies matching pairs of forecast and observed tracks, missed and false alarm tracks, and derives statistics using a multicategory contingency table methodology. The formalism includes track, intensity, and genesis.Two configurations of MPAS, a uniform 15-km mesh and a variable-resolution mesh transitioning from 60 km globally to 15 km over the eastern Pacific, are compared with each other and with the operational GFS. The two configurations of MPAS reveal highly similar forecast skill and biases through at least day 7. This result supports the effectiveness of TC prediction using variable resolution.Both MPAS and the GFS suffer from biases in predictions of genesis at longer time ranges; MPAS produces too many storms whereas the GFS produces too few. MPAS better discriminates hurricanes than does the GFS, but the false alarms in MPAS lower overall forecast skill in the medium range relative to GFS. The biases in MPAS forecasts are traced to errors in the parameterization of shallow convection south of the equator and the resulting erroneous invigoration of the ITCZ over the eastern North Pacific.

The Joint Polar Satellite System (JPSS) is a key contributor to the next-generation operational polar-orbiting satellite observing system. In the JPSS era, the complete polar-orbiting observing system will be composed of two satellites—in the midmorning (mid-AM) and afternoon (PM) orbits—each with thermodynamic sounding capabilities from both microwave and hyperspectral infrared instruments. JPSS will occupy the PM orbit, while the Meteorological Operational (MetOp) system, sponsored by the European Organisation for the Exploitation of Meteorological Satellites (EUMETSAT), will occupy the mid-AM orbit. While the current polar-orbiting satellite system has been thoroughly evaluated, information about its resilience and efficacy in the JPSS era is needed. A 7-month (August 2012–February 2013) observing system experiment (OSE) was run with the National Centers for Environmental Prediction (NCEP) Global Forecast System (GFS). Observations were selected from operational satellite data platforms to be representative of the polar-orbiting data in the JPSS era. Overall, removing data from the PM orbit produced inferior scores, with the impact greater in the Southern Hemisphere (SH) than in either the Northern Hemisphere (NH) or the tropics. For the entire 7 months, the time-mean 500-hPa geopotential height anomaly correlation (Z500AC) decreased by 0.005 and 0.013 in the NH and SH, respectively—both of which are statistically significant at the 95% level. Additionally, a detailed statistical analysis of the distribution of Z500AC skill scores is presented and compared with historical accuracy data. It was determined that eliminating PM orbit data resulted in a higher probability of producing low scores and a lower probability of producing high scores, counter to the trend in GFS forecast skill over the last 20 years.

The current constellation of environmental satellites is at risk of degrading due to several factors. This includes the following: 1) loss of secondary polar-orbiting satellites due to reaching their nominal lifetimes, 2) decrease in the density of extratropical radio-occultation (RO) observations due to a likely delayed launch of the Constellation Observing System for Meteorology, Ionosphere and Climate-2 (COSMIC-2) high inclination orbit constellation, and 3) the risk of losing afternoon polar-orbiting satellite coverage due to potential launch delays in the Joint Polar Satellite System (JPSS) programs. In this study, the impacts from these scenarios on the National Oceanic and Atmospheric Administration (NOAA) Global Forecast System skill are quantified. Performances for several metrics are assessed, but to encapsulate the results the authors introduce an overall forecast score combining metrics for all parameters, atmospheric levels, and forecast lead times. The first result suggests that removing secondary satellites results in significant degradation of the forecast. This is unexpected since it is generally assumed that secondary sensors contribute to system’s robustness but not necessarily to forecast performance. Second, losing the afternoon orbit on top of losing secondary satellites further degrades forecast performances by a significant margin. Finally, losing extratropical RO observations on top of losing secondary satellites also negatively impacts the forecast performances, but to a lesser degree. These results provide a benchmark that will allow for the assessment of the added value of projects being implemented at NOAA in support of mitigation strategies designed to alleviate the negative impacts associated with these data gaps, and additionally help NOAA to define requirements of the future global observing system architecture.

It is quite common that in a regional ensemble system the large-scale initial condition (IC) perturbations and the lateral boundary condition (LBC) perturbations are taken from a global ensemble prediction system (EPS). The choice of global EPS as a driving model can have a significant impact on the performance of the regional EPS. This study investigates the impact of large-scale IC/LBC perturbations obtained from different global EPSs on the forecast quality of a regional EPS. For this purpose several experiments are conducted where the Aire Limitée Adaption dynamique Développement International–Limited Area Ensemble Forecasting (ALADIN-LAEF) regional ensemble is forced by two of the world’s leading global ensembles, the European Centre for Medium-Range Weather Forecasts’ Ensemble Prediction System (ECMWF-EPS) and the Global Ensemble Forecasting System (GEFS) from the National Centers for Environmental Prediction (NCEP), which provide the IC and LBC perturbations. The investigation is carried out for a 51-day period during summer 2010 over central Europe. The results indicate that forcing of the regional ensemble with GEFS performs better for surface parameters, whereas at upper levels forcing with ECMWF-EPS is superior. Using perturbations from GEFS lead to a considerably higher spread in ALADIN-LAEF, which is beneficial near the surface where regional EPSs are usually underdispersive. At upper levels, forcing with GEFS leads to an overdispersion of ALADIN-LAEF as a result of the large spread of some parameters, where forcing ALADIN-LAEF with ECMWF-EPS provides statistically more reliable forecasts. The results indicate that the best global EPS might not always provide the best ICs and LBCs for a regional ensemble.

Forecast solar exposure fields produced by the Australian Bureau of Meteorology’s updated numerical weather prediction systems were validated against multiple sites for the 2012 calendar year. The updated systems are denoted as the Australian Community Climate and Earth-System Simulator (ACCESS) model and became operational in August 2010. The systems are based on the Met Office’s Unified Model and feature improved assimilation methods and radiation parameterizations that were expected to greatly improve forecasts of solar exposure several days in advance. In this study forecasts of global, direct, and diffuse exposure from the mesoscale model ACCESS-A were validated. Statistics were generated for all-sky and clear-sky conditions. Additionally, evaluation of the model’s forecast exposure through single-layer low clouds was conducted. Results show an improvement in global forecasts relative to the older operational model; however, forecasts of diffuse and direct exposure still suffer from large biases. These can be attributed to the choices of the asymmetry factor used in the two-stream approximation for incoming radiation, which determines scattering of the direct beam through clouds.

A global deterministic wave prediction system (GDWPS) is used to improve regional forecasts of waves off the Canadian coastline and help support the practice of safe marine activities in Canadian waters. The wave model has a grid spacing of ¼° with spectral resolution of 36 frequency bins and 36 directional bins. The wave model is driven with hourly 10-m winds generated by the operational global atmospheric prediction system. Ice conditions are updated every three hours using the ice concentration forecasts generated by the Global Ice–Ocean Prediction System. Wave forecasts are evaluated over two periods from 15 August to 31 October 2014 and from 15 December 2014 to 28 February 2015, as well as over select cases during the fall of 2012. The global system is shown to improve wave forecast skill over regions where forecasts were previously produced using limited-area models only. The usefulness of a global expansion is demonstrated for large swell events affecting the northeast Pacific. The first validation of a Canadian operational wave forecast system in the Arctic is presented. Improvements in the representation of forecast wave fields associated with tropical cyclones are also demonstrated. Finally, the GDWPS is shown to result in gains of at least 12 h of lead time.

Improving the initial conditions of short-range numerical weather prediction (NWP) models is one of the main goals of the meteorological community. Development of data assimilation and ensemble forecast systems is essential in any national weather service (NWS). In this sense, the local ensemble transform Kalman filter (LETKF) is a methodology that can satisfy both requirements in an efficient manner. The Weather Research and Forecasting (WRF) Model coupled with the LETKF, developed at the University of Maryland, College Park, have been implemented experimentally at the NWS of Argentina [Servicio Meteorológico Nacional (SMN)], but at a somewhat lower resolution (40 km) than the operational Global Forecast System (GFS) at that time (27 km). The purpose of this work is not to show that the system presented herein is better than the higher-resolution GFS, but that its performance is reasonably comparable, and to provide the basis for a continued improved development of an independent regional data assimilation and forecasting system. The WRF-LETKF system is tested during the spring of 2012, using the prepared or quality controlled data in Binary Universal Form for Representation of Meteorological Data (PREPBUFR) observations from the National Centers for Environmental Prediction (NCEP) and lateral boundary conditions from the GFS. To assess the effect of model error, a single-model LETKF system (LETKF-single) is compared with a multischeme implementation (LETKF-multi), which uses different boundary layer and cumulus convection schemes for the generation of the ensemble of forecasts. The performance of both experiments during the test period shows that the LETKF-multi usually outperforms the LETKF-single, evidencing the advantages of the use of the multischeme approach. Both data assimilation systems are slightly worse than the GFS in terms of the synoptic environment representation, as could be expected given their lower resolution. Results from a case study of a strong convective system suggest that the LETKF-multi improves the location of the most intense area of precipitation with respect to the LETKF-single, although both systems show an underestimation of the total accumulated precipitation. These preliminary results encourage continuing the development of an operational data assimilation system based on WRF-LETKF at the SMN. Fil: Dillon, María Eugenia. Ministerio de Defensa. Secretaria de Planeamiento. Servicio Meteorológico Nacional; Argentina. Universidad de Buenos Aires. Facultad de Ciencias Exactas y Naturales. Departamento de Ciencias de la Atmósfera y los Océanos; Argentina. Consejo Nacional de Investigaciones Científicas y Técnicas; Argentina Fil: Garcia Skabar, Yanina. Ministerio de Defensa. Secretaria de Planeamiento. Servicio Meteorológico Nacional; Argentina. Consejo Nacional de Investigaciones Científicas y Técnicas; Argentina Fil: Ruiz, Juan Jose. Consejo Nacional de Investigaciones Científicas y Técnicas. Oficina de Coordinación Administrativa Ciudad Universitaria. Centro de Investigaciones del Mar y la Atmósfera. Universidad de Buenos Aires. Facultad de Ciencias Exactas y Naturales. Centro de Investigaciones del Mar y la Atmósfera; Argentina Fil: Kalnay, Eugenia. University of Maryland; Estados Unidos Fil: Collini, Estela Angela. Ministerio de Defensa. Armada Argentina. Servicio de Hidrografía Naval; Argentina. Ministerio de Defensa. Secretaria de Planeamiento. Servicio Meteorológico Nacional; Argentina Fil: Echevarria, Pablo. Ministerio de Defensa. Secretaria de Planeamiento. Servicio Meteorológico Nacional; Argentina Fil: Saucedo, Marcos Adolfo. Consejo Nacional de Investigaciones Científicas y Técnicas. Oficina de Coordinación Administrativa Ciudad Universitaria. Centro de Investigaciones del Mar y la Atmósfera. Universidad de Buenos Aires. Facultad de Ciencias Exactas y Naturales. Centro de Investigaciones del Mar y la Atmósfera; Argentina Fil: Miyoshi, Takemasa. RIKEN Advanced Institute for Computational Science; Japón Fil: Kunii, Masaru. RIKEN Advanced Institute for Computational Science; Japón

Probabilistic severe weather forecasts for days 1 and 2 were produced using 30-member convection-allowing ensemble forecasts initialized by an ensemble Kalman filter data assimilation system during a 32-day period coinciding with the Mesoscale Predictability Experiment. The forecasts were generated by smoothing the locations where model output indicated extreme values of updraft helicity, a surrogate for rotating thunderstorms in model output. The day 1 surrogate severe probability forecasts (SSPFs) produced skillful and reliable predictions of severe weather during this period, after an appropriate calibration of the smoothing kernel. The ensemble SSPFs exceeded the skill of SSPFs derived from two benchmark deterministic forecasts, with the largest differences occurring on the mesoscale, while all SSPFs produced similar forecasts on synoptic scales. While the deterministic SSPFs often overforecasted high probabilities, the ensemble improved the reliability of these probabilities, at the expense of producing fewer high-probability values. For the day 2 period, the SSPFs provided competitive guidance compared to the day 1 forecasts, although additional smoothing was needed to produce the same level of skill, reducing the forecast sharpness. Results were similar using 10 ensemble members, suggesting value exists when running a smaller ensemble if computational resources are limited. Finally, the SSPFs were compared to severe weather risk areas identified in Storm Prediction Center (SPC) convective outlooks. The SSPF skill was comparable to the SPC outlook skill in identifying regions where severe weather would occur, although performance varied on a day-to-day basis.

Current methods of assimilation of precipitation into numerical weather prediction models are able to make the model precipitation become similar to the observed precipitation during the assimilation, but the model forecasts tend to return to their original solution after a few hours. To facilitate the precipitation assimilation, a logarithm transformation has been used in several past studies. Lien et al. proposed instead to assimilate precipitation using the local ensemble transform Kalman filter (LETKF) with a Gaussian transformation technique and succeeded in improving the model forecasts in perfect-model observing system simulation experiments (OSSEs). In this study, the method of Lien et al. is tested within a more realistic configuration: the TRMM Multisatellite Precipitation Analysis (TMPA) data are assimilated into a low-resolution version of the NCEP Global Forecast System (GFS). With guidance from a statistical study comparing the GFS model background precipitation and the TMPA data, some modifications of the assimilation methods proposed in Lien et al. are made, including 1) applying separate Gaussian transformations to model and to observational precipitation based on their own cumulative distribution functions; 2) adopting a quality control criterion based on the correlation between the long-term model and observed precipitation data at the observation location; and 3) proposing a new method to define the transformation of zero precipitation that takes into account the zero precipitation probability in the background ensemble rather than the climatology. With these modifications, the assimilation of the TMPA precipitation data improves both the analysis and 5-day model forecasts when compared with a control experiment assimilating only rawinsonde data.

The current operational eddy-diffusivity countergradient (EDCG) planetary boundary layer (PBL) scheme in the NCEP Global Forecast System (GFS) tends to underestimate the PBL growth in the convective boundary layer (CBL). To improve CBL growth, an eddy-diffusivity mass-flux (EDMF) PBL scheme is developed, where the nonlocal transport by large turbulent eddies is represented by a mass-flux (MF) scheme and the local transport by small eddies is represented by an eddy-diffusivity (ED) scheme. For the vertical momentum mixing, the MF scheme is modified to include the effect of the updraft-induced pressure gradient force. While the EDMF scheme displays better CBL growth than the EDCG scheme, it tends to overproduce the amount of low clouds and degrades wind vector forecasts over the tropical ocean where strongly unstable PBLs are rarely found. In order not to degrade the forecast skill in the tropics, a hybrid scheme is developed, where the EDMF scheme is applied only for the strongly unstable PBL, while the EDCG scheme is used for the weakly unstable PBL. Along with the hybrid EDMF scheme, the heating by turbulent kinetic energy (TKE) dissipation is parameterized to reduce an energy imbalance in the GFS. To enhance a too weak vertical turbulent mixing for weakly and moderately stable conditions, the current local scheme in the stable boundary layer (SBL) is modified to use an eddy-diffusivity profile method. The hybrid EDMF PBL scheme with TKE dissipative heating and modified SBL mixing led to significant improvements in some key medium-range weather forecast metrics and was operationally implemented into the NCEP GFS in January 2015.

Recently, two analog-based postprocessing methods were demonstrated to reduce the systematic and random errors from Weather Research and Forecasting (WRF) Model predictions of 10-m wind speed over the central United States. To test robustness and generality, and to gain a deeper understanding of postprocessing forecasts with analogs, this paper expands upon that work by applying both analog methods to surface stations evenly distributed across the conterminous United States over a 1-yr period. The Global Forecast System (GFS), North American Mesoscale Forecast System (NAM), and Rapid Update Cycle (RUC) forecasts for screen-height wind, temperature, and humidity are postprocessed with the two analog-based methods and with two time series–based methods—a running mean bias correction and an algorithm inspired by the Kalman filter. Forecasts are evaluated according to a range of metrics, including random and systematic error components; correlation; and by conditioning the error distributions on lead time, location, error magnitude, and day-to-day error variability.Results show that the analog methods are generally more effective than time series–based methods at reducing the random error component, leading to an overall reduction in root-mean-square error. Details among the methods differ and are elucidated upon in this study. The relative levels of random and systematic error in the raw forecasts determine, to a large extent, the effectiveness of each postprocessing method in reducing forecast errors. When the errors are dominated by random errors (e.g., where thunderstorms are common), the analog-based methods far outperform the time series–based methods. When the errors are strictly systematic (i.e., a bias), the analog methods lose their advantage over the time series methods. It is shown that slowly evolving systematic errors rarely dominate, so reducing the random error component is most effective at reducing the error magnitude. The results are shown to be valid for all seasons. The analog methods show similar performance to the operational model output statistics (MOS) while showing greater reduction of random errors at certain lead times.

Traditional pixel-versus-pixel forecast evaluation scores such as the critical success index (CSI) provide a simple way to compare the performances of different forecasts; however, they offer little information on how to improve a particular forecast. This paper strives to demonstrate what additional information an object-based forecast evaluation tool such as the Method for Object-Based Diagnostic Evaluation (MODE) can provide in terms of assessing numerical weather prediction models’ convective storm forecasts. Forecast storm attributes evaluated by MODE in this paper include storm size, intensity, orientation, aspect ratio, complexity, and number of storms. Three weeks of the High Resolution Rapid Refresh (HRRR) model’s precipitation forecasts during the summer of 2010 over the eastern two-thirds of the contiguous United States were evaluated as an example to demonstrate the methodology. It is found that the HRRR model was able to forecast convective storm characteristics rather well either as a function of time of day or as a function of storm size, although significant bias does exist, especially in terms of storm number and storm size. Another interesting finding is that the model’s ability of forecasting new storm initiation varies substantially by regions, probably as a result of its different skills in forecasting convection driven by different forcing mechanisms (i.e., diurnal heating vs synoptic-scale frontal systems).