SciELO - Scientific Electronic Library Online

 
vol.86 número209Comparison of alumina ball size distribution in two white cement grinding units using Swebrec functionA new scheme for multiple fault detection and isolation for rotational mechatronic systems, by means of analytical redundancy and adaptive filtering índice de autoresíndice de assuntospesquisa de artigos
Home Pagelista alfabética de periódicos  

Serviços Personalizados

Artigo

Indicadores

Links relacionados

  • Em processo de indexaçãoCitado por Google
  • Não possue artigos similaresSimilares em SciELO
  • Em processo de indexaçãoSimilares em Google

Compartilhar


DYNA

versão impressa ISSN 0012-7353

Dyna rev.fac.nac.minas vol.86 no.209 Medellín abr./jun. 2019

http://dx.doi.org/10.15446/dyna.v86n209.73774 

Artículos

Rainfall rate and rain attenuation contour maps for preliminary “Simon Bolivar” satellite links planning in Venezuela

Mapas de contorno de tasa de precipitación y atenuación por lluvias para planificación preliminar de enlaces del satélite “Simón Bolívar” en Venezuela

Angel Dario Pinto-Mangones a   , Nelson Alexander Pérez-García b   , Juan Manuel Torres-Tovio a   , Eduardo José Ramírez c   , Samir Oswaldo Castaño-Rivera d   , Jaime Velez-Zapata e   , John Dwiht Ferreira-Rodríguez f   , Leidy Marian Rujano-Molina c  

a Escuela de Ingeniería de Sistemas, Universidad del Sinú, Montería, Colombia. angelpinto@unisinu.edu.co, juantorrest@unisinu.edu.co

b Escuela de Ingeniería Eléctrica, Universidad de Los Andes, Mérida, Venezuela. perezn@ula.ve

c Postgrado en Telecomunicaciones, Universidad de Los Andes, Mérida, Venezuela. edramirez952@gmail.com, lmarian89@hotmail.com

d Departamento de Ingeniería de Sistemas y Telecomunicaciones, Universidad de Córdoba, Montería, Colombia. sacastano@correo.unicordoba.edu.co

e Departamento de Ciencias de la Computación y Electrónica, Universidad de la Costa, Barranquilla, Colombia. jvelez@cuc.edu.co

f Departamento de Física, Universidad de Los Andes, Mérida, Venezuela. fjohn@ula.ve

Abstract

Predicting precipitation rate and rainfall attenuation are key aspects in planning and dimensioning of wireless communication systems operating at frequencies above 10 GHz, such as satellite communication systems at Ku and Ka-bands. In this paper, contour maps of rainfall rate and rain attenuation are developed for the first time in Venezuela. It is based on 1-min rain rate statistics obtained from measurements carried out in Venezuela over a 30 years period with a higher integration time. The Rice-Holmberg model, refined Moupfouma-Martin model and Recommendation ITU-R P.837-7 for rain rate estimation were used, while Recommendation ITU-R P.618-13, Ramachandran-Kumar model, Yeo-Lee-Ong model and Rakshit-Adhikari-Maitra model were used for rain attenuation prediction in “Simón Bolívar” satellite links in Venezuela. The results of both types of maps represent a useful tool for preliminary planning of those links in the country, specifically, in Ku and Ka-bands.

Key words: contour maps; rain rate prediction; integration-time conversion; rain attenuation estimation; satellite links

Resumen

La tasa de precipitación y atenuación por lluvias son aspectos claves en la planificación y dimensionamiento de sistemas inalámbricos de comunicaciones que operan en frecuencias superiores a 10 GHz, por ejemplo, sistemas de comunicación vía satélite en las bandas Ku y Ka. En este artículo, se desarrollan por primera vez en Venezuela mapas de contorno de tasa de precipitación y atenuación por lluvias, en base a estadísticas de lluvia de 1-minuto obtenidas a partir de mediciones realizadas en Venezuela en un periodo de 30 años con alto tiempo de integración. Se usaron los modelos Rice-Holmberg, Moupfouma-Martin refinado y Recomendación ITU-R P.837-7, para estimar la precipitación, y los modelos Recomendación ITU-R P.618-13, Ramachandran-Kumar, Yeo-Lee-Ong y Rakshit-Adhikari-Maitra, para predecir la atenuación por lluvias para enlaces del satélite “Simón Bolívar” en Venezuela. Los resultados de ambos tipos de mapas representan una herramienta útil para la planificación preliminar de dichos enlaces en el país, específicamente en las bandas Ku y Ka.

Palabras-clave: mapas de contorno; predicción de tasa de precipitación; conversión de tiempos de integración; estimación de atenuación por lluvias; enlaces satelitales.

1. Introduction

In 2008, Venezuela launched its first satellite, “Simón Bolívar”, into space. It is a telecommunications satellite located in the geostationary orbit at 78o W longitude and operates at C, Ku and Ka-bands [1]; hence, it is necessary to consider the impact of rain events on satellite links performance, which deteriorates the carrier-to-noise ratio (C/N) in the uplink and downlink, to the point that the system may become unavailable in case of being undersized or involved in an investment higher than necessary, if the system is oversized [1].

Rain affects telecommunication systems operating in frequencies above 10 GHz. At these frequencies, raindrop sizes have similar magnitude that signal wavelength, and one part of the electromagnetic wave energy is absorbed and dispersed, resulting in attenuation at the received power [1]. If these effects are not properly estimated, it may result in undersized/oversized links.

Previous studies have shown several models rain attenuation prediction for satellite communication systems, for example, Recommendation ITU-R P.618-13 [2], Ramachandran-Kumar model [3], Yeo-Lee-Ong model [4] and Rakshit-Adhikari-Maitra model [5]. These models depend on rain rate, R, so accuracy in predicting rain attenuation is mainly related to whether or not rain measurements exist in the locality of interest. Here, it is important to emphasize that the rain rate cumulative distribution to be used in rain attenuation models should be constructed from measurements performed with of very short temporal resolution, i.e., 1-min integration time or less [1,6,7]. This means that in countries, like Venezuela, where such statistics are not available, it is necessary to implement conversion methodologies (e.g., Recommendation ITU-R P.837-7 Annex 2 [6]) or analytical formulations (e.g., Rice-Holmberg model [8], refined Moupfouma-Martin model [9] and Recommendation ITU-R P.837-7 Annex 1 [6]) to obtain cumulative distributions of precipitation rate suitable for satellite links planning, from statistics with longer integration times (e.g., hourly, daily and monthly), considering local climate information.

This study continues the research initiated by Rujano-Molina et al. [10]. It develops, for the first time in Venezuela, contour maps of rain rate for various exceedance probabilities of time, p%, based on 1-min rainfall rate statistics at several locations of the country. These were obtained through Rice-Holmberg model, refined Moupfouma-Martin model and Recommendation ITU-R P.837-7 Annex 1. In addition, rain attenuation contour maps are also generated by the widely recognized and accepted ITU model (Recommendation ITU-R P.618-13), as well as by some models developed for tropical climates, i.e., Ramachandran-Kumar model, Yeo-Lee-Ong model and Rakshit-Adhikari-Maitra model. Both maps were created and based on local rain rate and rain attenuation estimated by these models. The interpolation spatial method known as Spline [11] was used because it offers good results for less uniformly distributed data, such as those for Venezuela.

It is important to point out this work, with or without contour maps, follows other similar studies developed in other countries, for example, Colombia [12], Nigeria [13], South Africa [14], India [15], Bangladesh [16], Brazil [17], and Libya [18], among others.

2. Rainfall rate conversion models

In Rice-Holmberg (RH) model, data related with short-term precipitations was extrapolated to generate rain rate cumulative distributions based on measurements made at the meteorological stations. The probability for a rain rate, R, with a 1-min integration time is mainly based on average cumulative precipitation over an annual period, M, the highest monthly precipitation in the period, M m , and average number of days, U, with storm or convective rains within the same year [8]. The complementary cumulative distribution function (CCDF) of the rain rate, P(rR), is given by

It is worth highlighting that the RH model makes a good estimate of the rain rate distribution for percentages of link unavailability (p%), from 0.01% to 0.1% [1,12,19]. Additionally, the RH model allows to estimate rain rate cumulative distribution from measurements with high integration times (e.g. a week, month, etc.) [12,19].

Since the implementation of RH model requires the knowledge of data that is not always available from the meteorological services of a country, such as that related with convective rains, refined Moupfouma-Martin (RMM) model [9,14] can be implemented. This model is used for both tropical and subtropical regions and approximates low rain rates to a log-normal distribution, while high rainfall rate is approximated by gamma distribution. The CCDF of 1-min rain rate in RMM model is defined by [1,9,14]

where λ and ү are factors that depend on the climate zone. For example, for tropical and subtropical climates λ = 1.066 and ү = 0.214 [9].

Eq. (3)-(5), show that the rainfall rate probability, P(rR), in refined Moupfouma-Martin model also depends on rain rate, R 0.01 , exceeded at 0.01% of time in an average year. This can be estimated using the Chebil-Rahman model [20], based on a power law for rain rate calculation with a 1-min integration time by using only annual cumulative average precipitation M, in conjunction with regression parameters, α and β, as follows

where ( and ( are regression parameters equals to 12.2903 and 0.2973, respectively.

Finally, Recommendation ITU-R P.837-7 Annex 1 [6] is a global rain prediction method that considers the impact of monthly predictions of rainfall rate cumulative distribution. Input data involves the geographic coordinates (latitude, longitude) of the specific site, monthly rain rate in an annual period, time percentage of interest and monthly ground temperatures (referred 2 meters above the surface). If temperatures are not available, bi-linear interpolation from digital maps with long-term reliable data provided by [21] can be used. The Recommendation ITU-R P.837-7 Annex 1 is based on the following two assumptions [7]: a) the monthly rainfall rate statistics are approximated by log-normal distribution, whose scale parameter, ( i , does not depend on the locality and is equal to 1.26; b) the monthly mean rain rate, r i , is a function of the monthly mean ground temperature, T i . The annual rainfall rate probability of exceedance is obtained by [6,7]

where N i represents the number of days in each month.

In Eq. (7), P i (rR) corresponds to the monthly CCDF of rain rate for the ith month, defined as

where Q and erfc are the well-known Q-function and complementary error function, respectively.

In Eq. (8), represents the monthly probability of rain, given by

where MT i is monthly mean total rainfall. If is greater than 70, then = 70 and r i is calculated by

Eq. (10), shows that, if reliable long-term local data is not available, MT i at the desired location can be determined performing a bi-linear interpolation from digital maps of mean monthly value of rainfall, provided as an integral part of Recommendation ITU-R P.837-7.

3. Rain attenuation prediction models in satellite links

Recommendation ITU-R P.618-13 [2] shows a method to predict rain attenuation in satellite links considering some of the following input data: operating frequency, f, effective radius of the Earth, R E , mean elevation above sea level of the earth station, h s , earth station elevation angle, θ, earth station latitude, (, and rain rate to 0.01% of the time in an average year, R 0.01 . The method estimates rain for any percentage of time (p%) and consists of ten (10) steps that include, among others, calculation of h R (i.e., average annual rainfall height above the mean sea level), inclined path, L s (depends on h s , h R , R E and θ) and horizontal projection of the inclined path, L G (depends on L s and θ). The algorithm of the Recommendation ITUR P.618-13 also includes specific rain attenuation calculation for 0.01% of the time, ү 0.01 , horizontal reduction factor determined for 0.01% of the time, r 0.01 and adjustment factor in the vertical direction for 0.01% of the time,v 0.01 . Finally, rain attenuation, A p , exceeded during percentages of time, p(%), other than 0.01%, from 0.001% to 5% range, is calculated through interpolation in function of p, A 0.01 (attenuation exceeded for 0.01% of the time), earth station elevation angle, θ, and earth station latitude, (.

Ramachandran-Kumar (RK) model was developed on rain attenuation measurements performed on Ku-band satellite in Fiji, a tropical region [3]. The model considers the existence of a breakpoint, related to a change in rain structure, in the decrease of complementary cumulative distribution (CCD) of precipitation rate and shows two adjusted expressions of Recommendation ITU-R P.618-8 (published in 2003): one for lower rainfall rates at the breakpoint and another one for higher rainfall rates. According to the model, rain rate breakpoint takes place in a exceedance probability of time, p R , approximately equal to 0.01%, while for attenuation, the percentage of time, p A , in which the breakpoint occurs, is 2.14 times higher than p R . Therefore, the attenuation, A 0.01 , of original ITU-R model, becomes A B , that is the attenuation exceeded for 0.021% of time, calculated by the original version. In general, rainfall attenuation given by RK model includes an expression for (0.021 ≤ p <1) that incorporates a correction factor, C f , for link elevation angle, θ, between 40o and 60o (consider crossed slant path with more than one cell), and another expression used to estimate rain attenuation for p ≤ 0.021.

Yeo-Lee-Ong (YLO) model [4] was developed on measurements carried out in six tropical countries (Thailand, Indonesia, Malaysia, Singapore, Papua New Guinea and Brazil). The method is similar to that of Recommendation ITU-R P.618-11 (published in 2013), with the following variants: a) use of a single path adjustment factor instead of the horizontal and vertical adjustment factors; this factor depends on θ, R 0.01 , f, h s and rain height H = h R ; b) the regression parameters, including factor (, from the exponent in the interpolation expression to calculate the attenuation in original ITU-R model, are readjusted on the measurements of the six tropical regions mentioned above.

Finally, Rakshit, Adhikari and Maitra (RAM) [5] developed a model for rainfall decay parameter Г, of Simple Attenuation Model (SAM) [22], using rain attenuation measurements carried out on Kuband satellite in Calcutta, India, a tropical climate region. By applying the numerical method known as bisection, the model offers values of the aforementioned parameter depending on point rain rate, R p , for the percentage of the time, p, in span window of 2 mm/h. In the model, rain decay parameter decreases to 25 mm/h, does not change significantly between 25 mm/h and 40 mm/h, and increases at the 40 mm/h interval.

The procedure for calculating the rainfall attenuation for a 1-minute integration time in each of the aforementioned models is shown in Appendix.

4. Rain rate contour maps in Venezuela using local climatological data

Venezuela is located in the inter-tropical zone from 0o to 13o latitude North and 59o to 54o longitude West. This country has an interesting variety of climates, namely, tropical rainforest, tropical savanna, tundra, desert, steppes, and so on. In this work, monthly precipitation rate measurements made along 30 years in 35 of the Venezuela’s meteorological stations is used for generating rain rate contour maps of the country to generate rain rate contour maps of the country. The measurements were supplied by the Instituto Nacional de Meteorología e Hidrología (INAMEH) [23] and initially explored by Rujano-Molina et al. [10]. Table 1 summarizes the information about the stations.

Table 2 shows rain rates converted to 1-min integration time for exceedance probability 0.01% (input parameter considered in various rain attenuation models) and 0.5% (a typical percentage in planning of satellite communication operating at frequencies above 10 GHz), using RH and RMM models and Recommendation ITU-R P.837-7 Annex 1.

Contour maps were constructed using the Spline method in the ArcGIS software [24] to interpolate in space the rain rate obtained through the aforementioned models. Figs. 1 and 2 show the results obtained when the RMM model was applied for both percentages of link unavailability, 0.01% and 0.5%, respectively.

As shown, for both exceedance probabilities the highest rain rates occur at southern sector (Amazon state), northeastern zone (Monagas state), western area (Táchira state, Córdoba municipality) and the most northwestern region of the country (Zulia state, Catatumbo municipality). These rates are expected because these areas correspond to regions of the country with tropical rainforest climate (i.e., Amazon), tropical savanna climate (i.e., Monagas and Tachira) and a high number of thunderstorms in an average year (i.e., Catatumbo). The lowest rain rates, on the other hand, occur at northwestern region (Falcón state) and Margarita Island, having both regions tropical desert climate. Figs. 1 and 2 and Table 1, also show that for the percentage of time equal to 0.5% the values of the precipitation rate obtained are lower than the values of such rate for 0.01% of the time. This is consistent with the well-known fact that high rainfall intensities correspond to short-duration time intervals, i.e., lower percentages of time for which rain is exceeded [1].

Source: The Authors.

Figure 1 Rain rate contour map using RMM model for p = 0.01%. 

Source: The Authors.

Figure 2 Rain rate contour map using RMM model for p = 0.5%. 

Table 1 Precipitation Data in Venezuela. 

Source: The Authors.

Table 2 Rainfall Rate for 0.01% and 0.5% of time. 

Source: The Authors.

A similar result regarding the geographic behavior of rain rate in the country is observed in Figs. 3 and 4, if the RH model and Recommendation ITU-R P.837-7 are used, respectively, for 0.01% of time. However, the values varying between 40 mm/h and 136 mm/h (for RH model) and from 42 mm/h to 96 mm/h (for ITU-R P.837-7), are lower than those obtained with the RMM model for same percentage of time. Although the lack of measurements in Venezuela of rain rate for 1-min integration time not allow to indicate which rainfall conversion model implemented in this work shows the best performance in the estimation, the results obtained are quite valuable not only because they may be helpful when carrying out a future evaluation once when measurements could be available, but also because these results represent a very important input for the rain attenuation preliminary prediction in the country, for satellite and terrestrial links.

Source: The Authors.

Figure 3 Contour map for rain rate for p = 0.01%, using RH model. 

Source: The Authors.

Figure 4 Contour map for rain rate for p = 0.01%, using Recommendation ITU-R P.837-7. 

Source: The Authors.

Figure 5 Contour map for rain attenuation for p = 0.5%, using Recommendation ITU-R P.618-13 model, in Ku band. 

Source: The Authors.

Figure 6 Contour map for rain attenuation for p = 0.5%, using Recommendation ITU-R P.618-13 model, in Ka band. 

5. Rainfall attenuation contour maps in Venezuela

Using rain rate values estimated in section 4 and implementing Recommendation ITU-R P.618-13, Ramachandran-Kumar (RK), Yeo-Lee-Ong (YLO) and Rakshit-Adhikari-Maitra (RAM) models, rain attenuation contour maps in Venezuela werw obtained. Attenuation was calculated for all meteorological stations of the Table 1, considering the downlink in Ku and Ka-bands (at 12 GHz and 20 GHz, respectively) with horizontal polarization.

As an example, based on rain rate results from Recommendation ITU-R P.837-7, Figs. 5 and 6 illustrate rain attenuation contour maps for Ku and Ka bands, respectively, when ITU-R P.618-13 model to attenuation prediction, with exceedance probability of 0.5%, is used. It can be seen that, as expected, rain attenuation is greater, in all locations, for the Kaband. Furthermore, it is observed that the highest attenuations occur at northwestern zone (Zulia state), western region (Táchira state) and southern regions (Amazonas and Apure states), while the lowest values occur at western center zone (Lara state), Andes mountains area (Mérida state), northwestern region (Falcón state) and center region (Distrito Capital). Here it is important to emphasize that attenuation also depends on mean elevation above sea level of the earth station, elevation angle, among other factors.

Figs. 7 to 10 show rain attenuation contours maps for 0.01% of time, Ku band, using the four models considered in the study. The difference found in rain attenuation values expected for the country is observed. It is observed a difference in the rainfall attenuation values estimated for the country, by depending on the model used for the prediction. Again, it is important to point out that although the lack of measurements, this time, of rain attenuation in satellite links in Venezuela, it is not possible to determine which of the rainfall fade prediction models has the best performance, the results obtained are the preamble for the future real evaluation in the country of those models, or even other models not analyzed in this study, in order to reduce undersizing or oversizing of the system.

Source: The Authors.

Figure 7 Rain attenuation contour map for p = 0.01%, Ku band, using Recommendation ITU-R P.618-13. 

Source: The Authors.

Figure 8 Rain attenuation contour map for p = 0.01%, Ku band, using Ramachandran-Kumar model. 

Source: The Authors.

Figure 9 Rain attenuation contour map for p = 0.01%, Ku band, using Yeo-Lee-Ong model. 

Source: The Authors.

Figure 10 Rain attenuation contour map for p = 0.01%, Ku band, using Recommendation ITU-R P.618-13. 

Finally, by comparing Figs. 5 and 7, an increase in rain attenuation could be observed when the percentage of time unavailability decreases. This is another expected aspect that should be considered by system designers.

6. Conclusions

In this paper, contour maps of 1-min rain rates and rain attenuation for exceedance probability of 0.01% and 0.5% have been shown in detail for the first time in Venezuela. This was achieved by using the Rice-Holmberg model, refined Moupfouma and Martin model and Recommendation ITU-R P.837-7 Annex 1 in order to convert monthly rain rate statistics in the country to rainfall rates with 1-min integration time and using the Recommendation ITU-R P.618-13, Ramachandran-Kumar, Yeo-Lee-Ong and Rakshit-Adhikari-Maitra models to estimate rain attenuation. The maps generated may are useful for preliminary “Simon Bolivar” satellite links planning in the country, because they may provide information quickly about the rain rate and rain attenuation levels.

Nevertheless, measurements should be carried out in Venezuela in the future to evaluate performance of the models analyzed in this study, other models not included here, and also to develop rain rate prediction models and rain attenuation estimation models, adjusted to the climate characteristics of the country.

Acknowledgement

We are grateful to Instituto Nacional de Meteorología e Hidrología (INAMEH) from Venezuela for providing us with the rainfall data used in this research.

References

[1] Pérez-García, N.P., et al., Fenómenos atmosféricos que afectan a las comunicaciones vía satélite, in Perez-Garcia, N.A., Pinto-Mangones, A.D., Torres-Tovio, J.M. and Ramirez, E.J., Planificación y dimensionamiento de sistemas de comunicación vía satélite , Montería, Colombia Sello Editorial Corporación Universidad del Sinú, 2018, pp. 36-71. ISBN: 978-958-8533-52-8 [ Links ]

[2] International Telecommunication Union - Radiocommunication, Recommendation P.618-13: propagation data and prediction methods required for the design of earth-space telecommunication systems, Geneva, Switzerland, 2017. [ Links ]

[3] Ramachandran, V. and Kumar, V., Modified rain attenuation model for tropical regions for Ku-band signals. Int. J. Satell. Commun. Netw., 25(1), pp. 53-67, 2007. DOI: 10.1002/sat.846 [ Links ]

[4] Yeo, J.X., Lee, Y.H. and Ong, J.T., Rain attenuation prediction model for satellite communications in tropical regions. IEEE Trans. Antennas Propag., 62(11), pp. 5775-5781, 2014. DOI: 10.1109/TAP.2014.2356208 [ Links ]

[5] Rakshit, G., Adhikari, A. and Maitra, A., Modelling of rain decay parameter for attenuation estimation at a tropical location. Adv. Space Res., 59(12), pp. 2901-2908, 2017. DOI: 10.1016/j.asr.2017.03.028 [ Links ]

[6] International Telecommunication Union - Radiocommunication, Recommendation P.837-7: Characteristics of precipitation for propagation modelling, Geneva, Switzerland, 2017. [ Links ]

[7] Luini, L., Emiliani, L., Boulanger, X., Riva, C. and Jeannin N ., Rainfall rate prediction for propagation applications: Model performance at regional level over Ireland. IEEE Trans. Antennas Propag., 65(11), pp. 6185-6189, 2017. DOI: 10.1109/TAP.2017.2754448 [ Links ]

[8] Rice, P.L. and Holmberg, N.R., Cumulative time statistics of surface-point rainfall rates. IEEE Trans. Commun., 21(10), pp. 1131-1136, 1973.10. . DOI: 10.1109/TCOM.1973.1091546 [ Links ]

[9] Moupfouma, F. and Martin L. , Modelling of the rainfall rate cumulative distribution for the design of satellite and terrestrial communication systems. Int. J. Satell. Commun., 13(2), pp. 105-115, 1995. DOI: 10.1002/sat.4600130203 [ Links ]

[10] Rujano-Molina, L.M., Perez-Garcia, N.A. y Nariño-González T ., Distribuciones acumulativas de la tasa de lluvias con tiempo de integración de 1-minuto en Venezuela. Ingeniería Al Día [Online], 3(1), pp. 24-44, 2017 [date of reference: April 10th, 2018]. Available at: Available at: http://revista.unisinu.edu.co/revista/index.php/Links ]

[11] Vieux, B.E., Distributed hydrologic modeling using GIS, 3rd ed. Netherlands: Springer, 2016, pp. 54-61. ISBN: 978-94-024-0930-7 [ Links ]

[12] Emiliani, L.D., Agudelo, J., Gutierrez, E. and Fradique-Mendez C ., Development of rain-attenuation and rain-rate maps for satellite system design in the Ku and Ka bands in Colombia. IEEE Antennas Propag. Mag., 46(6), pp. 54-68, 2004. DOI: 10.1109/MAP.2004.1396736 [ Links ]

[13] Ojo, J.S., Ajewole, M.O. and Emiliani, L.D., One-minute rain-rate contour maps for microwave-communication-system planning in a tropical country: Nigeria. IEEE Antennas Propag. Mag., 51(5), pp. 82-88, 2009. DOI: 10.1109/MAP.2009.5432046 [ Links ]

[14] Ojo, J.S. and Owolawi, P.A., Development of one-minute rain-rate and rain-attenuation contour maps for satellite propagation system planning in a subtropical country: South Africa. Adv. Space Res., 54(8), pp. 1487-1501, 2014. DOI: 10.1016/j.asr.2014.06.028 [ Links ]

[15] Kestwal, M.C., Joshi, S. and Garia L.S ., Prediction of rain attenuation and impact of rain in wave propagation at microwave frequency for tropical region (Uttarakhand, India). Int. J. Microw. Sci. Technol., 2014, pp. 1-6. DOI: 10.1155/2014/958498 [ Links ]

[16] Imran, A.Z.M., Islam, M.T., Gafur, A. and Rabby, Y.W., Rain attenuation prediction analysis and contour map design over Bangladesh, Proceedings. of the 18th International Conference on Computer and Information Technology (ICCIT 2015), Dhaka, Bangladesh, pp. 208-212, 2015. DOI: 10.1109/ICCITechn.2015.7488069 [ Links ]

[17] Fialho, M.R.B., Silva, R.M.L., Pontes, M.S., Oliveira, C.H.R. and Gonsioroski, L.H., Mapas de contorno da taxa de chuva utilizando os históricos de dados acumulados de precipitação no Estado do Maranhão para modelagem de fenômenos radiometeorológicos, Anais do XXXIV Simpósio Brasileiro de Telecomunicações e Processamento de Sinais (SBrT 2016), Santarém, PA, Brasil, 2016, pp. 148-151. [ Links ]

[18] Islam, M.R., Budalal, A.A.H., Habaebi, M.H., Badron, K. and Ismail, A.F., Performance analysis of rain attenuation on earth-to-satellite microwave links design in Libya, Proceedings of the 6th International Conference on Mechatronics (ICOM 2017), Kuala Lumpur, Malaysia, 2017, pp. 1-6. [ Links ]

[19] Mandeep, J.S. and Hassan, S.I.S., Performances of existing rain rate models in equatorial region. J. Geophys. Res., 113(D11), pp. 1-7, 2008. DOI: 10.1029/2007JD009737 [ Links ]

[20] Chebil, J. and Rahman, T.A., Development of 1 min rain rate contour maps for microwave applications in Malaysia peninsula. Electron. Lett., 35(20), pp. 1172-1174, 1999. DOI: 10.1049/el:19991188 [ Links ]

[21] International Telecommunication Union - Radiocommunication Sector (ITU-R), P Series, Radiowave propagation, Recommendation P.1510-1: mean surface temperature, Geneva, Switzerland, 2017. [ Links ]

[22] Stutzman, W.L. and Yon, K.M., A simple rain attenuation model for earth-space radio links operating at 10-35 GHz. Radio Sci., 21(1), pp. 65-72, 1986. DOI: 10.1029/RS021i001p00065 [ Links ]

[23] Instituto Nacional de Meteorología e Hidrología [Online]. [date of reference: March 31th, 2018]. Available at: Available at: http://www.inameh.gob.ve/web/Links ]

[24] ArcGIS [Online]. [date of reference: June 23th, 2018]. Available at: Available at: https://www.esri.com/es-es/arcgis/about-arcgis/overviewLinks ]

[25] International Telecommunication Union - Radiocommunication, Recommendation P.839-4: rain height model for prediction methods, Geneva, Switzerland, 2013. [ Links ]

[26] International Telecommunication Union - Radiocommunication, Recommendation P.838-3: specific attenuation model for rain for use in prediction methods, Geneva, Switzerland, 2005. [ Links ]

How to cite: Pinto-Mangones, D, Pérez-García, N.A, Torres-Tovio, J.M, Ramírez E.J, Castaño-Rivera S.0, Velez-Zapata, J, Ferreira-Rodríguez, J.D. and Rujano-Molina L.M. Rainfall rate and rain attenuation contour maps for preliminary “Simon Bolivar” satellite links planning in Venezuela. DYNA, 86(209), pp. 30-39, April - June, 2019.

A.D. Pinto-Mangones, is a BSc. in Systems Eng. from the Universidad del Sinú, Monteria, Colombia, in 2004. He received the MSc degree in Telematics and PhD. degree in Management Science and Technology, both from URBE, Maracaibo, Venezuela, in 2015 and 2016, respectively. Since 2009, Dr Pinto-Mangones has been a full professor at the Departamento de Ingenieria de Sistemas, Universidad del Sinu, Colombia. ORCID: 0000-0002-6792-8095

N.A. Pérez-García, is a BSc. in Electrical Eng. from the Universidad de Los Andes (ULA), Mérida, Venezuela, in 1992. He obtained MSc. and PhD. degrees in Electrical Engineering (Applied Electromagnetism area), from the Pontificia Universidade Catolica do Rio de Janeiro (PUC/Rio), Brasil, in 2000 and 2003, respectively. Since 1994, Dr. Pérez-García he is a full professor of ULA, Mérida, Venezuela. ORCID: 0000-0002-6163-933X

J.M. Torres-Tovio, is a BSc. in Informatics Eng. from the Universidad Católica del Norte, in Medellin, Colombia, in 2005. He received the PhD degree in Management of Science and Technology from URBE, Maracaibo, Venezuela, in 2015. Since 2015. Dr. Torres-Tovio has been a full professor at the Universidad del Sinú, Monteria, Colombia. ORCID: 0000-0002-9368-1436

E.J. Ramírez, is an Electrical Engineer from the Universidad de Los Andes (ULA), Merida, Venezuela, in 2003. Currently, he is a candidate for MSc (Telecommunications) and PhD (Applied Science) degrees at ULA. ORCID: 0000-0002-3728-4196

S.O. Castaño-Rivera, is a Systems Engineer and holds a MSc degree in Open Software from Universidad Autónoma de Bucaramanga, Colombia, in 2009. He is currently a professor in Universidad de Córdoba, Montería, Colombia. ORCID: 0000-0003-3302-3439

J. Vélez-Zapata, is an Electronic Engineer from Universidad de Antioquia, Colombia, and obtained PhD degree in Management of Science and Technology, at URBE, Maracaibo, Venezuela, in 2017. Currently, he is professor at Corporación Universidad de la Costa, Barranquilla, Colombia. ORCID: 0000-0002-0704-2292

J.D. Ferreira-Rodríguez, is a BSc in Electrical Eng. from the Universidad de Los Andes (ULA), Merida, Venezuela, in 1992, and PhD. from the University of Manchester, United Kingdom, in 2003. Since 1993, Dr. Ferreira-Rodríguez has been a full professor and a researcher at ULA. He is currently the research leader of the Institutional Projects called “Tren Electromagnético” (TELMAG) and Cohete Sonda, at ULA. ORCID: 0000-0001-6195-3516

L.M. Rujano-Molina, is a BSc. in Telecommunications Eng. from the UNEFA, Merida, Venezuela, in 2012. Obtained Lawyer degree and MSc degree in Telecommunications, from the Universidad de Los Andes (ULA), Merida, Venezuela, in 2016 and 2017, respectively. ORCID: 0000-0002-7015-6559

7. Appendix: rain attenuation calculation according to ITU-R P.618-13, RK, YLO and RAM models

The rain attenuation estimation in satellite communication systems is based on the link geometry illustrated in Fig. 11, whose parameters have been defined in section 3.

7.1. Recommendation ITU-R P.618-13

Step 1

Source: The Authors.

Figure 11 Schematic diagram of the Earth-space path. 

where h 0 is the average height of the 0 oC isotherm. This height refers to the earth station site geographical coordinates and is given in Recommendation ITU-R P.839-4 [25].

Step 2

Step 3

Step 4

Obtain the rain rate R 0.01 , exceeded during 0.01% of time in an average year, for a 1-min integration time. This parameter should better be obtained from measurements of the rainfall rate corresponding to the location of the earth station. However, in the absence of such an information, the rain rate Recommendation ITU-R P.837-7 [6] database can be used.

Step 5

where k and α are the regression coefficients that depend on the operation frequency and wave polarization and are given by Recommendation ITU-R P.838-3 [26].

Step 6

Step 7

Step 8

Calculate the effective path length L E , as follows

Step 9

Step 10

7.2. Ramachandran-Kumar model

The first 4 steps of RK model are identical to those of Recommendation ITU-R P.618-13. Likewise, steps 5 to 7 also unlike that of model ITU-R, with the exception that ( 0.01 , r 0.01 and ( 0.01 , change to ( B , r B and ( B , where the subscript B refers to the breakpoint explained in section 3.

Step 8

Step 9

Step 10

Obtain estimated attenuation, in dB, to be exceeded for p(%) of time an average year.

For p ≤ 0.021

For 0.021 ≤ p < 1

where the parameter ( is given by Eq. (24).

7.3. Yeo-Lee-Ong model

The first 5 steps of the YLO model are the same as Recommendation ITU-R P.618-13.

Step 6

Step 7

Step 8

7.4. Rakshit-Adhikari-Maitra model

Step 1

In RAM model, rain height above the mean sea level is calculated as follows

where H 0 is the same elevation above sea level of the earth station, h s .

Step 2

In this case, L s is given only by

Step 3

Step 4

Calculate the attenuation to be exceeded for p(%) of time an average year.

For R p ≤ 10 m/h

For R p ≥ 10 m/h

In Eq. (38), the rainfall decay parameter (, depends on rainfall rate R p and is obtained from [5].

Recebido: 28 de Julho de 2018; Revisado: 04 de Fevereiro de 2019; Aceito: 25 de Fevereiro de 2019

Creative Commons License This is an open-access article distributed under the terms of the Creative Commons Attribution License