Four-dimensional Variational Data Assimilation System Research Articles

Impact of GNSS radio occultation bending angle data assimilation in YH4DVAR system

Global Navigation Satellite System (GNSS) radio occultation measurements have been assimilated into the four-dimensional variational data assimilation system (YH4DVAR) using a one-dimensional bending angle operator (GBAO) as a new type of observation. For the sake of verifying the impact of GNSS radio occultation (RO) measurements to the data assimilation system, three experiments have been conducted. The statistical results of the analysis error experiment and forecast skill experiment show that the GNSS RO measurements have an impact on the analysis system. The typhoon forecast experiment shows the impact on the important weather process. They all have a positive impact on the weather forecast. Lastly, we look forward to future work on the observation system simulation experiment (OSSE) to investigate the impact of GNSS RO measurements as a function of observation number, which is an effective method to estimate the saturation of the observation number.

Inferring regional sources and sinks of atmospheric CO<sub>2</sub> from GOSAT XCO<sub>2</sub> data

Abstract. We have examined the utility of retrieved column-averaged, dry-air mole fractions of CO2 (XCO2) from the Greenhouse Gases Observing Satellite (GOSAT) for quantifying monthly, regional flux estimates of CO2, using the GEOS-Chem four-dimensional variational (4D-Var) data assimilation system. We focused on assessing the potential impact of biases in the GOSAT CO2 data on the regional flux estimates. Using different screening and bias correction approaches, we selected three different subsets of the GOSAT XCO2 data for the 4D-Var inversion analyses, and found that the inferred global fluxes were consistent across the three XCO2 inversions. However, the GOSAT observational coverage was a challenge for the regional flux estimates. In the northern extratropics, the inversions were more sensitive to North American fluxes than to European and Asian fluxes due to the lack of observations over Eurasia in winter and over eastern and southern Asia in summer. The regional flux estimates were also sensitive to the treatment of the residual bias in the GOSAT XCO2 data. The largest differences obtained were for temperate North America and temperate South America, for which the largest spread between the inversions was 1.02 and 0.96 Pg C, respectively. In the case of temperate North America, one inversion suggested a strong source, whereas the second and third XCO2 inversions produced a weak and strong sink, respectively. Despite the discrepancies in the regional flux estimates between the three XCO2 inversions, the a posteriori CO2 distributions were in good agreement (with a mean difference between the three inversions of typically less than 0.5 ppm) with independent data from the Total Carbon Column Observing Network (TCCON), the surface flask network, and from the HIAPER Pole-to-Pole Observations (HIPPO) aircraft campaign. The discrepancy in the regional flux estimates from the different inversions, despite the agreement of the global flux estimates suggests the need for additional work to determine the minimum spatial scales at which we can reliably quantify the fluxes using GOSAT XCO2. The fact that the a posteriori CO2 from the different inversions were in good agreement with the independent data although the regional flux estimates differed significantly, suggests that innovative ways of exploiting existing data sets, and possibly additional observations, are needed to better evaluate the inferred regional flux estimates.

Wind extraction potential from 4D-Var assimilation of stratospheric O<sub>3</sub>, N<sub>2</sub>O, and H<sub>2</sub>O using a global shallow water model

Abstract. The wind extraction due to assimilation of stratospheric trace gas (tracer) data is examined using a 4D-Var (four-dimensional variational) data assimilation system based on the shallow water equations coupled to the tracer continuity equation. The procedure is outlined as follows. First, a nature run is created, simulating middle stratospheric winter conditions. Second, ozone (O3), nitrous oxide (N2O), and water vapor (H2O) (treated in this study as passive tracers) are initialized using Aura Microwave Limb Sounder (MLS) mixing ratios at 850 K potential temperature and are advected by the nature run winds. Third, the initial dynamical conditions are perturbed by using a 6 h offset. Fourth, simulated hourly tracer observations on the full model grid are assimilated with a 4D-Var system in which tracer and winds are coupled via the adjoint of the tracer continuity equation. Multiple assimilation experiments are performed by varying the amount of random observation error added to the simulated measurements. Finally, the wind extraction potential (WEP) is calculated as the reduction of the vector wind root mean square error (RMSE) relative to the maximum possible reduction. For a single 6 h assimilation cycle with the smallest observation error, WEP values are ~60% for all three tracers, while 10-day multi-cycle simulations result in WEP of ~90%, wind errors of ~0.3 m s−1, and height errors of ~13 m. There is therefore sufficient information in the tracer fields to almost completely constrain the dynamics, even without direct assimilation of dynamical information. When realistic observation error is added (based on MLS precisions at 10 hPa), the WEP after 10 days is 90% for O3, 87% for N2O, and 72% for H2O. O3 and N2O provide more wind information than H2O due to stronger background gradients relative to the MLS precisions. The RMSE for wind reach a minimum level of ~0.3–0.9 m s−1 for the MLS precisions, suggesting a limit to which realistic tracers could constrain the winds, given complete global cover age. With higher observation noise levels, the WEP values decrease, but the impact on the winds is still positive up to noise levels of 100% (relative to the global mean value) when compared to the case of no data assimilation.

The influence of erroneous background, beam-blocking and microphysical non-linearity on the application of a four-dimensional variational Doppler radar data assimilation system for quantitative precipitation forecasts

A series of observation system simulation experiments (OSSEs) and a real case study are conducted to investigate the application of the Doppler radar data assimilation technique for numerical model quantitative precipitation forecasts (QPFs). A four-dimensional variational Doppler radar analysis system (VDRAS) is adopted for all experiments. The first set of OSSEs demonstrates that when the background field contains the imperfect information predicted from a mesoscale model, the incorrect convective-scale perturbations in the background can result in spurious scattered precipitation. However, a smoothing procedure can be used to remove the fine structures from the primitive model output in order to avoid this over-prediction. Results from the second set of OSSEs indicate that the lack of low-elevation data owing to radar scan and/or beam blockage could significantly alter the retrieved low-level thermal and dynamical structures when a different number of data assimilation cycles is applied. These impacts could lower the rainfall forecast capability of the model. The third set of OSSEs shows that, when the rainwater is assimilated over a long assimilation window, the non-linearity embedded in the microphysical process could lead the minimization algorithm in a wrong direction, causing a further degradation of the rainfall prediction. However, using multiple short assimilation cycles produces better minimization and forecast results than those obtained with a single long cycle. A real case experiment based on data collected during Intensive Operation Period (IOP) #8 of the 2008 Southwest Monsoon Experiment (SoWMEX) is conducted to provide a verification of the conclusions obtained from OSSEs under a realistic framework.

Satellite Radiance Assimilation in the JMA Operational Mesoscale 4DVAR System

Abstract The direct radiance assimilation scheme used in the Japan Meteorological Agency (JMA) global analysis system is applied to the JMA mesoscale four-dimensional variational data assimilation (4DVAR) system with two modifications. First, the data-thinning distance is shortened, and, second, the atmospheric profiles are extrapolated from the mesoscale model top to the radiative transfer model top using the U.S. Standard Atmosphere lapse rate. Although the variational bias correction method is widely used in many numerical weather prediction centers for global radiance assimilations, a radiance bias correction method for regional models has not been established because of difficulties in estimating the biases within limited regions and times. This paper examined the use of the bias correction coefficients estimated in the global system for the mesoscale system when the radiance data were introduced instead of the retrievals. It was found that the profile extrapolation was necessary to reduce the biases. Moreover, the use of common bias coefficients enables the use of the radiance data in the same way as the global system. The radiance data assimilation experiments in the mesoscale system demonstrated considerable improvements to the tropospheric geopotential height forecasts and precipitation forecasts. The improvements resulted from the introduction of radiance data from multiple satellites into data-sparse regions and times. However, the major effect of the radiance assimilation on the precipitation forecasts was limited to weak precipitation areas over oceans; the effects on deep convective areas and over land were relatively small.

Impacts of 4DVAR Assimilation of Airborne Doppler Radar Observations on Numerical Simulations of the Genesis of Typhoon Nuri (2008)

AbstractThe Weather Research and Forecasting Model and its four-dimensional variational data assimilation (4DVAR) system are employed to examine the impact of airborne Doppler radar observations on predicting the genesis of Typhoon Nuri (2008). Electra Doppler Radar (ELDORA) airborne radar data, collected during the Office of Naval Research–sponsored Tropical Cyclone Structure 2008 field experiment, are used for data assimilation experiments. Two assimilation methods are evaluated and compared, namely, the direct assimilation of radar-measured radial velocity and the assimilation of three-dimensional wind analysis derived from the radar radial velocity. Results show that direct assimilation of radar radial velocity leads to better intensity forecasts, as this process enhances the development of convective systems and improves the inner-core structure of Nuri, whereas assimilation of the radar-retrieved wind analysis is more beneficial for tracking forecasts, as it results in improved environmental flows. The assimilation of both the radar-retrieved wind and the radial velocity can lead to better forecasts in both intensity and tracking, if the radial velocity observations are assimilated first and the retrieved winds are then assimilated in the same data assimilation window. In addition, experiments with and without radar data assimilation led to developing and nondeveloping disturbances in numerical simulations of Nuri’s genesis. The improved initial conditions and forecasts from the data assimilation imply that the enhanced midlevel vortex and moisture conditions are favorable for the development of deep convection in the center of the pouch and eventually contribute to Nuri’s genesis. The improved simulations of the convection and associated environmental conditions produce enhanced upper-level warming in the core region and lead to the drop in sea level pressure.

Radar Data Assimilation with WRF 4D-Var. Part II: Comparison with 3D-Var for a Squall Line over the U.S. Great Plains

Abstract The Weather Research and Forecasting Model (WRF) four-dimensional variational data assimilation (4D-Var) system described in Part I of this study is compared with its corresponding three-dimensional variational data assimilation (3D-Var) system using a Great Plains squall line observed during the International H2O Project. Two 3D-Var schemes are used in the comparison: a standard 3D-Var radar data assimilation (DA) that is the same as the 4D-Var except for the exclusion of the constraining dynamical model and an enhanced 3D-Var that includes a scheme to assimilate an estimated in-cloud humidity field. The comparison is made by verifying their skills in 0–6-h quantitative precipitation forecast (QPF) against stage-IV analysis, as well as in wind forecasts against radial velocity observations. The relative impacts of assimilating radial velocity and reflectivity on QPF are also compared between the 4D-Var and 3D-Var by conducting data-denial experiments. The results indicate that 4D-Var substantially improves the QPF skill over the standard 3D-Var for the entire 6-h forecast range and over the enhanced 3D-Var for most forecast hours. Radial velocity has a larger impact relative to reflectivity in 4D-Var than in 3D-Var in the first 3 h because of a quicker precipitation spinup. The analyses and forecasts from the 4D-Var and 3D-Var schemes are further compared by examining the meridional wind, horizontal convergence, low-level cold pool, and midlevel temperature perturbation, using analyses from the Variational Doppler Radar Analysis System (VDRAS) as references. The diagnoses of these fields suggest that the 4D-Var analyzes the low-level cold pool, its leading edge convergence, and midlevel latent heating in closer resemblance to the VDRAS analyses than the 3D-Var schemes.

Experiments of the WRF three-/four-dimensional variational (3/4DVAR) data assimilation in the forecasting of Antarctic cyclones

The three-/four-dimensional variational data assimilation systems (3/4DVAR) of the Weather Research and Forecasting (WRF) model were explored in the forecasting of two Antarctic synoptic cyclones, which had large influence on the Ross Sea/Ross Ice Shelf region in October 2007. A suite of variational data assimilation experiments, including regular 3DVAR, high-resolution 3DVAR, and 4DVAR experiments, were designed to evaluate their performances in weather analysis and forecasting in Antarctica. In general, both 4DVAR and high-resolution 3DVAR experiments showed better forecasting skill than regular 3DVAR experiments. High-resolution 3DVAR experiments were the most efficient in reducing the analysis errors of surface winds and temperature, and had the best performance during the first 24 h of forecasting. However, during the following forecast period, 4DVAR experiments showed either better or about comparable performance to high-resolution 3DVAR experiments. These results indicate that increasing the spatial resolution during 3DVAR is an economical approach to improving the weather analysis and forecasting over Antarctica. At the same time, the 4DVAR approach had a longer impact on forecasting than the high-resolution 3DVAR approach. Understandably, both of the variational assimilation approaches are promising techniques toward improving the regional analysis and forecasting over Antarctica.

Research on the implementation of GPS radio occultation one-dimensional bending angle operator in four-dimensional variational data assimilation system

GPS radio occultation data are implemented in the numerical weather prediction as a kind of very important observation data because of its excellent characteristics, which makes a remarkable positive effect in the system. In this article, we design an implement method and a quality control strategy of one-dimensional bending angle operator in four-dimensional variational data assimilation system, and we also describe the steps of implementing one-dimensional bending angle forward operator and its linear and adjoint operators in four-dimensional variational data assimilation system. The author chose Reger’s three term expression and its best average coefficients are chosen to calculate refractivity values. In this article, we introduce an observation error model of one-dimensional bending angle operator through the statistics of GPS radio occultation observation data of August, 2012. A comparison of experiment of the analysis field between in the cases with and without GPS radio occultation data is made. It is concluded that GPS radio occultation data plays a positive role in the global areas, especially in the low and middle latitude areas higher than 10 km. This conclusion is the foundation of the following-up GPS radio occultation data assimilation verification experiments.

Estimates of Analysis and Forecast Error Variances Derived from the Adjoint of 4D-Var

Abstract A method is presented in which the adjoint of a four-dimensional variational data assimilation system (4D-Var) was used to compute the expected analysis and forecast error variances of linear functions of the ocean state vector. The power and utility of the approach are demonstrated using the Regional Ocean Modeling System configured for the California Current system. Linear functions of the ocean state vector were considered in the form of indices that characterize various aspects of the coastal upwelling circulation. It was found that for configurations of 4D-Var typically used in ocean models, reliable estimates of the expected analysis error variances can be obtained both for variables that are observed and unobserved. In addition, the contribution of uncertainties in the model control variables to the forecast error variance was also quantified. One particularly powerful and illuminating aspect of the adjoint 4D-Var approach to the forecast problem is that the contribution of individual observations to the predictability of the circulation can be readily computed. An important finding of the work presented here is that despite the plethora of available satellite observations, the relatively modest fraction of in situ subsurface observations sometimes exerts a significant influence on the predictability of the coastal ocean. Independent checks of the analysis and forecast error variances are also presented, which provide a direct test of the hypotheses that underpin the prior error and observation error estimates used during 4D-Var.

Accounting for representativeness errors in the inversion of atmospheric constituent emissions: application to the retrieval of regional carbon monoxide fluxes

ABSTRACTA four-dimensional variational data assimilation system (4D-Var) is developed to retrieve carbon monoxide (CO) fluxes at regional scale, using an air quality network. The air quality stations that monitor CO are proximity stations located close to industrial, urban or traffic sources. The mismatch between the coarsely discretised Eulerian transport model and the observations, inferred to be mainly due to representativeness errors in this context, lead to a bias (average simulated concentrations minus observed concentrations) of the same order of magnitude as the concentrations. 4D-Var leads to a mild improvement in the bias because it does not adequately handle the representativeness issue. For this reason, a simple statistical subgrid model is introduced and is coupled to 4D-Var. In addition to CO fluxes, the optimisation seeks to jointly retrieve influence coefficients, which quantify each station's representativeness. The method leads to a much better representation of the CO concentration variability, with a significant improvement of statistical indicators. The resulting increase in the total inventory estimate is close to the one obtained from remote sensing data assimilation. This methodology and experiments suggest that information useful at coarse scales can be better extracted from atmospheric constituent observations strongly impacted by representativeness errors.

Design and implementation of four-dimensional variational data assimilation system constrained by the global spectral model

The mainstream data assimilation system in operation will still employ four-dimensional variational data assimilation(4DVAR) method in a long time of future. We develop a new 4DVAR system, i.e., YH4DVAR, using global spectral model as a constraint to impose a dynamic balance on the assimilation. The cost function of YH4DVAR consists of four terms: background, observations, digital filter, and bias correction, respectively. YH4DVAR employs the wavelet background error convariance, the multi-resolution incremental formulation, the tangent linear and adjoint models for dynamics core and physical processes, and ATOVS radiance data assimilation. Simultaneous spatial and spectral variations of horizontal and vertical covariance are achieved by dividing the control vector into several parts, each of which corresponds to a band of total spherical wavenumbers. Using NWP consisting of YH4DVAR and global spectral model, the anomaly correlation with the verifying analysis for geopotential height 8-day forecast on the 500 hPa isobaric surface at 12-month mean is above 60%. The formulation and implementation of YH4VAR are described in detail in this paper.

Hindcasting the January 2009 Arctic Sudden Stratospheric Warming with Unified Parameterization of Orographic Drag in NOGAPS. Part II: Short-Range Data-Assimilated Forecast and the Impacts of Calibrated Radiance Bias Correction

Abstract This study is Part II of the effort to improve the forecasting of sudden stratospheric warming (SSW) events by using a version of the Navy Operational Global Atmospheric Prediction System (NOGAPS) that covers the full stratosphere. In Part I, extended-range (3 week) hindcast experiments (without data assimilation) for the January 2009 Arctic major SSW were performed using NOGAPS with a unified orographic drag parameterization that consists of the schemes employed by Webster et al., as well as Kim and Arakawa and Kim and Doyle. Part I demonstrated that the model with upgraded middle-atmospheric orographic drag physics better forecasts the magnitude and evolution of the SSW and better simulates the trend of the Arctic Oscillation (AO) index. In this study (Part II), a series of 5-day hindcast experiments is performed with cycling data assimilation using the Naval Research Laboratory Atmospheric Variational Data Assimilation System-Accelerated Representer (NAVDAS-AR), a four-dimensional variational data assimilation (4DVAR) system. Further efforts are made to improve the hindcasting of SSW by improving the satellite radiance bias correction process that strongly affects the data assimilation. The innovation (observation minus background) limit is optimally determined to reduce the rejection of useful radiance data. It is found that when the innovation limit is properly set, both the analysis and forecast of the SSW event can be improved, and that the orographic drag helps improve the SSW forecast.

Four-Dimensional Assimilation of Multitime Wind Profiles over a Single Station and Numerical Simulation of a Mesoscale Convective System Observed during IHOP_2002

Abstract This study examines the impact of assimilating multitime wind profiles over a single station on the numerical simulation of a warm season mesoscale convective system over the region from the Kansas and Oklahoma border to the Texas Panhandle, observed 12–13 June 2002 during the International H2O Project (IHOP_2002). Wind profile observations, obtained from Goddard Lidar Observatory for Winds (GLOW) are assimilated into an advanced research version of the Weather Research and Forecasting (WRF) model using its four-dimensional variational data assimilation (4DVAR) system. Results indicate that the assimilation of high temporal and vertical resolution GLOW wind profiles has a significant influence on the numerical simulation of the convective initiation and evolution. Besides the wind fields, the structure of the moisture fields associated with the convective system is also improved. Data assimilation has also resulted in a more accurate prediction of the locations and timing of the convection initiations; as a consequence, the skill of quantitative precipitation forecasting is enhanced greatly. The positive impact of 4DVAR assimilation of multitime wind profiles over a single station on the mesoscale prediction in this study presents a successful procession of the traditional technique in time to space conversion. However, when the data from conventional networks are assimilated into the model with GLOW wind profiles, the data impact is not compatible with that from the assimilation of GLOW wind profiles only, implying the need for a high temporal and spatial resolution wind profile network in order to achieve reasonable mesoscale analysis and forecasting.

Rainfall assimilation using a new four-dimensional variational method: A single-point observation experiment

Accurate forecast of rainstorms associated with the mei-yu front has been an important issue for the Chinese economy and society. In July 1998 a heavy rainstorm hit the Yangzi River valley and received widespread attention from the public because it caused catastrophic damage in China. Several numerical studies have shown that many forecast models, including Pennsylvania State University National Center for Atmospheric Research’s fifth-generation mesoscale model (MM5), failed to simulate the heavy precipitation over the Yangzi River valley. This study demonstrates that with the optimal initial conditions from the dimension-reduced projection four-dimensional variational data assimilation (DRP-4DVar) system, MM5 can successfully reproduce these observed rainfall amounts and can capture many important mesoscale features, including the southwestward shear line and the low-level jet stream. The study also indicates that the failure of previous forecasts can be mainly attributed to the lack of mesoscale details in the initial conditions of the models.

A Cloud-Resolving 4DVAR Assimilation Experiment for a Local Heavy Rainfall Event in the Tokyo Metropolitan Area

Abstract A cloud-resolving nonhydrostatic four-dimensional variational data assimilation system (NHM-4DVAR) was modified to directly assimilate radar reflectivity and applied to a data assimilation experiment using actual observations of a heavy rainfall event. Modifications included development of an adjoint model of the warm rain process, extension of control variables, and development of an observation operator for radar reflectivity. The responses of the modified NHM-4DVAR were confirmed by single-observation assimilation experiments for an isolated deep convection, using pseudo-observations of rainwater at the initial and end times of the data assimilation window. The results showed that the intensity of convection could be adjusted by assimilating appropriate observations of rainwater near the convection and that undesirable convection could be suppressed by assimilating small or no reflectivity. An assimilation experiment using actual observations of a local heavy rainfall in the Tokyo, Japan, metropolitan area was conducted with a horizontal resolution of 2 km. Precipitable water vapor derived from global positioning system data was assimilated at 5-min intervals within 30-min assimilation windows, and surface and wind profiler data were assimilated at 10-min intervals. Doppler radial wind and radar-reflectivity data below the elevation angle of 5.4° were assimilated at 1-min intervals. The 4DVAR assimilation reproduced a line-shaped rainband with a shape and intensity consistent with the observation. Assimilation of radar-reflectivity data intensified the rainband and suppressed false convection. The simulated rainband lasted for 1 h in the extended forecast and then gradually decayed. Sustaining the low-level convergence produced by northerly winds in the western part of the rainband was key to prolonging the predictability of the convective system.

Optimizing global CO emission estimates using a four-dimensional variational data assimilation system and surface network observations

Abstract. We apply a four-dimensional variational (4D-VAR) data assimilation system to optimize carbon monoxide (CO) emissions for 2003 and 2004 and to reduce the uncertainty of emission estimates from individual sources using the chemistry transport model TM5. The system is designed to assimilate large (satellite) datasets, but in the current study only a limited amount of surface network observations from the National Oceanic and Atmospheric Administration Earth System Research Laboratory (NOAA/ESRL) Global Monitoring Division (GMD) is used to test the 4D-VAR system. By design, the system is capable to adjust the emissions in such a way that the posterior simulation reproduces background CO mixing ratios and large-scale pollution events at background stations. Uncertainty reduction up to 60 % in yearly emissions is observed over well-constrained regions and the inferred emissions compare well with recent studies for 2004. However, with the limited amount of data from the surface network, the system becomes data sparse resulting in a large solution space. Sensitivity studies have shown that model uncertainties (e.g., vertical distribution of biomass burning emissions and the OH field) and the prior inventories used, influence the inferred emission estimates. Also, since the observations only constrain total CO emissions, the 4D-VAR system has difficulties in separating anthropogenic and biogenic sources in particular. The inferred emissions are validated with NOAA aircraft data over North America and the agreement is significantly improved from the prior to posterior simulation. Validation with the Measurements Of Pollution In The Troposphere (MOPITT) instrument version 4 (V4) shows a slight improved agreement over the well-constrained Northern Hemisphere and in the tropics (except for the African continent). However, the model simulation with posterior emissions underestimates MOPITT CO total columns on the remote Southern Hemisphere (SH) by about 10 %. This is caused by a reduction in SH CO sources mainly due to surface stations on the high southern latitudes.

Data Assimilation of Cloud-Affected Radiances in a Cloud-Resolving Model

AbstractAssimilation of cloud-affected infrared radiances from the Geostationary Operational Environmental Satellite-8 (GOES-8) is performed using a four-dimensional variational data assimilation (4DVAR) system designated as the Regional Atmospheric Modeling Data Assimilation System (RAMDAS). A cloud mask is introduced in order to limit the assimilation to points that have the same type of cloud in the model and observations, increasing the linearity of the minimization problem. A series of experiments is performed to determine the sensitivity of the assimilation to factors such as the maximum-allowed residual in the assimilation, the magnitude of the background error decorrelation length for water variables, the length of the assimilation window, and the inclusion of other data such as ground-based data including data from the Atmospheric Emitted Radiance Interferometer (AERI), a microwave radiometer, radiosonde, and cloud radar. In addition, visible and near-infrared satellite data are included in a separate experiment. The assimilation results are validated using independent ground-based data. The introduction of the cloud mask where large residuals are allowed has the greatest positive impact on the assimilation. Extending the length of the assimilation window in conjunction with the use of the cloud mask results in a better-conditioned minimization, as well as a smoother response of the model state to the assimilation.

A Frequent-Updating Analysis System Based on Radar, Surface, and Mesoscale Model Data for the Beijing 2008 Forecast Demonstration Project

Abstract The Variational Doppler Radar Analysis System (VDRAS) was implemented in Beijing, China, and contributed to the Beijing 2008 Forecast Demonstration Project (B08FDP) in support of the Beijing Summer Olympics. VDRAS is a four-dimensional variational data assimilation system that produces frequently updated analyses using Doppler radar radial velocities and reflectivities, surface observations, and mesoscale model data. The system was tested in real time during the B08FDP pretrials in the summers of 2006 and 2007 and run during the Olympics to assist the 0–6-h convective weather nowcasting. This paper provides a description of the upgraded system and its Beijing implementation, an evaluation of the system performance using data collected during the pretrials, and its utility on convective weather nowcasting through two case studies. Verification of VDRAS wind against a wind profiler shows that the analyzed wind is reasonably accurate with a smaller RMS difference for 2006 than for 2007 due to better radar data coverage in 2006. The analyzed cold pools in three convective episodes are compared with surface observations at selected stations. The result shows good agreement between the analysis and the observations. The two case studies demonstrate the role that VDRAS could play in nowcasting convective initiation.

The structure of background-error covariance in a four-dimensional variational data assimilation system: Single-point experiment

A four dimensional variational data assimilation (4DVar) based on a dimension-reduced projection (DRP-4DVar) has been developed as a hybrid of the 4DVar and Ensemble Kalman filter (EnKF) concepts. Its good flow-dependent features are demonstrated in single-point experiments through comparisons with adjoint-based 4DVar and three-dimensional variational data (3DVar) assimilations using the fifth-generation Pennsylvania State University-National Center for Atmospheric Research Mesoscale Model (MM5). The results reveal that DRP-4DVar can reasonably generate a background error covariance matrix (simply B-matrix) during the assimilation window from an initial estimation using a number of initial condition dependent historical forecast samples. In contrast, flow-dependence in the B-matrix of MM5 4DVar is barely detectable. It is argued that use of diagonal estimation in the B-matrix of the MM5 4DVar method at the initial time leads to this failure. The experiments also show that the increments produced by DRP-4DVar are anisotropic and no longer symmetric with respect to observation location due to the effects of the weather trends captured in its B-matrix. This differs from the MM5 3DVar which does not consider the influence of heterogeneous forcing on the correlation structure of the B-matrix, a condition that is realistic for many situations. Thus, the MM5 3DVar assimilation could only present an isotropic and homogeneous structure in its increments.