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.     

Seasonal patterns of condensation and sublimation cycles in the cryptic and non-cryptic regions of the South Pole

The South Pole of Mars is characterized by an asymmetric residual ice cap composed of water ice and CO₂ ice. On the opposite side of the residual cap, there exists an area called cryptic region which is relatively free of ice during summer time. Many fan-shaped km-scale structures apparently caused by a wind-blown system of dust-laden gas jets occurred dozens degrees of Ls before the complete sublimation of the CO₂ frost layer. We have examined the seasonal cycles of condensation and sublimation in the cryptic and non-cryptic regions by using the topographic data from the MOLA/MGS measurements. Using the MOLA topography data collected over one Martian year (1999–2001), we have studied the temporal elevation change and the seasonal cycle of the carbon dioxide frost on the southern polar caps. We have produced mapping of the seasonal CO₂ frost thickness variation for seven Ls (30°, 60°, 90°, 120°, 150°, 180°, 210°, 240°, 270° and 330°). It is found that the time variations of the CO₂ frost thickness in these two regions are quite similar. The greatest thickness of the CO₂ frost layer is about 0.76–0.78 m in both places occurs at Ls = 150°.     

A feasibility study for the detection of the diurnal variation of tropospheric NO2 over Tokyo from a geostationary orbit

We have conducted a feasibility study for the geostationary monitoring of the diurnal variation of tropospheric NO₂ over Tokyo. Using NO₂ fields from a chemical transport model, synthetic spectra were created by a radiative transfer model, SCIATRAN, for summer and winter cases. We then performed a Differential Optical Absorption Spectroscopy (DOAS) analysis to retrieve NO₂ slant column densities (SCDs), and after converting SCDs into vertical column densities (VCDs), we estimated the precision of the retrieved VCDs. The simulation showed that signal-to-noise ratio (SNR) ? 500 is needed to detect the diurnal variation and that SNR ? 1000 is needed to observe the local minimum occurring in the early afternoon (LT13–14) in summer. In winter, the detection of the diurnal variation during LT08–15 needs SNR ? 500, and SNR ? 1000 is needed if early morning (LT07) and early evening (LT16) are included. The currently discussed sensor specification for the Japanese geostationary satellite project, GMAP-Asia, which has a horizontal resolution of 10 km and a temporal resolution of 1hr, has demonstrated the performance of a precision of several percent, which is approximately corresponding to SNR = 1000–2000 during daytime and SNR ? 500 in the morning and evening. We also discuss possible biases caused by the temperature dependence of the absorption cross section utilized in the DOAS retrieval, and the effect of uncertainties of surface albedo and clouds on the estimation of precisions.     

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.     

Variation of monthly mean foF2 and hmF2 over a mid latitude station during the period 1997–2006

Ionosonde data of a mid latitude station Novosibirsk (Geog. Lat. 54.6°N, Geog. Long. 83.2°E) has been analyzed for the years 1997–2006 that covers the major part of solar cycle 23. Our results show the presence of winter anomaly in the daytime F2 layer critical frequency during different phases of solar activity. Results also reveal a semiannual variation of foF₂ with two maxima and a minimum that always appears in summer. While the first maximum is in the spring equinox, the second one is found to shift from autumn to winter with the increase of solar activity. The maximum height of F2 layer during the daytime shows variation with the solar activity. It is higher during the higher activity periods and lower during the periods of low activity. Results of ionosonde observations have been compared with those obtained from IRI-2007 model and it is found that model reproduces gross features of foF₂ variation. However, the modeled hmF₂ variations during equinoxes are significantly different from the ones derived using the ionosonde data. The model also underestimates the hmF₂ values.     

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.     

Study of maximum electron density NmF2 at Karachi and Islamabad during solar minimum (1996) and solar maximum (2000) and its comparison with IRI

The monthly hourly medians of maximum electron density, N_mF₂, at two Pakistani ionospheric stations, Karachi and Islamabad, have been determined for solar minimum (1996) and solar maximum (2000) and compared with IRI predictions using the URSI coefficients. At night and pre-noon period the N_mF₂ values at both stations are almost equal during the 2 years. However, at post-noon the values at Karachi are considerably larger than those at Islamabad due to the equatorial or geomagnetic anomaly. Karachi (geomag. coord. 16.44°N, 139.08°E) lies near the region of the equatorial anomaly (+20 and −20 geomagnetic latitude), so most of the N_mF₂ values at Karachi are larger than those at Islamabad (geomag. coord. 24.46°N, 145.67°E). The maximum monthly values of N_mF₂ show a semi-annual variation at Karachi and Islamabad both during 1996 and 2000 as predicted by IRI.     

