segunda-feira, 21 de setembro de 2009

Article Review: CEDAR CREEK

Debele, B., Srinivasan, R., and Parlange, J. Y. (2008). "Coupling upland watershed and downstream waterbody hydrodynamic and water quality models (SWAT and CE-QUAL-W2) for better water resources management in complex river basins." Environmental Modeling and Assessment, 13(1), 135-153.

The article presents a modeling framework of water quantity and quality of upland watershed and a downstream reservoir. Two models have been coupled: SWAT (Soil & Water Assessment Tool) and CE-QUAL-W2. SWAT is a hydrological and hydraulic 1-D hydrodynamic and water quality model used to simulate upland watershed and narrower streams. CE-QUAL-W2 is a waterbody laterally averaged (2-D) hydrodynamic model, suited for narrow and deep waterbodies where lateral variation is minimum. This modeling approach has been tested in Cedar Creek watershed, Texas. The main land uses is agriculture, forest and urban areas. The watershed is drained by two main creeks (Cedar and Kings) and seven smaller tributaries which feed into the Cedar Creek Reservoir. The mean depth of the reservoir is 6.5 meters, a surface area of 34,300 acres and a volume of 566,000 acre-feet. About 25 % of total reservoir inflow is lost through evaporation and approximately 10 % is withdrawn for urban consumption (Fort Worth and nearby towns). Data from a period of 5 years, from 1997 to 2001, was used to calibrate the models. The calibration of the reservoir hydrodynamic showed good waterbalance match between the observed data and simulated results. The water quality calibration was performed with first temperature followed by Dissolved Oxygen. A second dataset of water quality parameters such as organic nitrogen, organic phosphorus, phosphate, chlorophyll a, total nitrogen, and total phosphorus concentration, were available and used in the calibration. Comparison among observed and simulated water quality parameters showed reasonable agreements. It was concluded that this framework can properly simulate seasonal and intraseasonal distributions of many water quality variables, with exception of Chl-a, TP, TN and NH3/NH4. The reasons for that could be by low quality of these input dataset and error propagation and compounding. In conclusion, the authors state that the integration of such models of spatially variable upland watershed hydrological and water quality models with downstream waterbody hydrodynamic and water quality models can help the water resources management, allowing the backtrack of the spatial sources of contaminants and support a good decision making for water quality standards in downstream waterbody.

sexta-feira, 28 de agosto de 2009

Summer Summary

This post is a summary of the main ideas collected during the research period of summer. In total, I read completely 16 articles, and another number of incompletely readings. The research focused mainly in two topics: my present research, which is stormwater sustainability, and optimization techniques, mainly focused in multi-objective genetic algorithm.

For the first and more important topic of my research in this summer, the readings helped me to support the main ideas of the hydrologic footprint residence. HFR is a new metric used to assess the impacts of urbanization in the hydrologic flow regime. Many articles served as references to characterize the impact of urbanization in the hydrological regime. The transformation of natural land covers as forest and grassland to impervious covers alters most of the hydrological processes such as interception, evapotranspiration, infiltration and overland runoff.

In the revised literature, it was identified the main hydrologic metrics used to assess the impact of urbanization. Most of the metrics are based in time series of flow discharges. The main metrics found are: flow frequency curve, flow duration curve, TQmean, which is fraction of year that flow rate exceeds the mean flow rate of the entire flow time series, TQ0.5yr, which is the fraction of year that the flow rate exceeds the historical half-year peak flow rate, peak flow, storm flow duration, time-to-peak, ratio of time-to-peak to storm flow event duration, event-mean rise rate, event-mean fall rate, duration that the flow rate exceeds 75% of the peak flow rate in an event, ratio of event-peak-duration to total event duration, and IHA (indices of hydrological alteration). Those metrics can be used to evaluate different alternatives of development, Best Management Practices (BMPs) and Low Impact Development (LID).

Besides that, an interesting discussion about the effectiveness of BMPs was addressed in some articles. The argument is BMPs introduces a sustained flow for most of the storms which can harm downstream communities. BMPs are design to attenuate bigger storms. Besides that, urbanization increases the frequency of flow events. The combination of more frequent events with sustained flow is being questioned by recent researches.

Another topics explored during this research are the links between the hydrological flow regime and heath of instream ecosystems. Poff et al (1997) summarizes the knowledge about formation of habitats and the natural hydrological regime. Scientific community agrees that stream is a master variable for instream environments because it regulates most of biological and physicochemical processes in rivers. The natural variability of low and high flows plays a crucial role in the entire life cycle of instream and riverine species. Watershed management has been limited to water quality and one aspect of water quantity, which is minimum flow.

