www.biogeosciences.net/5/421/2008/ © Author(s) 2008. This work is distributed under the Creative Commons Attribution 3.0 License.

Using a six year data set of eddy covariance flux measurements of sensible and latent heat, soil heat flux, net radiation, above-ground phytomass and meteorological driving forces energy partitioning was investigated at a temperate mountain grassland managed as a hay meadow in the Stubai Valley (Austria). The main findings of the study were: (i) Energy partitioning was dominated by latent heat, followed by sensible heat and the soil heat flux; (ii) When compared to standard environmental forcings, the amount of green plant matter, which due to three cuts varied considerably during the vegetation period, explained similar, and partially larger, fractions of the variability in energy partitioning; (iii) There were little, if any, indications of water stress effects on energy partitioning, despite reductions in soil water availability in combination with high evaporative demand, e.g. during the summer drought of 2003.


Introduction
The partitioning of surface net radiation, that is the net balance of incident and reflected/emitted radiative energy, into the fluxes of sensible and latent heat and heat storage controls the degree of coupling between the land and the atmosphere (Oke, 1987).Thereby, surface energy partitioning determines several properties of the planetary boundary layer, such as the surface temperature (Schneider and Eugster, 2005), boundary layer height (Pan and Mahrt, 1987), cloud development and convective precipitation (Pielke, 2001;Ray et al., 2003).Since the planetary boundary layer is where humans spend most of their lives, it is of interest to understand the variability of land surface energy partitioning (e.g.Wilson et al., 2002a), the causes for this variability (e.g.Gu et al., Correspondence to: A. Hammerle (Albin.Hammerle@uibk.ac.at) 2006), and also how global changes in climate and land use affect energy partitioning (e.g.Rosset et al., 2001;Beringer et al., 2005;Schneider and Eugster, 2005).In addition, there is a close link between ecosystem energy partitioning and the cycling of carbon and nutrients (Meyers, 2001;Wilson et al., 2003).
A convenient framework for analysing the surface energy partitioning is the energy balance theorem (e.g.Oke, 1987), which states that the net radiation (R net ) of a uniform ecosystem patch equals the sum of the fluxes of latent (LE), sensible (H ) heat and heat storage ( S), i.e.R net = LE + H + S. (1) Rearranging Eq. ( 1), expanding terms 1-3 and neglecting the contribution of heat storage in metabolic processes and above-ground phytomass to S but accounting for the soil heat flux (G; described in detail in Sect.2.3) yields: s − G = ρc p g aH (T S − T a ) + λg tV (e s − e a ) P a , (2) where K ↓ and L ↓ represent the incident short-and long-wave radiation fluxes, respectively; α stands for the albedo, ε s and T s for the long-wave emissivity and (bulk) surface temperature, and σ represents the Stefan-Boltzmann constant; ρ and c p are the density and specific heat of air, g aH is the aerodynamic conductance for heat, and T a stands for air temperature at some reference height; the latent heat of vaporisation is given by λ, g tV represents the total conductance for water vapour, e s and e a refer to the saturation vapour pressure of the surface and the actual air vapour pressure, respectively, and P a is the air pressure.
The four main variables under biotic control in Eq. ( 2) are G, α, g aH and g tV , the latter deriving from 1/(1/g sV + 1/g aV ), where g sV and g aV are the surface and aerodynamic conductance for water vapour, respectively.The soil heat flux, G, depends on the amount of radiation absorbed at the soil surface and is thus inversely related to the amount of above-ground plant area.The albedo, α, depends on the amount and spatial organisation of the above-ground phytoelements and the optical properties of the vegetation and the soil surface (Ross, 1981).The aerodynamic conductance for heat and water vapour, at a given wind speed and atmospheric stability, are influenced by the amount, spatial organisation and drag of the phytoelements (Raupach, 1992;Massman, 1997).The surface conductance, g sV , which lumps together contributions from the soil surface (Mahfouf and Noilhan, 1991) and the canopy, depends on the wetness of the surface soil, the amount of transpiring leaf area and the leaf stomatal conductance.The latter is governed in a complex fashion by the aerial environment surrounding the leaves, as well as via root-to-shoot feedback mechanisms by soil water availability (Larcher, 2001).
In summary, the amount of vegetation is a crucial biotic control on land surface energy partitioning and its influence has been studied for several ecosystems in the past (e.g.Iritz and Lindroth, 1996;Gholtz and Clark, 2002;Beringer et al., 2005;Amiro et al., 2006).In contrast to evergreen canopies which show little seasonal leaf area dynamics or natural deciduous ecosystems whose leaf area development is often closely correlated with environmental conditions, grasslands managed for hay or silage are ideal for studying the effects of changes in the amount of vegetation on energy partitioning.Because of the vegetation cuttings, which are followed by rapid plant re-growth, their canopies undergo multiple growing cycles within a single vegetation period (Wohlfahrt and Cernusca, 2002), which allows to disentangle the biotic and abiotic controls on energy partitioning under a wide range of environmental conditions.
The objective of the present paper is to study the energy partitioning of a temperate mountain grassland in the Stubai Valley (Austria) managed as a hay meadow.To this end we report six years of sensible and latent energy fluxes measured by means of the eddy covariance method and accompanying measurements of net radiation, soil heat flux, above-ground plant matter and meteorological driving forces.We hypothesise that due to rapid changes in the amount of above-ground vegetation, the amount of leaf area will play a major role in determining energy partitioning in this grassland.

