For the determination of lateral velocity or absorption inhomogeneities, methods such as the generalized matrix inversion and its damped versions, for example the stochastic inverse, are usually applied in seismology to travel-time or amplitude anomalies. These methods are not appropriate for the solution of very extensive systems of equations. Reconstruction techniques as developed for computer tomography are suitable for operations with extremely large numbers of equations and unknown parameters. In this paper solutions obtained with the BPT (Back Projection Technique), ART (Algebraic Reconstruction Technique) and SIRT (Simultaneous Iterative Reconstruction Technique) are compared with those obtained from a damped version of the generalized inverse method. Data of 2-D model-seismic experiments are presented for demonstration.
We introduce a new framework of numerical multiscale methods for advection-dominated problems motivated by climate sciences. Current numerical multiscale methods (MsFEM) work well on stationary elliptic problems but have difficulties when the model involves dominant lower order terms. Our idea to overcome the associated difficulties is a semi-Lagrangian based reconstruction of subgrid variability into a multiscale basis by solving many local inverse problems. Globally the method looks like a Eulerian method with multiscale stabilized basis. We show example runs in one and two dimensions and a comparison to standard methods to support our ideas and discuss possible extensions to other types of Galerkin methods, higher dimensions and nonlinear problems.
We implement the ADER-DG numerical method using the CUDA-C language to run the code in a Graphic Processing Unit (GPU). We focus on solving linear hyperbolic partial differential equations where the method can be expressed as a combination of precomputed matrix multiplications becoming a good candidate to be used on the GPU hardware. Moreover, the method is arbitrarily high-order involving intensive work on local data, a property that is also beneficial for the target hardware. We compare our GPU implementation against CPU versions of the same method observing similar convergence properties up to a threshold where the error remains fixed. This behaviour is in agreement with the CPU version but the threshold is larger that in the CPU case. We also observe a big difference when considering single and double precision where in the first case the threshold error is significantly larger. Finally, we did observe a speed up factor in computational time but this is relative to the specific test or benchmark problem
Computing forecasts of hazards, such as tsunamis, requires fast reaction times and high precision, which in turn demands for large computing facilities that are needed only in rare occasions. Cloud computing environments allow to configure largely scalable on-demand computing environments. In this study, we tested two of the major cloud computing environments for parallel scalability for relevant prototypical applications. These applications solve stationary and non-stationary partial differential equations by means of finite differences and finite elements. These test cases demonstrate the capacity of cloud computing environments to provide scalable computing power for typical tasks in geophysical applications. As a proof-of-concept example of an instant computing application for geohazards, we propose a workflow and prototypical implementation for tsunami forecasting in the cloud. We demonstrate that minimal on-site computing resources are necessary for such a forecasting environment. We conclude by outlining the additional steps necessary to implement an operational tsunami forecasting cloud service, considering availability and cost.
It has been attempted to observe the subsynoptic scale turbulent kinetic energy with high resolution upper-air wind profiles. After filtering out the characteristic large-scale profile, the remaining disturbance part of the observed profile is deemed subsynoptic scale in nature. It is then subject to a direct harmonic analysis for the portion of the atmosphere above 5,000 m in height to examine the distribution of variances over harmonic numbers.It is suggested that the major portion of the subsynoptic scale turbulent energy seems to be contained in the variances from the harmonic analysis of the isolated component of the wind profiles. Examination of obtained variance spectra indicates that they can be grouped in a number of characteristic types according to the concentration of kinetic energy in certain harmonic number ranges.It is suggested that this is related to the flow of kinetic energy from large to small scales or vice versa through subsynoptic scale ranges. Those subsynoptic scale components of kinetic energy in the upper and mid troposphere are also examined in relation to the Richardson number.
scite is a Brooklyn-based organization that helps researchers better discover and understand research articles through Smart Citations–citations that display the context of the citation and describe whether the article provides supporting or contrasting evidence. scite is used by students and researchers from around the world and is funded in part by the National Science Foundation and the National Institute on Drug Abuse of the National Institutes of Health.