Next Article in Journal
Signal Spectrum Analysis of Sediment Water Impact of Hydraulic Turbine Based on ICEEMDAN-Wavelet Threshold Denoising Strategy
Previous Article in Journal
Wastewater Treatment System Optimization for Sustainable Operation of the SHARON–Anammox Process under Varying Carbon/Nitrogen Loadings
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Mechanism Analysis and Demonstration of Effective Information Extraction in the System Differential Response Inversion Estimation Method

1
Jiangsu Maritime Institute, Nanjing 211199, China
2
College of Hydrology and Water Resources, Hohai University, Nanjing 210098, China
3
Beijing IWHR Technology Co., Ltd., Beijing 100038, China
4
Changtan Reservoir Affairs Center, Taizhou 318020, China
5
ShaoXing Tanpu Reservoir Co., Ltd., Shaoxing 312000, China
*
Author to whom correspondence should be addressed.
Water 2023, 15(22), 4016; https://doi.org/10.3390/w15224016
Submission received: 14 September 2023 / Revised: 23 October 2023 / Accepted: 27 October 2023 / Published: 19 November 2023
(This article belongs to the Section Hydrology)

Abstract

:
The system differential response method for inverse estimation has received much attention in the hydrology literature. However, its underlying mechanisms remain largely unexplored, highlighting the need for this study. This study proposes the relation degree coefficient (RC) concept describing the nonlinear relationships between different variables and demonstrates that the selective information extraction ability of the method was first demonstrated from a theoretical perspective. Synthetic cases were conducted to demonstrate the method performance under various variables to be estimated. The results show that the useful information is extracted from the relationship between the variables to be estimated and the observed discharge. In addition, there is a general trend that suggests that incorporating more variables into the inversion estimation can lead to enhanced estimation performance.

1. Problems of the Optimization Based on Objective Functions

