Abstract

Evapotranspiration (ET) is an important component of the wetland water budget. Water level declines in Poyang Lake, the largest freshwater lake in China, have caused concerns, especially during low water levels. However, how wetland ET and its partitioning respond to abnormally low water levels is unclear. In this study, wetland ET was estimated with MODIS data and meteorological data. The wetland ET partitioning and its relationship with abnormally low water levels were analyzed for 2000–2013. The results showed that the water evaporation rate () was larger than the land ET rate (); the / ranged from 0.77 to 0.99. When the water level was below 12.8 m, the ET partition ratio was larger than 1, which indicates that wetland ET comes from land surface ET more than water evaporation. The negative standardized water level index (SWI) was used to represent an abnormally low water level in the wetland. Although the monthly wetland ET decreased as the negative SWI decreased, was higher than the average under negative SWI conditions from September to December, when the water level decreased. The abnormally low water level induced more water loss from the land surface, especially when the water level decreased, which reduced the available water resources along the wetland shore.

1. Introduction

Evapotranspiration (ET) is the major factor that controls water balance in nearly all wetlands. Floodplain wetlands are hydrologically complex components of the landscape and are located at the transition between aquatic and terrestrial ecosystems [1, 2]. Water vapor fluxes, especially evapotranspiration, can be highly variable because of intermittent flooding and a diversity of plant communities, which may differ in their water use characteristics and responses to water level changes.

Changes in the proportion of water lost via ET from land surfaces versus evaporation () from water bodies relative to the total wetland ET affect ecological, biogeochemical, and hydrological cycles at multiple temporal and spatial scales in floodplain wetlands [3]. Many studies have highlighted the linkage between ET partitioning versus water level and vegetation cover in wetlands [48]. Changes in water levels for shallow lakes can transform an aquatic environment into an exposed environment, or vice versa. Land surface covers are typically reflected in the shifting dynamics, productivity, and spatial distribution of wetland communities [7, 9, 10]. Vegetation cover and soil conditions are also important factors that affect the rate of water loss [11]. Booth and Loheide [5] showed that changes in ET partitioning occurred because of changes in Leaf Area Index (LAI), plant water stress function, and available storage in the upper soil layer. Decreasing the LAI increased the partitioning of ET toward evaporation. Sánchez-Carrillo et al. [7] indicated that increased vegetation cover did not significantly reduce evapotranspiration but did increase the transpiration to evaporation ratio and the open water area, which increased the wetland evapotranspiration. Huxman et al. [12] identified the partitioning of ET as one of the most important ecohydrological challenges in understanding vegetation dynamics and succession in floodplain wetlands.

Poyang Lake is the largest freshwater lake in China [13, 14]. It is a typical floodplain wetland with dramatic seasonal water level variations that cause inundation areas to vary remarkably from several thousand square kilometers during the summer to less than one thousand square kilometers during the winter [1517]. Poyang Lake wetland is well known and was listed on the Ramsar Convention List of Wetlands of International Importance (the Ramsar Convention 2012) to protect migration bird species and habitats for rare and endangered species [1820]. Unfortunately, the water level and inundation area in Poyang Lake wetland have significantly decreased; moreover, the duration of extreme low water events has increased in the previous decade [13, 15, 21]. These changes have affected the water supply, water quality, wetland, and aquatic habitats [22].