The second set of articles is related to optimization algorithms. Two interesting articles present a new heuristic optimization technique called Harmony Search (HS). HS is inspired the way musician improvise harmony searching for maximize the aesthetic estimation of music (objective function). The decision variables are the pitches (sounds) musicians play. Several problems are tested in the articles, and HS presented better performance when compared to other optimization techniques such as Generalized Reduced Gradient (GRC), Evolutionary Programming (EP), Genetic Algorithm (GA), and Nonlinear Programming Gradient (NLPG). This better efficiency, however, should be looked with attention, since performance of methods can be highly dependent of parameters and the problems itself.

The last group of articles is focused in Multi-objective optimization techniques, more specifically genetic algorithms. The objective here is to understand the main operators used in different techniques. Two methods have been studied in more details: MOGA and NSGA-II. One article combined the use of MOGA and Constrained Differential Dynamic Programming in a surface reservoir allocation problem.

In summary, important goals have been achieved during this summer period. The idea of reading two papers a week, one specific of my research and other more general, to discuss in the research group is really efficient. My reading and writing skills have improved a lot in period. I hope we can continue to do that with, hopefully, more students contributing to the discussion.

End of Summer I - List of articles

The summer is over... here is the list of articles that I read in this period:

1. Bertrand-Krajewski, J.-L., Barraud, S., and Chocat, B. (2000). "Need for improved methodologies and measurements for sustainable management of urban water systems." Environmental Impact Assessment Review, 20(3), 323-331.
2. Deb, K., Pratap, A., Agarwal, S., and Meyarivan, T. (2002). "A fast and elitist multiobjective genetic algorithm: NSGA-II." Evolutionary Computation, IEEE Transactions on, 6(2), 182-197.
3. Derek B. Booth, James R. Karr, Sally Schauman, Christopher P. Konrad, Sarah A. Morley, Marit G. Larson, and Stephen J. Burges. (2004). "Reviving Urban Streams: Land Use, Hydrology, Biology, and Human Behavior." Journal of the American Water Resources Association, 40(5), 1351-1364.
4. Egderly, J. L., Roesner, L. A., Rohrer, C. A., and Gironas, J. A. "Quantifying Urban-Induced Flow Regime Alteration and Evaluating Mitigation Alternatives Using Mathematical Models and Hydrologic Metrics." 425.
5. Fan, C., and Li, J. (2004). "A Modelling Analysis of Urban Stormwater Flow Regimes and their Implication for Stream Erosion." Water Quality Research Journal of Canada, 39(4), 356-61.
6. Geem, Z. W. (2009). "Harmony search optimisation to the pump-included water distribution network design." Civil Engineering and Environmental Systems, 26(3), 211 - 221.
7. Getter, K. L., Rowe, D. B., and Andresen, J. A. (2007). "Quantifying the effect of slope on extensive green roof stormwater retention." Ecological Engineering, 31(4), 225-231.
8. Graham, P., Maclean, L., and Medina, D. (2004). "The Role of Water Balance Modelling in the Transition to Low Impact Development." Water Quality Research Journal of Canada, 39(4), 331-42.
9. Konak, A., Coit, D. W., and Smith, A. E. (2006). "Multi-objective optimization using genetic algorithms: A tutorial." Reliability Engineering & System Safety, 91(9), 992-1007.
10. McCuen, R. H. (1979). "Downstream Effects of Stormwater Management Basins." Journal of the Hydraulics Division, 105(11), 1343-1356.
11. Nehrke, S. M., and Roesner, L. A. (2004). "Effects of Design Practice for Flood Control and Best Management Practices on the Flow-Frequency Curve." Journal of Water Resources Planning and Management, 130(2), 131-139.
12. Poff, N. L., Allan, J. D., Bain, M. B., Karr, J. R., Prestegaard, K. L., Richter, B. D., Sparks, R. E., and Stromberg, J. C. (1997). "The Natural Flow Regime." BioScience, 47(11), 769-784.
13. Reddy, M., and Kumar, D. (2006). "Optimal Reservoir Operation Using Multi-Objective Evolutionary Algorithm." Water Resources Management, 20(6), 861-878.
14. Roesner, L. A., Bledsoe, B. P., and Brashear, R. W. (2001). "Are Best-Management-Practice Criteria Really Environmentally Friendly?" Journal of Water Resources Planning and Management, 127(3), 150-154.
15. Yang, C.-C., Chang, L.-C., Yeh, C.-H., and Chen, C.-S. (2007). "Multiobjective Planning of Surface Water Resources by Multiobjective Genetic Algorithm with Constrained Differential Dynamic Programming." Journal of Water Resources Planning and Management, 133(6), 499-508.
16. Zong Woo Geem, Joong Hoon Kim, and Loganathan, G. V. (2001). "A New Heuristic Optimization Algorithm: Harmony Search." SIMULATION, 76(2), 60-68.

