Studies of electron and ion temperatures at 500 km altitude during sunrise using Indian SROSS C2 satellite

Solar dependence of electron and ion temperatures (T_e and T_i) in the ionosphere is studied using RPA data onboard SROSS C2 at an altitude of ∼500 km and 77°E longitude during early morning hours (04:00–07:00 LT) for three solar activities: solar minimum, moderate and maximum during winter, summer and equinox months in 10°S–20°N geomagnetic latitude. In winter the morning overshoot phenomenon is observed around 06:00 LT (T_e enhances to ∼4000 K) during low-solar activity and to T_e ∼ 3800 K, during higher solar activity. In summer, it is observed around 05:30 LT, but the rate of T_e enhancement is higher during moderate solar activity (∼2700 K/hr) than the low-solar activity (∼1700 K/hr). During equinox, this phenomenon is delayed and is observed around 06:00 LT (∼4200 K) during all three activities.     

The asymmetrical features in electron density during extreme solar minimum

The variations of plasma density in topside ionosphere during 23rd/24th solar cycle minimum attract more attentions in recently years. In this analysis, we use the data of electron density (Ne) from DEMETER (Detection of Electromagnetic Emissions Transmitted from Earthquake Regions) satellite at the altitude of 660–710 km to investigate the solstitial and equinoctial asymmetry under geomagnetic coordinate system at LT (local time) 1030 and 2230 during 2005–2010, especially in solar minimum years of 2008–2009. The results reveal that ΔNe (December–June) is always positive over Southern Hemisphere and negative over northern part whatever at LT 1030 or 2230, only at 0–10°N the winter anomaly occurs with ΔNe (December–June) > 0, and its amplitude becomes smaller with the declining of solar flux from 2005 to 2009. The ΔNe between September and March is completely negative during 2005–2008, but in 2009, it turns to be positive at latitudes of 20°S–40°N at LT 1030 and 10°S–20°N at LT 2230. Furthermore, the solstitial and equinoctial asymmetry index (AI) are calculated and studied respectively, which all depends on local time, latitude and longitude. The notable differences occur at higher latitudes in solar minimum year of 2009 with those in 2005–2008. The equinoctial AI at LT 2230 is quite consistent with the variational trend of solar flux with the lowest absolute AI occurring in 2009, the extreme solar minimum, but the solstitial AI exhibits abnormal enhancement during 2008 and 2009 with bigger AI than those in 2005–2007. Compared with the neutral compositions at 500 km altitude, it illustrates that [O/N₂] and [O] play some roles in daytime and nighttime asymmetry of Ne at topside ionosphere.     

Variability of the ionosphere over Irkutsk at low solar activity

Statistical and spectral analyses are performed to investigate variations of two ionosphere F2 layer key parameters, the critical frequency (foF2) and the peak height (hmF2), that were measured over Irkutsk (52.5°N, 104.0°E) from December 2006 to January 2008 under solar minimum. The analyses showed that both parameters contain quasi-harmonic oscillations with periods of T_n = 24/n hours (n = 1–7), among which the diurnal (n = 1) and semidiurnal (n = 2) ones are the strongest. Seasonal variations are explored of mean and median values, spectrum, amplitude, and phase of the diurnal and semidiurnal components of foF2 and hmF2.     

Contribution of satellite laser ranging to combined gravity field models