Site description
The study site is located at a meadow in the vicinity of the village Neustift (47 • 07 N, 11 • 19 E) in the Stubai Valley (Austria).The site is situated at an elevation of 970 m a.s.l. in the middle of the flat valley bottom.The fetch is homogenous up to 300 m to the east and 900 m to the west of the instrument tower, the dominant day and night time wind directions, respectively.
The average annual temperature is 6.5 • C, average annual precipitation amounts to 852 mm.The vegetation has been classified as a Pastinaco-Arrhenatheretum and consists mainly of a few dominant graminoid (Dactylis glomerata, Festuca pratensis, Phleum pratensis, Trisetum flavescens) and forb (Ranunculus acris, Taraxacum officinale, Trifolium repens, Trifolium pratense, Carum carvi) species.The meadow is cut three times a year; fertilisation (solid manure) usually occurs in late autumn and occasionally after cutting (liquid manure).
The soil has been classified as a Fluvisol (FAO classification) and is approximately 1 m deep.Below a thin (0.001 m) organic layer, an A horizon, with an organic volume fraction of approximately 14%, extends down to 0.02 m, followed by the B horizon, which is best described as a (sandy) loam.Roots reach down to 0.5 m, but 80% of them are concentrated in the upper 0.13 m of the soil.

Eddy covariance (EC)
EC measurements at this site are operational since March 2001 and measurements continue as of this writing.Within this paper, data from March 2001 until December 2006 are presented.Sensible (H ) and latent (LE) heat fluxes were measured using the eddy covariance method (Baldocchi et al., 1988;Baldocchi, 2003) using the same instrumentation as and following the procedures of the EUROFLUX project (Aubinet et al., 2000).The three wind components and the speed of sound were measured by a three-dimensional sonic anemometer (R3IA, Gill Instruments, Lymington, UK).H 2 O mole fractions were measured by a closed-path infra-red gas analyser Lincoln,NE,USA).Air was sucked from the inlet, a distance of 0.1 m from the centre of the sensor volume of the sonic anemometer mounted at 3 m above ground, through a 4 m Teflon tube of 0.004 m inner diameter through a filter (Acro 50, Gelman, Ann Arbor, MI, USA) to the infra-red gas analyser at a flow rate of 9 l min −1 (N035ANE, KNF Neuberger, Freiburg, Germany).The infra-red gas analyser was operated in the absolute mode, flushing the reference cell with dry N 2 from a gas cylinder at 100 ml min −1 .Raw voltage signals of the H 2 O mole fractions were output at 10 Hz to the analogue input of the sonic, where they were synchronised with the sonic signals, which were measured at 20 Hz.All raw data were saved to a hard disc of a personal computer for postprocessing using the EdiSol software (University of Edinburgh).Half-hourly mean eddy fluxes of latent and sensible heat were calculated as the covariance between the turbulent departures from the mean of the vertical wind speed and the H 2 O mixing ratio and the sonic temperature (T s ), respectively, using the post-processing software EdiRe (University of Edinburgh).Means and turbulent departures there from were calculated by Reynolds (block) averaging.The tubeinduced time delay of the H 2 O signal (approx.1.2 s) was determined by optimising the correlation coefficient with the vertical wind velocity (McMillen, 1988) within a given time window.A three-axis co-ordinate rotation was performed aligning the co-ordinate system's vector basis with the mean wind streamlines (Kaimal and Finnigan, 1994).Frequency response corrections were applied to raw eddy fluxes accounting for low-pass (sensor separation, dynamic frequency response, scalar and vector path averaging, frequency response mismatch and the attenuation of concentration fluctuations down the sampling tube) and high-pass filtering following Moore (1986) and Aubinet et al. (2000).Experimentally derived frequency response correction factors, according to Aubinet et al. (2000Aubinet et al. ( , 2001)), were used to calibrate and assess the validity of the theoretical low-pass filtering correction method, as described in detail in Wohlfahrt et al. (2005).Finally, the sensible heat flux was corrected for the effects of air humidity following Schotanus et al. (1983).
Net fluxes of sensible and latent heat were calculated as the sum of the corrected vertical eddy term and the storage fluxes, the latter being estimated from the time-rate-ofchange of the H 2 O concentration and air temperature at the reference height, which in a previous comparison with a profiling system was found to be sufficiently accurate.Negative fluxes represent transport from the atmosphere towards the surface, positive ones the reverse.

