论文信息 - Numerical Efficiency in Monte Carlo Simulations—Case Study of a River Thermodynamic Model

Numerical Efficiency in Monte Carlo Simulations—Case Study of a River Thermodynamic Model

Trade-offs between precision of numerical solutions to deterministic models of the environment, and the number of model realizations achievable within a framework of Monte Carlo simulation, are investigated and discussed. A case study of a model of river thermodynamics is employed. It is shown that the tractability of Monte Carlo simulation relies on adaptation of the numerical solution time-step, giving results with a guaranteed error in the time domain as well as near-optimum speed of calibration under any chosen accuracy criteria. Time-step control is implemented using two adaptive Runge-Kutta methods: a second order scheme with first order error estimator, and an embedded fourth-fifth order scheme. In the case study, where the effects of sparse and imprecise data dominate the overall modeling error, both the schemes appear adequate. However, the higher order scheme is concluded to be generally more reliable and efficient, and has wide potential to improve the value of applying the Monte Carlo method to environmental simulation. The problem of reconciling spatial error with the specified temporal error is discussed.

[1] Keith Beven,et al. Bayesian estimation of uncertainty in land surface‐atmosphere flux predictions , 1997 .

[2] Hung Tao Shen,et al. A Mathematical Model for River Ice Processes , 1991 .

[3] R. Bras. Hydrology : an introduction to hydrologic science , 1990 .

[4] Kjell Gustafsson,et al. Control of Error and Convergence in ODE Solvers , 1992 .

[5] Peter Reichert,et al. On the usefulness of overparameterized ecological models , 1997 .

[6] Hung Tao Shen,et al. Simulation of Growth and Decay of River Ice Cover , 1984 .

[7] M. B. Carver. The Choice of Algorithms in Automated Method of Lines Solution of Partial Differential Equations , 1976 .

[8] Alan H. Karp,et al. A variable order Runge-Kutta method for initial value problems with rapidly varying right-hand sides , 1990, TOMS.

[9] Heinz G. Stefan,et al. Riverbed heat conduction prediction , 1994 .

[10] Howard S. Wheater,et al. Estimation and propagation of parametric uncertainty in environmental models , 2002 .

[11] M. D. McKay,et al. A comparison of three methods for selecting values of input variables in the analysis of output from a computer code , 2000 .

[12] E. Hairer,et al. Solving ordinary differential equations I (2nd revised. ed.): nonstiff problems , 1993 .

[13] Keith Beven,et al. The future of distributed models: model calibration and uncertainty prediction. , 1992 .

[14] S. Chapra. Surface Water-Quality Modeling , 1996 .