Seasonally steady planetary disturbances in the troposphere and stratosphere as seen in 30 years of NCEP reanalysis data

Zonal velocity and temperature daily global reanalysis data of 30 years are used to search seasonally steady planetary disturbances in the middle troposphere (400 hPa) and middle stratosphere (10 hPa). Significant wavenumber 1, 2 and 3 modes are found. Constant phase lines of zonal velocity 1 modes exhibit significant inclination angles with respect to the meridians. The winter hemisphere generally shows a more significant presence of structures. The Northern Hemisphere (NH) exhibits all over the year a larger amount of structures and more intense amplitudes than the Southern Hemisphere (SH). Middle latitudes exhibit the most significant cases and low latitudes the least significant ones. Longitudinally oriented land–sea transitions at ±$\pm$65° and −35° latitudes appear to play a significant role for the presence of steady planetary modes. The stratosphere exhibits a much simpler picture than the troposphere. Large scale structures with respectively NE–SW (NH) and NW–SE (SH) tilts in the observed temperature and zonal velocity constant phase lines recall the quasi-stationary Rossby wave trains that favor the poleward transport of angular momentum.     

Response of the extratropical middle atmosphere to the September 2002 major stratospheric sudden warming

The effects of a major stratospheric sudden warming (SSW) at extratropical latitudes have been investigated with wind and temperature observations over a Brazilian station, Cachoeira Paulista (22.7°S, 45°W) during September–October 2002. In response to the warming at polar latitudes a corresponding cooling at tropical and extratropical latitudes is prominent in the stratosphere. A conspicuous signature of latitudinal propagation of a planetary wave of zonal wavenumbers 1 and 2 from polar to low latitude has been observed during the warming period. The polar vortex which split into two parts of different size is found to travel considerably low latitude. Significant air mass mixing between low and high latitudes is caused by planetary wave breaking. The meridional wind exhibits oscillations of period 2–4 days during the warming period in the stratosphere. No wave feature is evident in the mesosphere during the warming period, although a 12–14 day periodicity is observed after 2 weeks of the warming event, indicating close resemblance to the results of other simultaneous investigations carried out from high latitude Antarctic stations. Convective activity over the present extratropical station diminishes remarkably during the warming period. This behavior is possibly due to destabilization and shift of equatorial convective active regions towards the opposite hemisphere in response to changes in the mean meridional circulation in concert with the SSW.     

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.     

Ozone mini-hole observation over the Balkan Peninsula in March 2005

Areas with dimensions of 1000–3000 km in which the total ozone content (TOC) decreases fast are called ozone mini-holes. They are generated mainly dynamically in two ways, either by poor-ozone air mass transport from the tropics to higher latitudes by planetary wave activity or, they are connected with strong adiabatic uplifting of the tropopause height. An ozone mini-hole, generated by the second mechanism, was observed over the Balkan Peninsula on 19/21 March 2005. In the middle of March, the polar vortex was strongly disturbed by Rossby waves, reaching up to the lower stratosphere. Warming episodes over a geographical area, covering the Barents Sea and the Polar Sea north from Central Siberia, displaced a polar vortex fragment extremely southwards. However, the vorticity was weak and the stratospheric temperatures did not reach low values, providing conditions for ozone chemical destruction via heterogenic reactions. At the same time, a Rossby wave ridge was located below the European polar fragment. In the period from 13 to 19 March, the thermal tropopause over Sofia was uplifted almost by 3 km. Ozone distributions observed by the SCIAMACHY instrument on 18–21 March show a fast TOC decrease westwards from Ireland, which was moving eastwards during the next days, increasing the area in which the ozone content decreased. On 20/21 March low ozone content was observed above the Stara Zagora (42°N, 25°E) ground-based station by means of the GASCOD instrument, using DOAS technique. The TOMS Earth probe instrument detected 237 DU over Sofia. This is a record low March value from the beginning of the TOMS instrument measurements in 1978. In March/April the ozone distribution was characterized by its mean annual maximum of 360 DU at 42°N.     

Establishment of atmospheric weighted mean temperature model in the polar regions

