Dynamic Coupling Analysis of Urbanization and Water Resource Utilization Systems in China

While urbanization brings economic and social benefits, it also causes water pollution and other environmental ecological problems. This paper provides a theoretical framework to quantitatively analyze the dynamic relationship between water resource utilization and the process of urbanization. Using data from Jiangsu province, we first construct indices to evaluate urbanization and water resource utilization. We then adopt an entropy model to examine the correlation between urbanization and water resource utilization. In addition, we introduce a dynamic coupling model to analyze and predict the coupling degree between urbanization and water resource utilization. Our analyses show that pairing with rising urbanization during 2002–2014, the overall index of water resource utilization in Jiangsu province has experienced a “decline -rise-decline” trend. Specifically, after the index of water resource utilization reached its lowest point in 2004, it gradually began to rise. Water resource utilization reached its highest value in 2010. The coupling degree between urbanization and water resource utilization was relatively low in 2002 and 2003 varying between −90◦ and 0◦. It has been rising since then. Out-of-sample forecasts indicate that the coupling degree will reach its highest value of 74.799◦ in 2016, then will start to gradually decline. Jiangsu province was chosen as our studied area because it is one of the selected pilot provinces for China’s economic reform and social development. The analysis of the relationship between provincial water resource utilization and urbanization is essential to the understanding of the dynamic relationship between these two systems. It also serves as an important input for developing national policies for sustainable urbanization and water resource management.


Introduction
Urbanization draws the rural population into cities as secondary and tertiary industries, and also develop into urban areas.As a result, the scale of urbanization expands over time [1,2].Urbanization also reflects changes in industrialization, technological modernization, and social changes in the country [3].In the 1800s, the world's urban population accounted for 2.4% of the total population.By 1950, it had increased to 29.2%.The urban population then surged to 42.6% by 1990, and rising further to 50% by 2000.The global urban population is expected to reach 62.5% by 2025 [4].Greater increasing urbanization increases the urban population, the number of cities, and also improves social wellbeing and promotes intensified uses of economic resources.Urban development is a mixture of a complex spatial morphological changes and social progress [5,6].While urbanization brings many social and economic benefits, it can also lead to serious ecological and environmental problems caused by industrialization and dense population concentrations.With rapid urbanization, water scarcity, water pollution, and issues of water usage efficiency have become prominent problems in some Chinese areas.Potential struggles between urbanization and water conservation have clearly become an important issue [7][8][9].Therefore, it is important to learn how to coordinate water resource utilization during urbanization process.
Urbanization and efficient water resource utilization are essential to economic development.Studies on this topic focus mainly on two aspects: first, how do water resources support and accelerate the urbanization process.Abundant water resources not only provide for the daily needs of city-dwellers, water is also an important input for various industrial and tertiary activities [10,11].Furthermore, abundant water resources support improved ecological, environmental, and climatological factors that are attractive to business firms and talented workers.This can increase both the speed and quality of urbanization [12].Lei and Chen [13] studied the importance of East Lake to the urbanization of Wuhan city using the water ecological carrying capacity (WECC) measure.Studying Barcelona and Madrid, March [14] considers how urbanization mobilizes water resources to keep pace with sustainable urban growth.
The second topic studied relates to changes in water resources accompanying urbanization.Urbanization is usually characterized by construction that increases impervious areas and reductions in agricultural and forest lands.This can lead to many well-known environmental problems, such as water shortages [15,16], water and air pollution [17], and the creation of urban heat islands [18].Studies have demonstrated that degradation of river structures has created severe problems related to river health, including the "urban stream syndrome".These problems can be particularly devastating to human development [19][20][21].On the other hand, some remarkable improvements in water quality accompanied by urbanization in highly developed watersheds.primarily due to the influence of comprehensive water management practices (particularly for COD, BOD, NH3-N, and TP) [22].Shrestha et al. [23] found that the Bagmati river water quality was so chemically and biologically degraded that it could not be used for any purpose due to rapid population growth and urbanization.This is particularly a problem during the region's dry season.
While many studies focus on the linkage between urbanization and water resource utilization, few have conducted systematic analyses of the interaction between the two.In fact, urbanization itself is an integral process of urban growth and water resource utilization.All elements in these two processes function together under both social and economic restraints.
The coupling effect takes a spiraling progression from a lower degree to a higher degree of coordination [24][25][26].In this study, we analyze the coupling effect between the evolution of urbanization and water resources utilization using data from Jiangsu province in Southeastern China.This province's economy is booming and is highly urbanized.Its GDP in 2014 was about $10 trillion, and it was 65.2% urbanized.In addition, Jiangsu's water system is unique: it has insufficient local water resources but abundant transit water [27].Using the coupling model, we analyze trends associated with coordinating development between urbanization and water utilization.We further forecast the future coupling degrees between the two systems and provide some policy recommendations for sustainable urbanization in Jiangsu province.By presenting a theoretical framework and quantitative analyses of the relationship between water resources and urbanization processes in this province, we hope to provide provincial and local governments with a better understanding of the interplay of resource management and urbanization.This will assist provincial governments in implementing effective policies for sustainable water resource management as urbanization progresses.

