Search NASA⌕ Search

SEARCH · Search NASA

Results for “Theoretical and computational geophysics”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

30 records · Page 2
All sourcesNASA NTRSDOE OSTI
NASA NTRS1994

Multifractal characterizations of nonstationarity and intermittency in geophysical fields: Observed, retrieved, or simulated

Geophysical data rarely show any smoothness at any scale, and this often makes comparison with theoretical model output difficult. However, highly fluctuating signals and fractal structures are typical of open dissipative systems with nonlinear dynamics, the focus of most geophysical research. High levels of variability are excited over a large range of scales by the combined actions of external forcing and internal instability. At very small scales we expect geophysical fields to be smooth, but these are rarely resolved with available instrumentation or simulation tools; nondifferentiable and even discontinuous models are therefore in order. We need methods of statistically analyzing geophysical data, whether measured in situ, remotely sensed or even generated by a computer model, that are adapted to these characteristics. An important preliminary task is to define statistically stationary features in generally nonstationary signals. We first discuss a simple criterion for stationarity in finite data streams that exhibit power law energy spectra and then, guided by developments in turbulence studies, we advocate the use of two ways of analyzing the scale dependence of statistical information: singular measures and qth order structure functions. In nonstationary situations, the approach based on singular measures seeks power law behavior in integrals over all possible scales of a nonnegative stationary field derived from the data, leading to a characterization of the intermittency in this (gradient-related) field. In contrast, the approach based on structure functions uses the signal itself, seeking power laws for the statistical moments of absolute increments over arbitrarily large scales, leading to a characterization of the prevailing nonstationarity in both quantitative and qualitative terms. We explain graphically, step by step, both multifractal statistics which are largely complementary to each other. The geometrical manifestations of nonstationarity and intermittency, 'roughness' and 'sparseness', respectively, are illustrated and the associated analytical (differentiability and continuity) properties are discussed. As an example, the two techniques are applied to a series of recent measurements of liquid water distributions inside marine stratocumulus decks; these are found to be multifractal over scales ranging from approximately 60 m to approximately 60 km. Finally, we define the 'mean multifractal plane' and show it to be a simple yet comprehensive tool with many applications including data intercomparison, (dynamical or stochastic) model and retrieval validations.

Davis, Anthony↗
NASA NTRS1993

Laboratory and theoretical models of planetary-scale instabilities and waves

Research work is proceeding in theoretical, numerical, and experimental geophysical fluid dynamics leading up to a reflight of the GFFC (Geophysical Fluid Flow Cell Experiment) on USML-2. The work is intended not only to generate ideas for future space experiments, but to provide fundamental results concerned with nonlinear and chaotic properties of thermal convection and baroclinic waves in terrestrial and planetary atmospheres. The major efforts are focussed on thermal convection in a rapidly rotating annulus relevant to Jovian atmospheric dynamics, and on the chaotic behavior of baroclinic waves relevant to the Earth's atmosphere. The approach, in preparation for USML-2, is primarily theoretical and numerical. Mechanistic process models are solved numerically in order to identify physical mechanisms that may be observed in the GFFC, and which are important in real geophysical applications. The results from numerical simulations of geophysical fluid flow (subject to rotation and stratification) are compared with previous GFFC experiments on Spacelab-3 and with existing and proposed terrestrial laboratory experiments of various types. Pattern recognition algorithms have been employed to generate low-dimensional descriptions of the highly nonlinear and turbulent numerical simulations. Such empirically truncated descriptions provide for simplified but robust physical interpretations of the dynamics, as well as yielding highly efficient computations of these chaotic flows.

Hart, John E.↗
NASA NTRS1990

Theoretical and experimental studies of baroclinic processes

The geophysical fluid flow simulation code, GEOSIM, is being used to study the phenomenon of vacillation in the baroclinic annulus. Having verified that the code predicts vacillation for the same points as the experiments, the work is aiming toward explaining the mechanics of vacillation and pointing out some of the sensitivities of the results to the numerical method. Researchers are finding that there is a structural change associated with amplitude vacillation, where the structural changes are in the vertical. The results disagree with the premise of Lindzen et al, that the vacillation is due to constructive and destructive interference of neutral modes with different phase speeds. The researchers are continuing to study the Spacelab 3 Geophysical Fluid Flow Cell (GFFC) results with horizontal temperature gradients and heating from below. GEOSIM has been used to compute a wide range of cases, and these are being compared with the observations. The computations and observations compare well, and the model is being used to extend the results beyond cases studied in the experiments and to study the mechanics and predictability of the flows. The study of fully nonlinear baroclinic instability using the GFFC apparatus is proceeding with the numerical code. While the first instability that occurs is of planetary scale, secondary instabilities consisting of small-scale, penetrative convection occurs where cold fluid flows over a warm surface. The simultaneous modeling of the planetary scale and the convective scale is possible because of the nonhydrostatic formulation of the model. Some of these results have been animated on the Stardent computer, which shows the explosive nature of the small-scale convection.

Miller, Tim↗
NASA NTRS2003

The Earth System Model

The Earth System Model is the natural evolution of current climate models and will be the ultimate embodiment of our geophysical understanding of the planet. These models are constructed from components - atmosphere, ocean, ice, land, chemistry, solid earth, etc. models and merged together through a coupling program which is responsible for the exchange of data from the components. Climate models and future earth system models will have standardized modules, and these standards are now being developed by the ESMF project funded by NASA. The Earth System Model will have a variety of uses beyond climate prediction. The model can be used to build climate data records making it the core of an assimilation system, and it can be used in OSSE experiments to evaluate. The computing and storage requirements for the ESM appear to be daunting. However, the Japanese ES theoretical computing capability is already within 20% of the minimum requirements needed for some 2010 climate model applications. Thus it seems very possible that a focused effort to build an Earth System Model will achieve succcss.

Schoeberl, Mark↗
DOE OSTI2022

Dynamic mode decomposition with core sketch

With the increase in collected data volumes, either from experimental measurements or high fidelity simulations, there is an ever-growing need to develop computationally efficient tools to process, analyze, and interpret these datasets. Modal analysis techniques have gained great interest due to their ability to identify patterns in the data and extract valuable information about the system being considered. Dynamic mode decomposition (DMD) relies on elements of the Koopman approximation theory to compute a set of modes, each associated with a fixed oscillation frequency and a decay/growth rate. Extracting these details from large datasets can be computationally expensive due to the need to implement singular value decomposition of the input data matrix. Sketching algorithms have become popular in numerical linear algebra where statistical theoretic approaches are utilized to reduce the cost of major operations. A sketch of a matrix is another matrix, which is significantly smaller, but still sufficiently approximates the original system. We put forth an efficient DMD framework, SketchyDMD, based on a core sketching algorithm that captures information about the range and corange (their mutual relationship) of input data. The proposed sketching-based framework can accelerate various portions of the DMD routines, compared to classical methods that operate directly on the raw input data. We conduct numerical experiments using the spherical shallow water equations as a prototypical model in the context of geophysical flows. In conclusion, we show that the proposed SketchyDMD is superior to existing randomized DMD methods that are based on capturing only the range of the input data.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

A Fast Solver for the Fractional Helmholtz Equation

The purpose of this paper is to study a Helmholtz problem with a spectral fractional Laplacian, instead of the standard Laplacian. Recently, it has been established that such a fractional Helmholtz problem better captures the underlying behavior in Geophysical Electromagnetics. In this work, we establish the well-posedness and regularity of this problem. We introduce a hybrid spectral-finite element approach to discretize it and show well-posedness of the discrete system. In addition, we derive a priori discretization error estimates. Finally, we introduce an efficient solver that scales aswell as the best possible solver for the classical integer-order Helmholtz equation. We conclude withseveral illustrative examples that confirm our theoretical findings.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Using Deep Learning to Derive Shear-Wave Velocity Models from Surface-Wave Dispersion Data

We present a new algorithm for derivations of 1D shear-wave velocity models from surface-wave dispersion data using convolutional neural networks (CNNs). The technique is applied for continental China and the plate boundary region in southern California. Different CNNs are designed for these two regions and are trained using theoretical Rayleigh-wave phase and group velocity images computed from reference 1D VS models. The methodology is tested with 3260 phase–group images for continental China and 4160 phase–group images for southern California. The conversions of these images to velocity profiles take ~23 s for continental China and ~30 s for southern California on a personal laptop with the NVIDIA GeForce GTX 1060 core and a memory of 6 GB. The results obtained by the CNNs show high correlation with previous studies using conventional methods. The effectiveness of the CNN technique makes this fast method an important alternative for deriving shear-wave velocity models from large datasets of surface-wave dispersion data.

Geochemistry & Geophysics↗
NASA NTRS1985

Theoretical and experimental studies in support of the geophysical fluid flow experiment

Meteorologists and astrophysicists interested in large scale planetary and solar circulations have come to recognize the importance of rotation and stratification in determining the character of these flows. In particular, the effect of latitude-dependent Coriolis force on nonlinear convection is thought to play a crucial role in such phenomena as differential rotation on the Sun, cloud band orientation on Jupiter, and the generation of magnetic fields in thermally driven dynamos. The continuous low-gravity environment of the orbiting space shuttle offers a unique opportunity to make laboratory studies of such large-scale thermally driven flows under the constraint imposed by rotation and sphericity. This is possible because polarization forces in a dielectric liquid, which are linearly dependent on fluid temperature, give rise to an effectively radial buoyancy force when a radial electrostatic field is imposed. The Geophysical Fluid Flow Cell (GFFC) is an implementation of this ideal in which fluid is contained between two rotating hemispheres that are differentially heated and stressed with a large a-c voltage. The experiment, to be flown on Spacelab III (currently set for launch April 29, 1985), will explore non-linear mode selection and high Rayleigh number turbulence in a rotating convecting spherical shell of liquid. Experiments will be carried out in a low driving parameter range where some limited numerical experimentation is currently feasible, as well as in a parameter range significantly beyond numerical computation for many years.

Hart, J.↗
DOE OSTI2021

Multichannel Analysis of Surface Waves Accelerated (MASWAccelerated): Software for efficient surface wave inversion using MPI and GPUs

Multichannel Analysis of Surface Waves (MASW) is a technique frequently used in geotechnical engineering and engineering geophysics to infer 1D layered models of seismic shear wave velocities in the top tens to hundreds of meters of the subsurface. We aim to accelerate MASW calculations by capitalizing on modern computer hardware available in the workstations of most engineers: multiple cores and graphics processing units (GPUs). We propose new parallel and GPU accelerated algorithms for computing 1D MASW inversion, and provide software implementations in C using Message Passing Interface (MPI) and CUDA. These algorithms take advantage of sparsity that arises in the problem, and the work balance between processes considers typical data trends. We compare our methods to an existing open source Matlab MASW tool. Our serial C implementation achieves a 2x speedup over the Matlab software, and we continue to see improvements by parallelizing the problem with MPI. Here we see nearly perfect strong and weak scaling for uniform data, and improve strong scaling for realistic data by repartitioning the problem to process mapping. By utilizing GPUs available on most modern workstations, we observe an additional 1.3x speedup over the serial C implementation on the first use of the method. We typically repeatedly evaluate theoretical dispersion curves as part of an optimization procedure, and on the GPU the kernel can be cached for faster reuse on later runs. We observe a 3.2x speedup on the cached GPU runs compared to the serial C runs. This work is the first open-source parallel or GPU-accelerated software tool for MASW imaging, and should enable geotechnical engineers to fully utilize all computer hardware at their disposal.

58 GEOSCIENCES↗
DOE OSTI2025

Investigating the Relationship Between Bolide Entry Angle and Apparent Direction of Infrasound Signal Arrivals

Infrasound sensing offers critical capabilities for detecting and geolocating bolide events globally. However, the observed back azimuths, directions from which infrasound signals arrive at stations, often differ from the theoretical expectations based on the bolide’s peak brightness location. For objects with shallow entry angles, which traverse longer atmospheric paths, acoustic energy may be emitted from multiple points along the trajectory, leading to substantial variability in back azimuth residuals. This study investigates how the entry angle of energetic bolides affects the back azimuth deviations, independent of extrinsic factors such as atmospheric propagation, station noise, and signal processing methodologies. A theoretical framework, the Bolide Infrasound Back-Azimuth EXplorer Model (BIBEX-M), was developed to compute predicted back azimuths solely from geometric considerations. The model quantifies how these residuals vary as a function of source-to-receiver distance, revealing that bolides entering at shallow angles, e.g., 10°, can produce average residuals of 20°, with deviations reaching up to 46° at distances below 1000 km, and remaining significant even at 5000 km (up to 8°). In contrast, bolides with steeper entry angles, e.g., > 60°, show smaller deviations, typically under 5° at 1000 km and diminishing to less than ~1° beyond 5000 km. These findings attest to the need for careful interpretation when evaluating signal detections and estimating bolide locations. This work is not only pertinent to bolides but also to other high-energy, extended-duration atmospheric phenomena such as space debris and reentry events, where similar geometric considerations can influence infrasound arrival directions.

Acoustics↗
DOE OSTI2022

Combining different 3-D global and regional seismic wave propagation solvers towards box tomography in the deep Earth

SUMMARY In previous publications, we presented a general framework, which we called ‘box tomography’, that allows the coupling of any two different numerical seismic wave propagation solvers, respectively outside and inside a target region, or ‘box’. The goal of such hybrid wavefield computations is to reduce the cost of computations in the context of full-waveform inversion for structure within the target region, when sources and/or receivers are located at large distances from the box. Previously, we had demonstrated this approach with sources and receivers outside the target region in a 2-D acoustic spherical earth model, and demonstrated and applied this methodology in the 3-D spherical elastic Earth in a continental scale inversion in which all stations were inside the target region. Here we extend the implementation of the approach to the case of a 3-D global elastic earth model in the case where both sources and stations are outside the box. We couple a global 3-D solver, SPECFEM3D_GLOBE, for the computation of the wavefield and Green’s functions in a reference 3-D model, with a regional 3-D solver, RegSEM, for the computation of the wavefield within the box, by means of time-reversal mirrors. We briefly review key theoretical aspects, showing in particular how only the displacement is needed to be stored at the boundary of the box. We provide details of the practical implementation, including the geometrical design of the mirrors, how we deal with different sizes of meshes in the two solvers, and how we address memory-saving through the use of B-spline compression of the recorded wavefield on the mirror. The proposed approach is numerically efficient but also versatile, since adapting it to other solvers is straightforward and does not require any changes in the solver codes themselves, as long as the displacement can be recovered at any point in time and space. We present benchmarks of the hybrid computations against direct computations of the wavefield between a source and an array of stations in a realistic geometry centred in the Yellowstone region, with and without a hypothetical plume within the ‘box’, and with a 1-D or a 3-D background model, down to a period of 20 s. The ultimate goal of this development is for applications in the context of imaging of remote target regions in the deep mantle, such as, for example, Ultra Low Velocity Zones.

Geochemistry & Geophysics↗
DOE OSTI2020

Integration of Soft Data Into Geostatistical Simulation of Categorical Variables

Uncertain or indirect “soft” data, such as geologic interpretation, driller’s logs, geophysical logs or imaging, offer potential constraints or “soft conditioning” to stochastic models of discrete categorical subsurface variables in hydrogeology such as hydrofacies. Previous bivariate geostatistical simulation algorithms have not fully addressed the impact of data uncertainty in formulation of the (co) kriging equations and the objective function in simulated annealing (or quenching). This paper introduces the geostatistical simulation code tsim-s, which accounts for categorical data uncertainty through a data “hardness” parameter. In generating geostatistical realizations with tsim-s, the uncertainty inherent to soft conditioning is factored into both 1) the data declustering and spatial correlation functions in cokriging and 2) the acceptance probability for change of category in simulated quenching. The degree or sensitivity to which soft data conditions a realization as a function of hardness can be quantified by mapping category probabilities derived from multiple realizations. In addition to point or borehole data, arrays of data (e.g., as derived from a depth-dependency function, probability map, or “prior realization”) can be used as soft conditioning. The tsim-s algorithm provides a theoretically sound and general framework for integrating datasets of variable location, resolution, and uncertainty into geostatistical simulation of categorical variables. A practical example shows how tsim-s is capable of generating a large-scale three-dimensional simulation including curvilinear features.

54 ENVIRONMENTAL SCIENCES↗