SENSITIVITY OF TSUNAMI HEIGHT DISTRIBUTION TO FAULT MODEL PARAMETERS: CASE STUDY BASED ON THE 2016 OFF FUKUSHIMA EARTHQUAKE

2021 ◽  
Vol 77 (2) ◽  
pp. I_187-I_192
Author(s):  
Yutaka HAYASHI
Keyword(s):  
Fault Model ◽  
Model Parameters ◽  
Height Distribution ◽  
Tsunami Height

scholarly journals Real-time automatic uncertainty estimation of coseismic single rectangular fault model using GNSS data

2021 ◽  
Vol 73 (1) ◽  
Author(s):  
Keitaro Ohno ◽  
Yusaku Ohta ◽  
Satoshi Kawamoto ◽  
Satoshi Abe ◽  
Ryota Hino ◽  
...  
Keyword(s):  
Real Time ◽  
Disaster Response ◽  
Order Approximation ◽  
Nonlinear Problems ◽  
Large Earthquake ◽  
Fault Model ◽  
Deformation Monitoring ◽  
Model Parameters ◽  
The Sea Of Japan ◽  
Plate Interface

AbstractRapid estimation of the coseismic fault model for medium-to-large-sized earthquakes is key for disaster response. To estimate the coseismic fault model for large earthquakes, the Geospatial Information Authority of Japan and Tohoku University have jointly developed a real-time GEONET analysis system for rapid deformation monitoring (REGARD). REGARD can estimate the single rectangular fault model and slip distribution along the assumed plate interface. The single rectangular fault model is useful as a first-order approximation of a medium-to-large earthquake. However, in its estimation, it is difficult to obtain accurate results for model parameters due to the strong effect of initial values. To solve this problem, this study proposes a new method to estimate the coseismic fault model and model uncertainties in real time based on the Bayesian inversion approach using the Markov Chain Monte Carlo (MCMC) method. The MCMC approach is computationally expensive and hyperparameters should be defined in advance via trial and error. The sampling efficiency was improved using a parallel tempering method, and an automatic definition method for hyperparameters was developed for real-time use. The calculation time was within 30 s for 1 × 106 samples using a typical single LINUX server, which can implement real-time analysis, similar to REGARD. The reliability of the developed method was evaluated using data from recent earthquakes (2016 Kumamoto and 2019 Yamagata-Oki earthquakes). Simulations of the earthquakes in the Sea of Japan were also conducted exhaustively. The results showed an advantage over the maximum likelihood approach with a priori information, which has initial value dependence in nonlinear problems. In terms of application to data with a small signal-to-noise ratio, the results suggest the possibility of using several conjugate fault models. There is a tradeoff between the fault area and slip amount, especially for offshore earthquakes, which means that quantification of the uncertainty enables us to evaluate the reliability of the fault model estimation results in real time.

scholarly journals Parameter stability and consistency in an alongshore-current model determined with Markov chain Monte Carlo

2008 ◽  
Vol 10 (2) ◽  
pp. 153-162 ◽  
Author(s):  
B. G. Ruessink
Keyword(s):  
Monte Carlo ◽  
Markov Chain ◽  
Markov Chain Monte Carlo ◽  
Parameter Estimates ◽  
Model Parameters ◽  
Parameter Stability ◽  
Small Uncertainty ◽  
Best Fit ◽  
Stability And Consistency

When a numerical model is to be used as a practical tool, its parameters should preferably be stable and consistent, that is, possess a small uncertainty and be time-invariant. Using data and predictions of alongshore mean currents flowing on a beach as a case study, this paper illustrates how parameter stability and consistency can be assessed using Markov chain Monte Carlo. Within a single calibration run, Markov chain Monte Carlo estimates the parameter posterior probability density function, its mode being the best-fit parameter set. Parameter stability is investigated by stepwise adding new data to a calibration run, while consistency is examined by calibrating the model on different datasets of equal length. The results for the present case study indicate that various tidal cycles with strong (say, >0.5 m/s) currents are required to obtain stable parameter estimates, and that the best-fit model parameters and the underlying posterior distribution are strongly time-varying. This inconsistent parameter behavior may reflect unresolved variability of the processes represented by the parameters, or may represent compensational behavior for temporal violations in specific model assumptions.

Dynamic Data-Driven Spatiotemporal System Behavior Prediction With Simulations and Sensor Measurement Data