Maximizing effective information utilization is the basis of model establishment, optimal parameter search, and accurate simulation and prediction [1,2,3,4]. The flow information at the outlet of a river basin can be seen as the information treasure trove of hydrological modeling research because it contains all the information that impacts the calculation results of the model. The inversion estimation method of system differential response proposed by our research group provides an effective method for the selective utilization of flow information. It has been widely verified in real-time flood forecasting error corrections, model parameter calibration, and error estimations of rainfall, runoff yield, and initial soil water content. To deeply understand and analyze the theory and mechanism of the system differential response method, this paper analyzes and demonstrates the basic theoretical issues of the extraction effect and selective extraction ability of the inversion estimation information of the system differential response.
Model parameter calibration is the earliest application of the system differential response inversion estimation method [4,5], which has attracted extensive attention [6,7,8,9,10,11]. Calibration methods of model parameters are conventionally based on the responses of the objective function and use the provided information to search for the optimal parameters. These methods have two essential problems [4]. The first problem is that the objective function, e.g., the error square operation, not only complicates the surface of the objective function but also changes the original optimal values and induces irrelevant local optimal values. We illustrate this problem through a parameter estimation example (i.e., estimating a from the following nonlinear function):
y = a x + a 2 sin ( π 2 x )
Let us take samples of (x = 1, y = 2) and (x = 2, y =2). For the first sample, parameter a has two possible local optimal solutions, a1 = 1 and a2 = −2. There is only one solution for the second sample, a1 = 1. Parameter a should satisfy both samples; hence, only a1 = 1 is the unique solution of Equation (1). We can use the objective function of the sum of squared errors of the two samples, as shown in Equation (2).
F ( a ) = ( 2 a a 2 ) 2 + ( 2 2 a ) 2
The local optimal solutions can be obtained by formatting a cubic equation by the derivation dF (a)/da = 0
6 a 3 a 2 - 2 a 3 = 0
In theory, the above equation has three local optimal solutions, and interestingly, −2 is not the local optimal solution of Equation (3). This shows that the square error operation not only induces uncorrelated local optimal solutions (from 2 to 3) but also changes the solution of the original function.
The second problem is that the sum of error squares of the objective function greatly loses the effective information in determining parameters. Similarly, we take Equation (1) as an example without losing generality. If the error squares are not aggregated, take each error square as the individual objective function, as follows:
F 1 ( a ) = ( 2 a a 2 ) 2 F 2 ( a ) = ( 2 2 a ) 2
Similarly, with the derivation operation dF1 (a)/da = 0 and dF2 (a)/da = 0, Equation (5) is obtained
( 2 a a 2 ) ( 1 + 2 a ) = 0 2 ( 2 2 a ) = 0
The local optimal solutions of Equation (5) are a1 = 1, a2 = −2, and a3 = −0.5. This shows that if the error squares are not summated, the original two local optimal solutions a1 and a2 can still be found. But after summation, a2 cannot be found, as shown above. This implies a significant loss of information on the aggregation operations commonly used in optimization issues. If it is a complex nonlinear system such as the watershed hydrological model, information loss problems, the change in local optimal values, and the increase in uncorrelated local optimal values can be more pronounced using the objective functions based on the sum of error squares. Therefore, as a step forward, the method of system differential response proposed by our research group extracts parameters directly from the surface of the parametric function. This method does not lose the amount of information and does not increase the local optimal value, and the search speed is very fast. For more information, see papers [4,5,6,7].
Flood forecasting error correction is another important application field of the method. Hydrologic models are generally designed to model the rainfall–runoff physical process which is generally considered to be nonlinear and time-varying. Conceptual rainfall–runoff models simplify and conceptualize these complex processes using a set of simple mathematical equations. Conceptual rain–runoff models are generally reported to be robust and reliable in the applications of flood forecasting. However, it is generally accepted that the satisfactory application of such a model can be hampered by many factors including errors in the inputs, inadequacy of the model which generally includes errors in the model structure, errors in the model parameters, and errors in the model initial conditions. Many studies have indicated that sequential data assimilation methods represented by particle filtering and the ensemble Kalman filter, as well as variational assimilation methods exemplified by the 4D-variational method, can effectively reduce the uncertainties of flood forecasting. However, the successful application of these algorithms heavily relies on statistical assumptions that align with real-world conditions. In practice, it is often challenging to provide precise statistical assumptions. For instance, estimating the background error covariance in variational methods is difficult. Assumptions are typically made based on empirical evidence. If the given statistical assumptions deviate significantly from the actual conditions, the effectiveness of these methods might be compromised, sometimes even leading to divergence in results.
The method is applied to estimate the errors existing in rainfall [12], runoff yield [13], and soil water content [14], and the results demonstrate that the method can significantly improve the forecasting accuracy. For example, the proposed method outperforms the autoregressive technique when applying the methods to improve the forecast performance of a lumped hydrological model [15]. As the driving force input of the hydrological model, the error in rainfall is critical to the model’s accuracy. There are two types of rainfall errors: the observation error caused by instruments and equipment of measuring points and the error caused by obtaining the areal rainfall from the measured point rainfall. In arid and semi-arid basins, the second type of error is often the main error source of the model [12,13,16]. The dynamic system inversion estimation method of differential response is used to dynamically estimate the rainfall error, which indirectly modifies the flood prediction results and greatly improves the flood prediction accuracy [12,13]. Later, many scholars put forward improved methods accordingly [14,15,16,17]. Rainfall is not only the main error source of the hydrological model but also the main source of error in the integrated water–sediment model [16]. The static system differential response inversion method is further used to estimate the total rainfall of each flood or rainfall error [16] to reduce the influence of rainfall error on the flow and sediment model simulation results. The error of calculated runoff is one of the main sources of model calculation error [13].
The inversion estimation method of system differential response is used in the estimation of runoff error [13,18,19], and a variety of improved methods have been derived [20,21,22,23,24], which has received much attention in the context of hydrological modeling [25,26,27,28,29,30]. Other applications include the inversion correction of initial soil water content error [14], the inversion correction of sediment yield error [31], and the error correction of comprehensive multi-factors [15]. However, these applications are mainly demonstrations of the effect of the method, and some theoretical problems of the method itself are rarely explored. In this paper, the mechanism of how the system differential response extracts effective information from the flow hydrograph is expounded and demonstrated theoretically.

2. The Inverse Estimation Method of System Differential Response and Theoretical Demonstration of Its Effective Information Extraction

2.1. Basic Method

A hydrological model can be expressed by the following expression:
Q t = f ( Φ ; t )
where Q is the model output; Φ is a vector composed of model parameters, inputs, intermediate variables, and initial values to be estimated; and t is a discrete-time auxiliary variable. According to the inversion estimation method of system differential response [4,5], the cyclic estimation Formula (7) can be used to calculate the Φ to be estimated:
Φ ( j + 1 ) = Φ ( j ) + ( S T S ) 1 S T Δ Q
where j is the number of estimation cycles, Δ Q T = [ Q 1 f ( Φ ( j ) ; 1 ) , Q 2 f ( Φ ( j ) ; 2 ) Q n f ( Φ ( j ) ; n ) ] is the deviation of the model simulation, and S is the system differential response matrix. Equation (7) can be resolved iteratively by minimizing the sum of squared deviations of the model simulation. That is, the iterative estimation makes Equation (8) satisfy Equation (9):
F ( j ) = t = 1 n ( Q t f ( Φ ( j ) ; t ) ) 2
F ( j ) F ( j + 1 ) j min
Therefore, the system differential response method is convergent and optimal.

