Calculation of the contribution rate of China’s hydraulic science and technology based on a feedforward neural network

Quantitative analysis of the contribution rate of China’s hydraulic science and technology and analysis of the underlying reasons behind changes provide an important foundation upon which the government can formulate water policies. This paper abandons the assumption of a scale economy and separates the changes of benefits brought about by the scale from scientific and technological progress, thus changing the C-D production function from linear to nonlinear. Based on a feedforward neural network, it calculates the coefficient of the output elasticity, the economic contribution rate of China’s hydraulic science and technology and the scale economies for each year using relevant data from 1981 to 2016. The results show that (1) the average contribution rate of capital investment from 1981 to 2016 was 47.3%, and the average contribution rate of labor from 1981 to 2016 was 9.1%. It is not obvious that the significant increase in the labor force has contributed to the growth of China’s water conservancy industry. (2) The average contribution rate of scale economies in 1981–2016 was 26.7%, and the contribution rate of scale economies is negatively correlated with the capital contribution rate. (3) The average contribution rate of China’s hydraulic science and technology was 43.6% from 1981 to 2016, and the average contribution rate of the total factor productivity after removing scale economies from 1981 to 2016 was 16.9%. During the period of the 6th Five-Year Plan(1981~1985), the contribution rate of water conservancy science and technology was relatively high. Since that time, it has remained at 40%. In recent years, as water conservancy reforms in key areas have made positive progress, scientific and technological progress has increased the growth of water conservancy benefits annually.


Introduction
In recent years, water resources in China have faced practical problems of flooding, droughts, water quality pollution, and soil loss [1], and the government of China has stated that the PLOS  improvement of the total factor CO 2 emission performance relationship through scientific and technological progress. The production function and water price [31,32] have also been studied, and it has been shown that different production functions should be selected to better reflect reality [33]. The current method for measuring the contribution rate of water conservancy science and technology in a monopoly industry is mainly based on the C-D production function and the Solow residual value method [34,35]. Kong [36] used solow model to calculate the contribution rate of water conservancy science and technology in China, and predicted the future development of hydraulic technology [37]. The values of the capital-output elasticity coefficient α and the labor output elasticity coefficient β in the model strongly influence the final result of the contribution rate of water conservancy science and technology [38]. At present, there are some disadvantages to estimating model parameters using the empirical value method [39] or the neoclassical growth theory, the empirical value method is subjective, and the neoclassical growth theory lacks a solid foundation in reality [40]. The neoclassical growth theory has two core assumptions: one is exogenous technology, and the other is the constant return of the scale of production. Therefore, if there is a scale economy effect, the result of the water conservancy science and technology contribution rate will be distorted [41]. In this paper, we abandon the second assumption of the neoclassic growth theory, and separate the changes in benefits brought about by the scale from scientific and technological progress, thus changing the model from linear to nonlinear. Based on the partial derivative of the feedforward neural network output function, the partial derivative is used to fit the production function.
There are some limitations in using a multilayer feedforward network. The number of hidden layer nodes in neural network is related to the complexity of mapping relation. The fewer the number of hidden layer nodes, the more difficult the network is to describe complex problems. However, too much number will lead to the increase of training time of neural network, and too many networks will overfit. Therefore, when using neural network to calculate the contribution rate of water conservancy science and technology in China, it is necessary to choose the appropriate network structure to get the ideal result. Another defect of neural networks is local optimum, Most of the current algorithms can be divided into two categories: gradient descent and numerical computation. It is difficult to compare which algorithm is better in the actual situation, and researchers are often required to conduct simulation experiments.
Generally speaking, neural network is more complex than traditional solow residual method, and the network and algorithm need more elaborate design. However, through the neural network algorithm, the production function can be made more flexible, and the effect of economies of scale on the overall return can be well calculated. We abandon the unreasonable assumptions in the model and make the model closer to reality. In this paper, we design an algorithm to overcome. From this algorithm, the output elasticity of production factors is calculated, and the contribution rate of water conservancy science and technology is estimated for each year. The deep-seated reasons behind the changes are analyzed to provide a reference for the government to use when formulating water conservancy policies and long-term control measures.
This paper was organized as follows: Section 2 described the data and calculation method in detail. Section 3 set the key parameters of the neural network and described calculation results. Section 4 compared the results to those using another production function, and then discussed the development direction of water conservancy technology in China, summarized the concept and calculation method of the total factor productivity, and presented the future prospects for research on this topic.

