Search

The Gridpoint Statistical Interpolation (GSI)-based four- and three-dimensional ensemble–variational (4DEnVar and 3DEnVar) methods are compared as a smoother and a filter, respectively, for rapidly changing storms using the convective-scale direct radar reflectivity data assimilation (DA) framework. Two sets of experiments with varying DA window lengths (WLs; 20, 40, 100, and 160 min) and radar observation intervals (RIs; 20 and 5 min) are conducted for the 5–6 May 2019 case. The RI determines the temporal resolution of ensemble perturbations for the smoother and the DA interval for the filter spanning the WL. For experiments with a 20-min RI, evaluations suggest that the filter and the smoother have comparable performance with a 20-min WL; however, extending the WL results in the outperformance of the filter over the smoother. Diagnostics reveal that the degradation of the smoother is attributed to the increased degree of nonlinearity and the issue of time-independent localization as the WL extends. Evaluations for experiments with different RIs under the same WL indicate that the outperformance of the filter over the smoother diminishes for most forecast hours at thresholds of 30 dBZ and above when shortening the RI. Diagnostics show that more frequent interruptions of the model introduce model imbalance for the filter, and the increased temporal resolution of ensemble perturbations enhances the degree of nonlinearity for the smoother. The impact of model imbalance on the filter overwhelms the enhanced nonlinearity on the smoother as the RI reduces. Significance Statement: The background uncertainties of rapidly changing storms suffer from fast error growth and high degrees of nonlinearities during the data assimilation (DA) period. Two variants of the ensemble-based DA method can account for such temporal evolution. The smoother uses background ensemble from multiple observation times over an assimilation period to estimate the propagation of statistics. The filter frequently calculates the statistics at multiple observation times over the same period. Current comparisons of the smoother and the filter were mostly performed using simple models; however, unknowns remain for convection-allowing forecasts with additional complexities. This study compares the filter and the smoother for the convective-scale analysis and prediction using a real-data study and finds that the comparison varies with the assimilation period and the observation interval. [ABSTRACT FROM AUTHOR]

The Advanced Baseline Imager (ABI) aboard the GOES-16 and GOES-17 satellites provides high-resolution observations of cloud structures that could be highly beneficial for convective-scale DA. However, only clear-air radiance observations are typically assimilated at operational centers due to a variety of problems associated with cloudy radiance data. As such, many questions remain about how to best assimilate all-sky radiance data, especially when using hybrid DA systems such as EnVar wherein a nonlinear observation operator can lead to cost function gradient imbalance and slow minimization. Here, we develop new methods for assimilating all-sky radiance observations in EnVar using the novel Rapid Refresh Forecasting System (RRFS) that utilizes the Finite-Volume Cubed-Sphere (FV3) model. We first modify the EnVar solver by directly including brightness temperature (Tb) as a state variable. This modification improves the balance of the cost function gradient and speeds up minimization. Including Tb as a state variable also improves the model fit to observations and increases forecast skill compared to utilizing a standard state vector configuration. We also evaluate the impact of assimilating ABI all-sky radiances in RRFS for a severe convective event in the central Great Plains. Assimilating the radiance observations results in better spinup of a tornadic supercell. These data also aid in suppressing spurious convection by reducing the snow hydrometeor content near the tropopause and weakening spurious anvil clouds. The all-sky radiance observations pair well with reflectivity observations that remove primarily liquid hydrometeors (i.e., rain) closer to the surface. Additionally, the benefits of assimilating the ABI observations continue into the forecast period, especially for localized convective events. [ABSTRACT FROM AUTHOR]

The current global operational four-dimensional ensemble-variational (4DEnVar) data assimilation (DA) system at NCEP adopts a background ensemble at a reduced resolution, which restricts the range of spatial scales that the ensemble background error covariance can resolve. A prior study developed a multiresolution ensemble 4DEnVar method and determined that this approach can provide a comparable forecast to an approach using solely high-resolution members, while substantially reducing the computational cost. This study further develops the multiresolution ensemble 4DEnVar approach to allow for a flexible number of low- and high-resolution ensemble members as well as varying localization length scales between the high- and low-resolution ensembles. Three 4DEnVar experiments with the same computational costs are compared. The first experiment has an 80-member high-resolution background ensemble with single-scale optimally tuned localization (SR-High). The second and third experiments utilize the multiresolution background ensembles. One has 130 low-resolution and 40 high-resolution members (MR170) while the other has 180 low-resolution members and 24 high-resolution members (MR204). Both multiresolution ensemble experiments utilize differing localization radii with ensemble resolution. Despite having the same costs, both MR170 and MR204 improves global forecasts and decreases tropical cyclone track errors for up to 5 days' lead time compared to SR-High. Improvements are most apparent in larger-scale features, such as jet streams and the environmental steering flow of tropical cyclones. Additionally, MR170 outperforms MR204 in terms of global and tropical cyclone track forecasts, demonstrating the value of both increasing sampling at large scales and retaining substantial information at small scales. [ABSTRACT FROM AUTHOR]

A valid time shifting (VTS) method is explored for the GSI-based ensemble variational (EnVar) system modified to directly assimilate radar reflectivity at convective scales. VTS is a cost-efficient method to increase ensemble size by including subensembles before and after the central analysis time. Additionally, VTS addresses common time and phase model error uncertainties within the ensemble. VTS is examined here for assimilating radar reflectivity in a continuous hourly analysis system for a case study of 1–2 May 2019. The VTS implementation is compared against a 36-member control experiment (ENS-36), to increase ensemble size (3 × 36 VTS), and as a cost-savings method (3 × 12 VTS), with time-shifting intervals τ between 15 and 120 min. The 3 × 36 VTS experiments increased the ensemble spread, with largest subjective benefits in early cycle analyses during convective development. The 3 × 12 VTS experiments captured analysis with similar accuracy as ENS-36 by the third hourly analysis. Control forecasts launched from hourly EnVar analyses show significant skill increases in 1-h precipitation over ENS-36 out to hour 12 for 3 × 36 VTS experiments, subjectively attributable to more accurate placement of the convective line. For 3 × 12 VTS, experiments with τ ≥ 60 min met and exceeded the skill of ENS-36 out to forecast hour 15, with VTS-3 × 12τ90 maximizing skill. Sensitivity results demonstrate preference to τ = 30–60 min for 3 × 36 VTS and 60–120 min for 3 × 12 VTS. The best 3 × 36 VTS experiments add a computational cost of 45%–67%, compared to the near tripling of costs when directly increasing ensemble size, while best 3 × 12 VTS experiments save about 24%–41% costs over ENS-36. Significance Statement: The purpose of this work is to study a valid time shifting method to improve the prediction of severe convective storm systems over the continental United States. This method improves ensemble-based radar reflectivity analyses by including ensemble member information at times before and after the analysis time, thereby increasing the ensemble size at just a fractional added computational cost. The results show the method can boost the accuracy of high-resolution convection prediction out to at least 12 h. This case study motivates future systematic testing in a real-time setting and potential implementation to enhance the U.S. operational ensemble-based convection-allowing forecast model and data assimilation system. [ABSTRACT FROM AUTHOR]

