Formosa Satellite Mission Research Articles

Convective tropopause over the tropics: Climatology, seasonality, and inter-annual variability inferred from long-term FORMOSAT-3/COSMIC-1 RO data

The present study aims to provide a comprehensive analysis of the convective tropopause height (COTH) and temperature (COT-T) over the tropics using long-term radio occultation (RO) measurements from the Formosa Satellite Mission 3–Constellation Observing System for Meteorology, Ionosphere and Climate (FORMOSAT-3/COSMIC) from July 2006 to December 2018. The study reveals that the COT-H (COT-T) exhibits a significant longitudinal structure with a maximum (minimum) over the Western Pacific and two secondary maxima (minima) over the American and African regions. The seasonal cycles of COT-H and COT-T are also prominent, with higher and colder values over the Western Pacific in boreal winter and over the Indian region in boreal summer. This seasonality in COT-H and COT-T is consistent with the seasonal process of tropical convective activity and is explained by outgoing longwave radiation (OLR) data. The correlation analysis shows that COT-H (COT-T) is negatively (positively) correlated with OLR, with the highest correlation found in the Indian region. Additionally, the study investigates the relationship between COT-H (COT-T) and upper tropospheric water vapor (UTWV) over the various regions in the tropics. Over the Indian region, the COT-H (COT-T) exhibits a positive (negative) correlation with UTWV, while over the Western Pacific, American, and African regions, it displays a negative (positive) correlation. Finally, the spatial distribution of the convective tropopause in El Niño and La Niña winters was examined. Results show that the COT-H (COT-T) is higher (lower) over the central and eastern Pacific Ocean and lower (higher) over the maritime continent during El Niño compared with La Niña.

Characterizing the tropospheric water vapor spatial variation and trend using 2007–2018 COSMIC radio occultation and ECMWF reanalysis data

Abstract. Atmospheric water vapor plays a crucial role in the global energy balance, hydrological cycle, and climate system. High-quality and consistent water vapor data from different sources are vital for weather prediction and climate research. This study assesses the consistency between the Formosa Satellite Mission 3–Constellation Observing System for Meteorology, Ionosphere, and Climate (FORMOSAT-3/COSMIC) radio occultation (RO) and European Centre for Medium-Range Weather Forecasts (ECMWF) Reanalysis Model 5 (ERA5) water vapor datasets. Comparisons are made across different atmospheric pressure levels (300, 500, and 850 hPa) from 2007 to 2018. Generally, the two datasets show good spatial and temporal agreement. COSMIC's global water vapor retrieval is slightly lower than ERA5's at 500 and 850 hPa, with distinct latitudinal differences between hemispheres. COSMIC exhibits global water vapor increasing trends of 3.47 ± 1.77 % per decade, 3.25 ± 1.25 % per decade, and 2.03 ± 0.65 % per decade at 300, 500, and 850 hPa, respectively. Significant regional variability in water vapor trends, encompassing notable increasing and decreasing patterns, is observable in tropical and subtropical regions. At 500 and 850 hPa, strong water vapor increasing trends are noted in the equatorial Pacific Ocean and the Laccadive Sea, while decreasing trends are evident in the Indo-Pacific Ocean region and the Arabian Sea. Over land, substantial increasing trends at 850 hPa are observed in the southern United States, contrasting with decreasing trends in southern Africa and Australia. The differences between the water vapor trends of COSMIC and ERA5 are primarily negative in the tropical regions at 850 hPa. However, the water vapor increasing trends at 850 hPa estimated from COSMIC are significantly higher than the ones derived from ERA5 data for two low-height stratocumulus-cloud-rich ocean regions west of Africa and South America. These regions with notable water vapor trend differences are located in the Intertropical Convergence Zone (ITCZ) area with frequent occurrences of convection, such as deep clouds. The difference in characterizing water vapor distribution between RO and ERA5 in deep cloud regions may cause such trend differences. The assessment of spatiotemporal variability in RO-derived water vapor and reanalysis of atmospheric water vapor data helps ensure the quality of these datasets for climate studies.

Using the Commercial GNSS RO Spire Data in the Neutral Atmosphere for Climate and Weather Prediction Studies

