Next Article in Journal
Source Analysis and Human Health Risk Assessment Based on Entropy Weight Method Modification of PM2.5 Heavy Metal in an Industrial Area in the Northeast of China
Next Article in Special Issue
A Comparison of HWRF Six-Hourly 4DEnVar and Hourly 3DEnVar Assimilation of Inner Core Tail Dopper Radar Observations for the Prediction of Hurricane Edouard (2014)
Previous Article in Journal
Effects of Excessive Equatorial Cold Tongue Bias on the Projections of Tropical Pacific Climate Change. Part II: The Extreme El Niño Frequency in CMIP5 Multi-Model Ensemble
Previous Article in Special Issue
Analog Ensemble Methods for Improving Satellite-Based Intensity Estimates of Tropical Cyclones
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Applications of Radar Data Assimilation with Hydrometeor Control Variables within the WRFDA on the Prediction of Landfalling Hurricane IKE (2008)

1
Key Laboratory of Meteorological Disaster, Ministry of Education (KLME)/Joint International Research Laboratory of Climate and Environment Change (ILCEC)/Collaborative Innovation Center on Forecast and Evaluation of Meteorological Disasters (CIC-FEMD), Nanjing University of Information Science & Technology, Nanjing 210044, China
2
Heavy Rain and Drought-Flood Disasters in Plateau and Basin Key Laboratory of Sichuan Province, Chengdu 610225, China
3
The Institute of Atmospheric Environment, China Meteorological Administration, Shenyang 110000, China
4
Shanghai Typhoon Institute, China Meteorological Administration, Shanghai 200030, China
5
Chongqing Meteorological Observatory, Chongqing 401147, China
6
Key Laboratory of Agrometeorological Support and Applied Technique of Henan Province, CMA, Zhengzhou 450003, China
*
Author to whom correspondence should be addressed.
Submission received: 27 April 2021 / Revised: 22 June 2021 / Accepted: 25 June 2021 / Published: 30 June 2021
(This article belongs to the Special Issue Tropical Cyclones: Observation and Prediction)

Abstract

:
The impact of assimilating radar radial velocity and reflectivity on the analyses and forecast of Hurricane IKE is investigated within the framework of the WRF (Weather Research and Forecasting) model and its three-dimensional variational (3DVar) data assimilation system, including the hydrometeor control variables. Hurricane IKE in the year 2008 was chosen as the study case. It was found that assimilating radar data is able to effectively improve the small-scale information of the hurricane vortex area in the model background. Radar data assimilation experiments yield significant cyclonic wind increments in the inner-core area of the hurricane, enhancing the intensity of the hurricane in the model background. On the other hand, by extending the traditional control variables to include the hydrometeor control variables, the assimilation of radar reflectivity can effectively adjust the water vapor and hydrometeors of the background, further improving the track and intensity forecast of the hurricane. The precipitation forecast skill is also enhanced to some extent with the radar data assimilation, especially with the extended hydrometeor control variables.

1. Introduction

Tropical cyclones (TC) have become one of the most important natural disasters in the world because of their great destructive power. The accurate forecast of TC track and intensity can greatly reduce unnecessary economic losses. In the past 20 years, some achievements have been made in the prediction of TC track, but there are still many works to be further done in terms of TC intensity prediction (Dong and Xue [1] (2013)). For a large number of TCs, most of their lifespan is in the open ocean, and their movement and development are often influenced by the large-scale steering flow and the underlying surface. Therefore, the improvement of TC prediction mainly depends on numerical weather prediction (NWP). With the development of the NWP model, the accuracy of the initial condition (IC) has become a vital factor in determining the forecast skills of the TC systems.
Generally, conventional observation data are relatively scarcely distributed in the open ocean with relatively low temporal resolution. Thus, it is highly essential to introduce unconventional observation data (such as radar observations, satellite observations, etc.) to improve the ICs of the NWP model. As one of the most outstanding observing systems, Doppler weather radar is able to provide multi-dimensional thermal and dynamic observations [2,3]. Radar radial velocity and reflectivity can provide information on wind field and hydrometeors, respectively. How to effectively use radar observations to initialize the NWP model becomes a notable issue.
In the past decades, experts and scholars have done a lot of research on radar data assimilation. Xiao et al. [4] (2007) developed a direct assimilation scheme of radar reflectivity using total hydrometeors as the control variable, which used the warm rain scheme to refine total hydrometeors into water vapor, cloud rain, and rain phase states. Gao et al. [5] (2012) used the background temperature field of the NWP model to divide the hydrometeors, introduced the ice phase hydrometeors, and improved the analysis and prediction of storms; Wang et al. [6] (2013) used an indirect radar reflectivity assimilation scheme to improve the prediction of local precipitation by introducing a new observation operator based on the z-qr (reflectivity-rain) equation to avoid the linearization error of the observation operator. Several studies have focused on TC initialization based on radar data assimilation, among which most researches mainly focused on radar radial velocity observations [7,8,9,10,11]. However, the assimilation of radar reflectivity data, especially those with multiple hydrometeor control variables, obviously need more detailed investigation.
The control variable is the key component of the variational data assimilation technique in terms of optimizing the meteorological state by minimizing the cost function. Generally, the state control variables of the NWP model, including [12]: (1) momentum related control variables (e.g., stream function ψ and velocity potential χ); (2) unbalanced temperature ( T u ); (3) unbalanced surface pressure ( P s , u ); and (4) humidity related control variables (e.g., pseudo relative humidity ( R H s )), are adjusted or expanded to better optimize the analysis. Sun et al. [13] (2016) used horizontal wind components as control variables for improving the convective scale analysis. Multivariable correlation methods are proposed by Michel et al. (2011) and Descombes et al. (2015), which consider hydrometeors for the control variables to improve the analysis [14,15]. Meng et al. [16] (2018) adopted the EnVar method to include hydrometeors variable to effectively improve the simulation of the convective available potential energy and humidity of a convective system, leading to improvement in the prediction for the spatial distribution and intensity of the precipitation. Zhao et al. [17] (2008) updated the hydrometers in the NWP model directly by applying general analytic relationship when assimilating the radar reflectivity in two storm cases with promising results. As far as the authors know, there is very few published work with radar reflectivity within the framework of WRFDA with UV control variables (Li et al., (2016) [18], Sun et al., (2016) [13]). The control variables related to hydrometeors have become one of the key technologies to effectively assimilate hydrometeor-related observations, such as radar reflectivity. However, up to now, there are not many researches on radar data assimilation technology based on the control variables of hydrometeors, and there are few researches on radar reflectivity assimilation for landing TCs in America. Therefore, this study will discuss the impact of assimilation of radar reflectivity on the analysis and prediction of Hurricane IKE (2008) by using control variables from hydrometeors based on the WRF model and its assimilation system.

