Bending Angle Data Research Articles

Evaluating the Impacts of COSMIC-2 GNSS RO Bending Angle Assimilation on Atlantic Hurricane Forecasts Using the HWRF Model

Abstract This study evaluates the impact of assimilating Global Navigation Satellite System (GNSS) radio occultation (RO) bending angles from Formosa Satellite Mission-7/Constellation Observing System for Meteorology, Ionosphere, and Climate-2 (COSMIC-2) receiver satellites on Hurricane Weather Research and Forecasting (HWRF) Model tropical cyclone (TC) forecasts. Launched in June 2019, the COSMIC-2 mission provides significantly higher tropics data coverage compared to its predecessor COSMIC constellation. GNSS RO measurements yield information about atmospheric pressure, temperature, and water vapor profiles. HWRF is cycled with and without COSMIC-2 bending angle data assimilation for six 2020 Atlantic hurricane cases. COSMIC-2 assimilation has little impact on HWRF track forecasts, consistent with HWRF’s design limiting cycled data assimilation impacts on surrounding large-scale flows; however, COSMIC-2 assimilation results in a statistically significant ∼8%–12% mean absolute forecast error reduction in minimum central sea level pressure for t = 36-, 54-, 60-, and 108–120-h lead times. Forecasts initialized from analyses assimilating COSMIC-2 observations also have a 1%–4% smaller 600–700-hPa specific humidity (SPFH) root-mean-squared deviation compared to radiosondes and dropwindsondes for most lead times. While not all HWRF intensity forecasts benefit from COSMIC-2 assimilation, a few show notable improvement. For example, assimilating two COSMIC-2 profiles within the inner core of developing Hurricane Hanna (2020) increases 800-hPa SPFH by up to 1 g kg−1 locally, helping to correct a dry bias. The forecast initialized from this analysis better captures Hanna’s observed intensification rate, likely because its moister inner core facilitates development of persistent deep convection near the TC center, where diabatic heating is more efficiently converted to cyclonic wind kinetic energy. Significance Statement Tropical cyclone (TC) intensification can be strongly sensitive to the lower-to-midtropospheric water vapor distribution near the storm. The COSMIC-2 GNSS radio occultation (RO) receiver satellite mission provides denser spatial coverage of atmospheric water vapor and temperature profiles over the tropics compared to other GNSS RO observation platforms. Herein, using six 2020 Atlantic TC cases, we evaluate the impacts of assimilating COSMIC-2 RO bending angles into a regional forecast model that already assimilates clear-sky satellite radiances. It is shown that COSMIC-2 assimilation yields a modest ∼10% intensity forecast skill improvement for several lead times, although more substantial intensity forecast improvement is found for a few forecasts where the COSMIC-2 observation assimilation helps correct a lower-to-midtropospheric water vapor bias.

A new method of retrieving atmospheric refractivity structure

ABSTRACT It is proposed to obtain information on atmospheric refractivity structure by measuring the angle of arrival (AoA) of radio signals routinely broadcast by commercial aircraft. The angle of arrival would be measured at hill-top sites using a simple two-element interferometer. Knowledge of the aircraft’s location (information conveniently contained within the broadcasts) and the AoA will enable the bending angle of the signals to be calculated. As measurable bending will only occur at grazing incidence, sources of signals either very close to the radio horizon, or at a similar height to the interferometer, are essential. The routine navigational data broadcasts from civil aircraft represent the ideal source. In areas of high air traffic density such as the UK, ~-bending angle measurements may be possible each day. Numerical weather prediction models routinely assimilate bending angles retrieved from GNSS radio occultation data, so it is anticipated that assimilation methods could be developed that are able to make good use of this new source of bending angle data. Sensitivity tests were performed to estimate the resolution of humidity retrievals assuming a target AoA accuracy of 0.01°. Simulated annealing was used to demonstrate the ability to retrieve relative humidity and mixing ratio vertical profiles using AoA measurements. It is shown that for observed AoA measurements with an accuracy of 0.01° it should be possible to retrieve relative humidity and mixing ratio vertical profiles with an accuracy of ~ and ~ g/kg respectively. An AoA accuracy of 0.01° should be achievable using hardware costing ~€10k, however further hardware development is still required.

Energy Harvesting and Sensing Integrated Woven Structure Kneepad Based on Triboelectric Nanogenerators

AbstractTriboelectric nanogenerators (TENGs) are extensively studied because of their great potential in energy harvesting and sensing. Traditional chemical batteries cannot meet the requirements of wearable electronic products because of their inherent limitations such as large weight, large size, limited endurance, and limited lifespan. To address the issue, the self‐powered energy harvester has attracted considerable interest among researchers. Here, the woven structure TENG consists of flexible silicone mixed with polyvinylidene fluoride as dielectric layers, silver fabric, and conductive tape as electrodes. The electrical output performance of the pattern structure processed by tip design is nearly three times higher. The whole TENG and flexible printed circuit board are encapsulated by polydimethylsiloxane film, reducing inconvenient, and discomfort sensation during wearing. The woven structure TENG can not only harvest energy to power wearable electronic devices but also be used to recognize human motion as the waveforms changing when people are in different states. Taking the knee bending angle data as an example to conduct machine learning, it is able to classify different states of human motion. This work extends the application of the woven TENG as an integrated device for sensing and energy harvesting.

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.

A Time Synchronization Method Between FBG Deformation Measurement System and Airborne DPOS Based on Waveform Matching

The airborne Distributed Position and Orientation System (DPOS) is an important means for multi-point motion measurement required for earth observation with multi-task sensing loads. The sub-nodes in DPOS rely on the main node to perform transfer alignment for them to get high-precision motion data, which requires accurate deformation data and the relative motion relationship between different nodes. Fiber Bragg Grating (FBG) is one of the main means of deformation measurement at present. However, the time synchronization error between FBG and DPOS directly restricts the accuracy of transfer alignment. For this, a high-precision time synchronization method based on waveform matching is presented. Firstly, the data correlation between FBG and DPOS is analyzed. Then the main time error factors, sampling frequency and time delay are determined. On this basis, the empirical mode decomposition and wavelet analysis are used to extract and match the vibration data from the FBGS and bending angle data calculated by DPOS, thus updating the time-varying sampling frequency and time delay of FBGS, and finally achieving high-precision time synchronization between FBGS and DPOS. Results of ground dynamic experiment show that the standard deviation and mean value of FBG time error is less than 4ms and 2ms respectively, and the deformation displacement error of synchronized FBG data is less than 1%, which provides a prerequisite for high-precision transfer alignment of DPOS.

Assimilation Impact of Early FORMOSAT-7/COSMIC-2 GNSS Radio Occultation Data with Taiwan’s CWB Global Forecast System

AbstractThe FORMOSAT-7/COSMIC-2 Global Navigation Satellite System (GNSS) Radio Occultation (RO) satellite constellation was launched on June 2019 as a successor of the FORMOSAT-3/COSMIC mission. The Central Weather Bureau (CWB) of Taiwan has received FORMOSAT-7/COSMIC-2 GNSS RO data in real time from Taiwan Analysis Center for COSMIC. With the global numerical prediction system at CWB, a parallel semi-operational experiment assimilating the FORMOSAT-7/COSMIC-2 bending angle data with all other operational observation data has been conducted to evaluate the impact of the FORMOSAT-7/COSMIC-2 data. The first seven-month results show that the quality of the early FORMOSAT-7/COSMIC-2 data has been satisfactory for assimilation. Consistent and significant positive impacts on global forecast skills have been observed since the start of the parallel experiment, with the most significant impact found in the tropical region, reflecting the low-inclination orbital design of the satellites. The impact of the FORMOSAT-7/COSMIC-2 RO data is also estimated using the Ensemble Forecast Sensitivity to Observation Impact (EFSOI) method, showing an average positive impact per observation similar to other existing GNSS RO datasets, while the total impact is impressive by virtue of its large amount. Sensitivity experiments suggest that the quality control processes built in the Gridpoint Statistical Interpolation (GSI) system for RO data work well to achieve a positive impact by the low-level FORMOSAT-7/COSMIC-2 RO data, while more effort on observation error tuning should be focused to obtain an optimal assimilation performance. This study demonstrates the usefulness of the FORMOSAT-7/COSMIC-2 RO data in global numerical weather prediction during the calibration/validation period and leads to the operational use of the data at CWB.

Experimental analysis on dynamic response of pre-cracked aluminum plate subjected to underwater explosion shock loadings

As the main defense component of the ship against the attack of external damage sources, the liquid cabin will inevitably suffer some damage like holes or cracks in the battle because of the explosion fragments or the penetration of high-speed projectiles. The research on the rest capacity and dynamic response of the cabin with different damage cracks due to the underwater explosions is of great significance to the improvement of the protection design of cabin structures and the development of the fluid-solid coupling theory. For the sake of simplicity, an aluminum thin-walled plate with central prefabricated crack was adopted to solve the problem on coupling loads induced by underwater explosion, and an underwater shock wave tube was used to substitute an underwater explosion load to generate the shock loading on thin-walled aluminum plate. The shock wave pressures in the underwater explosion tube and the dynamic deformation responses of aluminum thin-walled plates with central prefabricated crack were recorded by three pressure sensors and two high-speed cameras, respectively. By analyzing experimental results, the history characteristics of the dynamic increase of plastic bulge height and the crack propagation of the aluminum plates were obtained. The results show that the time of the crack growth and the bulge increasing are not synchronized in the dynamic response of thin-walled plate; in addition, the effects of the crack shape and size on the bulge and the crack growth in dynamic response of target plate can be attributed to the trade-off between the energy leakage and the bending stiffness of the target plate. Finally, according to the velocity and acceleration on the target and based on a plastic-rigid material model, we obtain that the target bending angle sine is linear with the radial distance, and the conclusion was well validated by experimental target bending angle data.

Assimilation of GPSRO Bending Angle Profiles into the Brazilian Global Atmospheric Model

The Global Positioning System (GPS) Radio Occultation (RO) technique allows valuable information to be obtained about the state of the atmosphere through vertical profiles obtained at various processing levels. From the point of view of data assimilation, there is a consensus that less processed data are preferable because of their lowest addition of uncertainties in the process. In the GPSRO context, bending angle data are better to assimilate than refractivity or atmospheric profiles; however, these data have not been properly explored by data assimilation at the CPTEC (acronym in Portuguese for Center for Weather Forecast and Climate Studies). In this study, the benefits and possible deficiencies of the CPTEC modeling system for this data source are investigated. Three numerical experiments were conducted, assimilating bending angles and refractivity profiles in the Gridpoint Statistical Interpolation (GSI) system coupled with the Brazilian Global Atmospheric Model (BAM). The results highlighted the need for further studies to explore the representation of meteorological systems at the higher levels of the BAM model. Nevertheless, more benefits were achieved using bending angle data compared with the results obtained assimilating refractivity profiles. The highest gain was in the data usage exploring 73.4% of the potential of the RO technique when bending angles are assimilated. Additionally, gains of 3.5% and 2.5% were found in the root mean square error values in the zonal and meridional wind components and geopotencial height at 250 hPa, respectively.

Determining Atmospheric Boundary Layer Height with the Numerical Differentiation Method Using Bending Angle Data from COSMIC

This paper presents a new method to estimate the height of the atmospheric boundary layer (ABL) by using COSMIC radio occultation bending angle (BA) data. Using the numerical differentiation method combined with the regularization technique, the first derivative of BA profiles is retrieved, and the height at which the first derivative of BA has the global minimum is defined to be the ABL height. To reflect the reliability of estimated ABL heights, the sharpness parameter is introduced, according to the relative minimum of the BA derivative. Then, it is applied to four months of COSMIC BA data (January, April, July, and October in 2008), and the ABL heights estimated are compared with two kinds of ABL heights from COSMIC products and with the heights determined by the finite difference method upon the refractivity data. For sharp ABL tops (large sharpness parameters), there is little difference between the ABL heights determined by different methods, i.e., the uncertainties are small; whereas, for non-sharp ABL tops (small sharpness parameters), big differences exist in the ABL heights obtained by different methods, which means large uncertainties for different methods. In addition, the new method can detect thin ABLs and provide a reference ABL height in the cases eliminated by other methods. Thus, the application of the numerical differentiation method combined with the regularization technique to COSMIC BA data is an appropriate choice and has further application value.

Comparison study of COSMIC RO dry-air climatologies based on average profile inversion

Abstract. Global Navigation Satellite System (GNSS) radio occultation (RO) data enable the retrieval of near-vertical profiles of atmospheric parameters like bending angle, refractivity, pressure, and temperature. The retrieval step from bending angle to refractivity, however, involves an Abel integral with an upper limit of infinity. RO data are practically limited to altitudes below about 80 km and the observed bending angle profiles show decreasing signal-to-noise ratio with increasing altitude. Some kind of high-altitude background data are therefore needed in order to perform this retrieval step (this approach is known as high-altitude initialization). Any bias in the background data will affect all RO data products beyond bending angle. A reduction of the influence of the background is therefore desirable – in particular for climate applications. Recently a new approach for the production of GNSS radio occultation climatologies has been proposed. The idea is to perform the averaging of individual profiles in bending angle space and then propagate the mean bending angle profiles through the Abel transform. Climatological products of refractivity, density, pressure, and temperature are directly retrieved from the mean bending angles. The averaging of a large number of profiles suppresses noise in the data, enabling observed bending angle data to be used up to 80 km without the need for a priori information. Some background information for the Abel integral is still necessary above 80 km. This work is a follow-up study, having the focus on the comparison of the average profile inversion climatologies (API) from the two processing centers WEGC and DMI, which study monthly COSMIC (Constellation Observing System for Meteorology, Ionosphere, and Climate) data from January to March 2011. The impact of different backgrounds above 80 km is tested, and different implementations of the Abel integral are investigated. Results are compared for the climatological products with ECMWF analyses, MIPAS, and SABER data. It is shown that different implementations of the Abel integral have little impact on the API climatologies. On the other hand, different extrapolations of the bending angle profile above 80 km play a key role in the resulting monthly mean refractivities above 35 km in altitude. Below that respective altitude the API climatologies show a good agreement between the two processing centers WEGC and DMI. Due to the downward propagation within the retrieval, effects of the high-altitude initialization lead to differences in dry-temperature climatologies down to 20 km in altitude. When applying an exponential extrapolation to the bending angles above 80 km at both centers, the dry-temperature climatologies agree among WEGC, DMI, ECMWF analysis, and MIPAS up to 35 km in altitude within ±0.5 K and up to 40 km in altitude within ±1 K. We conclude that the API retrieval is a valid approach up to the lower stratosphere. It is a computationally efficient alternative method for producing dry atmospheric RO climatologies.

The Impact of Vertical Resolution in the Assimilation of GPS Radio Occultation Data

Abstract The sensitivity of GPS radio occultation (GPSRO) bending angle assimilation to vertical resolution was studied within a global three-dimensional variational data assimilation (3DVAR) system. The sensitivity experiments were performed using different vertical resolutions of GPSRO data at 2 km, 1 km, 500 m, and 200 m. The assimilation of the higher vertical resolution GPSRO data showed better consistency in the analysis–forecast cycle in terms of the differences between GPSRO bending angle data and 6-h forecasts (O-F). This resulted in an improved analysis of the temperature, geopotential height, and wind in the mid-/upper-level troposphere by the hydrostatic response and the related model dynamics. It should be noted that the highest vertical resolution of the GPSRO data (200 m in this study) improved the forecasting skill level in terms of the root-mean-square error (against the European Centre for Medium-Range Weather Forecasts analysis) and the anomaly correlation of the geopotential height forecasting at 500 and 200 hPa in both the Northern and Southern Hemispheres. The benefits of assimilating higher vertical resolution GPSRO data were more pronounced in the upper-level troposphere, which was in agreement with previous studies using real GPSRO observations.

Evaluating the Impact of the COSMIC RO Bending Angle Data on Predicting the Heavy Precipitation Episode on 16 June 2008 during SoWMEX-IOP8

Abstract Global positioning system (GPS) radio occultation (RO) data have been broadly used in global and regional numerical weather predictions. Assimilation with the bending angle often performs better than refractivity, which is inverted from the bending angle under spherical assumption and is sometimes associated with negative biases at the lower troposphere; however, the bending angle operator also requires a higher model top as used in global models. This study furnishes the feasibility of bending-angle assimilation in the prediction of heavy precipitation systems with a regional model. The local RO operators for simulating bending angle and refractivity are implemented in the Weather Research and Forecasting (WRF)–local ensemble transform Kalman filter (LETKF) framework. The impacts of assimilating RO data from the Constellation Observing System for Meteorology Ionosphere and Climate (COSMIC) using both operators are evaluated on the prediction of a heavy precipitation episode during Southwest Monsoon Experiment intensive observing period 8 (SoWMEX-IOP8) in 2008. Results show that both the refractivity and bending angle provide a favorable condition for generating this heavy rainfall event. In comparison with the refractivity data, the advantage of assimilating the bending angle is identified in the midtroposphere for deepening of the moist layer that leads to a rainfall forecast closer to the observations.

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.

Systematic residual ionospheric errors in radio occultation data and a potential way to minimize them

Abstract. Radio occultation (RO) sensing is used to probe the earth's atmosphere in order to obtain information about its physical properties. With a main interest in the parameters of the neutral atmosphere, there is the need to perform a correction of the ionospheric contribution to the bending angle. Since this correction is an approximation to first order, there exists an ionospheric residual, which can be expected to be larger when the ionization is high (day versus night, high versus low solar activity). The ionospheric residual systematically affects the accuracy of the atmospheric parameters at low altitudes, at high altitudes (above 25–30 km) it even is an important error source. In climate applications this could lead to a time dependent bias which induces wrong trends in atmospheric parameters at high altitudes. The first goal of our work was to study and characterize this systematic residual error. In a second step we developed a simple correction method, based purely on observational data, to reduce this residual for large ensembles of RO profiles. In order to tackle this problem, we analyzed the bending angle bias of CHAMP and COSMIC RO data from 2001–2011. We could observe that the nighttime bending angle bias stays constant over the whole period of 11 yr, while the daytime bias increases from low to high solar activity. As a result, the difference between nighttime and daytime bias increases from about −0.05 μrad to −0.4 μrad. This behavior paves the way to correct the solar cycle dependent bias of daytime RO profiles. In order to test the newly developed correction method we performed a simulation study, which allowed to separate the influence of the ionosphere and the neutral atmosphere. Also in the simulated data we observed a similar increase in the bias in times from low to high solar activity. In this simulation we performed the climatological ionospheric correction of the bending angle data, by using the bending angle bias characteristics of a solar cycle as a correction factor. After the climatological ionospheric correction the bias of the simulated data improved significantly, not only in the bending angle but also in the retrieved temperature profiles.

Analysis of GPS radio occultation data from the FORMOSAT-3/COSMIC and Metop/GRAS missions at CDAAC

Abstract. This study investigates the noise level and mission-to-mission stability of Global Positioning System (GPS) radio occultation (RO) neutral atmospheric bending angle data at the UCAR COSMIC Data Analysis and Archive Center (CDAAC). Data are used from two independently developed RO instruments currently flying in orbit on the FORMOSAT-3/COSMIC (F3C) and Metop/GRAS (GNSS Receiver for Atmospheric Sounding) missions. The F3C 50 Hz RO data are post-processed with a single-difference excess atmospheric phase algorithm, and the Metop/GRAS 50 Hz closed loop and raw sampling (down-sampled from 1000 Hz to 50 Hz) data are processed with a zero-difference algorithm. The standard deviations of the F3C and Metop/GRAS bending angles from climatology between 60 and 80 km altitude from June–December 2009 are approximately 1.78 and 1.13 μrad, respectively. The F3C standard deviation reduces significantly to 1.44 μrad when single-difference processing uses GPS satellites on the same side of the spacecraft. The higher noise level for F3C bending angles can be explained by additional noise from the reference link phase data that are required with single-difference processing. The F3C and Metop/GRAS mean bending angles differences relative to climatology during the same six month period are statistically significant and have values of −0.05 and −0.02 μrad, respectively. A comparison of ~13 500 collocated F3C and Metop/GRAS bending angle profiles over this six month period shows a similar mean difference of ~0.02 ± 0.02 μrad between 30 and 60 km impact heights that is marginally significant. The observed mean difference between the F3C and Metop/GRAS bending angles of ~0.02–0.03 μrad is quite small and illustrates the high degree of re-produceability and mission independence of the GPS RO data at high altitudes. Collocated bending angles between two F3C satellites from early in the mission differ on average by up to 0.5% near the surface due to systematically lower signal-to-noise ratio for one of the satellites. Results from F3C and Metop/GRAS differences in the lower troposphere suggest the Metop/GRAS bending angles are negatively biased compared to F3C with a maximum of several percents near the surface in tropical regions. This bias is related to different tracking depths (deeper in F3C) and data gaps in Metop/GRAS which make it impossible to process the data from both missions in exactly the same way.

GRAS radio occultation on-board of Metop

The GRAS radio occultation instrument is flying on Metop-A and belongs to the EPS (EUMETSAT Polar System). GRAS observes GPS satellites in occultation. Within this work, validation of GRAS closed-loop bending angle data against co-located ECMWF profiles extracted from model fields and occultations from the COSMIC constellation of radio occultation instruments is shown. Results confirm the high data quality and robustness, where GRAS shows lower bending angle noise against ECMWF than COSMIC and in terms of occultations per day, one GRAS ≈ two COSMIC satellites. This is partly due to the operational setup of EPS. For the investigation we focus on two observation periods where updates in the ECMWF (March 2009) and COSMIC processing (October 2009) have improved the statistics further. Bending angles biases agree to within 0.5% against ECMWF and to within 0.1% against COSMIC after the updates for altitudes between 8 and 40 km. In addition, we also analyze the impact of the Metop orbit processing on the derived GRAS bending angle data, where different GPS and Metop orbit solutions are analyzed. Results show that a batch based orbit processing would improve in particular the bending angle bias behavior at higher altitudes. Requirements for the operational processing of GRAS data are briefly outlined, options to ease the use of other positioning system satellites in the near future are discussed. A simplified analysis on the observation of several of these systems, e.g. GPS and Galileo, from one platform shows that about 16% of occultations are found within 300 km, ±3 h, thus providing similar information. A constellation of 2 GRAS like instruments would have only about 10% close-by.

The impact of GPS radio occultation assimilation at the Met Office

AbstractMonitoring of global positioning system radio occultation (GPSRO) data was performed. By comparisons between coincident GRAS and FORMOSAT‐3/COSMIC occultations and coincident GPSRO refractivity and Vaisala RS92 radiosonde‐derived refractivity data, we show the consistency of GPSRO data from independent missions (and differing processing methods) and with independent insitu observation data. It was found that the bending angle and refractivity data from COSMIC is positively biased relative to that of GRAS by 0.2–0.5% in the stratosphere, which is where the GRAS refractivity agreed better with the RS92 radiosonde‐derived refractivity than COSMIC.Forecast impact experiments using GPSRO refractivity and bending angle measurements were conducted (using 1D operators). It is shown that refractivity data (from FORMOSAT‐3/COSMIC, GRAS, CHAMP and GRACE‐A) provided a large positive impact, with the greatest improvements in the Southern Hemisphere extratropics, where r.m.s. fit to radiosondes improved by up to 10% at 250 hPa. A refractivity to bending angle assimilation comparison was made. We found that bending angle assimilation (1D operator) generally provides a greater positive impact than refractivity, with a notable improvement in the Northern Hemisphere extratropics, where a ∼1% improvement in r.m.s. fit to radiosonde temperatures, geopotential heights and winds in the lower stratosphere was seen. © Crown Copyright 2010. Published by John Wiley & Sons, Ltd.

Quality Control, Error Analysis, and Impact Assessment of FORMOSAT-3/COSMIC in Numerical Weather Prediction

Following several years of experimentation with the GPS radio occultation technique, the 6-satellite FORMOSAT-3/COSMIC (F3C) mission was launched mid-2006 and has been collecting data since then. In this paper we present early findings of research performed at Meteo-France regarding the use of these data for assimilation in numerical weather prediction. Benefiting from the dense global coverage allowed by F3C refraction-induced observations, we first assess the quality of these data at four levels: bending angle, refractivity, refractivity lapse rate, and temperature. We compare them with calculations from Meteo-France numerical weather forecasts. Leaning from these various levels of data we devise quality control procedures that rely on the refractivity lapse rate. Applying a recent methodology developed in data assimilation we calculate observation bending angle error variances for our assimilation system. Using these new quality control procedures and observation error estimates we run an assimilation and forecast experiment with Meteo-France's operational global 4DVAR data assimilation system used as a reference. Our results indicate a very clear positive impact of the assimilation of F3C bending angle data in the Southern hemisphere for the prediction of geo potential heights and winds. We also observe an improvement in wind forecast skill in the Northern hemisphere, albeit such an improvement is smaller than in the Southern hemisphere.

Sensitivity of GPS occultation to the stratopause height

We scrutinize temperature profiles collected with radio occultation measurement for an imprint of the stratopause. In the retrieval step that integrates bending angle data to atmospheric refractivity, the falloff toward infinite altitude is constrained in a boundary condition with statistical optimization or extrapolation. We point to the coherence between temperature and density scale height as advocator for a fitting procedure that exploits data down to the stratopause. Data below have higher signal‐to‐noise ratio, but fitting them in costs additional parameters for the shifted sign in scale height gradient. On the basis of noise free simulation using a climatology covering all latitudes, seasons, and hours and on the basis of validation against data collected with weather balloons, laser imaging, and limb sounding, we find that adaptation to the fluctuating stratopause is crucial for the accuracy of the retrieved temperature. To facilitate investigation, the stratopause altitude was preset in the boundary condition according to climatological temperature. Biases between 10 and 50 hPa cancel out when the fitting interval is extended to 7.2 km above the stratopause, almost the average displacement of the maximum density scale height caused by nonlinearity in the temperature profile. This is a reminder that the stratopause should be read from bending angle, not temperature. Moreover, bias would be approaching zero at all levels if the a priori information for extrapolation had been lapse rate, not isothermal conditions. Keeping the model seed for temperature conversion to subsequent retrieval steps eliminates external information from the deconvolved refractivity. It will help argue for radio occultation as independent vehicle for climate monitoring.

Derivation of the Water Vapor Content from the GNSS Radio Occultations Observations

Abstract The present work investigates the possibility of retrieving humidity using the bending angle data obtained from radio occultation of GPS signals without additional external information. In particular, with the proposed approach, the dry pressure profiles are obtained by fitting the bending angles of the outer-troposphere layers (from h = h250K up to the stratopause) using the Hopfield dry atmosphere model. The ground pressure and temperature are the parameters of the model to be estimated. In the second step the humidity profiles are extracted by subtracting the contribution resulting from the dry atmosphere from the measured bending angles. Such derivation implies a complex mathematical treatment of the relationship between the bending angle and the refractivity, which is fully explained herein. Furthermore, the method was applied on Challenging Minisatellite Payload (CHAMP) profiles. The CHAMP profiles are achieved by applying heuristic retrieval algorithms based on the canonical transform. The algorithms are applied to minimize the negative refractivity bias that is observed for low-latitude GNSS RO. Thus, the results are shown and discussed in the second part of the paper. Finally, it is widely discussed how the proposed method is able to retrieve refractivity profiles without using the Abel inversion.