GPS Radio Occultation Research Articles

The Boundary Layer Height Obtaining by New Spline Numerical Regularization Method Based on Gps Radio Occultation Data: A Case Study of the Qinghai-Tibet Plateau

The Boundary Layer Height Obtaining by New Spline Numerical Regularization Method Based on Gps Radio Occultation Data: A Case Study of the Qinghai-Tibet Plateau

Global 3D Features of Error Variances of GPS Radio Occultation and Radiosonde Observations

Global Positioning System (GPS) radio occultation (RO) and radiosonde (RS) observations are two major types of observations assimilated in numerical weather prediction (NWP) systems. Observation error variances are required input that determines the weightings given to observations in data assimilation. This study estimates the error variances of global GPS RO refractivity and bending angle and RS temperature and humidity observations at 521 selected RS stations using the three-cornered hat method with additional ERA-Interim reanalysis and Global Forecast System forecast data available from 1 January 2016 to 31 August 2019. The global distributions, of both RO and RS observation error variances, are analyzed in terms of vertical and latitudinal variations. Error variances of RO refractivity and bending angle and RS specific humidity in the lower troposphere, such as at 850 hPa (3.5 km impact height for the bending angle), all increase with decreasing latitude. The error variances of RO refractivity and bending angle and RS specific humidity can reach about 30 N-unit2, 3 × 10−6 rad2, and 2 (g kg−1)2, respectively. There is also a good symmetry of the error variances of both RO refractivity and bending angle with respect to the equator between the Northern and Southern Hemispheres at all vertical levels. In this study, we provide the mean error variances of refractivity and bending angle in every 5°-latitude band between the equator and 60°N, as well as every interval of 10 hPa pressure or 0.2 km impact height. The RS temperature error variance distribution differs from those of refractivity, bending angle, and humidity, which, at low latitudes, are smaller (less than 1 K2) than those in the midlatitudes (more than 3 K2). In the midlatitudes, the RS temperature error variances in North America are larger than those in East Asia and Europe, which may arise from different radiosonde types among the above three regions.

E Layer Dominated Ionosphere Occurrences as a Function of Geophysical and Space Weather Conditions

At some locations, especially in the auroral regions, the ionization of the E layer can dominate over that of the F2 layer, which is called the E layer dominated ionosphere (ELDI). In the present work we investigate the spatiotemporal variation of the ELDI depending on the season, solar activity, geomagnetic activity, interplanetary magnetic field, convection electric field, and solar wind energy. We specify each distribution of ELDI events by the values of four parameters. In this regard, we compute the height, full width at half maximum, and position of a Gaussian function relative to a precomputed reference ellipse as parameters to describe the spatial distribution of ELDI events in geocentric latitude/longitude coordinates. To study the temporal variation of the ELDI events, we estimate the weighted mean local time of the distribution as the fourth parameter. The database used for our investigations contains more than 3.5 million vertical electron density profiles derived from ionospheric GPS radio occultation observations on board the COSMIC/FORMOSAT-3 (Constellation Observing System for Meteorology, Ionosphere, and Climate/Formosa Satellite Mission 3) mission, covering a period of almost 13 years. The analysis of observations representing changing geophysical conditions results in clear trends for all ELDI parameters. In this context, the mean local time varies mostly between 01:00 and 02:00 local time, while the probability of ELDI occurrence is increased in local winter and in the case of low solar activity. Likewise, an increase in the solar wind parameters increases the number of ELDI events and leads to an equatorward shift of their position. The relationships found in our investigations can serve as a basis for future modeling studies addressing ELDI occurrences as a function of selected geophysical quantities.

Estimating GPS radio occultation observation error standard deviations over China using the three‐cornered hat method

AbstractObservation error variances are required input parameters for assimilating GPS radio occultation (RO) data in numerical weather prediction systems. A three‐cornered hat method is applied to estimate RO observation error standard deviation (SD) of both refractivity and bending angle. RO data from the Global Navigation Satellite System Receiver for Atmospheric Sounding on board the MetOp‐A/B satellites from 1 January 2016 to 31 August 2019, along with radiosonde observations, the European Centre for Medium‐Range Weather Forecasts Interim reanalysis, and the Global Forecast System forecast analysis, are employed in this study. Taking the Guam radiosonde station as an example, we first introduce a detailed procedure of estimating GPS RO fractional error SDs, which shows the significant impacts of outliers on the resulting error SDs, especially in the lower and middle troposphere. Fractional bending angle error SDs near 48 radiosonde stations in China are then estimated for the same time period of more than 3 years. The RO fractional error SDs near the 48 radiosonde stations have a significant latitudinal dependence, larger at lower than higher latitudes. This characteristic feature of the RO bending angle errors agrees with that of the bending angle uncertainty information provided by the Constellation Observing System for Meteorology, Ionosphere, and Climate Data Analysis and Archival Center.

On the Seeding of Periodic Equatorial Plasma Bubbles by Gravity Waves Associated With Tropical Cyclone: A Case Study

AbstractUnderstanding the day‐to‐day variability of equatorial plasma bubble (EPB) irregularities has been a major challenge for many decades that makes the deterministic prediction of EPBs remain elusive. In the present study, we report a case of intense and periodic EPBs observed during 8 and 9 April 2013 by the 47‐MHz Equatorial Atmosphere Radar at Kototabang, Indonesia. The periodic EPBs separated by about 200–250 km were initiated before sunset. The presunset onset and development of these periodic EPBs were discussed in light of the gravity waves (GWs) excited in connection with the deep convection due to the tropical cyclone (TC) Victoria. The outgoing long‐wave radiation measurement by very high‐resolution radiometer (VHRR) onboard Indian meteorological satellite Kalpana‐1 shows the occurrence of deep convective activity during these days. The presence of upward propagating GWs from the deep convective region associated with TC Victoria was confirmed using the GPS radio occultation (GPS‐RO) observations in the troposphere and stratosphere. The GW signatures at ionospheric altitudes were also observed from the Ionosonde observations over magnetic equator, and medium‐scale (~300 km) GWs were observed from the GPS‐TEC (Total Electron Content) data near to the magnetic equator and cyclone center. From the GW parameters observed from GPS‐TEC and GPS‐RO, we surmise that the secondary GWs generated by the dissipation of primary GWs associated with TC Victoria could have served as a seeding source on the generation of periodic EPBs during these two consecutive days.

Marine Boundary Layer Height Obtained by New Numerical Regularization Method Based on GPS Radio Occultation Data.

The boundary layer height (BLH) determines the interface between the lower and the free atmosphere, and it is a key variable in numerical simulations and aerosol and environmental pollution studies. This article proposes a novel method in conjunction with numerical regularization to analyze the climate characteristics of the marine boundary layer height (MBLH) using 2007–2011 GPS-RO data from the COSMIC mission. The MBLH corresponds to the smallest gradient, which is calculated using the numerical regulation method where the regulation parameters are determined by the double-parameter model function method. The results reveal the relationship between the MBLH and ocean currents for the first time. A low MBLH is associated with cold seasons and seas where cold ocean currents prevail whereas a high MBLH is related to warm seasons and seas where warm currents prevail. This correlation was validated by comparing the obtained results with different occultation data including atmprf and echprf, which also showed that atmprf is more sensitive to convective cloud top capture. To test the credibility of the results, the standard deviation was used to express the MBLH confidence level. The results show that the standard deviation of the MBLH was highest in low latitudes and lowest in the middle and high latitudes. Furthermore, we analyzed the trends in interannual MBLH variations, which display significant seasonal variations and spatial distributions that correspond with the current and subsolar point. Finally, we conducted a case study in the South China Sea, and identified a distinctive seasonal change and interannual decline in MBLH.

Climatology of Planetary Boundary Layer Height-Controlling Meteorological Parameters Over the Korean Peninsula

Planetary boundary layer (PBL) height plays a significant role in climate modeling, weather forecasting, air quality prediction, and pollution transport processes. This study examined the climatology of PBL-associated meteorological parameters over the Korean peninsula and surrounding sea using data from the ERA5 dataset produced by the European Centre for Medium-range Weather Forecasts (ECMWF). The data covered the period from 2008 to 2017. The bulk Richardson number methodology was used to determine the PBL height (PBLH). The PBLH obtained from the ERA5 data agreed well with that derived from sounding and Global Positioning System Radio Occultation datasets. Significant diurnal and seasonal variability in PBLH was observed. The PBLH increases from morning to late afternoon, decreases in the evening, and is lowest at night. It is high in the summer, lower in spring and autumn, and lowest in winter. The variability of the PBLH with respect to temperature, relative humidity, surface pressure, wind speed, lower tropospheric stability, soil moisture, and surface fluxes was also examined. The growth of the PBLH was high in the spring and in southern regions due to the low soil moisture content of the surface. A high PBLH pattern is evident in high-elevation regions. Increasing trends of the surface temperature and accordingly PBLH were observed from 2008 to 2017.

A Statistical Analysis of Global Ionospheric E‐Layer Scintillation During 2007–2014

AbstractThis paper presents the global three‐dimensional morphology and seasonal variation of the scintillation index at E‐layer altitudes varying from 70 to 140 km. An attempt is conducted to analyze the scintillation distribution in different latitude bins and to explore solar and geomagnetic activity‐dependent characteristic of E‐layer scintillations. The scintillation index is derived from the signal‐to‐noise ratio intensity fluctuations of L1 channels in GPS radio occultation (RO) using the Constellation Observing System for Meteorology, Ionosphere, and Climate (COSMIC) during 2007–2014. The results show that the summer scintillation index has the greatest values, followed by the winter, autumn, and spring months. The summer scintillation is in close connection with the solar radiation intensity and sunshine duration, but it exhibits seasonal asymmetry whereby the summer scintillation in the northern hemisphere is more intense than in the southern hemisphere. Meanwhile, the scintillation in the autumn hemisphere is slightly more intense than in the spring hemisphere in approximately midlatitude regions during most daytimes. The conspicuous scintillation activity in the middle and low latitudes starts post‐sunrise and often persists until post‐midnight. Some observed results show that the greatest scintillation occurrence rate of the latitude profile occurs at approximately sunset and at an altitude range of 100–120 km in the middle‐low latitudes. Additionally, the enhanced solar activity possibly can result in the slight decrease of scintillation occurrence, and geomagnetic activity has slight effect on E‐layer scintillation in most cases.

Observations of Upper-Tropospheric Temperature Inversions in the Indian Monsoon and Their Links to Convectively Forced Quasi-Stationary Kelvin Waves

Abstract High-vertical-resolution temperature measurements from GPS radio occultation data show frequent upper-tropospheric inversions over the equatorial Indian Ocean during the summer monsoon season. Each year, around 30% of profiles in this region have temperature inversions near 15 km during the monsoon season, peaking during July–September. This work describes the space–time behavior of these inversions, and their links to transient deep convection. The Indian Ocean inversions occur episodically several times each summer, with a time scale of 1–2 weeks, and are quasi stationary or slowly eastward moving. Strong inversions are characterized by cold anomalies in the upper-troposphere (12–15 km), warm anomalies in the tropopause layer (16–18 km), and strong zonal wind anomalies that are coherent with temperature anomalies. Temperature and wind anomalies are centered over the equator and show a characteristic eastward phase tilt with height with a vertical wavelength near 5 km, consistent with a Kelvin wave structure. Composites of outgoing longwave radiation (OLR) show that strong inversions are linked to enhanced deep convection over the equatorial Indian Ocean, preceding the inversions by ~2–6 days. These characteristics suggest that the inversions are linked to convectively forced Kelvin waves, which are Doppler shifted by the easterly monsoonal winds such that they remain quasi stationary in the equatorial Indian Ocean. These large-scale waves influence circulation on the equatorial side of the Indian monsoon anticyclone; they may provide a positive feedback to the underlying convection, and are possibly linked with regions of shear-induced turbulence.

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.

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

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

Evaluation of the 15-year ROM SAF monthly mean GPS radio occultation climate data record

Abstract. We here present results from an evaluation of the Radio Occultation Meteorology Satellite Application Facility (ROM SAF) gridded monthly mean climate data record (CDR v1.0), based on Global Positioning System (GPS) radio occultation (RO) data from the CHAMP (CHAllenging Minisatellite Payload), GRACE (Gravity Recovery and Climate Experiment), COSMIC (Constellation Observing System for Meteorology, Ionosphere, and Climate), and Metop satellite missions. Systematic differences between RO missions, as well as differences of RO data relative to ERA-Interim reanalysis data, are quantified. The methods used to generate gridded monthly mean data are described, and the correction of monthly mean RO climatologies for sampling errors, which is essential for combining data from RO missions with different sampling characteristics, is evaluated. We find good overall agreement between the ROM SAF gridded monthly mean CDR and the ERA-Interim reanalysis, particularly in the 8–30 km height interval. Here, the differences largely reflect time-varying biases in ERA-Interim, suggesting that the RO data record has a better long-term stability than ERA-Interim. Above 30–40 km altitude, the differences are larger, particularly for the pre-COSMIC era. In the 8–30 km altitude region, the observational data record exhibits a high degree of internal consistency between the RO satellite missions, allowing us to combine data into multi-mission records. For global mean bending angle, the consistency is better than 0.04 %, for refractivity it is better than 0.05 %, and for global mean dry temperature the consistency is better than 0.15 K in this height interval. At altitudes up to 40 km, these numbers increase to 0.08 %, 0.11 %, and 0.50 K, respectively. The numbers can be up to a factor of 2 larger for certain latitude bands compared to global means. Below about 8 km, the RO mission differences are larger, reducing the possibilities to generate multi-mission data records. We also find that the residual sampling errors are about one-third of the original and that they include a component most likely related to diurnal or semi-diurnal cycles.

Improvement of Accuracy of Global Numerical Weather Prediction Using Refined Error Covariance Matrices

Abstract In data assimilation for NWP, accurate estimation of error covariance matrices (ECMs) and their use are essential to improve NWP accuracy. The objective of this study is to estimate ECMs of all observations and background variables using sampling statistics, and improve global NWP accuracy by using them. This study presents the first results of such all ECM refinement. ECM diagnostics combining multiple methods, and analysis and forecast cycle experiments were performed on the JMA global NWP system, where diagonal components of all ECMs and off-diagonal components of radiance observations are refined. The ECM diagnostic results are as follows: 1) the diagnosed error standard deviations (SDs) are generally much smaller than those of the JMA operational system (CNTL); 2) interchannel correlations in humidity-sensitive radiance errors are much larger than 0.2; and 3) horizontal correlation distances of AMSU-A are ~50 km, excluding channel 4. The experimental results include the following: 1) the diagnosed ECMs generally improve forecast accuracy over CNTL even without additional tunings; 2) the supplemental tuning parameter, which is the deflation factor (0.6 in SD) applied for the estimated ECMs of nonsatellite conventional data and GPS radio occultation data, statistically significantly improves forecast accuracy; 3) this value 0.6 is set as the same value as the ratio of the estimated background error SD to that in CNTL; 4) high-density assimilation (10 times) of AMSU-A is better than CNTL, not better than that with 5 times; and 5) ECMs estimated using boreal summer data can improve forecast accuracy in winter, which indicates their robustness.

Ducting and Biases of GPS Radio Occultation Bending Angle and Refractivity in the Moist Lower Troposphere

AbstractRadio occultation (RO) can provide high-vertical-resolution thermodynamic soundings of the planetary boundary layer (PBL). However, sharp moisture gradients and strong temperature inversion lead to large gradients in refractivity N and often cause ducting. Ducting results in systematically negative RO N biases resulting from a nonunique Abel inversion problem. Using 8 years (2006–13) of Constellation Observing System for Meteorology, Ionosphere and Climate (COSMIC) RO soundings and collocated European Centre for Medium-Range Weather Forecasts interim reanalysis (ERA-I) data, we confirm that the large lower-tropospheric negative N biases are mainly located in the subtropical eastern oceans and we quantify the contribution of ducting for the first time. The ducting-contributed N biases in the northeast Pacific Ocean (160°–110°W; 15°–45°N) are isolated from other sources of N biases using a two-step geometric-optics simulation. Negative bending angle biases in this region are also observed in COSMIC RO soundings. Both the negative refractivity and bending angle biases in COSMIC soundings mainly lie below ~2 km. Such bending angle biases introduce N biases that are in addition to those caused by ducting. Following the increasing PBL height from the southern California coast westward to Hawaii, centers of maxima bending angles and N biases tilt southwestward. In areas where ducting conditions prevail, ducting is the major cause of the RO N biases. Ducting-induced N biases with reference to ERA-I compose over 70% of the total negative N biases near the southern California coast, where strongest ducting conditions prevail, and decrease southwestward to less than 20% near Hawaii.

Consistency and structural uncertainty of multi-mission GPS radio occultation records

Abstract. Atmospheric climate monitoring requires observations of high quality that conform to the criteria of the Global Climate Observing System (GCOS). Radio occultation (RO) data based on Global Positioning System (GPS) signals are available since 2001 from several satellite missions with global coverage, high accuracy, and high vertical resolution in the troposphere and lower stratosphere. We assess the consistency and long-term stability of multi-satellite RO observations for use as climate data records. As a measure of long-term stability, we quantify the structural uncertainty of RO data products arising from different processing schemes. We analyze atmospheric variables from bending angle to temperature for four RO missions, CHAMP, Formosat-3/COSMIC, GRACE, and Metop, provided by five data centers. The comparisons are based on profile-to-profile differences aggregated to monthly medians. Structural uncertainty in trends is found to be lowest from 8 to 25 km of altitude globally for all inspected RO variables and missions. For temperature, it is < 0.05 K per decade in the global mean and < 0.1 K per decade at all latitudes. Above 25 km, the uncertainty increases for CHAMP, while data from the other missions – based on advanced receivers – are usable to higher altitudes for climate trend studies: dry temperature to 35 km, refractivity to 40 km, and bending angle to 50 km. Larger differences in RO data at high altitudes and latitudes are mainly due to different implementation choices in the retrievals. The intercomparison helped to further enhance the maturity of the RO record and confirms the climate quality of multi-satellite RO observations towards establishing a GCOS climate data record.

Occurrence climatology of equatorial plasma bubbles derived using FormoSat-3 ∕ COSMIC GPS radio occultation data

Abstract. The Global Positioning System – Radio Occultation (GPS-RO) observations from FormoSat-3 ∕ COSMIC are used to comprehend the global distribution of equatorial plasma bubbles which are characterized by depletion regions of plasma in the F region of the ionosphere. Plasma bubbles that cause intense scintillation of the radio signals are identified based on the S4 index derived from the 1 Hz raw signal-to-noise ratio measurements between 2007 and 2017. The analyses revealed that bubbles influenced by background plasma density occurred along the geomagnetic equator and had an occurrence peak around the dip equator during high solar activity. The peak shifted between the African and American sectors, depending on different solar conditions. Plasma bubbles usually developed around 19:00 local time (LT), with maximum occurrence around 21:00 LT during solar maximum and ∼22:00 LT during solar minimum. The occurrence of bubbles showed a strong dependence on longitudes, seasons, and solar cycle with the peak occurrence rate in the African sector around the March equinox during high solar activity, which is consistent with previous studies. The GPS-RO technique allows an extended analysis of the altitudinal distribution of global equatorial plasma bubbles obtained from high vertical resolution profiles, thus making it a convenient tool which could be further used with other techniques to provide a comprehensive view of such ionospheric irregularities.

Seasonal and Annual Changes of the Regional Tropical Belt in GPS-RO Measurements and Reanalysis Datasets

AbstractThe width of the tropical belt has been analyzed with a variety of metrics, often based on zonal-mean data from reanalyses. However, constraining the global and regional tropical width requires both a global spatial-resolving observational dataset and an appropriate metric to take advantage of such data. The tropical tropopause break is arguably such a metric. This study aims to evaluate the performance of different reanalyses and metrics with a focus on depicting regional tropical belt width. We choose four distinct tropopause-break metrics derived from global positioning system radio occultation (GPS-RO) satellite data and four modern reanalyses (ERA-Interim, MERRA-2, JRA-55, and CFSR). We show that reanalyses generally reproduce the regional tropical tropopause break to within 10° of that in GPS-RO data—but that the tropical width is somewhat sensitive (within 4°) to how data are averaged zonally, moderately sensitive (within 10°) to the dataset resolution, and more sensitive (20° over the Northern Hemisphere Atlantic Ocean during June–August) to the choice of metric. Reanalyses capture the poleward displacement of the tropical tropopause break over land and equatorward displacement over ocean during summertime, and the reverse during the wintertime. Reanalysis-based tropopause breaks are also generally well correlated with those from GPS-RO, although CFSR reproduces 14-yr trends much more closely than others (including ERA-Interim). However, it is hard to say which dataset is the best match of GPS-RO. We further find that the tropical tropopause break is representative of the subtropical jet latitude and the Northern Hemisphere edge of the Hadley circulation in terms of year-to-year variations.

Characterizing Extratropical Tropopause Bimodality and its Relationship to the Occurrence of Double Tropopauses Using COSMIC GPS Radio Occultation Observations

Lapse rate tropopause (LRT) heights in the extratropics have been shown to display a bimodal distribution, with one modal maxima above 15 km (typical of the tropical tropopause) and the other below 13 km (typical of the extratropical tropopause). The climatology of the tropopause is studied by characterizing tropopause bimodality and how it relates to the occurrence of double tropopauses (DTs). LRT heights are derived from Constellation Observing System for Meteorology, Ionosphere, and Climate (COSMIC) GPS Radio Occultation temperature profiles from 2006 to 2017. Tropopause bimodality occurs most frequently within a subtropical band (20°–40°) in both hemispheres. A distinct seasonality is observed as bimodality occurs most frequently in winter except for another local maximum along the northern edge of the Asian summer monsoon. The regions with a bimodal height distribution nearly overlap the regions that experience a high frequency of DTs. DTs occur most frequently in winter (50%–70% of the time) along the poleward edge of the bimodal band, and most LRT heights are within the extratropical mode (>80%), whereas DT occurrence decreases quickly toward the equatorward edge (<20%) along with fewer LRT heights in the extratropical mode (<50%). These results indicate that LRT height bimodality occurs along the equatorward edge due to the occurrences of double tropopauses, while the poleward edge is due to single tropopause profiles that are more tropical in nature.

Homogeneity of the Temperature Data Series from ERA5 and MERRA2 and Temperature Trends

The stratosphere and its dynamics are a very important part of atmospheric circulation. We need to analyze its climatology, as well as long-term trends. A long-term trend study needs homogenous datasets without significant artificial discontinuities. The analysis is based on the two newest released reanalyses, Modern Era-Retrospective Analysis (MERRA2) and European Center for Medium-Range Weather Forecast Reanalysis (ERA5). The aim of this study is to detect discontinuities in the temperature time series from the above reanalyses with the help of the Pettitt homogeneity test for pressure layers above 500 hPa up to 1 hPa in January and February, and show a comparison of temperature trends from the studied reanalyses and GPS radio occultation (GPS RO). We search for individual grid points where these discontinuities occur, and also for the years when they occur (geographical and temporal distribution). As expected, the study confirms better results for the Northern Hemisphere due to the denser data coverage. A high number of grid points with jumps on the Southern Hemisphere is found, especially at higher pressure levels (from 50 hPa). The spatial and vertical distribution of discontinuities is also presented. The vertical distribution reveals the reduction of the number of jumps around 10 hPa, especially for ERA5 reanalysis. The results show that ERA5 has significantly less jumps than MERRA2. We also study temperature trends from reanalyses and GPS RO and our analysis shows that the agreement between the reanalyses and observations are very good for the period 2006–2018.

Comparison of Advanced Technology Microwave Sounder Biases Estimated Using Radio Occultation and Hurricane Florence (2018) Captured by NOAA-20 and S-NPP

The second Advanced Technology Microwave Sounder (ATMS) was onboard the National Oceanic and Atmospheric Administration (NOAA)-20 satellite when launched on 18 November 2017. Using nearly six months of the earliest NOAA-20 observations, the biases of the ATMS instrument were compared between NOAA-20 and the Suomi National Polar-Orbiting Partnership (S-NPP) satellite. The biases of ATMS channels 8 to 13 were estimated from the differences between antenna temperature observations and model simulations generated from Meteorological Operational (MetOp)-A and MetOp-B satellites’ Global Positioning System (GPS) radio occultation (RO) temperature and water vapor profiles. It was found that the ATMS onboard the NOAA-20 satellite has generally larger cold biases in the brightness temperature measurements at channels 8 to 13 and small standard deviations. The observations from ATMS on both S-NPP and NOAA-20 are shown to demonstrate an ability to capture a less than 1-h temporal evolution of Hurricane Florence (2018) due to the fact that the S-NPP orbits closely follow those of NOAA-20.