High vertical resolution water vapour profiles in the upper troposphere and lower stratosphere retrieved from MAESTRO solar occultation spectra

An algorithm has been developed that retrieves water vapour profiles in the upper troposphere and lower stratosphere from optical depth spectra obtained by the Measurements of Aerosol Extinction in the Stratosphere and Troposphere Retrieved by Occultation (MAESTRO) instrument onboard the SCISAT satellite as part of the Atmospheric Chemistry Experiment (ACE) mission. The retrieval relies on ro-vibrational absorption of solar radiation by water vapour in the 926–970 nm range. During the iterative inversion process, the optical depth spectra are simulated at the spectral resolution and sampling frequency of MAESTRO using the correlated-k approximation. The Chahine inversion updates the water vapour volume mixing ratio (VMR), adjusting all retrieval layers simultaneously, to match the observed differential optical depth due to absorption by water vapour and ozone at each tangent height. This approach accounts for significant line saturation effects. Profiles are typically obtained from ∼22 km down to the cloud tops or to 5 km, with relative precision as small as 3% in the troposphere. In the lower stratosphere, the precision on water vapour VMR is ∼1.3 μmol/mol in an individual retrieval layer (∼1 km thick). The spectral capability of MAESTRO allows for the clear separation of extinction due to water vapour and aerosol, and for the fitting quality to be quantified and used to determine an altitude-dependent convergence criterion for the retrieval. In the middle troposphere, interhemispheric differences in water vapour VMR are driven by oceanic evaporation whereas in the upper troposphere, deep convection dominates and a strong seasonal cycle is observed at high latitudes.     

Validation of the operational IASI level 2 processor using AIRS and ECMWF data

The algorithms being implemented in EUMETSAT’s IASI Level 2 Product Processing Facility are validated with real case situations using AIRS data and comparing the retrieved atmospheric states with ECMWF analyses. The tests have been performed for clear-sky ocean scenes during daytime.

The Empirical Orthogonal Function (EOF) retrievals show very good performance, with retrieved atmospheric states standard deviations between 1 and 2 K in temperature and 10% and 20% in relative humidity when compared with ECMWF analysis in the troposphere. The EOF retrievals show relatively smooth profiles.

Results from an iterative retrieval show a standard deviation between 2 and 3 K in temperature and 10% and 30% in relative humidity when compared with ECMWF analyses in the troposphere. They tend to show meteorologically reasonable discontinuities in both temperature and relative humidity. This seems to be the reason why they do not compare as well with ECMWF analyses as the EOF retrievals do. Whether they are closer to reality or not will have to be tested with co-located radiosondes or similar more accurate data, which generally do not exhibit such smooth vertical profiles as ECMWF analyses do.     

Spaceborne lidar measurement accuracy: Simulation of aerosol,cloud, molecular density,and temperature retrievals

Spaceborne lidar measurements and retrievals are simulated using realistic errors in signal, conventional density information, atmospheric transmission, and lidar calibration. We find that by day, independent analysis of returns at wavelengths of 0.53 and 1.06 μm yields vertical profiles (0.1- to 1-km resolution) of tenuous clouds and boundary-layer, Saharan, and strong volcanic stratospheric aerosols to accuracies of 30% or better, provided particulate optical depth does not exceed ?0.3. By night all these constituents are retrieved, plus noctilucent clouds, mesospheric aerosols, and upper tropospheric/nonvolcanic stratospheric (UT/NVS) areosols. Molecular-density uncertainties are a dominant source of error for UT/NVS retrievals.To reduce these errors and also to provide density and temperature profiles, we developed a procedure that combines returns at 0.35 and 1.06 μm. This technique significantly improves UT/NVS aerosol retrieval accuracy and also yields useful density and temperature profiles there. Strong particulate contamination limits the technique to the cloud-free upper troposphere and above.     

Gravity wave propagation studies using the Indian MST radar observations

