Surface Layer Turbulent Characteristics over the Complex Terrain of the Loess Plateau Semiarid Region

Surface layer turbulence has an important influence on land-air interactions and pollutant dispersion, and studying the characteristics of surface layer turbulence in complex terrain can contribute to understanding land-air interactions, improving model surface layer parameterization, and enhancing pollution prediction capabilities. (e surface layer turbulence observations from the Semi-Arid Climate and Environment Observatory (SACOL) of Lanzhou University in 2008 were processed in this study. High-quality turbulence parameters were calculated, the statistical and transfer characteristics of turbulence were analyzed, and the formation of turbulence was assessed in terms of dynamics and thermodynamics. (e atmospheric stability in the semiarid region of the Loess Plateau is basically dominated by neutral/near-neutral and weakly unstable/weakly stable conditions; this pattern is significantly different from the preponderance of unstable and stable layers at other sites.(e turbulence intensity differs significantly in both the horizontal and vertical directions and basically shows the relationship Iu ≈ Iv > Iw. (e mean values of Iu and Iv are 0.42 and 0.40, respectively, and the mean value of Iw is 0.14, which is different from the general pattern of Iu > Iv > Iw, indicating that the turbulence intensity at SACOL is characterized by a large lateral wind contribution. (e dimensionless standard deviations of the nonneutral wind velocity components satisfy the “1/3rd law,” and the dimensionless standard deviations of u, v, and w components are 3.35, 2.98, and 1.26, respectively, in the semiarid Loess Plateau.(ese values are larger than those over flat terrain.(e contribution of the shear term to the formation of turbulence is greater than that of the buoyancy term, and the mean annual values of the shear term and the buoyancy in the kinetic energy budget equation in SACOL are 47.94×10 m·s 3 and 11.32×10 m·s , respectively.(e annual mean values of themomentum transfer coefficientCD and the heat transfer coefficient CH under near-neutral conditions are 8.54 × 10 − 3 and 2.52×10 , respectively.


Introduction
e study of the atmospheric boundary layer is a major focus in Earth system science, as the exchange of materials and energy within the layer has a significant impact on weather and climate processes [1]. e boundary layer is also a major conduit for ecosystem feedback to the atmosphere and is the region where many pollutants accumulate, including aerosols, with aerosol distribution and pollutant dispersion being directly influenced by turbulent activity within the boundary layer [2]. Sudden meteorological hazards such as dense fog and dust are also closely related to boundary layer processes. Due to surface forcing processes, such as frictional dragging, heat transfer, and evapotranspiration, turbulence is the dominant form of movement in the surface layer atmosphere and serves as a bridge and channel for the exchange of energy and materials between the surface and the atmosphere, such as heat and moisture. erefore, an in-depth study of the turbulence characteristics in the atmospheric boundary layer, especially in the surface layer, is important for understanding land-air interactions, improving the model's performance over complex terrain, improving air pollution prediction, and the efficient use of wind energy.
At present, a great deal of research has been carried out on the near-surface formation problem over horizontal uniform flat underlying surfaces [3][4][5][6][7]. e Monin-Obukhov similarity theory (MOST) has been validated and largely agreed upon in unstable conditions and flat terrain areas. However, the understanding of complex surface turbulence is still poor, and the need to characterize surface layer turbulence over complex terrain is becoming more urgent with the development of higher-resolution weather models and the increased accuracy of simulation results for boundary layer parameterization and turbulence closure schemes. Rotach et al. [8] and Fernando et al. [9] studied land-atmosphere interactions in complex terrain and pointed out that secondary circulation, such as valley winds, causes spatial inhomogeneity in turbulence. Hooke et al. [10] verified the applicability of MOST to complex terrain, noting that it is only applicable to horizontal wind velocities in a specific stability interval (− 0.5 to 0). Moraes et al. [11] stated that the normalized standard deviation of horizontal wind velocity increases with increasing stability, mainly due to the fact that u * decreases rapidly with increasing stability. For complex terrain, Martins et al. [12] observed the standard deviations of velocities and temperature fluctuations at the edge of a sharp cliff, suggesting that complex topography causes quasimesoscale motions, such as gravity waves and local circulation, which change the structure of the turbulence. de Franceschi et al. [13] fitted the empirical coefficient with along-valley and cross-valley winds for u, v, w, and θ, respectively, and confirmed the validity of MOST over both stability ranges and for all wind velocity components. Babic et al. [14] used the stationary datasets in Owens valley, in which nonstationarity in turbulence was removed to investigate how the stationarity criterion affects the fluxvariance similarity functions and systematic differences from previously near-neutral values of the parameters in the flux-variance similarity functions over flat terrain were found, indicating a larger anisotropy of the flow over complex than over flat terrain. Nadeau et al. [15] found for a steep slope that the scaling was successful for the fluxvariance similarity functions, but not for the flux-gradient similarity functions, especially during stable conditions. Disparities such as those presented above call for further measurement campaigns to be conducted over complex terrain, mainly to obtain a more definitive and conclusive range of the applicability of similarity theory. So, study on characteristics of surface layer turbulence in complex terrain as the Loess Plateau, the similarity of variance relationships in typical complex terrain, the thermal and dynamic contributions of turbulent kinetic energy, and transfer characteristics is demanded. e semiarid zone of the Loess Plateau, as a globally important and unique geographic region, is widely distributed and sensitive to climate change, making it a crucial area in the study of global change [16]. e eddy-correlation system of the Semi-Arid Climate and Environment Observatory of Lanzhou University (SACOL) was calibrated in 2008 with good continuity and quality of data. A series of previous pieces of literature have analyzed the turbulence and land-atmosphere interaction at SACOL, with this data in 2008, so the data quality is recognized as credible. For instance, Guanghul et al. [17] compared the effects of different processing methods (including the double rotation (DR), planar fit (PF), and fetch planar fit (FPF)) on the results and invested the sensitivity of their results to the processing, suggesting that DR is recommended in the complicated terrain for reducing calculation and improving the data quality. Liang et al. [18][19][20] explored the effects of low-level jet and mesoscale motion on turbulence, and they also discussed the applicability of similarity theory to SACOL, the energy closure problem, and the influence of topography on turbulence. However, further research is needed on the similarity in the variance characteristics of turbulence influenced by the complex topography of the Loess Plateau region and the dynamics and thermal contributions of turbulence formation. A systematic understanding of surface layer turbulence in the Loess Plateau region will help to provide a parameterization of the surface layer material and energy transfer processes on the complex underlying surface of the Loess Plateau region.
In this paper, the SACOL 2008 surface layer turbulence observation data are quality controlled and revised to calculate the quality of high-quality turbulence parameters (including atmospheric stability, turbulence intensity, dimensionless standard deviation of wind velocity components, turbulence momentum, sensible heat and latent heat flux, buoyancy and shear terms in the turbulence kinetic energy budget equation, and turbulence bulk transfer coefficient). en, the surface layer turbulence characteristics over complex terrain are analyzed, and the dynamic and thermodynamic function contributions to turbulence formation are evaluated. is provides a database for land-air interaction studies and model parameterization.

Observation Sites.
e Semi-Arid Climate and Environment Observatory of Lanzhou University (SACOL) is approximately 48 km southeast of downtown Lanzhou, China, on top of the Cuiying Mountains (35°57′46″N, 104°8′13″E) at an elevation of 1965.8 meters above sea level. On a large scale, the Cuiying Mountains lie in a southeastnorthwest-oriented basin with an average elevation of approximately 1750 m at the base of the basin and a relative height of approximately 200 m. On a smaller scale, SACOL is located on a loess mesa surrounded by typical landforms of the Loess Plateau, such as intricate gullies, beams, and ridges [19]. e surface cover is largely native vegetation [16], the flux observation field is largely flat (extending over 1000 m in the north-south direction and approximately 200 m in the east-west direction), the average plant height is 0.2 m, the averaged aerodynamic roughness length z 0 is 0.04 m, and the displacement distance is d � 0.13 m. SACOL is a good representation of the semiarid zone of the Loess Plateau. Since its establishment, SACOL has provided a large amount of valuable information for climate change research in the semiarid zone of the Loess Plateau and has filled the gap in long-term observations of land-air exchange in the semiarid region of the Loess Plateau summer. For more information about SACOL, please refer to Huang et al. [16]. e surface radiation observation system is installed at a height of 1.5 m and consists of a total diurnal radiometer (CM21, Kipp & Zonen) and a surface radiometer (CG4, Kipp & Zonen) with observation values for total solar radiation, reflected solar radiation, atmospheric counterradiation, and surface longwave radiation. e main observations, instrument types, and installation locations are listed in Table 1.
To ensure the accuracy and validity of the results, the EC data at 10 Hz are processed according to the recommendations of the Global Flux Observation Network (FLUX-NET) by eliminating the noise signal (diagnostic value diag_csat ≠ 0 for CSAT3) caused by inherent conditions, such as instrument power failure and malfunction, and environmental factors, such as rain, dust, and flying insects. en, despiking [22], DR of coordinates [4], humidity correction for ultrasound virtual temperature [23,24], and density correction [25] were performed in an orderly manner. Pulsation, variance, and covariance were calculated on a 30 minute averaging period for the quality-controlled data. It is worth mentioning that Beijing time (BJT) was used for all times; BJT is 8 h earlier than universal time and 1 h earlier than the local time of SACOL. e seasons are divided into spring (March-May), summer (June-August), autumn (September-November), and winter (December-February). e hours are divided into 07:30-19:00 for the day and 19: 30-07:00 for the night. e gradient term (zu/zz), required to calculate the tangential term in the turbulence kinetic energy budget equation, is calculated from the wind velocity at the heights of the boundary layer gradient towers of 2 m and 4 m. In addition, a reference height of 4 m is used to calculate the wind velocity and temperature difference between the reference height and the ground surface when calculating the bulk transfer coefficient.
Surface radiation data are used to calculate surface temperatures: ULR is surface upward longwave radiation, DLR is atmospheric downward longwave radiation, ε is the surfacespecific emissivity (ε � 0.95), and σ is the Stefan-Boltzmann

EC Method.
According to the Reynolds decomposition, the physical quantity s can be decomposed in the average duration t (we use 30 minutes here) into the sum of the average quantity s and the pulsating quantity, and the turbulent pulsating quantity can be obtained by removing the average quantity from the measured instantaneous value as follows: e variance s ′ s ′ and the covariance of the vertical wind speed and the physical quantities w ′ s ′ are averaged over the mean duration t to obtain the mean variance s ′ s ′ and the mean covariance w ′ s ′ ; thus, information on turbulence intensity and flux is demonstrated. e momentum flux τ, the sensible heat flux Hs, and the latent heat flux λE of turbulence in the vertical direction are where ρ a is the air density (kg·m − 3 ); w and u are the vertical and horizontal components of the wind velocity (m·s − 1 ), respectively; ρ v is the water vapor density (kg·kg − 1 ); E is the water vapor flux (kg·m − 2 ·s − 1 ); λ is the latent heat of vapor evaporation (J·kg − 1 ); and C p is the specific heat of air at constant pressure (J·kg − 1 K − 1 ).

Atmospheric Stability Parameter.
According to MOST, the properly scaled dimensionless characteristics of the turbulence at reference height z m are universal functions of a stability parameter, ζ � z m − d/L, where d is the displacement distance. It is generally considered that the displacement distance of the underlying surface under uniform dense vegetation coverage is d � 2/3h, h is the mean height of vegetation, and L is the Obukhov length:

Advances in Meteorology
In equation (4), is the friction velocity, κ is the von Karman constant (whose value interval is usually 0.35∼0.40, and we adopt 0.40 here), and g is the acceleration due to gravity. Figure 2 shows the distribution of atmospheric stability for SACOL 2008 throughout the year and shows that the surface layer atmospheric stability is dominated by neutral/ near-neutral and weakly unstable/stable conditions, with few cases of strong stability and strong instability. is differs from the classical stability characteristics over flat terrain [26], which is usually dominated by unstable and stable conditions. is is due to the fact that, under complex topographic conditions, there is strong turbulence-related mixing in the surface layers, which makes it difficult to maintain strong stratification, such as strong stability and strong instability, which has been confirmed by multiple observation experiments [27,28]. e daytime stratification is dominated by instability, with a mean value of − 0.27, while the nighttime stratification is dominated by stability, with a mean value of 0.19. Considering that there are significant seasonal and diurnal variations in atmospheric stability, the daytime and nighttime distributions of stability in different seasons are presented in Table 2. e surface layer stratification was mainly unstable in the daytime and mainly stable at night. e average stability in winter was 0.03, which was significantly higher than that in other seasons.

Turbulence Intensity.
Turbulence intensity, which is the ratio of the standard deviation of wind velocity to the mean wind velocity, is a commonly used measure of turbulence strength and is an important physical quantity for characterizing atmospheric turbulence motion: where I x is the turbulence intensity in three directions, x ′ 2 is the standard deviation of the corresponding wind velocity, and U � ������ is the average horizontal wind velocity over average duration t. Figure 3 shows the average diurnal variation in turbulence intensity, wind velocity, and stability throughout the year. e diurnal variation in turbulence intensity is clearly characterized by a single-peak pattern, with  e average wind velocity and stability show opposite trends to the intensity of turbulence. At midday, when the average wind velocity is minimal, the atmospheric stratification is predominantly unstable, and turbulence develops most intensively, at night, when the wind velocity is higher, the boundary layer is more clearly stratified, stable stratification predominates, turbulence is inhibited by the stratification, and the turbulence intensity is limited to a low level. e turbulence intensity varies considerably in both horizontal and vertical directions, basically as I u ≈ I v > I w , which is different from the general characteristic of turbulence intensity over most of the other land surfaces, such as grassland of the central United States [4], which always exhibit the pattern I u > I v > I w . e annual averages of I u , I v , and I w are 0.42, 0.40, and 0.14, respectively. e complex topographic relief in SACOL results in a local wind field that is subject to a variety of motions, such as secondary circulation, quasisteady motion, and other changes; consequently, the wind direction is complex and variable, with turbulent lateral wind contributing significantly.
Atmospheric stratification is divided into unstable (ζ < 0), stable (ζ > 0), and near-neutral (− 0.1 ≤ ζ ≤ 0.1) conditions to further investigate the turbulence characteristics under different stratification conditions. Table 3 shows the distribution of turbulence intensity in different seasons and under different stratification conditions. e turbulence intensity shows significant differences between different stability levels within the same season, and I u , I v , and I w decrease uniformly with increasing stability. It is worth noting that the season-related differences are far less than the stability-related differences.

Normalized Standard Deviation of the Velocity
Fluctuations. According to the Monin-Obukhov similarity theory (MOST), the dimensionless standard deviation of wind velocity in surface layer turbulence is a universal function of ζ [29,30], and there is a uniform form between them: where σ j is the standard deviation of the variable j, u * is the frictional velocity, ϕ j is a universal function, and C x1 and C x2 are the empirical coefficients. Under neutral conditions, ζ � 0, where turbulence is purely mechanical, the dimensionless standard deviations of the wind velocity components u, v, and w will be constants, independent of roughness and height [31]. e similarity of the variance relationship in the form of equation (6) is referred to below as the "1/3rd law." e empirical relations (equations (7) and (8)) summarized by Panofskyet al. [32] for flat and homogeneous surfaces are the most widely used but have different statistical results for many different surfaces and

Near-Neutral Stratification
(|ζ| ≤ 0.1). e dimensionless standard deviations of the three-dimensional velocity fluctuations σ u /u * , σ v /u * , and σ w /u * in SACOL are 3.35, 2.98, and 1.26, respectively, which are close to the results for an area of undulating terrain in a Florida state park [35] (3.20, 2.90, and 1.24) and greater than the average results (2.39 ± 0.03, 1.92 ± 0.05, and 1.25 ± 0.03) for several sites with flat terrain, including Roskilde (Denmark), Ladner (Canada), and Donaldson (Minnesota) [35]. e dimensionless standard deviation of the horizontal wind velocity is larger than that of the vertical wind velocity. σ w /u * in complex terrain differs little from that in flat terrain, while σ u /u * and σ v /u * are significantly larger than those in flat terrain. is is because horizontal fluctuations are mainly generated by relatively large quasihorizontal turbulences, which typically have diameters of several hundred meters or more and are less adaptable to the terrain than vertical pulsations, which are dominated by the contribution of relatively small turbulences [31], so the influence of topography on turbulence is more pronounced in the horizontal direction.

Nonneutral Stratification (Stable (ζ > 0) and Unstable
Conditions (ζ < 0)). Although the 1/3rd law of dimensionless standard deviation of wind velocity under unstable conditions has been widely demonstrated, it is debatable whether the dimensionless standard deviation of wind velocity under stable stratification is satisfied. Figure 4 shows the relationship between the dimensionless standard deviations of the SACOL wind velocity components u, v, and w with respect to the stability parameter ζ. e dimensionless standard deviation of the wind velocity in all three directions satisfies the 1/3rd law reasonably well when unstable, in the form of formulas (9)- (11). e standard deviation of the dimensionless wind velocity increasing with ζ in the stable boundary layer has been observed at many sites and has been the focus of discussion in this field of study. On the one  hand, under stable conditions, σ j /u * , (j � u, v, w) and u * are both small, so there is large uncertainty in the ratio between them. Additionally, there is also a serious uncertainty in u * itself, and a tendency to overestimate σ j /u * due to underestimation of u * can easily be observed [31]. In recent years, an increasing number of studies have found that, in a stable boundary layer, turbulent motion is influenced by a variety of unsteady motions that are significantly intermittent, leading to an enhancement in turbulent motion over short periods of time, resulting in an increase in σ j /u * [20,36]. Furthermore, considering that u * occurs in both ζ and σ j /u * , ζ is proportional to 1/u 3 * , and σ j /u * is proportional to 1/u * , resulting in a spurious 1/3rd order relationship, which is generally called the self-correlation relationship [37]. It has been shown that this self-correlation is more pronounced in the strong stable case (such as ζ > 1～2) and when the data are discrete [38,39].

Turbulent Momentum Flux and Heat
Flux. e transfer of momentum and heat between the Earth's surface and the atmosphere is a driver of climate from local to global levels and has an important influence on climate change [40]. e semiarid zone of the Loess Plateau has a unique terrain and vegetation type; in this case, a correct understanding of the characteristics of turbulent fluxes, such as momentum and heat, is useful for understanding the earth-air interaction processes in this region. Figure 5 displays the diurnal variation in the SACOL turbulent momentum. e momentum flux rises gradually during the day to above 0.08 kg·m − 1 ·s − 2 in the afternoon and decreases at night to a minimum of 0.03 kg·m − 1 ·s − 2 at approximately 6:00 p.m. Depending on the season, the momentum flux is more similar in spring and summer than in autumn and winter due to seasonal differences in the stability of the surface layer (Table 2). Due to surface-cooling effects, the nighttime stratification is usually stable, which has a restraining effect on turbulence generation, the momentum transfer is low, during the daytime, and heat bubbles caused by heterogeneous surface heating destabilize the atmosphere and promote turbulence formation, resulting in a larger momentum flux. e heat fluxes show a single-peak diurnal variation pattern in Figure 6, with the sensible and latent heat flux peaking at 148 W·m − 2 and 76 W·m − 2 at approximately 14:00 during the day. e latent heat fluxes are almost always positive, and the sensible heat fluxes are positive during the day and negative at night because the surface is exposed to solar radiation during the day; the ground temperature rises and transfers heat to the atmosphere, while at night, the ground temperature is lower than the air temperature due to ground cooling, resulting in downward heat transfer. e peak heat fluxes in the different seasons show that the heat fluxes are higher in spring and summer than in autumn and winter, with diurnal maximum sensible heat fluxes of approximately 200 W·m − 2 at midday in spring and summer and only approximately 100 W·m − 2 in autumn and winter. e latent heat flux in winter is much lower than the sensible heat flux, which is related to the lower atmospheric water vapor in winter.

Diurnal Variation in the Turbulent Kinetic Energy Budget
Term. Moeng and Wyngaard [41] pointed out that the energy needed to sustain turbulence development comes mainly from the shear and buoyancy terms of the turbulence kinetic energy equation (TKE). To simplify the calculations, the coordinate system direction u is taken to be in the same direction as the mean wind direction, the turbulence is assumed to be uniform, and subsidence is ignored. e following equation for the kinetic energy of turbulence is obtained: e left-hand side of equation (12) represents the local variation in the turbulent kinetic energy (ze/zt), the first term on the right-hand side is the buoyancy term (g/Tw ′ T ′ ), and the second term is the shear generation term (the shear term is simplified as − u ′ w ′ (zu/zz) when the u direction of the coordinate system coincides with the mean wind direction). zu/zz is obtained from the wind velocity data of the boundary layer gradient towers at 2 m and 4 m, and w ′ T ′ and u ′ w ′ are calculated using the EC system to reach the buoyancy term and the shear term of the turbulence kinetic energy budget equation. Figure 7 shows the diurnal variation in the turbulent kinetic energy budget. e buoyancy term is mainly influenced by thermal stratification, with a maximum positive value at noon and a smaller negative value at night, which is basically the same as the diurnal variation in the heat flux. e shear term is mainly produced by wind velocity and surface friction, with a minimum value at noon and a maximum value at night, which is close to the diurnal variation in wind velocity but lags behind the diurnal variation in wind velocity by approximately 4 hours. Table 4 gives the shear and buoyancy terms for turbulence in different seasons. e buoyancy term is larger in spring and summer and smaller in autumn and winter, and the shear term is also significantly smaller in winter than in other seasons. e shear term is also significantly smaller in winter than in other seasons, which is consistent with the seasonal variation in the stability condition. e annual mean of the shear term is greater than that of the buoyancy term, which is 47.94 × 10 − 4 m 2 ·s − 3 and 11.31 × m 2 ·s − 3 , respectively. Under the influence of complex surface friction, Advances in Meteorology the contribution of the shear term to turbulence formation in the semiarid region of the Loess Plateau is greater than that of the buoyancy term.

Bulk Transfer Coefficient of Turbulence.
e bulk transfer coefficient is an important parameter for describing the transfer of turbulence at the surface; however, the few studies of bulk transfer coefficients conducted on the land surface mainly focus on uniform underlying surfaces, and the understanding of bulk transfer coefficients on complex terrain is still inadequate, which makes it very difficult to parameterize the boundary layer in atmospheric models. e bulk transfer coefficients of turbulence in complex underlying surfaces in semiarid regions are of great importance for improving the prediction of numerical models and estimating momentum and heat fluxes. In circulation models, momentum transfer and heat transfer coefficients (C D and C H ) are often simply considered constant coefficients (C D is approximately 2 × 10 − 3 and C H is approximately 1.5 × 10 − 3 ), without considering the fact that differences in different areas and different underlying surfaces significantly affect the accuracy of the simulations.
Equations (13) and (14) are simplified definitions of the momentum and heat bulk transfer coefficient of turbulence. T s is the surface temperature, u and T are the average wind velocity and temperature at reference height (here, we select 4 m above ground as the reference height), respectively, and the covariances w ′ u ′ and w ′ T ′ are directly evaluated from the sonic anemometer measurements. Figure 8 shows the diurnal variation in the turbulent bulk transfer coefficients of momentum and heat and their distribution with stability. e turbulent bulk transfer of   Advances in Meteorology momentum and heat is highest at noon and smaller at night; as the stability increases, both C D and C H show a decreasing trend, and the bulk transfer coefficient of heat changes more markedly. Moreover, in all cases, the bulk transfer coefficient of momentum is larger than that of heat. Table 5        8.54 × 10 − 3 and 2.52 × 10 − 3 , which are not listed in Table 5, respectively. e values of C D and C H are much larger than those typically used in models, which should be attributed to the complex terrain in the Loess Plateau. Owing to the same small values used in numerical models in complex terrain as in flat areas, the heat and momentum transport in the complex terrain would be underestimated.

Conclusions
e EC and boundary layer gradient observatory tower data from the Semi-Arid Climate and Environment Observatory of Lanzhou University (SACOL) were processed with the EC method for quality control and data processing to obtain the statistical parameters of the surface layer turbulence in order to analyze the turbulence characteristics of the complex terrain in the semiarid zone of the Loess Plateau. e formation of turbulence in terms of dynamics and thermal conditions is discussed in this study. e main results are summarized as follows.
e frictional action of the complex surface increases turbulent mixing in the surface layer and weakens stratification, resulting in predominantly neutral/near-neutral and weakly unstable/weakly stable conditions in the semiarid zone of the Loess Plateau. e daytime stratification is predominantly unstable, with a mean value of − 0.27, while the nighttime stratification is predominantly stable, with a mean value of 0.19. Stability and turbulence intensity show opposite diurnal variations, reaching a minimum and maximum, respectively, at noon, with turbulence intensity decreasing significantly with increasing stability. e turbulence intensity is shown by I u ≈ I v > I w , which is different from the general pattern of I u > I v > I w . e average values of I u and I v are 0.42 and 0.40, respectively, and the average value of I w is 0.14, indicating that lateral wind makes a large contribution to turbulence. e dimensionless standard deviations of u, v, and w components of the wind velocity under nonneutral stratification satisfy the 1/3 law relatively well. e dimensionless standard deviations of the wind velocity are significantly larger in the horizontal direction than in the vertical direction, which is related to the larger turbulence scale in the horizontal direction. e dimensionless standard deviations of the three-dimensional wind velocity σ u /u * , σ v /u * , and σ w /u * at SACOL under near-neutral stratification are 3.35, 2.98, and 1.26, respectively. e influence of complex terrain is mainly observed in the horizontal direction, resulting in a significantly higher standard deviation in the horizontal wind velocity over complex terrain than over a flat and uniform surface. e turbulent kinetic energy budget term and the bulk transfer coefficient of turbulence have been calculated, and their diurnal variation patterns are given. e shear terms greater than the buoyancy term in the semiarid zone of the Loess Plateau are 47.94 × 10 − 4 m 2 ·s − 3 and 11.31 × 10 − 4 m 2 ·s − 3 , respectively, indicating that wind shear contributes more than buoyancy to turbulence formation in the semiarid zone of the Loess Plateau. e annual mean of the momentum transfer coefficient C D is larger than the heat transfer coefficient C H under near-neutral conditions, which have values of 8.54 × 10 − 3 and 2.52 × 10 − 3 , respectively.
Data Availability e topographic map in Figure 1(a) is drawn using ASTER GDEM v2 data. ASTER GDEM is a product of METI NASA and is available at http://srtm.csi.cgiar.org. e data used in this study are archived in the sharing data of SACOL (http:// climate.lzu.edu.cn/English/Home_Page.htm) with an application.

Conflicts of Interest
e authors declare that there are no conflicts of interest regarding the publication of this paper.