Recently, the NOAA has included GNSS (Global Navigation Satellite System) Radio Occultation (RO) data as one of the crucial long-term observables for weather and climate applications. To include more GNSS RO data in its numerical weather prediction systems, the NOAA Commercial Weather Data Pilot program (CWDP) started to explore the commercial RO data available on the market. After two rounds of pilot studies, the CWDP decided to award the first Indefinite Delivery Indefinite Quantity (IDIQ) contract to GeoOptics and Spire Incs. in 2020. This study examines the quality of Spire RO data products for weather and climate applications. Spire RO data collected from commercial CubeSats are carefully compared with data from Formosa Satellite Mission 7–Constellation Observing System for Meteorology, Ionosphere, and Climate-2 (COSMIC-2), the fifth-generation European Centre for Medium-Range Weather Forecasts (ECMWF) atmospheric reanalysis (ERA5), and high-quality radiosonde data. The results demonstrate that, despite their generally lower Signal-Noise-Ratio (SNR), Spire RO data show a pattern of lowest penetration height similar to that of COSMIC-2. The Spire and COSMIC-2 penetration heights are between 0.6 and 0.8 km altitude over tropical oceans. Although using different GNSS RO receivers, the precision of Spire STRATOS receivers is of the same quality as those of the COSMIC-2 TriG (Global Positioning System—GPS, GALILEO, and GLObal NAvigation Satellite System—GLONASS) RO Receiver System (TGRS) receivers. Furthermore, the Spire and COSMIC-2 retrieval accuracies are quite comparable. We validate the Spire temperature and water vapor profiles by comparing them with collocated radiosonde observation (RAOB) data. Generally, over the height region between 8 km and 16.5 km, the Spire temperature profiles match those from RS41 RAOB very well, with temperature biases of <0.02 K. Over the height range from 17.8 to 26.4 km, the temperature differences are ~−0.034 K, with RS41 RAOB being warmer. We also estimate the error covariance matrix for Spire, COSMIC-2, and KOMPSAT-5. The results show that the COSMIC-2 estimated error covariance values are slightly more significant than those from Spire over the oceans at the mid-latitudes (45°N–30°N and 30°S–45°S), which may be owing to COSMIC-2 SNR being relatively lower at those latitudinal zones.

Spire RO Thermal Profiles for Climate Studies: Initial Comparisons of the Measurements from Spire, NOAA-20 ATMS, Radiosonde, and COSMIC-2

Global Navigation Satellite System (GNSS) Radio Occultation (RO) data play an essential role in improving numerical weather prediction (NWP) and monitoring climate change. The NOAA Commercial RO Purchase Program (CDP) purchased RO data provided by Spire Global Inc. To ensure the data quality from Spire Global Inc. is consistent with other RO missions, we need to quantify their accuracy and retrieval uncertainty carefully. In this work, Spire Wet Profile (wet temperature profile) data from 7 September 2021 to 31 October 2022, processed by the University Corporation for Atmospheric Research (UCAR), and COSMIC-2 (Constellation Observing System for Meteorology, Ionosphere, and Climate-2/Formosa Satellite Mission 7) data are evaluated through comparison with NOAA-20 Advanced Technology Microwave Sounder (ATMS) microwave sounder measurements and collocated RS41 radiosonde measurements. Through the Community Radiative Transfer Model (CRTM) simulation, we convert the Spire and COSMIC-2 RO retrievals to ATMS brightness temperature (BT) at sounding channels CH07 to CH14 (temperature channels), with weighting function peak heights from 8 km to 35 km, and CH19 to CH22 (water vapor channels), with weighting function peak heights ranging from 3.2 km to 6.7 km, and compare the simulations with the collocated NOAA-20 ATMS measurements over ocean. Using ATMS observations as references, Spire and COSMIC-2 BTs agree well with ATMS within 0.07 K for CH07-14 and 0.20 K for CH19-22. The trends between Spire and COSMIC-2 are consistent within 0.07 K/year over the oceans for ATMS CH07-CH13 and CH19-22, indicating that Spire/COSMIC-2 wet profiles are, in general, compatible with each other over oceans. The RO retrievals and RS41 radiosonde observation (RAOB) comparison shows that above 0.2 km altitude, RS41 RAOB matches Spire/COSMIC-2 temperature profiles well with a temperature difference of <0.13 K, and the trends between Spire and COSMIC-2 are consistent within 0.08 K/year over land, indicating that Spire/COSMIC-2 wet profiles are overall compatible with each other through RS41 RAOB measurements over land. In addition, the consistency of Spire and COSMIC-2 based on different latitude intervals, local times, and signal-to-noise ratios (SNRs) through ATMS was evaluated. The results show that the performance of Spire is comparable to COSMIC-2, even though COSMIC-2 has a higher SNR. The high quality of RO profiles from Spire is expected to improve short- and medium-range global numerical weather predictions and help construct consistent climate temperature records.

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.

COSMIC-2 soundings impacts on a RO-based NOAA microwave satellite data quality monitoring system

The Constellation Observing System for Meteorology, Ionosphere and Climate-2/Formosa Satellite Mission 7 (COSMIC-2) Global Navigation Satellite System (GNSS) Radio Occultation (RO) constellation is the follow-on to the highly successful COSMIC-1 program. The GNSS RO atmospheric soundings have historically been used to generate Community Radiative Transfer Model simulated background (B) microwave (MW) brightness temperature data needed to monitor NOAA operational MW sounding instrument observed (O) antenna temperature (Ta) product quality. This study is motivated by the need to determine the impact of COSMIC-2 RO soundings on this critical long-term monitoring capability. This study is based on individual MW sensor O-B Ta bias (∆Ta) statistics and "double-difference" inter-sensor Ta bias (δTa) statistics for two time periods. Time Period 1 (TP1—May 1, 2017, to September 30, 2019) exclusively uses COSMIC-1 and Korea Multi-Purpose Satellite-5 soundings, while Time Period 2 (TP2—October 1, 2019, to December 31, 2020) expands the analysis with COSMIC-2 soundings. The TP1 and TP2 ∆Ta statistics comparisons indicate COSMIC-2 data population augmentation and latitudinal distribution impact the MW sounder performance monitoring tool transition from TP1 to TP2. COSMIC-2 competently supports long-term MW individual sensor and inter-sensor product monitoring for MW radiometer channels with weighting functions that peak between 8 and 30 km. Individual sensor ∆Ta monitoring for other microwave channels that depend on COSMIC-2 data below 8 km and from 30 km to 60 has limitations, but these limitations are shown not to inhibit critical δTa monitoring.

A Tomographic Method for the Reconstruction of the Plasmasphere Based on COSMIC/ FORMOSAT-3 Data

A tomographic method has been developed for reconstructing the topside ionospheric and plasmaspheric electron density distribution using total electron content (TEC) measurements from global positioning system (GPS) receivers aboard the constellation observing system for meteorology, ionosphere, and climate&#x002F;Formosa Satellite Mission 3 (COSMIC&#x002F;FORMOSAT-3). Since the COSMIC&#x002F;FORMOSAT-3 constellation has an orbit altitude of about 800 km, the integral TEC measurements obtained from the topside GPS navigation data are rather small and imposed relevant challenges to obtaining stable electron density reconstructions. However, the developed method can represent the natural variability of the plasma ambient in terms of latitude, altitude, solar activity, season, and local time when analyzing electron density reconstructions during 2008&#x2013;2013. The method employs independent spatial grids for satellite rising and setting geometries and imposes a set of constraints to stabilize the solution in the presence of noise and ill-conditioned geometry. We further consider background ionosphere and plasmasphere models for electron density initialization and filling data gaps. The quality assessment using TEC and <i>in-situ</i> electron density measurements has shown that the proposed method performs better than the background model, with improvements of about 26&#x0025; in TEC and 20&#x0025; in terms of electron density. Our investigation also reveals the necessity of more accurate background electron density representations and precise TEC measurements in order to have better plasmaspheric specifications at high altitudes.

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 &lt; 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 &lt;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.

Validation and Improvement of NeQuick Topside Ionospheric Formulation Using COSMIC/FORMOSAT‐3 Data

AbstractWe examine systematic differences between topside electron density measurements and different topside model formulations including α‐Chapman, NeQuick, and an improved NeQuick (NeQuick‐corr) topside formulation, recently proposed. The global topside electron density data set used, was extracted from Radio Occultation (RO) topside electron density profiles on board Low‐Earth Orbit satellites from the COSMIC/FORMOSAT‐3 (Constellation Observing System for Meteorology, Ionosphere, and Climate and Formosa Satellite) mission. By using RO topside electron density measurements collocated with digisonde stations, we ensure that our investigation is based on two independent data sets (RO and digisondes). A subset of these profiles, with matched (within 5%) peak RO‐digisonde characteristics (foF2 and hmF2) is also exploited. This subset is exploited to extend the investigation on the basis of a higher quality validation data set. The comparison demonstrates that α‐Chapman and NeQuick‐corr underestimate, whereas NeQuick overestimates COSMIC topside electron density observations. The main outcome of this study is the significant NeQuick topside representation improvement that can be achieved near the F region peak, if the key parameter g, which controls the change of scale height with respect to altitude, is optimized to a value of 0.15 (compared to a currently adopted value of 0.125). The NeQuick‐corr topside formulation using the optimized g value of 0.15 outperforms all other topside formulations.