The above problems have caused great concern for researchers. Some studies have focused on investigating recent low water levels and explaining their causes [1315, 23]. Liu et al. [13] revealed a regime shift in the lake size that began in 2006; the trigger for the recent lake level declines was principally attributed to low water levels in the Yangtze River and enhanced outflows from Poyang Lake. Zhang et al. [24] related the low water level to the effect of local catchment and lake-river interactions. Lai et al. [23] indicated that increasing the discharge ability of Poyang Lake into the Yangtze River at low water levels enhanced the drought risk in Poyang Lake. Other studies have examined the effect of recent water level declines on vegetation dynamics and hydrological processes in Poyang Lake wetland [2528]. Zhang et al. [26] analyzed the relationship between hydrological conditions and the distribution of vegetation communities, whereas Yu et al. [28] indicated that recent low water levels have accelerated positive succession in wetland vegetation. Zhao and Liu [27] examined the impact of water level fluctuations on wetland ET, which indicated that the annual ET decreased because of falling water levels. The effect of water level changes on vegetation and hydrological processes has been attributed to water allocation and hydrological gradients in the wetland. Water vapor loss along the shore of Poyang Lake wetland was the primary component of ET when the water level decreased, which affected the water budget and vegetation growth [27]. Some studies have indicated that the long-term exposure of the wetland shore promoted vegetation growth from September to December, during which the water level decreased [25]; this high biomass may increase the land surface ET. However, the long-term exposure of the land surface also causes soil water deficits [7], which typically decrease the land surface ET. How the wetland ET and its partitioning in Poyang Lake wetland respond to abnormally low water levels in different seasons remains unclear.

A better understanding of the response of wetland ET partitioning is essential for the rational utilization of water resources and wetland ecosystem protection. This study estimated wetland ET based on remote sensing data and meteorological data. The objectives were to () identify abnormally low water levels over the previous decade, () analyze variations in wetland ET and its partitioning, and () explore the relationship between ET partitioning and abnormally low water levels.

2. Study Area

Poyang Lake is located in the middle of the Yangtze River Basin (28.2°–30.0°N and 115.5°–116.5°E) (Figure 1). Poyang Lake wetland has a subtropical monsoon climate. The annual mean temperature is 17.1°C, and the annual precipitation is 1620 mm (1958–2013), of which 55% occurs from March to June. Poyang Lake receives inflows from five major rivers (i.e., Xiushui, Ganjiang, Fuhe, Xinjiang, and Raohe) within the basin and discharges into the Yangtze River. The dominant land covers include forest, agricultural fields, grasslands, bare lands, and water surfaces.

Under the Ramsar wetland conservation treaty, Poyang Lake belongs to a wetland type of seasonal or intermittent freshwater lake under inland wetland classifications. Figure 1 shows Poyang Lake wetland with the minimum inundated area. The water levels were measured at Xingzi using the Wusong reference elevation system. Poyang Lake undergoes very large seasonal water level variations, which alters the lake’s inundation area. During the flood season (from April to September), these five contributing tributaries are flooded because of concentrated rainfall, which results in a maximum lake surface area of over 3000 km2 [29]. During the dry season (from October to March), the surface area of the lake shrinks to less than 1000 km2 [13, 15].

3. Methods and Data Source

3.1. Method Overview

Poyang Lake wetland includes water and land surfaces, which have different evaporation or ET processes. Evaporation from the water surface () and ET from the land surface () were both estimated in this study. The wetland ET is partitioned into water surface evaporation and land surface ET using the water surface area percentage () in Poyang Lake wetland. The wetland ET was estimated usingwhere is the spatially averaged wetland ET rate, which is separated into a water surface evaporation component () and a land surface ET component (). Moreover, is the water evaporation rate (, mm/day/m2) multiplied by the water surface area proportion (, %), and is the land surface rate (, mm/day/m2) multiplied by the land surface area proportion (, %). To analyze wetland ET partition, the ET partition ratio is defined as the ratio of to :

Based on the Priestley-Taylor equation, ET can be expressed as follows:where ET is evapotranspiration or evaporation (W m−2), is the Priestley-Taylor parameter, is the net radiation (W m−2), is the soil heat flux (W m−2), is the slope of the saturation vapor pressure at the air temperature (), and is the psychrometric constant (hPa K−1).

The land surface ET was estimated by the -VI triangle method [30] from MODIS data. The modified Priestley-Taylor equation was used on meteorological data to estimate water body evaporation. The lake surface area percentage was extracted by the normalized difference water index (NDWI) from MODIS data. The details of these methods are described as follows.

3.2. Evapotranspiration Estimation for the Land Surface by the Triangle Method

