Assimilation Of Radio Occultation Data Research Articles

Detection of Superrefraction at the Top of the Atmospheric Boundary Layer from COSMIC-2 Radio Occultations

Abstract Superrefraction at the top of the atmospheric boundary layer introduces problems for assimilation of radio occultation data in weather models. A method of detection of superrefraction by spectral analysis of deep radio occultation signals introduced earlier has been tested using 2 years of COSMIC-2/FORMOSAT-7 radio occultation data. Our analysis shows a significant dependence of the probability of detection of superrefraction on the signal-to-noise ratio, which results in a certain sampling nonuniformity. Despite this nonuniformity, the results are consistent with the known global distribution of superrefraction (mainly over the subtropical oceans) and show some additional features and seasonal variations. Comparisons to the European Centre for Medium-Range Weather Forecasts analyses and limited set of radiosondes show reasonable agreement. Being an independent measurement, detection of superrefraction from deep radio occultation signals is complementary to its prediction by atmospheric models and thus should be useful for assimilation of radio occultation data in the atmospheric boundary layer.

Recent Impact of COSMIC-2 with Improved Radio Occultation Data Assimilation Algorithms

Abstract A Constellation Observing System for Meteorology, Ionosphere and Climate (COSMIC) follow-on constellation, COSMIC-2, was successfully launched into equatorial orbit on 24 June 2019. With an increased signal-to-noise ratio due to improved receivers and digital beam-steering antennas, COSMIC-2 is producing about 5000 high-quality radio occultation (RO) profiles daily over the tropics and subtropics. The initial evaluation of the impact of assimilating COSMIC-2 into NOAA’s Global Forecast System (GFS) showed mixed results, and adjustments to quality control procedures and observation error characteristics had to be made prior to the assimilation of this dataset in the operational configuration in May 2020. Additional changes in the GFS that followed this initial operational implementation resulted in a larger percentage of rejection (∼90%) of all RO observations, including COSMIC-2, in the mid- to lower troposphere. Since then, two software upgrades directly related to the assimilation of RO bending angle observations were developed. These improvements were aimed at optimizing the utilization of COSMIC-2 and other RO observations to improve global weather analyses and forecasts. The first upgrade was implemented operationally in September 2021 and the second one in November 2022. This study describes both RO software upgrades and evaluates the impact of COSMIC-2 with this most recently improved configuration. Specifically, we show that the assimilation of COSMIC-2 observations has a significant impact in improving temperature and winds in the tropics, though benefits also extend to the extratropical latitudes.

An Improved One-Dimensional Bending Angle Forward Operator for the Assimilation of Radio Occultation Profiles in the Lower Troposphere

Abstract Under very large vertical gradients of atmospheric refractivity, which are typical at the height of the planetary boundary layer, the assimilation of radio occultation (RO) observations into numerical weather prediction (NWP) models presents several serious challenges. In such conditions, the assimilation of RO bending angle profiles is an ill-posed problem, the uncertainty associated with the RO observations is higher, and the one-dimensional forward operator used to assimilate these observations has several theoretical deficiencies. As a result, a larger percentage of these RO observations are rejected at the NWP centers by existing quality control procedures, potentially limiting the benefits of this data type to improve weather forecasting in the lower troposphere. To address these problems, a new methodology that enables the assimilation of RO data to be extended to the lower moist troposphere has been developed. Challenges associated with larger atmospheric gradients of refractivity are partially overcome by a reformulation that has minimal effect at higher altitudes. As a first step toward this effort, this study presents both the theoretical development of this new methodology and a forecast impact assessment of it using the NCEP NWP system. Though using a conservative approach, benefits in the lower tropical troposphere are already noticeable. The encouraging results of this work open the potential for further exploitation and optimization of RO assimilation.

Potential Impacts of Radio Occultation Data Assimilation on Forecast Skill of Tropical Cyclone Formation in the Western North Pacific

AbstractThis study assesses the potential influence of global navigation satellite system (GNSS) radio occultation (RO) data assimilation on the forecast skill of tropical cyclone formation over the western North Pacific in September–October 2019 through a regional model. Data from the Constellation Observing System for Meteorology, Ionosphere, and Climate mission II are applied. The forecast skill considers the hits and misses for nine developing cases and the false alarms and correct negatives for 23 non‐developing cases. Forecasts assimilating GNSS RO data reduce the false alarm ratio by 20% and increase the accuracy rate by 19%, compared to forecasts without GNSS RO data. Assimilation of GNSS RO data increases mid‐level moisture around the disturbance centers at the initial time of the forecasts. It also increases low‐level vorticity for developing cases but decreases vorticity throughout most of the troposphere for non‐developing cases. These lead to improved forecast performance for tropical cyclone formation.

Assimilation of Radio Occultation Data Using Measurement-Based Observation Error Specification: Preliminary Results

Abstract The local spectral width (LSW) of a radio occultation (RO) observation in impact parameter representation is a useful parameter for providing information on the uncertainty associated with the RO bending angle measurement. The LSW can potentially be used to specify the bending angle observation error (BaOE) in the lower troposphere for each individual sounding. This study assesses the usefulness and limitations of LSW in representing BaOE for a global data assimilation system. A two-step scheme is proposed to derive profile-dependent BaOE from LSW. Since the LSW-based BaOE varies with each individual RO observation, it is here designated as a dynamic BaOE (DBOE) in contrast to the traditional statistics-based BaOE specification. A benchmark control run and two sensitivity experiments are conducted with continuous cycling data assimilation using the National Centers for Environmental Prediction (NCEP) Gridpoint Statistical Interpolation (GSI) and Global Forecast System (GFS). The usefulness and impact of the LSW-based DBOE are evaluated using radiosonde observations and global analyses. Results show that DBOE is able to improve the assimilation of RO data, leading to better forecast skill scores. Another experiment, in which the GSI statistical observation error of the benchmark run is replaced by the average of LSW-based DBOE, shows that the ability to assign larger weighting for high-quality observation and lower weighting for low-quality observation is the key factor for the success of the LSW-based DBOE.

Estimation of refractivity uncertainties and vertical error correlations in collocated radio occultations, radiosondes, and model forecasts

Abstract. Random uncertainties and vertical error correlations are estimated for three independent data sets. The three collocated data sets are (1) refractivity profiles of radio occultation measurements retrieved from the Metop-A and B and COSMIC-1 missions, (2) refractivity derived from GRUAN-processed RS92 sondes, and (3) refractivity profiles derived from ERA5 forecast fields. The analysis is performed using a generalization of the so-called three-cornered hat method to include off-diagonal elements such that full error covariance matrices can be calculated. The impacts from various sources of representativeness error on the uncertainty estimates are analysed. The estimated refractivity uncertainties of radio occultations, radiosondes, and model data are stated with reference to the vertical representation of refractivity in these data sets. The existing theoretical estimates of radio occultation uncertainty are confirmed in the middle and upper troposphere and lower stratosphere, and only little dependence on latitude is found in that region. In the lower troposphere, refractivity uncertainty decreases with latitude. These findings have implications for both retrieval of tropospheric humidity from radio occultations and for assimilation of radio occultation data in numerical weather prediction models and reanalyses.

Impacts of Radio Occultation Data on Typhoon Forecasts as Explored by the Global MPAS-GSI System

Global Navigation Satellite System (GNSS) radio occultation (RO) provides plentiful sounding profiles over regions lacking conventional observations. The Gridpoint Statistical Interpolation (GSI) hybrid system for assimilating RO data is integrated in this study with the Model for Prediction Across Scales–Atmosphere (MPAS) to improve tropical cyclone forecasts. After the MPAS-GSI assimilation cycles, dynamical vortex initialization (DVI) that may effectively spin up the initial inner typhoon vortex through cycled model integration is implemented to improve the initial analysis fit to the best track position as well as maximum wind or pressure intensity for Typhoon Nepartak (2016) that moved northwestward toward southern Taiwan. During the cycling assimilation, assimilation with RO data improves the temperature and moisture analysis, and largely reduces the forecast errors compared to those without RO data assimilation. The two RO operators that assimilate local bending angle or refractivity produce similar analyses, but the temperature and moisture increments from bending angle assimilation are slightly larger than those from refractivity assimilation. The MPAS forecasts at 60-15 km resolution show that the typhoon track prediction is improved with RO data, especially using bending angle data. The reduction in track deviations is explained by the wavenumber-one potential vorticity budget for several forecasts associated with the track deflection near southern Taiwan. Assimilation of RO data has fewer impacts on the typhoon intensity forecast compared to the DVI that largely improves the initial and thus forecasted intensity of the typhoon but at the cost of a slightly degraded track. Use of the enhanced 3 km resolution in the typhoon path also further improved the forecasts with and without the DVI. The feasible performance of the MPAS-GSI system with the RO data impact is also illustrated for Typhoon Mitag (2019), that passed around northern Taiwan.

Impact of Radio Occultation Data on the Prediction of Typhoon Haishen (2020) with WRFDA Hybrid Assimilation

FORMOSAT-7/COSMIC-2 (FS7/C2) satellite was successfully launched in June 2019. The satellite provides about 5000 radio occultation (RO) soundings daily over the tropical and partial subtropical regions. Such plentiful RO soundings with high accuracy and vertical resolution could be used to improve model initial analysis for typhoon prediction. In this study, assimilation experiments with and without the RO data were conducted with the WRFDA hybrid system for the prediction of Typhoon Haishen (2020). The experimental results show a remarkable improvement in typhoon track prediction with RO data assimilation, especially when using a nonlocal refractivity operator. Results in cycling DA and forecast are analyzed and verified for the RO data impact. Diagnostics of potential vorticity (PV) tendency budget helps explain the typhoon translation induced by different physical processes in the budget. The typhoon translation is essentially dominated by horizontal PV advection, but the track deviation can increase due to the vertical PV advection with opposite effects in the absence of RO data. Sensitivity experiments for different model initial times, physics schemes, and RO observation amounts show positive RO data impacts on typhoon prediction, mainly contributed from FS7. Complementary, an improved forecast of Typhoon Hagupit (2020) is also illustrated for the RO data impact.

Evaluation of Operational Monsoon Moisture Surveillance and Severe Weather Prediction Utilizing COSMIC-2/FORMOSAT-7 Radio Occultation Observations

Operational monsoon moisture surveillance and severe weather prediction is essential for timely water resource management and disaster risk reduction. For these purposes, this study suggests a moisture indicator using the COSMIC-2/FORMOSAT-7 radio occultation (RO) observations and evaluates numerical model experiments with RO data assimilation. The RO data quality is validated by a comparison between sampled RO profiles and nearby radiosonde profiles around Taiwan prior to the experiments. The suggested moisture indicator accurately monitors daily moisture variations in the South China Sea and the Bay of Bengal throughout the 2020 monsoon rainy season. For the numerical model experiments, the statistics of 152 moisture and rainfall forecasts for the 2020 Meiyu season in Taiwan show a neutral to slightly positive impact brought by RO data assimilation. A forecast sample with the most significant improvement reveals that both thermodynamic and dynamic fields are appropriately adjusted by model integration posterior to data assimilation. The statistics of 17 track forecasts for typhoon Hagupit (2020) also show the positive effect of RO data assimilation. A forecast sample reveals that the member with RO data assimilation simulates better typhoon structure and intensity than the member without, and the effect can be larger and faster via multi-cycle RO data assimilation.

An Impact Study of GNSS RO Data on the Prediction of Typhoon Nepartak (2016) Using a Multiresolution Global Model with 3D-Hybrid Data Assimilation

AbstractConventional soundings are rather limited over the western North Pacific and can be largely compensated by GNSS radio occultation (RO) data. We utilize the GSI hybrid assimilation system to assimilate RO data and the multiresolution global model (MPAS) to investigate the RO data impact on prediction of Typhoon Nepartak that passed over southern Taiwan in 2016. In this study, the performances of assimilation with local RO refractivity and bending angle operators are compared for the assimilation analysis and typhoon forecast. Assimilations with both RO data have shown similar and comparable temperature and moisture increments after cycling assimilation and largely reduce the RMSEs of the forecast without RO data assimilation at later times. The forecast results at 60–15-km resolution show that RO data assimilation largely improves the typhoon track prediction compared to that without RO data assimilation, and assimilation with bending angle has better performances than assimilation with refractivity, in particular for wind forecast. The improvement in the forecasted track is mainly due to the improved simulation for the translation of the typhoon. Diagnostics of wavenumber-1 potential vorticity (PV) tendency budget indicates that the northwestward typhoon translation dominated by PV horizontal advection is slowed down by the southward tendency induced by the stronger differential diabatic heating south of the typhoon center for bending-angle assimilation. Simulations with the enhanced resolution of 3 km in the region of the storm track show further improvements in both typhoon track and intensity prediction with RO data assimilation. Positive RO impacts on track prediction are also illustrated for two other typhoons using the MPAS-GSI system.

The Impact of GPS RO Data on the Prediction of Tropical Cyclogenesis Using a Nonlocal Observation Operator: An Initial Assessment

Abstract In this study, the impact of global positioning system (GPS) radio occultation (RO) data on the prediction of the genesis of 10 tropical cyclones over the western North Pacific Ocean is assessed. With the use of a nonlocal excess phase observation operator in cycling data assimilation, the probability of detection for tropical cyclogenesis is increased from 30% to 70% for the cases considered, all of which developed into typhoons. However, the probability of detection is only increased to 40% when a local observation operator is used, indicating that the observation operator can significantly influence the performance of RO data assimilation in capturing tropical cyclogenesis. A nonlocal excess phase operator, which considers the atmospheric horizontal gradients by integrating the refractivity along a ray path, gives superior performance over the local observation operator. Additional sensitivity experiments on 3 of the 10 typhoon cases show that the RO data in the vicinity of the incipient cyclones (within 500 km of the cyclone center) are most critical to successful cyclogenesis prediction. This reflects the fact that having good RO observations at the right time and place is critical for RO to have beneficial impacts on tropical cyclogenesis. Further analyses for Typhoon Nuri (2008) show that assimilation of RO data using the nonlocal operator leads to moistening of the lower and middle troposphere, organized convection, robust grid-scale vertical motions, and the development of midlevel relative vorticity, all of which are favorable for tropical cyclogenesis.

Assimilation of GNSS tomography products into the Weather Research and Forecasting model using radio occultation data assimilation operator

Abstract. From Global Navigation Satellite Systems (GNSS) signals, accurate and high-frequency atmospheric parameters can be determined in all-weather conditions. GNSS tomography is a technique that takes advantage of these parameters, especially of slant troposphere observations between GNSS receivers and satellites, traces these signals through a 3-D grid of voxels, and estimates by an inversion process the refractivity of the water vapour content within each voxel. In the last years, the GNSS tomography development focused on numerical methods to stabilize the solution, which has been achieved to a great extent. Currently, we are facing new challenges and possibilities in the application of GNSS tomography in numerical weather forecasting, the main research objective of this paper. In the first instance, refractivity fields were estimated using two different GNSS tomography models (TUW, WUELS), which cover the area of central Europe during the period of 29 May–14 June 2013, when heavy-precipitation events were observed. For both models, slant wet delays (SWDs) were calculated based on estimates of zenith total delay (ZTD) and horizontal gradients, provided for 88 GNSS sites by Geodetic Observatory Pecny (GOP). In total, three sets of SWD observations were tested (set0 without compensation for hydrostatic anisotropic effects, set1 with compensation of this effect, set2 cleaned by wet delays outside the inner voxel model), in order to assess the impact of different factors on the tomographic solution. The GNSS tomography outputs have been assimilated into the nested (12 and 36 km horizontal resolution) Weather Research and Forecasting (WRF) model, using its three-dimensional variational data assimilation (WRFDA 3D-Var) system, in particular, its radio occultation observation operator (GPSREF). As only total refractivity is assimilated in GPSREF, it was calculated as the sum of the hydrostatic part derived from the ALADIN-CZ model and the wet part from the GNSS tomography. We compared the results of the GNSS tomography data assimilation to the radiosonde (RS) observations. The validation shows the improvement in the weather forecasting of relative humidity (bias, standard deviation) and temperature (standard deviation) during heavy-precipitation events. Future improvements to the assimilation method are also discussed.

Sensitivity of Forward-Modeled Bending Angles to Vertical Interpolation of Refractivity for Radio Occultation Data Assimilation

Assimilation of radio occultation (RO) observations into numerical weather prediction (NWP) models has improved forecasts, where RO is typically one of the top five observational systems contributing to forecast accuracy. By measuring the phase delay of radio waves traversing Earth’s atmosphere between global positioning system (GPS) and low-Earth orbiting satellites, RO obtains quasi-vertical profiles of bending angles (BA) of the radio waves’ trajectories. BA are the RO observation most often assimilated into NWP models, but since they are not computed or analyzed in the models, they must be computed from model data using a forward model. First, model refractivity N is computed from variables specified on the model’s vertical levels, then using the Abel integral, BA are computed from N. The forward model requires vertical differentiation of N, and accurate differentiation requires vertical interpolation of N between model levels. The interpolation results in errors, which then propagate through the forward model to produce BA errors. In this study, we investigate the sensitivity of forward-modeled BA to five different methods of vertical interpolation of N between model levels to determine a method that minimizes interpolation errors. We use RO-observed N to isolate the interpolation errors and determine an accurate method that can be applied to any NWP model. Of the five methods investigated, the log-spline interpolation reduces N and BA errors the most, regardless of the vertical resolution of the model grid.

A Quality Control Procedure Based on Bending Angle Measurement Uncertainty for Radio Occultation Data Assimilation in the Tropical Lower Troposphere

AbstractThe fluctuation of radio occultation (RO) signals in the presence of refractivity irregularities in the moist lower troposphere results in uncertainties of retrieved bending angle and refractivity profiles. In this study the local spectral width (LSW) of RO signals, transformed to impact parameter representation, is used for the characterization of the uncertainty (random error) of retrieved bending angle and refractivity profiles. A large LSW has some correlation with the large mean difference (bias) of retrieved refractivity and bending angle from radiosondes and European Centre for Medium-Range Weather Forecasts analyses based on data from 2008 to 2014. An LSW-based quality control (QC) procedure is developed to eliminate low-quality (large random errors and biases) profiles from data assimilation. The LSW-based QC procedure is tested and evaluated in the assimilation of Constellation Observing System for Meteorology, Ionosphere and Climate RO data using the NCAR Data Assimilation Research Testbed and the Weather Research and Forecasting Model. Preliminary results, based on a 2-week data assimilation cycle, show that the LSW-based QC procedure improves water vapor analyses in the moist lower troposphere.

Inclusion of Global Navigation Satellite System radio occultation data into Center for Weather Forecast and Climate Studies Local Ensemble Transform Kalman Filter (LETKF) using the Radio Occultation Processing Package as an observation operator

ABSTRACTThe Local Ensemble Transform Kalman Filter (LETKF) is a system for data assimilation in the testing phase at the Center for Weather Forecast and Climate Studies of the Brazilian Institute for Space Research (CPTEC/INPE). This system introduces significant advantages over the traditional variational methodology in the inclusion of new observation systems because determining the observation operator does not involve inverse problems. The refractivity profile from the Global Navigation Satellite System (GNSS) radio occultation (RO) technique is an example of an additional source of unconventional observation with promising results over regions with low data coverage. To assimilate RO profiles using the LETKF system at CPTEC/INPE, the Radio Occultation Processing Package (ROPP) from Radio Occultation Meteorology Satellite Applications Facilities/European Organization for the Exploration of Meteorological Satellites (ROM SAF/EUMETSAT) was used as an observation operator. The goal of this work was to investigate the impact on the Atmospheric General Circulation Model of the CPTEC using ROPP as an operator of the assimilation of RO data into the LETKF system. The results showed that positive impacts were obtained in several evaluated variables, including humidity, particularly over South America, where the model used by the CPTEC produces better results.

Assimilation of global navigation satellite radio occultation observations in GRAPES: Operational implementation

This paper presents the design of an observation operator for assimilation of global navigation satellite system (GNSS) radio occultation (RO) refractivity and the related operational implementation strategy in the global GRAPES variational data assimilation system. A preliminary assessment of the RO data assimilation effect is performed. The results show that the RO data are one of the most important observation types in GRAPES, as they have a significant positive impact on the analysis and forecast at all ranges, especially in the Southern Hemisphere and the global stratosphere where in-situ measurements are lacking. The GRAPES model error cannot be controlled in the Southern Hemisphere without RO data being assimilated. In addition, it is found that the RO data play a key role in the stable running of the GRAPES global assimilation and forecast system. Even in a relatively simple global data assimilation experiment, in which only the conventional and RO data are assimilated, the system is able to run for more than nine months without drift compared with NCEP analyses. The analysis skills in both the Northern and Southern Hemispheres are still relatively comparable even after nine-month integration, especially in the stratosphere where the number of conventional observations decreases and RO observations with a uniform global coverage dominate gradually.

An Impact Assessment of GPS Radio Occultation Data on Prediction of a Rapidly Developing Cyclone over the Southern Ocean*

Abstract The impact of global positioning system (GPS) radio occultation (RO) data on an intense synoptic-scale storm that occurred over the Southern Ocean in December 2007 is evaluated, and a synoptic explanation of the assessed impact is offered. The impact is assessed by using the three-dimensional variational data assimilation scheme (3DVAR) of the Weather Research and Forecasting (WRF) Model Data Assimilation system (WRFDA), and by comparing two experiments: one with and the other without assimilating the refractivity data from four different RO missions. Verifications indicate significant positive impacts of the RO data in various measures and parameters as well as in the track and intensity of the Antarctic cyclone. The analysis of the atmospheric processes underlying the impact shows that the assimilation of the RO data yields substantial improvements in the large-scale circulations that in turn control the development of the Antarctic storm. For instance, the RO data enhanced the strength of a 500-hPa trough over the Southern Ocean and prevented the katabatic flow near the coast of East Antarctica from an overintensification. This greatly influenced two low pressure systems of a comparable intensity, which later merged together and evolved into the major storm. The dominance of one low over the other in the merger dramatically changed the track, intensity, and structure of the merged storm. The assimilation of GPS RO data swapped the dominant low, leading to a remarkable improvement in the subsequent storm’s prediction.

Observation, analysis, and modeling of deep radio occultation signals: Effects of tropospheric ducts and interfering signals

AbstractGPS radio occultation (RO) signals are sometimes observed very deep in the Earth's shadow. To investigate these phenomena, one of the FORMOSAT‐3/COSMIC RO receivers was set to track RO signals deep below the limb, down to a height of straight line −350 km on 5–6 October 2010. Analysis of the spectrograms revealed the existence of two types of signals below −200 km, RO signals induced by tropospheric propagation and interfering signals not transmitted by the occulted GPS. The RO signals induced by tropospheric propagation arrive from impact heights corresponding to inversion layers. Wave optics modeling of RO signals showed that deep signals exist when the refractivity gradient exceeds critical (super‐refraction). The existence of such signals is a diffractional phenomenon, which offers a new quality control parameter to identify occultations that may be affected by super‐refraction. This is important for RO data assimilation in weather models in the moist lower troposphere because assimilation of RO data affected by super‐refraction is an ill‐conditioned problem. Detection of the tropospheric ducts also may be useful for evaluation of radio wave propagation conditions. For infinitely horizontally extended ducts, the deep signals are extended in duration, have amplitudes of about 0.1%, and exist for only elevated ducts. For ducts of limited horizontal extension, the deep signals are shorter in duration, have amplitude of about 1%, and may exist for both elevated and surface ducts. The interfering signals were found in about half of occultations. Based on frequency modeling, in most cases, the interfering signal was identified with non‐occulted GPS. The disturbance of retrieved bending angle induced by an interfering signal from a non‐occulted GPS in a region of strong defocusing and significant spectral spread of RO signal was modeled and determined to be quite large, up to 10% of bending angle. However, the probability of occurrence of such interference (not estimated), generally, must be low. Removal of the interfering signal by high‐pass filtering was tested, but it introduced a disturbance of bending angle of about the same magnitude as the interfering signal itself.

Parallelization Strategies for the GPS Radio Occultation Data Assimilation with a Nonlocal Operator in the Weather Research and Forecasting Model

Abstract The nonlocal excess phase observation operator for assimilating the global positioning system (GPS) radio occultation (RO) sounding data has been proven by some research papers to produce significantly better analyses for numerical weather prediction (NWP) compared to the local refractivity observation operator. However, the high computational cost and the difficulties in parallelization associated with the nonlocal GPS RO operator deter its application in research and operational NWP practices. In this article, two strategies are designed and implemented in the data assimilation system for the Weather Research and Forecasting Model to demonstrate the capability of parallel assimilation of GPS RO profiles with the nonlocal excess phase observation operator. In particular, to solve the parallel load imbalance problem due to the uneven geographic distribution of the GPS RO observations, round-robin scheduling is adopted to distribute GPS RO observations among the processing cores to balance the workload. The wall clock time required to complete a five-iteration minimization on a demonstration Antarctic case with 106 GPS RO observations is reduced from more than 3.5 h with a single processing core to 2.5 min with 106 processing cores. These strategies present the possibility of application of the nonlocal GPS RO excess phase observation operator in operational data assimilation systems with a cutoff time limit.

Assimilation of GPS Radio Occultation Refractivity Data with WRF 3DVAR and Its Impact on the Prediction of a Heavy Rainfall Event

AbstractTo accommodate accurate analyses and forecasts of a heavy rainfall event over the Korean Peninsula, the authors assimilated the GPS radio occultation (RO) data by using the Weather Research and Forecasting Model (WRF) and its three-dimensional variational data assimilation system (3DVAR). The employed datasets are from the Constellation Observing System for Meteorology, Ionosphere, and Climate (COSMIC) and Challenging Minisatellite Payload (CHAMP) missions. The selected case was from late October 2006, which intensively hit the northeastern part of the Korean Peninsula with record-breaking rainfall. In this study, the local refractivity observation operator was used in assimilating GPS RO soundings. The results are more pronounced for the cycling assimilation of GPS RO data than for the one-time data assimilation. From all of the parameters investigated (temperature, geopotential height, specific humidity, and winds), the GPS RO soundings highly modified the moisture distribution in the lower troposphere and also changed the wind field via the model dynamics. For the heavy rainfall forecast, the quantitative accuracy of the precipitation forecast with the GPS RO data assimilation was in good agreement with observations in terms of the maximum rainfall amount and threat scores. The improved forecast in the experiment came from the exact positioning of the low pressure system and its consequent convergence near the rainfall area. When RO data and GPS precipitable water data were assimilated simultaneously, the moisture distribution changed horizontally and vertically such that it increased the amount of rainfall, and an accurate description of the convective system development was feasible.