Influence of solar-geomagnetic disturbances on SABER measurements of 4.3 μm emission and the retrieval of kinetic temperature and carbon dioxide

Thermospheric infrared radiance at 4.3 μm is susceptible to the influence of solar-geomagnetic disturbances. Ionization processes followed by ion-neutral chemical reactions lead to vibrationally excited NO⁺ (i.e., NO⁺(v)) and subsequent 4.3 μm emission in the ionospheric E-region. Large enhancements of nighttime 4.3 μm emission were observed by the TIMED/SABER instrument during the April 2002 and October–November 2003 solar storms. Global measurements of infrared 4.3 μm emission provide an excellent proxy to observe the nighttime E-region response to auroral dosing and to conduct a detailed study of E-region ion-neutral chemistry and energy transfer mechanisms. Furthermore, we find that photoionization processes followed by ion-neutral reactions during quiescent, daytime conditions increase the NO⁺ concentration enough to introduce biases in the TIMED/SABER operational processing of kinetic temperature and CO₂ data, with the largest effect at summer solstice. In this paper, we discuss solar storm enhancements of 4.3 μm emission observed from SABER and assess the impact of NO⁺(v) 4.3 μm emission on quiescent, daytime retrievals of Tk/CO₂ from the SABER instrument.     

Study of satellite retrieved CO2 and CH4 concentration over India

This paper reports a study of spatial and temporal variations of columnar averaged concentration of CO₂ and CH₄ over India using SCanning Imaging Absorption spectroMeter for Atmospheric CHartographY (SCIAMACHY) and Greenhouse gas Observing SATellite (GOSAT) data. Comparison of these data with the global view National Oceanic and Atmospheric Administration (NOAA) land data and also location specific flask data is made. The temporal variation in column averaged global CO₂ is similar to that over India and it is also similar to the NOAA surface flask data and global view. The variation in NOAA surface CH₄ is location dependent and its global view appears to vary seasonally in opposite phase with the column averaged CH₄ values from satellites, reflecting the limited comparability of surface and column averaged data. Over India the CO₂ maximum is in May and minimum in August/September while for CH₄ the maximum is in September and minimum in February/March. The seasonal variation of CH₄ over India is correlated with the eastern coastal rice cultivation.     

A study of L band scintillations during the initial phase of rising solar activity at an Indian low latitude station

The ionospheric scintillation and TEC (Total Electron Content) variations are studied using GPS (Global Positioning System) measurements at an Indian low latitude station Surat (21.16°N, 72.78°E; Geomagnetic: 12.90°N, 147.35°E), situated near the northern crest of the equatorial anomaly region. The results are presented for data collected during the initial phase of current rising solar activity (low to moderate solar activity) period between January 2009 and December 2011. The results show that within a total number of 656 night-time scintillation events, 340 events are observed with TEC depletions, Rate of change of TEC (ROT) fluctuations and enhancement of Rate of change of TEC Index (ROTI). A comparison of night-time scintillation events from the considered period reveal strong correlation amongst the duration of scintillation activity in S₄ index_, TEC depletion, ROT fluctuations and ROTI enhancement in the year 2011, followed by the year 2010 and least in 2009. The statistical analyses of scintillation activity with enhancement of ROTI also show that about 70–96% scintillation activity took place in equinox and winter months. Moreover, from a nocturnal variation in occurrence of scintillation with (S₄ ? 0.2) and enhancement of ROTI with (ROTI ? 0.5), a general trend of higher occurrence in pre-midnight hours of equinox and winter seasons is observed in both indices during the year 2011 and 2010, while no significant trend is observed in the year 2009. The results suggest the presence of F-region ionospheric irregularities with scale sizes of few kilometers and few hundred meters over Surat and are found to be influenced by solar and magnetic activity.     

Pre-storm behaviour of NmF2 and TEC (GPS) over equatorial and low latitude stations in the Indian sector

The pre-storm behavior of NmF₂ and TEC over an equatorial station, Trivandrum (8.47°N, 76.91°E, dip 0.6°S) and a low latitude station, Waltair (17.7°N,83.3°E, dip 20°N) has been studied for a total of 18 strong geomagnetic storms with DST ? −100 nT. The simultaneous measurements of GPS-TEC and NmF₂ over Trivandrum and Waltair during the period 2000–2005 have been considered for the present study. It is found that there is a substantial increase in NmF₂ and TEC before the onset of the storm over Waltair, while the increase is not present at Trivandrum. The origin of pre-storm enhancements in electron density still remains unresolved owing to several conditions in their potential sources and occurrence mechanisms. In the present study an attempt is made to identify the possible mechanisms responsible for such enhancements in electron density of the F-region.     

