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

: 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 signiﬁcant for propagating electromagnetic waves to achieve the required quality of service. In this survey, different slant link rain attenuation prediction models have been examined, classiﬁed, and analyzed, and various features like improve-ments, drawbacks, and particular aspects of these models have been tabulated. This survey provides various techniques for obtaining input data sets, including rain height, efﬁcient 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.


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  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, Table 1. Rain attenuation survey papers. [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.

Ref. Survey Concentrations
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-todate 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.

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.

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.

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].
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: where L sys is system loss at the receiver and transmitter; λ 0 is free-space wavelength (m), which is determined from the frequency as λ 0 = 0.3 f (GHz) ; G tx and G rx are transmitting and receiving antenna gains, respectively; k is the Boltzmann constant 1.38 × 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]: N(dBm) = 174(dBm/Hz) − 10 log 10 B(Hz) Equation (2) can equivalently be computed in unit decibels (dBm) as SNR(dB) = P r (dBm) − 10 log 10 B(Hz) + 174(dBm/Hz) (6) 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]: 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]: where A 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 rain by total rain attenuation at a rainfall rate of exceedance of A 0.01 .

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 minute 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].

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.

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.

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) where A(dB) is the total attenuation due to rain, and specific attenuation of the link is γ, then the "effective path length" is defined by L e f f . Table 4 shows a list of additional techniques to determine 'effective path length. Table 3. Different techniques to estimate rain height.
[36] 10 log T m is the medium temperature, T a , under clear sky, and T a is under rainy condition temperatures, and γ is the specific attenuation coefficient.
[40] 3.6693 · R 0.097 R is the rain rate. [41] ϕ is the latitude in degrees. Table 4. Effective path length estimation techniques. [35] L e f f = 13.367R −0.21 , where R is the rain rate.

Ref. Mathematical Procedure
[42] L e f f = 1 1+ L S cos θ L G · L S , where θ 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.
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 = H R −H s tan θ (km), and if R ≥ R 0 then L H = a(R) (km) (Figure 1), and 'N e f f ' is the effective number of rain cells given by N e f f (θ, R) = n 1 (θ) + n 2 (θ) · R, where n 1 (θ) = 2.38θ −1.45 + 0.82, and n 2 (θ) = −0.01θ 0.52 + 0.09 (θ 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 .

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: 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.  [49] Using the Rain rate Γ τ (R), where Γ τ (R) = P 1−min (R)/P τ−min (R), the model can convert rain between 1-and τ-min, where P 1−min is the probability that the mean rainfall rate is exceeded for a 1 min accumulation time, and P τ−min is the probability that the mean rainfall rate is exceeded for τ min accumulation time.

Ref. Conversion Techniques
[50] Using cumulative rainfall distributions, the model can convert rain between 1-and τ-min. [51] Using the conversion factor R 1 (P)/R 60 (P) = CF 60 , where (CF 60 = aP b + ce (d·P) ) the model can convert 60 min to a 1 min integration time rain rate, where CF 60 is the conversion factor, a, b, c, d are constants, and P is the probability of rainfall rate exceedance within 0.001% ≤ P ≤ 1%. It can also be used to convert between different rain rates [52].
[54] Using the equation R 1−min = (R τ−min ) d , where d = 0.987(τ) 0.061 the model can convert the rain rate from τ min to 1 min. [55] Using the equation, log[R 1−min (P)] = α r · log[R τ−min (P)], where α r : regression component determined statistically, and 5.0% ≤ P ≤ 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] , R: rain rate, u: empirical constant, a, b: depending on the location and integration time of the rain gauge.
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 β 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(P), and β is the thunderstorm ratio (β = 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 initial = 23.50390R 0.27896 and β 0.34162 R 0.01 initial = 3.69786R 0.46613 β 0.43482 where R is the rain rate in mm/h and β is the regional climatic parameter.

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.

Rain Cell Size
In [43], the horizontal length a(R) and the slant length b(R) (R < 31.6 mm/h case only) estimation technique is proposed as: where R is the rainfall rate and dimension parameters a(R) and b(R) are shown in Figure 1.  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.

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 (R), root mean square (RMS) of rainfall rate (R 2 rms ), 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 (ρ) 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 2 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: Table 6. Parameters of EXCELL model.

Monoaxial Rain Cell Biaxial Rain Cell
Peak rainfall rate (R M ) Cell radius/radii ρ 0 ρ x , ρ y 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 ≤ 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 2 and the other distinct parameters were similar to the EXCELL model (R M , ρ 0 ) and additional number k = 1, 2 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(ρ 0 |R M ).

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: where A p is the rain attenuation, R p is the point rainfall rate exceeded at p% of the time, k and α 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): where a 1 , a 2 , a 3 , a 4 are constants that can be obtained through the curve fitting technique, L S = (H R − H S )/sinθ is the slant path in the precipitation (km), θ 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 • . It requires the advection velocity v (m/s) of rain, the length of the satellite path ζ (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 is the melting layer top limit, and ϕ is the latitude of the Earth station (degrees), if not mentioned otherwise.

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 sizedistribution-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].

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.
(2) The rainfall rate varies along with the link, and the rainfall rate can be defined as: where R p is the point rainfall intensity, γ decay is the decay profile along the horizontal axis, and l ≤ L S is defined by Equation (17): where rain height, H R , is approximated by Equation (18): (3) Finally the rain attenuation is determined by Equation (19): Algorithm 1 shows the procedure for this model.

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.

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.
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 • to 40 • . The detailed steps of the model are given in Algorithm 3.

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.

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

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]
/* Valid for the frequency range 10 GHz ≤ f ≤ 100 GHz */ The distribution of rainfall rate projection on the Earth is defined as: where b = 1/ρ here ρ (km) is the same radius for which the rainfall rate is fixed R max /e.
(2) The attenuation along with the slant link: where r(p) is the rain rate adjustment factor, and R(p) is the rain rate exceeded for p% of an average year.
(3) The coefficients a i (i = 1, 2, ..., 6) are taken using the genetic algorithm (GA) and annealing algorithm; thus, exploiting the rain databank of the ITU-R for the slant link: The procedure to determine attenuation by this model is shown in Algorithm 5.

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.

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.

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. 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: where P 1 , and P 2 are given by Equations (25) and (26): The experimental results agreed with the ITU-R model [107] regarding correction factor-induced attenuation, especially at elevation angles less than 60 • . Algorithm 7 shows the procedure for this model.

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 • ≤ θ < 60 • .

Limitations
The correction factor that is proposed in this model has not been analytically validated. Invalid for tropical region  [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):

Algorithm 7: Breakpoint Model [29]
where L S and R e f f (R p , L S , θ) are given by Equations (28) and (29): R e f f (R p , L S , θ) = 1.74R 0.786+0.197/L S cos θ cos θ + 120 The working procedure of this model is shown in Algorithm 8.

Advantages
The presented results show that the model outperforms ITU-R P.618-12 [108]. 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. Elevation angle error/* Invalid for θ < 5 • elevation angles */

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.

Maseng-Bakken (M-K) Model
This model [90] is based on the log-normal distribution of rain and utilizes a 1dimensional 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): The symbols k and α 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 γ is log-normal distributed. Consequently, following the log-normal distribution, rainfall intensity can be written as: Λ(R m , σ R 2 ; R). The justification of the modeling of the rainfall intensity distribution is given in Equation (31): The conformity of mapping the distribution of rainfall intensity through γ can be accessed through determining the auto-correlation between the variables, as in Equation (32): After normalization and the removing the bias: where ρ(τ) = exp(−β|τ|) and the value of β can be determined as: where ∆t is the sampling time interval determined from the time series of γ. 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. 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.

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]
1 Prepare sub-databank with m-group parameters /* Using ARIMA (1,1,7) and k-NN */ 2 Initialize and define the population N /* Maximum population, N max = 2 m */ 3 Compute parameters of ARIMA /* Generate N groups of 1 × m arrays with 0/1 */ 4 Calculate the average of each "chromosome"; 2 A t means the real value of rain-induced attenuation, and ARIMA(A t−1 , ϕ, θ i ) represents the predicted value using the parameters based on one of N. */ 6 if f i > η then 7 Confirm parameters 8 else 9 Duplicate, Overlap and Mutate and then go to step 3: 10 Return rain attenuation time series;

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.

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

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.

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.

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.

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.

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

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.

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]
1 Choose the region and select parameters (P c ,R c , P D , R D , σ D , R 0.01 ) // Projected path rain cell and debris region 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.

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) An illustration of these parameters is given in Figure 2.
Along the horizontal axis, the speed of the wind is given by: Transforming Equation (45) into a space spectrum, and with further analysis, Equations (46) and (47), yield the attenuation for the elevation angles θ = 90 • and θ = 90 • , respectively.
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.

Advantages
It is the first model that comes with the practical melting and rain layer consideration and applying a 1 minute 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.

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 • isotherm height; β i : the ratio between the convective and the total rain amounts calculated over each successive 6 hour 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 hour 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).
Equation (48) is valid for the following conditions: initial density of the melting particles d 0 = 0.16 g/cm 3 ; particle size distribution (gamma with µ = 0); and frequency (10 GHz < f < 50 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.
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.

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.

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 ∆t to be 2 s. A typical fade slope diagram is shown in Figure 5.

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): where M is equal to 2m + 1 and ∆ is equal to 1 s.
The detailed steps of the Japan model are given in Algorithm 19.

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.

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 synthesizerbased technique [136]. Algorithm 20 shows the procedure for this model.

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.

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 • . 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 , ∆t): We then use F( f B , ∆t) to determine σ ζ : where s is the constant that depends on elevation angle and climate. We calculate the slope: The fade attenuation determination procedure by this model is presented in Algorithm 21.

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

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: After the scintillation is filtered out, the fade slope is calculated according to: where ∆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: Further, this assumption (Equation (64)) has been used to modify the ITU-R model: where Q is used as the "adjusting parameter" with the standard deviation corresponding to the measured data-matching capability.
The detailed steps of the model are given in Algorithm 22.

Algorithm 22: Dao model [99]
1 Apply MAF /* Equation (62)  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: 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: where η is the learning rate, i is the input, ∆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: 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.

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].  as α is dependent on frequency (not polarization), ∩ mentioned but no mathematical form is found, ∪ rain attenuation is used, cylindrical rain cell, the rain structure is composed of heavy and light rainfall components, frequency parameter is not directly used in the calculation of the attenuation, as the applicable frequency range is limited to 10-50 GHz, ⊕ not mentioned, s-constant that depends on elevation angle and climate.  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 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.

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.

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 ∈ 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].
where distance is expressed as: 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.

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. 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 conflicts of interest.

Abbreviations
The following abbreviations are used in this manuscript: Probability of the mean rainfall rate is exceeded for τ-min P 1−min Probability of the mean rainfall rate is exceeded for 1-min P t Transmitter power expressed in dBm) R Average rainfall rate r(p) Rain rate adjustment R(p) Rain rate exceeded for p% of an average year R 0 Boundary rain rate R M Local peak rainfall rate R max (p) Maximum rain rate for p% of an average year r p Path length reduction factor [Equation (13) (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