Description of Studied Area
Jiangsu province is located on the southeastern coast of China (see Figure 1).It consists of 13 cities, including Nanjing and Suzhou, covering an area of 102,600 km 2 , accounting for 1.06 percent of China's land mass.Jiangsu is 70,600 km 2 of plain with a water surface area of 17,300 km 2 .This province's climate is in a transition zone between subtropical and temperate.Separated by the Huaihe River and an irrigation canal, the southern part of Jiangsu is in the humid, subtropical, monsoon-climate area, while the northern part has a drier, semi-humid climate [28].Annual precipitation gradually increases from north to south, with average annual rainfall between 700-1100 mm.Annual average sunshine hours are between 2000-2400 h.The average temperature also gradually increases from north to south, but the difference is negligible [29].The annual average temperature is between 14-16 • C and annual average wind speed is between 2-4 m per second with relatively small regional differences.
The economy of Jiangsu province has been growing rapidly in the 21st century.GDP per capita was $1423 in 2000 growing to $13,328 by 2014.Rapid economic growth accelerated urbanization with the urban population increasing from 41.5% in 2000 to 65.2% in 2014.Employment in urban areas expanded from 6.7 million to 15.1 million between 2000 and 2014.
Water resources in Jiangsu province are unevenly distributed, both seasonally and locationally.Precipitation in the southern part is greater than the northern region with 60%-70% of the rain falling between June and September.Average annual surface water evaporation is about 950-1100 mm.Water evaporation, however, increases from southwest to the north.Most precipitation evaporates with the annual average runoff depth of about 246 mm, which is only about 2.25 percent of total rainfall.The province is located downstream of the Yangtze and Huaihe Rivers.About 1.03 trillion m 3 of water flows into the province annually from these rivers.This serves as the main water source for urban residents and industries.Since the province is flat and not suitable for constructing dams or large reservoirs, it is difficult to retain transitional and local surface water.Poor capacity in regulating or storing water has led to insufficient supplies of water [30,31] with seasonal and regional water shortages often occurring in dry years.
In 2014, the total provincial water usage was 48 billion m 3 .The industrial and agricultural water usage was 44.2 billion m 3 (92% of total water usage), domestic water usage was 3.58 billion m 3 (7.4% of total water usage), and the remaining 0.6% was urban ecological and environmental usage.Compared to water usage in 2013, industrial and agricultural water usage decreased by 3.9%, while domestic usage increased by 0.8%.One serious problem with provincial water resources is water pollution.In 2014, industrial waste water totaled 2.05 billion tons, and domestic waste water reached 3.96 billion tons.To address these issues, understanding of the dynamic relationship between water resource utilization and urbanization is essential as urbanization progresses.
province's climate is in a transition zone between subtropical and temperate.Separated by the Huaihe River and an irrigation canal, the southern part of Jiangsu is in the humid, subtropical, monsoon-climate area, while the northern part has a drier, semi-humid climate [28].Annual precipitation gradually increases from north to south, with average annual rainfall between 700-1100 mm.Annual average sunshine hours are between 2000-2400 h.The average temperature also gradually increases from north to south, but the difference is negligible [29].The annual average temperature is between 14-16 °C and annual average wind speed is between 2-4 m per second with relatively small regional differences.
The economy of Jiangsu province has been growing rapidly in the 21st century.GDP per capita was $1423 in 2000 growing to $13,328 by 2014.Rapid economic growth accelerated urbanization with the urban population increasing from 41.5% in 2000 to 65.2% in 2014.Employment in urban areas expanded from 6.7 million to 15.1 million between 2000 and 2014.
Water resources in Jiangsu province are unevenly distributed, both seasonally and locationally.Precipitation in the southern part is greater than the northern region with 60%-70% of the rain falling between June and September.Average annual surface water evaporation is about 950-1100 mm.Water evaporation, however, increases from southwest to the north.Most precipitation evaporates with the annual average runoff depth of about 246 mm, which is only about 2.25 percent of total rainfall.The province is located downstream of the Yangtze and Huaihe Rivers.About 1.03 trillion m 3 of water flows into the province annually from these rivers.This serves as the main water source for urban residents and industries.Since the province is flat and not suitable for constructing dams or large reservoirs, it is difficult to retain transitional and local surface water.Poor capacity in regulating or storing water has led to insufficient supplies of water [30,31] with seasonal and regional water shortages often occurring in dry years.
In 2014, the total provincial water usage was 48 billion m 3 .The industrial and agricultural water usage was 44.2 billion m 3 (92% of total water usage), domestic water usage was 3.58 billion m 3 (7.4% of total water usage), and the remaining 0.6% was urban ecological and environmental usage.Compared to water usage in 2013, industrial and agricultural water usage decreased by 3.9%, while domestic usage increased by 0.8%.One serious problem with provincial water resources is water pollution.In 2014, industrial waste water totaled 2.05 billion tons, and domestic waste water reached 3.96 billion tons.To address these issues, understanding of the dynamic relationship between water resource utilization and urbanization is essential as urbanization progresses.