The Indian MST radar facility at Gadanki (13.5°N, 79.2°E) has been utilised to study the propagation of gravity waves from the troposphere/lower stratosphere to the mesosphere and their interaction with the radar backscattered signal variations. The main objective is to correlate vertically propagating gravity waves derived from the tropospheric velocity fields with the dynamics of mesospheric scattering centres. The tropospheric wind velocities and signal strengths over the entire height range have been subjected to power spectral and wavelet analysis to determine the predominant wave periods/amplitudes and the coupling between the lower atmosphere and mesosphere. Results show that (a) the gravity waves are clearly detectable near tropopause heights, (b) while relatively higher period gravity waves (20–50 min) interact with mesospheric scattering centres, the lower period waves (<20 min) are absorbed in the troposphere itself, (c) the mesospheric scattering layers are affected by gravity waves of complementary periods.     

CRISTA-NF measurements of water vapor during the SCOUT-O3 Tropical Aircraft Campaign

The new remote sensing experiment CRISTA-NF (Cryogenic Infrared Spectrometers and Telescopes for the Atmosphere – New Frontiers) successfully participated in the SCOUT-O3 Tropical Aircraft Campaign in November and December 2005. CRISTA-NF operated aboard the high-altitude research aircraft M-55 Geophysica. Mid-infrared spectra (4–15 μm) were measured in the limb sounding geometry with high spatial resolution (250 m vertical sampling, 5–15 km along track sampling). Measurements were carried out during transfer flights between Oberpfaffenhofen, Germany, and Darwin, Australia, as well as during several local flights near Darwin. Water vapor volume mixing ratios in the upper troposphere and lower stratosphere were derived from the CRISTA-NF radiance measurements by utilizing a rapid radiative transfer forward model and the optimal estimation retrieval approach. CRISTA-NF water vapor measurements below the hygropause have a total retrieval error of 15–40% (i.e. root mean square of accuracy and precision). The systematic terms are dominating in the retrieval error budget. The contributions of a priori information to the retrieval results are less than 5–10%. The vertical resolution of the observations is about 250–500 m when permitted by instrument sampling. In this paper we present first results for three transfer flights of the campaign. Being generally in good agreement with corresponding ECMWF operational analyzes, the CRISTA-NF measurements show significantly higher variability and local structures in the upper tropospheric water vapor distributions.     

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 8 km and 17 km (or between 100 hPa and 340 hPa). 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.5 km and 10 km. 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.     

Decadal variability of tropical Pacific temperature in relation to solar cycles

We use the 8-year long satellite temperature data (2002–2010) from Atmospheric InfraRed Sounder (AIRS) and Atmospheric Microwave Sounding Unit (AMSU) on the Aqua satellite to identify temperature trends in the troposphere and low stratosphere over the Niño 3.4 region of the Tropical Pacific Ocean in the most recent 11-year solar cycle. Employing more extended sea surface temperature (SST) data for five solar cycles (1950–2009) in this region we show that the satellite trends reflect a typical decrease of the sea surface temperature (SST) in the Niño 3.4 region in the declining phase of the solar cycle. The magnitude of the SST decrease depends on the solar cycle and ranges between 0.07 K/yr and 0.27 K/yr for the last five solar cycles.     

Sodium lidar measurements of mesopause region temperatures at 23° S

A sodium lidar, capable of measuring temperature in the 80–100 km region, has been in operation at São José dos Campos (23° S, 46 W) since March 2007. Good quality data have been obtained for late autumn, winter and spring, but weather conditions make it extremely difficult to make measurements from mid-November to mid- February. We find the temperature structure to be strongly modulated by tides and gravity waves, but average profiles typically show a primary mesopause height close to 100 km with temperatures around 180 K, and a tendency for a secondary minimum of about 185 K to occur close to 90 km. Vertical temperature gradients greater than 50 K/km are sometimes seen even on profiles averaged over several hours. The strongest gradients are always positive and are frequently associated with strong gradients in sodium concentration. On the other hand, we frequently see rapid changes in the temperature profile, suggesting that models and non-local temperature measurements, as made by satellite radiometers, for example, are of little use in applications such as the analysis of gravity wave propagation seen in airglow images.     

