Eunseo Choi edited untitled.tex  almost 9 years ago

Commit id: 6a76a5447c6b1fc0aa464de8912db3025a60d9ba

deletions | additions      

       

\section{Calculating induced dynamic pressure changes in the magma chamber}  \subsection{Model Setup}  We %We  numerically estimate pressure changes induced in the magma chamber by a nuclear explosion. We choose as a numerical tool use  PyLith, ``a finite-element code for dynamic and quasistatic simulations of crustal deformation, primarily earthquakes and volcanoes''\cite{Aagaard_2013}. This code's versatility volcanoes''\cite{Aagaard_2013} to compute pressure changes induced  in meshing and , element types and types of analyses. the magma chamber by a nuclear explosion.  The region around Baekdu and the nuclear test site is represented by a 400 $\times$ 400 $\times$ 100 km domain. While the velocity structure is taken from the regional seismological constraints \cite{DUAN_2005}, we assume that the entire domain is a single layer of crust with Vs = 3.5 km/s and Vp = 6.25 km/s. A magma chamber has been detected beneath Baekdu from low shear wave velocity in the depth range of 10-16 km \cite{ZHANG_2002}. Based on this observation, we assume that the magma chamber is a sphere of 3 km radius, of which center is at 13 km depth and 142 km away to the northeast of the center of the domain, where we put the nuclear explosion site. The S and P wave velocities assigned to the magma chamber are 3.0 km/s and 5.7 km/s. To represent anelastic attenuation in the crustal materials, the isotropic linear Maxwell viscoelastic rheology is assumed for the crust. The constant viscosities of the crust and the magma chamber are 10$^{19}$ and 10$^{18}$ Pa$\cdot$s, respectively.