3D Galileo Reference Antenna Pattern for Space Service Volume Applications

Abstract

There is an increasing demand for navigation capability for space vehicles. The exploitation of the so-called Space Service Volume (SSV), and hence the extension of the Global Navigation Satellite System (GNSS) from terrestrial to space users, is currently considered a fundamental step. Knowledge of the constellation antenna pattern, including the side lobe signals, is the main input for assessing the expected GNSS signal availability and navigation performance, especially for high orbits. The best way to define and share this information with the final GNSS user is still an open question. This paper proposes a novel methodology for the definition of a high-fidelity and easy-to-use statistical model to represent GNSS constellation antenna patterns. The reconstruction procedure, based on antenna characterization techniques and statistical learning, is presented here through its successful implementation for the “Galileo Reference Antenna Pattern (GRAP)” model, which has been proposed as the reference model for the Galileo programme. The GRAP represents the expected Equivalent Isotropic Radiated Power (EIRP) variation for the Galileo FOC satellites, and it is obtained by processing the measurements retrieved during the characterization campaign performed on the Galileo FOC antennas. The mathematical background of the model is analyzed in depth in order to better assess the GRAP with respect to different objectives such as improved resolution, smoothness and proper representation of the antenna pattern statistical distribution. The analysis confirms the enhanced GRAP properties and envisages the possibility of extending the approach to other GNSSs. The discussion is complemented by a preliminary use case characterization of the Galileo performance in SSV. The accessibility, a novel indicator, is defined in order to represent in a quick and compact manner, the expected Galileo SSV quality for different altitudes and target mission requirements. The SSV characterization is performed to demonstrate how simply and effectively the GRAP model can be inserted into user analysis. The work creates the basis for an improved capability for assessing Galileo-based navigation in SSV according to the current knowledge of the antenna pattern.

1. Introduction

Use of global navigation satellite systems (GNSSs) in the Space Service Volume (SSV), and particularly for certain missions (MEO, GEO, HEO, Lunar, and beyond), requires knowledge about the GNSS transmission antenna patterns, or a representative model thereof, including side lobes [1]. Different studies of navigation in high orbit ([2,3]) rely on the knowledge of the received signal power, which is instrumental for a wide range of link budget analysis used for spaceborne receiver design, navigation performance assessment and implementation of filtering techniques of the GNSS observables [4,5,6,7]. The best way to retrieve and share this information with the final GNSS user is still an open research issue. The availability of the antenna pattern measurements at system level can be sparse and prone to limitation of the testing facilities. In fact, different pattern realizations can be experienced, even with the same antenna design. This generally corresponds to a burden for users, who still need to aggregate the data ([2,3]) in order to mitigate the discontinuities and measurement degradation effects. This paper proposes a novel approach aiming to achieve the same objectives by considering a more general representation, which preserves all the features of the constellation antenna pattern with improved statistical significance. This derives from a data-driven model reconstruction, providing antenna pattern expected mean and variance distributions. Such an approach is here implemented and adopted for the definition of the Galileo Reference Antenna Pattern (GRAP), for which metadata have recently been published in [8]. The proposed reconstructed model is an efficient and compact tool to provide the GNSS user and community with the most stable representation of the expected Galileo FOC transmitting antenna pattern, according to the knowledge matured at system level before the deployment of the service. The paper aims at providing a more in-depth presentation of the methodologies put in place for the model development, and showing how it can be easily exploited to perform a preliminary characterization of the expected Galileo SSV.

The first contribution is a detailed description of the GRAP model synthesis [8]. The procedure ingests the information retrieved during the antenna characterization campaign performed on the Galileo FOC satellite up to the highest available off-boresight angles, before their integration on the satellites. Instead of just pooling all observations to achieve a better statistical representation of the expected system performance, an ad-hoc multi-step procedure has been developed. The approach aims at optimally aggregating the available antenna pattern information using methodologies derived from antenna characterization techniques [9,10,11,12,13,14] and statistical learning [15]. The first step consists of a frequency resample of the power distribution. After that, a reconstruction algorithm based on representation of the pattern in a sparse spherical harmonic basis is considered. The algorithm exploits an elastic net regularization [16,17] in order to improve the smoothness of the 3D spatial distribution and obtain its representation in sparse spherical harmonic parameter space. The final step introduces the concept of representing the constellation antenna pattern as an hierarchical population model [18]. The reconstructed patterns from different data collections are addressed as different realizations of the same model and merged, in order to achieve a general and easy-to-use representation of the antenna for the final GNSS space user. One of the most relevant features of the proposed statistical process is the possibility to define the GRAP values and their correspondent confidence bounds, so the idea is to complement the antenna pattern metadata [8] with information representative of the uncertainty level associated with the different pattern realizations. This level of confidence can be reflected in the user’s analysis; hence, their conclusions inherit the capability to assess what will likely happen according to current knowledge of the pattern. For the sake of brevity and without loss of generality, all the results presented hereafter refer to the Galileo E1 carrier frequency, but extension to other Galileo frequencies can be easily derived by using the metadata provided in [8] and adjusting all frequency-dependent parameters.

The second contribution of this paper is to process the developed GRAP model into a preliminary use case characterization of the Galileo SSV performance. Such an assessment is based on a simplified SSV geometrical model allowing easy definition of a target user altitude, the GRAP region spanned, and the contribution to the overall received signal power. More importantly, the proposed SSV analysis introduces a novel compact Key Performance Indicator (KPI), namely the accessibility. The introduced KPI expresses the effective arc of a rising GNSS Space Vehicle (SV) that can be exploited for navigation purposes, considering a target received power or receiver sensitivity. The results show that this parameter, even though limited by the selected geometry and signal propagation assumptions, provides a new, competitive and compact performance index to assess antenna pattern-driven GNSS SSV peculiarities.

The paper is organized as follows. Section 2 defines the general constellation antenna pattern reconstruction procedure and its mathematical framework. The resulting 3D GRAP is presented in Section 3 and complemented with the analysis of properties deriving from the reconstruction procedure. In Section 4.1, the simplified SSV model and accessibility index are introduced. The information about the GRAP is injected into the SSV model and, in Section 4.2, a Galileo SSV characterization is provided in terms of accessibility charts. Section 5 assesses the general conclusions and identifies the way forward for the current and future improvement of the Galileo SSV.

2. Multi-Step 3D Reference Antenna Pattern Reconstruction

The GRAP model is the result of processing of the dataset made available by the Galileo programme [8] and retrieved from several testing campaigns performed on the Galileo FOC satellite’s isolated antennas, i.e., before their integration on the satellite.

The measurement dataset has been provided as Continuous Wave (CW) spherical scanning of the radiated field for a target antenna $i$. The pattern measurements are indicated as

$${\mathrm{G}}^i\left(\left.\mathsf{\theta },\mathsf{\phi },{\mathrm{f}}_{CW}\right|{\mathsf{\theta }}_1,{\mathsf{\phi }}_1,{\mathrm{f}}_1,{\mathsf{\theta }}_2,{\mathsf{\phi }}_2,{\mathrm{f}}_2,\dots {\mathsf{\theta }}_{\mathrm{n}},{\mathsf{\phi }}_{\mathrm{n}},{\mathrm{f}}_{\mathrm{n}}\right)$$