Initial Assessment of the COSMIC-2/FORMOSAT-7 Neutral Atmosphere Data Quality in NESDIS/STAR Using In Situ and Satellite Data

A COSMIC-1/FORMOSAT-3 (Constellation Observing System for Meteorology, Ionosphere, and Climate-1 and Formosa Satellite Mission 3) follow-on mission, COSMIC-2/FORMOSAT-7, had been successfully launched into low-inclination orbits on 25 June 2019. COSMIC-2 has a significantly increased Signal-to-Noise ratio (SNR) compared to other Radio Occultation (RO) missions. This study summarized the initial assessment of COSMIC-2 data quality conducted by the NOAA (National Oceanic and Atmospheric Administration) Center for Satellite Applications and Research (STAR). We use validated data from other RO missions to quantify the stability of COSMIC-2. In addition, we use the Vaisala RS41 radiosonde observations to assess the accuracy and uncertainty of the COSMIC-2 neutral atmospheric profiles. RS41 is currently the most accurate radiosonde observation system. The COSMIC-2 SNR ranges from 200 v/v to about 2800 v/v. To see if the high SNR COSMIC-2 signals lead to better retrieval results, we separate the COSMIC-2–RS41 comparisons into different SNR groups (i.e., 0–500 v/v group, 500–1000 v/v group, 1000–1500 v/v group, 1500–2000 v/v group, and &gt;2000 v/v group). In general, the COSMIC-2 data quality in terms of stability, precision, accuracy, and uncertainty of the accuracy is very compatible with those from COSMIC-1. Results show that the mean COSMIC-2–RS41 water vapor difference from surface to 5 km altitude for each SNR groups are equal to −1.34 g/kg (0–500 v/v), −1.17 g/kg (500–1000 v/v), −1.33 g/kg (1000–1500 v/v), −0.93 g/kg (1500–2000 v/v), and −1.52 g/kg (&gt;2000 v/v). Except for the &gt;2000 v/v group, the high SNR measurements from COSMIC-2 seem to improve the mean water vapor difference for the higher SNR group slightly (especially for the 1500–2000 v/v group) comparing with those from lower SNR groups.

Ionospheric electron density characteristics over Africa from FORMOSAT-3/COSMIC radio occultation

With the widespread availability of ground and space-based global navigation satellite system (GNSS) observables, continuous and long-term explorations of ionospheric variations have been made possible worldwide or on regional basis with improved accuracy. The Formosa Satellite Mission#3/Constellation Observing System for Meteorology, Ionosphere, and Climate (FORMOSAT-3/COSMIC) mission has a huge database of radio occultation (RO) soundings at regional and global scales with a high vertical resolution. Comparative studies between radio occultation, incoherent scatter radar and ionosonde observations indicate that COSMIC profiles agree well with ground measurements. The present paper investigates the ionospheric profiles over Africa using COSMIC data for the period from 2006 to 2017, representing almost a solar cycle year of study. The spatiotemporal variation of electron density confirms a hemispheric asymmetry among the equinoctial seasons and the solstice seasons during both low and moderate solar activity. Seasonal/winter anomaly manifestation is also clearly noticed in our observations with relatively high electron density during the winter solstice than the summer solstice. Moreover, the electron density over the region show apparent spatial and temporal variations identical to earlier ground-based ionospheric monitoring results over the African region. The outcomes from this study would strengthen the understanding of the ionospheric alterations and modelling activities in Africa, especially the areas with inadequate ground-based measuring instruments, hence, our results may complement the progress in global ionospheric modelling.

The COSMIC/FORMOSAT-3 Radio Occultation Mission after 12 Years: Accomplishments, Remaining Challenges, and Potential Impacts of COSMIC-2