In the framework of satellite-only gravity field modeling, satellite laser ranging (SLR) data is typically exploited to recover long-wavelength features. This contribution provides a detailed discussion of the SLR component of GOCO02S, the latest release of combined models within the GOCO series. Over a period of five years (January 2006 to December 2010), observations to LAGEOS-1, LAGEOS-2, Ajisai, Stella, and Starlette were analyzed. We conducted a series of closed-loop simulations and found that estimating monthly sets of spherical harmonic coefficients beyond degree five leads to exceedingly ill-posed normal equation systems. Therefore, we adopted degree five as the spectral resolution for real data analysis. We compared our monthly coefficient estimates of degree two with SLR and Gravity Recovery and Climate Experiment (GRACE) time series provided by the Center for Space Research (CSR) at Austin, Texas. Significant deviations in C₂₀ were noted between SLR and GRACE; the agreement is better for the non-zonal coefficients. Fitting sinusoids together with a linear trend to our C₂₀ time series yielded a rate of (−1.75 ± 0.6) × 10⁻¹¹/yr; this drift is equivalent to a geoid change from pole to equator of 0.35 ± 0.12 mm/yr or an apparent Greenland mass loss of 178.5 ± 61.2 km³/yr. The mean of all monthly solutions, averaged over the five-year period, served as input for the satellite-only model GOCO02S. The contribution of SLR to the combined gravity field model is highest for C₂₀, and hence is essential for the determination of the Earth’s oblateness.     

Latitudinal variation of the topside electron temperature at different levels of solar activity

A database of electron temperature (T_e) measurements comprising of most of the available satellite measurements in the topside ionosphere is used for studying the solar activity variations of the electron temperature T_e at different latitudes, altitudes, local times and seasons. The T_e data are grouped into three levels of solar activity (low, medium, high) at four altitude ranges, for day and night, and for equinox and solstices. We find that in general T_e changes with solar activity are small and comparable in magnitude with seasonal changes but much smaller than the changes with altitude, latitude, and from day to night. In all cases, except at low altitude during daytime, T_e increases with increasing solar activity. But this increase is not linear as assumed in most empirical T_e models but requires at least a parabolic approximation. At 550 km during daytime negative as well as positive correlation is found with solar activity. Our global data base allows to quantify the latitude range and seasonal conditions for which these correlations occur. A negative correlation with solar activity is found in the invdip latitude range from 20 to 55 degrees during equinox and from 20 degrees onward during winter. In the low latitude (20 to −20 degrees invdip) F-region there is almost no change with solar activity during solstice and a positive correlation during equinox. A positive correlation is also observed during summer from 30 degrees onward.     

Heliotropic dust rings for Earth climate engineering

This paper examines the concept of a Sun-pointing elliptical Earth ring comprised of dust grains to offset global warming. A new family of non-Keplerian periodic orbits, under the effects of solar radiation pressure and the Earth’s J₂ oblateness perturbation, is used to increase the lifetime of the passive cloud of particles and, thus, increase the efficiency of this geoengineering strategy. An analytical model is used to predict the orbit evolution of the dust ring due to solar-radiation pressure and the J₂ effect. The attenuation of the solar radiation can then be calculated from the ring model. In comparison to circular orbits, eccentric orbits yield a more stable environment for small grain sizes and therefore achieve higher efficiencies when the orbit decay of the material is considered. Moreover, the novel orbital dynamics experienced by high area-to-mass ratio objects, influenced by solar radiation pressure and the J₂ effect, ensure the ring will maintain a permanent heliotropic shape, with dust spending the largest portion of time on the Sun facing side of the orbit. It is envisaged that small dust grains can be released from a circular generator orbit with an initial impulse to enter an eccentric orbit with Sun-facing apogee. Finally, a lowest estimate of 1 × 10¹² kg of material is computed as the total mass required to offset the effects of global warming.     

Solar cycle effect on temporal and spatial variation of topside ion density measured by SROSS C2 and ROCSAT 1 over the Indian longitude sector

