Analysis of Climate Variability in a Time Series of Precipitation and Temperature Data: A Case Study in Cartagena de Indias, Colombia

: Anthropogenic climate change is a global trend, hitherto incontrovertible, causing immense social and economic damage. Although the this is evident at the global level, at the local level, there is still debate about the most appropriate analyses to support this fact. This debate is particularly relevant in developing countries, such as Colombia, where there is a signiﬁcant lack of data at the local level that require analysis and interpretation. Consequently, studies are often superﬁcially conducted to support climate change theory at the local level. However, such studies are then used to design hydraulic infrastructure, with potential catastrophic errors for human and environmental health. In this study, we sought evidence of climate change through an analysis of a series of data on temperature (maximum, mean and minimum), as well as total annual and maximum rainfall in 24 h registered at the Rafael Nuñez Airport station in the city of Cartagena, Colombia, from 1941 to 2015. The hypotheses of homogeneity, trend, stationarity and non-stationarity were analyzed. Problems of non-homogeneity and the presence of periodicity in the analyzed series were found, showing a trend and apparent non-stationarity in the original series. This could be associated with the effects of climate change. In this case, no correlation was found between temperatures and rainfall. Spectral analysis was performed for all series, and residual series were generated by extracting the harmonics of greatest signiﬁcance. It was found that the series data generated from the third harmonic are generally stationary and without trend. Therefore, the trend and non-stationarity of the original series are due to problems of non-homogeneity and periodicity in the series. In the results of the stationarity test conducted according to the Phillips–Perron criterion, all series were non-stationary. For the two additional criteria of stationarity tests, 40% were shown to be stationary, and 60% were non-stationary. Speciﬁcally, non-homogeneity problems and apparent trends associated with climate change could have negative implications for the design of drainage systems.


Introduction
The effects of climate change have been mainly associated with the gradual change in long-term average temperature. This hypothesis, which has been internationally proposed and discussed [1], impacts the magnitude of average and extreme rainfall, as well as the maximum intensities that are used in the design of engineering works. It has been argued that climate change has increased the change in the natural limits of the environment, creating more frequent extreme climate events.
Furthermore, it has been claimed that the effect of climate change on hydrological series ruins the stationarity hypothesis, which has been questioned in recent years by some authors, who stated that it has died and that it is impossible to revive it because anthropogenic effects on the climate are increasing [2]. The concern is so great that in 2010, an interinstitutional workshop involving entities that manage water in the United States

Location of the Study Area
The study zone chosen was the city of Cartagena, located in the central region of the Colombian Caribbean, at approximately 10 • 25 20" N, 75 • 32 25" W. Cartagena is a tourist city and has been declared a Historical and Cultural Heritage of Humanity by UNESCO (Figure 1). precipitation registers low volumes that do not exceed 1000 mm per year. The intra-annual regime is approximately bimodal, with a main rainy season in September, October and November and two dry seasons from December to April and in June-July. The average number of rainy days during the year is 50 to 100. Due to its physiognomy, without major mountainous accidents, the temperature distribution is very homogeneous, registering a little more than 28 °C. According to the Caldas-Lang indices for climatic classification, the climate presented is warm, semi-arid and semi-humid.

Methodology
For the study, the time series of the IDEAM station with the code name 14015020, located at the Rafael Núñez Airport in the city of Cartagena, Colombia, were selected. The series of total annual rainfall, maximum rainfall in 24 h, minimum annual temperature, mean annual temperature and maximum annual temperature were selected for the period between 1941 and 2015. For statistical analysis of the series, XLSTAT ® software was used.
The data were reviewed. There were some missing data, such as the maximum temperature and minimum temperature for the years 1941, 1942 and 1960, as well as the minimum temperature for the years 1951 and 1952. The missing data were filled using the principle of conservation of the mean and variance to calculate the descriptive statistics of all series. The initial trend of the series was analyzed with the original data, filling in the missing data. A homogeneity test was performed on all series, considering the criteria of According to the climatological atlas of Colombia from the Colombian Institute of Hydrology, Meteorology and Environmental Studies (IDEAM, for its initials in Spanish), precipitation registers low volumes that do not exceed 1000 mm per year. The intra-annual regime is approximately bimodal, with a main rainy season in September, October and November and two dry seasons from December to April and in June-July. The average number of rainy days during the year is 50 to 100. Due to its physiognomy, without major mountainous accidents, the temperature distribution is very homogeneous, registering a little more than 28 • C. According to the Caldas-Lang indices for climatic classification, the climate presented is warm, semi-arid and semi-humid.

Methodology
For the study, the time series of the IDEAM station with the code name 14015020, located at the Rafael Núñez Airport in the city of Cartagena, Colombia, were selected. The series of total annual rainfall, maximum rainfall in 24 h, minimum annual temperature, mean annual temperature and maximum annual temperature were selected for the period between 1941 and 2015. For statistical analysis of the series, XLSTAT ® software was used.
The data were reviewed. There were some missing data, such as the maximum temperature and minimum temperature for the years 1941, 1942 and 1960, as well as the minimum temperature for the years 1951 and 1952. The missing data were filled using the principle of conservation of the mean and variance to calculate the descriptive statistics of all series. The initial trend of the series was analyzed with the original data, filling in the missing data. A homogeneity test was performed on all series, considering the criteria of Pettitt [18,19], SNHT [20,21], Buishand [22] and Von Neuman [23]. Additionally, three stationarity tests were carried out: the Dickey-Fuller test [24], the Phillips-Perron test [25] and the KPSS test [26,27]. The series were corrected for non-homogeneity. The Mann-Kendall test was applied to analyze the trend [19,28,29]. A spectral analysis was performed on the series after correction for homogeneity to determine the significant periodicities of the series.
In addition, the eight harmonic components that showed the greatest significance in each of the series were selected. Residual series of analysis were generated by extracting the periodic functions with the two most significant frequencies from the original series, with three, four and eight significant frequencies. Finally, the stationarity and trend tests used in the original series were applied to the series adjusted for non-homogeneity and to the generated residual functions. In the literature, the Mann-Kendall test has been used to execute analyses of hydroclimatic trends and to detect seasonal and annual droughts in a basin, revealing significant warming throughout the basin [30]. On the other hand, the Mann-Kendall test has been used to detect the trend and change point in hydroclimatic variables during a specific period in a basin, demonstrating good performance [31].

Results
The data series taken from IDEAM are shown in Appendix A, with the original data presented in the series on the left. The values shaded in blue represent doubtful values, whereas those in gray are the missing data. The series in green were filled with random data that met the mean and variance.
In Table 1, the descriptive statistical parameters of the original series provided by IDEAM are presented. It can be seen that the coefficients of variation of the maximum rainfall in 24 h (P24) and the total annual rainfall (PT) are similar (0.40 and 0.46), and the temperatures (mean, maximum and minimum) are lower (0.02, 0.05 and 0.06). The temporal correlation coefficient shows an intermediate value for P24 and PT (0.39 and 0.45), a significant value for the minimum temperature (Tmin) (0.66) and a very low correlation with the mean and maximum temperature. In Figure 2, the IDEAM data are presented for the maximum precipitation in 24 h (P24) and the total annual precipitation (PT), plotted as a function of time, to which a linear function was fitted in order to visualize the trend. The results show low correlation coefficients and an increase in the mean values of 0.68 mm per year for maximum precipitation in 24 h (P24), as well as an increase of 8.85 mm per year for total annual precipitation (PT). Figure 3 presents the maximum, mean and minimum temperatures, in degrees Celsius, recorded at the Rafael Núñez Airport station. To visualize the possible trends, a linear function is presented. The series of maximum rainfall shows a very slight decreasing trend of 0.0068 degrees per year. The series of mean temperatures tends to be constant, with a slight decrease of 0.00004 degrees per year and very low correlations. The minimum temperature shows a significant increase of 0.0357 • C per year and a significant coefficient of determination (R2 = 0.4327).

Homogeneity Tests
To evaluate the consistency of the historical data, XLSTAT software was used, and the Pettitt, SNHT, Buishand and Von Neuman homogeneity tests were performed on the P24, PT, Tmax, Tmean and Tmin series. In all cases, the results of the Pettitt, SNHT and Buishand tests showed that the data are not homogeneous, with days on which there were changes in the data. In one case, the tests showed homogeneity (P24, Von Neumann) ( Table 2). The results of the Buishand homogeneity test for P24 and PT are shown in Figure 4. The mean P24 in the period between 1941 and 1968 was 73.075 mm, and for the period from 1969 to 2015, it was 107.451 mm. For the PT, in the period from 1941 to 1968, the mean value was 705.304 mm (µ1), and between 1969 and 2015, the mean value was 1084 mm (µ2).
To evaluate the consistency of the historical data, XLSTAT software was used, and the Pettitt, SNHT, Buishand and Von Neuman homogeneity tests were performed on the P24, PT, Tmax, Tmean and Tmin series. In all cases, the results of the Pettitt, SNHT and Buishand tests showed that the data are not homogeneous, with days on which there were changes in the data. In one case, the tests showed homogeneity (P24, Von Neumann) ( Table 2). The results of the Buishand homogeneity test for P24 and PT are shown in Figure 4. The mean P24 in the period between 1941 and 1968 was 73.075 mm, and for the period from 1969 to 2015, it was 107.451 mm. For the PT, in the period from 1941 to 1968, the mean value was 705.304 mm (μ1), and between 1969 and 2015, the mean value was 1084 mm (μ2).   The results of the Buishand homogeneity test are shown for Tmax and Tmean. The mean value of Tmax between 1941 and 1970 was 36.296 • C, and between 1971 and 2015, it was 35.016 • C. For the Tmean, in the period between 1941 and 1958, it showed a mean value of 28.150 • C (µ1), and for the period from 1959 to 2015, the mean value was 27.716 • C (µ2).
The results of the Buishand homogeneity test are shown for the minimum temperature (Tmin). The mean value of Tmin between 1941 and 1974 was 19.882 • C (µ1), and between 1975 and 2015, the mean value was 21.620 • C (µ2).
These results show the non-homogeneity of the series, in which the total annual rainfall (PT) and the maximum rainfall in 24 h (P24) increased significantly (mean value) around the year 1968. The maximum temperature shows a sharp decline around the year 1970.
Similarly, the average temperature shows a sharp decrease around 1958. However, the minimum temperature shows a sharp increase around the year 1974. These non-homogeneities are associated with changes in the location of the station or in instrumentation, among other aspects [32].

Tests of Stationarity and Trend
Using XLSTAT software, three tests of stationarity were performed (the Dickey-Fuller test, the Phillips-Perron test and the KPSS test) on the PT, P24, Tmax, Tmean and Tmin series. The results are shown in Table 3. According to the test results for the PT series, for two criteria (Dickey-Fuller and KPSS), the series was stationary, but for the Phillips-Perron criterion, the series was non-stationary. For the P24 series, two criteria showed results of nonstationarity (Dickey-Fuller and Phillips-Perron), and one showed stationarity (KPSS). For the maximum and mean temperature series, all the criteria showed non-stationarity results. For the minimum temperature series, two criteria showed non-stationarity (Dickey-Fuller and Phillips-Perron), and the KPSS criterion showed that the series was stationary. According to the Mann-Kendall criterion, the series of total precipitation (PT), maximum precipitation in 24 h (P24) and minimum temperature (Tmin) showed trends, whereas the series of maximum temperature (Tmax) and average temperature (Tmean) showed no trend.
According to the results of the stationarity test obtained with the Phillips-Perron criterion, all series were non-stationary. For the two additional stationarity tests, 40% were shown to be stationary and 60% were non-stationary. For the Mann-Kendall trend test, 60% had a trending series, and 40% were trendless.
It is interesting to compare some case studies in the literature with the present results. Viloria-Marimón et al. (2019) evaluated the trend of 19 registered series of maximum rainfall in 24 h from the department of Atlántico, Colombia. They found that 10 series showed increasing trends, and eight series showed decreasing trends. Additionally, they estimated the maximum rainfall in 24 h, applying the criteria of stationarity and non-stationarity regardless of whether the series were stationary or not. In a similar way, Tekleab et al. They found heterogeneous results, where the mean temperature had increased (0.3 • C in the rainy season and 0.6 • C in the dry season) but the time series of precipitation did not show statistically significant trends. Finally, Gonzalez-Alvarez et al. (2018), analyzing the 24 h precipitation series at the Rafael Núñez Airport station of Cartagena (the same series analyzed in this work), identified a trend over time. Therefore, they proposed the use of the frequency criterion for evaluating non-stationary and stationary conditions, showing that the non-stationary condition better represented the series. The estimates of the maximum precipitation in 24 h were higher (6-44%) than the estimates with stationarity. However, the series presented homogeneity problems (Figure 3), which generated the apparent trend.

Correction for Non-Homogeneity
Considering that all the analyzed series presented homogeneity problems, they were adjusted to coincide with the final section of the series. The series of adjusted temperatures were labeled as Tmax-R, -R, Tmin-R, P24-C and PT-C. To correct the series for non-homogeneity, the following equation was used (1): where X t is the value adjusted for non-homogeneity, x t is the initial value of the subseries used to adjust for non-homogeneity, X − 1 is the mean value of subseries 1, X − 2 is the mean value of subseries 2, S1(X) is the standard deviation of subseries 1 and S2(X) is the standard deviation of subseries 2. Table 4 shows how the temporal correlation of the minimum temperature decreases to almost zero after correction for homogeneity. The correlations of the maximum rainfall in 24 h and the total rainfall with the minimum temperature, after correcting for nonhomogeneity, decreased by more than 100%. A higher temporal correlation value appears for the maximum (0.26) and mean (0.35) temperatures after correcting for non-homogeneity. For the correlation of the maximum rainfall in 24 h and the total annual rainfall with the mean temperature, the two analyses showed similar results. The correlation of the maximum precipitation in 24 h and the total annual precipitation with the maximum temperature tended to zero after adjusting the series for non-homogeneity. Non-homogeneity alters the statistical parameters of the series and the cross-correlations between them, which can lead to biased interpretations.
In analyses of the series of total annual precipitation and maximum precipitation in 24 h, as well as maximum, average and minimum temperatures, at the Rafael Núñez Airport station in the city of Cartagena from 1941 to 2015, it can be seen that there is a significant correlation between maximum rainfall in 24 h and total annual precipitation (0.66).
Additionally, all the series presented non-homogeneity problems. The non-homogeneity of the minimum temperature series suggests a significant temporal correlation (0.66) and an intermediate correlation (0.33, 0.38) with the maximum precipitation in 24 h and with the total annual precipitation. These correlations tended to low values when non-homogeneity correction was applied.
The mean temperature showed an intermediate temporal correlation (0.35) and negative correlations with the maximum precipitation in 24 h (−0.24) and with the total annual precipitation (−0.28). Ultimately, non-homogeneity tests are one of the most important tests in climate studies. To obtain reliable results, meteorological data must be homogeneous [32].

Periodicity Analysis
To identify the periodic components of the series corrected for homogeneity, we subjected them to spectral analysis using XLSTAT software.
In Table 5 According to this criterion, the frequency of greatest significance is 3.75 years, followed by 37.5 years and, in descending order, 75, 4.69 and 6.25 years.
The eight frequencies with the highest intensities were selected for analysis. For each of the selected frequencies, its harmonic components were calculated, and the periodic functions of each frequency were generated. A residual function consisting of the original function minus the periodic functions selected from the spectral analysis was generated. The frequencies that showed the strongest signals in the spectral analysis for all the series were 75, 37.5, 18.75, 15.6. 25, 5, 4.69, 3.75, 3.41, 2.88 and 2.42 years.
By considering the weighted sum of the signal of each series, we found that the periodic frequencies of greatest significance, in descending order, were 3.75 years, followed by 37.5, 75, 4.69 and 6.25 years (the five largest). All the series analyzed from the second harmonic showed no trend. Table 5. Summary of the periodogram resulting from the spectral analysis, expressed as a percentage for each series analyzed (based on XLSTAT results). A significant presence of frequencies of 3.75 and 4.69 years was found, which could be considered to be part of the anomalies of the Pacific and tropical Atlantic. In the maximum rainfall in 24 h, the frequency of 2.42 years appears, which could be associated with the quasi-biennial oscillation of the zonal wind in the lower stratosphere, whereas the frequency of 18.75 is associated with the Pacific decadal oscillation. The 4-5-year signal associated with the tropical Pacific and Atlantic anomalies is marked in the minimum temperature series with frequencies of 3.95, 4.41 and 4.69 years. The frequency of 37.5 years is marked in the series for total annual precipitation, as well as maximum and mean temperature.

Periodicity (Years) P24-C (%) PT-C (%) Tmax-R (%) Tmean-R (%) Tmin-R (%)
In studies of the time series of aquifers [33] using piezometers, for example, cyclicities (11-22 years) related to solar cycles of 11 and 3.2 years associated with the North Atlantic oscillation (NAO) were found. Rodríguez and Llasat (1997) [34] analyzed the monthly precipitation series of Barcelona (Spain) between 1850 and 1991. They found a trend of increasing precipitation. A semiannual periodicity (0.5 years) and two periodicities between 6.0 and 7.2 years and 19 to 24 years with little amplitude and no real meaning were found. For the series of total annual precipitation, the frequencies with the greatest inten are, in decreasing order, 37.5, 3.75, 6.25, 2.88, 5, 25, 3.13 and 15 years. The Residual 8 se for total annual precipitation is the total annual precipitation series (PT-C) minus the ei periodic series with frequencies of 37.5, 3.75, 6.25, 2.88, 5, 25, 3.13 and 15 years. The ot residual series are constructed in a similar way.
Similar analyses were carried out for the other series. The results for the total ann precipitation and Residuals 8, 4, 3 and 2 are shown, with their respective trend lines a the linear function calculated for each series. The smallest slope of the residual serie represented by Residual 2, followed by Residual 3.
For the mean temperature, the predominant frequencies, in order of importance, 75, 37.5, 5, 12.5, 5.77, 3.57, 15 and 6.25 years. The slopes of the lines of the residual se are lower than the slope of the original series.
For maximum temperature, the frequencies of the greatest significance, descending order, are 37.5, 15, 75, 10.71, 2.78, 2.5, 3.57 and 25 years. The decrease in slopes of the residual series is appreciable compared with the original series.
For the minimum temperature, the frequencies of the greatest significa determined by spectral analysis are 4.69, 9.38, 3.41, 3.95, 6.82, 75, 4.41 and 2.03 years. T slopes of the original function and the residual functions remained at values of a sim magnitude. The minimum temperatures showed different frequency signals from other series in general terms.
In summary (Table 5), the total annual precipitation, the maximum temperature a the average temperature showed a significant frequency of 37.5 years. Maximum rain in 24 h and minimum temperature showed significant frequency signals of 4.69 and 3 years, respectively.
The slope of the trend line showed a general trend of a decrease compared with number of harmonics considered (Table 6). Only the maximum precipitation in 2 showed different behavior. Up to the third harmonic, all slopes decreased. For the series of total annual precipitation, the frequencies with the greatest intensity are, in decreasing order, 37.5, 3.75, 6.25, 2.88, 5, 25, 3.13 and 15 years. The Residual 8 series for total annual precipitation is the total annual precipitation series (PT-C) minus the eight periodic series with frequencies of 37.5, 3.75, 6.25, 2.88, 5, 25, 3.13 and 15 years. The other residual series are constructed in a similar way.
Similar analyses were carried out for the other series. The results for the total annual precipitation and Residuals 8, 4, 3 and 2 are shown, with their respective trend lines and the linear function calculated for each series. The smallest slope of the residual series is represented by Residual 2, followed by Residual 3.
For the mean temperature, the predominant frequencies, in order of importance, are 75, 37.5, 5, 12.5, 5.77, 3.57, 15 and 6.25 years. The slopes of the lines of the residual series are lower than the slope of the original series.
For maximum temperature, the frequencies of the greatest significance, in descending order, are 37.5, 15, 75, 10.71, 2.78, 2.5, 3.57 and 25 years. The decrease in the slopes of the residual series is appreciable compared with the original series.
For the minimum temperature, the frequencies of the greatest significance determined by spectral analysis are 4.69, 9.38, 3.41, 3.95, 6.82, 75, 4.41 and 2.03 years. The slopes of the original function and the residual functions remained at values of a similar magnitude. The minimum temperatures showed different frequency signals from the other series in general terms.
In summary (Table 5), the total annual precipitation, the maximum temperature and the average temperature showed a significant frequency of 37.5 years. Maximum rainfall in 24 h and minimum temperature showed significant frequency signals of 4.69 and 3.41 years, respectively.
The slope of the trend line showed a general trend of a decrease compared with the number of harmonics considered (Table 6). Only the maximum precipitation in 24 h showed different behavior. Up to the third harmonic, all slopes decreased.
In Appendix B, the results of the stationarity and trend tests are shown for total precipitation (PT); adjusted total precipitation (PT-C); and the corresponding residual functions with two harmonics (PTC-R2), three harmonics (PTC-R3), four harmonics (PTC-R4) and eight harmonics (PTC-R8). Likewise, stationarity and trend tests were carried out for the other series under study (P24, P24-C, P24-C-R2, P24-C-R3, P24-C-R4 and P24-C-R8 for total rainfall in 24 h and all temperature series (maximum, mean and minimum)), as shown in Appendix B. According to the results of the Phillips-Perron test, all the series are non-stationary. Therefore, this test was not sensitive in the case of our study and was considered in the analysis of the results.
According to the Dickey-Fuller test, the original series of total precipitation was stationary after adjustment for non-homogeneity. When the periodicities (two, three, four and eight harmonics) were extracted, stationary series were revealed. According to the KPSS test, all total precipitation series were stationary. The Mann-Kendall trend tests showed that the original series of total precipitation had a trend, but after adjustment for non-homogeneity, no trend was recorded for any of the residual series.
The original series of maximum precipitation in 24 h was non-stationary according to the Dickey-Fuller test. On the other hand, when it was adjusted for non-homogeneity, the result was stationary, and when the periodicities (two, three, four and eight harmonics) were extracted, the result showed a stationary series. According to the KPSS test, all 24 h precipitation series were stationary. The Mann-Kendall trend test showed that the original 24-h precipitation series had a trend, but after it was adjusted for non-homogeneity, no trend was recorded in any of the residual series.
Moreover, according to the Dickey-Fuller test for the original maximum temperature series adjusted for non-homogeneity and after two harmonics (Residual 2) had been extracted, a non-stationary series was found. The series corresponding to Residuals 3, 4 and 8 are stationary. The KPSS test gave non-stationary results both for the original maximum temperature series and the series adjusted for non-homogeneity. The series of Residuals 2, 3, 4 and 8 gave stationary results. The Mann-Kendall trend tests showed that for all series, no trends were recorded.
The results of the Dickey-Fuller test showed a non-stationary series both for the original mean temperature series, as well as the series adjusted for non-homogeneity. The series corresponding to Residuals 2, 3, 4 and 8 provided stationary results. According to the KPSS test, the original mean temperature series and the series adjusted for non-homogeneity showed non-stationarity results. The series of Residuals 2, 3, 4 and 8 showed stationary results. The Mann-Kendall trend tests showed that in the series of mean temperature adjusted for non-homogeneity, there was a trend, whereas in the other series, no trend was found.
Finally, according to the Dickey-Fuller test for the minimum temperature series adjusted for non-homogeneity, the result was a stationary series. The other series gave non-stationary results. According to the KPSS test, all the minimum temperature series gave results of stationarity. The series of Residuals 2, 3, 4 and 8 gave stationary results. The Mann-Kendall trend tests showed a trend in the original minimum temperature series, whereas in the other cases, no trend was shown.
In summary, all the series were insensitive to the Phillips-Perron criterion because all the series in any condition were non-stationary. Here, 60% of the original series (PT, P24, Tmax, Tmean and Tmin) were shown to be non-stationary, and 40% were stationary. In 60% of the cases, a trend was shown. When the series were adjusted for non-homogeneity (PT-C, P24-C, Tmax-R, Tmean-R and Tmin-R), 50% remained non-stationary, and the other 50% were stationary. In this case, only 20% showed a trend, revealing the strong impact of non-homogeneity on the trend.
After extraction of the two harmonics of the greatest significance from the series, the percentage of non-stationary series was 20%, 80% were stationary and no series showed a trend. After extraction of the three, four and eight harmonics of greatest significance from the series, 10% of the series displayed non-stationarity, and 90% were stationary. None of the series showed a trend (Appendix C).
In earlier work, Viloria-Marimón et al. (2019) [35] studied the maximum precipitation data in 19 pluviographs in a region of Colombia. They evaluated the trends of the series and found positive and negative trends. Additionally, they applied the concepts of frequency for stationary and non-stationary series, estimating the values of maximum precipitation in 24 h for different return periods. However, the homogeneity of the series and the periodicity contained in the records, which could affect the reliability of the results, were not considered [32]. When these adjustments are not applied, apparent trends can be obtained. In the case of negative trends, this could lead to underestimations of the design parameters, leaving communities exposed and vulnerable to the original problem. In the opposite case, an overestimation of the design parameters could lead to negative economic consequences for the designs [36].

Discussion
For the stationarity test, we applied several criteria. The results obtained with the Phillips-Perron criterion were all non-stationary series; when the two additional criteria were applied, the results obtained were 40% stationary and 60% non-stationary; and when the Mann-Kendall criterion was applied, 60% were found to be trending series, and 40% were trendless.
It is interesting to compare some case studies in the literature with the present results. As previously mentioned, the work by Viloria-Marimón et al. (2019) [35] evaluated the trend of 19 registered series of maximum rainfall in 24 h in the department of Atlántico, Colombia. They found that 10 series showed increasing trends, whereas eight series showed decreasing trends. Additionally, they estimated the maximum rainfall in 24 h by applying criteria of stationarity and non-stationarity, regardless of whether the series were stationary or not. In a similar way, Tekleab et al. (2013) [19] analyzed the hydroclimatic trends in the upper part of the Blue Nile basin in Ethiopia. They found heterogeneous results, where the mean temperature had increased (0.3 • C in the rainy season and 0.6 • C in the dry season) but the time series of precipitation did not show statistically significant trends. Finally, Gonzalez-Alvarez et al. (2018) [11] evaluated the 24 h precipitation series at the Rafael Núñez Airport station of Cartagena and identified a trend over time. Consequently, they proposed the use of the frequency criterion for non-stationary and stationary conditions, showing that the non-stationary condition better represents the series. Nevertheless, according to the results of this research, this is an apparent trend.
Therefore, non-homogeneity problems and apparent trends associated with climate change could have negative implications for the design of drainage systems. There is a growing tendency to claim that the theory of non-stationarity has died due to the effects of climate change. This has led to inadequate legislation in some countries, such as Colombia, forcing engineers to apply new and developing calculation hypotheses.
The main limitation of this study is the data, as they were obtained from a single station.

Conclusions
The objective of the present work was to analyze descriptive statistical parameters as basic elements for making decisions about the effects of climate change in an undeveloped area in Latin America and to realize a statistical analysis of a time series of the basic weather parameters in the city of Cartagena normally used to design drainage systems. Currently, the stationarity criterion is in force for calculating the hydrological design parameters of the Rafael Núnez Airport station in the city of Cartagena. In the case of this study, a non-homogeneity problem was highlighted, producing apparent trends possibly associated with climate change effects. Nevertheless, the studied series are non-stationary, and a new paradigm is most likely required for estimating the design parameters of engineering works.
The series are non-stationary due to the omission of homogeneity and a periodicity analysis of the original series. This could lead to misevaluated engineering design parameters, which could result in unpredictable negative consequences. Finally, the apparent trend shown by the original series is due to problems of non-homogeneity and the influence of periodicity (climate variability), whereby no evidence of climate change effects was found.

Acknowledgments:
The authors would like to acknowledge the IDEAM for providing the data on which this study was based.

Series Adjusted for Non-Homogeneity
Years P24-C PT-C Tmean-R Tmax-R Tmin-R