The atmosphere works as a safety barrier, without it; temperature would reach +75˚C during the day and −130˚C below 0 at night, it is made up from a gases mixture that surround the planet with thick concentric layers and variable densities. The atmosphere base is in contact with the land’s surface and its irregularities such as the continent’s protrudes  . The boundary layer or superficial layer is the lowest layer of the atmosphere, where the atmospheric turbulence is created by thermal or mechanic imbalances, produces an increase in air chaotic movements, and it facilitates the diffusion phenomenon, energy and matter transport   .
The atmosphere is not a closed system, it produces exchanges between various elements and has related properties according to their physical state and chemical composition. Internal procedures and external interactions maintain or modify the atmospheric dynamic    . Gases that made up the atmosphere, are well mixed, but it is not physically constant because it has some significant temperature and pressure variations, in relation with the height above sea level  , phenomenon that affects the Ecuadorian Highland, whose height fluctuates between 1000 m a.s.l and 6310 m a.s.l, with an average surrounding 2500 m a.s.l, and it is located in the Equatorial zone, a length around zero degrees.
Most of the atmospheric phenomenon occur in the troposphere, the atmospheric lowest layer, and particularly in the atmospheric boundary layer, superficial troposphere sub layer, where important dynamic processes are carried out. These dynamic processes are given because of the mechanical turbulence development, which is related to the friction between air and the solid surface and the convective turbulence, produced by the heat fluid among soil and air, the result of the soil heating during day hours  . Natural processes that occur in atmosphere are turbulent, which means, a nonlinear system, being a reason for the development several models that describe its dynamic  -  .
The main objective of this research was to study the atmospheric micro dynamic in the Ecuadorian plateau, for establishing the characteristics according to the produced effects because of the height and geographical position in the zone, by the utilization of proposed methods by Van Ulden and Hostlag in 1985  , which are adjusted to air physic properties in the zone, Haro A. et al. (2014), using meteorological data measured at Universidad Nacional de Chimborazo station, located in Riobamba with a geographic position 1˚38' and 1˚4' south latitude, 78˚39' y 78˚40' west longitude, highland at 2754 m a.s.l located in the Ecuadorian plateau.
2. Material and Methods
Using proposed methods by Van Ulden and Hostlag  , a study of superficial fluids has been developed in the boundary atmospheric surface. It is considered that when an air mass circulates over a surface and goes to another one with different characteristics, its physical properties change   .
The balance energy equation in simplified form is given by Equation (1):
where: net radiation flux, sensible heat flux, latent heat flux and stored heat flux in soil, Equation (2).
Where: net radiation fluid , sensitive heat fluid , latent heat fluid and heat fluid stored in the soil , Equation (2).
With , a constant equal to 0.1 in rural areas and 0.3 for urban areas, proposed by Doll D., Ching J. K. S. and Kaneshire J.  .
Sensitive heat is given by Equation (3)  .
where is a constant equal to 20 W∙m−2 and the parameter S is defined by Equation (4).
Net radiation can be parameterized according Equation (5), based on albedo A, solar elevation angle θ, cloudiness fraction N, flux radiation incident Q and air temperature T.
s = 5.67 × 10−8 W∙m−2∙K−4 it is the Stefan-Boltzmann constant.
c1 = 5.31 × 10−13 W∙m−2∙K−6 y c2 = 60 W∙m−2.
With h humidity, the c3 parameter is given by the formula:
a1 = 990 W∙m−2
a2 = −30 W∙m−2
b1 = 0.75
b2 = 3.4
With solar elevation, , it can be calculated according to latitude function, and longitude to the studied zone and the Julian day  .
The cloudiness is approximately determined with Equation (13).
The kinetic energy for the atmospheric turbulent is given by the Equation (9)    .
The left term represents the total energy fluid, while the first right term represents the thermal energy fluctuation, the second term represents the mechanical energy fluctuation, the third one the energy dissipation because of the friction effect and the fourth one, the term of the energy fluid that mixes the transport and pressure transport.
Obukhov’ lenght (L) is determined by Equation (9), where thermal and mechanical energy are balanced.
According to the similarity theory of the low atmosphere, , where friction velocity, with U(z) wind speed to a height z and z0 of roughness height that can be found in charts  , and K the constant Van Karmans, Cp specific air heat and ρ air density that adjusts to the height through Equation (12)  . Replacing in Equation (6) and considering Z = L, the Obukhov length, it is described as:
The parameter L is considered surface and can be calculated with meteorological measurements (wind speed, temperature, pressure, net radiation, humidity) equation that is used as an indicator of the dynamic state of the atmosphere  .
The dynamic state is evident when a plot of land is moved upward, if it does not suffer an alteration with height, it means, if the external forces of the surrounding plots of land do not act, so; the atmosphere is considered as neutral or adiabatic; if the forces act and the small air portion is moved upward, it is said that is non-stable, but if the external forces tend to make the plot of land to go back to the initial position, then the atmosphere is in stable conditions.
Quantitatively, stability can be measured with Obukhov’s length.
If L < 0 is unstable
If L > 0 is stable
If L ± ∞ is neutral or stable
Depending on the fluctuation degree of mechanic or thermic energy, L is negative when the thermic energy is positive (soil transfers heat to atmosphere) and positive when the opposite occurs, when this fluctuations ten to give a L large number (positive o negative) it is said that the is neutral or adiabatic   .
3. Results and Discussion Results and Discussion
This research seeks characterize the micro-meteorological parameters variation, as a function of physical and meteorological parameters, according to Ecuadorian Highlands regions. For this, a software in MATLAB has been developed, based on the purpose of Van Ulden and Hostlang, which determines some parameters such as: Obukhov’s length, superficial heat fluids, sensitive and latent, friction speed and net radiation, which characterize the micro dynamic boundary atmospheric zone layer.
Table 1 presents annual averages, minimums, and maximums, of the meteorological variables measured in the research zone, where the height effects in temperature and pressure can be noticed, they decrease considerably according to the sea level, in the position that is the research zone (zero latitude). The results indicate strong radiation intensities with low air temperatures, which will produce temperature changes, specifically during day hours.
Table 2 presents the average, minimum, and maximum from the micrometeorological variables studied, where small variations where observed among latent heat fluid, sensitive and superficial, with a major fluctuation in the sensitive heat fluid, which varies between −27.19 to 359,76 W/m2, the same occurs with Obukhov length which takes values among −6732.61 to 31020.98 m, which indicates the strong energy fluctuations that occur between Earth and atmosphere, variables such as net radiation suffer a decrease in the found values regarding to the incident radiation due to the reflected radiation and the friction speed is characterized by maintaining values that fluctuate a small range from 0 until 1.10 m/s.
Figure 1 represents the variation of the thermal energy and net average radiation, an hour from a year (2015) they have been determined according to the researches by Van Ulden y Hostlag, with an important presence of sensitive and superficial heat fluid, due to the height and geographical position, a strong interaction between earth surface and superficial atmospheric air is produced, which is represented in the high values of the sensitive energy.
Figure 2 represents the speed variation of the hour average friction, in the
Table 1. Average, maximum and minimum methereological parameters from the year 2015 in Riobamba city.
Data weather station of the air quality of the UNACH project.
Table 2. Average, maximum and minimum micrometeorological parameters from the year 2015 from Riobamba city.
Figure 1. Sensitive heat fluids, latent heat, surface heat and net radiation, average per hour of the year 2015.
Figure 2. Average friction speed hour of the year 2015.
research year, with a tendency to ascend at midday and to decrease at night, being a regular behavior through the year due to the geographical position.
Figure 3 reflects the Obukhov length performance, an average hour in the research year, which presents a regular performance in different hours, with a strong variation between 17:00 and 19:00 during the sunset, which produces a thermal reversal, modifying the dynamic atmospheric state.
The variation of the regular meteorological parameters is determined during the year, with significant variations through the day, whose height and geographical
Figure 3. Length of Obukov average hour of year 2015.
position effects are specifically noticed in the temperature with a 13.55˚C average and a pressure with a 731.54 mB average, these are small according to the results above sea level a.s.l in those latitudes, where micrometeorological variables calculated with the model, present regular performances throughout the year, similar to the meteorological conditions, variations are presented mostly during the day, with significant variations in thermal energy fluids, notice the height effects and geographical position with a high sensitive heat fluid. Besides, Obukhov’s length presents a regular performance through the year with significant variations between 17:00 and 19:00, where clearly a thermal reversal is presented which changes the atmospheric dynamic state.
 Doll, D., Ching, J.K.S. and Kaneshire, J. (1985) Parametrization of Surfaces Heating for Soli and Concrete Using Net Radiation Data. Boundary Layer Meteorology, 32, 351-372.
 Cohen, A., et al. (2015) A Review of Planetary Boundary Layer Parameterization Schemes and Their Sensitivity in Simulating Southeastern U.S. Cold Season Severe Weather Environments. Weather and Forecasting, 30, 591-612.
 Businger, J.A., et al. (1971) Flux-Profile Relationships in the Atmospheric Surface Layer. Journal of the Atmospheric Sciences, 28, 181-189.
 Rodríguez, D., et al. (2015) Variación de la estabilidad y altura de la capa de mezcla en la ciudad de pinar del río: Su relación con condiciones sinópticas. [Stability and Change in Mixing in City of Pinar Del Rio: It’s Relationship with Synoptic Conditions] Revista Brasileira de Meteorología, 30, 1-15.
 Marks, D., et al. (2007) Comparing Simulated and Measured Sensible and Latent Heat Fluxes over Snow under a Pine Canopy to Improve an Energy Balance Snowmelt Model. Journal of Hydrometeorology, 9, 1506-1522.