High resolution limb observations of clouds by the CRISTA-NF experiment during the SCOUT-O3 tropical aircraft campaign

The Cryogenic Infrared Spectrometers and Telescopes for the Atmosphere – New Frontiers (CRISTA-NF) experiment on board the Russian research aircraft Geophysica measures limb emission spectra with an unprecedented vertical and horizontal resolution in the 4–15 μm wavelength region. The IR spectra measured during the SCOUT-O3 Tropical Aircraft Campaign have been analysed with respect of cloud occurrence, cloud vertical and horizontal extent, cloud spatial structures and their utilisation for trace gas retrievals. In addition indicators for ice water content and optical thickness of the clouds have been adopted. These new kinds of measurements in the upper troposphere/lower stratosphere region are especially valuable for the design and development of future space borne high resolution limb sounders.     

Vertical wavenumber spectra of atmospheric ozone measured from ozonesonde observations

Vertical profiles of ozone have been measured at balloon altitudes. Our purpose is to examine the character of vertical wavenumber spectra of ozone fluctuations, to assess the possible roles of gravity wave field in ozone fluctuations, and to determine dominant vertical wavelengths of ozone spectra. Vertical wavenumber spectra of 12 ozone fluctuations obtained during June–August 2003 are presented. Results indicate that mean spectral slopes in the wavenumber range from 4.69 × 10⁻⁴ to 2.50 × 10⁻³ cyc/m are about −2.91 in the troposphere and −2.87 in the lower stratosphere, which is close to the slope of −3 predicted by current gravity wave saturation models. The consistency of the observed spectral slopes with the value of −3 predicted by current gravity wave saturation models suggests that the observed ozone fluctuations are due primarily to atmospheric gravity waves. At m = 1/(1000 m) the mean spectral amplitude is over 30 times larger in the lower stratosphere than in the troposphere. Mean vertical wavenumber spectra in area-preserving form reveal dominant vertical wavelengths of ∼2.6 km in the troposphere and ∼2.7 km in the lower stratosphere, which is consistent with the values varying between 1.5 and 3.0 km estimated from the velocity field and temperature field at these heights.     

Determination of K,Ar, Cl,S, Si and Al flare abundances from RESIK soft X-ray spectra

The RESIK is a high sensitivity, uncollimated bent crystal spectrometer which successfully operated aboard Russian CORONAS-F solar mission between 2001 and 2003. It measured for the first time in a systematic way solar soft X-ray spectra in the four wavelength channels from 3.3 Å to 6.1 Å. This range includes characteristic strong lines of H- and He-like ions of K, Ar, Cl, Si, S and Al in the respective spectral channels. A distinguishing feature of RESIK is its possibility of making reliable measurements of the continuum radiation in flares. Interpretation of line and the continuum intensities observed in vicinity of respective strong lines provides diagnostics of plasma temperature and absolute abundances of K, Ar, Cl, S, Si and Al in several flares. We analyzed the observed intensities of spectral lines and the nearby continuum using the CHIANTI v5.2 atomic data package. A specific, so-called “locally isothermal” approach has been used in this respect allowing us to make not only flare-averaged abundance estimates, but also to look into a possible variability of plasma composition during the course of flares.     

Temperature comparison between CHAMP radio occultation and TIMED/SABER measurements in the lower stratosphere

Global Positioning System (GPS) receiver on the CHAllenging Mini-satellite Payload (CHAMP) and the Sounding of the Atmosphere using Broadband Emission Radiometry (SABER) instrument, one of four on board the TIMED satellite, provide middle atmosphere temperature profiles by Radio Occultation (RO) and limb viewing infrared emission measurements, respectively. These temperature profiles retrieved by two different techniques in the stratosphere are compared with each other using more than 1300 correlative profiles in March, September and December 2005. The over-all mean differences averaged over 15 and 35 km are approximately −2 K and standard deviation is less than 3 K. Below 20 km of altitude, relatively small mean temperature differences ∼1 K are observed in wide latitudinal range except for June (during the SABER nighttime observation). In the middle to low latitudes, between 30°S and 30°N, the temperature difference increases with height from ∼0–1 K at 15 km, to ∼−4 K at 35 km of altitude. Large temperature differences about −4 to −6 K are observed between 60°S and 30°N and 31–35 km of altitude for all months and between 0° and 30°N below 16 km during June (nighttime).     

