# Urban Metabolic Analysis of a Food-Water-Energy System for Sustainable Resources Management

^{1}

^{2}

^{3}

^{*}

Next Article in Journal

Next Article in Special Issue

Next Article in Special Issue

Previous Article in Journal

Previous Article in Special Issue

Previous Article in Special Issue

Department of Bioenvironmental Systems Engineering, National Taiwan University, No. 1, Sec. 4, Roosevelt Road, Taipei 10617, Taiwan

Department of Urban Planning, National Cheng Kung University, No. 1, University Road, Tainan City 70101, Taiwan

Department of Agricultural Economics, National Taiwan University, No. 1, Sec. 4, Roosevelt Road, Taipei 10617, Taiwan

Author to whom correspondence should be addressed.

Received: 9 December 2018 / Revised: 24 December 2018 / Accepted: 25 December 2018 / Published: 30 December 2018

(This article belongs to the Special Issue Land management accounting inter-correlations among soil erosion risk, cropping suitability and environmental footprint)

Urban metabolism analyzes the supply and consumption of nutrition, material, energy, and other resources within cities. Food, water, and energy are critical resources for the human society and have complicated cooperative/competitive influences on each other. The management of interactive resources is essential for supply chain analysis. This research analyzes the food-water-energy system of urban metabolism for sustainable resources management. A system dynamics model is established to investigate the urban metabolism of food, water, and energy resources. This study conducts a case study of Shihmen Reservoir system, hydropower generation, paddy rice irrigation of Taoyuan and Shihmen Irrigation Associations, and water consumption in Taoyuan, New Taipei, and Hsinchu cities. The interactive influence of the food-water-energy nexus is quantified in this study; the uncertainty analysis of food, water, and energy nexus is presented.

This research aims to analyze food, water, and energy resources of urban metabolism for sustainable resources management. The interactions between the natural environment and humans in cities are complex systems. Urban metabolism is a concept simulating the production, distribution, intake, digestion, and absorption of nutrition, material, energy, and other resources within cities [1,2,3,4,5,6,7,8]. Food, water, and energy are scarce natural resources for human society. The cooperative and competitive relationships are important topics for food-water-energy resources management. The cooperative and competitive relationships of food, water, and energy attract lots of attention; they needed to be quantified [9,10,11,12,13,14,15,16,17].

Previous studies have investigated urban metabolism using life cycle, material flow, meta-analysis, input-output analysis, and simulation tools. The nexus of water, energy, material, carbon dioxide emission, ecological system, climate, and human society has been explored and discussed. Chen and Chen (2017) analyzed the coupling of carbon and energy flows in the cities of Beijing and Issaquah. The research used meta-analysis and nexus modeling to evaluate the energy-carbon nexus [18]. Chen et al. (2018) conducted multi-regional physical input-output analysis for the water-energy nexus in Hong Kong and South China [19]. De Stercke et al. (2018) investigated end-use interactions of the urban water-energy nexus to reduce pressure of water, electricity, and carbon dioxide emissions. The dynamic interactions between London’s water and energy systems was simulated [20]. Facchini et al. (2017) explored the energy metabolisms in 27 of the world’s megacities, including mobile and stationary energy consumption patterns, fuels used, and the end-use and electricity generation mix [21]. Martinez-Hernandez et al. (2017) analyzed water-energy-food and ecosystem interactions. The concomitant Nexus Simulation System (NexSym) is presented to assess a synergistic nexus system in a UK eco-town [22]. Ravalde and Keirstead (2017) compared performance metrics for multi-resource systems of urban metabolism [23]. Schlör et al. (2018) investigated interactions among nature and society of an energy-mineral-society nexus. The nexus is analyzed by the social life cycle assessment model with the scenario based on the Human Development Index [24]. Wang and Chen (2016) addressed the energy–water nexus of urban agglomeration based on multiregional input–output tables and ecological network analysis in the Beijing–Tianjin–Hebei region [25]. Yang et al. (2018) explored the environmental pressures in urban sectors using a energy-water-carbon nexus perspective. This study employed environmental input–output analysis to estimate sectoral embodied energy, water, and CO_{2} in Beijing and Shanghai [26]. Zhang et al. (2018) utilized statistics, equilibrium, econometric, ecological, life cycle, system dynamics, agent-based, and integrated index methods to explore the water-energy-food nexus [27].

Related studies of the urban metabolism involve large-scale, conceptual, life cycle, or simulation modeling analysis. This research aims to investigate the cooperative/competitive interaction of the water-energy-food nexus. Furthermore, quantitative analysis of the interactive influence of the food-water-energy nexus is presented in this study. A case study in a system with water supply, water consumption, hydropower generation, and paddy rice growing is also conducted. Figure 1 depicts the cooperative/competitive relationship in the food-water-energy nexus for urban metabolic analysis. The Integrated management and tradeoff of water, food, and energy is analyzed in this study.

