An updated model of atomic oxygen redline dayglow emission

A comprehensive model is developed using the updated rate coefficients and transition probabilities to study the redline dayglow emission of atomic oxygen. The solar EUV fluxes are obtained from the Solar Irradiance Platform (SIP), and incorporated into the model successfully. All possible production and loss mechanisms of O(¹D) are considered in the model. The neutral number densities and temperature are adopted from the NRLMSISE-00 model. The ion and electron densities, and electron temperature are adopted from the IRI-07 model. The model results are validated with the help of measurements as provided by the Wind Imaging Interferometer (WINDII) on board Upper Atmosphere Research Satellite (UARS). The present results are found in better agreement with the measurements in comparison with the earlier model. The measured volume emission rate profiles are reproduced quite well by the present model. The model results show that the updated rate coefficients and transition probabilities are quite consistent and may be used in the aeronomical studies.     

Effect of severe geomagnetic storm conditions on atomic oxygen greenline dayglow emission in mesosphere

Severe geomagnetic storms and their effects on the 557.7 nm dayglow emission are studied in mesosphere. This study is primarily based on photochemical model with the necessary input obtained from a combination of experimental observations and empirical models. The model results are presented for a low latitude station Tirunelveli (8.7°N, 77.8°E). The volume emission rates are calculated using MSISE-90 and NRLMSISE-00 neutral atmospheric models. A comparison is made between the results obtained from these two models. A positive correlation amongst volume emission rate (VER), O, O₂ number densities and Dst index has been found. The present results indicate that the variation in emission rate is more for MSISE-90 than in NRLMSISE-00 model. The maximum depletion in the VER of greenline dayglow emission is found to be about 30% at 96 km during the main phase of the one of the geomagnetic storms investigated in the case of MSISE-90 (which is strongest with Dst index −216 nT). The O₂ density decreases about 22% at 96 km during the main phase of the same geomagnetic storm.The NRLSMSISE-00 model does not show any appreciable change in the number density of O during any of the two events. The present study also shows that the altitude of peak emission rate is unaffected by the geomagnetic storms. The effect of geomagnetic storm on the greenline nightglow emission has also been studied. It is found that almost no correlation can be established between the Dst index and variations in the volume emission rates using the NRLMSISE-00 neutral model atmosphere. However, a positive correlation is found in the case of MSISE-90 and the maximum depletion in the case of nightglow is about 40% for one of the storms. The present study shows that there are significant differences between the results obtained using MSISE-90 and NRLMSISE-00.     

Comparison between African equatorial station ground-based inferred vertical E × B drift,Jicamarca direct measured drift,and IRI model

The Incoherent Scatter Radar measurement over Jicamarca, together with the IRI model-2007 measurements were compared with ground-based digisonde inferred E × B drift over Ilorin in the African region during year of solar minima (F10.7 = 81). Seasonally, Ilorin pre-reversal enhancement (PRE) had peak drift velocities of 7.2, 3.7 and 7.9 m/s for March equinox, September equinox and December solstice respectively, while Jicamarca drifts indicated 13.0, 10.5 and 5.2 m/s; as well as the IRI model with 14.3, 8.4 and 0.7 m/s in similar order. PRE value was insignificant during June solstice. The PRE magnitude of the IRI-model during the equinoxes is twice the value obtained at Ilorin. The daytime E × B drift peaked over Ilorin 1–2 h earlier than both the modeled and Jicamarca observations. This could be due to the difference in sunset time at the conjugate points corresponding to the altitude of the observation. During the evening time PRE, the respective correlation coefficients (R) for Vz–F10.7 relation over Jicamarca, Ilorin and the modeled observations are −0.5559, 0.4796 and −0.4979. Similarly, the Vz–Ap relation exhibit excellent anti-correlation coefficient (R = −0.8637) for the IRI-model, −0.4827 over Jicamarca and 0.3479 for Ilorin. Annual mean drift velocities over Jicamarca, Ilorin and IRI model measurements respectively are 10, 5.6 and 10 m/s for the peak PRE observation; 15, 16 and 21 m/s for the daytime pre-sunrise peak values; and −21, −9 and −16 m/s for the nighttime downward reversals. The root-mean square (RMS) deviation between IRI-model and the Ilorin drift between 2000 and 0500 h is 4.37, 2.03, 3.71 and 2.42 m/s for March equinox, June solstice, September equinox and December solstice respectively. For Jicamarca–Ilorin drift relation, RMS deviation is 5.48, 2.30, 3.47 and 1.27 m/s in the same order respectively. Annual hmF2 inferred drift over Ilorin during daytime is higher by a factor of ≈2 and 3 at Jicamarca and IRI model measurements respectively; and by a factor of ≈5 for both during the night-time period. The limitations in using hmF2 to infer drifts are discussed.     