对流层顶变化对上对流层/下平流层臭氧分布的影响

上对流层和下平流层(UT／LS)，位于8-25km高度之间，是大气中一个很特殊的区域．大部分的臭氧分布在下平流层，在下平流层臭氧的含量发生一个很小的变化，就会对气候和地面的紫外辐射产生很大的影响．而作为气象参数的对流层顶，是充分混合、缺乏臭氧的上对流层和层结稳定、臭氧丰富的下平流层之间的边界或过渡层，其变化对臭氧总量和分布有直接和明显的影响．本文使用二维模式模拟研究对流层顶变化对臭氧在UT／LS分布的影响．模拟结果表明对流层顶的季节变化对UT/LS的臭氧分布有明显的影响，臭氧的局地变化可以超过10％在冬季北半球中纬度对流层顶高度升高1km时，模式结果表明对臭氧分布的影响比较显著，局地变化可超过6％，但是对臭氧总量的影响较小，变化不超过5DU，小于观测资料统计分析的结果。     

The GFZ real-time GNSS precise positioning service system and its adaption for COMPASS

Motivated by the IGS real-time Pilot Project, GFZ has been developing its own real-time precise positioning service for various applications. An operational system at GFZ is now broadcasting real-time orbits, clocks, global ionospheric model, uncalibrated phase delays and regional atmospheric corrections for standard PPP, PPP with ambiguity fixing, single-frequency PPP and regional augmented PPP. To avoid developing various algorithms for different applications, we proposed a uniform algorithm and implemented it into our real-time software. In the new processing scheme, we employed un-differenced raw observations with atmospheric delays as parameters, which are properly constrained by real-time derived global ionospheric model or regional atmospheric corrections and by the empirical characteristics of the atmospheric delay variation in time and space. The positioning performance in terms of convergence time and ambiguity fixing depends mainly on the quality of the received atmospheric information and the spatial and temporal constraints. The un-differenced raw observation model can not only integrate PPP and NRTK into a seamless positioning service, but also syncretize these two techniques into a unique model and algorithm. Furthermore, it is suitable for both dual-frequency and sing-frequency receivers. Based on the real-time data streams from IGS, EUREF and SAPOS reference networks, we can provide services of global precise point positioning (PPP) with 5–10 cm accuracy, PPP with ambiguity-fixing of 2–5 cm accuracy, PPP using single-frequency receiver with accuracy of better than 50 cm and PPP with regional augmentation for instantaneous ambiguity resolution of 1–3 cm accuracy. We adapted the system for current COMPASS to provide PPP service. COMPASS observations from a regional network of nine stations are used for precise orbit determination and clock estimation in simulated real-time mode, the orbit and clock products are applied for real-time precise point positioning. The simulated real-time PPP service confirms that real-time positioning services of accuracy at dm-level and even cm-level is achievable with COMPASS only.     

Theoretical simulation of electron density and temperature distribution at Indian equatorial and low latitude ionosphere

The electron density and temperature distribution of the equatorial and low latitude ionosphere in the Indian sector has been investigated by simultaneously solving the continuity, momentum and energy balance equations of ion and electron flux along geomagnetic field lines from the Northern to the Southern hemisphere. Model algorithm is presented and results are compared with the electron density and electron temperature measured in situ by Indian SROSS C2 satellite at an altitude of ∼500 km within 31°S–34°N and 75 ± 10°E that covers the Indian sector during a period of low solar activity. Equatorial Ionization Anomaly (EIA) observed in electron density, morning and afternoon enhancements, equatorial trough in electron temperature have been simulated by the model within reasonable limits of accuracy besides reproducing other normal diurnal features of density and temperature.     

On the measurements of the moon’s infrared temperature and its relation to the phase angle

Radiometric measurements of the thermal radiation originating from the moon’s surface were obtained using an infrared detector operating at wavelengths between 8 and 14 μm. The measurements cover a full moon cycle. The variation of the moon’s temperature with the lunar phase angle was established. The lunar temperatures were 391 ± 2.0 K for the full moon, 240 ± 3.5 K for the first quarter, and 236 ± 3 K for the last quarter. For the rest of the phase angles, the lunar temperature varied between 170 and 380 K. Our results are comparable with those obtained previously at these phase angles. For the new moon phase, the obtained temperature was between 120 and 133 K. With the exception of the new moon phase, our measurements at all the phase angles were consistent with those obtained using Earth-based data and those obtained by the Diviner experiment and the Clementine spacecraft. At the new phase, our measurements were comparable with those obtained from the ground but were significantly higher than those obtained by the Diviner and Clementine data. We attribute this inconsistency to either the calibration curve of our detector, which does not perform well at very low temperatures, or to infrared emission from the atmosphere. A simple linear model to predict the lunar temperature as a function of the phase angle was proposed. The experimental errors that affect the measured temperatures are discussed.     

Sensor calibration impacts on dust detection based on MODIS and VIIRS thermal emissive bands

Dust detection using remotely sensed measurements has been one of the challenging problems encountered by atmospheric scientists. MODerate Resolution Imaging Spectroradiometer (MODIS) on the Terra (T) and Aqua (A) platforms have been a versatile sensor for well over 21 and 18 years respectively, and have been extremely useful in the retrieval of aerosol information over the entire globe. The MODIS radiances from the Level 1B in general are expected to be within 5% accuracy in the reflective wavelengths and within 1% in the thermal emissive wavelengths. In this paper, we evaluate the sensitivity of previously developed dust detection technique based on thermal emissive wavelengths, which correspond to MODIS bands 20, 29, 31, and 32 respectively. The Thermal Emissive Dust Index (TEDI) performed very comparably to the traditional Aerosol Optical Thickness (AOT) retrievals by MODIS reflective channels. Since the MODIS Thermal Emissive Bands (TEB) are well calibrated on-orbit using a BlackBody (BB) source, the calibration of these long wave infrared bands is quite robust. As A-MODIS continues to perform well beyond its designed lifetime of 6 years, the instrument has undergone various levels of degradation during its mission time. As a consequence, it is imperative to check the impacts of calibration on the higher-level retrievals. In this paper, we rigorously analyze the sensitivity of TEDI due to the impact of calibration by the afore-mentioned TEB. The perturbation of the dominant (linear) calibration term demonstrated the following: first, there was a correlation in the sensitivity of the TEDI due to the uncertainty in the linear calibration term. Based on a perturbation in the linear calibration term for all aforementioned bands over a range of ±5% yielded the TEDI sensitivity to vary from approximately ?3.2% to about ?3.6%. When considering the uncertainty in each individual band significant changes were observed. The least change was observed for the perturbation in the calibration of band 20 with the TEDI sensitivity and the largest sensitivity in TEDI was observed in the perturbation of band 31 calibration. Thus, in the case of TEDI, noticeable sensitivity due to calibration uncertainty was observed in bands 29, 31, and 32, reiterating the importance of the TEB calibration in these bands. Also, the dust detection scheme based on A-MODIS was successfully transferred to the follow-on sensors such as Suomi (SNPP) and NOAA 20 (N20) VIIRS. The results presented in this paper would be extremely helpful in understanding impacts of calibration on the higher-level products for both current and future missions based on the MODIS heritage. Finally, the work also identifies the importance of radiometric fidelity in maintaining the accuracy of the dust detection. Results presented will show drastic improvement of the Saharan dust detection after the reduction of the electronic crosstalk in the 8.5 µm channel of T-MODIS.     

Comparison of RO tropopause height based on different tropopause determination methods

The tropopause is an important boundary in the Earth’s atmosphere, and has been the subject of close attention from atmosphere and climate researchers. To monitor the global tropopause using radio occultation (RO) data, there are two primary methods, one is the widely used temperature lapse rate method, and the other is the bending angle covariance transform method which is unique to RO data. We use FengYun3-C (FY3C) and Meteorological Operational Satellite Program (MetOp) RO data and European Centre for Medium-Range Weather Forecasts (ECMWF) operational analysis data to determine differences in RO tropopause height calculated by these two methods. We compute biases of the RO lapse rate tropopause height (LRTH) and the RO bending angle tropopause height (BATH) relative to the ECMWF LRTH. The dependences of the tropopause height biases on tropopause height (TPH) retrieval method, latitude, season and RO mission are investigated. The results indicate that BATH show a consistent 0.8–1.2 km positive bias over the tropics and high latitude regions compared with LRTH, however, over 25° to 40° latitude in both hemisphere, BATH results are less stable. Furthermore, the mean bias between BATH and LRTH displayed a different symmetrical characteristic from 2017.12 to 2018.2 (DJF) compared to 2018.6–2018.8 (JJA). However, except for some bias over Antarctica, the mean value of both LRTH and BATH show a similar tropopause variation, indicating the consistency of both methods.     

Recent results on the Venus atmosphere from pioneer Venus radio occultations

Radio occultation measurements of the temperature structure of the Venus atmosphere have been obtained during seven occultation “seasons” extending from December 1978 to December 1983. Approximately 123 vertical profiles of temperature from about 40 km to about 85 km altitudes have been derived. Since these measurements cover latitudes from both poles to the equator, they have shown the latitudinal dependence of thermal structure. There is a smooth transition from the troposphere to the mesosphere at latitudes below about 45°, with the tropopause at about 56 km. The troposphere then rises to about 62 km in the “collar cloud” region between about 60° and 80° latitude, where a strong temperature inversion (up to 30 K) is present. In the polar areas, 80°–90°, the mesosphere becomes isothermal and there is no inversion. This latitudinal behavior is related to the persistent circulation pattern, in which a predominantly zonal retrograde motion at latitudes below 45° gradually changes to a circumpolar vortex at the “collar cloud” latitudes. Indeed, the radio occultation data have been used in a cyclostrophic balance model to derive zonal winds in the Venus atmosphere, which showed a mid-latitude (50°–55°) jet with a speed of about 120–140 ms^?1 at about 70 km altitude /1,2/. The observations obtained in 1983 and 1984 have shown that above the tropopause there is considerable temporal variability in the detailed thermal structure, suggesting that the persistent circulation pattern is subject to weather-like variability.     

Deriving a sea surface temperature record suitable for climate change research from the along-track scanning radiometers

We describe the approach to be adopted for a major new initiative to derive a homogeneous record of sea surface temperature for 1991–2007 from the observations of the series of three along-track scanning radiometers (ATSRs). This initiative is called (A)RC: (Advanced) ATSR Re-analysis for Climate. The main objectives are to reduce regional biases in retrieved sea surface temperature (SST) to less than 0.1 K for all global oceans, while creating a very homogenous record that is stable in time to within 0.05 K decade⁻¹, with maximum independence of the record from existing analyses of SST used in climate change research. If these stringent targets are achieved, this record will enable significantly improved estimates of surface temperature trends and variability of sufficient quality to advance questions of climate change attribution, climate sensitivity and historical reconstruction of surface temperature changes. The approach includes development of new, consistent estimators for SST for each of the ATSRs, and detailed analysis of overlap periods. Novel aspects of the approach include generation of multiple versions of the record using alternative channel sets and cloud detection techniques, to assess for the first time the effect of such choices. There will be extensive effort in quality control, validation and analysis of the impact on climate SST data sets. Evidence for the plausibility of the 0.1 K target for systematic error is reviewed, as is the need for alternative cloud screening methods in this context.