2.2. The Relation Degree Coefficient and the Information Extraction Mechanism

For the convenience of demonstration without losing generality, it is assumed that the model has only one factor P to be estimated. According to the inversion estimation method of system differential response (Equation (7)), P is estimated as Equation (10).
P ( j + 1 ) = P ( j ) + t = 1 n ( Q t f ( P ( j ) ; t ) ) f ( P ; t ) P P = P ( j ) t = 1 n ( f ( P ; t ) P P = P ( j ) ) 2
where f ( P ; t ) P P = P ( j ) t = 1 n , when P = P ( j ) , is the system differential response curve series, P (j) is known at step j, and P (j + 1) is a new value to be estimated from the system differential response. In analogy to the concept of the correlation coefficient in statistics, Equation (10) is transformed to Equation (11) to demonstrate the mechanism of information extraction of system differential response.
Δ P = t = 1 n ( Q t f ( P ( j ) ; t ) ) f ( P ; t ) P P = P ( j ) t = 1 n ( Q t f ( P ( j ) ; t ) ) 2 t = 1 n ( f ( P ; t ) P P = P ( j ) ) 2 t = 1 n ( Q t f ( P ( j ) ; t ) ) 2 t = 1 n ( f ( P ; t ) P P = P ( j ) ) 2
The relation degree coefficient (RC), remaining information (RI), and unit information (UI) are then defined as follows.
R C = t = 1 n ( Q t f ( P ( j ) ; t ) ) f ( P ; t ) P P = P ( j ) t = 1 n ( Q t f ( P ( j ) ; t ) ) 2 t = 1 n ( f ( P ; t ) P P = P ( j ) ) 2 R I = t = 1 n ( Q t f ( P ( j ) ; t ) ) 2 U I = t = 1 n ( f ( P ; t ) P P = P ( j ) ) 2
It can be seen from Equation (12) that if the two series ( Q t f ( P ( j ) ; t ) ) ; f ( P ; t ) P P = P ( j ) t = 1 n have a zero mean, RC is equivalent to the correlation coefficient. Therefore, RC can be understood as the coefficient reflecting the correlation degree between the series of differential response curves of the system and the series of deviations calculated by the model. Moreover, it can be proved mathematically that the range of the RC value is between −1 and 1, similar to that of the correlation coefficient. The residual information RI is expressed by overall differences between the observed flow and the flow calculated by the model and thus usually contains both effective and harmful information. The harmful information should not affect the inversion estimation result of P. Therefore, RC can be used to filter out harmful information. UI is the amount of effective information needed to estimate a unit variable, which is numerically equal to the sum of the system differential response curve squares. It can be seen that RC expresses the correlation between the differential response curve and the deviation series of the model calculation, suggesting that the method extracts the effective information needed for estimation through the correlation between the differential response curve and the residual information series. The differential response curve of the system varies with different variables to be updated, and the statistical characteristics of the associated residual information series also differ. This achieves the distinction and extraction of the necessary information for updating different factors.

3. Analysis and Demonstration of Information Extraction Effect

It is sometimes difficult to demonstrate the rationality of a new method from the mathematical perspective, and hence, synthetic cases are used for in-depth analysis in this study.

3.1. Synthetic Case

Because the true model structure, parameters, inputs, outputs, and all related errors are unknown for real-world applications, we conduct a synthetic case using a conceptual hydrological model, where all the true values are known (predefined). Therefore, the performances of the method can be directly verified and compared by evaluating the accuracy of the updated variables. In other words, in the synthetic case, we assume that the model structure, parameters, inputs, outputs, and all related errors are known. In other words, the model structure is set, the model parameters and input are given, and the model output is calculated according to the input. As such, the error-free synthetic case is constructed. Then, a series of random errors are generated and superimposed into the “true” value, according to the need for demonstration, to form a synthetic case with an error closer to reality. In this study, the random number series with uniform distribution ranges from 0 to 1. Then, the error series e t t = 1 n are calculated by Equation (13):
e t = 2 ( r t 0.5 ) e max , t = 1 , 2 , n
where e max is the maximum possible error proportional to the mean absolute error of the study variable.

3.2. Single-Factor Information Extraction under the Influence of Noise Intensity

