Hybrid Data Assimilation System Research Articles

Assimilating AMSU-A Radiance Observations with an Ensemble Four-Dimensional Variational (En4DVar) Hybrid Data Assimilation System

Many ensemble-based data assimilation (DA) methods use observation space localization to mitigate the sampling errors due to the insufficient ensemble members. Observation space localization is simpler and more timesaving than model space localization in implementation, but more difficult to directly assimilate satellite radiance observations, a kind of non-local observations. The vertical locations of radiance observations are undetermined and the transmission of observational information is thereby obstructed. To determine the vertical coordinates of radiance observations, a weighted average hypsometry is proposed. Using this hypsometry, AMSU-A radiance observations are directly assimilated with an ensemble four-dimensional variational (En4DVar) DA system. It consists of a four-dimensional ensemble-variational (4DEnVar) system providing ensemble covariance and a 4DVar system. Observing system simulation experiments show that the hypsometry alleviates the degradations in the late period of medium-range forecast in the Northern Extratropics that occur in the traditional peak-based hypsometry. It obviously improves the analysis qualities and forecast skills of the En4DVar system and its two components, especially in the Southern Extratropics, when incorporating AMSU-A radiance observations. The improvement in the En4DVar-initialized forecast is comparable to that in the 4DVar-initialized forecast in the Southern Extratropics and Tropics. It indicates that a proper hypsometry enables efficient extraction of useful information from AMSU-A radiance observations by 4DEnVar with observation space localization. Therefore, the 4DEnVar provides high-quality ensemble covariances for En4DVar.

The impacts of assimilating Aeolus horizontal line-of-sight winds on numerical predictions of Hurricane Ida (2021) and a mesoscale convective system over the Atlantic Ocean

Abstract. The Aeolus satellite, the first spaceborne wind lidar launched by the European Space Agency (ESA) on 22 August 2018, provides global measurements of horizontal line-of-sight (HLOS) winds. The assimilation of Aeolus HLOS winds has been proven to improve numerical weather predictions (NWPs). Still, its influences on forecasts of tropical cyclones (TCs) and tropical convective systems have yet to be examined in detail. This study investigates the impacts of assimilating Aeolus HLOS winds on the analysis and forecasts of Hurricane Ida (2021) and a mesoscale convective system (MCS) embedded in an African easterly wave (AEW) over the Atlantic Ocean (AO) with the mesoscale community Weather Research and Forecasting (WRF) model and the NCEP Gridpoint Statistical Interpolation (GSI)-based three-dimensional ensemble-variational (3DEnVAR) hybrid data assimilation (DA) system. Mie-cloudy and Rayleigh-clear winds are assimilated. The results for Ida (2021) show that assimilating Aeolus HLOS winds leads to better track predictions. The intensity forecasts are improved in some cases, even with limited coverage of Aeolus HLOS winds within the inner core region of Ida (2021). In addition, the structure of heavy precipitation associated with Ida (2021) is refined after the assimilation of Aeolus HLOS winds. Further diagnosis demonstrates that the improved intensity and precipitation forecasts result from enhanced divergence in the upper level of the troposphere after the assimilation of Aeolus HLOS winds. Additional results from the MCS associated with an AEW indicate that assimilating Aeolus HLOS winds enhances forecasts of its precipitation structure and the associated low-level divergence. Findings from this study suggest that the assimilation of Aeolus HLOS winds has the potential to improve forecasts for TCs and tropical convective systems.

The Use of Regional Data Assimilation to Improve Numerical Simulations of Diurnal Characteristics of Precipitation during an Active Madden–Julian Oscillation Event over the Maritime Continent

This study examines the impact of regional data assimilation on diurnal characteristics of precipitation and winds over the Maritime Continent (MC) using a set of cloud-permitting-scale (~3 km) numerical simulations with the mesoscale community Weather Research and Forecasting (WRF) model and the NCEP Gridpoint Statistical Interpolation (GSI)-based ensemble-3DVAR hybrid data assimilation system. Numerical experiments focus on January 2018, when a well-defined, active Madden–Julian Oscillation (MJO) propagated through the MC region. Available conventional and satellite data are assimilated. Results show that simulated precipitation with data assimilation generally agrees better with satellite-derived rainfall than the control simulation without data assimilation. Simulations with data assimilation also reproduce the diurnal cycle of precipitation better, especially for the timing of the precipitation peak. Data assimilation modulates the overstrong (overweak) diurnal forcing over the land (ocean) in the control simulation. The vertical phase shift of the thermodynamic environment, associated with the timing of vertical motion transition along with low-level water vapor supplies, results in maximum precipitation occurring later, especially over land. To further demonstrate the impact of data assimilation, an additional experiment assimilates NASA Cyclone Global Navigation Satellite System (CYGNSS)-derived ocean surface winds. The results indicate that the assimilation of CYGNSS data exhibits an evident impact on the diurnal variation of surface variables and a similar shift in the diurnal cycle of precipitation. Overall, this study highlights the importance of regional data assimilation in improving the representation of precipitation over the MC, paving the way for a better understanding of the interactions of local diurnal convective precipitation cycles with MJO.

A hybrid data assimilation system based on machine learning

In the earth sciences, numerical weather prediction (NWP) is the primary method of predicting future weather conditions, and its accuracy is affected by the initial conditions. Data assimilation (DA) can provide high-precision initial conditions for NWP. The hybrid 4DVar-EnKF is currently an advanced DA method used by many operational NWP centres. However, it has two major shortcomings: The complex development and maintenance of the tangent linear and adjoint models and the empirical combination of the results of 4DVar and EnKF. In this paper, a new hybrid DA method based on machine learning (HDA-ML) is presented to overcome these drawbacks. In the new method, the tangent linear and adjoint models in the 4DVar part of the hybrid algorithm can be easily obtained by using a bilinear neural network to replace the forecast model, and a CNN model is adopted to fuse the analysis of 4DVar and EnKF to adaptively obtain the optimal coefficient of combination rather than the empirical coefficient as in the traditional hybrid DA method. The hybrid DA methods are compared with the Lorenz-96 model using the true values as labels. The experimental results show that HDA-ML improves the assimilation performance and significantly reduces the time cost. Furthermore, using observations instead of the true values as labels in the training system is more realistic. The results show comparable assimilation performance to that in the experiments with the true values used as the labels. The experimental results show that the new method has great potential for application to operational NWP systems.

A Bias Correction Scheme with the Symmetric Cloud Proxy Variable and Its Influence on Assimilating All-Sky GOES-16 Brightness Temperatures

Abstract All-sky assimilation of brightness temperatures (BTs) from GOES-16 infrared water vapor channels is challenging, primarily because these channels are sensitive to cloud ice that causes large nonlinear errors in the forecast and forward models. Thus, bias correction (BC) for all-sky assimilation of GOES-16 BTs is vital. This study examines the impacts of different BC schemes, especially for a scheme with a quartic polynomial of cloud predictors (the ASRBC4 scheme), on the analysis and WRF Model forecasts of tropical cyclones when assimilating the all-sky GOES-16 channel-8 BTs using the NCEP GSI-based 3D ensemble–variational hybrid data assimilation (DA) system with variational BC (VarBC). Long-term statistics are performed during the NASA Convective Processes Experiment field campaign (2017). Results demonstrate that the ASRBC4 scheme effectively reduces the average of all-sky scaled observation-minus-backgrounds (OmBs) in a cloudy sky and alleviates their nonlinear conditional biases with respect to the symmetric cloud proxy variable, in contrast to the BC schemes without the cloud predictor or with a first-order cloud predictor. In addition, adopting the ASRBC4 scheme in DA decreases the positive temperature increments at 200 hPa and the accompanying midlevel cyclonic wind increments in the analysis of Tropical Storm (TS) Cindy (2017). Applying the ASRBC4 scheme also leads to better storm-track predictions for TS Cindy (2017) and Hurricane Laura (2022), compared to experiments with other BC schemes. Overall, this study highlights the importance of reducing nonlinear biases of OmBs in a cloudy sky for successful all-sky assimilation of BTs from GOES-16 infrared water vapor channels.

Hybrid ensemble-variational data assimilation in ABC-DA within a tropical framework