Short-term spinup for strong storms is a known difficulty for the operational Hurricane Weather Research and Forecasting (HWRF) Model after assimilating high-resolution inner-core observations. Our previous study associated this short-term intensity prediction issue with the incompatibility between the HWRF Model and the data assimilation (DA) analysis. While improving physics and resolution of the model was found to be helpful, this study focuses on further improving the intensity predictions through the four-dimensional incremental analysis update (4DIAU). In the traditional 4DIAU, increments are predetermined by subtracting background forecasts from analyses. Such predetermined increments implicitly require linear evolution assumption during the update, which are hardly valid for rapidly evolving hurricanes. To confirm the hypothesis, a corresponding 4D analysis nudging (4DAN) method, which uses online increments is first compared with the 4DIAU in an oscillation model. Then, variants of 4DIAU are proposed to improve its application for nonlinear systems. Next, 4DIAU, 4DAN and their proposed improvements are implemented into the HWRF 4DEnVar DA system and are investigated with Hurricane Patricia (2015). Results from both the oscillation model and HWRF Model show that 1) the predetermined increments in 4DIAU can be detrimental when there are discrepancies between the updated and background forecasts during a nonlinear evolution; 2) 4DAN can improve the performance of incremental update upon 4DIAU, but its improvements are limited by the overfiltering; 3) relocating initial background before the incremental update can improve the corresponding traditional methods; and 4) the feature-relative 4DIAU method improves the incremental update the most and produces the best track and intensity predictions for Patricia among all experiments. [ABSTRACT FROM AUTHOR]

A convective-scale static background-error covariance (BEC) matrix is further developed to include the capability of direct reflectivity assimilation and evaluated within the GSI-based three-dimensional variational (3DVar) and hybrid ensemble–variational (EnVar) methods. Specific developments are summarized as follows: 1) Control variables (CVs) are extended to include reflectivity, vertical velocity, and all hydrometeor types. Various horizontal momentum and moisture CV options are included. 2) Cross correlations between all CVs are established. 3) A storm intensity-dependent binning method is adopted to separately calculate static error matrices for clear-air and storms with varying intensities. The resultant static BEC matrices are simultaneously applied at proper locations guided by the observed reflectivity. 4) The EnVar is extended to adaptively incorporate static BECs based on the quality of ensemble covariances. Evaluation and examination of the new static BECs are first performed on the 8 May 2003 Oklahoma City supercell. Detailed diagnostics and 3DVar examinations suggest zonal/meridional winds and pseudo–relative humidity are selected as horizontal momentum and moisture CVs for direct reflectivity assimilation, respectively; inclusion of cross correlations favors spin up and maintains the analyzed storms; application of binning improves characteristics and persistence of the simulated storm. Relative to an experiment using the full ensemble BECs (Exp-PureEnVar), incorporating static BECs in hybrid EnVar reduces spinup time and better analyzes reflectivity distributions while the background ensemble is deficient in sampling errors. Compared to both pure 3DVar and Exp-PureEnVar, hybrid EnVar better predicts reflectivity distributions and better maintains a strong mesocyclone. Further examination through the 20 May 2013 Oklahoma supercells confirms these results and additionally demonstrates the effectiveness of adaptive hybridization. [ABSTRACT FROM AUTHOR]

Although numerous studies have demonstrated that increasing model spatial resolution in free forecasts can potentially improve tropical cyclone (TC) intensity forecasts, studies on the impact of model resolution during data assimilation (DA) on TC prediction are lacking. In this study, using the ensemble-variational DA system for the Hurricane Weather Research and Forecasting (HWRF) Model, we investigated the individual impact of increasing the model resolution of first guess (FG) and background ensemble (BE) forecasts during DA on initial analyses and subsequent forecasts of Hurricane Patricia (2015). The impacts were compared between horizontal and vertical resolutions and also between the tropical storm (TS) and hurricane assimilation during Patricia. The results show that increasing the horizontal or vertical resolution in FG has a larger impact than increasing the resolution in BE on improving the analyzed TC intensity and structure for the hurricane stage. The result is reversed for the TS stage. These results are attributed to the effectiveness of increasing the FG resolution in intensifying the background vortex for the hurricane stage relative to the TS stage. Increasing the BE resolution contributes to improving the analyzed intensity through the better-resolved background correlation structure for both the hurricane and TS stages. Increasing horizontal resolution has an overall larger effect than increasing vertical resolution in improving the analysis at the hurricane stage and their effects are close for the analysis at the TS stage. Additionally, the more accurately analyzed primary circulation, secondary circulation, and warm-core structures via the increased resolution in DA lead to improved TC intensity forecasts. [ABSTRACT FROM AUTHOR]

A scale-dependent localization (SDL) method was formulated and implemented in the Gridpoint Statistical Interpolation (GSI)-based four-dimensional ensemble-variational (4DEnVar) system for NCEP FV3-based Global Forecast System (GFS). SDL applies different localization to different scales of ensemble covariances, while performing a single-step simultaneous assimilation of all available observations. Two SDL variants with (SDL-Cross) and without (SDL-NoCross) considering cross-wave-band covariances were examined. The performance of two- and three-wave-band SDL experiments (W2 and W3, respectively) was evaluated through 1-month cycled data assimilation experiments. SDL improves global forecasts to 5 days over scale-invariant localization including the operationally tuned level-dependent scale-invariant localization (W1-Ope). The W3 SDL-Cross experiment shows more accurate tropical storm–track forecasts at shorter lead times than W1-Ope. Compared to the W2 SDL experiments, the W3 SDL counterparts applying tighter horizontal localization at medium-scale wave band generally show improved global forecasts below 100 hPa, but degraded global forecasts above 50 hPa. While the outperformance of the W3 SDL-NoCross experiment versus the W2 SDL-NoCross experiment below 100 hPa lasts for 5 days, that of the W3 SDL-Cross experiment versus the W2 SDL-Cross experiment lasts for 3 days. Due to local spatial averaging of ensemble covariances that may alleviate sampling error, the SDL-NoCross experiments show slightly better forecasts than the SDL-Cross experiments at shorter lead times. However, the SDL-Cross experiments outperform the SDL-NoCross experiments at longer lead times, likely from retention of more heterogeneity of ensemble covariances and resultant analyses with improved balance. Relative performance of tropical storm–track forecasts in the W2 and W3 SDL experiments are generally consistent with that of global forecasts. [ABSTRACT FROM AUTHOR]

