It has long been recognized that the least-squares estimation method of fitting the best straight line to data points having normally distributed errors yields identical results for the slope and intercept of the line as does the method of maximum likelihood estimation. We show that, contrary to previous understanding, these two methods also give identical results for the standard errors in slope and intercept, provided that the least-squares estimation expressions are evaluated at the least-squares-adjusted points rather than at the observed points as has been done traditionally. This unification of standard errors holds when both x and y observations are subject to correlated errors that vary from point to point. All known correct regression solutions in the literature, including various special cases, can be derived from the original York equations. We present a compact set of equations for the slope, intercept, and newly unified standard errors.

Purely numerical methods based on finite-element approximation of the acoustic or elastic wave equation are becoming increasingly popular for the generation of synthetic seismograms. We present formulas for the grid dispersion and stability criteria for some popular finite-element methods (FEM) for wave propagation, namely, classical and spectral FEM. We develop an approach based on a generalized eigenvalue formulation to analyze the dispersive behavior of these FEMs for acoustic or elastic wave propagation that overcomes difficulties caused by irregular node spacing within the element and the use of high-order polynomials, as is the case for spectral FEM. Analysis reveals that for spectral FEM of order four or greater, dispersion is less than 0.2% at four to five nodes per wavelength, and dispersion is not angle dependent. New results can be compared with grid-dispersion results of some classical finite-difference methods (FDM) used for acoustic or elastic wave propagation. Analysis reveals that FDM and classical FEM require a larger sampling ratio than a spectral FEM to obtain results with the same degree of accuracy. The staggered-grid FDM is an efficient scheme, but the dispersion is angle dependent with larger values along the grid axes. On the other hand, spectral FEM of order four or greater is isotropic with small dispersion, making it attractive for simulations with long propagation times.

S U M M A R YRecently, there has been an increased interest in applying the discontinuous Galerkin method (DGM) to wave propagation. In this work, we investigate the applicability of the interior penalty DGM to elastic wave propagation by analysing it's grid dispersion properties, with particular attention to the effect that different basis functions have on the numerical dispersion. We consider different types of basis functions that naturally yield a diagonal mass matrix. This is relevant to seismology because a diagonal mass matrix is tantamount to an explicit and efficient time marching scheme. We find that the Legendre basis functions that are traditionally used in the DGM introduce numerical dispersion and anisotropy. Furthermore, we find that using Lagrange basis functions along with the Gauss nodes has attractive advantages for numerical wave propagation.

S U M M A R YWe investigate the stability of some high-order finite element methods, namely the spectral element method and the interior-penalty discontinuous Galerkin method (IP-DGM), for acoustic or elastic wave propagation that have become increasingly popular in the recent past. We consider the Lax-Wendroff method (LWM) for time stepping and show that it allows for a larger time step than the classical leap-frog finite difference method, with higher-order accuracy. In particular the fourth-order LWM allows for a time step 73 per cent larger than that of the leap-frog method; the computational cost is approximately double per time step, but the larger time step partially compensates for this additional cost. Necessary, but not sufficient, stability conditions are given for the mentioned methods for orders up to 10 in space and time. The stability conditions for IP-DGM are approximately 20 and 60 per cent more restrictive than those for SEM in the acoustic and elastic cases, respectively.

No abstract