Due to the special geographical location and extreme climate environment, the polar regions (Antarctic and Arctic) have an important impact on global climate change. Atmospheric weighted mean temperature (Tm) is a crucial parameter in the retrieval of precipitable water vapor (PWV) from the zenith wet delay (ZWD) of ground-based Global Navigation Satellite System (GNSS) signal propagation. In this paper, the correlation between weighted mean temperature and surface temperature (Ts) is studied firstly. It is shown that the correlation coefficients between Tm and Ts are 0.93 in the Antarctic and 0.94 in the Arctic. The linear regression Tm model and quadratic function Tm model of the Antarctic and the Arctic are established respectively using the radiosonde profiles of 12 stations in the Antarctic and 58 stations in the Arctic from 2008 to 2015. The accuracies of the linear regression Tm model, the quadratic function Tm model and GPT2w Tm model which is a state-of-the-art global Tm model are verified using the radiosonde profiles from 2016 to 2018 in the Antarctic and Arctic. Root Mean Square (RMS) errors of the linear regression Tm model, the quadratic function Tm model and GPT2w Tm model in the Antarctic are 3.07 K, 2.87 K and 4.32 K respectively, and those in the Arctic are 3.53 K, 3.38 K and 4.82 K, which indicates that the quadratic function Tm model has a higher accuracy compared to linear regression Tm model, and the accuracies of the two regional Tm models are better than that of GPT2w Tm model in the polar regions. In order to better evaluate the accuracy of Tm in the PWV retrieval, the PWV values of radiosondes are used for comparisons as the reference value. The RMS errors of PWV derived from the two Tm models are similar for 1.28 mm in the Antarctic and 1 mm in the Arctic respectively. In addition, the spatial and temporal variation characteristics of Tm are analyzed in the polar regions by spectral analysis of Tm data using fast Fourier transform. The results show that the Tm has obvious seasonality and annual periodicity in the polar regions, and the maximum difference between warm season and cold season is about 63 K. After comparing and analyzing the influences of latitude, longitude and elevation on the Tm in the polar regions, it is found that latitude and elevation have a greater influence on the Tm than the longitude. As the latitude and elevation increase, the Tm decreases, and vice versa in the polar regions.     

Height and time variability of planetary wave activity

The height–season and year-to-year regularities of parameters of first and second spatial harmonics determine the structure of the stratosphere and mesosphere circulation and its variability. In the period 1992–2002 at heights 0–55 km, the amplitudes and phases of the first and second spatial harmonics in the field of temperature, geopotential height, zonal and meridional wind were calculated by the method of harmonic decomposition. Dispersion (standard or mean square deviation) of their day-to-day and year-to-year variations was calculated by their wavelength constants. Height and season patterns of variability have been estimated. The difference in height–longitude variability for wave numbers m = 1 and 2 has been discovered. At the same time, the intensity of wave disturbances for m = 1 is less than for m = 2 excluding the polar areas, where a significant variability appears at the heights 0–55 km. There is also a tendency for the intensity of year-to-year variations to decrease in comparison with day-to-day variations. In cold and warm periods the amplitude of perturbation waves with m = 2 both for day-to-day and year-to-year variations is greater than for waves with m = 1. Transient height areas in the interval of 20–30 km are more distinct for day-to-day variations of polar area.     

Global structure,seasonal and interannual variability of the eastward propagating tides seen in the SABER/TIMED temperatures (2002–2007)

The present paper is focused on the global spatial (altitude and latitude) structure, seasonal and interannual variability of the most stable in amplitude and phase eastward propagating diurnal and semidiurnal tides with zonal wavenumbers 2 and 3 derived from the SABER/TIMED temperatures for full 6 years (January 2002–December 2007). The tidal results are obtained by an analysis method where the tides (migrating and nonmigrating) and the planetary waves (zonally travelling, zonally symmetric and stationary) are simultaneously extracted from the satellite data. It has been found that the structures of the eastward propagating diurnal tides with zonal wavenumbers 3 and 2 change from antisymmetric with respect to the equator below ∼85 km height, to more symmetric above ∼95 km. The seasonal behavior of the DE3 is dominated by annual variation with maximum in August–September reaching average (2002–2007) amplitude of ∼15 K, while that of the DE2 by semiannual variation with solstice maxima and with average amplitude of ∼8 K. These tides revealed some interannual variability with a period of quasi-2 years. The seasonal behavior of the eastward propagating semidiurnal tide with zonal wavenumber 2 in the southern hemisphere (SH) is dominated by annual variation with maximum in the austral summer (November–January) while that in the northern hemisphere (NH) by semiannual variation with equinoctial maxima. The SE2 maximizes near 115 km height and at latitude of ∼30° reaching an average amplitude of ∼6 K. The seasonal behavior of the eastward propagating semidiurnal tide with zonal wavenumber 3 in both hemispheres indicates a main maximum during June solstice and a secondary one during December solstice. The tide maximizes near 110–115 km height and at a latitude of ∼30° reaching an average amplitude of ∼4.8 K in the SH and ∼4 K in the NH. The tidal structures of the two eastward propagating semidiurnal tides are predominantly antisymmetric about the equator.     

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.     

