Estimated errors in a global gravity wave climatology from GPS radio occultation temperature profiles

In a previous paper by Schmidt et al. (2008), from CHAllenging Minisatellite Payload (CHAMP) Global Positioning System (GPS) radio occultation data, a comparison was made between a Gaussian filter applied to the “complete” temperature profile and to its “separate” tropospheric and stratospheric height intervals, for gravity wave analyses. It was found that the separate filtering method considerably reduces a wave activity artificial enhancement near the tropopause, presumably due to the isolation process of the wave component. We now propose a simple approach to estimate the uncertainty in the calculation of the mean specific wave potential energy content, due exclusively to the filtering process of vertical temperature profiles, independently of the experimental origin of the data. The approach is developed through a statistical simulation, built up from the superposition of synthetic wave perturbations. These are adjusted by a recent gravity wave (GW) climatology and temperature profiles from reanalyses. A systematic overestimation of the mean specific wave potential energy content is detected and its variability with latitude, altitude, season and averaging height interval is highlighted.     

Uncertainties in gravity wave parameters,momentum fluxes,and flux divergences estimated from multi-layer measurements of mesospheric nightglow layers

Measurements of dynamic parameters of atmospheric gravity waves, mainly the vertical wavelength, the momentum flux and the momentum flux divergence, are affected by large uncertainties crudely documented in the scientific literature. By using methods of error analysis, we have quantified these uncertainties for frequently observed temporal and spatial wave scales. The results show uncertainties of ～10%, ～35%, and ～65%, at least, in the vertical wavelength, momentum flux, and flux divergence, respectively. The large uncertainties in the momentum flux and flux divergence are dominated by uncertainties in the Brunt-Väisälä frequency and in spatial separation of the nightglow layers, respectively. The measured uncertainties in fundamental wave parameters such as the wave amplitude, intrinsic period, horizontal wavelength, and wave orientation are ～10% or less and estimated directly from our nightglow image data set. Other key environmental quantities such as the scale height and the Brunt-Väisälä frequency, frequently considered as constants in gravity wave parameter estimations schemes, are actually quite variable, presenting uncertainties of ～4% and ～9%, respectively, according to the several solar activity and seasonal atmosphere scenarios from the NRLMSISE-00 model simulated here.     

Application of heterogeneous mesospheric ion-chemistry model to MST radar echoes over low latitude

MST radar studies at low latitude stations have documented regions in the mesosphere from where enhanced echoes (Low Latitude Mesospheric Echoes (LMEs)) are observed. Such echoes cannot, in general, be explained by considering the dynamical aspects (such as turbulence, winds, waves, etc.) of the region alone. Mesospheric dust/aerosols can enhance the radar echoes considerably and dust is known to exist at all heights and latitudes of the mesosphere. This study investigates the presence of dusty plasma in the mesosphere through the heterogeneous ion-chemistry of the region.Dust of meteoric origin is incorporated in the conventional ion chemistry scheme and the equilibrium height profiles of charged and neutral dust densities corresponding to effective dust sizes (radii) of 1, 10 and 30 nm are computed for the equatorial quiet daytime conditions.The model derived dust density profiles show structures with respect to dust size, height and season that are indicative of the possible role of mesospheric dust in the production/enhancement mechanisms of the LMEs observed over the equatorial station at Gadanki (13.5°N, 79.2°E), India.     

Comparative evaluation of five global gravity models over a part of the Bay of Bengal

