Observation and Estimation of Evapotranspiration from an Irrigated Green Roof in a Rain-Scarce Environment

While the rain-driven evapotranspiration (ET) process has been well-studied in the humid climate, the mixed irrigation and rain-driven ET process is less understood for green roof implementations in dry regions, where empirical observations and model parameterizations are lacking. This paper presents an effort of monitoring and simulating the ET process for an irrigated green roof in a rain-scarce environment. Annual ET rates for three weighing lysimeter test units with non-vegetated, sedums, and grass covers were 2.01, 2.52, and 2.69 mm d−1, respectively. Simulations based on the three Penman–Monteith equation-derived models achieved accuracy within the reported range of previous studies. Compared to the humid climate, the overestimation of high ET rates by existing models is expected to cause a larger error in dry environments, where the enhanced ET process caused by repeated irrigations overlapped with hot, dry conditions often occurs during summer. The studied sedum species did not show significantly lower ET rates than native species, and could not effectively take advantage of the deep moisture storage. Therefore, native species, instead of the shallow-rooted species commonly recommended in humid climates, might be a better choice for green roofs in rain-scarce environments.

Green roofs mainly rely on evapotranspiration (ET) to regenerate the retention capacity between rainfall-runoff events [22].The evaluation of other green roof sustainability objectives (like roof cooling and food production) also requires the accurate estimation of ET.In addition to its insulation [14,16,21], Water 2018, 10, 262 4 of 23 green roof ET rates, because of their previous successful applications in GI [24,33,62] and their capacity to account for the advection effect.The water balance of a weighing lysimeter is: where ET a is actual evapotranspiration (mm), P is precipitation (mm), Irrig is irrigation (mm), O is surface overflow (mm), Q is discharge from the underdrain (mm), and ∆S is the change in water storage (mm).The surface overflow was neglected in this study due to the highly pervious soil used and the lack of observed runoff.Infiltration is implicitly tracked by Q and ∆S.During the periods when precipitation, irrigation, surface overflow, and discharge were zero, ET was then determined as the mass change tracked by the weighing lysimeters: ET a = −∆S. ( Water 2018, 10, x FOR PEER REVIEW 4 of 22 where ETa is actual evapotranspiration (mm), P is precipitation (mm), Irrig is irrigation (mm), O is surface overflow (mm), Q is discharge from the underdrain (mm), and ∆S is the change in water storage (mm).The surface overflow was neglected in this study due to the highly pervious soil used and the lack of observed runoff.Infiltration is implicitly tracked by Q and ∆S.During the periods when precipitation, irrigation, surface overflow, and discharge were zero, ET was then determined as the mass change tracked by the weighing lysimeters: The lysimeter shells were rectangular acrylic containers (1.22 m by 0.61 m by 0.36 m) sealed with adhesives and silicone.To match the green roof, an underdrain, made of a bulkhead fitting and a valve, was installed at the bottom of the lysimeter on one side.Each lysimeter contained a replica of layers of the existing green roof at the site.The depth of the growth medium within the container was designed to be 254 mm, 17% less than the depth of the referenced green roof.A deeper lysimeter requires a more costly load cell with a greater loading capacity.Such a load cell, however, would have a lower resolution that affects the accuracy of observed ET rates.Therefore, compared to the referenced green roof, a slightly shallower soil depth was selected to achieve more accurate and economical ET observations.Three custom-built weighing lysimeters of the type described above were established to contain covers including non-vegetated growing medium (as the control), sedums (Red Carpet Stonecrop, Sedum spurium 'Red Carpet' and Russian Stonecrop, Sedum kamtschaticum), and grass (Blue Grama grass, Bouteloua gracilis) (Figure 2).The plants were transplanted into the lysimeters in September 2013, after several years of growth on site as part of the green roof installation.Rice Lake BenchMark HE weighing scales (Rice Lake, WI, capacity: 454 kg, resolution: 0.05 kg ~ 0.06 mm, tolerance: 0.08 kg ~ 0.10 mm) were adopted as the load cells managed by IQ355 indicators (Rice Lake, WI, USA).Buried at the surface (~10 mm deep) and bottom (~360 mm deep) of each lysimeter after lab calibrations, Decagon (Pullman, WA, USA) 5TM sensors were used to measure soil moisture (resolution: 0.08% VWC).Attached to and protected by a PVC pipe on the outside, a tipping bucket rain gauge (TE525, resolution: 0.25 mm, Campbell Scientific, Logan, UT, USA) was half buried underground to monitor The lysimeter shells were rectangular acrylic containers (1.22 m by 0.61 m by 0.36 m) sealed with adhesives and silicone.To match the green roof, an underdrain, made of a bulkhead fitting and a valve, was installed at the bottom of the lysimeter on one side.Each lysimeter contained a replica of layers of the existing green roof at the site.The depth of the growth medium within the container was designed to be 254 mm, 17% less than the depth of the referenced green roof.A deeper lysimeter requires a more costly load cell with a greater loading capacity.Such a load cell, however, would have a lower resolution that affects the accuracy of observed ET rates.Therefore, compared to the referenced green roof, a slightly shallower soil depth was selected to achieve more accurate and economical ET observations.Three custom-built weighing lysimeters of the type described above were established to contain covers including non-vegetated growing medium (as the control), sedums (Red Carpet Stonecrop, Sedum spurium 'Red Carpet' and Russian Stonecrop, Sedum kamtschaticum), and grass (Blue Grama grass, Bouteloua gracilis) (Figure 2).The plants were transplanted into the lysimeters in September 2013, after several years of growth on site as part of the green roof installation.Rice Lake BenchMark HE weighing scales (Rice Lake, WI, capacity: 454 kg, resolution: 0.05 kg ~0.06 mm, tolerance: 0.08 kg ~0.10 mm) Water 2018, 10, 262 5 of 23 were adopted as the load cells managed by IQ355 indicators (Rice Lake, WI, USA).Buried at the surface (~10 mm deep) and bottom (~360 mm deep) of each lysimeter after lab calibrations, Decagon (Pullman, WA, USA) 5TM sensors were used to measure soil moisture (resolution: 0.08% VWC).Attached to and protected by a PVC pipe on the outside, a tipping bucket rain gauge (TE525, resolution: 0.25 mm, Campbell Scientific, Logan, UT, USA) was half buried underground to monitor the timing of the underdrain discharge for each of the three lysimeters (Figure 3).The observation during the drainage periods could then be tracked and excluded from the analysis.Additional TE525 tipping bucket rain gauges were placed next to each of the lysimeters on the roof in order to measure the precipitation and irrigation.All data were recorded as 5-min averages by a Campbell Scientific CR1000 datalogger (Campbell Scientific, Logan, UT, USA) linked with a Campbell Scientific AM16/32B multiplexer (Campbell Scientific, Logan, UT, USA).The measured 5-min ET a rates were summed to daily values to be used in the following analyses and comparison.the timing of the underdrain discharge for each of the three lysimeters (Figure 3).The observation during the drainage periods could then be tracked and excluded from the analysis.Additional TE525 tipping bucket rain gauges were placed next to each of the lysimeters on the roof in order to measure the precipitation and irrigation.All data were recorded as 5-min averages by a Campbell Scientific CR1000 datalogger (Campbell Scientific, Logan, UT, USA) linked with a Campbell Scientific AM16/32B multiplexer (Campbell Scientific, Logan, UT, USA).The measured 5-min ETa rates were summed to daily values to be used in the following analyses and comparison.The micro-meteorological conditions were monitored by a weather station (Met.Tower in Figure 2).A Kipp & Zonen CNR-1 (Delft, The Netherlands) net radiometer was mounted 3 m above the roof level, which measured both incoming and outgoing components of the short and longwave radiation.The 2-m wind speed was monitored by a Campbell Scientific CSAT3 sonic anemometer (20 Hz) (Campbell Scientific, Logan, UT, USA).Four thermocouple probes (Omega HTTC36-T-18G-6, Stamford, CT, USA) were mounted at 0.12 m, 0.60 m, 1.50 m, and 2.50 m above the roof level to determine the vertical air temperature profile.Relative humidity was measured at 2 m using a Campbell Scientific HMP155A probe (Campbell Scientific, Logan, UT, USA).Meteorological data were recorded by a Campbell Scientific CR5000 datalogger (Campbell Scientific, Logan, UT, USA).The whole experiment was powered by directly connected solar panels.Standard operational the timing of the underdrain discharge for each of the three lysimeters (Figure 3).The observation during the drainage periods could then be tracked and excluded from the analysis.Additional TE525 tipping bucket rain gauges were placed next to each of the lysimeters on the roof in order to measure the precipitation and irrigation.All data were recorded as 5-min averages by a Campbell Scientific CR1000 datalogger (Campbell Scientific, Logan, UT, USA) linked with a Campbell Scientific AM16/32B multiplexer (Campbell Scientific, Logan, UT, USA).The measured 5-min ETa rates were summed to daily values to be used in the following analyses and comparison.The micro-meteorological conditions were monitored by a weather station (Met.Tower in Figure 2).A Kipp & Zonen CNR-1 (Delft, The Netherlands) net radiometer was mounted 3 m above the roof level, which measured both incoming and outgoing components of the short and longwave radiation.The 2-m wind speed was monitored by a Campbell Scientific CSAT3 sonic anemometer (20 Hz) (Campbell Scientific, Logan, UT, USA).Four thermocouple probes (Omega HTTC36-T-18G-6, Stamford, CT, USA) were mounted at 0.12 m, 0.60 m, 1.50 m, and 2.50 m above the roof level to determine the vertical air temperature profile.Relative humidity was measured at 2 m using a Campbell Scientific HMP155A probe (Campbell Scientific, Logan, UT, USA).Meteorological data were recorded by a Campbell Scientific CR5000 datalogger (Campbell Scientific, Logan, UT, USA).The whole experiment was powered by directly connected solar panels.Standard operational The micro-meteorological conditions were monitored by a weather station (Met.Tower in Figure 2).A Kipp & Zonen CNR-1 (Delft, The Netherlands) net radiometer was mounted 3 m above the roof level, which measured both incoming and outgoing components of the short and longwave radiation.The 2-m wind speed was monitored by a Campbell Scientific CSAT3 sonic anemometer  (20 Hz) (Campbell Scientific, Logan, UT, USA).Four thermocouple probes (Omega HTTC36-T-18G-6, Stamford, CT, USA) were mounted at 0.12 m, 0.60 m, 1.50 m, and 2.50 m above the roof level to determine the vertical air temperature profile.Relative humidity was measured at 2 m using a Campbell Scientific HMP155A probe (Campbell Scientific, Logan, UT, USA).Meteorological data were recorded by a Campbell Scientific CR5000 datalogger (Campbell Scientific, Logan, UT, USA).The whole experiment was powered by directly connected solar panels.Standard operational weather stations (ID: WBB and MTMET) from the nearby campus operated by Mesowest (http://mesowest.utah.edu/)were used to fill the gaps caused by power failures.A Decagon SC-1 Leaf Porometer (Pullman, WA, USA) was employed to measure the diurnal stomatal resistances of sedums and grass in the lysimeters on 17 June 2015, in order to corroborate ET estimates.The leaf area was estimated by branch sampling and image analysis via Easy Leaf Area, which uses the color ratio of image pixels to distinguish leaves from the background [63].

Evapotranspiration Estimation
Based on the Monin-Obukhov Similarity Theory and the surface energy balance, the Penman equation and its family, including the Penman-Monteith (P-M) equation, were derived for ET estimation.The P-M equation serves as the keystone of other simplified equations like the Priestley-Taylor (P-T) equation [64] and the Advection-Aridity (A-A) equation [65], or revised forms, including the Food and Agriculture Organization (FAO) of the United Nations' irrigation and drainage method (FAO56) [48] and the American Society of Civil Engineers standardized reference ET scheme [66].As the P-M equation and its descendants are commonly used for green roof ET estimation in the humid environment [32][33][34][35]42,43], the P-M equation and its two revisions are tested for the two vegetated covers.Non-vegetated cover, without the transpiration part (no crop coefficient and water stress coefficient), was excluded from the ET estimation.
The FAO56 ET scheme has been widely used as the baseline to estimate green roof ET [32,34,43].It can be calculated by the following workflow.The reference ET (ET o ) was first computed as follows [48]: where ET o is the reference ET corresponding to well-watered short grass (mm h −1 ), ∆ represents the slope of the saturation-vapor-pressure temperature curve (kPa • C −1 ), R n is the net radiation (MJ m −2 h −1 ), G is the soil heat flux (MJ m −2 h −1 ), γ is the psychrometric constant (kPa • C −1 ), T is the air temperature at 2 m ( • C), e s is saturation vapor pressure (kPa), e a is actual vapor pressure (kPa), and u 2 is the mean hourly wind speed at the 2-m height (m s −1 ).Daily summed reference ET can be further revised by a crop coefficient and a water stress coefficient to compute ET a (mm d −1 ) [48]: where K s is the water stress coefficient, K c is the crop coefficient (parameterized as monthly constants in this study), and the product of K c and ET o results in PET (mm d −1 ).Converted to the volume demand with the green roof area, PET can be directly used to estimate the irrigation cost ($): where TP j is the tiered water price ($/m 3 ) corresponding to the jth tier, and V j is the net water demand (m 3 ) of the jth tier during the growth season (April-October), calculated as the residue of subtracting the product of PET and area by the rainfall amounts.The water stress coefficient can be estimated as follows [48]: Water 2018, 10, 262 where θ is the moisture content (fraction) as the average of the measured surface and the bottom soil moisture contents of the lysimeter, θ wp is the wilting point (fraction), θ fc is the field capacity (fraction), and p is the average fraction of the available moisture that can be depleted from the root zone before moisture stress occurs.The field capacity (0.35) was determined via lab tests by oven-drying soil samples.There is no established way to measure the wilting point directly other than visual observation, which, however, is difficult to determine and could be inconsistent between observers and cases [29,33].As no visual wilting was clearly observed, the wilting point of this study was approximated based on the following empirical equations [67,68]: where θ wp is the wilting point (fraction), θ 1500t is the moisture content at a tension of 1500 kPa, S is the weight fraction of sand, C is the weight fraction of clay, and OM is the percent organic matter.
The fraction p can be estimated as follows [48]: where p FAO56 is the recommended value from Table 22 in the FAO56 report [48].The average of the ET a measured from the lysimeter observation was used.The p FAO56 value (0.4) corresponding to bluegrass was used to represent the studied grass.No p FAO56 value was published for sedum species.The fraction p represents the plant's capability of extracting moisture from the root zone.The documented p FAO56 value (0.5) of cucumber was assumed as that of the sedum, since both are succulent species.Once K s is determined, the daily K c can be back-calculated given the computed ET o and the measured ET a .The developed daily K c time series were summarized into the monthly averages.The Thornthwaite-Mather equation (T-M), as the second model tested here, has been often used to convert ET o to ET a for green roofs by only adjusting soil moisture variations [33,38]: The ET o Equation (3) implicitly fixes the surface resistance at 70 s m −1 , or at 50 s m −1 for daytime and 200 s m −1 for nighttime.Only representing the short, cool-season, well-watered reference grass, though, such values are presumed for green roofs as a routine [32,33,42,43].As the cornerstone of all P-M derived methods, the original P-M equation was also examined to estimate ET a after calibrating the surface resistance parameter [69]: where ρ a is the mean air density (kg m −3 ), c p is the specific heat of the air (MJ kg −1 • C −1 ), r a is the aerodynamic resistance (s m −1 ) which can be calculated based on FAO56 scheme [48], λ is the latent heat of vaporization, which can be calculated as λ = 2.501 − 0.00237 × T [70], but 2.45 MJ kg −1 was used in this study, γ is the psychrometric constant (kPa • C −1 ), and r s is surface resistance (s m −1 ).
Given the tracked ET a and the measured meteorological conditions, the hourly surface resistances could be back-calculated via Equation ( 11) [71,72].
As a short summary on the parameterizations, back-calculating the FAO-56 method generates the daily time series of K c , which were then reduced to monthly K c .Back-calculating for the P-M equation produces the hourly time series of r s , which were then reduced to the hourly r s averaged Water 2018, 10, 262 8 of 23 over the month and over the year.The former gives the hourly surface resistances averaged over the month, and the latter gives the hourly surface resistances averaged over the year.The FAO56 method then took the reduced monthly crop coefficients to estimate the ET a , the T-M method did not need to calibrate any parameters, and the P-M method used the reduced hourly surface resistances averaged over the month (Table 1).
Table 1.Parameters used for three evapotranspiration estimation methods.

Models
Note: r s refers to the surface resistance; K c refers to the crop coefficient; K s refers to the water stress coefficient; FAO56 refers to the Food and Agriculture Organization of the United Nations' irrigation and drainage method; T-M refers to the Thornthwaite-Mather equation; P-M refers to the Penman-Monteith equation.

Statistics
The Kruskal-Wallis Analysis of Variance (ANOVA) test and the Mood's Median test were used to determine the difference in the ET a observations between covers.Partial least squares (PLS) analysis was made to examine the influence of selected environmental variables (air temperature, relative humidity, wind speed, air pressure, incoming solar radiation, the sum of precipitation and irrigation, surface soil moisture, and bottom soil moisture) on the ET process.To conduct the PLS analyses on a daily basis, the solar radiation and the sum of precipitation and irrigation were aggregated to daily values, while a daily average was used to represent all other measured variables.For all three ET a models, their 5-min ET a estimates were summed to daily amounts to be compared with the ET a daily measurements.The variable importance in projection (VIP) scores summarize the influence of individual variables on the PLS model, and give a measure useful to select the independent variables that contribute the most to the dependent variable's variance [73].It is generally accepted in practice that variables having a VIP > 1.0 are highly influential, values between 0.8 < VIP < 1.0 indicate moderately influential variables, and variables with VIP < 0.8 are less important [73,74].The coefficients of determination (R 2 ), the root mean square error (RMSE), the Nash-Sutcliffe efficiency (NSE) [75], and the percent bias (PBIAS) were calculated to evaluate the goodness of fit of the three ET a models.Model simulation can be judged satisfactory if NSE > 0.65 or PBIAS < ±15% [43].The statistical analysis was conducted via OriginPro 2015.
To quantify the contributions of the uncertainty of the parameterization process to ET a estimates, the generalized likelihood uncertainty estimation (GLUE) [76] was used in Python 2.7.13 to conduct the uncertainty analysis for the FAO56 model.As only the FAO56 model includes the full suite of parameters of interest (K c , K s , and r s for this study), the TM and PM models were not included in the uncertainty analysis.A uniform distribution within a feasible range (Table 2) was assigned to the parameters, as no evidence existed for other statistical distributions [77,78].The maximum of K c was set according to the previous observation [79], and the threshold of r s was the maximum value observed in the field.The range of the water stress coefficient was based on the definition of the FAO56 method [48].A Monte Carlo simulation with 100,000 runs was performed to create a random parameter space for sedum and grass ET a estimates, respectively.Each run was driven by the available hourly meteorological data of 2014 and independently randomized parameters.As a formal likelihood measure could be more problematic due to the subjectivity of choosing the shape factor [80], a less formal likelihood measure (NSE) was used in the form: Water 2018, 10, 262 where L(θ i |Y) is the likelihood measure for the model with the ith parameter set conditioned on the observations Y, σ 2 i is the error variance for the ith model, σ 2 obs is the observed variance, and N is the shape factor.As higher N values have the effect of accentuating the weight given to the better simulations and narrowing the confidence limits, N = 1 was used [81].The simulated and measured daily ET a time series were compared to compute a yearly NSE as the indicator of the model efficiency.Runs with a negative NSE were excluded.Followed by Staudt et al. [77] and Prihodko et al. [78], the top 10% behavioral threshold for NSE was used to select runs for further analysis.The posterior likelihood can be updated by the application of Bayes equation [81]: where L(Y|θ i ) is the posterior likelihood for the simulation of Y given the ith parameter set and the scaling constant C, L o (θ i ) is the prior likelihood for the parameter set θ i .The scaling constant is used to scale the cumulative of the posterior likelihood to equal to unity.Based on the assumed uniform distribution, the prior likelihood for each parameter set was equal.Then the scaling constant could be calculated as [82]: Predictive uncertainty bounds were determined as the 5% and 95% quantiles of the derived cumulative probability distribution of the ranked model efficiency [83].The cumulative likelihood of each parameter was also plotted against the uniform distribution.The deviation from the uniform distribution indicates the parameter sensitivity, and the steepest slope in the cumulative likelihood curves reflects the optimal range of the parameter values [77].The two-sample Kolmogorov-Smirnov (K-S) test [84] was conducted to evaluate the difference between the generated distribution and the uniform distribution.The K-S value was used to rank the sensitivity of parameters [77,78].

Evapotranspiration Regime
The average observed ET rates are 2.01 ± 1.16 (standard deviation), 2.52 ± 1.79, and 2.69 ± 1.69 mm d −1 over the one-year study period for the non-vegetated, sedum, and grass covers, respectively (Figure 4).The p values of the Kruskal-Wallis ANOVA test and Mood's Median test for the daily ET time series are 2.08 × 10 −4 and 1.53 × 10 −6 , respectively, at a significance level of 0.05.Therefore, the observed ET rates are statistically different among the three surface covers.Compared to the non-vegetated roofs, vegetation could elevate ET rates by over 25% on average.
The average observed ET rates are 2.01 ± 1.16 (standard deviation), 2.52 ± 1.79, and 2.69 ± 1.69 mm d −1 over the one-year study period for the non-vegetated, sedum, and grass covers, respectively (Figure 4).The p values of the Kruskal-Wallis ANOVA test and Mood's Median test for the daily ET time series are 2.08 × 10 −4 and 1.53 × 10 −6 , respectively, at a significance level of 0.05.Therefore, the observed ET rates are statistically different among the three surface covers.Compared to the nonvegetated roofs, vegetation could elevate ET rates by over 25% on average.The observed ET rates of sedums and grass tend to be slightly higher than most studies conducted in rain-rich environments [23,29,33,[85][86][87].For example, Rezaei [85] observed ET rates of 1.68 mm/d for sedums and 1.06 mm/d for the non-vegetated medium with unlimited water supply in an indoor environment.Marasco et al. [86] measured via a dynamic chamber 1.71 mm/d on average for sedums in New York City, NY with 1078 mm annual rainfall, while DiGiovanni et al. [33] measured via lysimeters 2.27 mm/d for sedums in the Bronx, NY with an average annual rainfall of 1100 mm.In Auckland, New Zealand, Voyde et al. [22] reported 2.21 mm/d for sedums with an annual rainfall of 1240 mm via lysimeter measurements.With a lower water supply (a total of 938 mm water: 496 mm precipitation + 442 mm irrigation) in the studied sub-humid climate compared to the above-listed wetter environments, the values of the observed annual ET rates are not correspondingly smaller.In other words, the rain-scarce region does not necessarily lead to a lower yearly ET rate of the green roof compared to the rain-rich region.
Intuitively, the larger yearly ET rates seem to be caused by the added irrigation only.The maximums of the observed monthly ET averages of the sedum (4.07 mm d −1 ) and grass (4.93 mm d −1 ) are higher than the maximum daily rate of 3.4 mm/d from the grass and sedum covers after a 7-day dry period measured in a climate-controlled chamber [29].Yet, the observed ET rates during the summer with irrigation applied are not always significantly larger than the reported values in the rain-rich environments.The maximum monthly ET rates of the sedum and grass are lower than the maximum monthly average (7 mm d −1 ) reported by a multi-year study of sedums in Villanova, PA, USA with an annual rainfall of 817-1352 mm [42], but are close to the reported 4.94 mm d −1 in the rain-rich region [86].There is no consistent relationship between the observed summer ET rates and the reported values.But overall, comparable summer ET rates were found between the studied Water 2018, 10, 262 11 of 23 rain-scarce region and other rain-rich regions.Added irrigation appears to be able to support the energy-limited ET regime in the studied environment.
Four variables were found as the important factors influencing the ET process of the vegetated covers via the PLS analysis (Figure 5).For the sedum cover, the VIP scores for air temperature, surface moisture, relative humidity, and solar radiation are 1.79, 1.26, 1.22, and 0.95, respectively.For the grass cover, the VIP scores for air temperature, solar radiation, surface moisture, and relative humidity are 1.48, 1.33, 1.07, and 1.01, respectively.dry period measured in a climate-controlled chamber [29].Yet, the observed ET rates during the summer with irrigation applied are not always significantly larger than the reported values in the rain-rich environments.The maximum monthly ET rates of the sedum and grass are lower than the maximum monthly average (7 mm d −1 ) reported by a multi-year study of sedums in Villanova, PA, USA with an annual rainfall of 817-1352 mm [42], but are close to the reported 4.94 mm d −1 in the rain-rich region [86].There is no consistent relationship between the observed summer ET rates and the reported values.But overall, comparable summer ET rates were found between the studied rainscarce region and other rain-rich regions.Added irrigation appears to be able to support the energylimited ET regime in the studied environment.
Four variables were found as the important factors influencing the ET process of the vegetated covers via the PLS analysis (Figure 5).For the sedum cover, the VIP scores for air temperature, surface moisture, relative humidity, and solar radiation are 1.79, 1.26, 1.22, and 0.95, respectively.For the grass cover, the VIP scores for air temperature, solar radiation, surface moisture, and relative humidity are 1.48, 1.33, 1.07, and 1.01, respectively.The VIP scores mean that the climate factors play a large role in determining the variations of observed ET rates.Air temperature is discovered to be the most important factor determining ET rates for sedum and grass covers, as it is closely related to other meteorological variables (like vapor pressure deficit) and plant physiological responses (like stomatal conductance).The wind speed has the largest influence on evaporation rates for non-vegetated cover.This is likely a result of the different turbulent flow characteristics associated with bare soil and vegetative canopies.The VIP scores mean that the climate factors play a large role in determining the variations of observed ET rates.Air temperature is discovered to be the most important factor determining ET rates for sedum and grass covers, as it is closely related to other meteorological variables (like vapor pressure deficit) and plant physiological responses (like stomatal conductance).The wind speed has the largest influence on evaporation rates for non-vegetated cover.This is likely a result of the different turbulent flow characteristics associated with bare soil and vegetative canopies.
The important role of climate factors indicates that the ET process of the studied species in the rain-scarce region prevalently had an energy-limited regime due to a combined impact.The climate and irrigation worked together to promote the high ET rates.Applied irrigation could promote an energy-limited regime; within the regime, the ET process would be independent of the irrigation but taken over by climatic factors.With intensified atmospheric demands and prevailing heat advection, the hot, dry summers in the dry regions would have a larger climatic driving force than the summers in the wet regions, for the energy-limited ET regime given the same moisture supply.However, maintaining an energy-limited ET regime in a rain-scarce environment can be costly and may not make sense in terms of water conservation.Native species with stronger drought tolerance, therefore, are recommended as green roof species in a dry environment.

Design Considerations for Water Conservation
Via PLS analysis, solar radiation is another important contributor to ET rates for both studied plant species (Figure 5).Strategies like shading to reduce the duration and magnitude of sunshine might be helpful for suppressing PET rates.Selecting a rooftop that can be partially shaded by (higher part of or nearby) buildings becomes an important design factor for green roof investments in rain-scarce environments.A heterogeneous landscape design with higher native shrubs being able to shade less drought-resistant grass plants on the nearby roof surface is also recommended.In addition to climatic factors, the moisture storage of green roofs also has a larger influence on ET variations than wetting events (Figure 5).The storage created by the medium pores and drainage mats could delay the response of ET rates corresponding to the wetting events.Increasing the size of this moisture storage (by using, for example, a deeper medium or thicker drainage mats) could help green roofs collect more moisture from the wetting events for use during the drought periods.Therefore, shading and increasing the moisture storage would be recommended for green roof applications in dry regions.
Actual evapotranspiration rates of sedum and grass statistically differ from each other.From an analysis of the vertical moisture profiles, this ET difference is related to the diverse rooting depths between the two species.The surface moisture content of the sedum cover was consistently lower than that of the grass cover during April-July (Figure 6).During this period, however, sedums tended to have lower ET rates than the grass (Figure 4).Seemingly, during the early growing stage, sedums, as succulent species, tended to store moisture rather than releasing it.Since August, when the surface moisture of the sedum cover became close to or even larger than that of the grass, the ET rates of sedums became larger than those of the grass.No consistent relation, however, was found between the moisture content of the bottom layer and the observed ET rates.Even when the deep moisture of the sedum-covered lysimeter, for instance, was highest among all three covers in July, the sedum ET amount was still much lower than the grass ET.These indicate that the ET rates and the survival of the studied sedums are more restricted by the surface moisture.Therefore, though the extra soil depth in an intensive green roof could provide more moisture and heat storage for runoff reduction and roof cooling, sedums may not be able to fully take advantage of the moisture stored in the deep soil.From the perspective of water conservation, sedums may not be the ideal species for green roofs in the rain-scarce regions.However, as a succulent species, sedums need a moist surface soil layer to grow healthily.This has to be satisfied by repeated irrigation when precipitation is limited.Dying sedums were observed from another green roof on the University of Utah campus due to the lack of enough irrigation (according to personal correspondence with the campus irrigation technicians).This also partly explains why native species established more readily than exotic succulent species like sedums under extreme conditions [25], where the capacity of tapping extra deeper moisture becomes necessary.

Parameterization for Evapotranspiration Modeling
Calculated yearly averages of crop coefficients are 0.59 ± 0.29 (standard deviation) for sedum and 0.61 ± 0.25 for grass.The sedum Kc found by this study is generally close to and slightly larger than most previous studies conducted in moisture-rich environments.For example, 0.53 was reported for a well-watered sedum canopy measured by a weighing lysimeter of the same base size as this study in Durham, NH, USA [87].An annual average of 0.51 was reported for different sedums in College Park, MD, USA with a 1004 mm annual rainfall, which was estimated based on moisture contents of green roof platforms (1.31 m 2 ) larger than the area of this study [88], while 0.32-0.45was reported as the annual averages for three sedum types also in College Park, MD, USA [33].Similarly, 0.35-0.52 was observed for well-watered sedums calculated based on an energy balance approach for a 1000 m 2 green roof in Vicenza, Italy with an approximately 1087 mm annual rainfall [89].The sedum Kc of this study case is lower than some findings with smaller lysimeters but more moisture supplied.For example, Voyde [90] reported 0.85 for the well-watered sedum measured using 0.072 m 2

Parameterization for Evapotranspiration Modeling
Calculated yearly averages of crop coefficients are 0.59 ± 0.29 (standard deviation) for sedum and 0.61 ± 0.25 for grass.The sedum K c found by this study is generally close to and slightly larger than most previous studies conducted in moisture-rich environments.For example, 0.53 was reported for a well-watered sedum canopy measured by a weighing lysimeter of the same base size as this study in Durham, NH, USA [87].An annual average of 0.51 was reported for different sedums in College Park, MD, USA with a 1004 mm annual rainfall, which was estimated based on moisture contents of green roof platforms (1.31 m 2 ) larger than the area of this study [88], while 0.32-0.45was reported as the annual averages for three sedum types also in College Park, MD, USA [33].Similarly, 0.35-0.52 was observed for well-watered sedums calculated based on an energy balance approach for a 1000 m 2 green roof in Vicenza, Italy with an approximately 1087 mm annual rainfall [89].The sedum K c of this study case is lower than some findings with smaller lysimeters but more moisture supplied.For example, Voyde [90] reported 0.85 for the well-watered sedum measured using 0.072 m 2 weighing trays in Auckland, New Zealand; Schneider [79] observed 1.0-1.7 for sedums measured with a 0.21 m 2 weighing lysimeter in Villanova, PA; and Rezaei [84] observed 1.35 on average for sedums estimated from 0.56 m 2 indoor greenhouse trays.
In spite of the absolute difference, the crop coefficients of sedums estimated by this study are still within the range of the above-listed studies.As the crop coefficients adjust ET o mainly for the characteristics of species and climates, this reinforces that the applied irrigation effectively made up for the lack of rainfall in the studied environment.Notably, the studied sedum species have higher K c than grass for some months (Figure 7).So, if well-watered, the studied sedum species could consume more water than grass, as also observed by Wolf and Lundholm [91].Noticeably, both sedums and grass were observed to have higher K c than the native species [92], indicating the potential application of the native species in green roofs in dry environments given enough soil space for the longer and wider root systems.
The estimated surface resistances appear to vary widely among different plant covers, different months, and different hours (Figure 8).As no value larger than 2500 s m −1 was ever observed for both species during the diurnal field camp, any simulated surface resistances above this threshold were regarded as indicating the closure of stomata (i.e., negligible ET).Not representing the average resistance in the P-M equation, those infinite values were excluded from this study.The annual surface resistances of sedums and grass were estimated as 678 s m −1 and 651 s m −1 , respectively.The corresponding stomatal resistances of sedums and grass during the middle growing stage (as leaf areas keep varying) were estimated as 958 and 606 s m −1 , respectively.The estimated stomatal resistance of sedums is higher than previously reported values (250 s m −1 [90] and 750 s m −1 [79]).A diurnal measurement of stomatal resistance was further made on 17 June 2015 for the sedum and grass in the same lysimeters.Their values were converted to surface resistances after leaf areas were measured, and then compared with the hourly averages of estimated values of June 2014.In spite of one year's difference, the estimated surface resistance still showed a relatively good match to the observed values (Figure 8), validating this parameterization procedure.Notably, the estimated nighttime non-zero resistances in this study seem to be lower than the measured nighttime non-zero resistances.This discrepancy might be because the former calculated the average resistance for the entire plant, while the latter only represented the average resistance of a few sampled leaves.
measured, and then compared with the hourly averages of estimated values of June 2014.In spite of one year's difference, the estimated surface resistance still showed a relatively good match to the observed values (Figure 8), validating this parameterization procedure.Notably, the estimated nighttime non-zero resistances in this study seem to be lower than the measured nighttime non-zero resistances.This discrepancy might be because the former calculated the average resistance for the entire plant, while the latter only represented the average resistance of a few sampled leaves.measured, and then compared with the hourly averages of estimated values of June 2014.In spite of one year's difference, the estimated surface resistance still showed a relatively good match to the observed values (Figure 8), validating this parameterization procedure.Notably, the estimated nighttime non-zero resistances in this study seem to be lower than the measured nighttime non-zero resistances.This discrepancy might be because the former calculated the average resistance for the entire plant, while the latter only represented the average resistance of a few sampled leaves.It has been recognized that sedums could have a crassulacean acid metabolism (CAM) mechanism, which shuts down stomata during the daytime to save carbon dioxide and water, and opens stomata to resume metabolism during the cooler and wetter nights [22,71,93,94].Based on the estimated and measured surface resistances, however, no strong CAM was observed for the sedums in this study, as found by others as well [88,90].The night irrigation might prevent sedums switching from the C 3 mode, with which plants usually open stomata during the daytime, to CAM mode [95][96][97].

Comparison of Model Estimates
The selected three methods (FAO56, T-M, and P-M) were tested to simulate daily ET a values, which were then compared to the observed ET a time series (Figure 9).In spite of a slightly overall overestimation, all methods are close to the 1:1 line.The T-M method had the lowest R 2 for the grass (Table 3).This may be because the T-M method is highly sensitive to the input parameters like the wilting point and field capacity [33], which are often difficult to accurately quantify.The FAO56 method achieved the best fit among the three methods for both vegetated covers.The accuracy of the P-M method was similar to the FAO56 method for grass cover.With less tedious computational work, the P-M method might be used to replace the FAO56 method for estimating ET rates of the grass-covered green roof.Note: R 2 is the coefficient of determination, RMSE is the root mean square error (mm d −1 ), NSE is the Nash-Sutcliffe efficiency, and PBIAS is the percent bias; FAO56 refers to the Food and Agriculture Organization of the United Nations' irrigation and drainage method, T-M refers to the Thornthwaite-Mather equation, and P-M refers to the Penman-Monteith equation.
Water 2018, 10, x FOR PEER REVIEW 14 of 22 It has been recognized that sedums could have a crassulacean acid metabolism (CAM) mechanism, which shuts down stomata during the daytime to save carbon dioxide and water, and opens stomata to resume metabolism during the cooler and wetter nights [22,71,93,94].Based on the estimated and measured surface resistances, however, no strong CAM was observed for the sedums in this study, as found by others as well [88,90].The night irrigation might prevent sedums switching from the C3 mode, with which plants usually open stomata during the daytime, to CAM mode [95][96][97].

Comparison of Model Estimates
The selected three methods (FAO56, T-M, and P-M) were tested to simulate daily ETa values, which were then compared to the observed ETa time series (Figure 9).In spite of a slightly overall overestimation, all methods are close to the 1:1 line.The T-M method had the lowest R 2 for the grass (Table 3).This may be because the T-M method is highly sensitive to the input parameters like the wilting point and field capacity [33], which are often difficult to accurately quantify.The FAO56 method achieved the best fit among the three methods for both vegetated covers.The accuracy of the P-M method was similar to the FAO56 method for grass cover.With less tedious computational work, the P-M method might be used to replace the FAO56 method for estimating ET rates of the grass-covered green roof.Note: R 2 is the coefficient of determination, RMSE is the root mean square error (mm d −1 ), NSE is the Nash-Sutcliffe efficiency, and PBIAS is the percent bias; FAO56 refers to the Food and Agriculture Organization of the United Nations' irrigation and drainage method, T-M refers to the Thornthwaite-Mather equation, and P-M refers to the Penman-Monteith equation.Although the model performance statistics are not particularly good, they fall within the range of most previous studies.The corresponding R 2 for the sedum cover is slightly lower than the best estimate (R 2 = 0.51) of Starry et al. [34] who used the FAO56 scheme to estimate sedum ETa in a rainrich climate, but falls within their reported range (0.07 < R 2 < 0.51).The R 2 of sedum ETa simulation is also lower than the reported R 2 value of DiGiovanni et al. [33] using the T-M method for sedums in a rain-rich climate, but the linear fitting equation of that study fixed the intercept at zero which could Although the model performance statistics are not particularly good, they fall within the range of most previous studies.The corresponding R 2 for the sedum cover is slightly lower than the best estimate (R 2 = 0.51) of Starry et al. [34] who used the FAO56 scheme to estimate sedum ET a in a rain-rich climate, but falls within their reported range (0.07 < R 2 < 0.51).The R 2 of sedum ET a simulation is also lower than the reported R 2 value of DiGiovanni et al. [33] using the T-M method for sedums in a rain-rich climate, but the linear fitting equation of that study fixed the intercept at zero which could lead to a higher R 2 value.The RMSEs of simulated sedum ET a are close to Marasco et al. [35] (1.19 to 1.32 mm d −1 ), who used a modified P-T equation, an A-A equation, and a storage model to simulate ET a for sedum green roofs in a rain-rich environment.The NSEs and PBIASs of ET a simulations for the studied sedums generally fall within the range of Berretta et al. [43] (NSE: −0.15 to 0.98; PBIAS: −74.02 to 15.13) based on a study also using the FAO56 scheme and the Hargreaves equation for sedums in a rain-rich environment; and NSEs of FAO56 and T-M methods are higher than the reported value (0.31) of Starry et al. [34].Overall, the results of the tested three methods fall within the normal range of the aforementioned previous studies.The efficiency of the three methods generally follows the pattern of FAO56 > T-M > P-M for sedums, and FAO56 > P-M > T-M for grass.
All three methods generate a better estimation of grass ET a than sedums ET a .But, the three methods couldn't fully explain the ET variance for the sedum cover, especially for the large ET a rates (over about 5 mm d −1 ).The underestimation of ET a for high ET-rate days has been observed in previous studies that used a fixed annual K c [34,35].The temporal upscaling of the calibrated parameters (i.e., averaging the daily crop coefficients into monthly or yearly values) does not seem to be a satisfactory approach to capture the ET response of sedums to large wetting events.Yearly or monthly constant crop coefficients might not have enough temporal resolution for ET prediction for succulent plants.This issue may be less problematic in wet regions where mild weather with higher humidity usually does not support a high ET rate after rain events.The accuracy of ET a estimation, however, might be more severely affected in dry regions, where frequent irrigation events often occur simultaneously with or right before hot and dry weather that favors high ET rates.Improved (e.g., non-linear) models of K c and K s are recommended for dry regions that depend on irrigation.
Most of the observations (overall 78% for sedums and 73% for grass) fell out of the 90% uncertainty envelope especially during the summer period (Figure 10).The NSE values of the top 10% of the runs ranged from 0.05 to 0.26 for sedums, and from 0.3 to 0.67 for grass.The variation in the selected parameters alone could not explain the disagreement between the estimates and measurements.A similar challenge was also found by Poyatos et al. [83], who modeled ET based on a stomatal conductance model in the Mediterranean climate.This deviation of the observed data from the uncertainty envelope might further indicate the fundamental deficiency in the current model structure [98,99].On the other hand, the estimates also slightly fell out of the uncertainty envelope.This might be due to the lack of seasonality in the GLUE analysis.The sampled parameters with a fixed value were used constantly through a year for one GLUE run.Unlike the calibrated parameters with a finer temporal resolution, the sampled parameters without accounting for the seasonality might not be able to predict the ET a variations with a similar accuracy.As evidence, the observations tend to be closer to the estimates than to the uncertainty bounds (Figure 10).A stratified sampling for different seasons and multi-year observations as the inputs might be necessary for the future uncertainty analysis of ET models.All three methods generate a better estimation of grass ETa than sedums ETa.But, the three methods couldn't fully explain the ET variance for the sedum cover, especially for the large ETa rates (over about 5 mm d −1 ).The underestimation of ETa for high ET-rate days has been observed in previous studies that used a fixed annual Kc [34,35].The temporal upscaling of the calibrated parameters (i.e., averaging the daily crop coefficients into monthly or yearly values) does not seem to be a satisfactory approach to capture the ET response of sedums to large wetting events.Yearly or monthly constant crop coefficients might not have enough temporal resolution for ET prediction for succulent plants.This issue may be less problematic in wet regions where mild weather with higher humidity usually does not support a high ET rate after rain events.The accuracy of ETa estimation, however, might be more severely affected in dry regions, where frequent irrigation events often occur simultaneously with or right before hot and dry weather that favors high ET rates.Improved (e.g., non-linear) models of Kc and Ks are recommended for dry regions that depend on irrigation.
Most of the observations (overall 78% for sedums and 73% for grass) fell out of the 90% uncertainty envelope especially during the summer period (Figure 10).The NSE values of the top 10% of the runs ranged from 0.05 to 0.26 for sedums, and from 0.3 to 0.67 for grass.The variation in the selected parameters alone could not explain the disagreement between the estimates and measurements.A similar challenge was also found by Poyatos et al. [83], who modeled ET based on a stomatal conductance model in the Mediterranean climate.This deviation of the observed data from the uncertainty envelope might further indicate the fundamental deficiency in the current model structure [98,99].On the other hand, the estimates also slightly fell out of the uncertainty envelope.This might be due to the lack of seasonality in the GLUE analysis.The sampled parameters with a fixed value were used constantly through a year for one GLUE run.Unlike the calibrated parameters with a finer temporal resolution, the sampled parameters without accounting for the seasonality might not be able to predict the ETa variations with a similar accuracy.As evidence, the observations tend to be closer to the estimates than to the uncertainty bounds (Figure 10).A stratified sampling for different seasons and multi-year observations as the inputs might be necessary for the future uncertainty analysis of ET models.The cumulative likelihood distribution for the model parameters deviates from the uniform distribution (Figure 11), suggesting that the tested model is sensitive to the selected parameters [77].The current model is sensitive to all three tested parameters via the K-S test at the 0.05 level.While the K-S values of three parameters were found to be close to each other, the model tended to be more sensitive to rs for both sedum and grass (Table 4).

Irrigation Cost
Since it is a major concern associated with the long-term maintenance cost in water-limited regions, irrigation cost was estimated in this study.The total PET amounts of the growing season (April-October) were calculated as 785 mm and 768 mm for growing sedum and grass, respectively.An extra 174 m 3 and 168 m 3 water would be required for sedum and grass green roofs (322 m 2 ), respectively, after assuming that the rainfall (246 mm) could be completely used.The estimated irrigation amounts could only be considered as a minimum, due to the assumption of no water loss The cumulative likelihood distribution for the model parameters deviates from the uniform distribution (Figure 11), suggesting that the tested model is sensitive to the selected parameters [77].The current model is sensitive to all three tested parameters via the K-S test at the 0.05 level.While the K-S values of three parameters were found to be close to each other, the model tended to be more sensitive to r s for both sedum and grass (Table 4).The cumulative likelihood distribution for the model parameters deviates from the uniform distribution (Figure 11), suggesting that the tested model is sensitive to the selected parameters [77].The current model is sensitive to all three tested parameters via the K-S test at the 0.05 level.While the K-S values of three parameters were found to be close to each other, the model tended to be more sensitive to rs for both sedum and grass (Table 4).

Irrigation Cost
Since it is a major concern associated with the long-term maintenance cost in water-limited regions, irrigation cost was estimated in this study.The total PET amounts of the growing season (April-October) were calculated as 785 mm and 768 mm for growing sedum and grass, respectively.An extra 174 m 3 and 168 m 3 water would be required for sedum and grass green roofs (322 m 2 ), respectively, after assuming that the rainfall (246 mm) could be completely used.The estimated irrigation amounts could only be considered as a minimum, due to the assumption of no water loss

Irrigation Cost
Since it is a major concern associated with the long-term maintenance cost in water-limited regions, irrigation cost was estimated in this study.The total PET amounts of the growing season (April-October) were calculated as 785 mm and 768 mm for growing sedum and grass, respectively.An extra 174 m 3 and 168 m 3 water would be required for sedum and grass green roofs (322 m 2 ), respectively, after assuming that the rainfall (246 mm) could be completely used.The estimated irrigation amounts could only be considered as a minimum, due to the assumption of no water loss during the irrigation process.From the water conservation perspective, the non-vegetated roof might be a good alternative to the vegetated green roofs in a dry climate.Based on the tiered water price in Salt Lake City, Utah [100], the total irrigation cost for the studied green roof was $117 for the entire growing season of 2014, which approximately equals to $0.4/m 2 .This irrigation cost is much less than a longer-term estimate ($2.9/m 2 ) for an extensive green roof [41].The difference, however, is highly dependent on the local water price.
Of the various environmental benefits that can be provided by green roofs [17,57], only stormwater runoff reduction is considered in this study.The economic returns of runoff reduction could be reflected by the potential to waive the stormwater impact fee, which is based for this study on the policy of Salt Lake City, UT, USA [101].Given all of the harvested rainwater processed by the green roof, the avoided impact fee for the entire year would be $34.This economic return of using the green roof for runoff reduction would cover only 28% of the irrigation cost.The remaining 72% of irrigation cost would need to be offset by other potential benefits of the green roof, or the irrigation demand would need to be reduced to make projects balance cost and benefits in an economic sense.However, if the native species with the smaller K c [92] are available to grow in this study setting (i.e., assuming the same PET rates), they could live without any irrigation and harvest a minimum of 60 mm of rainwater for the next year's use (equal to 22 m 3 harvested rainwater for the studied green roof).Also, if a retrofitting is possible, harvesting of discharge water from air conditioning units and other nearby sources can be considered solutions.

Conclusions
This study aimed to examine the ET behavior and to parameterize three ET models (FAO56, T-M, and P-M) for green roofs in a rain-scarce environment.A field experiment based on weighing lysimeter units was conducted in Salt Lake City, UT, USA, in 2014.The annual average ET rates of the studied grass, sedums, and non-vegetated covers were observed as 2.69, 2.52, and 2.01 mm d −1 , respectively, which tend to be higher than other studies with more moisture-rich environments.Irrigation was able to support the energy-limited ET regime in the studied rain-scarce region.Given enough irrigation, the ET process is determined more by climate conditions than the wetting events during the dry, hot summers.In terms of water conservation, native species with higher drought tolerance are recommended for the rain-scarce environments with less irrigation demand.Shading and increasing the moisture storage capacity are recommended measures for reducing the irrigation demands for green roof plants with deep rooting systems.
The parameters of the ET models determined by the designed experiment generally fell in the range of other studies conducted in more humid environments.The yearly averages of crop coefficients for the sedum and grass covers were found to be 0.59 and 0.61, respectively, while the yearly averages of surface resistances were estimated as 678 and 651 s m −1 , accordingly.These values are within the reported ranges for rain-rich areas.The physiological parameters of the sedums were found to be close to those of the grass species in this study.This finding indicates that sedums in the studied climate might not necessarily have a significantly lower ET capacity than other grass species as commonly assumed in humid climates.Extra irrigation would be required for sedums to keep them growing healthily in the dry environments.Furthermore, the studied sedums and possibly other shallow-rooted plants could not effectively utilize the deep moisture storage to survive drought periods like native species.On the other hand, the challenge of depleting the deep moisture storage reflects a less satisfactory runoff reduction capability for severe events compared to the long-rooted plants.Therefore, unlike in a humid climate, sedums may not be the ideal green roof species for a drier climate.
The tested three ET models generated reasonable estimates compared to previous similar studies.The FAO56 method overall achieved the best fit.The temporal upscaling of critical parameters like the crop coefficient may lead to the underestimation of high ET rates for a green roof.Compared to humid regions, where the rainfall-driven ET process is prevalent, this issue may lead to more severe errors for ET estimation in drier regions where the overlapping of repeated irrigation events with hot, dry meteorological conditions would often lead to higher ET rates.Improved ET models, and especially the non-linear models with adjustable parameters (e.g., K c and K s ) which can simulate low and high ET ranges, are needed for dry regions with an irrigation demand.
The irrigation cost of the studied green roof was estimated at $0.4/m 2 per year.This cost, however, is in excess of the savings from the avoided stormwater impact fee.At least 72% of the annual irrigation cost needs to be offset by other environmental benefits unless the irrigation demand can be reduced.Native species with the reported smaller crop coefficients may not need any irrigation in the study setting and, therefore, can be recommended for green roof implementations in dry environments.

Figure 1 .
Figure 1.Monthly precipitation and irrigation amounts, and average 2-m air temperature at the study site in 2014.

Figure 1 .
Figure 1.Monthly precipitation and irrigation amounts, and average 2-m air temperature at the study site in 2014.

Figure 2 .
Figure 2. Photograph of the experimental setup.

Figure 3 .
Figure 3. Schematic (not to scale) of the drainage observation setup underneath the lysimeter.

Figure 2 .
Figure 2. Photograph of the experimental setup.

Figure 2 .
Figure 2. Photograph of the experimental setup.

Figure 3 .
Figure 3. Schematic (not to scale) of the drainage observation setup underneath the lysimeter.

Figure 3 .
Figure 3. Schematic (not to scale) of the drainage observation setup underneath the lysimeter.

Figure 4 .
Figure 4. Monthly averages of observed green roof evapotranspiration (ET) rates of 2014 (error bars represent standard deviations).

Figure 4 .
Figure 4. Monthly averages of observed green roof evapotranspiration (ET) rates of 2014 (error bars represent standard deviations).

Figure 5 .
Figure 5. Variable importance in projection (VIP) plots for partial least squares analysis.Note: the VIP value of the influential variables is larger than 0.8 (red dashed line).

Figure 5 .
Figure 5. Variable importance in projection (VIP) plots for partial least squares analysis.Note: the VIP value of the influential variables is larger than 0.8 (red dashed line).

Water 2018 ,
10, x FOR PEER REVIEW 12 of 22

Figure 6 .
Figure 6.Monthly averages of soil moisture at the bottom (left) and the surface (right) soil layers.Data are missing from the non-vegetated lysimeter at the bottom soil in February and March due to broken sensors.The point symbols represent monthly averages (AVG), while the lines represent the monthly standard deviations (STD).

Figure 6 .
Figure 6.Monthly averages of soil moisture at the bottom (left) and the surface (right) soil layers.Data are missing from the non-vegetated lysimeter at the bottom soil in February and March due to broken sensors.The point symbols represent monthly averages (AVG), while the lines represent the monthly standard deviations (STD).

Figure 7 .
Figure 7. Parameterized monthly crop coefficients of sedums and grass surface covers for the Food and Agriculture Organization of the United Nations' irrigation and drainage method (FAO56).

Figure 8 .
Figure 8. Parameterized hourly surface resistances averaged over the year (red solid lines) and averaged over every month (short dashed lines, esp.over June: long dashed lines), and measured hourly surface resistances (blue triangles) of the sedum (top) and grass covers (bottom).* Note: values larger than 2500 s m −1 were excluded.

Figure 7 .
Figure 7. Parameterized monthly crop coefficients of sedums and grass surface covers for the Food and Agriculture Organization of the United Nations' irrigation and drainage method (FAO56).

Figure 7 .
Figure 7. Parameterized monthly crop coefficients of sedums and grass surface covers for the Food and Agriculture Organization of the United Nations' irrigation and drainage method (FAO56).

Figure 8 .
Figure 8. Parameterized hourly surface resistances averaged over the year (red solid lines) and averaged over every month (short dashed lines, esp.over June: long dashed lines), and measured hourly surface resistances (blue triangles) of the sedum (top) and grass covers (bottom).* Note: values larger than 2500 s m −1 were excluded.

Figure 9 .
Figure 9. Evapotranspiration (ETa) estimates based on FAO-56, T-M, and P-M methods plotted against ETa observations for the sedum (left); and grass green roofs (right).Note: FAO56 refers to the Food and Agriculture Organization of the United Nations' irrigation and drainage method, T-M refers to the Thornthwaite-Mather equation, and P-M refers to the Penman-Monteith equation.

Figure 9 .
Figure 9. Evapotranspiration (ET a ) estimates based on FAO-56, T-M, and P-M methods plotted against ET a observations for the sedum (left); and grass green roofs (right).Note: FAO56 refers to the Food and Agriculture Organization of the United Nations' irrigation and drainage method, T-M refers to the Thornthwaite-Mather equation, and P-M refers to the Penman-Monteith equation.

Water 2018 ,
10,  x FOR PEER REVIEW 15 of 22 lead to a higher R 2 value.The RMSEs of simulated sedum ETa are close to Marasco et al.[35] (1.19 to 1.32 mm d −1 ), who used a modified P-T equation, an A-A equation, and a storage model to simulate ETa for sedum green roofs in a rain-rich environment.The NSEs and PBIASs of ETa simulations for the studied sedums generally fall within the range of Berretta et al.[43] (NSE: −0.15 to 0.98; PBIAS: −74.02 to 15.13) based on a study also using the FAO56 scheme and the Hargreaves equation for sedums in a rain-rich environment; and NSEs of FAO56 and T-M methods are higher than the reported value (0.31) of Starry et al.[34].Overall, the results of the tested three methods fall within the normal range of the aforementioned previous studies.The efficiency of the three methods generally follows the pattern of FAO56 > T-M > P-M for sedums, and FAO56 > P-M > T-M for grass.

Figure 10 .
Figure10.Observed and estimated daily evapotranspiration (ETa) with the (5% and 95% quantiles) uncertainty bounds (shaded area) for the sedum (upper); and grass green roofs (lower).Note: the Food and Agriculture Organization of the United Nations' irrigation and drainage method (FAO56) with the parameterized crop coefficients was used here to estimate ETa.

Figure 11 .
Figure 11.Prior and posterior cumulative likelihood distribution for the crop coefficient (Kc) (left); water stress coefficient (Ks) (middle); and surface resistance (rs) (right).Note: the Food and Agriculture Organization of the United Nations' irrigation and drainage method (FAO56) with the parameterized crop coefficients was used here to estimate ETa.

Figure 10 .
Figure10.Observed and estimated daily evapotranspiration (ET a ) with the (5% and 95% quantiles) uncertainty bounds (shaded area) for the sedum (upper); and grass green roofs (lower).Note: the Food and Agriculture Organization of the United Nations' irrigation and drainage method (FAO56) with the parameterized crop coefficients was used here to estimate ET a .

Water 2018 , 22 Figure 10 .
Figure10.Observed and estimated daily evapotranspiration (ETa) with the (5% and 95% quantiles) uncertainty bounds (shaded area) for the sedum (upper); and grass green roofs (lower).Note: the Food and Agriculture Organization of the United Nations' irrigation and drainage method (FAO56) with the parameterized crop coefficients was used here to estimate ETa.

Figure 11 .
Figure 11.Prior and posterior cumulative likelihood distribution for the crop coefficient (Kc) (left); water stress coefficient (Ks) (middle); and surface resistance (rs) (right).Note: the Food and Agriculture Organization of the United Nations' irrigation and drainage method (FAO56) with the parameterized crop coefficients was used here to estimate ETa.

Figure 11 .
Figure 11.Prior and posterior cumulative likelihood distribution for the crop coefficient (K c ) (left); water stress coefficient (K s ) (middle); and surface resistance (r s ) (right).Note: the Food and Agriculture Organization of the United Nations' irrigation and drainage method (FAO56) with the parameterized crop coefficients was used here to estimate ET a .

Table 2 .
Range of the prior uniform distribution sampled for the uncertainty analysis.

Table 3 .
Statistical performance measures of evapotranspiration (ET a ) estimations.

Table 4 .
Kolmogorov-Smirnov coefficients for the sensitive parameters at the 0.05 level.

Table 4 .
Kolmogorov-Smirnov coefficients for the sensitive parameters at the 0.05 level.

Table 4 .
Kolmogorov-Smirnov coefficients for the sensitive parameters at the 0.05 level.