Diurnal and Seasonal Variations of Surface Energy and CO 2 Fluxes over a Site in Western Tibetan Plateau

: Land surface process observations in the western Tibet Plateau (TP) are limited because of the abominable natural conditions. During the ﬁeld campaign of the Third Tibetan Plateau Atmospheric Scientiﬁc Experiment (TIPEX III), continuous measurements on the four radiation ﬂuxes (downward / upward short / long-wave radiations), three heat ﬂuxes (turbulent sensible / latent heat ﬂuxes and soil heat ﬂux) and also CO 2 ﬂux were collected from June 2015 through January 2017 at Shiquanhe (32.50 ◦ N, 80.08 ◦ E, 4279.3 m above sea level) in the western Tibetan Plateau. Diurnal and seasonal variation characteristics of these surface energy and CO 2 ﬂuxes were presented and analyzed in this study. Results show that (1) diurnal variations of the seven energy ﬂuxes were found with di ﬀ erent magnitudes, (2) seasonal variations appeared for the seven energy ﬂuxes with their maxima in summer and minima in winter, (3) diurnal and seasonal variations of respiration caused by the biological and chemical processes within the soil were found, and absorption (release) of CO 2 around 0.1 mg m − 2 s − 1 occurred at afternoon of summer (midnight of winter), but the absorption and release generally canceled out from a yearly perspective; and (4) the surface energy balance ratio went through both diurnal and seasonal cycles, and in summer months the slopes of the ﬁtting curve were above 0.6, but in winter months they were around 0.5. Comparing the results of the Shiquanhe site with the central and eastern TP sites, it was found that (1) they all generally had similar seasonal and diurnal variations of the ﬂuxes, (2) caused by the low rainfall quantity, latent heat ﬂux at Shiquanhe (daily daytime mean always less than 90 W m − 2 ) was distinctively smaller than at the central and eastern TP sites during the wet season (generally larger than 100 W m − 2 ), and (3) a ﬀ ected by various factors, the residual energy was comparatively larger at Shiquanhe, which led to a small surface energy balance ratio.

where is the measured flux value at r, is the source emission rate at r + r', R is the domain of integration, f is the transfer function (footprint function), and it can be expressed in terms of a crosswind-integrated footprint and a crosswind dispersion function. A detailed description of the model can be found in Kljun et al. [35]. The footprint contour lines of 50%, 70%, and 90% are shown in Figure 1c. It can be found that 90% of the measured fluxes come from within 300 m of the observation tower, and the surface is relatively flat within 70% of the source area. However, there are some buildings and roads located in 90% of the source area.

Fast Response Measurements
The eddy covariance instruments were installed at 5.45 m (sensors' center height) above the ground surface (Figure 1d), and they consisted of a three-dimensional ultrasonic anemometer, a gas analyzer, and a controller. The three-dimensional ultrasonic anemometer (CSAT3A, Campbell To estimate both the source area and the relative contribution of each surface elements to the measured fluxes, the footprint analysis is conducted by using the Kljun et al. model [35]: (1) where ϕ is the measured flux value at r, Q ϕ is the source emission rate at r + r', R is the domain of integration, f is the transfer function (footprint function), and it can be expressed in terms of a crosswind-integrated footprint and a crosswind dispersion function. A detailed description of the model can be found in Kljun et al. [35]. The footprint contour lines of 50%, 70%, and 90% are shown in Figure 1c. It can be found that 90% of the measured fluxes come from within 300 m of the observation tower, and the surface is relatively flat within 70% of the source area. However, there are some buildings and roads located in 90% of the source area.

Fast Response Measurements
The eddy covariance instruments were installed at 5.45 m (sensors' center height) above the ground surface (Figure 1d), and they consisted of a three-dimensional ultrasonic anemometer, a gas analyzer, and a controller. The three-dimensional ultrasonic anemometer (CSAT3A, Campbell Scientific Inc., USA) was used to measure wind components (i.e., u, v, and w) and sonic virtual air temperature. The separated, in situ, open-path, mid-infrared absorption gas analyzer (EC150, Campbell Scientific Inc., Logan, UT, USA) was used to simultaneously measure the absolute density of carbon dioxide and water vapor. A controller (EC100, Campbell Scientific Inc., Logan, UT, USA) was used to control the gas analyzer and sonic anemometer, and synchronously sample data from them. The controller was also used to control a temperature thermistor probe and a barometer to measure the ambient air temperature and barometric pressure, respectively. All signals of the sensors were recorded at a sampling rate of 10 Hz and averaged over 30-min periods by a data-logger (CR3000, Campbell Scientific Inc., Logan, UT, USA). These sensors were calibrated before installation. Quality control and turbulent flux calculation were carried out with EddyPro 6. 2.1 (LI-COR Inc., Lincoln, USA, 2017) for the eddy covariance data. The processing steps and flux corrections are listed as follows: (1) Remove the spikes of the raw time series with a criterion of X − 4σ < x(t) < X + 4σ , where x(t) denotes the measurement (i.e., horizontal wind speed u and v, vertical wind speed w, etc.), X is the mean over the data block and σ is the standard deviation [36], (2) skip the turbulent data block missing more than 10%, (3) use double rotation method for the sonic anemometer tilt correction [37], (4) use 30-min as the flux averaging interval to compute turbulent fluxes, (5) use the WPL (Webb, Pearman, and Leuning) method for density corrections [38], (6) use high/low-frequency spectral corrections to compensate for flux losses [39], (7) remove the flux data in rainy (snowy) or fog conditions as eddy-covariance systems do not work properly under such environment.

Slow Response Measurements
Other supporting data, such as radiation, air temperature, air pressure, relative humidity, soil moisture, soil temperature, and soil heat flux, were collected synchronously during the experiment. Upward/downward shortwave/longwave radiation components were measured with two pyranometers (CMP22, Kipp & Zonen B.V., Delft, The Netherlands) and two pyrgeometers (CGR4, Kipp & Zonen B.V., Delft, The Netherlands) mounted at the height of 1.5 m. Soil heat flux at a depth of 0.05 m was measured by a heat flux plate (HFP01SC, Hukseflux Thermal Sensors B.V., Delft, The Netherlands), while soil temperature and volumetric water content at five depths (0.05, 0.1, 0.2, 0.4, and 0.80 m) were measured by a thermocouple (109, Campbell Scientific Inc., Logan, UT, USA) and a water content reflectometer (CS 616, Campbell Scientific Inc., Logan, UT, USA), respectively. The ground surface temperature was measured by an infrared radiometer (IRR-P, Apogee Inc., Austin, USA) mounted at the height of 1.5 m AGL. A tipping bucket rain gauge (TE525M, Texas Electronics Inc., Dallas, TX, USA) mounted at the height of 70 cm was used to measure precipitation. All the data were sampled each minute, averaged (Precipitation was accumulated.) every 30 min, and recorded with a CF card by a data-logger (CR3000, Campbell Scientific Inc., Logan, UT, USA).

Methods
The surface energy balance over the Gobi canopy can be approximated by: Atmosphere 2020, 11, 260 5 of 18 where Rn is the net radiation (W m −2 ), H and LE are the sensible and latent heat fluxes (W m −2 ), G 0 is the soil heat flux at the ground surface (W m −2 ), Re (W m −2 ) is the residual energy (W m −2 ). Rn can be approximated with the four components of radiation fluxes: where DSR is the downward shortwave radiation (W m −2 ), DLR is the downward longwave radiation (W m −2 ), USR is the upward shortwave radiation (W m −2 ), and ULR is the upward longwave radiation (W m −2 ). Turbulent sensible and latent heat fluxes are calculated as: where ρ, C p , and L are the density of air (kg m −3 ), the specific heat of air (J kg −1 K −1 ), and the latent heat of vaporization (J kg −1 ), respectively. w', T', and q' are the fluctuations in the vertical wind component (m s −1 ), air temperature (K) and specific humidity (g Kg −1 ), respectively. Relevant calculation methods can be found in [40,41]. The soil heat flux at the surface, G 0 , is estimated with: where G 0.05 is the soil heat flux measured at depth of 0.05 m with a heat flow plate, C g is the volumetric heat capacity of the soil (J m −3 K −1 ), which can be derived from soil components [42], ∆z is the thickness of the soil layer above the plate (i.e., 0.05 m) and T is the mean soil temperature of the layer, ∂T/∂t is the change rate of mean soil temperature. Here, T = (T 0 + T 0.05 )/2, where T 0 is the ground surface temperature, and T 0.05 is the soil temperature measured at the depth of 0.05 m. Neglecting the heat contribution of air in the soil, the volumetric heat capacity is estimated by [43]: where η is the volumetric moisture content of soil and η s is the saturation value of sand soil moisture (0.395 m 3 m −3 ) ( [43], its Table A9), C s is the volumetric heat capacity of dry sand soil (1.280 × 106 J m −3 K −1 ) ( [43], its Table A7), and C w is the volumetric heat capacity for water (4.186 × 106 J m −3 K −1 ). The residual energy Re, also called absolute energy balance residual (AEBR) [44], is defined as the following formula: It involves various processes, such as soil respiration, heat storage in soil and canopy [45,46], and phase change of water in the soil [47][48][49].
To assess the surface energy balance closure, two criteria are often used: Re (AEBR) and the surface energy balance ratio (SEBR) [50]. SEBR, marked as ε, also known as the surface heating rate, is defined as the ratio of the turbulent heat fluxes (H + LE) to the available energy fluxes (Rn − G 0 ) [51]. Ideally, when the surface energy balance (i.e., partitioning of Rn into LE, H, and G 0 ) is perfectly closed, ε should be equal to one. To obtain SEBR, one method is to derive linear regression coefficients from the least square relationship between (H + LE) and (Rn − G 0 ), and the slope is taken as the SEBR. The other method is to calculate the ratio between the cumulative sums of (H + LE) and (Rn -G 0 ) for a specific period [51][52][53], as follows: Atmosphere 2020, 11, 260 6 of 18 CO 2 flux is calculated using the following equation [54,55]: where F CO 2 is CO 2 flux (mg m −2 s −1 ) and c' is the fluctuation in the CO 2 concentration (mg m −3 ) [56].

The Meteorological Background
The variations of the daily mean wind speed (WS in m s −1 ), air temperature (T air in K), specific humidity (q in g kg −1 ), air pressure (P in hPa), and daily accumulated precipitation (Prec. in mm) from June 2015 to January 2017 are shown in Figure 2. Generally, all the meteorological elements had obvious seasonal cycles. June to September is categorized as wet season in this study due to the relatively high specific humidity and frequent precipitation events, while October to May is defined as dry season.
Atmosphere 2020, 11, 260 6 of 17 relatively high specific humidity and frequent precipitation events, while October to May is defined as dry season. At this site, the wind speed was relatively low with its daily mean mostly below 3 m s −1 for the whole year, and the maximum daily mean wind speed in this observation period was 5.8 m s −1 , but the maximum hourly mean wind speed reached 8.9 m s −1 . South-westerlies were dominant at this site and showed little season variation. However, the seasonal variation of air temperature was significant. The mean air temperature of the entire experiment was 276.6 K. In 2016, the annual mean air temperature was 276.5 K, the highest (lowest) monthly mean air temperature was 288.3 K (263.5 K) appearing in August (January). During the experiment period, daily mean air temperature reached a maximum of 291.3 K in July 2016, and a minimum of 257.5 K in January 2017, while the maximum 30-min air temperature of 298.1 K occurred at 19:30, 12 July 2016, and the minimum 30-min air temperature (251.0 K) occurred at 9:00, 25 December 2015.
The specific humidity was relatively low, and the mean specific humidity of the entire experiment was 2.6 g kg −1 . The difference in specific humidity between wet and dry season was At this site, the wind speed was relatively low with its daily mean mostly below 3 m s −1 for the whole year, and the maximum daily mean wind speed in this observation period was 5.8 m s −1 , but the maximum hourly mean wind speed reached 8.9 m s −1 . South-westerlies were dominant at this site and showed little season variation. However, the seasonal variation of air temperature was significant. The mean air temperature of the entire experiment was 276.6 K. In 2016, the annual mean air temperature was 276.5 K, the highest (lowest) monthly mean air temperature was 288.3 K (263.5 K) appearing in August (January). During the experiment period, daily mean air temperature reached a maximum of 291.3 K in July 2016, and a minimum of 257.5 K in January 2017, while the maximum 30-min air The specific humidity was relatively low, and the mean specific humidity of the entire experiment was 2.6 g kg −1 . The difference in specific humidity between wet and dry season was obvious. In the year 2016, the average of specific humidity in the wet season was 5.3 g kg −1 , and in the dry season, it was 1.2 g kg −1 . In response to the precipitation, the 30-min specific humidity reached 12 g kg −1 in the wet season (not shown) but kept less than 2 g kg −1 in the dry season.
Almost all precipitation occurred in June-September in both 2015 and 2016. The total precipitation amounts from June to September were 119.9 and 58.9 mm in 2015 and 2016, respectively. The annual maximum daily precipitation of 20.4 mm occurred on 4 August 2015, and 12.8 mm occurred on 5 August 2016.
The 30-min air pressure ranged from 612.8 to 594.1 hPa with an average of 605.0 hPa in the observation period. Different from the plain area, the atmospheric pressure at Shiquanhe in the warm summer (June, July, and August) was slightly higher than in the cold winter (December, January, and February), exhibiting the typical alpine air pressure characteristics. In winter, some strong low-pressure center moved over this place and decreased the daily mean air pressure to 595.0 hPa.

Radiation Components
To investigate the diurnal variations of the radiation components, the monthly averaged diurnal variations of the four radiation components (DSR, USR, DLR, and ULR) and the net radiation (Rn) are given in  15, 142.76, 196.25, and 363.33 W m −2 ) in the Jun, Jun, August, and Jun (December, January, January, and January), respectively. The monthly variation phase of DLR was one to two months behind the other radiation components, because of the slow response of the temperature of the upper atmosphere to the solar altitude change.
Rn was calculated with Equation (2) by using the four radiation components (i.e., DSR, USR, DLR, and ULR). Generally, Rn was mainly affected by DSR since it was the largest of the four radiation components and also the driving force of the other three. Here, Rn varied roughly from -130 at night to 600 W m −2 at noon, heated the surface in the daytime while imposed a cooling effect on the surface after sunset. Similar results can also be found in the observations from Lhasa and several other stations in the central, Nyingchi and Litang in the Eastern TP, and Haibei in the northwestern TP [57][58][59]. In the year 2016, the yearly averaged Rn was 64.70 W m −2 , and the largest (smallest) monthly averaged daily maximum of Rn was 589.81 W m −2 (311.64 W m −2 ) at noon in June (December), while the largest (smallest) monthly averaged daily minimum of Rn was −95.99 W m −2 (−134.70 W m −2 ) at midnight in January (September).
The surface albedo was calculated as the ratio of USR to DSR. Affected by the diurnal variation of solar altitude, the surface albedo was also found changing diurnally between roughly 0.2 and 0.4, with the maximum in the morning and minimum at noon. Besides, caused by the seasonal variation of the soil water content, the surface albedo was slightly smaller in the wet season than in the dry season. In the year 2016, the averaged surface albedo was 0.21 in the wet season and was 0.25 in the dry season.
fluctuations of DLR (about 50 W m −2 ). In the year 2016, the monthly averaged daily maximum DSR, USR, DLR, and ULR showed maximum (minimum) values of 1078.24, 206.06, 321.57, and 564.0 W m −2 (632. 15, 142.76, 196.25, and 363.33 W m −2 ) in the Jun, Jun, August, and Jun (December, January, January, and January), respectively. The monthly variation phase of DLR was one to two months behind the other radiation components, because of the slow response of the temperature of the upper atmosphere to the solar altitude change. Rn was calculated with Equation (2) by using the four radiation components (i.e., DSR, USR, DLR, and ULR). Generally, Rn was mainly affected by DSR since it was the largest of the four  The diurnal variation of latent heat flux was small during the dry season and at night values within 5 W m −2 were generally found, which could be attributed to the following factors: (1) Soil water was frozen when temperature was below zero, (2) soil water content was low and little water for evaporation at the surface, (3) little vegetation and then little water transpiration. The inter-month variation of latent heat was obvious. During the wet season, (1) the soil was thawed, (2) precipitation happened (as shown in Figure 2   CO 2 is one of the most important greenhouse gases and attracts attention globally. Wohlfahrt [60], Xie [61], and Mamtimin Ali [62] found that desert could behave as carbon sink, because: (1) The abiotic process and air-soil CO 2 concentration gradients in surface boundary layer caused such desert carbon uptake [62], (2) alkaline soil on land absorbed CO 2 with an inorganic, non-biological process [61], (3) a small amount of precipitation could activate bio-autotrophic respiration in the soil and heterogeneous respiration of soil organisms, resulting in enhanced autotrophic respiration and abnormal breathing of soil organisms in the surface moist layer. Here, measurements of F CO 2 were also carried out in the observation period. Figure 4d showed that F CO 2 was generally positive in autumn (September, October, and November) and winter but showed relatively significant negative values in spring (March, April, and May) and summer. Large absorption (release) of CO 2 occurred in the afternoon of summer (midnight of winter) and the monthly mean daily maximum absorption (release) was around 0.1 mg m −2 s −1 in Jun (December) 2016. From a yearly perspective, these positive and negative values generally canceled out, and the yearly average of F CO 2 in 2016 was close to zero. In this period, the diurnal variation in late autumn or early winter and late spring or early summer was generally more significant compared to the other periods, while the largest (smallest) monthly averaged daily maximum of F CO 2 was 0.  Figure 5 shows the daily daytime (10:00-19:00 local time) means of the four radiation components (DSR, USR, DLR, and ULR), net radiation, and albedo. All the four radiation components exhibited significant seasonal variation due to the seasonal variation of solar altitude, and DSR, USR, and DLR showed some scatters caused by clouds. DSR, USR, and ULR generally shared the same phase of seasonal variation and had their maxima (minima) in June (December) with values of 885.52 W m −2 , 181.15 W m −2 and 530.65 W m −2 (178.39 W m −2 , 40.42 W m −2 , and 286.95 W m −2 ), respectively. However, the seasonal variation phase of DLR was behind those of DSR and USR for about one to two months, respectively. DLR showed maximum (minimum) in August (January) with values of 336.53 W m −2 (151.34 W m −2 ). Caused by the lag of DLR seasonal variation, seasonal variation of daily mean Rn was also slightly behind the phase of DSR and USR, with its maximum (minimum) value of 465.41 W m −2 (68.79 W m −2 ) shown in July (December). The surface was bare and had little vegetation all year round so that the seasonal variation of albedo was relatively invariant, and for most of the observation period, it fluctuated around 0.20 ± 0.03. Some scattered dots lower than 0.17 in the wet season were caused by the precipitation event. Skipping these scattered dots, the slightly higher albedo was shown around December than around July, which was caused by the drier (moister) soil in the winter (summer) season.

Radiation Components
Atmosphere 2020, 11,260 10 of 17 in the wet season were caused by the precipitation event. Skipping these scattered dots, the slightly higher albedo was shown around December than around July, which was caused by the drier (moister) soil in the winter (summer) season.   Figure 6 shows the seasonal variations of daily daytime means of H, LE, G 0 , and F CO 2 , all of which had remarkable seasonal variations. The maximum of the daily daytime mean H was shown in May and a little ahead of that of Rn, but the minimum was in December and consistent with that of Rn.

Heat and CO 2 Fluxes
In the year 2016, the maximum and minimum of the daily daytime mean H were 189.77 W m −2 and 32.70 W m −2 , respectively. The daily daytime averaged LE was mostly within 10 W m −2 in the dry season but could be larger than 30 W m −2 in the wet season. The maximum daily daytime mean LE of 88.48 W m −2 occurred in August 2015 associated with a precipitation event, which was much smaller than the observations at Lhasa in the central, Nyingchi, Wudaoliang, and Litang in the eastern, Haibei, Nagqu, and Tanggula in the northeastern TP, where the latent heat flux was generally larger than 100 W m −2 in wet season [24,33,[57][58][59]63,64]. The fluctuation of G 0 was well correlated with Rn. The daily daytime averaged G 0 was always positive showing that the land constantly absorbed heat in the daytime. In the year 2016, the maximum (minimum) daily daytime mean G 0 was in June (December) with a value of 185.28 W m −2 (26.65 W m −2 ).

Energy Partitioning and Energy Balance Closure
To further understand the energy partitioning, the monthly averaged diurnal variation of H/Rn, LE/Rn, G0/Rn, Bowen ratio (β = H/LE), Re and SEBR are shown in Figure 7. As H was relatively small in the nighttime, H/Rn was close to zero at night in all months. After sunrise, H/Rn started to increase. Further, the diurnal variation phase of H was about one hour behind that of Rn (not shown), H/Rn reached to near unity at the time of sunset. LE accounted for a very small amount of Rn at this site, and LE/Rn was mostly around zero in the dry season but fluctuated within ±0.25 in the wet season. Diurnal variation of LE/Rn from negative at night to positive in the daytime during the wet season was caused by the change of the sign of Rn. G0/Rn was generally positive throughout the whole period. During the daytime, G0/Rn had a quasi-sinusoidal fluctuation with a maximum of about 0.4-0.5 around noon, while during nighttime, G0/Rn was generally invariant with a value of about 0.8-0.9 indicating that surface energy loss through radiation was mostly compensated by G0. The absolute values of the Bowen ratio in this site could be extremely large (>50) in the dry season which was caused by the near-zero latent heat flux. While in Figure 7d, only the absolute values smaller than 50 are shown. In the wet season, diurnal quasi-sinusoidal fluctuations of β were found, with the maximum values shown around noon. Besides, the year-to-year difference was also presented and β Daily daytime mean F CO 2 was positive in late autumn and early winter but negative in late spring and early summer. The negative period was relatively longer and last about eight months (February 2016 to September 2016), while the positive period was short and last for about three months (October to December 2016). In the year 2016, the maximum (minimum) daily daytime mean F CO 2 was in December (May) with a value of 0.13 mg m −2 s −1 (−0.13 mg m −2 s −1 ).

Energy Partitioning and Energy Balance Closure
To further understand the energy partitioning, the monthly averaged diurnal variation of H/Rn, LE/Rn, G 0 /Rn, Bowen ratio (β = H/LE), Re and SEBR are shown in Figure 7. As H was relatively small in the nighttime, H/Rn was close to zero at night in all months. After sunrise, H/Rn started to increase. Further, the diurnal variation phase of H was about one hour behind that of Rn (not shown), H/Rn reached to near unity at the time of sunset. LE accounted for a very small amount of Rn at this site, and LE/Rn was mostly around zero in the dry season but fluctuated within ±0.25 in the wet season. Diurnal variation of LE/Rn from negative at night to positive in the daytime during the wet season was caused by the change of the sign of Rn. G 0 /Rn was generally positive throughout the whole period. During the daytime, G 0 /Rn had a quasi-sinusoidal fluctuation with a maximum of about 0.4-0.5 around noon, while during nighttime, G 0 /Rn was generally invariant with a value of about 0.8-0.9 indicating that surface energy loss through radiation was mostly compensated by G 0 . The absolute values of the Bowen ratio in this site could be extremely large (>50) in the dry season which was caused by the near-zero latent heat flux. While in Figure 7d, only the absolute values smaller than 50 are shown. In the wet season, diurnal quasi-sinusoidal fluctuations of β were found, with the maximum values shown around noon. Besides, the year-to-year difference was also presented and β in the wet season of 2016 was generally larger than that of 2015, which was related to higher rainfall quantity in 2015.
In the dry season, at night, H was close to zero so that β also fluctuated around zero, but during the daytime, its absolute values could reach more than 50 and its sign was determined by the direction of the near-zero latent heat flux. Respondent to the quasi-sinusoidal variation of DSR, the residual energy Re in the daytime had also a quasi-sinusoidal fluctuation with a maximum of about 100-150 W m −2 around noon and a minimum reaching −50 W m −2 at sunset, while at night it was generally invariant around −20 W m −2 . SEBR showed significant diurnal variation as well. It was small and could be even negative in the nighttime, but after sunrise, it started to increase and reaching unity at sunset. The mean SEBR in the total period was 0.57, which was similar to but a little smaller than the finding of Xin et al. (0.65 from July to September) [34]. Many influencing factors contributed to the large Re and small SEBR. For example, the uncertainties caused by H and LE instrumental errors. Xin et al. pointed out that the uncertainties of Re caused by H and LE instrumental errors were generally smaller during the nighttime than those during the daytime, and the relative instrumental error in measuring H was nearly 12.8% whereas that for LE was 12.5% with the highest value reaching 22.9% [34]. Xin et al. also indicated that the soil heat storage term represented a significant source of uncertainty (~10%) to close the surface energy budget, which was affected by the soil respiration, heat storage in canopy, phase change of water in the soil, energy loss due to advection, while all of them were neglected in this study due to the difficulties in measurement. Further, even if the measurement's uncertainties could be reduced to zero, and all energy exchange processes and losses were measured precisely, the residual energy still could not be close to zero because of the lagging of the diurnal wave phases of G 0 to those of H and LE [65][66][67]. In our study, in addition to the above reasons, the following factors could also contribute to the non-closure: (1) The effect of freezing and thawing on the calculation of G 0 , as the latent heat related to the water phase change in the soil was not considered. (2) The lower moisture content in the atmosphere. The moisture was very low in our study area, which might lead to significant uncertainty in the measurement of latent heat. (3) The complex composition of soil surface. We regarded the soil here as uniform sand to set the relevant soil parameters (e.g., saturated soil moisture and soil volume heat capacity). However, the underlying surface in the observation site was a mixture of sand and a small number of pebbles. (4) The influence of nearby buildings and human activities on the site.
As shown in Figure 8, SEBR went through a regular diurnal cycle in every month, starting from roughly the origin point at morning then to the left of 1:1 line, it crossed the 1:1 line around noon and to its right at afternoon, at the same time both turbulent heat fluxes (H + LE) and surface available energy (Rn − G 0 ) reached their maxima, then at night, it retreated to the origin point gradually. Seasonal variation was also presented, and in summer months, the slopes of the fitting curve were mostly above 0.6, but in winter months, they were around 0.5.
of nearby buildings and human activities on the site.
As shown in Figure 8, SEBR went through a regular diurnal cycle in every month, starting from roughly the origin point at morning then to the left of 1:1 line, it crossed the 1:1 line around noon and to its right at afternoon, at the same time both turbulent heat fluxes (H + LE) and surface available energy (Rn − G0) reached their maxima, then at night, it retreated to the origin point gradually.
Seasonal variation was also presented, and in summer months, the slopes of the fitting curve were mostly above 0.6, but in winter months, they were around 0.5.

Summary and Conclusions
To examine the land-atmosphere interaction and understand the surface energy partitioning and CO2 exchange over western TP, continuous measurements on radiation and turbulent heat fluxes, soil heat fluxes and CO2 flux were conducted from June 2015 through January 2017 at Shiquanhe (32.50° N, 80.08° E, 4279.3 m above sea level) in the western Tibetan Plateau.
The Shiquanhe site was dominated by weak south-westerlies (generally below 3 m s −1 ) and could

Summary and Conclusions
To examine the land-atmosphere interaction and understand the surface energy partitioning and CO 2 exchange over western TP, continuous measurements on radiation and turbulent heat fluxes, soil heat fluxes and CO 2 flux were conducted from June 2015 through January 2017 at Shiquanhe (32.50 • N, 80.08 • E, 4279.3 m above sea level) in the western Tibetan Plateau.
The Shiquanhe site was dominated by weak south-westerlies (generally below 3 m s −1 ) and could be divided into a wet (June to September) and a dry season (October to May). All the rainfall was in the wet season, and the 30-min specific humidity reached 12 g kg −1 in the wet season but kept less than 2 g kg −1 in the dry season. Diurnal variations of the four radiation components, three heat fluxes, and CO 2, flux were found. In 2016, the magnitudes of the monthly averaged diurnal variations for DSR, USR, and ULR were the largest in June For F CO 2 , diurnal variation of respiration caused by the biological and chemical processes within the soil was found, large absorption (release) of CO 2 around 0.1 mg m −2 s −1 occurred in the afternoon of summer (midnight of winter). However, from a yearly perspective, the absorption and release generally canceled out, and yearly average of F CO 2 in 2016 was close to zero.
Seasonal variations also appeared for the daily daytime (10:00-19:00 local time) mean radiation, turbulent, and soil heat fluxes, and the maxima were all in summer and minimum in winter. DSR, USR, and ULR generally shared the same phase of seasonal variation and had their maxima (minima) in June (December) with values of 885.52, 181.15, and 530.65 W m −2 (178. 39, 40.42, and 286.95 W m −2 ), respectively. While DLR was one to two months behind and had its maximum (minimum) in August (January) with values of 336.53 W m −2 (151.34 W m −2 ). Rn was also slightly behind the phase of DSR, with its maximum (minimum) value of 465.41 W m −2 (68.79 W m −2 ) shown in July (December). The daily daytime mean albedo had only slight seasonal variation and fluctuated within 0.20 ± 0.03, and its maximum was in winter associating with the drier soil. The maximum (minimum) daily daytime mean H was shown in May (December) with a value of 189.77 and 32.70 W m −2 in 2016, while the maximum LE larger than 30 W m −2 generally occurred in July or August associated with precipitation events. Daily daytime mean G 0 was always positive, and its yearly maximum (minimum) in 2016 occurred in June (December) with a value of 185.28 W m −2 (26.65 Wm −2 ). In the daytime, the soil absorbed (released) CO 2 around May (December) at a rate of 0.13 mg m −2 s −1 .
Diurnal and seasonal variation of H/Rn, LE/Rn, G 0 /Rn, Bowen ratio (β = H/LE), Re and SEBR are also presented in this study. H/Rn was relatively large and could reach unity at the time of sunset, while LE/Rn was mostly around zero in the dry season but fluctuated within ±0.25 in the wet season. G 0 /Rn was about 0.4-0.5 around noon and 0.8-0.9 at midnight. β had extremely large absolute values (>50) because of the near-zero of LE in the dry season, while in the wet season, it was within 10-50 at noon. Re in the daytime had a maximum of about 100-150 W m −2 around noon and a minimum reaching −50 W m −2 at sunset, and at night it was generally invariant around −20 W m −2 . The mean SEBR in the total period was 0.58. It went through both regular diurnal cycle and season variation, and in summer months, the slopes of the fitting curve were mostly above 0.6, but in winter months, they were around 0.5. The measurement errors of the fluxes, the effects of frozen soil on the calculation of G 0 , the lower moisture content in the atmosphere, the complex composition of the soil surface, and the influence of nearby buildings and human activities were the possible reasons to the small SEBR.
Comparing the results from the Shiquanhe site with the central and eastern TP sites [24,33,[57][58][59]63,64], it could be found that (1) they all generally had similar seasonal and diurnal variation of the fluxes, (2) caused by the low rainfall quantity, LE at Shiquanhe (daily daytime mean always less than 90 W m −2 ) was distinctively smaller than the central and eastern TP sites during the wet season, where the latent heat flux was generally larger than 100 W m −2 , and (3) affected by various factors, the residual energy was comparatively larger at Shiquanhe, which led to a small surface energy balance ratio.