Nocturnal convection is often initiated by mechanisms that cannot be easily observed within the large gaps between rawinsondes or by conventional surface networks. To improve forecasts of such events, we evaluate the systematic impact of assimilating a collocated network of high-frequency, ground-based thermodynamic and kinematic profilers collected as part of the 2015 Plains Elevated Convection At Night (PECAN) experiment. For 13 nocturnal convection initiation (CI) events, we find small but consistent improvements when assimilating thermodynamic observations collected by Atmospheric Emitted Radiance Interferometers (AERIs). Through midlevel cooling and moistening, assimilating the AERIs increases the fractions skill score (FSS) for both nocturnal CI and precipitation forecasts. The AERIs also improve various contingency metrics for CI forecasts. Assimilating composite kinematic datasets collected by Doppler lidars and radar wind profilers (RWPs) results in slight degradations to the forecast quality, including decreases in the FSS and traditional contingency metrics. The impacts from assimilating thermodynamic and kinematic profilers often counteract each other, such that we find little impact on the detection of CI when both are assimilated. However, assimilating both datasets improves various properties of the CI events that are successfully detected (timing, distance, shape, etc.). We also find large variability in the impact of assimilating these remote sensing profilers, likely due to the number of observing sites and the strength of the synoptic forcing for each case. We hypothesize that the lack of flow-dependent methods to diagnose observation errors likely contributes to degradations in forecast skill for many cases, especially when assimilating kinematic profilers. [ABSTRACT FROM AUTHOR]

The potential future installation of a multifunction phased-array radar (MPAR) network will provide capabilities of case-specific adaptive scanning. Knowing the impacts adaptive scanning may have on short-term forecasts will influence scanning strategy decision-making in hopes to produce the most optimal ensemble forecast while also benefiting human severe weather warning decision-making. An ensemble-based targeted observation algorithm is applied to an observing system simulation experiment (OSSE) where the impacts of synthetic idealized supercell radial velocity observations are estimated before the observations are "collected" and assimilated. The forecast metric of interest is the low-level rotation forecast metric (0–1-km updraft helicity), a surrogate for tornado prediction. It is found that the ensemble-based targeted observation approach can reasonably estimate the true error variance reduction when an effective method that treats sampling error is applied, the period of model forecast is associated with less degrees of nonlinearity, and the observation information content relative to the background forecast is larger. In some scenarios, a subset of a full-volume scan assimilation produces better forecasts than all observations within the full volume. Assimilating the full-volume scan increases the number of potential spurious correlations arising between the forecast metric and radial velocity observation induced state perturbations, which may degrade the forecast metric accuracy. [ABSTRACT FROM AUTHOR]

A gridpoint statistical interpolation (GSI)-based hybrid ensemble–variational (EnVar) scheme was extended for convective scales—including radar reflectivity assimilation—and implemented in real-time spring forecasting experiments. This study compares methods to address model error during the forecast under the context of multiscale initial condition error sampling provided by the EnVar system. A total of 10 retrospective cases were used to explore the optimal design of convection-allowing ensemble forecasts. In addition to single-model single-physics (SMSP) configurations, ensemble forecast experiments compared multimodel (MM) and multiphysics (MP) approaches. Stochastic physics was also applied to MP for further comparison. Neighborhood-based verification of precipitation and composite reflectivity showed each of these model error techniques to be superior to SMSP configurations. Comparisons of MM and MP approaches had mixed findings. The MM approach had better overall skill in heavy-precipitation forecasts; however, MP ensembles had better skill for light (2.54 mm) precipitation and reduced ensemble mean error of other diagnostic fields, particularly near the surface. The MM experiment had the largest spread in precipitation, and for most hours in other fields; however, rank histograms and spaghetti contours showed significant clustering of the ensemble distribution. MP plus stochastic physics was able to significantly increase spread with time to be competitive with MM by the end of the forecast. The results generally suggest that an MM approach is best for early forecast lead times up to 6–12 h, while a combination of MP and stochastic physics approaches is preferred for forecasts beyond 6–12 h. [ABSTRACT FROM AUTHOR]

An object-based probabilistic (OBPROB) forecasting framework is developed and applied, together with a more traditional neighborhood-based framework, to convection-permitting ensemble forecasts produced by the University of Oklahoma (OU) Multiscale data Assimilation and Predictability (MAP) laboratory during the 2017 and 2018 NOAA Hazardous Weather Testbed Spring Forecasting Experiments. Case studies from 2017 are used for parameter tuning and demonstration of methodology, while the 2018 ensemble forecasts are systematically verified. The 2017 case study demonstrates that the OBPROB forecast product can provide a unique tool to operational forecasters that includes convective-scale details such as storm mode and morphology, which are typically lost in neighborhood-based methods, while also providing quantitative ensemble probabilistic guidance about those details in a more easily interpretable format than the more commonly used paintball plots. The case study also demonstrates that objective verification metrics reveal different relative performance of the ensemble at different forecast lead times depending on the verification framework (i.e., object versus neighborhood) because of the different features emphasized by object- and neighborhood-based evaluations. Both frameworks are then used for a systematic evaluation of 26 forecasts from the spring of 2018. The OBPROB forecast verification as configured in this study shows less sensitivity to forecast lead time than the neighborhood forecasts. Both frameworks indicate a need for probabilistic calibration to improve ensemble reliability. However, lower ensemble discrimination for OBPROB than the neighborhood-based forecasts is also noted. [ABSTRACT FROM AUTHOR]

A multiresolution ensemble (MR-ENS) method is developed to resolve a wider range of scales of the background error covariance (BEC) in the hybrid four-dimensional ensemble–variational (4DEnVar) while saving computational costs. MR-ENS is implemented in the NCEP Global Forecast System (GFS) gridpoint statistical interpolation (GSI) hybrid 4DEnVar. MR-ENS generates analysis increment by incorporating high-resolution static BEC and flow-dependent ensemble BECs from both high and low resolutions. MR-ENS is compared with three 4DEnVar update approaches: 1) the single-resolution (SR)-Low approach where the analysis increments are generated from the ensemble BEC and the static BEC at the same low resolution; 2) the dual-resolution (DR) approach where the analysis increment is generated using the high-resolution static BEC and low-resolution ensemble BEC; and 3) the SR-High approach, which is the same as 1) except that all covariances are at high-resolution. Experiments show that MR-ENS improves global and tropical cyclone track forecasts compared to SR-Low and DR. Inclusion of the high-resolution ensemble leads to increased background ensemble spread, better fitting of the background to observations, increased effective ranks, more accurate ensemble error correlation, and increased power of analysis increment at small scales. The majority of the improvement of MR-ENS relative to SR-Low is due to the partial use of high-resolution background ensemble. Compared to SR-High, MR-ENS decreases the overall cost by about 40% and shows comparable global and tropical cyclone track forecast performances. Diagnostics show that particularly in the tropics, MR-ENS improves the analysis increment over a wide range of scales and increases the effective rank of the ensemble BEC to the degree comparable to SR-High. [ABSTRACT FROM AUTHOR]

Numerical weather prediction models often fail to correctly forecast convection initiation (CI) at night. To improve our understanding of such events, researchers collected a unique dataset of thermodynamic and kinematic remote sensing profilers as part of the Plains Elevated Convection at Night (PECAN) experiment. This study evaluates the impacts made to a nocturnal CI forecast on 26 June 2015 by assimilating a network of atmospheric emitted radiance interferometers (AERIs), Doppler lidars, radio wind profilers, high-frequency rawinsondes, and mobile surface observations using an advanced, ensemble-based data assimilation system. Relative to operational forecasts, assimilating the PECAN dataset improves the timing, location, and orientation of the CI event. Specifically, radio wind profilers and rawinsondes are shown to be the most impactful instrument by enhancing the moisture advection into the region of CI in the forecast. Assimilating thermodynamic profiles collected by the AERIs increases midlevel moisture and improves the ensemble probability of CI in the forecast. The impacts of assimilating the radio wind profilers, AERI retrievals, and rawinsondes remain large throughout forecasting the growth of the CI event into a mesoscale convective system. Assimilating Doppler lidar and surface data only slightly improves the CI forecast by enhancing the convergence along an outflow boundary that partially forces the nocturnal CI event. Our findings suggest that a mesoscale network of profiling and surface instruments has the potential to greatly improve short-term forecasts of nocturnal convection. [ABSTRACT FROM AUTHOR]

The dropsondes released during the Tropical Cyclone Intensity (TCI) field campaign provide high-resolution kinematic and thermodynamic measurements of tropical cyclones within the upper-level outflow and inner core. This study investigates the impact of these upper-level TCI dropsondes on analyses and prediction of Hurricane Patricia (2015) during its rapid intensification (RI) phase using an ensemble–variational data assimilation system. In the baseline experiment (BASE), both kinematic and thermodynamic observations of TCI dropsondes at all levels except the upper levels are assimilated. The upper-level wind and thermodynamic observations are assimilated in additional experiments to investigate their respective impacts. Compared to BASE, assimilating TCI upper-level wind observations improves the accuracy of outflow analyses verified against independent atmospheric motion vector (AMV) observations. It also strengthens the tangential and radial wind near the upper-level eyewall. The inertial stability within the upper-level eyewall is enhanced, and the maximum outflow is more aligned toward the inner core. Additionally, the analyses including the upper-level thermodynamic observations produce a warmer and drier core at high levels. Assimilating both upper-level kinematic and thermodynamic observations also improves the RI forecast. Compared to BASE, assimilating the upper-level wind induces more upright and inward-located eyewall convection, resulting in more latent heat release closer to the warm core. This process leads to stronger inner-core warming. Additionally, the initial warmer upper-level inner core produced by assimilating TCI thermodynamic observations also intensifies the convection and latent heat release within the eyewall, thus further contributing to the improved intensity forecasts. [ABSTRACT FROM AUTHOR]

Four case studies from the Plains Elevated Convection at Night (PECAN) field experiment are used to investigate the impacts of horizontal and vertical resolution, and vertical mixing parameterization, on predictions of bore structure and upscale impacts of bores on their mesoscale environment. The reduction of environmental convective inhibition (CIN) created by the bore is particularly emphasized. Simulations are run with horizontal grid spacings ranging from 250 to 1000 m, as well as 50 m for one case study, different vertical level configurations, and different closure models for the vertical turbulent mixing at 250-m horizontal resolution. The 11 July case study was evaluated in greatest detail because it was the best observed case and has been the focus of a previous study. For this case, it is found that 250-m grid spacing improves upon 1-km grid spacing, LES configuration provides further improvement, and enhanced low-level vertical resolution also provides further improvement in terms of qualitative agreement between simulated and observed bore structure. Reducing LES grid spacing further to 50 m provided very little additional advantage. Only the LES experiments properly resolved the upscale influence of reduced low-level CIN. Expanding on the 11 July case study, three other cases from PECAN with diverse observed bore structures were also evaluated. Similar to the 11 July case, enhancing the horizontal and vertical grid spacings, and using the LES closure model for vertical turbulent mixing, all contributed to improved simulations of both the bores themselves and the larger-scale modification of CIN to varying degrees on different cases. [ABSTRACT FROM AUTHOR]

Two methods for assimilating radar reflectivity into deterministic convection-allowing forecasts were compared: an operationally used, computationally less expensive cloud analysis (CA) scheme and a relatively more expensive, but rigorous, ensemble Kalman filter–variational hybrid method (EnVar). These methods were implemented in the Nonhydrostatic Multiscale Model on the B-grid and were tested on 10 cases featuring high-impact deep convective storms and heavy precipitation. A variety of traditional, neighborhood-based, and features-based verification metrics support that the EnVar produced superior free forecasts compared to the CA procedure, with statistically significant differences extending up to 9 h into the forecast. Despite being inferior, the CA scheme was able to provide benefit compared to not assimilating radar reflectivity at all, but limited to the first few forecast hours. While the EnVar is able to partially suppress spurious convection by assimilating 0-dBZ reflectivity observations directly, the CA is not designed to reduce or remove hydrometeors. As a result, the CA struggles more with suppression of spurious convection in the first-guess field, which resulted in high-frequency biases and poor forecast evolution, as illustrated in a few case studies. Additionally, while the EnVar uses flow-dependent ensemble covariances to update hydrometers, thermodynamic, and dynamic variables simultaneously when the reflectivity is assimilated, the CA relies on a radar reflectivity-derived latent heating rate that is applied during a separate digital filter initialization (DFI) procedure to introduce deep convective storms into the model, and the results of CA are shown to be sensitive to the window length used in the DFI. [ABSTRACT FROM AUTHOR]

Assimilating inner-core observations collected from recent field campaign programs such as Tropical Cyclone Intensity (TCI) and Intensity Forecasting Experiment (IFEX) together with the enhanced atmospheric motion vectors (AMVs) produce realistic three-dimensional (3D) analyses using the newly developed GSI-based, continuously cycled, dual-resolution hybrid ensemble–variational data assimilation (DA) system for the Hurricane Weather Research and Forecasting (HWRF) Model for Hurricane Patricia (2015). However, more persistent surface wind maximum spindown is found in the intensity forecast initialized from the realistic analyses produced by the DA system but not from the unrealistic initial conditions produced through vortex modification. Diagnostics in this study reveal that the spindown issue is likely attributed to the deficient HWRF Model physics that are unable to maintain the realistic 3D structures from the DA analysis. The horizontal diffusion is too strong to maintain the realistically observed vertical oscillation of radial wind near the eyewall region. The vertical diffusion profile cannot produce a sufficiently strong secondary circulation connecting the realistically elevated upper-level outflow produced in the DA analysis. Further investigations with different model physics parameterizations demonstrate that spindown can be alleviated by modifying model physics parameterizations. In particular, a modified turbulent mixing parameterization scheme together with a reduced horizontal diffusion is found to significantly alleviate the spindown issue and to improve the intensity forecast. Additional experiments show that the peak-simulated intensity and rapid intensification rate can be further improved by increasing the model resolution. But the model resolution is not as important as model physics in the spindown alleviation. [ABSTRACT FROM AUTHOR]

The Nationwide Network of Networks (NNoN) concept was introduced by the National Research Council to address the growing need for a national mesoscale observing system and the continued advancement toward accurate high-resolution numerical weather prediction. The research test bed known as the Dallas–Fort Worth (DFW) Urban Demonstration Network was created to experiment with many kinds of mesoscale observations that could be used in a data assimilation system. Many nonconventional observations, including Earth Networks and Citizen Weather Observer Program surface stations, are combined with conventional operational data to form the test bed network. A principal component of the NNoN effort is the quantification of observation impact from several different sources of information. In this study, the GSI-based EnKF system was used together with the WRF-ARW Model to examine impacts of observations assimilated for forecasting convection initiation (CI) in the 3 April 2014 hail storm case. Data denial experiments tested the impact of high-frequency (5 min) assimilation of nonconventional data on the timing and location of CI and subsequent storm evolution. Results showed nonconventional observations were necessary to capture details in the dryline structure causing localized enhanced convergence and leading to CI. Diagnosis of denial-minus-control fields showed the cumulative influence each observing network had on the resulting CI forecast. It was found that most of this impact came from the assimilation of thermodynamic observations in sensitive areas along the dryline gradient. Accurate metadata were found to be crucial toward the future application of nonconventional observations in high-resolution assimilation and forecast systems. [ABSTRACT FROM AUTHOR]

A prototype convection-allowing system using the Advanced Research version of the Weather Research and Forecasting (WRF-ARW) Model and employing an ensemble Kalman filter (EnKF) data assimilation technique has been developed and used during the spring 2016 and 2017 Hazardous Weather Testbeds. This system assimilates WSR-88D reflectivity and radial velocity, geostationary satellite cloud water path (CWP) retrievals, and available surface observations over a regional domain with a 3-km horizontal resolution at 15-min intervals, with 3-km initial conditions provided by an experimental High-Resolution Rapid Refresh ensemble (HRRR-e). However, no information on upper-level thermodynamic conditions in cloud-free regions is currently assimilated, as few timely observations exist. One potential solution is to also assimilate clear-sky satellite radiances, which provide information on mid- and upper-tropospheric temperature and moisture conditions. This research assimilates GOES-13 imager water vapor band (6.5 μm) radiances using the GSI-EnKF system to take advantage of the Community Radiative Transfer Model (CRTM) integration. Results using four cases from May 2016 showed that assimilating radiances generally had a neutral-to-positive impact on the model analysis, reducing humidity bias and/or errors at the appropriate model levels where verification observations were present. The effects on high-impact weather forecasts, as verified against forecast reflectivity and updraft helicity, were mixed. Three cases (9, 22, and 24 May) showed some improvement in skill, while the other (25 May) performed worse, despite the improved environment. This research represents the first step in designing a high-resolution ensemble data assimilation system to use GOES-16 Advanced Baseline Imager data, which provides additional water vapor bands and increased spatial and temporal resolution. [ABSTRACT FROM AUTHOR]

Multiscale ensemble-based data assimilation and forecasts were performed in real time during the Plains Elevated Convection At Night (PECAN) field experiment. A 20-member ensemble of forecasts at 4-km grid spacing was initialized daily at both 1300 and 1900 UTC, together with a deterministic forecast at 1-km grid spacing initialized at 1300 UTC. The configuration of the GSI-based data assimilation and forecast system was guided by results presented in Part I of this two-part study. The present paper describes the implementation of the real-time system and the extensive forecast products that were generated to support the unique interests of PECAN researchers. Subjective and objective verification of the real-time forecasts from 1 June through 15 July 2015 is conducted, with an emphasis on nocturnal mesoscale convective systems (MCSs), nocturnal convective initiation (CI), nocturnal low-level jets (LLJs), and bores on the nocturnal stable layer. Verification of nocturnal precipitation during overnight hours, a proxy for MCSs, shows both greater skill and spread for the 1300 UTC forecasts than the 1900 UTC forecasts. Verification against observed soundings reveals that the forecast LLJs systematically peak, veer, and dissipate several hours before the observations. Comparisons with bores that passed over an Atmospheric Emitted Radiance Interferometer reveal an ability to predict borelike features that is greatly improved at 1-km, compared with 4-km, grid spacing. Objective verification of forecast CI timing reveals strong sensitivity to the PBL scheme but an overall unbiased ensemble. [ABSTRACT FROM AUTHOR]

In this exploratory study, a series of perturbations to the land surface model (LSM) component of the Weather Research and Forecasting (WRF) Model was developed to investigate the sensitivity of forecasts of severe thunderstorms and heavy precipitation at 4-km grid spacing and whether such perturbations could improve ensemble forecasts at this scale. The perturbations (generated using a combination of perturbing fixed parameters and using separate schemes, one of which-Noah-MP-is new among the WRF modeling community) were applied to a 10-member ensemble including other mixed physics parameterizations and compared against an identically configured ensemble that did not include the LSM perturbations to determine their impact on probabilistic forecasts. A third ensemble using only the LSM perturbations was also configured. The results from 14 (in total) 36-h ensemble forecasts suggested the LSM perturbations resulted in systematic improvement in ensemble dispersion and error characteristics. Lower-tropospheric temperature, moisture, and wind fields were all improved, as were probabilistic precipitation forecasts. Biases were not systematically altered, although some outlier members are present. Examination of near-surface temperature and mixing ratio fields, surface energy fluxes, and soil fields revealed tendencies caused by certain perturbations. A case study featuring tornadic supercells illustrated the physical causes of some of these tendencies. The results of this study suggest LSM perturbations can sample a dimension of model error not yet sampled systematically in most ensembles and should be included in convection-allowing ensembles. [ABSTRACT FROM AUTHOR]

The impacts of multiscale flow-dependent initial condition (IC) perturbations for storm-scale ensemble forecasts of midlatitude convection are investigated using perfect-model observing system simulation experiments. Several diverse cases are used to quantitatively and qualitatively understand the impacts of different IC perturbations on ensemble forecast skill. Scale dependence of the results is assessed by evaluating 2-h storm-scale reflectivity forecasts separately from hourly accumulated mesoscale precipitation forecasts. Forecasts are initialized with different IC ensembles, including an ensemble of multiscale perturbations produced by a multiscale data assimilation system, mesoscale perturbations produced at a coarser resolution, and filtered multiscale perturbations. Mesoscale precipitation forecasts initialized with the multiscale perturbations are more skillful than the forecasts initialized with the mesoscale perturbations at several lead times. This multiscale advantage is due to greater consistency between the IC perturbations and IC uncertainty. This advantage also affects the short-term, smaller-scale forecasts. Reflectivity forecasts on very small scales and very short lead times are more skillful with the multiscale perturbations as a direct result of the smaller-scale IC perturbation energy. The small-scale IC perturbations also contribute to some improvements to the mesoscale precipitation forecasts after the ~5-h lead time. Altogether, these results suggest that the multiscale IC perturbations provided by ensemble data assimilation on the convection-permitting grid can improve storm-scale ensemble forecasts by improving the sampling of IC uncertainty, compared to downscaling of IC perturbations from a coarser-resolution ensemble. [ABSTRACT FROM AUTHOR]

The efficacy of a stochastic kinetic energy backscatter (SKEB) scheme to improve convection-allowing probabilistic forecasts was studied. While SKEB has been explored for coarse, convection-parameterizing models, studies of SKEB for convective scales are limited. Three ensembles were compared. The SKMP ensemble used mixed physics with the SKEB scheme, whereas the MP ensemble was configured identically but without using the SKEB scheme. The SK ensemble used the SKEB scheme with no physics diversity. The experiment covered May 2013 over the central United States on a 4-km Weather Research and Forecasting (WRF) Model domain. The SKEB scheme was successful in increasing the spread in all fields verified, especially mid- and upper-tropospheric fields. Additionally, the rmse of the ensemble mean was maintained or reduced, in some cases significantly. Rank histograms in the SKMP ensemble were flatter than those in the MP ensemble, indicating the SKEB scheme produces a less underdispersive forecast distribution. Some improvement was seen in probabilistic precipitation forecasts, particularly when examining Brier scores. Verification against surface observations agree with verification against Rapid Refresh (RAP) model analyses, showing that probabilistic forecasts for 2-m temperature, 2-m dewpoint, and 10-m winds were also improved using the SKEB scheme. The SK ensemble gave competitive forecasts for some fields. The SK ensemble had reduced spread compared to the MP ensemble at the surface due to the lack of physics diversity. These results suggest the potential utility of mixed physics plus the SKEB scheme in the design of convection-allowing ensemble forecasts. [ABSTRACT FROM AUTHOR]

A GSI-based data assimilation (DA) system, including three-dimensional variational assimilation (3DVar) and ensemble Kalman filter (EnKF), is extended to the multiscale assimilation of both meso- and synoptic-scale observation networks and convective-scale radar reflectivity and velocity observations. EnKF and 3DVar are systematically compared in this multiscale context to better understand the impacts of differences between the DA techniques on the analyses at multiple scales and the subsequent convective-scale precipitation forecasts. Averaged over 10 diverse cases, 8-h precipitation forecasts initialized using GSI-based EnKF are more skillful than those using GSI-based 3DVar, both with and without storm-scale radar DA. The advantage from radar DA persists for ~5 h using EnKF, but only ~1 h using 3DVar. A case study of an upscale growing MCS is also examined. The better EnKF-initialized forecast is attributed to more accurate analyses of both the mesoscale environment and the storm-scale features. The mesoscale location and structure of a warm front is more accurately analyzed using EnKF than 3DVar. Furthermore, storms in the EnKF multiscale analysis are maintained during the subsequent forecast period. However, storms in the 3DVar multiscale analysis are not maintained and generate excessive cold pools. Therefore, while the EnKF forecast with radar DA remains better than the forecast without radar DA throughout the forecast period, the 3DVar forecast quality is degraded by radar DA after the first hour. Diagnostics revealed that the inferior analysis at mesoscales and storm scales for the 3DVar is primarily attributed to the lack of flow dependence and cross-variable correlation, respectively, in the 3DVar static background error covariance. [ABSTRACT FROM AUTHOR]

Two approaches for accounting for errors in quantitative precipitation forecasts (QPFs) due to uncertainty in the microphysics (MP) parameterization in a convection-allowing ensemble are examined. They include mixed MP (MMP) composed mostly of double-moment schemes and perturbing parameters within the Weather Research and Forecasting single-moment 6-class microphysics scheme (WSM6) MP scheme (PPMP). Thirty-five cases of real-time storm-scale ensemble forecasts produced by the Center for Analysis and Prediction of Storms during the NOAA Hazardous Weather Testbed 2011 Spring Experiment were examined. The MMP ensemble had better fractions Brier scores (FBSs) for most lead times and thresholds, but the PPMP ensemble had better relative operating characteristic (ROC) scores for higher precipitation thresholds. The pooled ensemble formed by randomly drawing five members from the MMP and PPMP ensembles was no more skillful than the more accurate of the MMP and PPMP ensembles. Significant positive impact was found when the two were combined to form a larger ensemble. The QPF and the systematic behaviors of derived microphysical variables were also examined. The skill of the QPF among different members depended on the thresholds, verification metrics, and forecast lead times. The profiles of microphysics variables from the double-moment schemes contained more variation in the vertical than those from the single-moment members. Among the double-moment schemes, WDM6 produced the smallest raindrops and very large number concentrations. Among the PPMP members, the behaviors were found to be consistent with the prescribed intercept parameters. The perturbed intercept parameters used in the PPMP ensemble fell within the range of values retrieved from the double-moment schemes. [ABSTRACT FROM AUTHOR]

Multiscale convection-allowing precipitation forecast perturbations are examined for two forecasts and systematically over 34 forecasts out to 30-h lead time using Haar Wavelet decomposition. Two small-scale initial condition (IC) perturbation methods are compared to the larger-scale IC and physics perturbations in an experimental convection-allowing ensemble. For a precipitation forecast driven primarily by a synoptic-scale baroclinic disturbance, small-scale IC perturbations resulted in little precipitation forecast perturbation energy on medium and large scales, compared to larger-scale IC and physics (LGPH) perturbations after the first few forecast hours. However, for a case where forecast convection at the initial time grew upscale into a mesoscale convective system (MCS), small-scale IC and LGPH perturbations resulted in similar forecast perturbation energy on all scales after about 12 h. Small-scale IC perturbations added to LGPH increased total forecast perturbation energy for this case. Averaged over 34 forecasts, the small-scale IC perturbations had little impact on large forecast scales while LGPH accounted for about half of the error energy on such scales. The impact of small-scale IC perturbations was also less than, but comparable to, the impact of LGPH perturbations on medium scales. On small scales, the impact of small-scale IC perturbations was at least as large as the LGPH perturbations. The spatial structure of small-scale IC perturbations affected the evolution of forecast perturbations, especially at medium scales. There was little systematic impact of the small-scale IC perturbations when added to LGPH. These results motivate further studies on properly sampling multiscale IC errors. [ABSTRACT FROM AUTHOR]

Forecasts generated by the Center for Analysis and Prediction of Storms with 1- and 4-km grid spacing using the Advanced Research Weather Research and Forecasting Model (ARW-WRF; ARW1 and ARW4, respectively) for the 2009-11 NOAA Hazardous Weather Testbed Spring Experiments are compared and verified. Object-based measures, including average values of object attributes, the object-based threat score (OTS), and the median of maximum interest (MMI) are used for the verification. Verification was first performed against observations at scales resolvable by each forecast model and then performed at scales resolvable by both models by remapping ARW1 to the ARW4 grid (ARW1to4). Thirty-hour forecasts of 1-h accumulated precipitation initialized at 0000 UTC on 22, 36, and 33 days during the spring of 2009, 2010, and 2011, respectively, are evaluated over a domain covering most of the central and eastern United States. ARW1, ARW1to4, and ARW4 all significantly overforecasted the number of objects during diurnal convection maxima. The overforecasts by ARW1 and ARW1to4 were more pronounced than ARW4 during the first convection maximum at 1-h lead time. The average object area and aspect ratio were closer to observations for ARW1 and ARW1to4 than for ARW4. None of the models showed a significant advantage over the others for average orientation angle and centroid location. Increased accuracy for ARW1, compared to ARW4, was statistically significant for the MMI but not the OTS. However, ARW1to4 had similar MMI and OTS as ARW4 at most lead times. These results are consistent with subjective evaluations that the greatest impact of grid spacing is on the smallest resolvable objects. [ABSTRACT FROM AUTHOR]

Convection intensity and longevity is highly dependent on the surrounding environment. Ensemble sensitivity analysis (ESA), which quantitatively and qualitatively interprets impacts of initial conditions on forecasts, is applied to very short-term (1–2 h) convective-scale forecasts for three cases during the Mesoscale Predictability Experiment (MPEX) in 2013. The ESA technique reveals several dependencies of individual convective storm evolution on their nearby environments. The three MPEX cases are simulated using a previously verified 36-member convection-allowing model (Δx = 3 km) ensemble created via the Weather Research and Forecasting (WRF) Model. Radar and other conventional observations are assimilated using an ensemble adjustment Kalman filter. The three cases include a mesoscale convective system (MCS) and both nontornadic and tornadic supercells. Of the many ESAs applied in this study, one of the most notable is the positive sensitivity of supercell updraft helicity to increases in both storm inflow region deep and shallow vertical wind shear. This result suggests that larger values of vertical wind shear within the storm inflow yield higher values of storm updraft helicity. Results further show that the supercell storms quickly enhance the environmental vertical wind shear within the storm inflow region. Application of ESA shows that these storm-induced perturbations then affect further storm evolution, suggesting the presence of storm–environment feedback cycles where perturbations affect future mesocyclone strength. Overall, ESA can provide insight into convection dependencies on the near-storm environment. [ABSTRACT FROM AUTHOR]

This article summarizes research and risk reduction that will inform acquisition decisions regarding NOAA's future national operational weather radar network. A key alternative being evaluated is polarimetric phased-array radar (PAR). Research indicates PAR can plausibly achieve fast, adaptive volumetric scanning, with associated benefits for severe-weather warning performance. We assess these benefits using storm observations and analyses, observing system simulation experiments, and real radar-data assimilation studies. Changes in the number and/or locations of radars in the future network could improve coverage at low altitude. Analysis of benefits that might be so realized indicates the possibility for additional improvement in severe-weather and flash-flood warning performance, with associated reduction in casualties. Simulations are used to evaluate techniques for rapid volumetric scanning and assess data quality characteristics of PAR. Finally, we describe progress in developing methods to compensate for polarimetric variable estimate biases introduced by electronic beam-steering. A research-to-operations (R2O) strategy for the PAR alternative for the WSR-88D replacement network is presented. [ABSTRACT FROM AUTHOR]

The increasing maturity of the Warn-on-Forecast System (WoFS) coupled with the now operational GOES-16 satellite allows for the first time a comprehensive analysis of the relative impacts of assimilating GOES-16 all-sky 6.2-, 6.9-, and 7.3-μm channel radiances compared to other radar and satellite observations. The WoFS relies on cloud property retrievals such as cloud water path, which have been proven to increase forecast skill compared to only assimilating radar data and other conventional observations. The impacts of assimilating clear-sky radiances have also been explored and shown to provide useful information on midtropospheric moisture content in the near-storm environment. Assimilation of all-sky radiances adds a layer of complexity and is tested to determine its effectiveness across four events occurring in the spring and summer of 2019. Qualitative and object-based verification of severe weather and the near-storm environment are used to assess the impact of assimilating all-sky radiances compared to the current model configuration. We focus our study through the entire WoFS analysis and forecasting cycle (1900–0600 UTC, daily) so that the impacts throughout the evolution of convection from initiation to large upscale growth can be assessed. Overall, assimilating satellite data improves forecasts relative to radar-only assimilation experiments. The retrieval method with clear-sky radiances performs best overall, but assimilating all-sky radiances does have very positive impacts in certain conditions. In particular, all-sky radiance assimilation improved convective initiation forecast of severe storms in several instances. This work represents an initial attempt at assimilating all-sky radiances into the WoFS and additional research is ongoing to further improve forecast skill. [ABSTRACT FROM AUTHOR]

The 2016–18 NOAA Hazardous Weather Testbed (HWT) Spring Forecasting Experiments (SFE) featured the Community Leveraged Unified Ensemble (CLUE), a coordinated convection-allowing model (CAM) ensemble framework designed to provide empirical guidance for development of operational CAM systems. The 2017 CLUE included 81 members that all used 3-km horizontal grid spacing over the CONUS, enabling direct comparison of forecasts generated using different dynamical cores, physics schemes, and initialization procedures. This study uses forecasts from several of the 2017 CLUE members and one operational model to evaluate and compare CAM representation and next-day prediction of thunderstorms. The analysis utilizes existing techniques and novel, object-based techniques that distill important information about modeled and observed storms from many cases. The National Severe Storms Laboratory Multi-Radar Multi-Sensor product suite is used to verify model forecasts and climatologies of observed variables. Unobserved model fields are also examined to further illuminate important intermodel differences in storms and near-storm environments. No single model performed better than the others in all respects. However, there were many systematic intermodel and intercore differences in specific forecast metrics and model fields. Some of these differences can be confidently attributed to particular differences in model design. Model intercomparison studies similar to the one presented here are important to better understand the impacts of model and ensemble configurations on storm forecasts and to help optimize future operational CAM systems. [ABSTRACT FROM AUTHOR]

Due to lack of high spatial and temporal resolution boundary layer (BL) observations, the rapid changes in the near-storm environment are not well represented in current convective-scale numerical models. Better representation of the near-storm environment in model initial conditions will likely further improve the forecasts of severe convective weather. This study investigates the impact of assimilating high temporal resolution BL retrievals from two ground-based remote sensing instruments for short-term forecasts of a tornadic supercell event on 13 July 2015 during the Plains Elevated Convection At Night field campaign. The instruments are the Atmospheric Emitted Radiance Interferometer (AERI) that retrieves thermodynamic profiles and the Doppler lidar (DL) that measures horizontal wind profiles. Six sets of convective-scale ensemble data assimilation (DA) experiments are performed: two control experiments that assimilate conventional and WSR-88D radar observations using either relaxation-to-prior-spread (RTPS) or the adaptive inflation (AI) technique and four experiments similar to the control but that assimilate either DL or AERI or both observations in addition to all other observations that are in the control experiments. Results indicate a positive impact of AERI and DL observations in forecasting convective initiation (CI) and early evolution of the supercell storm. The experiment that employs the AI technique to assimilate BL observations in DA enhances the humidity in the near-storm environment and low-level convergence, which in turn helps forecasting CI. The forecast improvement is most pronounced during the first ~3 h. Results also indicate that the AERI observations have a larger impact compared to DL in predicting CI. [ABSTRACT FROM AUTHOR]

Sampling errors and model errors are major drawbacks from which ensemble Kalman filters suffer. Sampling errors arise because of the use of a limited ensemble size, while model errors are deficiencies in the dynamics and underlying parameterizations that may yield biases in the model's prediction. In this study, we propose a new time-adaptive posterior inflation algorithm in which the analyzed ensemble anomalies are locally inflated. The proposed inflation strategy is computationally efficient and is aimed at restoring enough spread in the analysis ensemble after assimilating the observations. The performance of this scheme is tested against the relaxation to prior spread (RTPS) and adaptive prior inflation. For this purpose, two model are used: the three-variable Lorenz 63 system and the Community Atmosphere Model (CAM). In CAM, global refractivity, temperature, and wind observations from several sources are incorporated to perform a set of assimilation experiments using the Data Assimilation Research Testbed (DART). The proposed scheme is shown to yield better quality forecasts than the RTPS. Assimilation results further suggest that when model errors are small, both prior and posterior inflation are able to mitigate sampling errors with a slight advantage to posterior inflation. When large model errors, such as wind and temperature biases, are present, prior inflation is shown to be more accurate than posterior inflation. Densely observed regions as in the Northern Hemisphere present numerous challenges to the posterior inflation algorithm. A compelling enhancement to the performance of the filter is achieved by combining both adaptive inflation schemes. [ABSTRACT FROM AUTHOR]

In this study, the ensemble of three-dimensional variational data assimilation (En3DVar) method for convective-scale weather is adopted and evaluated using an idealized supercell storm simulated by the Weather Research and Forecasting (WRF) Model. Synthetic radar radial velocity, reflectivity, satellitederived cloud water path (CWP), and total precipitable water (TPW) data are produced from the simulated supercell storm and then these data are assimilated into another WRF Model run that starts with no convection. Two types of experiments are performed. The first assimilates radar and satellite CWP data using a perfect storm environment. The second assimilates additional TPW data using a storm environment with dry bias. The first set of experiments indicates that incorporating CWP and radar data into the assimilation leads to a much faster initiation of supercell storms than found using radar data alone. Assimilating CWP data primarily improves the analyses of nonprecipitating hydrometeor variables. The results from the second set of experiments demonstrate the critical importance of the storm environment. When using the biased storm environment, assimilation of CWP and radar data enhances the analyses, but the forecast skill rapidly decreases over the subsequent 1-h forecast. Further experiments show that assimilating theTPWdata has a large impact on storm environment that is essential to the accuracy of the storm forecasts. In general, the combination of radar data and satellite data within the En3DVar results in better analyses and forecasts than when only radar data are used, especially for an imperfect storm environment. [ABSTRACT FROM AUTHOR]

A coupled ensemble square root filter-three-dimensional ensemble-variational hybrid (EnSRF-En3DVar) data assimilation (DA) system is developed for the operational Rapid Refresh (RAP) forecasting system. The En3DVar hybrid system employs the extended control variable method, and is built on the NCEP operational gridpoint statistical interpolation (GSI) three-dimensional variational data assimilation (3DVar) framework. It is coupled with an EnSRF system for RAP, which provides ensemble perturbations. Recursive filters (RF) are used to localize ensemble covariance in both horizontal and vertical within the En3DVar. The coupled En3DVar hybrid system is evaluated with 3-h cycles over a 9-day period with active convection. All conventional observations used by operational RAP are included. The En3DVar hybrid system is run at ⅓ of the operational RAP horizontal resolution or about 40-km grid spacing, and its performance is compared to parallel GSI 3DVar and EnSRF runs using the same datasets and resolution. Short-term forecasts initialized from the 3-hourly analyses are verified against sounding and surface observations. When using equally weighted static and ensemble background error covariances and 40 ensemble members, the En3DVar hybrid system outperforms the corresponding GSI 3DVar and EnSRF. When the recursive filter coefficients are tuned to achieve a similar height-dependent localization as in the EnSRF, the En3DVar results using pure ensemble covariance are close to EnSRF. Two-way coupling between EnSRF and En3DVar did not produce noticeable improvement over one-way coupling. Downscaled precipitation forecast skill on the 13-km RAP grid from the En3DVar hybrid is better than those from GSI 3DVar analyses. [ABSTRACT FROM AUTHOR]

A regional ensemble Kalman filter (EnKF) system is established for potential Rapid Refresh (RAP) operational application. The system borrows data processing and observation operators from the gridpoint statistical interpolation (GSI), and precalculates observation priors using the GSI. The ensemble square root Kalman filter (EnSRF) algorithm is used, which updates both the state vector and observation priors. All conventional observations that are used in the operational RAP GSI are assimilated. To minimize computational costs, the EnKF is run at ⅓ of the operational RAP resolution or about 40-km grid spacing, and its performance is compared to the GSI using the same datasets and resolution. Short-range (up to 18 h, the RAP forecast length) forecasts are verified against soundings, surface observations, and precipitation data. Experiments are run with 3-hourly assimilation cycles over a 9-day convectively active retrospective period from spring 2010. The EnKF performance was improved by extensive tuning, including the use of height-dependent covariance localization scales and adaptive covariance inflation. When multiple physics parameterization schemes are employed by the EnKF, forecast errors are further reduced, especially for relative humidity and temperature at the upper levels and for surface variables. The best EnKF configuration produces lower forecast errors than the parallel GSI run. Gilbert skill scores of precipitation forecasts on the 13-km RAP grid initialized from the 3-hourly EnKF analyses are consistently better than those from GSI analyses. [ABSTRACT FROM AUTHOR]