Several global gravity models (GGMs) are freely available in the public domain, which can be utilised to study the earth's gravity field in almost every part of the globe. The present study compared the free-air gravity anomalies calculated from the five GGMs EGM2008, EIGEN6C4, GECO, XGM2019e_2159, and SGG-UGM-2 archived by the International Centre for Global Earth Models (ICGEM) with respect to shipborne gravity in the Bay of Bengal. The average correlation coefficient and covariance are ～ 96 % and ～ 450mGal². The mean difference between the shipborne and the modelled gravity is ? 5 mGal. Relatively higher amplitude gravity differences observed at the continental-oceanic transition, the 85°E and Ninetyeast ridges, and the western basin are possibly due to high gradient, dominant density contrasts, and rugged topography. The average standard deviation and root-mean-square-error (RMSE) of the differences are ～ 6.5 mGal and ～ 7.5 mGal. A significantly lower standard deviation and RMSE found for the models generated at higher degree/order compared to lower degree/order is due to diminishing omission error of the GGMs with increasing degrees of truncation. The spectral analysis and coherence estimation of the modelled gravity demonstrate excellent correspondence for anomalies wider than ～ 25 km. The agreement between anomaly amplitudes and shapes and calculated statistics indicates that the accuracy and resolution of the modelled gravity data are certainly good enough for regional-scale studies aiming to map Moho topography and mantle structures. However, the delineation of shorter wavelength features originating from the shallow-level basement/sedimentary might be uncertain and requires further validations. The present study confirms that all five models show promising results in terms of their accuracy and resolution. Moreover, the SGG-UGM-2 and XGM2019e_2159 models compare favourably with respect to the GECO, EIGEN6C4 and EGM2008 models in the Bay of Bengal.     

Seasonal variation of gravity wave sources from satellite observation

Recent high-resolution satellite observations of gravity waves in the middle atmosphere have shown correlations with the strength of the stratospheric jet stream, surface topography, and tropical convection. Seasonal variations of wave-induced stratospheric radiance variances are often the manifestations of modulations of these sources and refractive influences. In this paper, we focus on the seasonal climatology of gravity waves observed by the UARS MLS, while also showing some new results from GPS and AMSU instruments. Our analysis is aided by MWFM modeling of mountain waves at high latitudes and CMAP precipitation indices in the subtropics to provide a clearer picture of global gravity wave dynamics.     

A numerical model characterizing internal gravity wave propagation into the upper atmosphere

A new two-dimensional, time-dependent and fully nonlinear model is developed to numerically simulate plane wave motions for internal gravity waves in a non-isothermal and windy atmosphere that accounts for the dissipation due to eddy and molecular processes. The atmosphere has been treated as a well mixed total gas with a constant mean molecular weight. The effects of Rayleigh friction and Newtonian cooling are applied near the upper boundary of the model to simulate the radiation conditions as well as to act as a sponge layer; lateral boundaries are periodic over a horizontal wavelength to simulate a horizontally infinite domain. The thermal excitation to initiate upward propagating waves is spatially localized in the troposphere and is a Gaussian function of time. A time-splitting technique is applied to the finite difference equations that are derived from the Navier–Stokes equations. The time integration for these highly coupled equations is performed using an explicit second order Lax–Wendroff scheme and an implicit Newton–Raphson scheme. The wave solutions derived from the model are found to be broadly agreeable with those derived from a Wentzel–Kramers–Brillouin theory.     

Ground-based measurement of strato–mesospheric CO by a FTIR spectrometer over Poker Flat,Alaska

Upper atmospheric CO above 24 km has been observed over Poker Flat (147°W, 65°N, altitude 0.61 km), Alaska using ground-based solar absorption infrared spectroscopy. This is the first reported detection of stratospheric–mesospheric CO using this method from the ground. The results clearly indicate that there is a seasonal variation of the CO profile with enhanced abundances in spring while remaining low from May onwards.The Poker Flat Research Range is one of the many measurement sites that constitute the Network for the Detection of Stratospheric Change (NDSC). The method used in this work, estimating the CO partial column abundances above the middle stratosphere, can be applied to spectra observed using FTIR spectrometry at many other NDSC sites. This suggests the availability of this established technique as a new method for CO measurements in the upper atmosphere.     

A comparison between CRISTA satellite data and Warner and McIntyre gravity wave parameterization scheme: Horizontal and vertical wavelength filtering of gravity wave momentum flux

An along-track analysis of temperature height profiles measured by the CRISTA satellite instrument was performed to obtain height profiles of absolute values of gravity wave momentum flux. Global distributions of momentum flux from CRISTA are compared to results of the Warner and McIntyre gravity wave parameterization scheme. Only a limited range of horizontal and vertical wavelengths is visible to the CRISTA instrument. Therefore the model results have to be filtered accordingly to make comparisons meaningful. We demonstrate the importance of such a wavelength filtering which significantly increases the correlation of observed and modeled momentum fluxes.     