Ancillary data
Supporting meteorological measurements of relevance to this study included albedo, net-(R net ) and global radiation (R g ) by separately measuring the up-and downward longand short-wave radiation components (CNR-1, Kipp & Zonen, Delft, The Netherlands), air temperature (T air ) and humidity (RH) at 2 m height and soil temperature (T soil ) at 0.05 m depth, measured by the means of a combined temperature/humidity sensor (RFT-2, UMS, Munich, Germany) and an averaging soil thermocouple probe (TCAV, Campbell Scientific, Logan, UT, USA), respectively, soil heat flux (G) measured by the means of heat flux plates (3 replicates at 0.05 m depth (HFP01, Hukseflux, Delft, The Netherlands)), corrected for the change in heat storage above that depth by using the calorimetric method (Sauer and Horten, 2005), soil water content (SWC) (ML2x, Delta-T Devices, Cambridge, UK) and precipitation (Precip) (52202, R. M. Young, Traverse City, MI, USA).Volumetric soil water content was converted to plant available water (PAW, %) by normalising between the water content at field capacity (100% PAW) and the wilting point (0% PAW), which were determined from water retention curve analysis (Hillel, 1980).
Because green, transpiring stems make up an appreciable fraction of the above-ground biomass (Wohlfahrt et al., 2001), we use the green area index (GAI, m 2 m −2 ), which comprises leaves and green stems, instead of the commonly used leaf area index for quantifying the amount of transpiring surface.The GAI was assessed (i) in a destructive fashion by clipping of square plots of 0.09 m 2 (3-5 replicates) and sub-sequent plant area determination (Li-3100, Li-Cor, Lincoln, NE, USA) and (ii) from measurements of maximum canopy height (h, m) which was related to destructively measured GAI using the following relationship: Continuous time series of the GAI were derived by fitting sigmoid and quadratic functions to measured data separately for each growing phase before and after the third cutting, respectively.