A superposed epoch study of the effects of solar wind stream interface events on the upper mesospheric and lower thermospheric temperature

The response of mesosphere and lower thermosphere (MLT) temperature to energetic particle precipitation over the Earth’s polar regions is not uniform due to complex phenomena within the MLT environment. Nevertheless, the modification of MLT temperatures may require an event-based study to be better observed. This work examines the influence of precipitation, triggered by solar wind stream interfaces (SI) event from 2002 to 2007, on polar MLT temperature. We first test the relationship between the ionospheric absorption measured by the SANAE IV (South African National Antarctic Expedition IV) riometer and the layer of energetic particle precipitation from POES (Polar Orbiting Environmental Satellites). The combined particle measurements from POES 15, 16, 17 and 18 were obtained close in time to the pass of the SABER (Sounding of the Atmosphere using Broadband Emission Radiometry) temperature retrieval. Here, a superposed epoch technique is described and implemented to obtain average temperature profiles during SI-triggered particle precipitation. The superposed epoch average shows no significant temperature decrease below 100 km prior to the onset of SI-triggered precipitation, whereas a clear superposed average temperature decrease is observed at 95 km after the SI impact. A case study of SI event also yields similar observations. Results indicate that cooling effects due to the production of mesospheric odd hydrogen might be major contributors to temperature decrease under compressed solar wind stream.     

The impact on tsunami detection from space using GNSS-reflectometry when combining GPS with GLONASS and Galileo

The devastating Sumatra tsunami in 2004 demonstrated the need for a tsunami early warning system in the Indian Ocean. Such a system has been installed within the German-Indonesian Tsunami Early Warning System (GITEWS) project. Tsunamis are a global phenomenon and for global observations satellites are predestined. Within the GITEWS project a feasibility study on a future tsunami detection system from space has therefore been carried out. The Global Navigation Satellite System Reflectometry (GNSS-R) is an innovative way of using GNSS signals for remote sensing. It uses ocean reflected GNSS signals for sea surface altimetry. With a dedicated Low Earth Orbit (LEO) constellation of satellites equipped with GNSS-R receivers, densely spaced sea surface height measurements could be established to detect tsunamis. Some general considerations on the geometry between LEO and GNSS are made in this simulation study. It exemplary analyzes the detection performance of a GNSS-R constellation at 900 km altitude and 60° inclination angle when applied to the Sumatra tsunami as it occurred in 2004. GPS is assumed as signal source and the combination with GLONASS and Galileo signals is investigated. It can be demonstrated, that the combination of GPS and Galileo is advantageous for constellations with few satellites while the combination with GLONASS is preferable for constellations with many satellites. If all three GNSS are combined, the best detection performance can be expected for all scenarios considered. In this case an 18 satellite constellation will detect the Sumatra tsunami within 17 min with certainty, while it takes 53 min if only GPS is considered.     

Estimation of the Love and Shida numbers: h2, l2 using SLR data for the low satellites

In this paper we present results for the global elastic parameters: Love number h₂ and Shida number l₂ derived from the analysis of Satellite Laser Ranging (SLR) data. SLR data for the two low satellites STELLA (H = 800 km) and STARLETTE (H = 810 km) observed during 2.5 years from January 3, 2005 until July 1, 2007 with 18 globally distributed ground stations were analyzed. The analysis was done separately for the two satellites. We do a sequential analysis and study the stability and convergence of the estimates as a function of length of the data set used.     

Global water vapour distribution in the upper troposphere and lower stratosphere during CRISTA 2

