The impact of microwave absorber and radome geometries on GNSS measurements of station coordinates and atmospheric water vapour

We have used microwave absorbing material in different geometries around ground-based Global Navigation Satellite System (GNSS) antennas in order to mitigate multipath effects on the estimates of station coordinates and atmospheric water vapour. The influence of a hemispheric radome – of the same type as in the Swedish GPS network SWEPOS – was also investigated. Two GNSS stations at the Onsala Space Observatory were used forming a 12 m baseline. GPS data from October 2008 to November 2009 were analyzed by the GIPSY/OASIS II software using the Precise Point Positioning (PPP) processing strategy for five different elevation cutoff angles from 5° to 25°. We found that the use of the absorbing material decreases the offset in the estimated vertical component of the baseline from ∼27 mm to ∼4 mm when the elevation cutoff angle varies from 5° to 20°. The horizontal components are much less affected. The corresponding offset in the estimates of the atmospheric Integrated Water Vapour (IWV) decreases from ∼1.6 kg/m² to ∼0.3 kg/m². Changes less than 5 mm in the offsets in the vertical component of the baseline are seen for all five elevation cutoff angle solutions when the antenna was covered by a hemispheric radome. Using the radome affects the IWV estimates less than 0.4 kg/m² for all different solutions. IWV comparisons between a Water Vapour Radiometer (WVR) and the GPS data give consistent results.     

Optical studies of the ionospheric irregularities over the brazilian region by nocturnal images of the OI 630 nm emission

This study is an extension of previous statistical studies (Sobral et al., 1990, 1991, 1999) of both the local time and latitude variations of the zonal drift velocities of ionospheric plasma depletions, over the Brazilian low latitude station Cachoeira Paulista — CP (22.54°S, 45.00°W). The past studies were based on OI 630 nm scanning photometer data and the present one is based on digital OI 630 nm airglow images obtained by an all-sky imager system. These data were gathered between October 1998 and October 1999, at CP. The present results show that, in general, the velocities clearly tended to decrease with local time. Such a decrease should be associated with decreasing intensity of the vertical component of the ambient electric field which, in turn can be accounted for by recombination. All zonal drifts obtained for the 18 nights were eastwards. During equinox, the velocities clearly tended to decrease with local time at lower rates as compared with spring and summer. The highest and lowest zonal drift velocities, from all three seasons considered here, were observed to be in the summer ≈180 ms⁻¹ at 21:45 LT, and in the spring ≈25 ms⁻¹ at 03:15 LT, respectively. Ionospheric plasma bubbles were detected out to the maximum extra-tropical geographical latitude of ≈28° S, which was the highest latitude position analyzed in this study.     

Coordinated in situ measurements of plasma irregularities and ground based scintillation observations at the crest of equatorial anomaly

First comparison of in situ density fluctuations measured by the DEMETER satellite with ground based GPS receiver measurements at the equatorial anomaly station Bhopal (geographic coordinates (23.2°N, 77.6°E); geomagnetic coordinates (14.29°N, 151.12°E)) for the low solar activity year 2005, are presented in this paper. Calculation of the diurnal maximum of the strength of the equatorial electrojet, which can serve as precursor to ionospheric scintillations in the anomaly region is also done. The Langmuir Probe experiment and Plasma Analyzer onboard DEMETER measure the electron and ion densities respectively. Irregularities in electron density distribution cause scintillations on transionospheric links and there exists a close relationship between an irregularity and scintillation. In 40% of the cases, DEMETER detects the irregularity structures (dNe/Ne ? 5% and dNi/Ni (O⁺) ? 5%) and GPS L band scintillations (S4 ? 0.2) are also observed around the same time, for the low solar activity period. It is found that maximum irregularity intensity is obtained in the geomagnetic latitude range of 10–20° for both electron density and ion density. As the GPS signals pass through this irregularity structure, scintillations are recorded by the GPS receiver installed at the equatorial anomaly station, Bhopal it is interesting to note that in situ density fluctuations observed on magnetic flux tubes that pass over Bhopal can be used as indicator of ionospheric scintillations at that site. Many cases of density fluctuations and associated scintillations have been observed during the descending low solar activity period. The percentage occurrence of density irregularities and scintillations shows good correspondence with diurnal maximum of the strength of electrojet, however this varies with different seasons with maximum correspondence in summer (up to 66%) followed by equinox (up to 50%) and winter (up to 46%). Also, there is a threshold value of EEJ strength to produce density irregularities ((dNe/Ne)_max ? 5%) and for moderate to strong scintillations (S4 ? 0.3) to occur. For winter this value is found to be ∼40 nT whereas for equinox and summer it is around 50 nT.     

Impact of geomagnetic variation over sub-auroral ionospheric region during high solar activity year 2014

The present work is an attempt to evaluate the impact of changing space weather condition over sub-auroral ionosphere during high solar activity year 2014. In view of this, the GPS based TEC along with Ionosonde data over Indian permanent scientific base “Maitri”, Antarctica (70°46′00″S, 11°43′56″E) has been utilized. The results suggested that the nature of ionospheric responses to the geomagnetic disturbances not only depended upon the status of high latitudinal electro-dynamic processes but also influenced by the seasonal variations. The results revel both negative and positive type of ionospheric response in a single year but during different seasons. The study suggested that the combination of equator-ward plasma transportation along with ionospheric compositional changes causes a negative ionospheric impact especially during summer and equinox seasons. However, the combination of pole-ward contraction of the oval region along with particle precipitation may lead to exhibit positive ionospheric response during the winter season. The plasma transportation direction has been validated with the help of convection boundary (HM boundary) deduced with the help of SuperDARN observations. The ground based ionosonde observations clearly provided the evidence of deep penetration of high energetic particles up to the E-layer heights which results a sudden and strong appearance of E-layer. The strengthening of E-layer is responsible for modification of auroral electrojet and field-aligned current system. Also, the sudden appearance of E-layer along with a decrease in F-layer electron density suggested the dominance of NO⁺ over O⁺ in a considered region under geomagnetic disturbed condition.     

Plasmaspheric electron content derived from GPS TEC and FORMOSAT-3/COSMIC measurements: Solar minimum condition

The plasmaspheric electron content (PEC) was estimated by comparison of GPS TEC observations and FORMOSAT-3/COSMIC radio occultation measurements at the extended solar minimum of cycle 23/24. Results are retrieved for different seasons (equinoxes and solstices) of the year 2009. COSMIC-derived electron density profiles were integrated up to the height of 700 km in order to retrieve estimates of ionospheric electron content (IEC). Global maps of monthly median values of COSMIC IEC were constructed by use of spherical harmonics expansion. The comparison between two independent measurements was performed by analysis of the global distribution of electron content estimates, as well as by selection specific points corresponded to mid-latitudes of Northern America, Europe, Asia and the Southern Hemisphere. The analysis found that both kinds of observations show rather similar diurnal behavior during all seasons, certainly with GPS TEC estimates larger than corresponded COSMIC IEC values. It was shown that during daytime both GPS TEC and COSMIC IEC values were generally lower at winter than in summer solstice practically over all specific points. The estimates of PEC (h > 700 km) were obtained as a difference between GPS TEC and COSMIC IEC values. Results of comparative study revealed that for mid-latitudinal points PEC estimates varied weakly with the time of a day and reached the value of several TECU for the condition of solar minimum. Percentage contribution of PEC to GPS TEC indicated the clear dependence from the time with maximal values (more than 50–60%) during night-time and lesser values (25–45%) during day-time.     

Preliminary results of fair-weather atmospheric electric field in the proximity of Main Boundary Thrust,Northern Pakistan

The atmospheric electric Potential Gradient (PG) variation of fair-weather days at Muzaffarabad (MZF, Northern Pakistan) station is presented for the period of January 2015–October 2017. The present investigations focus on the diurnal variation of atmospheric PG on the seasonal and annual time scale. The fair-weather seasonal and annual average PG variation has revealed two peaks, i.e. a primary maximum peak and a secondary maximum peak. The average maximum PG value found is ～410?Vm^?1. The results are then compared with the well-known standard oceanic Carnegie curve. The diurnal curve of PG is found to deviate from the Carnegie curve. The seasonal PG variation at MZF shows lower values during the summer (monsoon) and autumn (post-monsoon) as compared to the winter and spring (pre-monsoon) which could be due to local aerosol concentration. Global comparison of PG (%) values of the annual mean at different longitude, with MZF observatory, shows local morning and evening peaks depicting local sunrise and sunset effects on the PG. The overall outcomes will certainly contribute to further investigate the Global Electric Circuit (GEC).     

第23太阳活动周行星际磁场扇形结构对中低纬地磁场的影响

用特征向量分析法对第23太阳活动周天津静海磁场强度水平分量H的时均值进行研究,分析行星际磁场扇形结构的地磁效应（简称扇形效应）对中低纬地磁场H分量日变化的贡献.研究结果表明,中低纬扇形效应为3～11nT,在太阳活动高年扇形效应达到最大值（约11nT）,低年达到最小值（约3nT）.太阳活动高年扇形效应引起的地磁H分量值变化与太阳活动低年的情况相反,但是扇形效应在夏季对地磁H分量的影响较小.太阳活动高年扇形效应日均值的增减与上升年的相反,与下降年相同,夏季扇形效应平均增量最小且无规律.春、夏和秋三个季节的扇形效应最大值都出现在太阳活动高年,冬季的扇形效应在太阳活动峰年两年后才出现最大值（约11nT）.在太阳活动低年（或高年）,当扇形磁场背离（或指向）太阳时,夏季扇形效应白天引起地磁H分量增大（或减小）,夜间导致地磁H分量减小（或增大）,其他季节全天都会导致地磁H分量减小（或增大）.用特征向量推断行星际磁场扇形极性的符合率在春夏秋三个季节高达60%,在冬季为55%.      

太阳活动低年低纬地区VTEC 变化特性分析

利用福州台站(26.1°N, 119.3°E, 磁纬14.4°N)电离层闪烁与TEC监测仪2006-2010年的观测数据, 对该地区垂直总电子含量(VTEC)进行时间变化特性分析. 结果表明, 春秋冬三季的VTEC平均最高值出现在06:00UT, 夏季出现在08:00UT, 所有季节的平均最低值均出现在21:00UT; VTEC变化存在季节异常和弱冬季异常, 春秋季节高, 冬夏季节低, 夏季VTEC比冬季低且最大值出现时间延迟; VTEC在2006-2009年呈现下降的变化趋势, 2010年开始增强, 年际变化与太阳活动及地磁活动变化趋势具有较好的对应关系; VTEC变化与太阳活动存在很好的相关性, 相关系数达到0.5以上, 地磁活动则显示了弱相关的特性; F_10.7与VTEC的相关性随着每天Kp指数总值Σ_kp的增大而减小.      

Stability of the seasonal variations in diurnal and semidiurnal components of mid-latitude F2 layer parameters

We report work utilizing 15-min resolution ionospheric data obtained with DPS-4 digisonde in 2003–2011 to study the seasonal variations in amplitudes and phases of the most powerful spectral components of the F2 layer critical frequency (foF2) and peak height (hmF2) fluctuations over Irkutsk (52.5°N, 104.0°E). We show that fluctuations of both parameters contain quasi-harmonic components with periods of T_n = 24/n h (n = 1–7). The number of distinct spectral peaks varies from 3 in summer to 7 in winter. Amplitude and phase characteristics of the diurnal (n = 1) and semidiurnal (n = 2) components is studied using the data sets extracted from the original data sets with band-pass filter. It has been found that the amplitudes of diurnal/semidiurnal foF2 and diurnal hmF2 components are maximum in winter and minimum in summer. Amplitudes of the diurnal components vary gradually; those of the foF2 semidiurnal one, abruptly, thus forming a narrow winter maximum in November–January. The phase (local time of maximum) of the diurnal foF2 component increases gradually by 4–6 h from winter to summer. The phase of the semidiurnal foF2 component is nearly stable in winter/summer and sharply decreases (increases) by 2–3 h near the spring (autumn) equinox. The phase of the diurnal component of hmF2 (local time of minimum) varies slightly between 1130 and 1300 LT; that of the semidiurnal one decreases (increases) by 4–6 h from January to March (from September to November). The results obtained show that the main features of seasonal variations in the diurnal and semidiurnal components of the mid-latitude F2 layer parameters recur consistently during the solar activity growth and decline phases.     

Time variations of the total electron content in the Southeast Asian equatorial ionization anomaly for the period 2006–2011

This study presents the time variations of the total electron content in the South East Asian equatorial ionization anomaly. The time variation of the TEC is analyzed through the period 2006–2011 by using a latitudinal chain of GPS stations extending in the northern and southern hemisphere. The data shows that the shape of the diurnal variation of the TEC depends on the latitude: a plateau is observed at the stations near the equator and a Gaussian at the station distant from the equator. We observe a semiannual pattern in all the stations with maxima at equinox. In both hemispheres, the amplitude of the crest is larger in spring than autumn from 2006 to 2008 and smaller in spring than in autumn from 2009 to 2011. We also observe an asymmetry between the amplitude and the position of the two crests of ionization. There is a very high level of correlation between the amplitude of the TEC at the two crests and the sunspot number: ∼0.88. During the deep solar minimum 2008–2009, the amplitude of crests of ionization becomes small during several months in summer and winter. The results show that both crests move significantly equatorward in winter than other seasons and there is a tendency for both crests to appear earlier in winter and later in summer.     

Assessment of improvement of the IRI model over Ethiopia for the modeling of the variability of TEC during the period 2013–2016

This paper discusses the monthly and seasonal variation of the total electron content (TEC) and the improvement of performance of the IRI model in estimating TEC over Ethiopia during the solar maximum (2013–2016) phase employing as reference the GPS derived TEC data inferred from four GPS receivers installed in different regions of Ethiopia; Assosa (geog 10.05°N, 34.55°E, Geom. 7.01°N), Ambo (8.97°N, 37.86°E, Geom. 5.42°N), Nazret (8.57°N, 39.29°E, Geom. 4.81°N) and Arba Minch (6.06°N, 37.56°E, Geom. 2.62°N). The results reveal that, in the years 2013–2016, the highest peak GPS-derived diurnal VTEC is observed in the March equinox in 2015 over Arba Minch station. Moreover, both the arithmetic mean GPS-derived and modelled VTEC values, generally, show maximum and minimum values in the equinoctial and June solstice months, respectively in 2014–2015. However, in 2013, the minimum and maximum arithmetic mean GPS-derived values are observed in the March equinox and December solstice, respectively. The results also show that, even though overestimation of the modelled VTEC has been observed on most of the hours, all versions of the model are generally good to estimate both the monthly and seasonal diurnal hourly VTEC values, especially in the early morning hours (00:00–03:00?UT or 03:00–06:00?LT). However, it has also been shown that the IRI 2007 and IRI 2012 versions generally perform best in matching the diurnal GPS derived TEC values as compared to that of the IRI 2016 version. In addition, the IRI 2012 version with IRI2001 option for the topside electron density shows the highest overestimation of the VTEC as compared to the other options. None of the versions of the IRI model are proved to be able to capture the effects of geomagnetic storms.     

Low solar activity variability and IRI 2007 predictability of equatorial Africa GPS TEC

Diurnal, seasonal and latitudinal variations of Vertical Total Electron Content (VTEC) over the equatorial region of the African continent and a comparison with IRI-2007 derived TEC (IRI-TEC), using all three options (namely; NeQuick, IRI01-corr and IRI-2001), are presented in this paper. The variability and comparison are presented for 2009, a year of low solar activity, using data from thirteen Global Positioning System (GPS) receivers. VTEC values were grouped into four seasons namely March Equinox (February, March, April), June Solstice (May, June, July), September Equinox (August, September, October), and December Solstice (November, December, January). VTEC generally increases from 06h00 LT and reaches its maximum value at approximately 15h00–17h00 LT during all seasons and at all locations. The NeQuick and IRI01-corr options of the IRI model predict reasonably well the observed diurnal and seasonal variation patterns of VTEC values. However, the IRI-2001 option gave a relatively poor prediction when compared with the other options. The post-midnight and post-sunset deviations between modeled and observed VTEC could arise because NmF2 or the shape of the electron density profile, or both, are not well predicted by the model; hence some improvements are still required in order to obtain improved predictions of TEC over the equatorial region of the Africa sector.     

Behavior of the equivalent slab thickness over three European stations

The total electron content (TEC) derived from the global positioning system (GPS) and the F2-layer peak electron density obtained from Digisonde data have been used to study the diurnal, seasonal and solar activity variations of the ionospheric equivalent slab thickness (τ) over three European stations located at Pruhonice (50.0°N, 15.0°E), Ebro (40.8°N, 0.5°E) and El Arenosillo (37.1°N, 353.3°E). The diurnal variation of the τ is characterized by daytime values lower than nighttime ones for all seasons at low solar activity while daytime values larger than nighttime characterizes the diurnal variation for summer at high solar activity. A double peak is noticeable at dusk and at dawn, better expressed for winter at low solar activity. The seasonal variations of τ depend on local time and solar activity, the daytime values of τ increases from winter to summer whereas nighttime values of τ show the opposite. The effect of the solar activity on τ depends on local time and season, there being very sensitive for winter nighttime values of τ. The results of this study are compared with those presented by other authors.     

Variations of the O(D) and O(S) peak volume emission rates without direct solar effects

Airglow volume emission rates of the O(¹D) red line at 630.0 nm and the O(¹S) green line at 557.7 nm were measured by the Wind Imaging Interferometer (WINDII) on the Upper Atmospheric Research Satellite (UARS) during 1991–1997. Focus of this study is on the peak volume emission rates of the two airglows after removing the direct solar effect, which are referred to as the ‘dark’ peak emission rates. The main results are as follows. For the red line emission, at low and mid-latitudes the daytime variation does not have a clear pattern except an enhancement at dusk; during nighttime there is an enhancement in the equatorial region at 20–03 h, which has a semiannual variation with maxima at equinoxes; at solstices the daytime O(¹D) dark emission rate is stronger in winter than in summer. For both the green line E-and F-layers the distribution of the dark peak volume emission rate is symmetric about noon in all seasons, symmetric about the equator at equinoxes, and stronger in summer than in winter. The O(¹S) E-layer is profoundly affected by tides. For the first time the diurnal and semidiurnal amplitudes for the emission rates are derived using 24-h zero-sun data. The amplitude of the diurnal tide can be as large as 20% of the mean peak volume emission rate, and has maxima at the equator and about 40°N/S, and minima at about 20°N/S. The daily diurnal maximum is at noon at the equator but at midnight at 40°N/S. There is a clear semiannual variation of the diurnal amplitude in the equatorial region with maxima at equinoxes. The amplitude of the semidiurnal tide is mostly less than 10% of the mean peak volume emission rate with maximum amplitudes at noon and midnight. There is an annual variation of the semidiurnal amplitude at mid-latitudes peaking in summer. Aurorae appear in all three emission layers day and night. The green aurorae are brighter than the red aurorae, and the green E-layer aurorae are 2–3 times stronger than the F-layer aurorae. The green aurora has a clear midday gap in the F-layer and an afternoon gap in the E-layer. The red aurorae are particularly strong in the so-called cusp region at equinoxes.     

Transient variations of vertical total electron content over some African stations from 2002 to 2012

This paper presents the vertical total electron content vTEC variations for three African stations, located at mid-low and equatorial latitudes, and operating since more than 10 years. The vTEC of the middle latitude GPS station in Alexandria, Egypt (31.2167°N; 29.9667°E, geographic) is compared to the vTEC of two others GPS stations: the first one in Rabat/Morocco (33.9981°N; 353.1457°E, geographic), and the second in Libreville/Gabon (0.3539°N; 9.6721°E, geographic). Our results discussed the diurnal, seasonal, and solar cycle dependences of vTEC at the local ionospheric conditions, during different phases of solar cycle in the light of the classification of Legrand and Simon. The vTEC over Alexandria exhibits the well-known equinoctial asymmetry which changes with the phases of the solar cycle; the spring vTEC is larger than that of autumn during the maximum, decreasing and minimum phases of solar cycle 23. During the increasing phase of solar cycle 24, it is the contrary. The diurnal variation of the vTEC presents multiple maxima during the equinox from 2005 to 2008 and during the summer solstice from 2006 to 2012. A nighttime vTEC enhancement and winter anomaly are also observed. During the deep solar minimum (2006–2009) the diurnal variation of the vTEC observed over Alexandria is similar to the diurnal variation observed during quiet magnetic period at equatorial latitudes. We observed also that the amplitude of vTEC at Libreville is larger than the amplitude of vTEC observed at Alexandria and Rabat, indeed Libreville is near the southern crest of the Equatorial Ionization anomaly. Finally, the correlation coefficient between vTEC and the sunspot number Rz is high and changes with solar cycle phases.     

Statistical ionospheric E layer properties measured with the Cyprus Digisonde and comparisons with IRI predictions

This paper reports the diurnal, seasonal, and long term variability of the E layer critical frequency (foE) and peak height (hmE) derived from Digisonde measurements from 2009 to 2016 at the low-middle latitude European station of Nicosia, Cyprus (geographical coordinates: 35°N, 33°E, geomagnetic lat. 29.38°N, I = 51.7°). Manually scaled monthly median values of foE and hmE are compared with IRI-2012 predictions with a view to assess the predictability of IRI. Results show that in general, IRI slightly overestimates foE values both at low and high solar activity. At low solar activity, overestimations are mostly limited to 0.25?MHz (equivalent electron density, 0.775?×?10³?el/m^?3) but can go as high as 0.5?MHz (equivalent electron density, 3.1?×?10³?el/m^?3, during noon) around equinox. In some months, underestimations, though sporadic in nature, up to 0.25?MHz are noted (mostly during sunrise and sunset). At high solar activity, a similar pattern of over-/underestimation is evident. During the entire period of study, over-/under estimations are mostly limited to 0.25?MHz. In very few cases, these exceed 0.25?MHz but are limited to 0.5?MHz. Analysis of hmE reveals that: (1) hmE remains almost constant during ±2 to ±4?h around local noon, (2) hmE values are higher in winter than in spring, summer and autumn, (3) there are two maxima near sunrise and sunset with a noontime minimum in between. During the entire period of study, significant differences between observed hmE and the IRI predictions have been noted. IRI fails to predict hmE and outputs a constant value of 110?km, which is higher than most of the observed values. Over- and under estimations range from 3 to 13?km and from 0 to 3?km respectively.     

Analysis of electron content variations over Japan during solar minimum: Observations and modeling

Comparative analysis of GPS TEC data and FORMOSAT-3/COSMIC radio occultation measurements was carried out for Japan region during period of the extremely prolonged solar minimum of cycle 23/24. COSMIC data for different seasons corresponded to equinox and solstices of the years 2007–2009 were analyzed. All selected electron density profiles were integrated up to the height of 700 km (altitude of COSMIC satellites), the monthly median estimates of Ionospheric Electron Content (IEC) were retrieved with use of spherical harmonics expansion. Monthly medians of TEC values were calculated from diurnal variations of GPS TEC estimates during considered month. Joint analysis of GPS TEC and COSMIC data allows us to extract and estimate electron content corresponded to the ionosphere (its bottom and topside parts) and the plasmasphere (h > 700 km) for different seasons of 2007–2009. Percentage contribution of ECpl to GPS TEC indicates the clear dependence from the time and varies from a minimum of about 25–50% during day-time to the value of 50–75% at night-time. Contribution of both bottom-side and topside IEC has minimal values during winter season in compare with summer season (for both day- and night-time). On average bottom-side IEC contributes about 5–10% of GPS TEC during night and about 20–27% during day-time. Topside IEC contributes about 15–20% of GPS TEC during night and about 35–40% during day-time. The obtained results were compared with TEC, IEC and ECpl estimates retrieved by Standard Plasmasphere–Ionosphere Model that has the plasmasphere extension up to 20,000 km (GPS orbit).     

Ionospheric behavior of the F2 peak parameters foF2 and hmF2 at Hainan and comparisons with IRI model predictions

Monthly median values of foF2, hmF2 and M(3000)F2 parameters, with hourly time interval resolution for the diurnal variation, obtained with DPS-4 digisonde observations at Hainan (19.4°N, 109.0°E) are used to study the low latitude ionospheric variation behavior. The observational results are compared with the International Reference Ionospheric Model (IRI) predictions. The time period coverage of the data used for the present study is from March 2002 to February 2005. Our present study showed that: (1) In general, IRI predictions using CCIR and URSI coefficients follow well the diurnal and seasonal variation patterns of the experimental values of foF2. However, CCIR foF2 and URSI foF2 IRI predictions systematically underestimate the observed results during most time period of the day, with the percentage difference ΔfoF2 (%) values changing between about −5% and −25%, whereas for a few hours around pre-sunrise, the IRI predictions generally overestimate the observational ones with ΔfoF2 (%) sometimes reaching as large as ∼30%. The agreement between the IRI results and the observational ones is better for the year 2002 than for the other years. The best agreement between the IRI results and the observational ones is obtained in summer when using URSI coefficients, with the seasonal average values of ΔfoF2 (%) being within the limits of ±10%. (2) In general, the IRI predicted hmF2 values using CCIR M(3000)F2 option shows a poor agreement with the observational results. However, when using the measured M(3000)F2 as input, the diurnal variation pattern of hmF2 given by IRI2001 has a much better agreement with the observational one with the detailed fine structures including the pre-sunrise and post-sunset peaks reproduced reasonably well. The agreement between the IRI predicted hmF2 values using CCIR M(30,000)F2 option and the observational ones is worst for the afternoon to post-midnight hours for the high solar activity year 2002. During daytime hours the agreement between the hmF2 values obtained with CCIR M(30,000)F2 option and the observational ones is best for summer season. The discrepancy between the observational hmF2 and that obtained with CCIR M(30,000)F2 option stem from the CCIR M(3000)F2 model, which does not produce the small scale structures observed in the measured M(3000)F2.     

Semi-annual variations of chlorophyll concentration on the Eastern Tropical Pacific coast of Mexico

An analysis of the main remotly sensed oceanographic variables was conducted to characterize the dominant drivers that modulate the spatial-temporal variability of coastal phytoplankton biomass in the northern limit of the Eastern Tropical Pacific. The phytoplankton biomass was analyzed using monthly average satellite chlorophyll (Chlo) concentration data from MODIS sensors for the period 2003–2017. Oceanographic conditions were analyzed by considering (i) Sea Surface Temperature (SST) high-resolution data from the GHRSST project, (ii) wind stress calculated with data from the zonal and meridional components of the CCMP product, (iii) climatological precipitation, and (iv) climatological river flow. Chlo time series and spatial variability were analyzed using Hovmöller diagrams and Empirical Orthogonal Functions (EOF), respectively. A strong semi-annual signal in Chlo concentration along the coast was observed: the first peak occurs in winter-spring (5.0 mg·m⁻³) and a second one in summer-autumn (6.5 mg·m⁻³). A high year-round average of Chlo concentration (3.0–15.0 mg·m⁻³) was maintained by a 10 km wide along-shore coast, with maximum values spatially associated with river mouths, of which Santiago river registered the highest Chlo values (20.0 mg·m⁻³). Surface oceanographic conditions showed a marked annual cycle with warmer conditions (30–31 °C) from July to October and colder ones (23–24 °C) from December to April; during the cold period, seasonal wind stress stimulated coastal upwelling. The EOF showed that 70% of Chlo concentration variability was controlled by the semi-annual pattern, which responded to coastal upwelling conditions during cold period, while rivers outflows influenced high Chlo concentration during the warm period. These results highlighted the importance for land-ocean interface to sustain coastal ecosystems' biological production and the major role of watersheds as sources of nutrients to maintain high biological production during warm periods in transitional tropical-subtropical zones.     

Variability and climatology of precipitable water vapor from 12-year GPS observations in Taiwan

Because of global warming, global sea levels have risen, the frequency of drought in Taiwan is much more frequent in winter and spring, and rainfall tends to concentrate in summer. The probability of disaster-type weather has also increased significantly. Estimating precipitable water vapor (PWV) through GPS signals, related studies and analyses of weather conditions, and the effective use of meteorological forecasts have been valued by many meteorological research organizations and officials. In this study, PWV data from 2006 to 2017 and rainfall data were used for long-term harmonic analysis. PWV data calculated by ECMWF (ECMWF-PWV) and PWV data calculated by GPS (GPS-PWV) were subjected to regression analysis to verify the reliability of the GPS-PWV data. The research results show that GPS-PWV and ECMWF-PWV have extremely high correlations; however, the climatic characteristics of some regions and the high spatial resolution of GPS-PWV are able to accurately calculate the high topographic relief of small areas. It is judged that the GPS-PWV is more accurate than the ECMWF-PWV. It is worth noting that the PWV trend of the regions during the 6-year-before period has not changed very much, but the rainfall trend has changed obviously. Except for the eastern region, most of the regions show a decreasing trend year by year. More long-term observations are still needed to prove whether this phenomenon relates to global warming. Long-term rainfall analysis showed that the topography blocked water vapor to the western, southern, and mountainous regions, making them distinctly wet or dry. The harmonic curve showed great consistency with the peaks of PWV and rainfall. However, in the northern and eastern parts of the windward side, the time when maximum rainfall occurred each year may be one month later than the time when the maximum PWV value occurred each year. The reason for this difference is likely to be a decrease in the number of autumn typhoons, resulting in a nearly one-month difference in PWV peaks and rainfall peaks. Finally, we analyzed the linear trend of GPS-PWV and temperature for all regions in Taiwan, and found that annual increasing rate of GPS-PWV and temperature of all regions are within 0.4–0.5 mm/year and 0.04–0.11 C°/year, respectively.