Solar flux effects on equatorial ionization anomaly and total electron content over Brazil: Observational results versus IRI representations

The F layer critical frequency (foF2) as measured by Digisondes in the equatorial and low latitude locations in Brazil is analyzed to investigate the seasonal and solar flux controls of the intensity of the equatorial ionization anomaly (EIA) in the equinoctial month of March. The analysis also included the total electron content (TEC) as measured by a GPS receiver operated at the EIA crest location. The foF2 data set covered a period of large solar flux variation from 1996 to 2003, while the GPS TEC data was for a period in 2002–2003 when the solar flux parameter F10.7 underwent large variations, permitting in both cases an examination of the solar flux effects on these parameters. The seasonal variation pattern in TEC shows a maximum in equinoctial months and a minimum in June solstice, with similar variations for foF2. The solar flux dependence of the TEC is a maximum during equinoxes, especially for post-sunset TEC values at times when the latitudinal distribution is controlled by the equatorial evening plasma fountain processes. Significant variations with local time are found in the degree of solar flux dependence for both the TEC and EIA. The EIA intensity shows large dependence on F10.7 during post-sunset to midnight hours. These results are discussed in comparison with their corresponding IRI representations.     

The post sunset equatorial F- region zonal drift variability and its linkage with equatorial spread F onset and duration over Indian longitudes

The seasonal and solar activity variation of the post sunset F- region zonal plasma drift, at the magnetic equatorial region over Indian longitudes is analyzed using the Republic of China Satellite-1 data from January 2000 to April 2004. The post sunset F- region zonal drifts are observed to be higher in the years of high solar activity in comparison with low solar activity, while seasonally the drifts are minimum in summer with much higher values in other seasons. The seasonal and solar activity variations of zonal plasma drift are attributed to the corresponding variations in the neutral winds. The dependences of the F region peak vertical drift on the zonal plasma drift at 18.5 IST (Indian Standard Time) and the time difference of the conjugate points sunset times, are quantitatively analyzed. Further an integrated parameter (incorporating the above mentioned two independent factors), which is able to predict the peak vertical drift and growth rate of Rayleigh Taylor instability is proposed. The other major outcome of the study is the successful prediction of the Equatorial Spread F (ESF) onset time and duration using the new integrated parameter at 18.5 IST. ESF irregularities and associated scintillations adversely affect communication and navigation systems. Hence, the present methodology for the prediction of the characteristics of these nocturnal irregularities becomes relevant.     

Longitudinal differences in the equatorial spread F characteristics between Vietnam and Brazil

We present the results of a comparative study of the equatorial spread F (ESF) and the F layer critical parameter, the base height of the F layer bottomside (h′F) over the two equatorial sites, Ho Chi Minh City – HCM (dip latitude: 2.9°N) in Vietnam and Sao Luis – SL (dip latitude: ∼2°S) in Brazil. The study utilizes simultaneous data collected by a CADI at HCM and a digisonde at SL during the year 2002 with the monthly mean solar 10.7 cm flux (F10.7) varying from ∼120 to ∼185. This study focuses on the quiet time seasonal behavior of the F layer parameters in the two widely separated longitude sectors, and addresses the question as to what can we learn from such comparative studies with respect to the ambient ionospheric and thermospheric parameters that are believed to control the ESF generation and hence its longitudinal occurrence pattern. The observed differences/similarities in the diurnal and seasonal patterns of the F Layer height vis-à-vis the ESF occurrences are evaluated in terms of the known longitudinal differences in the F layer heights, thermospheric meridional winds and the geomagnetic peculiarities of the two sites.     

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).     

Total electron content at low latitudes and its comparison with the IRI90

We used total electron content (TEC) data measured by Faraday rotation technique over Cachoeira Paulista (22.5°S, 45°W), in Brazil, to study the TEC variations with the solar flux at 10.7 cm (F10.7) and to compare the results with the IRI90 predictions. The data were divided into summer, equinox and winter. During the analysed period F10.7 varied from 66 up to 330. Our data shows that the observed TEC at 1600 LT (around the diurnal maximum) and at 0500 LT (around the diurnal minimum) increases with F10.7 until saturation is reached which occurs at F10.7≈210 to 220 for equinox and summer, and at F10.7≈180 for winter months. Comparison with the IRI90 predictions shows that IRI overestimates the TEC at 0500 LT for all solar flux values. At 1600 LT, IRI overestimates the observed TEC for low solar flux but underestimates it for high solar flux values.     

