Bending Angle Profiles Research Articles

First Results of Airborne GNSS Radio Occultation Sounding From Airbus Commercial Aircraft

AbstractThe lack of high vertical resolution atmospheric thermodynamic structure observations inside or near major weather events impedes our understanding of physical processes and their predictability in numerical weather prediction (NWP) models. Airborne Global Navigation Satellite System (GNSS) radio occultation (airborne radio occultation [ARO]) has proven to be a viable remote sensing option to offer dense soundings near flight tracks. The global fleet of commercial aircraft already equipped with GNSS receivers could be leveraged to produce an unprecedented number of ARO soundings along global flight paths. Eleven cases of atmospheric bending angle and refractivity profiles were successfully retrieved and compared with the colocated European Center for Medium‐Range Weather Forecasting global reanalysis data. Good quality measurements are obtained with median refractivity differences less than 1% in the middle and upper troposphere, between 5.5 and 11.5 km. Given the use of aircraft data (e.g., Aircraft Meteorological DAta Relay) for data assimilation, incorporating ARO profiles would be a valuable addition, further enhancing the accuracy of aviation and weather forecasts.

УСВОЕНИЕ ДАННЫХ РАДИОЗАТМЕННЫХ НАБЛЮДЕНИЙ В ЧИСЛЕННОМ ПРОГНОЗЕ ПОГОДЫ

The paper discusses the principles of the assimilation of radio occultation observations into numerical weather prediction models. The most widely used approach is assimilation of the bending angle profiles obtained from the observations by the wave optics methods. The observation operators are constructed from the geometric optics equations and can be three-dimensional (3D), taking into account the actual field of meteorological parameters, or one-dimensional (1D), using only one vertical profile. In the approximation of weak gradients, the 3D operator is reduced to the two-dimensional (2D) one, which takes into account the section of the fields of meteorological parameters by the vertical plane of the occultation event. Another approach is the assimilation of the retrieved atmospheric refractivity profiles by means of a simplest 1D or more accurate 2D operator. Some examples of statistical comparisons of the profiles of temperature, pressure, and specific humidity obtained by variational assimilation with the GFS reanalysis data are given.

Understanding Ionospheric and Geomagnetic Effects on Residual Biases in Radio Occultation Data for Stratospheric Climate Monitoring

AbstractHigher‐order residual ionospheric errors in Global Navigation Satellite System (GNSS) radio occultation (RO) data can induce a systematic residual ionospheric bias (RIB) in RO bending angles, which can impact stratospheric climate monitoring. The main RIB causes are the ray path splitting of dual‐frequency GNSS signals, the electron density distribution along the ray path and at the RO receiver location, and the geomagnetic field. In this study, we investigate the ionospheric and geomagnetic effects on RO‐retrieved stratospheric bending angle and temperature profiles by inspecting the kappa and bi‐local RIB correction methods, the current state‐of‐the‐art methods, using multiple RO satellite mission data in different ionospheric and geomagnetic conditions. We find that, globally, the layer mean/median RIBs of the kappa and bi‐local correction methods exhibit similar behaviors for different RO missions: the estimated bending angle RIBs reach about −0.025/−0.01 and −0.024/−0.008 μrad in the upper (40–45 km) and lower (30–35 km) stratosphere, respectively, while the temperature RIBs reach about −1.0/−0.3 and −0.2/−0.1 K in these layers. However, in the equatorial day time region, the RIB statistics of the mission results diverge. Both the kappa and bi‐local RIBs increase in magnitude with increasing ionization and geomagnetic field strength but show no increase with an increasing degree of ionospheric asymmetry. Overall, the more refined bi‐local method has higher capacity than the comparatively simple kappa method to represent the variability of the ionospheric and geomagnetic conditions that affect RO events. This is important for regional‐scale studies, where the geomagnetic term can be of key relevance.

A Comparison of Atmospheric Boundary Layer Height Determination Methods Using GNSS Radio Occultation Data