2. Methodology

2.1. WRFDA 3DVar Data Assimilation System

Data assimilation aims to obtain optimized initial conditions of the NWP model by analyzing multi-source observations [19,20]. The typical three-dimensional variational DA technique in WRFDA specific formula is as follows:
J = ( x x b ) T B 1 ( x x b ) + ( y o H ( x ) ) T R 1 ( y o H ( x ) )
where vectors x and xb represent the atmospheric state and background state [21], y0 represents the observation vector. H is the observation operator. B and R are background error covariance matrix, and observation error covariance matrix, respectively.
In this study, the NMC (National Meteorological Center) method [22] is used to estimate the background error covariance. The differences in the forecasts valid at the same time but with different forecast lengths (24 h forecast and 12 h forecast) are taken as the approximation of the forecast error. The forecasts are launched twice a day from 1 September to 20 September 2008 to provide a large number of prediction samples for the estimation of background error covariance. According to the differences between samples, estimate the covariance matrix of background error [18], as follows:
B ( x 24 x 12 ) ( x 24 x 12 ) T ¯
Among them, x 24 and x 12 are 24 h and 12 h forecasts, respectively. It should be pointed out that both prediction fields have errors and are expressed as follows:
x 24 = x truth + ε 24 + b 24
x 12 = x truth + ε 12 + b 12
Here, x truth represents the real value. b 24 and b 12 are the biases in each prediction. ε 24 and ε 12 are random errors. Suppose that the prediction of 24 h and 12 h is unbiased, or the deviation is constant in time, because b 24 = b 12 , the prediction difference is expressed as:
x diff = ε 24 ε 12
B matrix can be expressed by formula (6):
B ( x diff ) ( x diff ) T ¯ = ( ε 24 ε 12 ) ( ε 24 ε 12 ) T ¯ ,
In the data assimilation system, due to the large dimension of the background error covariance matrix, it is difficult to conduct the minimum. Generally, the B matrix is decomposed into B = U U T , and the following preprocessing process can be formed through decomposition U v = x x b . δ x = x x b stands for the analysis increment and v represents the control variable. One approximation is applied as:
y o H ( x ) y o H ( x b ) H ( x x b ) ,
where H is the linearized observation operator of H. The final format of the cost function of 3DVAR is:
J = 1 2 v T v + 1 2 ( d H U v ) T R 1 ( d H U v ) ,
where d = y o H ( x b ) is the innovation.
There are two types of momentum control variables in WRFDA-3DVar. The default option in WRFDA-3DVar is the CV5 scheme. The control variables include the stream function ψ, the unbalanced velocity potential xu, the unbalanced temperature Tu, and the unbalanced surface pressure P su and pseudo relative humidity RHs. Another option is the CV7 scheme in which velocity component (U, V) is used as a dynamic control variable. Other control variables include total temperature (T), pseudo relative humidity (RHS), and total surface pressure (PS). Previous studies have discussed the influence of the two different momentum control variables in assimilating radar radial velocity [9,12,13,23].
In order to make full use of the radar data, control variables are extended to include cloud and precipitation information. The control variables used in this study can be designed as follows:
CV5: stream function (ψ), unbalanced velocity potential (Chi_u), unbalanced surface pressure (Ps_u), unbalanced temperature (T_u), pseudo relative humidity (RHs).
CV7: zonal wind (U), meridional wind (V), surface pressure (Ps), temperature (T), pseudo relative humidity (RHs).
CV8: zonal wind (U), meridional wind (V), surface pressure (Ps), temperature (T), pseudo relative humidity (RHs), cloud water (qc), rain water (qr), ice water (qi), snow water (qs), graupel (qg).

