deletions | additions       diff --git a/untitled.tex b/untitled.tex  index d04c693..a53ef2e 100644  --- a/untitled.tex  +++ b/untitled.tex 
...
We numerically estimate pressure changes induced in the magma chamber by a nuclear explosion.  We choose as a numerical tool PyLith, ``a finite-element code for dynamic and quasistatic simulations of crustal deformation, primarily earthquakes and volcanoes''\cite{Aagaard_2013}. This code's versatility in meshing and , element types and types of analyses.  The region aroundMt.  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. The stresses are initially zero everywhere in the domain and the gravitational body force is not considered. These conditions are suitable to this study's purpose of estimating only dynamic stress changes. Absorbing boundary conditions are applied to all the boundaries except the top boundary, which is a free surface. 
...