Abstract

To accurately estimate the magnitude and seasonal dynamics of evapotranspiration (ET) over an important a swamp meadow in the Fenghuoshan permafrost region, we employed the Food and Agriculture Organization- (FAO-) Penman-Monteith (P-M) model. The model was also used to investigate changes in the crop coefficient (), which was calculated as the ratio of the measured actual ET ( from the eddy covariance (EC) system) to the reference ET ( from the P-M model). The results indicated a reference ET of 900 mm/year from the swamp meadow ecosystem, which was significantly higher than the actual ET (426 mm/year). The reference ET peaked from April to July, while the actual ET was primarily in growing season. The value of exhibited significant seasonal variations within the range 0.3–1.0 with a mean of 0.55 during the growing season. The daily showed a linear increase with and and a linear decrease with the VPD. With respect to the biotic factors, the biomass exhibited a significant positive correlation with . Thus, a daily model is developed as a function of the VPD, , , and biomass.

1. Introduction

The hydrologic balance of terrestrial ecosystems is an important determinant of ecosystem structure, function, and productivity [1]. Evapotranspiration (ET), which is the second largest water flux in the terrestrial hydrologic cycle, plays an important role in the maintenance of the water and energy balance of the ground surface [2]. In addition, ET processes are closely related to vegetation conditions, the ecophysiological processes of plants, soil environments, and micrometeorological characteristics. Thus, ET is a pivotal water exchange process in the soil-plant-atmosphere continuum (SPAC). Therefore, accurate ET estimation is significant not only to the regulation and management of hydrologic cycles in ecosystems but also to scientific decisions regarding local ecological construction and production activities in agriculture and animal husbandry [35].

Domestic and international scholars have investigated land surface ET for more than 200 years and have developed a number of fitting models and observational methods for ecosystem evapotranspiration [69]. Currently, the eddy covariance (EC) system is the most extensively used and sophisticated micrometeorological approach. Ecosystem ET may be continuously monitored with the EC system without damage to the soil or vegetation. The Penman equation, which utilizes conventional meteorological data, is the most influential fitting equation for predicting ET. This equation was proposed in 1948 when Penman proposed a formula for calculating ET from a water surface, which considers radiant energy, air saturation deficit, wind speed, and other ET-influencing factors. Based on prior research, Monteith and Unsworth modified the Penman equation and proposed the Penman-Monteith (P-M) model, which is a canopy ET model that has been extensively applied [7]. In 1990s, the Food and Agriculture Organization (FAO) of the United Nations amended the P-M model for the estimation of actual ET from farmland and grassland and recommended this modified model as a standard method for calculating ET [10].

In the P-M model recommended by the FAO (the FAO-P-M model), a crop coefficient () has been introduced to correct reference ET values and to accurately estimate actual ET values in specific ecosystems. Studies have reported that is related to various biotic factors, such as crop type and growth stage [5, 11]. For homogeneous vegetation, which consists of plants in a particular growth stage, may be regarded as a constant. For example, in the FAO-56 method, values are assigned based on a meadow’s growth stage; thus, values for the initial, mid, and late season growth stages are 0.4, 1.05, and 0.85, respectively [10]. Recent studies have demonstrated that exhibits specific variations and is affected by radiant energy, moisture levels, and other environmental factors [12, 13]. The value is the key to using the FAO-P-M model to accurately estimate the actual ET of an ecosystem. Thus, the study of the patterns and characteristics of changes in in natural meadowland can improve the accuracy and simplicity of the calculation of the ecological water demand of meadowlands and provide a theoretical foundation for the grazing production of these lands.

The Sanjiangyuan Region (i.e., the source of the Yangtze, Yellow, and Mekong rivers, well known as the “water tower of Asia”) is located in the hinterlands of the Qinghai-Tibet Plateau, which plays a pivotal role in the global hydrologic cycle and the global water balance. Swamp meadow ecosystem, which is one of the most extensively distributed types of vegetation in the Sanjiangyuan Region, covers an area of approximately 127.63 km2 [14]. The swamp meadow is also a unique natural landscape and one of the most important grassland resources for grazing [15]. Few studies have explored the ET characteristics of meadow ecosystems [16]; as a result, studies on applicable models for accurately and conveniently evaluating ET in the swamp meadow ecosystem also remain scarce.

This study aimed to achieve the following objectives: (1) use the FAO-P-M model and eddy covariance system to explore the dynamics of the actual ET and reference ET changes in the swamp meadow ecosystem of the Sanjiangyuan Region and (2) derive a suitable -based model of the swamp meadow ecosystem by determining how values vary with changes in meadow climate and vegetation and by establishing the relationships between and the factors (including both the environmental and biotic factors).

2. Material and Methods

2.1. Study Site

The Fenghuoshan permafrost region of the Sanjiangyuan Region (92°50′–93°3′ E and 34°40′–34°48′ N), located in the upper watershed of the Zuomo Xikong River, a tributary of the Yangtze River, with typical swamp meadow ecosystem, was chosen as the system of the study. The region covers a total area of 127.63 km2 at elevations ranging from 4680 to 5360 m, belonging to an arid climate area and without glacial or multiyear snowpack. The mean annual (1973 to 2005) air temperature, the highest air temperature, the lowest air temperature, precipitation, evaporation, and relative humidity are −5.2°C, 23.2°C, −37.7°C, 290.9 mm, 1316.9 mm, and 57%, respectively. The mean annual ground temperature ranges from −1.5°C to 4.0°C and the main frozen soil depth ranges from 50 m to 120 m. The permafrost table varies between 0.8 and 2.5 m depth. The average annual duration of sunshine is 2462.7 h, and the total radiation () received per year ranges from 6000 to 7000 MJ m−2. The swamp meadow features high vegetation cover composed of short and densely distributed plants.

The vegetation primarily consists of meadow species, such as Stipa aliena, Kobresia tibetica, Festuca sp., Carex atrofusca, Leontopodium leontopetaloides, and other alpine plants. The meadow land, which is primarily used for meadow (January to May and September to December), serves as an important grazing resource in the plateau region. The aboveground biomass of the swamp meadow increases in early May and peaks in August, with a multiyear average biomass of approximately 350 g/m2. The study area receives a yearly precipitation of approximately 600 mm. This precipitation is primarily concentrated between May and September, and changes in the soil water content (SWC) are subject to precipitation-related effects.

2.2. Measurements
2.2.1. Flux Measurements

EC observation system was placed in the center of the study area. This area features flat terrain, which enables unobscured observation and provides a sufficiently large “fetch” to satisfy the required physical conditions for eddy and meteorological observations. ET was measured using a H2O/CO2 infrared gas analyzer (Li-7500, Li-Cor, USA) and a three-dimensional sonic anemometer (CSAT3, CSI, USA) with a sensor placed 2.2 m above the ground. A data logger (CR5000, Campbell Inc.) was employed to continuously record water vapor flux data and output average values at 15 min intervals. The actual ET () was measured by the EC system in this study.

2.2.2. Meteorological Measurements

Environmental variables were continuously measured at this site. Air temperature (), humidity, and actual vapor pressures were measured at 110 cm and 220 cm (HMP45C, CSI); soil temperature was measured at 0, 5, 10, 20, 30, and 50 cm depths (thermocouple); solar radiation () and net radiation () were measured at 150 cm (CNR-1, Kipp and Zonen, Netherlands); a horizontal wind-speed sensor (014A and 034A-L, CSI) was attached at 110 cm and 220 cm to measure horizontal wind speeds; precipitation was measured at 70 cm (TE525MM, CSI); SWC was measured at 5, 20, and 50 cm depths (TDR soil moisture sensor; CS615, CSI); and the soil heat flux was set at a depth of 2 cm (HFT-3, CSI). The signals were sampled at 10 Hz, and 15 min mean data were logged by the data logger (CR23X, CSI).

2.2.3. Observations of Vegetation Data

Aboveground biomass in the study area was measured semimonthly during the growing season. A harvesting method was adopted to obtain these measurements. In particular, five 0.5 m × 0.5 m quadrants were randomly selected; within each quadrant, plants were harvested by cutting the plants at ground level and loading the plants into sampling bags. Each sample was numbered, rapidly transported to a laboratory, and dried at 65°C until a constant weight was obtained for each sample. These weights were subsequently converted into g/m2.

2.3. The FAO-P-M Model
2.3.1. The Calculation of Reference ET

The reference crop comprises a well-managed short grass of uniform height (8–15 cm) that completely covers the ground and grows lushly across open ground without experiencing water stress. The reference crop ET refers to the ET of the reference crop under these conditions [10]. In this study, the standardized FAO-P-M model was used to calculate the actual ET in the swamp meadow ecosystem. This model includes two components—reference ET and —as indicated in the following equation:

In this equation, is the actual ET of an ecosystem in the FAO-P-M model, is the crop coefficient, and is the reference ET. The reference surface consists of well-watered grass with a height of 12 cm and a fixed surface resistance of 70 s/m. The applicable equation is expressed aswhere is the slope of the saturation vapor pressure curve at the examined temperature; is the net radiation; is the soil heat flux; is the psychrometric constant; is the mean daily temperature at a height of 2 m; is the wind speed at a height of 2 m; and and are the saturation vapor pressure and actual vapor pressure, respectively. and are calculated aswhere is the specific heat at constant pressure, which has a fixed value of 1.013 × 10−3 MJ/kg°C; is the atmospheric pressure; 0.662 is the ratio of the molecular weight of water vapor to the molecular weight of dry air; and is the latent heat of vaporization, which represents the energy per unit volume of water required to convert water into steam under ambient temperature and pressure conditions. The study meadow area is located on a plateau zone at an average elevation of 3250 m. Using (4), the atmospheric pressure () can be calculated from this average elevation ( m):

2.3.2. The Crop Coefficient ()

The crop coefficient was obtained from the ratio of the reference ET in 2004 to the ET obtained from EC system observations in 2004 and fitted with environmental factors to obtain an empirical equation, which was eventually validated by 2005 data. In the FAO-56 approach, the growth of meadow plants was divided into different growth stages [10]. The study area of the current investigation features an alpine climate; in this area, plant growth typically begins in late April, and the growing season extends from May to September.

2.4. Evaluation Methods

In this study, the slope, linear correlation coefficient (), relative root-mean-squared error (RRMSE), index of agreement (IA), and coefficient of determination (CD) were utilized to examine the extent of the differences between the observed and calculated values of ET and to statistically analyze the accuracy of the adopted fitting approach from multiple perspectives. The slope and were calculated by Origin 8.0. RRMSE, IA, and CD values were computed as [13, 17]where indicates the th observed value, represents the th simulated or estimated value, represents the number of samples, and is the mean observed value. RRMSE is a measure of the relative magnitude of residuals; smaller values of RRMSE indicate better model calculation results. The IA can be used to evaluate the correlation between the observed values and the simulated values; the closer the IA is to 1, the more closely the fitted values match the observed values. The CD can be used to measure the dispersion of the simulated values from the mean observed values; a CD greater than 0.8 indicates satisfactory simulation results [18].

3. Results and Discussion

3.1. Seasonal Variations in Reference ET () and Actual ET

As indicated in Figure 1, the reference ET derived from the P-M model was significantly higher than the actual ET measured by the EC system. The reference ET () refers to the reference capacity of moisture diffusion from the ecosystem to the atmosphere under specific environmental conditions; thus, represents the maximum ET. The annual reference ET of the examined swamp meadow was 900.2 mm, which was significantly higher than the actual ET of 425.8 mm. These parameters also exhibited significant seasonal variations. began to significantly increase in March and attained peak levels from April to July; however, it began to decline in August. The total during the growing season (from May to September) was 487.1 mm, which represented 54% of the annual cumulative . Conversely, seasonal variations in actual ET () in the swamp meadow lagged behind variations in reference ET. The was concentrated in the growing season; between May and September the was 337.1 mm, which represented 79% of the annual cumulative .

According to its definition, the reference ET () was only affected by climatic parameters [10]. As shown in (2), the net radiation (), vapor pressure deficit (VPD, ), air temperature (), and wind speed () determined the reference ET. The annual variation of these environmental factors was shown in Figure 2. As the maximum of , VPD, and occurred during April to June (Figure 2), reached its maximum value at a similar period (March to June). However, with the canopy development in growing season, the effect of biotic factors will exert increasing influence on the actual ET. Therefore, with the rapidly increasing in growing season, the difference between and became smaller (Figure 1).

Moreover, the reference ET in growing season (487.1 mm) was lower than the value of precipitation (579.7 mm). The ratio of to was 84% during growing season. This result was similar to a nonirrigated pasture [19]; however, it was much lower than some grassland ecosystems [9, 20]. In these grasslands, the reference ET was significantly higher than the precipitation in growing season.

Comparing to these grasslands [9, 20], the lower may be caused by the unique climate in this plateau. The net radiation () in the Qinghai-Tibetan Plateau was much lower than that for lowland grasslands [21], despite the high incident solar radiation (), due to the fact that the net long-wave radiation in alpine region is much higher than that for lowland regions [22]. In this swamp meadow, the was much lower than the , and the average and values were 7.7 and 16.6 mol/(m2 d), respectively (Figure 2(a)). The low energy available for water evaporation may limit the reference ET to some extent. Moreover, VPD and were meteorological variables for control of water vapour exchange between the atmosphere and vegetation [23]. Gu et al. [24] indicated that low VPD and due to the frequent precipitation and the altitude were characteristics of the climate in the swamp meadow in Northeast Qinghai-Tibetan Plateau. Figure 2(b) showed the variation of VPD and in this study site. The maximum values of daily mean VPD and were only 1.5 kPa and 14.5°C, which were greatly lower than many other grasslands with the maximum VPD ranging from about 2 to 5 kPa and with the maximum ranging from about 20 to 30°C [2527]. The low VPD and might imply the weak driving power and are considered as the factors to restrict the reference ET in growing season.

3.2. The Seasonal Variation in

Due to the difference between reference ET and actual ET, the correct determination of , which is the ratio of daily ET to , is important for the accurate estimation of the actual ET. In this study, exhibited a gradual rise from mid-April to mid-June and remained at a high level in July and August (Figure 3). The value of rapidly declined in September and was less than 0.2 by the end of October. During the growing season, which extended from May to September, exhibited a mean value of 0.55 and fluctuated within the range of 0.30–0.92. This range is consistent with the range of values reported in FAO-56 (0.30–1.05) [10] but is significantly higher than the range of values observed for a typical steppe (0.32–0.68) [28] and temperate desert steppe (0.02–0.50) [13] of Inner Mongolia. These results indicated that the for the swamp meadow features more adequate moisture and better plant growth conditions than arid and semiarid steppe ecosystems.

3.3. Factors That Impact
3.3.1. The Effect of Environmental Factors on

Previous studies have focused on values related to local climate conditions [12, 13]. In this meadow, the Pearson product-moment correlation coefficients for the relationships between and its main environmental factors were computed on a daily scale (Table 1). This computation showed that the VPD and were critical for controlling . also showed a significant relation to , whereas and SWC at a depth of 5 cm were not significant at the 95% confidence level.

The responses of to changes in VPD, , and are illustrated in Figure 4. To reduce or offset the errors associated with the values, VPD, , and were grouped into bins with the following criteria: 0.2 kPa for VPD, 1 MJ/m2 d for , and 2°C for . increased linearly in response to an increase in and and in response to a decrease in VPD (Figure 4). L. Zhou and G. S. Zhou [29] also obtained a similar result from a study of a reed marsh. However, our result was different from a study in a temperate desert steppe [13], which demonstrated that the soil water content (SWC) was the most important factor for determining . In this study, the root-layer (0–10 cm) SWC fluctuated within the range 0.2–0.6 m3/m3 (Figure 1(c)), which was much better than the SWC at desert steppe in Inner Mongolia (0.05–0.15 m3/m3). The higher soil moist in this swamp meadow might weaken the SWC impact on .

Daily for the swamp meadow can be expressed by an empirical equation from the correlation and regression analyses of the relationships between and its statistically significant environmental variables:

3.3.2. The Effects of Biotic Factors on

In addition to environmental factors, the plant species composition, growth conditions, and other biotic factors can significantly affect [3, 30, 31]. Biomass directly reflects the growth status of a biological community. In this study, the sampled plants began to grow at the end of April, which caused an increase in the biomass that primarily began in early May. The biomass rapidly increased during June and August and attained a peak in late August (with a maximum value of 353.1 g/m2 in the study year). After August, the biomass rapidly declined as plants gradually died (Figure 5(a)). The linear regression analysis of the relationship between and the biomass demonstrated a significant positive correlation between and the biomass; however, the significance of this correlation changed during different growth stages. As indicated in Figure 5(b), the biomass increased and rapidly increased from 0.45 to 0.93 from May to middle August. After biomass attained its peak, began to decline and decreased to 0.39 in late October, while no significant decrease in was observed from September to October, because only decreased from 0.54 to 0.39 during this period.

The following linear equations describe how increases with biomass:

Considering both biotic and environmental factors, daily actual ET can be calculated as follows:

3.4. Test for Modeling Evapotranspiration

Equation (9) was validated using the EC data collected during the 2005 growing season. The model performed well during the growing season. The simulated daily ET (ETmod) followed the seasonal trend of the EC measured ET () (Figure 6). The slope (=0.95) of the regression line between the measured and simulated ET passes through the origin close to 1, with , a RRMSE of 0.23 mm/d, and an IA of 0.92.

In this study, the abiotic and biotic factors enabled a better estimation of evapotranspiration. In Table 2, the slope, , IA, and CD of ETmod1 are greater and the RRMES is smaller than those of ETmod2. These results indicated that relative to the fitting method that accounts for meteorological variables, the fitting method that comprehensively considers meteorological factors and vegetation conditions can be used to derive simulated ET results that are closer to the actual observed ET results (Table 2).

4. Conclusions

The FAO-Penman-Monteith model was used to simulate ET over the swamp meadow on the Sanjiangyuan Region. The reference ET was estimated by the FAO-P-M model and showed difference from the actual ET measured by the EC system. The crop coefficients () were used to estimate . The average value of during growing season was 0.55 (ranging from 0.30 to 0.92). VPD, , and explained the majority of the daily variation in , which showed a linear increase with an increase in and and a linear decrease with an increase in VPD. We also considered the impact of biomass on , which showed a significant positive effect on . A daily empirical model driven by VPD, , , and biomass was developed to estimate daily actual ET using the FAO-56 approach and the P-M model for reference ET. This ET model was validated against 2005 growing season data for this meadow and demonstrated a suitable and consistent performance between the simulated and measured ET.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgments

This study was supported by the National Natural Science Foundation of China (31070433) and the Japan-China Research Cooperative Program (2010DFA31290).