Abstract

In this study, the characteristics of energy flux partitioning and evapotranspiration of P. euphratica forests were examined in the extreme arid region of Northwest China. Energy balance closure of the ecosystem was approximately 72% ( + LE = 0.72 ∗ , , ), where is the net radiation, is the soil heat flux, is the sensible heat flux, and LE is the latent heat flux. LE was the main term of energy consumption at annual time scale because of higher value in the growing season. The ratios of the latent (LE) and sensible () heat fluxes to net radiation were 0.47 and 0.28 throughout the year, respectively. Moreover, the yearly evapotranspiration of P. euphratica forests was 744 mm year−1. And the mean daily ET was 5.09 mm·d−1 in the vibrant growing season. In particular, a small spike in the actual evapotranspiration distribution occurred during the soil ablation period due to the higher temperature and sufficient soil moisture associated with soil thawing. This period is accompanied by a series of physical processes, such as moisture transfer and heat exchange.

1. Introduction

The Ejina oasis is located at the interface between the Alxa highland and desert in the lower reach of Heihe River Basin. This area exhibits a dry climate, harsh natural conditions, and a fragile ecology, and the region is highly vulnerable to climate change. Thus, the sound and steady progress of the Ejina oasis will aid in the prevention of further desertification and help maintain ecological security for the entire Heihe River Basin. Water is the key ecological factor in the inland river basin-oasis ecological processes. In semiarid ecosystems, evapotranspiration (ET) generally represents the greatest flux of water out of an ecosystem, and the relative partitioning of ET between plant and soil surfaces is controlled by water supply and available energy for ET, as well as the path resistance and capacitance along the water travel path [1]. Furthermore, the partitioning between ecosystem latent (LE) and sensible () heat fluxes is critical in determining the hydrological cycle, boundary layer development, weather, and climate [2]. A full understanding of the variability of partitioning across different climates and ecosystems and the mechanisms for this variability are necessary. Therefore, an in-depth study of energy flux partitioning and evapotranspiration in a desert riparian P. euphratica forest ecosystem is needed to increase the efficiency of water distribution and management in the Heihe River Basin and thereby aid the rational development of the oasis.

Natural vegetation in arid areas plays an essential role in the conservation of biodiversity and reduction of desertification [3]. The natural forest of the study area is dominated by P. euphratica. This is a light-demanding species that grows along rivers and represents a component of the natural desert forest ecosystems in China and Middle Eastern countries [4]. The growth of this species is dependent on the groundwater table, which is determined by the amount of water transfer in hyperarid regions. Of course, increased delivery of water downstream is beneficial. However, water shortage and overexploitation of water as well as conflicts between water demand and supply exist over the entire extent of the Heihe River Basin. The amount and time of ecological water conveyance need to be ensured according to the ecohydrological requirements. Thus, understanding the processes of energy flux partitioning and evapotranspiration of P. euphratica forests in different growth stages and their response to anthropogenic and natural disturbances (e.g., drought events, rain pulses, and irrigation) is of foremost importance to ecosystem restoration and estimation of catchment water balance.

The eddy covariance method, a micrometeorological technique, has been widely used to study exchanges of carbon dioxide, water vapor, and energy between the biosphere and atmosphere. Many of these studies have concentrated on the mass and energy exchange between forest ecosystems and atmosphere, such as tropical forests [5, 6], temperate forests [79], and boreal forests [1012]. In addition, several studies were conducted in endemic forest ecosystems, for instance, bamboo forests [13], olive trees [1416], riparian woodland [17, 18], and mangrove forest [19]. Nevertheless, few studies have used the eddy covariance method to measure the energy and water exchange between desert riparian P. euphratica forest and the atmosphere in hyperarid regions. For example, Yuan et al. (2015) reported that evapotranspiration controlled by vegetation phenological characteristics was mainly partitioned to plant transpiration and that groundwater was the main water source as well as the key factor controlling the water cycle process of the desert riparian forest in the lower Tarim River Basin [18]. However, the characteristics of the energy flux partitioning and evapotranspiration of the P. euphratica forest are still unclear for different phenophases and areas. These topics require intensive study. The objectives of this research are as follows: to study the characteristics and differences of diurnal variation in energy flux partitioning at different phenological phases in a desert riparian P. euphratica forest; to analyze the patterns of seasonal variation in evapotranspiration and crop coefficients in the P. euphratica forest; to determine the contribution of latent (LE) and sensible () heat fluxes to the available energy of the P. euphratica forest in 2014; to identify the evapotranspiration at different phenological phase and total yearly evapotranspiration.

2. Materials and Methods

2.1. Site Description

The study site is dominated by P. euphratica and is located in the lower basin of the Heihe River (101°7′25′′E, 41°59′34′′N) within the continental arid zone (Figure 1), where terrain was relatively regular with highly homogenous land surface cover. The values of the roughness length, zero-displacement height, and site’s slope were 1.38 m, 7.9 m, and 0.51 degrees, respectively. The study of the energy flux partitioning and evapotranspiration was conducted throughout the year of 2014. This period was then divided into six subperiods for further investigation. The climate at the site is a typical temperate continental climate, and it can also be classified as an extremely arid climate. The yearly potential evapotranspiration is approximately 2062 mm, while the rainfall is less than 35 mm. According to the data provided by the China Meteorological Data Sharing Service System (Climatic Data Center, National Meteorological Information Center, China Meteorological Administration, Beijing, China) for 1960–2012, the annual mean temperature is 9°C and the annual mean rainfall is 33.24 mm, with 60% of the rainfall occurring during June–October. The soil is derived from fluvial sediments mixed with gray-brown desert deposits [20]. The soil exhibits high salinity as well as low organic matter and nutrient content. The riparian forest is composed of broad-leaved P. euphratica and extensive thickets of T. ramosissima.

Because of overexploitation of the water resources in the middle reaches of Heihe River beginning in the 1960s, water resources discharged to lower reaches were decreased, which led to degradation of the surrounding ecology [21]. In order to restore the deteriorated ecohydrology in the lower reaches of the Heihe River, the ecological water conveyance project (EWCP) was established in 2002, and water resources have since been under strict artificial controls. As a result of this artificial water management, the groundwater level has increased, the vegetation has experienced a favorable change, and most plants can now acquire water access groundwater [22, 23]. The growth and distribution of vegetation were directly affected by groundwater. Besides, the level and the seasonal variation of evapotranspiration were important factors that impacted the amount and time of water conveyance. The time of water conveyance was September 15, 2014.

2.2. Energy Fluxes and Measurement of Environmental Variables

Fluxes of energy and water between the forest and the atmosphere were measured during January 1, 2014, to December 31, 2014, using the eddy covariance (EC) method. The EC system was installed on terrain that was relatively regular with highly homogenous land surface cover. The instruments included a three-dimensional sonic anemometer (CSAT3, Campbell Scientific, Inc., Logan, UT, USA) to measure wind speed and air temperature and an open-path CO2/H2O infrared gas analyzer (LI-7500, LI-COR Inc., Lincoln, NE, USA) to measure the water vapor density and carbon dioxide fluctuations, which was installed at 22 m above the ground. Solar radiation components were measured at 20 m above the ground by a 4-Component Net Radiation Sensor (CNR4, Kipp & Zonen, Delft, NL). The turbulence fluxes were recorded by a CR3000 data logger (Campbell Scientific, Logan, UT, USA) at a 10 Hz sampling frequency, and 30-minute block average fluxes were determined. Radiation data was collected with interval set to 10 minutes. In addition, air temperature and relative humidity were measured at 10 m and 28 m above the ground. Three soil heat flux plates were buried at 6 cm below the soil surface to measure soil heat flux. All these meteorological data were recorded by a CR1000 (Campbell Scientific, Logan, UT, USA).

2.3. Data Processes and Gap Filling

The half-hour sensible () and latent (LE) heat fluxes were calculated by Eddypro software following the standard quality control procedures including despiking, 2D coordinate rotation, time lag removal, frequency response correction using model spectra and transfer functions [24], and air density correction [25].

Data gaps (not exceeding two hours) were filled using linear interpolation. Large gaps in energy flux data were replaced with average values calculated using lookup table methods with assorted meteorological conditions [26, 27]. Then, mean diurnal variations (hourly averages binned over periods of 15 days) were used to fill residual gaps.

2.4. Calculations and Statistical Analysis

Surface hourly evapotranspiration of the desert riparian P. euphratica forest was calculated as follows:where LE is the latent heat flux (MJ·m−2·h−1); is the latent heat of vaporization of water (2.45 kJ·g−1); is water density (1 g·cm−3).

The FAO Penman-Monteith method to estimate daily reference evapotranspiration (mm·d−1) can be derived:where ET0 is reference evapotranspiration [mm·day−1]; is net radiation at the crop surface [MJ·m−2·day−1]; is soil heat flux density [MJ·m−2·day−1]; is mean daily air temperature at 2 m height [°C]; is wind speed at 2 m height [m·s−1]; is saturation vapor pressure [kPa]; is actual vapor pressure [kPa]; Δ is slope vapor pressure curve [kPa·°C−1]; and is the psychrometric constant [kPa·°C−1].

Thus, the crop coefficient of P. euphratica can be calculated as follows:where ET and ET0 are the actual and reference evapotranspiration, respectively.

The evaporative fraction (EF) is a derived variable to describe the partition of available energy into the sensible and latent heat fluxes. The EF was nearly constant during daytime and could be substituted by the midday (10:00–16:00) EF [28, 29]. The EF was calculated by the following equation [30]:where LE is latent heat fluxes and is the sensible heat fluxes.

3. Results and Discussion

3.1. Meteorological Conditions

In order to study the water and energy exchange between the P. euphratica forest and the atmosphere, it is necessary to understand the seasonality of local key environmental variables. Seasonal variations in air temperature (), photosynthetic active radiation (PAR), soil temperature (), potential evapotranspiration (EVP), vapor pressure deficit (VPD), soil volumetric water content (), and precipitation are presented in Figure 2.

The maximum and minimum value of PAR were 1109.98 μmol·m−2·s−1 (July 3) and 127.83 μmol·m−2·s−1 (October 3), respectively. The PAR showed significant day-to-day fluctuations due to frequent cloud cover, except in winter. In addition, the minimum and maximum daily average air temperature () at 28 m above the ground were −14.32 (January 8) and 31.29°C (July 16), respectively. Mean annual was 10.83°C with winter average of −7.11°C and summer average of 24.98°C. The mean annual VPD was 1.35 kPa, and VPD winter and summer averages were 0.27 and 2.52 kPa, respectively. The obvious seasonal variation of soil temperature () at 0, 2, and 4 cm depths was similar to that of . Annual precipitation was 16.4 mm in 2014, and most of the precipitation during 2014 was concentrated in the vibrant growing season. The seasonal variation of was small, and the value of was less than 4% throughout the year excluding precipitation days. In other words, the top-soil (upper 5 cm) was dry on nonprecipitation days, and the fluctuation of top-soil volumetric water content was dependent on the precipitation.

3.2. Seasonal and Diurnal Patterns in Energy Balance Components
3.2.1. Diurnal Patterns in Energy Balance Components at Each Stage

Because of the distinct differences found in the energy partitioning pattern and plant phenology, we divided the year into six periods: (a) the soil-frozen period (days of year (DOY) 323–64); (b) the soil ablation period (DOY 65–130); (c) the germination, flowering, and leaf expansion period (DOY 131–169); (d) the vibrant growing season (DOY 170–253); (e) the leaf senescence period including coloration and defoliation stages (DOY 254–288); and (f) the early soil-frozen period (DOY 289–322). The values of energy balance components and environmental variables were distinct for the different periods (Table 1).

Figure 3 shows the diurnal variation in the ensemble half-hourly means of , LE, , and during each period at the study site. The daily maximum of during each period was 380.72 W·m−2, 515.79 W·m−2, 576.28 W·m−2, 632.63 W·m−2, 547.05 W·m−2, and 407.53 W·m−2 for the first to the sixth stages, respectively. The magnitude and scope of variation of LE were distinct in each period. LE showed evident variance at each period with values ranging from −2.96 to 24.96 W m−2, 0.33 to 50.78 W m−2, 3.94 to 106.05 W m−2, 5.25 to 160.68 W m−2, 17.43 to 218.13 W m−2, and 22.02 to 360.52 W m−2 in the soil-frozen period, the early soil-frozen period, the soil ablation period, the leaf senescence period, the germination, flowering, and leaf expansion period, and the vibrant growing season, respectively. In addition, the magnitude and scope of were largest in the soil ablation period, and the value of ranged from −34.11 to 309.77 W·m−2. The value of ranged from −35.94 to 174.20 W·m−2, −34.11 to 309.77 W·m−2, −43.11 to 268.24 W·m−2, −41.89 to 157.64 W·m−2, −38.66 to 209.22 W·m−2, and −45.04 to 189.51 W·m−2 for the first to the sixth stage, respectively. Moreover, the maximum value of latent heat flux was approximately 360 W·m−2 in the vibrant growing season, while the maximum sensible heat flux of 310 W·m−2 occurred in the soil ablation period. The sensible heat flux at night was negative in each period and the result was consistent with other findings [31]. In contrast, the latent heat flux was positive at night during the growing season. Si et al. (2015) reported that significant stomatal opening persists during the night in P. euphratica forests, leading to the occurrence of substantial water losses, and nighttime sap flow was observed during the whole growing season [20]. This indicated that tree transpiration may have occurred during the night and can explain why the latent heat flux was positive in our study. The midday half-hourly value of was higher than that of LE in each period, except the vibrant growing season. During the daytime, the value of was higher than LE in the nongrowing season, was lower than LE in the vibrant growing season, and was lower than LE before or after the midday in the pregrowing season and senescence periods.

3.2.2. Seasonal Variation in Energy Balance Components

The seasonal variation of daily half-hourly means of , LE, , and is shown in Figure 4. exhibited significant seasonal and day-to-day variation, ranging from 12 W·m−2 (November 26) to 261 W·m−2 (July 3). was mostly converted into LE in the growing season (DOY 131–288). The seasonal variation of LE was distinctly related to plant phenology. And exhibited reverse course with the LE seasonal course in the growing season. As a result of water regulation on the Heihe River in mid-September, the groundwater level and soil water content increased (Figure 5). Soil freezing and thawing occurred in winter and spring, respectively. Freezing and thawing of soil moisture are accompanied by a series of physical processes, such as moisture transfer and heat exchange. With the start of soil thawing, and LE increased in the soil ablation period. LE continued to maintain a rapid increase and tended to decrease during the beginning of the vegetation growth in the germination, flowering, and leaf expansion period. In the vibrant growing season, LE sustained higher values ranging from approximately 150 to 200 W·m−2 and exhibited lower values ranging from approximately −20 to 60 W·m−2. Additionally, LE began to decrease and started to increase in the leaf senescence period. Over the whole soil-freezing period, ranged from approximately 0 to 40 W m−2 and LE ranged from approximately 0 to 10 W·m−2. For the soil heat flux (), there was a small yearly cycle in ranging from approximately −15 to 15 W·m−2. Studying evaporative fraction (EF) plays an important role in interpreting the components of energy budget and evapotranspiration (ET) [28]. Seasonal variation of was similar to the variation of LE. EF was distinct during different periods, and the mean value of EF was 0.14, 0.25, 0.49, 0.75, 0.44, and 0.24 from the first to sixth stage, respectively (Table 1).

3.3. Energy Balance and Evapotranspiration
3.3.1. Energy Balance Closure

Meteorological data (net radiation, ) and soil heat flux () are combined with the EC data (sensible and latent heat fluxes, and LE) to analyze the energy balance closure. and LE were obtained from EC data, was measured by a Net Radiation Sensor and was measured at 6 cm below the soil surface by three soil heat flux plates. The result of energy balance closure was determined by a linear regression method using 30 min averaged fluxes. The energy closure of the ecosystem was approximately 72% (, , and ) over the whole study period. Since the ratio of to net radiation () was less than 0.4% in the different phenological periods, the energy storage term () was neglected in the calculation process. Wilson et al. (2002) comprehensively evaluated energy balance closure across 22 sites and 50 site-years. This assessment indicated [32] a general lack of closure and the slope values ranged from 0.55 to 0.99 with a mean of . Thus, the surface energy fluxes ( + LE) are underestimated by approximately 27% relative to estimates of available energy (-), and the closure value was reasonable for the P. euphratica forest ecosystem. The energy balance closure in our study (72%) agreed with an eddy covariance based estimate of riparian forest (Populus fremontii S. Watson) in the western US, which was 75% [33], but was lower than the value estimated of riparian forest (Tamarix ramosissima) in other desert regions of Northwest China, which was 84% [18].

There is a general lack of the energy balance closure at eddy covariance sites, which could be attributed to several hypotheses [32], including () a systematic bias in instrumentation, () neglected energy sinks, () the loss of low or high frequency contributions to the turbulent flux, () neglected advection of scalars, and (5) sampling errors associated with different measurement source areas for terms in energy balance components. In addition, the imbalance of energy could also be caused by neglected energy storage term and biased soil heat flux in this paper. And it was difficult to isolate the contribution of these factors to the lack of energy balance. Heat flux plates have been used widely to measure soil heat flux for their simplicity and durability but suffer from errors such as heat flow distortion and soil-plate contact resistance [34]. What is more, measurement source areas of soil heat flux and other terms in energy balance components do not match. So it is important to improve and calibrate soil heat flux measurement for reducing the imbalance of energy in surface flux measurements.

3.3.2. Seasonal Variation of Evapotranspiration

The seasonal variation of actual evapotranspiration (ET) and reference evapotranspiration (ET0) of the P. euphratica forest is shown in Figure 6 and Table 1. ET and ET0 exhibited similar seasonal course. The value of ET is much smaller than ET0 in the soil-freezing season (Figure 7; : 0.12), and averages of ET and ET0 were 0.17 and 1.42 mm·d−1, respectively. In comparison, was high (: 0.75) in the vibrant growing season, and averages of ET and ET0 were 5.09 and 6.81 mm·d−1, respectively. The daily ET over the P. euphratica forest observation site exhibited distinct variation at different stages of growth. ET began to increase rapidly from approximately 1 mm·d−1 in early May to 5 mm·d−1 in mid-June during the germination, flowering, and leaf expansion of the P. euphratica forest. ET then fluctuated between approximately 4 mm·d−1 and 7 mm·d−1 according to the weather and reached its peak (8.38 mm·d−1) on July 7. Afterwards, ET dropped rapidly until the end of November, from 3 mm·d−1 to 0.6 mm·d−1. However, a small spike in the ET distribution occurred during the ablation period due to the higher temperature and increased water level associated with soil thawing. During the growing season (DOY 131–288), ET was 628.6 mm·year−1; otherwise ET was 744.4 mm·year−1 over an entire year. Plant transpiration accounted for most of the ET with weak soli evaporation for the following reason. ET was much less than ET0 in the nongrowing season when plant transpiration stopped, indicating that soil evaporation was weak.

The mean daily ET was 5.09 mm·d−1 over P. euphratica stand during the vibrant growing season in our study. This was similar to results in other desert regions of Northwest China, with a mean value of 4.52 mm·d−1 [18]. And cumulative ET was 744.4 mm·year−1 over the whole year. The magnitude of the annual cumulative ET is less than that observed over the riparian forest (Populus fremontii S. Watson) in the US, where yearly ET was 1095 mm [33, 35]. These differences may be caused by two reasons. First, there was the difference in vegetation coverage at two study sites. The stand density index of the P. euphratica stand observed in our study (148 trees per hectare) was less than what is observed in the western US (430 trees per hectare). Second, the soil evaporation was different at two study sites. Our study was conducted in the hyperarid region, where precipitation was less than 35 mm. In contrast, precipitation over the riparian forest (Populus fremontii S. Watson) in the US (463 mm) was more, so that the soil evaporation was more strong.

4. Conclusions

Energy balance partitioning and evapotranspiration of a P. euphratica forest in an extreme arid region were analyzed using EC method. Several conclusions can be drawn from this study. In summary, we found that the seasonal processes of energy balance partitioning and evapotranspiration of a P. euphratica forest were distinctly related to plant phenology in an extreme arid region. The energy balance closure was approximately 72% over the study period and agreed with results estimated from other eddy covariance measurements of riparian forest. The closure status was higher in the vibrant growing seasons than that in the soil-frozen period. Over the whole year, the ratios of the latent (LE), sensible (), and soil () heat fluxes to net radiation were 0.47, 0.28, and 0.01, respectively. Moreover, characteristics of the energy partitioning were distinct at different stages. In the vibrant growing season, LE sustained higher values ranging from approximately 150 to 200 W·m−2 and exhibited lower values ranging from approximately −20 to 60 W·m−2. Otherwise, was mostly converted into and LE was weak in the nongrowing season. In addition, the yearly evapotranspiration of P. euphratica forest was 744 mm·year−1 with a mean of 5.09 mm·d−1 in the vibrant growing season. And ET was 629 mm·year−1 during the growing season (DOY 131–288). ET over riparian forest from our study site was consistent with that from other desert regions’ stand of Northwest China but was less than that from field in the western US. The difference may be caused by distinct soil evaporation and vegetation coverage. In conclusion, features of the energy and water exchanges over desert riparian P. euphratica forest were basically consistent with general results of other riparian forest.

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this paper.

Acknowledgments

This work was funded by the National Natural Science Foundation of China (no. 31370466, no. 91125002, no. 41401033, and no. 31370467).