In this study, the -VI triangle method [30] was applied to estimate ET over Poyang Lake wetland. The advantages of the -VI triangle method for ET estimation are that () the complex parameterization of aerodynamic resistance and uncertainty from the replacement of aerodynamic temperature is bypassed, () no ground-based near-surface measurements are required, and () the ET estimation is independent of land-cover types [31]. The Priestley-Taylor parameter, which was proposed by Jiang and Islam [32] and is based on the Priestley-Taylor equation and the existence of a physically meaningful relationship between remotely sensed surface temperatures and the vegetation index, is expressed as in the following equation:where is the Priestly-Taylor coefficient of 1.26 [33]. is the land surface temperature (K). is the lowest temperature in the full vegetation cover, which forms the wet edge in the triangle space of versus NDVI. is the highest temperature at a given NDVI and can be obtained by extrapolating the dry edge to intersect with NDVI = 0. The two boundaries of the /NDVI feature space constitute limiting conditions for the surface fluxes. Three steps are thus required to estimate as follows: () establish the boundaries of the triangle, () interpolate along the dry edge, and () linearly interpolate between the highest and lowest for a given value of NDVI. Then, the wet and dry edges of the triangular space are determined by using the algorithm that was described by Tang et al. [31]. Finally, the value for each pixel in the study area is calculated from the triangular space.

The net radiation () is the sum of the shortwave net radiation and longwave net radiation. The algorithm that was proposed by Bisht et al. [34] was applied to estimate the daily entirely from MODIS products. Then, the daily ET was estimated from the daily net radiation and near-noon instantaneous evaporative flux (EF) ratio. EF was defined as the ratio of ET to available energy () [35]. Missing ET values in the time series from cloud cover in satellite images were reconstructed based on the potential ET (based on the Penman-Monteith equation [31]) and the estimated ET on a clear day [36]. Additional detailed information about soil heat flux estimation and daily ET estimation is available in Zhao and Liu [27]. The calibration and validation of ET estimation in Poyang Lake wetland were also shown in Zhao and Liu [27]. The estimated land surface ET by the triangle method showed good agreement with the observed ET by Lysimeter in Poyang Lake wetland, that is, a regression slope of 0.97 and a correlation coefficient () of 0.49. At the annual scale, the relative error of the estimated ET was 5.4%, which suggests a satisfactory accuracy.

3.3. Evaporation Estimation for Water Bodies by the Modified Priestley-Taylor Equation

The Priestley-Taylor (P-T) coefficient, which is based on the Priestley-Taylor equation (2), was assumed to be constant (i.e., 1.26) for water evaporation estimation studies [3739]. However, research has shown that the P-T coefficient for evaporation over an open water body is not constant because of its association with ambient atmospheric vapor pressure deficits [40]. In this study, the modified Priestley-Taylor equation was developed to estimate water body evaporation through revising the P-T coefficient equation. According to bulk aerodynamic methods and the evaporation observations [41, 42], the P-T coefficient can be calculated as a function of the wind velocity and vapor pressure deficit as follows:where is the wind velocity at a height of 2 m (m s−1) and VPD is the saturation vapor pressure deficit (kPa), which is equivalent to , where is the saturation vapor pressure at a given air temperature and is the actual vapor pressure. These parameters can be obtained via in situ meteorological measurements. Compared with land surfaces, the latent heat flux from water surfaces has a high proportion of net radiation because of the large heat capacity of water. The heat flux that is transferred in water was determined from profile measurements of the water temperature [43]. The meteorological parameters for estimating the water evaporation in this study were obtained from in situ measurements.

3.4. Water Surface Area Extraction

The lake’s surface area was extracted by using the normalized difference water index (NDWI) [44]. The NDWI can be expressed as follows:where NIR is the surface reflectance in the NIR band and Green is the reflectance in the green band. The MODIS product MOD09 provides surface reflectance in the green and NIR bands. The water surface area was delineated by using a thresholding segmentation algorithm that was based on the NDWI [45, 46]. The regression function for the water surface area proportion () and water level () at the Xingzi station was found to be (). Based on this function, the daily lake surface area percentage could be calculated for 2000 to 2013.

3.5. Water Level Index

