The contribution of local gravimetric geoid models to the calibration of satellite altimetry data and an outlook of the latest GOCE GGM performance in Gavdos

The use of geoid heights has been one of the available methodologies utilized for the independent calibration/validation of altimeters on-board satellites. This methodology has been employed for long in the Gavdos dedicated cal/val facility (Crete, Greece), where calibration results for the Jason satellites have been estimated, both for ascending and descending passes. The present work gives a detailed overview of the methodology followed in order to estimate a high-resolution and accuracy gravimetric geoid model for the wider Gavdos area, in support of the on-going calibration work. To estimate the geoid model, the well-known remove-compute-restore method is used while residual geoid heights are estimated through least-squares collocation so that associated errors are determined as well. It is found that the estimated formal geoid errors from LSC along passes 018 and 109 of Jason satellites, used for the bias estimation, range between ±0.8–1.6 cm. The so-derived geoid heights are employed in the determination of the Jason-2 altimeter bias for all available cycles (cycles 1-114, spanning the period from July 2008 to August 2011) together with the RioMed DOT model. From the results acquired the Jason-2 bias has been estimated to be +196.1 ± 3.2 mm for pass 109 and +161.9 ± 5.1 mm for pass 018. Within the same frame, the GOCE/GRACE-based geopotential model GOCO02s has been used to estimate the mean dynamic ocean topography and the steady-state circulation in the area around Gavdos. The so-derived DOT model was used to estimate the Jason-2 bias in an effort to evaluate the performance of satellite-only geoid models and investigate whether their spatial resolution and accuracy provides some improvement w.r.t. traditional local gravimetric geoids. From the results acquired with geoid heights from GOCO02s, the estimated Jason-2 bias deviates significantly from that of the local gravimetric model, which can be attributed to a possible mean offset and the low resolution of GOCE-based GGMs. On the other hand, when the newly estimated GOCE-based DOT was employed with geoid heights from the local gravimetric geoid model, the Jason-2 bias has been estimated to be +185.1 ± 3.2 mm for pass 109 and +130.2 ± 5.1 mm for pass 018.     

Solar and interplanetary origins of the November 2004 superstorms

During the first half of November 2004, many solar flares and coronal mass ejections (CMEs) were associated with solar active region (AR) 10696. This paper attempts to identify the solar and interplanetary origins of two superstorms which occurred on 8 and 10 November with peak intensities of Dst = −373 nT and −289 nT, respectively. Southward interplanetary magnetic fields within a magnetic cloud (MC), and a sheath + MC were the causes of these two superstorms, respectively. Two different CME propagation models [Gopalswamy, N., Yashiro, S., Kaiser, M.L. et al. Predicting the 1-AU arrival times of coronal mass ejections. J. Geophys. Res. 106, 29207–29219, 2001; Gopalswamy, N.S., Lara, A., Manoharan, P.K. et al. An empirical model to predict the 1-AU arrival of interplanetary shocks. Adv. Space Res. 36, 2289–2294, 2005] were employed to attempt to identify the solar sources. It is found that the models identify several potential CMEs as possible sources for each of the superstorms. The two Gopalswamy et al. models give the possible sources for the first superstorm as CMEs on 2330 UT 4 November 2004 or on 1454 UT 5 November 2004. For the second superstorm, the possible solar source was a CME that on 0754 UT 5 November 2004 or one that occurred on 1206 UT 5 November 2004. We note that other propagation models sometimes agree and other times disagree with the above results. It is concluded that during high solar/interplanetary activity intervals such as this one, the exact solar source is difficult to identify. More refined propagation models are needed.     

Measuring atmospheric density using GPS–LEO tracking data

We present a method to estimate the total neutral atmospheric density from precise orbit determination of Low Earth Orbit (LEO) satellites. We derive the total atmospheric density by determining the drag force acting on the LEOs through centimeter-level reduced-dynamic precise orbit determination (POD) using onboard Global Positioning System (GPS) tracking data. The precision of the estimated drag accelerations is assessed using various metrics, including differences between estimated along-track accelerations from consecutive 30-h POD solutions which overlap by 6 h, comparison of the resulting accelerations with accelerometer measurements, and comparison against an existing atmospheric density model, DTM-2000. We apply the method to GPS tracking data from CHAMP, GRACE, SAC-C, Jason-2, TerraSAR-X and COSMIC satellites, spanning 12 years (2001–2012) and covering orbital heights from 400 km to 1300 km. Errors in the estimates, including those introduced by deficiencies in other modeled forces (such as solar radiation pressure and Earth radiation pressure), are evaluated and the signal and noise levels for each satellite are analyzed. The estimated density data from CHAMP, GRACE, SAC-C and TerraSAR-X are identified as having high signal and low noise levels. These data all have high correlations with anominal atmospheric density model and show common features in relative residuals with respect to the nominal model in related parameter space. On the contrary, the estimated density data from COSMIC and Jason-2 show errors larger than the actual signal at corresponding altitudes thus having little practical value for this study. The results demonstrate that this method is applicable to data from a variety of missions and can provide useful total neutral density measurements for atmospheric study up to altitude as high as 715 km, with precision and resolution between those derived from traditional special orbital perturbation analysis and those obtained from onboard accelerometers.     

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.     

The hydroxyl radical as an indicator of SEP fluxes in the high-latitude terrestrial atmosphere

The low background values at nighttime of the mesospheric hydroxyl (OH) radical make it easier to single out the atmospheric response to the external solar forcing in Polar Regions. Because of the short lifetime of HO_x, it is possible to follow the trails of Solar Energetic Particle (SEP) events in the terrestrial atmosphere, as shown by Storini and Damiani (2008). The sensitivity of this indicator makes discernible not only extreme particle events with a flux peak of several thousand pfu [1 pfu = 1 particle/(cm² s sr)] at energies >10 MeV, but also those with lower flux up to about 300 pfu. Using data from the Microwave Limb Sounder (MLS) on board the EOS AURA satellite, we examined the correlation of OH abundance vs. solar proton flux for almost all the identified SEP events spanning from November 2004 to December 2006 (later on no more SEP events occurred during Solar Cycle no. 23). The channels at energies greater than 5 MeV and 10 MeV showed the best correlation values (r ∼ 0.90–0.95) at altitudes around 65–75 km whereas, as expected, the most energetic channels were most highly correlated at lower altitudes. Therefore, it is reasonably possible to estimate the solar proton flux from values of mesospheric OH (and viceversa) and it could be useful in studying periods with gaps in the records of solar particles.     

Solar particle precipitation into the polar atmosphere and their dependence on hemisphere and local time

The precipitation of solar energetic particles, protons as well as electrons, at high latitudes is commonly assumed to be homogeneous across both polar caps. Using Low-Earth Orbit POES (Polar Orbiting Environmental Satellites) we determine particle penetration ratios into the polar atmosphere for protons ranging from about 0.1 MeV to 500 MeV and for electrons spanning about one order of magnitude in energy with a maximum of 0.3 MeV. Based on power law fits for the POES spectrum we show, that for energies interesting for middle and lower atmosphere chemistry, particle flux over the poles is comparable in magnitude to flux at the geostationary orbit or at L1 in interplanetary space. The time period under study are the solar energetic particle (SEP) event series of October/November 2003 and January 2005.     

Height-latitude structure of stationary planetary waves in the stratosphere and lower mesosphere

Daily UK Met Office stratospheric assimilated data for the Northern and Southern Hemispheres, accumulated for the period from 2004 to 2012 and pressure range of 1000–0.1 hPa, are used in this paper. The paper presents and thoroughly discusses spatial–temporal distributions of stationary planetary wave (SPW) amplitudes and phases, calculated on the basis of geopotential height, temperature, zonal and meridional wind data for zonal wave numbers 1 and 2 (SPW1 and SPW2). The climatological planetary wave amplitudes and phases are calculated by extracting waves from three types of data: daily, monthly mean and climatological monthly mean. It has been established that magnitude of amplitudes and height-latitude distribution of amplitudes of SPW1 and SPW2 depend on data processing method for all parameters. It has also been established that height-latitude distribution amplitudes and phases significantly differ for geopotential height, temperature, zonal and meridional wind and depend on wave number and hemisphere. However, height-latitude distributions of phases are little different from each other for the used methods of data processing.     

Influence of different atmospheric models and seeing effects on the nonlinear calibration

This paper describes a technique intended to improve calibration of vector magnetograms achieved by the spectral line Fei 5324.19 Å which was first described in detail by Hagyard and Kineke [Hagyard, M.J., Kineke, J.I., Improved method for calibrating filter vector magnetographs, Solar Phys., 158, 11–28, 1995.]. The technique attempts to compensate for nonlinear behavior of the measured linear and circular polarizations of the line as functions of inclination angle and field strength. This method depends on different atmospheric models. To the three different atmospheric models, the VAL-C quiet photosphere model, the Ding and Fang penumbra model and Allen umbral model, the circular polarization signal is directly linearly proportional to the longitudinal field strength in the range of 800 G, 1000 G and 1300 G, respectively. For the transverse field, however, this linear relationship holds over much smaller ranges which are 200 G, 300 G and 200 G. The improved method improves the accuracy in calibrating the polarization signals over a wide range of field strengths and inclination angles. With the seeing effects decreasing in the magnetic observation, the longitudinal fields and the transverse ones are underestimated. The differences in the zero-iterative longitudinal and transverse field are relatively minor, but differences in the multi-iterative longitudinal and transverse field can be as high as 260 G and 500 G. The improved calibration method is more sensitive to the seeing effects variation than the former calibration method.     

Measurement of mesopause temperature from hydroxyl nightglow at Kolhapur (16.8°N, 74.2°E), India

The paper reports the nightglow observations of hydroxyl (8–3), (7–2) and (6–2) Meinel band carried out at a low latitude station Kolhapur (16.8°N, 74.2°E, dip latitude 10.6°N), India during November 2002 to May 2005 with the objective of investigating mesopause dynamics based on derived OH rotational temperature. Overall, 132 nights of quality data were collected using filter-tilting photometer and an all-sky scanning photometer. The mean mesopause temperature observed at Kolhapur is 195 ± 11, 196 ± 9 and 195 ± 7 K from OH (8–3), (7–2) and (6–2) band emissions, respectively, using transition probabilities given by Langhoff et al. [Langhoff, S.R., Werner, H.J., Rosmus, P. Theoretical transition probabilities for the OH Meinel system. Journal of Molecular Spectroscopy 118, 507–529, 1986]. Small wave-like variations (periodicities ∼ few hours) existing over long period variations in derived temperatures are also present. A steady decrease of emission intensities from evening to dawn hours has been observed in approximately 59% of nights. No significant change of nightly mean temperatures has been noted. Furthermore, about 62% of observed nightly mean temperatures lie within one error bar of MSISE-90 model predictions.     

Response of the mesopause temperatures to solar activity over Yakutia in 1999–2013

OH(6-2) rotational temperature trends and solar cycle effects are studied. Observations were carried out at the Maimaga station (63.04°N, 129.51°E) for the period August 1999 to March 2013. Measurements were conducted with an infrared spectrograph. Temperatures were determined from intensity ratios in the P branch of the OH band. The monthly average residuals of temperature after the subtraction of the mean seasonal variation were used for a search for the solar component of temperature response. The dependence of temperatures on solar activity has been investigated using the Ottawa 10.7 cm flux as a proxy. A linear regression fitting on residual temperatures yields a solar cycle coefficient of 4.24 ± 1.39 K/100 solar flux units (SFU). The cross-correlation analyses showed that changes of the residual temperature follow changes of solar activity with a quasi-two year delay (25 months). The temperature response at the delay of 25 months reaches 7 K/100 SFU. The possible reason of the observed delay can be an influence of quasi-biennial oscillations (QBO) of the atmosphere on the relation of temperature and solar activity. The value of the temperature trend after the subtraction of seasonal and solar components is not statistically significant.     

Effect of the troposphere on surface neutron counter measurements

Surface neutron counter data are often used as a proxy for atmospheric ionisation from cosmic rays in studies of extraterrestrial effects on climate. Neutron counter instrumentation was developed in the 1950s and relationships between neutron counts, ionisation and meteorological conditions were investigated thoroughly using the techniques available at the time; the analysis can now be extended using modern data. Whilst surface neutron counts are shown to be a good proxy for ionisation rate, the usual meteorological correction applied to surface neutron measurements, using surface atmospheric pressure, does not completely compensate for tropospheric effects on neutron data. Residual correlations remain between neutron counts, atmospheric pressure and geopotential height, obtained from meteorological reanalysis data. These correlations may be caused by variations in the height and temperature of the atmospheric layer at ∼100 hPa. This is where the primary cosmic rays interact with atmospheric air, producing a cascade of secondary ionising particles.     

Response of low latitude D-region ionosphere to the total solar eclipse of 22 July 2009 deduced from ELF/VLF analysis

Response of the D-region of the ionosphere to the total solar eclipse of 22 July 2009 at low latitude, Varanasi (Geog. lat., 25.27° N; Geog. long., 82.98° E; Geomag. lat. = 14° 55’ N) was investigated using ELF/VLF radio signal. Tweeks, a naturally occurring VLF signal and radio signals from various VLF navigational transmitters are first time used simultaneously to study the effect of total solar eclipse (TSE). Tweeks occurrence is a nighttime phenomena but the obscuration of solar disc during TSE in early morning leads to tweek occurrence. The changes in D-region ionospheric VLF reflection heights (h) and electron density (n_e: 22.6–24.6 cm⁻³) during eclipse have been estimated from tweek analysis. The reflection height increased from ∼89 km from the first occurrence of tweek to about ∼93 km at the totality and then decreased to ∼88 km at the end of the eclipse, suggesting significant increase in tweek reflection height of about 5.5 km during the eclipse. The reflection heights at the time of totality during TSE are found to be less by 2–3 km as compared to the usual nighttime tweek reflection heights. This is due to partial nighttime condition created by TSE. A significant increase of 3 dB in the strength of the amplitude of VLF signal of 22.2 kHz transmitted from JJI-Japan is observed around the time of the total solar eclipse (TSE) as compared to a normal day. The modeled electron density height profile of the lower ionosphere depicts linear variation in the electron density with respect to solar radiation as observed by tweek analysis also. These low latitude ionospheric perturbations on the eclipse day are discussed and compared with other normal days.     

Connection between winter mesopause region temperatures and diurnal LF reflection height variations measured at Collm

Scale height, H, estimates are calculated from the decrease/increase of ionospheric virtual reflection heights of low-frequency (LF) radio waves at oblique incidence in suitably defined morning intervals around sunrise during winter months. The day-to-day variations of H qualitatively agree with daily mean temperature variations around 90 km from meteor radar measurements. Since mesospheric long-period temperature variations are generally accepted to be the signature of atmospheric planetary waves, this shows that LF reflection height measurements can be used for monitoring the dynamics of the upper middle atmosphere. The long-term variations of monthly mean H estimates have also been analysed. There is no significant trend, which is in agreement with other measurements of mesopause region temperature trends.     

Effects on winter circulation of short and long term solar wind changes

Indices of the North Atlantic Oscillation and the Arctic Oscillation show correlations on the day-to-day timescale with the solar wind speed (SWS). Minima in the indices were found on days of SWS minima during years of high stratospheric aerosol loading. The spatial distribution of surface pressure changes during 1963–2011 with day-to-day changes in SWS shows a pattern resembling the NAO. Such a pattern was noted for year-to-year variations by Boberg and Lundstedt (2002), who compared NAO variations with the geo-effective solar wind electric field (the monthly average SWS multiplied by the average southward component, i.e., negative Bz component, of the interplanetary magnetic field). The spatial distribution of the correlations of geopotential height changes in the troposphere and stratosphere with the SWS; the geo-effective electric field (SWS∗Bz); and the solar 10.7 cm flux suggests that solar wind inputs connected to the troposphere via the global electric circuit, together with solar ultraviolet irradiance acting on the stratosphere, affect regional atmospheric dynamics.     

Physics results with the ARGO-YBJ experiment

ARGO-YBJ is a multipurpose experiment consisting in a dense sampling air shower array with 93% sensitive area located at very high altitude. The apparatus is in stable data taking since November 2007 at the YangBaJing Cosmic Ray Laboratory (Tibet, PR China, 4300 m a.s.l., 606 g/cm²). In this paper we report the main results in Gamma-Ray Astronomy and Cosmic Ray Physics after about 3 years of operation.     

Is there an abnormal enhancement of atmospheric aerosol before the 2008 Wenchuan earthquake?

There are extensive reports of ionospheric disturbances before the great 2008 Wenchuan earthquake, which are possibly explained by seismogenic electric field hypotheses linked with the aerosols injected in atmosphere. This paper attempts to investigate the possible change of atmospheric aerosol optical depth (AOD) associated with this earthquake by using MODIS data from both Terra and Aqua satellites. The result shows a clear enhancement of AOD along the Longmenshan faults 7 days before the quake, which is 1 day and 4 days earlier than the reported negative and positive ionospheric disturbances, respectively, and is 1 day earlier than or quasi-synchronism with other reported atmospheric anomalies including air temperature, outgoing longwave radiation and relative humidity. Particularly, the spatial distribution of AOD enhancement is very local and it is correlated well with the active faults and surface ruptures. We suggest that this unique enhancement could be associated with the Lithosphere–Atmosphere–Ionosphere coupling process during the preparation of the Wenchuan earthquake.     

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.     

Ozone trends in the mid-latitude stratopause region based on microwave measurements at Lindau (51.66° N, 10.13° E), the ozone reference model,and model calculations

We compared 8 years of ozone measurements taken at Lindau (51.66° N, 10.13° E) at altitudes between 40 and 60 km using the microwave technique with the CIRA ozone reference model that was established 20 years ago (Keating et al., 1990). We observed a remarkable decrease in ozone density in the stratopause region (i.e., an altitude of 50 km), but the decrease in ozone density in the middle mesosphere (i.e., up to 60 km in altitude) is slight. Likewise, we observed only a moderate decrease in the atmospheric region below the stratopause. Other studies have found the strongest ozone decrease at 40 km and a more moderate decrease at 50 km, which is somewhat in contradiction to our results. This decrease in ozone density also strongly depends on the season. Similar results showed model calculations using the GCM COMMA-IAP when considering the increase in methane. In the lower mesosphere/stratopause region, the strongest impact on the concentration of odd oxygen (i.e., O₃ and O) was observed due to a catalytic cycle that destroys odd oxygen, including atomic oxygen and hydrogen radicals. The hydrogen radicals mainly result from an increase in water vapor with the growing anthropogenic release of methane. The finding suggesting that the stratopause region is apparently attacked more strongly by the water vapor increase has been interpreted in terms of the action of this catalytic cycle, which is most effective near the stratopause and amplified by a positive feedback between the ozone column density and the ozone dissociation rate, thereby chemically influencing the ozone density. However, the rising carbon dioxide concentration cools the middle atmosphere, thereby damping the ozone decline by hydrogen radicals.     

Analysis of precipitable water vapor from GPS measurements in Chengdu region: Distribution and evolution characteristics in autumn

The rainfall process of Chengdu region in autumn has obvious regional features. Especially, the night-time rain rate of this region in this season is very high in China. Studying the spatial distribution and temporal variation of regional atmospheric precipitable water vapor (PWV) is important for our understanding of water vapor related processes, such as rainfall, evaporation, convective activity, among others in this area. Since GPS detection technology has the unique characteristics, such as all-weather, high accuracy, high spatial and temporal resolution as well as low cost, tracking and monitoring techniques on water vapor has achieved rapid developments in recent years. With GPS–PWV data at 30-min interval gathered from six GPS observational stations in Chengdu region in two autumns (September 2007–December 2007 and September 2008–December 2008), it is revealed that negative correlations exist between seasonally averaged value of GPS–PWV as well as its variation amplitude and local terrain altitude. The variation of PWV in the upper atmosphere of this region results from the water vapor variation from surface to 850 hPa. With the help of Fast Fourier Transform (FFT), it is found that the autumn PWV in Chengdu region has a multi-scale feature, which includes a seasonal cycle, 22.5 days period (quasi-tri-weekly oscillation). The variation of the GPS–PWV is related to periodical change in the transmitting of the water vapor caused by zonal and meridional wind strengths’ change and to the East Asian monsoon system. According to seasonal variation characteristics, we concluded that the middle October is the critical turning point in PWV content. On a shorter time scale, the relationship between autumn PWV and ground meteorological elements was obtained using the composite analysis approach.     

Electric discharge in the Martian atmosphere,Paschen curves and implications for future missions

Electric discharge between two electrically charged surfaces occurs at a well-defined, gas-dependent combination of atmospheric pressure and the distance between those two surfaces, as described by Paschen’s law. The understanding of when the discharge will occur in the conditions present on Mars is essential for designing space-flight hardware that will operate on the Martian surface as well as understanding electrical discharge processes occurring in the Martian atmosphere. Here, we present experimentally measured Paschen curves for a gas mixture representative of the Martian atmosphere and compare our results to breakdown voltages of carbon dioxide, nitrogen, and helium as measured with our system and from the literature. We will discuss possible implications for instrument development as well as implications for processes in the Martian atmosphere. The DC voltage at which electric discharge occurred between two stainless steel spheres was measured at pressures from 10⁻² to 100 torr in all gases. We measured a minimum voltage for discharge in the Mars ambient atmosphere of 410 ± 10 V at 0.3 torr cm. As an application, the breakdown properties of space-qualified, electrical wires to be used in the Sample Analysis at Mars (SAM) instrument suite on the Mars Science Laboratory (MSL) were studied.