Simultaneous perturbation stochastic approximation for tidal models
KAUST DepartmentPhysical Sciences and Engineering (PSE) Division
Earth Science and Engineering Program
Permanent link to this recordhttp://hdl.handle.net/10754/334621
MetadataShow full item record
AbstractThe Dutch continental shelf model (DCSM) is a shallow sea model of entire continental shelf which is used operationally in the Netherlands to forecast the storm surges in the North Sea. The forecasts are necessary to support the decision of the timely closure of the moveable storm surge barriers to protect the land. In this study, an automated model calibration method, simultaneous perturbation stochastic approximation (SPSA) is implemented for tidal calibration of the DCSM. The method uses objective function evaluations to obtain the gradient approximations. The gradient approximation for the central difference method uses only two objective function evaluation independent of the number of parameters being optimized. The calibration parameter in this study is the model bathymetry. A number of calibration experiments is performed. The effectiveness of the algorithm is evaluated in terms of the accuracy of the final results as well as the computational costs required to produce these results. In doing so, comparison is made with a traditional steepest descent method and also with a newly developed proper orthogonal decompositionbased calibration method. The main findings are: (1) The SPSA method gives comparable results to steepest descent method with little computational cost. (2) The SPSA method with little computational cost can be used to estimate large number of parameters.
CitationAltaf MU, Heemink AW, Verlaan M, Hoteit I (2011) Simultaneous perturbation stochastic approximation for tidal models. Ocean Dynamics 61: 1093-1105. doi:10.1007/s10236-011-0387-6.
The following license files are associated with this item:
Except where otherwise noted, this item's license is described as This is an Open Access article distributed under the terms of the Creative Commons-Attribution Noncommercial License, which permits unrestricted noncommercial use, distribution, and reproduction in any medium, provided the original author(s) and source are credited.
Showing items related by title, author, creator and subject.
Control of modeling error in calibration and validation processes for predictive stochastic modelsOden, J. Tinsley; Prudhomme, Serge (International Journal for Numerical Methods in Engineering, Wiley-Blackwell, 2010-09-28) [Article]The idea of adaptive control of modeling error is expanded to include ideas of statistical calibration, validation, and uncertainty quantification. © 2010 John Wiley & Sons, Ltd.
Clustered iterative stochastic ensemble method for multi-modal calibration of subsurface flow modelsElsheikh, Ahmed H.; Wheeler, Mary Fanett; Hoteit, Ibrahim (Journal of Hydrology, Elsevier, 2013-05) [Article]A novel multi-modal parameter estimation algorithm is introduced. Parameter estimation is an ill-posed inverse problem that might admit many different solutions. This is attributed to the limited amount of measured data used to constrain the inverse problem. The proposed multi-modal model calibration algorithm uses an iterative stochastic ensemble method (ISEM) for parameter estimation. ISEM employs an ensemble of directional derivatives within a Gauss-Newton iteration for nonlinear parameter estimation. ISEM is augmented with a clustering step based on k-means algorithm to form sub-ensembles. These sub-ensembles are used to explore different parts of the search space. Clusters are updated at regular intervals of the algorithm to allow merging of close clusters approaching the same local minima. Numerical testing demonstrates the potential of the proposed algorithm in dealing with multi-modal nonlinear parameter estimation for subsurface flow models. © 2013 Elsevier B.V.
Boosting iterative stochastic ensemble method for nonlinear calibration of subsurface flow modelsElsheikh, Ahmed H.; Tavakoli, Reza; Wheeler, Mary Fanett; Hoteit, Ibrahim (Computer Methods in Applied Mechanics and Engineering, Elsevier BV, 2013-06) [Article]A novel parameter estimation algorithm is proposed. The inverse problem is formulated as a sequential data integration problem in which Gaussian process regression (GPR) is used to integrate the prior knowledge (static data). The search space is further parameterized using Karhunen-Loève expansion to build a set of basis functions that spans the search space. Optimal weights of the reduced basis functions are estimated by an iterative stochastic ensemble method (ISEM). ISEM employs directional derivatives within a Gauss-Newton iteration for efficient gradient estimation. The resulting update equation relies on the inverse of the output covariance matrix which is rank deficient.In the proposed algorithm we use an iterative regularization based on the ℓ2 Boosting algorithm. ℓ2 Boosting iteratively fits the residual and the amount of regularization is controlled by the number of iterations. A termination criteria based on Akaike information criterion (AIC) is utilized. This regularization method is very attractive in terms of performance and simplicity of implementation. The proposed algorithm combining ISEM and ℓ2 Boosting is evaluated on several nonlinear subsurface flow parameter estimation problems. The efficiency of the proposed algorithm is demonstrated by the small size of utilized ensembles and in terms of error convergence rates. © 2013 Elsevier B.V.