Variability of foE in the equatorial ionosphere with solar activity

This research examined the variability of f_oE in the equatorial ionosphere with solar activity within the equatorial ionospheric anomaly region. Ionosonde data recorded at Ouagadougou (lat. 12.4°N, long. 1.5°W and magnetic dip 1.43°N) were engaged to study the transient variations of the critical frequency of the E-layer (f_oE) and its dependence on solar activity. The study revealed that f_oE increases with the increase in solar intensity of the sun. The variability of the f_oE decreases with increases in the solar activity. The maximum value of the f_oE is at local noon when the ionosphere is stable; the variability at this local time is minimal. The minimum value of the f_oE is at sunrise and sunset, at this period on local time the equatorial ionosphere recorded its maxima variability. Irrespective of the degree of solar activity, f_oE is observed to be maximum in June solstice, followed by the equinoxes and minimum in December solstice. Equinoctial asymmetry occurred in the variation of the relative standard deviation of f_oE with maximum in September/March equinox for low/high solar activity.     

Variations of total electron content over high latitude region during the ascending phase of 24th solar cycle

The solar cycle variation and seasonal changes significantly affects the ionization process of earth’s ionosphere and required to be monitored in real time basis for regional level refinement of existing models. In view of this, the present study has been carried out by using the ionospheric Total Electron Content (TEC) data observed with the help of Global Ionospheric Scintillation and TEC monitoring (GISTM) system installed at Indian Antarctic Research Station, “Maitri” [70°46′00″S 11°43′56″E] during the ascending phase of 24th solar cycle. The daily values of solar extreme ultraviolet (EUV) flux (0.1–50?nm wavelength), 10.7?cm radio flux F_10.7 and Sunspot number (SSN) has been taken as a proxy to represent the solar cycle variation to correlate with TEC. The linear regression results revels better correlation of TEC with EUV flux rather than F_10.7 and SSN. Also, the EUV and TEC show better agreement during summer as compared to winter and equinox period. Correlation between TEC and EUV appears significantly noticeable during ten internationally defined quiet days of each month (stable background geophysical condition) as compared to the overall days (2010–2014). Further, saturation effect has been observed on TEC values during the solar maxima year 2014. The saturation effects are more prominent during the night hours of winter and equinox season due to transportation losses manifested by the equator-ward direction of meridional wind.     

An evaluation of the IRI-2007 storm time model at low latitude stations

This paper discusses the ability of the International Reference Ionosphere IRI-2007 storm time model to predict foF2 ionospheric parameter during geomagnetic storm periods. Experimental data (based on availability) from two low latitude stations: Vanimo (geographic coordinates, 2.7 °S, 141.3 °E, magnetic coordinates, 12.3 °S, 212.50 °E) and Darwin (geographic coordinates, 12.45 °S, 130.95 °E, magnetic coordinates, 22.9 °S, 202.7 °E) during nine storms that occurred in 2000 (Rz₁₂ = 119), 2001(Rz₁₂ = 111) and 2003 (Rz₁₂ = 64) are compared with those obtained by the IRI-2007 storm model. The results obtained show that the percentage deviation between the experimental and IRI predicted foF2 values during these storm periods is as high as 100% during the main and recovery phases. Based on the values of “relative deviation module mean” (RDMM) obtained (i.e. between 0.08 and 0.60), it is observed that there is a reasonable to poor agreement between measured foF2 values and the IRI-storm model prediction values during main and recovery phases of the storms under investigation. As a result, in addition to other studies that have been carried out from different sectors, more studies are required to be carried out. This will enable IRI community to improve on the present performance of the model. In general the IRI-storm model predictions follow normal trend of the foF2 measured values but does not reproduce well the measured values.     

Equatorial Ionospheric Anomaly (EIA) and comparison with IRI model during descending phase of solar activity (2005–2009)

