Next Article in Journal
Effects of the Concept of the Dualist Division of Lakes in Polish Law on Shaping Lake Water Policy from the Hydrological Perspective
Previous Article in Journal
A Multi-Criteria Assessment of Water Supply in Ugandan Refugee Settlements
Previous Article in Special Issue
Evaporation from (Blue-)Green Roofs: Assessing the Benefits of a Storage and Capillary Irrigation System Based on Measurements and Modeling
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Monitoring and Modeling the Long-Term Rainfall-Runoff Response of the Jacob K. Javits Center Green Roof

1
CAEE Department, College of Engineering, Drexel University, Philadelphia, PA 19104, USA
2
Department of Civil Engineering, Albert Nerken School of Engineering, The Cooper Union, New York, NY 10003, USA
3
Arcadis of New York, Inc., Long Island City, NY 13057, USA
*
Author to whom correspondence should be addressed.
Water 2018, 10(11), 1494; https://doi.org/10.3390/w10111494
Submission received: 15 June 2018 / Revised: 16 August 2018 / Accepted: 25 September 2018 / Published: 23 October 2018
(This article belongs to the Special Issue Hydrological Performance of Green Roofs)

Abstract

:
Drainage from the 27,316-m2 Jacob K. Javits Convention Center (JJCC) green roof was investigated in the field to quantify the system’s long-term rainfall-runoff response. The JJCC hosts one of the largest extensive green roofs in the United States. Utilizing four years of rooftop monitoring data collected using a weather station, custom designed and built drainage systems, three Parshall flumes equipped with pressure transducers, and weighing lysimeters, this study quantified the 25.4-mm-deep green roof’s ability to decrease the volume and peak rate of runoff. With parameters derived from the site, the Environmental Protection Agency Stormwater Management Model (EPA-SWMM) predicted event total runoff volume and event peak runoff rates to within +10% to −20% and +25% to −15% of the observations, respectively. The analysis further indicated that approximately 55% of the cumulative precipitation that fell on the JJCC extensive green roof during the monitoring period (warm weather months, June 2014–November 2017) was captured and retained. The average percent retained on an event-basis was 77%, and average event runoff coefficient was 0.7, implying a substantial reduction in the volume and rate of runoff generated from the roof compared to the pre-green roof condition, when most, if not all, of the precipitated water would have immediately resulted in runoff. Our research suggests that, on average, 96% of rainfall events 6.35 mm or less were retained within the green roof, whereas 27% of the total event volume was retained for events greater than 12.7 mm in depth. A sensitivity analysis suggests if the substrate depth were increased, better stormwater capture performance would be achieved, but only up 127 mm, whereas increased precipitation coupled with warmer temperatures as a result of climate change could decrease the performance by up to 5%, regardless of substrate depth. An equivalency analysis suggested that even shallow green roofs can significantly reduce the required stormwater detention volume that New York City requires on new development. This particular green roof appears to be more than 18 times as cost-effective as a subsurface cistern would be for managing an equivalent volume of stormwater in Midtown Manhattan.

1. Introduction

Built even in ancient times [1,2], green roofs have emerged as one of the most effective methods of distributed storm water control in contemporary cities, particularly where storm water runoff can cause flooding, sewer surcharges and overflows, and pollutant loading to receiving waters [3,4,5,6]. In dense cities with limited undeveloped space, one key advantage of green roofs over other types of green infrastructure (GI) is that they can be retrofitted onto existing buildings [7]. As the pilot phase of GI implementation in many US cities comes to an end, street-level runoff capture sites are increasingly scarce, making rooftops, which typically account for 20–25% of impervious cover in American cities [8,9], one of the next big opportunities for decentralized urban storm water management. Precipitation incident to green roofs is intercepted, detained, and evapotranspired, reducing the volume of roof runoff while also delaying and lowering the peak rate of discharge [10,11,12]. A review of some of the key hydrologic processes at work in green roofs was provided by Stovin et al. [13].
Green roofs can be classified based on the thickness of the substrate layer as either intensive (substrate thickness > 150 mm) or extensive (substrate thickness < 150 mm) [10,14,15]. Intensive green roofs with thicker soil layers can accommodate larger vegetation resulting in higher retention rates, but often require building reinforcements or retrofits to support the additional weight, which also increases the initial costs of installation. Extensive green roofs have much thinner substrate layers and can be installed at a lower cost, often without the need for reinforcements [7,16].
The hydrologic performance of green roofs has been extensively studied in both monitoring and modeling studies. Fifty of the more than 1200 Web of Science papers referencing green roofs investigated their hydrology [10,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25]. Two decades of early European (mostly German) green roof monitoring was reviewed by Mentens et al. [14]. A decade of literature from field and controlled experiments was reviewed by Berndtsson [10] and later updated by Li and Babcock [19] and Carson et al. [20]. The most recent summary of studies on the hydrologic performance of green roofs was published in 2016 by Cipolla et al. [21]. These studies generally suggest that the performance of individual systems varies widely due to differences in local climate and microclimate [22]; incident precipitation amounts, intensities, and durations [14]; the type [23] and depth of the growing medium [14,23]; and the presence [23] and type [24,19] of vegetation. Seasonal differences have also been observed, with more runoff leaving the systems during winter [25,14].
Many computer simulations of green roof hydrology have also been attempted. Modeling approaches have advanced in number and sophistication since 2007, when Elliot and Trowsdale [26] reported that only one model (e.g., Water Balance Model) could be used to model green roof hydrology explicitly. Recent modeling efforts represent green roof storm water attenuation using runoff coefficients or curve numbers [27,28,29,30] and use more elaborate physically-based models such as HYDRUS-1D (Prague, Czech Republic) [31,32,33] to represent saturated and unsaturated flow patterns.
Given its extensive use in contemporary urban storm water planning, the Environmental Protection Agency Storm Water Management Model (EPA-SWMM) has also been used to simulate green roof hydrology. In early versions of the software, green roofs could only be represented in ad hoc ways, for example, using storage nodes [27]. SWMM 5.1 and subsequent versions, include the Low Impact Development Green Roof (LID-GR) module to explicitly represent the properties of a green roof’s various layers: surface, media, and drainage layer. Despite the widespread use of SWMM’s LID-GR in planning studies throughout the world [34,35,36,37], few attempts to validate it with actual monitoring data have been attempted. Burszta-Adamiak and Mrowiec [38] used another SWMM LID control (the bio-retention control or LID-BR) to perform single-event simulations from a small-scale (~3 m2) system, with mixed success. Calibration of the model to match the observed runoff volume greatly overestimated the peak rates of runoff, whereas calibration of the peak rates underestimated the volume. Though both the peak rate and total volume of runoff were underestimated, Palla and Gnecco [39] had greater success. Their single-event simulations were also based on observations from small-scale systems (~2.5 m2) but achieved Nash Sutcliff Errors of >0.78 using SWMM’s LID-GR.
In continuous simulations, additional provisions are required to adequately represent evapotranspiration between rain events. Cipolla et al. [21] used the LID-BR in SWMM to model a full-scale (120 m2) green roof in a one-year-long continuous simulation. Daily potential evapotranspiration was estimated using the Hargreaves equation based on local air temperature and radiation, modified to account for the sedum found at the site, using monthly soil recovery coefficients developed by the authors. Good agreement with experimental observations was found. Krebs et al. [40] used SWMM’s LID-GR control to perform continuous simulations of small-scale (2 m2) for growing seasons of two different years, including inter-event periods. SWMM’s Potential Evapotranspiration (PET) values were corrected using a coefficient to account for the influence of vegetation on evapotranspiration (ET) rates. Good agreement with the observations was found, though the model results were highly sensitive to the soil and drainage layer parameters and the ET coefficient, and hydrograph recession limbs were not well simulated. The authors suggested that crop coefficients for shorter periods of time should be attempted.
The goals of this paper were: (1) to present long-term rainfall-runoff observations from three 186 m2 drainage areas of the Jacob K. Javits Convention Center (JJCC) green roof (40.7578° N, 74.0027° W) in New York City; (2) to further investigate the accuracy with which SWMM’s LID-GR control can be used to represent long-term rainfall-runoff relationships in continuous simulations when PET values are adjusted for moisture availability; and (3) to use the calibrated model to explore the storm water capture implications of different depths and configurations of green roofs in New York City. The research is novel for several reasons. First, the size of this green roof (currently the second largest in the United States) and its many internal roof drains required installation of a unique monitoring system to direct runoff from three 186 m2 custom-designed rooftop drainage areas into Parshall flumes for runoff measurement. Second, to the best of our knowledge, no prior research has attempted to validate SWMM’s LID-GR control in continuous simulations of runoff from a large-scale green roof with ET rates adjusted for soil moisture availability. Finally, the study is novel because it is situated in New York City, which has a $1.4 billion USD plan [41] to capture 25 mm of precipitation over 10% of the city’s impervious surfaces using various forms of GI, including green roofs. As opportunities for street-level GI implementation become fewer, the city is interested in cost-effective strategies for managing stormwater from its more than 56 million m2 of rooftops [41]. It has been estimated that by greening all rooftop space in New York City air temperature could be reduced by an average 1.6 °C, annual stormwater reduced by 38 billion liters, and energy cost savings amounting to roughly $130 million USD each year achieved [12,41]. Research documenting the hydrologic benefits of green roofs in New York City’s unique urban climate is a first step toward evaluating whether such goals are achievable.

