Search Results
You are looking at 1 - 10 of 14 items for
- Author or Editor: Zhiquan Liu x
- Refine by Access: All Content x
Abstract
Analyses with 20-km horizontal grid spacing were produced from parallel continuously cycling three-dimensional variational (3DVAR), ensemble square root Kalman filter (EnSRF), and “hybrid” variational–ensemble data assimilation (DA) systems between 0000 UTC 6 May and 0000 UTC 21 June 2011 over a domain spanning the contiguous United States. Beginning 9 May, the 0000 UTC analyses initialized 36-h Weather Research and Forecasting Model (WRF) forecasts containing a large convection-permitting 4-km nest. These 4-km 3DVAR-, EnSRF-, and hybrid-initialized forecasts were compared to benchmark WRF forecasts initialized by interpolating 0000 UTC Global Forecast System (GFS) analyses onto the computational domain.
While important differences regarding mean state characteristics of the 20-km DA systems were noted, verification efforts focused on the 4-km precipitation forecasts. The 3DVAR-, hybrid-, and EnSRF-initialized 4-km precipitation forecasts performed similarly regarding general precipitation characteristics, such as timing of the diurnal cycle, and all three forecast sets had high precipitation biases at heavier rainfall rates. However, meaningful differences emerged regarding precipitation placement as quantified by the fractions skill score. For most forecast hours, the hybrid-initialized 4-km precipitation forecasts were better than the EnSRF-, 3DVAR-, and GFS-initialized forecasts, and the improvement was often statistically significant at the 95th percentile. These results demonstrate the potential of limited-area continuously cycling hybrid DA configurations and suggest additional hybrid development is warranted.
Abstract
Analyses with 20-km horizontal grid spacing were produced from parallel continuously cycling three-dimensional variational (3DVAR), ensemble square root Kalman filter (EnSRF), and “hybrid” variational–ensemble data assimilation (DA) systems between 0000 UTC 6 May and 0000 UTC 21 June 2011 over a domain spanning the contiguous United States. Beginning 9 May, the 0000 UTC analyses initialized 36-h Weather Research and Forecasting Model (WRF) forecasts containing a large convection-permitting 4-km nest. These 4-km 3DVAR-, EnSRF-, and hybrid-initialized forecasts were compared to benchmark WRF forecasts initialized by interpolating 0000 UTC Global Forecast System (GFS) analyses onto the computational domain.
While important differences regarding mean state characteristics of the 20-km DA systems were noted, verification efforts focused on the 4-km precipitation forecasts. The 3DVAR-, hybrid-, and EnSRF-initialized 4-km precipitation forecasts performed similarly regarding general precipitation characteristics, such as timing of the diurnal cycle, and all three forecast sets had high precipitation biases at heavier rainfall rates. However, meaningful differences emerged regarding precipitation placement as quantified by the fractions skill score. For most forecast hours, the hybrid-initialized 4-km precipitation forecasts were better than the EnSRF-, 3DVAR-, and GFS-initialized forecasts, and the improvement was often statistically significant at the 95th percentile. These results demonstrate the potential of limited-area continuously cycling hybrid DA configurations and suggest additional hybrid development is warranted.
Abstract
The impact of assimilating radiance observations from the Advanced Microwave Sounding Unit-A (AMSU-A) on forecasts of several tropical cyclones (TCs) was studied using the Weather Research and Forecasting Model (WRF) and a limited-area ensemble Kalman filter (EnKF). Analysis/forecast cycling experiments with and without AMSU-A radiance assimilation were performed over the Atlantic Ocean for the period 11 August–13 September 2008, when five named storms formed. For convenience, the radiance forward operators and bias-correction coefficients, along with the majority of quality-control decisions, were computed by a separate, preexisting variational assimilation system. The bias-correction coefficients were obtained from 3-month offline statistics and fixed during the EnKF analysis cycles. The vertical location of each radiance observation, which is required for covariance localization in the EnKF, was taken to be the level at which the AMSU-A channels’ weighting functions peaked.
Deterministic 72-h WRF forecasts initialized from the ensemble-mean analyses were evaluated with a focus on TC prediction. Assimilating AMSU-A radiances produced better depictions of the environmental fields when compared to reanalyses and dropwindsonde observations. Radiance assimilation also resulted in substantial improvement of TC track and intensity forecasts with track-error reduction up to 16% for forecasts beyond 36 h. Additionally, assimilating both radiances and satellite winds gave markedly more benefit for TC track forecasts than solely assimilating radiances.
Abstract
The impact of assimilating radiance observations from the Advanced Microwave Sounding Unit-A (AMSU-A) on forecasts of several tropical cyclones (TCs) was studied using the Weather Research and Forecasting Model (WRF) and a limited-area ensemble Kalman filter (EnKF). Analysis/forecast cycling experiments with and without AMSU-A radiance assimilation were performed over the Atlantic Ocean for the period 11 August–13 September 2008, when five named storms formed. For convenience, the radiance forward operators and bias-correction coefficients, along with the majority of quality-control decisions, were computed by a separate, preexisting variational assimilation system. The bias-correction coefficients were obtained from 3-month offline statistics and fixed during the EnKF analysis cycles. The vertical location of each radiance observation, which is required for covariance localization in the EnKF, was taken to be the level at which the AMSU-A channels’ weighting functions peaked.
Deterministic 72-h WRF forecasts initialized from the ensemble-mean analyses were evaluated with a focus on TC prediction. Assimilating AMSU-A radiances produced better depictions of the environmental fields when compared to reanalyses and dropwindsonde observations. Radiance assimilation also resulted in substantial improvement of TC track and intensity forecasts with track-error reduction up to 16% for forecasts beyond 36 h. Additionally, assimilating both radiances and satellite winds gave markedly more benefit for TC track forecasts than solely assimilating radiances.
Abstract
Dual-resolution (DR) hybrid variational-ensemble analysis capability was implemented within the community Weather Research and Forecasting (WRF) Model data assimilation (DA) system, which is designed for limited-area applications. The DR hybrid system combines a high-resolution (HR) background, flow-dependent background error covariances (BECs) derived from a low-resolution ensemble, and observations to produce a deterministic HR analysis. As DR systems do not require HR ensembles, they are computationally cheaper than single-resolution (SR) hybrid configurations, where the background and ensemble have equal resolutions.
Single-observation tests were performed to document some characteristics of limited-area DR hybrid analyses. Additionally, the DR hybrid system was evaluated within a continuously cycling framework, where new DR hybrid analyses were produced every 6 h over ~3.5 weeks. In the DR configuration presented here, the deterministic backgrounds and analyses had 15-km horizontal grid spacing, but the 32-member WRF Model–based ensembles providing flow-dependent BECs for the hybrid had 45-km horizontal grid spacing. The DR hybrid analyses initialized 72-h WRF Model forecasts that were compared to forecasts initialized by an SR hybrid system where both the ensemble and background had 15-km horizontal grid spacing. The SR and DR hybrid systems were coupled to an ensemble adjustment Kalman filter that updated ensembles each DA cycle.
On average, forecasts initialized from 15-km DR and SR hybrid analyses were not statistically significantly different, although tropical cyclone track forecast errors favored the SR-initialized forecasts. Although additional studies over longer time periods and at finer grid spacing are needed to further understand sensitivity to ensemble perturbation resolution, these results suggest users should carefully consider whether SR hybrid systems are worth the extra cost.
Abstract
Dual-resolution (DR) hybrid variational-ensemble analysis capability was implemented within the community Weather Research and Forecasting (WRF) Model data assimilation (DA) system, which is designed for limited-area applications. The DR hybrid system combines a high-resolution (HR) background, flow-dependent background error covariances (BECs) derived from a low-resolution ensemble, and observations to produce a deterministic HR analysis. As DR systems do not require HR ensembles, they are computationally cheaper than single-resolution (SR) hybrid configurations, where the background and ensemble have equal resolutions.
Single-observation tests were performed to document some characteristics of limited-area DR hybrid analyses. Additionally, the DR hybrid system was evaluated within a continuously cycling framework, where new DR hybrid analyses were produced every 6 h over ~3.5 weeks. In the DR configuration presented here, the deterministic backgrounds and analyses had 15-km horizontal grid spacing, but the 32-member WRF Model–based ensembles providing flow-dependent BECs for the hybrid had 45-km horizontal grid spacing. The DR hybrid analyses initialized 72-h WRF Model forecasts that were compared to forecasts initialized by an SR hybrid system where both the ensemble and background had 15-km horizontal grid spacing. The SR and DR hybrid systems were coupled to an ensemble adjustment Kalman filter that updated ensembles each DA cycle.
On average, forecasts initialized from 15-km DR and SR hybrid analyses were not statistically significantly different, although tropical cyclone track forecast errors favored the SR-initialized forecasts. Although additional studies over longer time periods and at finer grid spacing are needed to further understand sensitivity to ensemble perturbation resolution, these results suggest users should carefully consider whether SR hybrid systems are worth the extra cost.
Abstract
A probability matching (PM) product using the ensemble maximum (EnMax) as the basis for spatial reassignment was developed. This PM product was called the PM max and its localized version was called the local PM (LPM) max. Both products were generated from a 10-member ensemble with 3-km horizontal grid spacing and evaluated over 364 36-h forecasts in terms of the fractions skill score. Performances of the PM max and LPM max were compared to those of the traditional PM mean and LPM mean, which both used the ensemble mean (EnMean) as the basis for spatial reassignment. Compared to observations, the PM max typically outperformed the PM mean for precipitation rates ≥5 mm h−1; this improvement was related to the EnMax, which had better spatial placement than the EnMean for heavy precipitation. However, the PM mean produced better forecasts than the PM max for lighter precipitation. It appears that the global reassignment used to produce the PM max was responsible for its poorer performance relative to the PM mean at light precipitation rates, as the LPM max was more skillful than the LPM mean at all thresholds. These results suggest promise for PM products based on the EnMax, especially for rare events and ensembles with insufficient spread.
Abstract
A probability matching (PM) product using the ensemble maximum (EnMax) as the basis for spatial reassignment was developed. This PM product was called the PM max and its localized version was called the local PM (LPM) max. Both products were generated from a 10-member ensemble with 3-km horizontal grid spacing and evaluated over 364 36-h forecasts in terms of the fractions skill score. Performances of the PM max and LPM max were compared to those of the traditional PM mean and LPM mean, which both used the ensemble mean (EnMean) as the basis for spatial reassignment. Compared to observations, the PM max typically outperformed the PM mean for precipitation rates ≥5 mm h−1; this improvement was related to the EnMax, which had better spatial placement than the EnMean for heavy precipitation. However, the PM mean produced better forecasts than the PM max for lighter precipitation. It appears that the global reassignment used to produce the PM max was responsible for its poorer performance relative to the PM mean at light precipitation rates, as the LPM max was more skillful than the LPM mean at all thresholds. These results suggest promise for PM products based on the EnMax, especially for rare events and ensembles with insufficient spread.
Abstract
Two parallel experiments were designed to evaluate whether assimilating microwave radiances with a cyclic, limited-area ensemble adjustment Kalman filter (EAKF) could improve track, intensity, and precipitation forecasts of Typhoon Morakot (2009). The experiments were configured identically, except that one assimilated microwave radiances and the other did not. Both experiments produced EAKF analyses every 6 h between 1800 UTC 3 August and 1200 UTC 9 August 2009, and the mean analyses initialized 72-h Weather Research and Forecasting model forecasts. Examination of individual forecasts and average error statistics revealed that assimilating microwave radiances ultimately resulted in better intensity forecasts compared to when radiances were withheld. However, radiance assimilation did not substantially impact track forecasts, and the impact on precipitation forecasts was mixed. Overall, net positive results suggest that assimilating microwave radiances with a limited-area EAKF system is beneficial for tropical cyclone prediction, but additional studies are needed.
Abstract
Two parallel experiments were designed to evaluate whether assimilating microwave radiances with a cyclic, limited-area ensemble adjustment Kalman filter (EAKF) could improve track, intensity, and precipitation forecasts of Typhoon Morakot (2009). The experiments were configured identically, except that one assimilated microwave radiances and the other did not. Both experiments produced EAKF analyses every 6 h between 1800 UTC 3 August and 1200 UTC 9 August 2009, and the mean analyses initialized 72-h Weather Research and Forecasting model forecasts. Examination of individual forecasts and average error statistics revealed that assimilating microwave radiances ultimately resulted in better intensity forecasts compared to when radiances were withheld. However, radiance assimilation did not substantially impact track forecasts, and the impact on precipitation forecasts was mixed. Overall, net positive results suggest that assimilating microwave radiances with a limited-area EAKF system is beneficial for tropical cyclone prediction, but additional studies are needed.
Abstract
To improve the wind speed forecasts at turbine locations and at hub height, this study develops the WRFDA system to assimilate the wind speed observations measured on the nacelle of turbines (hereafter referred as turbine wind speed observations) with both 3DVAR and 4DVAR algorithms. Results exhibit that the developed data assimilation (DA) system helps in greatly improving the analysis and the forecast of wind turbine speed. Among three experiments with no cycling DA, with 2-h cycling DA, and with 4-h cycling DA, the last experiment generates the best analysis, improving the averaged forecasts (from T + 9 to T + 24) of wind speed over all wind farms by 32.5% in the bias and 6.3% in the RMSE. After processing the turbine wind speed observations into superobs, even bigger improvements are revealed when validating against either the original turbine wind speed observations or the superobs. Taken the results validated against the superobs as an example, the bias and RMSE of the forecasts (from T + 9 to T + 24) averaged over all wind farms are reduced by 38.8% and 12.0%, respectively. Compared to the best-performed 3DVAR experiment (4-h cycling and superobs), the experiment following the same DA strategy but using 4DVAR algorithm exhibits further improvements, especially for the averaged bias in the forecasts of all wind farms, and the changing amount in the forecasts of the enhanced wind farms. Compared to the control experiment, the 4DVAR experiment reduces the bias and RMSE in the forecasts (from T + 9 to T + 24) by 54.6% (0.66 m s−1) and 12.7% (0.34 m s−1).
Abstract
To improve the wind speed forecasts at turbine locations and at hub height, this study develops the WRFDA system to assimilate the wind speed observations measured on the nacelle of turbines (hereafter referred as turbine wind speed observations) with both 3DVAR and 4DVAR algorithms. Results exhibit that the developed data assimilation (DA) system helps in greatly improving the analysis and the forecast of wind turbine speed. Among three experiments with no cycling DA, with 2-h cycling DA, and with 4-h cycling DA, the last experiment generates the best analysis, improving the averaged forecasts (from T + 9 to T + 24) of wind speed over all wind farms by 32.5% in the bias and 6.3% in the RMSE. After processing the turbine wind speed observations into superobs, even bigger improvements are revealed when validating against either the original turbine wind speed observations or the superobs. Taken the results validated against the superobs as an example, the bias and RMSE of the forecasts (from T + 9 to T + 24) averaged over all wind farms are reduced by 38.8% and 12.0%, respectively. Compared to the best-performed 3DVAR experiment (4-h cycling and superobs), the experiment following the same DA strategy but using 4DVAR algorithm exhibits further improvements, especially for the averaged bias in the forecasts of all wind farms, and the changing amount in the forecasts of the enhanced wind farms. Compared to the control experiment, the 4DVAR experiment reduces the bias and RMSE in the forecasts (from T + 9 to T + 24) by 54.6% (0.66 m s−1) and 12.7% (0.34 m s−1).
Abstract
The microphysical parameterization scheme employed in four-dimensional variational data assimilation (4D-Var) plays an important role in the assimilation of humidity and cloud-sensitive observations. In this study, a newly developed full-hydrometeor assimilation scheme, integrating warm-rain and cold-cloud processes, has been implemented in the Weather Research and Forecasting (WRF) 4D-Var system. This scheme is based on the WRF single-moment 6-class microphysics scheme (WSM6). Its primary objective is to directly assimilate radar reflectivity observations, with the goal of evaluating its effects on model initialization and subsequent forecasting performance. Four assimilation experiments were conducted to assess the performance of the full-hydrometeor assimilation scheme against the warm-rain assimilation scheme. These experiments also investigated reflectivity assimilation using both indirect and direct methods. We found that the nonlinearity of the radar operator in the two direct reflectivity assimilation experiments requires more iterations for cost function reduction than in the indirect assimilation method. The hydrometeor fields were reasonably analyzed using the full-hydrometeor assimilation scheme, particularly improving the simulation of ice-phase hydrometeors and reflectivity above the melting layer. The assimilation of radar reflectivity led to improvements in short-term (0–3 h) precipitation forecasting with the full-hydrometeor assimilation scheme. Assimilation experiments across multiple case studies reaffirmed that assimilating radar reflectivity observations with the full-hydrometeor assimilation scheme accelerated model spinup and yielded enhancements in 0–3-h total accumulate precipitation forecasts for a range of precipitation thresholds.
Abstract
The microphysical parameterization scheme employed in four-dimensional variational data assimilation (4D-Var) plays an important role in the assimilation of humidity and cloud-sensitive observations. In this study, a newly developed full-hydrometeor assimilation scheme, integrating warm-rain and cold-cloud processes, has been implemented in the Weather Research and Forecasting (WRF) 4D-Var system. This scheme is based on the WRF single-moment 6-class microphysics scheme (WSM6). Its primary objective is to directly assimilate radar reflectivity observations, with the goal of evaluating its effects on model initialization and subsequent forecasting performance. Four assimilation experiments were conducted to assess the performance of the full-hydrometeor assimilation scheme against the warm-rain assimilation scheme. These experiments also investigated reflectivity assimilation using both indirect and direct methods. We found that the nonlinearity of the radar operator in the two direct reflectivity assimilation experiments requires more iterations for cost function reduction than in the indirect assimilation method. The hydrometeor fields were reasonably analyzed using the full-hydrometeor assimilation scheme, particularly improving the simulation of ice-phase hydrometeors and reflectivity above the melting layer. The assimilation of radar reflectivity led to improvements in short-term (0–3 h) precipitation forecasting with the full-hydrometeor assimilation scheme. Assimilation experiments across multiple case studies reaffirmed that assimilating radar reflectivity observations with the full-hydrometeor assimilation scheme accelerated model spinup and yielded enhancements in 0–3-h total accumulate precipitation forecasts for a range of precipitation thresholds.
Abstract
A variational bias correction (VarBC) scheme is developed and tested using regional Weather Research and Forecasting Model Data Assimilation (WRFDA) to correct systematic errors in aircraft-based measurements of temperature produced by the Tropospheric Airborne Meteorological Data Reporting (TAMDAR) system. Various bias models were investigated, using one or all of aircraft height tendency, Mach number, temperature tendency, and the observed temperature as predictors. These variables were expected to account for the representation of some well-known error sources contributing to uncertainties in TAMDAR temperature measurements. The parameters corresponding to these predictors were evolved in the model for a two-week period to generate initial estimates according to each unique aircraft tail number. Sensitivity experiments were then conducted for another one-month period. Finally, a case study using VarBC of a cold front precipitation event is examined. The implementation of VarBC reduces biases in TAMDAR temperature innovations. Even when using a bias model containing a single predictor, such as height tendency or Mach number, the VarBC produces positive impacts on analyses and short-range forecasts of temperature with smaller standard deviations and biases than the control run. Additionally, by employing a multiple-predictor bias model, which describes the statistical relations between innovations and predictors, and uses coefficients to control the evolution of components in the bias model with respect to their reference values, VarBC further reduces the average error of analyses and short-range forecasts with respect to observations. The potential impacts of VarBC on precipitation forecasts were evaluated, and the VarBC is able to indirectly improve the prediction of precipitation location by reducing the forecast error for wind-related synoptic circulation leading to precipitation.
Abstract
A variational bias correction (VarBC) scheme is developed and tested using regional Weather Research and Forecasting Model Data Assimilation (WRFDA) to correct systematic errors in aircraft-based measurements of temperature produced by the Tropospheric Airborne Meteorological Data Reporting (TAMDAR) system. Various bias models were investigated, using one or all of aircraft height tendency, Mach number, temperature tendency, and the observed temperature as predictors. These variables were expected to account for the representation of some well-known error sources contributing to uncertainties in TAMDAR temperature measurements. The parameters corresponding to these predictors were evolved in the model for a two-week period to generate initial estimates according to each unique aircraft tail number. Sensitivity experiments were then conducted for another one-month period. Finally, a case study using VarBC of a cold front precipitation event is examined. The implementation of VarBC reduces biases in TAMDAR temperature innovations. Even when using a bias model containing a single predictor, such as height tendency or Mach number, the VarBC produces positive impacts on analyses and short-range forecasts of temperature with smaller standard deviations and biases than the control run. Additionally, by employing a multiple-predictor bias model, which describes the statistical relations between innovations and predictors, and uses coefficients to control the evolution of components in the bias model with respect to their reference values, VarBC further reduces the average error of analyses and short-range forecasts with respect to observations. The potential impacts of VarBC on precipitation forecasts were evaluated, and the VarBC is able to indirectly improve the prediction of precipitation location by reducing the forecast error for wind-related synoptic circulation leading to precipitation.
Abstract
This study examines the impact of assimilating Microwave Humidity Sounder (MHS) radiances in a limited-area ensemble Kalman filter (EnKF) data assimilation system. Two experiments spanning 11 August–13 September 2008 were run over a domain featuring the Atlantic basin using a 6-h full cycling analysis and forecast system. Deterministic 72-h forecasts were initialized at 0000 and 1200 UTC for a comparison of forecast impact. The two experiments were configured identically with the exception of the inclusion of the MHS radiances (AMHS) in the second to isolate the impacts of the MHS radiance data. The results were verified against several sources, and statistical significance tests indicate the most notable differences are in the midlevel moisture fields. Both configurations were characterized by high moisture biases when compared to the European Centre for Medium-Range Weather Forecasts interim reanalysis (ERA-Interim, also known as ERA-I) specific humidity fields, as well as precipitable water vapor from an observationally based product. However, the AMHS experiment has midlevel moisture fields closer to the ERA-I and observation datasets. When reducing the verification domain to focus on the subtropical and easterly wave regions of the North Atlantic Ocean, larger improvements in midlevel moisture at nearly all lead times is seen in the AMHS simulation. Finally, when considering tropical cyclone forecasts, the AMHS configuration shows improvement in intensity forecasts at several lead times as well as improvements at early to intermediate lead times for minimum sea level pressure forecasts.
Abstract
This study examines the impact of assimilating Microwave Humidity Sounder (MHS) radiances in a limited-area ensemble Kalman filter (EnKF) data assimilation system. Two experiments spanning 11 August–13 September 2008 were run over a domain featuring the Atlantic basin using a 6-h full cycling analysis and forecast system. Deterministic 72-h forecasts were initialized at 0000 and 1200 UTC for a comparison of forecast impact. The two experiments were configured identically with the exception of the inclusion of the MHS radiances (AMHS) in the second to isolate the impacts of the MHS radiance data. The results were verified against several sources, and statistical significance tests indicate the most notable differences are in the midlevel moisture fields. Both configurations were characterized by high moisture biases when compared to the European Centre for Medium-Range Weather Forecasts interim reanalysis (ERA-Interim, also known as ERA-I) specific humidity fields, as well as precipitable water vapor from an observationally based product. However, the AMHS experiment has midlevel moisture fields closer to the ERA-I and observation datasets. When reducing the verification domain to focus on the subtropical and easterly wave regions of the North Atlantic Ocean, larger improvements in midlevel moisture at nearly all lead times is seen in the AMHS simulation. Finally, when considering tropical cyclone forecasts, the AMHS configuration shows improvement in intensity forecasts at several lead times as well as improvements at early to intermediate lead times for minimum sea level pressure forecasts.
Abstract
The Weather Research and Forecasting Model (WRF) “hybrid” variational-ensemble data assimilation (DA) algorithm was used to initialize WRF model forecasts of three tropical cyclones (TCs). The hybrid-initialized forecasts were compared to forecasts initialized by WRF's three-dimensional variational (3DVAR) DA system. An ensemble adjustment Kalman filter (EAKF) updated a 32-member WRF-based ensemble system that provided flow-dependent background error covariances for the hybrid. The 3DVAR, hybrid, and EAKF configurations cycled continuously for ~3.5 weeks and produced new analyses every 6 h that initialized 72-h WRF forecasts with 45-km horizontal grid spacing. Additionally, the impact of employing a TC relocation technique and using multiple outer loops (OLs) in the 3DVAR and hybrid minimizations were explored.
Model output was compared to conventional, dropwindsonde, and TC “best track” observations. On average, the hybrid produced superior forecasts compared to 3DVAR when only one OL was used during minimization. However, when three OLs were employed, 3DVAR forecasts were dramatically improved but the mean hybrid performance changed little. Additionally, incorporation of TC relocation within the cycling systems further improved the mean 3DVAR-initialized forecasts but the average hybrid-initialized forecasts were nearly unchanged.
Abstract
The Weather Research and Forecasting Model (WRF) “hybrid” variational-ensemble data assimilation (DA) algorithm was used to initialize WRF model forecasts of three tropical cyclones (TCs). The hybrid-initialized forecasts were compared to forecasts initialized by WRF's three-dimensional variational (3DVAR) DA system. An ensemble adjustment Kalman filter (EAKF) updated a 32-member WRF-based ensemble system that provided flow-dependent background error covariances for the hybrid. The 3DVAR, hybrid, and EAKF configurations cycled continuously for ~3.5 weeks and produced new analyses every 6 h that initialized 72-h WRF forecasts with 45-km horizontal grid spacing. Additionally, the impact of employing a TC relocation technique and using multiple outer loops (OLs) in the 3DVAR and hybrid minimizations were explored.
Model output was compared to conventional, dropwindsonde, and TC “best track” observations. On average, the hybrid produced superior forecasts compared to 3DVAR when only one OL was used during minimization. However, when three OLs were employed, 3DVAR forecasts were dramatically improved but the mean hybrid performance changed little. Additionally, incorporation of TC relocation within the cycling systems further improved the mean 3DVAR-initialized forecasts but the average hybrid-initialized forecasts were nearly unchanged.