2018 ◽  
Author(s):  
Xiangxue Zhao ◽  
Shapour Azarm ◽  
Balakumar Balachandran
Keyword(s):  
Dynamical System ◽  
Real Time ◽  
Measurement Data ◽  
Data Driven ◽  
Sensor Data ◽  
Model Parameters ◽  
High Fidelity ◽  
System Behavior ◽  
Sensor Measurement

Online prediction of dynamical system behavior based on a combination of simulation data and sensor measurement data has numerous applications. Examples include predicting safe flight configurations, forecasting storms and wildfire spread, estimating railway track and pipeline health conditions. In such applications, high-fidelity simulations may be used to accurately predict a system’s dynamical behavior offline (“non-real time”). However, due to the computational expense, these simulations have limited usage for online (“real-time”) prediction of a system’s behavior. To remedy this, one possible approach is to allocate a significant portion of the computational effort to obtain data through offline simulations. The obtained offline data can then be combined with online sensor measurements for online estimation of the system’s behavior with comparable accuracy as the off-line, high-fidelity simulation. The main contribution of this paper is in the construction of a fast data-driven spatiotemporal prediction framework that can be used to estimate general parametric dynamical system behavior. This is achieved through three steps. First, high-order singular value decomposition is applied to map high-dimensional offline simulation datasets into a subspace. Second, Gaussian processes are constructed to approximate model parameters in the subspace. Finally, reduced-order particle filtering is used to assimilate sparsely located sensor data to further improve the prediction. The effectiveness of the proposed approach is demonstrated through a case study. In this case study, aeroelastic response data obtained for an aircraft through simulations is integrated with measurement data obtained from a few sparsely located sensors. Through this case study, the authors show that along with dynamic enhancement of the state estimates, one can also realize a reduction in uncertainty of the estimates.

Hierarchical Bayesian Parameter Estimation for Modeling and Analysis of User Affective Influence

2013 ◽  
Author(s):  
Feng Zhou ◽  
Jianxin (Roger) Jiao
Keyword(s):  
Decision Making ◽  
Interior Design ◽  
Computational Models ◽  
Human Perception ◽  
Model Parameters ◽  
Hierarchical Bayesian ◽  
Affective States ◽  
Modeling And Analysis ◽  
Human Decision

Traditional user experience (UX) models are mostly qualitative in terms of its measurement and structure. This paper proposes a quantitative UX model based on cumulative prospect theory. It takes a decision making perspective between two alternative design profiles. However, affective elements are well-known to have influence on human decision making, the prevailing computational models for analyzing and simulating human perception on UX are mainly cognition-based models. In order to incorporate both affective and cognitive factors in the decision making process, we manipulate the parameters involved in the cumulative prospect model to show the affective influence. Specifically, three different affective states are induced to shape the model parameters. A hierarchical Bayesian model with a technique called Markov chain Monte Carlo is used to estimate the parameters. A case study of aircraft cabin interior design is illustrated to show the proposed methodology.

scholarly journals SHIMMER (1.0): a novel mathematical model for microbial and biogeochemical dynamics in glacier forefield ecosystems

2015 ◽  
Vol 8 (10) ◽  
pp. 3441-3470 ◽  
Author(s):  
J. A. Bradley ◽  
A. M. Anesio ◽  
J. S. Singarayer ◽  
M. R. Heath ◽  
S. Arndt
Keyword(s):  
Organic Matter ◽  
Microbial Biomass ◽  
Microbial Growth ◽  
Model Development ◽  
Biogeochemical Cycling ◽  
Biogeochemical Model ◽  
Model Parameters ◽  
Glacier Forefield ◽  
Unconstrained Model

Abstract. SHIMMER (Soil biogeocHemIcal Model for Microbial Ecosystem Response) is a new numerical modelling framework designed to simulate microbial dynamics and biogeochemical cycling during initial ecosystem development in glacier forefield soils. However, it is also transferable to other extreme ecosystem types (such as desert soils or the surface of glaciers). The rationale for model development arises from decades of empirical observations in glacier forefields, and enables a quantitative and process focussed approach. Here, we provide a detailed description of SHIMMER, test its performance in two case study forefields: the Damma Glacier (Switzerland) and the Athabasca Glacier (Canada) and analyse sensitivity to identify the most sensitive and unconstrained model parameters. Results show that the accumulation of microbial biomass is highly dependent on variation in microbial growth and death rate constants, Q10 values, the active fraction of microbial biomass and the reactivity of organic matter. The model correctly predicts the rapid accumulation of microbial biomass observed during the initial stages of succession in the forefields of both the case study systems. Primary production is responsible for the initial build-up of labile substrate that subsequently supports heterotrophic growth. However, allochthonous contributions of organic matter, and nitrogen fixation, are important in sustaining this productivity. The development and application of SHIMMER also highlights aspects of these systems that require further empirical research: quantifying nutrient budgets and biogeochemical rates, exploring seasonality and microbial growth and cell death. This will lead to increased understanding of how glacier forefields contribute to global biogeochemical cycling and climate under future ice retreat.

scholarly journals Evaluating performances of simplified physically based models for landslide susceptibility

2015 ◽  
Vol 12 (12) ◽  
pp. 13217-13256 ◽  
Author(s):  
G. Formetta ◽  
G. Capparelli ◽  
P. Versace
Keyword(s):  
Landslide Susceptibility ◽  
Goodness Of Fit ◽  
Shallow Landslides ◽  
Model Parameters ◽  
Test Case ◽  
Susceptibility Analysis ◽  
Physically Based ◽  
Private And Public ◽  
Physically Based Models

Abstract. Rainfall induced shallow landslides cause loss of life and significant damages involving private and public properties, transportation system, etc. Prediction of shallow landslides susceptible locations is a complex task that involves many disciplines: hydrology, geotechnical science, geomorphology, and statistics. Usually to accomplish this task two main approaches are used: statistical or physically based model. Reliable models' applications involve: automatic parameters calibration, objective quantification of the quality of susceptibility maps, model sensitivity analysis. This paper presents a methodology to systemically and objectively calibrate, verify and compare different models and different models performances indicators in order to individuate and eventually select the models whose behaviors are more reliable for a certain case study. The procedure was implemented in package of models for landslide susceptibility analysis and integrated in the NewAge-JGrass hydrological model. The package includes three simplified physically based models for landslides susceptibility analysis (M1, M2, and M3) and a component for models verifications. It computes eight goodness of fit indices by comparing pixel-by-pixel model results and measurements data. Moreover, the package integration in NewAge-JGrass allows the use of other components such as geographic information system tools to manage inputs-output processes, and automatic calibration algorithms to estimate model parameters. The system was applied for a case study in Calabria (Italy) along the Salerno-Reggio Calabria highway, between Cosenza and Altilia municipality. The analysis provided that among all the optimized indices and all the three models, the optimization of the index distance to perfect classification in the receiver operating characteristic plane (D2PC) coupled with model M3 is the best modeling solution for our test case.

scholarly journals Introductory Models of COVID-19 in the United States

2021 ◽  
Vol 2 (3) ◽  
pp. 74-98
Author(s):  
Peter Hugo Nelson
Keyword(s):  
United States ◽  
The United States ◽  
Transition Function ◽  
Model Parameters ◽  
Microsoft Excel ◽  
New Model ◽  
Social Distancing ◽  
Unlimited Growth ◽  
Least Squares Techniques

ABSTRACT Students develop and test simple kinetic models of the spread of coronavirus disease 2019 (COVID-19) caused by the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) virus. Microsoft Excel is used as the modeling platform because it is nonthreatening to students and it is widely available. Students develop finite difference models and implement them in the cells of preformatted spreadsheets following a guided inquiry pedagogy that introduces new model parameters in a scaffolded step-by-step manner. That approach allows students to investigate the implications of new model parameters in a systematic way. Students fit the resulting models to reported cases per day data for the United States using least squares techniques with Excel's Solver. Using their own spreadsheets, students discover for themselves that the initial exponential growth of COVID-19 can be explained by a simplified unlimited growth model and by the susceptible-infected-recovered (SIR) model. They also discover that the effects of social distancing can be modeled using a Gaussian transition function for the infection rate coefficient and that the summer surge was caused by prematurely relaxing social distancing and then reimposing stricter social distancing. Students then model the effect of vaccinations and validate the resulting susceptible-infected-recovered-vaccinated (SIRV) model by showing that it successfully predicts the reported cases per day data from Thanksgiving through the holiday period up to 14 February 2021. The same SIRV model is then extended and successfully fits the fourth peak up to 1 June 2021, caused by further relaxation of social distancing measures. Finally, students extend the model up to the present day (27 August 2021) and successfully account for the appearance of the delta variant of the SARS-CoV-2 virus. The fitted model also predicts that the delta variant peak will be comparatively short, and the cases per day data should begin to fall off in early September 2021, counter to current expectations. This case study makes an excellent capstone experience for students interested in scientific modeling.

Modeling and Control of a Complex Interacting Process

2011 ◽  
Vol 403-408 ◽  
pp. 3758-3762
Author(s):  
Subhajit Patra ◽  
Prabirkumar Saha
Keyword(s):  
Storage System ◽  
Linear Quadratic Regulator ◽  
Model Parameters ◽  
Linear Quadratic ◽  
Modeling And Control ◽  
Dynamic Matrix ◽  
Liquid Storage ◽  
Efficient Control ◽  
And Control

In this paper, two efficient control algorithms are discussed viz., Linear Quadratic Regulator (LQR) and Dynamic Matrix Controller (DMC) and their applicability has been demonstrated through case study with a complex interacting process viz., a laboratory based four tank liquid storage system. The process has Two Input Two Output (TITO) structure and is available for experimental study. A mathematical model of the process has been developed using first principles. Model parameters have been estimated through the experimentation results. The performance of the controllers (LQR and DMC) has been compared to that of industrially more accepted PID controller.

Cosmolian: High-resolution Integrated OSL and CN modelling Program for Sand Transportation in Eolian Realms (HIPSTER)

2021 ◽  
Author(s):  
Shlomy Vainer ◽  
Yoav Ben Dor
Keyword(s):  
Landscape Evolution ◽  
Sand Dunes ◽  
Integrative Approach ◽  
Model Parameters ◽  
Arid Environments ◽  
Aeolian Transport ◽  
Migration Patterns ◽  
Aeolian Processes ◽  
Simpson Desert

<p>The extensivity of sand dunes in continental interiors makes the understating of their morphodynamical properties valuable for palaeoenvironmental reconstructions and the interpretation of landscape evolution. Nevertheless, the study of aeolian landscape development at the million-years timescale is hampered by the complex interaction of factors determining dune migration and the inherently self-destructive nature of their chronostratigraphy, thus limiting the applicability of traditional luminescence-based dating methods for configuring processes beyond ~300 Ka. In this study, we present a standalone program that simulates aeolian transport based on luminescence-derived chronologies coupled with numerical modelling of cosmogenic nuclides accumulation. This integrative approach reveals ancient phases of sand irruption and provides a data-based scheme facilitating the morphodynamical study of aeolian processes over multiple timescales. We present a case study of the program application by analyzing data from the Australian Simpson Desert, unfolding several phases of aeolian vitality since the late Pliocene. The synchronicity of the results with drastic changes in environmental settings exemplifies the applicability of process-based modelling in constructing a timeframe of key landscape evolution events in arid environments by studying aeolian landforms. Finally, the relationships between model parameters used to determine environmental settings on sand migration patterns make the program a powerful tool to further investigating triggers and mechanisms of aeolian processes.</p>

scholarly journals Interpreting the Estimates from the Full VECH Model with Asymmetry: The Case of US and Canadian Equity Prices

2017 ◽  
Vol 6 (4) ◽  
pp. 236
Author(s):  
Chikashi Tsuji
Keyword(s):  
Parameter Estimates ◽  
Model Parameters ◽  
Model Structure ◽  
Equity Prices ◽  
Related Research ◽  
Asymmetric Effects ◽  
Equity Index ◽  
Autoregressive Conditional Heteroscedasticity ◽  
Full Vector

This paper attempts to derive careful interpretation of the parameter estimates from one of the multivariate generalized autoregressive conditional heteroscedasticity (MGARCH) models, the full vector-half (VECH) model with asymmetric effects. We also consider and interpret the parameter estimates from a case study of US and Canadian equity index returns by applying this model. More specifically, we firstly inspect the model formula and derive general interpretation of the model parameters. We consider this is particularly useful for understanding not only the full VECH model structure but also similar MGARCH models. After the general considerations, we also interpret the case results that are derived from our application of the full VECH model to US and Canadian equity index returns. We consider that these concrete illustrations are also very helpful for future related research.

Sign in / Sign up

Export Citation Format

Share Document