The Coupling Model
There is a dynamic coupling relationship between urbanization and water resource utilization.This relationship can c considered as a non-equilibrium, self-regulating, and dynamically-oscillating system [32].This study applies the idea of systematic evolution to establish a dynamic coupling model between urbanization and water resource utilization.The progression of urbanization and water resource utilization is a nonlinear process and the evolution equation is: Since the nonlinear system's stability is mainly determined by the characteristics of the root function, we can expand Equation (1) at t = 0 as a Taylor series expansion given the system's stability.By omitting the higher-order terms E(x 1 , x 2 , x n ), a linear function of the system can be approximated as: The evolution function for urbanization and water resource utilization can be expressed as: where x, y represent the elements of urbanization and water utilization respectively; and α, β are their respective weights.Given the complex relationships between urbanization and water resource utilization, we transform the two systems into one comprehensive system.According to system structure theory, the evolution functions can be written as: where A and B represent the dynamic evolution trend of urbanization and water resource utilization respectively, and V A and V B are the evolution speeds under various conditions for the two systems.The function of the evolution speed for the comprehensive/overall system, V, can be expressed as V = F(V A , V B ).The path of function V can be used to analyze the dynamic coupling effect between urbanization and water utilization systems.
Considering the evolution path of function V as the result of the interactions between V A and V B , the coordinated system can be used to analyze the evolution path of V A and V B as shown in Figure 2. The angle between V A and V B , α, represents the coupling degree of the two systems, and α satisfies tanα = V A /V B .By calculating α, α = arctan(V A /V B ), the evolution trend of the overall system can be further depicted.The dynamic coupling degree between urbanization and water resources utilization is clearly presented in Figure 2.
Based on the values of the coupling degree between urbanization and water utilization, there are four stages in the coupling process shown in Panel A in Figure 2: quadrant I, coordinated symbiosis stage; quadrant II, coordinated development stage; quadrant III, limited development stage; and quadrant IV, spiral progressing stage.
In quadrant I, −90 • < α < 0 • , the two systems are in stage I: coordinated symbiosis.During this stage, the two systems develop independently.The urbanization process is relatively slow and water resource constraints are relatively weak.There is relatively low integration between the two systems.In quadrant II, 0° < α < 90°, the two systems are in a coordinated development stage.With the acceleration of urbanization, problems with water shortages and pollution start to develop.Conflicts emerge between the two systems as urbanization progresses.
In quadrant III, 90° < α < 180°, the two systems are in a limited development stage.The rapid growth of urbanization further intensifies the conflicts between urbanization and water resources.With the accelerating urbanization, the demand for water resources surges sharply due to the growing demographic and economic factors.This causes serious ecological and environmental pollution and promotes misuse of water resources.This further strains the relationship between urbanization and water utilization.At this stage, the overall system may evolve into two possible scenarios: (1) conflicts between urbanization and water utilization cannot be resolved.As water shortage problems increase, the whole system collapses and civilization and its economy regress; and (2) people attempt to confront the water crisis caused by urbanization.The community can adopt effective measures to increase the local water supplies with government stringently regulating industrial and water-consuming activities to optimize water demand [33].With continued adjustments and regulations, the internal elements of the comprehensive plan can be optimally harmonized.Then the coupling relationship between urbanization and water utilization can evolve favorably until the coupling relationship reaches the optimal symbiosis as shown in Panel B of Figure 2.
In quadrant IV, −180° < α < −90°, the coupling degree is in the rising spiral stage.As urbanization progresses, industrial and residential water usage improve and are optimized.A water-conserving society is established and enhanced.The coupling relationship between urbanization and water utilization is stabilized, and the overall system is in the rising spiral stage.
Although the four quadrants represent the entire interaction process of the two systems, it does not mean every provincial/local interaction between water resources and urbanization would experience all four quadrants.Areas/regions experiencing slow urban growth tend to stay in quadrants I and II for long periods of time.This suggests that slow urban growth does not present a threat to water resources.For rapid urbanization areas, the conflicts between urbanization and water resources may not be resolved, with the two systems in a "downward spiral" remaining in quadrant III.If the conflicts stay unresolved, it may ultimately lead to total destruction of the eco-system with the city becoming uninhabitable.Another possible development in quadrant III is that local government and residents adopt effective measures to resolve the conflicts.Then the two systems progress into quadrant IV-the coupling relationship between urbanization and water utilization is stabilized and optimized.