segunda-feira, 27 de julho de 2009

Harmony Search

Article Review

Geem, Z.W., Kim, J.H., Loganathan, G.V. (2001) “A New Heuristic Optimization Algorithm: Harmony Search”. 76(2): 60-68.

This paper presents the algorithm Harmony Search (HS), a new heuristic optimization algorithm inspired in the way musician improvise harmonies. The analogy with music is: musical performance seeks a harmony with maximum aesthetic estimation, as the optimization algorithm seeks the best solution. The aesthetic estimation represents the objective function and the pitches of instruments represent decision variables. The decision variables achieve better solution by a process of practice. After each iteration, the solution is evaluated. If the harmony is better than the previous ones, it is included in the Harmony Memory. The process is repeated until the stopping criteria is satisfied.

The harmony memory is initialized with random values. A new harmony is generated based on the solutions already present in the harmony search. If the new harmony is better than the ones in the harmony memory, it is included on it, and the worst harmony is excluded. In order to increase the decision space and find the global optimum, a parameter called Harmony Memory Considering Rate (HMCR), ranging from 0 to 1, is used to randomly generate new solutions. HMCR is uniformly generated. If it is above the current value of the HMCR, the algorithm selects a new variable with certain probability and introduces the new random value.
Another operator is used to escape local optima, similar to crossover in the genetic algorithm (GA). The Pitch Adjusting Rate (PAR) is used to select a variable with certain probability and shifts the value with its neighbors. The article uses the example of the solution {1,3,4,6,7,9}. The value {6} can be moved to neighboring {4} or {7}.

The article illustrates the harmony search with 3 problems found in the literature. The first one is the Traveling Salesman Problem (TSP), which tries to find the shortest path for a salesman that has to visit a certain number of cities once. The second problem is a 2 dimension constrained minimization problem. The performance of HS is compared with others methods like GA, Generalized Reduced Gradient (GRC) and Evolutionary Programming (EP). HS had better performance compared to the other methods. The last problem is the design of a pipeline network for water supply consisted of 32 nodes, 34 links and 3 loops. The results were compared with GA and Nonlinear Programming Gradient (NLPG). The HS produced the solution with the least cost related to the other methods.

Discussion

The HS algorithm is a very interesting methodology. Particularly for me that play the guitar, was very curious to see how real life phenomena can be invently used to develop new techniques to solve older and new problems. This paper was also interesting because I had the opportunity to attend the talk of one of the authors in the 2009 US-Korea Conference on Science, Technology and Entrepreneurship, last July 16-19, 2009 in Raleigh, North Carolina.
During our weekly discussion meeting, we tried to compared the differences and similarities between HS and GA. As a heuristic method, HS is similar to GA. One aspect that seems important is that HS creates a new solution based in all existing solutions presented in the harmony memory, while GA creates a new solution only with two previous solutions (parents).

sexta-feira, 26 de junho de 2009

Downstream Effects of Stormwater Management Basins

Article Review

McCuen, R.H. (1979) “Downstream effects of Stormwater Management Basins” Journal of the Hydraulics Division, Vol. 105, No. 11.

Prof. McCuen addresses the problem of hydrological alteration due urbanization. The use of impervious materials decreases natural storage, increasing peak rates, decreasing time-to-peak, and increasing the volumes of runoff and sediments. To mitigate these detrimental effects storm water management (SWM) is required. However, according to the paper, many SWM strategies fail in mitigate the urbanization effects because “they only attempt to limit peak flow rates and sediments volumes” to predevelopment conditions. This is in general attempted by the use of detention basins. This strategy, however, is not sufficient. This problem lead to main question of the article: “Are SWM basins effectively meeting the intent of stormwater management?”

“Urban development causes a decrease in natural storage, including infiltration, interception, and depression storage”. SWM basin is an attempt to restore the loss of natural storage in the watershed. This volume can be artificially restored, but the timing characteristics introduced by detention basins are significantly different from the predevelopment conditions. In natural conditions, the major part of the storage volume is filled during the storm event, while in a detention facility this storage is rapidly filled during the risen portion of the hydrograph.