The accurate determination of Atmospheric Boundary Layer Height (ABLH) is crucial in various atmospheric studies and practical applications. In this study, we present a comprehensive comparative analysis of five distinct methods for estimating ABLH using Global Navigation Satellite System (GNSS) Radio Occultation (RO) data. These methods encompass the use of bending angle and refractivity profiles, namely Minimum Gradient methods of the Bending Angle (MGBA) and Refractivity (MGR) profiles, breaking point, Wavelet Covariance Transform (WCT), and Double-Parameter Model Function (DPMF). GNSS-RO data from COSMIC-2 and Spire are used. To establish robust validation, radiosonde data are employed as a reference, ensuring the reliability of our findings. The results reveal notable variations in the performances of these ABLH estimation methods. Specifically, the MGBA, MGR, breaking point, and DPMF methods exhibit strong correlations with the reference data. Conversely, the WCT method displays weaker correlations, higher biases, and elevated root-mean-square-errors, suggesting limitations in capturing the true ABLH. Furthermore, we remove outlier screening to facilitate a comparison of the differences among the five methods. The WCT and DPMF methods can detect strong variations in the profiles near the Earth’s surface and consider them as ABLH. However, these variations are caused by errors. The MGBA method emerges as a reliable and stable option, while the WCT and DPMF methods should be used with caution due to the lower quality of the GNSS-RO profiles near the Earth’s surface.

Assessment of GRAS Ionospheric Measurements for Ionospheric Model Assimilation

We conducted a study to assess the GNSS (Global Navigation Satellite System) Receiver for Atmospheric Sounding (GRAS) ionospheric data quality by processing Radio Occultation (RO) observations of ionospheric products. The main objective of the study is to validate ionospheric data generated at EUMETSAT, such as ionospheric bending angle profiles, amplitude and phase scintillations, topside Total Electron Content (TEC) from MetOp-A GRAS instrument as well as generating and validating new ionospheric products derived from GRAS RO observations such as the TEC, rate of TEC and vertical electron density profiles. The assessment is conducted by comparing and evaluating the systematic differences between similar products from other Low Earth Orbit (LEO) satellite missions or from ground-based ionospheric measurements. The study confirms that the GNSS topside and RO observations recorded by the GRAS instrument onboard MetOp satellites are of good quality and are a valuable source of data for ionospheric research.

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.

Study of Ionospheric Bending Angle and Scintillation Profiles Derived by GNSS Radio-Occultation with MetOp-A Satellite

In this work, a dedicated campaign by MetOp-A satellite is conducted to monitor the ionosphere based on radio-occultation (RO) measurements provided by the onboard GNSS (Global Navigation Satellite System) Receiver for Atmospheric Sounding (GRAS). The main goal is to analyze the capabilities of the collected data to represent the bending angle and scintillation profiles of the ionosphere. We compare the MetOp-A products with those generated by other RO missions and explore the spatial/temporal distributions sensed by the MetOp-A campaign. Validation of dual frequency bending angles at the RO tangent points, S4 index, and Rate of the Total electron content Index (ROTI) is performed against independent products from Fengyun-3D and FORMOSAT-7/COSMIC-2 satellites. Our main findings constitute the following: (1) bending angle profiles from MetOp-A agree well with Fengyun-3D measurements; (2) bending angle distributions show a typical S-shape variation along the altitudes; (3) signatures of the sporadic E-layer and equatorial ionization anomaly crests are observed by the bending angles; (4) sharp transitions are observed in the bending angle profiles above ~200 km due to the transition of the daytime/nighttime in addition to the transition of the bottom-side/top-side; and (5) sporadic E-layer signatures are observed in the S4 index distributions by MetOp-A and FORMOSAT-7/COSMIC-2, with expected differences in magnitudes between the GPS (Global Positioning System) L1 and L2 frequencies.

Processing and Validation of the STAR COSMIC-2 Temperature and Water Vapor Profiles in the Neutral Atmosphere