Calculation of the annual water conservancy benefits in China.
The scope of the water conservancy benefits includes six major benefits: the irrigation benefit, hydropower benefit, flood control benefit, water supply benefit, soil and water conservation benefit, and water logging control benefit.
The irrigation benefit is mainly reflected in the increase in the yield and output of agricultural products due to irrigation. The increased agricultural output is the result of a combination of irrigation and agricultural measures and is deduced by the apportion coefficient of the water conservancy irrigation benefit.
The direct financial benefit of hydropower is obtained by multiplying the electric energy production (after deducting factory power consumption and line losses) by the electricity price. It is calculated by multiplying the hydroelectric energy production in China by the national average electricity price for water conservancy systems annually.
The benefits of new flood protection projects that reduce flood damage include the sum of the direct economic losses of various sectors of the national economy, including industry, agriculture, commerce, transportation, construction, and property in inundated areas. The flood control benefit is obtained by multiplying the area experiencing disaster reduction due to the flood prevention project by the total loss in terms of the Mu unit price.
The water supply benefit refers to the urban water supply benefits from water conservancy projects. After utilizing the input-output coefficient obtained from the weighted average of the water supply investment and water supply benefits over the years, the water supply benefits can be calculated according to the national annual investment in water supply projects.
The soil and water conservation benefits mainly include direct income from agriculture, forestry, animal husbandry, and fishing in addition to other social and economic benefits.
Waterlogging control benefits refer to the increased production of primary crops after building a waterlogging control project. This benefit generally reflects the national waterlogging control benefit, which is calculated according to the annual waterlogging control area and the waterlogging control benefit per acre.
The retail price index is selected for analysis, the economic benefits of China's water conservancy system are shown in Table 1 Table 2.
Labor input estimation: The contribution of labor to the water conservancy system, excluding labor by those engaged in work with earth and stone, should be attributed to the labor input. According to the relevant data from the "China Water Conservancy Statistical Yearbook" for past years, the annual labor input of the national water conservancy system is calculated by taking the average between the beginning and ending dates of the year, as shown in Table 2.

Production function selection
Among the many production functions, the most commonly used production functions are the Cobb-Douglas production function, linear production function, Leontief production function [42], constant elasticity of substitution production function [43] and translog production function [44], as shown in Table 3. The main differences between the production functions lie in the assumption of elasticity of substitution, and the values of the elasticity of substitution of these production functions are shown in Table 4.
According to the characteristics of each production function, we can obtain some conclusions: 1. The linear production function describes a production model with a constant return on scale, which can be used in economic analysis, including constant economies of scale. In Table 2. Capital and labor input to the water conservancy system over time.

Year Fixed asset investment price index
Fixed assets formation (current prices) this paper, we abandon the assumption that the economies of scale of the CD production function are constant.
2. The constant elasticity of the substitution production function (CES) improves the elasticity of substitution of the Cobb-Douglas production function. For different research objects or different sample intervals of the same research object, the elasticity of substitution is different because of the difference in sample values. The substitution elasticity of the constant elasticity of the substitution production function is 1 1þr , while the substitution elasticity of the Cobb-Douglas production function is 1, and this makes the constant elasticity of the substitution production function more realistic than that of Cobb-Douglas production function.
3. The elasticity of substitution of the Leontief production function is zero, and this shows that there is no substitution between input factors.
4. The translog production function can be converted to any form of second-order Taylor approximation of the production function, and it can be used to test whether the elasticity of substitution is constant. The translog production function has good mathematical properties, but it requires complete panel data, so the transcendental logarithm function is not used for comparison in this paper.
Based on the data in this article, we used a more general C-D production function that considers scale economy effects for analysis, and we also used CES production functions to compare the results of the C-D production function. The methods of CES production functions are shown in the S1 File. Table 3. Commonly used form of the production function.