(1)

where ${\mathrm{G}}^i\in \left[0, 1\right]$ represents the radiated field normalized power in linear scale for the i-th antenna (Figure 1), measured at the different co-elevation ${\theta }_j\in [35,90°]$ and azimuth ${\phi }_k\in [0,360°]$ scanning points for a fixed CW frequency ${\mathrm{f}}_{\mathrm{C}\mathrm{W}}$.

The reference antenna pattern derivation is here intended as a data processing approach aiming to achieve the following three objectives:

(1)

Provide to the user a representation of the pattern applicable to target Galileo signals, hence the need to report the measurement at different CW to the Galileo central frequencies.

(2)

Overcome the facility sampling limitations and achieve a high-resolution 3D antenna pattern with smooth distribution, filtering out discontinuities and measurement degradation.

(3)

Handle the uncertainty introduced by different realizations of the same antenna during different SVs testing sessions by using a statistical representation in terms of expected pattern values and bounds.

The GRAP synthesis procedure has been split into three different steps. In the first step (Step-1 in Figure 2), a Power Spectral Density (PSD) Weighting is applied to the CW antenna measurements in order to resample the power distribution of the modulated signal in the target frequency band:

$${\mathrm{G}}_f^i\left(\left.\mathsf{\theta },\mathsf{\phi }\right|{\mathsf{\theta }}_1,{\mathsf{\phi }}_1,{\mathrm{f}}_1,{\mathsf{\theta }}_2,{\mathsf{\phi }}_1,\dots {\mathsf{\theta }}_{\mathrm{j}},{\mathsf{\phi }}_{\mathrm{k}}\right)$$

(2)

A different GRAP model can then be derived for the selected Galileo central frequencies considering $f=\{E1,E5a,E5b,E6\}$ [8]. Step-2 (Figure 2 ) introduces an improved 3D spatial model reconstruction to overcome scanning facility limitations and obtain an estimation of a smooth model ($m$) of the i-th antenna pattern $G_{f,m}^i\left(\mathsf{\theta },\mathsf{\phi }\right)$ working closer to the expected physical behavior. Such an estimation can be obtained together with corresponding information about pattern spatial uncertainty, expressed by its model variance distribution ($\sigma )$, which we indicate as ${\mathrm{G}}_{f,\sigma }^i\left(\mathsf{\theta },\mathsf{\phi }\right)$. The 3D patterns obtained are then used to feed the reference GRAP constellation model for a target central frequency ${\mathrm{G}}_{f,m}^R\left(\mathsf{\theta },\mathsf{\phi }\right)$ considering a hierarchical approach [18]. Specifically, the different radiation pattern realizations ${\mathrm{G}}_{f,m}^i\left(\mathsf{\theta },\mathsf{\phi }\right)$ relying on a specific space vehicle or test campaign session are processed as observations nested within different subjects at the bottom level, and belonging to the same constellation population at the top level. This defines a regression model allowing generalization of the antenna behavior (features) on top of its specific subject (single antenna patterns) or within subject realizations (measurements for different cuts and angles). The resulting ${\mathrm{G}}_{f,m}^R\left(\mathsf{\theta },\mathsf{\phi }\right)$ allows extension of the model to a more general representation compatible with possible realizations of the pattern family. The residual errors are still considered through estimation of the corresponding variance distribution ${\mathrm{G}}_{f,\sigma }^R\left(\mathsf{\theta },\mathsf{\phi }\right)$. The aggregation procedure can be generally described by (3).

$$\left\{{\mathrm{G}}_{f,m}^1\left(\mathsf{\theta },\mathsf{\phi }\right),{\mathrm{G}}_{f,m}^2\left(\mathsf{\theta },\mathsf{\phi }\right),\dots ,{\mathrm{G}}_{f,m}^n\left(\mathsf{\theta },\mathsf{\phi }\right),{\mathrm{G}}_{f,\sigma }^1\left(\mathsf{\theta },\mathsf{\phi }\right),{\mathrm{G}}_{f,\sigma }^2\left(\mathsf{\theta },\mathsf{\phi }\right),\dots ,{\mathrm{G}}_{f,\sigma }^n\left(\mathsf{\theta },\mathsf{\phi }\right)\right\}\to \left\{{\mathrm{G}}_{f,m}^R\left(\mathsf{\theta },\mathsf{\phi }\right),{\mathrm{G}}_{f,\sigma }^R\left(\mathsf{\theta },\mathsf{\phi }\right)\right\}$$

(3)

Such a procedure corresponds to a specific statistical process, the mathematical implementation of which is described in the next sections.

2.1. Frequency Synthesis Step

The $M$ available radiated field linear normalized power CW measurements at elevation and azimuth scanning points ${\mathsf{\theta }}_{\mathrm{j}},{\mathsf{\phi }}_{\mathrm{k}}$ are combined as per the following:

$${\left.{\mathrm{G}}_f^i\right|}_{{\mathsf{\theta }}_{\mathrm{j}},{\mathsf{\phi }}_{\mathrm{k}}}=\sum _{m=1}^M{W}_{{∆f}_m}^{PSD}{\mathrm{G}(\mathsf{\theta }}_{\mathrm{j}},{\mathsf{\phi }}_{\mathrm{k}})$$

(4)

The normalized weighting $W_{{∆f}_n}^{PSD}$ corresponds to the PSD of the E5a, E5b, and E6-B/C signals provided at intervals ∆f of 2.5 MHz. Uniform weighting is applied for the E1-B/C signal given the narrowband nature of this signal compared to the CW sampling interval.

The distribution of $W_{{∆f}_n}^{PSD}$ is represented with respect to frequency in Figure 3. The E5AltBOC combination is not considered hereafter, since the processing of the E5a and E5b patterns have been carried out separately [8]. The procedure moves the different CW samples to the target discrete spatial distribution ${\left.{\mathrm{G}}_f^i\right|}_{{\mathsf{\theta }}_{\mathrm{j}},{\mathsf{\phi }}_{\mathrm{k}}}$, that can then be used to feed the reconstruction procedure to obtain the specific subject ${\mathrm{G}}_{f,m}^i\left(\mathsf{\theta },\mathsf{\phi }\right)$ pattern realization for the selected frequency $f$.

2.2. A 3D Optimal Spherical Harmonic Based Spatial Reconstruction

The antenna spherical field characterization generally leads to a complex analysis, the results of which can also be far from the most effective realization. The common approach is to characterize the pattern by using a proper spherical harmonic expansion. As per [9,11] and [12,13], the reconstructed pattern can be expressed by using the following spherical harmonic $SH\left(\mathsf{\theta },\mathsf{\phi }\right)$ expansion up to the order N

$$SH\left(\mathsf{\theta },\mathsf{\phi }\right)=\sum _{l=0}^N\sum _{m=-l}^l{Y}_{lm}\left(\mathsf{\theta },\mathsf{\phi }\right){\beta }_{lm}$$

(5)

with