A second effect introduced by detention basin is that the increased volume of direct runoff is released at a higher rate when compared to natural storage. This cause increase in the duration of bankfull flows.

Due the reasons explained so far, the work examines the hypothesis “that, when changes in both timing characteristics and flow volumes are not corrected by SWM, they may have an adverse effect on flow conditions downstream from the SWM basin.”

A 2.12 square mile watershed was simulated using the Soil Conservation Service TR-20 hydrological model. The predominant soils are silt loams, characterized as deep, well-drained, and moderately erodible. Simulations indicate that peak flow increases from 11.7 to 59.3 cfs (factor of five!!!) from pre-development to urbanized condition of one catchment for the 10 yr-storm. In terms of volume, there is an increase of 100%. The time-to-peak decrease 7 minutes. Introducing a detention basin, the peak flow is reduced to 13.08 cfs, but the time-to-peak is 36 minutes greater than that which existed prior to development. A compared analysis is done for the 100 yr-storm. The effects are the same, but the differences are much lower.

Looking further downstream in the watershed, the computation indicates that for the 100 yr-storm the SWM basin increases the peak flow compared to the post development with no control. This is true for approximately 5 miles downstream the storage facility. The reason for that is changes in timing characteristics affects the runoff conditions downstream.

After the hydrological analysis, the author has approached the bedload transport. Although its is recognized the difficulty in estimating bedload discharges, a theoretical model can be used to assess different scenarios and trace comparisons between them. A model developed by Goncharov was used to predict the sediment discharge as a function of mean flow velocity, mean sediment diameter, the critical velocity of sediment, and the mean depth. The previous discussed scenarios were evaluated to generate estimative of sediment yield. The results indicate that even with the SWM basin there is a significant increase in total volume of sediment. This is due the sustain flow introduced by the basin, which induce higher velocities than the critical velocity over a longer period.

The author finalizes the article with some considerations related with SWM policies and designs. He considers three aspects important:

(1) “The recognition of a comprehensive interpretation to the definition of the intent of SWM;

(2) A policy statement that reflects this interpretation and provides the guidelines that are required to institute proper design methods;

(3) Design methods and construction practices that are representative of the state-of-art and the fully consider hydrologic changes that go along with the urbanization process.”

He suggests that alternatives like grass-lined swales, rooftops detention, and porous pavements are SWM measures that better represents the natural storage found in predevelopment conditions, both in volume and temporal characteristics. A second alternative is to make change in the basin designs, although that may cause conflicts between flood control and pollution/sediment control, for example.

quinta-feira, 25 de junho de 2009

Flow duration curves

Article Review

Fan, C., Li, James (2004) “A Modelling Analysis of Urban Stormwater Flow Regimes and their Implication for Stream Erosion” Water Quality Re. Journal Canada, Vol. 39, No. 4.

The objective of this paper is to analyze flow-duration characteristics of a stream, with and without control measures applied during and after urbanization and its potential for erosion.

“Flow duration curve is the complement of the cumulative distribution function of stream flow. It represents the relationship between the magnitude and frequency of daily, weekly, monthly or even longer period of time of stream flow for a particular watershed.” The flow duration curve is the relationship between exccedance probability and discharge.

According to the authors, stream erosion is a continuous process of removal and deposition of sediments. Urbanization can alters drastically this balance due alteration of the flow regime. Since measuring both discharge of solid and liquid is something difficult to perform, especially in the planning process of development, the hypothesis of this work is that hydrological models can be used to estimate the flow regimes alteration by the flow-duration curve and evaluate the implication on stream erosion.

A period of 3 months of measure data were used to calibrate the model ISWMS (Integrated Storm Water Management Simulator). The watershed was discretized in 9 sub catchments. The calibrated parameters were base flow and horton’s equation for infiltration. Twenty-seven years of hourly rainfall data was used in the long-term simulations.

Different degrees of development were tested by changing the percent of imperviousness (20, 30 and 40 %), of all catchment, with exception of catchment #5 which represents a wetland area. It was noted that urbanized duration curves are shifted from the original curve as more imperviousness is used, indicating that for the same discharge, the frequency is higher.

To analyze the impact of flow control measures, three scenarios were tested: (1) an extended detention basin was added to each subcatchment except subcatchment 5 (the wetland); (2) two extended detention basins were added, one more upstream and other downstream; and (3) an extended detention basin was added downstream of all the subcatchments. The basins were sized considering the criteria of capturing 5 mm of rainfall and the outlet orifices were sized for 24 and 48 hours.

