A Survey of Rain Fade Models for Earth–Space Telecommunication Links—Taxonomy, Methods, and Comparative Study

Abstract

Satellite communication is a promising transmission technique to implement 5G and beyond networks. Attenuation due to rain begins at a frequency of 10 GHz in temperate regions. However, some research indicates that such attenuation effects start from 5–7 GHz, especially in tropical regions. Therefore, modeling rain attenuation is significant for propagating electromagnetic waves to achieve the required quality of service. In this survey, different slant link rain attenuation prediction models have been examined, classified, and analyzed, and various features like improvements, drawbacks, and particular aspects of these models have been tabulated. This survey provides various techniques for obtaining input data sets, including rain height, efficient trajectory length measurement techniques, and rainfall rate conversion procedures. No survey of the Earth–space link models for rain attenuation is available to the best of our knowledge. In this study, 23 rain attenuation models have been investigated. For easy readability and conciseness, the details of each model have not been included. The comparative analysis will assist in propagation modeling and planning the link budget of slant links.

1. Introduction

Satellite communication is playing major role in the back-haul data network. In addition, satellite communication links can be used to enhance the existing telecommunication infrastructure in 5G and beyond networks [1,2,3], and satellite communication also has enormous applications for unmanned aerial vehicles (UAVs) and the Internet of Things (IoT) [4]. Because of radio congestion and broader bandwidth requirements, satellites switch to higher-frequency (33–75 GHz and 75–110 GHz) bands. However, atmospheric disturbances wreak havoc on these frequencies, causing attenuation, scintillation, and depolarization, lowering service efficiency [5]. Rain is one of the significant factors that creates attenuation on the propagation of electromagnetic waves. This effect has directed the interest of researchers in reducing the effect of rain on radio waves by controlling the transmitted power. Thus, multiple studies have been conducted on this globally. The research on rain attenuation is used to predict rain attenuation in various geographical areas over a wide range of frequency bands, especially for frequency bands over 10 GHz, and determine a suitable model that can predict attenuation. To develop such a model, the factors that impact attenuation first need to be determined. Some of the factors are related to the infrastructure setup (latitude, longitude, antenna size, antenna height, elevation angle, operating frequency) and some of them are related to the rain events such as wind flow, humidity, wind direction, and temperature.

It may be possible to evaluate the attenuation of atmospheric elements such as wind flow and the effect of the direction of wind flow on propagated signal attenuation separately, but it is conventional to consider attenuation due to other factors such as temperature, wind flow, and wind direction, along with rain attenuation. Even though rain attenuation issues were noticed in the middle of the last century, the problem has not yet been efficiently addressed. Besides the necessity to transfer a massive volume of data due to the fourth industrialization, the lower frequency bands are getting exhausted. This has led researchers to develop rain attenuation models that can work with higher frequency bands such as Q/V, W, and E-bands. For example, the application of a rain attenuation model for slant link, where the operating frequency is at 75 GHz, was assessed in [6,7]. Several models of rain attenuation have been proposed in the literature, and scientists are aiming to refine existing models for local climate conditions [8,9,10,11,12,13].

Precise rain attenuation determinations are critical for preparing the link budget, maintaining communication efficiency, and the architecture of the system. In addition, following probable rain attenuation determined by the rain attenuation model, the over- or underestimation of the power of the transmitter can be avoided. Radio-frequency engineers must comply with the permissible power transmission specifications in compliance with spectrum management rules for each frequency band. If the engineers do not follow such power regulations, the transmitted signal power may interfere with another frequency band, e.g., neighbor channel, which may disturb the neighboring telecommunication equipment. Thus, by deploying a rain attenuation model in Earth–space telecommunication, such disturbance can also be avoided. In the literature, there is currently a lack of survey papers addressing the issues and algorithms of slant link fade prediction models, which inspired us to write this survey paper. In Table 1, some of the existing survey papers are presented, which are not adequate for covering the slant link fade prediction models.

Previously, we surveyed the recent literature in relation to terrestrial link rain attenuation models [17], which is also clearly different compared to the current study, as our previous terrestrial link survey was on Earth-to-Earth links only. Rain affects the propagation of electromagnetic waves both on the Earth–space and terrestrial links. The rain attenuation model for the Earth–space link is used to manage the rain attenuation on the Earth–satellite telecommunication link management. In comparison, the rain attenuation model is used to manage the rain attenuation on the Earth-to-Earth links. Thus, these two types of link are different although they have a resemblance. Even ITU-R proposed a separate model for rain attenuation for the Earth–space (ITU-R P.618 series) and terrestrial links (ITU-R P.530 series). Further differences between the Earth–space and terrestrial links are given in Table 2.

This paper presents findings such as classification of well-known rain attenuation models, supporting the frequency range, applicability to Earth–satellite links, and improvements or limitations of the models.

This study extensively investigates the well-known and essential features of up-to-date models, drawbacks, and unique features. We have also qualitatively compared the slant link rain attenuation models. The significant contributions of this study are as follows:

• To the best of our knowledge, slant link rain fade models have not been classified broadly. In this work, we have prepared a taxonomy of the rain fade models for Earth–space links.

• A brief overview of each of the models selected is presented. In addition, we have provided algorithms of different models.

• Quantitative and qualitative characteristics of different models are organized in tables to review comparative studies.

• We noticed that each model was improved inherently and criticized the prototype by finding the inconveniences, and the specific characteristics are listed.

• Finally, the open research issues are summarized.

The remainder of this paper is organized as follows: Section 2 describes related parameters of Earth–space rain attenuation models. In Section 3, a taxonomy of existing models is presented. In Section 4, core insights of the models, algorithms, advantages, and limitations are discussed. In Section 5, comparative studies among the models are presented. In Section 6, open issues and research challenges are summarized and discussed. Finally, concluding remarks are presented in Section 7.

2. Background Study

The most significant impacts on the propagation of slant microwave and millimetric wave links are hydrometeors such as thunder, hail, rain, ice, and snow. Among these parameters, rainfall attenuation is the most important phenomenon. The raindrops act as dielectric particle loss, which causes dispersion and depolarization.

2.1. Rain Attenuation Parameters

Developing a rain fade model involves mathematical analysis of rain attenuation phenomena by reasoning and cause-based interaction. A broad range of variables that can impact the rain attenuation is given in [15,18,19,20,21,22,23,24,25,26]. The parameters that affect the rain attenuation are rainfall rate, frequency, path length in the precipitation, temperature, humidity, wind velocity, wind direction, visibility, polarization, raindrop radius, raindrop size, latitude, elevation angle, station height, raindrop cross-section, depolarization loss, and refractive index of the air.

2.2. Earth–Space Link Budget

For rain attenuation studies and measurement, link budgeting is concerned with the losses and gains at the receiver antenna. The link budget does not differ much from terrestrial links except that free space loss (FSL) is much higher due to the considerable distance between the transponder and the Earth station [23,24,25,27].

(1)$P_r=P_t+G_{tx}+G_{rx}-FSL-L$

where $P_t$ is the transmitter power. This is the microwave carrier output power, expressed in $dBm$ in Equation (1); L is losses due to the presence of atmospheric gases, vegetation, buildings, clouds, and fog; $G_{tx}$ is the transmitter antenna gain, $G_{rx}$ is the receiver antenna gain, FSL is the free space path loss in dB, and $P_r$ is the received signal level in dBm. Compared to terrestrial links discussed in [17], FSL is higher, e.g., above 200 dB in Ka/Q bands [25]. The performance of a communication system is estimated based on the achievable signal-to-noise ratio (SNR) at the receiver. The term SNR (in dB) refers to the estimation of signal strength as a function of signal degradation and background noise. This SNR can be expressed as:

(2)$SNR=\frac{P_t{G}_{tx}{G}_{rx}}{kTB{L}_{sys}}{\left(\frac{4\pi l}{{\lambda }_0}\right)}^2$

where $L_{sys}$ is system loss at the receiver and transmitter; ${\lambda }_0$ is free-space wavelength (m), which is determined from the frequency as $\left({\lambda }_0=\frac{0.3}{f\left(\mathrm{GHz}\right)}\right)$; $G_{tx}$ and $G_{rx}$ are transmitting and receiving antenna gains, respectively; k is the Boltzmann constant $1.38\times {10}^{-23}$ J/K; B is the bandwidth in (Hz); l is the link distance (m); and T is the noise temperature (K) of the system, which is assumed to be 290 K.

The receiver antenna noise figure can be computed as [28]:

(3)$N\left(dBm\right)=10{log}_{10}\left(kTB\right)$

(4)$N\left(dBm\right)=174(dBm/Hz)-10{log}_{10}B\left(Hz\right)$

Equation (2) can equivalently be computed in unit decibels (dBm) as

(5)$\begin{array}{c}\hfill SNR\left(dB\right)=P_t\left(dBm\right)+G_{tx}\left(dB\right)+G_{rx}\left(dB\right)-FSL\left(dB\right)-L_{sys}\left(dB\right)\\ \hfill +174(dBm/Hz)-10{log}_{10}B\left(Hz\right)\end{array}$

(6)$SNR\left(dB\right)=P_r\left(dBm\right)-10{log}_{10}B\left(Hz\right)+174(dBm/Hz)$

where $P_r$ is determined from the link profile given parameters. In real-world radio links, for the communications system to acquire the least tolerable quality of service (usually termed as minimum SNR), the received power level often needs to be better than the threshold level. When the received signal magnitude is below the threshold value, a network outage will occur. As a result, the radio link performance evaluation is done by comparing the percentage of outage time with the total time duration. This design allowance for received power is termed as the link margin and defined as the difference between the design value of received power and the minimum threshold value of received power [25,26]:

(7)$LM=P_r\left(dBm\right)-P_{r,min}\left(dBm\right)$

where LM denotes the link margin measured in $dB$. LM, which is used to account for fading effects, is sometimes referred to as fade margin.

Clear Air SNR calculation: The received signal for clear-sky conditions is determined using Equation (1) and the SNR of this link in a clear-sky environment can be determined using (6).

Rainy time SNR calculation: Here the received signal is additionally degraded by the total rain attenuation amount. Thus, SNR at a rainy time is calculated as [23]:

(8)$SNR=(P_r-A_{rain})-10{log}_{10}B\left(Hz\right)+174\left(dBm\right)$

where $A_{\mathrm{rain}}$ is the attenuation induced by rain. If the SNR at a rainfall rate of exceedance of $0.01\%$ is required, it can be calculated by replacing $A_{\mathrm{rain}}$ by total rain attenuation at a rainfall rate of exceedance of $A_{0.01}$.

2.3. Rain Attenuation Anomalies: Breakpoint

At a breakpoint, the monotonous behavior of the time exceedance of rain changes abruptly in tropical regions. The breakpoint is more visible if the integration is considered over a smaller time duration; however, according to [29], a 1 min rainfall rate integration is the optimum time to clarify the rain structure to avoid noise. The main reason for modifying the model is considering the non-monotonous behavior of rain attenuation in tropical areas. It is found that there is a breakpoint of rainfall rate attenuation while moving from the temperate region to the tropical area, where the integration time is comparatively short. For such existing breakpoints, many of the temperate region models do not fit well in the tropical regions, and based on this, breakpoint idea-based models have been proposed in [30,31].

2.4. Slant Path in Rain Cell

The path length that remains in the rain cell structure suffers rain attenuation. There may be some other meteorological effects on the extended portion of the slant link, but there exists no rain attenuation effect. As a result, the slant path length that falls within the rain cell must be determined. However, since the length exists in the space, it is impossible to determine precisely, rather than considering the estimated length. The rain height estimation helps to estimate the slant path length in the precipitation.

2.4.1. Rain Height

In many rain attenuation models [32,33,34], the rain height is used to determine the slant link length. However, the rain height is not uniform in every location of the world, so it is crucial to estimate the rain height in a specific geographical area. For a low rainfall rate, the minimum rain height is about 4 km, although there are variations in different models. However, some of the models estimate rain heights with elevation angles and rainfall rate as well. Table 3 lists most of the estimation techniques for rain height determination.

2.4.2. Effective Slant Path