The most concerning problem in the demonstration of the information extraction method is whether the extraction method is effective and how its effect is affected by the error. So, in this paper, we take the unit hydrograph confluence model as an example to analyze the error inversion estimation of a flood runoff. The unit hydrographs are shown in Figure 1. Suppose there are two runoff generation periods with 10 mm and 12 mm runoff depths in the synthetic case. One hundred thousand sets of error series generated by Equation (13) are added to the process of runoff and flow in each period. The calculation results are shown in Table 1. In the table, V1 and V2 are the mean square error of the estimated runoff. δ1 and δ2 are the relative errors of the estimated runoff yield. α represents the proportional coefficient between the maximum possible error and the average discharge, as shown in Equation (14).
e max = α Q ¯
where Q ¯ is the average flow. We can see from Equation (14) that α also reflects the ratio of the absolute mean value of the error series to the average discharge.
According to the results in Table 1, the mean square error and the relative error of the runoff estimation are both at a small level, less than 0.048 and 3.8%, respectively, although the error proportion coefficient relative to the average flow is 0.2, which is quite a large scale for the average flow error. This indicates that the inversion estimation method of system differential response has a good effect when there is a single factor error.

3.3. Multi-Factor Information Extraction Affected by Noise Intensity

Multi-factor information extraction is more complex than that of a single factor, and the biggest problem is whether it can separate extraction when there are errors from multiple factors. For the convenience of the graphic display of the extraction effect, this paper takes the unit hydrograph routing model as an example to analyze and demonstrate the mechanism. Then, the general effects of the method under multi-factor errors are verified through the final performance of improving runoff simulations.

3.3.1. Feature Analysis of Distinguishing and Extracting Information for Two-Factor Inversion Estimation

Figure 1 shows the unit hydrographs (i.e., system differential response curve) of the surface and underground runoff components (denoted as “Us” and “Ug”, respectively). Figure 2 shows the total residual information of discharge and the error series (denoted as “DQ” and “e”, respectively). It can be seen from the figure that the residual information with rapid changes in the early stage is positive, while negative information with relatively gradual changes is observed in the late stage. Figure 3 shows the effective information process of surface runoff and underground runoff extracted by the inversion estimation method of system differential response. These graphs show the information features extracted by the system differential response inversion method, that is, the updated surface and underground runoff extracted from the information process with the same shape as the surface and underground runoff differential response curve through the correlation between the surface and underground runoff differential response curve and the residual flow process information.

3.3.2. Demonstration of Multi-Factor Information Extraction Effect

In the watershed hydrological system, rainfall (P), evaporation (E), initial soil water content (W), and runoff yield (R) are the main factors affecting the final flow calculation performance. Although conducting theoretical and visual analyses like single-factor and two-factor analyses is impossible, the effect of system differential response inversion estimation can still be comprehensively verified using a synthetic case. We constructed the synthetic case by applying the Xinanjiang hydrological model in the Shuangjiangxi basin. The Xinanjiang (XAJ) rainfall–runoff model is a conceptual rainfall–runoff model proposed by Zhao [32] and has been extensively used in most humid and semi-humid regions in China. At the heart of its architecture is the proposition that a watershed’s soil moisture storage capacity distribution can be systematically represented by a specific curve. A typical structure of the XAJ model consists of a runoff generation module, a runoff separation module, a flow concentration module, and a routing module.
Rainfall and evaporation observation data from 32 historical floods are used to calculate the synthetic flood process. Error series of rainfall, evaporation, initial soil water, and discharge are randomly generated similarly as in the previous section. The same statistical index of flow is used to compare the effect of error inversion estimation of different factors. The calculation procedure is shown in Figure 4, and the results are in Table 2. In the table, R2 is the coefficient of determination, and RE is the relative effective coefficient. An R2 value of 1 indicates a perfect match of the model output to the observed data, while an R2 value of 0 means that the performance of the model is only as accurate as the mean of the observations. RE represents the degree of reduction of the model calculation error using the inverse estimation method of system differential response compared to the original simulations. P, E, W, and R in brackets after RE are the estimated rainfall, evaporation, initial soil water, and runoff, respectively. The R2 and RE calculation formulas are shown in Equations (15) and (16), respectively.
R 2 = 1 t = 1 n ( Q t Q C t ) 2 t = 1 n ( Q t Q ¯ ) 2
R E = 1 t = 1 n ( Q t Q C U i ) 2 t = 1 n ( Q t Q C i ) 2
where Q, QC, and QCU represent the observed flow, model calculated flow, and estimated flow, respectively.
In Table 2, findings derived from 32 distinct flood events are detailed. The results of R2 underscore the significant reduction in forecasting errors achieved by implementing distinct inversion estimates for each flood event. The average REs of 32 floods are 0.832, 0.726, 0.695, 0.904, and 0.835 for single-factor P, e, W, R, and multi-factor PEW, respectively. The most effective among them is the inversion estimation of runoff, followed by the simultaneous inversion estimation of P, E, and W, followed by the single-factor estimation of P, E, and W. The reason is that most of the error factors, including the errors from rainfall, evaporation, initial soil moisture, and the structure of the runoff model, are embedded in runoff errors and thus are considered indirectly when we estimate the runoff errors. However, the simultaneous estimation of P, E, and W only considers the errors of P, E, and W without considering those of the model structure. The single-factor estimation considers a single error factor, resulting in a lower RE. The results show that the more error factors there are involved in the inversion estimation, the more effective information contained in the flow hydrograph can be extracted by the system differential response, which makes the inversion estimation better.