Thermospheric winds over Abuja during solar minimum period

Monthly variations of averaged nighttime thermospheric winds have been investigated over Abuja, Nigeria (Geographic: 9.06°N, 7.5°E; Geomagnetic: 1.60°S). The reports are based on Fabry-Perot interferometer measurements of Doppler shifts and Doppler broadening of the 630.0 nm spectral emission. The results were obtained during a period of weak solar activity with the solar flux (F10.7) typically below 70 solar flux units. Inspection of the average monthly thermospheric winds from October 2017 to December 2017 found December meridional winds to be more equatorward than the October and November winds. Zonal winds are eastward with pre-midnight maximum speeds going above 100 m/s. Compared to Jicamarca zonal winds in the Peruvian sector for the same month of October, the magnitude of maximum Abuja zonal wind speed is weaker. We compare the observed diurnal variation with the recently updated Horizontal wind model (HWM 14). Most of the observational features of thermospheric wind diurnal variation are captured in the model variation. The HWM14 generally showed good agreement with the Abuja October and November zonal wind observations but overestimates the December meridional winds. Expected longer period analysis of the results from Abuja will stimulate a better understanding of wind climatology over the West African sector.     

Equatorial F-region vertical ion drifts during quiet solar maximum

Median values of ionosonde h′F data acquired at Ibadan (Geographic:7.4°N, 3.9°E, Magnetic: dip 6°S, and magnetic declination, 3°W), Nigeria, West Africa, have been used to determine vertical ion drift (electric field) characteristics in the postsunset ionosphere in the African region during a time of high solar activity (average F_10.7 −208). The database spans from January and December 1958 during the era of International Geophysical Year (IGY) for geomagnetic quiet conditions. Bimonthly averaged diurnal variations patterns are very similar, but differ significantly in magnitude and in the evening reversal times. Also, monthly variations of F-region vertical ion drift reversal times inferred from the time of h′F maximum indicates early reversal during equinoxes and December solstice months except for the month of April. Late reversal is observed during the June solstice months. The equatorial evening prereversal enhancement in vertical ion drift (V_zp) occurs largely near 1900 LT with typical values 20–45 m/s. Comparison of Ibadan ionosonde V_zp with the values of prereversal peak velocity reported for Jicamarca (South America), Kodaikanal (India), and Scherliess and Fejer global model show considerable disparity. The changes of postsunset peak in virtual height of F-layer (h′F_P) with prereversal velocity peak V_zp are anti-correlated. Investigation of solar effects on monthly values of V_zp and h′F_P revealed that these parameters are independent of monthly averaged solar flux intensity during quiet-time sunspot maximum conditions.     

Equatorial anomaly according to the Interkosmos-19 data and IRI model: A comparison

The comparison of the IRI model with the foF2 distribution in the equatorial anomaly region obtained by topside sounding onboard the Interkosmos-19 satellite has been carried out. The global distribution of foF2 in terms of LT-maps was constructed by averaging Intercosmos-19 data for summer, winter, and equinox. These maps, in fact, represent an empirical model of the equatorial anomaly for high solar activity F10.7 ~ 200. The comparison is carried out for the latitudinal foF2 profiles in the characteristic longitudinal sectors of 30, 90, 210, 270, and 330°, as well as for the longitudinal variations in foF2 over the equator. The largest difference between the models (up to 60%) for any season was found in the Pacific longitudinal sector of 210°, where there are a few ground-based sounding stations. Considerable discrepancies, however, are sometimes observed in the longitudinal sectors, where there are many ground-based stations, for example, in the European or Indian sector. The discrepancies reach their maximum at 00 LT, since a decay of the equatorial anomaly begins before midnight in the IRI model and after midnight according to the Interkosmos-19 data. The discrepancies are also large in the morning at 06 LT, since in the IRI model, the foF2 growth begins long before sunrise. In the longitudinal variations in foF2 over the equator at noon, according to the satellite data, four harmonics are distinguished in the June solstice and at the equinox, and three harmonics in the December solstice, while in the IRI model only two and one harmonics respectively are revealed. In diurnal variations in foF2 and, accordingly, in the equatorial anomaly intensity, the IRI model does not adequately reproduce even the main, evening extremum.     