Production function name Production function form
Cobb-Douglas production function Constant elasticity of substitution production function(CES) b nm lnx n lnx m Þ https://doi.org/10.1371/journal.pone.0222091.t003 Table 4. Commonly used forms of production functions.

Production function name Elasticity of substitution
Cobb-Douglas production function 1 The C-D production function and Solow residual value method assume that the sum of the output elasticities of production factors is 1. If there is no scale economy, this assumption will attribute the changes in benefits brought about by the scale change to scientific and technological progress. The C-D production function and Solow residual method do not consider the impacts of scale economies, which will lead to errors in the results of the contribution rate of scientific and technological progress.
This paper does not assume scale economies. The production function: Let capital productivity P 1 = Y/K and labor productivity P 2 = Y/L, Convert Eq (1) to: In Eq (2), Change Eq (1) to a differential form: ð3Þ α t and β t are the output elasticities of capital and labor in year t. The first item in Eq (3) is the Solow residual value produced after eliminating the scale economy. The second item is the output change (s t ) brought about by the scale economy effect. The sum of the first two items constitutes the total factor productivity (TFP). The third item is the output changes brought about by the change in production factors. Through this method, the increased output brought by scale economies can be removed from the TFP.

The C-D production function calculation method
2.3.1 The challenges of the C-D production function calculation in a multilayer feedforward network. Next, the effect of the scale economy on the output growth must be calculated based on the Eq (3). When using multivariate regression analysis, which requires a longer time series, it is easy to miss changes in capital, labor, and technology for the current year. Additionally, in this paper, if α and β are constants, the output change (s t ) brought about by the scale economy effect will be equal to zero, which will be contrary to the original intention of this paper. Therefore, ordinary least squares(OLS) cannot be used.
A multilayer feedforward network can describe extremely complex nonlinear systems and handle the random characteristics of the original data. A multilayer feedforward network is a universal approximator that can automatically simulate a system when there is sufficient training time.
In the structural design of a feedforward network, the node scale should be reduced as much as possible to reduce the system complexity and learning time while still meeting the requirements. The number of nodes is related to the complexity of the object to be expressed in the model. If the number is too small, the model struggles to describe real problems. If the number is too large, however, the degree of nonlinearity becomes higher, the training time increases, and the model becomes prone to the "overlearning" phenomenon. As a result, the ability of the model to generalize deteriorates. Therefore, the principle used to select the number of nodes in this paper is based on the correct reflection of the relationship between the input and output, and fewer hidden layer nodes are used to make the network structure as simple as possible.
The multilayer feedforward network uses a single-track multilayer structure in which each layer contains several nodes. The nodes in the same layer are not connected and the information between the layers is transmitted in only one direction. A single-layer feedforward network can only distinguish linearly separable modes while a multilayer feedforward network can be used for any classification problem. In addition, a multilayer feedforward network can be used as a general function approximator. A two-layer feedforward network can approximate any complex continuous function as long as there are enough nodes in its hidden layer with an s-shaped activation function of the hidden layer and a linear activation function of the output node. Therefore, this paper selects a two-layer feedforward network with a hidden layer as the research method.
In the research process, we find that there are two main problems of the BP algorithm: the slow convergence speed and the local minimum point of the function. To resolve the slow convergence speed, a network pruning algorithm is used to delete redundant input nodes. At the same time, based on the sensitivity analysis, we analyze the influence of parameters on the system, the influence of small disturbances in the parameters on the network output, and the importance of the network parameters. When the function has a local minimum point, the adaptive variable step length algorithm is selected, which means that the learning compensation adjusts with the variation in the error surface.