4. Conclusions

The system differential response method for inverse estimation has received much attention in the hydrology literature, while its underlying mechanisms remain largely unexplored. To this end, this paper analyzes the method’s formation mechanism of selective extraction from a theoretical perspective and proposes a new concept (i.e., relation degree coefficient) to analyze the correlation between nonlinear functional variables. Then, we analyze the mechanism of selective extraction from multiple perspectives and construct synthetic cases to demonstrate the effect of the selective extraction of inversion estimation information of the system differential response. The major conclusions of this paper are as follows:
(1)
The inversion estimation method of system differential response can selectively extract effective information, and its mechanism is based on the correlation between the system differential response curve of the factors to be estimated and the information contained in the flow hydrograph.
(2)
The more error factors considered in the inversion estimation, the more effective information contained in the flow hydrograph can be extracted by the system differential response, which makes the inversion estimation better.
In this study, synthetic cases were utilized to verify the effectiveness of the proposed method. Therefore, an intuitive idea for future research is to employ real-world cases to further substantiate and augment our findings. The method should be more broadly tested by expanding case studies across different hydrological models and study areas. For example, one can utilize the method to enhance the forecast performance of the HyMOD hydrological model using the CAMELS dataset [33,34]. Additionally, this study analyzed the mechanism of information extraction through the correction of model inputs and model variables. Therefore, future work is needed to include more factors, such as the model parameters and model structure. While the proposed method has already proved its usefulness in the hydrologic community, it also holds potential as a robust estimation tool in other domains, such as target tracking and assessing environmental sustainability [35,36]. For instance, the proposed method can be used to estimate the states and parameters of a maneuvering target by augmenting the state vector with model parameters.

Author Contributions

Y.C.: Conceptualization, Methodology, Software, Validation, Formal analysis, Data curation, Writing—original draft, Writing—review and editing, Visualization, Funding acquisition. K.L.: Writing—review and editing, Supervision. S.J.: Writing—review and editing, Visualization. Y.S.: Visualization. H.C.: Editing. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

The datasets generated and/or analyzed during the current study are available upon request to the corresponding author.

Conflicts of Interest

The authors declare that they have no known competing financial interest or personal relationships that could have appeared to influence the work reported in this paper.

References

  1. Huang, Y.X.; Wang, Q.Z.; Liang, Z.M.; Deng, X.D. Research advances on real-time correction methods for flood forecasting. South-North Water Transf. Water Sci. Tecnol. 2021, 19, 12–35. [Google Scholar] [CrossRef]
  2. Liang, Z.M.; Huang, Y.X.; Hu, Y.M.; Li, B.Q.; Wang, J. The entire-process correction approach for flood forecasting. South-North Water Transf. Water Sci. Tecnol. 2020, 18, 1–10+17. [Google Scholar] [CrossRef]
  3. Zhang, Y.; Liang, Z.M.; Chen, Z.N.; Zhu, Y.; Wang, J. Study on Flood Forecast and Real-time Correction in Upper Reaches of Dadu River Basin. Water Power 2020, 46, 13–15+21. [Google Scholar]
  4. Bao, W.M.; Zhang, X.Q.; Zhao, L.P. Parameter estimation method based on parameter function surface. Sci. China Technol. Sci. 2013, 56, 1485–1498. [Google Scholar] [CrossRef]
  5. Bao, W.M.; Li, Q.; Qu, S.M. Efficient Calibration Technique under Irregular Response Surface. J. Hydrol. Eng. 2013, 18, 1140–1147. [Google Scholar] [CrossRef]
  6. Li, Q.; Bao, W.; Liang, G.; Qian, J.; Chen, X. Effective Optimization Technique for a Nonlinear Rainfall-Runoff Model. J. Hydrol. Eng. 2014, 19, 1312–1319. [Google Scholar] [CrossRef]
  7. Bao, W.M.; Zhao, L.P. Application of Linearized Calibration Method for Vertically Mixed Runoff Model Parameters. J. Hydrol. Eng. 2014, 19, 04014007. [Google Scholar] [CrossRef]
  8. Baojian, L.I.; Cheng, C.T. Monthly discharge forecasting using wavelet neural networks with extreme learning machine. Sci. China Technol. Sci. 2014, 57, 2441–2452. [Google Scholar]
  9. Yan, R.; Huang, J.; Wang, Y.; Gao, J.; Qi, L. Modeling the combined impact of future climate and land use changes on streamflow of Xinjiang Basin, China. Nord. Hydrol. 2016, 47, 356–372. [Google Scholar]
  10. Liu, X.; Xu, X.; Yu, M.; Lu, J. Hydrological Drought Forecasting and Assessment Based on the Standardized Stream Index in the Southwest China. Procedia Eng. 2016, 154, 733–737. [Google Scholar] [CrossRef]
  11. Patsalides, M.; Efthymiou, V.; Anastassiou, C.; Georghiou, G.E. Transient PV system models for power quality studies. In Proceedings of the Qatar Foundation Annual Research Conference, ARC16, Doha, Qatar, 22–23 March 2016. [Google Scholar]
  12. Si, W.; Bao, W.; Gupta, H.V. Updating real-time flood forecasts via the dynamic system response curve method. Water Resour. Res. 2015, 51, 5128–5144. [Google Scholar] [CrossRef]
  13. Weimin, B.; Wei, S.; Simin, Q. Flow Updating in Real-Time Flood Forecasting Based on Runoff Correction by a Dynamic System Response Curve. J. Hydrol. Eng. 2014, 19, 747–756. [Google Scholar] [CrossRef]
  14. Sun, Y.; Bao, W.; Jiang, P.; Wang, X.; He, C.; Zhang, Q.; Wang, J. Development of dynamic system response curve method for estimating initial conditions of conceptual hydrological models. J. Hydroinform. 2018, 20, 1387–1400. [Google Scholar] [CrossRef]
  15. Sun, Y.; Bao, W.; Jiang, P.; Ji, X.; Gao, S.; Xu, Y.; Zhang, Q.; Si, W. Development of multivariable dynamic system response curve method for real-time flood forecasting correction. Water Resour. Res. 2018, 54, 4730–4749. [Google Scholar] [CrossRef]
  16. Weimin, B.; Yuwei, G.; Wei, S.; Lu, H.; Jinli, L.; Quanfu, L. Application of rainfall dynamic system response curve method for streamflow and sediment simulation in loess region. J. Lake Siences 2020, 2, 528–538. [Google Scholar] [CrossRef]
  17. Si, W.; Zhong, H.; Jiang, P.; Bao, W.; Shi, P.; Qu, S. A dynamic information extraction method for areal mean rainfall error and its application in basins of different scales for flood forecasting. Stoch. Environ. Res. Risk Assess. 2021, 35, 255–270. [Google Scholar] [CrossRef]
  18. Weimin, B.; Wei, S.; Guohua, S.; Xiaoqin, Z.; Qian, L. Runoff error updating based on unit hydrograph inversion. Adv. Water Sci. 2012, 23, 317–322. [Google Scholar]
  19. Li, Q.; Li, C.; Yu, H.; Qian, J.; Hu, L.; Ge, H. System response curve correction method of runoff error for real-time flood forecast. Hydrol. Res. 2020, 51, 1312–1331. [Google Scholar] [CrossRef]
  20. Liu, K.X.; Zhang, X.Q.; Bao, W.M.; Zhao, L.; Li, J. A system response correction method with runoff error smooth matrix. J. Hydraul. Eng. 2015, 46, 960–966. [Google Scholar] [CrossRef]
  21. Zhang, X.Q.; Liu, K.X.; Bao, W.M.; Li, J.J.; Lai, S.Z. Runoff error proportionality coefficient correction method based on system response. Adv. Water Sci. 2014, 25, 789–796. [Google Scholar] [CrossRef]
  22. Zhang, X.Q.; Bao, W.M.; Sun, Y.Q. Enhancing the hydrologic system differential response method for flood forecasting correction. J. Hydrol. 2021, 592, 125793. [Google Scholar] [CrossRef]
  23. Zhang, X.; Bao, W.; Yuan, F. Spatial runoff updating based on the hydrologic system differential response for flood forecasting. J. Hydroinform. 2020, 22, 1573–1587. [Google Scholar] [CrossRef]
  24. Sun, Y.; Bao, W.; Jiang, P.; Si, W.; Zhou, J.; Zhang, Q. Development of a regularized dynamic system response curve for real-time flood forecasting correction. Water 2018, 10, 450. [Google Scholar] [CrossRef]
  25. Lv, N.; Liang, X.; Chen, C.; Zhou, Y.; Li, J.; Wei, H.; Wang, H. A long Short-Term memory cyclic model with mutual information for hydrology forecasting: A Case study in the xixian basin. Adv. Water Resour. 2020, 141, 103622. [Google Scholar] [CrossRef]
  26. Yang, R.; Hou, B.; Xiao, W.; Liang, C.; Zhang, X.; Li, B.; Yu, H. The applicability of real-time flood forecasting correction techniques coupled with the Muskingum method. Hydrol. Res. 2020, 51, 17–29. [Google Scholar] [CrossRef]
  27. Farfán, J.F.; Palacios, K.; Ulloa, J.; Avilés, A. A hybrid neural network-based technique to improve the flow forecasting of physical and data-driven models: Methodology and case studies in Andean watersheds. J. Hydrol. Reg. Stud. 2020, 27, 100652. [Google Scholar] [CrossRef]
  28. Zhang, X.; Liu, P.; Cheng, L.; Liu, Z.; Zhao, Y. A back-fitting algorithm to improve real-time flood forecasting. J. Hydrol. 2018, 562, 140–150. [Google Scholar] [CrossRef]
  29. Ocio, D.; Le Vine, N.; Westerberg, I.; Pappenberger, F.; Buytaert, W. The role of rating curve uncertainty in real-time flood forecasting. Water Resour. Res. 2017, 53, 4197–4213. [Google Scholar] [CrossRef]
  30. Liu, Z.; Guo, S.; Zhang, H.; Liu, D.; Yang, G. Comparative study of three updating procedures for real-time flood forecasting. Water Resour. Manag. 2016, 30, 2111–2126. [Google Scholar] [CrossRef]
  31. Hou, L.; Bao, W.; Si, W.; Jiang, P.; Shi, P.; Qu, S.; Ye, F. Sediment Yield Error Correction by Dynamic System Response Curve Method in Real-Time Flood Forecasting. J. Hydroinform. 2020, 203, 28–46. [Google Scholar] [CrossRef]
  32. Zhao, R.J. The Xinanjiang model applied in China. J. Hydrol. 1992, 135, 371–381. [Google Scholar]
  33. Moradkhani, H.; Hsu, K.; Gupta, H.; Sorooshian, S. Uncertainty assessment of hydrologic model states and parameters: Sequential data assimilation using the particle filter. Water Resour. Res. 2005, 41, W05012. [Google Scholar] [CrossRef]
  34. Addor, N.; Newman, A.J.; Mizukami, N.; Clark, M.P. The CAMELS data set: Catchment attributes and meteorology for large-sample studies. Hydrology and Earth System Sciences 2017, 21, 5293–5313. [Google Scholar] [CrossRef]
  35. Spreafico, C.; Landi, D.; Russo, D. A new method of patent analysis to support prospective life cycle assessment of eco-design solutions. Sustain. Prod. Consum. 2023, 38, 241–251. [Google Scholar] [CrossRef]
  36. Farrell, W.J. Interacting multiple model filter for tactical ballistic missile tracking. IEEE Trans. Aerosp. Electron. Syst. 2008, 44, 418–426. [Google Scholar] [CrossRef]
