Meteorological Temperature and Humidity Prediction from Fourier-Statistical Analysis of Hourly Data

The temperature readings for all the 365 days and the 24 hours may be fitted through a 3 × 3 matrix (the so-called T-matrix). The mean square deviation between this fit and the actual meteorological measurements is smaller than three degrees Celsius. Four entries of this (nonsymmetric) matrix may be fixed by other means, leaving only five independent components. However, the same method applied to the humidity measurements produces a larger mean square deviation. A strong stochastical connection is found between the T-temperature matrix and the U-humidity matrix. The computer program, in C, may be used to adjust a (2M + 1) × (2m + 1) matrix simply by changing the arguments at the command line and has been tested with m and M ranging from zero to 11 (eleven) (more than 24 readings per day are necessary for larger values of m). The physical meaning of these constants is given only in the case m = M = 1. Our results have also been connected to fundamental cosmological properties: Earth’s orbit, the ecliptic angle, and the latitude of Querétaro (or whatever geographical location is chosen). A separate program calculates the angular position of the Sun as measured in the sky of Querétaro, to determine the length of the day or the mean value of the solar cosine. This work introduces several new variables which happen to be stochastically connected.


Introduction
Statisticians must face problems at three different levels. e first one concerns the connection between meaningful statistical data and probability distributions. e second, first documented and discussed by Simon de Laplace [1] in his "Essai philosophique sur les probabilities," is related somehow with "corrupted" data, due to some particular interest in influencing the facts. Examples may be found in medicine at all levels, whether in the accounts of pharmacological substances or reports on weight and shape of the body control, and even in surgery. Politics has also made contributions in this second level. e third is related with "accidentally" corrupted data: the operators do not endeavor to make careful readings because, on the average, the final results must be meaningful anyway. Currently, the solar radiation prediction proposes machine learning algorithms to predict the hourly solar irradiance [2], but the historical data are fundamental for a correct operation [3]. Additionally, the daily basis in spatial-temporal analysis uses single datasets that provide data for a certain region, valid for a certain time period [4]. Consequently, new and efficient methodologies are required, to be able to achieve the inference for nontrivial models within an affordable time [5], including embedded, distributed, and parallel computing of inferential processes and anomalies from global or local distributions [6]. Recently, the wind speed prediction using Wind Speed and Turbulence Intensity-based Recursive uses artificial neural network [7] but requires that the maximum time-resolution of turbulence intensity should be longer than the prediction interval [8]. Likewise, the temperature prediction uses the artificial neural network model and it was able to predict the indoor temperature [9], but the majority of previous studies applying either ambient or ground have tended to emphasize the structural improvement of individual forecasting models without considering the periodicity of data [10]. In this approach, annual average soil temperatures are determined by air temperature, solar radiance, wind speed, and relative humidity [11]. e objective of this paper is to provide the theoretical framework to meaningfully analyze data on temperature and humidity in general. us, the present work was carried out choosing the case of Querétaro as a particular example, but this model can be used in other sites, i.e., their approach is universal and may be applied to other sites as such. erefore, this scheme for Fourier-statistical analysis of hourly data considers daily temperature variations as superposition of annual average temperatures prediction (long-term climate impact) and daily temperature amplitude predictions (short-term climate impact) for various climates on agricultural production [9] and for meteorological temperature and humidity prediction on environmental-scale applications. e results of this work are useful for establishing a solar energy research program, for the buildings and road construction industries, for agriculture, and for all those people whose work depends on weather.

Statistical Fit.
e fitting principle is very simple; any place on the Earth is subjected to two types of physical conditions according to their predictability: the systematic type, as the flight of the Sun through the sky, which can be described to a high degree of accuracy, and the stochastic type related to physical conditions varying randomly day after day, such as the horizontal air transport (wind) coming from the oceans or from the polar circle and the precipitations, including phenomena such as the storms or tornados. e resultant weather depends thus on such superposition. e daily systematic changes should produce systematic variations in the weather variables, such as the temperature or the humidity, following an approximate periodical behavior. e objective of this paper is mainly related to the systematic part so that when we speak of "variables," we really mean the causal part. In Section 3, we shall propose an operational definition of this "causal part" concerning the meteorological measurements of temperature and humidity. e Sun's flight will be taken as causal. Since the variables are periodic (more precisely, approximately periodic, since the Sun's trajectory along the sky is slightly different each 24 hr period) in time, they can be expanded in a Fourier series whose m-th partial sum is as follows: where is the angular frequency in radians, f is the frequency, and b − a is the time period. We expect that for a small number m, our m-order Fourier expansion, either on temperature or on humidity, can reproduce accurately the meteorological measurements. Our very first step was to write a computer program to calculate numerically, from the raw data, the Fourier coefficients A n and B n , defined by the usual Fourier integrals as follows: where the time period (b − a) is 24 hours. e program runs fast, and it is easy to obtain all the (2m + 1) coefficients. A similar analysis was made by Carson [12] in his pioneering work. We realized that a large amount of coefficients were necessary to obtain a good fit and decided to try a different method: to fit a "trigonometric polynomial" of a fixed degree m by the least-squares method, that is, to choose the coefficients A n and B n so that the squared distance between Φ m (t i ) and the measured variable T i at time t i , added over all the measurements, is a minimum. is total distance σ is given by where n i represents the number of measurements made to obtain the experimental pair (t i , T i ) and N is the total number of data. Normally, n i � 1, but frequently, on holidays, the temperature or humidity measurement is lacking so that n i � 0. If we are working with averaged numbers, over a period of a month, for example, n i is the number of experimental points used to establish the average value of T i . e definition of σ is the standard deviation between the experimental data T i and the "expected" Φ m (t i ). e following equations result from the condition of minimum and are represented by  n i y i sin jωx i , for j � 0, . . . , m.

(7)
From equations (6) and (7), we can constitute a system of 2m + 1 equations with 2m + 1 as the unknown variable and it may be written in the matrix form as follows: where the matrix E is defined by e inhomogeneous vector F is described by And the vector C, corresponding to the Fourier coefficients to be fitted, is represented by It is the solutions C to equation (8), obtained through Cramer's rule, that will be tabulated in this work along with the corresponding value of the mean square deviation σ given by equation (5). ere is an important difference between the least-squares approach, based on equations (6) and (7) and the first one, which employs the Fourier coefficients given by equations (2)- (4). Suppose that we have adjusted a trigonometric polynomial of degree m using the Fourier coefficients of the equations (2)-(4). If now a polynomial of degree m + 1 is desired, all we have to do is to calculate the A m + 1 and B m + 1 coefficients using equations (2)-(4). On the contrary, if we have already fitted an m-degree trigonometric polynomial using equations (4) and (5) and a polynomial of degree m + 1 is desired, all the (2m + 1) + 1 coefficients must be calculated. A word of caution should be mentioned concerning the least-squares approach. e numbers of hourly readings averaged or not are always, at the most, 24. If the number 2m + 1 of coefficients fitted equals the number of readings, the trigonometric polynomial, although passing through all the experimental points, shows unacceptably strong variations of temperature or humidity. Smooth and acceptable variations are obtained only with small m numbers. Anyway, our physical interpretation will be restricted to m � 1.

Statistical Considerations.
e hourly readings of temperature and humidity are sometimes subjected to important fluctuations along the day. Furthermore, as the readings of consecutive days, corresponding to the same time, are compared, more fluctuations may be revealed. It is a usual practice, however, to calculate monthly averages for each one of the 24 hrs. is procedure leads to a smooth variation, along the day, of both temperature and humidity, but has a drawback that only twelve averages are obtained for the whole year, one for each month, a fact that prevents us from obtaining a detailed variation along the whole year. Instead, we have defined an l-average, for l is a nonnegative integer. Let i stand for a i-th day of the year, 1 ≤ i ≤ 365. e laverage of temperature at time t for the i-th day, T (t, i, l), is the average of the temperature readings at time t corresponding to the days i − l, i − l + 1,. . .,i + l. Since our data are limited to the weather data of 1994, T (t, 365, 1) is the average calculated over those days whose numbers are 364, 365, and 366, that is, days with numbers 364, 365, and 1, the integer reduction modulo 365. Our programs receive directly the number l in the command line and their running time is not sensitive to the value chosen for l. We have found that l � 5 produces already smooth curves. e program may certainly be run with l � 0. For practical reasons, we shall not present our results concerning the Fourier coefficients A n and B n as they appear in equation (1), but in terms of C n and δ n (in degrees) defined by