2.2. Radar Observation Operator

According to Tong and Xue (2005) [24], the observation operator of radar radial wind is calculated based on the wind field as:
V r = u x x i r i + v y y i r i + ( w v T ) z z i r i
where (x, y, z) denotes the radar location elevation, (xi, yi, zi) is the location for the radar observation, and (u, v, w) is the three-dimensional wind field. ri is the distance between the observations and the radar location. vT (m/s) is the terminal speeds. The effects of atmospheric refraction and earth curvature on radar-beam height and slope angle is neglected [25].
For reflectivity factor (RF) data, it should be noted that the z-qr (reflectivity-rain) relationship [26] is used as the observation operator in WRFDA for direct assimilation. Usually, the humidity control variable is expanded to the total hydrometeor variable with the warm rain scheme without any ice phase process considered. Thus, there are obvious defects in terms of describing the microphysical process of the deep convection system. In addition, WRFDA adopts the incremental method in the minimization based on the linearization of the forward observation operator. The linearization of the RF operator is rather sensitive to the humidity field in the model background field. When the model background field is dry, a large analysis error is introduced. Therefore, indirect assimilation of RF data is popularly used based on the assumption that the water vapor in the cloud is saturated (Wang et al., 2013) [27]. To be specific, when the reflectivity factor is greater than 30 dBZ, the relative humidity in the cloud is considered as 100%. In this way, the pseudo observation q v o can be obtained based on the pressure and temperature at that time. For this pseudo observation q v o , its nonlinear observation operator is:
q v = r h × q v s ,
where qv stands for the water vapor mixing ratio (unit, kg/kg), rh is the relative humidity, and qvs is the saturated water vapor mixing ratio, respectively. The forward observation operator of reflectivity factor is based on the interaction of rain mixing ratio, snow mixing ratio, and graupel mixing ratio. A more advanced water vapor assimilation technique to alleviate the potential overprediction of precipitation is beyond the scope of the current study (refer to Xu et al. (2010) [28]).

3. Case Introduction and Experimental Setup

3.1. Hurricane IKE

On 1 September 2008, Hurricane IKE developed from a tropical disturbance in the west of Cape Verde. Later, its intensity gradually increased to Category 4 hurricane in the central Atlantic on September 4 by moving westward. On September 9, IKE made landfall in eastern Cuba and was weakened prior to entering the Gulf of Mexico. On 13 September, it increased by the time of the final landfall in Galveston, Texas, as an extratropical storm, causing considerable damage inland.
The environment had a great influence on the evolution of TC. Figure 1 is the circulation from FNL analysis data. As can be seen from Figure 1a at 850 hPa, the TC center was affected by the southerly wind flow, which is able to supply sufficient moisture. It is noted that the water vapor flux near the TC center reached 40 g/(cm·hPa·s). At 500 hPa in Figure 1b, it can be seen that IKE was embedded in the southwest periphery of the subtropical. The easterly jet on the south side of the subtropical made the high-level TC outflow stronger, leading the TC to move westward (Wang (2011) [29]).

3.2. Model and Experimental Design

The WRF model is a commonly used NWP model in weather research and forecast, which mainly serves the meteorological application of space scale from tens of meters to thousands of kilometers. In this study, the ARW (the advanced research WRF) version for research users is used. In this study, the 4.2 version of the WRF model is used to simulate Hurricane IKE, and the center of simulation domain is set at 30° N, −95° W. The resolution of the domain is 4 km (Figure 2) following Dong and Xue [1] (2013), and Li et al., [30] (2015). The number of grid points of the domain are 501 × 501, respectively, and the number of vertical layers is 53 with 10 hPa as the model top. To provide initial and boundary conditions for all experiments, 0.5° GFS (Global Forecast System) data with a 6-hour interval are used. The parameterization of the numerical experiment includes: (1) RRTM (rapid radiative transfer model) long wave radiation scheme, (2) WDM6 (WRF Double-Moment 6-class) Microphysics scheme, (3) Dudhia short wave radiation scheme, (4) Noah land surface model land surface scheme, (5) MM5 similarity near ground surface scheme, and (6) YSU (Yonsei University) boundary layer scheme.
The radar data used in this study are from the radars of Weather Surveillance Radar 88 Doppler (WSR-88D) Houston, Texas (KHGX), and Lake Charles, Louisiana (KLCH). Before the data assimilation, the 88D2ARPS module is used to perform the quality control of the radar observations. Finally, the processed data are interpolated to the grid of WRF model for assimilation. In this study, the observation error of radial wind is set as 2.5 m/s, and the observation error of reflectivity is 5 dBZ. The details of the radar observations can be found in Shen et al. (2017) [7] for reference.
In this study, three experiments are performed. The first experiment is the control experiment named as Conv, which only assimilates conventional observations from Global Telecommunications System by using the UV control variable scheme; the second experiment conducted the radar data assimilation with the same control variable scheme as in the first experiment; the third experiment further includes the control variables of hydrometeors in the radar data assimilation named as DA_Radar_CVH.
In the control experiment, a 6 h forecast is first carried out from 1800 UTC on 12 September 2008, before GTS conventional observations are assimilated from 0000 UTC on 13 September 2008 to 0600 UTC on 13 September 2008 every 30 min. At the end of DA, an 18 h deterministic forecast is launched at 0600 UTC on 13 September 2008. Figure 3 shows the distribution of conventional observation data of GTS in two different analysis times.
Similarly, two sets of radar data assimilation experiments are conducted with the same background used in the DA_Conv experiment at 0000 UTC 13 September. Then, the radar radial velocity and reflectivity were assimilated from 0000 UTC 13 September to 0600 UTC 13 September every 30 min. With the final analysis, a deterministic 18 h forecast is carried out from 0600 UTC 13 September to 0000 UTC 14 September (Figure 4).

4. Results

4.1. Wind Analysis

Figure 5 shows the wind increment at 700 hPa of the radar DA experiment at the first assimilation time (i.e., 0000 UTC 13 September). The corresponding radar radial velocity from the observation can be referred to in Shen et al. [8] (2018). It can be seen that the inner-core area of TC is totally covered by radar observations, and the wind increment is generated in both the TC core area and the outer environment field. It should be noted that the pattern of wind increment is cyclonic circulation, and the wind increment at the center of TC is most obvious. The results suggested that the intensity of TC simulated by DA_Conv experiment is still weak after assimilating the GTS observation data, which is commonly found in the simulation of the IKE case without any radar data assimilation [1,8]. On the other hand, the direct assimilation of radar radial wind observations is able to enhance the intensity of TC by adjusting the dynamic field of Hurricane IKE.
In order to objectively check whether the radar data are properly assimilated, the scatter plots of the radar radial velocity in the background and the analysis versus the observed radar radial velocity are shown (Figure 6). It can be seen that the scatters of the radar radial velocity of the background and the observations are obviously divergent from the diagonal. It means that the wind speed simulated by the model background has significant bias from the observation. It also can be seen that the assimilation of radar radial wind data is able to adjust the background field significantly in Figure 6b to fit the observation with reduced root mean square error (RMSE) and standard deviation.
The influence of radar data assimilation is also evaluated by calculating the RMSE of the simulated radial velocity against the observed radial velocity. The RMSE variations from each background and analysis for the 13 DA cycles are illustrated in Figure 7. It can be found that after radar data assimilation, the error of the analysis field at each assimilation time is significantly reduced. The RMSE error of the first analysis cycle (0000 UTC on 13 September 2008) is reduced from 9.4 m/s to 2.8 m/s, while the RMSE of the second analysis cycle is 2.6 m/s. Generally, the RMSE of each analysis in subsequent cycles is lower than 2.4 m/s. The magnitude of the RMSE is consistent with the radar observation error of 2.5 m/s that is applied in our current study.

4.2. Increment for the Hydrometeors

Figure 8 shows the increment of the hydrometeors water path of DA_Radar experiment (Figure 8a,c,e,g) and DA_Radar_CVH experiment (Figure 8b,d,f,h) at 0300 UTC on 13 September 2008. In the Radar experiment, increments of these hydrometeors are created in the procedure of radar data assimilation, since the observation operators of the indirect assimilation include the variables of rain, snow water, and hail water. Compared with the Radar experiment, the increment of hydrometeors in the DA_Radar_CVH experiment is more obvious, especially near the core of the TC, due to the use of rain water, snow water, and hail water variables as the inputs in the observation operator as well as the inclusion of all hydrometeors in the control variables.
Figure 9 shows the RMS profile of analysis increment of the model variables and hydrometeors variables at 0000 UTC on 13 September, 0300 UTC on 13 September, 0600 UTC on 13 September. It can be seen that at 0000 UTC on 13 September, the increment of hydrometeors was less than 0.06 g/kg. With the increase in assimilation times, the increment of hydrometeors gradually increased for most hydrometeors. At 0600 UTC on 13 September, the maximum of the hydrometeor increment reached 0.08 g/kg, indicating that the radar data assimilation has made significant modifications in the hydrometeors. In addition, it is worth noting that the maximum increment of qr existed mainly lower than the 19th level, while the maximum increment of qs was located higher than the 23th level, and the maximum increment of qg was at about the 20th level.