Forecasting Model
An exponential smoothing model is often used in forecasting time series models.It typically includes three smoothing methods: single exponential smoothing, double exponential smoothing, In quadrant II, 0 • < α < 90 • , the two systems are in a coordinated development stage.With the acceleration of urbanization, problems with water shortages and pollution start to develop.Conflicts emerge between the two systems as urbanization progresses.
In quadrant III, 90 • < α < 180 • , the two systems are in a limited development stage.The rapid growth of urbanization further intensifies the conflicts between urbanization and water resources.With the accelerating urbanization, the demand for water resources surges sharply due to the growing demographic and economic factors.This causes serious ecological and environmental pollution and promotes misuse of water resources.This further strains the relationship between urbanization and water utilization.At this stage, the overall system may evolve into two possible scenarios: (1) conflicts between urbanization and water utilization cannot be resolved.As water shortage problems increase, the whole system collapses and civilization and its economy regress; and (2) people attempt to confront the water crisis caused by urbanization.The community can adopt effective measures to increase the local water supplies with government stringently regulating industrial and water-consuming activities to optimize water demand [33].With continued adjustments and regulations, the internal elements of the comprehensive plan can be optimally harmonized.Then the coupling relationship between urbanization and water utilization can evolve favorably until the coupling relationship reaches the optimal symbiosis as shown in Panel B of Figure 2.
In quadrant IV, −180 • < α < −90 • , the coupling degree is in the rising spiral stage.As urbanization progresses, industrial and residential water usage improve and are optimized.A water-conserving society is established and enhanced.The coupling relationship between urbanization and water utilization is stabilized, and the overall system is in the rising spiral stage.
Although the four quadrants represent the entire interaction process of the two systems, it does not mean every provincial/local interaction between water resources and urbanization would experience all four quadrants.Areas/regions experiencing slow urban growth tend to stay in quadrants I and II for long periods of time.This suggests that slow urban growth does not present a threat to water resources.For rapid urbanization areas, the conflicts between urbanization and water resources may not be resolved, with the two systems in a "downward spiral" remaining in quadrant III.If the conflicts stay unresolved, it may ultimately lead to total destruction of the eco-system with the city becoming uninhabitable.Another possible development in quadrant III is that local government and residents adopt effective measures to resolve the conflicts.Then the two systems progress into quadrant IV-the coupling relationship between urbanization and water utilization is stabilized and optimized.

Forecasting Model
An exponential smoothing model is often used in forecasting time series models.It typically includes three smoothing methods: single exponential smoothing, double exponential smoothing, and triple exponential smoothing.Studies show that using triple exponential smoothing ensures superior forecasts for trended time series data [34,35], that is, the forecasts are closer to the actuals using triple exponential smoothing.Thus, we choose the triple exponential smoothing method to forecast the degrees of coupling.In triple exponential smoothing, we need to transform the original data.We first calculate the single exponential smoothing value, the double exponential smoothing value, and then the triple exponential smoothing value.Forecasts are then estimated.Suppose time series data are: X 1 , X 2 , X 3 , . . ., X n , the calculation process is as follows: where α is the smoothing factor, and 0 • < α < 1 • .The smoothed statistic, S t , is a simple weighted average of the observations X t and the previous smoothed statistic S t−1 .S t (1) is the single exponential smoothing value, S t (2) is the double exponential smoothing value, and S t (3) is the triple exponential smoothing value.The forecasting model with triple exponential smoothing is as follows: where Y t+T is the forecasted value, t is the base year, T is the forecasted period, a t , b t , and c t are the smoothing coefficients, which are calculated as: 2.3.Index System for Urbanization and Water Resources Utilization

Construction of the Urbanization Index System
Urbanization is an evolutionary process that transforms economic production and lifestyle.Studies [36][37][38][39][40] have shown that in constructing an urbanization evaluation index, four aspects should be considered: demographic, economic, social, and environmental factors.As urbanization progresses, urban population and city size continue to expand and the economic engine shifts from primary to secondary and tertiary industries.The concentration of resources and the changes in industrial structure bring rapid economic growth [36,37].Therefore, we include the urban to total population ratio, employment ratio of tertiary industries, and urban employment ratio to account for the demographic aspect.
Urbanization improves the industrial structure and resource utilization efficiency; as a result, the local economy and household incomes grow rapidly.We adopt GDP per capita and disposable income per urban household as proxies for economic growth and the ratio of value-added by tertiary industries to reflect changes in economic structure.
The growth of the urban population and changes in economic structure promote improvement in social welfare.Urban residents desire generous transfer programs/safety net, better education services, and better medical services, hoping to improve quality of life [38,39].Thus, we adopt four measures to reflect quality of social welfare.Social benefit expenses mainly include pensions for disabled and killed-in-action soldiers, living allowance for veterans, natural disaster relief funds, and reconstruction subsidies for post-extraordinary natural disasters.Social benefit expenses reflect the extent of financial assistance provided to the society's most vulnerable groups by regional governments.Therefore, we include these to reflect the quality of social welfare provided by the government.Employment in the health sector (in 10,000 s) reflects the quality of medical services.The number of college students is a proxy for the quality of education services.An area with a large number of college students tends to have more universities and research institutions, which can easily lead to technological innovations and developments to provide higher quality of social welfare.We also adopt the Engel coefficient of urban households as another indicator for quality of social welfare.A small Engel coefficient indicates a relatively low level of expenditure on food and beverages by urban households.This means the expenditures on social and other activities would be relatively high, which can improve the quality of social welfare.
Urbanization brings economic and social benefits to urban dwellers, but it also carries costs, including environmental pollution.Urbanization promotes rapid growth in population and industrial output.As a result of growing population and industrial output, air pollution and industrial solid waste become serious environmental problems.In the calculation of our urbanization index, we include two environmental factors reflecting air pollution and industrial waste: sulfur dioxide emissions and comprehensive utilization rate of industrial solid waste [22,23].Other factors, such as the ratio of urban green area to housing area, and green area per capita, are typical indicators for a positive eco-environment.
We construct our urbanization evaluation index using data from Jiangsu province.The detailed aspects of the four factors of our urbanization index are shown in Table 1.The determinations of weights of the index system are calculated using the entropy method.This method assigns weights based on raw data and are not influenced by human factors [54].
The method first conducts an information entropy calculation using standardized index data.
The formula is: Y ij , and Y ij is the standardized index data.The weight of entropy of the jth indicator could be defined as: , in which 0 ≤ W j ≤ 1, and m ∑ j=1 W j = 1.Calculated weights of urbanization index are shown in Table 1.

Construction of the Water Resources Utilization Index System
Water is an essential natural resource.Water shortages combined with inefficient uses of water often lead to lower water tables, degraded forests, disappearing lakes and rivers, and water pollution problems.This environmental degradation adversely affects development [55,56].In Jiangsu province, water supply differs by location and season.Thus, the distribution, protection, and utilization of water resources can have a significant impact on regional socio-economic development [57].Water resource utilization is a complex economic and engineering problem.We constructed a water resource utilization evaluation index for Jiangsu province based on three dimensions: volume of water resources, water utilization efficiency, and water environmental management.We include four factors to measure the volume of water resources.Three water utilization efficiency measures are included: overall, industrial, and agricultural water utilization efficiency.Four water environmental management measures are selected: industrial waste water, residential wastewater, ratio of sewage control to total pollution control investment, and the industrial wastewater compliance rate.The weight of each index layer is calculated using the same entropy method described above.Table 2 shows the factors for each aspect and their respective weights.

Data Standardization
The data for evaluating the water resources utilization index was not available until 2002, therefore, we collected data for the period of 2002-2014.
Since the properties of the two indices are not the same, their dimensions and magnitudes also differ.We use the min-max standardizing method to normalize the data.The formulas are as follows: Y ij = max x j − x ij max x j − min x j (15) where X ij is the initial value of the jth factor at year i.Y ij represents the corresponding standardized data.For positive data, we use Equation ( 14) to normalize it; for negative data, Equation ( 15) is used.The values of normalized data are between 0 and 1.

The Index of the Water Resource Utilization System for Jiangsu Province
Table 2 provides the weights for calculating the water resources utilization index.According to Table 2, the two most important factors affecting the water resources utilization index in Jiangsu are water usage per 10,000 GDP (with a weight of 0.177) and the total volume of production/industrial water (with a weight of 0.124).This means that the water usage per 10,000 GDP and production water demand are the two most important factors in evaluating water resource utilization in Jiangsu province.The computed indices are shown in Figure 3.
According to Figure 3, the index of urbanization is rising throughout the sample period, while the index of water resources utilization in Jiangsu province experienced a pattern of decline-increase-decline.The index of water utilization reached its lowest point in 2004 (0.299).It then gradually improved and reached its highest value in 2010 (0.614) after which it started to decline.There are three possible explanations for the high value in 2010: first, the total volume of production/industrial water reached the highest value, 51.5 billion m 3 .Second, after the "2007 Oxygenation Crisis of Taihu Lake", the provincial government significantly increased investment in managing water pollution and improving water resources utilization.As a result, the ratio of sewage to total pollution control investment reached a high value of 0.4, it then gradually declined to 0.219 in 2014.Third, industrial wastewater discharge compliance rate was the highest for the sample period, 98.03%.This rate started to decline after 2010.where Xij is the initial value of the jth factor at year i.Yij represents the corresponding standardized data.For positive data, we use Equation ( 14) to normalize it; for negative data, Equation ( 15) is used.The values of normalized data are between 0 and 1.

The Index of the Water Resource Utilization System for Jiangsu Province
Table 2 provides the weights for calculating the water resources utilization index.According to Table 2, the two most important factors affecting the water resources utilization index in Jiangsu are water usage per 10,000 GDP (with a weight of 0.177) and the total volume of production/industrial water (with a weight of 0.124).This means that the water usage per 10,000 GDP and production water demand are the two most important factors in evaluating water resource utilization in Jiangsu province.The computed indices are shown in Figure 3.
According to Figure 3, the index of urbanization is rising throughout the sample period, while the index of water resources utilization in Jiangsu province experienced a pattern of decline-increase-decline.The index of water utilization reached its lowest point in 2004 (0.299).It then gradually improved and reached its highest value in 2010 (0.614) after which it started to decline.There are three possible explanations for the high value in 2010: first, the total volume of production/industrial water reached the highest value, 51.5 billion m 3 .Second, after the "2007 Oxygenation Crisis of Taihu Lake", the provincial government significantly increased investment in managing water pollution and improving water resources utilization.As a result, the ratio of sewage to total pollution control investment reached a high value of 0.4, it then gradually declined to 0.219 in 2014.Third, industrial wastewater discharge compliance rate was the highest for the sample period, 98.03%.This rate started to decline after 2010.The pattern of the water utilization index can be approximated with a nonlinear fitting process as shown in Equation (10).We find that the cubic function has the best fit as shown in Equation ( 16) and Figure 4.The summary statistic for the regression is as follow: R 2 = 0.891, F value = 24.439, the t statistics for the coefficients are 5.666, −2.181, −3.585, and 3.302, respectively; all are significant at the 5% level.The pattern of the water utilization index can be approximated with a nonlinear fitting process as shown in Equation (10).We find that the cubic function has the best fit as shown in Equation ( 16) and Figure 4.The summary statistic for the regression is as follow: R 2 = 0.891, F value = 24.439, the t statistics for the coefficients are 5.666, −2.181, −3.585, and 3.302, respectively; all are significant at the 5% level.