2. Materials and Methods

2.1. Site Description

The Jacob K. Javits Convention Center’s 27,316 m2 extensive green roof (JGR) was completed in the spring of 2014, stretching across both the northern and southern wings of the building (Figure 1A). In section, the design consists of a pre-vegetated sedum mat installed on top of a 25.4 mm growing medium, a retention fleece layer, a drainage layer, and a root barrier, as described in Alvizuri et al. [42]. The XF301 pre-vegetated sedum mat (XeroFlor, New York, NY, USA) was grown and obtained from various nursery farms in upstate New York, South Carolina, and the mid-west United States. The sedum mat is an integrated unit of plant materials and growing medium meshed onto a geotextile combing these components.
Per the manufacturer, the XeroTerr (XeroFlor, New York, NY, USA) growing medium is a low-organic/high-mineral content growing mix composed of composted plant materials and lightweight porous aggregate. Its mineral fraction includes heat-treated clay, slate, or shale and geologically-formed pumice. Its organic content, estimated at 7%, is formed from compost that is aged sufficiently to prevent unwanted microbial growth, weed seed survival, or re-sprouting of plant litter debris. Its grain size distribution can be represented by the percent of a typical batch passing the #4, #8, #16, #30, #50, and #100 sieves. These values are 75–95%, 35–80%, 15–50%, 10–25%, ≤10%, and ≤5%, respectively.
The XF157 (XeroFlor, New York, NY, USA) water retention fleece is a fabric composed of recycled fiber materials that provide a tough robust layer for plant roots to latch on and grow. The fleece also serves as a filter layer to screen water that may carry minerals from the growing medium layer or from the precipitation. The XF108H (XeroFlor, New York, NY, USA) drain mat is a geotextile layer attached to a mesh of polymeric fibers (2 mm thick) designed to create a clear space for excess water to drain out from the topmost layers.
Finally, the XF112 (XeroFlor, New York, NY, USA) root barrier is a nearly water-impermeable layer designed to deter the plant roots from establishing themselves any further down into the actual roof structure.
Runoff generated on the JGR flows into a series internal roof drains, spaced uniformly throughout the roof. Each roof drain conveys runoff generated by its associated tributary area to internal drain pipes, eventually discharging to the local combined sewer.

2.2. Overall Description of Monitoring Setup

A large number of sensors were installed on the JGR as part of a much broader investigation of its environmental performance, including a study of its ability to provide thermal buffering [42]. Only those sensors used in the rainfall-runoff study conducted on the south roof are described here.

2.2.1. Weather Station, Drainage Channels, Flumes, and Weighing Lysimeters

A weather station positioned on the south roof logs climatic data at five-minute intervals on a Campbell Scientific, Inc. (Logan, UT, USA) CR1000 data logger. The station is equipped with a Texas Electronics, Inc. (Dallas, TX, USA) Series 525-Rainfall Sensor and other sensors to monitor temperature, net-radiation, relative humidity, wind speed, and wind direction, as summarized in Table 1.
During installation of the JGR, custom-built aluminum drainage channels were affixed to the structural roof surface. The purpose of these channels was to clearly delineate three 186 m2 tributary plots for runoff measurement (Figure 1B). Each monitored tributary plot has a 27.4-m-long transverse and 2.4-m-long longitudinal section, forming a “T” in the plan view (Figure 1C). Each transverse section was assembled from 2.4-long built-up aluminum box components with 11 boxes constituting the transverse section and one box for the longitudinal section (Figure 1E). This design allows the drainage channels to conform to the irregular slope of the roof top. The transverse section of each channel has a 30.5 × 2.86 cm cross section, whereas the longitudinal section has a 22.2 × 5.7 cm cross section (Figure 2). These built-up boxes were constructed of aluminum 3003-H14 sheets and aluminum 6061-T6 sheets, 16 mm and 3.2 mm thick, respectively. The box sections are reinforced with aluminum 6063-T52 U-channels and aluminum Zee channel members. The Zee channel members are arranged with an offset to provide increased structural reinforcement of the box while allowing rain runoff to flow unimpeded within the box. The U channels are used to form the downstream edge of the boxes and are sealed tightly.
Parshall flumes (Figure 3) equipped with pressure transducers were used to measure discharge from the end of each of the three channel systems (Point B in Figure 1C). A Parshall flume imposes a critical flow condition within its throat so that the open channel discharge through it can be unequivocally obtained from measurements of the upstream water level. In this case, the upstream water level was measured continuously, also at 5 min intervals, with Campbell Scientific (Logan, UT, USA) CS451 pressure transducers, installed in a 51-mm-diameter polyvinyl chloride (PVC) well attached to the side of the Parshall flume. The pressure transducer was placed at the bottom of this well and continuously recorded the depth of water in the well to the data logger. The water level depths were converted to flow through the flume using a stage discharge curve obtained from the manufacturer:
  Q = 0.107 · H 1.55  
where Q is flow through the flume in m3/s and H is the upstream water level in meters.
The particular Parshall flume model selected for this study ((01-007) from Global Water (College Station, TX, USA), a division of Xylem (Rye Brook, NY, USA)) can measure flow rates between 2 × 10−5 and 1.92 × 10−2 m3/s. This range was considered appropriate after estimating the peak flow rates expected from each 186 m2 plot under New York City’s 5-year, 15 cm/h design storm, which would yield 5.5 × 10−3 m3/s over each flume’s tributary catchment area. The flume is constructed of a durable, polyester resin, and fiberglass, has a 5.08 cm (2”) throat width, and is 77.5 cm long, 34.8 cm high, and 21.4 cm wide.
Three weighing lysimeters were installed in each of the flume’s three tributary areas. The weighing lysimeters consist of a square piece of green roof matrix enclosed in a 0.372 m2 metal box. The metal box rests on a custom 0.372 m2 Rice Lake Roughdeck (Rice Lake, WI, USA) mild steel platform scale equipped with four mild steel load cells. The maximum capacity of the scale is 227 kg with a 0.02 kg resolution. The lysimeter system extends 0.19 m above the bare roof, corresponding to a total extension of 0.089 m above the green roof. Although physically isolated from the green roof, each lysimeter is surrounded by four sloped transition sections of the green roof to minimize thermal and heat boundary effects. The metal box is tilted on the weighing scale at a similar slope as the green roof to allow drainage from the lysimeter during a rain event.
The entire rainfall and runoff monitoring system was fully operational by June 2014. Rainfall and runoff observations were recorded for four monitoring seasons, roughly extending from July 2014 (when the monitoring setup was complete) through the end of August 2017. Monitoring was not conducted during the winter months to avoid frost damage to the pressure transducers.

2.2.2. Data Analysis: Event Discretization and Percent Retained

The continuous precipitation time series was discretized into individual storm events using a four-hour inter-event period, following the approach described in Yu et al. [43]. The four-hour inter-event period is the standard duration used in drainage studies in New York City. The individual rain events were then further classified into three categories based on the total rainfall: low rainfall (less than 6.35 mm), medium rainfall (6.35 to 12.7 mm), and heavy rainfall (greater than 12.7 mm), per the American Society of Civil Engineers (ASCE) Hydrology Handbook [44].
Runoff events from each of the flumes were associated with each of the storm events by visual inspection. For each event, the percent retained (%R) was computed as:
  % R = [ 1 D i s c h a r g e   V o l u m e R a i n f a l l   V o l u m e ] × 100  
where the discharge volume is the sum of the incremental discharge computations made for every non-zero water level measurement in the flumes following the initiation of rainfall. Because each of the tributary areas were equivalent in size, data from all three flumes were pooled.

2.3. Rainfall-Runoff Modeling Approach

A model of the green roof was created using EPA-SWMM [45], a dynamic hydrologic-hydraulic model that is used for long-term simulation of runoff quantity from impervious surfaces in urban areas. SWMM is an industry standard model, widely used in storm water planning in New York City and throughout the world. Recent versions of SWMM can be used to predict the performance of certain types of low impact development (LID) controls, including green roofs. In this application, and in contrast to several other recent attempts to model green roofs in SWMM [38,21], a single 186 m2 subcatchment area was 100% occupied by a LID-GR (not bioretention) control to represent one of the three identical tributary areas. Subcatchment properties used in SWMM are shown in Table 2.
SWMM requires two climate time series inputs: precipitation and ET. In this study, precipitation data were sourced from the onsite weather station. In continuous simulations, ET input data is especially critical because it establishes the initial moisture regime for each simulated event. Researchers [39] have identified ET as a critical determinant of simulated storm water capture performance. Because actual ET rates in the Northeast U.S. are typically constrained by moisture availability and not energy, a two-step process was developed to generate a time series of Actual Evapotranspiration (AET) time series for use in the SWMM simulations. The first step involved computing a PET time series by inserting the monitored rooftop air temperature (T) and solar radiation (Ra) values into the Hargreaves equation [32]:
  P E T ( t ) = 0.0023 × R a ( t ) × ( T m a x T m i n ) 0.5 × ( T a v g + 17.8 )  
where Tmax (°C) is the maximum daily air temperature, Tmin (°C) is the minimum daily air temperature, and Ra (MJ m2/day) is the extra-terrestrial solar radiation.
Next, an intermediate SWMM run was performed using the PET values, the goal of which was solely to export a predicted soil moisture (θ) time series for the green roof growing media. These soil moisture values were then used to derive AET values using the Thornthwaite-Mather equation [48]:
  A E T ( t ) = P E T ( t ) × ( ( ( t ) W P ) / ( F C W P ) )  
where θ is the soil moisture content, WP is the wilting point (Table 3), and FC is the field capacity (Table 3). The final model predictions were performed with this AET time series uploaded into SWMM’s climate editor.
Other model properties were derived through a sensitivity analysis using the Comprehensive R Archive Network (CRAN) Sensitivity Package in R [49]. The sensitivity analysis varied the flow width, depression storage, and roughness coefficients on a trial-and-error basis to maximize the number of predicted values of total volume and peak flow for each event were within ranges: +20% and −10%, and +25% and −15% of the observed values, respectively. These ranges are typically considered acceptable in urban drainage studies [50]. The initially assumed, and final calibrated, parameter values presented are in Table 3.

2.4. Model Validation

In addition to the parameter calibration described above, model validation was accomplished by performing a water budget analysis on each event. The purpose of this water budgeting exercise was to confirm that the volume of precipitation falling on the roof could be fully accounted for in the measured discharge, the change in soil moisture of the green roof (as represented by the increase in lysimeter mass over the event), and the evapotranspiration per Equation (5):
  P = R O + Δ S + E T  
where P is precipitation per event, RO is runoff per event, and Δ S is change in lysimeter mass, more extensively described in Smalls-Mantey et al. [52]. Because conditions were saturated during the rain event, ET in this case was computed using the Penman-Monteith equation [36].
All values were calculated in volume per area of green roof bay (m3 per green roof bay). A standard error of estimate equation was used to quantify the accuracy of the data collected. This model measures the level of scatter of N data points about a regression line between two variables X and Y via:
  S y x = ( Y Y e s t ) 2 N  
where Yest is the observed precipitation volume, and Y is the sum of the observed storm runoff, ET, and the storage change in the substrate. Y and Yest were plotted such that the 45° diagonal of the graph represents values that satisfy the water balance equation. Points that lie elsewhere indicate a discrepancy between Y and Yest and will have non-zero water balance.

2.4.1. Statistical Analysis

The SWMM model, calibrated to the observations, was then used to investigate which physical and climatic conditions were most strongly correlated with the roof’s storm water capture performance. New model runs were performed by incrementally changing individual parameter values. A binary logistic regression analysis of the impact of four independent variables on both annual and event-based performance was carried out using IBM SPSS [53]. Logistic models are well suited for predictor variables that are a mix of categorical and continuous variables, or for variables that are not normally distributed. Binary logistic models are used for dependent variables that are dichotomous. In this analysis, the dependent variable was coded as a binary variable indicating whether there was an increase or a decrease in green roof performance. The variables defined are as follows.
Green Roof Performance (GRP) (dependent variable): defined as the change (increase or decrease) in the percent of runoff retained by the green roof.
Climatic Conditions (CC): precipitation and temperature data were varied in the model based on climate projections to represent different climate change scenarios. Projected climate datasets were obtained using the Climate Adjuster Tool (CAT) in SWMM [54]. CAT allows for the incorporation of future climate change projection into a SWMM model by adding a set of monthly adjustment factors to each climatic time series. SWMM allows for the user to manually insert monthly adjustment factors to multiply to any climate data series; however, the CAT provides location-specific adjustment factors based on the World Climate Research Programme (WCRP) Coupled Model Intercomparison Project Phase 3 (CMIP3) archive. For this simulation, near-term (2020–2049) warm/wet climate projection adjustment factors were used, reflecting the climate model that represents the lowest annual temperature and the highest annual rainfall.
Design Parameters (DP): the depth of the soil layer was varied in the model within the thickness range for extensive green roofs (30–150 mm). Since SWMM allows users to add the depth of a fleece layer, its thickness is assumed to be included in the soil depth.

2.4.2. Runoff Coefficients and Green Infrastructure Equivalency

The validated model presented an opportunity to quantify the effective runoff coefficient for green roofs of different depths. Runoff coefficients were calculated from the observed peak flow rates for the actual JGR and for hypothetical green roofs with other depths obtained from the SWMM model by modifying the substrate depth property of the LID-GR control. For each of the observed and simulated events, the runoff coefficients were computed using the Rational Equation:
  C = A i / Q  
where C is the runoff coefficient (dimensionless), Q is peak runoff flow rate (cm/s), A is area (ha), and i is rainfall intensity (mm/h).
Individual C values were developed for each event and box plots used to represent the distribution of C values obtained for each green roof configuration (e.g., actual and simulated).
Because it is an existing building, the Javits Center technically is not required to manage any storm water onsite. However, a new building of this size would be required to detain storm water onsite following a procedure enforced by the New York City Department of Environmental Protection (NYCDEP) [55]. The observed and predicted data sets offered a unique opportunity to quantify the percent reduction in required detention volume that could be achieved by green roof installation. This analysis was performed for different green roof depths and roof areas by following the procedure outlined in the NYCDEP guidelines [55]. A series of empirical equations are provided that require the analysis be performed in English units. The first step is to calculate the duration of the storm in minutes with a 10-year return frequency requiring the maximum detention volume with a variable outflow, TV:
  T V = 0.27 / ( C A / Q D R R ) 0.5 15  
