Office of Water Environmental Protection Agency August 2010 E PA 800-R-10-001 Climate Change Vulnerability Assessments: A Review of Water Utility Practices ------- Acknowledgements This report was prepared for EPA by Stratus Consulting Inc., Boulder, CO, under subcontract to Abt Associates, Inc., under EPA Contract EP-C-07-023, Work Assignment 2-13. Karen Metchis served as Technical Project Officer and principal reviewer of this document. EPA wishes to acknowledge the principle authors, Jason M. Vogel and Joel B. Smith for their diligence and expertise in preparing this report. Further thanks go to EPA ORD project officer, Jill Neal, and ORD scientists Jeffrey Yang, Ph.D., and Thomas Johnson, Ph.D., for their comments and suggestions in helping shape this report. This report represents one product resulting from the 'First National Expert Workshop on Water Utilities and Climate Change," held January 20-21, 2009, in Washington, D.C. Acknowledgement also goes to Elizabeth Corr, Robert Cantilli, EPA Office of Water, and Jeff Yang, Ph.D., ORD, for their collaboration in formulating the workshop. Notice The U.S. Environmental Protection Agency, through its Office of Water and Office of Research and Development, funded and managed the development of this report. Any opinions expressed in this report are those of the author(s) and do not necessarily reflect the views of the Agency, therefore, no official endorsement should be inferred. Any mention of trade names or commercial products does not constitute endorsement or recommendation for use. ------- Contents 1. EXECUTIVE SUMMARY 1 2. INTRODUCTION 2 3. APPROACHES TO ASSESSING CLIMATE CHANGE VULNERABILITY 4 3.1 Top-down Modeling Assessments 5 3.1.1 Scenario analyses 6 3.1.2 Sensitivity analyses 6 3.1.3 Paleoclimate or historic analyses 7 3.2 Bottom-up Threshold Analyses 7 3.3 Policy Environment as a Constraint 8 4. SOURCES OF CLIMATE INFORMATION 9 4.1 The Instrumental Record 10 4.2 Paleoclimate Data 10 4.3 Literature Reviews 11 4.4 Climate Projections 12 4.4.1 Climate scenarios 12 4.4.2 Global climate models 14 4.4.3 Downscaling 16 4.4.4 Combining methods: Climate projections with paleoclimate data 18 5. MODELING CHANGES IN WATER RESOURCES 20 5.1 Water Quantity 20 5.2 Water Demand 21 6. SUMMARY 23 7. RECOMMENDATIONS FOR FURTHER STUDY 25 8. REFERENCES 26 ------- 1. Executive Summary Climate change poses a variety of challenges for water management, and there is a need to develop methods for understanding and managing risk. While much has been written about the projected impacts of climate change at the continental or regional scale, scientists are quick to caution decision makers about using projections based on global circulation models (GCMs) for local decision making. This 'uncertainty' about specific impacts on local systems has raised concern about the ability of water resource managers to plan for climatic and hydrological changes at the local scale, and has spurred recent activity to develop methods for understanding vulnerabilities, including how to downscale climate models. This study examines and documents the steps taken by some of the leading utilities in an attempt to identify the emergent characteristics of water utility climate change vulnerability assessments. By examining the approaches taken and articulating the steps, information, and judgments needed for such decision making, we hope to contribute to the collaborative problem solving among the user and research communities who are working to further refine and validate such procedures. The study describes the activities of eight water utilities who have conducted climate vulnerability assessments: East Bay Municipal Utility District, City of Boulder Utilities Division, Denver Water, Massachusetts Water Resources Authority, New York City Department of Environmental Protection, Portland Water Bureau, Lower Colorado River Authority, and Seattle Public Utilities. The following general observations can be made. First, while the intention was to evaluate both drinking water and wastewater utilities in the systems evaluated, most vulnerability assessment work focused on water quantity and water demand and in only one case on water quality. Second, while utility managers typically possess expertise about their systems, their hydrological and management models, and the local hydrology, they have more limited access to climate change information. Utilities are able to examine impacts on their operations due to changes in climate-sensitive variables (e.g., flow, demand), however, in many cases, in order to understand the changes in climate that affect those variables, utilities are engaging outside climate expertise. Third, in each case study presented in this report, utilities are compensating for the 'uncertainty' about climate change by evaluating a range of scenarios or models, and in most cases the resulting information is used to test the robustness of existing decision making or planning against different future, plausible scenarios. Climate change is a complex issue, and more work is needed to establish reliable practices for incorporating climate change into water utility decisions and planning. At the same time, utilities appear to have benefited from their various efforts to understand their potential vulnerabilities and to evaluate long term planning options. Despite the uncertainty, there are reasonable and prudent steps utilities and other water managers can take to better understand and manage climate risk. 1 ------- 2. Introduction For many decades, water utilities have proactively assessed the ability of their systems to provide a reliable supply of drinking water or adequate wastewater services under assumptions of population change, municipal expansion (e.g., annexation), technological innovation, and changes in regulations. These were the foreseeable challenges utilities addressed to ensure that ratepayers received high quality water services at minimal cost. Historically, utilities assumed stationarity of climate in their water resource planning - the idea that natural systems fluctuate within an unchanging envelope of variability. Recently, however, this assumption has been challenged, forcing managers to rethink future water resources planning and management with regard to climate change (Milly et al., 2008). In response, over the last decade, water utilities have expanded their risk assessment efforts to address changing climate conditions. Generally speaking, vulnerability assessments done by water utilities are grounded in a thorough understanding of their water system. On a practical level, this means that vulnerability studies use a variety of tools and models that have been developed and refined to reflect the local hydrology, climate, infrastructure, operations, and demands that are unique to a particular water utility. Consequently, no single tool currently exists that can comprehensively address the vulnerability of diverse water systems to climate change. For example, most utilities have operational or management models that estimate the day-to-day operation of their water or wastewater systems. By necessity, these operational models are either unique to a particular utility (e.g., a proprietary system model) or in some cases, customized versions of a tool used by many utilities (e.g., the Water Evaluation and Planning, or WEAP, system model). The purpose of this report is to illustrate the approaches used by water utilities to assess their vulnerability to climate change. This is a review of best practices in this emerging effort across the industry for the purpose of informing utilities considering engaging in this issue about the various methods used by their peers. The report does not judge or evaluate the efforts of utilities or the merits of different vulnerability assessment methods, but describes the efforts of eight utilities as broad approaches, tools, and methods worthy of consideration by other utilities depending upon their needs, available resources, and other factors. In this report we examined the ways in which several water utilities evaluated their potential vulnerability to climate change. Utilities included in this study were identified from published studies on water resources and climate change, including six Water Research Foundation reports (Strange et al., 2009; WRF, 2010a, 201 Ob, 2010c, 2010d, 201 Of), two reports by the Water Utility Climate Alliance (Barsugli, et al., 2009; Means et al., 2010), and participants in the ongoing Joint Front Range Climate Change Vulnerability Study (WRF, 2010e). We also included utilities selected from the list of participants in two recent workshops focused on water resources and climate change: the National Drinking Water Advisory Council's Climate Ready Water Utilities Working Group and U.S. Environmental Protection Agency's (EPA's) First National Expert and Stakeholder Workshop on Water Infrastructure Sustainability and 2 ------- Adaptation to Climate Change (U.S. EPA, 2009).l In total, we identified 50 water utilities as potential candidates for in-depth case studies. Initial review, however, revealed that many of the selected utilities were interested in tracking the climate change issue and learning what others had done, but had not yet engaged in their own climate change vulnerability assessments. If this was the case, the utility was eliminated from the initial list for further study. We identified five utilities that had published reports on their own climate change vulnerability work, and another five utilities that appeared to have conducted vulnerability work but had not published the modeling results. Because of time and resource limitations, the study was limited to eight of these 10 utilities. The two utilities not included were Inland Empire Utilities Agency and King County Wastewater Treatment Division. The eight utilities included in the assessment included: > East Bay Municipal Utility District (EBMUD) n A water supply and wastewater utility serving 1.3 million customers. > City of Boulder Utilities Division ° A water supply and wastewater utility serving 113,000 customers. > Denver Water n A water supply utility serving 1.3 million customers. > Massachusetts Water Resources Authority (MWRA) ° A water supply and wastewater utility serving 2.2 million customers. > New York City Department of Environmental Protection (NYCDEP) n A water supply and wastewater utility serving 9.2 million customers. > Portland Water Bureau (PWB) n A water supply utility serving 860,000 customers. > Lower Colorado River Authority (LCRA) n A conservation and reclamation district that manages water supply along a 600- mile stretch of the Colorado River in Texas, operates six dams, and helps communities plan and coordinate their water and wastewater needs. > Seattle Public Utilities (SPU) n A water supply utility serving 1.35 million customers. Finally, we studied these eight utilities' climate change vulnerability analyses in depth to identify tools and approaches used by those water utilities in their assessments. The remainder of this report synthesizes the insights from the analysis of these eight utilities in the following sections: approaches to assessing climate change vulnerability, sources of climate information, modeling changes in water resources, summary, and recommendations for further study. 1. A forthcoming report from the Water Research Foundation (WRF, 2010d) offers case study information relevant to this report, but the final draft was not publicly available for inclusion in this report. ------- 3. Approaches to Assessing Climate Change Vulnerability Most water managers engaged in climate vulnerability analyses have a strong technical understanding of their water systems, including local hydrology, historical operating conditions, and standard operational practices.2 However, they typically are not climatologists and have limited experience assessing risks from climate change. As a result, a water utility can either engage outside expertise to do a sophisticated exploration of the implications of climate change or they can make do with the information and capabilities readily available. Utilities that choose not to engage outside experts are generally limited to analyses based on their internal understanding of their water system and/or sensitivity analyses to assess what hypothetical climate changes would mean for their water system. In contrast, utilities that engage outside expertise can generate climate projections specific to their watershed. This more resource- intensive approach enables computationally sophisticated vulnerability analyses. Approaches to assessing climate change risks are generally classified as either "top-down" modeling assessments or "bottom-up" threshold analyses (e.g., Miller and Yates, 2006; Freas et al., 2008; Stratus Consulting and MWH Global, 2009). This categorization is useful, and, generally speaking, utilities choose one of these two approaches to initiate their climate vulnerability assessment. However, it should be noted that these two approaches are not mutually exclusive and utilities often incorporate elements of both, either in series or in parallel. The choice of one approach over another, however, has significant implications for study design and resource needs, making them useful as categorical descriptions. The following utilities initially used a top-down assessment: > City of Boulder Utilities Division > Denver Water > MWRA > NYCDEP > PWB > SPU The following utilities initially used a bottom-up assessment: > EBMUD > LCRA In general, top-down assessments are model and data driven. They often are more time and resource intensive than bottom-up assessments, sometimes requiring expertise beyond the 2. Furthermore, most utilities take land use changes, population projections, and economic development into account in infrastructure design and operational planning. Climate, however, often gets lumped in with hydrological considerations and does not rise to the level of an independent planning consideration. ------- capacity of many utilities. In general, bottom-up assessments are driven by knowledge of the utility system itself. Such assessments often focus on defining critical system thresholds or specific decisions that may be sensitive to climate change. This can often be done qualitatively, and the necessary expertise often resides within the utility. As suggested above, these two approaches are not mutually exclusive. 3.1 Top-down Modeling Assessments Six utilities in our study initiated their vulnerability assessments with a top-down approach: Boulder Utilities (Smith et al., 2009), Denver Water (Laurna Kaatz, personal communication, April 8, 2009), MWRA (Stephen Estes-Smargiassi, personal communication, February 26, 2010), NYCDEP (2008; Major et al., 2007), PWB (Palmer and Hahn, 2002), and SPU (Palmer, 2007). This method projects future climate conditions specific to the watershed of concern and is the most computationally and resource intensive of the analysis approaches used by utilities. Because of this, the top-down modeling assessments done by the utilities examined for this report were often completed by the water utility in conjunction with an academic institution or consulting firm with expertise in climate modeling. For example, SPU and PWB partnered with the University of Washington's Climate Impacts Group (UW-CIG), NYCDEP partnered with the City University of New York Institute for Sustainable Cities (CISC) and Columbia University's Center for Climate Systems Research (CCSR), LCRA used CH2M Hill as a consultant, and the City of Boulder Utilities Division engaged Stratus Consulting and Hydrosphere (now AMEC). Top-down assessments can involve relying on historic hydrology, using paleoclimate records, engaging in literature reviews, or they can be done using climate projections. All these options are discussed in more detail in Section 3 (Sources of Climate Information). By far, the most common form of top-down assessment involves generating downscaled climate projections for a utility's local watershed using global climate model (GCM) output. To do this, the analysis generally includes assumptions about future emission levels (i.e., global socioeconomic development scenarios), the effects of those emission levels on global climate (e.g., temperature and precipitation), the translation of global climate effects to the region or watershed, and the specification of climate changes on hydrology (e.g., streamflow), as well as operational, management, and demand models. All of the utilities we studied in depth assessed a number of climate variables; however, the most common climate variables of interest were those that fed into their watershed-specific hydrology models or system-specific management or operations models. Most often, this included temperature and precipitation data or streamflow data at the spatial and temporal scales of each individual utility's hydrology, management, operations, or other models.3 Therefore, to make meaningful comparisons in this study, we focused our investigation of the different modeling approaches to develop projections of temperature and precipitation. It is worth noting, however, 3. This conclusion was also reached in a study of the Water Utility Climate Alliance utilities (Barsugli et al., 2009). 5 ------- that some utilities also investigated issues of special concern for their water system. For example, NYCDEP was concerned about source water quality, especially due to heavy precipitation events and resulting turbidity, and MWRA was concerned about likely impacts on the safe yield of its sources. Top-down modeling has three important subcategories: scenario analyses, sensitivity analyses, and paleoclimate data or historic climate observations to define temperature and precipitation patterns for water system planning purposes. Each of these subcategories is described in detail in the following sections. 3.1.1 Scenario analyses The scenario approach to a top-down assessment begins by defining at least one, but more commonly two or three plausible greenhouse gas (GHG) emission futures that merit consideration by a water utility. A scenario analysis then works step by step through GCMs, downscaling those results,4 and finally using those results in water system specific hydrology, demand, operational, and/or management models as appropriate to investigate the vulnerability of a water system to each scenario. Section 3.4 (Climate Projections) discusses these individual steps in greater detail. In many cases, scenarios are developed to capture a wide range of plausible futures by incorporating an extensive range of models, emissions scenarios, and projected demands. As an example, the SPU study applied the scenario analysis technique and used GCMs and socioeconomic projections to assess a "middle of the road," a "warmer and wetter," and a "warmer and drier" scenario. MWRA's approach used the outcome of all GCMs to estimate the probability distributions of likely changes in temperature and precipitation and then utilized those distributions in the statistical downscaling. This methodology, developed at NCAR, avoids using a scenario tied to a single GCM. 3.1.2 Sensitivity analyses Another form of the top-down approach uses a sensitivity analysis, incorporating the use of incremental changes in climate such as 1°, 2°, and 3°C (1.8°, 3.6°, and 5.4°F) annual temperature increases combined with +/- 0%, 10%, or 20% annual changes in precipitation.5 PWB used this type of sensitivity analysis to bracket a plausible range of climate-altered future hydrology. But utilities can also engage in a sensitivity analysis, as was done by EBMUD, with a bottom-up 4. The term "downscaling" has come to mean the use of higher resolution regional climate models (RCMs; dynamic downscaling) or complex statistical approaches (statistical downscaling), but in this case can also refer to simple downscaling, which typically involves adding the temperature increase in a GCM grid to observed temperatures for stations in the grid box and multiplying the percentage change in precipitation by the observed precipitation record. 5. Sometimes sensitivity analyses use arbitrary increments, but they can also be informed by the scientific literature. Sensitivity analyses can use different changes in temperature and precipitation by season, but this approach was not used in the vulnerability analyses investigated for this report. ------- approach to avoid the computationally and resource intensive steps of a scenario analysis. A sensitivity analysis does not require using GHG emissions projections, using climate model projections, or downscaling climate model output. 3.1.3 Paleoclimate or historic analyses A third top-down approach incorporates paleoclimate studies or historic climate observations to define temperature and precipitation patterns for water system planning purposes. The method often defines a worst case scenario (e.g., three consecutive years of the drought of record) or examines the water system effects of droughts in the paleorecord outside of observed variability. While the City of Boulder Utilities Division, Denver Water, and EBMUD used paleoclimate studies to complement their climate change vulnerability studies, none of them relied exclusively on paleoclimate or conservative extrapolations from historic data.6 3.2 Bottom-up Threshold Analyses The bottom-up approach to climate change vulnerability assessments is grounded in knowledge of the water system itself. Using this approach, a qualitative system assessment is done to determine which system components are potentially vulnerable to change. The results of the assessment can be used to focus further study on specific impacts of concern. EBMUD, which used this approach, stated, "In a 'Bottom-Up' approach, the most critical vulnerabilities of the District's water supply system are identified, the causes of those vulnerabilities are articulated, and then steps are taken to better address and solve the vulnerability in the face of climatic uncertainty" (EBMUD, 2009, p. 4-16). They noted their concept of a bottom-up approach is adapted from the AwwaRF (now the Water Research Foundation) publication, Climate Change and Water Resources: A Primer for Municipal Water Providers (Miller and Yates, 2006). More specifically, the EBMUD research approach involved a portfolio evaluation that first identified potential portfolio components (e.g., new reservoirs, expanded reservoir storage, increased conservation, conjunctive use, water reclamation, desalination, interbasin transfers), screened those components for technical, environmental, and economic feasibility, and then constructed alternate portfolios of multiple components that could meet projected demands (e.g., increased conservation and conjunctive use, or water reclamation and interbasin transfers). Then, a preliminary portfolio analysis was conducted using a combination of the WEAP system model and the district's EBMUDsim model - known collectively as the "W-E model." Portfolios that performed poorly under current hydrological conditions were eliminated and the remaining portfolios were subjected to detailed analyses under anticipated climate change conditions using 6. The lack of evidence for sole reliance on these types of analyses in the cases investigated for this report is largely a product of the case selection methodology, which focused on utilities with a published record of sophisticated vulnerability work. 7 ------- the W-E model. By screening out portfolios that perform poorly under current conditions, EBMUD implicitly limited options to those that could be implemented under current conditions without significant reduction of reliability, but which presumably would improve the system under climate change relative to the current system. By comparing the climate change scenarios to a baseline scenario, the sensitivity of the EBMUD water system to each component was assessed to identify critical vulnerabilities and identify portfolios that addressed those vulnerabilities. LCRA also used a bottom-up approach, and in CH2M Hill's report for the authority, the threshold approach was defined in this way: "This threshold approach identifies system components that are dependent on the status of climate variables (e.g., precipitation, temperature, etc.) and the overall system risk to climate change, resulting in a preliminary risk assessment based on the professional judgment of experts who know the system and the planning area" (CH2M Hill, 2008, p. 3-2). This approach is a qualitative or semi-quantitative analysis that consisted of (1) identifying the climate variables of importance and exploring the sensitivity of LCRA to these variables; (2) determining water system responses to a range of potential climate changes; (3) assessing the vulnerability of LCRA to climate change impacts; (4) assessing system performance according to the uncertainty associated with climate change factors driving LCRA vulnerability; and (5) evaluating overall system risk and identifying areas in need of further analysis. 3.3 Policy Environment as a Constraint A key constraining factor affecting any utility's ability to assess climate change vulnerability is the level of support of the local community about climate change, in particular utility ratepayers, boards of directors, and locally elected officials. In many circumstances, it is not feasible for a utility to engage in a "climate change" vulnerability study because there is no political support for the effort. Therefore, the policy environment directly impacts a utility's approach to climate change vulnerability assessment. For example, New York City Mayor Michael Bloomberg has taken climate change seriously during his tenure, providing political support for climate change adaptation efforts. As a result, NYCDEP has been able to do sophisticated analyses and make long-term investments to assess potential risks. Many water resource managers in less supportive political environments, however, do not have the political support needed to take on climate change directly or explicitly and must use historic or paleoclimate data rather than downscaled climate model data to justify operational changes or investments that require a rate increase. If support does not exist for climate studies, a water manager might elect to do a sensitivity analysis as a lower-profile way to understand the potential vulnerabilities of their water system to changes in climate conditions - climate change or climate variability.7 7. Water utilities might also select these methods for other reasons, such as a lack of funding or because uncertainty in the projections does not yield enough value added to justify the effort of an intensive climate change vulnerability assessment. ------- Vulnerability to changing environmental conditions is often incorporated into contemporary utility planning efforts even if those efforts are not billed as climate change vulnerability assessments. Indeed, during research for this report, we discovered a utility that reported taking sea level rise into account in their infrastructure design as far back as 1989, but without advertising the analysis as "climate change." In summary, the policy environment in which a utility operates may dictate what approaches to climate vulnerability assessment are acceptable, how such studies are conducted, what data sources are used, and other important aspects of vulnerability assessment design. 4. Sources of Climate Information Regardless of how a climate change vulnerability assessment is framed, climate information that operates at the spatial and temporal scale of a utility's hydrology, planning, and/or operational models is almost always used in vulnerability assessments.8 Most commonly, temperature and precipitation estimates at daily, weekly, or monthly timesteps are input into utility-specific hydrology, demand, operations, and/or management models to assess the effects of different climate conditions on water supply or wastewater services.9 Temperature and precipitation estimates can be derived from a number of sources, including the instrumental record, paleoclimate studies of variability beyond the historical record, literature reviews, and downscaling GCM projections to utility-specific watershed(s) to generate projected future climate conditions. These sources of climate information have different associated logistical and resource constraints and each is explored in more depth below. The key criterion for selecting a source of climate information is that the approach and the data can be justified, sometimes only to a technical audience, but increasingly to ratepayers, boards of directors, and elected officials as well. Utilities need credible climate information in order to consider climate change impacts on their water systems. Historically, most utilities that address climate change issues have relied on partnerships with academic and research institutions to help them select appropriate climate data, however, many utilities are increasingly turning to private consulting companies to help them with this task. Small- and medium-sized utilities, however, often lack the financial resources needed to engage consultants to provide technical guidance on climate change vulnerability studies and their ability to partner with academic institutions is not clear. The scalability of the climate information approaches described below must be carefully considered for potential 8. Note that climate information includes more than GCM output, including historic weather station data, stream gauge data, paleoclimate reconstructions, and more. However, when GCM output is used, it typically has to be downscaled or assumptions made to reduce the temporal and spatial resolution to be usable in utility hydrology or operational models. 9. Climate models often provide monthly average changes that must then be combined with the observed record at a daily timestep to generate projected daily changes in climate. ------- applicability beyond the relatively large and highly-resourced water utilities that have led the movement to consider climate change explicitly in their infrastructure and operational planning. 4.1 The Instrumental Record The instrumental record is often used to define utility-specific models for demand, hydrology, operations, and management. These models can be simple statistical relationships or complex physical models of local hydrology that incorporate vegetation cover, aspect, slope angle, and soil type in addition to precipitation and temperature. Regardless, the instrumental record forms a fundamental basis for most utility vulnerability assessments in one form or another. Utility vulnerability assessments have historically used the drought of record or flood of record to assess the vulnerability of a water system by applying the most extreme climate conditions on record to contemporary circumstances. In many cases, this type of analysis can be useful for identifying current vulnerabilities of a water system to current climate. ° Although none of the utilities discussed in this report relied exclusively on this assessment methodology, many utilities still routinely rely on this method due to its familiarity and ease of use.11 It is important to note that using only the instrumental record implicitly retains the assumption of climate stationarity - fluctuation within an unchanging envelope of variability - as observed typically within the last century or less.12 4.2 Paleoclimate Data Some researchers have begun developing proxy records for precipitation and streamflow based on tree rings and geological evidence.13 These researchers have found evidence of dry periods prior to the 20th century that were more persistent and perhaps more intense than what is observed in the instrumental record (e.g., Woodhouse and Lukas, 2006a). Because natural climate variability can often be quite large compared to changes projected from climate change, some utilities have used paleoclimate data as a complement to their climate vulnerability assessments. 10. With a changing climate, there should be more caution about interpreting the historic record. What may have been a 100-year drought or flood calculated from the historic record may shift significantly with projected climate change, or even by adding the paleorecord to historic data. 11. For example, both Contra Costa Water District and Santa Clara Valley Water District in California make conservative assumptions about "worst-case scenarios" for drought planning instead of engaging in climate change-specific studies as described in this report. 12. The exception to this would be to identify trends in the historic record and extrapolate them into the future. None of the utilities surveyed carried out such an approach. 13. The scientific field of dendrohydrology involves using tree-ring qualities, such as the width of the annual growth rings, to estimate pre-instrumental streamflow of a specific river. 10 ------- For example, EBMUD used paleoclimate reconstructions of overflow on the Sacramento River to examine the resilience of their drought planning sequence (Meko and Woodhouse, 2005). Their use of paleoclimate data was largely qualitative, but the analysis demonstrated EBMUD's drought planning sequence would likely function adequately under the drought conditions evidenced in the paleorecord. Denver Water has also used 375-year tree ring reconstructions of the Colorado and Platte River basins in conjunction with their water supply simulation model - PACSM - to estimate the frequency and severity of drought within Denver Water's collection system (Woodhouse and Lukas, 2006a). Boulder Utilities used a 300-year tree ring reconstruction of hydrology for similar purposes (Hydrosphere, 2003). Tree-ring reconstructions hold much promise because of an excellent network of tree-ring data, especially (but not limited to) the western United States (TreeFlow, 2010). The use of paleoclimate reconstructions, however, is limited because tree-ring reconstructions are not available for all areas, they can be difficult to translate into the modeling environment for water management, good relationships with hydrology have only been established for the growing season - giving an incomplete record of annual hydrology, and the instrumental records necessary for high-quality calibration of tree-ring data can be difficult to obtain. For these reasons, and because of the specialized skill set necessary to conduct such studies, paleoclimate data may not be feasible for many water utilities (Garrick and Jacobs, 2005). 4.3 Literature Reviews Literature reviews can generate information for projected future climate conditions for many utilities at low cost. The most common sources for initiating climate literature reviews are the Intergovernmental Panel on Climate Change Fourth Assessment Report (IPCC, 2007) and the U.S. Global Change Research Program's report Global Climate Change Impacts in the United States (Karl et al., 2009). The National Oceanic and Atmospheric Administration (NOAA) is currently developing a climate information web portal (http://www.climate.gov), and NOAA's Regional Integrated Sciences and Assessments (RISA) program offices often provide regional literature and data resources (http://www.climate.noaa.gov/cpo_pa/risa/). With the dramatic increase of climate research across the United States, regional specification of GCM output can often be found in the scientific peer-reviewed literature or on web portals. For example, EBMUD decided not to expend resources on developing its own downscaled climate projections because sufficient data could be obtained in works of the IPCC (2007), the U.S. National Research Council (NRC, 2002), and the U.S. Geological Survey (Dettinger, 2005). Using these data, EBMUD ran a sensitivity analysis on their water system using the following factors: > Increased customer demand resulting from a 4°C (7.2°F) increase in air temperature between 1980 and 2040; > Changes in the timing of runoff in the Mokelumne River corresponding to 2°, 3°, and 4°C (3.6°, 5.4°, and 7.2° F) increases in air temperature; and 11 ------- > Reductions in Mokelumne River runoff corresponding to a 10% and 20% decrease in precipitation. EBMUD combined these three factors to derive seven scenarios that represented an adequate range of conditions to analyze the potential sensitivity of the EBMUD water system to climate change. 4.4 Climate Projections The development of watershed-specific projected future climate conditions is the most computationally and resource intensive approach to gathering climate information. Data obtained by this method were employed by many utilities, including the City of Boulder Utilities Division (Smith et al, 2009), Denver Water (Laurna Kaatz, personal communication, April 8, 2009), SPU (Palmer, 2007), MWRA (Stephen Estes-Smargiassi, personal communication, February 26, 2010), NYCDEP (2008), PWB (Palmer and Hahn, 2002), and LCRA (CH2M Hill, 2008). Climate projections consist of three main elements: the use of socioeconomic scenarios of GHG emissions, the selection of GCMs, and the region-specific downscaling of climate model output. Each of these elements is described in greater detail below. Also described below is the innovative combination of climate projections with other data sources for analytical purposes. It is important to note that the complexity and resource intensity of making climate projections for water resources is being rapidly reduced. For example, the U.S. Bureau of Reclamation now provides a downscaled dataset of statistically downscaled monthly climate projections from 1950 to 2099 that consists of 112 climate projections based on 16 climate models and three GHG emissions scenarios downscaled to a spatial resolution of 1/8 degree across the coterminous United States (Maurer et al, 2007). The National Center for Atmospheric Research (NCAR) is also developing a dataset of GCMs combined with different regional climate models (RCMs) (UCAR, 2007). As these datasets and tools become more widespread, the time and resource threshold to engage in climate projections likely will be lowered. CH2M Hill's assessment of the LCRA system used the Bureau of Reclamation dataset as the foundation for its analysis. 4.4.1 Climate scenarios Utility climate scenarios generally consist of two distinct elements: GHG emissions scenarios and defined time periods. Probabilities have not been assigned to the emissions scenarios (they are all assumed equally likely to occur) and, as a result, probabilities have not been assigned to climate change scenarios.14 GHG emissions scenarios have been developed to identify plausible future global energy use, economic growth, land-use change, population growth, technological innovation, and other 14. NCAR has developed probabilities of regional changes in temperature and precipitation based on analysis of GCM output (Tebaldi et al., 2005, 2006), but these probabilities were not used by the utilities examined in this study. 12 ------- factors that could affect future emissions of GHGs. Most utilities in our analysis used GHG emissions scenarios defined in the 2000 IPCC Special Report on Emissions Scenarios (SRES; see Box I).15 In general, water utilities have selected two or three GHG emissions scenarios to provide substantially different future GHG emission levels. Using multiple GHG emissions scenarios avoids putting too much emphasis on any single projection of the future. For example, the SPU study used emissions scenarios A2 and Bl; NYCDEP used emissions scenarios A2, A1B, and Bl; and LCRA used emissions scenarios A2 and Bl. While, by definition, these scenarios are equally plausible, water utilities have generally selected both a low emissions/environmentally friendly scenario (i.e., Bl) and a high emissions/development-oriented scenario (i.e., A2) to ensure representation of a wide range of plausible futures. Box 1. IPCC SRES scenarios. > The Al storyline and scenario family describes a future world of very rapid economic growth, global population that peaks in mid-century and declines thereafter, and the rapid introduction of new and more efficient technologies. Major underlying themes are convergence among regions, capacity building, and increased cultural and social interactions, with a substantial reduction in regional differences in per capita income. The Al scenario family develops into three groups that describe alternative directions of technological change in the energy system. The three Al groups are distinguished by their technological emphasis: fossil intensive (A1FI), non-fossil energy sources (AIT), or a balance across all sources (A1B). > The A2 storyline and scenario family describes a very heterogeneous world. The underlying theme is self-reliance and preservation of local identities. Fertility patterns across regions converge very slowly, which results in continuously increasing global population. Economic development is primarily regionally oriented and per capita economic growth and technological change are more fragmented and slower than in other storylines. > The Bl storyline and scenario family describes a convergent world with the same global population that peaks in midcentury and declines thereafter, as in the Al storyline, but with rapid changes in economic structures toward a service and information economy, with reductions in material intensity, and the introduction of clean and resource-efficient technologies. The emphasis is on global solutions to economic, social, and environmental sustainability, including improved equity, but without additional climate initiatives. > The B2 storyline and scenario family describes a world in which the emphasis is on local solutions to economic, social, and environmental sustainability. It is a world with continuously increasing global population at a rate lower than A2, intermediate levels of economic development, and less rapid and more diverse technological change than in the Bl and Al storylines. While the scenario is also oriented toward environmental protection and social equity, it focuses on local and regional levels. Source: Nakicenovic et al., 2000 (p. 4-5). 15. IPCC is now moving toward a set of scenarios based on radiative forcing and not tied to specific socioeconomic scenarios (Moss et al., 2010). 13 ------- Furthermore, nearly all utilities utilizing climate projections looked at two or three future time periods, often defined by relevance to their planning processes. Typically, projections were made for approximately 30 years for operations or strategic planning and approximately 50 years for infrastructure and capital improvement plans. For example, the SPU study used three time periods (2025, 2050, and 2075); Portland used two time periods (2025 and 2045); and LCRA used two time periods (2050 and 2080). Each time period generally represented an average of 10 of more years of climate data to smooth out year-to-year variation in model output. 4.4.2 Global climate models After selecting appropriate scenarios, the water utilities generally identified a subset of GCMs to use in their study. GCMs are generally selected from the World Climate Research Programme's Coupled Model Intercomparison Project (CMIP3).16 CMIP3 contains data contributed by 23 of the world's leading climate models (see Table 1). Table 1. GCMs in the World Climate Research Programme's CMIP Phase Three Model IPCC First designation published Sponsor BCC-CM1 BCCR-BCM2.0 CCSM3 CGCM3.1(T47) CGCM2.1(T63) CNRM-CM3 CSIRO-MK3.0 ECHAM5/MPI-OM ECHO-G FGOALS-gl.O GFDL-CM2.0 GFDL-CM2.1 GISS-AOM 2005 2005 2005 2005 2005 2004 2001 2005 1999 2004 2005 2005 2004 Beijing Climate Center Bjerknes Centre for Climate Research National Center for Atmospheric Research (NCAR) Canadian Centre for Climate Modeling and Analysis Canadian Centre for Climate Modeling and Analysis Meteo-France/Centre National de Recherches Meteorologiques (CNRM) Commonwealth Scientific and Industrial Research Organisation (CSIRO) Atmospheric Research Max Plank Institute for Meteorology Meteorological Institute of the University of Bonn, Meteorological Research Institute of the Korea Meteorological Administration (KMA), and Model and Data Group National Key Laboratory of Numerical Modeling for Atmospheric Sciences and Geophysical Fluid Dynamics (LASG)/mstitute of Atmospheric Physics U.S. Department of Commerce/NCAA Geophysical Fluid Dynamics Laboratory (GFDL) NOAA GFDL National Aeronautics and Space Administration (NASA) Goddard Institute for Space Studies (GISS) Country China Norway USA Canada Canada France Australia Germany Germany/ Korea China USA USA USA 16. Note that CMIP is in the process of updating model outputs in anticipation of the Fifth Assessment Report of the IPCC. The next iteration of CMIP will be called CMIP5 - skipping "CMIP4" - in order to align the numbering with the IPCC reports. 14 ------- Table 1. GCMs in the World Climate Research Programme's CMIP Phase Three (cont.) Model IPCC First designation published Sponsor GISS-EH GISS-ER INM-CM3.0 IPSL-CM4 MIROC3.2 (high resolution) MIROC3.2 (medium resolution) MRI-CGCM2.3.2 PCM UKMO-HadCM3 UKMO-HadGEMl 2004 2004 2004 2005 2004 2004 2003 1998 1997 2004 NASA/GISS NASA/GISS Institute for Numerical Mathematics Institut Pierre Simon Laplace Center for Climate System Research (University of Tokyo), National Institute for Environmental Studies, and Frontier Research Center for Global Change (JAMSTEC) Center for Climate System Research (University of Tokyo), National Institute for Environmental Studies, and Frontier Research Center for Global Change (JAMSTEC) Meteorological Research Institute NCAR Hadley Centre for Climate Prediction and Research/Met Office Hadley Centre for Climate Prediction and Research/Met Office Country USA USA Russia France Japan Japan Japan USA UK UK This climate model output, however, can vary dramatically by region. Consequently, for this and other practical considerations, most water utilities select a subset of climate models to use in their study. Three common methods for selecting models are described below. One method for selecting climate models is hindcasting, which refers to using climate models to project climate over years for which instrumental data are available. The ability of various models to replicate observed climate is often used as a means of selecting a subset of climate models for use in the study under the justification that those models more accurately model the climate of that watershed or region. For example, in the first phase of its work, the NYCDEP Climate Change Task Force selected three GCMs namely the National Center for Atmospheric Research (NCAR), Goddard Institute of Space Studies (GISS) and European Center Hamburg Model (ECHAM) (NYCDEP, 2008). In the second phase of work, NYCDEP evaluated the -20 GCMs for various combinations of five meteorological variables (precipitation, average, maximum and minimum temperatures and wind speed), over four seasons using probability based skill scores. Since no single model performs best for all variables and seasons they arrived at the most objective way to choose a subset of models to identify the models with the highest mean skill scores (averaged across all variables). A second method for selecting climate models is to purposefully select models that project a range of climate conditions. This was the strategy employed in the SPU study (Polebitski et al., 2007a). The three GCM-scenario combinations selected by the UW-CIG were the GISS Model ER/B1 combination, the MPIECHAM5/A2 combination, and the Institut Pierre Simon Laplace (IPSL) CM4/A2 combination. These models have performed well in other studies when replicating the temperature and precipitation trends of the Pacific Northwest during the 20th century (Mote et al., 2005). The MPI ECHAM5/A2 model represents a middle-of-the-road 15 ------- scenario with moderate warming and precipitation increase. The IPSL-CM4/A2 model scenario is significantly wetter and warmer, and the GISS-ER/B1 scenario is slightly drier and warmer. The LCRA study (CH2M Hill, 2008) used a third method for selecting climate models that is effectively a hybrid of the two methods described above. Initially, six models were selected for further assessment based on expert judgment of the leading models for California selected by researchers at the Scripps Institute of Oceanography: NCAR CCSM3.0, NCARPCM1, GFDL CM2.1, CNRM CMS, CCSRMICRO3.2 (medres), and MPIECHAM5.17 Projected changes in temperature and precipitation from all six models under two emissions scenarios (A2 and Bl) and two time periods (2050 - averaged 2036-2065; and 2080 - averaged 2066-2099) were plotted onto a grid representing the 10th and 90th percentile of all 112 downscaled model simulations from the U.S. Bureau of Reclamation. Scenarios that plotted consistently outside of this range were discarded, resulting in the selection of two climate models for the LCRA study that were considered to provide a reasonable range of potential climate change effects in the Colorado River Basin. 4.4.3 Downscaling Finally, all climate projections must be downscaled to the spatial and temporal resolution to be useful inputs into water system models. This is generally done using one of three downscaling methodologies: simple, statistical, or dynamic downscaling (see Box 2). For example, the PWB study used simple downscaling. It started with relatively coarse spatial resolution model results - on the order of three degrees or approximately 300-km grids (Palmer and Hahn, 2002). UW downscaled the data to one-degree or approximately 100-km grids using the Symap algorithm. However, even at one degree, the model results were still coarse. Consequently, the climate change signal for each model and time period was calculated as the difference between average monthly temperature and precipitation from a control run simulating current climate and future climate model predictions. This yielded delta values by month for 2025 and 2045 that indicated the modeled change in temperature and precipitation. These deltas were then applied to an observed climate dataset used by UW as input into its Distributed Hydrology, Soil-Vegetation Model (DHSVM) to characterize the hydrology of PWB's Bull Run watershed. This represents a form of simple downscaling because it used coarse GCM data and related them to the observed record to define a dataset at the spatial and temporal scale of the utility's models. 17. NCAR CCSM3.0 and GFDL CM2.1 were ultimately chosen for use in the LCRA study. 16 ------- Box 2. Downscaling methodologies Simple downscaling involves adding projected changes in temperature and precipitation from GCMs to weather observations for a historical control period. This preserves observed spatial differences and locally observed seasonality of weather, and allows for the selection of different historic control periods for different purposes (e.g., a drought year or a year of heavy flooding). However, this method does not account for differences in climate change within a grid box. Statistical downscaling is a method by which climate change can be projected for a particular location or small geographic area. Variables found in GCMs such as pressure are correlated with observations such as observed temperature or precipitation. The correlation uses GCMs' simulation of current climate. This correlation is used to project changes in the predicted variable(s) based on the estimated changes in GCMs. Statistical downscaling is computationally easier than dynamic downscaling, but does not allow for changes in the relationship between the predictor and predicted variables. Dynamic downscaling uses RCMs to give higher resolution projections than GCMs. These models are like GCMs, but simulate only a portion of the globe and can therefore have much smaller grid boxes than GCMs and can incorporate sub-GCM-grid scale variables (e.g., mountain ranges, lakes). RCMs are "nested" within GCMs using "boundary conditions" from GCMs to drive them. RCMs typically have resolutions of 50 km or less. RCMs have the advantage of simulating dynamic relationships between climate variables, but require extensive computing power. A key limitation of all three downscaling methods is that they depend on GCMs. Errors in GCMs are typically not corrected by downscaling techniques, although each technique will supply sub-GCM-grid scale variations in climate. NYCDEP used a similar change factor approach, using the difference in daily GCM control and future scenarios to calculate change factors that were subsequently applied to local meteorological records used to drive watershed and reservoir models. In the first phase of NYCDEP climate change work, monthly change factors (additive or multiplicative) were calculated from pooled daily data associated with each month in each scenario. In the second phase of the work the change factor methodology was evaluated, and an improved method was developed that calculated 25 additive change factors for 25 equally spaced bins of the frequency distribution associated with each month's pooled data. City of Boulder Utilities Division also used simple downscaling, but combined GCM output with paleoclimate reconstructions. The LCRA study used a different technique for downscaling. The selection of GCMs for the LCRA study ultimately led to eight datasets - two models (NCAR CCSM3.0 and GFDL CM2.1), two scenarios (A2 and Bl), and two time periods (2050 - averaged 2036-2065; and 2080 - averaged 2066-2099). These eight datasets were statistically downscaled from the model grid scale to regional or watershed scales by retrieving the appropriate temporal-spatial data in the archived U.S. Bureau of Reclamation downscaled dataset (1/8 degree; 12-km grids).18 A suite 18. The Bureau of Reclamation developed an archive of statistically downscaled and bias-corrected climate projections for the contiguous 48 states at 1/8 degree resolution, which is what was used in the LCRA study (see http://gdo-dcp.ucllnl.org/downscaled cmip3_projections/dcplnterface.html). It should be noted that it is possible to statistically downscale GCM projections using a variety of techniques independent of the U.S. Bureau of Reclamation dataset. The SPU study provides an example of a different technique. 17 ------- of data processing tools and geographic information systems were applied to these data to delineate the datasets for the State of Texas and the Colorado River Basin. Simulated future climate was analyzed by comparing the 2050 and 2080 periods with a reference period of 1970-1999 (CH2M Hill, 2008). This use of sophisticated statistical techniques to downscale GCM output to local watersheds is known as statistical downscaling. The SPU study conducted its own method of statistical downscaling (Polebitski et al., 2007b). First, the SPU study did a bias correction of GCM simulation of hindcast climate using historic data. Those bias corrections were then applied to coarse (GCM grid scales) future climate data. The bias-corrected coarse grid was then related to a regional grid (1/8 degree) using scaling factors (additive factors for temperature and multiplicative factors for precipitation). The regional, bias-corrected future climate data are then related to specific weather station locations and again bias-corrected to yield a monthly time-series at the station level. These station level changes can be combined with historic data to yield scenarios of daily data that preserve historic climate variability. For its statistical downscaling, the MWRA used a non-parametric method for generating the requisite daily weather sequences, specifically the K-nearest neighbor (K-NN) bootstrapping technique (Yates et al., 2003). This method samples the historical record in order to generate a large set of individual weather sequences that replicates the statistical characteristics of local weather but are consistent with the range of GCM-derived temperature and precipitation trends. Using probability density functions derived from 21 GCMs following the procedure outlined in Tebaldi et al. (2005), the K-NN algorithm was used to develop individual sequences of weather variables for the key weather station locations that are used by the system simulation model (WEAP) of MWRA's supply system. The final downscaling methodology is dynamic downscaling. Dynamic downscaling uses RCMs to translate GCM output into higher-resolution regional climate data. Dynamic downscaling captures the effects of mesoscale features such as narrow mountain ranges, complex land/waterbody interactions, or variations in land use and land cover in a way that GCM currently cannot. Until recently, the computing power necessary for dynamic downscaling to provide useful data was prohibitive and the results were not high-resolution enough to improve upon statistical downscaling approaches. Current research into RCMs and dynamic downscaling shows promise. For example, UW-CIG is engaged in research to develop a fully functional RCM for the Pacific Northwest. NYCDEP also worked with Columbia University's Center for Climate Systems Research (CCSR) and the City University of New York to establish a dynamic downscaling strategy for future assessments. 4.4.4 Combining methods: Climate projections with paleoclimate data The study conducted for the City of Boulder combined climate model output with paleoclimate data (Smith et al., 2009). To our knowledge, this is the first study to do so in the United States. The study matched paleoclimate reconstructions of streamflow in Boulder Creek with years during the observed climate record that had similar streamflow. The climate in the similar years 18 ------- was used as a proxy for the climate of the paleorecord. In other words, the reconstruction of streamflow in a paleo year (e.g., 1635) was then matched to a year in the observed record with a similar annual streamflow. A best match between reconstructed streamflow and years in the observed record did not always exist. Therefore, the technique used a random process to select years from the observed record with similar characteristics. The model was run many times to yield different combinations of the observed record and mimic the reconstruction. This introduces many combinations of reconstructions, which allows for more variability in results. The Boulder study used a simple downscaling of GCM data. The changes in temperature and precipitation were added to or multiplied by temperature and precipitation, respectively, in the proxy record. This yielded a set of climate change scenarios that mimicked the year-to-year variability of the paleorecord, but with long-term average changes in climate from the GCMs imposed on the set. The reconstruction of streamflow for Boulder Creek utilized a 437-year record from 1566 to 2002 (Woodhouse and Lukas, 2006b). The reconstruction technique estimates 65% (i.e., has an o R of 0.65) of the observed streamflow variance. The study used output from three GCMs under three emissions scenarios for the decades of the 2030s and 2070s. The selected emissions scenarios (A2, A1B, and Bl) represented a wide range of future conditions. As with SPU, the GCMs were selected to capture a wide range of potential changes in regional climate. While all climate models project that the central Rocky Mountains and adjacent plains will become warmer, the models disagree as to whether precipitation will increase or decrease. To overcome this, the study used a relatively wet model (the Canadian Climate Model), a relatively dry model (GFDL Version 0), and a model roughly in the middle of the range of model output (GFDL Version 1). 19 ------- 5. Modeling Changes in Water Resources Water utilities use a variety of models and analytical techniques for various purposes. As a general organizational distinction, it simplifies things to examine how a utility approaches water quantity and water demand separately.19 5.1 Water Quantity Water utilities may use one or more water management model(s) for day-to-day operations and long-term water planning. Many, but not all, also use a hydrologic model to translate temperature and precipitation data into streamflow, reservoir storage, evaporative loss, and other variables input directly into water system management models. The specific constellation of models used by a utility often determines what climate information it needs to carry out a climate change vulnerability assessment. PWB, for example, input temperature and precipitation data into the DHSVM to generate averaged annual hydrographs. DHSVM is a physically based hydrology model that characterized the entire Bull Run watershed into 150-m grid cells with grid-specific data on soil and vegetation type, soil depth, vegetation height, and surface elevation and slope. DHSVM represents the more detailed end of utility hydrology models that have been empirically characterized for a number of parameters in addition to temperature and precipitation. PWB applied the UW-CIG developed climate altered streamflows through DHSVM and then input the data into PWB's Supply and Transmission Model (STM), an operational model of PWB's terminal water storage and groundwater resources. STM operates at a daily timestep and simulates the flow of water throughout the water transmission system using seasonally varying rule curves for the reservoirs, as well as modeling water releases for hydropower production and instream flows. Groundwater operations are coordinated with reservoir operations to enable the investigation of a variety of operating alternatives using variables such as length of drawdown period, amount of groundwater pumped during drawdown, minimum storage during drawdown, and water used during drawdown. As a part of a Water Research Foundation project, the PWB used the STM data sets to look at how the WEAP platform worked in comparison with the outputs of the STM, however, the WEAP model was not used for the initial climate change study done by UW-CIG. For purposes of this analysis it is important to note the use of multiple operational models for different purposes by a single utility. 19. Some utilities in this study, such as NYCDEP, were also concerned about climate change effects on water quality. However, generally speaking, water quality involves a greater variety of models (e.g., reservoir specific eutrophication models) and/or micro-climate processes that often occur at a higher resolution than can be captured by current GCMs or RCMs (e.g., intense precipitation events). The scope of this report did not allow us to fully explore water quality modeling. 20 ------- In order to use SPU's systems model, climate change parameters needed to be converted to inflow. UW-CIG used the DHSVM hydrology model to produce climate-altered hydrologic datasets based on climate model output for use in the Conjunctive Use Evaluation (CUE). SPU then used the results from DHSVM in its CUE systems model - a weekly timestep simulation model of the Cedar and Tolt River systems - for water supply planning. CUE is used for 90 calculating and evaluating the firm yield and reliability of Seattle's water supply system and potential future water supply projects. CUE results indicated that yield decreased under all climate change scenarios for all time periods. SPU also ran several planning scenarios through CUE to determine whether available supply could be increased to compensate for anticipated shortfalls. The use of river gauge data as input into CUE was a policy decision for SPU, because it also has the capability to use the Seattle Forecast Model (SEAFM), a proprietary hydrology model calibrated to the SPU watersheds and used in operational forecasting and operations planning. UW-CIG chose to use DHSVM instead of SEAFM (or its successor models) for its climate change vulnerability analysis because UW-CIG developed DHSVM and it was available to all watersheds in the broader UW-CIG led regional study in which SPU participated. NYCDEP has a very complex modeling environment. It currently does analysis using the Variable Source Loading Function (VSLF) or the SWAT watershed models. These can in turn provide inputs to one of two reservoir models: a one-dimensional reservoir eutrophication model; or a two-dimensional reservoir turbidity transport model (CEQual W2). The OASIS reservoir system operation model requires measured or simulated water inputs to the reservoirs and simulates the storage and transfer of water between the 19 reservoirs comprising the NYC water supply system. The watershed models, take daily temperature and precipitation data to generate streamflow and evapotranspiration, as well as a number of water quality parameters. Outputs of the watershed models can be used to drive both the OASIS system model and reservoir water quality models. A knowledge of reservoir operations is an important determinate of reservoir water quality, and a needed input to the reservoir models. For future climate simulations, where historical operations are not known, OASIS simulation are needed to specify reservoir operation scenarios associated with a given climate scenario. All of the above models were in use at NYCDEP to evaluate the impacts of changes in watershed management, land-use, and reservoir operation policies. This integrated modeling system was easily adapted for evaluation of the effects of climate change, once credible future scenarios of the meteorological variables needed to drive the models were available. 5.2 Water Demand Changes in climate can also have an important effect on water demand. Demand means different things to different utilities depending upon their customer base and the timeframe of their assessment. Several utilities evaluated in this report provided interesting examples of different 20. Firm yield is based on the 98% reliability standard. 21 ------- approaches or concerns to project changing water demand. EBMUD was focused primarily on exterior water use by urban and suburban water customers, while SPU and PWB were focused on the balance between population growth and climate altered demand over time. According to EBMUD, customer demands are projected to vary predominantly based on temperature change. While indoor water use is not expected to change significantly with global warming, outdoor water use could change dramatically. To account for this, 2040 customer demands were re-normalized using a projected temperature increase of 2.15°C (3.87°F) between 2005 and 2040 with no change in precipitation. EBMUD's analysis indicated that a 20% reduction in precipitation had little influence on overall customer demands compared to the projected temperature change. The demand model suggested a 3.6% increase in customer demands by the year 2040 (EBMUD, 2009). The SPU study examined the effect of climate change on water demand using a dual approach of regression analysis and forecast modeling. First, SPU performed a regression analysis of peak season consumption for 1982-2007 using monthly consumption data, maximum temperature, and rainfall at SeaTac Airport for May through September. This relationship was assumed to hold in the future. SPU had already developed a demand forecasting model for its 2007 Water System Plan, which forecasted non-climate altered demand change over time. Under this model, demand was forecasted to decrease below historic levels through 2050, but increase above historic levels by 2075. Applying the results of the regression analysis to these forecasts adjusts demand slightly upward due to climate change in 2025 and 2050. But in 2075, the climate- induced increase becomes more significant on top of the additional forecasted upswing from the non-climate altered demand model. In 2002, PWB examined the effects of climate-altered streamflow on system operations from a demand perspective using their STM. The evaluation investigated the climate and population growth impact on demand and supply separately and then jointly to discern the discrete effects of climate change. The process for calculating the climate impact on demand was derived from the Joint Institute on the Study of Atmosphere and Oceans (JISAO) report, Impacts of Climate Variability and Change in the Pacific Northwest (JISAO-CIG, 1999). This study used seven featured years and the ECHAM4 climate scenario for detailed analysis. Growth in demand not considering changes in climate decreased the average minimum storage dramatically. This was exacerbated (approximately doubled) by the impacts of climate on hydrology and demand. PWB also utilized their water demand econometric model to estimate near- and long-term water demand that established a relationship between total water demand and selected economic, demographic, seasonal, and day-to-day weather variables. This econometric model is used as a forecasting tool by projecting changes in the economic and demographic variables. In order to gauge the effects of climate change on demand, past weather years with patterns that approximate climate projections are selected to forecast future demands. For example, in projecting demands in 2050, PWB uses 1948 as the wettest, or best-case scenario (corresponding to a 3.67°F/2.04°C increase in average peak season temperature), 1980 as a middle-of-the-road scenario (corresponding to a 4.83°F/2.68°C increase in average peak season temperature), and 22 ------- 1991 as the driest, or worst-case scenario (corresponding to a 6.72°F/3.73°C increase in average peak season temperature). 6. Summary We investigated assessments of climate change vulnerability conducted by eight utilities to determine the emergent characteristics of water utility climate change vulnerability assessments (see Table 2). We found that two standard approaches were used for the risk assessments: the top-down modeling assessment and the bottom-up threshold analysis. While both methods produce important risk information, the selection of an assessment model is often determined by available fiscal and technical resources as well as the policy environment faced by the utility. These two approaches are not mutually exclusive, and are often run in series or in parallel to provide a robust vulnerability assessment. Notably, top-down assessments come in a variety of types, including scenario analyses, sensitivity analyses, and paleoclimate or historic analyses. Sources of climate data used by the utilities included the instrumental record, paleoclimate data, literature reviews, and climate projections from GCMs. While information from each of these sources played an important role in one or more utility vulnerability assessments, climate projections offered the most computationally demanding approach. Climate projections required the selection of GHG emissions scenarios and defined time periods, the selection of a subset of GCMs (typically by hindcasting, selecting a broad range of outputs, or eliminating outliers), and downscaling the GCM data to the spatial and temporal scale necessary for input into utility hydrology, management, or operations models through simple, statistical, or dynamical downscaling techniques. Finally, the climate information was run through utility-specific models to determine the effects of projected climate conditions on water supply and demand. We found that on the supply side, the models typically used included both hydrology and operational or management models. Demand modeling by utilities ranged from service-area-specific correlations between customer demand and temperature/precipitation to utility-specific demand models developed for long- range planning. 23 ------- Table 2. Aspects of utility vulnerability assessment efforts Boulder Denver Water MWRA Assessment approach Top-down (paleo) Top-down (sensitivity, scenario in process) Top-down NYCDEP Top-down PWB SPU LCRA (scenario) Top-down (sensitivity) Top-down (scenario) Bottom-up (scenario) EBMUD Bottom-up a. b. (sensitivity) The 1% annual Climate model selection Purposeful selection for range Current study: Purposeful selection for range All models (Yates/NCAR project) Probability based skill score approach Purposeful selection for range Purposeful selection for range Hindcast and purposeful selection for range Conducted sensitivity analysis GHG emissions scenarios A2,A1B, Bl Current study: A2, A1B,B1 A1B A2,A1B, Bl 1% annual increase in C02a 2007 study: A2, Bl,2002 study: A2 A2,B1 N/A increase in C02 scenario was Because EBMUD chose to use a applicable. Downscaling Simple Current study: statistical Statistical Delta change factor methodology by month and frequency distribution, statistical Statistical Statistical Statistical (Bureau of Reclamation dataset) N/A Time periods 2030, 2070 In process study: 2040 and 2070 Weekly projections through 2060 2046-2065 and 2081- 2100 2025, 2045 2007 study: 2000, 2025, 2050, 2075 2002 study: 2000, 2020, 2040 2050, 2080 N/A System models Boulder Creek Model PACSM, ESP WEAP OASIS, CEQual W2, 1-D Eutrophication STM, WEAP Conjunctive Use Evaluation (CUE) WAMWRAP, EBMUDSim, WEAP often used prior to development of the SRES sensitivity analysis, the climate modeling Runoff models CLIRUN Sacramento Soil Moisture coupled with Anderson Snow- 17, WEAP abed watershed model VSLF, SWAT DHSVM DHSVM Variable Infiltration Capacity N/A scenarios. ; aspects of this table are not 24 ------- 7. Recommendations for Further Study This report describes a range of practices used by eight water utilities to evaluate the effects of climate change on water quantity and water demand using hydrological, management and operational models. Climate change is a complex issue and will require ongoing work to establish reliable practices for incorporating climate change into water utility decisions and planning. Additional questions to develop this line of inquiry further include: • How are decision makers and rate payers responding to these types of analyses? • Would the results vary if different methods were used by the same utility? • How do the different methods compare in effectiveness, over time? • What other decisions, utility models, etc. need climate vulnerability assessments? • How are vulnerability assessments being conducted for understanding climate impacts on managing: o Water quality? o Intense precipitation and storms (stormwater, floods, wind)? o Sea level rise? As other water resource managers gain more experience with climate vulnerability assessments, it would benefit the water resource management field to continue to document lessons learned in an effort to create, over time, a solid foundation of acceptable industry practices. 25 ------- 8. References Barsugli, J., C. Anderson, J.B. Smith, and J.M. Vogel. 2009. Options for Improving Climate Modeling to Assist Water Utility Planning for Climate Change. Final Report. Prepared for Water Utility Climate Alliance. December. CH2M Hill. 2008. Climate Change Study: Report on Evaluation Methods and Climate Scenarios. Final Report. Lower Colorado River Authority and San Antonio Water System. November. Dettinger, M.D. 2005. From climate-change spaghetti to climate-change distributions for 21st Century California. San Francisco Estuary and Watershed Science 3(1). Available: http://repositories.cdlib.org/jmie/sfews/vol3/issl/art4. EBMUD. 2009. WSMP 2040: Water Supply Management Program 2040 Plan. East Bay Municipal Utility District. October. Freas, K., B. Bailey, A. Munevar, and S. Butler. 2008. Incorporating climate change in water planning. Journal of the American Water Works Association 100:6. June. Garrick, D. and K. Jacobs. 2005. Tree-ring records inform water management decisions: Workshop strives for better integration of past with present. Southwest Climate Outlook. May. Hydrosphere. 2003. Simulation of Hypothetical Climate Change Scenarios Using the Boulder Creek Water shed Model. Hydrosphere Resource Consultants (now AMEC), Boulder, CO. IPCC. 2007. Climate Change 2007: The Physical Science Basis. Contribution of Working Group I to the Fourth Assessment Report of the Intergovernmental Panel on Climate Change, S. Solomon, D. Qin, M. Manning, Z. Chen, M. Marquis, K.B. Averyt, M. Tignor and H.L. Miller (eds.). Cambridge University Press, Cambridge, UK and New York. JISAO-CIG. 1999. Impacts of Climate Variability and Change: Pacific Northwest. A Report of the Pacific Northwest Regional Assessment Group for the US Global Change Research Program. Joint Institute for the Study of Atmosphere and Oceans - Climate Impacts Group. Office of Global Programs at the National Oceanographic and Atmospheric Administration. Seattle, WA. November. Karl, T.R., J.M. Melillo, and T.C. Peterson. 2009. Global Climate Change Impacts in the United States. Cambridge University Press, New York. Major, D.C., R. Horton, and C. Rosenzweig. 2007. Integrated Modeling Project for Water Quantity and Quality: Guidelines for Project Final Report. Climate Impact Assessment of Environmental Infrastructure Systems. New York City Department of Environmental Protection. June 29. Maurer, E.P., L. Brekke, T. Pruitt, and P.B. Duffy. 2007. Fine-resolution climate projections enhance regional climate change impact studies. Eos Trans. AGU88(4iy.504. 26 ------- Means III, E.G., M.C. Laugier, and J.A. Daw. 2010. Evaluating Decision Support Planning Methods for Incorporating Climate Change Uncertainties into Water Planning. Final Report. Prepared for Water Utility Climate Alliance. January. Meko, D.M. and C.A. Woodhouse. 2005. Tree-ring footprint of joint hydrologic drought in Sacramento and Upper Colorado River basins, western USA. Journal of Hydrology 308:196-213. Miller, K. and D. Yates. 2006. Climate Change and Water Resources: A Primer for Municipal Water Providers. Jointly sponsored by Awwa Research Foundation, Denver, CO and University Corporation for Atmospheric Research, Boulder, CO. Awwa Research Foundation, American Water Works Association, and IWA Publishing. Milly, P.C.D., J. Betancourt, M. Falkenmark, R.M. Hirsch, Z.W. Kundzewicz, D.P. Lettenmaier, and R.J. Stouffer. 2008. Stationarity is dead: Whither water management? Science 319:573-574. Moss, R.H., J.A. Edmonds, K.A. Hibbard, M.R. Manning, S.K. Rose, D.P. van Vuuren, T.R. Carter, S. Emori, M. Kainuma, T. Kram, G.A. Meehl, J.F.B. Mitchell, N. Nakicenovic, K. Riahi, S.J. Smith, R.J. Stouffer, A.M. Thomson, J.P. Weyant, and T.J. Wilbanks. 2010. The next generation of scenarios for climate change research and assessment. Nature 463:747-756. Mote, P.W., E. Salathe, and C. Peacock. 2005. Scenarios of Future Climate for the Pacific Northwest. University of Washington: Climate Impacts Group. Nakicenovic, N., J. Alcamo, G. Davis, B. de Vries, J. Fenhann, S. Gaffm, K. Gregory, A. Grubler, T.Y. Jung, T. Kram, E.L. La Rovere, L. Michaelis, S. Mori, T. Morita, W. Pepper, H. Pitcher, L. Price, K. Raihi, A. Roehrl, H.-H. Rogner, A. Sankovski, M. Schlesinger, P. Shukla, S. Smith, R. Swart, S. van Rooijen, N. Victor, and Z. Dadi. 2000. Emissions Scenarios. A Special Report of Working Group III of the Intergovernmental Panel on Climate Change. Cambridge University Press, New York. UCAR. 2007. North American Regional Climate Change Assessment Program. Available: http://www.narccap.ucar.edu/. Accessed May 28, 2010. NRC. 2002. Abrupt Climate Change: Inevitable Surprises. National Academies Press: Washington, D.C. Available: http://www.nap.edu/openbook.php7isbn = 0309074347. NYCDEP. 2008. The New York City Department of Environmental Protection Climate Change Program Assessment and Action Plan: Report 1. May. Palmer, R.N. 2007. Final Report of the Climate Change Technical Committee. A report prepared by the Climate Change Technical Subcommittee of the Regional Water Supply Planning Process, Seattle, WA. Palmer, R.N. and M. Hahn. 2002. The Impacts of Climate Change on Portland's Water Supply: An Investigation of Potential Hydrologic and Management Impacts on the Bull Run System. Portland Water Bureau. January. 27 ------- Polebitski, A., L. Traynham, and R.N. Palmer. 2007a. Technical Memorandum #4: Approach for Developing Climate Impacted Meteorological Data and its Quality Assurance/Quality Control. A report prepared by the Climate Change Technical Subcommittee of the Regional Water Supply Planning Process, Seattle, WA. Polebitski, A., M.W. Wiley, and R.N. Palmer. 2007b. Technical Memorandum #2: Methodology for Downscaling Meteorological Data for Evaluating Climate Change. A report prepared by the Climate Change Technical Subcommittee of the Regional Water Supply Planning Process, Seattle, WA. Smith, J.B., K. Strzepek, L. Rozaklis, C. Ellinghouse, and K.C. Hallett. 2009. The Potential Consequences of Climate Change for Boulder Colorado's Water Supplies. Prepared for NOAA's Climate Program Office. February. Strange, E.M., D.R. Lane, and C.N. Herrick. 2009. Utility Guidance for Mitigating Catastrophic Vegetation Change in Watersheds. Project #4009. Water Research Foundation, Denver, CO. Stratus Consulting and MWH Global. 2009. Implications of Climate Change for Adaptation by Wastewater andStormwater Agencies. Water Environment Research Foundation Report No. CC2R08. Tebaldi, C., J.M. Arblaster, K. Hayhoe, and G.M. Meehl. 2006. Going to extremes: An intercomparison of model-simulated historical and future change extremes. Climatic Change 79:185- 211. Tebaldi, C., R.L. Smith, D. Nychka, and L.O. Mearns. 2005. Quantifying uncertainty in projection of regional climate change: ABayesian approach to analysis of multi-model ensembles. Journal of Climate 18:1524-1540. TreeFlow. 2010. Streamflow reconstructions from treerings website. Available: http://treeflow.info/. Accessed March 30, 2010. U.S. EPA. 2009. Proceedings of the First National Expert and Stakeholder Workshop on Water Infrastructure Sustainability and Adaptation to Climate Change. Publication No. EPA 600/R-09/010. Available: http://www.epa.gov/nrmrl/wswrd/wqm/wrap/workshop.html. Woodhouse, C.A. and JJ. Lukas. 2006a. Multi-century tree-ring reconstructions of Colorado Streamflow for water resource planning. Climatic Change. DOT: 10.1007/sl0584-006-9055-0. Woodhouse, C. and J. Lukas. 2006b. Streamflow Reconstructions for Boulder Creek, South Boulder Creek, and the Colorado River. Report to City of Boulder and Hydrosphere Resource Consultants. AMEC Consulting, Boulder, CO. WRF. 2010a. Analysis of Changes in Water Use Under Regional Climate Change Scenarios. Project #4263. Project in progress; completion year 2013. Water Research Foundation, Denver, CO. WRF. 2010b. Estimating Salinity Effects Due to Climate Change on the Georgia and South Carolina Coasts. Project #4285. Project Description; completion year 2012. Water Research Foundation, Denver, CO. 28 ------- WRF. 2010c. Evaluating Effects of Climate Change on Water Utility Planning Criteria and Design Standards. Project #4154. Project Description and Abstract; completion year 2010. Water Research Foundation, Denver, CO. WRF. 2010d. Incorporating Climate Change Information in Water Utility Planning: A Collaborative, Decision Analytic Approach. Project #3132. Project Description and Abstract; completion year 2010. Water Research Foundation, Denver, CO. WRF. 2010e. Project Snapshot: Joint Front Range Climate Change Vulnerability Study #4205. Available: http ://www. waterresearchfoundation. org/research/TopicsAndProj ects/proj ectSnapshot. aspx?pn = 4205. Accessed March 30, 2010. WRF. 201 Of. Water Utilities and Climate Change: A Research Workshop on Effective System Adaptation. Project #4228. Project Description; completion year 2010. Water Research Foundation, Denver, CO. Yates, D., S. Gangopadhyay, B. Rajagopalan, and K. Strzepek (2003). A technique for generating regional climate scenarios using a nearest neighbor algorithm, Water Resources Research, 39(7), 1199. 29 ------- |