1. Introduction
China’s topography is high in the west and low in the east, presenting a step-like distribution, and the typical terrain in the central and western regions are mountains and plateau basins. According to the internationally accepted altitude classification standard above the sea level: 1500–3500 m is high altitude; 3500–5500 m is ultra-high altitude; and above 5500 m is extremely high altitude [
1]. Most of the mountainous areas in central and western China belong to high-altitude areas where wind characteristics of the surface layer are quite different from the vast flat areas in eastern China [
2]. It is difficult to observe and evaluate the wind power in these mountainous areas, both for technology and for work processing. Therefore, in order to explore the influence of complex terrain on wind characteristics of the surface layer and to better develop and utilize wind energy resources of high-altitude regions in central and western China, it is necessary to carry out research to discover the wind characteristics of the surface layer over some typical terrains.
The stability characteristics of the atmospheric boundary layer have been focused on for a long time by scholars all over the world and studies on the algorithm of atmospheric stability have been reported [
3,
4,
5,
6,
7]. Jin et al. [
8] analyzed the applicability of different atmospheric stability classification methods by using observational data from four 100 mmeteorological towers in Urumqi. It has been found that the temperature-difference-wind-speed method is the most suitable computing method for atmospheric stability–after comparing the temperature-difference method, temperature-difference-wind-speed method, wind-speed-ratio method, Richardson number method, and the overall Richardson number method. Ameya et al. [
9] compared the diurnal and seasonal variation characteristics of atmospheric stability by calculating the Obukhov length based on the observation data of meteorological masts from two offshore wind farms in the Dutch North Sea. Ameya et al. [
10] compared the mean wind speed and atmospheric stability distribution by calculating the Obukhov length based on four test stations: Høvsøre, Egmond aan Zee-OWEZ, Östergarnsholm, and Hurghada. Gryning et al. [
11] illustrated the frequency distribution of atmospheric stability for rural, residential, and urban homogeneous areas by employing the data measured from three-dimensional sonic anemometers in Denmark and Germany. Ameya et al. [
12] calculated the frequency occurrence of the wind speed in different atmospheric stability at two offshore Danish sites to estimate the fatigue damage of wind turbines.
In addition, the study of the eddy covariance fluxes is a popular field for atmospheric boundary layer research [
13,
14,
15,
16,
17,
18,
19,
20]. Chen et al. [
21] compared the diurnal and seasonal variations of heat and momentum flux by using the turbulent data from January 2007 to December 2008, which were collected from the Semi-arid Climate & Environment Observatory of Lanzhou University (SACOL) and the Tongyu Observatory of Jilin. He et al. [
22] revealed the impact of pollutant concentration on turbulent fluxes by studying the diurnal variation of fluxes during a heavy PM2.5 pollution period based on data from a meteorological gradient tower in Beijing. Wang et al. [
23] penetrated the basic characteristics of the land surface energy budget, seasonal and diurnal variations of moisture, and heat flux over the Loess Plateau by employing the field observations from SACOL. Sun et al. [
24] studied the seasonal and daily variations of sensible heat flux and latent heat flux over a lake in the Badain Jaran desert to explore surface energy exchange between the atmosphere and water bodies.
Different studies have been published on the atmospheric stability and the variation characteristics of eddy covariance fluxes. However, little attention has been paid to the influence of different complex underlying surfaces on wind characteristics in the surface layer of high-altitude regions in central and western China. To fill this gap, the diurnal and monthly variations of atmospheric stability were contrasted in two typical topographies, the Qiaodi Village in Sichuan (in western China, site 1) and the Nanhua Mountain in Shanxi (in central China, site 2), according to the Obukhov length calculated by the eddy covariance data. The energy exchange process between complex underlying surfaces and the atmospheric boundary layer can be reflected to a certain extent. The references for the wind resource assessment of high-altitude areas in central and western China can be provided by investigating the diurnal variation differences of the turbulent fluxes at the two sites.
2. Field Test
Site 1, located in the alpine mixed forests of the western Sichuan Plateau, is surrounded by relatively steep hillsides (102.087° E, 27.628° N, 2786 m above sea level). The underlying surface is mainly covered by shrubs, weeds, and pines about 5 m tall (
Figure 1a). Site 2 is located on a gentle slope near the ridge (111.991° E, 39.272° N, 1800 m above sea level, 500 m above the foothill village). The local terrain is relatively flat, which is a typical Loess Plateau landform. The underlying surface is covered by sparse trees and weeds usually less than 30 cm (
Figure 1b).
Figure 2 shows the eddy covariance data at the two sites.
The local land use categories of the Weather Research and Forecasting (WRF) Model near the two sites are shown in
Figure 3. The land use categories are shown in
Table 1. The land use description of the mixed forest near site 1 is completely different from the land use description of the grassland near site 2.
The data of the two sites are measured by the eddy covariance system at 10 m above the ground, consisting of a three-axis Sonic Anemometer (CSAT3B, Campbell Scientific Inc., Logan, UT, US) with the sampling frequency of 10 Hz. The data files with a 0.1 s time resolution contain three wind components for the , , and directions, air temperature, CO2 concentration, water vapor concentration, and air pressure.
The observation period at site 1 is from 9 January to 20 March 2020, while the observation period at site 2 is from 9 January to 20 March 2021.
4. Results and Discussion
To study the diurnal variation characteristics of atmospheric stability on the complex terrain of the Yalong River Basin,
is calculated as shown in Equation (5) based on the observed data from site 1 (
Figure 4). The atmospheric stability categories are shown in
Table 3. As shown in
Figure 4, the relative proportion of the neutral and the unstable boundary layer changes significantly at 8:00 (Beijing Time, BJT, the same as below) and 18:00. The dominant boundary layer transforms from the neutral boundary layer to the unstable boundary layer at 8:00, and it is the opposite at 18:00. From 18:00 to 8:00 the next day, the boundary layer is dominated by a neutral boundary layer with a certain proportion of a weak-stable and stable boundary layer. From 8:00 to 18:00, the boundary layer is dominated by a weak-unstable and unstable boundary layer. The proportion of the unstable boundary layer is more than 90% and reaches a peak of the day at 12:00.
With the same method as site 1, the diurnal variation characteristics of atmospheric stability on the complex terrain of the Loess Plateau is plotted (
Figure 5). It can be seen that the relative proportion of the stable and unstable boundary layer changes significantly at 7:00 and 18:00. The dominant boundary layer transforms from the stable boundary layer to the unstable boundary layer at 7:00 (1 h earlier than that at site 1), and it is the opposite at 18:00. From 18:00 to 7:00 the next day, the boundary layer is dominated by the stable boundary layer, with a certain proportion of the weak-stable boundary layer. From 7:00 to 18:00, the unstable boundary layer is dominated with a certain proportion of the weak-unstable boundary layer. Comparing the changes of the atmospheric boundary layer at the two sites, it can be derived that the nighttime dominant boundary layer at site 1 is the neutral boundary layer, while at site 2 it is the stable boundary layer.
Correspondingly, the distribution of atmospheric stability in different periods at site 1 is analyzed in
Table 4. During the study period, the surface layer is dominated by the neutral boundary layer, with a proportion of more than 50%, followed by the unstable boundary layer, accounting for more than 30%. The occurrence probability of the stable boundary layer is the lowest, only about 10%. At site 2, the stable boundary layer is the main, accounting for nearly 40%, followed by the unstable boundary layer, accounting for about 30%. The proportion of the weak-unstable boundary layer is the lowest, less than 6% (
Table 5). Compared with previous results, the diurnal variation characteristics and the distribution of atmospheric stability at the two sites are quite different from those in Xilinhot, Jiangxi, and the Netherlands [
25]. This may be caused by the differences in the topography and the underlying surface between the two sites. Site 1 is located in the hinterland of the western Sichuan Plateau and on the east bank of the Yalong River. It is a complex mountainous terrain of the plateau and the underlying surface is mixed forest. The topographic effect has an impact on the anemometer tower, thus, the neutral boundary layer dominated during the nighttime at site 1. Site 2 is located on the Loess Plateau. Although the overall altitude is high, the local underlying surface is a flat and gentle slope, thus the diurnal variation characteristics of the atmospheric stability at site 2 is essentially consistent with the classical similarity theory.
Accordingly, the average horizontal wind speed at 10 m above the ground at different atmospheric stabilities during each period at the two sites can be calculated (
Table 5 and
Table 6). Compared with site 2, the horizontal wind speed at site 1 is significantly smaller at the same atmospheric stability. As evident from
Table 6, the horizontal wind speed in the neutral (unstable) boundary layer is significantly higher (lower) than that at the other atmospheric stability at site 1. As shown in
Table 7, the horizontal wind speed in the neutral and weak-unstable boundary layers are clearly higher than that in other boundary layers, and in the stable boundary layer it is the smallest. This result is consistent with the opinion of Gong et al. [
25], and it is because wind shear is the main source of turbulent energy in the neutral boundary layer and buoyancy is negligible. While in the non-neutral boundary layer, the dynamic and thermodynamic factors act simultaneously. However, the difference of the stability variations and the average horizontal wind speed at 10 m above the ground of the two sites are likely associated with the discrepancy of the momentum and heat flux transfer caused by the diverse underlying surfaces. Therefore, the variations of the momentum and heat flux at the two sites are investigated.
The greater the absolute value of the momentum flux is, the greater the kinetic energy in the surface layer is, and the stronger the vertical upward transport of the turbulence is.
Figure 6 shows the variation of the momentum flux at the two sites. On the whole, the diurnal variations of the momentum flux at the two sites are all unimodal. In different boundary layers, the momentum flux is the largest in the neutral, followed by in the weak-unstable, and the smallest in the stable, which is in accordance with the conclusion obtained above. Compared with site 2, the average momentum flux at site 1 is evidently larger, and the diurnal range of the momentum flux at site 1 (0.379 N.m
−2) is greater than that at site 2 (0.297 N.m
−2). Compared with the values of Yuzhong in Gansu (0.08 N.m
−2) and Tongyu in Jilin (0.25 N.m
−2) [
21], the diurnal range of the momentum flux at site 1 is clearly higher, which may be mainly caused by the complex underlying surface (alpine mixed forests of the western Sichuan Plateau), and at site 2 it is similar with that in Tongyu.
Figure 7 shows the diurnal variation of the sensible heat flux at the two sites. The diurnal variations of the sensible heat flux at the two sites are all unimodal and values of the sensible heat flux change remarkably at 8:00 and 18:00. Due to enhanced solar radiation at 8:00, a sharp vertical temperature gradient is generated in the surface layer to promote the development of turbulence. At this time, the sensible heat flux begins to change from the negative to the positive, implying that the ground cooling effect during the nighttime converts to the ground heating effect. After that, the sensible heat flux continues to increase and reaches a peak value between 13:00 and 15:00. Hereafter, the sensible heat flux gradually decreases. Between 18:00 and 19:00, the sensible heat flux turns to negative again. It means that the sun sets at this time and the solar radiation received by the surface layer is weakened, thus the vertical temperature gradient is reduced and the energy provided by thermal driving for the development of turbulence is reduced. From 19:00 to 7:00 the next day, the sensible heat flux remains near zero. In the different stability, values of the sensible heat flux are in the order: the weak-unstable boundary layer > the unstable boundary layer > the neutral boundary layer > the stable boundary layer > the weak-stable boundary layer, illustrating the dominant role of thermal driving in the unstable boundary layer. Compared with site 2, the peak time at site 1 (14:00) is 1 h behind site 2 (13:00). The peak value of site 1 (446.239 W.m
−2) is two and a half times as much as that of site 2 (180.035 W.m
−2), and both are higher than those of Yuzhong, Gansu (140.3 W.m
−2) and Jilin Tongyu (157.3 W.m
−2) [
21]. It is mainly because the altitude of site 1 is higher than that of site 2, there is a difference in solar radiation intensity, resulting in a large discrepancy in the sensible heat flux between the two sites.
Figure 8 shows the diurnal variation of the latent heat flux in the surface layer at the two sites. The diurnal variations of the latent heat flux at the two sites are all unimodal. In a different stability, the latent heat flux in the unstable and weak-unstable boundary layers are the largest, followed by the neutral boundary layer, and in the stable and weak-stable boundary layers are the smallest, maintaining below 10 W.m
−2. Comparing the diurnal variation of the latent heat flux at the two sites, it can be seen that the peak value at site 1 (79.902 W.m
−2) appears at 14:00, while the peak value at site 2 (57.082 W.m
−2) appears at 13:00. Compared with those values of Gansu Yuzhong in the summer (195.4 W.m
−2) and Jilin Tongyu in the summer (211.4 W.m
−2) [
21], the peak values of the latent heat flux at the two sites are significantly low. This is because our study periods of the two sites are relatively dry, resulting in a lower latent heat flux than Yuzhong and Tongyu.