The ionospheric variability at equatorial and low latitude region is known to be extreme as compared to mid latitude region. In this study the ionospheric total electron content (TEC), is derived by analyzing dual frequency Global Positioning System (GPS) data recorded at two stations separated by 325 km near the Indian equatorial anomaly region, Varanasi (Geog latitude 25°, 16^/ N, longitude 82°, 59^/ E, Geomagnetic latitude 16°, 08^/ N) and Kanpur (Geog latitude 26°, 18^/ N, longitude 80°, 12^/ E, Geomagnetic latitude 17°, 18^/ N). Specifically, we studied monthly, seasonal and annual variations as well as solar and geomagnetic effects on the equatorial ionospheric anomaly (EIA) during the descending phase of solar activity from 2005 to 2009. It is found that the maximum TEC (EIA) near equatorial anomaly crest yield their maximum values during the equinox months and their minimum values during the summer. Using monthly averaged peak magnitude of TEC, a clear semi-annual variation is seen with two maxima occurring in both spring and autumn. Results also showed the presence of winter anomaly or seasonal anomaly in the EIA crest throughout the period 2005–2009 only except during the deep solar minimum year 2007–2008. The correlation analysis indicate that the variation of EIA crest is more affected by solar activity compared to geomagnetic activity with maximum dependence on the solar EUV flux, which is attributed to direct link of EUV flux on the formation of ionosphere and main agent of the ionization. The statistical mean occurrence of EIA crest in TEC during the year from 2005 to 2009 is found to around 12:54 LT hour and at 21.12° N geographic latitude. The crest of EIA shifts towards lower latitudes and the rate of shift of the crest latitude during this period is found to be 0.87° N/per year. The comparison between IRI models with observation during this period has been made and comparison is poor with increasing solar activity with maximum difference during the year 2005.     

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.     

Variability study of ionospheric total electron content at crest of equatorial anomaly in China from 1997 to 2007

The variation of TEC data at Wuhan station (geographic coordinate: 30.5°N, 114.4°E; geomagnetic coordinate: 19.2°N, 183.8°E) at crest of equatorial anomaly in China from January 1997 to December 2007 were analyzed. Variability with solar activity, annual, semiannual, diurnal and seasonal variation were also analyzed. The MSIS00 model and ISR model were used to analyze the possible mechanisms of the variabilities found in the results. The TEC data in 1997 and 2001 deduced from another crest station Xiamen (geographic coordinate: 24.4°N, 118.1°E; geomagnetic coordinate: 13.2°N, 187.4°E) were used to contrast. Analysis results show that long-term variations of TEC at Xiamen station are mainly controlled by the variations of solar activities. Typical diurnal variation behaves as a minimum of the TEC in the pre-dawn hours around 05:00–06:00LT and a maximum on the afternoon hours around 13:00–15:00LT. Some features like the semiannual anomaly and winter anomaly in TEC have been reported. The anomaly may be the result of common action of the electric field over the magnetic equatorial and the [O/N₂] at the crest station.     

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.     

NO2 field observations in the stratosphere during Map/Globus campaign in 1983 and 1985

NO₂ has been measured from balloons, plane and ground during Map/Globus Campaign in Europe in 1983. The results have shown an agreement between in situ and remote sensing techniques in the lower stratosphere, but some divergence at higher altitude. Part of the differences might be due to unknown experimental systematic errors and part to real atmospheric variations as the observations were not performed on the same day. In order to improve the knowledge of the instruments performances and to evaluate the NO_x stratospheric budget a new field campaign has been set up in Europe in 1985. Simultaneous in situ and balloons, ground and satellites remote sensing observations of NO and NO₂ were obtained, covering a 24 H period.     

Prediction of clear-sky biologically effective erythematic radiation (EER) From global solar radiation (250–2800 nm) at Cairo,Egypt

The hourly and daily measured clear-sky global solar radiation (G) and biologically important effective erythematic radiation (EER) incident on a horizontal surface at Cairo, Egypt (latitude 30° 05′ N & Longitude 31° 15′ E), during the period from January 1995 to December 2005 are used in this paper. The relationship between daily integrated totals of EER and the daily totals of broadband global solar radiation (250–2800 nm) is established. The temporal variability of the percentage ratio of the total daily erythema to total daily broadband solar global irradiation (EER/G) is determined. The monthly and the seasonal averages of the extraterrestrial UVB solar radiation, Mesurad and estímated UVB solar radiation and clearness index K_tUVB of UVB radiation are discussed. The average monthly mean variation of slant ozone (Z) and UVB transmission (K_tUVB) at the present work are found. The two variables show an opposite seasonal behavior, and the average monthly of slant ozone column and UVB transmission values shows the relationship between them in a clearer way than those of daily values. The estimated values of UVB solar radiation a good agreement with the measured values of the UVB solar radiation, the difference between the estimated and measured values of UVB solar radiation varies from 1.2% to 2.8%. The effect of the annual cycles of solar zenith angle (SZA) and total column ozone (TCO) on the ratios (EER/G) are presented and the correction factors are determined for removal of the ozone cycle. The seasonal variability of EER/G is also discussed. The effect of the annual cycles of solar zenith angle (SZA) and total column ozone (TCO) on the ratios (EER/G) is presented and the correction factors are determined for removal of the ozone cycle.