Variation in the ionospheric propagating factor M(3000)F2 at Ouagadougou,Burkina Faso

We use hourly monthly median values of propagation factor M(3000)F2 data observed at Ouagadougou Ionospheric Observatory (geographic12.4°N, 1.5°W; 5.9^o dip), Burkina Faso (West Africa) during the years Januar1987–December1988 (average F10.7 < 130 × 10⁻²² W/m²/Hz, representative of low solar flux conditions) and for January 1989–December1990 (average F10.7 ? 130 × 10⁻²² W/m²/Hz, representative of high solar epoch) for magnetically quiet conditions to describe local time, seasonal and solar cycle variations of equatorial ionospheric propagation factor M(3000)F2 in the African region. We show that that seasonal trend between solar maximum and solar minimum curves display simple patterns for all seasons and exhibits reasonable disparity with root mean square error (RMSE) of about 0.31, 0.29 and 0.26 for December solstice, June solstice and equinox, respectively. Variability Σ defined by the percentage ratio of the absolute standard deviation to the mean indicates significant dissimilarity for the two solar flux levels. Solar maximum day (10–14 LT) and night (22–02 LT) values show considerable variations than the solar minimum day and night values. We compare our observations with those of the IRI 2007 to validate the prediction capacity of the empirical model. We find that the IRI model tends to underestimate and overestimate the observed values of M(3000)F2, in particular, during June solstice season. There are large discrepancies, mainly during high solar flux equinox and December solstice between dawn and local midnight. On the other hand, IRI provides a slightly better predictions for M(3000)F2 between 0900 and 1500 LT during equinox low and high solar activity and equinox high sunspot number. Our data are of great importance in the area of short-wave telecommunication and ionospheric modeling.     

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.     

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.     

Low latitude nighttime ionospheric vertical E × B drifts at African region

The nighttime vertical E × B drifts velocities of the F2-region were inferred from the hourly hmF2 values obtained from ionosonde data over an African equatorial station, Ilorin (8.50^oN, 4.68^oE; dip lat. 2.95^o) during period of low solar activity. For each season, the plasma drift Vz is characterized by an evening upward enhancement, then by a downward reversal at 1900 LT till around 0000 LT, except for June solstice. This was explained using the Rayleigh–Taylor (R-T) instability mechanism. The occasional drift differences in Vz obtained by inferred and direct measurement over Ilorin and Jicamarca, respectively are reflective of the importance of chemistry and divergent transport system due to both the E region electric and magnetic fields instead of simple motions. The pre-reversal enhancement (PRE) magnitude is higher during the equinoctial months than the solsticial months over Jicamarca, highest during December solstice and the equinoctial months over Ilorin, suggesting the dominance of higher E × B fountain during equinoxes at both stations. The lowest PRE magnitude was in June solstice. The appearance of post-noon peak in NmF2 around 1700 LT is highest during the equinoctial months and lowest during the solsticial period. A general sharp drop in NmF2 around 1800 LT is distinct immediately after sunset, lowest during June solstice and highest in March equinox. Our result suggests that between 0930 and 2100 LT, the general theory that vertical drifts obtained by digisonde measurements only match the E × B drift if the F layer is higher than 300 km is reliable, but does not hold for the nighttime period of 2200–0600 LT under condition of solar minima. Hence, the condition may not be sufficient for the representation of vertical plasma drift at nighttime during solar minima. This assertion may still be tentative, as more equatorial stations needed to be studied for better confirmation.     

Temporal and spatial variation of equatorial ionization anomaly by using multistation ionosonde data for the 19th solar cycle over the Indian region

The paper deals with the study of temporal and spatial variation of equatorial ionization anomaly (EIA) phenomenon along with its dependence on solar activity and season during the 19th solar cycle by using seven Indian ionosonde stations. Present study is an attempt to carry out the comprehensive study of EIA by using the limited number of ground based instruments. This has been achieved by performing the Gaussian fitting over the latitudinal distribution of F2-region critical frequency (foF2) data. Results reveal that the phenomenon of EIA has a strong dependence on solar activity and seasons. The EIA crest exhibits the feature of latitudinal shifting and expansion with increasing solar activity. It is found out that the effect of solar cycle and seasons on EIA is local time dependent. The observations were also compared with the IRI-2001 model predictions and results reveal that the model values are in general agreement with the observed values with some discrepancies, particularly during the high solar activity period and morning sector. The results have been discussed in the light of relative contribution from transequatorial interhemispheric neutral wind and strength of equatorial fountain process during different local time, season and solar activity levels. Furthermore, an attempt is made to parameterize the location and foF2 of the EIA crest by using the regression analysis. These results can be used to predict the latitudinal position and foF2 of the EIA crest for any given 12-month running average sunspot number (R12).     

Characteristics of the equatorial F-region vertical plasma drift observed during post-sunset and pre-sunrise hours

Post-sunset and pre-sunrise vertical plasma drifts at the equatorial F-region have been investigated using the HF Doppler radar and ionosonde observations. Observed vertical plasma drift features during the sunrise are found to complement that observed during the evening. The post-sunset vertical plasma drift is characterized by an upward enhancement, a pre-reversal enhancement and a reversal in the drift direction. Similarly, the pre-sunrise plasma drift is characterized by a sudden downward excursion followed by an upward turning. The wavelet analysis of the plasma drift shows the presence of fluctuations in the period range 4–32 min and the short period fluctuations are attributed to the atmospheric gravity waves.     

Parameterized Regional Ionospheric Model and a comparison of its results with experimental data and IRI representations

We describe a Parameterized Regional Ionospheric Model (PARIM) to calculate the spatial and temporal variations of the ionospheric electron density/plasma frequency over the Brazilian sector. The ionospheric plasma frequency values as calculated from an enhanced Sheffield University Plasmasphere–Ionosphere Model (SUPIM) were used to construct the model. PARIM is a time-independent 3D regional model (altitude, longitude/local time, latitude) used to reproduce SUPIM plasma frequencies for geomagnetic quiet condition, for any day of the year and for low to moderately high solar activity. The procedure to obtain the modeled representation uses finite Fourier series so that all plasma frequency dependencies can be represented by Fourier coefficients. PARIM presents very good results, except for the F region peak height (hmF2) near the geomagnetic equator during times of occurrence of the F₃ layer. The plasma frequency calculated by IRI from E region to bottomside of the F region present latitudinal discontinuities during morning and evening times for both solar minimum and solar maximum conditions. Both the results of PARIM and the IRI for the E region peak density show excellent agreement with the observational values obtained during the conjugate point equatorial experiment (COPEX) campaign. The IRI representations significantly underestimate the foF2 and hmF2 compared to the observational results over the COPEX sites, mainly during the evening–nighttime period.     

Solar activity variations of equatorial spread F occurrence and sustenance during different seasons over Indian longitudes: Empirical model and causative mechanisms

A comprehensive analysis using nearly two decades of ionosonde data is carried out on the seasonal and solar cycle variations of Equatorial Spread F (ESF) irregularities over magnetic equatorial location Trivandrum (8.5°N, 77°E). The corresponding Rayleigh Taylor (RT) instability growth rates (γ) are also estimated. A seasonal pattern of ESF occurrence and the corresponding γ is established for low solar (LSA), medium solar (MSA) and high solar (HSA) activity periods. For LSA, it is seen that the γ maximizes during post sunset time with comparable magnitudes for autumnal equinox (AE), vernal equinox (VE) and winter solstice (WS), while for summer solstice (SS) it maximizes in the post-midnight period. Concurrent responses are seen in the ESF occurrence pattern. For MSA, γ maximizes during post-sunset for VE followed by WS and AE while SS maximises during post-midnight period. The ESF occurrence for MSA is highest for VE (80%), followed by AE (70%), WS (60%) and SS (50%). In case of HSA, maximum γ occurs for VE followed by AE, WS and SS. The concurrent ESF occurrence maximizes for VE and AE (90%), WS and SS at 70%.The solar cycle variation of γ is examined. γ shows a linear variation with F10.7?cm flux. Further, ESF percentage occurrence and duration show an exponential and linear variation respectively with γ. An empirical model on the solar activity dependence of ESF occurrence and sustenance time over Indian longitudes is arrived at using the database spanning two solar cycles for the first time.     

Performance evaluation of IRI-2007 at equatorial latitudes and its Matlab version for GNSS applications

International Reference Ionosphere (IRI) model is the widely used empirical model for ionospheric predictions, especially TEC which is an important parameter for radio navigation and communication. The Fortran based IRI-2007 does not support real-time interactive visualization and debugging. Therefore, the source code is converted into Matlab and is validated for the purposes of this study. This facilitates easy representation of results and for near real-time implementation of IRI in the applications including spacecraft launching, now casting, pseudolite based navigation systems etc. In addition, the vertical delay results over the equatorial region derived from IRI and GPS data of three IGS stations namely Libreville (Garbon, Africa), Brasilia (Brazil, South America) and Hyderabad (India, Asia) are compared. As the IRI model does not account for plasmasphere TEC, the vertical delays are underestimated compared to vertical delays of GPS signals. Therefore, the model should be modified accordingly for precise TEC estimation.     

Characterization of GPS-TEC over African equatorial ionization anomaly (EIA) region during 2009–2016

This study characterizes total electron content (TEC) measured by Global Positioning System (GPS) over African equatorial ionization anomaly (EIA) region for 2009–2016 period during both quiet geomagnetic conditions (Kp?≤?1) and normal conditions (1?>?Kp?≤?4). GPS-TEC data from four equatorial/low-latitude stations, namely, Addis Ababa (ADIS: 9.04°N, 38.77°E, mag. lat: 0.2°N) [Ethiopia]; Yamoussoukro (YKRO: 6.87°N, 5.24°W, mag. lat: 2.6°S) [Ivory Coast]; Malindi (MAL2; 3.00°S, 40.19°E, mag. lat: 12.4°S) [Kenya] and Libreville (NKLG; 0.35°N, 9.67°W, mag. lat: 13.5°S) [Gabon] were used for this study. Interesting features like noontime TEC bite-out, winter anomaly during the ascending and maximum phases of solar cycle 24, diurnal and seasonal variations with solar activity have been observed and investigated in this study. The day-to-day variations exhibited ionospheric TEC asymmetry on an annual scale. TEC observed at equatorial stations (EIA-trough) and EIA-crest reach maximum values between ～1300–1600 LT and ～1300–1600 LT, respectively. About 76% of the high TEC values were recorded in equinoctial months while the June solstice predominantly exhibited low TEC values. Yearly, the estimated TEC values increases or decreases with solar activity, with 2014 having the highest TEC value. Solar activity dependence of TEC within the EIA zone reveals that both F10.7?cm index and EUV flux (24–36?nm) gives a stronger correlation with TEC than Sunspot Number (SSN). A slightly higher degree of dependence is on EUV flux with the mean highest correlation coefficient (R) value of 0.70, 0.83, 0.82 and 0.88 for quiet geomagnetic conditions (Kp?≤?1) at stations ADIS, MAL2, NKLG, and YKRO, respectively. The correlation results for the entire period consequently reveals that SSN and solar flux F10.7?cm index might not be an ideal index as a proxy for EUV flux as well as to measure the variability of TEC strength within the EIA zone. The estimated TEC along the EIA crest (MAL2 and NKLG) exhibited double-hump maximum, as well as post-sunset peaks (night time enhancement of TEC) between ～2100 and 2300 LT. EIA formation was prominent during evening/post-noon hours.     

Equatorial total electron content (TEC) at low and high solar activity

Total electron content (TEC) derived from ionosonde data recorded at the station of Korhogo (Lat = 9.33°N, Long = 5.43°W, Dip = 0.67°S) are compared to the International Reference Ionosphere (IRI) model predicted TEC for high (1999) and low (1994) solar activity conditions. The results show that the model represents the diurnal variation of the TEC as well as a solar activity and seasonal dependence. This variation is closer to that of the ionosonde-inferred TEC at high solar activity. However, at low solar activity the IRI overestimates the ionosonde-inferred TEC. The relative deviation ΔTEC is more prominent in the equinoctial seasons during nighttime hours where it is as high as 70%. At daytime hours, the relative deviation is estimated to 0–30%.