On the Correlation of Solar Activity and Troposphere on the GNSS/EGNOS Integrity. Fuzzy Logic Approach

: There is a noticeable number of studies carried out on both the operational parameters of Global Navigation Satellite System (GNSS) and the satellite signal itself. Researchers look for, inter alia, proven sources of errors affecting the integrity of the satellite signal because this parameter determines the approval of the system’s operational use. It also seems of key importance that the atmospheric conditions, in any area of satellite signal usage, should not be underestimated due to their extensive impact. As the ionospheric refraction seriously limits the operational use of the satellite navigation signal, in this article, the authors attempted to quantify the effect of solar activity (expressed by sunspots) on the signal integrity using fuzzy logic. Fuzzy reasoning is used when information is inaccurate or incomplete and necessitates making decisions under conditions of uncertainty. Thanks to fuzzy sets, there are no obstacles to characterize the degree of intensity of a given phenomenon. In order to look at the problem more broadly, attention was also paid to the tropospheric conditions, and it was veriﬁed whether, against the background of cloudiness, precipitation, humidity, pressure and temperature, solar activity affects the integrity to the greatest extent. The integrity measurements from the EGNOS system (PRN120 and PRN126) collected at the monitoring station in Warsaw, Poland in 2014 were used.


Introduction
An interest in satellite systems applicable to various engineering fields, especially to navigation, has been observed for a long time. Indeed, navigation systems based on GNSS (Global Navigation Satellite System) have their practical application in many areas, e.g., geodesy, transport, and robotics. It may be written that these systems have an unparalleled advantage in navigational technologies due to their high-precision delivery of location in terms of position, time, and velocity on any object or person all over the world. An example of an exemplification of this feature may be Real-Time Precise Point Positioning (RT-PPP) since it can provide centimeter-to-decimeter-level positioning accuracy by using a single receiver anywhere in the world [1]. The wide application of satellite navigation can be found in transport, where it can be used both as object (vehicle) navigation, but also with toll-collection systems on highways, parking-space monitoring systems, or as a vehicle fleet-management subsystem [2,3]. The sources of satellite navigation application should be sought in maritime transport, in which GNSS applications have been supporting for years the Automatic Identification System (AIS) [4].
However, it turned out that the use of a satellite navigation system entails the need to correct positioning errors caused by the atmosphere or by signal reflection [5]. The article 1. striving to minimize, and if possible eliminate, errors in the propagation of the satellite signal [13]; 2.
studying the influence of atmospheric conditions on the satellite signal [16,22,23].
Not surprisingly, model analysis is predominantly used in these trends, and also nowadays, artificial neural networks are used.
In this article, the authors attempt to quantify the effect of solar activity (expressed by sunspots) on satellite navigation signal integrity, using fuzzy logic. However, in order to look at the problem more broadly, attention was also paid to other atmospheric (tropospheric) conditions, and in this case, unlike other research in the literature [24][25][26][27], it was verified whether, against the background of cloudiness, precipitation, humidity, pressure and temperature, solar activity affects the integrity to the greatest extent. Moreover, research usually deals with accuracy rather than signal integrity. For this purpose, the integrity measurements from the EGNOS (PRN120 and PRN126) system collected at the monitoring station in Warsaw, Poland, in 2014 were used. The research contribution of this article focuses on several aspects. Firstly, studies of the effects of the solar activity on satellite signal integrity were conducted. Secondly, the use of a fuzzy approach was considered. Thirdly, result verification with a statistical analysis that shows the direct correlations between variables was undertaken. Finally, tropospheric factors were added to the research, including not only humidity, pressure and temperature, but also cloudiness and precipitation, which is not found in the traditional approach. The article has the following structure: introduction and review of state-of-the-art research (1), wider presentation of the research problem (2), materials and method devoted to fuzzy logic (3), results and discussion (4), verification (5), conclusions (6), literature.

Research Problem
The operational limitation of satellite systems usage can be caused by various errors that are inherent in the operation of these systems. In general, errors in satellite systems can be divided into four groups: signal propagation errors (ionospheric errors, tropospheric errors, multipath errors); relativistic errors; system operation errors (satellite ephemeris errors, satellite clock errors); receiver errors, including DOP (Dilution of Precision) [1].
On the way from the satellite to the receiver antenna, the satellite signal passes through the atmosphere, which affects the signal velocity. From the radio wave propagation perspective, the Earth's atmosphere can be divided into the neutral part-troposphere, and the ionized part-ionosphere. The state of the ionosphere, in the context of satellite systems, is determined by Total Electron Content (TEC): the total number of electrons concentrated between a satellite and a receiver along a column with a cross-section equal to 1 m 2 . TEC is directly related to the ionospheric delay of satellite signals (1). The delay results from a change in the group and phase velocity of the wave and depends on the frequency of the signal. Since the ionosphere is a dispersive medium for GNSS frequencies, the ranging codes are delayed and the carrier phases are advanced [28,29].
where: ∆ Ion is the GNSS signal delay; f is the signal frequency. The TEC level is primarily driven by the solar radiation, mainly ultraviolet and x-rays. Consequently, the state of the ionosphere is driven by solar activity, which is manifested by changes on the surface of the sun and its atmosphere. The level of solar activity is related to the changes in the distribution and number of sunspots (SSN). A sunspot is a darker area on sun surface, which is a source of increased electromagnetic radiation. The solar sunspot cycle is approximately 11 years, and in one cycle, the sun changes its activity from a minimum to the next minimum [30]-see Figure 1.
The study of TEC, its modeling and prediction has become a popular problem in modern science in the context of progress in satellite technology. Due to the proven significant influence of TEC on the quality of the received satellite signal, much attention has been focused on TEC prediction, e.g., by using artificial neural networks [31][32][33][34]. Another important field of study is radio occultation technique, which is used both for studying ionosphere, space weather, and neutral atmosphere [35,36].  The troposphere, on the other hand, is the lowest and the thinnest layer of the atmosphere. It extends from the Earth's surface to a height of about 40 km [1]. All the most important phenomena that shape the weather and climate on Earth take place in the troposphere. For GNSS frequencies the troposphere is a non-dispersive medium, so both code measurements and carrier phases are delayed, regardless of the signal frequency. Tropospheric error is caused by both the dry and wet parts of the troposphere. The dry (hydrostatic) part consists of atmospheric pressure, air density and air temperature (2). The wet (non-hydrostatic) part of the troposphere is related to the water vapor content (3). This content is due to weather patterns in the troposphere, for example, the presence of clouds and the precipitation introduce almost 99% of moisture in the lowest part of the atmosphere. The water vapor content, thus, also depends on altitude and atmospheric phenomena (storms, fronts) [38].
where: and is the tropospheric delay for the dry and wet component, respectively; is the elevation angle at the observing site; is the atmospheric pressure; is the temperature; is the partial pressure of water vapor. While integrity-according to the International Civil Aviation Organization (ICAO), is defined as a measure of confidence in the correctness of information provided by the system, it includes the ability of the system to deliver appropriate warnings (alerts) to the user at a predetermined time, which contain information on when the system should not be used. In the context of integrity modeling, in order to determine if the position error is acceptable, an alert limit is set so that it is possible to reflect the maximum allowable position error that will not jeopardize the ongoing operation using a satellite signal (e.g., aircraft take-off). It should be remembered that navigation using satellite systems, and thereby, the signal, is provided simultaneously to many objects over a large area-often one or several continents. The impact of the satellite system integrity loss could be much greater than with conventional navigation methods. Sufficiently, information provided about the loss of signal integrity should result in the resignation from the use of satellite navigation or the cessation of operations (e.g., take-off or landing of an aircraft, moves of The troposphere, on the other hand, is the lowest and the thinnest layer of the atmosphere. It extends from the Earth's surface to a height of about 40 km [1]. All the most important phenomena that shape the weather and climate on Earth take place in the troposphere. For GNSS frequencies the troposphere is a non-dispersive medium, so both code measurements and carrier phases are delayed, regardless of the signal frequency. Tropospheric error is caused by both the dry and wet parts of the troposphere. The dry (hydrostatic) part consists of atmospheric pressure, air density and air temperature (2). The wet (non-hydrostatic) part of the troposphere is related to the water vapor content (3). This content is due to weather patterns in the troposphere, for example, the presence of clouds and the precipitation introduce almost 99% of moisture in the lowest part of the atmosphere. The water vapor content, thus, also depends on altitude and atmospheric phenomena (storms, fronts) [38].  is the tropospheric delay for the dry and wet component, respectively; E is the elevation angle at the observing site; p is the atmospheric pressure; T is the temperature; e is the partial pressure of water vapor.
While integrity-according to the International Civil Aviation Organization (ICAO), is defined as a measure of confidence in the correctness of information provided by the system, it includes the ability of the system to deliver appropriate warnings (alerts) to the user at a predetermined time, which contain information on when the system should not be used. In the context of integrity modeling, in order to determine if the position error is acceptable, an alert limit is set so that it is possible to reflect the maximum allowable position error that will not jeopardize the ongoing operation using a satellite signal (e.g., aircraft take-off). It should be remembered that navigation using satellite systems, and thereby, the signal, is provided simultaneously to many objects over a large area-often one or several continents. The impact of the satellite system integrity loss could be much greater than with conventional navigation methods. Sufficiently, information provided about the loss of signal integrity should result in the resignation from the use of satellite navigation or the cessation of operations (e.g., take-off or landing of an aircraft, moves of ship or robots). Moreover, it is an individual and unique feature of satellite navigation to adapt the navigation capability to the changing constellation of satellites. The impact of changes in the space segment may be increased by an additional fault in the ground segment, for example, damage of one of the components [39].
GNSS augmentation systems are also susceptible to the influence of atmospheric effects. Rapid changes in the ionosphere influence signal delays that cannot be corrected in real time. Such cases are rare, and the likelihood of their occurrence varies by region, but cannot be excluded. As a result, errors that arise may distort information and should be taken into account at the system design stage. Determining the requirements for the parameters of the satellite signal and the implementation of augmentation systems are the key processes for introducing satellite navigation [40,41]. In Europe, the turning point in this activity was the launch of the SBAS-EGNOS (Satellite Based Augmentation System-Geostationary Navigation Overlay Service) system [15,42]. It also reaching the milestone in Poland in implementing aviation procedures based on GNSS navigation.
The EGNOS system is a European GPS and GLONASS augmentation system in the fields of air, maritime and road transport. It was designed by the ETG (European Tripartite Group) group, which includes: the European Space Agency ESA, the European Commission EC and the European Organization for the Safety of Air Navigation EUROCON-TROL. The aim of EGNOS is to monitor the integrity of GPS and GLONASS and to increase their accuracy by introducing data corrections. The principle of operation is based on the reception of GNSS signals by ground reference stations [19]. The EGNOS terrestrial segment includes, among others, RIMS (Ranging and Integrity Monitoring Stations) stations. RIMS can be analyzed in terms of the channels constituting them (A, B, and C). Data from channel A is used to calculate data needed for EGNOS messages. Data from the B channel is used to verify the message calculated from the A channel. Data from the C channel is to detect errors in the signal provided by the GNSS satellites. The task of RIMS is to collect data from GNSS satellites and verify them. By dividing RIMS stations into three channels, the EGNOS data waveform is presented in two separate but related cycles, which, according to the minimum fault warning requirements, should take a total duration of no more than 6 s. The first cycle is the processing cycle. It uses only channel A. The course of this cycle consists of measuring the pseudorange to the satellites through the antennas and receivers of the A channel of the RIMS station, then transmission of the received data to the MCC (Mission Control Center) cell: CPF (Central Processing Unit), then generating corrections, integrity information and the entire navigation message that is sent by the NLES (Navigation Land Earth Stations) stations to the EGNOS geostationary satellites. The second cycle is the checking one. The B and C channels are used in here. Data collected by RIMS antennas and receivers from EGNOS satellites and other navigation satellite systems are sent to the MCC-CPF cell. In the CPF, these data are checked and verified for accuracy, correctness, truthfulness and integrity of information. In case of detection of an error in emitted information, the so-called alarm flag informing system recipients of an error or failure is sent.
Integrity for EGNOS is presented in the form of a Stanford diagram. The diagram shows boundaries of the relationship between PE (Position Error), PL (Protection Level) and AL (Alert Limit). The PE level should always be lower than AL. A favorable relationship for integrity also occurs when PE < PL. It can, therefore, be concluded that data integrity will be ensured when PE < PL < AL. The relationship PL < PE < AL may also be considered relatively positive, however, when PL < PE, then information distortions may occur [19,43]. The distribution of the Stanford diagram is shown in Figure 2.
Owing to the measurements of the EGNOS satellite system signal carried out by, among others, the Polish Air Navigation Services Agency, it is possible to analyze this signal on the basis of its daily parameters. The agency uses a specialized receiver (Septentrio PolaRx3) and appropriate software (PEGASUS 4.7.5.) that can decode the EGNOS data. PEGASUS, developed for research purposes, determines the difference between the position identified by the EGNOS system and the known position of the EGNOS receiver. It is then possible to obtain a positioning error and a number of other important quantities. This software was developed to test the performance of SBAS and GBAS (Ground Based Augmentation System) systems in the first place. For each of the geostationary satellites, PEGASUS automatically generates a report once a day. The integrity is interpreted by the number of misleading information (MI) and the horizontal (HPL-Horizontal Protection Level and HPE-Horizontal Position Error) and vertical (VPL-Vertical Protection Level and VPE-Vertical Position Error) statistics in terms of HPL and VPL standard deviations (std), mean values, and 50%, 95% and 99% significance levels [2]. Graphical visualization of integrity is presented via Stanford diagrams. Owing to the measurements of the EGNOS satellite system signal carried out by, among others, the Polish Air Navigation Services Agency, it is possible to analyze this signal on the basis of its daily parameters. The agency uses a specialized receiver (Septentrio PolaRx3) and appropriate software (PEGASUS 4.7.5.) that can decode the EGNOS data. PEGASUS, developed for research purposes, determines the difference between the position identified by the EGNOS system and the known position of the EGNOS receiver. It is then possible to obtain a positioning error and a number of other important quantities. This software was developed to test the performance of SBAS and GBAS (Ground Based Augmentation System) systems in the first place. For each of the geostationary satellites, PEGASUS automatically generates a report once a day. The integrity is interpreted by the number of misleading information (MI) and the horizontal (HPL-Horizontal Protection Level and HPE-Horizontal Position Error) and vertical (VPL-Vertical Protection Level and VPE-Vertical Position Error) statistics in terms of HPL and VPL standard deviations (std), mean values, and 50%, 95% and 99% significance levels [2]. Graphical visualization of integrity is presented via Stanford diagrams.
The analysis and interpretation of PEGASUS reports are particularly important for the purposes of this article. As a result of the agreement and official cooperation with the Polish Air Navigation Services Agency in the conducted research, the data from 2014 generated for the Warsaw monitoring station for the EGNOS PRN120 and PRN126 satellites operating in this period were used. Note that the last solar maximum took place in 2014.

Materials and Methods
The analysis of phenomena related to the satellite signal is not possible without model analysis. Experimenting on the actual operation can lead to dramatic consequences. Therefore, numerous attempts are made to study the signal and the processes taking place in satellite systems during operation, which have a potential impact on its disturbance.
Fuzzy sets are a suitable tool for assessing the influence of solar activity on the integrity of a satellite signal. A fuzzy set is a mathematical object with a defined membership function that takes values from the range [0,1]. The fuzzy set theory was introduced by Lotfi A. Zadeh in 1965 as an extension of the classical set theory [45]. This reasoning is used when information is inaccurate or incomplete, which necessitates a decision under The analysis and interpretation of PEGASUS reports are particularly important for the purposes of this article. As a result of the agreement and official cooperation with the Polish Air Navigation Services Agency in the conducted research, the data from 2014 generated for the Warsaw monitoring station for the EGNOS PRN120 and PRN126 satellites operating in this period were used. Note that the last solar maximum took place in 2014.

Materials and Methods
The analysis of phenomena related to the satellite signal is not possible without model analysis. Experimenting on the actual operation can lead to dramatic consequences. Therefore, numerous attempts are made to study the signal and the processes taking place in satellite systems during operation, which have a potential impact on its disturbance.
Fuzzy sets are a suitable tool for assessing the influence of solar activity on the integrity of a satellite signal. A fuzzy set is a mathematical object with a defined membership function that takes values from the range [0,1]. The fuzzy set theory was introduced by Lotfi A. Zadeh in 1965 as an extension of the classical set theory [45]. This reasoning is used when information is inaccurate or incomplete, which necessitates a decision under uncertainty. Knowledge about the phenomenon is subjective and is expressed through the opinion or descriptive research of an expert. It is informally referred to as linguistic uncertainty. Their big advantage is the possibility of graphical presentation of research results, modeling uncertainty and examining the influence of many factors on the phenomenon. There are numerous references to the use of fuzzy sets in the literature for navigation, e.g., [46][47][48]. For the sake of order, the process of building a linguistic fuzzy model is presented below. Consequently, a universal, simplified course of action was created, based on the available literature on the subject [45,[49][50][51], enabling the construction of a comprehensive fuzzy model. In the process of its building, it is possible to distinguish several stages. Stage 1. Definition of a fuzzy set. Fuzzy set A in the space of considerations X = {x}, meaning set A in X is defined as a set of pairs: where: 1] is the degree of belonging of an element x ∈ X in fuzzy set A.
The membership function µ A : X → [0, 1] to each element x ∈ X assigns its degree of belonging to the fuzzy set A. Three cases of belonging can be distinguished: Stage 2. Fuzzification. The selection of the membership function shape takes place before defining the linguistic area of considerations. The mathematical notation of Gauss function is presented below (5). Gaussian function is one of the commonly used membership functions where x is the center and σ is the width of the Gaussian curve. Also, this curve was chosen for researches carried out in this paper.
In Table 1, the first two stages: fuzzy set definition, determination of linguistic variables and membership function identification are presented. Table 1 refers to the research problem presented in this paper where inputs are identified by sunspots, number of daily sunspot observations and standard deviation of daily sunspot number. The outputs are identified by integrity of PRN120 and PRN126. Even though we do not expect the number of sunspot observations and related standards deviation to have any impact on the signal integrity, we added these parameters to test/validate our approach.   Each curve represents the degree of belonging of a given linguistic variable to a function. The x-axis (horizontal) specifies a range of values for each variable, e.g., the number of sunspots ranges from 40 to 220. The y-axis values, in turn, range from 0 to 1 and represent the severity of the phenomenon, e.g., the value of 100 sunspots is small with probability 0.5, medium with probability 0.7, and large with probability 0.  Each curve represents the degree of belonging of a given linguistic variable to a function. The x-axis (horizontal) specifies a range of values for each variable, e.g., the number of sunspots ranges from 40 to 220. The y-axis values, in turn, range from 0 to 1 and represent the severity of the phenomenon, e.g., the value of 100 sunspots is small with probability 0.5, medium with probability 0.7, and large with probability 0.
Stage 3. Inference. The inference is served by the base of rules, which is created after the process of assigning the membership function to linguistic variables. Two types of architectures can be used in this step: Mamdani or Takagi-Sugeno [48,49]. The architecture of Mamdani was chosen in this research. Rules may be written as follows: where: , -linguistic variables; , -values of linguistic variables (low, moderate, high); ∆ -control variable (output variable). There are different sources of creating the base of rules. Among them, the following may be distinguished: knowledge obtained from experts, scientific experience, literature sources, other research tools. The rule base, according to the selected architecture, is used for inference, which consists of determining the fuzzy value of the output variable, knowing the set of rules representing the relationships between the input and output variables and some fuzzy values of the input variables.
Stage 4. Defuzzification. Having obtained a fuzzy value of the output variable, it should be defuzzified, that is, the corresponding strict (non-fuzzy) value should be found. where: x, y-linguistic variables; PM, PB-values of linguistic variables (low, moderate, high); ∆u-control variable (output variable).
There are different sources of creating the base of rules. Among them, the following may be distinguished: knowledge obtained from experts, scientific experience, literature sources, other research tools. The rule base, according to the selected architecture, is used for inference, which consists of determining the fuzzy value of the output variable, knowing the set of rules representing the relationships between the input and output variables and some fuzzy values of the input variables.
Stage 4. Defuzzification. Having obtained a fuzzy value of the output variable, it should be defuzzified, that is, the corresponding strict (non-fuzzy) value should be found. Below, a mathematical formula representing a center of gravity method is presented (7). It is one of the most commonly used for defuzzification process. It is assumed that there is a fuzzy set A identified in X = {x 1 , x 2 , . . . , x n }. A non-fuzzy value is searched for a ∈ [x 1 , x n ], which best represents A.
Shall it be assumed that a non-fuzzy value is a ∈ X = {x 1 , . . . , x n }, for which µ A (x) takes the maximum value.
Building a model and conducting research is subject to appropriate data. In order to build the fuzzy model, measurements from year 2014 and the Warsaw station for PRN120 and PRN 126 satellites were used together with data on solar activity, including the number of sunspots, the number of observations, and SSN standard deviation from The Space Research Center of the Polish Academy of Sciences. In total, 181 rules of inference were distinguished.

Results and Discussion
As an initial result of the conducted research, it is interesting to see that there is some coincidence between the number of the SSN observations and the integrity of the satellite signal. In case of both, the PRN120 and PRN126 satellites, a large number of SSN observations correlate positively with the signal integrity ( Figure 5).
Below, a mathematical formula representing a center of gravity method is presented (7). It is one of the most commonly used for defuzzification process. It is assumed that there is a fuzzy set identified in = { 1 , 2 , … , }. A non-fuzzy value is searched for ∈ [ 1 , ], which best represents .
Shall it be assumed that a non-fuzzy value is ∈ = { 1 , … , }, for which ( ) takes the maximum value.
Building a model and conducting research is subject to appropriate data. In order to build the fuzzy model, measurements from year 2014 and the Warsaw station for PRN120 and PRN 126 satellites were used together with data on solar activity, including the number of sunspots, the number of observations, and SSN standard deviation from The Space Research Center of the Polish Academy of Sciences. In total, 181 rules of inference were distinguished.

Results and Discussion
As an initial result of the conducted research, it is interesting to see that there is some coincidence between the number of the SSN observations and the integrity of the satellite signal. In case of both, the PRN120 and PRN126 satellites, a large number of SSN observations correlate positively with the signal integrity ( Figure 5).
It can be also observed that changes in the standard deviation of SSN measurement correlate with the signal integrity ( Figure 6). Of course, this comes from obvious correlation between the number of measurements carried out, and the related standard deviation of the measured parameter.  It can be also observed that changes in the standard deviation of SSN measurement correlate with the signal integrity ( Figure 6). Of course, this comes from obvious correlation between the number of measurements carried out, and the related standard deviation of the measured parameter.
More interesting results come from SSN analysis. The correlation between the number of sunspots and signal integrity for both EGNOS satellites is presented in Figure 7. Interestingly, the signal in the analyzed period obtained the highest integrity when the solar activity, measured by the number of sunspots, reached an average value in the range of 101-160. It is also interesting that low SSN is correlated with low signal integrity. However, low solar activity is connected with more frequent occurrence of the coronal holes, which, in turn, are sources of high-speed solar wind that affect Earth's magnetosphere [30]. And this results in adverse ionospheric effects that influence signal integrity. On the other hand, high SSN is connected with a high TEC level and greater chances for intense geomagnetic storms, which, clearly, have negative effects [22]. This shows that SSN is not necessarily the best metric to be used in the context of GNSS signal integrity. Therefore, our future studies will include analysis of more directly related factors, such as TEC and geomagnetic indices. More interesting results come from SSN analysis. The correlation between the number of sunspots and signal integrity for both EGNOS satellites is presented in Figure  7. Interestingly, the signal in the analyzed period obtained the highest integrity when the solar activity, measured by the number of sunspots, reached an average value in the range of 101-160. It is also interesting that low SSN is correlated with low signal integrity. However, low solar activity is connected with more frequent occurrence of the coronal holes, which, in turn, are sources of high-speed solar wind that affect Earth's magnetosphere [30]. And this results in adverse ionospheric effects that influence signal integrity. On the other hand, high SSN is connected with a high TEC level and greater chances for intense geomagnetic storms, which, clearly, have negative effects [22]. This shows that SSN is not necessarily the best metric to be used in the context of GNSS signal integrity. Therefore, our future studies will include analysis of more directly related factors, such as TEC and geomagnetic indices.  More interesting results come from SSN analysis. The correlation between the number of sunspots and signal integrity for both EGNOS satellites is presented in Figure 7. Interestingly, the signal in the analyzed period obtained the highest integrity when the solar activity, measured by the number of sunspots, reached an average value in the range of 101-160. It is also interesting that low SSN is correlated with low signal integrity. However, low solar activity is connected with more frequent occurrence of the coronal holes, which, in turn, are sources of high-speed solar wind that affect Earth's magnetosphere [30]. And this results in adverse ionospheric effects that influence signal integrity. On the other hand, high SSN is connected with a high TEC level and greater chances for intense geomagnetic storms, which, clearly, have negative effects [22]. This shows that SSN is not necessarily the best metric to be used in the context of GNSS signal integrity. Therefore, our future studies will include analysis of more directly related factors, such as TEC and geomagnetic indices.  The last stage of the research is the process of defuzzification, i.e., searching for the non-fuzzy values of the variables that give the best result (highest integrity)- Figure 8. The combination from Figure 8 was used as an example for the first 20 rules. A relatively high integrity value was achieved under certain input conditions: x 1 =127, x 2 = 9.88, x 3 =20.6. (a) (b) re 6. The SSN standard deviation (x2) and the integrity of PRN120 (a) and PRN126 (b) (Own study in MATLAB Fuzzy c Toolbox).
More interesting results come from SSN analysis. The correlation between the number of sunspots and signal integrity for both EGNOS satellites is presented in Figure 7. Interestingly, the signal in the analyzed period obtained the highest integrity when the solar activity, measured by the number of sunspots, reached an average value in the range of 101-160. It is also interesting that low SSN is correlated with low signal integrity. However, low solar activity is connected with more frequent occurrence of the coronal holes, which, in turn, are sources of high-speed solar wind that affect Earth's magnetosphere [30]. And this results in adverse ionospheric effects that influence signal integrity. On the other hand, high SSN is connected with a high TEC level and greater chances for intense geomagnetic storms, which, clearly, have negative effects [22]. This shows that SSN is not necessarily the best metric to be used in the context of GNSS signal integrity. Therefore, our future studies will include analysis of more directly related factors, such as TEC and geomagnetic indices.   The observations of behavior of the fuzzy model led to interesting conclusions, resulting from the research-the highest integrity of the satellite signal was observed for the "medium" value of the number of sunspots (in the range of 101-160). There is also surprisingly high correlation between measurements of SSN (number and std) and the signal integrity.

Materials and Method
The linguistic fuzzy model presented in the previous point must be verified. This verification is necessary to confirm, in a documented manner, that the models actually lead to the presented results. Moreover, this stage is particularly important when the studied phenomenon, and especially the research results, may have an impact on the safety of navigation of people, vehicles or robots. At the same time, in order not to generalize the results, the authors proposed a statistical verification of the direct influence of solar activity (defined by the number of sunspots) on the integrity of the satellite signal against the background of other atmospheric factors. Five variables were selected: cloudiness, relative air humidity, precipitation, temperature and atmospheric pressure.
The verification was performed using Hellwig's integral capacity method, also known as the optimal predicant selection method [52]. This is a formal technique for selecting the best explanatory variables. This technique is considered to be optimal because, as a result of verification, explanatory (input) variables with the largest impact on the explained (output) variable may be determined. Furthermore, the rules of Hellwig's verification method are as follows: • it allows to eliminate the explanatory variables that do not affect the explained variable; • explanatory (input) variables should be weakly related to each other, which prevents duplication of information provided by them; • a large difference between the number of observations and the number of explanatory variables positively affects the accuracy of the predicate selection.
At the source of Hellwig's method lies the mathematical statistics and the correlation calculus. Its essence and, at the same time, the result reflect in a combination of explanatory variables where the integral information capacity is the largest. The entire calculation procedure can be summarized in a few steps: 1. the concept of the information carrier, which is the explanatory variable (input variable), is introduced; 2. the number of all possible combinations of explanatory variables is determined according to the formula: where m is the number of explanatory variables; 3. all types of combinations between explanatory variables are determined; 4. the individual capacity of information h kj is calculated according to the formula: where: k-in h is the number of combinations k = 1, 2, . . . , K; j-in h specifying the number of the explanatory variables distinguished in a given combination remaining in dependencies with the explained variable; r j -correlation coefficient of the j-th explanatory variable with the explained variable; r ij -correlation coefficient of i-th and j-th explanatory variable (i, j = 1, 2, . . . , m); 6. the optimal combination of explanatory variables H k opt. -the one for which the integral information capacity is the largest is selected: For the purposes of calculations, the values of new variables were obtained in the tested time-series. The summary in Table 2 presents a measure of each variable and range of measurement values from the given period and station. It is worth emphasizing that the selection of these variables for comparison was not accidental. In the literature, the role of the tropospheric factor in shaping the satellite signal path is repeatedly emphasized. For example, the tropospheric water vapor content is a derivative parameter from zenith wet delay (ZWD) and can be calculated as a function of the refractive index, depending on temperature, pressure and relative humidity [25][26][27]. Selbesoglu [24] claims that the troposphere affects GNSS signals due to the fast variability of the refraction index. Determination of the wet delay is more difficult than the hydrostatic delay due to much faster changes in the water vapor in the troposphere.

Results and Discussion
As presented in Table 2, the following explanatory variables will take part in the calculations: x 5 -atmospheric pressure; • x 6 -solar activity.
The explained variables will be an integrated integrity from the satellite signal (y). Therefore, the number of possible combinations of explained variables according to the Formula (8) is: All types of combinations between explanatory variables are shown in Table 3. It is possible to calculate the individual capacities of the h kj information carrier after prior determination of correlation coefficients. Firstly, the correlation coefficients r j , j-th explanatory variable with the explained variable. Secondly, the coefficients of correlation r ij , i-th and j-th explanatory variable. The coefficient determining the level of linear where: Specialists [53,54] have interpreted its range, but it should not be treated too closely. For the needs of these studies, the below assumptions were taken. In this, if r = 0, variables are generally described as correlated. And moreover: The absolute value of the correlation coefficient indicates the strength of statistical interdependence of variables, understood as a degree of determination of one variable value by the value of the second variable. It is worth mentioning that when the regression lines have a positive slope, this is unambiguous with a positive correlation of variables. On the other hand, when the regression lines have a negative slope, this indicates that the variables are negatively correlated and the correlation coefficient takes a negative value. At the same time, it should be noted that: and the other coefficients can be considered in the same way. In Table 4 comparison of correlation coefficients for r j (explanatory variable with the explained variable) is shown. In Table 5 comparison of correlation coefficients for r ij (i-th and j-th explanatory variable) is shown. The next stage of Hellwig's method is the calculation of the capacity of the individual information carrier h kj . A large number of combinations between explanatory variables generates a complex process of determining all individual capacities. For order, the calculation example below shows the individual capacities of the two selected combinations K 17 , K 53 : h 173 = r 2 y 1 x 3 1 + |r x 3 x 5 | = 0.059 (16) h 175 = r 2 y 1 x 5 1 + |r x 3 x 5 | = 0.021 (17) h 532 = r 2 y 1 x 2 1 + |r x 2 x 3 | + | r x 2 x 4 | + |r x 2 x 6 | = 0.0 (18) h 533 = r 2 h 534 = r 2 y 1 x 4 1 + |r x 2 x 4 | + |r x 3 x 4 | + |r x 4 x 6 | = 0.008 (20) h 536 = r 2 y 1 x 6 1 + |r x 2 x 6 | + |r x 3 x 6 | + |r x 4 x 6 | = 0.041 (21) The fifth step of the presented method is to calculate the integral capacities of the information carrier H k for individual combinations of variables. For previously selected combinations, the formula for determining the value of H k is presented below. This formula refers to the combination of K 17 and takes the form: And for K 53 presents as follows: H 53 = h 532 + h 533 + h 534 + h 536 = 0.095 (23) The last, but the key, step of the method is to select the optimal combination of explanatory variables H k opt. . The list of H k opt. values for each combination is shown in Table 6. And, as a result of the calculations carried out, the maximum value of H k opt. was determined for K 40 (signed with yellow in Table 6).
H k opt. = 0.123 (24) The largest value of the integral information capacity was related to the combination K 40 , in which x 3 , x 5 , x 6 explanatory variables performed. Variable x 3 corresponds to the precipitation, while x 5 corresponds to pressure. Variable x 6 is responsible for solar activity. This suggests that even though the wet delay is relatively small (10-20 cm), its fast changes associated with passages of weather-fronts may have a negative effect on signal integrity.

Conclusions
The purpose of this article was to attempt to quantify the effect of solar activity (expressed by sunspots) and tropospheric parameters on the satellite navigation signal integrity using fuzzy logic. Due to the fuzzy model, the highest integrity of the satellite signal was observed for the "medium" value of the number of sunspots (in the range 101-160). Integrity seemed to be very sensitive to changes in solar activity manifested by a significant decrease or a significant increase in the number of sunspots. What is more, attention was also paid to other atmospheric (tropospheric) conditions and, in this case, it was verified whether, against the background of cloudiness, precipitation, humidity, pressure and temperature, solar activity affects the integrity to the greatest extent. Interestingly, despite the fact that the tropospheric water vapor content (which is a derivative parameter from ZWD) and its components are considered in the context of the impact on the GNSS signal, it was actually the precipitation which played an important