<?xml version="1.0" encoding="UTF-8"?><!DOCTYPE article  PUBLIC "-//NLM//DTD Journal Publishing DTD v3.0 20080202//EN" "http://dtd.nlm.nih.gov/publishing/3.0/journalpublishing3.dtd"><article xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink" dtd-version="3.0" xml:lang="en" article-type="research article"><front><journal-meta><journal-id journal-id-type="publisher-id">JWARP</journal-id><journal-title-group><journal-title>Journal of Water Resource and Protection</journal-title></journal-title-group><issn pub-type="epub">1945-3094</issn><publisher><publisher-name>Scientific Research Publishing</publisher-name></publisher></journal-meta><article-meta><article-id pub-id-type="doi">10.4236/jwarp.2020.124016</article-id><article-id pub-id-type="publisher-id">JWARP-99361</article-id><article-categories><subj-group subj-group-type="heading"><subject>Articles</subject></subj-group><subj-group subj-group-type="Discipline-v2"><subject>Earth&amp;Environmental Sciences</subject></subj-group></article-categories><title-group><article-title>
 
 
  Implications of the New Operational Rules for Cantareira Water System: Re-Reading the 2014-2016 Water Crisis
 
</article-title></title-group><contrib-group><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>Karinne</surname><given-names>Reis Deusdará-Leal</given-names></name><xref ref-type="aff" rid="aff1"><sup>1</sup></xref><xref ref-type="corresp" rid="cor1"><sup>*</sup></xref></contrib><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>Luz</surname><given-names>Adriana Cuartas</given-names></name><xref ref-type="aff" rid="aff1"><sup>1</sup></xref></contrib><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>Rong</surname><given-names>Zhang</given-names></name><xref ref-type="aff" rid="aff2"><sup>2</sup></xref></contrib><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>Guilherme</surname><given-names>S. Mohor</given-names></name><xref ref-type="aff" rid="aff3"><sup>3</sup></xref></contrib><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>Luíz</surname><given-names>Valério de Castro Carvalho</given-names></name><xref ref-type="aff" rid="aff4"><sup>4</sup></xref></contrib><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>Carlos</surname><given-names>Afonso Nobre</given-names></name><xref ref-type="aff" rid="aff5"><sup>5</sup></xref></contrib><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>Eduardo</surname><given-names>Mario Mendiondo</given-names></name><xref ref-type="aff" rid="aff6"><sup>6</sup></xref></contrib><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>Elisângela</surname><given-names>Broedel</given-names></name><xref ref-type="aff" rid="aff1"><sup>1</sup></xref></contrib><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>Marcelo</surname><given-names>Enrique Seluchi</given-names></name><xref ref-type="aff" rid="aff1"><sup>1</sup></xref></contrib><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>Regina</surname><given-names>Célia dos Santos Alvalá</given-names></name><xref ref-type="aff" rid="aff1"><sup>1</sup></xref></contrib></contrib-group><aff id="aff6"><addr-line>Sao Carlos School of Engineering, University of Sao Paulo, USP, Sao Carlos, Brazil</addr-line></aff><aff id="aff3"><addr-line>University of Potsdam, Potsdam, Germany</addr-line></aff><aff id="aff4"><addr-line>Universidade de Taubaté, Taubaté, Brazil</addr-line></aff><aff id="aff2"><addr-line>Department of Hydrology and Water Resources, Nanjing Hydraulic Research Institute (NHRI) under the Ministry of Water 
Resources, The Ministry of Transport and National Energy Administration, Nanjing, China</addr-line></aff><aff id="aff5"><addr-line>Institute of Advanced Studies, University of Sao Paulo, Sao José dos Campos, Brazil</addr-line></aff><aff id="aff1"><addr-line>Brazilian Center for Monitoring and Early Warning of Natural Disaster-Cemaden, Sao José dos Campos, Brazil</addr-line></aff><pub-date pub-type="epub"><day>07</day><month>04</month><year>2020</year></pub-date><volume>12</volume><issue>04</issue><fpage>261</fpage><lpage>274</lpage><history><date date-type="received"><day>28,</day>	<month>February</month>	<year>2020</year></date><date date-type="rev-recd"><day>5,</day>	<month>April</month>	<year>2020</year>	</date><date date-type="accepted"><day>8,</day>	<month>April</month>	<year>2020</year></date></history><permissions><copyright-statement>&#169; Copyright  2014 by authors and Scientific Research Publishing Inc. </copyright-statement><copyright-year>2014</copyright-year><license><license-p>This work is licensed under the Creative Commons Attribution International License (CC BY). http://creativecommons.org/licenses/by/4.0/</license-p></license></permissions><abstract><p>
 
 
  During the 2014-2016 water shortage crisis, the Metropolitan Area of Sao Paulo (MASP) water supply system extracted pumping volume from the Cantareira System. Before the crisis, between 1984 and 2013, the reservoir’s average water extraction flow was 29.6 m
  <sup>3</sup>
  &#183;s
  <sup>-1</sup>. During the period of pumping volume usage, the average extraction flow was 16.2 m
  <sup>3</sup>
  &#183;s
  <sup>-1</sup>. Following the crisis, two new mitigation policies were implemented: a water extraction Resolution (in 2017) and a Resolution for water reallocation from another basin (in 2018). This study provides a novel investigation of the Cantareira System water crisis by assessing the mitigation policies impacts on storage level dynamics. The system storage level was evaluated using the reservoir simulation module of PDM-Cemaden hydrological model, assuming that the new policies had already been implemented prior to the crisis. A control simulation was run with observed in- and out-flow and operationally-practiced extraction flow. The storage level dynamics impacts were evaluated under 4 water mitigation policies scenarios varying the policies implementation starting date, the extraction flow range and including the water reallocation variable. Results showed that pumping volume would only need extraction during a short period (Scenarios I, III and IV), and considering the water reallocation, pumping volume extraction would not have been necessary (Scenario II). Although the pumping volume would still have been extracted during a short period, water shortage impact would have been lessened, had the policies been already implemented before the crisis. The water mitigation policies implementation supports the reservoirs storage management but does not guarantee that MASP water demand is fully met. Therefore, in order to effectively improve water security, further policies and practices to reduce water demand and enhance supply should be considered.
 
</p></abstract><kwd-group><kwd>Cantareira System</kwd><kwd> Drought</kwd><kwd> Sao Paulo</kwd><kwd> Brazil</kwd><kwd> Dead Storage</kwd><kwd> Water Crisis</kwd></kwd-group></article-meta></front><body><sec id="s1"><title>1. Introduction</title><p>Drought is a complex natural hazard that impacts ecosystems and human activities in several ways, mainly associated with hydrological impacts. Hydrological droughts (HD) can be defined as the low streamflow, low inflow levels, and reduced available groundwater [<xref ref-type="bibr" rid="scirp.99361-ref1">1</xref>]. The often unnoticed onset and slow development of HD converges to devastating impacts on agriculture, hydroelectric power generation, water supply, public health, navigation, and recreation. In recent years, many HD events occurred around the world. For example, the metropolitan area of Barcelona faced a severe drought that affected various economic sectors between 2007 and 2008. Important communication campaigns were put into place which led to a significant water demand reduction and to the setup of mechanisms for public participation for future water management. Although HD has been gaining attention in the last decades, one further step is to integrate science, management, and policy. A strong interface between policymakers and scientists is necessary to ensure that research better addresses the drought impacts management.</p><p>Prolonged dry and hot weather, and consequently less than normal water availability, have historically challenged Northeast Brazil, but in the last decade, the Southeast has also been experiencing similar issues. Between 2014 and 2016, the Southeast and Central-West Regions of Brazil experienced an unprecedented HD [<xref ref-type="bibr" rid="scirp.99361-ref2">2</xref>] [<xref ref-type="bibr" rid="scirp.99361-ref3">3</xref>] [<xref ref-type="bibr" rid="scirp.99361-ref4">4</xref>]. The unfavorable climatological condition, combined with political management challenges, seriously impacted water supply [<xref ref-type="bibr" rid="scirp.99361-ref5">5</xref>] [<xref ref-type="bibr" rid="scirp.99361-ref6">6</xref>], mainly to the Metropolitan Area of S&#227;o Paulo (MASP), hydropower generation, agriculture [<xref ref-type="bibr" rid="scirp.99361-ref7">7</xref>] [<xref ref-type="bibr" rid="scirp.99361-ref8">8</xref>] as well as the water quality [<xref ref-type="bibr" rid="scirp.99361-ref9">9</xref>]. For example, in the upper stream of the S&#227;o Francisco river basin, an important basin that connects Northeast and Southeast Regions, the reservoir recorded a historical minimum storage level (2.6%) in November 2014, challenging the hydropower generation, water supply and irrigation [<xref ref-type="bibr" rid="scirp.99361-ref10">10</xref>]. In the United States, as well as in Brazil, the recurrent HD in California has challenged human activities along the history. The 2012-2016 drought in California was deep, prolonged and unusually warm. Interties, water stored in aquifers and reservoirs, water market agreements, and long-term water conservation avoided, with few exceptions, emergency conditions in the cities. In 1976-1977, the California drought event included the driest year on the record and deep urban water conservation was achieved as a result of significant water demand management programs [<xref ref-type="bibr" rid="scirp.99361-ref11">11</xref>].</p><p>The Cantareira water supply system, located in Southeast Brazil (<xref ref-type="fig" rid="fig1">Figure 1</xref>), is used to be responsible for supplying water for 8.8 million people [<xref ref-type="bibr" rid="scirp.99361-ref12">12</xref>], corresponding to 46% of MASP (or 65% of Sao Paulo city demand). This number was reduced to 5.5 million people during the 2014-2016 HD, when it received support from other water systems, and returned to attending 7.5 million people after the drought. Cantareira System pumping volume was extracted, between July 2014 and December 2015, after depletion of its 982 hm<sup>3</sup> active storage. Pumping volume refers to reservoir negative water levels which cannot be extracted by gravity, but only pumped out, demanding more energy for the water extraction [<xref ref-type="bibr" rid="scirp.99361-ref13">13</xref>]. In May 2014, the first pumping volume quota was set operational, adding 182.5 hm<sup>3</sup>, and in October 2014 the second quota was activated, adding 105 hm<sup>3</sup>, reaching a total storage capacity of 1284.5 hm<sup>3</sup>. At the beginning of January 2016, both pumping volume quotas were recovered, because of the rainy season above the historical mean (SOND 2015 registered 28% above the Long-Term Mean—LTM) and the emergency practices implemented. At the end of the rainy season, in March 2016, the reservoir reached 36% of its active capacity.</p><p>The S&#227;o Paulo State Sanitation Company—SABESP, a mixed-economy company responsible for water supply, sewage collection and treatment in the MASP, which holds the Cantareira System water permit, adopted a set of practices to reduce consumption and preserve the Cantareira reservoirs storage. Among the practices implemented, stand out: 1) tariff discount for reduced consumption; 2) additional use of other MASP water systems [<xref ref-type="bibr" rid="scirp.99361-ref14">14</xref>], namely: the Upper Tiet&#234; System (capacity of 520 hm<sup>3</sup>), the Guarapiranga System (capacity of 171 hm<sup>3</sup>), the Rio Grande System (capacity of 112 hm<sup>3</sup>), the Upper Cotia System</p><p>(capacity of 16.5 hm<sup>3</sup>), and the Rio Claro System (capacity of 16 hm<sup>3</sup>); and 3) pressure decrease in the water distribution network to leakage reductions [<xref ref-type="bibr" rid="scirp.99361-ref15">15</xref>]. Notwithstanding the existence of other MAPS water supply systems, the Cantareira System is the main and largest source of water.</p><p>After the 2014-2016 HD, in order to mitigate the vulnerability of the Cantareira System, a water extraction rule and water reallocation from another basin were implemented. In 2017, the Water and Energy Resources Department of S&#227;o Paulo state (DAEE) and the National Water Agency (ANA) implemented a water extraction rule controlled by the hydrologically-conditioned reservoir storage level, via the ANA/DAEE resolution N˚925 [<xref ref-type="bibr" rid="scirp.99361-ref16">16</xref>]. Moreover, in 2018, water reallocation was implemented from the Jaguari reservoir, located in the Paraiba do Sul River Basin (Jaguari-PS), to Atibainha, a Cantareira System reservoir, guided by the resolution ANA N˚1931, issued in 2017 [<xref ref-type="bibr" rid="scirp.99361-ref17">17</xref>]. Notwithstanding the implementation of these mitigation policies, at the end of the 2018 dry season, the Cantareira System storage situation was worse than in 2017, reaching 33.5% of the active volume, which once again alarmed the authorities about the capacity to MASP water supply [<xref ref-type="bibr" rid="scirp.99361-ref18">18</xref>] [<xref ref-type="bibr" rid="scirp.99361-ref19">19</xref>] [<xref ref-type="bibr" rid="scirp.99361-ref20">20</xref>] [<xref ref-type="bibr" rid="scirp.99361-ref21">21</xref>].</p><p>Therefore, this study provides a novel investigation of the Cantareira System water crisis by assessing the mitigation policies impacts on storage level dynamics. The system storage level was evaluated using the Reservoir Simulation module of Probability Distributed Model (PDM-Cemaden) hydrological model, assuming that the new policies had already been implemented prior to the crisis. Therefore, this paper addresses three questions: 1) could the new water mitigation policies have avoided the pumping volume extraction, had they been implemented in the beginning of 2014 or 2) some years prior to the crisis? 3) are these rules effective to ensure the MASP water security?</p></sec><sec id="s2"><title>2. Material and Methods</title><sec id="s2_1"><title>2.1. The Cantareira Water Supply System</title><p>The Cantareira System (<xref ref-type="fig" rid="fig1">Figure 1</xref>) comprises four water producing reservoirs: Jaguari-Jacare&#237;, Cachoeira, Atibainha (these three constitute the “Equivalent System”), and Paiva Castro. The first three reservoirs capture and divert water through tunnels and channels, from some tributaries of the Piracicaba river basin to the Juqueri river basin, up to the Paiva Castro reservoir [<xref ref-type="bibr" rid="scirp.99361-ref12">12</xref>]. Finally, the water is pumped from Paiva Castro to the &#193;guas Claras reservoir, in the lift station Santa In&#234;s—ESI. This reservoir allows a continuous flow of 33 m<sup>3</sup>∙s<sup>−1</sup> to be maintained for 3 hours if the ESI stops. The water is finally transferred by gravity to the Guara&#250; water treatment plant (ETA) before distribution in the network. The total drainage area is approximately 2300 km<sup>2</sup>. The donating reservoirs upstream of the Piracicaba river basin also contribute with downstream reservoirs to supply cities within the Piracicaba, Capivari, and Jundia&#237; river basins (PCJ basins), which are managed together with Cantareira System.</p><p>In this study, the Cantareira System is considered as the equivalent system plus the Paiva Castro reservoir. The inflow is therefore the total of the reservoir’s inflow, the outflow is the total of reservoir’s outflow and the water extraction is the flow to the lift station Santa In&#234;s, used for the MASP water supply. <xref ref-type="fig" rid="fig2">Figure 2</xref> shows the time series of the observed in- and out-flow, the operationally-practiced extraction flow and also the storage level.</p><p>The climatological precipitation (1983-2013) over the Cantareira System is 1543 mm, with a markedly rainy season (between October and March), responsible for 74% of the total rainfall. The rainfall for the period from October 2013 to September 2014 was 39% below the climatological rainfall (<xref ref-type="fig" rid="fig3">Figure 3</xref>), which lead to inflow much lower than the LTM, and also below the minimum records [<xref ref-type="bibr" rid="scirp.99361-ref2">2</xref>]. The annual LTM inflow is 42.8 m<sup>3</sup>∙s<sup>−1</sup>. The average inflow between October 2013 and March 2014 was 19 m<sup>3</sup>∙s<sup>−1</sup>, 36% of the LTM for the period. More information about the meteorological conditions that lead to this precipitation deficit is presented in [<xref ref-type="bibr" rid="scirp.99361-ref2">2</xref>] and [<xref ref-type="bibr" rid="scirp.99361-ref4">4</xref>].</p></sec><sec id="s2_2"><title>2.2. Cantareira System Operational Rules</title><p>Before the crisis, between 1984 and 2013, the average extraction flow was 29.6 m<sup>3</sup>∙s<sup>−1</sup>, ranging between 18.6 and 32.8 m<sup>3</sup>∙s<sup>−1</sup>, and the average outflow was 10.6 m<sup>3</sup>∙s<sup>−1</sup>, varying between 0.8 and 91.7 m<sup>3</sup>∙s<sup>−1</sup>. This highest outflow value was</p><p>observed in 2010, during an extreme event when Cantareira System reservoirs were not able to support the large inflow. During the period of volume pumping utilization, from July 2014 to December 2015, the average extraction flow was 16.2 m<sup>3</sup>∙s<sup>−1</sup>, varying between 13.4 m<sup>3</sup>∙s<sup>−1</sup> and 21.8 m<sup>3</sup>∙s<sup>−1</sup>, and the average outflow was 2.2 m<sup>3</sup>∙s<sup>−1</sup>, varying between 0.5 and 4.2 m<sup>3</sup>∙s<sup>−1</sup> (<xref ref-type="fig" rid="fig2">Figure 2</xref>).</p><p>Since May 2017, ANA and DAEE implement a new water extraction rule via the ANA/DAEE resolution N˚925 [<xref ref-type="bibr" rid="scirp.99361-ref16">16</xref>]. This resolution defines storage levels operational ranges and their respective maximum extraction flow, as shown in <xref ref-type="table" rid="table1">Table 1</xref>, ranging from a minimum of 15.5 m<sup>3</sup>∙s<sup>−1</sup> to a maximum of 33.0 m<sup>3</sup>∙s<sup>−1</sup>, replacing a previous rule in which the maximum extraction had been set in 33 m<sup>3</sup>∙s<sup>−1</sup> for any storage condition. An evaluation on the water allocation rules in the Cantareira System, its drivers and how it has reflected in the regional sharing of the benefits can be found in [<xref ref-type="bibr" rid="scirp.99361-ref22">22</xref>].</p><p>In order to increase water supply, ANA implemented the water reallocation between the Cantareira System and the Jaguari reservoir from the Para&#237;ba do Sul River Basin (PS) via resolution ANA N˚1931 [<xref ref-type="bibr" rid="scirp.99361-ref17">17</xref>], published in October 2017, and effectively started in April 2018. The water reallocation can operate both ways, although this process has just been favored Cantareira System. According to this resolution, the reallocation flow annual mean must be up to 5.13 m<sup>3</sup>∙s<sup>−1</sup>, when the Cantareira storage level is lower than 60%, and the Jaguari-PS reservoir water level is ranging from 603 to 623 m [<xref ref-type="bibr" rid="scirp.99361-ref17">17</xref>].</p></sec><sec id="s2_3"><title>2.3. Control Simulation of the Storage Level Dynamic</title><p>The System storage level was simulated by a Reservoir Simulation module coupled to Probability Distributed Model—PDM implemented at the Brazilian National Center for Monitoring and Early Warning of Natural Disasters (CEMADEN), PDM-CEMADEN [<xref ref-type="bibr" rid="scirp.99361-ref23">23</xref>]. The PDM is a conceptual hydrological model that transforms rainfall and potential evapotranspiration into streamflow and actual evapotranspiration. A detailed description of the PDM-Cemaden model can be found in [<xref ref-type="bibr" rid="scirp.99361-ref23">23</xref>] [<xref ref-type="bibr" rid="scirp.99361-ref24">24</xref>] [<xref ref-type="bibr" rid="scirp.99361-ref25">25</xref>]. Although this model is generally implemented to simulate the streamflow, in this study just the Storage Simulation</p><table-wrap id="table1" ><label><xref ref-type="table" rid="table1">Table 1</xref></label><caption><title> Cantareira System (CS) reservoirs operational ranges and their respective maximum extraction flows for the MASP water supply (resolution ANA/DAEE N˚925) and the water reallocation mean flow from Para&#237;ba do Sul river basin (PS) (resolution ANA N˚1931)</title></caption><table><tbody><thead><tr><th align="center" valign="middle" >Storage Level (%)</th><th align="center" valign="middle" >Condition</th><th align="center" valign="middle" >Maximum extraction flow (m<sup>3</sup>∙s<sup>−1</sup>)</th><th align="center" valign="middle" >Water reallocation mean flow PS-CS (m<sup>3</sup>∙s<sup>−1</sup>)</th></tr></thead><tr><td align="center" valign="middle" >&gt;60%</td><td align="center" valign="middle" >Normal</td><td align="center" valign="middle" >33.0</td><td align="center" valign="middle" >-</td></tr><tr><td align="center" valign="middle" >60%&lt; and &gt;40%</td><td align="center" valign="middle" >Attention</td><td align="center" valign="middle" >31.0</td><td align="center" valign="middle" >5.13</td></tr><tr><td align="center" valign="middle" >40%&lt; and &gt;30%</td><td align="center" valign="middle" >Alert</td><td align="center" valign="middle" >27.0</td><td align="center" valign="middle" >5.13</td></tr><tr><td align="center" valign="middle" >30%&lt; and &gt;20%</td><td align="center" valign="middle" >Restriction</td><td align="center" valign="middle" >23.0</td><td align="center" valign="middle" >5.13</td></tr><tr><td align="center" valign="middle" >&lt;20%</td><td align="center" valign="middle" >Special</td><td align="center" valign="middle" >15.5</td><td align="center" valign="middle" >5.13</td></tr></tbody></table></table-wrap><p>module for the superficial reservoir was used. Therefore, the observed in- and out-flow were used for the control simulation and also for all scenarios’ simulations (<xref ref-type="fig" rid="fig4">Figure 4</xref>). The module was run starting in 2004 for the control simulation and used the operationally-practiced extraction flow data (<xref ref-type="fig" rid="fig4">Figure 4</xref>(a)). Flow data were obtained from SABESP website (http://mananciais.sabesp.com.br/Home).</p></sec><sec id="s2_4"><title>2.4. Storage Level Dynamic Scenarios Simulation</title><p>In order to assess the Cantareira System storage level dynamics impacts, this paper addressed four scenarios based on water mitigation policies, varying the policies implementation starting date, the extraction flow range and including the water reallocation variable, as described in <xref ref-type="table" rid="table2">Table 2</xref>.</p><p>The scenario I considers water extraction following the Resolution ANA/DAEE N˚925. Scenario II considers, in addition to Scenario I, the water reallocation from Jaguari-PS, following the Resolution ANA N˚1931. Therefore, Scenario II</p><table-wrap id="table2" ><label><xref ref-type="table" rid="table2">Table 2</xref></label><caption><title> Description of the four simulation based on water mitigation policies scenarios adressed in this paper</title></caption><table><tbody><thead><tr><th align="center" valign="middle" >Scenario</th><th align="center" valign="middle" >Description</th><th align="center" valign="middle" >Time period of policies implementation</th></tr></thead><tr><td align="center" valign="middle" >I</td><td align="center" valign="middle" >Maximum water extraction for storage level operational ranges from the resolution ANA/DAEE N˚925 (shown in <xref ref-type="table" rid="table1">Table 1</xref>)</td><td align="center" valign="middle" >2014-2018</td></tr><tr><td align="center" valign="middle" >II</td><td align="center" valign="middle" >Scenario (I) including the water reallocation flow, according to the resolution ANA N˚1931 (shown in <xref ref-type="table" rid="table1">Table 1</xref>)</td><td align="center" valign="middle" >2014-2018</td></tr><tr><td align="center" valign="middle" >III</td><td align="center" valign="middle" >Same extraction rule as Scenario I, but assuming an earlier policy implementation</td><td align="center" valign="middle" >2005-2018</td></tr><tr><td align="center" valign="middle" >IV</td><td align="center" valign="middle" >Intermediate water extraction flow for storagelevel operational ranges from the Resolution ANA/DAEE N˚925: Normal range: 32 m<sup>3</sup>∙s<sup>−1</sup>; Attention range: 29 m<sup>3</sup>∙s<sup>−1</sup>; Alert range: 25 m<sup>3</sup>∙s<sup>−1</sup>; Restriction range: 19.5 m<sup>3</sup>∙s<sup>−1</sup>; Special range: 15.5 m<sup>3</sup>∙s<sup>−1</sup>.</td><td align="center" valign="middle" >2005-2018</td></tr></tbody></table></table-wrap><p>is more optimistic than Scenario I. These scenario simulations start the policies implementation in March 2014, at the end of the rainy season in the Southeast Brazil and when concern with the critical reservoir storage was raised. In March 2014, the reservoir storage was only 17% of the active volume, and the rainfall, from Oct 2013 to Feb 2014, accumulated a climatological rainfall deficit of 61%. Both conditions could be enough to initiate a water use restriction. The water use restriction management in real-time during an HD is challenging, such that making decision to restrict reservoir water extraction in advance would be difficult.</p><p>Scenarios III and IV were simulated for a longer policy implementation time period, starting one year after the permit renewal (2004). Scenario III used the same water extraction rule as Scenario I, although under different initial conditions. Since the Resolution ANA/DAEE N˚925 was implemented, in 2017, operationally-practiced water extraction has been below the maximum allowed, therefore, a new scenario (Scenario IV) was elaborated to assess this impact on storage level dynamics. Scenario IV used intermediate extraction flow for each storage level operational range from the Resolution, as described in <xref ref-type="table" rid="table2">Table 2</xref>.</p><p>In scenarios simulations, pumping volume utilization just would begin when Cantareira System storage level is below 0%, although it starts prior to this level in operationally-practice.</p></sec></sec><sec id="s3"><title>3. Results</title><sec id="s3_1"><title>3.1. Control Simulation</title><p>The control simulation performance indicators (Nash-Sutcliff = 0.98; Mean Absolute Error = 137.8; Root Mean Square Error = 7.49) showed the great ability of the PDM-CEMADEN model in simulating the storage level dynamics.</p></sec><sec id="s3_2"><title>3.2. Scenarios Simulation</title><p>The Cantareira System 4-year simulations of storage level and extraction flow under Scenarios I and II are shown, respectively, in <xref ref-type="fig" rid="fig5">Figure 5</xref> and in <xref ref-type="fig" rid="fig6">Figure 6</xref>. In Scenario I, the model estimated that, during the 2014-2016 water crisis, pumping volume would have been extracted during two shorter periods than operationally-practiced (observed), lasting 239 and 141 days, respectively, interspersed by 34 days of active volume utilization (around June, 2015), totalizing 380 days of pumping volume utilization (<xref ref-type="table" rid="table3">Table 3</xref>). Regarding pumping volume storage, the minimum estimated level would have been −11% (−110 hm&#179;), thus the pumping volume first quota would have been utilized, and the total pumped water volume would have reached 509 hm&#179;, equivalent to 52% of active volume (<xref ref-type="table" rid="table3">Table 3</xref>).</p><p>Cantareira System level storage under Scenario II showed that pumping volume utilization would not have been necessary during the 2014-2016 water crisis. The minimum estimated storage level would have been reached in January 2015, representing 3.9% (38 hm<sup>3</sup>) of the active storage. This more optimistic storage dynamics is attributed to water reallocation, which would contribute to 243 hm<sup>3</sup>, equivalent to 25% of the active volume, during the period between July 2014 to December 2015.</p><table-wrap id="table3" ><label><xref ref-type="table" rid="table3">Table 3</xref></label><caption><title> Cantareira System pumping volume utilization data: operationally-practiced and based on 4 water mitigation policies scenarios</title></caption><table><tbody><thead><tr><th align="center" valign="middle"  rowspan="2"  >Scenario</th><th align="center" valign="middle"  colspan="4"  >Pumping volume utilization</th></tr></thead><tr><td align="center" valign="middle" >Period</td><td align="center" valign="middle" >Days</td><td align="center" valign="middle" >Extraction flow mean (m<sup>3</sup>∙s<sup>−1</sup>)</td><td align="center" valign="middle" >Total extracted volume hm<sup>3</sup> and (%)<sup>a</sup></td></tr><tr><td align="center" valign="middle" >Practiced</td><td align="center" valign="middle" >Jul 12<sup>th</sup>, 2014 to Dec 30<sup>th</sup>, 2015</td><td align="center" valign="middle" >537</td><td align="center" valign="middle" >16.2</td><td align="center" valign="middle" >745 (76%)</td></tr><tr><td align="center" valign="middle" >I</td><td align="center" valign="middle" >Oct 06<sup>th</sup>, 2014-Jun 01<sup>st</sup>, 2015 and Jul 06<sup>th</sup>, 2015-Nov 23<sup>th</sup>, 2015</td><td align="center" valign="middle" >380</td><td align="center" valign="middle" >15.5</td><td align="center" valign="middle" >509 (52%)</td></tr><tr><td align="center" valign="middle" >II</td><td align="center" valign="middle" >Not used</td><td align="center" valign="middle" >0</td><td align="center" valign="middle" >-</td><td align="center" valign="middle" >0 (0%)</td></tr><tr><td align="center" valign="middle" >III</td><td align="center" valign="middle" >Oct 23<sup>th</sup>, 2014-Mar 01<sup>st</sup>, 2015</td><td align="center" valign="middle" >130</td><td align="center" valign="middle" >15.5</td><td align="center" valign="middle" >174 (18%)</td></tr><tr><td align="center" valign="middle" >IV</td><td align="center" valign="middle" >Jan 22<sup>th</sup>, 2015-Feb 07<sup>th</sup>, 2015</td><td align="center" valign="middle" >17</td><td align="center" valign="middle" >15.5</td><td align="center" valign="middle" >23 (2%)</td></tr></tbody></table></table-wrap><p><sup>a</sup>In relation to the active volume, i.e. 982 hm<sup>3</sup>.</p><p>The Cantareira System 13-year simulations under Scenarios III and IV are shown in <xref ref-type="fig" rid="fig7">Figure 7</xref>. In Scenario III, the model estimated that pumping volume would have been extracted during 130 days, a number of days lower than the</p><p>operationally-practiced. The minimum estimated storage level would have been −6% (−61 hm<sup>3</sup>), reached in February, 2015, therefore, just the pumping volume first quota would have been utilized. Total pumped water volume would have reached 174 hm<sup>3</sup>, equivalent to 18% of the active volume. In Scenario IV, the storage level simulated showed that pumping volume would still have been extracted, although just during 17 days. Regarding pumping volume utilization, the estimated minimum storage level would have been −0.4% (−4 hm<sup>3</sup>), reached in February, 2015 and the total pumped volume would reach 23 hm<sup>3</sup>, equivalent to 2% of the active volume.</p></sec></sec><sec id="s4"><title>4. Discussion</title><p>The scenario simulations pointed out that extraction flow mean would be 15.5 m<sup>3</sup>∙s<sup>−1</sup> (except for Scenario II) during the period July 2014 to December 2015, a lower value than real practiced (16.2 m<sup>3</sup>∙s<sup>−1</sup>) (<xref ref-type="table" rid="table3">Table 3</xref>). Regarding this time period, the simulated storage level under Scenario II would have been higher than under other scenarios and also than observed storage level, which would have allowed a higher extraction flow mean (17.2 m<sup>3</sup>∙s<sup>−1</sup>). However, under all scenarios simulated, it would have not been possible extracting the total demand observed before the water crisis (29.6 m<sup>3</sup>∙s<sup>−1</sup>).</p><p>After the 2014-2016 water crisis, all four simulations showed that storage level would have been lower than the observed (<xref ref-type="fig" rid="fig5">Figure 5</xref> and <xref ref-type="fig" rid="fig7">Figure 7</xref>) which could be explained by higher extraction flow mean estimated. Successively, higher extraction from Cantareira System could have lessened impact the MASP water supply.</p><p>Scenarios I and III followed the same extraction rule as established on Resolution ANA/DAEE N˚925, however, their initial storage level was 17% in 2014 and 47% in 2005, respectively, leading to different storage levels dynamics during the 2014-2016 water crisis. Pumping volume would have been less utilized (250 days; 335 hm<sup>3</sup>) had the Resolution been implemented years (Scenario III) earlier compared to months (Scenario I) prior to the crisis. The early implementation of this policy could have reduced the water system vulnerability.</p><p>Simulations results showed that Cantareira System total storage (active and pumping volume) would not be enough to meet the total MASP water demand observed prior to the crisis, even considering the water reallocation from another basin. Therefore, other water systems located in the MASP should have been utilized in order to supply 29.6 m<sup>3</sup>∙s<sup>−1</sup>. The Alto Tiet&#234; System utilization increased during the 2014-2016 water crisis, which also led this system into a crisis. Wherefore, requiring more inflow from this system could not be possible.</p><p>The water reallocation over longer distances probably increases the financial resources utilization and energy’s needs. For that reason, the sustainability of the water supply project for crisis management should be carefully analyzed. In addition, prolonged and severe droughts usually affect large territorial areas, which challenges the regional water management, as occurred in the HD addressed in this paper. Several reservoirs in the Southeastern Brazil faced an unfavored condition during 2014-2016 [<xref ref-type="bibr" rid="scirp.99361-ref8">8</xref>] [<xref ref-type="bibr" rid="scirp.99361-ref13">13</xref>] [<xref ref-type="bibr" rid="scirp.99361-ref26">26</xref>]. The Jaguari-PS reservoir recorded the historical minimum value of 3% of active volume in December 2014 [<xref ref-type="bibr" rid="scirp.99361-ref27">27</xref>], and the water reallocation to Cantareira System was questionable. Therefore, the assumption that water problems there will finally be resolved with the construction of more infrastructure could fail due to growing human needs and climate change [<xref ref-type="bibr" rid="scirp.99361-ref28">28</xref>].</p><p>Water systems management can be influenced by various resources: financial, natural and adaptative capacity. Appropriate operational practices can substantially reduce these resources impacts on water systems under water scarcity and drought [<xref ref-type="bibr" rid="scirp.99361-ref29">29</xref>]. Each hypothetical scenario in this study presents alternatives for the Cantareira System operation in order to support the management of possible future water scarcity.</p></sec><sec id="s5"><title>5. Conclusions</title><p>In 2014-2016 Cantareira System faced an unprecedented HD that overwhelmed the MASP water supply. This study assessed Cantareira System mitigation policies (implemented after the crisis) impacts on the system storage level dynamics. Impacts were evaluated under four water mitigation policies scenarios by varying policies implementation starting date, extraction flow magnitude range and including the water reallocation variable.</p><p>Answering the three addressed questions in this paper, 1) the storage level simulation under Scenario I pointed out that Cantareira System pumping volume would have been utilized to supply the MASP water system, which had the extraction rule (Resolution ANA/DAEE N˚925) been implemented in March 2014. Moreover, Scenario II results showed that the Cantareira System active volume would have been able to supply the MASP water system, provided that surplus input from another basin (Resolution ANA N˚1931) had been used since March 2014. 2) Scenario III and IV results showed that pumping volume would still have been utilized during a period of time shorter than operationally practiced, had the policies already been implemented some years prior to the crisis. 3) The water mitigation policies implementation supports the water supply via reservoirs storage management. Further, despite water reallocation from Scenario II suffices to pumping volume saving, the MASP water demand prior to the crisis is not fully met. Considering the observed negative impacts of precipitation and inflow in recent years and S&#227;o Paulo’s economic contribution to the country’s Gross Domestic Product (GDP), a strategic planning approach to water security is necessary, via the mismatches evaluation between water supply and demand. The proposed flow magnitude reduction addressed in Scenario IV could support policy makers to improve the Resolution ANA/DAEE N˚925 aiming to further decrease the Cantareira System vulnerability.</p><p>Lastly, policies and practices could shift from water-supply to demand-management approaches contributing to water security improvement. An implemented practice example was the temporarily incentive to reduce consumption applied during the crisis. Actions should also address reducing leakages from water supply distribution network and enhance sewage collection and treatment and also the public participation for future water management.</p></sec><sec id="s6"><title>Acknowledgements</title><p>The authors are very grateful to Diogo Amore, hydrological scientist, for text and English revising. The authors are grateful for the finantial support from CEMADEN Institutional Capacity Building Program, (CNPq Process N˚454809/2015-8); from CAPES Pro Alertas CEPED/USP (88887.091743/2014-01), from the Consolidation of research and technological development work, Process N˚401897/2014-1, and from the National Institute of Science and Technology for Climate Change Phase 2 under: CNPq Grant 465501/2014-1, FAPESP Grant 2014/50848-9. The authors gratefully acknowledge the S&#227;o Paulo State Basic Sanitation Company (Sabesp) and the National Water Agency (ANA) for providing the data required in this study.</p></sec><sec id="s7"><title>Conflicts of Interest</title><p>The authors declare no conflicts of interest regarding the publication of this paper.</p></sec><sec id="s8"><title>Data Availability Statement</title><p>The data that support the findings of this study are available from the corresponding author upon reasonable request.</p></sec><sec id="s9"><title>Cite this paper</title><p>Deusdar&#225;-Leal, K.R., Cuartas, L.A., Zhang, R., Mohor, G.S., de Castro Carvalho, L.V., Nobre, C.A., Mendiondo, E.M., Broedel, E., Seluchi, M.E. and dos Santos Alval&#225;, R.C. (2020) Implications of the New Operational Rules for Cantareira Water System: Re-Reading the 2014-2016 Water Crisis. Journal of Water Resource and Protection, 12, 261-274. https://doi.org/10.4236/jwarp.2020.124016</p></sec></body><back><ref-list><title>References</title><ref id="scirp.99361-ref1"><label>1</label><mixed-citation publication-type="other" xlink:type="simple">Van Loon, A.F. (2015) Hydrological Drought Explained. Wiley Interdisciplinary Reviews: Water, 2, 359-392. https://doi.org/10.1002/wat2.1085</mixed-citation></ref><ref id="scirp.99361-ref2"><label>2</label><mixed-citation publication-type="other" xlink:type="simple">Nobre, C.A., Marengo, J.A., Seluchi, M.E., Cuartas, A. and Alves, L.M. (2016) Some Characteristics and Impacts of the Drought and Water Crisis in Southeastern Brazil during 2014 and 2015. Journal of Water Resource and Protection, 8, 252-262. https://doi.org/10.4236/jwarp.2016.82022</mixed-citation></ref><ref id="scirp.99361-ref3"><label>3</label><mixed-citation publication-type="other" xlink:type="simple">Coelho, C.A.S., de Oliveira, C.P., Ambrizzi, T., Reboita, M.S., Carpenedo, C.B., Campos, J.L.P.S., et al. (2016) The 2014 Southeast Brazil Austral Summer Drought: Regional Scale Mechanisms and Teleconnections. Climate Dynamics, 46, 3737-3752. https://doi.org/10.1007/s00382-015-2800-1</mixed-citation></ref><ref id="scirp.99361-ref4"><label>4</label><mixed-citation publication-type="other" xlink:type="simple">Coelho, C.A.S., Cardoso, D.H.F. and Firpo, M.A.F. (2016) Precipitation Diagnostics of an Exceptionally Dry Event in Sao Paulo, Brazil. Theoretical and Applied Climatology, 125, 769-784. https://doi.org/10.1007/s00704-015-1540-9</mixed-citation></ref><ref id="scirp.99361-ref5"><label>5</label><mixed-citation publication-type="other" xlink:type="simple">Millington, N. (2018) Producing Water Scarcity in Sao Paulo, Brazil: The 2014-2015 Water Crisis and the Binding Politics of Infrastructure. Political Geography, 65, 26-34. https://doi.org/10.1016/j.polgeo.2018.04.007</mixed-citation></ref><ref id="scirp.99361-ref6"><label>6</label><mixed-citation publication-type="other" xlink:type="simple">Otto, F.E.L., Coelho, C.A.S., King, A., de Perez, E.C., Wada, Y., van Oldenborgh, G.J., et al. (2015) Factors Other than Climate Change, Main Drivers of 2014/15 Water Shortage in Southeast Brazil. Bulletin of the American Meteorological Society, 96, S35-S40. https://doi.org/10.1175/BAMS-D-15-00120.1</mixed-citation></ref><ref id="scirp.99361-ref7"><label>7</label><mixed-citation publication-type="other" xlink:type="simple">Dobrovolski, R. and Rattis, L. (2015) Water Collapse in Brazil: The Danger of Relying on What You Neglect. Natureza e Conservacao, 13, 80-83. https://doi.org/10.1016/j.ncon.2015.03.006</mixed-citation></ref><ref id="scirp.99361-ref8"><label>8</label><mixed-citation publication-type="other" xlink:type="simple">Hanbury, S. (2015) Brazil’s Worst Drought in History Prompts Protests and Blackouts. The Guardian. https://www.theguardian.com/world/2015/jan/23/brazil-worst-drought-history</mixed-citation></ref><ref id="scirp.99361-ref9"><label>9</label><mixed-citation publication-type="other" xlink:type="simple">Taffarello, D., Samprogna Mohor, G., do Carmo Calijuri, M. and Mendiondo, E.M. (2016) Field Investigations of the 2013-14 Drought through Quali-Quantitative Freshwater Monitoring at the Headwaters of the Cantareira System, Brazil. Water International, 41, 776-800. https://doi.org/10.1080/02508060.2016.1188352</mixed-citation></ref><ref id="scirp.99361-ref10"><label>10</label><mixed-citation publication-type="other" xlink:type="simple">Cunha, A.P.M.A., Zeri, M., Deusdará Leal, K., Costa, L., Cuartas, L.A., Tomasella, J., et al. (2019) Extreme Drought Events over Brazil from 2011 to 2019. Atmosphere, 10, 20. https://doi.org/10.3390/atmos10110642</mixed-citation></ref><ref id="scirp.99361-ref11"><label>11</label><mixed-citation publication-type="book" xlink:type="simple">Lund, J. and Medellín-Azuara, J. (2018) California: Water Security from Infrastructure, Institutions, and the Global Economy. In: Council, W.W., Ed., Global Water Security: Water Resources Development and Management, Springer, Singapore, 267-279. https://doi.org/10.1007/978-981-10-7913-9_11</mixed-citation></ref><ref id="scirp.99361-ref12"><label>12</label><mixed-citation publication-type="other" xlink:type="simple">Whately, M. and Cunha, P. (2006) Cantareira 2006: Um olhar sobre o maior manancial de água da Regiao Metropolitana de Sao Paulo. Resultados do Diagnóstico Socioambiental Participativo do Sistema Cantareira. Instituto socioambiental, Sao Paulo. (In Portuguese)</mixed-citation></ref><ref id="scirp.99361-ref13"><label>13</label><mixed-citation publication-type="other" xlink:type="simple">Kelman, J. (2015) Water Supply to the Two Largest Brazilian Metropolitan Regions. Aquatic Procedia, 5, 13-21. https://doi.org/10.1016/j.aqpro.2015.10.004</mixed-citation></ref><ref id="scirp.99361-ref14"><label>14</label><mixed-citation publication-type="other" xlink:type="simple">Sabesp (2014) Sistema Cantareira—Plano de Contingência II. Acoes Contingenciais e Resultados. 75. (In Portuguese)</mixed-citation></ref><ref id="scirp.99361-ref15"><label>15</label><mixed-citation publication-type="other" xlink:type="simple">Leite, F., Italiani, R. and Burgarelli, R. (2015) O Passado e o Futuro do Cantareira. Estadao. (In Portuguese) https://infograficos.estadao.com.br/especiais/passado-futuro-cantareira</mixed-citation></ref><ref id="scirp.99361-ref16"><label>16</label><mixed-citation publication-type="other" xlink:type="simple">ANA and DAEE (2017) Resolucao Conjunta ANA/DAEE N&amp;deg;925, de 29 de Maio de 2017. (In Portuguese) http://arquivos.ana.gov.br/resolucoes/2017/925-2017.pdf</mixed-citation></ref><ref id="scirp.99361-ref17"><label>17</label><mixed-citation publication-type="other" xlink:type="simple">ANA (2017) Resolucao N&amp;deg;1.931, de 30 de Outubro de 2017. (In Portuguese) http://www.ceivap.org.br/resolucoes/ana/2017/1931-2017.pdf</mixed-citation></ref><ref id="scirp.99361-ref18"><label>18</label><mixed-citation publication-type="other" xlink:type="simple">Vieira, B.M. (2018) Haverá novo racionamento em SP? Veja o que dizem especialistas e a Sabesp. G1. (In Portuguese) https://g1.globo.com/sp/sao-paulo/noticia/2018/08/02/havera-novo-racionamento-em-sp-veja-o-que-dizem-especialistas-e-a-sabesp.ghtml</mixed-citation></ref><ref id="scirp.99361-ref19"><label>19</label><mixed-citation publication-type="other" xlink:type="simple">Vieira, B.M. (2018) Volume do Sistema Cantareira tem ligeira recuperacao, mas situacao ainda é preocupante. G1. (In Portuguese) https://g1.globo.com/sp/sao-paulo/noticia/2018/11/26/volume-do-sistema-cantareira-tem-ligeira-recuperacao-mas-situacao-ainda-e-preocupante-diz-especialista.ghtml</mixed-citation></ref><ref id="scirp.99361-ref20"><label>20</label><mixed-citation publication-type="other" xlink:type="simple">Cunha, A.P., Marengo, J.A., Alvalá, R.C., Deusdará-Leal, K.R., Cuartas, L.A., Seluchi, M.E., et al. (2019) Secas e Seus Impactos no Brasil em 2018. Sao José dos Campos-SP. (In Portuguese) http://www.cemaden.gov.br/wp-content/uploads/2019/01/Boletim_Anual_SECAS_f.pdf</mixed-citation></ref><ref id="scirp.99361-ref21"><label>21</label><mixed-citation publication-type="other" xlink:type="simple">Cemaden (2019) Situacao Atual e Projecao Hidrológica. Sao José dos Campos-SP. (In Portuguese)https://www.cemaden.gov.br/categoria/monitoramento/monitoramento-hidrologico/relatoriocantareira</mixed-citation></ref><ref id="scirp.99361-ref22"><label>22</label><mixed-citation publication-type="other" xlink:type="simple">Leao, R.S. and De Stefano, L. (2019) Making Concrete Flexible: Adapting the Operating Rules of the Cantareira Water System (Sao Paulo, Brazil). Water Security, 7, Article ID: 100032. https://doi.org/10.1016/j.wasec.2019.100032</mixed-citation></ref><ref id="scirp.99361-ref23"><label>23</label><mixed-citation publication-type="other" xlink:type="simple">Zhang, R., Cuartas, L.A., Carvalho, L.V.C., Leal, K.R.D., Mendiondo, E.M., Abe, N., et al. (2018) Season-Based Rainfall-Runoff Modelling Using the Probability-Distributed Model (PDM) for Large Basins in Southeastern Brazil. Hydrological Processes, 32, 2217-2230. https://doi.org/10.1002/hyp.13154</mixed-citation></ref><ref id="scirp.99361-ref24"><label>24</label><mixed-citation publication-type="other" xlink:type="simple">Moore, R.J. (2007) The PDM Rainfall-Runoff Model. Hydrology and Earth System Sciences, 11, 483-499. https://doi.org/10.5194/hess-11-483-2007</mixed-citation></ref><ref id="scirp.99361-ref25"><label>25</label><mixed-citation publication-type="other" xlink:type="simple">Moore, R.J. and Bell, V.A. (2002) Incorporation of Groundwater Losses and Well Level Data in Rainfall-Runoff Models Illustrated Using the PDM. Hydrology and Earth System Sciences, 6, 25-38. https://doi.org/10.5194/hess-6-25-2002</mixed-citation></ref><ref id="scirp.99361-ref26"><label>26</label><mixed-citation publication-type="other" xlink:type="simple">Escobar, H. (2015) Drought Triggers Alarms in Brazil’s Biggest Metropolis. Science, 347, 812-812. https://doi.org/10.1126/science.347.6224.812</mixed-citation></ref><ref id="scirp.99361-ref27"><label>27</label><mixed-citation publication-type="other" xlink:type="simple">ANA (2019) Boletim de monitoramento dos reservatórios do sistema hidráulico do rio Paraíba do Sul. https://www.ana.gov.br/sala-de-situacao/paraiba-do-sul/boletins/mensal/05-boletim_monitoramento_reservatorios_pb_sul_2019_maio.pdf/view</mixed-citation></ref><ref id="scirp.99361-ref28"><label>28</label><mixed-citation publication-type="other" xlink:type="simple">Vorosmarty, C.J., Green, P., Salisbury, J. and Lammers, R.B. (2000) Global Water Resources: Vulnerability from Climate Change and Population Growth Contemporary Population Relative to Demand per Discharge. Science, 289, 284-288. https://doi.org/10.1126/science.289.5477.284</mixed-citation></ref><ref id="scirp.99361-ref29"><label>29</label><mixed-citation publication-type="other" xlink:type="simple">Martin-Ortega, J., González-Eguino, M. and Markandya, A. (2012) The Costs of Drought: The 2007/2008 Case of Barcelona. Water Policy, 14, 539-560. https://doi.org/10.2166/wp.2011.121</mixed-citation></ref></ref-list></back></article>