In this research, system dynamics models are established to analyze the urban metabolic food-water-energy systems. System dynamics was developed in 1950; the simulation method is widely applied in industrial, engineering, natural science, and public policy fields [28,29,30,31,32,33,34,35,36]. Furthermore, system dynamics can be used to examine optimal state trajectories and tradeoffs of systems by solving dynamic differential system equations [37,38,39,40]. Applications of system dynamics and optimization include mechanical, aerospace, chemical engineering, economics, management sciences, and other related fields [38,39,40,41,42].

The system dynamics models are used to examine cooperative/competitive behaviors, including the water supply, water consumption, hydropower generation, and food crops growing in food-water-energy systems. The urban metabolism simulation model is established on NetLogo platform. NetLogo system is a popular and powerful system dynamics simulation and agent-based interactive simulation platform. NetLogo models utilize system dynamics, patch, turtle, plot, and monitor agents to analyze the system performance and agent interactions. In our analysis, NetLogo uses patches and turtle agents to simulate water, paddy rice, hydropower production, and demand/supply interactions. Furthermore, monitor agents collect and display dynamics changes of food, water, and energy nexus systems.

Compared with related studies, the significant contribution of this research includes the establishment of a framework for analyzing food-water-energy nexus, quantifying competition/cooperation among resources, and integrated resources management. In addition, this research constructs system dynamics models for participants in food-water-energy systems. Accordingly, the system dynamics model is established on the NetLogo platform and then equilibrium solutions of competition/cooperation are determined. Moreover, a case study of competition and cooperation among resources is conducted and discussed. This paper is organized as follows. Section 2 formulates system dynamics models to analyze the urban metabolic food-water-energy system. Section 3 presents results and discussion of a case study in Taiwan. The conclusion is addressed in Section 4.

This section begins with formulating a system dynamic model to investigate the food, water, and energy nexus. The model is capable of analyzing complex interactive component behaviors and nonlinear system performance over time. The system dynamics model of the food-water-energy system is established as follows. The water resources network model is formulated including dams, water flow, water consumption, and hydropower production. In the model, water inflow, loss, and spillage at node $\mathrm{i}$ are indicated as ${\mathrm{IN}}_{\mathrm{i}}^{}$, ${\mathrm{LOSS}}_{\mathrm{i}}^{}$, and ${\mathrm{spill}}_{\mathrm{i}}^{}$. The storage capacity of node $\mathrm{i}$ is ${\mathrm{STORAGE}}_{\mathrm{i}}^{}$. ${\mathrm{flow}}_{\mathrm{i},\mathrm{j}}^{}$ represents water flow delivered from node $\mathrm{i}$ to node $\mathrm{j}$, and ${\mathrm{supply}}_{\mathrm{i},\mathrm{k}}^{}$ indicates water supply from node $\mathrm{i}$ to demand node $\mathrm{k}$.

The model is a system dynamics simulation model with multiple objectives. Three objective functions are formulated, including food production maximization and energy production maximization. Water resources supply and demand are considered in the constraints. The objective functions are subject to constraints of water mass balance, food yield, hydropower generation, bioenergy production, reservoir operation, and non-negativity equations.

Assuming $\mathrm{R}{1}_{\mathrm{p}}^{}$ to be the food production rate of crop $\mathrm{p}$, Equation (1) calculates the food production of crop $\mathrm{p}$. $\mathrm{R}{2}_{\mathrm{p}}^{}$ represents the bioenergy production rate and total bioenergy production is estimated in Equation (2). Equation (3) calculates total hydropower generation of the water resources system. Equation (4) establishes a differential equation of water dynamic change, $\frac{{\mathrm{ds}}_{\mathrm{i}}^{}\left(t\right)}{\mathrm{dt}}$, considering total upstream inflow ($\sum}_{\mathrm{j}\in \mathrm{up}\left(\mathrm{i}\right)}^{}\left(\mathrm{R}{3}_{\mathrm{j},\mathrm{i}}^{}\times {\mathrm{flow}}_{\mathrm{j},\mathrm{i}}^{}\right)$), total downstream outflow ($\sum}_{\mathrm{j}\in \mathrm{down}\left(\mathrm{i}\right)}^{}{\mathrm{flow}}_{\mathrm{i},\mathrm{j}}^{$), loss, supply, and spill at node $\mathrm{i}$. Equation (5) decide crop growing plans of farm $\mathrm{I}$; ${\mathrm{x}}_{\mathrm{i},\mathrm{p}}^{}$ is tillage size of crop $\mathrm{p}$ with water supply from node $\mathrm{i}$ and ${\mathrm{FARM}}_{\mathrm{i}}^{}$ is farm size with water supply from node $\mathrm{i}$. In Equation (6), agricultural water demand in node $\mathrm{i}$ at time $\mathrm{t}$ is $\sum}_{\mathrm{p}}^{}\left(\mathrm{R}{4}_{\mathrm{p}}^{}\times {\mathrm{x}}_{\mathrm{i},\mathrm{p}}^{}\left(\mathrm{t}\right)\right)$ where water demand of crop $\mathrm{p}$ is $\mathrm{R}{4}_{\mathrm{p}}^{}$. Hydropower production rate at node $\mathrm{i}$ is $\mathrm{R}{5}_{\mathrm{i}}^{}$ and hydropower generation is calculated in Equation (7). The water supply of ${\mathrm{supply}}_{\mathrm{i},\mathrm{k}}^{}\left(\mathrm{t}\right)$ for residential, industrial, and agricultural demand of ${\mathrm{waterdmd}}_{\mathrm{i},\mathrm{k}}^{}\left(\mathrm{t}\right)$ in node $\mathrm{i}$ at time $\mathrm{t}$ are constructed in Equation (8). The capacity constraint of dam storage ${\mathrm{s}}_{\mathrm{i}}^{}\left(\mathrm{t}\right)$ is formulated in Equation (9). Non-negativity constraints are established in Equations (10) and (11). The multi-objective system dynamics model in Equations (1)–(11) is applied to simulate and analyze the tradeoff among food, renewable energy production, and water supply.

$${\mathrm{food}}_{\mathrm{p}}^{}={{\displaystyle \int}}_{\mathrm{t}}^{}{{\displaystyle \sum}}_{\mathrm{i}}^{}\left(\left(\mathrm{R}{1}_{\mathrm{p}}^{}\right)\left({\mathrm{x}}_{\mathrm{i},\mathrm{p}}^{}\left(\mathrm{t}\right)\right)\right)\mathrm{dt}$$

$$\mathrm{biofuel}={{\displaystyle \int}}_{\mathrm{t}}^{}{{\displaystyle \sum}}_{\mathrm{i},\mathrm{p}}^{}\left(\left(\mathrm{R}{2}_{\mathrm{p}}^{}\right)\left({\mathrm{x}}_{\mathrm{i},\mathrm{p}}^{}\left(\mathrm{t}\right)\right)\right)\mathrm{dt}$$

$$\mathrm{hydropower}={{\displaystyle \int}}_{\mathrm{t}}^{}{{\displaystyle \sum}}_{\mathrm{i}}^{}{\mathrm{hydropower}}_{\mathrm{i}}^{}\left(\mathrm{t}\right)\mathrm{dt}$$

$$\frac{{\mathrm{ds}}_{\mathrm{i}}^{}\left(\mathrm{t}\right)}{\mathrm{dt}}={\mathrm{IN}}_{\mathrm{i}}^{}+{{\displaystyle \sum}}_{\mathrm{j}\in \mathrm{up}\left(\mathrm{i}\right)}^{}\left(\left(\mathrm{R}{3}_{\mathrm{j},\mathrm{i}}^{}\right)\left({\mathrm{flow}}_{\mathrm{j},\mathrm{i}}^{}\right)\right)-{\mathrm{LOSS}}_{\mathrm{i}}^{}-{{\displaystyle \sum}}_{\mathrm{j}\in \mathrm{down}\left(\mathrm{i}\right)}^{}{\mathrm{flow}}_{\mathrm{i},\mathrm{j}}^{}-{{\displaystyle \sum}}_{\mathrm{k}}^{}{\mathrm{supply}}_{\mathrm{i},\mathrm{k}}^{}-{\mathrm{spill}}_{\mathrm{i}}^{}\text{}\forall \text{}\mathrm{i},\text{}\mathrm{t}$$

$${{\displaystyle \sum}}_{\mathrm{p}}^{}{\mathrm{x}}_{\mathrm{i},\mathrm{p}}^{}\left(\mathrm{t}\right)\le {\mathrm{FARM}}_{\mathrm{i}}^{}\text{}\forall \text{}\mathrm{i},\text{}\mathrm{t}$$

$${\mathrm{waterdmd}}_{\mathrm{i},\mathrm{k}}^{}\left(\mathrm{t}\right)={{\displaystyle \sum}}_{\mathrm{p}}^{}\left(\left(\mathrm{R}{4}_{\mathrm{p}}^{}\right)\left({\mathrm{x}}_{\mathrm{i},\mathrm{p}}^{}\left(\mathrm{t}\right)\right)\right)\text{}\forall \text{}\mathrm{i},\text{}\mathrm{k}=\mathrm{AG},\text{}\mathrm{t}$$

$${\mathrm{hydropower}}_{\mathrm{i}}^{}\left(\mathrm{t}\right)=\mathrm{R}{5}_{\mathrm{i}}^{}\times \left({{\displaystyle \sum}}_{\mathrm{j}\in \mathrm{down}\left(\mathrm{i}\right)}^{}{\mathrm{flow}}_{\mathrm{i},\mathrm{j}}^{}+{{\displaystyle \sum}}_{\mathrm{k}}^{}{\mathrm{supply}}_{\mathrm{i},\mathrm{k}}^{}+{\mathrm{spill}}_{\mathrm{i}}^{}\right)\text{}\forall \text{}\mathrm{i},\text{}\mathrm{t}$$

$${\mathrm{supply}}_{\mathrm{i},\mathrm{k}}^{}\left(\mathrm{t}\right)\ge {\mathrm{waterdmd}}_{\mathrm{i},\mathrm{k}}^{}\left(\mathrm{t}\right)\text{}\forall \text{}\mathrm{i},\text{}\mathrm{k},\text{}\mathrm{t}$$

$${\mathrm{s}}_{\mathrm{i}}^{}\left(\mathrm{t}\right)\le {\mathrm{STORAGE}}_{\mathrm{i}}^{}\text{}\forall \text{}\mathrm{i},\text{}\mathrm{t}$$

$$\mathrm{biofuel},{\text{}\mathrm{flow}}_{\mathrm{i},\mathrm{j}}^{},{\text{}\mathrm{food}}_{\mathrm{p}}^{},{\text{}\mathrm{hydropower}}_{\mathrm{i}}^{}\ge 0\text{}\forall \text{}\mathrm{i},\text{}\mathrm{j},\text{}\mathrm{p},\text{}\mathrm{t}$$

$${\mathrm{s}}_{\mathrm{i}}^{},{\text{}\mathrm{spill}}_{\mathrm{i}}^{},{\text{}\mathrm{supply}}_{\mathrm{i},\mathrm{k}}^{},{\text{}\mathrm{waterdmd}}_{\mathrm{i},\mathrm{k}}^{},{\text{}\mathrm{x}}_{\mathrm{i},\mathrm{p}}^{}\ge 0\text{}\forall \text{}\mathrm{i},\text{}\mathrm{j},\text{}\mathrm{k},\text{}\mathrm{p},\text{}\mathrm{t}$$

The food-water-energy system dynamics model is a time-dependent simulation problem. The study applies the system dynamics simulation method to analyze the model in Equations (1)–(11). The method simulates control and state variables of the system under different uncertain scenarios.

This research analyzes the system dynamics model of the urban metabolic food-water-energy system in Equations (1)–(11). A system dynamics model is established by using NetLogo to simulate the components interaction and system performance. Different scenarios are simulated and compared by the NetLogo system dynamics model. NetLogo is a simulation modeling platform and is able to conduct system dynamics and agent-based simulations. NetLogo is an open source software widely applied for simulations in physics, chemistry, biology, economics, social science and other related research fields. The model simulates agent-based interaction by using stationary Patch agents and mobile Turtle agents. Additionally, Link agents are used for connecting multiple agents, and Observer agents are set to record individual behavior and system performance.

The research conducts a case study of the urban metabolic food-water-energy system in the Shihmen Reservoir system. The food-water-energy nexus system includes the Shihmen Reservoir water resources system, hydropower generation system, paddy rice land, and household water consumption in Taoyuan, New Taipei, and Hsinchu cities (see Figure 2). Accordingly, water supply, water demand, hydropower generation, and paddy rice production are considered and simulated in the system dynamics modeling.

This research constructs a system dynamics model. The model includes system dynamics simulation and agent-based models. Figure 3 presents the system dynamics simulation of the food-water-energy system for urban metabolic analysis. The system dynamics model simulates time-varying water demand, water supply, food production, hydropower generation, and related resources interactions. In this research, an agent-based model is established to examine the interactive behavior and system performance of food-water-energy. The model displays a multi-agent based NetLogo model with patch, turtle, plot, and monitor agents to analyze the agent interactions. In the multi-agent based NetLogo model, patches are created to simulate water, paddy rice, and hydropower agents. Turtle agents are established to analyze water demand and supply, food production and consumption, and hydropower generation. In addition, monitors and graphs of food, water, and energy are established to collect and display dynamics changes of the systems.

In the simulation, farm land for paddy rice production is about 36,000 hectares in Taoyuan; farm land includes approximately 24,000 hectares of the Taoyuan Irrigation Association and 12,000 hectares of the Shihmen Irrigation Association. Rice production rate is about 5 ton/hectare. For the water supply system, the Shihmen Reservoir is a multi-purpose reservoir located in the upstream of Dahan River in northern Taiwan (Figure 2). It is the third largest reservoir in Taiwan with a storage capacity of 210 million m^{3} and a catchment area of 760 km^{2}. The reservoir is confronted with a serious sediment accumulation problem. One-third of the active storage capacity has been taken over by sediments.

In the model, the Shihmen Reservoir supplies water to a population of 3.5 million people in Taoyuan, New Taipei, and Hsinchu cities. The Shihmen Reservoir delivers water to irrigation land for agricultural water consumption of the Taoyuan Irrigation Association and the Shihmen Irrigation Association. Agricultural water demand, on average, accounts for 45% of the total water supply of the Shihmen Reservoir while the rest goes to the domestic (household and industrial) water demand. Furthermore, the Shihmen Reservoir has two hydropower generators, and each of them has a capacity of 45 MW. The hydropower system has a maximal water flow of 140 m^{3}/s and total installed capacity of 90 MW. The average annual hydropower generation is about 230,000 MWh. The maximal rice production is estimated at 180,000 tons per season for 36,000 hectares of paddy rice land of the Taoyuan Irrigation Association and the Shihmen Irrigation Association. The average water consumption of paddy rice irrigation is approximately 350 million m^{3}.

In addition, uncertainty analysis of the water inflow of reservoir, water demand, and paddy rice production are simulated using a Monte Carlo simulation in the model. Table 1 demonstrates the simulation results under various scenarios. The percentage of uncertainty ranges from 5 to 40%. The total water inflow/outflow, average water storage and average generated hydropower under each scenario are shown with the standard deviation. In the results, the interactive influence between water storage and hydropower under uncertain scenarios are also presented. In the study, Scenarios 1 to 3 evaluate the 10, 20, and 30% uncertainty of water inflow of the Shihmen Reservoir and demand. Figure 4, Figure 5 and Figure 6 present the dynamic water storage of the Shihmen Reservoir correspondingly.

Food, water, and energy are essential and scarce resources. The interactive relationship and allocation of food, water, and energy resources are crucial topics. This research analyzes the integrated food-water-energy systems for urban metabolism. The contribution of this study is to establish a framework for quantifying the cooperation/competition of multiple resources and formulating the system dynamics models of the urban metabolic food-water-energy system. The system dynamics model is capable of simulating different uncertain scenarios for uncertain analysis. A case study of food-water-energy system is conducted consisting of the Shihmen Reservoir water resources system, hydropower generation, paddy rice production of the Taoyuan and Shihmen Irrigation Associations, and domestic water supply in Taoyuan, New Taipei, and Hsinchu cities. The results quantify the interactive influence between food, water, and energy resources. In addition, stochastic simulation of water inflow shows uncertain water supply, paddy rice growing, and hydropower production systems. Future studies include simulation of uncertain climate scenarios, human behavior, hydrological process, energy availability, and energy load of nexus systems. Interactive agent-based modeling with multi-objective in food-water-energy systems are also of interest.

Conceptualization, M.-C.H., C.F., T.H. and Y.-H.C.; Data curation, C.-F.W. and Y.-H.C.; Funding acquisition, M.-C.H. and C.F.; Investigation, M.-C.H., C.F. and C.-F.W.; Methodology, M.-C.H., C.F., T.H., C.-F.W. and Y.-H.C.; Supervision, T.H.; Writing—original draft, M.-C.H., C.-F.W. and Y.-H.C.; Writing—review & editing, T.H.

This research was funded by National Taiwan University under grant NTU-CCP-106R891007 and Ministry of Science and Technology of Taiwan under grant MOST-106-2627-M-002-027.

The authors would like to thank the anonymous referees and editors for their thoughtful comments and suggestions. The authors are responsible for the accuracy of the information presented in this paper and for all opinions expressed herein.

The authors declare no conflict of interest.

In this paper, lowercase letters are used to indicate indices and decision variables, and uppercase letters indicate given parameters. Indices, given parameters, and variables are listed below and followed by the definitions.

Indices | |

$\mathrm{i},\text{}\mathrm{j}$ | node of water network, $\mathrm{i},\text{}\mathrm{j}=1,\text{}2,\text{}3,\text{}\dots ,\text{}\mathrm{I}$ |

$\mathrm{k}$ | residential, industrial, and agricultural water demand, $\mathrm{i}=\mathrm{RES},\text{}\mathrm{IND},\text{}\mathrm{AG}$ |

$\mathrm{p}$ | crop, $\mathrm{p}=1,\text{}2,\text{}3,\text{}\dots ,\text{}\mathrm{P}$ |

$\mathrm{t}$ | time, $\mathrm{t}=1,\text{}2,\text{}3,\text{}\dots ,\text{}\mathrm{T}$ |

Parameters | |

${\mathrm{FARM}}_{\mathrm{i}}^{}$ | farm size with water supply from node $\mathrm{i}$ |

${\mathrm{IN}}_{\mathrm{i}}^{}$ | water inflow at node $\mathrm{i}$ |

${\mathrm{LOSS}}_{\mathrm{i}}^{}$ | water loss at node $\mathrm{i}$ |

$\mathrm{R}{1}_{\mathrm{p}}^{}$ | production rate of crop $\mathrm{p}$ |

$\mathrm{R}{2}_{\mathrm{p}}^{}$ | biofuel production rate of crop $\mathrm{p}$ |

$\mathrm{R}{3}_{\mathrm{j},\mathrm{i}}^{}$ | water loss rate from node $\mathrm{j}$ to node $\mathrm{i}$ |

$\mathrm{R}{4}_{\mathrm{p}}^{}$ | water demand rate of crop $\mathrm{p}$ |

$\mathrm{R}{5}_{\mathrm{i}}^{}$ | hydropower production rate of node $\mathrm{i}$ |

${\mathrm{STORAGE}}_{\mathrm{i}}^{}$ | water storage capacity at node $\mathrm{i}$ |

Variables | |

$\mathrm{biofuel}$ | biofuel production |

${\mathrm{flow}}_{\mathrm{i},\mathrm{j}}^{}$ | water flow from node $\mathrm{i}$ to node $\mathrm{j}$ |

${\mathrm{food}}_{\mathrm{p}}^{}$ | food production of crop $\mathrm{p}$ |

${\mathrm{hydropower}}_{\mathrm{i}}^{}$ | hydropower generation at node $\mathrm{i}$ |

${\mathrm{s}}_{\mathrm{i}}^{}$ | water storage at node $\mathrm{i}$ |

${\mathrm{spill}}_{\mathrm{i}}^{}$ | water spillage at node $\mathrm{i}$ |

${\mathrm{supply}}_{\mathrm{i},\mathrm{k}}^{}$ | water supply of $\mathrm{i}$ to demand $\mathrm{k}$ |

${\mathrm{waterdmd}}_{\mathrm{i},\mathrm{k}}^{}$ | water demand $\mathrm{k}$ supplied from node $\mathrm{i}$ |

${\mathrm{x}}_{\mathrm{i},\mathrm{p}}^{}$ | tillage size of crop $\mathrm{p}$ with water supply from node $\mathrm{i}$ |

- Lei, K.; Liu, L.; Lou, I. An evaluation of the urban metabolism of Macao from 2003 to 2013. Resour. Conserv. Recycl.
**2018**, 128, 479–488. [Google Scholar] [CrossRef] - Ramaswami, A.; Boyer, D.; Nagpure, A.S.; Fang, A.; Bogra, S.; Bakshi, B.; Cohen, E.; Rao-Ghorpade, A. An urban systems framework to assess the trans-boundary food-energy-water nexus: Implementation in Delhi, India. Environ. Res. Lett.
**2017**, 12, 025008. [Google Scholar] [CrossRef] - Newman, P.W. Sustainability and cities: Extending the metabolism model. Landsc. Urban Plan.
**1999**, 44, 219–226. [Google Scholar] [CrossRef] - Kennedy, C.; Pincetl, S.; Bunje, P. The study of urban metabolism and its applications to urban planning and design. Environ. Pollut.
**2011**, 159, 1965–1973. [Google Scholar] [CrossRef] [PubMed] - Pincetl, S.; Bunje, P.; Holmes, T. An expanded urban metabolism method: Toward a systems approach for assessing urban energy processes and causes. Landsc. Urban Plan.
**2012**, 107, 193–202. [Google Scholar] [CrossRef] - Feng, Y.Y.; Chen, S.Q.; Zhang, L.X. System dynamics modeling for urban energy consumption and CO
_{2}emissions: A case study of Beijing, China. Ecol. Model.**2013**, 252, 44–52. [Google Scholar] [CrossRef] - Zhang, Y. Urban metabolism: A review of research methodologies. Environ. Pollut.
**2013**, 178, 463–473. [Google Scholar] [CrossRef] [PubMed] - Rosado, L.; Niza, S.; Ferrão, P. A material flow accounting case study of the Lisbon metropolitan area using the urban metabolism analyst model. J. Ind. Ecol.
**2014**, 18, 84–101. [Google Scholar] [CrossRef] - Bekchanov, M.; Lamers, J.P.A. The Effect of Energy Constraints on Water Allocation Decisions: The Elaboration and Application of a System-Wide Economic-Water-Energy Model (SEWEM). Water
**2016**, 8, 253. [Google Scholar] [CrossRef] - Chang, Y.; Li, G.; Yao, Y.; Zhang, L.; Yu, C. Quantifying the Water-Energy-Food Nexus: Current Status and Trends. Energies
**2016**, 9, 65. [Google Scholar] [CrossRef] - Russo, T.; Alfredo, K.; Fisher, J. Sustainable Water Management in Urban, Agricultural, and Natural Systems. Water
**2014**, 6, 3934–3956. [Google Scholar] [CrossRef][Green Version] - Blekking, J.; Tuholske, C.; Evans, T. Adaptive Governance and Market Heterogeneity: An Institutional Analysis of an Urban Food System in Sub-Saharan Africa. Sustainability
**2017**, 9, 2191. [Google Scholar] [CrossRef] - Liu, L.; Feng, C.; Zhang, H.; Zhang, X. Game Analysis and Simulation of the River Basin Sustainable Development Strategy Integrating Water Emission Trading. Sustainability
**2015**, 7, 4952–4972. [Google Scholar] [CrossRef][Green Version] - Tracy, C.R. A model of the dynamic exchanges of water and energy between a terrestrial amphibian and its environment. Ecol. Monogr.
**1976**, 46, 293–326. [Google Scholar] [CrossRef] - Lee, C.L.; Huang, S.L.; Chan, S.L. Synthesis and spatial dynamics of socio-economic metabolism and land use change of Taipei Metropolitan Region. Ecol. Model.
**2009**, 220, 2940–2959. [Google Scholar] [CrossRef] - Bazilian, M.; Rogner, H.; Howells, M.; Hermann, S.; Arent, D.; Gielen, D.; Steduto, P.; Mueller, A.; Kmoor, P.; Tol, R.S.J.; et al. Considering the energy, water and food nexus: Towards an integrated modelling approach. Energy Policy
**2011**, 39, 7896–7906. [Google Scholar] [CrossRef] - Calise, F.; Cipollina, A.; d’Accadia, M.D.; Piacentino, A. A novel renewable polygeneration system for a small Mediterranean volcanic island for the combined production of energy and water: Dynamic simulation and economic assessment. Appl. Energy
**2014**, 135, 675–693. [Google Scholar] [CrossRef] - Chen, S.; Chen, B. Coupling of carbon and energy flows in cities: A meta-analysis and nexus Modelling. Appl. Energy
**2017**, 194, 774–783. [Google Scholar] [CrossRef] - Chen, P.-C.; Alvarado, V.; Hsu, S.-C. Water energy nexus in city and hinterlands: Multi-regional physical input-output analysis for Hong Kong and South China. Appl. Energy
**2008**, 225, 986–997. [Google Scholar] [CrossRef] - De Stercke, S.; Mijic, A.; Buytaert, W.; Chaturvedi, V. Modelling the dynamic interactions between London’s water and energy systems from an end-use perspective. Appl. Energy
**2018**, 230, 615–626. [Google Scholar] [CrossRef] - Facchini, A.; Kennedy, C.; Stewart, I.; Mele, R. The energy metabolism of megacities. Appl. Energy
**2017**, 186, 86–95. [Google Scholar] [CrossRef] - Martinez-Hernandez, E.; Leach, M.; Yang, A. Understanding water-energy-food and ecosystem interactions using the nexus simulation tool NexSym. Appl. Energy
**2017**, 206, 1009–1021. [Google Scholar] [CrossRef] - Ravalde, T.; Keirstead, J. Comparing performance metrics for multi-resource systems: The case of urban metabolism. J. Clean. Prod.
**2017**, 163, S241–S253. [Google Scholar] [CrossRef] - Schlör, H.; Venghaus, S.; Zapp, P.; Marx, J.; Schreiber, A.; Hake, J.-F. The energy-mineral-society nexus—A social LCA model. Appl. Energy
**2018**, 228, 999–1008. [Google Scholar] [CrossRef] - Wang, S.; Chen, B. Energy–water nexus of urban agglomeration based on multiregional input–output tables and ecological network analysis: A case study of the Beijing–Tianjin–Hebei region. Appl. Energy
**2016**, 178, 773–783. [Google Scholar] [CrossRef] - Yang, X.; Wang, Y.; Sun, M.; Wang, R.; Zheng, P. Exploring the environmental pressures in urban sectors: An energy-watercarbon nexus perspective. Appl. Energy
**2018**, 228, 2298–2307. [Google Scholar] [CrossRef] - Zhang, C.; Chen, X.; Li, Y.; Ding, W.; Fu, G. Water-energy-food nexus: Concepts, questions and methodologies. J. Clean. Prod.
**2018**, 195, 625–639. [Google Scholar] [CrossRef] - Carbonell, G.; Ramos, C.; Pablos, M.V.; Ortiz, J.A.; Tarazona, J.V. A system dynamic model for the assessment of different exposure routes in aquatic ecosystems. Sci. Total Environ.
**2000**, 247, 107–118. [Google Scholar] [CrossRef] - Badescu, V. Dynamic model of a complex system including PV cells, electric battery, electrical motor and water pump. Energy
**2003**, 28, 1165–1181. [Google Scholar] [CrossRef] - Kato, T. Simulation of water quality with the application of system dynamics model for population and land-use changes. Paddy Water Environ.
**2005**, 3, 103–109. [Google Scholar] [CrossRef] - Churkina, G. Modeling the carbon cycle of urban systems. Ecol. Model.
**2008**, 216, 107–113. [Google Scholar] [CrossRef] - Jin, W.; Xu, L.; Yang, Z. Modeling a policy making framework for urban sustainability: Incorporating system dynamics into the Ecological Footprint. Ecol. Econ.
**2009**, 68, 2938–2949. [Google Scholar] [CrossRef] - Liu, T.M.; Tung, C.P.; Ke, K.Y.; Chuang, L.H.; Lin, C.Y. Application and development of a decision-support system for assessing water shortage and allocation with climate change. Paddy Water Environ.
**2009**, 7, 301–311. [Google Scholar] [CrossRef] - Luo, G.; Yin, C.; Chen, X.; Xu, W.; Lu, L. Combining system dynamic model and CLUE-S model to improve land use scenario analyses at regional scale: A case study of Sangong watershed in Xinjiang, China. Ecol. Complex.
**2010**, 7, 198–207. [Google Scholar] [CrossRef] - Chen, S.; Chen, B. Network environ perspective for urban metabolism and carbon emissions: A case study of Vienna, Austria. Environ. Sci. Technol.
**2012**, 46, 4498–4506. [Google Scholar] [CrossRef] - Li, Y.H.; Chen, P.Y.; Lo, W.H.; Tung, C.P. Integrated water resources system dynamics modeling and indicators for sustainable rural community. Paddy Water Environ.
**2015**, 13, 29–41. [Google Scholar] [CrossRef] - Chalabi, Z.S.; Bailey, B.J.; Wilkinson, D.J. A real-time optimal control algorithm for greenhouse heating. Comput. Electron. Agric.
**1996**, 15, 1–13. [Google Scholar] [CrossRef] - Brouwer, D.R.; Jansen, J.D. Dynamic optimization of water flooding with smart wells using optimal control theory. In Proceeding of the European Petroleum Conference, Aberdeen, UK, 29–31 October 2002. [Google Scholar]
- Sohngen, B.; Mendelsohn, R. An optimal control model of forest carbon sequestration. Am. J. Agric. Econ.
**2003**, 85, 448–457. [Google Scholar] [CrossRef] - López-Ibáñez, M.; Prasad, T.D.; Paechter, B. Ant colony optimization for optimal control of pumps in water distribution networks. J. Water Res. Plan. Manag.
**2008**, 134, 337–346. [Google Scholar] [CrossRef] - Zhang, Q.; Maeda, S.; Kawachi, T. Stochastic multiobjective optimization model for allocating irrigation water to paddy fields. Paddy Water Environ.
**2007**, 5, 93–99. [Google Scholar] [CrossRef] - Chono, S.; Maeda, S.; Kawachi, T.; Imagawa, C.; Buma, N.; Takeuchi, J. Optimization model for cropping-plan placement in paddy fields considering agricultural profit and nitrogen load management in Japan. Paddy Water Environ.
**2012**, 10, 113–120. [Google Scholar] [CrossRef]

Uncertainty | Water Inflow | Water Outflow | Water Storage | Hydropower |
---|---|---|---|---|

Total (Standard Deviation) | Total (Standard Deviation) | Mean (Standard Deviation) | Mean (Standard Deviation) | |

million m^{3} | million m^{3} | million m^{3} | MWh | |

5% | 776.13 | 773.61 | 124.82 | 322,338 |

(4.92) | (15.60) | (15.57) | (6499) | |

10% | 785.45 | 789.23 | 133.36 | 328,846 |

(5.35) | (17.72) | (19.81) | (7384) | |

15% | 817.02 | 787.63 | 150.95 | 328,181 |

(7.10) | (19.17) | (15.30) | (7988) | |

20% | 764.54 | 754.18 | 135.54 | 314,241 |

(5.82) | (14.16) | (13.06) | (5900) | |

25% | 759.71 | 820.12 | 59.62 | 341,717 |

(12.27) | (22.75) | (36.85) | (9479) | |

30% | 763.26 | 776.79 | 108.59 | 323,662 |

(12.55) | (16.74) | (21.53) | (6975) | |

35% | 759.14 | 759.88 | 143.83 | 316,616 |

(16.49) | (18.18) | (19.40) | (7576) | |

40% | 834.72 | 754.82 | 121.33 | 314,510 |

(15.68) | (24.33) | (39.35) | (10,136) |

© 2018 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).