This paper revisits the implementation of surrogate modeling (metamodeling) techniques within seismic risk assessment, for applications that the seismic hazard is described through stochastic ground motion models. Emphasis is placed on how to efficiently address the aleatoric uncertainty in the ground motions, stemming in this case from the stochastic sequence utilized within the excitation model. Previous work has accommodated this uncertainty by approximating the statistics of the engineering demand parameters (EDPs), something that required a large number of replication simulations (for different stochastic sequences) for each training point that was used to inform the metamodel calibration. Using kriging (Gaussian Process regression) as surrogate model, an alternative formulation is discussed here, aiming to minimize the replications for each training point. This is achieved by approximating directly the EDP distribution. It is shown that accommodating heteroscedastic behavior with respect to the aleatoric uncertainty is absolutely critical for achieving an accurate approximation, and two different approaches are presented for establishing this objective. The first approach adopts a stochastic kriging formulation, utilizing a small number of replications for judicially selected inputs, leveraging a secondary surrogate model over the latter inputs to address the heteroscedastic behavior. The second approach uses no replications, establishing a heteroscedastic nugget formulation to accommodate the EDP distribution estimation. A functional relationship is introduced between the nugget and the excitation intensity features to approximate the heteroscedastic behavior. This relationship is explicitly optimized during the metamodel calibration.

During landfalling tropical storms, predictions of the expected storm surge are critical for guiding evacuation and emergency response/preparedness decisions, both at regional and national levels. Forecast errors related to storm track, intensity, and size impact these predictions and, thus, should be explicitly accounted for. The Probabilistic tropical storm Surge (P-Surge) model is the established approach from the National Weather Service (NWS) to achieve this objective. Historical forecast errors are utilized to specify probability distribution functions for different storm features, quantifying, ultimately, the uncertainty in the National Hurricane Center advisories. Surge statistics are estimated by using the predictions across a storm ensemble generated by sampling features from the aforementioned probability distribution functions. P-Surge relies, currently, on a full factorial sampling scheme to create this storm ensemble, combining representative values for each of the storm features. This work investigates an alternative formulation that can be viewed as a seamless extension to the current NHC framework, adopting a quasi-Monte Carlo (QMC) sampling implementation with ultimate goal to reduce the computational burden and provide surge predictions with the same degree of statistical reliability, while using a smaller number of sample storms. The definition of forecast errors adopted here directly follows published NWS practices, while different uncertainty levels are considered in the examined case studies, in order to offer a comprehensive validation. This validation, considering different historical storms, clearly demonstrates the advantages QMC can offer.

This paper investigates the impact of the different sources of excitation variability within the stochastic kriging framework recently developed by the authors for the estimation of the distribution of engineering demand parameters (EDPs) in applications that the seismic hazard is described through stochastic ground motion models. For a given seismic event, described by seismicity characteristics such as the moment magnitude and the distance to source, one can distinguish two type of uncertainties in the excitation description: (g.i) the stochastic sequence utilized within the excitation model; (g.ii) the predictive models relating the ground motion parameters, such as significant duration, arias intensity, or frequency properties of seismic waves, to the seismicity characteristics. The original formulation of the authors treats (g.i) as aleatoric uncertainty and (g.ii) as parametric uncertainty, including the latter in the model parameters for the risk characterization, representing the metamodel input. Implementation establishes the EDP distribution approximation by utilizing a database of EDP estimates for a set of different input parameters, considering replications of these estimates for different stochastic sequences for some small part of this database. The database with replications is first leveraged to approximate the heteroscedas-2466

No abstract

No abstract