Advances in Meteorology
A n · cos(nω · y) + B n · sin(nω · y) � C n · cos nω − δ n . (12) at is, according to Our Fourier polynomials are thus given by Table 1, because of space limitations, contains only the Fourier coefficients for the 15th day of each month with l � 15 and m � 0, 1, 2, and 3, with the corresponding values of the mean square deviation σ calculated from the year extension of equation (5). e σ values should be compared with the yearly average of temperature <T> � 18.3 C. ere are several variables which measure the atmosphere's water content. ese are, for example, the two temperatures from the psychrometer (actual and dew temperatures), the humidity H deduced from the previous temperatures and connected with the wetting power of the atmosphere, the water-vapor density ρ e , possibly the simplest to "understand," ranging from 0.7 to 17.1 g/m 3 , and the watervapor partial pressure p e whose values ranged from 0.95 to 23.85 mb (1 mb � 100 Pa) during 1994 in Querétaro. ese last two variables are connected by the Clapeyron's equation of state for ideal gases: p e � ρ e RT/m * , where T is the absolute temperature and m * � 18 g is the molecular weight of water. e last variable, for reasons to be given below, seems to be the most useful variable. Table 2 gives the Fourier coefficients for the water pressure p e and for m � 0, 1, 2, and 3, with the corresponding values of the mean square deviation σ. e σ values should be compared with the yearly average of water-vapor pressure <p e > � 11.8 mb.

Fourier Fitting the Fourier Coefficients.
After fitting a trigonometric polynomial of degree m, equation (14), to the meteorological data of each one of the 365 days, we obtain a matrix of size 365 × (2m + 1). Each one of its columns gives the daily variation along the year of the corresponding coefficient. e systematic part should repeat year after year, so that we can fit, for each one of the 2m + 1 columns, the 365 data by another trigonometric polynomial of degree M, obtaining 2M + 1 coefficients. e complete fit involves then a matrix (2M + 1) × (2m + 1). ere is no problem in running the program on the data with large values of m and M, let us say, m � 10 and M � 20. In this case, the fit is better and the mean quadratic deviation decreases, but the physical meaning of each coefficient becomes darker: as bigger values of m and M are chosen, not only the systematic, but also the stochastic variations of weather are fitted. e limit on the number m mentioned at the end of Section 2 is not applicable to number M, since in this case, we have 365 data to be fitted. Tables 3 and 4 give these fitted matrices to the temperature and water-vapor pressure, for m � 0, 1, 2 and M � 0, 1, 2. e meaning of Tables 3 and 4 is as follows: let us suppose the Fourier expansion of the water-vapor pressure for the k-th day of the year corresponding to m � 1 is desired, that is described by where ω � 360/24 * 60 * 60, since the phase is given in degrees and the time t in seconds. According to Table 4, we where Ω � 360/365. Table 4 also indicates that the mean square deviation between equation (15), obtained through equations (16)-(18), for k � 1, 2, 3,. . ., 365, and the whole set of measured p e values is σ � 1.1984.

ermodynamical Interpretation.
In this section, the analysis of our statistical results concerning temperature and humidity, both as daily and yearly functions, is studied in connection to the motion of the Sun along the sky of Querétaro. As we have previously stated, we are interested in the causal part of our variables. It is known [13] that the temperature of the atmosphere decreases as a function of altitude up to the tropopause where the temperature has a minimum. We are interested in a vertical column of atmosphere, of unit cross section, lying above the surface and whose upper bound is approximately at an altitude of some 20 km, slightly above the tropopause, which is located at a height of about 17 km in the tropics [13]. A small fraction of the atmosphere lies above 20 km since the pressure at such height is 0.06 atm [14]. e first law of thermodynamics "for any infinitesimal change" according to Planck [15] is as follows: where Q stands for the heat added to the system and W for the work done on the system. As a first approximation, we assume that, on average, equal quantities of air, in the same thermodynamic state, enter to and exit from our system so that the system is closed. We also assume that W � 0 since the work done on the system involves the pressure at the top, high above the tropopause, where the pressure is very low. We will also include in our system a layer of the ground below the surface, whose thickness is to be determined, on which no work is performed, either.

Heat Input
Balance. e heat added to our system has its main source from the radiated heat by the Sun while going through the sky of Querétaro. e heat radiated by our system to the outer space, due to its temperature, the socalled "radiation at long wavelengths," is a negative contribution. We shall assume in this approximation that the thermal loss may be completely accounted for by the Stefan-Boltzmann radiation law. e results will show to what extent is this hypothesis valid. e solar radiation intensity at the top of the atmosphere is given by the socalled solar constant S co � 1360 W/m 2 .
is is an average value, for not only the solar activity is continuously changing [16] (neglected in this work), but also the distance between the Earth and the Sun varies, and with it, the intensity, which  15th month   e C coefficients are given in°C and the phase δ in degrees; m is the number of harmonics, and σ is the mean square deviation computed from equation (3).
Advances in Meteorology goes as the square of the inverse distance. In the following S c � S c (t), t is the time along the year. Let us write θ 0 (t) for the zenith angle, that is, the angle between the zenith and the unit vector u 0 (t) locating the position of the Sun in the sky of Querétaro at time t. e total power from the Sun is given by where α ∈ <0, 1>, the so-called albedo, represents the fraction from the solar radiation intensity which is reflected back by the Earth's atmosphere. However, 1 − α should not be confused with the absorptivity [17] defined for a system submitted to black body radiation, mainly because the solar radiation at the Earth is not isotropic. It is known that the atmosphere absorbs a small amount from the solar energy, being essentially transparent, so that the main absorption occurs at the ground [18,19]. e clouds, on the other side [19], increase the reflectivity of the atmosphere and, in this way, increase the value of our α parameter.

Earth's Translational and Rotational
Motion. e Earth moves along an elliptical orbit whose eccentricity is e � 0.01678 [20], with a variable orbital speed, complying with Kepler's law of equal areas in equal times. However, this ellipse is not fixed but rotates slowly: its big axis, called apsis, advances 11 seconds of arc each year. e direction of this apsis is given by the directed vector from the Sun to the perihelion of the Earth's orbit. e Earth rotates around the poles axis with a constant angular velocity [21]. e direction of this axis may be given with respect to the orbit of the Earth so that two angles are necessary. e first angle, ϑ, is measured from this direction to that perpendicular to the orbit. According to Baize [20], we have ϑ � 23°26″40.62″″, whose time variation gives rise to the nutation of the Earth (neglected in this work). e second angle, ϕ, connected with the projection of this direction along the plane of the orbit-the line of solstices-whose time variation produces the so-called precession of the equinoxes, is specified by recording the angle between the winter solstice and the perihelion. According to Baize [20], ϕ � 10°20″15.18″″ and ϕ � 50.26″″ per year. e town of Querétaro extends from 20°30″ to 20°54″ of north latitude and from 100°17″ to 100°36″ of west longitude. We have taken average values for these angles, 20.7 for the (north) latitude and 100.44 for the (west) longitude. e longitude is used in this work only to connect the civil time, used to record the meteorological variables, with the solar time.

Heat Output.
e emitted power, that is, the thermal loss per unit time, may be written as We are not assuming the ground to behave as a gray body [16,22,23] so that 1 − β is not the corresponding emissivity. If the temperature of the ground is T, it certainly emits thermal radiation whose intensity is given by σT 4 , but it also receives thermal radiation from the atmosphere, whose temperature changes continuously with height. Such an exchange of heat, being interior to our system (ground plus atmosphere), gives a null contribution. e approximation behind equation (21) is that the temperature T of the ground is a parameter identifying the thermodynamic state of our system. If β stands for the probability of an emitted photon to be reabsorbed, 1 − β should represent the probability of it going to the outer space. Equation (21) is an average of all these thermal losses from our system to the outer space. Since water vapor absorbs strong radiation at those wavelengths associated with the Earth's temperature [13,18,19,[24][25][26][27], our β, connected with the probability of being absorbed by the atmosphere, is also related with the atmosphere's absorptivity. An important difference between α and β is that the first one is connected with short wave length coming from the Sun, corresponding to the Sun's temperature of 5985 K [24], while the second with the long wave length radiation from the Earth. Clear days should be connected with small values of α and β while the cloudy ones with the opposite. e total heat input per unit time is given by the sum of equations (20) and (21), obtained as follows:  (22) corresponding to a very long period of time should vanish since, on average, our system radiates as much energy as it receives. Failure to fulfill this requirement should warm up or cool down the Earth, which is contrary to the constant values of temperature averages. α and β are functions of the time t but will be, in the sequel, substituted each by the appropriate average. As a first approximation, we shall suppose that such time-integral, extended to the time duration of any day, vanishes: e time-integral of equation (23), corresponding to a very long period of time should vanish since, on average, our system radiates as much energy as it receives. Failure to fulfill this requirement should warm up or cool down the Earth, which is contrary to the constant values of temperature averages. α and β are functions of the time t but will be, in the sequel, substituted each by the appropriate average. As a first approximation, we shall suppose that such time-integral, extended to the time duration of any day, vanishes: where τ � 86400 s is the time duration of a day. It follows where the above equation follows from a Fourier expansion of S c (t) cos[θ 0 (t)] and of T 4 (t) in the form given by equation (15): C κ 0 is the C 0 coefficient for (S c (t)/S co ) cos [θ 0 (t)] and C R 0 (R for "radiation"), that for T 4 (t). e ratio c given by equation (25) is shown in Figure 1. e most striking fact is that α and β are approximately equal since this ratio varies around c � 1. ere is a sharp maximum around day number 349, with c � 1.38, and a flat minimum around day number 173 with c � 0.90. Since we expect α and β to increase with the water-vapor content in the atmosphere, which has its maximum value around day number 170 and its minimum value at the beginning or end of the year, we see that α increases more rapidly than β. e C κ 0 and C R 0 coefficients may be approximated by the corresponding daily averages of cos[θ 0 (t)] and of T 4 (t). According to equation (25), c depends on the ratio of two quantities which may be obtained easily. Figure 1 states clearly that the ratio, in the sense of equation (25), of the recorded temperatures (numerator) to the available solar energy (denominator) is large in winter or small in summer. ese facts are consistent with a positive (negative) heat input in winter (summer). However, the main wind transport in winter and summer is in the east-west direction, and the wind component along the south-north direction is small and alternates along the day. To account for these differences, we must accept that the air, coming from equal latitudes, is warmer in winter and colder in summer.
In Figure 1, the plot of equation (25) asserts that temperatures are low in summer, and this is true because of the rains and of the clouds covering the sky, thus preventing the temperature to increase. Our meteorological data concerning the clouds are scarce: For each month, the number n o of overcast days, n h of half-overcast days, and n c of clear days are tabulated. For Tables 1 and 2, we plotted n c (i) + (3/4)n h (i) + (1/2)n o and obtained a figure similar to Figure 1, where we have chosen the ratios 1, 3/4, and 1/2 arbitrarily. is fact may be considered in the opposite sense: c − 1 could be used as an objective measure of cloudiness. We shall return to this problem as follows.

Internal Energy U s of the Ground.
e solar heating of the ground rises its temperature, rising also the temperature below the surface. e heat propagates first downwards, warming the lower layers, and then upwards, during the night. e temperature of the ground as a function of time and depth may be denoted by T(y, t), where the ground corresponds to the negative values of y and the boundary between the ground and the atmosphere to y � 0. Joseph Fourier discovered "Fourier analysis" precisely working with heat problems in 1807 [28], and his book "La eorie Analytique de la Chaleur" includes interesting geophysical problems as the discussion of "the terrestrial temperatures" or the penetration depth of the temperature changes. e French Academy rejected for 17 years the publication of Fourier's work in the Memoires of Academie des Sciences, and finally in 1824, Fourier, already Secretary of the Academy, published in those Memoires [29] his work, in the same form as had been subjected in 1811. To determine T (y, t), we assume that the ground temperature satisfies the Fourier law of temperature diffusion [28,30] and is represented by where χ is the thermal conductivity of the ground, ρ s is its density, and c vs is its constant-volume specific heat. We suppose the surface temperature of the ground to vary as a Fourier harmonic as e temperature distribution should also satisfy the condition as follows: e solution to equation (26) satisfying equations (27) and (28) is given by

T(y, t) � T p + T n e y/y n
where from equation (29), we obtain Since heat equation (26) is a linear differential equation, for a surface temperature given by a Fourier series as in equation (1), with terms corresponding to frequencies ω, 2ω, 3ω,. . ., the temperature distribution in the ground should be given by the sum of the corresponding terms, having each one the form of equation (29). Let us stand for the specific internal energy of the ground. From the first law of thermodynamics, equation (19), we have We assume the ground to be rigid so that the work done on it is neglectable. e added heat Q may be written in terms of the constant-volume specific heat c vs to have dU s � c vs dT , (32) where dT is the infinitesimal temperature change. Assuming c vs to be constant along the small interval of ground temperature changes, the energy U s is given by U s � c vs T, where we have placed the integration additive constant equal to zero. e total energy U (n) s of the ground, corresponding to the n-th harmonic, may be written as follows: where we have used equations (29)-(32) to calculate the integral. If the surface temperature has several harmonics, the internal energy of the ground will be a sum of terms having the form of equation (33), one for each harmonic. Carson [12] has verified the validity of the heat equation concerning the space-time distribution of temperatures, equation (26), and the Fourier decomposition of the solution, although he has pointed out that small deviations were observed near the ground surface. We assume that these deviations are neglectable.

Results and Discussion
4.1. ermal State of the Atmosphere. According to [24], the chemical composition of the atmospheric air is approximately fixed for the first 100 km, that is, there is a fixed proportion of chemical species. e more abundant ones are nitrogen (N 2 ), oxygen (O 2 ), and argon (Ar). ere are also gases with varying proportions such as the water vapor (H 2 O) and carbon dioxide (CO 2 ). From the thermodynamic point of view, this means that the thermodynamic state at each point in the atmosphere may be determined by giving two thermodynamic variables, which we chose to be the absolute temperature T and the pressure p (water also requires the proportion of vapor, liquid, and solid).
Our first attempt to obtain the internal energy of the atmosphere depended on the barometric equation (33) and the adiabatic convection model for the temperature of the atmosphere [13] which establishes a connection between pressure and temperature. From this, a given ground surface temperature was sufficient to calculate the vertical profiles of pressure and temperature. e internal energy is obtained from a simple calculation. However, to follow the daily variations of temperature, neglecting even the radiated heat, the equations required the Sun to deliver 6 times more energy than ours. We realized that our model submitted the atmosphere to a uniform temperature shift as an instantaneous response at all heights. For this reason, we sought a temperature distribution propagating upwards, whose amplitude is a decreasing function of height. Let T a (y, t) be the temperature of the atmosphere at altitude y and time t. We assume that its average T p (y) with respect to time is the mean temperature as function of height reported by meteorologists and that − zT/zy corresponds to the so-called lapses rate. We may thus write T a (y, t) as where T(y, t) is connected with the daily oscillation of temperature produced by the Sun and transmitted through the atmosphere mainly by convection. e temperature distribution given by equation (34) should reduce at the surface of the ground to equation (27) and satisfy a limiting condition similar to that imposed by equation (28) and is represented by T(y, t) � T p (y), for large positive values of y. (35) e temperature changes at the surface of the ground should propagate upwards, to the atmosphere, and, although having pretty complicated details because of the convectionlike transport, should have, in the average, a simple spacetime behavior. For this reason, we assume that T a (y, t) may be written as: 8 Advances in Meteorology is temperature distribution, having the same form as equation (29), satisfies the imposed conditions. However, in the preceding case, the y n was connected with known thermal properties by means of equation (30). ere is no such connection in this case, although we still assume that Y n � Y 1 /√n. e physical interpretation of Y n is simple: it is the height in the atmosphere where the temperature oscillation equals the oscillation at the ground surface times exp(− 1). If equation (31) proves to be useful, Y n becomes an interesting parameter since it measures the thickness of the lower atmosphere which is submitted to the daily temperature variations. Concerning the atmospheric pressure we shall assume that its logarithm is a linear function of height.
is assumption is not only consistent with the table published by Petterssen [13] or that by Battan [18], taken from reference [14], but also physically acceptable as follows: We have obtained the slope µ � 1.445 × 10 − 4 m − 1 by fitting a least-squares straight line to the data published by Battan [28]. e equations (36) and (37) allow us to calculate the internal energy U a of the atmosphere. e pressure changes along the column are given by the barometric equation defined by dp dy Since the largest atmospheric pressures are close to 1 atm, we can take Clapeyron's equation of state for ideal gases defined by where p is the pressure, ρ is the density, R is the gas constant, T is the temperature and m * is the molar mass of the gas. Substituting the density ρ from equation (38) in the barometric equation (39), we obtain Comparing equations (37) and (40), we obtain µ as follows:

Internal
Energy U a of the Atmosphere. Since the energy per unit volume is ρc v T, the total energy U a in the atmosphere is given by is integral may be evaluated readily after substituting p(y) from equation (37) and T a (y, t) given by equation (36). e result is Assume that the argument of angle is tan(Δn) � 1/(1 + µY n ). e total energy of the atmosphere, in this first approximation, is the sum of an average energy, given by the first term of equation (43) and a perturbation produced by the daily solar heating, given by the second. e time derivative f n (t) of the sum of equations (33) and (43) is Finally, if the surface temperature of the ground is given by a Fourier expansion, written as equation (14), the first law of thermodynamics for our system may be written by equating equation (23) to the sum of the corresponding terms given by equation (44); the result is where it involves three unknown parameters: α, β, and Y 1 . We shall not try to write equation (45) for a choice of three different times t1, t2, and t3 and solve the system since the results would depend strongly in that choice. Rather, we consider the right hand side of equation (45) as a given function g(t, α, β) of α and β, as stressed by equation (23) and choose the values of these parameters so that the distance between g and f is defined by e equation (46) is a minimum. Here, τ is the time duration of a day and is given by τ � b − a, according to equation (1).
is definition allows the use of Fourier methods to solve this problem. In fact, let the mathematical rule associated 〈f | g〉, the scalar product between any two integrals functions, be defined by It follows that d 2 (f, g) � 〈f − g | f − g〉. Our object is to determine α and β so that for a given value of Y 1 , the squared distance given by equation (46) is a minimum. Equating to zero the derivatives of d 2 with respect to α and β, we obtain Since zg/zα � − S c (t) * cos θ Θ (t) and zg/zβ � σT 4 (t), the system of equation (48) may be written as where k � S c (t)/S co . e determination of α and β depends on the calculation of the five different scalar products appearing in equation (49), a very simple problem if we have Fourier expansions, in the form given by equation (1), for S c (t) cos[θ Θ (t)], T 4 (t), and f(t). In fact, if Φ(1) and Φ(2) are Fourier expansions of two different functions, their scalar product 〈Φ(1) | Φ(2)〉 is given by After writing equation (44) in the form of equation (1), it follows that the Fourier expansion of f(t) depends on that of T(t). To obtain expansion of T 4 , the temperature entries originally given in°C are transformed into absolute temperatures and squared twice, and from these data, the Fourier coefficients were determined. Finally, S c (t) cos[θ Θ (t)] is calculated from our program. Equation (50) enables us to obtain α and β for each day of the year. Since we do not know the value of Y 1 , we should run our programs for different values and then choose the one which seems to give the best results. It is true that Y 1 is a parameter in equation (46) and that we could also obtain its value by equating the derivative of d 2 with respect to Y 1 to zero. However, the Y 1 should be chosen so that the resulting albedo α values are acceptable. We wrote our program to determine α and β and also the distance d(α, β) given by equation (46) for a given value of Y 1 . From the tabulated values, a second program determined the year average of these three variables, determining also their maxima and minima and the day of their occurrence. ese two programs were iterated for 1000 different values of Y 1 , ranging from 10 to 20,000 m. e list of Y 1 values was chosen so that the ratio of one Y 1 value to the preceding one was a constant, that is, Y 1 (n) � const × q n . Concerning y 1 , defined through equation (30), we took the values of ρ s c vs and χ s /ρ s c vs published by Petterssen [13]: ρ s c vs � 0.4 Cal/K·cm 2 and χ s / ρ s c vs � 0.01 cm 2 /s, which correspond to y 1 � 16.6 cm. e resulting values of 1 − α(i) and 1 − β(i), where i is the day number, are increasing functions of Y 1 . eir behavior as a function of the i, for the different tested values of Y 1 , is qualitatively the same: β fluctuates around a fixed value, without showing a clear fixed tendency, as shown in Figure 2.
In Figure 2, while α, on the contrary, begins the year with high values, decreasing steadily and reaching a minimum by the middle of the year, when it begins to increase up to the end of the year. If we accept that β fluctuates around a fixed value, then the large values of α obtained in winter are independent of our thermodynamic model developed in previous sections since they were already apparent in Figure 1.
e ratio σ(1 − α)/(1 − α ′ ), of the root mean square deviation of 1 − α to average, ranges from 0.21 to 0.22 when the height Y 1 ranges from 700 to 3000 m. A particularly striking feature is that α ≈ ß for approximately 200 days, a fact consistent with the previous figure where we plotted their ratio c(i). As we analyzed the results corresponding to the list of values for Y 1 , we remarked that neither the distance d(α, β) given by equation (46) nor the mean square deviation σ(α) were sensitive functions to choose a good value of Y 1 . e criterion finally adopted is as follows: the total fraction of energy reflected by the Earth, the albedo, is, in a first approximation, the sum of that reflected by the ground plus that by the clouds. e first one, for a north latitude of 20°according to Kung et al. [31], is α g � 0.158. e clouds contribution is variable and ranges, according to Houghton, from 0.05 to 0.85 [16,33] although he gives average values from 0.3 to 0.7. We have chosen Y 1 so that the smallest albedo is α � 0.158. is choice leads us to Y 1 � 908 m. e following results were obtained with Y 1 � 900 m: 1 − α � � 0.536 and σ(1 − α) � 0.115, while the corresponding values for β are as follows: 1 − β � 0.493 and σ(1 − β) � 0.061.
at the root means square deviation of β is smaller than that of α confirms the lack of a clear tendency of its yearly variation, fluctuating around an average value, as compared with the previously noted yearly behavior of α with next assumptions listed as follows: (1) e thickness of the lower layer of the atmosphere submitted to the daily heating is Y 1 � 900 m. e meaning of this value is as fllows: if the Fourier C 1 coefficient, giving the lowest-order daily variation of temperature is six degrees, that is, if C 1 � 6°C, then the Fourier coefficient C 1 (y), according to equation (36), for the coefficient of daily temperature variation at y � Y 1 m is ≈ C 1 × exp(− 1) 2°.
(2) For the ground, this distance is 16.6 cm, with a similar meaning. (3) e daily energy change of the atmosphere is ♣ (4) e corresponding daily change for the ground is ♣ (5) e largest heat input, possibly due to the wind transport, is approximately ♣ (6) e mean quadratic distance between f(t) and Here, α, the albedo, is connected with the clouds and β, related to the probability of reabsorption, is also connected with the water amount in the atmosphere. Up to this point, we have chosen α and β as constants along the day, but they should depend on the water amount of the atmosphere. e partial pressure p e of the water vapor is, at least in principle, related to the amount of water lying above us. α and β should increase with increasing p e . We are interested now in determining new α ′ and β such that α � α ′ p e (t), β � β ′ p e (t). (51) We can determine p e (t) from the temperature and humidity readings and then fit Fourier polynomials [24][25][26][27][28] and determine again the value of these parameters by least squares. e plot of α ′ and β shows again that α ′ and β take essentially the same values for the summer days but that α ′ takes very small values, and even negative ones, in winter.

Conclusions
is paper presents a tool of Fourier-statistical analysis that has been used in a meteorological behavior for understanding the predicted values of temperature and humidity concerning the hourly data. Although we have defined c as a consequence of our equation (24), we can still keep that definition even if the supposed equilibrium hypothesis is false. However, our results are consistent with χ s /ρ s c vs � 0.0022 cm 2 /s and so they are closer to those reported by Carson [12]. Carson reports 40 different results concerning χ s /ρ s c vs and points that 6 of them are connected with results which deviate from the others. e average of his other 34 measurements is χ s /ρ s c vs � 0.53 cm 2 /s [33][34][35][36][37][38]. It is expected that this work will be useful for people who are interested in studying and researching meteorological temperature and humidity prediction from Fourier-statistical analysis of hourly data. In the performed study, the Fourier analysis can theoretically and practically analyze the data processed, where the length of time-years or just a year's-data would determine the wanted coefficients. However, from just a year's data, we can reach the average trend and isolate the stochastic part of the variables, i.e., the Fourier-statistical analysis is not considered for not being periodic and for the following events: the environmental disturbances induced by men, the natural disasters, the climate change, the ecological processes (the transformation of natural habitats into agricultural and urban land), and biodiversity loss, among others. is study shows that meteorological historical sequences may be employed for the training of neural networks with some predictability capabilities [7][8][9][10]. Furthermore, the key features of meteorological temperature and humidity behavior are investigated by this technique.

Nomenclature
A n : Fourier coefficients B n : Fourier coefficients b − a: Daily hours C j : Fourier coefficients c vs : Constant volume specific heat dp/dy: Barometric equation du: Infinitesimal change thermodynamics e: Elliptical orbit whose eccentricity is along the movement of the Earth f: Frequency F j : Inhomogeneous vector M: Trigonometric polynomial of degree m * : Molecular weight of water N: Total number of data n c : Clear days n h : Half-overcast days n i : Number of measurements made to obtain the experimental pair (t i , T i ) p: Clapeyron's equation p e (t): Fourier expansion of the water-vapor pressure Q: Heat added to the system R: Radiation S co : Solar constant t: Time in seconds T: Ground temperature 〈T〉: Yearly average of temperature T(t, i, l): l-average of temperature at time t for the i-th day u o : Locating the position of the Sun in the sky W: Work done on the system Y n : ickness of the lower atmosphere submitted to the daily temperature variations α: Fraction from the solar radiation which is reflected back by the Earth's atmosphere β: Probability of an emitted photon to be reabsorbed σ: Mean square deviation Φ m (t): Sun's trajectory along the sky is slightly different each 24 hr period θ 0 (t): Zenith angle ρ e : Water-vapor density Advances in Meteorology 11 ρ s : Sample density τ: Time duration of a day ω: Angular frequency σ: Total distance χ: ermal conductivity of the ground.
Data Availability e computer programs and data are available upon request.

Conflicts of Interest
e authors declare that they have no conflicts of interest.