The global navigation satellite system (GNSS) radio occultation (RO) is becoming an essential component of National Oceanic and Atmospheric Administration (NOAA) observation systems. The constellation observing system for meteorology, ionosphere, and climate (COSMIC) 2 mission and the Formosa satellite mission 7, a COSMIC follow-on mission, is now the NOAA’s backbone RO mission. The NOAA’s dedicated GNSS RO SAtellite processing and science Application Center (RO-SAAC) was established at the Center for Satellite Applications and Research (STAR). To better quantify how the observation uncertainty from clock error and geometry determination may propagate to bending angle and refractivity profiles, STAR has developed the GNSS RO data processing and validation system. This study describes the COSMIC-2 neutral atmospheric temperature and moisture profile inversion algorithms at STAR. We used RS41 and ERA5, and UCAR 1D-Var products (wetPrf2) to validate the accuracy and uncertainty of the STAR 1D-Var thermal profiles. The STAR-RS41 temperature differences are less than a few tenths of 1 K from 8 km to 30 km altitude with a standard deviation (std) of 1.5–2 K. The mean STAR-RS41 water vapor specific humidity difference and the standard deviation are −0.35 g/kg and 1.2 g/kg, respectively. We also used the 1D-Var-derived temperature and water vapor profiles to compute the simulated brightness temperature (BTs) for advanced technology microwave sounder (ATMS) and cross-track infrared sounder (CrIS) channels and compared them to the collocated ATMS and CrIS measurements. The BT differences of STAR COSMIC-2-simulated BTs relative to SNPP ATMS are less than 0.1 K over all ATMS channels.

Using GNSS Radio Occultation Data to Monitor Tropical Atmospheric Anomalies during the January–February 2009 Sudden Stratospheric Warming Event

We used Global Navigation Satellite System (GNSS) radio occultation (RO) temperature, density, and bending angle profiles to monitor tropical atmospheric anomalies during the January–February 2009 sudden stratospheric warming (SSW) event on a daily basis. We constructed RO anomaly profiles (tropical mean (30°S–30°N)) and gridded mean anomalies, as well as tropopause height and temperature anomalies. Based on the anomalies, we investigated the response time and region of the tropical atmosphere to SSW. It was found that the GNSS RO data were robust in monitoring tropical atmospheric anomalies during SSW. The tropical stratosphere revealed cooling simultaneously with polar stratospheric warming, although the magnitudes of the maximum tropical mean anomalies were 6–7 times smaller than the polar mean. Altitude variations showed that tropical stratospheric anomalies were largest within 35–40 km, which were 5 km higher than those in the polar region. On the onset day of 23 January, temperature anomalies over 0–30°N were mostly more than −5 K, which were larger than those of −2 K detected over the 0–30°S band, and the largest anomalies were detected over northern Africa with values more than −10 K. RO density and bending angle anomalies responded to SSW in a similar way as temperature but were 20 km higher. Following cooling, the tropical upper stratosphere and lower mesosphere revealed visible warming, with anomalies more than 10 K in the sector of 15°S–15°N. Tropopause anomalies revealed the largest variations over 20°N–30°N, further confirming that the extratropical region of the northern hemisphere is a key region for the dynamical coupling between the polar and tropical regions. Tropopause height anomalies had clear increase trends from 16 January to 8 February, with anomalies of the 20°N–30°N band that were −2 km on Jan 16 and increased to −0.5 km on Feb 6 with a variation of 1.5 km, while variations in other bands were within 0.5 km. Tropopause temperature anomalies had clear decrease trends over the same period, with anomalies at 20°N–30°N of 4 K on 16 January and decreasing to about −1 K on 8 February, while anomalies in other bands showed variations within 3 K.

Comparison of different machine learning approaches for tropospheric profiling based on COSMIC-2 data