Data analysis
Half-hourly flux data were screened for validity by removal of time periods with (i) the H 2 O and the temperature signal outside a specific range, (ii) the coefficient of variation for H 2 O concentration, temperature and pressure within the IRGA outside a specific range, (iii) the third rotation angle exceeding ±10 • (McMillen, 1988), (iv) the stationarity tests for latent and sensible heat fluxes exceeding 60% (Foken and Wichura, 1996), (v) the deviation of the integral similarity characteristics larger than 60% (Foken and Wichura, 1996) and (vi) the maximum of the footprint function (Hsieh et al., 2000) outside the boundaries of the meadow (Novick et al., 2004).For calculating midday average values of energy fluxes, the diurnal course of these fluxes was quantified using the diurnal centroid method (Wilson et al., 2003).The centroids for the latent and sensible heat flux were found to occur between 12:15 and 13:30 CET and between 10:45 and 12:45 CET, respectively.Based on these findings, the period for calculating the mean midday metrics of β (Bowenratio, H/LE), LE/R net , H/R net , G/R net and α was chosen from 10:00 CET until 14:00 CET.In the case of the albedo (α), this procedure had the additional benefit of minimizing the effect of low solar angles (Li et al., 2006).
Closure of the energy balance, i.e. matching the sum of LE and H to the difference between R net and G, is often considered to be an independent approach to check the overall performance of the EC system (Wilson et al., 2002b).The gap in energy balance closure at the study site in Neustift is on average 21% (range: 32% (2003)-14% ( 2001)) and is thus within the range reported for 22 FLUXNET sites in Wilson et al. (2002b).Presuming that both turbulent fluxes are underestimated by similar relative fractions (Wilson et al., 2002a), energy balance closure was achieved by increasing the turbulent fluxes by the mean imbalance, thereby preserving the Bowen-ratio (Twine et al., 2000).

Climatic conditions
As shown in Fig. 1

Canopy development and albedo
The amount of GAI ranged from close to nil after snow melt, up to 7.6 m 2 m −2 on the verge of cutting (Fig. 2).The meadow was cut three times a year, the cuts taking place between 2 and 16 June, 24 July and 12 August, and 21 September and 28 October.Cutting reduced the GAI down to 0.5-2.0m 2 m −2 (Fig. 2).
During periods of snow cover almost all shortwave radiation was reflected by the white surface and midday means of the albedo (α) varied between 0.8-1.0 (Fig. 2).During the snow-free (vegetation) period α varied with GAI in an asymptotic fashion (α=0.1765GAI 0.1295 , R 2 =0.40, p<0.001), 90% of the maximum vegetation period α (0.23; midday mean) being reached at a GAI of 3.4 m 2 m −2 (Figs. 2  and 3).Effects of changing solar angles over the course of the vegetation period were not evident (data not shown).

Seasonal variation in energy partitioning
During most of the winter time, energy fluxes were minimal -R net and H were slightly negative, while G and LE undulated around zero (Fig. 4).With the initiation of snowmelt, R net quickly increased and the mean monthly diurnal courses ranged from −70 W m −2 up to 520 W m −2 during the vegetation period, closely following global radiation (R g , R net =0.68 R g -40.82,R 2 =0.96, p<0.001).Due to the low canopy height and the absent snow cover, G increased rapidly after snowmelt, consuming up to 25% of R net and Biogeosciences, 5, 421-431, 2008 reaching maximum monthly averaged values of 135 W m −2 in March/April (Fig. 4).During the vegetation period, H ranged from −50 W m −2 up to 170 W m −2 (90% of data between 10% and 44% of R net (midday means)) and LE ranged from −10 W m −2 to 440 W m −2 (90% of data between 35% and 73% of R net (midday means)) as shown in Figs. 4 and 5.
Overall LE was the dominating energy flux over the course of the vegetation period (54% of R net (midday means)) and midday means of β amounted on average to 0.55, with 90% of the values ranging between 0.04 and 1.06 (vegetation period only).Because H increased faster than LE in spring and sustained high values longer than LE in autumn (Figs. 4  and 5), midday means of β exhibited a concave shape with β being one or higher during spring and autumn and close to 0.15 in summer (Fig. 5).Deviations from these general patterns were caused by the cutting events, which led to sustained increases in G (e.g. in 2006; Fig. 5), but also increased the fraction of H at the expense of LE for a few days after cutting (e.g.first cut 2001 and 2002, second cut 2006; Fig. 5).

Biotic controls on energy partitioning
The dependence of energy partitioning on GAI is shown in p<0.001),G decreased (R 2 =0.32, p<0.001), while H decreased initially and then increased again at high GAI values (R 2 =0.10, p<0.001;Fig. 6).The canopy drag coefficient (C d =u 2 * u −2 ), and thus the aerodynamic conductance (g a ), increased within increasing GAI (C d =0.0058e 0.1087GAI , R 2 =0.14, p<0.001).Table 1.Coefficients of determination (R 2 ) of linear regressions of the fractions of latent heat (LE), sensible heat (H ) and soil heat flux (G; dependent variables) against solar radiation (R g ), air temperature (T air ), friction velocity (u * ), plant available water (PAW) and vapour pressure deficit (VPD; independent variables).Values in parentheses show the coefficients of determination after the influence of GAI, as shown in Fig. 6, has been removed.

Abiotic controls on energy partitioning
Among the environmental influence factors considered (Table 1), air temperature (T air ) and vapour pressure deficit (VPD) explained around 40 and 30% of the variability in LE/R net and H/R net , respectively (p<0.001), but were highly correlated with each other (R 2 =0.82, p<0.001).Removing the influence of the GAI on LE/R net (as depicted in Fig. 6) reduced the correlation between T air and LE/R net (p<0.001), while little change was observed for VPD (p<0.001,Table 1).None of the environmental parameters considered explained more than 16% of the variability in G/R net , which reduced to 6% after removing the influence of GAI (Table 1).The ratio of H to LE, β, decreased roughly exponentially with increasing VPD, irrespective of soil water availability (Fig. 7).

Discussion
We have quantified the fluxes of sensible and latent heat, the soil heat flux, net radiation, above-ground phytomass and the major environmental driving forces at a temperate mountain grassland in the Stubai Valley (Austria) over a period of six years.The study site, managed as a hay meadow, is cut three times per year, causing significant temporal dynamics in vegetation cover.Based on these rapid changes in the amount of vegetation, we hypothesised that the partitioning of net radiation into the fluxes of sensible and latent heat and the soil heat flux would be mainly under biotic control (through the amount of aboveground vegetation) and that abiotic (environmental) controls would be of secondary importance for energy partitioning.
The amount of green area (GAI) was equally, and in the case of G even more, important for explaining the observed variability in energy partitioning as the investigated environmental driving forces (Fig. 6, Table 1): The fraction of R net consumed in LE increased with increasing GAI in an asymptotic fashion from initially 37% up to 61%.Increases in LE with increasing leaf areas have been observed for crops (Iritz and Lindroth, 1996) and grasslands ( Rosset et al., 1997;Li et al., 2006) and are generally attributed to increases in surface conductance, which comprises contributions from both the soil surface and the canopy (i.e.under dry conditions mostly stomatal) conductance.The saturation of LE with respect to GAI at high values of GAI has two probable causes: first, the contribution of the soil surface to total ecosystem evapotranspiration decreases with increasing leaf area (Schulze et al., 1996;Li et al., 2006;Saigusa et al., 1998;Wilson and Baldocchi, 2000).According to Rutter (1975) soil evaporation makes up around 10% of total evapotranspiration in reasonably closed canopies, while Körner (1977) reports soil evaporation contributions of up to 30% for low-stature Alpine grassland ecosystems with low canopy cover.Second, the fraction of shaded leaves, which receive only diffuse light, increases with increasing LAI (Ross, 1981).In contrast to sunlit leaves, the photosynthesis of shaded leaves is usually limited by light and not by the partial pressure of CO 2 (Wohlfahrt et al., 1999).CO 2 enters leaves through the same pathway water vapour is lost, the stomata (Larcher, 2001) and in order to avoid excessive water loss, shaded leaves thus generally reduce stomatal conductance and thus transpiration (Mooney et al., 1983;Bunce, 2000;Matsumoto et al., 2005).An alternative explanation of the saturation-type response of LE/R net with GAI would be that the high evapotranspiration rates lead to an increase in vapour pressure and thus a decrease of the evaporative demand, which in turn feeds back on evapotranspiration (e.g.Freedman et al., 2001).However, as shown in Fig. 7, high evapotranspiration rates (i.e.low bowen-ratios) occur during times with high VPD.Apparently, the evaporated water is transported away from the surface by efficient turbulent mixing, preventing low VPDs and their above-mentioned negative feedback on LE.
The fraction of net radiation partitioned into G decreased from 25% at low GAI values to 9% at high GAIs (Fig. 6), a trend which is in accordance with findings from various ecosystems (Song et al., 1997;Huang and Lyons, 1995;Shen et al., 2004;Li et al., 2005).According to Yang et al. (1999) there are several explanations for the reduction of G: the increasing vegetation cover reduces solar energy received by the ground surface, leading to lower amplitudes of the temperature wave in the soil and thereby reducing G. Another cause quoted by Yang et al. (1999) is heat storage in the canopy layer, which we did not account for in our energy balance assessment, which increases with increasing aboveground phytomass.According to Iritz and Lindroth (1996) this may lead to a phase shift of the peak of G towards the afternoon with increasing ground cover, a finding corroborated in the present study (data not shown).Yang et al. (1999) also discuss the role of within-canopy thermal radiation: sunlit leaves at the top of the canopy receive more radiation than shaded leaves close to the ground, which remain cooler and emit less thermal radiation to the soil.In addition, the energy storage term in the soil layer above the heat flux plates also decreases with increasing LAI (Iritz and Lindroth, 1996).Furthermore, the boundary layer resistance for heat increases with decreasing wind speed (Monteith and Unsworth, 1990).Wind speed decreases within the canopy and is thus lower under a vegetated surface than over bare soil.Hence the soil surface conductance for heat and water vapour decrease with increasing GAI, reducing the exchange of sensible and latent heat at the soil surface.Other factors explaining some of the scatter around the correlation between G and GAI are the SWC (Table 1), which is positively correlated with G ( Santanello and Friedl, 2003).
The fraction of R net used for H initially decreased with increasing GAI up to around GAI=4−5 m 2 m −2 ; at higher values the trend reversed (Fig. 6).H/R net decreasing with GAI is somewhat unexpected, as the drag coefficient, which is directly related to the aerodynamic conductance, g a , increased with increasing GAI, a finding supported by several other studies (e.g.Beringer et al., 2005;Hall, 2002).According to Eq. ( 2), an increase in g a should yield an increase in H , provided that the difference between surface and air temperature remains constant.We thus conclude that the decrease in H with increasing GAI reflects a decrease of the surface relative to air temperature caused by the vigorously transpiring canopy and that H represents some sort of remainder term in the energy balance.
The albedo plays an important role for energy partitioning, as it determines the amount of energy available to the ecosystem.The portion of reflected short wave radiation can be influenced by many factors, such as changes in land cover types (Forster et al., 2007), phenology (Song, 1999), canopy structure (Rosset et al., 2001), tissue chemistry, atmospheric conditions, solar and viewing geometry (Asner, 1998;Oguntunde and van de Giesen, 2004).Beside the presence of snow cover, the amount of GAI was found to be of dominating influence for the albedo during the vegetation period.The albedo increased with GAI in a saturation type fashion, resulting in maximum values of 0.23, which is in accordance with other grasslands studies (e.g.Rosset et al., 1997;Li et al., 2006;Song, 1999).The dependency of albedo on GAI can be explained by two factors.(i) The steep increase of the albedo at low GAI values results from the decreasing contribution of the soil to the canopy reflectance (Rosset et al., 2001;Asner, 1998;Oguntunde and Van de Giesen, 2004).According to Goudrian (1977) the effect of the soil on the canopy reflectance is practically negligible at LAI>2 m 2 m −2 .(ii) There is an increase in the canopy reflectance in the near-infra-red (NIR) band with increasing GAI, which diminishes in importance with increasing GAI, whereas reflectance in the visible region is sensitive to LAI only at very low values (LAI<1.0m 2 m −2 ; Asner, 1998), resulting in a saturation-type relationship.Scatter around this trend can be further explained by factors like changes in sky turbidity, leaf angle distribution (Asner, 1998) and phenology (Rosset et al., 2001).The decrease of albedo at low GAI (Fig. 3) contrasts other studies (e.g.Li et al., 2000), where albedo has been found to increase as canopy cover decreases.The likely cause for this difference is soil reflectance, which apparently is relatively low (i.e. the soil appears dark) at our site when the GAI is small, that is during early spring when soils are still water-saturated (Monteith and Unsworth, 1990) after snow melt (Fig. 1).
The dominating abiotic factor, controlling the energy availability, as well as energy partitioning was the presence/absence of a snow cover.Beside low air temperatures during winter, also the high reflectance of the white surface (Fig. 2), led to slightly negative R net values during these times of the year (Fig. 4).Accordingly small were G, H and LE; the latter two, in addition, were inherently noisy due to the combination of a stable stratification of the atmosphere and low wind speeds.
The two abiotic factors most frequently quoted for affecting energy partitioning are soil moisture, and thus the ability of plants to take up water, as well as the vapour pressure deficit (VPD), a measure of the evaporative demand of the atmosphere (Gu et al., 2006;Li et al., 2007).Low soil moisture and/or high VPD lead to stomatal closure (Larcher, 2001) and thus a reduction in the latent heat flux, which most often is compensated by increases in the flux of sensible heat (Valentini et al., 1995;Grant et al., 2006).Despite high VPDs and reductions in soil moisture close to the wilting point (Fig. 1), in particular during 2003 when severe water stress was diagnosed for several forest and agricultural ecosystems in Europe (Granier et al., 2007), no indication of major shifts in energy partitioning was found in the present study: The Bowen-ratio (β, H/LE) decreased with increasing VPD irrespective of soil water availability, which indicates that LE was not restricted by low soil water availability and high evaporative demand (Fig. 7, Table 1).
At present we can offer two explanations for this lack of water stress symptoms despite very low soil water contents: First, the 40% PAW threshold as an indicator of water stress, proposed and tested by Granier et al. (1999Granier et al. ( , 2007) ) for trees, may simply not be applicable to grassland ecosystems, which differ physiologically and structurally considerably from forest ecosystems.Second, while 80% of the roots are confined to the 0-0.12 m soil layer, plants may continue taking up water from deeper in the soil when the upper soil dries out, causing our soil water content measurements at 0.05 m soil depth to underestimate true soil water availability (Miller et al., 2007).Supporting evidence for this hypothesis derives from soil water content measurements at 0.2 m soil depth during 2006, which showed that even when the soil water content at 0.05 m soil depth dropped below 5%, it remained at about 40% at 0.2 m depth.Accordingly, Granier's 40% PAW threshold for the onset of water stress at 0.05 m relates to PAW=55% at 0.2 m.

Fig. 1 .
Fig. 1.Seasonal and inter-annual variability (March 2001-December 2006) of daily sums of global radiation (R g ) and precipitation (Precip; vertical bars) and daily averages of air temperature, vapour pressure deficit (T air and VPD, at 2 m above ground) and plant available water (PAW; closed symbols; at 0.05 m soil depth).

Fig. 2 .
Fig. 2. Seasonal and inter-annual variability (March 2001-December 2006) of green area index (GAI; measured data (open circles) and sigmoid/quadratic fits (lines)) and albedo (midday means).Snow cover periods and cutting events are indicated by grey horizontal and vertical bars respectively.

Fig. 4 .
Fig. 4. Mean monthly diurnal courses of net radiation (R net ), ground heat flux (G), latent (LE) and sensible heat flux (H ).The latter two were increased by the mean imbalance of 21% preserving the Bowen-ratio.Shaded areas refer to ±1 standard deviation.

Fig. 5 .
Fig. 5. Seasonal and inter-annual variability (March 2001-December 2006) of midday means of latent heat flux (LE), sensible heat flux (H ) and ground heat flux (G) expressed relative to net radiation (R net ) and Bowen-ratio (β).Ratios of H/R net and LE/R net were increased in order to meet energy balance closure, preserving the Bowen-ratio.Vertical grey bars indicate cutting events.
Ratios of H/R net and LE/R net were increased in order to meet energy balance closure, preserving the Bowen-ratio.Open symbols refer to midday means, closed symbols to bin averaged data (1 m 2 m −2 bins).Lines show best fits to the raw data.Error bars refer to ±1 standard deviation of bin-averaged data.