Water vapour measurements during the second mission of the CRyogenic Infrared Spectrometers and Telescopes for the Atmosphere (CRISTA) instrument are presented in the altitude regime 8–20 km. Mixing ratios are shown on isentropic surfaces (300–500 K) as global zonal means and as averages in 60° longitude sectors. Transports are indicated to occur preferentially on isentropic surfaces in the northern hemisphere, but not in the tropics and in the south. The hygropause is found shifted away from the equator towards the winter hemisphere. The amount of shift is longitude dependent and can be as great as 20–30°. The water content of the hygropause area shows considerable zonal asymmetries.     

Simulation software for the spiral trajectory of our Moon

At 4.56 Ga, the accretion of the slowly rotating Solar Nebula led to the formation of Sun and its Planets in the plane of disc of accretion. Moon was formed by accretion from a circumterrestrial disk of debris generated by the glancing angle impact of the young Earth by a Mars size planetary embryo at about 4.5 Ga at a distance of 15,000 km. The Moon since then has migrated to the present position of 384,400 km from the center of the Earth. In course of this outward migration it has slowed down the spin rate of Earth and caused the lengthening of diurnal day length from 5 h initially to 24 h presently. The basic mechanics of Earth–Moon System has been worked out and theoretical determination of lengthening of day curve is carried out. This theoretical lengthening of day curve is compared with the observed lengthening of day curve based on paleobotanical evidences, ancient tidalites and Australian Banded Iron Formation. There is a remarkable correspondence between the two curves except for intermittent deviations due to geographical and geophysical factors. Based on the theoretical curve of lengthening of day, an empirical formula for the lunar orbital radius expansion is determined. Based on this empirical formula, simulation software is developed that gives the correct evolution of the semi-major axis (a) of our Moon for any time span from the inception to the time chosen under study. For mathematical simplicity the system is considered to be a two body rotating system throughout its evolutionary history of 4.5 Gyrs. This simulation draws the Moon’s spiral trajectory from its inception to any subsequent epoch. The terminal epoch is an input to the simulation software to arrive at the spiral trajectory of the Moon from the inception to the given epoch. The basic mechanics of Earth–Moon System and this simulation can be generalized to lay the foundation of simulation software for any Planet–Satellite pair or any Sun-Planet pair in our Solar System or Star-Planet pair in any Extra-Solar System. The basic dynamics has been found to be valid for Star–Planet pair also. So this Simulation Methodology can as well be applied to study the migratory evolution of Gas Giants also.     

An update global model of hmF2 from values estimated from ionosonde and COSMIC/FORMOSAT-3 radio occultation

In this paper, we present our recent work on developing an updated global model of the ionospheric F2 peak height hmF2 parameter by combining data from the Constellation Observing System for Meteorology, Ionosphere and Climate (COSMIC/FORMOSAT-3) radio occultation (RO) measurements and from the extended global ionosonde stations. In particular, 10 Chinese ionosonde stations’ data are newly introduced into this study. The modeling technique used is based on a two-layer empirical orthogonal function (EOF) expansion. Global distributions of hmF2 maps calculated using the newly constructed global model and the one provided by the International Reference Ionosphere model (IRI-ITU-R) are compared with the global distributions of hmF2 obtained by the COSMIC RO measurements and quantitative statistical analysis of the differences between the model results and those of the COSMIC RO measurements is made for the low (2008) and high (2012) solar activity years. The obtained average root-mean-square differences (RMSEs) for our model are 27.7 km (11.1%) and 31.0 km (9.8%), respectively for the years 2008 and 2012, whereas those for the IRI-ITU-R model are 39.9 km (16.9%) and 35.0 km (11.6%), respectively. Comparison of the results calculated both by our model and the IRI-ITU-R model with the digisonde observation is also made. The comparisons show that the newly constructed global hmF2 model can reproduce reasonably well the observations and perform better than IRI-ITU-R model.     

Observation of the lunar topography by the laser altimeter LALT on board Japanese lunar explorer SELENE

The SELENE Laser Altimeter (LALT) is designed to map the Moon’s topography and will be launched in summer 2007. LALT incorporates Q-switched Cr doped Nd:YAG laser (1064 nm) with an output energy of 100 mJ and 1 Hz repetition frequency for about one year mission period. The laser pulse travels to the Moon’s surface and reflections from the surface are detected by a silicon avalanche photo-diode. The ranging distance is 50–150 km with about 5 m accuracy. Several corrections for accurate ranging data are investigated. The flight hardware has been qualified and passed all the integration tests. A principal goal of the LALT instrument is to obtain a much more detailed lunar topographic map which is superior in global coverage, measurement accuracy and number of data points to previous observations and models. The overall science objectives of LALT are (1) determination of lunar global figure, (2) internal structure and surface processes, (3) exploration of the lunar pole regions, and (4) reduction of lunar occultation data.     

Study of seasonal variation of winds in upper stratosphere over Hyderabad

A predictability of the stratospheric zonal winds above 38 km during the turnaround is an essential parameter for planning of the high-altitude scientific balloon flights. This information is more relevant in the case of Hyderabad balloon facility which is closer to equator and has much more unstable wind reversal patterns which appears to have changed enormously during the last decade probably in correlation with the global warming. With a majority of our flights reaching the altitudes of 38–42 km and the requirement of long float durations, a prior knowledge of wind pattern during the summer and winter turnaround seasons is highly desirable. Furthermore, the flight operation corridor for balloon flights from Hyderabad is limited to 400 km and though in the west direction there are flat lands, in all other three directions, the landscape is dotted by water bodies, reserve forests and hilly terrain, and therefore need of such a data is essential. In order to establish the climatology of the stratospheric winds and study their inter-annual variability over Hyderabad for the turnaround periods, we have made a detailed analysis of the United Kingdom Meteorological office data between 2000 and 2007, to derive average wind parameters (magnitude, direction) at different ceiling altitudes above 38 km. These results can be used only as general trend of stratospheric wind and should not be the limitation of the UKMO Data.     

Zonal mean temperature, pressure, zonal wind and geopotential height as functions of latitude

The new zonal mean COSPAR International Reference Atmosphere (CIRA-86) of temperature, zonal wind, and geopotential/geometric height is presented. This data can be used as a function of altitude or pressure and has nearly pole-to-pole coverage (80°S-80°N) extending from the ground to approximately 120 km. Data sources and methods of computation are described; in general, hydrostatic and thermal wind balance are maintained at all levels and latitudes. As shown by a series of cross sectional plots, the new CIRA accurately reproduces most of the characteristic features of the atmosphere such as the equatorial wind and the general structure of the tropopause, stratopause, and mesopause.     

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).     

Impact of GPS antenna phase center variations on precise orbits of the GOCE satellite

The first European Space Agency Earth explorer core mission GOCE (Gravity field and steady-state Ocean Circulation Explorer) has been launched on March 17, 2009. The 12-channel dual-frequency Global Positioning System receiver delivers 1 Hz data and provides the basis for precise orbit determination (POD) on the few cm-level for such a very low orbiting satellite (254.9 km). As a member of the European GOCE Gravity Consortium, which is responsible for the GOCE High-level Processing Facility (HPF), the Astronomical Institute of the University of Bern (AIUB) provides the Precise Science Orbit (PSO) product for the GOCE satellite. The mission requirement for 1-dimensional POD accuracy is 2 cm. The use of in-flight determined antenna phase center variations (PCVs) is necessary to meet this requirement. The PCVs are determined from 154 days of data and the magnitude is up to 3-4 cm. The impact of the PCVs on the orbit determination is significant. The cross-track direction benefits most of the PCVs. The improvement is clearly seen in the orbit overlap analysis and in the validation with independent Satellite Laser Ranging (SLR) measurements. It is the first time that SLR could validate the cross-track component of a LEO orbit.     

Temperature variability in the tropical mesosphere during the northern hemisphere winter

Temperature observations at 20–90 km height and 5°N–15°N during the period of December 1992–March 1993 from the WINDII and MLS experiments on the UARS satellite are analysed together with MF radar winds and UKMO assimilated fields of temperature and zonal and meridional winds. The correlation between the different datasets at the tropics and zonal mean wind data at mid latitudes is examined for period February–March 1993, when series of stratospheric warming events were observed at middle and high latitudes. Wavelet analysis is applied to investigate coupling between stationary and travelling planetary waves in the stratosphere and the upper mesosphere. Planetary waves m = 1 with periods of 4–7 days, 8–12 days and 13–18 days are found to dominate the period. Westward 7- and 16–18 day waves at the tropics appear enhanced by stationary planetary waves during sudden stratospheric warming events.