Precise and reliable information on the tropospheric temperature and water vapor profiles play a key role in weather and climate studies. Among the sensors supporting the observations of the troposphere, one can distinguish the Global Navigation Satellite System Radio Occultation (RO) technique, which provides accurate and high-quality meteorological profiles. However, external knowledge about temperature is essential to estimate other physical atmospheric parameters. To overcome this constraint, I trained and evaluated four different machine learning models comprising Artificial Neural Network (ANN) and Random Forest regression algorithms, where no auxiliary meteorological data is needed. To develop the models, I employed 150,000 globally distributed (45°S–45°N) RO profiles between October 2019 and December 2020. Input vectors consisted of bending angle or refractivity profiles from the Formosa Satellite-7/Constellation Observing System for Meteorology, Ionosphere, and Climate-2 mission together with the month, hour, and latitude of the RO event. While temperature, pressure, and water vapor profiles derived from the modern ERA5 reanalysis and interpolated to the RO location served as the models’ targets. Evaluation on the testing data set revealed a good agreement between all model outputs and ERA5 targets, where slightly better statistics were noted for ANN and refractivity inputs. Vertically averaged root mean square error (RMSE) did not exceed 1.7 K for the temperature and reached around 1.4 hPa and 0.45 hPa for the total and water vapor pressures. Additional validation with 477 co-located radiosonde observations and the operational one-dimensional variational product showed slightly larger discrepancies with the mean RMSE of around 1.9 K, 1.9 hPa, and 0.5 hPa for the temperature, pressure, and water vapor, respectively.Graphical

Comparisons of Different Methods to Determine Starting Altitudes for Dry Air Atmosphere by GNSS-RO Data

Boundary profile evaluation (BPV) is an approach proposed in order to estimate water vapor content in the atmosphere. It exploits radio occultation (RO) observations of the signals emitted by the satellites of global navigation systems (GNSS) which are eclipsing (rising) as viewed by a low earth orbit satellite (LEO). BPV requires, as a preliminary step, the estimation of the dry background atmosphere model of refractivity (i.e., obtained from bending angle profiles) to be subtracted from the real observations in order to extract water vapor profiles. The determination of the lowest layer of the atmosphere over which the concentration of water vapor can be deemed negligible is particularly crucial for a correct application of the BPV method. In this study, we have applied three methods to set the starting altitudes for the dry air layers of the atmosphere: (1) by air temperature below some threshold values (for example, 250 K); (2) by “smooth” bending angle profiles in ROs; (3) by saturated water vapor pressure. These methods were tested with thermodynamic and bending angle profiles from 912 radiosonde excursions colocated with RO observations. For every dry air starting altitude we determined the best estimator from each of the three methods. In particular, by comparing those estimators with the quantiles and momenta of the dry air starting altitude distributions, we achieved improvements of up to 50% of the humidity profiles.

Simultaneous Radio Occultation Predictions for Inter-Satellite Comparison of Bending Angle Profiles from COSMIC-2 and GeoOptics

The Global Navigation Satellite System (GNSS) radio occultation (RO) is a remote sensing technique that uses International System of Units (SI) traceable GNSS signals for atmospheric limb soundings. The RO bending angle/sounding profiles are needed for assimilation in Numerical Weather Prediction (NWP) models, weather, climate, and space weather applications. Evaluating these RO data to ensure the high data quality for these applications is becoming more and more critical. This study presents a method for predicting radio occultation events, from which simultaneous radio occultation (SRO) for a pair of low-Earth-orbit (LEO) satellites on the limb to the same GNSS satellite can be obtained. The SRO method complements the Simultaneous Nadir Overpass (SNO) method (for nadir viewing satellite instruments), which has been widely used to inter-calibrate LEO to LEO and LEO to geosynchronous-equatorial-orbit (GEO) satellites. Unlike the SNO method, the SRO method involves three satellites: a GNSS and two LEO satellites with RO receivers. The SRO method allows for the direct comparison of bending angles when the simultaneous RO measurements for two LEO satellites receiving the same GNSS signal pass through approximately the same atmosphere within minutes in time and within less than 200 km of distance from each other. The prediction method can also be applied to radiosonde overpass prediction, and coordinate radiosonde launches for inter-comparisons between RO and radiosonde profiles. The main advantage of the SRO comparisons of bending angles is the significantly reduced uncertainties due to the much shorter time and smaller atmospheric path differences than traditional RO comparisons. To demonstrate the usefulness of this method, we present a comparison of the Constellation Observing System for Meteorology, Ionosphere, and Climate-2 (COSMIC-2) and GeoOpitcs RO profiles using SRO data for two time periods: Commercial Weather Data (CWD) data delivery order-1 (DO-1): 15 December 2020–15 January 2021 and CWD DO-2: 17 March 2021–31 August 2021. The results show good agreement in the bending angles between the COSMIC-2 RO measurements and those from GeoOptics, although systematic biases are also found in the inter-comparisons. Instrument and processing algorithm performances for the signal-to-noise ratio (SNR), penetration height, and bending angle retrieval uncertainty are also characterized. Given the efficiency of this method and the many RO measurements that are publicly and commercially available as well as the expansion of receiver capabilities to all GNSS systems, it is expected that this method can be used to validate/inter-calibrate GNSS RO measurements from different missions.

Inverting COSMIC-2 Phase Data to Bending Angle and Refractivity Profiles Using the Full Spectrum Inversion Method

The radio occultation technique provides stable atmospheric measurements that can work as a benchmark for calibrating and validating satellite-sounding data. Launched on 25 June 2019, the Constellation Observing System for Meteorology, Ionosphere, and Climate 2 and Formosa Satellite Mission 7 (COSMIC-2/FORMOSAT-7) are expected to produce about 5000 high-quality RO observations daily over the tropics and subtropics. COSMIC-2 constellation consists of 6 Low Earth Orbit (LEO) satellites in 24° inclination orbits at 720 km altitude and distributed mainly between 45°N to 45°S. The COSMIC-2 observations have uniform temporal coverage between 30°N to 30°S. This paper presents an independent inversion algorithm to invert COSMIC-2 geometry and phase data to bending angle and refractivity. We also investigate the quality of Global Navigation Satellite System (GNSS) and LEO position vectors derived from the UCAR COSMIC Data Analysis and Archive Center (CDAAC). The GNSS and LEO position vectors are stable with LEO position variations < 1.4 mm/s. The signal-to-noise ratio (SNR) on the L1 band ranges from 300–2600 v/v with a mean of 1600 v/v. The inversion algorithm developed at NOAA Center for Satellite Applications and Research (STAR) uses the Full Spectrum Inversion (FSI) method to invert COSMIC-2 geometry and phase data to bending angle and refractivity profiles. The STAR COSMIC-2 bending angle and refractivity profiles are compared with in situ radiosonde, the current COSMIC-2 products derived from CDAAC, and the collocated European Center for Medium-Range Weather Forecasts (ECMWF) climate reanalysis data ERA5. The mean bias at 8–40 km altitude among the UCAR, ERA5, and STAR is <0.1% for both bending and refractivity, with a standard deviation in the range of 1.4–2.3 and 0.9–1.1% for bending angles refractivity, respectively. In the lowest 2 km, the RO bias relative to ERA-5 shows a strong latitudinal and SNR dependence.

Using A Sliding Window Phase Matching Method for Imaging of GNSS Radio Occultation Signals

Global Navigation Satellite System Radio Occultation (GNSS-RO) is a technique used to sound the atmosphere and derive vertical profiles of refractivity. Signals from GNSS satellites are received in a low-Earth orbit, and they are then processed to produce bending angle profiles, from which meteorological parameters can be retrieved. Generating two-dimensional images in the form of spectrograms from GNSS-RO signals is commonly done to, for instance, investigate reflections or estimate signal quality in the lower troposphere. This is typically implemented using, e.g., the Short-Time Fourier Transform (STFT) to produce a time-frequency representation that is subsequently transformed to bending angle (BA) and impact height (IH) coordinates by non-linear mapping. In this paper, we propose an alternative method based on a straightforward extension of the Phase Matching (PM) operator to produce two-dimensional spectral images in the BA-IH domain by applying a sliding window. This Sliding Window Phase Matching (SWPM) method generates the spectral amplitude on an arbitrary grid in BA and IH, e.g., along the coordinate axes. To illustrate, we show both SWPM and STFT methods applied to operational MetOp-A data. For SWPM we use a constant window in the BA-dimension, whereas for STFT we use a conventional constant time window. We show that the SWPM method produces the same result as STFT when the same window length is used for both methods. The sample points in impact parameter and bending angle are those generated by and the main advantage is that SWPM offers the user a convenient way to freely sample the BA-IH space. The cost for this is processing time that is somewhat longer than implementations based on the Fast Fourier Transform, such as the STFT method.

The boundary layer height obtained by the spline numerical differentiation method using COSMIC GPS radio occultation data: A case study of the Qinghai-Tibet Plateau

We apply a numerical differentiation method for gradients of bending angle (BA) vertical profiles to determine the atmospheric boundary layer (ABL) height. The basic idea to obtain profiles of BA vertical gradient is to minimize an objective functional, with the Tikhonov regularization term being the second order derivative of the function to be determined; the regularization parameter is chosen as the standard deviation of the COSMIC Radio Occultation's BA data. The optimal BA profile is in the form of cubic spline function, and its derivative is calculated, and the height with the minimum derivative is defined as the ABL height. In order to illustrate the effectiveness of this method, two numerical experiments and comparison with radio sonde data are carried out, which show that with the addition of Gaussian error, the numerical regularization method combined with cubic spline function can adequately maintain more stable and accurate results. Then, the ABL heights of the Qinghai-Tibet Plateau are obtained by the numerical differentiation method based on COSMIC radio occultation data from 2007 to 2011 (January, April, July, and October). The overall results show that the ABL height on the Qinghai-Tibet Plateau is closely related to the terrain height and the seasonal change. By comparing the results with the numerical differentiation method and “zbalmax”, it is obvious that both techniques show good consistency of spatial distribution and temporal evolution.

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.

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.

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.

Simulation of liquid water and ice contributions to bending angle profiles in the radio occultation technique

Inversion of radio occultation profiles usually relies on the Abel transform in a spherically symmetric atmosphere. The main contribution to the refractive index comes from gaseous constituents such as pressure, temperature and water vapor. However, the impact of other particles in the troposphere, neglected in the retrieval, can become important during specific weather events causing the commonly applied assumptions to be invalid. In this study, Global Forecast System model is used to simulate the effect of cloud water on bending angle profiles. The information of solid and liquid fractions was assessed from the cloud mixing ratio based on a temperature-dependent separation. The signal was propagated between GPS satellite and low-Earth orbiter through tangent point placed in the vicinity of cloud water. The assumption of a spherically symmetric distribution for cloud refractivity fields is validated based on 2-dimensional simulations with multiple phase screen method. Using one year of forecast data for 2016, the most statistically probable examples were selected and analyzed as well as the uppermost impact of cloud significance was determined in terms of bending angle errors. Typical refractivity values induced by liquid clouds are within 2 ppm and result in single-spike profile structures. Corresponding fractional errors of the bending angle can reach −4% with respect to retrievals with zero clouds contributions. In the presence of unusual weather phenomena, simulations in a horizontally inhomogeneous atmosphere suggest a large dependency on propagation paths that should be considered in the modeling of cloud impact. Bending angles can be affected by −10% fractional error, which inverted with Abel transform corresponds to refractivity error of −2%. The highest fractional difference induced by the ice water is generally below −1% in the bending angle and −0.5% in the refractivity. The most significant cloud content was observed in the first 6 km above the Earth’s surface, whereas ice clouds can affect radio occultation profiles up to the top height of 16 km.