The objective of the research consisted in demonstrating that the long-term meteorological prediction for December 10 is possible by means of the methodology of the Regressive Objective Regression (ROR) with an acceptable margin of error for eight meteorological variables in the Yabu station, Santa Clara, Villa Clara, Cuba. The daily data base of this meteorological station was used, where the following variables were taken into account: extreme temperatures, extreme humidity and its mean value, precipitation, wind force and cloudiness corresponding to the period 1977-2021. Each variable, which explains its behavior according to these variables, was obtained 15, 13, 10 and 8 years in advance. The calculation of the mean error with respect to the persistence forecast in temperatures, wind strength and cloudiness, as well as the persistence model was better with respect to humidity, which makes it possible to have valuable long-term information on the climate in a given locality, which translates into better decision-making in the different aspects of the economy and society, impacted by climate variability. It is concluded that these models allow long-term weather forecasting, opening a new possibility of forecasting, so that meteorological chaos can be overcome if this form of forecasting is used; besides being the first time that an ROR model is applied to weather forecasting processes for a specific day so many years in advance.
In different fields of study, cumulative growth models over time play an important role, and many researchers have contributed to the knowledge of the developed models [1-2]. Among the most common nonlinear models are: Gompertz, Weibull, negative exponential, Richard's model, logistic, monomolecular, Brody, Mitcherlich, von Betalanffy, S-shaped, among others [2-3]. There are about 77 known equations referring to sigmoidal growth models, which are used in epidemics, bioassays, agriculture, engineering, tree diameter, forest height distribution, etc. [1,2,4,5]. Among the most widely used models are those of Gompertz, Richard and Weibull [2,3]. These models have been implemented in COVID-19, a new disease; where it was shown for COVID-19 cases in Iraq, that ROR modeling gave better results than nonlinear models, so for meteorology we will use this linear ROR model [6-7].
In this study we developed a linear ROR model of eight temporal variables, these were: maximum and minimum extreme temperatures, daily precipitation, maximum and minimum humidity, daily mean, wind strength and cloudiness, which we will compare with the persistence model, which is a good model obtained for meteorological prediction, and we will show how it is possible to predict it on a given day well in advance which shows that chaos theory can be defeated if this way of prediction is used [7-8].
In the case of this research we have data from the Yabu weather station (located in the city of Santa Clara, Villa Clara province, Cuba) for that day over the years, and it is with this data that we demonstrate that long-term weather prediction for this particular day (December 10) is also possible with an acceptable margin of error, which shows that chaos can be overcome if this form of prediction is used.
The objective of the research was to demonstrate that the long-term meteorological prediction for December 10 is possible by means of the methodology of the Regressive Objective Regression (ROR) with an acceptable margin of error for eight meteorological variables in the Yabu station, Santa Clara, Villa Clara, Cuba, thus overcoming the Chaos theory.
For this research, data from the Yabu weather station (located in the city of Santa Clara, Villa Clara province, Cuba) were used (Figure 1), for this day (December 10) over the years, in order to demonstrate that long-term weather forecasting for this particular day is also possible with an acceptable margin of error.
The data used were eight meteorological variables (extreme maximum and minimum temperatures, daily precipitation, maximum and minimum humidity, average daily humidity, wind strength and cloudiness). The forecast was made using the methodology of Regressive Objective Regression (ROR) that has been implemented in different variables, such as viruses circulating in the province of Villa Clara, the population dynamics of culicidae, the entomofauna of culicidae and copepods [9-11], and particularly in COVID-19 in Cuba [6,12]. A long-term forecast of up to 15 years was made, from December 10, 2020 to December 10, 2035.
In the linear ROR methodology, first the dichotomous variables DS, DI and NoC must be created, where: NoC: Number of base cases (its coefficient in the model represents the trend of the series). DS = 1, if NoC is odd; DI = 0, if NoC is even, and vice versa. DS represents a saw tooth function and DI this same function, but in inverted form, so that the variable to be modeled is trapped between these parameters and a large amount of variance is explained, a detailed explanation of which can be found in Fimia et al. [10].
General Description of Climate
According to Köppen's classification (modified), in most of Cuba the predominant climate is of the warm tropical type, with a rainy season in the summer. In general, it is quite accepted to express that Cuba's climate is tropical, seasonally humid, with maritime influence and semi-continental features, where Villa Clara province also responds to these characteristics.

Figure 1: Geographical Distribution of Meteorological Stations
Due to its geographical position, the city of Santa Clara is located at a latitude very close to the Tropic of Cancer, which conditions the reception of high values of solar radiation throughout the year, determining the warm character of its climate. In addition, it is located in the limit between the tropical and extratropical circulation zones, receiving the influence of both in a seasonal way. In the season from approximately November to April, the meteorological and climatic variations become more noticeable, with abrupt changes in the daily weather, associated with the passage of frontal systems, the anticyclonic influence of continental origin and the extratropical low pressure centers. From May to October, on the other hand, there are few variations in the weather, with the more or less marked influence of the North Atlantic anticyclone. The most important changes are related to the presence of disturbances in the tropical circulation (easterly waves and tropical cyclones).
The analysis and processing of the data was carried out with the help of the IBM SPSS statistical package, version 19.
Table 1 shows the model parameters, where the coefficients of the explained variance R are high, while the model errors were small. Temperatures depend on lags in 15 years the minimum, and 13 years the maximum, precipitation in eight years, minimum humidity in ten years, wind force (FF) in 13 years and cloudiness in 15 years.
Table 1: Parameters of the Models of the Different Time Variables in Yabu, Cuba
| Variables | R | Error of the model | Saw Tooth (DS) | Inverted Saw Tooth (DI) | Tend | Delay Coeficient | Fisher F |
T Min | 98.8 | 3.09 | 10.73 | 8.86 |
| 15 0.472 | 335 |
T max | 99.7 | 2.41 | 51.08 | 49.93 |
| 13 -0.833 | 1350 |
Precip | 100 | 0.0 | -8 E(-16) | -1.5E(-15) | 2.98E(-17) | 8 -5.8E(-17) | 7.51i |
HRx | 100 | 1.92 | 96.44 | 96.11 | 0.041 |
| 36861 |
Hrn | 98.7 | 9.80 | 77.99 | 73.77 | -0.054 | 10 -0.333 | 277 |
Hrm | 99.7 | 6.16 | 81.57 | 79.62 | 0.070 |
| 2490 |
FF | 90.2 | 2.91 | 2.353 | 2.586 | 0.049 | 13 0.261 | 27.2 |
Cloudiness | 92.4 | 1.55 | 2.98 | 2.634 |
| 15 0.496 | 32.2 |
The analysis of variance of the model was highly significant, with Fisher's F of 1350.6 significant at 100 % (Table 2).
Table 2: Results of the Analysis of Variance
Model | Sum os Squares | Df | Cuadratic mean | F | Sig. | |
1 | Regression | 23618.784 | 3 | 7872.928 | 1350.618 | 0.000c |
Residuals | 163.216 | 28 | 5.829 |
|
| |
Total | 23782.000d | 31 |
|
|
| |
ANOVA a, b, a. Dependent variable: Tmax, b. Linear regression through the origin, c. Predictors: Lag13Tmax, DS, DI, d. This total sum of squares is not corrected for the constant because the constant is zero for the regression through the origin.
In table 3, it can be seen that all variables are significant, but the temperature trend was not significant, and therefore is not shown. The linear model depends on the data returned in 13 years (Lag13 Tmax).
Table 3: Results of the Linear Model for the Variables Included in the Study
| Model | Unstandarizaded Coefficient | Standardized Coefficient | t | Sig. | ||
B | Standard error | Beta | ||||
1 | DS | 51.079 | 6.010 | 1.283 | 8.499 | 0.000 |
DI | 49.939 | 5.930 | 1.295 | 8.422 | 0.000 | |
Lag13Tmax | -0.833 | 0.216 | -0.831 | -3.864 | 0.001 | |
Coefficient a, b, a. Dependent variable: Tmax, b. Linear regression through the origin

Figure 2: Maximum Temperature Forecast According to ROR Linear Model for the Yabu, Cuba Meteorological Station
In figure 2, the modeling results show that the predicted values coincide with the actual values of the variable, where the errors are very small, close to zero. Since no trend is observed in the extreme temperatures, we can state that on December 10 throughout history there is no cooling or warming trend, which does not mean that for any other specific day in the daily data series there may be a significant trend.
In the 13 year forecast, it was obtained that December 10 will present days with high and low maximum temperatures. This is the first time that an ROR model is applied to the long-term weather forecast in Cuba 13 years in advance for a specific day, as it had already been done for climatic variables such as sand dumping on Varadero beach [13], also for electricity consumption in Villa Clara province [14]; in addition, other works also model the number of people in the long term with cerebrovascular accidents [7,12] as well as the modeling of atmospheric pressure and its impact on mosquito population density well in advance [10,11]. Another long-term model was carried out in Sancti Spiritus, Cuba, where a forecast of daily variables is made one year in advance using ROR regression [14], as well as in the prediction of earthquakes [1516,], and of different meteorological disturbances, also of epidemic outbreaks in different geographical latitudes, where global warming of the planet has a high share of responsibility, and therefore, modeling is an essential tool, and even human cloning [16-23].
Table 4 shows the results of the errors taking an independent sample of 11 cases, from 2010 to 2021 to see how the model behaves in an independent sample of 25 % of the cases. As can be seen, the errors are small for all variables, as are the standard deviations of the variables.
Table 4: Descriptive Statistics
| Parameters | N | Minimum | Maximum | Media | Standard desviation |
eTmin.mi | 11 | -7.82 | 3.89 | -1.1261 | 4.12549 |
eTmax.mi | 11 | -5.10 | 2.70 | -.3086 | 3.03874 |
er24h.mi | 11 | .00 | .00 | .0000 | .00000 |
eFFmed.mi | 11 | -4.23 | 4.83 | .1066 | 2.84084 |
eHrmax.mi | 11 | -4.65 | 2.51 | -.1206 | 1.92130 |
eHrmin.mi | 11 | -8.74 | 12.10 | .5386 | 7.09060 |
eHrmed.mi | 11 | -8.28 | 6.97 | .3082 | 5.00155 |
eNmed.mi | 10 | -1.87 | 1.81 | -.1527 | 1.34801 |
The forecast results for the remaining meteorological variables in Yabu, Cuba are shown below (Figures 3-9).

Figure 3: Extreme Minimum Temperature

Figure 4: Daily Precipitation

Figure 5: Daily Wind Force

Figure 6: Maximum Daily Relative Humidity

Figure 7: Minimum Relative Humidity

Figure 8: Mean Relative Humidity

Figure 9: Daily Cloudiness on December 10 in Yabu, Cuba
Finally, the rate of improvement of one model over another was calculated as follows, which is nothing else than the modeling ABILITY. Here a slight modification was made to the Wilks (1987) formula, adding the mean instead of the mean square error (MSE) according to the original formula, obtaining that the model established is the persistence model and the model to be tested is the ROR, this was done in the independent sample of 11 cases (Table 5).
Table 5: Results of Mean Error in an Independent Sample of 11 Cases
| Variable | Error ROR | Error Persistency | Improvment (%) |
Temperature Minimum (T Min) | -1.1261 | -1.3175 | 14.5 |
Temperature Máximum (T max) | -0.3086 | -0.7757 | 60.2 |
Precipitation | 0.000 | 0.000 | - |
Maximum Relative Humidity (HRx) | -0.1206 | -0.1105 | * |
Mínimum Relative Humidity (Hrn) | 0.5386 | 0.5250 | * |
Mean Relative Humidity (Hrm) | 0.3082 | 0.2002 | * |
Mean velocity of wind (FF) | 0.1066 | 0.1786 | 40.31 |
Cloudiness | -0.1527 | -0.3736 | 59.12 |
The improvement ranged from 14.5 % for minimum temperature to 60.2 % for maximum temperature, so we can state that in this case a linear ROR model outperforms a persistence model for most of the variables studied, except for humidities where persistence outperformed the ROR model by a narrow margin.
We can affirm in this case that the linear ROR model outperformed the persistence model in most of the variables studied, except in the case of humidity, where persistence outperformed the ROR model by a narrow margin, thus overcoming the chaos in weather forecasting, all of which is very positive, since it is the first time that a ROR model has been applied to weather forecasting processes for a specific day 8, 10, 13 and 15 years in advance.
Acknowledgement
We would like to thank the Villa Clara Cuba Provincial Meteorological Center for providing the computing resources without which this work would not have been possible.
De Bacco, C. et al. “A physical model for efficient ranking in networks.” Science Advances, vol. 4, 2018, eaar8260.
Dagogo, J., Nduka, E.C. and Ogoke, U.P. “Comparative analysis of Richards', Gompertz and Weibull models.” IOSR Journal of Mathematics, vol. 16, 2020, pp. 15–20.
Ban, G.A. “Statistical modeling of the novel COVID-19 epidemic in Iraq.” Epidemiological Methods, vol. 10, 2021, 20200025.
Wang, B., Chen, G. and Liu, F. “Diversity of the Madden–Julian oscillation.” Science Advances, vol. 5, 2019, eaax0220.
Yan, Z. et al. “Eastward shift and extension of ENSO-induced tropical precipitation anomalies under global warming.” Science Advances, vol. 6, 2019, eaax4177.
Osés, L.C. et al. “Comparison of lineal ROR vs nonlinear Weibull model for COVID-19 in Iraq.” Himalayan Journal of Applied Medical Sciences and Research, vol. 2, 2021, pp. 88–96.
Osés, R.R. et al. “Age prediction for COVID-19 suspects and contacts in Villa Clara Province, Cuba.” EC Veterinary Science, vol. 6, 2021, pp. 41–51.
Osés, R.R. et al. “Modelation of the quantity of deceaseds and revenues for stroke in Sagua La Grande, Villa Clara, Cuba.” The Biologist (Lima), vol. 17, 2019.
Sánchez, Á.M.L. et al. “Objective regressive regression beyond white noise for viruses circulating in Villa Clara Province, Cuba.” The Biologist (Lima), vol. 15, 2017.
Fimia, D.R. et al. “The entomofauna of culicidae and copepods approached from biological control alternatives to mathematical modeling in two central provinces of Cuba.” Anales de la Academia de Ciencias de Cuba, vol. 10, 2020, pp. 1–11.
Fimia, D.R. et al. “Mathematical modeling of the effect of atmospheric pressure on the population density of mosquitoes (Diptera: Culicidae) in Villa Clara, Cuba.” Revista de la Facultad de Medicina, vol. 68, 2020, pp. 541–549.
Osés, R.R. et al. “Forecast of new and deceased cases of COVID-19 in Cuba with an advance of 105 days.” Acta Scientific Veterinary Sciences, vol. 3, 2021, pp. 31–36.
Osés, R.R. and Pedraza, M.A. “Mathematical modeling and forecast 11 years in advance of the dumping of sand in Varadero Beach, Cuba.” Journal of Engineering and Applied Sciences Technology, vol. 3, 2021, pp. 1–3.
Osés, R.R. et al. “Study of the provincial electricity consumption of Villa Clara and its forecast 2019–2023 in Cuba.” Revista ECOSOLAR, vol. 65, 2019, pp. 32–43.
Osés, R.R. et al. “Prediction of latitude and longitude of earthquakes at global level using the regressive objective regression method.” Advances in Theoretical & Computational Physics, vol. 1, 2018, pp. 1–5.
Osés, R.R. et al. “Long term forecast of meteorological variables in Sancti Spiritus, Cuba.” Applied Ecology and Environmental Sciences, vol. 2, 2014, pp. 37–42.
Rubin, D.B. “Bayesian inference for causal effects: the role of randomization.” The Annals of Statistics, vol. 6, 1978, pp. 34–58.
Wilks, S.D. “Statistical methods in the atmospheric sciences.” International Geophysics Series, vol. 59, 1987, pp. 255–256.
Rael, V.C. “Yes to human cloning: eternal life thanks to science.” NOVA Distribution Corporation, 2001.
Rubin, D.B. “Bayesian data analysis.” Chapman and Hall/CRC, 2013.
Gelbrecht, M., Boers, M. and Kurths, N. “Phase coherence between precipitation in South America and Rossby waves.” Science Advances, vol. 4, 2018, pp. 1–9.
Chu, J.E., Timmermann, A. and Lee, J.Y. “North American April tornado occurrences linked to global sea surface temperature anomalies.” Science Advances, vol. 5, 2019, eaax0220.
Baker, R.E. et al. “Susceptible supply limits the role of climate in the early SARS-CoV-2 pandemic.” Science, vol. 369, 2020, pp. 315–319.
Reed, K.A. et al. “Forecasted attribution of the human influence on Hurricane Florence.” Science Advances, vol. 6, 2020, eaaw9253.
Zhang, G. et al. “Tropical cyclone motion in a changing climate.” Science Advances, vol. 6, 2020, eaaz7610.