Abstract. Hybrid ensemble-variational data assimilation (DA) methods have gained significant traction in recent years. These methods aim to alleviate the limitations and maximise the advantages offered by ensemble or variational methods. Most existing hybrid applications focus on the mid-latitudinal context; almost none have explored its benefits in the tropical context. In this article, hybrid ensemble-variational DA is introduced to a tropical configuration of a simplified non-hydrostatic convective-scale fluid dynamics model (the ABC model, named after its three key parameters: the pure gravity wave frequency A, the controller of the acoustic wave speed B, and the constant of proportionality between pressure and density perturbations C), and its existing variational framework, the ABC-DA system. The hybrid ensemble-variational DA algorithm is developed based on the alpha control variable approach, often used in numerical weather prediction. Aspects of the algorithm such as localisation (used to mitigate sampling error caused by finite ensemble sizes) and weighting parameters (used to weight the ensemble and climatological contributions to the background error covariance matrix) are implemented. To produce the flow-dependent error modes (ensemble perturbations) for the ensemble-variational DA algorithm, an ensemble system is also designed for the ABC model which is run alongside the hybrid DA system. A random field perturbations method is used to generate an initial ensemble which is then propagated using the ensemble bred vectors method. This setup allows the ensemble to be centred on the hybrid control analysis. Visualisation software has been developed to focus on the diagnosis of the ensemble system. To demonstrate the hybrid ensemble-variational DA in the ABC-DA system, sensitivity tests using observing system simulation experiments are conducted within a tropical framework. A 30-member ensemble was used to generate the error modes for the experiments. In general, the best performing configuration (with respect to the “truth”) for the hybrid ensemble-variational DA system used an 80%/20% weighting on the ensemble-derived/climatological background error covariance matrix contributions. For the horizontal wind variables though, full weight on the ensemble-derived background error covariance matrix (100%/0%) resulted in the smallest cycle-averaged analysis root mean square errors, mainly due to large errors in the meridional wind field when contributions from the climatological background error covariance matrix were involved, possibly related to a sub-optimal background error covariance model. The ensemble bred vectors method propagated a healthy-looking DA-centred ensemble without bimodalities or evidence of filter collapse. The ensemble was under-dispersive for some variables but for others, the ensemble spread approximately matched the corresponding root mean square errors. Reducing the number of ensemble members led to slightly larger errors across all variables due to the introduction of larger sampling errors into the system.

A 4DEnVar‐Based Ensemble Four‐Dimensional Variational (En4DVar) Hybrid Data Assimilation System for Global NWP: System Description and Primary Tests

AbstractThis study developed an ensemble four‐dimensional variational (En4DVar) hybrid data assimilation system. Different from most of the available En4DVar systems that adopted ensemble Kalman Filter class or ensemble data assimilation approaches to produce ensemble covariances for their hybrid background error covariances (BECs), it used a four‐dimensional ensemble variational (4DEnVar) system to obtain the ensemble covariance. The localization scheme for 4DEnVar applied orthogonal functions to decompose the correlation matrix so that it was implemented easily and rapidly. In terms of analysis quality and forecast skill, the En4DVar system was evaluated in the single‐point observation experiments and observing system simulation experiments (OSSEs) with sounding and cloud‐derived wind observations, using its standalone four‐dimensional variational (4DVar) and 4DEnVar components as references. The single‐point observation experiments visually verified the explicit flow‐dependent characteristic of the BEC due to the introduction of the ensemble covariance from the 4DEnVar system. The OSSE‐based sensitivity experiments revealed different contributions of the weight for the ensemble covariance in the En4DVar system to the forecasts in the Northern and Southern Extratropics and Tropics. A much higher weight for the ensemble covariance in a properly inflated hybrid covariance helped En4DVar produce the most reasonable analysis. The forecast initialized by En4DVar is overall better than by 4DVar and 4DEnVar, although the quality of En4DVar analysis is between those of 4DVar and 4DEnVar ensemble mean analyses. It indicates that the flow‐dependent ensemble covariance provided by 4DEnVar dominantly contributes to the improvements in the En4DVar‐initialized forecast, with certain but necessary constraint from the balanced climatological covariance.

Impact of Ensemble-Variational Data Assimilation in Heavy Rain Forecast over Brazilian Northeast

The Brazilian Northeast (BNE) is located in the tropical region of Brazil. It is bounded by the Atlantic Ocean, and its climate and vegetation are strongly affected by continental plateaus. The plateaus keep the humid air masses to the east and are responsible for the rain episodes, and at the west side the northeastern hinterland and dry air masses are observed. This work is a case study that aims to evaluate the impact of updating the model initial condition using the 3DEnVar (Three-Dimensional Ensemble Variational) system in heavy rain episodes associated with Mesoscale Convective Systems (MCS). The results were compared to 3DVar (Three-Dimensional Variational) and EnSRF (Ensemble Square Root Filter) systems and with no data assimilation. The study enclosed two MCS cases occurring on 14 and 24 January 2017. For that purpose, the RMS (Regional Modeling System) version 3.0.0, maintained by the Center for Weather Forecasting and Climate Studies (CPTEC), used two components: the Weather Research and Forecasting (WRF) mesoscale model and the GSI (Gridpoint Statistical Interpolation) data assimilation system. Currently, the RMS provides the WRF initial conditions by using 3DVar data assimilation methodology. The 3DVar uses a climatological covariance matrix to minimize model errors. In this work, the 3DEnVar updates the RMS climatological covariance matrix through the forecast members based on the errors of the day. This work evaluated the improvements in the detection and estimation of 24 h accumulated precipitation in MCS events. The statistic index RMSE (Root Mean Square Error) showed that the hybrid data assimilation system (3DEnVar) performed better in reproducing the precipitation in the MCS occurred on 14 January 2017. On 24 January 2017, the EnSRF was the best system for improving the WRF forecast. In general, the BIAS showed that the WRF initialized with different initial conditions overestimated the 24 h accumulated precipitation. Therefore, the viability of using a hybrid system may depend on the hybrid algorithm that can modify the weights attributed to the EnSRF and 3DVar matrix in the GSI over the assimilation cycles.

Comparison of the performance of HYBRID ETKF-3DVAR and 3DVAR data assimilation scheme on the forecast of tropical cyclones formed over the Bay of Bengal

This study compares the performance of hybrid ensemble transform Kalman filter – three dimensional variational data assimilation (HYBRID) system and three dimensional variational (3DVAR) data assimilation system in Weather Research and Forecasting Model (WRF) in simulating tropical cyclones (TC) formed over the Bay of Bengal. An Ensemble Transform Kalman Filter (ETKF) system updates the ensemble system that provides flow-evolving background error covariance for HYBRID data assimilation system. Results indicate that use of flow-evolving ensemble error covariance in 3DVAR system has systematically reduced the TC position and intensity errors in the analysis; however, adding more weights to the ensemble error covariance term in 3DVAR cost function has not made any significant impact. The 3DVAR analysis depicts a stronger TC vortex with a well pronounced warm core structure as compared to that in HYBRID analysis. The forecasts from HYBRID analysis outperform that from 3DVAR in reducing TC track forecast error. The relative improvement in TC landfall position is 43% and 49% for variously configured HYBRID experiments. The forecasts initiated from HYBRID analysis has higher skill in quantitative precipitation forecasts during TC landfall compared to 3DVAR, which may be attributed to improved track prediction in the HYBRID experiments.

Evaluation of the Four-Dimensional Ensemble-Variational Hybrid Data Assimilation with Self-Consistent Regional Background Error Covariance for Improved Hurricane Intensity Forecasts

The feasibility of a hurricane initialization framework based on the Gridpoint Statistical Interpolation (GSI)-based four-dimensional ensemble-variational (GSI-4DEnVar) hybrid data assimilation system for the Hurricane Weather Research and Forecasting model (HWRF) model is evaluated in this study. The system considers the temporal evolution of error covariances via the use of four-dimensional ensemble perturbations that are provided by high-resolution, self-consistent HWRF ensemble forecasts. It is different from the configuration of the GSI-based three-dimensional ensemble-variational (GSI-3DEnVar) hybrid data assimilation system, similar to that used in the operational HWRF, which employs background error covariances provided by coarser-resolution global ensembles from the National Centers for Environmental Prediction (NCEP) Global Forecast System (GFS) ensemble Kalman filtering data assimilation system. In addition, our proposed initialization framework discards the empirical intensity correction in the vortex initialization package that is employed by the GSI-3DEnVar initialization framework in operational HWRF. Data assimilation and numerical simulation experiments for Hurricanes Joaquin (2015), Patricia (2015), and Matthew (2016) are conducted during their intensity changes. The impacts of two initialization frameworks on the HWRF analyses and forecasts are compared. It is found that GSI-4DEnVar leads to a reduction in track, minimum sea level pressure (MSLP), and maximum surface wind (MSW) forecast errors in all of the HWRF simulations, compared with the GSI-3DEnVar initialization framework. With assimilating high-resolution observations within the hurricane inner-core region, GSI-4DEnVar can produce the initial hurricane intensity reasonably well without the empirical vortex intensity correction. Further diagnoses with Hurricane Joaquin indicate that GSI-4DEnVar can significantly alleviate the imbalances in the initial conditions and enhance the performance of the data assimilation and subsequent hurricane intensity and precipitation forecasts.

All‐sky microwave humidity sounder assimilation in the Korean Integrated Model forecast system

AbstractThe Korean Integrated Model (KIM) and hybrid data assimilation system were extended to assimilate all‐sky radiance from microwave satellite sensors. Initially, a radiative transfer model for the TIROS Operational Vertical Sounder (TOVS), called RTTOV‐SCATT (version 11.3), was implemented to assimilate the Microwave Humidity Sounder (MHS) 183 GHz channels over the ocean. Cloud and precipitation parameters are not directly assimilated into our system, but temperature and humidity profiles are improved in the all‐sky assimilation. In the cycled analysis and forecast experiments, an assimilation of the MHS in a cloudy region shows globally substantial benefits: the coverage of the MHS radiance data is increased by 23–28% in the all‐sky assimilation. It is demonstrated that RTTOV‐SCATT well describes the brightness temperature (TB) on selected cloudy pixels over the Tropics and the Southern Hemisphere, with potential for linking the TB innovation to the increment for the hydrometeor variable. As a result, a reduction of 1.11% in the specific humidity root‐mean‐square error occurs in the 6‐hr forecast field as compared to the European Centre for Medium‐Range Weather Forecasts (ECMWF)'s Integrated Forecasting System (IFS) analysis. Even though the all‐sky MHS assimilation's impact on the Northern Hemisphere midlatitudes is not remarkable, it nonetheless produces an improved humidity analysis increment for a heavy rainfall case over East Asia.

Assimilation of Himawari-8 Rapid-Scan Atmospheric Motion Vectors on Tropical Cyclone in HWRF System

This study investigates the assimilation impact of rapid-scan (RS) atmospheric motion vectors (AMVs) derived from the geostationary satellite Himawari-8 on tropical cyclone (TC) forecasts. Forecast experiments for three TCs in 2016 in the western North Pacific basin are performed using the National Centers for Environmental Prediction (NCEP) operational Hurricane Weather Research and Forecasting Model (HWRF). An ensemble-variational hybrid data assimilation system is used as an initialization. The results show that the assimilation of RS-AMVs can improve the track forecast skill, while the weak bias or slow intensification bias increases at the shorter forecast lead time. A vortex initialization in HWRF has a substantial impact on TC structure, but it has neutral impacts on the track and intensity forecasts. A thinning of AMVs mitigates the weak bias caused by RS-AMV assimilation, resulting in the reduction of intensity error. However, it degrades the track forecast skill for a longer lead time. A decomposition of the TC steering flows demonstrated that the change in TC-induced flow was a primary factor for reducing the track forecast error, and the change in environmental flow has less impact on the track forecast. The investigation of the structural change from the assimilation of RS-AMV revealed that the following two factors are likely related to the intensity forecast degradation: (1) an increase of inertial stability outside the radius of maximum wind (RMW), which weakens the boundary layer inflow; and (2) a drying around and outside the RMW.

Impact of Model Bias Correction on a Hybrid Data Assimilation System

Hybrid data assimilation combines a conventional 3-D or 4-D variational system with background error covariance (BEC) generated from ensemble forecast systems. In order to achieve better BEC, three perturbation schemes, namely, the random combination of multiple physical paramterization schemes (referred to as MP), the MP plus stochastical perturbation on physical process tendencies (MP-SPPT), and the unified perturbation of stochastic physics with bias correction (UPSB, proposed by the authors of this paper in a previous work), were first used in a regional ensemble model, i.e., the Global and Regional Assimilation and Prediction System-Regional Ensemble Prediction System (GRAPES-REPS), and the BECs thus obtained were compared for 7-day ensemble forecasts. The results show that UPSB, which is in fact an MP-SPPT but with the systematic model bias removed, has a better consistency, i.e., the ratio between root-mean-square error (RMSE) and ensemble spread is much closer to 1, especially at low model levels, compared to the other two schemes. Moreover, the BEC derived from UPSB captured more reasonable distributions of forecast errors. Second, performance of a hybrid data assimilation system (the GRAPES-MESO hybrid En-3DVar) was evaluated by using the BECs from the three perturbation schemes for 7-day hybrid data assimilation forecasts, and thus disclosing the effect of the model bias correction (assuming that the random stocastical features are in general offset in the three perturbation schemes) on the hybrid system forecasts. A covariance weight of 0.8 was prescribed, and this value was determined through sensitivity experiments. The forecast results from the hybrid data assimilation system show that UPSB reduced the false correlation between distant points. The quality of analysis fields of the UPSB scheme shows visible improvement, i.e., the analysis fields produced by UPSB have much smaller RMSEs than those of the other two schemes, at all vertical model levels. The quality of the hybrid data assimilation forecast fields was also improved by this scheme. Furthermore, the improvement was much greater in the early stage of the assimilation cycle than in the late stage. Generally, the quality of the hybrid data assimilation of GRAPES-MESO hybrid En-3DVar could be efficiently improved by the model bias correction in the UPSB scheme.

Improving Hurricane Analyses and Predictions with TCI, IFEX Field Campaign Observations, and CIMSS AMVs Using the Advanced Hybrid Data Assimilation System for HWRF. Part II: Observation Impacts on the Analysis and Prediction of Patricia (2015)

Abstract Diverse observations, such as the High Definition Sounding System (HDSS) dropsonde observations from the Tropical Cyclone Intensity (TCI) program, the Tail Doppler Radar (TDR), Stepped Frequency Microwave Radiometer (SFMR), and flight-level observations from the Intensity Forecasting Experiment (IFEX) program, and the atmospheric motion vectors (AMVs) from the Cooperative Institute for Meteorological Satellite Studies (CIMSS) simultaneously depicted the three-dimensional (3D) structure of Hurricane Patricia (2015). Experiments are conducted to understand the relative impacts of each of these observation types on Patricia’s analysis and prediction using the Gridpoint Statistical Interpolation (GSI)-based ensemble-variational data assimilation system for the Hurricane Weather Research and Forecasting (HWRF) Model. In comparing the impacts of assimilating each dataset individually, results suggest that 1) the assimilation of 3D observations produces better TC structure analysis than the assimilation of two-dimensional (2D) observations; 2) the analysis from assimilating observations collected from platforms that only sample momentum fields produces a less improved forecast with either short-lived impacts or slower intensity spinup as compared to the forecast produced after assimilating observations collected from platforms that sample both momentum and thermal fields; and 3) the structure forecast tends to benefit more from the assimilation of inner-core observations than the corresponding intensity forecast, which implies better verification metrics are needed for future TC forecast evaluation.

Improving Hurricane Analyses and Predictions with TCI, IFEX Field Campaign Observations, and CIMSS AMVs Using the Advanced Hybrid Data Assimilation System for HWRF. Part I: What is Missing to Capture the Rapid Intensification of Hurricane Patricia (2015) when HWRF is already Initialized with a More Realistic Analysis?

Abstract Assimilating inner-core observations collected from recent field campaign programs such as Tropical Cyclone Intensity (TCI) and Intensity Forecasting Experiment (IFEX) together with the enhanced atmospheric motion vectors (AMVs) produce realistic three-dimensional (3D) analyses using the newly developed GSI-based, continuously cycled, dual-resolution hybrid ensemble–variational data assimilation (DA) system for the Hurricane Weather Research and Forecasting (HWRF) Model for Hurricane Patricia (2015). However, more persistent surface wind maximum spindown is found in the intensity forecast initialized from the realistic analyses produced by the DA system but not from the unrealistic initial conditions produced through vortex modification. Diagnostics in this study reveal that the spindown issue is likely attributed to the deficient HWRF Model physics that are unable to maintain the realistic 3D structures from the DA analysis. The horizontal diffusion is too strong to maintain the realistically observed vertical oscillation of radial wind near the eyewall region. The vertical diffusion profile cannot produce a sufficiently strong secondary circulation connecting the realistically elevated upper-level outflow produced in the DA analysis. Further investigations with different model physics parameterizations demonstrate that spindown can be alleviated by modifying model physics parameterizations. In particular, a modified turbulent mixing parameterization scheme together with a reduced horizontal diffusion is found to significantly alleviate the spindown issue and to improve the intensity forecast. Additional experiments show that the peak-simulated intensity and rapid intensification rate can be further improved by increasing the model resolution. But the model resolution is not as important as model physics in the spindown alleviation.

Development of an Operational Hybrid Data Assimilation System at KIAPS

This study introduces the operational data assimilation (DA) system at the Korea Institute of Atmospheric Prediction Systems (KIAPS) to the numerical weather prediction community. Its development history and performance are addressed with experimental illustrations and the authors’ previously published studies. Milestones in skill improvements include the initial operational implementation of three-dimensional variational data assimilation (3DVar), the ingestion of additional satellite observations, and changing the DA scheme to a hybrid four-dimensional ensemble-variational DA using forecasts from an ensemble based on the local ensemble transform Kalman filter (LETKF). In the hybrid system, determining the relative contribution of the ensemble-based covariance to the resultant analysis is crucial, particularly for moisture variables including a variety of horizontal scale spectra. Modifications to the humidity control variable, partial rather than full recentering of the ensemble for humidity further improves moisture analysis, and the inclusion of more radiance observations with higher-level peaking channels have significant impacts on stratosphere temperature and wind performance. Recent update of the operational hybrid DA system relative to the previous 3DVar system is described for detailed improvements with interpretation.

A Prototype Regional GSI-based EnKF-Variational Hybrid Data Assimilation System for the Rapid Refresh Forecasting System: Dual-Resolution Implementation and Testing Results

A dual-resolution (DR) version of a regional ensemble Kalman filter (EnKF)-3D ensemble variational (3DEnVar) coupled hybrid data assimilation system is implemented as a prototype for the operational Rapid Refresh forecasting system. The DR 3DEnVar system combines a high-resolution (HR) deterministic background forecast with lower-resolution (LR) EnKF ensemble perturbations used for flow-dependent background error covariance to produce a HR analysis. The computational cost is substantially reduced by running the ensemble forecasts and EnKF analyses at LR. The DR 3DEnVar system is tested with 3-h cycles over a 9-day period using a 40/∼13-km grid spacing combination. The HR forecasts from the DR hybrid analyses are compared with forecasts launched from HR Gridpoint Statistical Interpolation (GSI) 3D variational (3DVar) analyses, and single LR hybrid analyses interpolated to the HR grid. With the DR 3DEnVar system, a 90% weight for the ensemble covariance yields the lowest forecast errors and the DR hybrid system clearly outperforms the HR GSI 3DVar. Humidity and wind forecasts are also better than those launched from interpolated LR hybrid analyses, but the temperature forecasts are slightly worse. The humidity forecasts are improved most. For precipitation forecasts, the DR 3DEnVar always outperforms HR GSI 3DVar. It also outperforms the LR 3DEnVar, except for the initial forecast period and lower thresholds.

Impact of Representing Model Error in a Hybrid Ensemble-Variational Data Assimilation System for Track Forecast of Tropical Cyclones over the Bay of Bengal

Uncertainties in the numerical weather prediction models are generally not well-represented in ensemble-based data assimilation (DA) systems. The performance of an ensemble-based DA system becomes suboptimal, if the sources of error are undersampled in the forecast system. The present study examines the effect of accounting for model error treatments in the hybrid ensemble transform Kalman filter—three-dimensional variational (3DVAR) DA system (hybrid) in the track forecast of two tropical cyclones viz. Hudhud and Thane, formed over the Bay of Bengal, using Advanced Research Weather Research and Forecasting (ARW-WRF) model. We investigated the effect of two types of model error treatment schemes and their combination on the hybrid DA system; (i) multiphysics approach, which uses different combination of cumulus, microphysics and planetary boundary layer schemes, (ii) stochastic kinetic energy backscatter (SKEB) scheme, which perturbs the horizontal wind and potential temperature tendencies, (iii) a combination of both multiphysics and SKEB scheme. Substantial improvements are noticed in the track positions of both the cyclones, when flow-dependent ensemble covariance is used in 3DVAR framework. Explicit model error representation is found to be beneficial in treating the underdispersive ensembles. Among the model error schemes used in this study, a combination of multiphysics and SKEB schemes has outperformed the other two schemes with improved track forecast for both the tropical cyclones.

Influence of the Self-Consistent Regional Ensemble Background Error Covariance on Hurricane Inner-Core Data Assimilation with the GSI-Based Hybrid System for HWRF

Abstract An initial vortex spindown, or strong adjustment to the structure and intensity of a hurricane’s initial vortex, presents a significant problem in hurricane forecasting, as with the NCEP Hurricane Weather Research and Forecasting Model (HWRF), because it can cause significantly degraded intensity forecasts. In this study, the influence of the self-consistent regional ensemble background error covariance on assimilating hurricane inner-core tail Doppler radar (TDR) observations in HWRF is examined with the NCEP gridpoint statistical interpolation (GSI)-based ensemble–three-dimensional variational (3DVAR) hybrid data assimilation system. It is found that the resolution of the background error covariance term, coming from the ensemble forecasts, has notable influence on the assimilation of hurricane inner-core observations and subsequent forecasting results. Specifically, the use of ensemble forecasting at high-resolution native grids results in significant reduction of the vortex spindown problem and thus leads to improved hurricane intensity forecasting. Further diagnoses are conducted to examine the spindown problem with a gradient wind balance. It is found that artificial vortex initialization, performed before data assimilation, can cause strong supergradient winds or imbalance in the vortex inner-core region. Assimilation of hurricane inner-core TDR data can significantly mitigate this imbalance by reducing the supergradient effects. Compared with the use of a global ensemble background error term, application of the self-consistent regional ensemble background covariance to inner-core data assimilation leads to better representation of the mesoscale hurricane inner-core structures. It can also result in more realistic vortex structures in data assimilation even when the observational data are unevenly distributed.

Mesoscale Hybrid Data Assimilation System based on JMA Nonhydrostatic Model

This paper discusses the benefits of using a hybrid ensemble Kalman filter and four-dimensional variational (4D-Var) data assimilation (DA) system rather than a 4D-Var system employing the National Meteorological Center (NMC, now known as NCEP) method (4D-Var-Bnmc) to predict severe weather events. An adjoint-based 4D-Var system was employed with a background error covariance matrix [Formula: see text] constructed from the NMC method and perturbations in a local ensemble transform Kalman filter system. The DA systems are based on the Japan Meteorological Agency’s nonhydrostatic model. To reduce the sampling noise, three types of implementation (the spatial localization, spectral localization, and neighboring ensemble approaches) were tested. The assimilation of a pseudosingle observation of sea level pressure located at a tropical cyclone (TC) center yielded analysis increments physically consistent with what is expected of a mature TC in the hybrid systems at the beginning of the assimilation window, whereas analogous experiments performed using the 4D-Var-Bnmc system did not. At the end, the structures of the 4D-Var-based increments became similar to one another, while the analysis increment by the 4D-Var-Bnmc system was broad in the horizontal direction. Realistic DA experiments showed that all of the hybrid systems provided initial conditions that yielded more accurate TC track and intensity forecasts than those achievable by the 4D-Var-Bnmc system. The hybrid systems also yielded some statistically significant improvements in forecasting local heavy rainfall events in terms of fraction skill scores when a 160 km × 160 km window size was used. The overall skills of the hybrid systems were relatively independent of the choice of implementation.