$$Y_{lm}\left(\mathsf{\theta },\mathsf{\phi }\right)=\left\{\begin{array}{c}{\left(-1\right)}^m\sqrt{2}\sqrt{\frac{2l+1}{4\pi }\frac{\left(l-\left|m\right|\right)!}{\left(l+\left|m\right|\right)!}}{P}_l^{\left|m\right|}\left(cos\theta \right)\sin \left(\left|m\right|\mathsf{\phi }\right) \mathrm{i}\mathrm{f} \mathrm{m}<0\\ \sqrt{\frac{2l+1}{4\pi }}{P}_l^{\left|m\right|}\left(cos\theta \right) \mathrm{i}\mathrm{f} \mathrm{m}=0\\ {\left(-1\right)}^m\sqrt{2}\sqrt{\frac{2l+1}{4\pi }\frac{\left(l-m\right)!}{\left(l+m\right)!}}{P}_l^m\left(cos\theta \right)\cos \left(\left|m\right|\mathsf{\phi }\right) \mathrm{i}\mathrm{f} \mathrm{m}<0\end{array}\right.$$

(6)

where ${\beta }_{lm}$ indicates the ${\left(N+1\right)}^2$ expansion coefficients and $P_l^m$ are the corresponding Legendre polynomials. The reconstruction problem can then be expressed in matrix form as

$$G_{j,k}^i=\mathit{Y}\left({\mathsf{\theta }}_{\mathrm{j}},{\mathsf{\phi }}_{\mathrm{k}}\right){\mathit{\beta }}^i+{\nu }_{j,k}^i$$

(7)

where $\mathit{Y}$ is the spherical harmonic matrix operator computed by evaluating the spherical harmonic functions defined in (6) in the available observation points. The value $G_{j,k}^i$ indicates the i-th antenna measurements ${\left.G_f^i\right|}_{{\theta }_j,{\phi }_k}$ with spatial error ${\nu }_{j,k}^i$.

The standard least square solution can be poorly conditioned as far as higher harmonic order N is considered. Popular regularization strategies and basis pursuit favors physically admissible spherical fields with smooth variations, and allows better exploitation of the inner sparsity of the selected spherical harmonic feature space [9]. The application of a Sparse Estimator can select the relevant antenna model components [11] according to the data sampling and the level of noise. This work considers an elastic-net [15] variant, also proposed in [16], hence the target pattern parameter set can be computed as

$$\widehat{{\mathsf{\beta }}_{\mathit{i}}}=\underset{\mathsf{\beta }}{\mathit{min}}\sum _{\mathrm{j}=1}^{\mathrm{n}}{\left\{G_{j,k}^i-{\mathit{Y}}^{\mathrm{*}}\left({\mathsf{\theta }}_{\mathrm{j}},{\mathsf{\phi }}_{\mathrm{k}}\right){\mathit{\beta }}_{\mathit{i}}\right\}}^{\mathrm{T}}\left\{G_{j,k}^i-{\mathit{Y}}^{\mathrm{*}}\left({\mathsf{\theta }}_{\mathrm{j}},{\mathsf{\phi }}_{\mathrm{k}}\right){\mathit{\beta }}_{\mathit{i}}\right\}+\mathrm{n}\mathsf{\lambda }\left[{(1-\mathsf{\alpha })‖{\mathit{\beta }}_{\mathit{i}}‖}_2^2+{\mathsf{\alpha }‖{\mathit{\beta }}_{\mathit{i}}‖}_1\right]$$

(8)

where the loss function has been slightly modified considering

$${\mathit{Y}}^{\mathbf{*}}\left({\mathsf{\theta }}_{\mathrm{j}},{\mathsf{\phi }}_{\mathrm{k}}\right)={\Gamma }^{-1/2}\mathit{Y}\left({\mathsf{\theta }}_{\mathrm{j}},{\mathsf{\phi }}_{\mathrm{k}}\right)$$

(9)

where $\Gamma$ is a weighting variance matrix iteratively estimated from the regression residuals in order to approximate the normalized linear pattern log-normal distribution [19].

A balance between smoothing and shrinkage has been considered (α = 0.5), whereas the $\lambda$ key hyper-parameter has been found by using a Bayesian minimization of the k-fold cross-validation error (with k = 10) [15]. The cross-validation splits the dataset of the pattern into training and validation subsets. The optimal lambda ${\lambda }_{min}$ will be one that minimize the loss function computed by solving the regression on the training set and evaluating model predictions against the validation data. For the antenna pattern problem, this corresponds to learning the harmonic content of a subset of data and optimizing the capability of the model to predict another validation section (i.e., angular sector). The maximum order considered for $SH\left(\mathsf{\theta },\mathsf{\phi }\right)$ is $N=60$, so the maximum feature space $\widehat{{\mathit{\beta }}_{\mathit{i}}}$ considers 3721 components. The hyper-parameter is derived once per pattern family, i.e., E1 frequency patterns, and then included to solve (8). The learning curve for the E1 frequency pattern family is shown in Figure 4: the process achieves a minimum ${\lambda }_{min}$ for which the model does not learn much more by adding non-zero ${\beta }_{lm}$ coefficients. The curve can be assumed as a reference example of the optimization process for the other frequency cases $E5a,E5b,E6$ as well, which converge to their specific hyper-parameter values.

The target regressor problem in (8) aims at estimating not only the reconstruction coefficient $\widehat{{\mathit{\beta }}_{\mathit{i}}}$, providing the 3D reconstruction:

$$G_{f,m}^i\left(\theta ,\phi \right)=\mathit{Y}\left(\theta ,\phi \right)\widehat{{\mathit{\beta }}_{\mathit{i}}}$$

(10)

but also the associated covariance matrix ${\Sigma }_{\widehat{{\mathit{\beta }}_{\mathit{i}}}}$ allowing approximation of the expected spatial uncertainty in different pattern regions by using the SH operator:

$$G_{f,\sigma }^i\left(\theta ,\phi \right)=\mathit{Y}\left(\theta ,\phi \right){\Sigma }_{\widehat{{\mathit{\beta }}_{\mathit{i}}}}{\mathit{Y}\left(\theta ,\phi \right)}^{\mathit{T}}$$

(11)

A bootstrap approach [15] has been considered here to better estimate ${\Sigma }_{\widehat{{\mathit{\beta }}_{\mathit{i}}}}$ for the regressor exploited in (8). Residual resampling performed one hundred times is sufficient to converge to a good approximation, in line with the small deviation experienced from the data.

2.3. Hierarchical Constellation Model

The constellation pattern hierarchical model considers two representation layers:

• The individual models of the specific antenna pattern ${\mathrm{G}}_{f,m}^i\left(\mathsf{\theta },\mathsf{\phi }\right)$.
• The population average model or the constellation pattern ${\mathrm{G}}_{f,m}^R\left(\mathsf{\theta },\mathsf{\phi }\right)$.

The method of obtaining the individual estimate of ${\mathrm{G}}_{f,m }^i$ is the reconstruction procedure defined in the previous section, with the possibility of mapping the problem from the spatial domain to the parameter space ${\mathit{\beta }}_{\mathit{i}}$. In such a space, the population model can be expressed by a simplified normally distributed population ${\mathsf{\beta }}_{\mathit{i}}~N\left({\mathsf{\beta }}_{\mathit{\mu }},{\mathsf{\Sigma }}_{{\mathsf{\beta }}_{\mathit{i}}}\right)$, where individual specific realizations are expressed as

$${\mathsf{\beta }}_{\mathit{i}}={\mathsf{\beta }}_{\mathit{\mu }}+{\mathit{b}}_{\mathit{i}} E\left({\mathit{b}}_{\mathit{i}}\right)={\mathsf{\mu }}_{\mathbf{b}}=0 Var\left({\mathit{b}}_{\mathit{i}}\right)=\mathbf{D}$$

(12)

which is a linear combination of the population grand mean coefficients ${\mathsf{\beta }}_{\mathit{\mu }}$ and ${\mathit{b}}_{\mathit{i}}$, referred to as between-subject random effects [18]. $\mathbf{D}$ defines the between-subject associated covariance, and ${\mathsf{\mu }}_{\mathbf{b}}$ is their mean. Considering that the first reconstruction step provides only an estimate of parameters, (12) can be rewritten as

$$\widehat{{\mathsf{\beta }}_{\mathit{i}}}\approx {\mathsf{\beta }}_{\mathit{\mu }}+{\mathit{b}}_{\mathit{i}}+{\mathit{\epsilon }}_{\mathit{i}}^{\mathit{*}} Var\left({\mathit{\epsilon }}_{\mathit{i}}^{\mathit{*}}\right)={\mathsf{\Sigma }}_{\widehat{{\mathsf{\beta }}_{\mathit{i}}}}$$

(13)

where ${\mathit{\epsilon }}_{\mathit{i}}^{\mathit{*}}$ is the reconstruction error and ${\mathsf{\Sigma }}_{\widehat{{\mathsf{\beta }}_{\mathit{i}}}}$ is the uncertainty associated with the reconstruction procedure for the subject i. Equation (13) allows us to estimate the population average parameters ${\mathsf{\beta }}_{\mathit{\mu }}$ by the following weighted least square solution:

$$\widehat{{\mathsf{\beta }}_{\mathit{\mu }}}={\left({\mathsf{{\rm I}}}^{\mathit{T}}{\left(\mathbf{V}\right)}^{-1}\mathbf{I}\right)}^{-1}{\mathsf{{\rm I}}}^{\mathit{T}}{\left(\mathit{V}\right)}^{-1}\widehat{{\mathit{\beta }}_{\mathrm{1,2},\dots L}}{\Sigma }_{{\mathit{\beta }}_{\mathit{\mu }}}={\left({\mathsf{{\rm I}}}^{\mathit{T}}{\left(\mathbf{V}\right)}^{-1}\mathbf{I}\right)}^{-1}$$

(14)

Here, $\widehat{{\mathit{\beta }}_{1,2,\dots \mathit{L}}}$ is the concatenation of dimension LNx1 of vector parameters estimated for all antenna patterns $\mathrm{1,2},\dots L$, $\mathbf{I}$ is a block diagonal matrix with unit matrix blocks $1_{\mathit{N}\mathit{x}\mathit{N}}$, and $\mathit{V}$ is the block diagonal matrix with element ${\mathbf{V}}_i=\mathbf{D}+{\mathsf{\Sigma }}_{\widehat{{\mathsf{\beta }}_{\mathit{i}}}}$. This solution assumes that covariances $\mathbf{D}$ and ${\mathsf{\Sigma }}_{\widehat{{\mathsf{\beta }}_{\mathit{i}}}}$ are known. ${\mathsf{\Sigma }}_{\widehat{{\mathsf{\beta }}_{\mathit{i}}}}$ is provided by the first regression step, whereas $\mathbf{D}$ is computed according to the sample variance estimator:

$$\frac{\Delta {\mathsf{\beta }}_{\mathit{i}}\Delta {{\mathsf{\beta }}_{\mathit{j}}}^T}{N-1}={\left(N-1\right)}^{-1}\sum _{i=1}^N{\left(\widehat{{\mathsf{\beta }}_{\mathit{i}}}-\widehat{{\mathsf{\beta }}_{\mathit{\mu }}}\right)\left(\widehat{{\mathsf{\beta }}_{\mathit{j}}}-\widehat{{\mathsf{\beta }}_{\mathit{\mu }}}\right)}^T i\ne j$$

(15)

and then refined by recursive iterations. Once the population mean $\widehat{{\mathsf{\beta }}_{\mathit{\mu }}}$ and its asymptotic covariance ${\mathsf{\Sigma }}_{\widehat{{\mathsf{\beta }}_{\mathit{\mu }}}}$ are determined, it is possible to define the final population model by averaging (13) with respect to single subject realization, hence:

$$E\left(\widehat{{\mathsf{\beta }}_{\mathit{i}}}\right)=\widehat{{\mathsf{\beta }}_{\mathit{\mu }}} Var\left(\widehat{{\mathsf{\beta }}_{\mathit{i}}}\right)={\mathsf{\Sigma }}_{\widehat{{\mathsf{\beta }}_{\mathit{\mu }}}}+\mathbf{D}+{\mathsf{\Sigma }}_{\widehat{{\mathsf{\beta }}_{\mathit{i}}}}$$

(16)

with ${\mathsf{\Sigma }}_{\widehat{\mathsf{\beta }}}={\left(\mathrm{N}\right)}^{-1}\sum _{\mathrm{i}=1}^{\mathrm{N}}{\mathsf{\Sigma }}_{\widehat{{\mathsf{\beta }}_{\mathit{i}}}}$. The constellation model GRAP can be finally defined as per Section 2.2, in terms of its mean ${\mathrm{G}}_{m,\sigma }^R\left(\mathsf{\theta },\mathsf{\phi }\right)$ and variance ${\mathrm{G}}_{f,\sigma }^R\left(\mathsf{\theta },\mathsf{\phi }\right)$ by simply mapping equation (16) through a spherical harmonic operator:

$$\begin{array}{c}{\mathrm{G}}_{f,m}^R\left(\mathsf{\theta },\mathsf{\phi }\right)=\mathbf{Y}\left(\mathsf{\theta },\mathsf{\phi }\right)\widehat{{\mathsf{\beta }}_{\mathit{\mu }}}\\ {\mathrm{G}}_{f,\sigma }^R\left(\mathsf{\theta },\mathsf{\phi }\right)=\mathbf{Y}\left(\mathsf{\theta },\mathsf{\phi }\right)\left({\mathsf{\Sigma }}_{\widehat{{\mathsf{\beta }}_{\mathit{\mu }}}}+\mathbf{D}+{\mathsf{\Sigma }}_{\widehat{\mathsf{\beta }}}\right){\mathbf{Y}\left(\mathsf{\theta },\mathsf{\phi }\right)}^{\mathit{T}}\end{array}$$

(17)

3. Model Derivation Results

In this section, the GRAP for the E1 frequency is analyzed with respect to the implemented procedure.

The representation in the parameter space is provided in Figure 5. The upper plot provides the E1 GRAP coefficients $\widehat{{\mathsf{\beta }}_{\mathit{\mu }}}$, confirming the good capability of the regression process to emphasize or shrink toward zero the different harmonic components according to their relevance. The bottom plot shows the different covariance contributions defined in (16) ${\mathsf{\Sigma }}_{\widehat{{\mathsf{\beta }}_{\mathit{\mu }}}},\mathbf{D},{\mathsf{\Sigma }}_{\widehat{\mathsf{\beta }}}$ in terms of corresponding standard deviations in the coefficient domain. The three contributions respectively represent the mean model estimation variance (yellow), the between-subject variance (blue), and the specific subject mean variance (red). The proposed hierarchical aggregation in the parameter space allows us to isolate the principal antenna pattern components $\widehat{{\mathsf{\beta }}_{\mathit{\mu }}}$ from the realizations of specific testing campaigns or physical antenna detours, and move them into different variance contributions. This is in line with objective (3) defined in Section 2.

The $\mathit{Y}$ operator allows us to generate the ${\mathrm{G}}_{f,m}^R\left(\mathsf{\theta },\mathsf{\phi }\right)$ and ${\mathrm{G}}_{f,\sigma }^R\left(\mathsf{\theta },\mathsf{\phi }\right)$ for a target point $\left(\mathsf{\theta },\mathsf{\phi }\right)$ according to (17). This possibility fits with objective (2) in Section 2, hence the model can cope with improved resolution with respect to that provided by measured samples.

The reconstructed GRAP information has been made available to the final user in [8], evaluating the pattern function on a 1 deg resolution grid $\mathsf{\theta }\mathrm{x}\mathsf{\phi }$ = [0, 90] × [0, 360], which simultaneously guarantees enough accuracy via standard interpolation methods and provides an efficient lookup table for user implementation.

According to the reference document [8], the GRAP metadata are represented on such a grid in terms of expected Equivalent Isotropic Radiated Power (EIRP) variation with respect to the azimuth and co-elevation angles. This representation can be easily derived from the estimated pattern ${\mathrm{G}}_{m,\sigma }^R\left(\mathsf{\theta },\mathsf{\phi }\right)$, considering

$${EIRP}_{f,GRAP,dBW}^R\left(\mathsf{\theta },\mathsf{\phi }\right)={Pt}_{f,dBW}+{\mathrm{G}}_{f,m,dB}^R\left(\mathsf{\theta },\mathsf{\phi }\right)$$

(18)

where ${\mathrm{G}}_{f,m,dB}^R\left(\mathsf{\theta },\mathsf{\phi }\right)$ is the dB conversion of ${\mathrm{G}}_{f,m}^R\left(\mathsf{\theta },\mathsf{\phi }\right)$ considering the transformation provided in Appendix A, and ${Pt}_{f,dBW}$ is the in-band transmitted power. As per [8], the considered ${Pt}_{f,dBW}$ offset has been computed in order to ensure that the minimum user received power at ground level does correspond to the Galileo commitments for the different signals E5a, E5b, E6, and E1 defined in [20].

The final 3D plot of the ${EIRP}_{E1,m,dBW}^R\left(\mathsf{\theta },\mathsf{\phi }\right)$ and the corresponding polar plot (i.e., with respect to azimuth and radial component) are provided in Figure 6. Polar plot distribution is particularly relevant since it shows that the final user information considers the azimuth asymmetries of the antenna pattern. Those details can be considered as constellation antenna design specific features emerging from the different realizations.

One of the most relevant properties of the reconstruction based on spherical harmonic expansion is to filter out spatial discontinuities and noise effect in the data, moving the distribution closer to its expected physical behavior. This can be appreciated computing the 2D Laplacian of the averaged raw spherical scanning distribution for the target E1 frequency with respect to the developed 3D reconstructed pattern.

Such a comparison is provided in Figure 7, where upper and lower distributions represent 2D Laplacian applied to the raw data sample mean and the reconstructed pattern, respectively. In the latter case, it can be appreciated that the secondary lobes variation will follow a smooth transition from one angular position to the other. This completes GRAP model objective (2) by providing the advantage on the user side of mitigating unexpected discontinuities during GNSS receiver tracking simulations.

According to (17), from the parameter space reconstruction it is also possible to reconstruct the pattern variance ${\mathrm{G}}_{f,\sigma }^R\left(\mathsf{\theta },\mathsf{\phi }\right)$. As per [8], this information can be effectively exploited to define a model confidence bound, which is expressed as user EIRP to represent 95% of the expected variation.

Therefore, the upper ${\left.{EIRP}_{CI,f,dBW}\left(\mathsf{\theta },\mathsf{\phi }\right)\right|}_{+}$ and lower ${\left.{EIRP}_{CI,f,dBW}\left(\mathsf{\theta },\mathsf{\phi }\right)\right|}_{(-)}$ bounds have been defined by using the following expressions:

$${\left.{EIRP}_{CI,f,dBW}\left(\mathsf{\theta },\mathsf{\phi }\right)\right|}_{(+/-)}={EIRP}_{f,GRAP,dBW}^R\left(\mathsf{\theta },\mathsf{\phi }\right)+{2\sigma }_{f,dB}^R\left(\mathsf{\theta },\mathsf{\phi }\right)$$

(19)

Here, ${\sigma }_{f,dB}^R\left(\mathsf{\theta },\mathsf{\phi }\right)$ is derived from the variance matrix ${\mathrm{G}}_{f,\sigma }^R\left(\mathsf{\theta },\mathsf{\phi }\right)$ and converted into dB according to the formulas provided in Appendix A.

Figure 8 shows the final 3D representation of the estimated uncertainty ${2\sigma }_{f,dB}^R\left(\mathsf{\theta },\mathsf{\phi }\right)$. As expected, the uncertainty radially increases towards the extreme ends of the distribution and preserves azimuth diversity due to antenna-specific features.

Evidence that the defined bounds are able to asymptotically represent the data excursion is provided in Figure 9. It can be recognized as a good match with the measurement samples, which fall within the upper and lower bounds for different co-elevation and azimuth cuts.

Figure 9 also provides in the bottom left plot the${3\sigma }_{f,dB}^R\left(\mathsf{\theta },\mathsf{\phi }\right)$ bound in order to appreciate that the model can be generally used to extend the statistical test coverage. The bound slightly overestimates the variability at higher off-boresight angles ($\mathsf{\theta }$> 60), reflecting the complexity of predicting the pattern behavior in the regions that experience very small values and high spatial variation. However, the effect applies only to very small EIRP values, so in terms of power distribution its effect has been considered acceptable, also taking into account that it provides a conservative estimation. It is worth noting that the work performed on the definition of bounds is mainly devoted to complementing the model with compact information for the user about the expected excursion in a specific angular region.

In Figure 10 (left-hand side), the polar plot for ${\sigma }_{f,dB}^R\left(\mathsf{\theta },\mathsf{\phi }\right)<3\mathrm{d}\mathrm{B}$ is provided. A few areas, mainly related to the highest off-boresight angles and around nulls, do not cope with this requirement (white areas), confirming the expected high stability of Galileo FOC constellation patterns across the different realizations. Figure 10 (right-hand side) also provides ${\left.{EIRP}_{CI,f,dBW}\left(\mathsf{\theta },\mathsf{\phi }\right)\right|}_{(-)}>0\mathrm{d}\mathrm{B}$, which can be considered representative of the Galileo E1 EIRP strength. It can be seen that Galileo can guarantee a minimum positive EIRP up to 55 degrees off boresight angles (about 0.6 radii in polar plot), excluding a few localized regions (white areas).

Those kinds of features can be further investigated by using the GRAP metadata provided in [8] for the other Galileo signal carrier frequencies. In the next section, a simplified SSV geometrical analysis is introduced to better connect the GRAP to the final Galileo SSV link budgets and address its exploitation for a target space user.

4. Galileo Pattern Driven Space Service Volume Analysis

4.1. SSV Geometrical Model and Accessibility Index Definition

The Galileo Antenna Reference Pattern information can be adopted for mission feasibility analysis and spaceborne receiver design. This section proposes a first-order analytical procedure allowing preliminary characterization, with respect to a target mission altitude, of the effect of the GRAP EIRP on the expected SSV performance.

The approach is based on a simplified service geometry represented in Figure 11.

Let us consider:

• A target spacecraft rx receiving the Galileo signal and crossing the plane defined by the orbital plane of a reference tx Galileo satellite. Such an event defines a simplified geometry condition where the rx user can be assumed at fixed position ${\mathbf{R}}^{rx}$ and the transmitting Galileo is located at its rising position ${\mathbf{R}}^{tx}$ along its orbit.
• An rx Earth-centred 2D reference frame lying in the Galileo orbital plane and considering y-axis aligned to ${\widehat{\mathbf{R}}}^{rx}$. With this assumption, the 2D rx position vector is ${\mathbf{R}}^{rx}=\left\{0,{\mathrm{r}}_{\mathrm{y}}\right\}$ and depends on the user altitude considering only ${\mathrm{r}}_{\mathrm{y}}={\mathrm{r}}_{\mathrm{E}\mathrm{a}\mathrm{r}\mathrm{t}\mathrm{h}}+{\mathrm{h}}_{rx}$. In Figure 11, a GEO satellite use case is shown, so ${\mathrm{r}}_{\mathrm{y}}$ would be ${\mathrm{r}}_{\mathrm{E}\mathrm{a}\mathrm{r}\mathrm{t}\mathrm{h}}+{\mathrm{h}}_{\mathrm{G}\mathrm{E}\mathrm{O}}$.
• Different crossing events i can be defined so far as the transmitting MEO rising position can be placed at any point in its orbit (circularly approximated). This is represented by different MEO SV positions drawn in Figure 11. Such a sequence of Galileo tx positions can be defined in the RXRF reference frame as ${\mathbf{R}}_i^{tx}=\{{\mathrm{r}}_{\mathrm{x},i}^{tx},{\mathrm{r}}_{\mathrm{y},i}^{tx}\}$ with $\left|{\mathbf{R}}_i^{tx}\right|={\mathrm{r}}_{\mathrm{E}\mathrm{a}\mathrm{r}\mathrm{t}\mathrm{h}}+{\mathrm{h}}_{\mathrm{M}\mathrm{E}\mathrm{O}}$. It is assumed $\widehat{{\mathbf{R}}_i^{tx}}\equiv -\widehat{{\mathbf{n}}_i}$, where ${\mathbf{n}}_i$ is the transmitting antenna boresight nadir-pointing direction. Attitude corrections, as per yaw steering, are outside the scope of the simplified geometry.
• The symmetry of the problem allows us to consider the subset of events where ${\mathrm{r}}_{\mathrm{x},i}^{tx}>0$.

The proposed geometry focuses on the upper SSV layer, where target space vehicle altitudes are higher than the MEO constellation ones, so it is assumed ${\mathrm{h}}_{\mathrm{r}\mathrm{x}}>{\mathrm{h}}_{\mathrm{t}\mathrm{x}}$. This is in line with the standard expected SSV receiving conditions, considering receiver antennas pointing toward the Earth and the GNSS signals generated by satellites rising above the Earth disk. Now, consider the correspondingly defined angles:

${\theta }_i$ the anomaly angle sequence spanned by ${\mathbf{R}}_i^{tx}$ with ${\theta }_i\in \left[0,\pi \right]$, stated ${\mathrm{r}}_{\mathrm{x},i}^{tx}>0$;

${\alpha }_i$ angle defined by ${\mathbf{R}}^{tx}$ and ${\mathbf{R}}^{rx}$ according to the geometrical relation

$${\alpha }_i{=\cos }^{-1}\frac{{\left|{\mathbf{R}}^{tx}\right|}^2+{\left|{\mathbf{R}}^{tx-rx}\right|}^2-{\left|{\mathbf{R}}^{rx}\right|}^2}{2\left|{\mathbf{R}}^{tx}\right|\left|{\mathbf{R}}^{tx-rx}\right|}$$

(20)

where ${\mathbf{R}}^{tx-rx}={\mathbf{R}}^{tx}-{\mathbf{R}}^{rx}$. The proposed domain $\{{\alpha }_i,{\theta }_i\}$ is particularly suitable for linking the pattern region spanned by different relative $rx$ and $tx$ positions considering that each ${\alpha }_i$ corresponds to the antenna pattern off-boresight angles experienced by the GNSS SV for a target transmitting position ${\theta }_i$ (Figure 11, left side pattern representation). It is clear that, according to the altitude of the user ${\mathrm{h}}_{\mathrm{r}\mathrm{x}}$, the sequence ${\theta }_i$ of a rising satellite can cover the whole antenna pattern off-boresight range [0, π/2] without spanning the whole ${\theta }_i$ interval $[0,\pi ]$. Therefore, it is possible to define that a rising transmitter position $\theta$ can be considered compatible only for $\theta : \alpha \le \pi /2$. It is expected that as the altitude increases from GEO to HEO (and above) the $\theta$ compatible interval shrinks according to (20). The interval ${\theta }_i$ is also restricted to [${\theta }_{Earth}$,π] considering the Earth shadowing angle provided by ${\theta }_{Earth}={\sin }^{-1}\left(\left|{\mathbf{R}}_{\mathrm{E}\mathrm{a}\mathrm{r}\mathrm{t}\mathrm{h}}\right|/\left|{\mathbf{R}}^{rx}\right|\right)$.

The compatible set [${\theta }_{Earth}$,${\theta }_{\alpha =\pi /2, h=h_{rx}}$] can be intended as a measure of the geometrical compatibility of the Galileo transmitting positions with respect to the SSV end-user. Figure 12 represents the distribution ${\alpha }_i\left({\theta }_i,h_{rx}\right)$ of the spanned off-boresight angles on a 2D grid. The ${\theta }_i$ step size is 1.25 deg and ${ h}_{rx}$ is considered equally spaced in log-scale and spanning the interval $\left[28\times {10}^6, 390\times {10}^6\right]$ meters, which covers SSV altitudes up to lunar orbits. This $\left\{{\theta }_i,h_{rx}\right\}$ grid will be assumed as a reference to represent all the distributions hereafter considered. In Figure 12, the two competitive effects affecting the ${\theta }_i$ set are clearly highlighted as far as the altitude increase: (1) the reduction of the Earth disk obscuration; (2) the reduction of maximum compatible ${\theta }_i$ corresponding to the maximum ${\alpha }_i=\pi /2$.

Considering this domain $\left\{{\theta }_i,h_{rx}\right\}$, the SSV Received Signal Power ${\mathrm{C}}_{SSV}^R$ can be easily defined to map the experienced off-boresight angles through the GRAP EIRP model and considering the correspondent path loss [21]:

$${\mathrm{C}}_{SSV}^R\left({\theta }_i,\mathsf{\phi },h_{rx}\right)={\mathrm{E}\mathrm{I}\mathrm{R}\mathrm{P}}_{\mathsf{\phi }\mathrm{*}}^R\left({\alpha }_i,\mathsf{\phi }\right)-20log\left(\frac{4\pi \left|{\mathbf{R}}_{h_{rx}}^{tx-rx}\right|}{\lambda }\right)$$

(21)

where$\lambda$ is the E1 carrier signal wavelength, and (21) is intended as a function of the receiver altitude $h_{rx}$ through ${\mathbf{R}}_{h_{rx}}^{tx-rx}$ and ${\theta }_i$ according to the representation defined. Equation (21) can be used to express the system performance for different EIRP azimuth angles spanning $\mathsf{\phi }\in [0,2\pi ]$. In addition, assuming that for the simplified geometry, different $\mathsf{\phi }$ can be equivalently experienced, it is possible to build up synthetic indexes like:

$$\begin{array}{c}{\mathrm{C}}_{SSV,min}^R\left({\theta }_i,h_{rx}\right)=\underset{\mathsf{\phi }\in \left[\mathrm{0,2}\pi \right]}{\min }{\mathrm{C}}^R\left({\theta }_i,\mathsf{\phi },h_{rx}\right)\\ {\mathrm{C}}_{SSV,ave}^R\left({\theta }_i,h_{rx}\right)=\frac{1}{2\pi }{\displaystyle \sum _{\mathsf{\phi }\in \left[\mathrm{0,2}\pi \right]}}{\mathrm{C}}^R\left({\theta }_i,\mathsf{\phi },h_{rx}\right)\end{array}$$

(22)

which define the worst case and average SSV received signal power, respectively. The minimum and average are computed considering the 1 deg resolution discrete grid defined for $\mathsf{\phi }\in \left[\mathrm{0,2}\pi \right]$. The so-defined ${\mathrm{C}}_{SSV}^R\left({\theta }_i,\mathsf{\phi },h_{rx}\right)$ can also be exploited to introduce a novel compact index, namely Galileo SSV accessibility.

This index is intended as a measure of the fraction of the rising Galileo orbit arc that the user can potentially exploit to access an SV signal during a crossing event for a target minimum received signal power threshold ${\mathrm{C}}_{Min}^R$. In fact, considering the crossing event as uniformly distributed with respect to the current position of the rising Galileo tx, such an indicator can be expressed as a percentage $A_{\%}$, considering the ratio between the effective arc ${\theta }_{eff}$ (spanning signal power values ${{\mathrm{C}}^R>\mathrm{C}}_{Min}^R$) and the maximum allowed arc ${\theta }_{tot}$

$${\left.A_{\%}\right|}_C=\frac{{\left.{\theta }_{eff}\right|}_{{{\mathbf{C}}^{\mathit{R}}>\mathbf{C}}_{Min}^{\mathit{R}}}}{{\theta }_{tot}}=\frac{{\left.{\theta }_{eff}\right|}_{{{\mathbf{C}}^{\mathit{R}}>\mathbf{C}}_{Min}^{\mathit{R}}}}{\pi } { with \theta }_{eff}(h_{rx})={\theta }_{max}-{\theta }_{min }$$

(23)

where ${\theta }_{max}$ and ${\theta }_{min }$are respectively the maximum angle ${\theta }_i$ satisfying ${{\mathrm{C}}^R>\mathrm{C}}_{Min}^R$, and the minimum ${\theta }_i$ geometrically compatible at altitude $h_{rx}$. ${\theta }_{tot}=\pi$ can be easily retrieved by the geometry considered herein (Figure 11). For ${\theta }_{max}={\theta }_{\alpha =\pi /2}$, the Equation (23) defines the geometrical accessibility, hence one that can be achieved removing the limitation on the minimum received signal power.

Moving from service to user side, the introduction of a high-gain antenna is generally considered in order to leverage SSV navigation limits. Assuming the geometry in Figure 11, and a receiving antenna pattern equally distributed in azimuth, the receiving antenna gain can be evaluated by

$${\mathrm{G}}_{Rx}\left({\beta }_i\right) \mathrm{with} {\beta }_i{=\cos }^{-1}\frac{{\left|{\mathbf{R}}^{rx}\right|}^2+{\left|{\mathbf{R}}^{tx-rx}\right|}^2-{\left|{\mathbf{R}}^{tx}\right|}^2}{2\left|{\mathbf{R}}^{tx}\right|\left|{\mathbf{R}}^{tx-rx}\right|}$$

(24)

where ${\beta }_i$is the receiving antenna off-boresight angle corresponding to a transmitting position ${\theta }_i$.

As shown in Figure 13 (top), the usage of a high-gain antenna is also justified by considering that the spanned elevation $\frac{\pi }{2}-{\beta }_i$ converges to very high values as the altitude increases. The restricted receiving working region is shown in the bottom of Figure 13 for a representative high-gain antenna model, which will be considered as reference in the rest of the analysis.

It is now possible to introduce in the link budget assessment [21] the receiver sensitivity expressed in terms of carrier to noise ratio C/N0 [dB/Hz]. Such a value, without considering implementation losses [21], can be easily computed from the available signal carrier power ${\mathrm{C}}^R\left({\theta }_i,\mathsf{\phi },h_{rx}\right)$, adding the amplification effect of ${\mathrm{G}}_{Rx}\left({\beta }_i\right)$ from (24) and subtracting the receiver equivalent thermal noise $N_0$

$${C/N_0}_{SSV}^R\left({\theta }_i,\mathsf{\phi },h_{rx}\right)={\mathrm{C}}^R\left({\theta }_i,\mathsf{\phi },h_{rx}\right)+{\mathrm{G}}_{Rx}\left({\left.{\beta }_i\right|}_{{\alpha }_i}\right)-10log10\left(K{T}_{Sys}\right)$$

(25)

where K is the Boltzmann constant and $T_{Sys}$ is the equivalent system temperature in Kelvin [21]. From (25), synthetic performance indexes, i.e., expected average or worst-case conditions, can be derived as per (22).

Calculating the full link budget allows us to also express the accessibility in terms of receiver sensitivity, i.e., the fraction of the rising Galileo orbit arc providing a carrier-to-noise ratio above a predefined threshold ${\left.C/{\mathit{N}}_0\right|}_{Min}$:

$${\left.A_{\%}\right|}_{C/{\mathit{N}}_0}=\frac{{\left.{\theta }_{eff}\right|}_{C/{\mathit{N}}_0>{\left.C/{\mathit{N}}_0\right|}_{\mathit{M}\mathit{i}\mathit{n}}}}{\pi }{ with \theta }_{eff}\left(h_{rx}\right)={\theta }_{max}-{\theta }_{min }$$

(26)

It is worth bearing in mind that $C/N_0$ accessibility based on (25) is intended as a simplified version of the link budget, which is sufficient for a first-order characterization of the SSV general performance. Starting from (25), specific user evaluation can customize the process by considering more complex signal propagation and implementation losses. The receiving antenna pattern can deviate from the assumed reference (Figure 13), and it can also be a function of the azimuth. Those limits do not jeopardize the generality of the results provided in the next section and the assessment of the expected E1 Galileo SSV volume peculiarities.

4.2. SSV Galileo E1 Characterization Results

The first significant variable to be considered is the SSV available signal power.

The average distribution expressed by (22) is represented in Figure 14. The effect of the GRAP pattern variation for high off-boresight angles (side-lobes) can be clearly indicated by correlating the spanned off-boresight angles in Figure 12 with the corresponding signal power in Figure 14. It is also possible to see that as the altitude increases, the system service converges into three main regions. They are identified by different levels of power, approximately ${\mathrm{C}}_{SSV,ave}^R$ > 185 dBW, −205 dBW ${<\mathrm{C}}_{SSV,ave}^R$ < −185 dBW, ${\mathrm{C}}_{SSV,ave}^R$ < −205. According to (25), similar behavior on the user side is expected by computing ${\mathrm{C}/\mathrm{N}0}_{SSV,ave}^R\left({\theta }_i,h_{rx}\right)$. Such a distribution is provided in Figure 15, where the asymptotic separation of the three regions is even sharper.

The distributions provided can be easily expressed in terms of accessibility, which aggregates the information in a more compact index $A_{\%}$. Considering Equations (23) and (26) for different thresholds of the ${\mathrm{C}}_{SSV,ave}^R\left({\theta }_i,h_{rx}\right)$ and ${\mathrm{C}/\mathrm{N}0}_{SSV,ave}^R\left({\theta }_i,h_{rx}\right)$, it is possible to define the average accessibility charts describing the expected improvements when moving from one threshold to the other.

Those charts are represented in Figure 16 and Figure 17, respectively, considering the interval ${\mathrm{C}}_{min}^R\in \left[-220,-155\right]$ with a step size of 5 dBW for the received signal power and an interval of ${\mathrm{C}/\mathrm{N}0}_{min}^R\in [0, 35]$ with a step size of 5 dB/Hz for the carrier-to-noise ratio. The charts should be interpreted with respect to system and user requirements as follows. The system can guarantee more than 20% average accessibility up to 200 thousand kilometers for a user whose receiver design can cope with a received power level of −195 dBW (cyan curve in Figure 16); considering a fixed GNSS assembly (receiver plus receiving antenna), and a sensitivity of 25 dB/Hz (cyan curve in Figure 17), the user can benefit from at least 10% accessibility from HEO to the Moon.

It can be noted that the accessibility index possesses a great capability for representing the expected non-linear SSV performance variation. The second region, identifying side lobes access, is very well represented since the same signal power or receiver sensitivity increment corresponds to larger variations of the percentage, hence a larger increment of the effective Galileo orbit arc that can be exploited for navigation. It should also be noted that the sensitivity ranges 10–25 dB/Hz, accessing the second region at very high altitudes, is in line with current receiver solutions proposed for navigation in SSV up to lunar orbits [1].

Accessibility also has the advantage of aggregating the information in a scalar index, so it allows better visualization of E1 GRAP driven SSV with respect to azimuth variation. In fact, it is possible to compute (23) and (26) as functions of the azimuth through (21) and (25), hence providing azimuth-dependent accessibility distributions.

For the sake of brevity, only ${\left.A_{\%}\right|}_{C/{\mathit{N}}_0}$ variation with respect to azimuth is provided in Figure 18.

As the sensitivity increases, the accessibility increases, but more importantly several azimuth intervals become active in the sense of accessing a longer Galileo satellite arc. In the azimuth-dependent accessibility plots (Figure 18), sixteen beacons can be clearly recognized, including eight which are much stronger than the rest. For more than 20 dB/Hz, such beacons support satisfactory accessibility up to lunar altitude. The azimuth-dependent results complement the average accessibility, revealing more details about azimuth intervals where service can outperform the average value. The two metrics can be used for optimization of future antennas, i.e., maximizing accessibility in certain pattern regions with the constraint of achieving a target average.

5. Conclusions

This paper has introduced a novel methodology for the definition of GNSS constellation antenna patterns with the objective of shortening and simplifying the information transfer from the system to the final SSV user. The proposed procedure has been successfully implemented for the processing of on-ground isolated antenna measurements of Galileo FOC antenna patterns, leading to the definition of a Galileo Reference Antenna Pattern model, recently published in [8]. The solution introduced aims not only at regulating how the system and user share the pattern information, but also setting the basis for a standard procedure that can potentially be adopted by any GNSS system. The approach directly provides a 3D model at user level, which filters out discontinuities and separates what is best considered uncertain due to test-specific realizations and measurement degradation effects. Providing a statistical distribution in the form of confidence bounds, the system informs the user community about the regions or the specific points where performance can deviate from nominal values. The mathematical details of the antenna pattern reconstruction procedure have been presented herein to better understand the GRAP model properties with respect to the objectives of improved resolution, smoothness, and more general representation of the expected constellation EIRP. The approach and its statistical framework could also be considered a fundamental step towards updating GNSS antenna data processing to the state of the art of antenna characterization techniques and statistical learning. The statistical model developed can be considered as a-priori knowledge that can be exploited in the future for new data monitoring, prediction, and integration.

The additional results based on simplified geometrical analysis demonstrate how the GRAP model can be exploited for final Galileo SSV characterization. A new metric has also been introduced, accessibility, which can potentially also be adopted for other GNSS systems to better assess their SSV peculiarities. This index expresses the fraction of the rising GNSS SV orbit arc that the user can potentially exploit to access an SV signal for a target minimum received power. The approach does not pretend to replace standard assessments, which directly derive the GNSS SVs availability and navigation performance alongside the simulated space user trajectory. Accessibility has been proposed as a more compact and faster tool allowing the user to get general information about service quality for a target mission altitude and system or user requirements. In fact, the analysis presented for the Galileo E1 identifies different Galileo SSV operative regions, which reflect different ways to span the GRAP antenna patterns at different altitudes. The transition from one region to another is mainly driven by the capability of the user to access antenna pattern side-lobes for the target mission altitude, hence accessibility can be considered as an effective tool for the preliminary design of the on-board spaceborne GNSS system. It is also not excluded that, even limited by its assumptions, this compact and scalar index can be used as a key performance indicator for the design or adjustment of next-generation GNSS constellation antenna patterns. Small modifications in antenna pattern can correspond to relevant changes in accessibility. As such, this could represent a valid starting point for an optimization-based approach able to leverage SSV criticalities.

Future work will fit within the ongoing international effort for the definition of an Interoperable Global Navigation Satellite Systems Space Service Volume [1]. In the next steps, the GRAP model will be exploited to perform additional Galileo and multi-constellation SSV analysis with the possibility of fostering and better designing next-generation missions in high orbit.