The Index of Urbanization in Jiangsu Province
Among the 14 factors affecting the urbanization index in Table 1, the two most factors are urban employment ratio and employment in the health sector, with weights of 0.185 and 0.093 respectively.This indicates that the urban employment share and urban health environment are most important in computing the urbanization index for Jiangsu province.
The computed index of urbanization is depicted in Figure 3.According to this index, urbanization in Jiangsu province has increased markedly, resulting from rapid economic growth and increased industrial production.In 2002, the index of urbanization was only 0.044, but had reached 0.978 by 2013.The urbanization process of Jiangsu province has been much faster than that of other provinces in China.One of the major drivers for this has been provincial government policies.On 1 May 2003, Jiangsu implemented a "unified household registration" policy.It abolished separate registrations for agricultural residents and non-agricultural residents.The provincial government also relaxed regulations for migrant workers.Anyone with a legal local residency or a steady job in the province was allowed to apply for urban household registration.These policies drew many talented workers to the province, thus speeding up the provincial urbanization after 2003.In addition, the provincial government introduced the "urban planning template in Jiangsu Province (2012-2030)" in 2012 to help develop an affluent society and to promote sustainable urbanization.As a result the urbanization index jumped from 0.705 in 2011 to 0.923 in 2013.
Studies have shown that many regional urbanization trends can be approximated with the S-curve [84,85].We adopt a logistic curve to fit the pattern of urbanization index.The results show the logistic curve is good fit (R 2 = 0.93, F value = 146.444)as shown in Equation ( 17) and Figure 5.The t values for the two coefficients are 5.490 and 43.579 respectively, and both are significant at 1% level.Equation ( 17) can be transformed into Equation (18), which is used to compute the forecasted index of urbanization: ln 2 9 .6 2 ln 0 .6 2 7 1 ( )

The Index of Urbanization in Jiangsu Province
Among the 14 factors affecting the urbanization index in Table 1, the two most are urban employment ratio and employment in the health sector, with weights of 0.185 and 0.093 respectively.This indicates that the urban employment share and urban health environment are most important in computing the urbanization index for Jiangsu province.
The computed index of urbanization is depicted in Figure 3.According to this index, urbanization in Jiangsu province has increased markedly, resulting from rapid economic growth and increased industrial production.In 2002, the index of urbanization was only 0.044, but had reached 0.978 by 2013.The urbanization process of Jiangsu province has been much faster than that of other provinces in China.One of the major drivers for this has been provincial government policies.On 1 May 2003, Jiangsu implemented a "unified household registration" policy.It abolished separate registrations for agricultural residents and non-agricultural residents.The provincial government also relaxed regulations for migrant workers.Anyone with a legal local residency or a steady job in the province was allowed to apply for urban household registration.These policies drew many talented workers to the province, thus speeding up the provincial urbanization after 2003.In addition, the provincial government introduced the "urban planning template in Jiangsu Province (2012-2030)" in 2012 to help develop an affluent society and to promote sustainable urbanization.As a result the urbanization index jumped from 0.705 in 2011 to 0.923 in 2013.
Studies have shown that many regional urbanization trends can be approximated with the S-curve [84,85].We adopt a logistic curve to fit the pattern of urbanization index.The results show the logistic curve is good fit (R 2 = 0.93, F value = 146.444)as shown in Equation ( 17) and Figure 5.The t values for the two coefficients are 5.490 and 43.579 respectively, and both are significant at 1% level.Equation ( 17) can be transformed into Equation (18), which is used to compute the forecasted index of urbanization:

Coupling Analysis
Taking the derivative with respect to time, t, in Equations ( 17) and ( 18), we can obtain evolution equations for urbanization and water resources utilization: 2 ( ) / 0.108 0.054 0.003 (1 ) Based on Equations ( 19) and ( 20), we calculate the values of VA and VB.The variable α and tanα can be calculated using formula: α = arctan(VA/VB).The computed annual values of VA, VB, α, and tanα are shown in Table 3.  6 show that the pattern of the coupling process between urbanization and water resource utilization includes two phases.The first phase covers 2002-2003, the coupling degree was in quadrant I, −90° < α < 0°, and the two systems are in the coordinated symbiosis stage.At this stage, the level of water resource utilization declined slightly, while urbanization had just started.The interaction between the two systems was not very pronounced.
The second phase covers the years from 2004-2014, when the coupling degree was in quadrant II, 0° < α < 90°, the coordinated development stage.The coupling degree itself was rising but at a relatively slow rate.To better analyze this second phase, we further divide it into two stages.During

Coupling Analysis
Taking the derivative with respect to time, t, in Equations ( 17) and ( 18), we can obtain evolution equations for urbanization and water resources utilization: 627 × e ln29.62+ln0.627×t (1 + e ln29.62+ln0.627×t ) Based on Equations ( 19) and ( 20), we calculate the values of V A and V B .The variable α and tanα can be calculated using formula: α = arctan(V A /V B ).The computed annual values of V A , V B , α, and tanα are shown in Table 3.  6 show that the pattern of the coupling process between urbanization and water resource utilization includes two phases.The first phase covers 2002-2003, the coupling degree was in quadrant I, −90 • < α < 0 • , and the two systems are in the coordinated symbiosis stage.At this stage, the level of water resource utilization declined slightly, while urbanization had just started.The interaction between the two systems was not very pronounced.
The second phase covers the years from 2004-2014, when the coupling degree was in quadrant II, 0 • < α < 90 • , the coordinated development stage.The coupling degree itself was rising but at a relatively slow rate.To better analyze this second phase, we further divide it into two stages.During the earlier stage, years from 2004-2007, the coupling degree rose from 28.63 • to 45.218 • in which the speed of urbanization far exceeded the evolution rate of water utilization.This resulted in water shortages and water pollution problems.Water resource problems became so bad that they started to threaten the urbanization process.The "2007 Oxygenation Crisis of Taihu Lake" is a notorious example of water crisis caused by rapid urbanization.The "Oxygenation Crisis of Taihu Lake" occurred in May of 2007.The environmental pollution of Taihu Lake was caused by increasing urban residential waste water and industrial waste water (two factors in water utilization index).The severe pollution led to the cyanobacteria outbreak in the lake.Seriously polluted Taihu Lake water caused 70% of the city water pollution in Wuxi city.Sheng [86] estimated that this crisis affected the drinking water supply of two million residents in Eastern China and resulted in direct economic losses of more than $30 billion.
The second stage covers 2008-2014, the coupling degree further rose from 46.609 • to 72.143 • .During this stage, urbanization progressed at a slower pace.In addition, following the Oxygenation Crisis in 2007, people became more environmentally conscientious and the provincial government rigorously implemented an industrial restructuring program and closed many chemical plants.Both contributed to the increase in the degree of coupling between urbanization and water resource utilization.

Forecasting the Coupling Degree
The forecast accuracy of the exponential smoothing method depends on the value of α [87].Thus, we compare forecasted coupling degrees using three values of α.Table 4 shows the predicted coupling degrees and corresponding squared errors of three values of α.Given the values of α, the corresponding mean squared errors (MSE) are calculated as follows: The mean squared errors indicate that the gap between forecasts and actuals is minimized when the smoothing factor α = 0.9.Therefore we adopt smoothing factor α = 0.9.
Taking year 2014 as the first period t, we calculate the smoothing coefficient a t , b t , and c t according to Equations ( 11)-(13).Based on Equation (10) and the values of a t , b t , and c t , the forecast equation can be written as: where T represents years.The predicted coupling degrees using Equation ( 21) for the period of 2015 to 2020 are shown in Table 5.  (forecasted coupling degrees for years 2015-2020).The forecasted coupling degrees between urbanization and water resources utilization for Jiangsu province exhibit an "increase-decrease" pattern.The coupling degree is predicted to reach its highest value 74.799 • in 2016, then forecast to decline.As shown in Figure 6, the pattern of coupling degrees consists of two turning points: one in 2004 and the other for the forecast value in 2016.Prior to 2004, the coupling degree between urbanization and water resources utilization increased rapidly and stayed in quadrant I. Since 2004, the coupling degree increased gradually and is predicated to peak in 2016.All predicted coupling degrees stay in quadrant II, which means the two systems are expected to stay in the coordinated development stage during the forecast period.Thus, the demand for water resources will increase to facilitate the rapid urban, economic, and social advancements.These developments will eventually cause major problems with water shortages and water pollution.Such water crises will further hinder sustainable urbanization and aggravate conflicts between urbanization and efficient water utilization.However, as people become more environmentally conscientious, they will try to optimize the tradeoff between urbanization and water resource management.Unless such changes occur, the challenge of allowing both to achieve optimal levels will remain.Figure 6 plots data from Table 3 (the actual coupling degrees for years 2002-2014) and Table 5 (forecasted coupling degrees for years 2015-2020).The forecasted coupling degrees between urbanization and water resources utilization for Jiangsu province exhibit an "increase-decrease" pattern.The coupling degree is predicted to reach its highest value 74.799° in 2016, then forecast to decline.As shown in Figure 6, the pattern of coupling degrees consists of two turning points: one in 2004 and the other for the forecast value in 2016.Prior to 2004, the coupling degree between urbanization and water resources utilization increased rapidly and stayed in quadrant I. Since 2004, the coupling degree increased gradually and is predicated to peak in 2016.All predicted coupling degrees stay in quadrant II, which means the two systems are expected to stay in the coordinated development stage during the forecast period.Thus, the demand for water resources will increase to facilitate the rapid urban, economic, and social advancements.These developments will eventually cause major problems with water shortages and water pollution.Such water crises will further hinder sustainable urbanization and aggravate conflicts between urbanization and efficient water utilization.However, as people become more environmentally conscientious, they will try to optimize the tradeoff between urbanization and water resource management.Unless such changes occur, the challenge of allowing both to achieve optimal levels will remain.

Conclusions
This paper provides a theoretical framework analyzing the dynamic relationship between water resources utilization and the process of urbanization.We empirically examine this dynamic relationship using data from Jiangsu province in China.The results serve as an important input for developing sustainable water resource management policies during the urbanization process.
The coupling process between urbanization and water resource utilization not only reflects the complex and nonlinear interactions among various factors of two index systems, it also exhibits evolutionary phases and patterns.Based on our constructed evaluation index system, we analyze the coupling degrees of urbanization and water resources utilization in Jiangsu province.We also forecast the coupling degree of the two systems for the next six years using a cubic exponential smoothing model.Regarding the urbanization process, we find that industrial restructuring and economic development in Jiangsu province sped up the progress of urbanization.In 2002 the level of

Conclusions
This paper provides a theoretical framework analyzing the dynamic relationship between water resources utilization and the process of urbanization.We empirically examine this dynamic relationship using data from Jiangsu province in China.The results serve as an important input for developing sustainable water resource management policies during the urbanization process.
The coupling process between urbanization and water resource utilization not only reflects the complex and nonlinear interactions among various factors of two index systems, it also exhibits evolutionary phases and patterns.Based on our constructed evaluation index system, we analyze the coupling degrees of urbanization and water resources utilization in Jiangsu province.We also forecast the coupling degree of the two systems for the next six years using a cubic exponential smoothing model.Regarding the urbanization process, we find that industrial restructuring and economic development in Jiangsu province sped up the progress of urbanization.In 2002 the level of urbanization was only 0.044, but it reached 0.978 by 2013.Our analyses show the coupling degrees between urbanization and water utilization exhibit a pattern of "decline-rise-decline" for the sample period of 2002-2014.The coupling degree reached its lowest point in 2004 with the value of 0.299.Then it began to rise gradually and reached its highest value in 2010.We further divided the coupling effect between urbanization and water resource utilization into two phases.During the first phase, from 2002-2003, the two systems were at a low coordinated stage and the coupling degree stayed in quadrant I, where −90 • < α < 0 • .The second phase, 2004-2014, was characterized by the rising coupling degree and evident conflicts between the two systems.The coupling degree during phase two stayed in quadrant II, 0 • < α < 90 • .In addition, we forecast the coupling degree between these two systems.Our forecasts indicate that the coupling degree will reach its peak in 2016, and after that it starts to gradually decline.We believe that policy-makers could take our results into consideration and focus on optimizing the advancement of urbanization and sustainable management of water resources to achieve a coordinated development stage for urbanization and water resource utilization.
Possible future research could include using the theoretical framework of this study to analyze coupling degrees between urbanization and water resources utilization of other Chinese provinces.The results of such analyses would help in designing policies that coordinate urbanization and sustainable water utilization.Additional work incorporating more factors in computing the water utilization index and urbanization index would improve the quality of indices and study results.For example, adding factors to reflect the quality of urban life for the urbanization index and incorporating factors, such as ecological/environmental water usage, in computing the water utilization index.

Figure 1 .
Figure 1.Area map of Jiangsu province, China.

Figure 1 .
Figure 1.Area map of Jiangsu province, China.

Figure 2 .
Figure 2. The stages of coupling process of urbanization and water resource system.(A) The evolution path of VA; (B) The evolution path of VB.

Figure 2 .
Figure 2. The stages of coupling process of urbanization and water resource system.(A) The evolution path of V A ; (B) The evolution path of V B .

Figure 3 .
Figure 3.The indices of urbanization and water resources utilization system.

Figure 3 .
Figure 3.The indices of urbanization and water resources utilization system.

Figure 4 .
Figure 4. Cubic fitting curve for water resource utilization in Jiangsu.

Figure 4 .
Figure 4. Cubic fitting curve for water resource utilization in Jiangsu.

Figure 5 .
Figure 5. Fitting Curve of the urbanization process in Jiangsu.

Figure 5 .
Figure 5. Fitting Curve of the urbanization process in Jiangsu.

Figure 6
Figure 6 plots data from Table3(the actual coupling degrees for years 2002-2014) and Table5(forecasted coupling degrees for years 2015-2020).The forecasted coupling degrees between urbanization and water resources utilization for Jiangsu province exhibit an "increase-decrease" pattern.The coupling degree is predicted to reach its highest value 74.799 • in 2016, then forecast to decline.As shown in Figure6, the pattern of coupling degrees consists of two turning points: one in 2004 and the other for the forecast value in 2016.Prior to 2004, the coupling degree between urbanization and water resources utilization increased rapidly and stayed in quadrant I. Since 2004, the coupling degree increased gradually and is predicated to peak in 2016.All predicted coupling degrees stay in quadrant II, which means the two systems are expected to stay in the coordinated development stage during the forecast period.Thus, the demand for water resources will increase to facilitate the rapid urban, economic, and social advancements.These developments will eventually cause major problems with water shortages and water pollution.Such water crises will further hinder sustainable urbanization and aggravate conflicts between urbanization and efficient water utilization.However, as people become more environmentally conscientious, they will try to optimize the tradeoff between urbanization and water resource management.Unless such changes occur, the challenge of allowing both to achieve optimal levels will remain.

Table 1 .
Urbanization evaluation index system for Jiangsu province.

Table 2 .
Evaluation index system of water resource utilization.

Table 3 .
The coupling statistics between urbanization and water resources utilization for Jiangsu province: 2002-2014.

Table 3 .
The coupling statistics between urbanization and water resources utilization for Jiangsu province: 2002-2014.

Table 3 and
Figure
Numbers in parenthesis are squared forecasting errors.