- NSF-PAR ID:
- Date Published:
- Journal Name:
- Geophysical Journal International
- Page Range / eLocation ID:
- 586 to 626
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
The analyses of interior penalty discontinuous Galerkin methods of any order k for solving elliptic and parabolic problems with Dirac line sources are presented. For the steady state case, we prove convergence of the method by deriving a priori error estimates in the L 2 norm and in weighted energy norms. In addition, we prove almost optimal local error estimates in the energy norm for any approximation order. Further, almost optimal local error estimates in the L 2 norm are obtained for the case of piecewise linear approximations whereas suboptimal error bounds in the L 2 norm are shown for any polynomial degree. For the time-dependent case, convergence of semi-discrete and of backward Euler fully discrete scheme is established by proving error estimates in L 2 in time and in space. Numerical results for the elliptic problem are added to support the theoretical results.more » « less
We introduce a new finite-element (FE) based computational framework to solve forward and inverse elastic deformation problems for earthquake faulting via the adjoint method. Based on two advanced computational libraries, FEniCS and hIPPYlib for the forward and inverse problems, respectively, this framework is flexible, transparent and easily extensible. We represent a fault discontinuity through a mixed FE elasticity formulation, which approximates the stress with higher order accuracy and exposes the prescribed slip explicitly in the variational form without using conventional split node and decomposition discrete approaches. This also allows the first order optimality condition, that is the vanishing of the gradient, to be expressed in continuous form, which leads to consistent discretizations of all field variables, including the slip. We show comparisons with the standard, pure displacement formulation and a model containing an in-plane mode II crack, whose slip is prescribed via the split node technique. We demonstrate the potential of this new computational framework by performing a linear coseismic slip inversion through adjoint-based optimization methods, without requiring computation of elastic Green’s functions. Specifically, we consider a penalized least squares formulation, which in a Bayesian setting—under the assumption of Gaussian noise and prior—reflects the negative log of the posterior distribution. The comparison of the inversion results with a standard, linear inverse theory approach based on Okada’s solutions shows analogous results. Preliminary uncertainties are estimated via eigenvalue analysis of the Hessian of the penalized least squares objective function. Our implementation is fully open-source and Jupyter notebooks to reproduce our results are provided. The extension to a fully Bayesian framework for detailed uncertainty quantification and non-linear inversions, including for heterogeneous media earthquake problems, will be analysed in a forthcoming paper.
We generated dense, high‐resolution 3‐D ground displacement maps for the 2016 MW7.8 Kaikōura, New Zealand earthquake—the most geometrically and kinematically complex rupture yet recorded—from stereo WorldView optical satellite imagery using a new methodology that combines subpixel image correlation with a ray‐tracing approach. Our analysis reveals fundamental new details of near‐field displacement patterns, which cannot easily be obtained through other methods. From our detailed correlation maps, we measured fault slip in 3‐D along 19 faults at 500‐m spacing. Minimum resolvable horizontal slip is ~0.1 m, and vertical is ~0.5 m. Net slip measurements range from <1 to ~12 m. System‐level kinematic analysis shows that slip on faults north of the Hope fault was oriented primarily subparallel to the Pacific‐Australian plate motion direction. In contrast, slip on faults to the south was primarily at high angle to the plate motion and secondarily parallel to plate motion. Fault kinematics are in some locations consistent with long‐term uplift patterns, but inconsistent in others. Deformation within the Seaward Kaikōura Range may indicate an attempt by the plate boundary fault system to geometrically simplify. Comparison of published field measurements along the Kekerengu fault with our correlation‐derived measurements reveals that ~36% of surface displacement was accommodated as distributed off‐fault deformation when considering only field measurements of discrete slip. Comparatively, field measurements that project previously linear features (e.g., fence lines) into the fault over apertures >5–100 m capture nearly all (~90%) of the surface deformation.
An H(div)-conforming finite element method for the Biot’s consolidation mo- del is developed, with displacements and fluid velocity approximated by elements from BDM_k space. The use of H(div)-conforming elements for flow variables ensures the local mass conservation. In the H(div)-conforming approximation of displacement, the tan- gential components are discretised in the interior penalty discontinuous Galerkin frame- work,and the normal components across the element interfaces are continuous. Having introduced a spatial discretisation, we develop a semi-discrete scheme and a fully dis- crete scheme,prove their unique solvability and establish optimal error estimates for each variable.more » « less
Abstract An interior penalty discontinuous Galerkin method is devised to approximate minimizers of a linear folding model by discontinuous isoparametric finite element functions that account for an approximation of a folding arc. The numerical analysis of the discrete model includes an a priori error estimate in case of an accurate representation of the folding curve by the isoparametric mesh. Additional estimates show that geometric consistency errors may be controlled separately if the folding arc is approximated by piecewise polynomial curves. Various numerical experiments are carried out to validate the a priori error estimate for the folding model.more » « less