4.3. Precipitation Forecasts

The accurate forecast of TC precipitation is of great significance to the flood warning of inland coastal areas. Figure 10 shows the 3 h accumulated precipitation of all the experiments from 0600 UTC, 13 September to 0900 UTC, 13 September. The observed large precipitation center is located southeast of the State of Texas. All the experiments are able to simulate the distribution of the rainfall to some extent, while the precipitation center in the northeast of TC core is largely underestimated. It can be seen that the predicted precipitation intensity in the DA_Conv experiment in the southeast of the TC core is weaker than the observation, while the intensity in two radar data assimilation experiments is stronger than the observation.
Furtherly, the TS (Threat Score) of 3 h accumulated precipitation in each experiment is evaluated for different thresholds (Figure 11). Two radar data assimilation experiments obviously yield higher TS scores than DA_Conv experiment does for all thresholds. For thresholds less than 25mm, the score of DA_Radar_CVH experiment is much higher than those of the Radar experiment. However, the score of the DA_Radar_CVH experiment is similar to that of the Radar experiment for a 25 mm threshold.

4.4. Track Forecasts

In this study, an 18 h deterministic forecast was carried out at the end of the analysis. As shown in Figure 12, the DA_Conv experiment had a significant westward bias for the track, whereas the tracks from the two radar DA experiments match better with the best track. Compared to the Radar experiment, the DA_Radar_CVH experiment yielded track fitting better with the best track, especially after 6 h. The track error of the DA_Conv experiment at the initial time was less than 14 km, while the errors of the two radar data assimilation experiments were reduced to about 4.5 km at the initial time. After 12 h of prediction, it can be found that the tracks from the three experiments all shifted to the northwest to some extent, which is consistent with the best track. At 0000 UTC on 14 September, the track errors of the DA_Conv, DA_Radar, and DA_Radar_CVH experiments were 92.8 km, 70.5 km, and 50.7 km, respectively, which are their maximum errors. Overall, the mean value of track error of DA_Radar_CVH was the smallest of all the three experiments.

5. Summary

In this study, the impact of assimilating radar radial velocity and reflectivity on the analysis and forecast of the Hurricane IKE is evaluated with hydrometeors considered as control variables. The WRF model and its assimilation system WRFDA are applied for the simulation of the strong Hurricane IKE in 2008. The main results and conclusions are as follows:
The TC intensity simulated by the model is weak without radar data assimilation. After assimilating the radar radial velocity data, a more obvious cyclonic circulation is generated in the TC area of the model analysis. Generally, the assimilation of radar radial velocity data can effectively enhance the dynamic field of the TC. The hydrometeors are added to the control variables to improve the analysis of the cloud and precipitation. Outstanding increments for all hydrometeors exist after considering the all-type hydrometeors in the control variables. Compared with the control experiment, the forecast skills of the precipitation and track from the two radar data assimilation experiments are enhanced. The forecasts of the precipitation and track are further improved with the use of the hydrometeor control variables.
Although the results are encouraging, there are areas for further improvements and explorations. Further studies with nested domains with more advanced data assimilation methods (such as 4DVar and ensemble techniques) are needed for better understanding of the impacts of including additional control variables from hydrometers. Our effort in this study represents a step in that direction. It is also worth pointing out that the impact of hydrometeors on the TCs in different regions and different stages of the development of TCs need to be further investigated. For better utilization of all types of observations (such as infrared radiance) that detect hydrometeors, it is necessary to explore the different contributions of various observations in the analytical procedure.

Author Contributions

Data curation, D.Z.; Formal analysis, J.M., H.L. and D.X.; Investigation, F.S. and H.L.; Methodology, F.S.; Project administration, Y.G.; Resources, J.M.; Software, H.L.; Validation, D.X.; Visualization, A.S.; Writing—original draft, A.S. and L.S.; Writing—review & editing, F.S., J.M., H.L., D.X., A.S., D.Z., Y.G. and L.S. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Chinese National Natural Science Foundation of China (G41805070), the Chinese National Key R&D Program of China (2018YFC1506404), the Chinese National Natural Science Foundation of China (G41805016), the Chinese National Key R&D Program of China (2018YFC1506603), the research project of Heavy Rain and Drought–Flood Disasters in Plateau and Basin Key Laboratory of Sichuan Province in China (SZKT201901, SZKT201904), the research project of the Institute of Atmospheric Environment, China Meteorological Administration, Shenyang in China (2020SYIAE07, 2020SYIAE02), and the Joint Open Project of KLME & CIC-FEMD, NUIST (KLME201906).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Acknowledgments

We acknowledge the support of the High Performance Computing Center of Nanjing University of Information Science & Technology.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Dong, J.; Xue, M. Assimilation of radial velocity and reflectivity data from coastal WSR-88D radars using an ensemble Kalman filter for the analysis and forecast of landfalling hurricane Ike (2008). Q. J. R. Meteorol. Soc. 2013, 139, 467–487. [Google Scholar] [CrossRef] [Green Version]
  2. Sokol, Z. Assimilation of extrapolated radar reflectivity into a NWP model and its impact on a precipitation forecast at high resolution. Atmos. Res. 2011, 100, 201–212. [Google Scholar] [CrossRef]
  3. Courtier, P. The ECMWF implementation of three-dimensional variational assimilation (3D-Var). I: Formulation. Q. J. R. Meteorol. Soc. 1998, 124, 1783–1807. [Google Scholar] [CrossRef]
  4. Xiao, Q.; Kuo, Y.; Sun, J.; Xiao, Q.; Kuo, Y.H.; Sun, J.; Lee, W.C.; Barker, D.M.; Lim, E. An approach of radar reflectivity data assimilation and its assessment with the inland QPF of TC Rusa (2002) at landfall. J. Appl. Meteorol. Climatol. 2007, 46, 14–22. [Google Scholar] [CrossRef] [Green Version]
  5. Gao, J.; Stensrud, D. Assimilation of reflectivity data in a convective-scale, cycled 3DVAR framework with hydrometeor classification. J. Atmos. Sci. 2012, 69, 1054–1065. [Google Scholar] [CrossRef]
  6. Wang, H.; Sun, J.; Fan, S.; Huang, X.-Y. Indirect assimilation of radar reflectivity with WRF 3D-Var and its impact on prediction of four summertime convective events. J. Appl. Meteorol. Climatol. 2013, 52, 889–902. [Google Scholar] [CrossRef]
  7. Shen, F.; Xue, M.; Min, J. A comparison of limited-area 3DVAR and ETKF-En3DVAR data assimilation using radar observations at convective scale for the prediction of TC Saomai (2006). Meteorol. Appl. 2017, 24, 628–641. [Google Scholar] [CrossRef] [Green Version]
  8. Shen, F.; Xu, D.; Xue, M.; Min, J. A comparison between EDA-EnVar and ETKF-EnVar data assimilation techniques using radar observations at convective scales through a case study of Hurricane Ike (2008). Meteorol. Atmos. Phys. 2018, 130, 649–666. [Google Scholar] [CrossRef]
  9. Shen, F.; Xu, D.; Min, J. Effect of momentum control variables on assimilating radar observations for the analysis and forecast for TC Chanthu (2010). Atmos. Res. 2019, 230, 104622. [Google Scholar] [CrossRef]
  10. Shen, F.; Xu, D.; Min, J.; Chu, Z.; Xia, X. Assimilation of radar radial velocity data with the WRF hybrid 4DEnVar system for the prediction of hurricane Ike (2008). Atmos. Res. 2020, 234, 104771. [Google Scholar] [CrossRef]
  11. Xu, D.; Shen, F.; Min, J. Effect of background error tuning on assimilating radar radial velocity observations for the forecast of hurricane tracks and intensities. Meteorol. Appl. 2020, 27, e1280. [Google Scholar] [CrossRef] [Green Version]
  12. Derber, J.; Bouttier, F. A reformulation of the background error covariance in the ECMWF global data assimilation system. Tellus A 1999, 51, 195–221. [Google Scholar] [CrossRef]
  13. Sun, J.; Wang, H.; Tong, W. Comparison of the impacts of momentum control variables on high-resolution variational data assimilation and precipitation forecasting. Mon. Weather Rev. 2016, 144, 149–169. [Google Scholar] [CrossRef]
  14. Michel, Y.; Auligné, T.; Montmerle, T. Heterogeneous convective-scale background error covariances with the inclusion of hydrometeor variables. Mon. Wea. Rev. 2011, 139, 2994–3015. [Google Scholar] [CrossRef]
  15. Descombes, G.; Auligné, T.; Vandenberghe, F.; Barker, D.; Barre, J. Generalized background error covariance matrix model (GEN_BE v2.0). Geosci. Model Dev. 2015, 8, 669–696. [Google Scholar] [CrossRef] [Green Version]
  16. Meng, D.; Chen, Y.; Wang, H.; Gao, Y.; Potthast, R.; Wang, Y. The evaluation of EnVar method including hydrometeors analysis variables for assimilating cloud liquid/ice water path on prediction of rainfall events. Atmos. Res. 2019, 219, 1–12. [Google Scholar] [CrossRef]
  17. Zhao, Q.; Cook, J.; Xu, Q.; Harasti, P. Improving short-term storm predictions by assimilating both radar radial-wind and reflectivity observations. Weather Forecast. 2008, 23, 373–391. [Google Scholar] [CrossRef]
  18. Li, X.; Zeng, M.; Wang, Y.; Wang, W.; Wu, H.; Mei, H. Evaluation of two momentum control variable schemes and their impact on the variational assimilation of radar wind data: Case study of a squall line. Adv. Atmos. Sci. 2016, 33, 1143–1157. [Google Scholar] [CrossRef]
  19. Barker, D.; Huang, X.-Y.; Liu, Z.; Auligné, T.; Zhang, X.; Rugg, S.; Ajjaji, R.; Bourgeois, A.; Bray, J.; Chen, Y.; et al. The weather research and forecasting model’s community variational/ensemble data assimilation system: WRFDA. Bull. Am. Meteorol. Soc. 2012, 93, 831–843. [Google Scholar] [CrossRef] [Green Version]
  20. Ide, K.; Courtier, P.; Ghil, M.; Lorenc, A.C. Unified notation for data assimilation: Operational sequential and variational. Meteorol. Soc. 1997, 75, 181–189. [Google Scholar] [CrossRef] [Green Version]
  21. Wang, H.; Huang, X.; Sun, J.; Xu, D.; Zhang, M.; Fan, S.; Zhong, J. Inhomogeneous background error modeling for WRF-Var using the NMC method. Appl. Meteor. Clim. 2014, 53, 2287–2309. [Google Scholar] [CrossRef] [Green Version]
  22. Parrish, D.F.; Derber, J.C. The national meteorological center’s spectral statistical-interpolation analysis system. Mon. Weather Rev. 1992, 120, 1747–1763. [Google Scholar] [CrossRef]
  23. Xu, Q. On the choice of momentum control variables and covariance modeling for mesoscale data assimilation. J. Atmos. Sci. 2019, 76, 89–111. [Google Scholar] [CrossRef]
  24. Tong, M.; Xue, M. Ensemble Kalman filter assimilation of Doppler radar data with a compressible nonhydrostatic model: OSS experiments. Mon. Weather Rev. 2005, 133, 1789–1807. [Google Scholar] [CrossRef] [Green Version]
  25. Xu, Q.; Wei, L. Prognostic equation for radar radial velocity derived by considering atmospheric refraction and earth curvature. J. Atmos. Sci. 2013, 70, 3328–3338. [Google Scholar] [CrossRef]
  26. Sun, J.; Crook, N.A. Dynamical and microphysical retrieval from Doppler radar observations using a cloud model and its adjoint. Part I: Model development and simulated data experiments. Atmos. Sci. 1997, 54, 1642–1661. [Google Scholar] [CrossRef]
  27. Wang, H.; Sun, J.; Zhang, X.; Huang, X.-Y.; Auligné, T. Radar data assimilation with WRF 4D-Var. Part I: System development and preliminary testing. Mon. Weather Rev. 2013, 141, 2224–2244. [Google Scholar] [CrossRef]
  28. Xu, Q.; Wei, L.; Gu, W.; Gong, J.; Zhao, Q. A 3.5-dimensional variational method for Doppler radar data assimilation and its application to phased-array radar observations. Adv. Meteorol. 2010, 2010, 797265. [Google Scholar] [CrossRef] [Green Version]
  29. Wang, X. Application of the WRF Hybrid ETKF–3DVAR Data Assimilation System for Hurricane Track Forecasts. Weather Forecast. 2011, 26, 868–884. [Google Scholar] [CrossRef]
  30. Li, X.; Ming, J.; Xue, M.; Wang, Y.; Zhao, K. Implementation of a dynamic equation constraint based on the steady state momentum equations within the WRF hybrid ensemble-3DVar data assimilation system and test with radar T-TREC wind assimilation for tropical Cyclone Chanthu (2010). J. Geophys. Res. Atmos. 2015, 120, 4017–4039. [Google Scholar] [CrossRef] [Green Version]
