scholarly journals Improving Large-Domain Convection-Allowing Forecasts with High-Resolution Analyses and Ensemble Data Assimilation

2016 ◽  
Vol 144 (5) ◽  
pp. 1777-1803 ◽  
Author(s):  
Craig S. Schwartz

Analyses with 20-km horizontal grid spacing were produced from continuously cycling three-dimensional variational (3DVAR), ensemble square root Kalman filter (EnSRF), and “hybrid” variational–ensemble data assimilation (DA) systems over a domain spanning the conterminous United States. These analyses initialized 36-h Weather Research and Forecasting Model forecasts containing a large convection-allowing 4-km nested domain, where downscaled 20-km 3DVAR, EnSRF, and hybrid analyses initialized the 4-km forecasts. Overall, hybrid analyses initialized the best 4-km precipitation forecasts. Furthermore, whether 4-km precipitation forecasts could be improved by initializing them with true 4-km analyses was assessed. As it was computationally infeasible to produce 4-km continuously cycling ensembles over the large 4-km domain, several “dual-resolution” hybrid DA configurations were adopted where 4-km backgrounds were combined with 20-km ensembles to produce 4-km hybrid analyses. Additionally, 4-km 3DVAR analyses were produced. In both hybrid and 3DVAR frameworks, initializing 4-km forecasts with true 4-km analyses, rather than downscaled 20-km analyses, yielded superior precipitation forecasts over the first 12 h. Differences between forecasts initialized from 4-km and downscaled 20-km hybrid analyses were smaller for 18–36-h forecasts, but there were occasionally meaningful differences. Continuously cycling the 4-km backgrounds and using static background error covariances with larger horizontal length scales in the hybrid led to better forecasts. All hybrid-initialized forecasts, including those initialized from downscaled 20-km analyses, were more skillful than forecasts initialized from 4-km 3DVAR analyses, suggesting the analysis method was more important than analysis resolution.

2013 ◽  
Vol 141 (12) ◽  
pp. 4350-4372 ◽  
Author(s):  
Craig S. Schwartz ◽  
Zhiquan Liu ◽  
Xiang-Yu Huang ◽  
Ying-Hwa Kuo ◽  
Chin-Tzu Fong

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.


2014 ◽  
Vol 142 (2) ◽  
pp. 716-738 ◽  
Author(s):  
Craig S. Schwartz ◽  
Zhiquan Liu

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.


2011 ◽  
Vol 2011 ◽  
pp. 1-10 ◽  
Author(s):  
Chien-Ben Chou ◽  
Huei-Ping Huang

This work assesses the effects of assimilating atmospheric infrared sounder (AIRS) observations on typhoon prediction using the three-dimensional variational data assimilation (3DVAR) and forecasting system of the weather research and forecasting (WRF) model. Two major parameters in the data assimilation scheme, the spatial decorrelation scale and the magnitude of the covariance matrix of the background error, are varied in forecast experiments for the track of typhoon Sinlaku over the Western Pacific. The results show that within a wide parameter range, the inclusion of the AIRS observation improves the prediction. Outside this range, notably when the decorrelation scale of the background error is set to a large value, forcing the assimilation of AIRS data leads to degradation of the forecast. This illustrates how the impact of satellite data on the forecast depends on the adjustable parameters for data assimilation. The parameter-sweeping framework is potentially useful for improving operational typhoon prediction.


2014 ◽  
Vol 142 (6) ◽  
pp. 2309-2320 ◽  
Author(s):  
Erika L. Navarro ◽  
Gregory J. Hakim

Abstract A significant challenge for tropical cyclone ensemble data assimilation is that storm-scale observations tend to make analyses that are more asymmetric than the prior forecasts. Compromised structure and intensity, such as an increase of amplitude across the azimuthal Fourier spectrum, are a routine property of ensemble-based analyses, even with accurate position observations and frequent assimilation. Storm dynamics in subsequent forecasts evolve these states toward axisymmetry, creating difficulty in distinguishing between model-induced and actual storm asymmetries for predictability studies and forecasting. To address this issue, a novel algorithm using a storm-centered approach is proposed. The method is designed for use with existing ensemble filters with little or no modification, facilitating its adoption and maintenance. The algorithm consists of 1) an analysis of the environment using conventional coordinates, 2) a storm-centered analysis using storm-relative coordinates, and 3) a merged analysis that combines the large-scale and storm-scale fields together at an updated storm location. This algorithm is evaluated in two sets of observing system simulation experiments (OSSEs): first, no-cycling tests of the update step for idealized three-dimensional storms in radiative–convective equilibrium; second, full cycling tests of data assimilation applied to a shallow-water model for a field of interacting vortices. Results are compared against a control experiment based on a conventional ensemble Kalman filter (EnKF) scheme as well as an alternative EnKF scheme proposed by Lawson and Hansen. The storm-relative method yields vortices that are more symmetric and exhibit finer inner-core structure than either approach, with errors reduced by an order of magnitude over a control case with prior spread consistent with the National Hurricane Center (NHC)’s mean 5-yr forecast track error at 12 h. Azimuthal Fourier error spectra exhibit much-reduced noise associated with data assimilation as compared to both the control and the Lawson and Hansen approach. An assessment of free-surface height tendency of model forecasts after the merge step reveals a balanced trend between the storm-centered and conventional approaches, with storm-centered values more closely resembling the reference state.


2013 ◽  
Vol 141 (4) ◽  
pp. 1263-1284 ◽  
Author(s):  
Glen S. Romine ◽  
Craig S. Schwartz ◽  
Chris Snyder ◽  
Jeff L. Anderson ◽  
Morris L. Weisman

Abstract During the spring 2011 season, a real-time continuously cycled ensemble data assimilation system using the Advanced Research version of the Weather Research and Forecasting Model (WRF) coupled with the Data Assimilation Research Testbed toolkit provided initial and boundary conditions for deterministic convection-permitting forecasts, also using WRF, over the eastern two-thirds of the conterminous United States (CONUS). In this study the authors evaluate the mesoscale assimilation system and the convection-permitting forecasts, at 15- and 3-km grid spacing, respectively. Experiments employing different physics options within the continuously cycled ensemble data assimilation system are shown to lead to differences in the mean mesoscale analysis characteristics. Convection-permitting forecasts with a fixed model configuration are initialized from these physics-varied analyses, as well as control runs from 0.5° Global Forecast System (GFS) analysis. Systematic bias in the analysis background influences the analysis fit to observations, and when this analysis initializes convection-permitting forecasts, the forecast skill is degraded as bias in the analysis background increases. Moreover, differences in mean error characteristics associated with each physical parameterization suite lead to unique errors of spatial, temporal, and intensity aspects of convection-permitting rainfall forecasts. Observation bias by platform type is also shown to impact the analysis quality.


Author(s):  
Y. Hu ◽  
M. Zhang ◽  
Y. Liang ◽  
L. Ye ◽  
D. Zhao ◽  
...  

<p><strong>Abstract.</strong> Background error covariance (BEC) plays a key role in a variational data assimilation system. It determines variable analysis increments by spreading information from observation points. In order to test the influence of BEC on the GSI data assimilation and prediction of aerosol in Beijing-Tianjin-Hebei, a regional BEC is calculated using one month series of numerical forecast fields of November 2017 based on the National Meteorological Center (NMC) method, and compared with the global BEC.The results show that the standard deviation of stream function of the regional BEC is larger than that of the global BEC. And the horizontal length-scale of the regional BEC is smaller than that of the global BEC, white the vertical length-scale of the regional BEC is similar with that of the global BEC. The increments of the assimilation experiment with the regional BEC present more small scale information than that with the global BEC. The forecast skill of the experiment with the regional BEC is higher than that with the global BEC in the stations of Beijing, Tianjin, Chengde and Taiyuan, and the average root-mean-square errors (RMSE) reduces by over 13.4%.</p>


2013 ◽  
Vol 141 (2) ◽  
pp. 754-772 ◽  
Author(s):  
Sara Q. Zhang ◽  
Milija Zupanski ◽  
Arthur Y. Hou ◽  
Xin Lin ◽  
Samson H. Cheung

Abstract Assimilation of remotely sensed precipitation observations into numerical weather prediction models can improve precipitation forecasts and extend prediction capabilities in hydrological applications. This paper presents a new regional ensemble data assimilation system that assimilates precipitation-affected microwave radiances into the Weather Research and Forecasting Model (WRF). To meet the challenges in satellite data assimilation involving cloud and precipitation processes, hydrometeors produced by the cloud-resolving model are included as control variables and ensemble forecasts are used to estimate flow-dependent background error covariance. Two assimilation experiments have been conducted using precipitation-affected radiances from passive microwave sensors: one for a tropical storm after landfall and the other for a heavy rain event in the southeastern United States. The experiments examined the propagation of information in observed radiances via flow-dependent background error auto- and cross covariance, as well as the error statistics of observational radiance. The results show that ensemble assimilation of precipitation-affected radiances improves the quality of precipitation analyses in terms of spatial distribution and intensity in accumulated surface rainfall, as verified by independent ground-based precipitation observations.


2020 ◽  
Vol 10 (24) ◽  
pp. 9010
Author(s):  
Sujeong Lim ◽  
Myung-Seo Koo ◽  
In-Hyuk Kwon ◽  
Seon Ki Park

Ensemble data assimilation systems generally suffer from underestimated background error covariance that leads to a filter divergence problem—the analysis diverges from the natural state by ignoring the observation influence due to the diminished estimation of model uncertainty. To alleviate this problem, we have developed and implemented the stochastically perturbed hybrid physical–dynamical tendencies to the local ensemble transform Kalman filter in a global numerical weather prediction model—the Korean Integrated Model (KIM). This approach accounts for the model errors associated with computational representations of underlying partial differential equations and the imperfect physical parameterizations. The new stochastic perturbation hybrid tendencies scheme generally improved the background error covariances in regions where the ensemble spread was not sufficiently expressed by the control experiment that used an additive inflation and the relaxation to prior spread method.


Sign in / Sign up

Export Citation Format

Share Document