- Open Access
Effects of changes in precipitation on energy and water balance in a Eurasian meadow steppe
Ecological Processes volume 8, Article number: 17 (2019)
Water such as precipitation is the most critical environment driver of ecosystem processes and functions in semi-arid regions. Frequency and intensity of drought and transient waterlogging are expected to increase in the meadow steppe in northeastern China. Using a 4-year dataset of eddy covariance flux measurements, ground measurements of biomass, phenology, and meteorological conditions, we investigated the changes in energy fluxes at multiple temporal scales and under different precipitation regimes.
The meadow steppe was latent heat (LE) dominated when soil water content was > 0.3 m3 m−3, but switched to sensible heat (H) dominated status when soil water content fell below 0.3 m3 m−3. LE dominated the energy exchange of the meadow grasslands on a yearly basis. Intensive precipitation had a profound impact on water-energy balance that could reduce the damages of drought by elevating deep soil moisture. The influence of LE on waterlogging depended on timing, with increased LE at the beginning of growing season and decreased LE after waterlogging. Spring and summer droughts resulted in different energy partitioning between latent and sensible heat energies, with spring drought dramatically decreased the LE fraction due to the change in water. In contrast, summer drought had little impact on LE due to the sufficient water input from large precipitation events at the beginning of the growing season.
There existed great seasonal and interannual variabilities in energy balance and partitioning in the meadow steppe over the 4-year study period, which were strongly influenced by changes in precipitation. The water loss through latent heat was more sensitive to spring drought than to summer drought, while summer drought had negligible impact on LE. Waterlogging contributed to LE by enhancing its values during and after the waterlogged periods at the beginning of the growing season in a dry year, but lowering its value after the waterlogged periods in growing season.
Water is the most influential resource in ecosystems and societies in arid and semi-arid areas (Chen et al. 2013; Feng et al. 2016; Jaeger et al. 2017; Qi et al. 2017). In drylands, the spatio-temporal distribution and availability of water resources have been the major focuses of terrestrial ecosystem studies at local, regional and global levels (Groisman et al. 2017; Monier et al. 2017; Soja and Groisman 2018). In terrestrial ecosystems, hydrological cycling is largely motivated through energy transfer because of the close coupling of water and heat. Precipitation, via either rainfall or snow, is the largest flux term of water budget and can lead to changes in soil moisture and latent heat (LE), sensible heat (H), soil heat flux (G), and heat storage (S) by influencing plant transpiration and soil evaporation (Rodrigues et al. 2013). Knowledge about coupled ecosystem water cycling and energy balance can not only be applied to explain the integrated water-energy exchange between ground and atmosphere but also be used to address water scarcity and maintain ecosystem functions.
Rapid changes in land use and land cover and global warming have produced an increasing incidence and intensity of climatic extremes such as frequent and heavy precipitation events (Eade et al. 2012; IPCC 2013; Rey et al. 2017) and more intensive and prolonged drought. Changes in precipitation regimes contribute to higher uncertainty in ecosystem water-energy balance that requires further investigation of biomes and habitat. For instance, Minderlein and Menzel (2015) reported direct influences of precipitation on soil moisture, evapotranspiration (ET), and energy balance of shrub-grass ecosystems in the semi-arid areas of northern Mongolia. Nevertheless, there remains limited knowledge and consensus on the relationships between changes in precipitation and water-energy balance in grassland ecosystems.
Climate change, together with human disturbances, has had a large influence on ecological processes in the Eurasian Steppes, such as carbon assimilation and emission and energy and hydrological cycling. The meadow steppe appears to be the most sensitive ecosystem type in the grassland biome because of its relatively high biodiversity and high demands for water (Wang et al. 2007). Using a long-term research site of a meadow steppe in Northeast China, we explore the changes over time in water and energy fluxes in response to changes in climate, especially precipitation. Here, we focus particularly on the changes of major energy flux terms and the overall balances under different water availability. Based on 4 years of eddy covariance measurements, we first emphasize the distribution of energy and water fluxes and the possible ecological and climatological controls for a Songnen meadow steppe. Precipitation and its extremes, especially drought and transient waterlogging events, are examined for their roles in regulating the energy balance. Two major scientific questions are as follows: (1) What are the roles of shifting precipitation on water-energy balance of a meadow steppe? (2) What are the differences in how each energy flux term responds to spring/summer drought and transient waterlogging?
Materials and methods
This study was conducted at the Songnen Grassland Ecology Field Station of the Northeast Normal University (NENU) in Changling, Jilin, China (123° 30′ E, 44° 35′ N, 171 m a.s.l.). The site is representative of a temperate meadow steppe dominated by perennial grasses, including Leymus chinensis and Phragmitis communis. The study area has a temperate, semi-arid continental monsoonal climate, with a cold-dry spring and a warm-wet summer. The mean annual temperature is 5 °C, with the maximum and minimum temperatures being 39.2 °C and − 33.9 °C, respectively. The frost-free period is approximately 130 days to 165 days. The precipitation varied greatly within and between years, with an annual mean of 350 mm, 80% of which falls between June and August. The annual pan evaporation varies from 1200 to 1600 mm, which is approximately three to four times higher than the mean annual precipitation. The growing season runs from May through September. Main soil types in this area are alkaline soil, chernozem, and saline soil.
To evaluate the energy and water exchanges between the land surface and the atmosphere, an open-path eddy-covariance (EC) flux tower was deployed in 2007 as part of the US-China Carbon Consortium (USCCC). A three-dimensional ultrasonic anemometer (CSAT3, Campbell Scientific Inc. (CSI)) and a CO2/H2O infrared gas analyzer (LI-7500, LI-COR) were mounted 2 m above the ground. Fluxes were calculated as the mean covariance of vertical wind speed fluctuation and the scalar fluctuation. Downward fluxes were indicated as negative and upward as positive. All raw data was sampled at 10 Hz and logged using a CR3000 data-logger (CSI). Footprint analysis using the flux source area model (FSAM) suggested a footprint in the prevalent wind direction (150–240°) extended to approximately 6.5 m to 77.8 m during unstable conditions (Monin-Obukhov length (L) < 0) and to about 11.3 m to 229.9 m during stable conditions (L > 0) at night, representing 90% of the total flux.
Meteorological measurements were taken starting in May 2007. Net radiation (Rn) was measured using a four-component net radiometer (CNR-1, Kipp & Zonen) at 2 m above the ground. Photosynthetically active radiation (PAR), precipitation (P), wind speed (WS), and wind direction (WD) were measured by a quantum sensor (LI190SB, LI-COR), a tipping-bucket rain gauge (TE525MM, CSI), and a propeller anemometer (034B-L Met One Windset, CSI), respectively. Air temperature (Ta) and relative humidity (Rh) were measured at 2, 4, and 6 m with HMP45C probes (Vaisala). Soil temperature profiles (Ts) were measured at depths of 0.05, 0.1, and 0.3 m below the ground by thermistors (107-L, CSI). Soil heat flux (G) was measured at a depth of 0.05 m in three separate locations (HFT-3, REBS), and soil water content (SWC) was measured at a depth of 0.1 m by a time-domain reflectometry probe (CS616, CSI), which was buried horizontally in the soil. Additionally, soil water potential (SWP) was estimated using watermark probes at 0.1 and 0.3 m (257-L, CSI). All micrometeorological data were recorded and averaged or summed over a 30-min interval by a separate CR3000 data-logger.
Biometric measurements of aboveground biomass (AGB), vegetation height (Vt), and leaf area index (LAI) were conducted monthly during the growing season (May to September). Within a radius of 200 m around the EC tower, we placed 12 0.5 m2 quadrats for measuring Vt and LAI by a plant canopy analyzer (LAI-2000, LI-COR), while AGB was destructively sampled in 12 additional 0.5 m2 quadrats.
Data processing and gap filling
The EC flux data was processed using EdiRe. The planar fit rotation and Webb–Pearman–Leaning (WPL) density correction were applied to the flux calculations (Webb et al. 1980). Data from stable nocturnal periods were excluded, specifically when the friction velocity (u*) was < 0.15 m s−1. Anomalous or spurious data due to sensor malfunction, sensor maintenance, precipitation events, IRGA calibration, and power failure were also rejected. Consequently, 25.4%, 20.9%, 27.1%, and 30.6% of the data obtained during growing season in 2007, 2008, 2009, and 2010, respectively, were discarded. The introduced data gaps were filled with following the methods: for gaps of less than 2 h, linear interpolation was used by averaging the fluxes before and after the gaps; for larger data gaps, we used both the empirical relationships and look-up tables (Falge et al. 2001).
Calculation of soil heat storage
The Rn was partitioned into convective H for heating the atmosphere, LE for evaporating water from plants and the soil, G for heating the soil, and S for heat storage in the soil. Because the canopy and air heat storage are expected to be negligible in short canopies with minimal biomass, the canopy heat storage was not included in our energy balance algorithm:
where soil heat storage (S) was calculated as follows:
where Ts is the average soil temperature (K) above heat flux plates, t is time (Δt = 1800 s), d is a depth between the heat flux plate and soil surface, and Cs is the soil heat capacity calculated from Shao et al. (2014):
where ρb is soil bulk density, ρw is the density of water, Cd and Cw are the specific heat capacity of dry mineral soil (Cd = 890 J kg−1 K− 1) and soil water (Cw = 4190 J kg−1 K−1), respectively, and θv is volumetric SWC (%). The heat storage of air and organic matter in the soil were neglected due to their small quantities.
The energy balance ratio (EBR) was calculated as a ratio of the turbulent fluxes (LE + H) and the radiative energy fluxes (Rn–G–S):
Canopy conductance (gc) was calculated using the inverted form of the Penman–Monteith equation (Monteith and Unsworth 1990):
where ρ is air density, Cp is the specific heat of the air, Δ is the change of saturation vapor pressure with temperature, γ is the psychrometric constant, D is the vapor press deficit of air, and ga is the aerodynamic conductance of the air layer between the canopy and the flux measurement height. ga was calculated flowing Monteith and Unsworth (1990) as follows:
where u is wind speed and u* is friction velocity.
To calculate the Priestley-Taylor coefficient (LE/LEeq), the equilibrium LE flux (LEeq) was determined using (Priestley and Taylor 1972):
The LEeq is dependent on the Rn and temperature. Lower and higher values indicate evaporation rates that are lower and higher than the equilibrium rate, respectively (Wilson et al. 2002b). The sensitivity of evapotranspiration to stomatal control and the degree of aerodynamic coupling between vegetation and the atmosphere was expressed by the decoupling factor (Ω), obtained by the equation (Jarvis and McNaughton 1986; Runkle et al. 2014):
The value of Ω ranges from 0 to 1. When it approaches 0, the vegetation and the atmosphere are fully aerodynamically coupled and the evapotranspiration proceeds at rates imposed by vapor pressure deficit; when it approaches 1, the vegetation and the atmosphere are completely aerodynamically decoupled and the evapotranspiration is controlled by Rn (Goldberg and Bernhofer 2001).
The mean Ta in 2007–2010 was 8.23, 6.83, 5.64, and 7.91 °C, respectively (Table 1). For the growing season, the mean Ta in 2007 (21.26 °C) and 2010 (20.67 °C) were higher than that in 2008 (19.76 °C) and 2009 (19.88 °C). The maximum Ta occurred in July of 2007 and 2008, August of 2009, and June of 2010, while the minimum Ta occurred in January in all 4 years. For the non-growing seasons (October–April) of 2007–2008, 2008–2009, 2009–2010, and 2010–2011, the mean Ta was − 2.82, − 3.38, − 6.69, and − 5.62 °C, respectively. The Ts closely followed the changes of Ta (Fig. 1a). However, there was no significant difference between Ts (21.15 °C) and Ta (21.26 °C) in the growing season of 2007.
Precipitation amounts and distribution were significantly different during the 4 years. Annual precipitation from 2007 to 2010 was 207.9, 384.2, 281.5, and 265.6 mm, respectively. These years can be classified as dry (2007, 2009, and 2010) or wet years (2008) based on to the precipitation trending in the Songnen meadow steppe for the past 50 years (Fig. 2). The driest period was from DOY181 to DOY279 in 2007, with the first rain of > 3 mm. Precipitation or snowfall throughout the initial 6 months was minimal. In 2008, the meadow steppe received 85% more precipitation and experienced 18 more single precipitation pulses than in 2007. The first precipitation in 2008 (March 22nd) was 3 months earlier than in 2007, and the rainy season therein lasted longer and ended on November 15th. In 2009 and 2010, precipitation started in mid-April, and it ended in mid-October and mid-November, respectively. Total precipitation of March–May in dry year 2007 accounted for only 8.8% and 5.1% of that in 2009 and 2010, respectively, indicating that the meadow steppe experienced a serious spring drought in 2007. On the other hand, June–August precipitation was low in 2010 (90.7 mm) compared to that in 2007 (170 mm) and 2008 (231.9 mm). In the meadow steppe, transient waterlogging occurs usually after a single precipitation of > 30 mm, and the waterlogging condition often lasts for 2–3 days. Six transient waterlogging events were detected during the study period: DOY181, 190, 212 in 2007; DOY187, 214 in 2008; DOY125 in 2010.
Seasonal variation of SWC closely followed the changes in precipitation. The temporal change of SWC was clearly asymmetrical within any of the 4 years (Fig. 1e). A long dry period in 2007 resulted in severe water deficit in the spring and early summer when the SWC fell below 0.2 m3 m−3 in the upper soil layer. In contrast, frequent precipitation events during the summer of 2008 resulted in a consistently higher SWC. The highest SWC occurred in the early spring of 2009 and 2010, which declined sharply afterwards. SWC decreased to its minimum before July in 2010, followed by a summer drought. The SWC recovered in June of 2009 after a series of continuous precipitation events. A large precipitation event combined with several isolated precipitation events led to a temporary SWC peak during the 2010 summer compared to the previous summer. Interestingly, SWC at 10 cm depth did not drop below 0.4 m3 m−3 during the water-stressed periods in 2007 and 2009 (Fig. 2).
A unimodal pattern of growing season AGB and LAI showed an overall higher AGB and earlier peak in 2008 than in the dry years (Fig. 3). Peak value of AGB in August 2007 was ~ 68% of that in July and August of 2009 and 2010. LAI was highest in 2008 and lowest in 2007.
Energy balance closure
The energy balance enclosures using half-hourly and daily values produced reasonable results for all 4 years (Table 2). The slope of the enclosure was < 1 for all years, with an average value (± SD) of 0.75 (± 0.06), intercept of 10.48–21.17 W m−2, and R2 value of 0.84–0.98. For the growing season, the slope values were similar in 2008, 2009, and 2010. The average slope value increased to 0.85 (± 0.09), ranging from 0.74 to 0.93, when using the daily sums instead of the half-hourly data to calculate the regression coefficients of [LE + H] against [Rn–G–S].
The EBR provides an overall evaluation of energy closure at long temporal scale by averaging over random errors in the flux measurement. The mean values of EBR for the four growing seasons were 0.94 (2007), 0.83 (2008), 0.85 (2009), and 0.92 (2010), respectively, and increased to 1.00, 0.97, 0.91, and 0.93 at annual scale for the four consecutive years (Table 1).
Intra- and interannual variations in energy fluxes
The seasonal changes of daily Rn, LE, H, and G were similar (Fig. 1b–d). Rn showed maximum values of 15.67–20.86 MJ m−2 day−1 in July and the minimum values of from − 1.03–3.56 MJ m−2 day−1 in December among the 4 years. Seasonal changes in LE matched well with those of Rn, which gradually increased from May, peaked in July–August, and gradually decreased until soil was frozen in late October through February. During the spring drought in 2007, LE in March–May (3.51 MJ m−2 day−1) was 70% of the normal values (4.65–5.08 MJ m−2 day−1). LE during the peak growing season of 2010 (July–August) was 7.31 MJ m−2 day−1 and remained lower during 2007–2009 (~ 6.4 MJ m−2 day−1) despite the summer drought of 2010. Surprisingly, H showed a bimodal seasonal change, which increased at the beginning of a year, peaked in May, and began to decrease in June even though Rn continued its high values. Interestingly, H started to increase again in August and reached the second peak in September. As expected, G values were low, especially during the wintertime. At the beginning of March, G switched from negative to positive and continued to increase, reaching a maximum in May for 2007 (3.39 MJ m−2 day−1), 2009 (2.16 MJ m−2 day−1), and 2010 (2.15 MJ m−2 day−1), but in June of 2008 (2.10 MJ m−2 day−1).
Monthly LE value was higher than H value during the growing season in 2008 (Fig. 4). In dry years, the duration of higher LE over H appeared short (2 months) in 2007 (spring drought), mediate (3 months) in 2009, and long (4 months) in 2010, with different starting days in July, June, and May, respectively. For accumulated G, however, a positive G (22.72 MJ m−2 month−1) was observed in March 2008, which was − 2.00 MJ m−2 month−1 in 2009 and − 7.84 MJ m−2 month−1 in 2010. More importantly, G in May (67.91 MJ m−2 month−1) and June (132.91 MJ m−2 month−1) of 2007 was about twice the values of 2008 (35.21/71.93 MJ m−2 month−1), 2009 (39.85/63.93 MJ m−2 month−1), and 2010 (41.44/76.21 MJ m−2 month−1).
Sensible heat (H) was a major factor of Rn at annual scale, but LE was more influential during the growing seasons, with the mean LE/Rn ratio 0.46 to 0.53. LE/Rn values remained low during the non-growing season, with mean values in 0.12, 0.18, and 0.19 for 2007, 2008, and 2009, respectively (Table 1). The peak LE/Rn in 2007 arrived half a month later than that in 2010 (Fig. 5). Throughout the spring drought period in 2007, LE/Rn remained low until late summer. In contrast, LE/Rn in 2010 remained relatively high from the beginning of the growing season, when soil moisture was high due to the large precipitation input. The seasonal changes of H/Rn seemed to be a mirror image of LE/Rn in all 4 years. The lowest value of H/Rn was observed during the growing season when LE/Rn was the highest, while H/Rn began increasing in late September. The annual mean G/Rn was higher in 2007 (0.171) than in other years (0.0045–0.055) (Table 1). Although there was a higher G/Rn at the beginning of growing season (May–June) in 2007, no significant difference was found among LEs of the 4 years despite a severe summer drought in July–August of 2010.
Biophysical controls of energy partitioning
The Bowen ratio (β, H/LE) during the growing season varied from 0.21 to 4.63 (Fig. 6b). Due to a faster increase in H than that in LE during the spring, midday means of β over the growing season exhibited a concave shape, with β being > 1 during spring and autumn, but ~ 1 in summer. Over the growing season, midday means of β amounted to 1.39, 0.95, 1.29, and 1.04 for 2007–2010, respectively. The low LE and high H during the non-growing season resulted in a dramatic increase in the mean daytime β. Mean daily β varied from 0.02 to 12.85 over the 4-year study period, with the largest value found in the non-growing season of 2007 (DOY310), and the lowest in the growing season of 2010 (DOY233). The average of daily β during the spring drought period was 1.97, which was higher than that in the summer drought period (1.17). The gc value ranged from 0.02 to 19.67 mm s−1 during the growing season, with a peak in July of 2008, 2009, and 2010, but not in 2007, when the maximum gc was delayed to August (Fig. 6c). In 2007, LAI maintained a higher level after July while gc declined during the same time. During the growing season, the Ω value amounted a peak of 0.73 (DOY229), 0.81 (DOY214), 0.75 (DOY204), and 0.56 (DOY213), respectively, for the 4 years. Its seasonal changes appeared less pronounced compared to other measures of the fluxes, with the change in 2007 very different from that in other 3 years (Fig. 6).
Impact of precipitation intensity on water and energy balance
In this study, P and ET appeared to be more or less balanced in the wet year. Differing from previous reports on P vs ET in neighboring typical steppes on the Mongolia Plateau (Chen et al. 2009; Foken et al. 2006), the meadow steppe in our study showed an overall higher cumulative ET than cumulative P in dry years (2007, 2009, and 2010), suggesting there were other sources of water to enrich the soil (Fig. 7). The mean annual ET (409.5 mm) (Table 1) in the meadow steppe lies in the higher end of ET values (163~481 mm) for grasslands worldwide (Aires et al. 2008; Krishnan et al. 2012; Li et al. 2007; Ryu et al. 2008; Zha et al. 2010). However, ET (334.5 mm) in 2007—a spring drought—fell below the mean value. The ratio of ET/P varied between 1.05 and 1.71 for all study years (Table 1), and was higher than those in the Mediterranean grasslands (0.4–0.87) (Aires et al. 2008; Ryu et al. 2008), American semi-arid prairies (0.72–1.20) (Krishnan et al. 2012), northern temperate grasslands (0.98) (Hao et al. 2007; Wever et al. 2002), and alpine Kobresia meadows (0.60) (Li et al. 2013), but similar to those in the alpine wetland meadows (1.27) (Hu et al. 2009).
Intensity of precipitation events had a crucial effect on water-energy balance and ecosystem recovery after drought (Cabral et al. 2010; Heisler-White et al. 2008; Knapp et al. 2008). Differences in seasonal ET/P in the meadow steppe suggested that non-growing season P from snow would supplement water resources required for the following growing season. The relative contribution of temporally isolated non-growing season precipitation events to water availability in the next growing season is considered a critical determinant of ecological processes (Eamus et al. 2001; Ma et al. 2013). Previous studies have indicated that freezing and thawing have important influences on ET and surface energy fluxes (Chen et al. 2011; Gu et al. 2005; Yao et al. 2008), by increasing SWC and ET from melted snow and frozen soil (Zhang et al. 2014; Zhang et al. 2003). Although cumulative P was lower than ET at annual scale, the monthly P could sometimes be higher than ET in the meadow steppe during the growing season. For example, heavy and frequent precipitation events in July of 2008 brought ~ 132.66 mm precipitation and 90.8 mm ET at our site, yielding a peak soil water content of 0.91. Similarly, although the maximum soil water content was equivalent to that of 2008, P was 86.5% lower and ET was 13.8% higher following the consecutive precipitation events in July 2009 compared to those in July 2008 (Figs. 1 and 7), suggesting a proportion of precipitation infiltrated to recharge groundwater in 2008. Heavy precipitation events in May 2010 (P = 105.1 mm, ET = 58.4 mm) alleviated summer drought in the growing season through water infiltration into deep soil under similar maximum soil water content (0.89 mm) in May 2009 (P = 14.5 mm, ET = 45.6 mm). Therefore, AGB between 2009 and the summer drought year of 2010 remained similar, likely because both the low infiltration rate in the alkali-saline soil of the meadow steppe and the limited precipitation input from short and sporadic showers could substantially reduce soil moisture. However, low ET/P in July 2007 (0.74), in July (0.68), and August (0.68) of 2008 and in May 2010 (0.56) (Fig. 7) indicated that infiltration might have happened after occasional and large precipitation events (Knapp et al. 2008; Thomey et al. 2011).
In our study, LE was apparently influenced by the process of temporary waterlogging caused by precipitation > 30 mm in a day (Fig. 8). LE increased because of low vegetation cover and high groundwater level at the beginning of the growing season, which is consistent with the findings of spring drought influence on evapotranspiration and water use efficiency in the same steppe (Dong et al. 2011). However, LE declined in response to the 3-day waterlogging at the growing season peak of 2007 and 2008, which could possibly be attributed to stomatal closure and reduction in both leaf photosynthesis and root respiration. Clearly, the severity of waterlogging effects on LE depend on the developmental stage of the grass.
Influence of drought on energy balance and partitioning
The Songnen meadow steppe exhibited the great intra- and interannual variability of energy partitioning into H and LE, which also differed greatly between the dry and wet years. The meadow steppe was LE-dominated when soil water content was > 0.3 m3 m−3, but switched to H-dominated status when soil water content fell below 0.3 m3 m−3 (Fig. 9). The timing of the shift from an H-dominated to a LE-dominated energy balance depended on the magnitude of the first precipitation event (> 30 mm day−1) or accumulated precipitation within five consecutive days over 30 mm, which enabled grass seed germination with SWC exceeding 0.3 m3 m−3. Similarly, the switch from a LE-dominated to an H-dominated ecosystem happened in concert with the start of biomass senescence when SWC fell below 0.3 m3 m−3. However, LE was dominant in energy balance even though SWC dropped to 0.3 m3 m−3 in the growing season of 2010 (since DOY200), which likely was due to the lag effect of a large precipitation event at early stage. In 2008–2010, the change in energy partitioning occurred in May–June. However, the spring drought led to a shift from H to LE occurring in mid-July in 2007. The increase in SWC and vegetation growth during the growing season altered solar heating of the surface via a decrease in albedo that likely was driven by the dark soil and vegetation surfaces after the precipitation events (Krishnan et al. 2012; Thompson et al. 2004).
Interannual climate, especially the occurrence, duration, and extent of drought, can significantly reduce LE in this site. On the one hand, limited precipitation could cause a decrease in soil moisture and limit soil water recharge, thus decreasing LE. On the other hand, droughts with high temperature and radiation would result in a reduction in transpiration with lower stomatal conductance of foliage and canopy (Chen et al. 2009). For instance, spring drought in 2007 resulted in a higher LE over H in the beginning of growing season (May–June), but the summer drought of 2010 did not change the interannual relation of higher LE over H in June–August. The growth of plants was triggered by the monsoon precipitation for all grasses. Delayed and low LE in spring drought year appeared to be a direct consequence of slow grass germination and late starting of transpiration. The main reason for higher LE during summer drought meanwhile might be related to the soil water availability. Several large precipitation pulses at the beginning of the growing season could have replenished the soil profile to allow grass to reach a biomass peak despite the summer drought in 2010. The results suggested that LE was more sensitive to spring drought than to summer drought, and there was hardly a negative consequence of summer drought to LE in the meadow steppe.
Impacts of biotic and climate on LE were further examined for other flux terms and key parameters, including LE/LEeq, gc, and Ω. Evapotranspiration (i.e., latent heat) is mainly composed of plant transpiration (Tr), soil evaporation (Es), and the negligible proportion of canopy intercepted precipitation (I) in grasslands that can be inferred from LE/LEeq, gc, and Ω. LE/LEeq is conventionally used to estimate daily evapotranspiration (Wilson et al. 2002a). LE/LEeq of spring drought (0.53) and summer drought (0.61) in this study were typical events for the semi-arid grassland LE/LEeq (0.6–0.7) (Meyers 2001), indicating that the restriction of the water supply for LE over the meadow steppe was less compared to other grasslands. Decreased gc is usually attributed to low LAI during the growing season (Baldocchi et al. 2004) and decreasing stem-root hydraulic conductivity under conditions of low soil water content (Rodrigues et al. 2013). According to Dong et al. (2011), Es was significantly higher than Tr under spring drought conditions, during which Es released more than 85% of available soil water in May 2007. During the spring drought period, the low gc did not reflect the low transpiration rate. On the contrary, higher gc during summer drought may have been influenced by high LAI in 2007. Ω evaluated the response of LE to air humidity and net radiation changes (Jarvis and Mcnaughton 1986). Ω of growing seasons was higher in summer drought (0.34) than in spring drought (0.24), suggesting that LE was highly sensitive to decreasing air humidity. Moreover, similar Ω (0.33) in summer drought years and wet years further proved the significant effect of spring drought in increasing sensitivity of LE to low air humidity. In summary, although LE was not limited by water stress in the meadow steppe (4.58 MJ m−2 day−1 in spring drought and 5.43 MJ m−2 day−1 in summer drought), spring drought had an even stronger effect on LE than summer drought did.
The proportion of individual surface energy balance component is a result of complex, long-term interactions between biogeochemical cycling, disturbance and climate, as well as short-term interactions of plant physiology and the dynamics of the atmospheric boundary layer (Aires et al. 2008; Baldocchi et al. 2004; Shao et al. 2012; Sun et al. 2010). The annual values of LE/Rn in our study (0.28–0.40) are comparable to the 0.31–0.35 that were reported for an annual grassland (Baldocchi et al. 2004), but lower than those in a tallgrass prairie (0.48–0.58) (Burba and Verma 2005) and in a Mediterranean C4 grassland (0.37–0.55) (Jongen et al. 2011). For the growing season, LE/Rn (0.49) was similar to that in an alpine meadow (Li et al. 2013), and higher than in temperate semi-arid grasslands in North America (Krishnan et al. 2012) and C3/C4 grassland in Southern Portugal (Aires et al. 2008), but less than the values reported in the Mongolian steppe (Li et al. 2007). When compared among the 4 years, the spring drought in 2007 caused the highest H/Rn and G/Rn but the lowest LE/Rn. The short-term water deficit at the initial stage of growing season could be partially responsible for the low LAI and AGB values. The cumulative available energy was 639.8 MJ m−2 in the extended drought period of 2010 (July–August), 6.6% lower than that of the same month in 2007. However, there was no obvious decrease in LE/Rn (Fig. 5), suggesting that the magnitude of LE was not constrained by the availability of energy but rather by the vegetation cover at dry period. As a result of earlier water inputs from precipitation, grass germinated earlier in the season, and the growing period was lengthened to maintain a stable and high AGB.
Soil heat flux (G) played an important role in energy partitioning at our site, which is similar to previous findings at other typical steppes (e.g., Shao et al. 2008). Soil heat flux tended to increase with Rn and H. A low canopy cover in the spring drought year jointly led to the highest G values (0.356 MJ m−2 day−1). Other than for the spring drought year, G/Rn was < 0.1 at the meadow steppe, which was lower than those reported in other grasslands (Aires et al. 2008) but similar to those at the Mongolian steppe and temperate semi-arid prairies (Krishnan et al. 2012; Li et al. 2007). More importantly, soil heat flux had wider diel variations in the dry year of 2007 than in other years (Fig. 10). The large seasonal fluctuations also reflected the seasonal variation in soil thermal conductivity, which were regulated by precipitation, soil moisture, and seasonal changes in vegetation coverage that directly determine value Rn (Heusinkveld et al. 2004; Rodrigues et al. 2013). A higher G was expected in a wet year (2008) because soil thermal conductivity increased with soil moisture (Hillel et al. 2005). However, the shading effect of vegetation cover likely offset the positive effect of soil thermal conductivity and resulted in lower G (0.69 MJ m−2 day−1). Finally, a higher S was found in spring drought year. A critical take-home message is that, if S is omitted from our analysis, G value at the topsoil may underestimate the surface heat flux by > 50% at our site (Ochsner et al. 2007; Shao et al. 2008).
The lack of energy closure is a common feature in most eddy flux sites. According to the FLUXNET study (Wilson et al. 2002a), the EBR for all vegetation types of site-years ranged from 0.53 to 0.99, with a mean of 0.79 ± 0.01. For the ChinaFLUX sites, the annual EBR ranged from 0.58 to 1.00, with a mean of 0.83 (Chen et al. 2011). In this study, the mean energy imbalance (0.89 ± 0.05) at meadow steppe is generally within, or higher than, the acceptable range reported within the FLUXNET and ChinaFLUX sites, which was higher compared with typical steppe 0.83 ± 0.07 (Chen et al. 2009).
Changes in precipitation resulted in frequent drought and waterlogging in meadow steppe of northeastern China. We measured the energy and water fluxes over a meadow steppe in Songnen Plain over a 4-year study period and found that spring and summer droughts resulted in a rapid shift in energy partitioning due to the change of water availability. However, spring drought produced more severe effects than summer drought. The rise of the groundwater level in our study region usually followed a large precipitation event or a series of continuous precipitation pulses that have profound influence on ecosystem water balance. Waterlogging influenced LE by enhancing its values during and after the waterlogged periods at the beginning of growing season in a dry year, but lowering them after the waterlogged periods in growing season. Disentangling the roles of different water table heights in controlling the energy balance, particularly LE remains a challenging research topic and should be prioritized in future studies.
Energy balance ratio
- Es :
- G :
Soil heat flux
- g c :
- H :
Sensible heat flux
Leaf area index
- LEeq :
Equilibrium LE flux
- P :
Photosynthetically active radiation
- R h :
- R n :
Soil water content
Soil water potential
- T a :
- T r :
- T s :
Soil temperature profiles
Vapor pressure deficit
Aires LM, Pio CA, Pereira JS (2008) The effect of drought on energy and water vapour exchange above a mediterranean C3/C4 grassland in Southern Portugal. Agric For Meteorol 148(4):565–579
Baldocchi DD, Xu L, Kiang N (2004) How plant functional-type, weather, seasonal drought, and soil physical properties alter water and energy fluxes of an oak–grass savanna and an annual grassland. Agric For Meteorol 123(1–2):13–39
Burba GG, Verma SB (2005) Seasonal and interannual variability in evapotranspiration of native tallgrass prairie and cultivated wheat ecosystems. Agric For Meteorol 135(1–4):190–201
Cabral OMR, Rocha HR, Gash JHC, Ligo MAV, Freitas HC, Tatsch JD (2010) The energy and water balance of a Eucalyptus plantation in southeast Brazil. J Hydrol 388(3–4):208–216
Chen J, Wan S, Henebry G, Qi J, Gutman G, Sun G, Kappas M (2013) Dryland East Asia: land dynamics amid social and climate change. Higher Education Press, Beijing
Chen N, Guan D, Jin C, Wang A, Wu J, Yuan F (2011) Influences of snow event on energy balance over temperate meadow in dormant season based on eddy covariance measurements. J Hydrol 399(1–2):100–107
Chen S, Chen J, Lin G, Zhang W, Miao H, Wei L, Huang J, Han X (2009) Energy balance and partition in Inner Mongolia steppe ecosystems with different land use types. Agric For Meteorol 149(11):1800–1809
Dong G, Guo J, Chen J, Sun G, Gao S, Hu L, Wang Y (2011) Effects of spring drought on carbon sequestration, evapotranspiration and water use efficiency in the songnen meadow steppe in Northeast China. Ecohydrology 4(2):211–224
Eade R, Hamilton E, Smith DM, Graham RJ, Scaife AA (2012) Forecasting the number of extreme daily events out to a decade ahead. J Geophys Res Atmos. https://doi.org/10.1029/2012JD018015
Eamus D, Hutley LB, O'Grady AP (2001) Daily and seasonal patterns of carbon and water fluxes above a north Australian savanna. Tree Physiol 21(12–13):977–988
Falge E, Baldocchi D, Olson R, Anthoni P, Aubinet M, Bernhofer C, Burba G, Ceulemans R, Clement R, Dolman H, Granier A, Gross P, Grunwald T, Hollinger D, Jensen NO, Katul G, Keronen P, Kowalski A, Lai CT, Law BE, Meyers T, Moncrieff H, Moors E, Munger JW, Pilegaard K, Rannik U, Rebmann C, Suyker A, Tenhunen J, Tu K, Verma S, Vesala T, Wilson K, Wofsy S (2001) Gap filling strategies for long term energy flux data sets. Agric For Meteorol 107(1):71–77
Feng X, Fu B, Piao S, Wang S, Ciais P, Zeng Z, Lü Y, Zeng Y, Li Y, Jiang X (2016) Revegetation in China’s Loess Plateau is approaching sustainable water resource limits. Nat Clim Chang 6(11):1019–1022
Foken T, Wimmer F, Mauder M, Thomas C, Liebethal C (2006) Some aspects of the energy balance closure problem. Atmos Chem Phys 6(12):4395–4402
Goldberg V, Bernhofer C (2001) Quantifying the coupling degree between land surface and the atmospheric boundary layer with the coupled vegetation-atmosphere model HIRVAC. Ann Geophys 19(5):581–587
Groisman P, Shugart H, Kicklighter D, Henebry G, Tchebakova N, Maksyutov S, Monier E, Gutman G, Gulev S, Qi J, Prishchepov A, Kukavskaya E, Porfiriev B, Shiklomanov A, Loboda T, Shiklomanov N, Nghiem S, Bergen K, Albrechtová J, Chen J, Shahgedanova M, Shvidenko A, Speranskaya N, Soja A, Kd B, Bulygina O, McCarty J, Zhuang Q, Zolina O (2017) Northern Eurasia Future Initiative (NEFI): facing the challenges and pathways of global change in the twenty-first century. Prog Earth Planetary Sci 4(1):41
Gu S, Tang Y, Cui X, Kato T, Du M, Li Y, Zhao X (2005) Energy exchange between the atmosphere and a meadow ecosystem on the Qinghai–Tibetan Plateau. Agric For Meteorol 129(3–4):175–185
Hao Y, Wang Y, Huang X, Cui X, Zhou X, Wang S, Niu H, Jiang G (2007) Seasonal and interannual variation in water vapor and energy exchange over a typical steppe in Inner Mongolia, China. Agric For Meteorol 146(1–2):57–69
Heisler-White J, Knapp A, Kelly E (2008) Increasing precipitation event size increases aboveground net primary productivity in a semi-arid grassland. Oecologia 158(1):129–140
Heusinkveld BG, Jacobs AFG, Holtslag AAM, Berkowicz SM (2004) Surface energy balance closure in an arid region: role of soil heat flux. Agric For Meteorol 122(1–2):21–37
Hillel D, Hatfield JH, Powlson DS, Rosenzweig C, Scow KM, Singer MJ, Sparks DL (2005) Thermal properties and processes. In: H Daniel (ed) Encyclopedia of Soils in the Environment. Elsevier, Oxford, pp 156–163
Hu Z, Yu G, Zhou Y, Sun X, Li Y, Shi P, Wang Y, Song X, Zheng Z, Zhang L, Li S (2009) Partitioning of evapotranspiration and its controls in four grassland ecosystems: application of a two-source model. Agric For Meteorol 149(9):1410–1420
IPCC (2013) Summary for policymakers. In: Climate change 2013: the physical science basis. Cambridge University Press, New York, pp 1–30
Jaeger WK, Amos A, Bigelow DP, Chang H, Conklin DR, Haggerty R, Langpap C, Moore K, Mote PW, Nolin AW (2017) Finding water scarcity amid abundance using human-natural system models. Proc Natl Acad Sci U S A 114(45):11884–11889
Jarvis PG, McNaughton KG (1986) Stomatal control of transpiration: scaling up from leaf to region. Adv Ecol Res 15:1–49
Jongen M, Pereira JS, Aires LMI, Pio CA (2011) The effects of drought and timing of precipitation on the inter-annual variation in ecosystem-atmosphere exchange in a Mediterranean grassland. Agric For Meteorol 151(5):595–606
Knapp AK, Beier C, Briske DD, Classen AT, Luo Y, Reichstein M, Smith MD, Smith SD, Bell JE, Fay PA (2008) Consequences of more extreme precipitation regimes for terrestrial ecosystems. BioScience 58(9):811–821
Krishnan P, Meyers TP, Scott RL, Kennedy L, Heuer M (2012) Energy exchange and evapotranspiration over two temperate semi-arid grasslands in North America. Agric For Meteorol 153:31–44
Li J, Jiang S, Wang B, Jiang W-w, Tang Y-h, Du M-y, Gu S (2013) Evapotranspiration and its energy exchange in alpine meadow ecosystem on the Qinghai-Tibetan Plateau. J Integr Agric 12(8):1396–1401
Li S-G, Asanuma J, Kotani A, Davaa G, Oyunbaatar D (2007) Evapotranspiration from a Mongolian steppe under grazing and its environmental constraints. J Hydrol 333(1):133–143
Ma X, Huete A, Yu Q, Coupe NR, Davies K, Broich M, Ratana P, Beringer J, Hutley LB, Cleverly J, Boulain N, Eamus D (2013) Spatial patterns and temporal dynamics in savanna vegetation phenology across the North Australian Tropical Transect. Remote Sens Environ 139:97–115
Meyers TP (2001) A comparison of summertime water and CO2 fluxes over rangeland for well watered and drought conditions. Agric For Meteorol 106(3):205–214
Minderlein S, Menzel L (2015) Evapotranspiration and energy balance dynamics of a semi-arid mountainous steppe and shrubland site in Northern Mongolia. Environ Earth Sci 73(2):1–17
Monier E, Kicklighter DW, Sokolov AP, Zhuang Q, Sokolik IN, Lawford R, Kappas M, Paltsev SV, Groisman PY (2017) A review of and perspectives on global change modeling for Northern Eurasia. Environ Res Lett 12(8):083001
Monteith JL, Unsworth MH (1990) Principles of environmental physics, 2nd edn. Edward Arnold, New York
Ochsner TE, Sauer TJ, Horton R (2007) Soil heat storage measurements in energy balance studies. Agron J 99(1):311
Priestley CHB, Taylor RJ (1972) On the assessment of surface heat flux and evaporation using large-scale parameters. Mon Weather Rev 100(2):81–92
Qi J, Xin X, John R, Groisman P, Chen J (2017) Understanding livestock production and sustainability of grassland ecosystems in the Asian Dryland Belt. Ecol Process 6(1):22
Rey A, Oyonarte C, Morán-López T, Raimundo J, Pegoraro E (2017) Changes in soil moisture predict soil carbon losses upon rewetting in a perennial semiarid steppe in SE Spain. Geoderma 287:135–146
Rodrigues TR, de Paulo SR, Novais JWZ, Curado LFA, Nogueira JS, de Oliveira RG, Lobo FA, Vourlitis GL (2013) Temporal patterns of energy balance for a Brazilian tropical savanna under contrasting seasonal conditions. Int J Atmos Sci https://doi.org/10.1155/2013/326010
Runkle BRK, Wille C, Gažovič M, Wilmking M, Kutzbach L (2014) The surface energy balance and its drivers in a boreal peatland fen of northwestern Russia. J Hydrol 511:359–373
Ryu Y, Baldocchi DD, Ma S, Hehn T (2008) Interannual variability of evapotranspiration and energy exchange over an annual grassland in California. J Geophys Res Atmos 113(D9):D09104
Shao C, Chen J, Li L, Xu W, Chen S, Gwen T, Xu J, Zhang W (2008) Spatial variability in soil heat flux at three Inner Mongolia steppe ecosystems. Agric For Meteorol 148(10):1433–1443
Shao C, Chen J, Li L, Zhang L (2012) Ecosystem responses to mowing manipulations in an arid Inner Mongolia steppe: an energy perspective. J Arid Environ 82:1–10
Shao C, Li L, Dong G, Chen J (2014) Spatial variation of net radiation and its contribution to energy balance closures in grassland ecosystems. Ecol Process 3(1):7
Soja AJ, Groisman PY (2018) Earth science and the integral climatic and socio-economic drivers of change across northern Eurasia: the NEESPI legacy and future direction. Environ Res Lett 13(4):040401
Sun G, Noormets A, Gavazzi MJ, McNulty SG, Chen J, Domec JC, King JS, Amatya DM, Skaggs RW (2010) Energy and water balance of two contrasting loblolly pine plantations on the lower coastal plain of North Carolina, USA. For Ecol Manag 259(7):1299–1310
Thomey ML, Collins SL, Vargas R, Johnson JE, Brown RF, Natvig DO, Friggens MT (2011) Effect of precipitation variability on net primary production and soil respiration in a Chihuahuan Desert grassland. Glob Chang Biol 17(4):1505–1515
Thompson CC, Beringer J, Chapin FS, Mcguire AD (2004) Structural complexity and land-surface energy exchange along a gradient from arctic tundra to boreal forest. J Veg Sci 15(3):397–406
Wang Y, Zhou G, Wang Y (2007) Modeling responses of the meadow steppe dominated by Leymus chinensis to climate change. Clim Chang 82(3–4):437–452
Webb EK, Pearman GI, Leuning R (1980) Corrections of flux measurements for density effects due to heat and water vapor transfer. Q J Roy Meteor Soc 106(447):85–100
Wever LA, Flanagan LB, Carlson PJ (2002) Seasonal and interannual variation in evapotranspiration, energy balance and surface conductance in a northern temperate grassland. Agric For Meteorol 112(1):31–49
Wilson KB, Baldocchi DD, Aubinet M, Berbigier P, Bernhofer C, Dolman H, Falge E, Field C, Goldstein A, Granier A, Grelle A, Halldor T, Hollinger D, Katul G, Law BE, Lindroth A, Meyers T, Moncrieff J, Monson R, Oechel W, Tenhunen J, Valentini R, Verma S, Vesala T, Wofsy S (2002a) Energy partitioning between latent and sensible heat flux during the warm season at FLUXNET sites. Water Resour Res 38(12):1294–1305
Wilson KB, Goldstein AH, Falge E, Aubinet M, Baldocchi DD, Berbigier P, Bernhofer C, Ceulemans R, Dolman H, Field CB (2002b) Energy balance closure at FLUXNET sites. Agric For Meteorol 113(1):223–243
Yao J, Zhao L, Ding Y, Gu L, Jiao K, Qiao Y, Wang Y (2008) The surface energy budget and evapotranspiration in the Tanggula region on the Tibetan Plateau. Cold Reg Sci Technol 52(3):326–340
Zha T, Barr AG, van der Kamp G, Black TA, McCaughey JH, Flanagan LB (2010) Interannual variation of evapotranspiration from forest and grassland ecosystems in western Canada in relation to drought. Agric For Meteorol 150(11):1476–1484
Zhang S-Y, Li X-Y, Ma Y-J, Zhao G-Q, Li L, Chen J, Jiang Z-Y, Huang Y-M (2014) Interannual and seasonal variability in evapotranspiration and energy partitioning over the alpine riparian shrub Myricaria squamosa Desv. on Qinghai–Tibet plateau. Cold Reg Sci Technol 102:8–20
Zhang Y, Ohata T, Kadota T (2003) Land-surface hydrological processes in the permafrost region of the eastern Tibetan Plateau. J Hydrol 283(1):41–56
This study was supported by the Major State Research Development Program of China (2016YFC0500600, 2017YFE0104500), Natural Science Foundation of China (31800512, 31870466), and the US-China Carbon Consortium (USCCC). We are grateful to all members who worked at the Changling station. Sincere thanks go to Ge Sun, Jingfeng Xiao, and Haiqiang Guo for their useful comments and valuable suggestions on early versions of the manuscript. We also greatly appreciate the careful reviews by the anonymous reviewers and editing of Kristine Blakeslee.
Availability of data and materials
The datasets generated and/or analyzed during the current study are available from the corresponding author upon request.
Ethics approval and consent to participate
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Chen, J., Shao, C., Jiang, S. et al. Effects of changes in precipitation on energy and water balance in a Eurasian meadow steppe. Ecol Process 8, 17 (2019). https://doi.org/10.1186/s13717-019-0170-z
- Energy balance
- Energy partition
- Precipitation change
- Eddy covariance