Figure 1. (a) Geopotential height (contour, unit: dagpm), wind field (barb, unit: m/s), and water vapor flux (shaded, unit: g/(cm·hPa·s)) at 850 hPa, (b) Height field (contour, unit: dagpm) and wind field (barb, unit: m/s) at 500 hPa at 0000 UTC on 13 September 2008.
Figure 1. (a) Geopotential height (contour, unit: dagpm), wind field (barb, unit: m/s), and water vapor flux (shaded, unit: g/(cm·hPa·s)) at 850 hPa, (b) Height field (contour, unit: dagpm) and wind field (barb, unit: m/s) at 500 hPa at 0000 UTC on 13 September 2008.
Atmosphere 12 00853 g001
Figure 2. The WRF simulation domain overlapped with the Houston, Texas (KHGX), and Lake Charles, Louisiana (KLCH), Weather Surveillance Radar 88 Doppler (WSR-88D) locations.
Figure 2. The WRF simulation domain overlapped with the Houston, Texas (KHGX), and Lake Charles, Louisiana (KLCH), Weather Surveillance Radar 88 Doppler (WSR-88D) locations.
Atmosphere 12 00853 g002
Figure 3. Distribution of conventional data (a) at 0000 UTC on 13 September; (b) at 0600 UTC on 13 September. SATOB (satellite atmospheric motion vectors), AIREP (aircraft data), SOUND (radiosondes).
Figure 3. Distribution of conventional data (a) at 0000 UTC on 13 September; (b) at 0600 UTC on 13 September. SATOB (satellite atmospheric motion vectors), AIREP (aircraft data), SOUND (radiosondes).
Atmosphere 12 00853 g003
Figure 4. Flow chart of the experiment with data assimilation conducted every 30 min.
Figure 4. Flow chart of the experiment with data assimilation conducted every 30 min.
Atmosphere 12 00853 g004
Figure 5. The 700 hPa wind analysis increments (m/s) at 0000 UTC 13 Sep 2008 from experiment DA_Radar.
Figure 5. The 700 hPa wind analysis increments (m/s) at 0000 UTC 13 Sep 2008 from experiment DA_Radar.
Atmosphere 12 00853 g005
Figure 6. Scatter of radar radial velocity in (a) observation and background, (b) in observation and analysis for experiment DA_Radar at 0000 UTC 13 September 2008.
Figure 6. Scatter of radar radial velocity in (a) observation and background, (b) in observation and analysis for experiment DA_Radar at 0000 UTC 13 September 2008.
Atmosphere 12 00853 g006
Figure 7. RMSE before and after the data assimilation from 0000 UTC on 13 September to 0600 UTC on September 13 for the 13 DA cycles for experiment DA_Radar_CVH.
Figure 7. RMSE before and after the data assimilation from 0000 UTC on 13 September to 0600 UTC on September 13 for the 13 DA cycles for experiment DA_Radar_CVH.
Atmosphere 12 00853 g007
Figure 8. DA_Radar (a,c,e,g) and DA_Radar_CVH (b,d,f,h) water vapor increment (unit: kg/m2) of qv (line 1), qr (line 2), qg (line 3), qs (line 4) at 0300 UTC on 13 September.
Figure 8. DA_Radar (a,c,e,g) and DA_Radar_CVH (b,d,f,h) water vapor increment (unit: kg/m2) of qv (line 1), qr (line 2), qg (line 3), qs (line 4) at 0300 UTC on 13 September.
Atmosphere 12 00853 g008
Figure 9. RMS profiles of the analysis increment for each variable at (a) 0000 UTC, 13 September (b) 0300 UTC, 13 September (c) 0600 UTC, 13 September (U, V in m/s, t in K, and other hydrometeors in 10 g/kg) for experiment DA_Radar_CVH.
Figure 9. RMS profiles of the analysis increment for each variable at (a) 0000 UTC, 13 September (b) 0300 UTC, 13 September (c) 0600 UTC, 13 September (U, V in m/s, t in K, and other hydrometeors in 10 g/kg) for experiment DA_Radar_CVH.
Atmosphere 12 00853 g009
Figure 10. The 3 h accumulated precipitation in (a) observation, (b) DA_CONV experiment, (c) DA_Radar experiment, (d) DA_Radar_CVH experiment from 0600 UTC, 13 September to 0900 UTC, 13 September.
Figure 10. The 3 h accumulated precipitation in (a) observation, (b) DA_CONV experiment, (c) DA_Radar experiment, (d) DA_Radar_CVH experiment from 0600 UTC, 13 September to 0900 UTC, 13 September.
Atmosphere 12 00853 g010
Figure 11. TS score of 3 h accumulated precipitation from 0600 UTC, 13 September to 0900 UTC, 13 September.
Figure 11. TS score of 3 h accumulated precipitation from 0600 UTC, 13 September to 0900 UTC, 13 September.
Atmosphere 12 00853 g011
Figure 12. (a) track, (b) track error from 0600 UTC, 13 September to 0000 UTC, 14 September.
Figure 12. (a) track, (b) track error from 0600 UTC, 13 September to 0000 UTC, 14 September.
Atmosphere 12 00853 g012
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Shen, F.; Min, J.; Li, H.; Xu, D.; Shu, A.; Zhai, D.; Guo, Y.; Song, L. Applications of Radar Data Assimilation with Hydrometeor Control Variables within the WRFDA on the Prediction of Landfalling Hurricane IKE (2008). Atmosphere 2021, 12, 853. https://0-doi-org.brum.beds.ac.uk/10.3390/atmos12070853

AMA Style

Shen F, Min J, Li H, Xu D, Shu A, Zhai D, Guo Y, Song L. Applications of Radar Data Assimilation with Hydrometeor Control Variables within the WRFDA on the Prediction of Landfalling Hurricane IKE (2008). Atmosphere. 2021; 12(7):853. https://0-doi-org.brum.beds.ac.uk/10.3390/atmos12070853

Chicago/Turabian Style

Shen, Feifei, Jinzhong Min, Hong Li, Dongmei Xu, Aiqing Shu, Danhua Zhai, Yakai Guo, and Lixin Song. 2021. "Applications of Radar Data Assimilation with Hydrometeor Control Variables within the WRFDA on the Prediction of Landfalling Hurricane IKE (2008)" Atmosphere 12, no. 7: 853. https://0-doi-org.brum.beds.ac.uk/10.3390/atmos12070853

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop