Ensemble Kalman Filter Assimilation of Simulated HIWRAP Doppler Velocity Data in a Hurricane

2013 ◽  
Vol 141 (8) ◽  
pp. 2683-2704 ◽  
Author(s):  
Jason A. Sippel ◽  
Scott A. Braun ◽  
Fuqing Zhang ◽  
Yonghui Weng

Abstract This study utilizes ensemble Kalman filter (EnKF) observing system simulation experiments (OSSEs) to analyze the potential impact of assimilating radial velocity observations of hurricanes from the High-altitude Imaging Wind and Rain Airborne Profiler (HIWRAP). HIWRAP is a new Doppler radar mounted on the NASA Global Hawk unmanned airborne system that flies at roughly 19-km altitude and has the benefit of a 25–30-h flight duration, which is 2–3 times that of conventional aircraft. This research is intended as a proof-of-concept study for future assimilation of real HIWRAP data. The most important result from this research is that HIWRAP data can potentially improve hurricane analyses and prediction. For example, by the end of a 12-h assimilation period, the analysis error is much lower than that in deterministic forecasts. As a result, subsequent forecasts initialized with the EnKF analyses also improve. Furthermore, analyses and forecasts clearly benefit more from a 12-h assimilation period than for shorter periods, which highlights a benefit of the Global Hawk's potentially long on-station times.

2014 ◽  
Vol 142 (12) ◽  
pp. 4559-4580 ◽  
Author(s):  
Jason A. Sippel ◽  
Fuqing Zhang ◽  
Yonghui Weng ◽  
Lin Tian ◽  
Gerald M. Heymsfield ◽  
...  

Abstract This study utilizes an ensemble Kalman filter (EnKF) to assess the impact of assimilating observations of Hurricane Karl from the High-Altitude Imaging Wind and Rain Airborne Profiler (HIWRAP). HIWRAP is a new Doppler radar on board the NASA Global Hawk unmanned airborne system, which has the benefit of a 24–26-h flight duration, or about 2–3 times that of a conventional aircraft. The first HIWRAP observations were taken during NASA’s Genesis and Rapid Intensification Processes (GRIP) experiment in 2010. Observations considered here are Doppler velocity (Vr) and Doppler-derived velocity–azimuth display (VAD) wind profiles (VWPs). Karl is the only hurricane to date for which HIWRAP data are available. Assimilation of either Vr or VWPs has a significant positive impact on the EnKF analyses and forecasts of Hurricane Karl. Analyses are able to accurately estimate Karl’s observed location, maximum intensity, size, precipitation distribution, and vertical structure. In addition, forecasts initialized from the EnKF analyses are much more accurate than a forecast without assimilation. The forecasts initialized from VWP-assimilating analyses perform slightly better than those initialized from Vr-assimilating analyses, and the latter are less accurate than EnKF-initialized forecasts from a recent proof-of-concept study with simulated data. Likely causes for this discrepancy include the quality and coverage of the HIWRAP data collected from Karl and the presence of model error in this real-data situation. The advantages of assimilating VWP data likely include the ability to simultaneously constrain both components of the horizontal wind and to circumvent reliance upon vertical velocity error covariance.


Author(s):  
Frédéric Fabry

Abstract In the ensemble Kalman filter (EnKF), the covariance localization radius is usually small when assimilating radar observations because of high density of the radar observations. This makes the region away from precipitation difficult to correct if no other observations are available, as there is no reason to correct the background. To correct errors away from the innovating radar observations, a multiscale localization (MLoc) method adapted to dense observations like those from radar is proposed. In this method, different scales are corrected successively by using the same reflectivity observations, but with different degree of smoothing and localization radius at each step. In the context of observing system simulation experiments, single and multiple assimilation experiments are conducted with the MLoc method. Results show that the MLoc assimilation updates areas that are away from the innovative observations and improves on average the analysis and forecast quality in single cycle and cycling assimilation experiments. The forecast gains are maintained until the end of the forecast period, illustrating the benefits of correcting different scales.


Icarus ◽  
2010 ◽  
Vol 209 (2) ◽  
pp. 470-481 ◽  
Author(s):  
Matthew J. Hoffman ◽  
Steven J. Greybush ◽  
R. John Wilson ◽  
Gyorgyi Gyarmati ◽  
Ross N. Hoffman ◽  
...  

2011 ◽  
Vol 139 (11) ◽  
pp. 3389-3404 ◽  
Author(s):  
Thomas Milewski ◽  
Michel S. Bourqui

Abstract A new stratospheric chemical–dynamical data assimilation system was developed, based upon an ensemble Kalman filter coupled with a Chemistry–Climate Model [i.e., the intermediate-complexity general circulation model Fast Stratospheric Ozone Chemistry (IGCM-FASTOC)], with the aim to explore the potential of chemical–dynamical coupling in stratospheric data assimilation. The system is introduced here in a context of a perfect-model, Observing System Simulation Experiment. The system is found to be sensitive to localization parameters, and in the case of temperature (ozone), assimilation yields its best performance with horizontal and vertical decorrelation lengths of 14 000 km (5600 km) and 70 km (14 km). With these localization parameters, the observation space background-error covariance matrix is underinflated by only 5.9% (overinflated by 2.1%) and the observation-error covariance matrix by only 1.6% (0.5%), which makes artificial inflation unnecessary. Using optimal localization parameters, the skills of the system in constraining the ensemble-average analysis error with respect to the true state is tested when assimilating synthetic Michelson Interferometer for Passive Atmospheric Sounding (MIPAS) retrievals of temperature alone and ozone alone. It is found that in most cases background-error covariances produced from ensemble statistics are able to usefully propagate information from the observed variable to other ones. Chemical–dynamical covariances, and in particular ozone–wind covariances, are essential in constraining the dynamical fields when assimilating ozone only, as the radiation in the stratosphere is too slow to transfer ozone analysis increments to the temperature field over the 24-h forecast window. Conversely, when assimilating temperature, the chemical–dynamical covariances are also found to help constrain the ozone field, though to a much lower extent. The uncertainty in forecast/analysis, as defined by the variability in the ensemble, is large compared to the analysis error, which likely indicates some amount of noise in the covariance terms, while also reducing the risk of filter divergence.


2010 ◽  
Vol 138 (2) ◽  
pp. 517-538 ◽  
Author(s):  
Nusrat Yussouf ◽  
David J. Stensrud

Abstract The conventional Weather Surveillance Radar-1988 Doppler (WSR-88D) scans a given weather phenomenon in approximately 5 min, and past results suggest that it takes 30–60 min to establish a storm into a model assimilating these data using an ensemble Kalman filter (EnKF) data assimilation technique. Severe-weather events, however, can develop and evolve very rapidly. Therefore, assimilating observations for a 30–60-min period prior to the availability of accurate analyses may not be feasible in an operational setting. A shorter assimilation period also is desired if forecasts are produced to increase the warning lead time. With the advent of the emerging phased-array radar (PAR) technology, it is now possible to scan the same weather phenomenon in less than 1 min. Therefore, it is of interest to see if the faster scanning rate of PAR can yield improvements in storm-scale analyses and forecasts from assimilating over a shorter period of time. Observing system simulation experiments are conducted to evaluate the ability to quickly initialize a storm into a numerical model using PAR data in place of WSR-88D data. Synthetic PAR and WSR-88D observations of a splitting supercell storm are created from a storm-scale model run using a realistic volume-averaging technique in native radar coordinates. These synthetic reflectivity and radial velocity observations are assimilated into the same storm-scale model over a 15-min period using an EnKF data assimilation technique followed by a 50-min ensemble forecast. Results indicate that assimilating PAR observations at 1-min intervals over a short 15-min period yields significantly better analyses and ensemble forecasts than those produced using WSR-88D observations. Additional experiments are conducted in which the adaptive scanning capability of PAR is utilized for thunderstorms that are either very close to or far away from the radar location. Results show that the adaptive scanning capability improves the analyses and forecasts when compared with the nonadaptive PAR data. These results highlight the potential for flexible rapid-scanning PAR observations to help to quickly and accurately initialize storms into numerical models yielding improved storm-scale analyses and very short range forecasts.


2017 ◽  
Vol 24 (3) ◽  
pp. 329-341 ◽  
Author(s):  
Guocan Wu ◽  
Xiaogu Zheng

Abstract. The ensemble Kalman filter (EnKF) is a widely used ensemble-based assimilation method, which estimates the forecast error covariance matrix using a Monte Carlo approach that involves an ensemble of short-term forecasts. While the accuracy of the forecast error covariance matrix is crucial for achieving accurate forecasts, the estimate given by the EnKF needs to be improved using inflation techniques. Otherwise, the sampling covariance matrix of perturbed forecast states will underestimate the true forecast error covariance matrix because of the limited ensemble size and large model errors, which may eventually result in the divergence of the filter. In this study, the forecast error covariance inflation factor is estimated using a generalized cross-validation technique. The improved EnKF assimilation scheme is tested on the atmosphere-like Lorenz-96 model with spatially correlated observations, and is shown to reduce the analysis error and increase its sensitivity to the observations.


2006 ◽  
Vol 134 (2) ◽  
pp. 618-637 ◽  
Author(s):  
Martin Charron ◽  
P. L. Houtekamer ◽  
Peter Bartello

Abstract The ensemble Kalman filter (EnKF) developed at the Meteorological Research Branch of Canada is used in the context of synthetic radial wind data assimilation at the mesoscale. A dry Boussinesq model with periodic boundary conditions is employed to provide a control run, as well as two ensembles of first guesses. Synthetic data, which are interpolated from the control run, are assimilated and simulate Doppler radar wind measurements. Nine “radars” with a range of 120 km are placed evenly on the horizontal 1000 km × 1000 km domain. These radars measure the radial wind with assumed Gaussian error statistics at each grid point within their range provided that there is sufficient upward motion (a proxy for precipitation). These data of radial winds are assimilated every 30 min and the assimilation period extends over 4 days. Results show that the EnKF technique with 2 × 50 members performed well in terms of reducing the analysis error for horizontal winds and temperature (even though temperature is not an observed variable) over a period of 4 days. However the analyzed vertical velocity shows an initial degradation. During the first 2 days of the assimilation period, the analysis error of the vertical velocity is greater when assimilating radar observations than when scoring forecasts initialized at t = 0 without assimilating any data. The type of assimilated data as well as the localization of the impact of the observations is thought to be the cause of this degradation of the analyzed vertical velocity. External gravity modes are present in the increments when localization is performed. This degradation can be eliminated by filtering the external gravity modes of the analysis increments. A similar set of experiments is realized in which the model dissipation coefficient is reduced by a factor of 10. This shows the level of sensitivity of the results to the kinetic energy power spectrum, and that the quality of the analyzed vertical wind is worse when dissipation is small.


Sign in / Sign up

Export Citation Format

Share Document