Xiang-Yu Huang

Four-Dimensional Variational Data Assimilation for WRF: Formulation and Preliminary Results

Abstract The Weather Research and Forecasting (WRF) model–based variational data assimilation system (WRF-Var) has been extended from three- to four-dimensional variational data assimilation (WRF 4D-Var) to meet the increasing demand for improving initial model states in multiscale numerical simulations and forecasts. The initial goals of this development include operational applications and support to the research community. The formulation of WRF 4D-Var is described in this paper. WRF 4D-Var uses the WRF model as a constraint to impose a dynamic balance on the assimilation. It is shown to implicitly evolve the background error covariance and to produce the flow-dependent nature of the analysis increments. Preliminary results from real-data 4D-Var experiments in a quasi-operational setting are presented and the potential of WRF 4D-Var in research and operational applications are demonstrated. A wider distribution of the system to the research community will further develop its capabilities and to encoura...

The Weather Research and Forecasting Model's Community Variational/Ensemble Data Assimilation System: WRFDA

Data assimilation is the process by which observations are combined with short-range NWP model output to produce an analysis of the state of the atmosphere at a specified time. Since its inception in the late 1990s, the multiagency Weather Research and Forecasting (WRF) model effort has had a strong data assimilation component, dedicating two working groups to the subject. This article documents the history of the WRF data assimilation effort, and discusses the challenges associated with balancing academic, research, and operational data assimilation requirements in the context of the WRF effort to date. The WRF Models Community Variational/Ensemble Data Assimilation System (WRFDA) has evolved over the past 10 years, and has resulted in over 30 refereed publications to date, as well as implementation in a wide range of real-time and operational NWP systems. This paper provides an overview of the scientific capabilities of WRFDA, and together with results from sample operation implementations at the U.S. ...

Radar Data Assimilation with WRF 4D-Var. Part I: System Development and Preliminary Testing

AbstractThe major goal of this two-part study is to assimilate radar data into the high-resolution Advanced Research Weather Research and Forecasting Model (ARW-WRF) for the improvement of short-term quantitative precipitation forecasting (QPF) using a four-dimensional variational data assimilation (4D-Var) technique. In Part I the development of a radar data assimilation scheme within the WRF 4D-Var system (WRF 4D-Var) and the preliminary testing of the scheme are described. In Part II the performance of the enhanced WRF 4D-Var system is examined by comparing it with the three-dimensional variational data assimilation system (WRF 3D-Var) for a convective system over the U.S. Great Plains. The WRF 4D-Var radar data assimilation system has been developed with the existing framework of an incremental formulation. The new development for radar data assimilation includes the tangent-linear and adjoint models of a Kessler warm-rain microphysics scheme and the new control variables of cloud water, rainwater, an...

Intercomparison of an Ensemble Kalman Filter with Three- and Four-Dimensional Variational Data Assimilation Methods in a Limited-Area Model over the Month of June 2003

Abstract This study compares the performance of an ensemble Kalman filter (EnKF) with both the three-dimensional and four-dimensional variational data assimilation (3DVar and 4DVar) methods of the Weather Research and Forecasting (WRF) model over the contiguous United States in a warm-season month (June) of 2003. The data assimilated every 6 h include conventional sounding and surface observations as well as data from wind profilers, ships and aircraft, and the cloud-tracked winds from satellites. The performances of these methods are evaluated through verifying the 12- to 72-h forecasts initialized twice daily from the analysis of each method against the standard sounding observations. It is found that 4DVar has consistently smaller error than that of 3DVar for winds and temperature at all forecast lead times except at 60 and 72 h when their forecast errors become comparable in amplitude, while the two schemes have similar performance in moisture at all lead times. The forecast error of the EnKF is compa...

Indirect Assimilation of Radar Reflectivity with WRF 3D-Var and Its Impact on Prediction of Four Summertime Convective Events

AbstractAn indirect radar reflectivity assimilation scheme has been developed within the Weather Research and Forecasting model three-dimensional data assimilation system (WRF 3D-Var). This scheme, instead of assimilating radar reflectivity directly, assimilates retrieved rainwater and estimated in-cloud water vapor. An analysis is provided to show that the assimilation of the retrieved rainwater avoids the linearization error of the Z–qr (reflectivity–rainwater) equation. A new observation operator is introduced to assimilate the estimated in-cloud water vapor. The performance of the scheme is demonstrated by assimilating reflectivity observations into the Rapid Update Cycle data assimilation and forecast system operating at Beijing Meteorology Bureau. Four heavy-rain-producing convective cases that occurred during summer 2009 in Beijing, China, are studied using the newly developed system. Results show that on average the assimilation of reflectivity significantly improves the short-term precipitation f...

Impact of assimilating microwave radiances with a limited-area ensemble data assimilation system on forecasts of Typhoon Morakot

AbstractTwo parallel experiments were designed to evaluate whether assimilating microwave radiances with a cyclic, limited-area ensemble adjustment Kalman filter (EAKF) could improve track, intensity, and precipitation forecasts of Typhoon Morakot (2009). The experiments were configured identically, except that one assimilated microwave radiances and the other did not. Both experiments produced EAKF analyses every 6 h between 1800 UTC 3 August and 1200 UTC 9 August 2009, and the mean analyses initialized 72-h Weather Research and Forecasting model forecasts. Examination of individual forecasts and average error statistics revealed that assimilating microwave radiances ultimately resulted in better intensity forecasts compared to when radiances were withheld. However, radiance assimilation did not substantially impact track forecasts, and the impact on precipitation forecasts was mixed. Overall, net positive results suggest that assimilating microwave radiances with a limited-area EAKF system is beneficial...

Comparing Limited-Area 3DVAR and Hybrid Variational-Ensemble Data Assimilation Methods for Typhoon Track Forecasts: Sensitivity to Outer Loops and Vortex Relocation

AbstractThe 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 WRFs 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 dur...

Application of an Adiabatic WRF Adjoint to the Investigation of the May 2004 McMurdo, Antarctica, Severe Wind Event

Abstract The tangent linear and adjoint of an adiabatic version of the Weather Research and Forecasting (WRF) Model with its Advanced Research WRF (ARW) dynamic core have been developed. The source-to-source automatic differentiation tool [i.e., the Transformation of Algorithm (TAF) in FORTRAN] was used in the development. Tangent linear and adjoint checks of the developed adiabatic WRF adjoint modeling system (WAMS) were conducted, and all necessary correctness verification procedures were passed. As the first application, the adiabatic WAMS was used to study the adjoint sensitivity of a severe windstorm in Antarctica. Linearity tests indicated that an adjoint-based sensitivity study with the Antarctic Mesoscale Prediction System (AMPS) 90-km domain configuration for the windstorm is valid up to 24 h. The adjoint-based sensitivity calculation with adiabatic WAMS identified sensitive regions for the improvement of the 24-h forecast of the windstorm. It is indicated that the windstorm forecast largely reli...

Development of the upgraded tangent linear and adjoint of the weather research and forecasting (WRF) Model

The authors propose a new technique for parallelizations of tangent linear and adjoint codes, which were applied in the redevelopment for the Weather Research and Forecasting (WRF) model with its Advanced ResearchWRFdynamiccoreusingtheautomaticdifferentiation engine.Thetangentlinearandadjointcodes of the WRF model (WRFPLUS) now have the following improvements: A complete check interface ensures that developers write accurate tangent linear and adjoint codes with ease and efficiency. A new technique based on the nature of duality that existed among message passing interface communication routines was adopted to parallelize the WRFPLUS model. The registry in the WRF model was extended to automatically generate the tangent linear and adjoint codes of the required communication operations. This approach dramatically speeds up the software development cycle of the parallel tangent linear and adjoint codes and leads to improved parallel efficiency. Module interfaces were constructed for coupling tangent linear and adjoint codes of the WRF model with applications such as four-dimensional variational data assimilation, forecast sensitivity to observation, and others.

Bridging Research to Operations Transitions: Status and Plans of Community GSI

AbstractWith a goal of improving operational numerical weather prediction (NWP), the Developmental Testbed Center (DTC) has been working with operational centers, including, among others, the National Centers for Environmental Prediction (NCEP), National Oceanic and Atmospheric Administration (NOAA), National Aeronautics and Space Administration (NASA), and the U.S. Air Force, to support numerical models/systems and their research, perform objective testing and evaluation of NWP methods, and facilitate research-to-operations transitions. This article introduces the first attempt of the DTC in the data assimilation area to help achieve this goal. Since 2009, the DTC, NCEP’s Environmental Modeling Center (EMC), and other developers have made significant progress in transitioning the operational Gridpoint Statistical Interpolation (GSI) data assimilation system into a community-based code management framework. Currently, GSI is provided to the public with user support and is open for contributions from inter...