1) Sensitivity Analysis
Sensitivity analysis is also called perturbation analysis. It is an analytical method used to determine the sensitive factors and their degree of influence. The sensitivity function S (y i , x j ) indicates how the change in the system input x j affects the system output y i . There are four calculation methods, as follows: Normalized sensitivity: Seminormalized sensitivity by y i : Seminormalized sensitivity by x i : Unnormalized sensitivity:

2) Configuration Optimizing Theory
Sensitivity analysis is used to define the sensitivity of nodes, and then the relative total sensitivity is defined. The relative total sensitivity index is used to characterize the influence degree of the parameters on the model output. The input node, hidden layer node, and connection weight are selected using the relative total sensitivity index, and then the parameters with little or no effect can be deleted.
Definition 1 For mode u, x ðuÞ j represents the input of the j th input node, and y ðuÞ i represents the output of the i th output node. The sensitivity of the network output function F(x, w) with respect to the input node can be defined as: AS ij is the relative sensitivity index of output i to input j. The weight of each output variable in the system is not exactly the same: AS ij must be weighted. λ ij is the weight of AS ij , and satisfies ∑ j λ ij = 1, λ ij � 0, i = 1, 2, . . ., n.
In the training process, some nodes have a strong network correlation but represent small proportions. The average sensitivity AS ij calculation method may delete nodes. Further, by using the average calculation method, training will be reduced, but a stronger correlation will be suppressed. Therefore, we need to calculate the variance of the node to be deleted s 2 j ¼ Varððts ð1Þ j ; ts ð2Þ j ; . . . ; ts ðsÞ j Þ T Þ to make further judgments.

3) Sensitivity Calculation
To carry out sensitivity analysis, @y ðuÞ i @x j needs to be calculated. Take the derivative of Eqs (1)~(3): The training set instance x (u) can calculate the partial derivative of the network output function F(x, w) with respect to the input x. If u = l, 2, . . ., s, s partial derivative matrixes can be obtained by repeated calculations, and the elements in the matrix will be partial derivatives that are shaped like formula (10). If expressed by a matrix symbol, Eq (10) can be written as: J ðuÞ x is the Jacobian matrix of the network output function F(x, w) at x (u) ; W 1 is the connection weight matrix of the first layer; W 2 is the second-level connection weight matrix; The sensitivity of the hidden layer node, relative sensitivity, and partial derivative can be similarly obtained.

4) Calculation of the output elasticity coefficient
The econometric model of a multilayer feedforward network can be expressed as: It can be obtained from Eq (12) . . . ; x ðtÞ m Þ ¼ A t K a L b , where A t , α, and β are parameters that change with time.

5) Key parameters of production function fitting
The input data are normalized to speed up the network convergence. Three datasets are selected as the test set, and the remaining 33 datasets are used as the training set. The initial network structure can be taken as n-l-m, and when the RTSI j of a hidden layer node is smaller than a ¼ 1 10l (where l is the number of hidden layer nodes), the node should be deleted. The training goal is to achieve a mean square error (MSE) of 10 −6 , which can ensure the robustness of the model. The code can be found in S2 File.

Parameter settings of the neural network
This paper establishes a multilayer feedforward network model. After training MATLAB's neural network toolbox, the output function is obtained in the form of Q(t) = A t K α L β .
The input data are normalized to speed up the network convergence. Three datasets are selected as the test set, and the remaining 33 datasets are used as the training set. The training goal is to achieve a mean squared error (MSE) of 10 −6 . The maximum number of cycles is 5000.
The initial network structure is taken as 5-10-1, and when the RTSI j of a hidden layer node is smaller than a ¼ 1 10l (where l is the number of hidden layer nodes), the node should be deleted. The network structure optimization algorithm is applied to obtain the simplest network structure 3-4-1. The MSEs of the training set and the test set are 1.0264×10 −6 and 1.0057×10 −6 . The connection weight matrix and threshold are W 1 = [-1.0245-2.9913] W 2 = [-1.6944-9.2483] b 1 = -0.1367 b 2 = -6.9678.
When the multilayer feedforward network fits the production function, because the input and output data are normalized, the partial derivative needs to be multiplied by a corresponding coefficient to obtain the estimated value of the partial derivative of the potential function. Then the output elasticities α and β of the production factors are found, as shown in the S3 File.
To determine the contribution rate of scientific and technological progress, the following indicators need to be calculated.
First, the average annual growth rate y, k, l should be calculated. Because the scope of flood disaster reduction and the extent of disasters in water conservancy projects are also random variables affected by climatic factors, benefits such as flood control and water logging control show greater volatility interannually. Therefore, the horizontal method is used to calculate the output. Taking the output as an example: where Y t is the output efficiency of the water system during period t. According to Eq (3), the contribution rate of the Solow residual value to the total output growth rate EA t is: The contribution of the scale economy to the total output growth ES t is: The contribution of the capital input to the total output EK t is: The contribution rate of the labor input to the total output EL t is: The calculation results are shown in the S3 File.

The growth of benefits of water conservancy according to the input factor
The contributions of different input factors to the economic benefits of water conservancy over the years show different trends. The water conservancy industry is a traditional infrastructure industry. The average contribution rate of capital input from 1981 to 2016 was 47.3%, and the average contribution rate of labor input from 1981 to 2016 was 9.1%.
During the process of reform and opening up, the scale of capital input was relatively small, but as time passed, the water conservancy investment gradually increased. From 1981 to 1993, the contribution rate of capital investment rose from 14.0% to 60.8%, and then it declined slowly. After the outbreak of the international financial crisis in 2008, the central government of China implemented a proactive fiscal policy and increased water conservancy infrastructure construction. As a result, the contribution rate of capital investment rose again and then declined, and the capital contribution rate was 37.1% in 2016.
As shown in Fig 1, the contribution rate of the labor input remained relatively stable between 1981 and 2016. It is not obvious that the significant increase in the labor force contributed to the growth of China's water conservancy industry because laborers in China are not productive. Taking a long-term perspective, China's economy must enter an intensive growth period and policy makers should pay attention to the utilization of human resources.

The growth of the benefits of water conservancy based on the scale economy
The average contribution rate of the scale economy from 1981 to 2016 was 26.7%. As shown in Fig 2, the scale economy contribution rate has a negative correlation with the capital contribution rate.
From 1981 to 2007, the scale economy contribution rate was relatively stable. After the reform and opening up, the invisible hand of the market economy gradually played a role, and the scale economy contribution rate has increased; thus, the capital investment can bring multiple investment returns, and there has been less water conservancy industry duplicate construction during this period. After the implementation of the proactive fiscal policy in 2008, the central government of China implemented a large amount of investment in water conservancy construction, which led to a certain amount of duplication of construction, and the scale of the economic contribution declined and then returned to normal levels.

Scientific and technological progress in the growth of water conservancy
The average rate of the scientific and technical contribution to water conservancy from 1981 to 2016 was 43.6%. After excluding the economies of scale, the average contribution rate of the TFP was 16.9%.
During the period of the 6th Five-Year Plan (1981~1985), although the scale of capital investment was small, the national scientific and technological research plan within the 6th Five-Year Plan heavily subsidized water conservancy. The implementation of the "Study of the South-to-North Water Diversion Project," "Technology development of large-scale hydropower stations", "Studies on the Development and Utilization of Water Resources and Comprehensive Evaluation in East China" and other projects quickly compensated for technical problems in traditional infrastructure construction and became the main impetus promoting the growth of water conservancy benefits in the early stages of reform and opening up.
Since then, the contribution rate of water science and technology has remained at 40%. With the implementation of strict water resource management systems and the construction of large-scale water conservancy projects, water conservancy reforms in key areas have made positive progress. Scientific and technological progress has increased the growth of water conservancy benefits annually. In 2015, the contribution rate of science and technology was 53.6%, compared with 55.7% in 2016.
As shown in Fig 3, Solow's residual value has been decreasing since 1981. The growth of water conservancy emphasizes the formation of physical capital. If the capital allocation is not reasonable, the high capital accumulation rate will compensate for the economic operation inefficiency. The Chinese economy has shown extensive features over a period of time, but in the last ten years, the TFP has again increased, and economic growth has shifted from the extensive mode to the intensive mode.

Comparison of the results of the CES production function
There is no research on the contribution rate of water conservancy science and technology in China, so we choose the calculation results of the CES production function for a comparative analysis.
As shown in the S3 File, ρ is approximately equal to zero, and this proves that selecting the Cobb-Douglas production function can basically meet the research needs. Then we compare the contribution rates of the capital, labor, total factor productivity, economies of scale and Solow residual value as calculated by the Cobb-Douglas production function and the constant elasticity of substitution(CES) production function.
As shown in Fig 4, the capital contribution rate of the CES production function is relatively stable compared with that of the Cobb-Douglas production function. However, the capital contribution rate of the CES production function from 2008 to 2011 was higher than that of Cobb-Douglas production function. The labor contribution rate of the CES production function is consistent with that of the Cobb-Douglas production function. The total factor productivity of the CES production function has the same trend as that of the Cobb-Douglas production function, but it is lower than that of the Cobb-Douglas production function from 2008 to 2011. The contribution rates of capital and the total factor productivity as calculated by the CES production function are still inversely proportional.
As shown in Fig 5, the scale economy contribution rates and Solow residual values calculated by the CES production function and Cobb-Douglas production function are quite different in some years. The scale economy contribution rate of the CES production function was negative from 2008 to 2012. Although the scale economy contribution rate of the Cobb-Douglas production function also declined in 2008, the decline was not as significant as the scale economy efficiency of the CES production function. Because the CES production function has a scale economy coefficient m, it has an advantages over the Cobb-Douglas production function in calculating the contribution rate of the scale economy. From 1981 to 2007, the scale economy contribution rates calculated by the Cobb-Douglas production function and the CES production function were roughly the same. Because the two production functions have different ways to calculate the scale economy contribution rate, the two Solow residual values have some deviations, and this results in a certain deviation of the two Solow residuals calculated by the two production functions. Especially after 2008, the Solow residual value of the CES production function fluctuated greatly, which is quite different from the case for that of the Cobb-Douglas production function.

Prospects for the development of water conservancy science and technology in China
The development of water conservancy science and engineering technology come from the needs of the survival and development of human society. Its goal is to understand the laws of nature and to manage and allocate water resources artificially through a variety of technical and engineering measures. However, water security in China is still facing a serious situation, which is highlighted in the serious shortage of water resources, water pollution, the threat of water disasters, water ecological degradation and other aspects. At the same time, superstrong earthquake, excessive flooding huge geological disasters and other factors threaten the long-term safe operation of hydropower stations [45]. These problems endanger the downstream people's lives and property safety. In terms of science and technology, there are still many major problems that need to be solved. 1) Water resources security [46]: the prediction of the supply and demand of water resources affected by climate change in the future; the inter-basin spatial and temporal allocation of water resources; and the strictest water resources management technology. 2) Watershed water, sediment and environmental ecology [47]: the prediction of future river water, sediment, and ecological environment changes under the conditions of industrial and agricultural development, urbanization and largescale water conservancy project construction; the interaction of the fluxes of river water and sediment with the ecological environment; the hydrodynamics of river systems; the balance and regulation of geomorphology and biodiversity. 3) Hydropower energy development and long-term safe operation guarantee [48]: the effects of extreme disaster factors such as super earthquakes, excessive flooding, and complex geological conditions on the hazard chain risk of the high dam junction group of hydropower stations. 4) Flood and drought disaster prevention [49]: the disaster mechanism and risk control of river flood, urban floods, mountain mudslides, storm surges, and drought under the influence of extreme meteorological conditions and human activities.
The first, second, and third scientific and technological problems require the construction of a large number of water conservancy facilities such as dams, water diversion channels and so on. The fourth scientific question is inclined to the study of earth system science and the water circulation mechanisms. As shown in Fig 5, the efficiency of the scale economy in China declined significantly in 2008. As shown in Fig 2, there is a certain inverse relationship between capital and scale economies. As shown in Fig 6, the contribution rate of capital and the science and technology of water conservancy in China are basically the same, and in recent years, the contribution rate of science and technology has risen slowly. The policy of the Chinese government has gradually shifted from building water conservancy to harmonious coexistence between man and nature. This policy will minimize the construction of large water conservancy projects. It is expected that the future investment in water resources will gradually turn into the mechanical study of the water cycle, in which technology will play a further leading role.

The concept of the total factor productivity
The concept of the total factor productivity in this paper is Solow residual value, which uses the output growth rate after deducting each input factor growth rate to measure the total factor productivity. There is another definition of the total factor productivity, which refers to the ability to achieve the maximum output under a given input or the ability to minimize the input under a given output. This scenario is considered an example of Pareto optimality [50]. Under this definition of the total factor productivity, we use the DEA [51] and other methods [52] for the calculation. The calculation method and sphere of application are summarized in Table 5.
The calculation method of the DEA or SFA is often used to compare the total factor productivities of different departments, which requires complete panel data. The C-D production function is used to calculate the total factor productivity of the entirety. The two types of calculation methods have different definitions of the total factor productivity. However, both methods can reflect the efficiency of the research object's use of input factors. In follow-up study, we can combine a non-parametric method (DEA, SFA) with a parametric method, and unify the definition of the total factor productivity. At the same time, we can consider pollution emission from different section to further expand the comprehensiveness of the total factor productivity.

Conclusions
This paper establishes a C-D production function that considers scale economy effects and uses a feedforward neural network to fit the potential production function with historical data. It calculates the coefficient of the output elasticity and then determines China's TFP from 1981 to 2016. The conclusions are as follows: 1. The contributions of different input factors to the economic output of water conservancy over the years have shown different trends. The average contribution rate of the capital input from 1981 to 2016 was 47.3%, and the average contribution rate of the labor input was 9.1%. During the process of reform and opening up, capital investment in water conservancy increased from 14.0% to 60.8% between 1981 and 1993 and then declined slowly. After the international financial crisis in 2008, China implemented a proactive fiscal policy and increased water conservancy infrastructure construction, leading the contribution rate of capital investment to rise again and then fall. The contribution rate of capital in 2016 is 37.1%. It is not obvious that the significant increase in the labor force has contributed to the growth of China's water conservancy industry, because the laborers in China are not productive.
2. The average contribution rate of the scale economy from 1981 to 2016 was 26.7%, and the scale economy contribution rate was negatively correlated with the capital contribution rate. From 1981 to 2007, the contribution rate of the economies of scale was relatively stable. After the implementation of the proactive fiscal policy in 2008, the contribution rate of economies of scale decreased, and the efficiency of capital utilization fell further before returning to normal levels.
3. The average rate of the scientific and technical contribution of water conservancy from 1981 to 2016 was 43.6%. After excluding economies of scale, the average rate of the total factor contribution from 1981 to 2016 was 16.9%, and the contribution rate of water conservancy technology in 2016 was 55.7%. During the period of the 6th Five-Year Plan (1981~1985), the contribution rate of water conservancy science and technology is relatively high. Since that time, it has remained at 40%. In recent years, as water conservancy reforms in key areas have made positive progress, scientific and technological progress has gradually increased the growth of water conservancy benefits.
Supporting information S1 File. The CES production function.