A Doppler radar data assimilation system is developed based on ensemble Kalman filter (EnKF) method and tested with simulated radar data from a supercell storm. As a first implementation, we assume the forward models are perfect and radar data are sampled at the analysis grid points. A general purpose nonhydrostatic compressible model is used with the inclusion of complex multi-class ice microphysics. New aspects compared to previous studies include the demonstration of the ability of EnKF method in retrieving multiple microphysical species associated with a multi-class ice microphysics scheme, and in accurately retrieving the wind and thermodynamic variables. Also new are the inclusion of reflectivity observations and the determination of the relative role of radial velocity and reflectivity data as well as their spatial coverage in recovering the full flow and cloud fields. In general, the system is able to reestablish the model storm extremely well after a number of assimilation cycles, and best results are obtained when both radial velocity and reflectivity data, including reflectivity information outside precipitation regions, are used. Significant positive impact of the reflectivity assimilation

The possibility of estimating fundamental parameters common in single-moment ice microphysics schemes using radar observations is investigated for a model-simulated supercell storm by examining parameter sensitivity and identifiability. These parameters include the intercept parameters for rain, snow, and hail/graupel, and the bulk densities of snow and hail/graupel. These parameters are closely involved in the definition of drop/particle size distributions of microphysical species but often assume highly uncertain specified values. The sensitivity of model forecast within data assimilation cycles to the parameter values, and the issue of solution uniqueness of the estimation problem, are examined. The ensemble square root filter (EnSRF) is employed for model state estimation. Sensitivity experiments show that the errors in the microphysical parameters have a larger impact on model microphysical fields than on wind fields; radar reflectivity observations are therefore preferred over those of radial velocity for microphysical parameter estimation. The model response time to errors in individual parameters are also investigated. The results suggest that radar data should be used at about 5-min intervals for parameter estimation. The response functions calculated from ensemble mean forecasts for all five individual parameters show concave shapes, with unique

In Part I of this two-part work, the feasibility of using an ensemble Kalman filter (EnKF) for mesoscale and regional-scale data assimilation through various observing system simulation experiments was demon-strated assuming a perfect forecast model for a winter snowstorm event that occurred on 24–26 January 2000. The current study seeks to explore the performance of the EnKF for the same event in the presence of significant model errors due to physical parameterizations by assimilating synthetic sounding and surface observations with typical temporal and spatial resolutions. The EnKF performance with imperfect models is also examined for a warm-season mesoscale convective vortex (MCV) event that occurred on 10–13 June 2003. The significance of model error in both warm- and cold-season events is demonstrated when the use of different cumulus parameterization schemes within different ensembles results in significantly different forecasts in terms of both ensemble mean and spread. Nevertheless, the EnKF performed reasonably well in most experiments with the imperfect model assumption (though its performance can sometimes be significantly degraded). As in Part I, where the perfect model assumption was utilized, most analysis error reduction comes from larger scales. Results show that using a combination of different physical parameter-ization schemes in the ensemble forecast can significantly improve filter performance. A multischeme

Abstract not found

A hybrid analysis scheme is compared with an ensemble square root filter (EnSRF) analysis scheme in the presence of model errors as a follow-up to a previous perfect-model comparison. In the hybrid scheme, the ensemble perturbations are updated by the ensemble transform Kalman filter (ETKF) and the ensemble mean is updated with a hybrid ensemble and static background-error covariance. The experiments were conducted with a two-layer primitive equation model. The true state was a T127 simulation. Data assimilation experiments were conducted at T31 resolution (3168 complex spectral coefficients), assimilating imperfect observations drawn from the T127 nature run. By design, the magnitude of the truncation error was large, which provided a test on the ability of both schemes to deal with model error. Additive noise was used to parameterize model errors in the background ensemble for both schemes. In the first set of experiments, additive noise was drawn from a large inventory of historical forecast errors; in the second set of experiments, additive noise was drawn from a large inventory of differences between forecasts and analyses. The static covariance was computed correspondingly from the two inventories. The hybrid analysis was statistically significantly more accurate than the EnSRF analysis. The improvement of the hybrid over the EnSRF was smaller when differences of forecasts and analyses were used to form the random noise and the static co-variance. The EnSRF analysis was more sensitive to the size of the ensemble than the hybrid. A series of tests was conducted to understand why the EnSRF performed worse than the hybrid. It was shown that the inferior performance of the EnSRF was likely due to the sampling error in the estimation of the model-error co-variance in the mean update and the less-balanced EnSRF initial conditions resulting from the extra locali-zations used in the EnSRF. 1.

A hybrid parallel scheme for the ensemble square root filter (EnSRF) suitable for parallel assimilation of multiscale observations, including those from dense observational networks such as those of radar, is de-veloped based on the domain decomposition strategy. The scheme handles internode communication through a message passing interface (MPI) and the communication within shared-memory nodes via Open Multi-processing (OpenMP) threads. It also supports pure MPI and pure OpenMPmodes. The parallel framework can accommodate high-volume remote-sensed radar (or satellite) observations as well as conventional ob-servations that usually have larger covariance localization radii. The performance of the parallel algorithm has been tested with simulated and real radar data. The parallel program shows good scalability in pure MPI and hybrid MPI–OpenMP modes, while pure OpenMP runs exhibit limited scalability on a symmetric shared-memory system. It is found that inMPImode, better parallel performance is achieved with domain decomposition configurations in which the leading dimension of the state variable arrays is larger, because this configuration allows for more efficient memory access. Given a fixed amount of computing resources, the hybrid parallel mode is preferred to pure MPI mode on super-computers with nodes containing shared-memory cores. The overall performance is also affected by factors such as the cache size, memory bandwidth, and the networking topology. Tests with a real data case with a large number of radars confirm that the parallel data assimilation can be done on amulticore supercomputer with a significant speedup compared to the serial data assimilation algorithm. 1.

One of the long-standing problems in meteorology is the estimation of the three-dimensional wind and

A framework for Observing System Simulation Experiments (OSSEs) based on the ensemble square root Kalman filter (EnSRF) technique for assimilating data from more than one radar network is described. The system is tested by assimilating simulated radial velocity and reflectivity data from a Weather Surveillance Radar-1988 Doppler (WSR-88D) radar and a network of four low-cost radars planned for the Oklahoma test bed by the new National Science Foundation (NSF) Engineering Research Center for Collaborative Adaptive Sensing of the Atmosphere (CASA). Such networks are meant to adaptively probe the lower atmosphere that is often missed by the existing WSR-88D radar network, so as to improve the detection of low-level hazardous weather events and to provide more complete data for the initialization of numerical weather prediction models. Different from earlier OSSE work with ensemble Kalman filters, the radar data are sampled on the radar elevation levels and a more realistic forward operator based on the Gaussian power-gain function is used. A stretched vertical grid with high vertical resolution near the ground allows for a better examination of the impact of low-level data. Furthermore, the impacts of storm propagation and higher-volume scan frequen-cies up to one volume scan per minute on the quality of analysis are examined, using a domain of a sufficient

A regional ensemble Kalman filter (EnKF) system is established for potential Rapid Refresh (RAP) op-erational 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 compu-tational costs, the EnKF is run at 1/3 of the operational RAP resolution or about 40-km grid spacing, and its performance is compared to theGSI 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. Ex-periments are run with 3-hourly assimilation cycles over a 9-day convectively active retrospective period from spring 2010. TheEnKFperformancewas 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-kmRAP grid

Warnings about convective-scale hazards are currently based on observations, but the time has come to develop warning methods in which numerical model forecasts play a much larger role. The National Oceanic and Atmospheric Administration’s (NOAA’s) National Weather Service (NWS) issues warnings when there is a threat to life and property from weather events. A warning is an urgent call for the public to take action when a hazardous weather or hydrologic event is occurring,