The first analysis compared the use of detention basin in the outlet of the watershed (scenario 3). The results indicate that both 24 hrs and 48 hrs extended basins control the runoff greater than 0.1 m³/s. For lower flow than 0.05 m³/s, 24 hrs basins performed slightly better than 48 hrs basins, and the opposite for flows varying from 0.1 to 0.05 m³/s.

When other locations are analyzed, it was found that regional control (scenario 2) provide better control performance than scenario 3 for flows greater than 0.07 m³/s. Local control (scenario 1) generally performs best when the flow is less than 0.09 m³/s. The authors conclude that the optimal location of detention basis is dependent on the flow regime.

In summary, the paper concludes that urbanization can change significantly the hydrological regime of a watershed and that extended detention basins with 48 hours of detention time provide better control of flow-duration characteristics than 24 hours of detention for higher flows. For lower flows, the 24 hour basin perform better. This make the authors suggest that the detention basins should be sized based on the range of flows which determine the stability of a stream. Finally, the paper states that if a watershed is urbanized beyond 40 % of imperviousness, it may be impossible to restore the flow-duration curve to its original position.

Discussion

This paper seems valid in order to illustrate the impact of urbanization in the flow regime. It has used as a tool of analysis the flow duration curve, analyses that can be found in more literature (see last post). Some comments however I would like to do:

(1) The paper doesn’t specify the area of the watershed that is been analyzed. The findings of the work may not be extrapolate to other watershed with different scales;
(2) The period of calibration is very small (just 3 months). A period of verification of the modeling was not performed;
(3) The limit of 40 % of imperviousness indicating that it may be impossible to restore the flow-duration curve seems very specific for this particular watershed. Different hydrological characteristics may turn this number lower or even higher;
(4) The results shown in the work prove good control over the high flow portion of the flow duration curve. For lower flows (and high frequent event) none of the scenarios seems to replicate the predevelopment conditions (see last post);
(5) The original hypothesis of the work was partially proved. The hydrological alterations were demonstrated, but no further analysis was done to evaluate the implication of the alterations on stream erosion.

quinta-feira, 11 de junho de 2009

Are Best-Management-Practice really environmentally friendly?

Article Review

Roesner, L.A., Bledsoe, B.P., Brashear, R.W. (2001) “Are Best-Management-Practice really environmentally friendly?” Journal of Water Resources Planning and Management, v. 127, no. 3 (May/June).

This article examines the argument that best-management-pratice (BMP) doesn’t not protect urban streams properly. The definition of BMP is strategies to control of urban runoff, usually store structures as detention basins. The authors defend that the BMP is not the problem in itself, but rather the design guidance for BMP outlet flow control that doesn’t take into account the geomorphologic characteristics of the receiver stream.

The authors present a description of the impacts of urbanization in the hydrology of urban watershed. Some numbers illustrate the changes: the peak flow can increase by 2 to more than 10 in terms of magnitude. Also frequency of two-year peak runoff rate in the predevelopment state occurs three to eight times per year after urbanization. The effects of this changes cause severe stream channel erosion and downstream flooding. In terms of water quality, the effects of urban runoff also harm aquatic habitats.

Considering the increase of runoff, many municipalities require big storms to be controlled, in order that the post-development peak flow does not exceed the pre-development discharge. Common recurrence times are 10 and 25 years. Few locals control the 2-year storm. The article states that while detention basins reduce flooding, they are not effective in reducing erosion. One reason is due the sustained flow that detention ponds induce downstream. The outlet flow is prolonged, subjecting the stream channel for a longer period of time. The second reason is the basin only attenuates the flow of the design storm. Attenuation of different storms is coincidental. Considering that urbanization increases the frequency of pre-development peakflow and basins are designed to regulate storms higher than the 10 or 25 year-storms, most storms will pass through the structure unregulated, subjecting the downstream channel to erosive velocities more frequently.

For correcting the problem, the writers believe that “integrating the design of storage and outlet controls for flood control and BMPs it is possible to better control wet weather discharges to receiving waters over the entire range of flows that the facility experiences.” The question that stormwater managers should ask is: “given the potential flow of sediment and water out of the watershed, what are our options for achieving desired social and biological endpoints?”

The article proposes a two step approach to link biological indices to watershed parameters. The first step is to develop a relationship between stream geomorphologic characteristics and biologic indicators. The second step is to develop a relationship between watershed hydrogeometric parameters and stream geomorphologic parameters. It is believed that with this approach is possible to design holistic runoff control facilities into urban development so both flood protection and maintenance of ecosystem can be achieved.