We investigated the diurnal, seasonal and latitudinal variations of ion density Ni over the Indian low and equatorial topside ionosphere within 17.5°S to 17.5°N magnetic latitudes by combining the data from SROSS C2 and ROCSAT 1 for the 9 year period from 1995 to 2003 during solar cycle 23. The diurnal maximum density is found in the local noon or in the afternoon hours and the minimum occurs in the pre sunrise hours. The density is higher during the equinoxes as compared to that in the June and December solstice. The local time spread of the daytime maximum ion density increases with increase in solar activity. A north south asymmetry with higher ion density over northern hemisphere in the June solstice and over southern hemisphere in December solstice has been observed in moderate and high solar activity years. The crest to crest distance increases with increase in solar flux. Ion density bears a nonlinear relationship with F_10.7 cm solar flux and EUV flux in general. The density increases linearly with solar flux up to ∼150 sfu (1 sfu = 10⁻²²Wm⁻²Hz⁻¹) and EUV flux up to ∼50 units (10⁹ photons cm⁻² s⁻¹). But beyond this the density saturates. Inverse saturation and linear relationship have been observed in some season or latitude also. Inter-comparison of the three solar activity indices F_10.7 cm flux, EUV flux and F_10.7P (= (F_10.7 + F_10.7A)/2, where F_10.7A is the 81 day running average value of F_10.7) shows that the ion density correlates better with F_10.7P and F_10.7 cm fluxes. The annual average daytime total ion density from 1995 to 2003 follows a hysteresis loop as the solar cycle reverses. The ion density at 500 km over the Indian longitude sector as obtained by the international reference ionosphere is in general lower than the measured densities during moderate and high solar activity years. In low solar activity years the model densities are equal or higher than measured densities. The IRI EIA peaks are symmetric (±10°) in equinox while densities are higher at 10°N in June solstice and at 10°S in the December solstice. The model density follows F_10.7 linearly up to about F_10.7 > ∼150 sfu and then saturates.     

Validation of B2bot in the NeQuick model during high solar activity at Hainan station

NeQuick ionospheric electron density model, which has been developed to version 2, produces the full electron density profile in the ionosphere. Each part of the profile is modeled using Epstein layer formalism. Simple empirical relations are used to compute the thicknesses of each layer. In order to validate the B2_bot parameter in the NeQuick model during high solar activity, we use the data at Hainan, China (109.1°E, 19.5°N; Geomagnetic coordinates: 178.95°E, 8.1°N), measured with DPS-4, and study the diurnal and seasonal variations of B2_bot, ΔB2 (B2_best − B2_NeQuick 2) and the seasonal median values of B2_best/B2_NeQuick 2 at that region. The results show that, (1) The differences between B2_best and B2_NeQuick 2 have diurnal and seasonal variations. (2) The diurnal variations of B2_NeQuick 2 are smaller than those of B2_best. (3) Generally, except for early morning the experimental values are properly reproduced. (4) Generally, during morning the NeQuick model has an underestimation. The magnitude of underestimation varies with LT and season.     

Heating heavy ions in the polar corona by collisionless shocks: A one-dimensional simulation

Recently a new model for explaining the observations of preferential heating of heavy ions in the polar solar corona was proposed ( and ). In that model the ion energization mechanism is the ion reflection off supercritical quasi-perpendicular collisionless shocks in the corona and the subsequent acceleration by the motional electric field E = −V × B/c. The mechanism of heavy ion reflection is based on ion gyration in the magnetic overshoot of the shock. The acceleration due to the motional electric field is perpendicular to the magnetic field, giving rise to large temperature anisotropy with T_⊥ ? T_∥, in agreement with SoHO observations. Such a model is tested here by means of a one dimensional test particle simulation where ions are launched toward electric and magnetic profiles representing the shock transition. We study the dynamics of O⁵⁺, as representative of coronal heavy ions for Alfvénic Mach numbers of 2–4, as appropriate to solar corona. It is found that O⁵⁺ ions are easily reflected and gain more than mass proportional energy with respect to protons.     

Comparing IRI and a regional model with ionosonde measurements in Pakistan

We have used the technique of expansion in Empirical Orthogonal Functions (EOFs) to develop regional models of the critical frequencies of E and F₂ layers (f_oE, f_oF₂), peak height (h_mF₂), and semi-thickness of F₂ layer (Y_mF₂) over Pakistan. In the present study levels of solar activity specified by Smoothed Sunspot Number (R) from 10 to 200 are taken into account. The magnetic dip angle for the model ranges from 30° to 60°. We have compared the regional model and the International Reference Ionosphere (IRI) with measurements of three ionosondes in Pakistan. The model parameters f_oE and f_oF₂ are found overall comparable to the observed hourly median values during daytime at Karachi (geographic latitude = 24.95°N, longitude = 67.13°E, magnetic inclination = 37°), Multan (30.18°N, 71.48°E, 45°) and Islamabad (33.75°N, 73.13°E, 51.5°) during the years 1988, 1996 and 2000. For h_mF₂ the computed values by regional and IRI model for the year 1995 are found close to each other. However, for Y_mF₂the results are better during daytime as compared to nighttime.     

Occurrence of sporadic-E layer over the ionospheric station of Rome: Analysis of data for thirty-two years

Hourly systematic measurements of the highest frequency reflected by the sporadic-E layer (foEs) recorded at the Rome ionospheric observatory (Italy, 41.8° N, 12.5° E), were considered during the period January 1976–December 2007, to calculate the percentage of occurrence of sporadic-E layer with frequencies foEs greater than a given threshold value f_T, P(foEs > f_T).     

Quantifying the relationship between the measurement precision and specifications of a UV/visible sensor on a geostationary satellite

To investigate the feasibility of new satellite observations, including air quality (AQ) observations from geostationary (GEO) orbit, it is essential to link the measurement precision (ε) with sensor specifications in advance. The present study attempts to formulate the linkage between ε and specifications of a UV/visible sensor (signal-to-noise ratio (SNR), full width at half maximum (FWHM) of the slit function, and sampling ratio (SR)) on a GEO satellite. A sophisticated radiative transfer model (JACOSPAR) is used to calculate synthetic radiance spectra that would be measured by a UV/visible sensor observing the atmosphere over Tokyo (35.7°N, 139.7°E) from GEO orbit at 120°E longitude. The spectra, modified according to given sensor specifications, are analyzed by the differential optical absorption spectroscopy technique to estimate the ε for slant column densities of O₃ and NO₂. We find clear relationships: for example, the ε of the O₃ slant column density (molecules cm⁻²) and SNR at 330 nm are linked by the equation log(ε) = −1.06 · log(SNR) + 20.71 in the UV region, and the ε of the NO₂ slant column density and SNR at 450 nm are linked by log(ε) = −0.98 · log(SNR) + 18.00, at a FWHM = 0.6 nm (for the Gaussian slit function) and SR = 4. The relationships are mostly independent of other specifications (e.g., horizontal and temporal resolutions), as they affect ε primarily through SNR, providing constraints in determining the optimal SNR (and alternatively FWHM and SR) for similar UV/visible sensors dedicated for AQ studies.     

Equatorial vertical E × B drift velocities inferred from ionosonde measurements over Ouagadougou and the IRI-2007 vertical ion drift model

F-region vertical plasma drift velocities were deduced from the hourly hmF2 values acquired from ionogram data over a near dip equatorial station Ouagadougou (12.4°N, 358.5°E, dip angle 5.9°N) in Africa. Our results are compared against the global empirical model of Scherliess and Fejer (1999) incorporated in the IRI model (IRI-2007) for 1600 to 0800 LT from 1 year of data during sunspot maximum year of 1989 (yearly average solar flux intensity, F10.7 = 192) corresponding to the peak phase of solar cycle 22, under magnetically quiet conditions. The drifts are entirely downward between 2000 and 0500 LT bin for both techniques and the root mean square error (RMSE) between the modeled and the ionosonde vertical plasma drifts during these periods is 3.80, 4.37, and 4.74 m/s for June solstice, December solstice and equinox, respectively. Ouagadougou average vertical drifts show evening prereversal enhancement (PRE) velocity peaks (V_ZP) of about 16, 14, and 17 m/s in June solstice, December solstice, and equinox, respectively, at 1900–2000 LT; whereas global empirical model average drifts indicate V_ZP of approximately 33 m/s (June solstice), 29 m/s (December solstice), and 50 m/s (equinox) at 1800 LT. We find very weak and positive correlation (+0.10376) between modeled V_ZP versus F10.7, while ionosonde V_ZP against F10.7 gives worst and opposite correlation (−0.05799). The results also show that modeled V_ZP–Ap indicates good and positive correlation (+0.64289), but ionosonde V_ZP–Ap exhibits poor and negative correlation (−0.22477).     

Variation of ionospheric total electron content in Taiwan region of the equatorial anomaly from 1994 to 2003

The ionospheric total electron content (TEC) in the northern hemispheric equatorial ionospheric anomaly (EIA) region is studied by analyzing dual-frequency signals of the Global Position System (GPS) acquired from a chain of nine observational sites clustered around Taiwan (21.9–26.2°N, 118.4–112.6°E). In this study, we present results from a statistical study of seasonal and geomagnetic effects on the EIA during solar cycle 23: 1994–2003. It is found that TEC at equatorial anomaly crests yield their maximum values during the vernal and autumnal months and their minimum values during the summer (except 1998). Using monthly averaged I_c (magnitude of TEC at the northern anomaly crest), semi-annual variations is seen clearly with two maxima occurring in both spring and autumn. In addition, I_c is found to be greater in winter than in summer. Statistically monthly values of I_c were poorly correlated with the monthly Dst index (r = −0.22) but were well correlated with the solar emission F_10.7 index (r = 0.87) for the entire database for the period during 1994–2003. In contrast, monthly values of I_c were correlated better with Dst (r ? 0.72) than with F_10.7 (r ? 0.56) in every year during the low solar activity period (1994–1997). It suggests that the effect of solar activity on I_c is a longer term (years), whereas the effect of geomagnetic activity on I_c is a shorter term (months).     

Correction of B2bot for NeQuick during low solar activity at Hainan station

This work is a continuation of the previous article and it focuses on low solar activity and modeling effort. NeQuick model uses Epstein layer formalism to model each part of the profile. We study the diurnal and seasonal variations of B2_bot, ΔB2 (B2_best − B2_NeQuick2) and R (B2_best/B2_{NeQuick 2}) at Hainan station during low solar activity. The results show it is possible to improve the B2_bot parameter of the NeQuick model at that region during low solar activity. Then, we use a function ?(t) with LT in different seasons to correct the B2_bot formula of NeQuick 2. The correction shows that (1) By the correction formula, the B2_bot of NeQuick is improved. The maximum standard deviation is improved for 9 km. (2) The correction formula is more effective in summer than in equinox and winter and performs better during early morning hours than during the rest of the day.     

Validation of space-born Moderate Resolution Imaging Spectroradiometer remote sensors aerosol products using application of ground-based Multi-wavelength Radiometer

The measurements of aerosol optical properties were carried out during April 2006 to March 2011 over Mohal (31.9°N, 77.12°E) in the northwestern Indian Himalaya, using the application of ground-based Multi-wavelength Radiometer (MWR) and space-born Moderate Resolution Imaging Spectroradiometer (MODIS) remote sensors. The average (±standard deviation) values of aerosol optical depth (AOD) at 500 nm, Ångström exponent and turbidity coefficient during the entire measurement period were 0.25 ± 0.09, 1.15 ± 0.42 and 0.12 ± 0.06 respectively. About 86% AOD values retrieved from MODIS remote sensor were found within an uncertainty limit (Δτ_pλ = ±0.05 ± 0.15τ_pλ). In general, the MWR derived AOD values were higher than that of MODIS retrieval with absolute difference ∼0.02. During the entire period of measurement space-born MODIS remote sensor and ground-based MWR observation showed good correspondence with significant correlation coefficient ∼0.78 and root mean square difference ∼0.06. For daily observations the relative difference between these two estimates stood less than 9%. However, satellite-based and ground-based observation showed good correspondence, but further efforts still needed to eliminate systematic errors in the existing MODIS algorithm.     

Validation of B0 and B1 in the IRI 2001 model at low solar activity for Ilorin an equatorial station

A new set of data obtained at low solar activity from Ilorin, Nigeria (geog. latitude 8.5°N, geog longitude, 4.6°E, dip 4.1°S) is used to validate the IRI 2001 model at low solar activity. The results show in general a good agreement between model and observed B0 at night but an over estimation during daytime. The overestimation is greatest during the morning period (0600LT–1000LT). The model prediction for B1 is fairly good at night and during the day. A dependence of B0 on solar zenith angle χ is observed during the daytime. A formulation of the form B0 = A[cos(χ)ⁿ] is therefore proposed. Values of the constants n and A were determined for the period of low solar activity for this station.     

Polarization pattern of low and mid-frequency magnetic pulsations in the polar cap: A comprehensive analysis at Terra Nova Bay (Antarctica)

The polarization pattern of ULF pulsations (f ≈ 1–100 mHz) at Terra Nova Bay (Antarctica, CGM λ ∼ 80°) has been determined for the entire 2003, soon after the solar maximum. A comparison with the results of previous investigations, conducted at the same station close to the solar minimum (1994–96), allows to focus common elements and major differences among different frequency bands which persist through the entire solar cycle. Basically, between f ∼ 1.5 and 5 mHz, the day can be divided into four sectors with alternate polarizations. The local time and latitudinal dependence of the observed pattern can be tentatively interpreted in terms of a latitude of resonant field lines reaching λ ∼ 80° in the noon sector; on the other hand, resonance effects of lower latitude field lines can be clearly identified also far from the noon meridian when the station moves into the deep polar cap. Moreover, in the morning sector the resonance region would extend to lower latitudes than in the evening sector. The proposed profile of the resonant region can interpret also the results obtained at other cusp/auroral stations and appears consistent with that one inferred in the northern hemisphere at smaller latitudes. The resonance region progressively shifts toward lower latitude with increasing frequency; correspondingly, the four-sector pattern progressively disappears at TNB. Above f ∼ 20 mHz, the experimental observations might suggest an additional contribution from Sunward propagating waves, possibly via the magnetotail lobes.     

Changes in sea-level pressure over South Korea associated with high-speed solar wind events

We explore a possibility that the daily sea-level pressure (SLP) over South Korea responds to the high-speed solar wind event. This is of interest in two aspects: first, if there is a statistical association this can be another piece of evidence showing that various meteorological observables indeed respond to variations in the interplanetary environment. Second, this can be a very crucial observational constraint since most models proposed so far are expected to preferentially work in higher latitude regions than the low latitude region studied here. We have examined daily solar wind speed V, daily SLP difference ΔSLP, and daily log(BV²) using the superposed epoch analysis in which the key date is set such that the daily solar wind speed exceeds 800 km s⁻¹. We find that the daily ΔSLP averaged out of 12 events reaches its peak at day +1 and gradually decreases back to its normal level. The amount of positive deviation of ΔSLP is +2.5 hPa. The duration of deviation is a few days. We also find that ΔSLP is well correlated with both the speed of solar wind and log(BV²). The obtained linear correlation coefficients and chance probabilities with one-day lag for two cases are r ? 0.81 with P > 99.9%, and r ? 0.84 with P > 99.9%, respectively. We conclude by briefly discussing future direction to pursue.     

Interpretation of the SOHO/UVCS observations of two CME-driven shocks

We report on the analysis of two fast CME-driven shocks observed with the UltraViolet Coronagraph Spectrometer (UVCS) on board the Solar and Heliospheric Observatory (SOHO). The first event, detected on 2002 March 22 at 4.1 R_⊙ with the UVCS slit placed in correspondence with the flank of the expanding CME surface, represents the highest UV detection of a shock obtained so far with the UVCS instrument in the corona. The second one, detected on 2002 July 23 at 1.6 R_⊙ with the UVCS slit placed in correspondence with the front of the expanding CME surface, shows an anomalous deficiency of ion heating with respect to what observed in previous CME/shocks observed by UVCS, possibly reflecting the effect of different coronal plasma conditions over the solar cycle. From the two different sets of observations we derived an estimate for the shock compression ratio X, which turns out to be X = 2.4 ± 0.2 and X = 2.2 ± 0.1, respectively, for the first and second event. Comparison between the two events provides complementary perspectives on the dynamical evolution of CME-driven shocks.