Occurrence rate of SAR arcs during the 23rd solar activity cycle

The occurrence rate of SAR arcs during 1997–2007 has been analyzed based on the photometric observations at the Yakutsk meridian (Maimaga station, corrected geomagnetic coordinates: 57°N, 200°E). SAR arcs appeared in 114 cases (∼500 h) during ∼370 nights of observations (∼3170 h). The occurrence frequency of SAR arcs increases to 27% during the growth phase of solar activity and has a clearly defined maximum at a decline of cycle 23. The SAR arc registration probability corresponds to the variations in geomagnetic activity in this solar cycle. The dates, intervals of UT, and geomagnetic latitudes of SAR arc observations at the Yakutsk meridian are presented.     

On the behavior of thermospheric O1D 630.0 nm dayglow emission during counter electrojet events

This study presents the response of thermospheric O¹D 630.0 nm dayglow emission to the variability associated with equatorial Counter Electrojet (CEJ) events. The analysis based on the data from a meridian scanning Dayglow Photometer, Digital Ionosonde and Proton Precession Magnetometer over Trivandrum (8.5°N, 77°E, 0.5°dip lat.), indicates that the O¹D 630.0 nm emission behave distinctly different during the CEJ events compared to that on normal days. It has been observed that O¹D 630.0 nm emission shows enhancement during the negative excursion of the ΔH, followed by an unusual depletion during the peak CEJ time. The observed variability was found to be more pronounced in a latitudinal region of ±3° centered at around the dip equator. In addition, the emission intensities also exhibit the presence of enhanced short period oscillations of periodicity 20–30 min during the CEJ events. Analysis of the data from the collocated ionosonde revealed that the F-region electron density showed enhancement during the early phase of the CEJ and a decrease during the peak CEJ. Further, the simulation studies using a Quasi 2 dimensional ionospheric model showed that the modified plasma fountain during the CEJ can alter the plasma density at the emission centroid. The study reveals a strong dynamical coupling between the E and F-region of the dip equatorial ionosphere.     

Large geomagnetic storms of extreme solar event periods in solar cycle 23

During extreme solar events such as big flares or/and energetic coronal mass ejections (CMEs) high energy particles are accelerated by the shocks formed in front of fast interplanetary coronal mass ejections (ICMEs). The ICMEs (and their sheaths) also give rise to large geomagnetic storms which have significant effects on the Earth’s environment and human life. Around 14 solar cosmic ray ground level enhancement (GLE) events in solar cycle 23 we examined the cosmic ray variation, solar wind speed, ions density, interplanetary magnetic field, and geomagnetic disturbance storm time index (D_st). We found that all but one of GLEs are always followed by a geomagnetic storm with D_st  −50 nT within 1–5 days later. Most(10/14) geomagnetic storms have D_st index  −100  nT therefore generally belong to strong geomagnetic storms. This suggests that GLE event prediction of geomagnetic storms is 93% for moderate storms and 71% for large storms when geomagnetic storms preceded by GLEs. All D_st depressions are associated with cosmic ray decreases which occur nearly simultaneously with geomagnetic storms. We also investigated the interplanetary plasma features. Most geomagnetic storm correspond significant periods of southward B_z and in close to 80% of the cases that the B_z was first northward then turning southward after storm sudden commencement (SSC). Plasma flow speed, ion number density and interplanetary plasma temperature near 1 AU also have a peak at interplanetary shock arrival. Solar cause and energetic particle signatures of large geomagnetic storms and a possible prediction scheme are discussed.     

Variations of the solar acoustic high-degree mode frequencies over solar cycle 23

Using full-disk observations obtained with the Michelson Doppler Imager (MDI) on board the Solar and Heliospheric Observatory (SOHO) spacecraft, we present variations of the solar acoustic mode frequencies caused by the solar activity cycle. High-degree (100 < ? < 900) solar acoustic modes were analyzed using global helioseismology analysis techniques over most of solar cycle 23. We followed the methodology described in details in [Korzennik, S.G., Rabello-Soares, M.C., Schou, J. On the determination of Michelson Doppler Imager high-degree mode frequencies. ApJ 602, 481–515, 2004] to infer unbiased estimates of high-degree mode parameters ([see also Rabello-Soares, M.C., Korzennik, S.G., Schou, J. High-degree mode frequencies: changes with solar cycle. ESA SP-624, 2006]). We have removed most of the known instrumental and observational effects that affect specifically high-degree modes. We show that the high-degree changes are in good agreement with the medium-degree results, except for years when the instrument was highly defocused. We analyzed and discuss the effect of defocusing on high-degree estimation. Our results for high-degree modes confirm that the frequency shift scaled by the relative mode inertia is a function of frequency and it is independent of degree.     

Study of the 28 October 2003 and 20 January 2005 solar flares by means of 2.223 MeV gamma-emission line

By the data on intensity-time profiles of the neutron capture line of 2.223 MeV we have studied some characteristics of two solar flares, 28 October 2003 and 20 January 2005 (INTEGRAL and CORONAS-F observations, respectively). The SINP code was applied making allowance for the main processes of neutron interactions and deceleration in the solar plasma, character of neutron source, losses of neutrons and density model of the solar atmosphere. Comparison of the computed time profiles of 2.223 MeV line with observed ones for the flare of 28 October 2003 confirms the results obtained earlier for three other flares. Namely, the effect of density enhancement (EDE) in the sub-flare region, as well as the variations (hardening) of accelerated particle spectrum in the course of the event have been confirmed. The usual modeling procedure by the SINP code, however, seems to be inapplicable to the event of 20 January 2005. Possible causes of density enhancements during some flares and peculiarities of the 20 January 2005 flare are discussed.     

The use of MODIS 250 m bands to improve the MODIS 1 km ocean color atmospheric correction algorithm in turbid water

“Clear water” is a scale-dependent concept, so it is more likely to successfully find the “clear water” from images with smaller scale than that with larger scale data. In this study, an optimal spectral relationship of moderate-resolution imaging spectroradiometer (MODIS) 250 m and 1 km resolution data at near-infrared bands (OSRLM) is constructed for converting pseudo “clear water” reflectance at 859 nm to those at 748 and 869 nm. According to scale effects, the satellite-observed pseudo “clear water” reflectance is greater than 5.18%, larger than that derived from OSRLM model. An atmospheric correction model for MODIS 1km data using pseudo “clear water” reflectance of MODIS 250 m data (ACMM) was developed for improving the performance of traditional “clear water” atmospheric correction model (CWAC). The model validation results indicate that ACMM model has a better performance than CWAC model. By comparison, the uncertainty decreases by 19.18% in the use of ACMM model over CWAC model for deriving water-leaving reflectance in Taihu Lake, China. This uncertainty is significantly reduced in water-leaving reflectance estimation due to partial removal of scale effects on “clear water”. These findings imply that satellite-derived aerosol scattering contribution at smaller scale usually has a better performance than that at larger scale.     

Evolutionary features of polar and non-polar coronal holes during solar cycle 24 and the rising phase of cycle 25

This paper presents the results of the analysis of the evolution of coronal holes (CHs) on the Sun during the period May 13, 2010 – March 20, 2022, covering Solar Cycle 24. Our study uses images in the extreme-ultraviolet iron line (Fe XII 193 Å) obtained with the Atmospheric Imager Assembly of the Solar Dynamics Observatory (AIA/SDO). To localize CHs and determine their areas, we used the Heliophysics Event Knowledgebase (HEK). We separate the CHs into polar and non-polar and study the evolutionary features of each group. During this period, an asymmetry between the Northern (N) and Southern (S) Hemispheres (N-S or hemispheric asymmetry) is detected both in the solar activity (SA) indices and in the localization of the maximum areas of the polar and non-polar CHs. It is shown that the hemispheric asymmetry of the areas of polar and non-polar CHs varies significantly over time and that the nature of these changes is clearly related to the SA cycle. We find that for most of the period, the polar CHs were predominated generated in the S- hemisphere while the non-polar CHs were dominant in the N- hemisphere. It is found that the maximum and minimum of the hemispheric imbalance in the areas of non-polar CHs are close in time to the maximum and minimum of the asymmetry of the SA indices (the number and areas of sunspots). The maximum hemispheric imbalance of the polar CH areas is observed at the maximum of Cycle 24, and the minimum imbalance is found at the cycle minimum. These results confirm our assumption that these two types of CHs are of a different nature and that the non-polar CHs, like sunspots, are elements of the general magnetic activity.     

Observations of non-migrating tides and ionospheric perturbations of O(D) airglow by ISUAL instrument

In this paper, we present the spatial variations of O(¹D) airglow observed by the ISUAL (Imager of Sprites and Upper Atmospheric Lightning) instrument on board the FORMOSAT-2 satellite. With a CCD camera and a 630 nm filter, ISUAL can measure global atmospheric emissions lying between the heights of 80 and 300 km. In days of 3–6 September 2008 and 25–27 February 2009, ISUAL has measured the emissions of O(¹D) airglow with results showing strong longitudinal peak-3 and peak-4 structures. The Lomb-Scargle analyses for these two cases show periods of longitudes of 120° and 90° supporting the DE2 and DE3 non-migrating tides. The 630 nm emissions are enhanced in equatorial regions and are lying along the equator. Over Africa its intensity can sometimes increase up to 80% relative to other longitudes. The perturbation is so strong that non-migrating tides are erased. A case of bimodal distribution with strong emissions at latitudes in equator and mid-latitude in geographic coordinates was observed.     

Statistical properties of the IMF clock angle in the solar wind with northward and southward interplanetary magnetic field based on ACE observation from 1998 to 2009: Dependence on the temporal scale of the solar wind

Utilizing ACE satellite observations from 1998 to 2009, we performed the elaborate study on the properties of the clock angle θ_CA (arctan(By/Bz) (?90° to 90°) of the interplanetary magnetic field (IMF) in the solar wind at 1?AU. The solar wind with northward IMF (NW-IMF) and southward IMF (SW-IMF) are analyzed, independently. Statistical analysis shows that the solar wind with SW-IMF and NW-IMF has similar properties in general, including their durations, the IMF Bz and By components, and the IMF θ_CA. Then, the solar wind with NW-IMF (SW-IMF) is classified into five different temporal scales according to the duration of the NW-IMF (SW-IMF), i.e., very-short wind of 10–30?min, short-scale wind of 0.5–1?h, moderate-scale wind of 1–3?h, long-scale wind of 3–5?h, and super-long wind >5?h. Our analysis reveals that the IMF θ_CA has a distinct decrease with increase of the temporal scale of the solar wind. Next, the solar wind is classified into two groups, i.e., the high-speed solar wind (>450?km/s) and the low-speed solar wind (<450?km/s). Our analysis indicates that the IMF θ_CA depends highly on the solar wind speed. Statistically, high-speed solar wind tends to have larger IMF θ_CA than low-speed solar wind. The evolutions of the solar wind and IMF with the solar activity are further studied, revealing no clear solar variation of the IMF θ_CA. Finally, we analyze the monthly variation of the IMF θ_CA. Superposed epoch result strongly suggests the seasonal variation of the IMF θ_CA.     

Reverberation time lags in the high luminosity quasar PG 1247+267

Virial estimates of the black hole mass in the center of AGNs, derived from single-epoch observations of luminosity and emission line widths, are now available for several thousands of objects at all redshifts and luminosities, so that studies of the cosmological evolution of the AGN mass function are becoming possible. These estimates are based on the empirical luminosity-size relation measured through reverberation mapping at low redshifts and luminosities. For this reason the spectrophotometric monitoring of 4 luminous quasars was started in 2003, with the 1.8 m telescope of the Asiago Observatory. To measure the reverberation time lag we adopt a method, recently introduced by Zu et al. (2011), which takes advantage of the available statistical information on variability autocorrelation and makes simple assumptions on the continuum-emission line transfer function. We present the estimate of the broad line region size for the quasar PG 1247+267, which is the most luminous object with reverberation measures to date. Virial products computed for C iv and C iii] emission lines provide consistent results. A possible flattening of the size-luminosity relation for C iv at high luminosity is suggested.     

Evidence for the lowering of the centroid of daytime thermospheric O(D) 630.0 nm emission over the magnetic equator: First results

It is shown in this paper for the first time that the intensity of the daytime thermospheric O(¹D) 630.0 nm airglow as measured by the ground-based dayglow photometer over Trivandrum (8.5°N; 77°E; dip lat. 0.5°N), a geomagnetic dip equatorial station, exhibit a direct correlation with the electron density at 180 km. This altitude is about ∼40 km lower than the believed centroid of the O(¹D) 630.0 nm dayglow emission i.e. 220 km. This observation is contrary to the understanding of the behavior of O(¹D) 630.0 nm dayglow over equatorial/low latitudes. Over these latitudes, the variations of the measured intensity of O(¹D) 630.0 nm dayglow are known to be associated with the changes in the electron density at altitudes around 220 km, the centroid of this emission. In this context, the present results indicating the lowering of the peak altitude of O(¹D) 630.0 nm emission from ∼220 to ∼180 km over the dip equator is new. Recent results on solar XUV flux indicate that this could be an important parameter that controls the O(¹D) 630.0 nm dayglow excitation rates through modulations in the neutral and ionic composition in lower thermosphere-ionosphere region. However, the lowering of the centroid of O(¹D) 630.0 nm emission, as shown in this study, has been ascribed primarily to the fountain effect associated with the equatorial ionization anomaly.     

Long-term changes in indices of geomagnetic activity at the auroral station Sodankylä

Here we compare the traditional analog measure of geomagnetic activity, Ak, with the more recent digital indices of IHV and Ah based on hourly mean data, and their derivatives at the auroral station Sodankylä. By this selection of indices we study the effects of (i) analog vs. digital technique, and (ii) full local-time vs. local night-time coverage on quantifying local geomagnetic activity. We find that all other indices are stronger than Ak during the low-activity cycles 15–16 suggesting an excess of very low scalings in Ak at this time. The full-day indices consistently depict stronger correlation with the interplanetary magnetic field strength, while the night-time indices have higher correlation with solar wind speed. The Ak index correlates better with the digital indices of full-day coverage than with any night-time index. However, Ak depicts somewhat higher activity levels than the digital full-day indices in the declining phase of the solar cycle, indicating that, due to their different sampling rates, the latter indices are less sensitive to high-frequency variations driven by the Alfvén waves in high-speed streams. On the other hand, the night-time indices have an even stronger response to solar wind speed than Ak. The results strongly indicate that at auroral latitudes, geomagnetic indices with different local time coverage reflect different current systems, which, by an appropriate choice of indices, allows studying the century-scale dynamics of these currents separately.     

VHF radio response of the near Earth space during solar activity growth in October, 2003

The analysis of observations of very high frequency radio noise intensity at the middle latitude on a frequency f = 500 MHz from 14th till 26th of October, 2003 is presented. These data are compared with the solar radio bursts in the range of frequencies 1–14 MHz registered by RAD2 receiver of the WAVES device installed on board the WIND spacecraft.     

Relativistic solar protons in the ground level event of 23 February 1956: New study

In this study we applied again to the outstanding solar particle event of 23 February 1956, the largest one in the entire history of observations of solar cosmic rays. Due to significant improvement of the analysis/modeling techniques and new understanding of physical processes in the solar atmosphere and interplanetary space, a possibility arises to interpret the old data in the light of modern concept of multiple particle acceleration at/near the Sun. In our new analysis the data of available then neutron monitors and muon telescopes are used. The technique of the analysis includes: (a) calculation of asymptotic cones of ground-based detectors; (b) modeling of cosmic ray detector responses at variable parameters of the flux of solar relativistic protons; (c) determination of primary solar proton parameters outside magnetosphere by comparison of computed responses with observations. Certain evidence was obtained that the flux of relativistic solar protons consisted of two distinct components: prompt and delayed ones. The prompt component with exponential energy spectrum caused a giant impulse-like increase at a number of European cosmic ray stations. The delayed component had a power-law spectrum and was a cause of gradual increase at cosmic ray stations in the North American region. A numerical simulation of the proton acceleration in the vicinity of the magnetic reconnection region brings to the proton spectrum with exponential dependence on energy. This agrees with observational data for the prompt component. It is also shown that the huge increase in ∼5000% on neutron monitors was due to the prompt component only with the exponential proton spectrum. The power-law spectrum of comparable intensity gave considerably smaller effect.     

An investigation of the solar cycle impact on the lower thermosphere O(S) nightglow emission as observed by WINDII/UARS

The yearly variation of the integrated emission rate of the O(¹S) nightglow in the lower thermosphere is studied and the solar cycle impact is examined from the observations of the Wind Imaging Interferometer (WINDII) operated on the Upper Atmosphere Research Satellite (UARS). More than 300,000 volume emission rate profiles of the O(¹S) nightglow observed by WINDII for 40°S–40°N latitudes during November 1991–August 1997 over half of a solar cycle are utilized. These profiles are vertically integrated for the altitude range of 80–100 km and the equivalent column integrated emission rates are then zonally averaged for bins with 10° latitude and 3 month intervals. It is found that for each latitude the O(¹S) nightglow emission rate appears to increase with increasing solar F10.7 cm flux, following a linear relationship. This characterizes the solar cycle impact on the O(¹S) nightglow, while the solar influence is modulated by a seasonal variation. Based on these variations, an empirical formula is derived for predicting the three-month averages of the O(¹S) nightglow integrated emission rate. The standard error of the estimated values from the formula is smaller than 30 Rayleigh.     

Estimation of Galactic Cosmic Ray exposure inside and outside the Earth’s magnetosphere during the recent solar minimum between solar cycles 23 and 24

The evidently low solar activity observed between solar cycles 23 and 24 during the years 2008–2010 led to a substantial increase in the Galactic Cosmic Ray (GCR) intensity in comparison with preceding solar minima. As the GCRs consist of highly-ionizing charged particles having the potential to cause biological damage, they are a subject of concern for manned missions to space. With the enhanced particle fluxes observed between 2008 and 2010, it is reasonable to assume that the radiation exposure from GCR must have also increased to unusually high levels. In this paper, the GCR exposure outside and inside the Earth’s magnetosphere is numerically calculated for time periods starting from 1970 to the end of 2011 in order to investigate the increase in dose levels during the years 2008–2010 in comparison with the last three solar minima. The dose rates were calculated in a water sphere, used as a surrogate for the human body, either unshielded or surrounded by aluminium shielding of 0.3, 10 or 40 g/cm².     

Study of atmospheric air pollutants during the partial solar eclipse on 15 January 2010 over Udaipur: A semi-arid location in Western India

The paper describes behavior of surface ozone, its precursor gases, BC along with TOCC, TWVC, AOT1020 nm as well as UV and IR radiation intensities observed during the partial solar eclipse of 15th January, 2010 over Udaipur, where 52% solar disc is obscured due to the moon’s shadow. During the beginning to main eclipse phase, the deviation values of several air pollutants concentrations from eclipse to control day values vary in a small range from −9 to −2 ppb in case of surface ozone and −180 to −80 ppb for CO. The corresponding change in the values of BC observed from −3.3 to −.5 μg/m³. No significant change is found in NO₂, NO or in ratio of NO₂/NO values during the partial eclipse time. TOCC values decrease from 3 to 5 DU along with a reduction in UV radiation intensity from 20 to 35% from starting to the main eclipse phase. The AOT1020 nm values are found to increase from .2 to 1.0 along with a reduction in IR radiation intensity order of 50%. However, TWVC values decrease from .22 to .1 cm during the eclipse hours. The low level of dilution in surface ozone in eclipse period may be attributed with change in local atmospheric boundary layer dynamic conditions or limited air pollutants dispersion, in term of decreases in planetary boundary layer height, wind speed and hence ventilation coefficient in the same eclipse hours. Thus, present studies support the argument for the leading roles of photochemical reactions with its precursor gases under presence of solar radiation in surface ozone variability. Other possible controlling factors are advection of air pollutants from the polluted region as evident from backward wind trajectories and altering the local meteorological conditions.