AbstractLaunched in 2006, the Formosa Satellite Mission 3–Constellation Observing System for Meteorology, Ionosphere and Climate (FORMOSAT-3/COSMIC) was the first constellation of microsatellites carrying global positioning system (GPS) radio occultation (RO) receivers. Radio occultation is an active remote sensing technique that provides valuable information on the vertical variations of electron density in the ionosphere, and temperature, pressure, and water vapor in the stratosphere and troposphere. COSMIC has demonstrated the great value of RO data in ionosphere, climate, and meteorological research and operational weather forecasting. However, there are still challenges using RO data, particularly in the moist lower troposphere and upper stratosphere. A COSMIC follow-on constellation, COSMIC-2, was launched into equatorial orbit in 2019. With increased signal-to-noise ratio (SNR) from improved receivers and digital beam steering antennas, COSMIC-2 will produce at least 5,000 high-quality RO profiles daily in the tropics and subtropics. In this paper, we summarize 1) recent (since 2011 when the last review was published) contributions of COSMIC and other RO observations to weather, climate, and space weather science; 2) the remaining challenges in RO applications; and 3) potential contributions to research and operations of COSMIC-2.

Elliptical Structures of Gravity Waves Produced by Typhoon Soudelor in 2015 near Taiwan

Tropical cyclones (TCs) are complex sources of atmospheric gravity waves (GWs). In this study, the Weather Research and Forecasting Model was used to model TC Soudelor (2015) and the induced elliptical structures of GWs in the upper troposphere (UT) and lower stratosphere (LS) prior to its landfall over Taiwan. Conventional, spectral and wavelet analyses exhibit dominant GWs with horizontal and vertical wavelengths, and periods of 16–700 km, 1.5–5 km, and 1–20 h, respectively. The wave number one (WN1) wind asymmetry generated mesoscale inertia GWs with dominant horizontal wavelengths of 100–300 km, vertical wavelengths of 1.5–2.5 km (3.5 km) and westward (eastward) propagation at the rear of the TC in the UT (LS). It was also revealed to be an active source of GWs. The two warm anomalies of the TC core induced two quasi-diurnal GWs and an intermediate GW mode with a 10-h period. The time evolution of dominant periods could be indicative of changes in TC dynamics. The FormoSat-3/COSMIC (Formosa Satellite Mission-3/Constellation Observing System for Meteorology, Ionosphere, and Climate) dataset confirmed the presence of GWs with dominant vertical wavelengths of about 3.5 km in the UT and LS.

The Marine Boundary Layer Height over the Western North Pacific Based on GPS Radio Occultation, Island Soundings, and Numerical Models.

This paper estimates marine boundary layer height (MBLH) over the western North Pacific (WNP) based on Global Positioning System Radio Occultation (GPS-RO) profiles from the Formosa Satellite Mission 3 (FORMOSAT-3)/Constellation Observing System for Meteorology, Ionosphere, and Climate (COSMIC) satellites, island soundings, and numerical models. The seasonally-averaged MBLHs computed from nine years (2007–2015) of GPS-RO data are inter-compared with those obtained from sounding observations at 15 island stations and from the European Centre for Medium-Range Weather Forecasts (ECMWF) Reanalysis (ERA-Interim) and National Centers for Environmental Prediction Global Forecast System (NCEP GFS) data over the WNP from 2012 to 2015. It is found that the MBLH using nine years of GPS-RO data is smoother and more consistent with that obtained from sounding observations than is the MBLH using four years of GPS-RO data in a previous study. In winter, higher MBLHs are found around the subtropical latitudes and over oceans east of Japan, which are approximately located within the paths of the North Equatorial Current and the Kuroshio Current. The MBLH is also significantly higher in winter than in summer over the WNP. The above MBLH pattern is generally similar to those obtained from the analysis data of the ERA-Interim and NCEP GFS, but the heights are about 200 m higher. The verification with soundings suggests that the ERA-Interim has a better MBLH estimation than the NCEP GFS. Thus, the MBLH distributions obtained from both the nine-year GPS-RO and the ERA-Interim data can represent well the climatological MBLH over the WNP, but the heights should be adjusted about 30 m lower for the former and ~200 m higher for the latter. A positive correlation between the MBLH and the instability of the lower atmosphere exists over large near-shore areas of the WNP, where cold air can move over warm oceans from the land in winter, resulting in an increase in lower-atmospheric instability and providing favorable conditions for convection to yield a higher MBLH. During summer, the lower-atmospheric instability becomes smaller and the MBLH is thus lower over near-shore oceans.

Comparison of the Arctic upper-air temperatures from radiosonde and radio occultation observations

The air temperature is one of the most important parameters used for monitoring the Arctic climate change. The constellation observing system for meteorology, ionosphere, and climate and Formosa Satellite Mission 3 (COSMIC/FORMOSAT-3) radio occultation (RO) “wet” temperature product (i.e., “wetPrf”) is used to analyze the Arctic air temperature profiles at 925–200 hPa in 2007–2012. The “wet” temperatures are further compared with radiosonde (RS) observations. The results from the spatially and temporally synchronized RS and COSMIC observations show that their temperatures agree well with each other, especially at 400 hPa. Comparisons of seasonal temperatures and anomalies from the COSMIC and homogenized RS observations suggest that the limited number of COSMIC observations during the spatial matchup may be insufficient to describe the smallscale spatial structure of temperature variations. Furthermore, comparisons of the seasonal temperature anomalies from the RS and 5°×5° gridded COSMIC observations at 400 hPa during the sea ice minimum (SIM) of 2007 and 2012 are also made. The results reveal that similar Arctic temperature variation patterns can be obtained from both RS and COSMIC observations over the land area, while extra information can be further provided from the densely distributed COSMIC observations. Therefore, despite COSMIC observations being unsuitable to describe the Arctic temperatures in the lowest level, they provide a complementary data source to study the Arctic upper-air temperature variations and related climate change.

Quality aspects of the Wegener Center multi-satellite GPS radio occultation record OPSv5.6

Abstract. The demand for high-quality atmospheric data records, which are applicable in climate studies, is undisputed. Using such records requires knowledge of the quality and the specific characteristics of all contained data sources. The latest version of the Wegener Center (WEGC) multi-satellite Global Positioning System (GPS) radio occultation (RO) record, OPSv5.6, provides globally distributed upper-air satellite data of high quality, usable for climate and other high-accuracy applications. The GPS RO technique has been deployed in several satellite missions since 2001. Consistency among data from these missions is essential to create a homogeneous long-term multi-satellite climate record. In order to enable a qualified usage of the WEGC OPSv5.6 data set we performed a detailed analysis of satellite-dependent quality aspects from 2001 to 2017. We present the impact of the OPSv5.6 quality control on the processed data and reveal time-dependent and satellite-specific quality characteristics. The highest quality data are found for MetOp (Meteorological Operational satellite) and GRACE (Gravity Recovery and Climate Experiment). Data from FORMOSAT-3/COSMIC (Formosa Satellite mission-3/Constellation Observing System for Meteorology, Ionosphere, and Climate) are also of high quality. However, comparatively large day-to-day variations and satellite-dependent irregularities need to be taken into account when using these data. We validate the consistency among the various satellite missions by calculating monthly mean temperature deviations from the multi-satellite mean, including a correction for the different sampling characteristics. The results are highly consistent in the altitude range from 8 to 25 km, with mean temperature deviations less than 0.1 K. At higher altitudes the OPSv5.6 RO temperature record is increasingly influenced by the characteristics of the bending angle initialization, with the amount of impact depending on the receiver quality.

Electron density extrapolation above F2 peak by the linear Vary‐Chap model supporting new Global Navigation Satellite Systems‐LEO occultation missions

AbstractThe new radio‐occultation (RO) instrument on board the future EUMETSAT Polar System‐Second Generation (EPS‐SG) satellites, flying at a height of 820 km, is primarily focusing on neutral atmospheric profiling. It will also provide an opportunity for RO ionospheric sounding, but only below impact heights of 500 km, in order to guarantee a full data gathering of the neutral part. This will leave a gap of 320 km, which impedes the application of the direct inversion techniques to retrieve the electron density profile. To overcome this challenge, we have looked for new ways (accurate and simple) of extrapolating the electron density (also applicable to other low‐Earth orbiting, LEO, missions like CHAMP): a new Vary‐Chap Extrapolation Technique (VCET). VCET is based on the scale height behavior, linearly dependent on the altitude above hmF2. This allows extrapolating the electron density profile for impact heights above its peak height (this is the case for EPS‐SG), up to the satellite orbital height. VCET has been assessed with more than 3700 complete electron density profiles obtained in four representative scenarios of the Constellation Observing System for Meteorology, Ionosphere, and Climate (COSMIC) in the United States and the Formosa Satellite Mission 3 (FORMOSAT‐3) in Taiwan, in solar maximum and minimum conditions, and geomagnetically disturbed conditions, by applying an updated Improved Abel Transform Inversion technique to dual‐frequency GPS measurements. It is shown that VCET performs much better than other classical Chapman models, with 60% of occultations showing relative extrapolation errors below 20%, in contrast with conventional Chapman model extrapolation approaches with 10% or less of the profiles with relative error below 20%.

On the relationship between the QBO/ENSO and atmospheric temperature using COSMIC radio occultation data

In this paper, the spatial patterns and vertical structure of atmospheric temperature anomalies, in both the tropics and the extratropical latitudes, associated with the El Niño-Southern Oscillation (ENSO) and quasi-biennial oscillation (QBO) in the upper troposphere and stratosphere are investigated using global positioning system (GPS) radio occultation (RO) measurements from the Constellation Observing System for Meteorology, Ionosphere and Climate (COSMIC) Formosa Satellite Mission 3 mission from July 2006 to February 2014. We find that negative correlations between the atmospheric temperature in the tropics and ENSO are observed at 17–30km in the lower stratosphere at a lag of 1–4 months and at a lead of 1 month. Out-of-phase temperature variation is observed in the troposphere over the mid-latitude band and in-phase behaviour is observed in the lower stratosphere. Interestingly, we also find that there is a significant negative correlation at a lag of 1–3 months from 32km to 40km in the mid-latitude region of the Northern Hemisphere. The atmospheric temperature variations over mid-latitude regions in both hemispheres are closely related to the QBO. There are also two narrow zones over the subtropical jet zone where the QBO signals are strong in both hemispheres, approximately parallel to the equator. Finally, we develop a new robust index to describe the strength of the ENSO and QBO signal.

An Impact Study of GPS Radio Occultation Observations on Frontal Rainfall Prediction with a Local Bending Angle Operator

Abstract The impact of global positioning system (GPS) radio occultation (RO) soundings on the prediction of severe mei-yu frontal rainfall near Taiwan in June 2012 was investigated in this study using a developed local bending angle (LBA) operator. Two operators for local refractivity (REF) and nonlocal refractivity [excess phase (EPH)] were also used for comparisons. The devised LBA simplifies the calculation of the Abel transform in inverting model local refractivity without a loss of accuracy. These operators have been implemented into the three-dimensional variational data assimilation system of the Weather Research and Forecasting (WRF) Model to assimilate GPS RO soundings available from the Formosa Satellite Mission 3/Constellation Observing Systems for Meteorology, Ionosphere and Climate (FORMOSAT-3/COSMIC). The RO data are found to be beneficial to the WRF forecast of local severe rainfall in Taiwan. Characteristics of assimilation performance and innovation for the three operators are discussed. Both of the local operators performing assimilation at observation levels appear to produce mostly larger positive moisture increments than do the current nonlocal operators performing assimilation on the mean height of each model vertical level. As the information of the initial increments has propagated farther south with the frontal flow, the simulation for LBA shows better prediction of rainfall peaks in Taiwan on the second day than both REF and EPH, with a maximum improvement of about 25%. The positive impact of the RO data results partially from several RO observations near Mongolia and north China. This study provides an intercomparison among the three RO operators, and shows the feasibility of regional assimilation with LBA.

Characteristics of the seasonal variation of the global tropopause revealed by COSMIC/GPS data

Using the Global Navigation Satellite System (GNSS) radio occultation observations from Formosa Satellite mission-3/Constellation Observing System for Meteorology, Ionosphere, and Climate (FORMOSAT-3/COSMIC) from 2007 to 2012, the climatological characteristics of the global tropopause was studied, with the following features identified. The overall results generally agree with previous studies. The tropopause has an obvious zonal structure, with more zonal characteristics in the Southern Hemisphere than the Northern Hemisphere. The vertical shape of the tropopause is sharp in the tropics and broad in the sub-tropical latitudes, with the sharpest latitudinal gradient in the mid-latitudes of both hemispheres. The global tropopause exists in a large range between 8km and 17km (or between 100hPa and 340hPa). The highest tropopause is over the South Asian monsoon regions for the entire year. The spatial structure of the tropopause in the polar region is of concentric structure, with an altitude between 7.5km and 10km. It is more symmetric in the Antarctic than the Arctic. Differing from other places, the height of the tropopause in the Antarctic is higher in winter as opposed to summer. The tropopause has distinct seasonal variability, especially in polar regions.