Figure 1. The unit curves of flow routing. Us is the unit hydrograph of surface runoff, and Ug is the unit hydrograph of underground runoff.
Figure 1. The unit curves of flow routing. Us is the unit hydrograph of surface runoff, and Ug is the unit hydrograph of underground runoff.
Water 15 04016 g001
Figure 2. Model-calculated flow deviations and error series. DQ represents the total residual information process of the flood discharge, and e represents the error series.
Figure 2. Model-calculated flow deviations and error series. DQ represents the total residual information process of the flood discharge, and e represents the error series.
Water 15 04016 g002
Figure 3. The process line of extracted effective information. IRS and IRG represent the effective information-extracting process of surface runoff and underground runoff extracted by the system differential response inversion estimation method.
Figure 3. The process line of extracted effective information. IRS and IRG represent the effective information-extracting process of surface runoff and underground runoff extracted by the system differential response inversion estimation method.
Water 15 04016 g003
Figure 4. Calculation flow chart.
Figure 4. Calculation flow chart.
Water 15 04016 g004
Table 1. Results of the estimated runoff obtained by the system differential response method. V1 and V2 are the mean square errors of estimated runoff. δ1 and δ2 are the relative errors of the estimated runoff yield. α represents the proportional coefficient between the maximum possible error and the average discharge.
Table 1. Results of the estimated runoff obtained by the system differential response method. V1 and V2 are the mean square errors of estimated runoff. δ1 and δ2 are the relative errors of the estimated runoff yield. α represents the proportional coefficient between the maximum possible error and the average discharge.
αV1V2δ1δ2
00000
0.010.0020.0020.0020.001
0.020.0040.0040.0040.003
0.030.0070.0060.0060.005
0.040.0100.0080.0080.006
0.050.0130.0100.0100.008
0.060.0150.0130.0120.010
0.070.0170.0140.0140.011
0.080.0200.0170.0160.013
0.090.0230.0190.0180.015
0.100.0240.0200.0200.016
0.110.0270.0230.0220.018
0.120.0310.0260.0250.021
0.130.0320.0270.0260.021
0.140.0360.0300.0280.024
0.150.0390.0330.0310.026
0.160.0400.0340.0320.027
0.170.0420.0350.0330.028
0.180.0450.0380.0360.030
0.190.0490.0410.0390.032
0.200.0480.0400.0380.032
Table 2. Multi-factor inversion estimation performance. Flood code represents a unique identifier assigned to each analyzed flood event. R2 represents the coefficient of determination indicating the reliability of the model fit. RE(P), RE(E), and RE(W) represent the relative error attributed to discrepancies in precipitation estimation, evaporation estimation, and initial soil water content estimation, respectively. RE(PEW) represents the combined relative error for factors P, E, and W, showing their compounded impact.
Table 2. Multi-factor inversion estimation performance. Flood code represents a unique identifier assigned to each analyzed flood event. R2 represents the coefficient of determination indicating the reliability of the model fit. RE(P), RE(E), and RE(W) represent the relative error attributed to discrepancies in precipitation estimation, evaporation estimation, and initial soil water content estimation, respectively. RE(PEW) represents the combined relative error for factors P, E, and W, showing their compounded impact.
Flood CodeR2RE(P)RE(E)RE(W)RE(PEW)RE(R)
310605170.9940.4200.0000.2710.5810.736
310706130.9630.9070.8990.9010.9380.968
310709170.9860.7860.5410.5180.8270.891
310710070.9510.9320.8410.8280.9360.96
310806090.9850.9190.8780.8530.9350.963
310806100.9850.6120.6250.6290.6600.710
310908080.9680.9530.8930.9600.9200.967
310908120.9900.640.5540.4600.7900.797
311004100.9790.8520.8110.8540.5020.899
311106030.9500.9310.2880.6270.8500.952
311106110.9700.8520.6840.5480.8920.896
311106170.9730.8480.7370.7290.8720.944
311108060.9440.8270.9390.9370.9290.991
311108250.9700.890.8550.9090.6670.964
311206170.9830.7820.7810.6790.6940.839
311207140.9680.8310.5520.7690.8360.922
311208020.9540.8960.7370.8580.8330.945
311208070.9730.9360.9160.6900.9070.952
311304290.9670.840.7590.7390.9320.931
311306060.9840.7230.640.5570.3680.857
311306260.9740.8970.7980.8280.9200.939
311310050.9710.9280.8990.8300.9300.947
311405250.9870.6290.6630.5970.7780.633
311406200.9620.8770.7900.7610.9100.914
311408180.9860.8730.8530.2420.8950.933
311409180.9760.8540.7450.5750.9080.909
311507040.9810.8420.7750.6520.8760.887
311507100.9700.8180.6530.4920.8800.888
311509280.9750.8240.5170.7880.9250.941
311605270.9660.8750.8500.4010.9550.962
311606280.9390.9360.9160.9210.9450.957
311706110.9620.8860.8430.8280.930.947
Mean0.9720.8320.7260.6950.8350.904
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Chen, Y.; Liu, K.; Jiang, S.; Sun, Y.; Chen, H. Mechanism Analysis and Demonstration of Effective Information Extraction in the System Differential Response Inversion Estimation Method. Water 2023, 15, 4016. https://doi.org/10.3390/w15224016

AMA Style

Chen Y, Liu K, Jiang S, Sun Y, Chen H. Mechanism Analysis and Demonstration of Effective Information Extraction in the System Differential Response Inversion Estimation Method. Water. 2023; 15(22):4016. https://doi.org/10.3390/w15224016

Chicago/Turabian Style

Chen, Yang, Kexin Liu, Sijun Jiang, Yiqun Sun, and Hui Chen. 2023. "Mechanism Analysis and Demonstration of Effective Information Extraction in the System Differential Response Inversion Estimation Method" Water 15, no. 22: 4016. https://doi.org/10.3390/w15224016

APA Style

Chen, Y., Liu, K., Jiang, S., Sun, Y., & Chen, H. (2023). Mechanism Analysis and Demonstration of Effective Information Extraction in the System Differential Response Inversion Estimation Method. Water, 15(22), 4016. https://doi.org/10.3390/w15224016

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop