Computation of free oscillations of the earth

Raymond P. Buland; F. Gilbert

doi:10.1016/0021-9991(84)90141-4

Home/Archive/USGS

USGS · 70010302

Computation of free oscillations of the earth

Raymond P. Buland·F. Gilbert

View original source ↗Export metadata ↓

Abstract

Although free oscillations of the Earth may be computed by many different methods, numerous practical considerations have led us to use a Rayleigh-Ritz formulation with piecewise cubic Hermite spline basis functions. By treating the resulting banded matrix equation as a generalized algebraic eigenvalue problem, we are able to achieve great accuracy and generality and a high degree of automation at a reasonable cost.

Explore related subjects

Journal of Computational Physics Article

Keep this discovery

Explore connections, maps & timelines

BibTeX RIS

Raymond P. Buland, F. Gilbert. 1984. Computation of free oscillations of the earth. https://doi.org/10.1016/0021-9991(84)90141-4

Cite the original work for its findings. Save a collection to share your selection of sources.

KEEP EXPLORING

Related USGS reports

USGS2025

Total uncertainty quantification in inverse solutions with deep learning surrogate models

We propose an approximate Bayesian method for quantifying the total uncertainty in inverse partial differential equation (PDE) solutions obtained with machine learning surrogate models, including operator learning models. The proposed method accounts for uncertainty in the observations, PDE, and surrogate models. First, we use the surrogate model to formulate a minimization problem in the reduced space for the maximum a posteriori (MAP) inverse solution. Then, we randomize the MAP objective function and obtain samples of the posterior distribution by minimizing different realizations of the objective function. We test the proposed framework by comparing it with the iterative ensemble smoother and deep ensembling methods for a nonlinear diffusion equation with an unknown space-dependent diffusion coefficient. Among other applications, this equation describes the flow of groundwater in an unconfined aquifer. Depending on the training dataset and ensemble sizes, the proposed method provides similar or more descriptive posteriors of the parameters and states than the iterative ensemble smoother method. Deep ensembling underestimates uncertainty and provides less-informative posteriors than the other two methods. Our results show that, despite inherent uncertainty, surrogate models can be used for parameter and state estimation as an alternative to the inverse methods relying on (more accurate) numerical PDE solvers.

Journal of Computational Physics↗

USGS2020

Dimensional effects of inter-phase mass transfer on attenuation of structurally trapped gaseous carbon dioxide in shallow aquifers

Based on experimental evidence and using mathematical modeling, inter-phase mass transfer processes of CO 2 exsolving from and dissolving into water in heterogeneous porous media are investigated under two fundamentally different flow conditions: in a quasi one dimensional vertical column and in a two-dimensional tank with a lateral background water flow, both at laboratory scale. In both cases, the CO 2 dissolved in water under a given overpressure is injected for a certain period at the bottom of the tank, exsolves, and migrates upwards. A layer of fine sand is present in the tanks designed to mimic geological scenarios of accumulation and trapping of exsolved CO 2 in shallow aquifers. Then, clean water is injected and the accumulated CO 2 is dissolved back into the flowing water. The study aims to point out the differences in the mass transfer processes between the quasi-1D and 2D cases using a mathematical model of two-phase compositional flow in heterogeneous porous media calibrated to the experimental datasets, and expose strategies that should be explored in future research. Additionally, temperature variations observed during the 2D experiments allow for analysis of isothermal versus non-isothermal effects on the processes of multiphase CO 2 evolution. The mathematical model is discretized and solved using the mixed hybrid finite element method in 2D that allows for the simulation of both advection- and diffusion-dominated processes accurately.

Journal of Computational Physics↗

USGS2008

Ocean forecasting in terrain-following coordinates: Formulation and skill assessment of the Regional Ocean Modeling System

Systematic improvements in algorithmic design of regional ocean circulation models have led to significant enhancement in simulation ability across a wide range of space/time scales and marine system types. As an example, we briefly review the Regional Ocean Modeling System, a member of a general class of three-dimensional, free-surface, terrain-following numerical models. Noteworthy characteristics of the ROMS computational kernel include: consistent temporal averaging of the barotropic mode to guarantee both exact conservation and constancy preservation properties for tracers; redefined barotropic pressure-gradient terms to account for local variations in the density field; vertical interpolation performed using conservative parabolic splines; and higher-order, quasi-monotone advection algorithms. Examples of quantitative skill assessment are shown for a tidally driven estuary, an ice-covered high-latitude sea, a wind- and buoyancy-forced continental shelf, and a mid-latitude ocean basin. The combination of moderate-order spatial approximations, enhanced conservation properties, and quasi-monotone advection produces both more robust and accurate, and less diffusive, solutions than those produced in earlier terrain-following ocean models. Together with advanced methods of data assimilation and novel observing system technologies, these capabilities constitute the necessary ingredients for multi-purpose regional ocean prediction systems.

Journal of Computational Physics↗