The standardized water level index (SWI) was used to quantify the magnitude of water level deviations from the mean value over multiple time scales. McKee et al. [47] defined the SWI as follows:where is the normalized monthly water level of year and month , is the multiyear mean of monthly water level for month , and is the standard deviation of monthly water level for month . Negative SWI values represent below-average water levels, whereas positive values represent above-average water levels. The magnitude of the SWI indicates the extent of the water level anomaly from the average.

3.6. Remote Sensing Data

The Terra-MODIS Collection 5 data products that were used in this study included MODIS geolocation (MOD03), the atmospheric profile product (MOD07), the surface reflectance (MOD09), the land surface temperature/emissivity product (MOD11_L2), and the albedo product (MOD43). MOD03 contains solar zenith and azimuth angles and satellite zenith and azimuth angles. MOD07 provides air temperatures and dew point temperatures. MOD09 provides surface reflectance in seven reflective bands and includes values that were corrected for radiometric and atmospheric effects. MOD11_L2 contributes 1 km land surface temperatures and surface emissivity in bands 31 and 32. The MOD43B3 product provides clear-sky observations at 1 km spatial resolution and albedo. To satisfy the spatial resolution and maximize the remote sensing information, all the parameters from the MODIS products were resampled to a 250 m resolution by using a nearest-neighbor algorithm. Because of weather and cloud cover, remote sensing images were limited to available data in the visible band, near-infrared band, and thermal infrared band. Based on cloud detection, images with less than 15% clouds were selected to estimate the regional ET. Many of the 469 images satisfied the above condition in Poyang Lake Basin for 2000–2013. Pixel gaps in the images from cloud cover were gap-filled by using the nearest-neighbor method.

3.7. Evaporation and Meteorological Data

Water evaporation was measured with an E-601 evaporation pan at the Tangyin station (Figure 1). The P-T coefficient in (4) was determined and calibrated based on the evaporation measurements from 2000 to 2012. The pan evaporation data in 2013 were used to validate the estimated evaporation. Meteorological data, including the daily air temperature, relative humidity, solar radiation, precipitation, and 2 m wind speed measurements, from Boyang (29.0°N and 116.7°E) station are available from the China Meteorological Data Sharing Service System at http://data.cma.cn/ for 1958–2013. Water level data from Xiangzi station were obtained from the Hydrological Bureau of the Yangtze River Water Resources Commission for 1960–2013.

Based on meteorological stations around Poyang Lake, the coefficients of variation (CV) (; SD is standard deviation) of air temperature, relative humidity, and wind speed were calculated, which were 6%, 7%, and 25%, respectively. The spatial variability of net radiation in Poyang Lake was neglected due to mixing well in water body. According to sensitivity analysis of the modified Priestley-Taylor equation, the net radiation, air temperature, wind speed, and relative humidity can explain the evaporation variation of 58%, 18%, 13%, and 11%, respectively. In conclusion, the effect of spatial variability of meteorological parameters on evaporation estimation was less than 3%.

4. Results

4.1. Validation

Good agreement exists between the estimated evaporation and the pan evaporation in monthly scale over Poyang Lake ( of 0.67 and RMSE of 18.7 mm month−1; Figure 2). Though evaporation that was estimated by the modified Priestley-Taylor equation was lower than the measurements at high value, which is an obviously improved estimation accuracy compared with the evaporation from the Priestley-Taylor equation, ( of 0.23 and RMSE of 31.4 mm d−1).

4.2. Monthly Variation in Water Level Changes and SWI

Clear seasonal and interannual variations occur with monthly precipitation from Poyang Lake catchment and the water level at Xingzi site for 2000–2013 (Figure 3(a)). Based on seasonal variations in the water level, Poyang Lake wetland can be divided into four distinct periods as follows: a dry season (low water level from January to March), a rising period (increasing water level from April to June), a flood season (high water level from July to September), and a retreating period (decreasing water level from October to December) [48]. Notably, the monthly water level variation was not completely consistent with the observed precipitation. Precipitation was concentrated from March to June, decreasing rapidly thereafter; however, the water level remained high from July to September because of the high water level in Yangtze River and blocked outflows from Poyang Lake [14, 29].

The standardized water level index (SWI) indicates the extent that the water level is above or below average (Figure 3(b)). Abnormally low water levels are identified (negative SWI values) for 42 months (13 months during the retreating period, 12 months during the flood season, 10 months during the rising period, and 7 months during the dry season). The maximum negative SWI was −3.35, which occurred in May 2011. When considering both the duration and magnitude of the negative SWI conditions, typical drought years were confirmed in 2004 (from January to April, lasting 4 months), 2006 (from July to December, lasting 6 months), 2007 (January to June, lasting 6 months), 2009 (January to February and November to December, 4 months), 2011 (February to October, lasting 9 months), and 2013 (August to December, lasting 5 months).

4.3. Variation in Wetland ET and Its Partitioning

Poyang Lake wetland lost water from both evaporation from the water surface () and ET from the land surface (). Seasonal variations in and had a single peak (Figure 4(a)); the maximum ET typically occurred in July, and the minimum occurred in January. The monthly ranged between 33.9 mm and 183.9 mm from 2000 to 2013, which was higher than the (33.0 mm to 154.9 mm). The / was lower during the summer and higher during the winter (Figure 4(b)), ranging from 0.77 to 0.99 (average of 0.87).

The ET partition ratio indicates the partitioning of into water surface and land surface ET over Poyang Lake wetland. When the ET partition ratio is larger than 1, more wetland ET came from the land surface. The average monthly ET partition ratio ranged from 0.11 to 4.8; the ratio was less than 1 from May to September and exceeded 1 in the other months (Figure 5(a)). The ET partition ratio during the dry season was larger than that during the retreating season. The annual mean ET partition ratio ranged from 0.53 to 1.55 (Figure 5(b)); the high ratio values occurred in 2004, 2006, 2007, 2009, 2011, and 2013, corresponding to the drought years in the wetland (as indicated by the SWI values). In those years, the ET partition ratio in the months with abnormally low water levels was substantially above average. Drought years were divided into three categories based on the time of year that a negative monthly SWI value occurred as follows: () the dry and rising seasons (January to May), that is, in 2004 and 2007 (Figure 6(a)); () the flood and retreating seasons (September to December), that is, in 2006 and 2013 (Figure 6(b)); and () the dry and retreating seasons, that is, in 2009 and 2011 (Figure 6(c)). The magnitude and length of time of the high ratio values indicate the intensity and duration of the abnormally low water levels over Poyang Lake wetland.

4.4. Response in Wetland ET to Abnormally Low Water Levels

The magnitude, timing, and duration of low water levels affect wetland ET and its partitioning. The / reflects the effect of land surface changes that are caused by water level changes. The / was not correlated with the water level during the study period (Figure 7(a)). The data were scattered from 0.77 to 0.99 when the water level was below 12 m, whereas the / tended to be constant (0.9) when the water level exceeded 12 m. The same phenomenon was observed in the relationship between the / and the SWI (Figure 7(b)). Notably, the / during the retreating season was always larger than those during the dry and rising seasons during drought years under low water levels and negative SWI values. The response of the land surface ET and wetland ET anomaly with respect to SWI during typical drought months indicates that the land surface ET anomaly was primarily greater than 0 during the retreating season (September to December; Figure 8) and exhibited a negative SWI, whereas the opposite occurred during the dry and rising seasons (January to May) (Figure 8(a)). The anomaly in the wetland ET increased with increasing SWI during the retreating, dry, and rising seasons (Figure 8(b)), which indicates that the more severe the drought, the lower the wetland ET compared with the average.

5. Discussion

The ET in Poyang Lake wetland was separated into evaporation from the water body and ET from the land surface. The annual ET partition ratio ranged from 0.57 to 1.55 and the low water level corresponded to a high ET partition ratio. This relationship has also been observed in other wetlands [7]. For a semiarid wetland in central Spain, evaporation from open water exhibited twice the evapotranspiration rate during high water level periods [7]. As a result, a lower water level (more land surface exposed) can increase the land surface ET’s contribution to the wetland ET and decrease the wetland ET.