The effective path length is the length of a conceptual path derived from radio data by dividing overall attenuation by actual attenuation exceeded over the same percentage of the time. The total attenuation is determined as [35]:

(9)$A\left(dB\right)=\gamma (dB/\mathrm{km})\times L_{eff}\left(\mathrm{km}\right)$

where $A\left(dB\right)$ is the total attenuation due to rain, and specific attenuation of the link is $\gamma$, then the “effective path length” is defined by $L_{eff}$. Table 4 shows a list of additional techniques to determine ‘effective path length.

2.5. Rainfall Rate Conversion

The rainfall rate is the main parameter to determine the rain attenuation. The rainfall rate is measured in mm per time unit. Depending on the sensor device used to record the rainfall rate, the integration time unit may be 10 s, 20 s, to 1 min, 5 min, or even 1 h. The rain attenuation model may need to convert the rainfall rate over different integration times. The International Telecommunication Union-Radiocommunication sector (ITU-R) P.618 [45] recommends applying a 1 min integration time. If local rain data are available in other integration times, the rainfall rate needs to be converted. Table 5 lists most of the rainfall rate conversion techniques, as per the literature. Instead of a mathematical conversion procedure, the rainfall rate “contour map” can be used to obtain the rainfall rate [40,46,47,48].

Even if the rainfall rate is not readily available, or there are missing data in the collected records, there also exist synthetic techniques to generate rain attenuation time series. Mathematically, Equation (10) defines the rainfall rate in a typical rainfall rate acquisition system:

(10)$R_p=\frac{R_{tip}\times N}{T}$

where $R_p$ is rainfall rate, $R_{tip}$ is rainfall per tip (mm), N is number of tips, and T is the time gap between consecutive tips.

2.6. Rain Type

In total, two types of rain, known as stratiform and convective precipitation, exist. The convective rain is heavy; thus, the propagated radio waves attenuated by convective rain are heavily attenuated. Conversely, stratiform rain is comparatively light and comparatively less attenuated. Therefore considering the effects of different rain types can yield better-predicted attenuation. For example, the rain attenuation model in [62] considers stratiform and convective rainfall to predict rain attenuation. The convective type is heavy rain that corresponds to a high rainfall rate and stratiform rain corresponds to light rain and hence low rainfall rate. The structure of the rain cell is also incorporated with the rainfall rate, which is shown in Figure 1.

2.7. Rain Cell Size

In [43], the horizontal length $a\left(R\right)$ and the slant length $b\left(R\right)$ ($R<31.6$ mm/h case only) estimation technique is proposed as:

(11)$\begin{array}{c}\hfill a\left(R\right)=\frac{1}{0.056+0.012R}\end{array}$

(12)$\begin{array}{c}\hfill b\left(R\right)=-0.183-0.088R\end{array}$

where R is the rainfall rate and dimension parameters $a\left(R\right)$ and $b\left(R\right)$ are shown in Figure 1.

2.8. Rainfall Rate Missing Data

2.8.1. Temporal Missing: Time Series (TS)

The rain attenuation prediction requires long-term local rainfall rate information. However, in all parts of the world, such rainfall rate datasets are not available. In addition, even though the rainfall rate or rain attenuation datasets are recorded for a specific experimental setup, some parts of the dataset might be missing. To address these issues, a rainfall rate or attenuation time series needs to be generated. Fortunately, there are techniques available such as the synthetic storm technique (SST) [63] and enhanced synthetic storm technique (ESST) [64], and 2-D SST [65]. Rain attenuation prediction models such as the global synthetic storm technique (GSST) model [66] and the ESST model [64] use such types of time series generation procedures. Through applying such a procedure, as stated earlier, the missing datasets of any local climate can be generated.

2.8.2. Temporal Missing: Generation by Synthetic Means

The synthetic way to generate the rainfall rate of a locality needs local statistical distribution of the point rainfall intensity that feeds to the model, known as the exponential cell (EXCELL) model [67]. This model offers two types of rain cell structure, called the monoaxial model and the biaxial model. In a monoaxial model, the cell is presumed to be circular, and in a biaxial model, the model is assumed to be elliptical, and the monoaxial or biaxial property is determined by the orientation of central moments of inertia. In the EXCELL model, a typical rain cell can be defined by the following parameters: cell area (A), cumulative rain (Q), average rainfall rate ($\overline{R}$), root mean square (RMS) of rainfall rate ($R_{rms}^2$), rain cell gravity center ($x_0,y_0$), and the central moments of inertia ($I_x,I_y$). These parameters can be estimated from the local peak rainfall rate ($R_M$) and rain cell radius ($\rho$) in the ($x,y$) plane through which the rainfall rate decreases at the rate of $e^{-1}$ (note the exponential form). It is notable that the simple rain attenuation model (SAM) [32] and the LU model [68] inherently include distance as a negative exponential factor to determine the spatial rainfall rate. This hypothesis induces error if the rain cell area is smaller than 5 km² and the rainfall rate is smaller than 5 mm/h ($R_1$). To describe the rain cell model for the EXCELL model, the parameters that must be considered are as per Table 6:

A variant of the EXCELL model is the HYCELL model [69], where the rainfall rate is considered a result of an exponentially decreasing rainfall rate ($e^{-1}$), normally distributed (Gaussian process) rainfall rate, and as a finite linear combination of these two types of rainfall distribution. These two approaches of determining the rainfall rate are applicable separately depending on a specified rainfall rate ($R_2$) that acts as boundary value defining the application of the Gaussian process and the exponential rate to be used. Thus, if the rainfall rate satisfies $R_1\le R<R_2$, the rainfall rate distribution is generated through the exponential parameter, and if it is greater than $R_2$ the Gaussian process parameter is applied. The EXCELL and HYCELL models contribute to developing a single rain cell that may not cover a long link. Consequently, an extended version of the EXCELL model covering a long link was proposed in the literature, which is known as the MultiEXCELL model (a kind of propagation-oriented precipitation model based on a cellular rain formulation) [70]. In the MultiEXCELL model, instead of the biaxial elliptical approach, the monoaxial approach was considered. The minimum rain cell area was set 4 km² and the other distinct parameters were similar to the EXCELL model ($R_M$, ${\rho }_0$) and additional number $k=1,2\mathrm{or}3$. Compared with the EXCELL model, the MultiEXCELL model is additionally able to compute the synthetic rain cells’ probability of occurrence and to consider the convective and stratiform rain dominance. Additionally it considers the distribution along the radii conditioned to $R_M$ as $p\left({\rho }_0\right|R_M)$.

2.9. Effective Rainfall Rate

An inconsistency was found between the rainfall rate and link length reduction factor in [42]. As a result, in Equation (13), “path length reduction factor” was defined as:

(13)$r_p=\frac{A_p}{k{R}_p^{\alpha }·d}$

where $A_p$ is the rain attenuation, $R_p$ is the point rainfall rate exceeded at $p\%$ of the time, k and $\alpha$ are the specific attenuation coefficients, and d is the path length. Based on the developed idea, the effective rainfall rate concept was applied to the rain attenuation model. Proceeding with the same approach as per Equation (13), effective rainfall rate was derived as in Equation (14):

(14)$R_{eff}=R_{eff}\left(cos\theta +a_1·R^{a_2+a_3/L_{s}cos\theta }·L_S^{a_4}·sin\theta \right)$

where $a_1$, $a_2$, $a_3$, $a_4$ are constants that can be obtained through the curve fitting technique, $L_S=\left(H_R-H_S\right)/sin\theta$ is the slant path in the precipitation (km), $\theta$ is the elevation angle (see Figure 2). This Figure 2 presents a general Earth-space link scenario that will be used as the reference, from now and on-wards. With the SST technique, rainfall time series can be converted into rain attenuation time series, which is a useful tool in case of the absence of beacon measured attenuation data regardless of polarization and range of frequency bands but restricted to within elevations above 10${}^{\circ }$. It requires the advection velocity v (m/s) of rain, the length of the satellite path $\zeta$ (km), and the location of the ground station, which is indicated by $x_0$. In the next section, the meaning of these symbols will be that $L_G$ is the length of horizontal projection, $H_S$ is height above the mean sea level of the Earth station (km), $H_R$ is rain height, $H_R^{\prime }$ is the melting layer top limit, and $\phi$ is the latitude of the Earth station (degrees), if not mentioned otherwise.

2.10. Raindrop Size Distribution

Raindrop size distribution is used to generate synthetic rainfall rates using SST techniques [71,72] and the spatial distribution of rainfall rates. In [73], a raindrop size-distribution-based rain attenuation model was proposed. In the literature, different techniques are proposed to measure the raindrop size; two-dimensional video disdrometer [74], momentum disdrometer [75], and X-band polarimetric radar [76] are available, but there exist accuracy issues among these sensor devices as some of the devices were tested in [77]. Real measured accuracy comparison of rain distribution measured by the weather data and the DSD-based data was proposed in [78]. Raindrop size also affects the polarization property because propagating the radio wave through the rain creates depolarization of the propagated signal [79]. The drop size is one of the factors that affect the polarization. In the literature there are various types of raindrop size distribution, such as the normalized gamma model [80], seasonal variations in raindrop size [81], clouds and raindrop size distribution [82], vertical raindrop size distribution [83], exponential distribution [84], gamma distribution [85], log-normal distribution [86], and Weibull raindrop-size distribution [87].

3. Models to Predict Slant Link Attenuation

Apart from a formulation of the rain attenuation scheme, existing rain fade models for Earth–space links can be grouped into five categories. These include the empirical, physical, statistical, fade slope, and learning-based models.

• Empirical models: The model depends on experimental data findings rather than mathematically describable input–output relationships.

• Physical models: In the physical model, there exists some physical resemblance between the formulated rain attenuation model and the physical structure of rain.

• Statistical models: This type of model is built on the long-term data of rain attenuation, rainfall rate, and related atmospheric parameter statistical analysis.

• Fade slope models: In the fade slope model, a change in rain attenuation is determined from the fluctuations of measured experimental rain attenuation over time. These results can subsequently be used to forecast the attenuation of rain.

• Learning-based models: The learning-based rain attenuation is new in the domain of knowledge. Long-term rain attenuation and huge datasets of related parameters are used as the input to a learning network (i.e., artificial neural network) to train, and later this trained network (e.g., optimized weights) can be used to predict rain attenuation.

A complete taxonomy for the slant link rain attenuation models is presented in Figure 3. In addition, a typical illustration of an Earth–space link with related geometric parameters is shown in Figure 2.

4. Review: Slant Links Rain Fade Models

4.1. Empirical Model

4.1.1. SAM

The SAM model [32] has evolved through some previously published articles [101]. Here, the point rainfall rate at the ground is used. The model yields the attenuation in time series. However, it is not validated with heavy rainfall or with rainfall data from various locations worldwide. The quick steps of the SAM model are the following set of equations:

• (1). Calculate specific attenuation:

  (15)${\gamma }_{\mathrm{decay}}\left(R\right)=k{R}^{\alpha }(\mathrm{dB}/\mathrm{km})$

If the rainfall rate is unavailable, use the International Radio Consultative Committee (CCIR) recommended rainfall rate (former CCIR is now merged with ITU-R).

• (2). The rainfall rate varies along with the link, and the rainfall rate can be defined as:

  (16)$R\left(l\right)=\left\{\begin{array}{cc}{R}_p& R_p\le 10\mathrm{mm}/\mathrm{h}\\ R_p{e}^{-{\gamma }_{\mathrm{decay}}ln(R_p/10)lcos\theta }& R_p\ge 10\mathrm{mm}/\mathrm{h}\end{array}\right.$

  where $R_p$ is the point rainfall intensity, ${\gamma }_{\mathrm{decay}}$ is the decay profile along the horizontal axis, and $l\le L_S$ is defined by Equation (17):

  (17)$L_S=\frac{H_R-H_S}{sin\theta }$

  where rain height, $H_R$, is approximated by Equation (18):

  (18)$H_R=\left\{\begin{array}{cc}4.8\hfill & \mathrm{if}\left|\phi \right|\le {30}^{\circ }{\mathrm{R}}_p\le 10\mathrm{mm}/\mathrm{h}\hfill \\ 7.8-0.1\left|\phi \right|\hfill & \mathrm{if}\left|\phi \right|\ge {30}^{\circ }{\mathrm{R}}_p\ge 10\mathrm{mm}/\mathrm{h}\hfill \end{array}\right.$

• (3). Finally the rain attenuation is determined by Equation (19):

  (19)$A\left(R_0\right)=\left\{\begin{array}{cc}{R}_p^{\alpha }{L}_S\hfill & \mathrm{if}{R}_p\le 10\mathrm{mm}/\mathrm{h}\hfill \\ \frac{\left(k{R}_p^{\alpha }\right)(1-exp\left[-\alpha {\gamma }_{\mathrm{decay}}ln\left(R_p/10\right)L_{S}cos\theta \right])}{\alpha {\gamma }_{\mathrm{decay}}ln\left(R_p/10\right)cos\theta }\hfill & \mathrm{if}{R}_p\ge 10\mathrm{mm}/\mathrm{h}\hfill \end{array}\right.$

Algorithm 1 shows the procedure for this model.

Algorithm 1: SAM Model [32]

Advantages

It proposes a useful rainfall rate model that can yield the spatially distributed rainfall rate. The spatial distribution of rainfall concept later played a significant role in determining rainfall from microwave link attenuation [102] and real-time frequency scaling applications [103]. It extensively used slant-path attenuation, which combines the individual features of convective and stratiform types of rainfall.

Limitations

The usable frequency band is limited to 10–35 GHz.

4.1.2. Singapore Model

This model [33] uses the cumulative rainfall rate distribution modeled using the power-law relationship of 0.01% average year rain attenuation. The method of this model is shown in Algorithm 2.

Advantages

This model introduces the “path reduction factor” by modifying the CCIR model, analyzing the tropical climatic region in Singapore, which corrects the attenuation in the tropical region.

Limitations

The model underestimates the computed values down to outage times of 0.001% of a mean year. At a low elevation angle, the incoming signal beam intercepts more than one rain cell [31]. The model is suitable only within the frequency range of 4–6 GHz, and which is validated by the authors; more analysis is required for higher frequencies. However, the model does not consider rainfall distributed over the cell in a non-uniform pattern, as in a real case, non-uniformly distributed rainfall does not normally exist.

Algorithm 2: Singapore Model [33]

1. if$0\le \phi <{36}^{\circ }$then   2.   3. else   4.   5. if$\theta \ge 5^{\circ }$then   6.   7. else   8.   9. Compute $L_G⟵L_S·cos\theta$/*$\phantom{(}$$L_G$ in km $\phantom{(}$ */ 10. Compute ${RF}_{0.01}\leftarrow 1/(1+{L_G}^{0.9}/35exp(-0.015R_{0.01}))$ /*$\phantom{(}$ RF: reduction factor $\phantom{(}$ */ 11. Compute (k, $\alpha$) /*$\phantom{(}$ CCIR [104]$\phantom{(}$ */ 12. Compute $A_{0.01}\leftarrow k{\left(R_{0.01}\right)}^{\alpha }{L}_S·{RF}_{0.01}$ 13. Compute $\gamma \leftarrow k{R^{\alpha }}_{0.01}$ 14. Compute $A_p\leftarrow A_{0.01}0.12p^{-(0.632+0.082log(p\left)\right)}$ 15. Return $A_p$

4.1.3. Garcia Lopez Model

The Garcia Lopez Model [34] is based on measurement data from satellite links in Europe, Japan, Australia, and the USA. The model has four coefficient constants, $c_1,c_2,c_3,c_4$, dependent on the geographical area to support the regional climate. However, the physical significance of these coefficients is unknown, and these coefficients do not have any physical meaning. The model is suitable for predictions with elevation angles ranging from 10${}^{\circ }$ to 40${}^{\circ }$. The detailed steps of the model are given in Algorithm 3.

Algorithm 3: Garcia Lopez Model [34]

1. Calculate k and $\alpha$ /*$\phantom{(}$ CCIR [104]$\phantom{(}$ */   2. Calculate specific attenuation $\gamma \leftarrow k{R^{\alpha }}_{0.01}$   3. if $\phi <{36}^{\circ }$ then   4.   5. else   6.   7. if $\theta \ge 5^{\circ }$ then   8.   9. else 10. 11. Calculate, $A_p\leftarrow \frac{\gamma L_S}{[c_1+\left\{L_S(c_{2}R+c_3{L}_S+c_4)\right\}/e]}$/*$\phantom{(}$$c_1,c_2,c_3,c_4$ are constants $\phantom{(}$ */ 12. Return $A_p$

Advantages

The model can predict rain attenuation at the receiving ground station end by finding four coefficients (a, b, c, and d), which makes it possible to apply the model anywhere in the world.

Limitations

The model suffers from a problem in satellite links having a low elevation angle, with low rain intensities.

4.1.4. GSST Model

The GSST model [66] is based on the SST technique and produces all the necessary rain attenuation statistics at any frequency, elevation angle, or polarization. The derivation idea of this standalone attenuation formula is given in [105]. The algorithmic rule of this model is shown in Algorithm 4, where the significance of related parameters can be checked in Figure 2.

Advantages

The model can faithfully convert the rainfall rate into rain attenuation only from the probability distribution of rain. It can also produce rain attenuation time series at an arbitrary polarization and frequency for every slant path over an approximately 10${}^{\circ }$ elevation angle.

Limitations

It considers “isotropy” of rainfall, which does not consider the “denser inner side” and comparatively “lighter outer side” of the rain structure. This model’s main drawback is that the SST technique needs to consider the velocity of wind, which is not known accurately, and the model has not been directly verified outside the temperate region.

4.1.5. LU Model

The LU Model [68] is based on the correction of the rainfall rate modification factor’s distribution and an exponential rain cell structure. The comparative results show that it offers the lowest RMS and standard deviation (STD) over the full range of the percentage of the time.

Algorithm 4: GSST Model [66]

• (1). The distribution of rainfall rate projection on the Earth is defined as:

  (20)$R(x,p)=R_{max}\left(p\right)e^{-b\left|x-L_G\right|}$

  where $b=1/\rho$ here $\rho$ (km) is the same radius for which the rainfall rate is fixed $R_{max}/e$.

• (2). The attenuation along with the slant link:

  (21)$A\left(p\right)=k{\left[R\left(p\right)r\left(p\right)\right]}^{\alpha }{L}_S$

  where $r\left(p\right)$ is the rain rate adjustment factor, and $R\left(p\right)$ is the rain rate exceeded for $p\%$ of an average year.

• (3). The coefficients $a_i(i=1,2,\dots ,6)$ are taken using the genetic algorithm (GA) and annealing algorithm; thus, exploiting the rain databank of the ITU-R for the slant link:

  (22)$r\left(p\right)=3.78(1-\frac{0.85p^{0.065}}{1+0.12L_S})R{\left(p\right)}^{-0.56+\frac{1.51}{L_S}}$

The procedure to determine attenuation by this model is shown in Algorithm 5.

Algorithm 5: LU Model [68]

//$\phantom{(}$ Assumption: abrupt change more than 1 dB is neglected $\phantom{(}$ 1. Compute $L_S$ 2. Calculate $L_G$ 3. Calculate $R(x,p)$ /*$\phantom{(}$ Equation (20)$\phantom{(}$ */ 4. Determine $r\left(p\right)$ /*$\phantom{(}$ Equation (22)$\phantom{(}$ */ 5. Compute $A\left(p\right)\leftarrow k{\left[R\left(p\right)r\left(p\right)\right]}^{\alpha }{L}_S$ 6. Return $A_p$

Advantages

This model introduces the concept of “rainfall rate correction factor”, and it shows monotonic attenuation behaviors with an elevation angle compared to other rain attenuation models. The unusual behavior of rain attenuation (such as a nonmonotonic change in attenuation with respect to elevation angle and percentage of time or rainfall rate) was improved.

Limitations

The monotonically changing behaviors of rain attenuation are a valuable property of the LU Model, but applicability is limited by low latitudes and low elevation angles.

4.1.6. Karasawa Model

In the Karasawa model [88], the rain area size parameter acts as a function of the rainfall rate for 0.01% of the time to minimize the prediction error. The model considers a “single volume” for falling rain and does not account for the two types of rainfall, convective and stratified rain. The model’s drawbacks are that it applies only to mid-latitude regions, at relatively high elevation angles, and supports frequency bands between 10–20 GHz. The detailed steps of the Karasawa model are given in Algorithm 6.

Algorithm 6: Karasawa Model [88]

1. Calculate k and $\alpha$   2. Calculate specific attenuation $\gamma \leftarrow k{R^{\alpha }}_{0.01}$   3. if $\phi <{36}^{\circ }$ then   4.   5. else   6.   7. if $\theta <{10}^{\circ }$ then   8. /*$\phantom{(}$ CCIR 564-2 [106]$\phantom{(}$ */   9. else 10. 11. Calculate, $L_0$←$aexp(-b{R}_{0.01})$ /*$\phantom{(}$ a, b are constants.$\phantom{(}$ */ 12. Calculate, $A_{0.01}\leftarrow \gamma ·L_S$ 13. Calculate for p%, $A_p\leftarrow 0.12A_{0.01}{p}^{-(0.546+0.043logp)}$ /*$\phantom{(}$${}_{0.001\le p\le 1}$$\phantom{(}$ */ 14. Return $A_p$

Advantages

The model uses vertical and horizontal path length “reduction factors.” It also includes a rain cell size parameter that depends on the rainfall rate (0.01%). This parameter can reduce the prediction error (e.g., eliminate the weakness of the former CCIR methods). Another beneficial aspect of this model is that, except for rainfall rate, the model considers the “probabilities of occurrence” and “mean rainfalls” for both cell and debris rain.

Limitations

This model considers uniform rainfall intensity (mm/h); in real life, rainfall is not normally uniform in a rain cell. Another drawback is that the model considers a single volume for falling rain and does not account for the difference between the cell and debris rain.

4.1.7. Breakpoint Model

This model [29] is a modification of the ITU-R model [107] and considers “breakpoints” in the rainfall rates in tropical and sub-tropical climates.

The main steps of determining attenuation with the Breakpoint model are as follows:

• (1). Calculate the “percentage of exceedance” at the breakpoint:

  (23)$p_{R_B}of0.01\%=p_{A_B}of0.021\%$

• (2). Compute:

  (24)$A_p=\left\{\begin{array}{cc}{A}_B{\left(\frac{p-0.011}{0.01}\right)}^{P_1},\hfill & 0.021\le p<1\\ A_B{\left(\frac{p}{0.021}\right)}^{P_2},\hfill & p<0.021\end{array}\right.$

  where $P_1$, and $P_2$ are given by Equations (25) and (26):

  (25)$\begin{array}{c}\hfill P_1=-(0.655+0.033ln(p-0.011)\\ \hfill -0.045ln\left(A_B\right)-\beta (0.989-p)sin\theta )\end{array}$

  (26)$\begin{array}{c}\hfill P_2=-0.5(0.655+0.033ln\left(p^2\right)\\ \hfill -0.045ln\left(A_B\right)-\beta (0.989-p)sin\theta )\end{array}$

The experimental results agreed with the ITU-R model [107] regarding correction factor-induced attenuation, especially at elevation angles less than 60${}^{\circ }$. Algorithm 7 shows the procedure for this model.

Algorithm 7: Breakpoint Model [29]

Advantages

This model determined the “breakpoint” in the tropical region and introduced a correction factor to the rain attenuation model with an elevation angle range $40{}^{\circ }\le \theta <60{}^{\circ }$.

Limitations

The correction factor that is proposed in this model has not been analytically validated.

4.1.8. Unified Model

This [89] model uses a “complete rain distribution” based on several non-linear rain regressions using several attenuations and rainfall databanks. This model proposes an “effective path length” determination procedure for the terrestrial and slant links in a common way; hence it is called “unified.” It was first designed for the terrestrial links and later added to the slant communication links. The rain attenuation is given by Equation (27):

(27)$A_p=k{\left[R_{eff}(R_p,L_S,\theta )\right]}^{\alpha }\frac{L_S}{1+L_{S}cos\theta /L_0}$

where $L_S$ and $R_{eff}(R_p,L_S,\theta )$ are given by Equations (28) and (29):

(28)$L_S=(H_R-H_S)/sin\theta$

(29)$R_{eff}(R_p,L_S,\theta )=1.74R^{0.786+0.197/L_{S}cos\theta }\left(cos\theta +\frac{120}{{L_S}^{2.88}}{R}^{-0.186}sin\theta \right)$

The working procedure of this model is shown in Algorithm 8.

Algorithm 8: Unified Model [89]

1. Choose rainfall rate, $R_p$ 2. Compute $L_S$ /*$\phantom{(}$ Equation (28)$\phantom{(}$ */ 3. Compute ${\mathrm{R}}_{\mathrm{eff}}$ /*$\phantom{(}$ Equation (29)$\phantom{(}$ */ 4. Compute $A_p$ /*$\phantom{(}$ Equation (27)$\phantom{(}$ */ 5. Return $A_p$

Advantages

The presented results show that the model outperforms ITU-R P.618-12 [108].

Limitations

The model has not been verified with measured attenuation data.

4.1.9. YJ.X. Yeo, Y.H. Lee and J.T. Ong (YLO) Model

This model [12] is based on “extensive satellite link measurement data” and uses the ITU-R P.838-3 prediction model [109]. The main area of concern for its application is the tropical region. The model has been verified up to a frequency of 30 GHz. The rain attenuation determination procedure by this model is presented in Algorithm 9.

Algorithm 9: YLO Model [12]

Advantages

It has a path length adjustment factor to refine the rain attenuation.

Limitations

It used only experimental data, which is of a short duration compared to long-term term data. Thus there is scope to verify the model with a long-term measured dataset.

4.2. Statistical Models

4.2.1. Maseng–Bakken (M-K) Model

This model [90] is based on the log-normal distribution of rain and utilizes a 1-dimensional Gaussian Markov process (generated by a memoryless non-linear device for transforming attenuation and rain intensity, Figure 4). With the assumption that the attenuation rate is constant throughout the path, rainfall intensity is given by Equation (30):

(30)$\gamma =k{R}^{\alpha }L$

The symbols k and $\alpha$ denote radio wave propagation parameters that depend on the frequency and polarization, and L is the length of the link. It was examined that the $\gamma$ is log-normal distributed. Consequently, following the log-normal distribution, rainfall intensity can be written as: $\Lambda (R_m,{{\sigma }_R}^2;R)$. The justification of the modeling of the rainfall intensity distribution is given in Equation (31):

(31)$\Lambda ({({\gamma }_m/k·L)}^{1/\alpha }m{\sigma }_{\gamma }^2/{\alpha }^2;R)=\Lambda (R_m,{{\sigma }_R}^2;R)$

The conformity of mapping the distribution of rainfall intensity through $\gamma$ can be accessed through determining the auto-correlation between the variables, as in Equation (32):

(32)$R_{\gamma }\left(\tau \right)={\gamma }_m^{2}exp\left({\sigma }_{\alpha }^2(1+\rho \left(\tau \right))\right)$

After normalization and the removing the bias:

(33)$k_{\gamma }\left(\tau \right)=\frac{exp\left({\sigma }_{\gamma }^2\rho \left(\tau \right)\right)-1}{exp\left({\sigma }_{\gamma }^2\right)-1}$

where $\rho \left(\tau \right)=exp(-\beta |\tau \left|\right)$ and the value of $\beta$ can be determined as:

(34)$\beta =E\left\{{(\gamma (t+\Delta t)-\gamma \left(t\right))}^2/\gamma \left(t\right)\right\}/2{\sigma }_{\alpha }^2\Delta t$

where $\Delta t$ is the sampling time interval determined from the time series of $\gamma$. A constraint of the M-K model is that the supporting frequency band only goes up to 35 GHz. This model has not been validated with heavy rainfall or with rainfall data from various locations around the world. A modified version of the M-K model, which introduced an offset parameter called the “calibration factor” after the memoryless non-linear device, was reported in [110]. The sequential steps of the model are given in Algorithm 10.

Algorithm 10: M-K Model [90]

1. Get the “LPF” filtered output; 2. Pass this filtered “signal” through the memoryless non-linear device; 3. Compute Equation (33); 4. Return attenuation;

Advantages

The model is straightforward, which allows simulation of communication link functionalities.

Limitations

The model has not been verified with diverse weather conditions such as low to high rainfall or rain type.

4.2.2. Modified Genetic Algorithm-ARIMA (MARIMA) Model

This model [91] combines GA concepts with the autoregressive integrated moving average (ARIMA) filter. The detailed steps of the MARIMA model are given in Algorithm 11. The validated results of the MARIMA model show a good agreement between measured and predicted attenuation.

Algorithm 11: MARIMA model [91]

Advantages

The model can determine fade at millimeter-wave bands dynamically, and also, it is an eight-parameter-based model where each value of parameter depends on a scheme (based on inter-parameter dependency).

Limitations

The local developing procedure is not easy to fit with the model; consequently, the parameters of this model have not been developed worldwide.

4.2.3. ITU-R P.618 Model

In the ITU-R P.618 model [45] the rain attenuation calculation procedure is complicated as several factors need to be estimated before determining the rain attenuation. The factors that need to be calculated are effective path length, vertical adjustment factor, estimation of the rain height, and calculation of the slant path length and its horizontal projection. A detailed procedure of the ITU-R P.618 model is given in Algorithm 12.

Algorithm 12: ITU-R P.618 Model [45]

Advantages

ITU-R P.618 showed good performance in a tropical Malaysian experiment [8,111], Italy [112], and a site diversity study in Japan showed good performance [113].

Limitations

It is reported that ITU-R P.618 overestimates the CCDF of rain attenuation with version P.618-11 in France [114].

4.2.4. PEARP Model

The PEARP model [92] is based on the probability of exceeding the attenuation threshold level. The steps of the PEARP model are shown in Algorithm 13. Validated information for the PEARP model was not found. The probabilistic weather forecasts could be useful for transmitting data maintaining a quality link in the temperate region. Hence, it can help to maximize the economic value of transmitted data at a high frequency. However, such information for the tropical region is not available yet. Algorithm 13 shows the procedure for this model.

Algorithm 13: PEARP model [92]

1. Get the time series attenuation data 2. Pass the data through LPF 3. Find fade slope $\zeta$ 4. Calculate PDF ($\zeta /A$) /* PDF: probability density function */ 5. Extract (${\sigma }_{\zeta }$, $\zeta$) 6. Find rain attenuation ${\sigma }_{\zeta }\left(A\right)$

Advantages

Useful for setting several types of rainfall rate based on the distinct threshold level. The model allows getting the benefit of ensemble weather forecasts using satellite communications.

Limitations

This model has low horizontal and temporal resolutions.

4.2.5. Das Model

The concept of the Das model [93,115] is based on considering the Gaussian distribution of the conditional occurrence of rain attenuation with a specific value of attenuation occurring before. The validation demonstrates a good agreement between measured and predicted attenuation using this model. However, it has not been compared with other relevant models. The working procedure of this model is shown in Algorithm 14.

Algorithm 14: Das Model [93]

Advantages

The model showed good performance for predicting the time series (of rain attenuation) at Ku bands over the Earth–space link. Another exemplary aspect of the Das model is that it can also be applied if there are some missing measured data.

Limitations

The Gaussian distribution was used in this model for predicting the attenuation time series, and the distribution was not fitted well with measured data. Another drawback is that if more than a 1 dB change in input power is detected, the procedure simply changes the current value by ±1, which fails to track more than a 1 dB change in the previous input.

4.3. Physical Models

4.3.1. Bryant Model

The Bryant model [37] comprises rainfall rate, Earth station height above sea level, and elevation angle at the Earth station. Equation (40) determines the rain attenuation, whereas Equations (35)–(39) are used to determine the required parameter values that are used to determine rain attenuation.

(35)$\begin{array}{cc}\hfill H_R& =4.5+0.0005R^{1.65}\left[\mathrm{km}\right]\hfill \end{array}$

(36)$\begin{array}{cc}\hfill D_m& =170R^{-1.2}\left[\mathrm{km}\right]\hfill \end{array}$

(37)$\begin{array}{cc}\hfill k_n& =e^{0.007R}\hfill \end{array}$

(38)$\xi =\left\{\begin{array}{cc}(1/\sqrt{2})e^{sin\theta }& \theta \le {55}^{\circ }\\ 1.1tan\theta & \theta >{55}^{\circ }\end{array}\right.$

(39)$L_S=\left\{\begin{array}{cc}\frac{H_R-H_S}{sin\theta }\left[\mathrm{km}\right]& \theta \ge 5^{\circ }\\ \frac{2(H_R-H_S)}{\sqrt{{sin}^2\theta +2(H_R-H_S)/R_E}+sin\theta }\left[\mathrm{km}\right]& \theta <5^{\circ }\end{array}\right.$

(40)$A=1.57D_m{k}_n\gamma \frac{L_S}{\xi L_G+D},\mathrm{for}3\mathrm{mm}/\mathrm{h}$

where $R_e$ is the equivalent radius of the Earth, $L_G$ is the projection of the slant path in the horizontal direction, and D is the single-cell diameter. The model has been validated using the ITU-R databank. The rain attenuation determination procedure by this model is presented in Algorithm 15.

Algorithm 15: Bryant Model [37]

Advantages

The model proposes rainfall rate distribution in two classes based on a breakpoint (a fixed rainfall rate), where rainfall rates below and above this breakpoint is uniform and non-uniformly distributed.

Limitations

For available rainfall data, it uses full rainfall rate distribution (e.g., effective rain cell) and argues a fixed breakpoint rainfall rate that is almost independent of the local climate. However, precise conclusions have not been made.

4.3.2. Crane TC

This model [94] has been developed primarily for Western Europe and the USA, with difficulty determining precipitation parameters for the weak and robust precipitation cell, including the probability of occurrence and average precipitation. These large and weak rain cells are often referred to as “debris” and the “cell.” For both satellite and terrestrial connections, the model has been tested. The procedure to determine attenuation by this model is shown in Algorithm 16.

Algorithm 16: Crane TC [94]

Advantages

It contains in the signal path the non-uniform definition of a heavy and light rain zone. The model also contains the standard statistics necessary for designing the versatile diversity scheme, the interference statistics due to rainfall, and attenuation statistics due to the percentage of the time.

Limitations

The required statistical parameters are available for the USA and part of the Europe region only. In addition, this method requires solving many complex equations.

4.3.3. Physical-Mathematical (PM) Model

The PM model [95] is based on experimental 1 min time series data, and it is modeled as two layers of precipitation phenomena and the cell-based rain structure. It is based on the CCIR model, along with the application of the SST. The model was tested experimentally with a long-term frequency of 11.6 GHz and an SST velocity speed of 10.6 m/s at three Italian cities; it showed −10.6% absolute error, 7.6% STD, and 13% RMS error, with an outage probability range of 0.005%. The detailed procedure of the PM model is as follows:

(41)$\begin{array}{cc}\hfill A\left(x_0\right)& =k_A{\int }_0^{L_A}{R}^{{\alpha }_A}(x_0+\Delta x_0,\zeta )d\zeta +k_B{r}^{{\alpha }_B}{\int }_{L_A}^{L_B}{R}^{{\alpha }_B}(x_0,\zeta )d\zeta \hfill \end{array}$

An illustration of these parameters is given in Figure 2.

(42)$\begin{array}{cc}\hfill L_A& =\frac{H_A-H_S}{sin\theta }\hfill \end{array}$

(43)$\begin{array}{cc}\hfill L_B& =\frac{H_B-H_S}{sin\theta }\hfill \end{array}$

(44)$\begin{array}{cc}\hfill \Delta L& =L_B-L_A=h/sin\left(\theta \right)\hfill \end{array}$

Along the horizontal axis, the speed of the wind is given by:

(45)$\begin{array}{cc}\hfill v\left(\theta \right)& =v/cos\left(\theta \right)\hfill \end{array}$

Transforming Equation (45) into a space spectrum, and with further analysis, Equations (46) and (47), yield the attenuation for the elevation angles $\theta \ne {90}^{\circ }$ and $\theta ={90}^{\circ }$, respectively.

(46)$\begin{array}{cc}\hfill \underset{v\to \infty }{lim}A\left(t\right)& =k_A{R}^{{\alpha }_B}\left(t\right)L_A+r^{{\alpha }_B}{a}_B{R}^{b_B}\left(t\right)(L_{\mathbf{B}}-L_A),for\theta \ne {90}^{\circ }\hfill \end{array}$

(47)$\begin{array}{cc}\hfill A\left(t\right)& =k_A{R}^{{\alpha }_B}\left(t\right)\left(H_A-H_S\right)+r^{{\alpha }_B}{k}_B{R}^{{\alpha }_B}\left(t\right)h,for\theta ={90}^{\circ }\hfill \end{array}$

The result reveals that the SST technique demonstrates the smallest RMS and average error rate compared to [32,34,116,117,118,119,120,121]. The detailed steps of the PM model are given in Algorithm 17.

Algorithm 17: PM Model [95]

Advantages

It is the first model that comes with the practical melting and rain layer consideration and applying a 1 min rainfall rate (to distance). The method uses the SST, which can also be used to calculate fade duration and its changes.

Limitations

The model needs a 1 min rainfall rate and needs a storm translation speed, which may not be readily available worldwide.

4.3.4. SCEXCELL Model

The SCEXCELL model [62] is based on the rain attenuation statistics generated from the Earth–space links. This model uses the synthetic EXCELL cell to obtain a finite value using a lower factor and modifying the rain profile. Furthermore, the model identifies different types of precipitating clouds, such as cumulus- and stratus-type clouds. The lower value of the rainfall rate is used to separate the yearly rain probability into stratiform and convective distribution, where it is assumed that for a stratiform rainfall, a rate of 10 mm/his rarely exceeded. The procedure for assessing attenuation due to rain involves the derivation of the probability of stratiform and convective rainfall from the yearly rainfall rate based on the discrimination method (steps are shown in Algorithm 18) [122].

The convective and stratiform rain heights were calculated using the local climatological parameters ($h_i$: the monthly mean values of the 0${}^{\circ }$ isotherm height; ${\beta }_i$: the ratio between the convective and the total rain amounts calculated over each successive 6 h period, as estimated by the European Center for Medium-Range Weather Forecasts (ECMWF) forecast process for the ECMWF Re-Analysis-15 (ERA-15) databank; and $p_i$: monthly mean values of the 6 h rainy periods’ probability) from the ERA-15 databank and the average equivalent bright-band height (BBH, a layer where maximum radar reflectivity happens), as defined by Equation (48).

(48)$BBH\left(f\right)=4.58e^{-0.0675f}+0.51\left[\mathrm{km}\right]$

Equation (48) is valid for the following conditions: initial density of the melting particles $\left(d_0=0.16\mathrm{g}/\mathrm{c}{\mathrm{m}}^3\right)$; particle size distribution (gamma with $\mu =0$); and frequency $(10\mathrm{GHz}<f<50\mathrm{GHz})$. Due to the concept of non-ideal exponentially shaped cells, the model considers the plateau contribution to the attenuation due to stratiform and convective types of rain differently.

(49)$\begin{array}{cc}\hfill A{\left(P\right)}_{str}& =A_C{\left(P\right)}_{str}+A_P{\left(P\right)}_{str}\hfill \end{array}$

(50)$\begin{array}{cc}\hfill A{\left(P\right)}_{cnv}& =A_C{\left(P\right)}_{cnv}+A_P{\left(P\right)}_{cnv}\hfill \end{array}$

(51)$\begin{array}{cc}\hfill P{\left(A\right)}_{tot}& =P{\left(A\right)}_{str}+P{\left(A\right)}_{cnv}\hfill \end{array}$

Although the model does not outperform the existing ITU-R model (it has approximately 2.4% error), it does consider supporting site diversity [123], the conversion of rainfall rate from 1 min to 1 h [124], and interference due to hydrometer scattering [125]. The rain attenuation determination procedure by this model is shown in Algorithm 18.

Algorithm 18: SCEXCELL Model [62]

Advantages

The model fuses the benefits of the EXCELL and HYCELL models. It facilitates the simulation of long microwave terrestrial links (and, consequently, extended terrestrial microwave networks) to evaluate site diversity systems.

Limitations

For long-range rain databank development, e.g., 150 km coverage, clutter values may appear, which could create an error.

4.4. Fade Slope Model

Fade mitigation techniques require proper design and implementation by introducing the first- and second-order statistics such as fade slope, duration of fades, and duration between fades [29]. The time interval between two crossings is defined as the time interval between the same attenuated threshold, and the device length is defined as the time interval between two crossings beneath the same attenuating threshold [126]. Various estimation studies of fade slope properties in temperate climatic conditions are reported in [127,128]. Since rain characteristics differ significantly over tropical regions, widely used ITU-R models are insufficient for predicting rain attenuation and modeling related propagation impairments [129]. Similar to the rain attenuation model, different variants of rain fade slope models are proposed for tropical climatic conditions [130]. By definition, the fade slope can be defined in [130] as Equation (52). The same technique is reported in [131] with consideration of the $\Delta t$ to be 2 s. A typical fade slope diagram is shown in Figure 5.

(52)$\zeta \left(t\right)=\frac{A(t+0.5\Delta t)-A(t-0.5\Delta t)}{\Delta t}$

4.4.1. Japan Model

In the Japan model [96], signal-level variations on Ku-band low-elevation Earth–space paths are examined. The significant signal fades due to rain, and tropospheric scintillations, sometimes coincide [132]. If the effect of scintillation becomes dominant, it needs to filter out, especially for low-margin systems operating at low elevation angles. To obtain the low-frequency components (X_low), a moving average filter (MAF) was used (53) to filter out the high-frequency components (X_high) that results from scintillation effects, according to the Equation (54):

(53)$\begin{array}{c}\hfill X_{\mathrm{low}}\left(t\right)=\frac{1}{M}\sum _{j=-m}^{m}X(t+j\Delta t)\end{array}$

where M is equal to $2m+1$ and $\Delta$ is equal to 1 s.

(54)$X_{\mathrm{high}}\left(t\right)=X\left(t\right)-X_{\mathrm{low}}\left(t\right)$

The detailed steps of the Japan model are given in Algorithm 19.

Algorithm 19: Japan Model [96]

/* Choose proper sampling time to avoid aliasing effect on $X\left(t\right)$ */1  Find Xlow /* Apply Equation (53) */2  Find Xhigh /* Apply Equation (54) */3  Apply FFT /* Choose number of FFT points for better visualization */4  Generate time domain diagram /* To observe the fluctuations */

Advantages

It provides a mechanism to identify the attenuation due to rain and the scintillation effect. The results show that rain attenuation and signal fade could be treated independently.

Limitations

The technique has been tested only at low elevation angle Earth stations and for the Ku-band only.

4.4.2. Das Fade Model

The Das Fade model [133] was implemented with a modification of the Van de Kamp (VDK) model [127]. This work also blended the fade slope concept to predict the time series attenuation series instead of customarily using velocity [64,97,134,135] or the synthesizer-based technique [136]. Algorithm 20 shows the procedure for this model.

(55)$\begin{array}{cc}\hfill \xi \left(i\right)& =(A\left(i\right)-A(i-1))/T_s(\mathrm{dB}/\mathrm{s})\hfill \end{array}$

(56)$\begin{array}{cc}\hfill P\left(\frac{\xi }{A}\right)& =2/\left(\pi {\sigma }_{\xi }{\left(1+{\left(\xi /{\sigma }_{\xi }\right)}^2\right)}^2\left({(\mathrm{dB}/\mathrm{s})}^{-1}\right)\right)\hfill \end{array}$

(57)$\begin{array}{cc}\hfill {\sigma }_{\xi }\left(A\right)& =-0.00031{\mathrm{A}}^3+0.0038{\mathrm{A}}^2-0.0047\mathrm{A}+0.04\hfill \end{array}$

(58)$\begin{array}{cc}\hfill A\left(t+T_s\right)& =A\left(t\right)+\xi \left(t\right)x{T}_s\hfill \end{array}$

Algorithm 20: Das Fade slope model [133]

1  Apply LPF and MAF;2  Find the fade slope;3  Calculate PDF /* PDF: probability density function */4  Extract statistical coefficient (${\sigma }_{\zeta }$) /* ${\sigma }_{\zeta }$: standard deviation of $\zeta$ */5  Fit polynomial to fit attenuation with ${\sigma }_{\zeta }$6  Use $\zeta$ to calculate rain attenuation for time series data prediction ($A_p$)7  Return $A_p$

Advantages

The enhanced fade-slope model updates the slope from the modified VDK model and can predict the time series of rain attenuation of the Ku-band signal during rain events for tropical locations.

Limitations

According to the results, there exists a 10% error compared to the measured value for predicting attenuation by this model.

4.4.3. ITU-R P.1623

The ITU-R fade slope model [98] is valid for the frequency range 10–30 GHz and applicable with elevation angles of 10–${50}^{\circ }$. ITU-R P.1623 defines the rate of change of fade by Equation (59). The distribution of the fade slope depends on the attenuation signal strength, the smaller time interval, and the specifications of moving-average filters. To compute the fade slope, wepass the attenuation time series through a low-pass filter (LPF). To calculate $F(f_B,\Delta t)$:

(59)$F\left(f_B,\Delta t\right)=\sqrt{\frac{2{\pi }^2}{{\left(1/f_B^{2.3}+{(2\Delta t)}^{2.3}\right)}^{1/2.3}}}$

We then use $F(f_B,\Delta t)$ to determine ${\sigma }_{\zeta }$:

(60)${\sigma }_{\zeta }=s·F\left(f_B,\Delta t\right)·A\left(\frac{\mathrm{dB}}{\mathrm{s}}\right)$

where s is the constant that depends on elevation angle and climate. We calculate the slope:

(61)$p(\zeta \mid A)=\frac{2}{\pi {\sigma }_{\zeta }{\left(1+{\left(\frac{\zeta }{{\sigma }_{\zeta }}\right)}^2\right)}^2}{\left(\frac{\mathrm{dB}}{\mathrm{s}}\right)}^{-1}$

The fade attenuation determination procedure by this model is presented in Algorithm 21.

Algorithm 21: ITU-R P.1623 model [98]

1  Calculate F /*$\phantom{(}$ Equation (59)$\phantom{(}$ */ 2  Calculate ${\sigma }_{\zeta }$ /*$\phantom{(}$ Equation (60)$\phantom{(}$ */ 3  Calculate $P\left(\zeta \right|A)$ /*$\phantom{(}$ Equation (61)$\phantom{(}$ */ Return A

Advantages

The ITU-R P.1623 model provides the fundamental idea for the fade slope model. Many other models have been developed using the ITU-R P.1623 model.

Limitations

The limitations of this model are that the applicable frequencies are limited to 10–30 GHz, and the elevation angles from 10–50${}^{\circ }$.

4.4.4. Dao Model

This model [99] for the fade slope was developed based on the VDK model, and a modification of the ITU-R P.1623 [98] model was proposed due to various parameters playing a dominant role in the fade slope. The datasets that were used were collected through an experimental setup in Kualalumpur, Malaysia. Details about the experimental setup description are available in [99]. The scintillation is filtered out from the attenuation through a moving average filter, whose cutoff frequency is defined by:

(62)$f_B=\frac{0.445}{t_a}$

After the scintillation is filtered out, the fade slope is calculated according to:

(63)$\zeta \left(t\right)=\frac{A\left(t+\frac{1}{2}\Delta t\right)-A\left(t-\frac{1}{2}\Delta t\right)}{\Delta t}\frac{\mathrm{dB}}{\mathrm{s}}$

where $\Delta t$ is the time interval length and t is sample number. It is noted that the climatic parameters need to be defined for the tropical regions as:

(64)$s=\frac{{\sigma }_{\zeta }}{F\left(f_B,\Delta t\right)A}$

Further, this assumption (Equation (64)) has been used to modify the ITU-R model:

(65)${\sigma }_{\zeta }=s(\Delta t)·F\left(f_B,\Delta t\right)·A·\mathbb{Q}\left(\frac{\mathrm{dB}}{\mathrm{s}}\right)$

where $\mathbb{Q}$ is used as the “adjusting parameter” with the standard deviation corresponding to the measured data-matching capability.

(66)$\mathbb{Q}=0.0035A^3-0.0868A^2+0.567A+0.1792$

The detailed steps of the model are given in Algorithm 22.

Algorithm 22: Dao model [99]

1. Apply MAF /*$\phantom{(}$ Equation (62)$\phantom{(}$ */ 2. Generate RN /*$\phantom{(}$ RN: Random number, between 0.1 to 60 s $\phantom{(}$ */ 3. Find $\mathbb{Q}$ /*$\phantom{(}$ Equation (66)$\phantom{(}$ */ 4. Estimate attenuation A 5. Return A

Advantages

The model proposed a fade slope determining technique with measured data from the tropical environment based on the nonlinear behavior of fade slope distribution.

Limitations

The model was tested only in Malaysia at 10.982 GHz. It has not been tested in other tropical climate conditions.

4.5. Learning-Based Model

4.5.1. Ahuna Model

The basis of this model [100] is the backpropagation neural network (BPNN), with the SST and ITU-R models. The BPNN network is composed of three input layers (I), three hidden layers (H), and one output layer (O), where a sigmoid function is used as the activation function, and the whole structure can be considered a “black box.” While backpropagating the signal through the network, the derivative of the error can be defined as:

(67)$\frac{\partial E_{\mathrm{total}}}{\partial O_a}=-\left(O_t-O_a\right)$

where $O_a$ is the actual output and $O_t$ is the desired output. In this study, the gradient descent algorithm was used to update the weights of the neurons as:

(68)$\Delta w_i=\eta \frac{\partial E}{\partial O_i}i=1,2,\dots ,I$

where $\eta$ is the learning rate, i is the input, $\Delta w_i$ is the weight change on the i-th input and $O_i$ is the output contributed by the i-th input. Finally, the minimum error gradient was achieved with a correlation of actual to predicted outputs at 91%, the optimized weighting and bias matrices were obtained as follows:

(69)$\begin{array}{cc}\hfill w_{\mathrm{HI}}& =\left[\begin{array}{ccc}-0.1974& -1.1467& -7.8568\\ -0.0115& 0.0086& -0.0476\\ 4.0797& 6.1838& 6.0543\end{array}\right]\hfill \end{array}$

(70)$\begin{array}{cc}\hfill w_{\mathrm{OH}}& =\left[\begin{array}{ccc}-78.4395\hfill & -32.0736\hfill & 31.939\hfill \end{array}\right]\hfill \end{array}$

(71)$\begin{array}{cc}\hfill b_{\mathrm{H}}& =\left[\begin{array}{c}-11.4485\\ 0.4845\\ -17.3805\end{array}\right]\hfill \end{array}$

(72)$\begin{array}{cc}\hfill b_0& =[-32.1027]\hfill \end{array}$

While millimeter-wave frequency information is unavailable from the models, the regional climatic parameter is simply the rainfall rate obtained from local rainfall data used to train the BPNN networks. The model is not validated with the standard DBSG3 or CCIR rain databanks.

Advantages

According to the validation outcome, the model can predict the rain attenuation at the accepted accuracy level.

Limitations

This model used rainfall data to predict future rainfall rate and hence the rain attenuation. Nevertheless, the rain attenuation does not depend on the rainfall rate only. There is research scope to add frequency, polarization, and rainfall distribution to enhance the model. However, in the BPNN system description, the name of the three input parameters was not documented.

5. Comparative Study of Slant Link Rain Fade Models

In the previous sections, we discussed the rain fade model prediction process. However, we may be interested in finding out some of the qualitative and quantitative characteristics of these models. Rain attenuation models for telecommunication links that support a wide variety of frequency spectra are a positive force. However, it is recommended to use a maximum 40 GHz frequency, as the next higher frequency band (above 40 GHz) may be used for astronomical and space communication purposes [137]. A 3–40 GHz frequency range is recommended for use in the telecommunication purposed Earth–space link to avoid possible interference with the deep-space communication link [138]. Table 7 presents the properties of slant link rain fade models such as whether they support terrestrial links; whether they consider rain structure including cloud, rainfall rate, frequency, elevation angle, rain layer height; and whether the models take into account the melting layer length within the effective path length. Table 8 shows the frequency spectrum supporting capabilities, the polarization used to evaluate the model, whether the model has a regional climatic parameter, whether the model considers rainfall rate distribution in a non-homogeneous means, and the recorded geographic area that suits the model properties. Table 9 shows whether these models were validated and the sources of rain attenuation related to measured databanks, constraints, and the complexity levels of the models in terms of low, medium, and high levels are classified. To validate the models, different validation tools such as the error figure and goodness-of-fit function, which is also known as relative error probability function, were mostly used [17,68,139].Figure 6 shows the frequency range supporting capacity of the studied models. As we can see, the maximum frequency range supporting model is the unified model [89], the second highest is the Crane TC model [94], and the third top model is the GSST model [66]. The lowest frequency band supporting models from lowest to highest bands are Singapore [33], the Breakpoint model [29], the LU model [68], and the Karasawa model [88].

6. Future Research Scope

6.1. Scaling Improvement

In satellite communication, scaling is the property that a model can be scaled for different input values. A model supporting a scaling property can be used for parameters with different frequencies, polarization techniques, elevation angle, and path length scaling. A model developed in a temperate climate region cannot be directly applied to a tropical area unless modified somehow. The proper scaling mechanism regarding these parameters is yet to be devised. The ITU-R model is the most significant model that has been moderated for tropical regions by many researchers [29,142,143,144,145,146]. However, these models do not support the scaling properties. There is research scope for adding various scaling properties to these models.

6.2. Non-Uniformity of Isothermal Heights

Accurate estimation of the slant path length in the rain structure up to the cloud is essential. Unfortunately, this part exists in the vertical space and is hard to estimate. ITU-R [38,41] has recommended that rain height is to be used only where another measured estimation is not available. Including ITU-R suggested techniques, many other rain height determination techniques are presented in Table 3. However, comparisons between these techniques are unavailable for different climatic zones. In this regard, there is research scope in a particular climatic region to compare rain height determination techniques or propose a better rain height estimation procedure.

6.3. Spatial Rainfall Distribution Along Slant Link

Using spatial interpolation techniques, the unknown rainfall rate values for the locations of unavailable data are estimated. The distance between each known point (with spatial coordinate ($x_i,y_i,z_i$) and the corresponding attribute is $P_i$), where $i\in N$ and the unknown point (with a spatial location of ($x_0,y_0,z_0$)) is $d_i(x,y,z)$. IDW is used to estimate the property value of each unknown point. Equation (73) represents the estimation procedure through the IDW method, where m is the power in the inverse distance power law [147].

(73)$P_R=\sum _{i=1}^N\frac{P_i}{{\left[d_i(x,y,z)\right]}^m}/\sum _{i=1}^N\frac{1}{{\left[d_i(x,y,z)\right]}^m}$

where distance is expressed as:

(74)$d_i=\sqrt[p]{{\left|x_0-x_i\right|}^p+{\left|y_0-y_i\right|}^p+{\left|z_0-z_i\right|}^p}$

The value p depends on the way the distance is considered. For example, in Euclidean distance, the value of p is given by 2. The IDW technique has been tested for a terrestrial link [148,149], but it has not been tested for rainfall estimation in slant links. Some other methods include modified inverse distance weighting [150,151], correlation coefficient weighting [150], multiple linear regression [152], and artificial neural networks [153], which can be examined for estimating the missing rainfall rate data. Now, for slant link rainfall estimation, research may be carried out based on these approaches to find suitable rainfall distribution along with the slant link.

7. Conclusions

In this review we have considered the well-known and most recent comprehensive surveys of the rain attenuation prediction models for the Earth–space link. The existing models have been classified as statistical, physical, empirical, learning-based, and fade slope models according to the model development and formulations. For these models, different aspects such as rain regions, rain structure, rainfall rate, elevation angle, rain layer height, melting layer height, frequency range, and polarization were considered. In this survey, we have revised more than 23 rain attenuation models for satellite links addressing their advantages and limitations. Accordingly, some models are good in the specific environments for which they have been developed and may not be well-performing in other geographical locations. Therefore, it is imperative to determine the attenuation parameters experimentally for different climate areas. We hope that this survey will inspire researchers to create a precise rain fade model for the slant link, either regionally or worldwide. The comparative analysis will assist people employed with slant relation architecture, budget making, and spreadsheets.

Author Contributions

M.A.S. anticipated, reviewed the related literature, evaluated, and outlined the rain fade models for Earth–space links. D.-Y.C. played a major role in coordinating research. M.A.S. conscripted the paper and subsequently modified and justified it by D.-Y.C. F.D.D. contributed to polishing the manuscript. All authors have read and agreed to the published version of the manuscript.

Funding

The BrainKorea21Four Program supported this research through the National Research Foundation of Korea (NRF) funded by the Ministry of Education (4299990114316). Additionally, this research was supported by the Basic Science Research Program through the National Research Foundation of Korea (NRF) funded by the Ministry of Education (2019R1F1A1058128).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

The writers would like to thank the editor and anonymous reviewers for their useful comments for improving the quality of this paper.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

The following abbreviations are used in this manuscript:

2-D SST	2-Dimentional SST
ARIMA	Autoregressive integrated moving average
BBH	Bright band height
BPNN	Backpropagation neural network
CCDF	Complementary cumulative distribution function
CCIR	Comité Consultatif International des Radiocommunications
CDF	Cumulative distribution function
DBSG3	Study group 3 databanks
DSD	Drop size distribution
ECMWF	European Center for Medium-Range Weather Forecasts
EMPI	Empirical
ERA-15	ECMWF Re-Analysis-15
ESST	Enhanced synthetic storm technique
EXCELL	EXponential CELL
FADE	Fade slope
FFT	Fast Fourier transform
FSL	Free space path loss
GA	Genetic algorithm
GE23	General Electric 23
GEO	Geostationary
GRV	Gaussian random variable
GSST	Global synthetic storm technique
HYCELL	Hybrid CELL
IDW	Inverse distance weighting
INTELSAT	International Telecommunications Satellite Organization
IoT	Internet of Things
ITU-R	International Telecommunication Union-Radiocommunication sector
LB	Learning-based
LEO	Low-Earth orbit
LPF	Low-pass filter
MAF	Moving average filter
MARIMA	Modified genetic algorithm-ARIMA
MEO	Medium-Earth orbit
M-K	Maseng–Bakken
MultiEXCELL	Multi-EXCELL
NM	Not mentioned
PDF	Probability density function
PEARP	Prévision d’Ensemble ARPEGE
PHY	Physical
PLF	Path length factor
PM	Physical-mathematical
POR	Pacific Ocean Region
PRF	Path reduction factor
RF	Reduction factor
RHCP	Right-hand circular polarization
RMS	Root mean square
RMSE	Root mean square error
RN	Random number
RSL	Received signal level
RV	Random variable
SAM	Simple rain attenuation model
SCEXCELL	Stratiform convective exponential CELL model
SNR	Signal-to-noise ratio
SST	Synthetic storm technique
STAT	Statistical
STD	Standard deviation
TC	Two component
TS	Time series
UAV	Unmanned aerial vehicle
VDK	Van de Kamp
WINDS	Wideband InterNet-working engineering test and Demonstration Satellite
Meanings of Used Symbols
$\Delta x_0$	A shift along horizontal axis due to the presence of layer B (Figure 2)
${\Gamma }_{\tau }\left(R\right)$	Rainfall rate conversion factor
${\gamma }_{decay}$	Decay profile along horizontal axis [Equation (16)]
$\gamma$	Specific attenuation
$\Lambda (\cdots )$	Log-normal distribution
${\lambda }_0$	Free-space wavelength (m) $\left({\lambda }_0=\frac{0.3}{f\left(GHz\right)}\right)$
$\rho$	Rain cell radius
$\theta$	Elevation angle in degrees of the earth station
$\phi$	Latitude of the earth station (degrees)
$\zeta$	Slant path
$v_{0.01}$	Vertical adjustment factor (Algorithm 7)
$a\left(R\right)$	Horizontal length of rain cell
$a,b$	Coefficient defined in “characteristic length” (Algorithm 6)
$A_B$	Breakpoint attenuation [Equation (23)]
$A_p$	Rain attenuation [Equation (13)]
$A_s$	Slant path attenuation
$b\left(R\right)$	Slant length of rain cell
B	Bandwidth (Hz)
${CF}_{60}$	Rainfall rate conversion factor
D	Rain cell diameter [Equation (36)]
$G_{rx}$	Receive antenna gain
$G_{tx}$	Transmit antenna gain
$H_R^{\prime }$	Top of melting level
$H_R$	Rain height
$H_S$	Height above mean sea level of the earth station (km)
$H_0$	$0^{\circ }$ isotherm height
$I_x,I_y$	Central moments of inertia
k	Boltzmann constant $=1.38\times {10}^{-23}$  J / K
$L_0$	Characteristic length
$L_{avg}$	Average long term slant path in the precipitation
$L_{eff}$	Effective path length
$L_G$	Horizontal projection length in the precipitation
$L_{sys}$	System loss at the receiver and transmitter
$L_H$	Projected path-length (Table 4)
l	Link distance (m)
L	Losses due to the presence of atmospheric gases, clouds, and fogs
$N_{eff}$	Effective number of rain cells (Table 4)
N	Number of tips [Equation (10)]
$P_{\tau -min}$	Probability of the mean rainfall rate is exceeded for $\tau$-min
$P_{1-min}$	Probability of the mean rainfall rate is exceeded for 1-min
$P_t$	Transmitter power expressed in dBm)
$\overline{R}$	Average rainfall rate
$r\left(p\right)$	Rain rate adjustment
$R\left(p\right)$	Rain rate exceeded for p% of an average year
$R_0$	Boundary rain rate
$R_M$	Local peak rainfall rate
$R_{max}\left(p\right)$	Maximum rain rate for p% of an average year
$r_p$	Path length reduction factor [Equation (13)]
$R_p$	Point rainfall rate [Equation (10)]
$R_p$	Point rainfall rate exceeded at $p\%$ of the time [Equation (13)]
${R_{rms}}^2$	Root mean square (RMS) of rainfall rate
$R_{tip}$	Rainfall per tip (mm) [Equation (10)]
$R_1$, $R_2$	Rainfall rate
T	Noise temperature (K) of the system which is assumed to be 290 K
T	Time gap in consecutive tips [Equation (10)]
$T_a^{\prime \prime }$	Under rainy conditions temperature
$T_a^{\prime }$	Under clear sky temperature
${TR}_T$	Maximum rainfall in mm for time interval T-min
u	Empirical constant (Table 5)
v (m/s)	Advection velocity of rain cells
$x_0$	Location of the ground station
$x_0,y_0$	Rain cell gravity center
$Z_{rain}$	Reflected signal in rainy condition
$Z_{th}$	Reflected signal in clear sky condition

Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Figures and Tables

Figure 1. Rain cell size for different rainfall rates: (a) R < 31.6 mm/h , (b) R 31.6 mm/h [43].

Figure 2. Illustration of two-layer structured rain attenuation model for the slant link.

Figure 3. Taxonomy of the rain fade models for Earth–space links [88,89,90,91,92,93,94,95,96,97,98,99,100].

Figure 4. Memoryless non-linear device used in M-K model, where β is the dynamic parameter of the model, which signifies the rate of change of attenuation; g is the gain; X(t) is the low-pass filtered output; and γ(t) is the attenuation.

Figure 5. Illustration of the fade slope and fade duration in the received signal level (RSL) time series [98].

Figure 6. Each“vertical bar” represents the supported frequency range of the studied models. The vertical axis represents frequency in GHz and the horizontal axis represents the models.

Table 1

Rain attenuation survey papers.

Ref.	Survey Concentrations
[14]	The results show that rain attenuation in horizontal polarization is slightly larger than that in vertical polarization. A differential coupling-based compensation was recommended to reduce the rain attenuation effects on polarization in the depolarization processes.
[15]	In this survey, the artificial-neural-network-based models were classified based on input parameters. In addition, the accuracy of the artificial-neural-network-based models was accessed through a comparative study.
[16]	Rain attenuation in a satellite system is more severe than in a terrestrial system, and for a good model, it should be implemented into the model through local climatic elements.

Table 2

Similarity and contrast between Earth–space to terrestrial links.

Factors	Earth–Space	Terrestrial
Rain	✓	✓
Cloud, sky noise, scintillation, latitude, Earth radius	✓	✗
Hop length	Varies according to orbit: LEO: ≈500–1000 km MEO: ≈8000–12,000 km GEO: ≈36,870 km	Varies according to necessity: Typical: <100 km Rare: Maximum 200 km
Path length reduction	Effective slant path	Effective terrestrial path

Table 3

Different techniques to estimate rain height.

Ref.	Mechanism	Symbol Meanings
[35]	$\left\{\begin{array}{ccc}{H}_0\hfill & \left(\mathrm{km}\right)\hfill & R\le 10\mathrm{mm}/\mathrm{h}\hfill \\ H_0+log(R/10)\hfill & \left(\mathrm{km}\right)\hfill & R>10\mathrm{mm}/\mathrm{h}\hfill \end{array}\right.$	$H_0$ is the $0^{\circ }$ isotherm height, and R is the rain rate.
[36]	$10log\left(\frac{T_m-T_a^{\prime }}{T_m-T_a^{\prime \prime }}\right)·{\gamma }^{-1}$	$T_m$ is the medium temperature, $T_a^{\prime }$, under clear sky, and $T_a^{″}$ is under rainy condition temperatures, and $\gamma$ is the specific attenuation coefficient.
[37]	$4.5+0.0005R^{1.65}\left(\mathrm{km}\right)$	R is the rain rate in the range 30–70 mm/h.
[38]	$H_0+0.36\left(\mathrm{km}\right)$, also using rain height map	$H_0$ is the $0^{\circ }$ isotherm height.
[32]	$\left\{\begin{array}{cc}4.8,& \left|\phi \right|\le {30}^{\circ }\mathrm{and}R\le 10\mathrm{mm}/\mathrm{h}\hfill \\ 7.8-0.1\left|\phi \right|,& \left|\phi \right|\ge {30}^{\circ }\mathrm{and}R\ge 10\mathrm{mm}/\mathrm{h}\hfill \end{array}\right.$	$\phi$ is the latitude in degrees.
[33]	$\left\{\begin{array}{cc}3+0.028\phi & 0\le \phi <{36}^{\circ }\\ 4-0.075(\phi -{36}^{\circ })& \phi \ge {36}^{\circ }\end{array}\right.$	$\phi$ is the latitude in degrees.
[34]	$\left\{\begin{array}{cc}4& 0<\left|\phi \right|<{36}^{\circ }\\ 4-0.075\left(\left|\phi \right|-{36}^{\circ }\right)& \left|\phi \right|\ge {36}^{\circ }\end{array}\right.$	$\phi$ is the latitude in degrees.
[39]	$14.114·R^{-0.6402}$	R is the rain rate.
[40]	$3.6693·R^{0.097}$	R is the rain rate.
[41]	$\left\{\begin{array}{cccc}5-0.075(\phi -{23}^{\circ })\hfill & \mathrm{for}\phi >{23}^{\circ }\hfill & & \\ 5\hfill & \mathrm{for}0\le \phi \le {23}^{\circ }\hfill & & \\ 5\hfill & \mathrm{for}0\ge \phi \ge -{21}^{\circ }\hfill & & \\ 5+0.1(\phi +{21}^{\circ })\hfill & \mathrm{for}-{71}^{\circ }\le \phi <-{21}^{\circ }\hfill & & \\ 0\hfill & \mathrm{for}\phi <-{71}^{\circ }\hfill & \end{array}\right.$	$\phi$ is the latitude in degrees.

Table 4

Effective path length estimation techniques.

Ref.	Mathematical Procedure
white [35]	$L_{eff}=13.367R^{-0.21}$, where R is the rain rate.
[42]	$L_{eff}=\frac{1}{1+\frac{L_{S}cos\theta }{L_G}}·L_S$, where $\theta$ is the elevation angle, $H_S$ is the station height above mean sea level of the Earth, and $L_G$ is the horizontal projection length in the precipitation.
[43]	$L_{eff}=\frac{L_H{N}_{eff}}{cos\theta }\left(\mathrm{km}\right)$ where $L_H$ is the projected path-length, which is defined by $R_0$, a boundary rain rate. If $R<R_0$ then $L_H=\frac{H_R-H_s}{tan\theta }\left(\mathrm{km}\right)$, and if $R\ge R_0$ then $L_H=a\left(R\right)\left(\mathrm{km}\right)$ (Figure 1), and ‘$N_{eff}$’ is the effective number of rain cells given by $N_{eff}\left(\theta ,R\right)=n_1\left(\theta \right)+n_2\left(\theta \right)·R$, where $n_1\left(\theta \right)=2.38{\theta }^{-1.45}+0.82$, and $n_2\left(\theta \right)=-0.01{\theta }^{0.52}+0.09$ ($\theta$ is elevation angle).
[44]	The estimation of the rain path can be determined by calibrating the reflected signal $Z_{th}$ in clear-sky conditions and in rainy conditions $Z_{rain}$. Once rainfall events are identified the first step of the procedure involves the determination of the volume of the rain cell with $Z_{rain}>Z_{th}$.

Table 5

Rain rate conversion techniques.

Ref.	Conversion Techniques
[49]	Using the Rain rate ${\Gamma }_{\tau }\left(R\right)$, where ${\Gamma }_{\tau }\left(R\right)=P_{1-\min }\left(R\right)/P_{\tau -\min }\left(R\right)$, the model can convert rain between 1- and $\tau$-min, where $P_{1-\min }$ is the probability that the mean rainfall rate is exceeded for a 1 min accumulation time, and $P_{\tau -\min }$ is the probability that the mean rainfall rate is exceeded for $\tau$ min accumulation time.
[50]	Using cumulative rainfall distributions, the model can convert rain between 1- and $\tau$-min.
[51]	Using the conversion factor $R_1\left(P\right)/R_{60}\left(P\right)={CF}_{60}$, where ($C{F}_{60}=a{P}^b+c{e}^{\left(d·P\right)}$) the model can convert 60 min to a 1 min integration time rain rate, where $C{F}_{60}$ is the conversion factor, a, b, c, d are constants, and P is the probability of rainfall rate exceedance within $0.001\%\le P\le 1\%$. It can also be used to convert between different rain rates [52].
[53]	Using the equation ($P_1=1.71·P_{\tau }·{10}^{\left(-0.242·exp(\tau /2428)\right)}$), $\tau$ min to 1 min conversion.
[54]	Using the equation $R_{1-\min }={\left(R_{\tau -\min }\right)}^d$, where $d=0.987{\left(\tau \right)}^{0.061}$ the model can convert the rain rate from $\tau$ min to 1 min.
[55]	Using the equation, $log\left[R_{1-\min }\left(P\right)]={\alpha }_r·log[R_{\tau -\min }\left(P\right)\right]$, where ${\alpha }_r$: regression component determined statistically, and $5.0\%\le P\le 0.01\%$.
[56]	In the absence of accurate long-term rainfall data, this method takes the percent of time exceedance, latitude, and longitude as input to generate a 1 min integration time rain rate.
[57]	The model simply calculates the rain rate $R=60·{TR}_T/T$, where ${TR}_T$: maximum rainfall in mm for time interval T min.
[58]	$P_1\left(R\right)=\left(a_1/a_t\right)R^{\left(b-b_1\right)}exp\left[\left(u_t-u_1\right)R\right]P_t\left(R\right)$, R: rain rate, u: empirical constant, a, b: depending on the location and integration time of the rain gauge.
[59]	$R_{0.01\%}=\frac{ln\left(0.03\beta \frac{M}{T}\right)}{0.03}$ where M is the average annual rainfall (mm), and T is the number of hours in the year during which rain rates exceed R (mm/h).
[60]	$R_p=a_p{M}^{b_p}{\beta }^{c_p}$ where p is the percentage of the time, M is the total annual rainfall accumulation, $a_p$, $b_p$ and $c_p$ are constants, which can be derived from $x=log\left(P\right)$, and $\beta$ is the thunderstorm ratio ($\beta =M_1/M_1+M_2$), where $M_1$ is the rain rate due to a thunderstorm and $M_2$ is the rain rate without a thunderstorm.
[61]	$R_{0.001\mathrm{initial}}=23.50390R^{0.27896}\mathrm{and}{\beta }^{0.34162}$$R_{0.01\mathrm{initial}}=3.69786R^{0.46613}{\beta }^{0.43482}$ where R is the rain rate in mm/h and $\beta$ is the regional climatic parameter.

Table 6

Parameters of EXCELL model.

Parameter	Monoaxial Rain Cell	Biaxial Rain Cell
Peak rainfall rate ($R_M$)	✓	✓
Cell radius/radii	${\rho }_0$	${\rho }_x,{\rho }_y$

Table 7

Characteristics table.

Tpe	Ref.	Slant or Both	Rain Structure	Rainfall Rate	Frequency	Angle(Earth Station)	Rain Height	Melting Layer Height	Time Series
EMPI	[32]	Slant	✓	✓	✓	✓	✗	✓	✓
	[33]	Slant	✗	✓	✓ ${}^{‡}$	✓${}^{⋏}$	✓	✗	✗
	[34]	Slant	✓	✓	✓	✓	✓	✗	✗
	[66]	Slant	✓	✓	✓	✓ ${}^{\diamond }$	✓	✓	✓
	[68]	Slant	✓ ${}^{◃}$	✗	✗	✗	✓	✓	✗
	[88]	Slant	✓	✓	✗	✓	✓ ${}^{▹}$	✗	✗
	[29]	Slant	✗	✓	✓ ${}^{‡}$	✓	✓	✗	✗
	[89]	Both	✗ ${}^{⋉}$	✓	✓	✓ ${}^{⋊}$	✗	✗	✗
	[12]	Slant	✗	✓	✓ ${}^{‡}$	✓ ${}^{‡}$	✓	✗ ${}^{⋉}$	✗
STAT	[90]	Slant	✗	✓	✓ ${}^{⊺}$	✗	✗	✗	✗
	[91]	Slant	✗	✓	✓	✓	✓	✗	✓
	[45]	Slant	✗	✓	✓	✓	✓	✗	✗
	[92]	Slant	✗	✓	✓ ${}^{\cap }$	✓ ${}^{\cap }$	✗	✗	✗
	[93]	Slant	✗	✗ ${}^{\cup }$	✗	✗	✗	✗	✓
PHY	[37]	Slant	✓	✓	✓	✓	✓	✗	✗
	[94]	Both	— ${}^{⋗}$	✓	— ${}^{\sqcap }$	✓	✓	✗	✗
	[95]	Both	✓ ${}^{◃}$	✓	✗	✓	✓	✓	✓
	[62]	Slant	✓	✓	✓	✓	✓	✓	✗
FADE	[96]	Slant	✗	✗	✓	✓	✗	✗	✗
	[97]	Slant	✗	✗	✗	✗	✗	✗	✓
	[98]	Slant	✗	✗	✓ ${}^{\bigsqcup }$	✓	✗	✗	✗
	[99]	Slant	✗	✗	— ${}^{\oplus }$	✓ ${}^{\ominus }$	✗	✗	✗
LB	[100]	Slant	✓	✓	✓	✓	✓	✓	✓

EMPI: empirical, STAT: statistical, PHY: physical, FADE: fade slope, LB: learning-based, ${}^{‡}$ k and $\alpha$, ${}^{⋏}$ 5–50${}^{\circ }$, ${}^{\diamond }$ 20–60${}^{\circ }$, ${}^{◃}$ ground layer and melting layer, ${}^{▹}$ exponential rain cell, ${}^{⋉}$ only rain height, no melting layer, ${}^{⋊}$ 6.50 72.8${}^{\circ }$. ${}^{⊺}$ as $\alpha$ is dependent on frequency (not polarization), ${}^{\cap }$ mentioned but no mathematical form is found, ${}^{\cup }$ rain attenuation is used, ${}^{⋖}$ cylindrical rain cell, ${}^{⋗}$ the rain structure is composed of heavy and light rainfall components, ${}^{\sqcap }$ frequency parameter is not directly used in the calculation of the attenuation, ${}^{\bigsqcup }$ as the applicable frequency range is limited to 10–50 GHz, ${}^{\oplus }$ not mentioned, ${}^{\ominus }$ s-constant that depends on elevation angle and climate.

Table 8

Characteristics table.

Type	Ref.	Frequency Range	Polarization	Regional Climatic Parameter	Spatial Friendly	Reported Region that Suits
EMPI	[32]	10–35 GHz	✓	✓	✓	The USA, Europe, and Japan
	[33]	4–6 GHz	NM ${}^{\nparallel }$	Rainfall rate	✗	Singapore
	[34]	NM	✓	4 coefficient constants called a,b,c and d depend on geographic area	NM	Best suited areas are in Europe, the USA, Japan, and Australia
	[66]	10–100 GHz	✓Circular, horizontal and vertical	Rainfall rate	✗	Temperate region
	[68]	Above 11 GHz	✗	Local rainfall rate	No information available	No practical information; test information is limited within DBSG3 databank
	[88]	10–20 GHz	NM	Rainfall rate, the probabilities of occurrence and mean rainfalls, for cell and debris	The model considers a single volume rainfalling area and does not consider convective (cell) and stratified rain (debris)	CCIR rain zone
	[29]	10.7–13 GHz (Ku band)	✓	Rainfall rate	NM	Fiji
	[89]	Terrestrial: 7–137 GHz	✓k and $\alpha$	Rainfall rate	✓	Worldwide (ITU-R databank)
	[12]	10–30 GHz	Both	Rainfall rate	✗	Tropical
STAT	[90]	NM	✓	Rainfall rate	✗	Temperate region
	[91]	10–50 GHz	✓	Climate characteristics, link elevation angle	Has spatial dynamics of rainfall rate parameter due to Van de Kamp [127]	Xi’an, China
	[45]	Up to 55 GHz	✓	Rainfall rate	✗	Worldwide
	[92]	NM (only 20 GHz is mentioned)	NM	NM	Information NM	France
	[93]	NM	✓Horizontal	Measured rain attenuation	NM	Kolkata, India
PHY	[37]	Reported frequencies are 15, 30, 39.6 GHz	NM	✓	Rainfall rate is not a function of distance	Temperate and tropical
	[94]	1–100 GHz	✗	Cell and debris parameter	No spatial correlation function is included	Temperate region
	[95]	11.6 GHz (Range NM)	✓Circular, Linear	Rainfall rate, wind velocity	✓(1-year rainfall rate PDF was used)	Italy
	[62]	10–50 GHz	✓	Rainfall rate; rain height; melting layer height, etc.	✓Vertical and horizontal	Worldwide
FADE	[96]	Ku-band	RHCP, Linear polarization	✗	✗	Japan (low-elevation Earth–space path)
	[97]	NM	Horizontal	✓	✗	Tropical area
	[98]	10–50 GHz	NM	✓ ${}^{\bowtie }$	✗	Designed for global perspective
	[99]	10.982 GHz (Ku-band)	Vertical	s-parameter ${}^{⊡}$	✗	Malaysia
LB	[100]	NM	NM	Local rainfall data: it needs to train the BPNN networks	NM	Butare, Rwanda

EMPI: empirical, STAT: statistical, PHY: physical, FADE: fade slope, LB: learning-based rain attenuation model, ${}^{\nparallel }$ not mentioned, ${}^{\bowtie }$ elevation angle, ${}^{⊡}$ s-parameter definition is given in Equation (60).

Table 9

Characteristics table.

Type	Ref.	Used/Validation Databank	Remarks about Validation	Complexity Level	Constraints
EMPI	[32]	62 experiments in the USA, Japan and Europe	Shows better prediction with minimum 2 years rainfall datasets	low	The operating frequency range is limited to 10–35 GHz
	[33]	INTELSAT POR experiment (Singapore)	The results have been validated with CCIR model and measured attenuation	low	The recommended frequency is limited to 4–6 GHz
	[34]	77 satellite links placed in Europe, the USA, Japan, and Australia (CCIR data bank)	Validated by the author	low	The model is suitable for prediction with an elevation angle ranging from 10${}^{\circ }$ to 40${}^{\circ }$
	[66]	ITU-R databank	ITU-R databank is used to examine the behavior of m parameter at different frequencies as well as different sites	low	Not spatial friendly
	[68]	DBSG3 databank	The RMS value show that the new model provides the smallest rms and STD in all percentages of time except at 0.001%	high	Application area is limited by low latitudes (between 36${}^{\circ }$ South and 36${}^{\circ }$ North) and low elevation angles (25${}^{\circ }$<)
	[88]	CCIR rain zone with different elevation angle	Validated by the author	medium ∔	10–20 GHz
	[29]	Used 6 tropical city’s rain databanks from ITU-R	The experimental result agreed about correction factor induced attenuation especially at elevation angle less than 60${}^{\circ }$ compared to ITU-R model	medium	In this model it is assumed that rain is uniformly distributed inside a rain cell
	[89]	ITU-R databank	The model gives decent results with the correctness of its terrestrial and slant direction	medium	The model was not checked with actual data for attenuation
	[12]	✗	Verified with the the beacon signals from WINDS and GE23 satellites (2009 to 2012)	low	Up to 30 GHz frequency is examined
STAT	[90]	✗	Not validated with heavy rainfall or with rainfall data around the world for different sites	low	The model is applicable only during the rain periods because no transition is included in the model to switch from rainy to clear sky conditions
	[91]	Validated based on the measured rain attenuation data at Xi’an, China	The validation is a good agreement between measured and predicted attenuation using this model, but has not been compared with other important rain attenuation models	high	The size of the area is not defined for a specific set of parameters
	[92]	✗	Not validated	—	The probabilistic weather forecasts could be beneficial to maximize the economic value accounting the transmitted data for higher frequencies (say 50 GHz).
	[93]	✗: used measured attenuation data	The validation shows a good agreement between measured and predicted attenuation using this model, but has not been compared with other important rain attenuation models	medium	It needs to derive mean and standard deviation of the Gaussian distribution for a different geographic area and find the coefficients of second-order polynomials from real measured attenuation
PHY	[37]	✓	Validated at Lae, Papua New Guinea	low	—
	[94]	Satellite: validated through CCIR procedure. Terrestrial: 35 terrestrial paths around the world	Both satellite and terrestrial links have been verified	high ${}^{†}$	The probabilities of incidence and mean precipitation for cells and debris are difficult to determine
	[95]	Experimental dataset	The experimental results show that the method show less error probability.	medium	Needs 1 min rainfall rate
	[62]	DBSG3 databank	Although the model does not outperform the existing ITU-R model (approximately 2.4% error), it supports a few additional facilities (e.g., site diversity) and takes account of the interference due to hydrometer scattering	high	The correctness is limited by the local values of input parameters like the melting layer and the rain plateau value, which might not be available everywhere
FADE	[96]	Experimental dataset	Validated with	low	—
	[97]	✗	Validated in Kolkata, India, with experimental set-up	low	To use fade margin data it needs to remove tropospheric scintillation
	[99]	Only experimental dataset	The resulted output agrees with the measured standard deviation of attenuation	low	The method was tested only at 10.982 GHz
	[98]	NM	Validated in [140,141]	low	Applicable for limited elevation angle 5–60${}^{\circ }$
LB	[100]	Durban, South Africa	The model was not validated with standard DSDB3 or CCIR rain databanks	—	The model was not tested with established well-known rain databanks like DBSG3 or CCIR

EMPI: empirical, STAT: statistical, ∔ it needs to solve at least seven equations. PHY: physical, FADE: fade slope, LB: learning-based rain attenuation model, ^† it needs to solve ten equations.