Observations of wave activity in the ionosphere over South Africa in geomagnetically quiet and disturbed periods

The present paper deals with observations of wave activity in the period range 1–60 min at ionospheric heights over the Western Cape, South Africa from May 2010 to July 2010. The study is based on the Doppler type sounding of the ionosphere. The Doppler frequency shift measurements are supplemented with measurements of collocated Digisonde DPS-4D at SANSA Space Sciences, Hermanus. Nine geomagnetically quiet days and nine geomagnetically active days were included in the study. Waves of periods 4–30 min were observed during the daytime independent of the level of geomagnetic activity. Amplitudes of 10–30 min waves always increased between 14:00 and 16:15 UT (16:00–18:15 LT). Secondary maxima were observed between 06:00 and 07:00 UT (08:00–09:00 LT). The maximum wave amplitudes occurred close to the time of passage of the solar terminator in the studied region which is known to act as a source of gravity waves.     

The amphibian egg as a model system for analyzing gravity effects.

Amphibian eggs provide several advantageous features as a model system for analyzing the effects of gravity on single cells. Those features include large size, readily tracked intracellular inclusions, and ease of experimental manipulation. Employing novel gravity orientation as a tool, a substantial data base is being developed. That information is being used to construct a 3-D model of the frog (Xenopus laevis) egg. Internal cytoplasmic organization (rather than surface features) are being emphasized. Several cytoplasmic compartments (domains) have been elucidated, and their behavior in inverted eggs monitored. They have been incorporated into the model, and serve as a point of departure for further inquiry and speculation.     

Observations of magnetospheric convection from low altitudes

Some of the observations and interpretive models that have provided a substantial increase in our knowledge of magnetospheric and ionospheric convection are discussed. While a two-cell convection pattern may be generally consistent with many ionospheric measurements, it is now clear that some significant departures from such a pattern must be considered. We can now specify more accurately the number of convection cells and their shape as well as the electrostatic potential distribution within the cells. All these factors can be shown to be sensitive functions of the interplanetary magnetic field (IMF). Interpretation of these findings in terms of the interaction of the earth's magnetosphere with the interplanetary medium has led to detailed consideration of the location of magnetic merging regions and the magnetic field topology of the outer magnetosphere. In addition, the relative importance of merging, viscous interaction and ionospheric processes in providing the driving force for convection has been considered. In general, the bulk of the driving force is magnetic reconnection; however, viscous processes play a significant role in times of northward interplanetary magnetic fields, and thermospheric drag may contribute to the maintainence of a convection pattern for several hours after such a northward turning.     

Empirical model of the TEC response to the geomagnetic activity over the North American region

The paper presents an empirical model of the total electron content (TEC) response to the geomagnetic activity described by the K_p-index. The model is built on the basis of TEC measurements covering the region of North America (50°W–150°W, 10°N–60°N) for the period of time between October 2004 and December 2009. By using a 2D (latitude-time) cross-correlation analysis it is found that the ionospheric response to the geomagnetic activity over the considered geographic region and at low solar activity revealed both positive and negative phases of response. The both phases of the ionospheric response have different duration and time delay with respect to the geomagnetic storm. It was found that these two parameters of the ionospheric response depend on the season and geographical latitude. The presence of two phases, positive and negative, of the ionospheric response imposed the implementation of two different time delay constants in order to properly describe the two different delayed reactions. The seasonal dependence of the TEC response to geomagnetic storms is characterized by predominantly positive response in winter with a short (usually ∼5–6 h) time delay as well as mainly negative response in summer with a long (larger than 15 h) time delay. While the TEC response in March and October is more close to the winter one the response in April and September is similar to the summer one.     

First observation of upper mesospheric semi annual oscillations using ground based airglow measurements from Indian low latitudes

We summarize two years of Mesosphere Lower Thermosphere Photometer (MLTP) operation of mesospheric OH and O₂ emission monitoring. The deduced mesospheric OH and O₂ temperatures show large variability. Nightly temperature variations over Gadanki (13.5°N, 79.2°E) are dominated by the short period wave features, while tidal amplitudes are relatively small. Our measurements are the first to report a long period seasonal variation at two upper mesospheric altitudes simultaneously over the Indian sector. Our observations reveal the presence of a dominant semi-annual oscillation (∼6 months periodicity) together with a shorter period (∼2.5  months periodicity) oscillation in both OH and O₂ data.     

Solar activity dependence of total electron content derived from GPS observations over Mbarara

Vertical total electron content (VTEC) observed at Mbarara (geographic co-ordinates: 0.60°S, 30.74°E; geomagnetic coordinates: 10.22°S, 102.36°E), Uganda, for the period 2001–2009 have been used to study the diurnal, seasonal and solar activity variations. The daily values of the 10.7 cm radio flux (F_10.7) and sunspot number (R) were used to represent Solar Extreme Ultraviolet Variability (EUV). VTEC is generally higher during high solar activity period for all the seasons and increases from 0600 h LT and reaches its maximum value within 1400 h–1500 h LT. All analysed linear and quadratic fits demonstrate positive VTEC-F_10.7 and positive VTEC-R correlation, with all fits at 0000 h and 1400 h LT being significant with a confidence level of 95% when both linear and quadratic models are used. All the fits at 0600 h LT are insignificant with a confidence level of 95%. Generally, over Mbarara, quadratic fit shows that VTEC saturates during all seasons for F_10.7 more than 200 units and R more than 150 units. The result of this study can be used to improve the International Reference Ionosphere (IRI) prediction of TEC around the equatorial region of the African sector.     

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.     

Analysis of global TEC prediction performance of IRI-PLAS model

The International Reference Ionosphere (IRI) empirical model provides valuable data for many fields including space and navigation applications. Since the IRI model gives the ionospheric parameters in the altitude range from 50?km to 2000?km, researchers focused on the IRI-PLAS model which is the plasmasphere extension of the IRI model. In this study, Total Electron Content (TEC) prediction performance of the IRI-PLAS model was examined at a global scale using the location of globally distributed 9 IGS stations. Besides the long term (01.01.2015–31.12.2015) behavior of the model, TEC predictions during the equinox and solstice days of 2014–2017 were also tested. IRI-PLAS-TEC values were examined in comparison with GPS-TEC data. Hourly interval of yearly profile exhibits that when the geomagnetic and solar active days are ignored, differences between IRI-PLAS-TEC and GPS-TEC are rather small (～2–3 TECU) at stations in the northern hemisphere, generally ～4–5 TECU level at the southern hemisphere stations and reaching above 10 TECU for few hours. While the IRI-PLAS-TEC generally overestimates the GPS-TEC at southern hemisphere stations during quiet days, the model-derived TEC underestimates GPS-TEC during solar active days. IRI-PLAS-TEC and GPS-TEC values exhibit similar trend for the equinoxes 21 March and 23 September which refer equivalent conditions.     

A gravity wave analysis near to the Andes Range from GPS radio occultation data and mesoscale numerical simulations: Two case studies

Global maps of potential wave energy per unit mass, recently performed with the Global Positioning System (GPS) Radio Occultation (RO) technique and different satellite missions (CHAMP and SAC-C since 2001, GRACE and COSMIC since 2006) revealed in Argentina, at the eastern side of the highest Andes Mountains, a considerable wave activity (WA) in comparison with other extra-tropical regions. The main gravity wave (GW) sources in this natural laboratory are deep convection (mainly during late Spring and Summer), topographic forcing and geostrophic adjustment.     

Quality assessment of sub-Nyquist recovery from future gravity satellite missions

Drawing on experience from Gravity Recovery and Climate Experiment (GRACE) data analysis, the scientific challenges were already identified in several studies. Any future mission should focus on improvement in both precision and resolution in space and time. For future gravity missions which use high quality sensors, aliasing of high frequency time-variable geophysical signals to the lower frequency signals is one of the most serious problems. The aliasing problem and the spatio-temporal resolution are mainly restricted by two sampling theorems describing the space-time sampling of satellite missions: (i) a Heisenberg-like uncertainty theorem which states that the product of spatial resolution and time resolution is constant, and (ii) the Colombo–Nyquist rule (CNR), which requires the number of satellite revolutions in a repeat period to be at least twice a given maximum spherical harmonic degree. The CNR holds under the assumption of equal ground-track spacing, and limits the spatial resolution of the gravity solution.