During drought years, the land surface ET () increased during the retreating season and decreased during the dry and rising seasons over Poyang Lake wetland. The land surface ET includes evaporation from soil and transpiration from vegetation. Abnormally low water levels occurred during the retreating seasons during 2006, 2009, 2011, and 2013, which increased the exposed land area and exposure time. In those years, the normalized difference vegetation indexes (NDVIs) during the retreating season (from September to December) were substantially higher than mean NDVI from 2000 to 2013 (Figure 9). Ye et al. [25] also indicated that a low water level during autumn improved vegetation growth over the wetland shore region. The high NDVI and biomass for vegetation increased vegetation transpiration [4, 5], which is the main reason for the higher land surface ET during the retreating seasons in drought years.

When the abnormally low water levels occurred during the dry and rising seasons in 2004, 2007, 2009, and 2011, the wetland shore exposure persisted for a few months from autumn to spring. The NDVI from January to May during those drought years was close to the multiyear average value (Figure 10(a)), whereas the ratio of the accumulated precipitation to ET from September to May was remarkably lower than mean ratio from 2000 to 2013 (Figure 10(b)). The ratios of the accumulated precipitation to ET before April during 2004, 2007, and 2009 were less than 1, which caused soil water loss. The soil moisture deficit in the wetland shore limited evaporation from the soil and vegetation sprouting and growth [3, 49]. Therefore, the lower land surface ET rates during the dry and rising seasons in drought years were mainly attributed to reduced evaporation from the soil surface.

6. Conclusions

Severely low water levels have occurred frequently over Poyang Lake wetland over the previous decade. In this study, wetland ET was separated into land ET and water evaporation in Poyang Lake wetland. The land ET was estimated by using the -VI triangle method with MODIS data, and water evaporation was estimated with the Priestley-Taylor method by using meteorological data. The water evaporation rate () was 1.01–1.30 times larger than the land ET rate (). When considering the inundated area in Poyang Lake wetland, the ET partition ratio was larger than 1 when the water level was below 12.8 m, and the water area percentage was less than 46%. The monthly wetland ET decreased as the negative standardized water level index (SWI) decreased. However, the was higher than the average under abnormally low water levels during the retreating season, and the opposite results occurred during the dry and rising seasons.

The abnormally low water level decreased the wetland ET over Poyang Lake wetland; however, this relationship does not suggest more water resources in wetland shore. The low water level, in terms of the reduced inundation area and increased exposure of the land surface, increased the ET partition ratio, which indicates that more water loss came from the wetland shore, decreasing the available water in the soil of the exposed land.

Nomenclature

Abbreviations and Symbols
:Water surface area proportion (%)
:Water surface evaporation component,
:Land surface ET component,
CV:Coefficient of variation
:Evaporation (mm/time/area)
:Actual vapor pressure (kPa)
:Saturation vapor pressure (kPa)
ET:Evapotranspiration (mm/time/area)
ET anomaly:Difference of ET to mean ET (mm/time/area)
ET partition ratio:Ratio of to ,
:ET from the land surface (mm/time/area)
:Wetland ET (mm/time/area)
:Evaporation from the water surface (mm/time/area)
:Soil heat flux (W m−2)
Green:Reflectance in the green band
NDVI:Normalized difference vegetation index
NDWI:Normalized difference water index
NIR:Reflectance in the NIR band
:Net radiation (W m−2)
SD:Standard deviation
SWI:Standardized water level index
:Air temperature (K)
:Highest temperature at a given NDVI (K)
:Lowest temperature in the full vegetation cover (K)
:Land surface temperature (K)
:Wind velocity at a height of 2 m (m s−1)
VI:Vegetation index
VPD:Saturation vapor pressure deficit (kPa)
:Water level (m)
:Normalized monthly water level of year and month (m)
:Multiyear mean of monthly mean water level for month (m)
:Psychrometric constant (kPa K−1)
:Slope of the vapor pressure curve (kPa K−1)
:Standard deviation of monthly water level for month
:Priestley-Taylor coefficient
:Priestly-Taylor coefficient of 1.26.

Competing Interests

The authors declare that they have no competing interests.

Acknowledgments

This work was jointly supported by the 973 Project (2012CB417003) and the National Natural Science Foundation of China (no. 41430855 and no. 41471022).