where C is the runoff coefficient (dimensionless), A is roof area (ft2), and QDRR is detention facility maximum release rate (assumed to be 0.25 cfs [40]).
Next, the maximum required detention facility volume is computed:
  V V = [ 0.19 C A / ( t V + 15 ) 40 Q D R R ] t V  
where VV is the maximum required detention volume (ft3). All other parameters are as defined previously.
The reduction in tank volume computes the percent reduction in VV between a roof with C of 0.95 (membrane roof), and a roof with a mean C value corresponding to different depths of green roof. The maximum storage depth is then calculated as follows:
  S D R = 1930 [ Q D R R / ( 3 2 ) ] 2 + ( 3 / 24 )  
The floor area of a required detention tank is equal to VV/SDR.

3. Results

All the raw data collected on the Javits Center by the research team can be accessed on a web portal (http://vistadv.cae.drexel.edu/vdv/index.html), with a periodically updated password that can be provided on request by contacting the authors.

3.1. Monitoring Results

A total of 118 pooled rain events were recorded including 31 events in 2014, 32 events in 2015, 24 events in 2016, and 31 events in 2017 (Table 4). Categorizing events by intensity, the study period included 53 light events, 27 medium events, and 38 heavy events (Table 5). The overall average percent retention for all events was 77.3%, whereas the mean annual event retention varied from 75.4% to 79.3% (Table 4). Approximately 55% of the total cumulative precipitation that fell on the JGR during the monitoring period (June 2014–November 2017, warm weather months) was captured and retained. Annually, event performance increased over time at a rate of 3.82% retention per year, though these increases were not found to be statistically significant through a T-test of means.
As part of the analysis, the percent retained values of each event were regressed against their Julian date, as well as against the maximum and minimum daily air temperature recorded on the day of each event. Although other researchers [14,25] found a seasonal relationship with green roof storm water retention, no such relationship was found here. The percent retained values were also regressed against the event precipitation depth, duration, and mean intensity, also without any significant statistical relationship evident.

Sample Hydrograph

To demonstrate the accuracy with which the SWMM LID-GR simulation matched the observed flows, Figure 4 shows the observed and predicted runoff and cumulative precipitation for Event 83, a heavy event that occurred on 27 October 2016. This event’s antecedent dry period was 18 days—the longest of all events included in this study. The storm resulted in 25.6 mm of precipitation, equivalent to 4.8 m3 over the 186 m2 tributary area. It resulted in 2.85 m3 of observed discharge. The observed retention rate of the green roof for this storm was calculated, using Equation (2), to be 41%. The SWMM model predicted 3.35 m3 of discharge, or a retention rate of 30%.

3.2. SWMM Model Validation

Tornado plots showing the role that each of the independent variables tested in the sensitivity analysis had on total and peak volume are shown in Figure 5 and Figure 6, respectively. The flow width had the greatest impact on both performance metrics.
The observed and modeled runoff volume and peak flowrate for each event were plotted against each other to determine whether the model accurately predicted the observations within an industry standard validation envelope: +20% to −10% for total volume, and +25% to −15% for peak flow [50]. The peak flow envelope is wider than the total volume envelope since model calibration generally identifies parameters best suited to the overall hydrograph, and not solely the timestep corresponding to peak flow. Figure 7 and Figure 8 show the validation of total volume and peak flow, respectively, for all events. Comparing observed and modeled runoff volume, most events fell roughly within the acceptable validation range with few exceptions (Figure 7). Peak runoff values generally fell within the acceptable range, with the model overpredicting the peak, resulting in the data points skewing toward the y-axis in Figure 8.

3.3. Water Balance

As discussed in the Methods section, a statistical standard error was used in concert with a water balance equation to validate the recorded runoff response and compare it to the onsite precipitation data. The black 45° diagonal line in Figure 9 represents the water balance as defined in Equation (5), whereas the red lines represent the standard error offsets of ±0.73 m3 from the 45° black line from Equation (6). A total of 80% of the events fell within the red lines. Events that fell outside the region outlined by the red lines strongly suggest the presence of discrepancies between the flume, and/or rain gage data. Note that using the Penman-Monteith approach, the maximum intra-storm ET that could have occurred for any of the events described in this report is 0.18 m3 per plot, or less than 3% of the total contribution to the water balance. Because this volume represents a small fraction of the overall water balance, it was ultimately neglected in the water balance equation. Outlier rain events should be examined with caution.
Figure 10 depicts the semi-log relationship between the change in lysimeter mass, an indicator of substrate water storage measured in kg, and the antecedent dry period (ADP) time in hours for heavy storms with ADP ≥ 24 h. A 24-h period was selected because the Natural Recourses Conservation Service (NRCS) data states that a heavy storm was one with a total precipitation of 12.7 mm (0.5”) or greater in one day. This relationship was also studied for light and medium storms that have precipitation values of less than 12.7 mm and little to no runoff. In Figure 10, the R2 of the black, semi-log best fit line is 0.594, with 11 of the 19 points falling within the standard deviation of 0.34 kg (indicated with red lines). The equation of the semi-log best fit line is:
  L = 0.4705 · ln A 1.093  
where L is the change in lysimeter mass in kg and A is the ADP in hours. With this equation, the change in lysimeter mass (e.g., the retention) can be predicted using the antecedent time, and the event discharge can be determined using the water balance in Equation (5). This computation was performed for over 100 storms and 80% of the events fell within the standard deviation.

3.4. Statistical Analysis

The results of the linear regression revealed that two of the co-variates examined (CC and DP) were significant predictors of the performance of the green roof (GRP). The depth of the substrate (DP) revealed the strongest relationship with performance: each 25.4 mm (1 inch) incremental increase in the depth of the substrate resulted in an increase in the likelihood that annual performance for the green roof will improve by an average of 4.4% (Exp(β) = 0.596, p < 0.001). The climate change adjusted scenario simulating warm/wet future climate resulted in decreased green roof performance by 5.4% (Exp(β) = 0.764, p < 0.005).
The effects of each of the two co-variates that proved statistically significant predictors of green roof performance are visualized in Figure 11, where the average percent retained of all events in each scenario is plotted against different media depths. The change in performance with different green roof depths is examined with the base model used in this analysis, then once again under climate change adjusted conditions. A simulated decrease in the depth of the vegetated medium from 25.4 mm (1 inch) to 12.7 mm (0.5 inch), as could occur, for example, through wind erosion of the growing medium (observed in some portions of the JGR), resulted in a statistically significant increase in the likelihood that performance will decrease. Retention for both scenarios presented in Figure 11 decreased by an average of 9% when the depth decreased to 12.7 mm.

3.5. Green Infrastructure Equivelancy

Quantification of the runoff coefficient associated with different media depths (25.4 to 127 mm) produced values ranging from 0.24 to 0.75, consistent with the literature [56]. Figure 12 includes boxplots of the event C values all events occurring in the simulations with green roofs of different depths. The first box represents the event runoff coefficient for the observed data for the JGR; the overall average event runoff coefficient for the 25.4 mm (1 inch) green roof was 0.7.
Replacing the green roof with an equivalent detention facility managing the same amount of runoff as the JGR would require a large volume cistern that would require a significant area of potentially developable floor space. For context, real estate in this section of Manhattan is currently selling at $1 billion USD per city block (Alan Steel, personal communication with Javits Center CEO). Figure 13 shows the percent reduction in detention facility volume that would be required of NYC buildings to retain storm water equivalent to green roofs of varying depths (12.7 to 127 mm) and varying green roof area (1250 to 30,000 m2). The size of the required detention facility to replace a green roof, yet manage the same amount of runoff, ranges from 11 to 1350 m3, depending on the area and depth of the growing medium. The detention facility volume required to replace the 25.4 mm (1 inch) deep, 27,000 m2 JGR was calculated to be 974 m3. At an assumed effective storage depth for the facility of 1.6 m, a detention facility with approximately 608 m2 of floor area would be required. A rough estimate of the cost to construct such a facility in Midtown Manhattan is ~$2.5 million USD (assuming approximately $4100 USD per m2 of subgrade foundation and cistern waterproofing).

4. Discussion

Despite the large area of the JGR, the rainfall retention values and runoff coefficients monitored on the three study plots are in the same range as reported for other green roofs [10,14,57]. Overall, 55% of the cumulative precipitation that fell on the green roof during the monitoring period was retained. The mean percent of event precipitation retained ranged from 75.4% to 79.3% over the four years of observations. Carson et al. [20] ran tests on three roof models at a slope of 2% and found that vegetated roofs retained 60.6% of rainfall, media-only roofs retained 50.4%, and a gravel ballast roof retained 27.2% of rainfall. For rain events resulting in less than 10 mm of precipitation, Simmons et al. [58] found similar retention percentages to those presented here, with values varying from 26% to 88% depending on the rain size and intensity.
This modeling approach builds on previous efforts to model green roofs in SWMM [21,40], by comparing the results of a four-year continuous simulation performed using SWMM’s LID-GR control (not its bioretention control) with observations for a large-scale system (186 m2). In contrast to previous work, we adjusted SWMM’s ET time series to account for the regularly occurring moisture limited conditions found in this climate. This adjustment was accomplished using a unique two-step process and SWMM’s own soil moisture accounting routine.
Although this approach accounts for moisture limitations, it did not incorporate any vegetation-specific crop coefficients, as previously attempted by others [40]. A potential next step would be to integrate both strategies. During high moisture conditions, ET rates obtained using the rooftop weighing lysimeters could be compared to reference ET measurements obtained using the ASCE standardized reference evapotranspiration equation to derive site specific crop coefficients. Next, these crop coefficients could be used to adjust the full time series of reference ET values obtained from the weather station. By running the model with this ET time series, interim soil moisture measurements could be exported from SWMM’s LID report and used to adjust the ET values for moisture limitations. It would be interesting to compare the results of our approach as part of the present effort to this proposed new approach to determine whether further improvements in model accuracy could be achieved.
Additional research could validate SWMM’s simulated soil moisture time series with substrate moisture content observations obtained using soil sensors installed on the Javits green roof. The sensors would first need to be calibrated to the actual substrate and rooftop data loggers, since differences in the logger’s excitation frequencies can alter the sensor calibration curve. To validate SWMM’s predicted soil moisture routine over various seasons and moisture conditions, sensor calibrations would need to be repeated periodically because these relationships have a tendency to vary over time (author observation based on work at other locations). Although soil sensors were installed on the JGR prior to the reporting period described here, sensor calibration activities were not frequent enough to ensure data quality.
The regression analysis suggests that the simulated depth of the green roof substrate layer played a significant role in predicted storm water retention rates over the four-year continuous simulation. Predicted retention increased significantly with each initial 25.4 mm (1 inch) increase in substrate depth, with less dramatic improvements in performance when the depth increased above 101.6 mm. However, these findings need to be interpreted in the context of other research regarding the role of substrate depth. In their review of observations in a suite of European green roofs, Mentens [14] found no relationship between seasonal runoff and substrate depth during the cold and transitional periods (November to April). However, they found that every 10 mm increase in substrate depth resulted in 2.5 mm less runoff (p < 0.05) during the warm seasons. ANOVA indicated that green roofs with substrates <50 mm and >150 mm had statistically significant differences in warm season runoff. In their study of small-scale (62 cm2) Midwest U.S. green roofs, Morgan et al. [23] found no significant difference in retention over an 18-month study for green roofs ranging from 5 to 20 cm in depth. During a second study, they found a small and statistically insignificant difference in runoff retained from a five-cm green roof, when compared to the deeper roofs.
Additional research might focus on the variable role of substrate depth under different climatic conditions. Deeper substrates theoretically can store more water. However, during cold seasons (when ET is low) in climates with regular and frequent precipitation, the substrate is more likely to be wetter at the initiation of rain. If the moisture content is near or at its field capacity, runoff from the green roof surface could be expected to occur through saturation excess, regardless of substrate depth or precipitation depth or intensity. This phenomenon, which could not be investigated as part of this study since no data collection was performed during the winter, would serve to diminish the importance of substrate depth as a design parameter. Runoff could also occur regardless of substrate depth in regions subject to short bursts of high intensity precipitation through the process of Horton overland flow. Specifically, if the rainfall intensity exceeds the infiltration capacity of the substrate surface, the excess precipitation would pond and/or flow overland across the green roof surface and into a roof drain. Additional research on the importance of substrate depth would ideally measure (and not model) runoff from side-by-side, full-scale green roofs constructed with different substrate depths over the course of both warm and cold seasons. This study would ideally be reproduced in regions with different precipitation characteristics in an effort to disentangle the various independent variables.
In this study, climate change simulations (assuming warm/wet conditions) resulted in significant, (~5%), decreases in green roof storm water retention, regardless of substrate depth. Future mean event retention was reduced to between 58% and 86% (depending on substrate depth). In the Northeast U.S., climate change is expected to cause both increased precipitation intensity and higher temperatures, trends that were evident in the SWMM-CAT files. SWMM’s predicted net decrease in green roof retention under these conditions suggests that, even though higher temperatures might increase evapotranspiration, more rapidly removing water from the green roof substrate in between rain events, the more intense precipitation will generate more runoff (e.g., Horton overland flow) due to the substrate’s inability to infiltrate it before it runs off the green roof surface. Additional research would test this hypothesis by measuring the green roof substrate’s infiltration capacity, for example with a ring infiltrometer, and comparing this value to the intensity of future climate-changed precipitation in this region.

5. Conclusions

This study investigated the rainfall-runoff response of the Javits Center green roof, currently the second largest green roof in the U.S. A novel monitoring system, installed during construction of the green roof, generated four years of fine temporal resolution discharge measurements from three 186 m2 rooftop tributary areas. Because the monitoring system was designed and installed in parallel to the green roof, we accurately delineated the three tributary areas using aluminum channels that were affixed to the structural roof surface. Pressure transducers installed in Parshall flumes yielded continuous measurements of discharge over the course of 118 events over four years. Storm water retention varied with event depth. The greatest retention rates (96% on average) were observed for light rainfall events (<6.35 mm). The average retention for medium rainfall events was 81.2%, whereas only 27% of heavy events were retained on average. These values, as well as the observed 55% of cumulative precipitation during the monitoring period, were consistent with the findings of other researchers working in other geographic regions with green roofs with other designs, who reported values between 65% and 85% for extensive green roofs in a review of several publications on water retention for intensive and extensive green roofs [10,14]. No relationships between event performance and Julian date, air temperature, mean event intensity, or mean event duration were observed, though no winter observations were used in this study. Along with the flume discharge measurements, the weather station and weighing lysimeter enabled performing a full water balance on each of the observed events. Between the increased moisture storage measured with the lysimeters and the discharge measured with the flumes, all of the incident precipitation accounted for 80% of all the events. Empirical relationships suggested that the green roof’s ability to store moisture is closely related with the duration of the event’s antecedent dry period.
With calibration, the SWMM model with LID-GR control and simulated AET values was able to predict the total volume and peak rate of discharge for most events to within the acceptable validation envelope in the four-year continuous simulation. The model tended to overestimate the peak rates, but not by much. Overall, this modeling exercise suggests that usage of the SWMM LID-GR control with AET values can yield reasonable estimates of green roof discharge and retention, at least during the non-winter months in a temperate climate. Following the work of Krebs et al. [40], future work further refining the AET values with crop coefficients to account for specific vegetation characteristics is proposed.
The monitoring and modeling activities presented a unique opportunity to explore the ability of green roofs with different substrate depths to offset the storm water detention requirements of buildings with different roof areas in New York City. The SWMM model suggests that in New York City climate, storm water retention increases with increasing substrate depth, but with diminishing returns after about 100 mm. If the depth of the JGR substrate were halved, for example as a result of wind erosion, average event retention could decrease by about 9%.
Curves were developed to estimate the percent reduction in the required detention facility volume that could be achieved by installing green roofs of different depths on buildings of different areas. The Javits Center is not required to detain its storm water, since it predates this particular regulation. However, a new building of equivalent roof area would require a cistern of up to 974 m3. If this facility has an effective depth of 1.6 m, its floor area would be roughly 608 m2. With commercial real estate in the area surrounding the JJCC valued at $1 billion USD per city block (or roughly $50,000 USD/m2 assuming a city block in Manhattan is roughly 22,000 m2), the lost real estate value of the 608 m2 footprint could be $30 million USD. The estimated construction cost of a cistern with this footprint is $2.5 million USD, or more than eighteen times the average cost to construct a green roof that would retain the same quantity of stormwater, estimated at $100–350 USD per square meter. With these general assumptions, the initial investment in construction of the JGR appears easily justified.
The observations, modeling, and analysis of the JGR suggests that for some, building green roofs represent a win-win opportunity. By utilizing the building’s rooftop space to mitigate incident rainfall, building owners can maximize developable floor space inside the building, while minimizing the cost to comply with current stormwater management requirements in New York City. Simultaneously, this same strategy can help prepare the city for climate change, while enhancing the city with the other co-benefits of this recently rediscovered approach to urban storm water management.

Author Contributions

Conceptualization and Methodology, F.A.M. and J.C.; Data Collection, D.F., Y.E., J.C. and F.A.M.; Drainage Design, Y.E. and J.C.; Literature Review, S.W., N.A. and J.C.; Model Validation, N.A., T.C. and F.A.M.; Water Balance Analysis, J.C. and D.F.; Writing—Original Draft Preparation, J.C., N.A.; Writing—Review & Editing, N.A., J.C. and F.A.M.

Funding

This research was funded by the Jacob K. Javits Center, with supplemental funding from the Consortium for Climate Risks in the Urban Northeast.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Dunnett, N.; Kingsbury, N. Planting Green Roofs and Living Walls; Timber Press: Portland, OR, USA, 2004. [Google Scholar]
  2. Osmundson, T. Roof Gardens, History, Design, Construction; Norton & Company: New York, NY, USA, 1999. [Google Scholar]
  3. White, R. Building the Ecological City; Woodhead Publication: Cambridge, UK, 2002. [Google Scholar]
  4. Niezgoda, S.L.; Johnson, P.A. Improving the urban stream restoration effort: Identifying critical form and processes relationships. Environ. Manag. 2005, 35, 579–592. [Google Scholar] [CrossRef] [PubMed]
  5. Sutherland, A.B.; Meyer, J.L.; Gardiner, E.P. Effects of land cover on sediment regime and fish assemblage structure in four southern Appalachian streams. Freshw. Biol. 2002, 47, 1791–1805. [Google Scholar] [CrossRef]
  6. U.S. EPA. EPA’S Report on the Environment; U.S. Environmental Protection Agency: Washington, DC, USA, 2008.
  7. Carter, T.; Jackson, R. Vegetated roofs for stormwater management at multiple spatial scales. Landsc. Urban Plan. 2007, 80, 84–94. [Google Scholar] [CrossRef]
  8. Akbari, H.; Rose, L.S. Urban surfaces and heat island mitigation potentials. J. Hum.-Environ. Syst. 2008, 11, 85–101. [Google Scholar] [CrossRef]
  9. Santamouris, M. Cooling the cities—A review of reflective and green roof mitigation technologies to fight heat island and improve comfort in urban environments. Sol. Energy 2014, 103, 682–703. [Google Scholar] [CrossRef]
  10. Berndtsson, J.C. Green Roof performance towards management of runoff water quantity and quality: A review. Ecol. Eng. 2010, 36, 351–360. [Google Scholar] [CrossRef]
  11. Bengtsson, L.; Grahn, L.; Olsson, J. Hydrological function of a thin extensive green roof in south Sweden. Hydrol. Res. 2005, 36, 259–268. [Google Scholar] [CrossRef]
  12. Rosenzweig, C.; Strzepek, K.M.; Major, D.C.; Iglesias, A.; Yates, D.N.; McCluskey, A.; Hillel, D. Water resources for agriculture in a changing climate: International case studies. Glob. Environ. Chang. A 2004, 14, 345–360. [Google Scholar] [CrossRef]
  13. Stovin, V.; Vesuviano, G.; Kasmin, H. The hydrological performance of a green roof test bed under UK climatic conditions. J. Hydrol. 2012, 414, 148–161. [Google Scholar] [CrossRef]
  14. Mentens, J.; Raes, D.; Hermy, M. Green roof as a tool for solving the rainwater runoff problem in urbanized 21st century? Landsc. Urban Plan. 2006, 77, 217–226. [Google Scholar] [CrossRef]
  15. Stovin, V.; Poe, S.; Berretta, C. A modelling study of long term green roof retention performance. J. Environ. Manag. 2013, 131, 206–215. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  16. Emilsson, T.; Rolf, K. Comparison of establishment methods for extensive green roofs in southern Sweden. Urban For. Urban Green. 2005, 3, 103–111. [Google Scholar] [CrossRef]
  17. Peck, S.W.; Kuhn, M. Design Guidelines for Green Roofs; Ontario Association of Architects & Canada Mortgage and Housing Corporation: Toronto, ON, Canada, 2003. [Google Scholar]
  18. Deville, S.; Menonb, M.; Jia, X.; Reed, G.; Stovin, V. The impact of green roof ageing on substrate characteristics and hydrological performance. J. Hydrol. 2017, 547, 332–344. [Google Scholar] [CrossRef]
  19. Li, Y.; Babcock, R. Green roof hydrology performance and modeling: A review. Water Sci. Technol. 2014, 69, 727–738. [Google Scholar] [CrossRef] [PubMed]
  20. Carson, T.; Keeley, M.; Marasco, D.E.; McGillis, W.; Culligan, P.J. Assessing methods for predicting green roof rainfall capture: A comparison between full-scale observations and four hydrological models. Urban Water J. 2017, 14, 589–603. [Google Scholar] [CrossRef]
  21. Cipolla, S.; Maglionico, M.; Stojkov, I. A long term hydrological modeling of an extensive green roof by means of SWIMM. Ecol. Eng. 2016, 95, 876–887. [Google Scholar] [CrossRef]
  22. Carson, T.B.; Marasco, D.E.; Culligan, P.J.; McGilis, W.R. Hydrological performance of extensive green roofs in New York City: Observations and multi-year modeling of three full-scale systems. Environ. Res. Lett. 2013, 8, 024036. [Google Scholar] [CrossRef]
  23. Morgan, S.; Celik, S.; Relziaff, W. Green Roof Storm-Water Runoff Quantity and Quality. J. Environ. Eng. 2013, 139, 471–478. [Google Scholar] [CrossRef]
  24. Nagase, A.; Dunnett, N. Amount of water from different vegetation types on extensive green roofs: Effects of plant species, diversity and plant structure. Landsc. Urban Plan. 2012, 104, 356–363. [Google Scholar] [CrossRef]
  25. Villarreal, E.L.; Bengtsoon, L. Response of a sedum green roof to individual rain event. Ecol. Eng. 2005, 25, 1–7. [Google Scholar] [CrossRef]
  26. Elliott, A.H.; Trowsdale, S.A. A review of models for low impact urban stormwater drainage. Environ. Model. Softw. 2007, 22, 394–405. [Google Scholar] [CrossRef]
  27. Alfredo, K.; Montalto, F.; Goldstein, A. Observed and Modeled Performances of Prototype Green Roof Test Plots Subjected to Simulated Low- and High-Intensity Precipitations in a Laboratory Experiment. J. Hydrol. Eng. 2010, 15, 444–457. [Google Scholar] [CrossRef]
  28. Fassman-Beck, E.; Hunt, W.; Berghage, R.; Carpenter, D.; Kurtz, T.; Stovin, V.; Wadzuk, B. Curve number and runoff coefficients for extensive living roofs. J. Hydrol. Eng. 2016, 21, 04015073. [Google Scholar] [CrossRef]
  29. Getter, K.L.; Rowe, D.B.; Andresen, J.A. Quantifying the effect of slope on extensive green roof stormwater retention. Ecol. Eng. 2007, 31, 225–231. [Google Scholar] [CrossRef]
  30. Hill, J.; Drake, J.; Sleep, B.; Margolis, L. Influences of four extensive green roof design variables on stormwater hydrology. J. Hydrol. Eng. 2017, 22, 04017019. [Google Scholar] [CrossRef]
  31. Hilten, R.N.; Lawrence, T.M.; Tollner, E.W. Modeling stormwater runoff from green roofs with HYDRUS-1D. J. Hydrol. 2008, 358, 288–293. [Google Scholar] [CrossRef]
  32. Carbone, M.; Brunetti, G.; Piro, P. Modelling the Hydraulic behavior of growing media with the explicit finite volume solution. Water 2015, 7, 568–591. [Google Scholar] [CrossRef]
  33. Soulis, K.X.; Valiantzas, J.D.; Ntoulas, N.; Kargas, G.; Nektarios, P.A. Simulation of green roof runoff under different substrate depths and vegetation covers by coupling a simple conceptual and a physically baed hydrological model. J. Environ. Manag. 2017, 200, 434–445. [Google Scholar] [CrossRef] [PubMed]
  34. Yang, Y.; May, T.F. Rapid Assessment of Hydrologic Performance of Low Impact Development Practices under Design Storms. J. Am. Water Resour. Assoc. 2018, 54, 613–630. [Google Scholar] [CrossRef]
  35. Mora-Melia, D.; Lopez-Aburto, C.S.; Ballestereros-Perez, P.; Muñoz-Velasco, P. Viability of Green Roofs as a Flood Mitigation Element in the Central Region of Chile. Sustainability 2018, 20, 1130. [Google Scholar] [CrossRef]
  36. Li, J.; Deng, C.; Li, Y.; Li, Y.; Song, J. Comprehensive Benefit Evaluation System for Low-Impact Development of Urban Stormwater Management Measures. Water Resour. Manag. 2017, 31, 4745–4758. [Google Scholar] [CrossRef]
  37. Urrestarazu Vincent, S.; Radhakrishnan, M.; Hayde, L.; Pathirana, A. Enhancing the Economic Value of Large Investments in Sustainable Drainage Systems (SuDS) through Inclusion of Ecosystems Services Benefits. Water 2017, 9, 841. [Google Scholar] [CrossRef]
  38. Burszta-Adamiak, E.; Mrowiec, M. Modelling of green roofs’ hydrologic performance using EPA’s SWMM. Water Sci. Technol. 2013, 61, 36–42. [Google Scholar] [CrossRef] [PubMed]
  39. Palla, A.; Gnecco, I. Hydrologic modeling of low impact development systems at the urban catchment scale. J. Hydrol. 2015, 528, 361–368. [Google Scholar] [CrossRef]
  40. Krebs, G.; Kuoppamaki, K.; Kokkonem, T.; Koivusalo, H. Simulation of green roof test bed runoff. Hydrol. Process. 2016, 30, 250–262. [Google Scholar] [CrossRef]
  41. NYDEP. Green Infrastructure Plan: Annual Report; New York City Department of Environmental Protection Police: New York, NY, USA, 2016.
  42. Alvizuri, J.; Cataldo, J.; Smalls-Mantey, L.A.; Montalto, F.A. Green roof thermal buffering: Insights derived from fixed and portable monitoring equipment. Energy Build. 2017, 151, 455–468. [Google Scholar] [CrossRef]
  43. Yu, Z.; Miller, S.; Montalto, F.A.; Lall, U. The bridge between precipitation and temperature—Pressure change events: Modeling future non-stationary precipitation. J. Hydrol. 2018, 562, 346–357. [Google Scholar] [CrossRef]
  44. ASCE. Hydrology Handbook, 2nd ed.; American Society of Civil Engineers: New York, NY, USA, 2000. [Google Scholar]
  45. Rossman, L.; Huber, W. Storm Water Management Model Reference Manual Volume I, Hydrology; EPA/600/R-15/162A; US EPA Office of Research and Development: Washington, DC, USA, 2015.
  46. McCuen, R.H.; Johnson, P.A.; Ragan, R.M. Hydrology. FHWA-SA-96-067; Federal Highway Administration: Washington, DC, USA, 1996. [Google Scholar]
  47. American Society of Civil Engineers. Design & Construction of Urban Stormwater Management Systems; ASCE: New York, NY, USA, 1992. [Google Scholar]
  48. Thornthwaite, C.W.; Mather, J.R. Instructions and tables for computing potential evapotranspiration and the water balance. Publ. Climatol. 1957, 10, 185–311. [Google Scholar]
  49. R Development Core Team. R: A Language and Environment Forstatistical Computing; R Foundation for Statistical Computing: Vienna, Austria, 2008; ISBN 3-900051-07-0. Available online: http://www.R-project.org (accessed on 15 January 2018).
  50. Wastewater Planning Users Group (WaPUG). Code of Practice for the Hydraulic Modeling of Sewer Systems, Version 3.001; WaPUG: London, UK, 2002. [Google Scholar]
  51. Rawls, W.J.; Brakensiek, D.L.; Miller, N. Green-ampt infiltration parameters from soils data. J. Hydraul. Eng. 1983, 109, 1316. [Google Scholar] [CrossRef]
  52. Smalls-Mantey, L. The Potential Role of Green Infrastructure in the Mitigation of the Urban Heat Island. Ph.D. Thesis, Drexel University, Philadelphia, PA, USA, 2017. Available online: http://hdl.handle.net/1860/idea:7596 (accessed on 22 October 2018).
  53. International Business Machines Corporation (IBM Corp.). IBM SPSS Statistics for Windows, Version 25.0; IBM Corp.: Armonk, NY, USA, 2017. [Google Scholar]
  54. Rossman, L. SWMM-CAT User’s Guide; EPA/600/R-14/428; US EPA Office of Research and Development: Washington, DC, USA, 2014.
  55. The New York City Department of Environmental Protection Police (NYCDEP). Criteria for Detention Facility Design; Bureau of Water and Sewer Operations: New York, NY, USA, 2012.
  56. Mann, G. Gründächer als lebende Systeme ohne Normzwang: Retentionsverhalten begrünter Dächer in Abhängigkeit der Niederschlagsregion. Dach + Grün 2002, 11, 30–33. [Google Scholar]
  57. DeNardo, J.; Jarrett, A.; Manbeck, H.; Beattle, D.; Berghoge, R. Stormwater mitigation and surface temperature reduction by green roofs. Trans. ASAE 2005, 48, 1491–1496. [Google Scholar] [CrossRef]
  58. Simmons, M.T.; Gardiner, B.; Windhager, S.; Tinsley, J. Green roofs are not created equal: The hydrologic and thermal performance of six different extensive green roofs and reflective and non-reflective roofs in a sub-tropical climate. Urban Ecosyst. 2008, 11, 339–348. [Google Scholar] [CrossRef]
Figure 1. (A) aerial view of the Javits Green Roof (JGR); (B) the three 186 m2 tributary areas on the south roof used for runoff collection; (C) flow pathways in each of the three tributary areas; (D) view of the end of the tributary area (point B in 1C), where it discharges into the flume; and (E) the flume and channel system under construction.
Figure 1. (A) aerial view of the Javits Green Roof (JGR); (B) the three 186 m2 tributary areas on the south roof used for runoff collection; (C) flow pathways in each of the three tributary areas; (D) view of the end of the tributary area (point B in 1C), where it discharges into the flume; and (E) the flume and channel system under construction.
Water 10 01494 g001
Figure 2. (a) Cross section of longitudinal section. (b) Cross section of transverse section.
Figure 2. (a) Cross section of longitudinal section. (b) Cross section of transverse section.
Water 10 01494 g002
Figure 3. Parshall Flume dimensions in cm (not to scale).
Figure 3. Parshall Flume dimensions in cm (not to scale).
Water 10 01494 g003
Figure 4. Runoff hydrograph of event 83 on 27 October 2016.
Figure 4. Runoff hydrograph of event 83 on 27 October 2016.
Water 10 01494 g004
Figure 5. Sensitivity analysis tornado plot of model parameters and their effect on total runoff volume.
Figure 5. Sensitivity analysis tornado plot of model parameters and their effect on total runoff volume.
Water 10 01494 g005
Figure 6. Sensitivity analysis tornado plot of model parameters and their effect on peak runoff rate.
Figure 6. Sensitivity analysis tornado plot of model parameters and their effect on peak runoff rate.
Water 10 01494 g006
Figure 7. Model vs. observed total runoff validation envelope (+20% to −10%).
Figure 7. Model vs. observed total runoff validation envelope (+20% to −10%).
Water 10 01494 g007
Figure 8. Model vs. observed peak runoff validation envelope (+25% to −15%).
Figure 8. Model vs. observed peak runoff validation envelope (+25% to −15%).
Water 10 01494 g008
Figure 9. Water balance analysis (2014–2017).
Figure 9. Water balance analysis (2014–2017).
Water 10 01494 g009
Figure 10. Relationship between lysimeter mass and antecedent dry period.
Figure 10. Relationship between lysimeter mass and antecedent dry period.
Water 10 01494 g010
Figure 11. Modeled green roof performance under varying design depths (12.7–127 mm).
Figure 11. Modeled green roof performance under varying design depths (12.7–127 mm).
Water 10 01494 g011
Figure 12. Variation in runoff coefficient (C) for varying green roof depths and areas. Boxes represent 25th and 75th percentiles, whiskers represent 5th and 95th percentiles.
Figure 12. Variation in runoff coefficient (C) for varying green roof depths and areas. Boxes represent 25th and 75th percentiles, whiskers represent 5th and 95th percentiles.
Water 10 01494 g012
Figure 13. Required detention facility volume to replace different sized green roofs.
Figure 13. Required detention facility volume to replace different sized green roofs.
Water 10 01494 g013
Table 1. Weather station equipment and specifications.
Table 1. Weather station equipment and specifications.
EquipmentManufacturerParameter MeasuredSpecificationsActive SinceRecording Interval
Young Wind Sentry Anemometer-Model 3002Campbell Scientific, Inc. (Logan, UT, USA)Wind speed, directionSpeed: ±0.5 m/s, Direction: ±5°August 20135 min
Series Rainfall Sensor-Model 525Texas Electronics, Inc. (Dallas, TX, USA)Total Rainfall1.0% at 10 mm/h or lessAugust 20135 min
Hukesflux 4 component net radiation sensor-Model NR01Campbell Scientific, Inc. (Logan, UT, USA)Incoming Solar/Shortwave, Incoming Longwave, Outgoing Solar/Shortwave, Outgoing Longwave±10%August 2013Average 5 min
Temperature and Relative Humidity Sensor-Model CS215Campbell Scientific, Inc. (Logan, UT, USA)Air temperature, Relative Humidity±0.3 °C, ±4%August 2013Average 5 min
Pressure TransducerCampbell Scientific, Inc. (Logan, UT, USA)Water level±0.3 °C, full scale rangeAugust 2013Average 5 min
Table 2. Subcatchment properties used in the Environmental Protection Agency (EPA) Stormwater Management Model.
Table 2. Subcatchment properties used in the Environmental Protection Agency (EPA) Stormwater Management Model.
PropertyDescriptionInitial ValueCalibrated ValueSource for Initial Value Selection
Subcatchment Area (m2)The area of the subcatchment drainage in acres-186
Width (m)Average width of overland flow path1213.62
% SlopeAverage surface slop (%)1.01.5
% ImperviousPercent of impervious area-100
N-ImperviousMannings N for impervious area0.0110.013[46]
N-PerviousMannings N for pervious area0.0300.035[46]
Dstore-Impervious (mm)Depth of depression storage on impervious area0.11.27[47]
Dstore-Pervious (mm)Depth of depression storage on pervious area0.12.17[47]
Table 3. Green roof Low Impact Development (LID) control properties defined in SWWM.
Table 3. Green roof Low Impact Development (LID) control properties defined in SWWM.
LayerPropertyInitial ValueCalibrated ValueSource for Initial Value Selection
Surface LayerArea of unit (m2)-186-
% of subcatchment occupied-100-
Surface width (m)1213.62-
% of impervious area treated-100-
Vegetation volume fraction-0.9-
Surface roughness (Manning’s n)0.0300.035[46]
Surface slope (%)1.01.5-
Soil LayerThickness (mm)-25.4-
Porosity-0.4139-
Field Capacity-0.371-
Wilting Point-0.12-
Conductivity (mm/hr)-381[51]
Drainage LayerThickness (mm)-76.2-
Void Fraction-0.75-
Roughness (Manning’s n)-0.05[46]
Table 4. Observed annual Javits Green Roof (JGR) performance.
Table 4. Observed annual Javits Green Roof (JGR) performance.
YearNo. EventsTotal Precipitation (mm)Precipitation Volume (m3)Total Observed Runoff (m3)Average Event % RetentionStandard DeviationCoefficient of Skewness
201431317.459.414.276.0%15.5−1.6
201532268.249.912.375.4%4.9−2.9
201624290.153.911.478.9%10.2−2.7
201731293.854.711.379.3%11.3−1.6
Table 5. Observed JGR storm retention.
Table 5. Observed JGR storm retention.
Type of EventNo. EventsPrecipitation Volume (m3)Total Observed Runoff (m3)Average Event % RetentionSDCoefficient of Skewness
Light (<6.35 mm)5332.01.196.4%0.1−2.0
Medium (6.35–12.7 mm)2742.78.081.2%0.02−0.4
Heavy (>12.7 mm)38167.0122.627.0%0.020.6

Share and Cite

MDPI and ACS Style

Abualfaraj, N.; Cataldo, J.; Elborolosy, Y.; Fagan, D.; Woerdeman, S.; Carson, T.; Montalto, F.A. Monitoring and Modeling the Long-Term Rainfall-Runoff Response of the Jacob K. Javits Center Green Roof. Water 2018, 10, 1494. https://doi.org/10.3390/w10111494

AMA Style

Abualfaraj N, Cataldo J, Elborolosy Y, Fagan D, Woerdeman S, Carson T, Montalto FA. Monitoring and Modeling the Long-Term Rainfall-Runoff Response of the Jacob K. Javits Center Green Roof. Water. 2018; 10(11):1494. https://doi.org/10.3390/w10111494

Chicago/Turabian Style

Abualfaraj, Noura, Joseph Cataldo, Yara Elborolosy, Daniel Fagan, Sloane Woerdeman, Tyler Carson, and Franco A. Montalto. 2018. "Monitoring and Modeling the Long-Term Rainfall-Runoff Response of the Jacob K. Javits Center Green Roof" Water 10, no. 11: 1494. https://doi.org/10.3390/w10111494

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop