Convective-Scale Data Assimilation and Precipitation Prediction with a Local Ensemble Transform Kalman Filter Radar Assimilation System Over Complex Terrain: A Thorough Investigation with the Heavy Rainfall in Taiwan on 16 June 2008

2021 ◽  
pp. 543-579
Author(s):  
Shu-Chih Yang ◽  
Hsiang-Wen Cheng ◽  
Pin-Ying Wu ◽  
Zih-Mao Huang ◽  
Chih-Chien Tsai
2017 ◽  
Author(s):  
Yun Liu ◽  
Eugenia Kalnay ◽  
Ning Zeng ◽  
Ghassem Asrar ◽  
Zhaohui Chen ◽  
...  

Abstract. We developed a Carbon data assimilation system to estimate the surface carbon fluxes using the Local Ensemble Transform Kalman Filter and atmospheric transfer model of GEOS-Chem driven by the MERRA-1 reanalysis of the meteorological fields based on the Goddard Earth Observing System Model, Version 5 (GEOS-5). This assimilation system is inspired by the method of Kang et al. [2011, 2012], who estimated the surface carbon fluxes in an Observing System Simulation Experiment (OSSE) mode, as evolving parameters in the assimilation of the atmospheric CO2, using a short assimilation window of 6 hours. They included the assimilation of the standard meteorological variables, so that the ensemble provided a measure of the uncertainty in the CO2 transport. After introducing new techniques such as variable localization, and increased observation weights near the surface, they obtained accurate carbon fluxes at grid point resolution. We developed a new version of the LETKF related to the Running-in-Place (RIP) method used to accelerate the spin-up of EnKF data assimilation [Kalnay and Yang, 2010; Wang et al., 2013, Yang et al., 2014]. Like RIP, the new assimilation system uses the no-cost smoothing algorithm for the LETKF [Kalnay et al., 2007b], which allows shifting at no cost the Kalman Filter solution forward or backward within an assimilation window. In the new scheme a long observation window (e.g., 7-days or longer) is used to create an LETKF ensemble at 7-days. Then, the RIP smoother is used to obtain an accurate final analysis at 1-day. This analysis has the advantage of being based on a short assimilation window, which makes it more accurate, and of having been exposed to the future 7-days observations, which accelerates the spin up. The assimilation and observation windows are then shifted forward by one day, and the process is repeated. This reduces significantly the analysis error, suggesting that this method could be used in other data assimilation problems.


2008 ◽  
Vol 60 (1) ◽  
pp. 113-130 ◽  
Author(s):  
Istvan Szunyogh ◽  
Eric J. Kostelich ◽  
Gyorgyi Gyarmati ◽  
Eugenia Kalnay ◽  
Brian R. Hunt ◽  
...  

2017 ◽  
Author(s):  
Yun Liu ◽  
Eugenia Kalnay ◽  
Ning Zeng ◽  
Ghassem Asrar ◽  
Zhaohui Chen ◽  
...  

Abstract. We developed a Carbon data assimilation system to estimate the surface carbon fluxes using the Local Ensemble Transform Kalman Filter and atmospheric transfer model of GEOS-Chem driven by the MERRA-1 reanalysis of the meteorological fields based on the Goddard Earth Observing System Model, Version 5 (GEOS-5). This assimilation system is inspired by the method of Kang et al. (2011, 2012), who estimated the surface carbon fluxes in an Observing System Simulation Experiment (OSSE) mode, as evolving parameters in the assimilation of the atmospheric CO2, using a short assimilation window of 6 hours. They included the assimilation of the standard meteorological variables, so that the ensemble provided a measure of the uncertainty in the CO2 transport. After introducing new techniques such as “variable localization”, and increased observation weights near the surface, they obtained accurate carbon fluxes at grid point resolution. We developed a new version of the LETKF related to the “Running-in-Place” (RIP) method used to accelerate the spin-up of EnKF data assimilation (Kalnay and Yang, 2010; Wang et al., 2013, Yang et al., 2014). Like RIP, the new assimilation system uses the “no-cost smoothing” algorithm for the LETKF (Kalnay et al., 2007b), which allows shifting at no cost the Kalman Filter solution forward or backward within an assimilation window. In the new scheme a long “observation window” (e.g., 7-days or longer) is used to create an LETKF ensemble at 7-days. Then, the RIP smoother is used to obtain an accurate final analysis at 1-day. This analysis has the advantage of being based on a short assimilation window, which makes it more accurate, and of having been exposed to the future 7-days observations, which accelerates the spin up. The assimilation and observation windows are then shifted forward by one day, and the process is repeated. This reduces significantly the analysis error, suggesting that this method could be used in other data assimilation problems. The newly developed assimilation method can be used with other Earth system models, especially for greater use of observations in conjunction with models.


2020 ◽  
Vol 148 (9) ◽  
pp. 3631-3652
Author(s):  
Pin-Ying Wu ◽  
Shu-Chih Yang ◽  
Chih-Chien Tsai ◽  
Hsiang-Wen Cheng

ABSTRACT Sampling error stems from the use of ensemble-based data assimilation (EDA) with a limited ensemble size and can result in spurious background error covariances, leading to false analysis corrections. The WRF-LETKF radar assimilation system (WLRAS) is performed separately with 256 and 40 members to investigate the characteristics of convective-scale sampling errors in the EDA and its impact on precipitation prediction based on a heavy rainfall event on 16 June 2008. The results suggest that the sampling errors for this event are sensitive to the relationships between the simulated observations and model variables, the intensity of reflectivity, and how the prevailing wind projects to the radial wind in the areas that the radar cannot resolve U or V wind. The sampling errors lead to an underprediction of heavy rainfall when the horizontal localization radius is inadequately large, but this can be mitigated when a more accurate moisture condition is provided. In addition, being able to use a larger vertical localization plays an important role in providing necessary adjustments for representing the vertical thermodynamical structure of convection, which further improves precipitation prediction. A strategy mitigating the impact of sampling errors associated with the limitation of radial wind measurement by inflating the observation error over sensitive areas can bring benefits to precipitation prediction.


2017 ◽  
Vol 145 (11) ◽  
pp. 4575-4592 ◽  
Author(s):  
Craig H. Bishop ◽  
Jeffrey S. Whitaker ◽  
Lili Lei

To ameliorate suboptimality in ensemble data assimilation, methods have been introduced that involve expanding the ensemble size. Such expansions can incorporate model space covariance localization and/or estimates of climatological or model error covariances. Model space covariance localization in the vertical overcomes problematic aspects of ensemble-based satellite data assimilation. In the case of the ensemble transform Kalman filter (ETKF), the expanded ensemble size associated with vertical covariance localization would also enable the simultaneous update of entire vertical columns of model variables from hyperspectral and multispectral satellite sounders. However, if the original formulation of the ETKF were applied to an expanded ensemble, it would produce an analysis ensemble that was the same size as the expanded forecast ensemble. This article describes a variation on the ETKF called the gain ETKF (GETKF) that takes advantage of covariances from the expanded ensemble, while producing an analysis ensemble that has the required size of the unexpanded forecast ensemble. The approach also yields an inflation factor that depends on the localization length scale that causes the GETKF to perform differently to an ensemble square root filter (EnSRF) using the same expanded ensemble. Experimentation described herein shows that the GETKF outperforms a range of alternative ETKF-based solutions to the aforementioned problems. In cycling data assimilation experiments with a newly developed storm-track version of the Lorenz-96 model, the GETKF analysis root-mean-square error (RMSE) matches the EnSRF RMSE at shorter than optimal localization length scales but is superior in that it yields smaller RMSEs for longer localization length scales.


Sign in / Sign up

Export Citation Format

Share Document