Hybrid Galerkin numerical modelling of elastodynamics and compressible Navier–Stokes couplings: applications to seismo-gravito acoustic waves

We introduce a hybrid Galerkin modelling tool for the nonlinear acoustic and gravity wave propagation in planetary atmospheres coupled through topography to a solid medium. We rely on a 2-D spectral-element technique to model linear visco-elastic solid media and couple it to a discontinuous Galerkin method for the atmosphere modelled by the fully nonlinear Navier–Stokes equations. Significant benefits of such a method are, first, its versatility because it handles both acoustic and gravity waves in the same simulation, second, it enables one to observe nonlinear effects as convection or wave-breaking and, finally, it allows one to study the impact of ground-atmosphere coupling for waves propagating from seismic sources. Simulations are performed for 2-D isothermal atmosphere models with complex wind and viscosity profiles. We validate the computations by comparing them to finite-difference solutions, already validated in a previous paper. Specific benchmark validation cases are considered for both acoustic and gravity waves subject to viscosity variations, wind duct and nonlinear wave breaking. We apply this tool to study acoustic and gravity waves generated by a strong seismic source and its nonlinear breaking in the upper atmosphere

Validation of a fast semi-analytic method for surface-wave propagation in layered media

Green’s functions provide an efficient way to model surface-wave propagation and estimate physical quantities for near-surface processes. Several surface-wave Green’s function approximations (far-field, no mode conversions and no higher mode surface waves) have been employed for numerous applications such as estimating sediment flux in rivers, determining the properties of landslides, identifying the seismic signature of debris flows or to study seismic noise through cross-correlations. Based on those approximations, simple empirical scalings exist to derive phase velocities and amplitudes for pure power-law velocity structures providing an exact relationship between the velocity model and the Green’s functions. However, no quantitative estimates of the accuracy of these simple scalings have been reported for impulsive sources in complex velocity structures. In this paper, we address this gap by comparing the theoretical predictions to high-order numerical solutions for the vertical component of the wavefield. The Green’s functions computation shows that attenuation-induced dispersion of phase and group velocity plays an important role and should be carefully taken into account to correctly describe how surface-wave amplitudes decay with distance. The comparisons confirm the general reliability of the semi-analytic model for power-law and realistic shear velocity structures to describe fundamental-mode Rayleigh waves in terms of characteristic frequencies, amplitudes and envelopes. At short distances from the source, and for large near-surface velocity gradients or high Q values, the low-frequency energy can be dominated by higher mode surface waves that can be captured by introducing additional higher mode Rayleigh-wave power-law scalings. We also find that the energy spectral density for realistic shear-velocity models close to piecewise power-law models can be accurately modelled using the same non-dimensional scalings. The frequency range of validity of each power-law scaling can be derived from the corresponding phase velocities. Finally, highly discontinuous near-surface velocity profiles can also be approximated by a combination of power-law scalings. Analytical Green’s functions derived from the non-dimensionalization provide a good estimate of the amplitude and variations of the energy distribution, although the predictions are quite poor around the frequency bounds of each power-law scaling

La planta hotelera en la ciudad de Cajamarca su oferta crecimiento anual entre los años 2012-2016

Desde tiempos remotos junto al desarrollo y crecimiento de las sociedades, la industria hotelera ha ido evolucionando a nivel mundial, la planta hotelera en la ciudad de Cajamarca no ha sido ajena a estos cambios, es por ello que este trabajo de investigación pretende a través de un estudio minucioso crear una base de datos actualizada anualmente respecto a la oferta hotelera existente en Cajamarca, para analizar el crecimiento anual durante los años 2012 y 2016, a fin de definir la magnitud y composición de la planta hotelera. Es así que se hizo un recorrido anualmente de una manera exhaustiva por todos los sectores de la ciudad de Cajamarca con la finalidad de recolectar dicha información. Una vez finalizado el recorrido, la recolección de datos fue analizada y organizada obteniendo el siguiente resultado: el crecimiento de la planta hotelera en la ciudad de Cajamarca entre los años 2012 y 2016 aumentó en 20%, en cuanto a su oferta la clase hostales tuvo mayor crecimiento en los últimos años respecto a las demás clases y la categoría dos estrellas de la clase hostales también tuvo un ligero aumento. Su oferta se ha expandido a lo largo de la ciudad en la mayoría de sus sectores, siendo los sectores 6 “Chontapaccha” y el sector 13 “San Martín” los que tuvieron mayor crecimiento, el sector 8 “La Merced” ubicado en el centro histórico contiene la mayor cantidad de establecimientos de hospedaje con y sin categoría respecto a los demás sectores, en la ciudad para el año 2016 existen 297 establecimientos en funcionamiento los cuales están clasificados en: hotel, hostal, hospedaje y casa hospedaje. Adicionalmente, se elaboró el mapa hotelero de la ciudad de Cajamarca 2016, que contiene todos los establecimientos de hospedaje inventariados durante el tiempo de duración de la investigación.Tesi

Broadband transmission losses and time dispersion maps from time-domain numerical simulations in ocean acoustics

In this letter, a procedure for the calculation of transmission loss maps from numerical simulations in the time domain is presented. It can be generalized to arbitrary time sequences and to elastic media and provides an insight into how energy spreads into a complex configuration. In addition, time dispersion maps can be generated. These maps provide additional information on how energy is distributed over time. Transmission loss and time dispersion maps are generated at a negligible additional computational cost. To illustrate the type of transmission loss maps that can be produced by the time-domain method, the problem of the classical two-dimensional upslope wedge with a fluid bottom is addressed. The results obtained are compared to those obtained previously based on a parabolic equation. Then, for the same configuration, maps for an elastic bottom and maps for non-monochromatic signals are computed

Finite-Difference Modeling of Acoustic and Gravity Wave Propagation in Mars Atmosphere: Application to Infrasounds Emitted by Meteor Impacts

The propagation of acoustic and gravity waves in planetary atmospheres is strongly dependent on both wind conditions and attenuation properties. This study presents a finite-difference modeling tool tailored for acoustic-gravity wave applications that takes into account the effect of background winds, attenuation phenomena (including relaxation effects specific to carbon dioxide atmospheres) and wave amplification by exponential density decrease with height. The simulation tool is implemented in 2D Cartesian coordinates and first validated by comparison with analytical solutions for benchmark problems. It is then applied to surface explosions simulating meteor impacts on Mars in various Martian atmospheric conditions inferred from global climate models. The acoustic wave travel times are validated by comparison with 2D ray tracing in a windy atmosphere. Our simulations predict that acoustic waves generated by impacts can refract back to the surface on wind ducts at high altitude. In addition, due to the strong nighttime near-surface temperature gradient on Mars, the acoustic waves are trapped in a waveguide close to the surface, which allows a night-side detection of impacts at large distances in Mars plains. Such theoretical predictions are directly applicable to future measurements by the INSIGHT NASA Discovery mission

Numerical Simulation of the Atmospheric Signature of Artificial and Natural Seismic Events

The mechanical coupling between solid planets and their atmospheres enables seismically induced acoustic waves to propagate in the atmosphere. We numerically simulate this coupled system for two application cases: active seismic experiments (ASEs) and passive seismic experiments. A recent ASE (Krishnamoorthy et al., 2018, https://doi.org/10.1002/2018GL077481) observed the infrasonic signals produced by a seismic hammer. To measure the sensitivity of observations to seismic parameters, we attempt to reproduce the results from this experiment at short range by considering a realistic unconsolidated subsurface and an idealized rock-solid subsurface. At long range, we investigate the influence of the source by using two focal mechanisms. We found surface waves generate an infrasonic plane head wave in the ASE case of the rock-solid material. For the passive seismic experiments, the amplitude of atmospheric infrasound generated by seismic surface waves is investigated in detail. Despite some limitations, the simulations suggest that balloon measurement of seismically induced infrasound might help to constrain ground properties

Atmospheric Science with InSight

International audienceIn November 2018, for the first time a dedicated geophysical station, the InSight lander, will be deployed on the surface of Mars. Along with the two main geophysical packages, the Seismic Experiment for Interior Structure (SEIS) and the Heat-Flow and Physical Properties Package (HP3), the InSight lander holds a highly sensitive pressure sensor (PS) and the Temperature and Winds for InSight (TWINS) instrument, both of which (along with the InSight FluxGate (IFG) Magnetometer) form the Auxiliary Sensor Payload Suite (APSS). Associated with the RADiometer (RAD) instrument which will measure the surface brightness temperature, and the Instrument Deployment Camera (IDC) which will be used to quantify atmospheric opacity, this will make InSight capable to act as a meteorological station at the surface of Mars. While probing the internal structure of Mars is the primary scientific goal of the mission, atmospheric science remains a key science objective for InSight. InSight has the potential to provide a more continuous and higher-frequency record of pressure, air temperature and winds at the surface of Mars than previous in situ missions. In the paper, key results from multiscale meteorological modeling, from Global Climate Models to Large-Eddy Simulations, are described as a reference for future studies based on the InSight measurements during operations. We summarize the capabilities of InSight for atmospheric observations, from profiling during Entry, Descent and Landing to surface measurements (pressure, temperature, winds, angular momentum), and the plans for how InSight’s sensors will be used during operations, as well as possible synergies with orbital observations. In a dedicated section, we describe the seismic impact of atmospheric phenomena (from the point of view of both “noise” to be decorrelated from the seismic signal and “signal” to provide information on atmospheric processes). We discuss in this framework Planetary Boundary Layer turbulence, with a focus on convective vortices and dust devils, gravity waves (with idealized modeling), and large-scale circulations. Our paper also presents possible new, exploratory, studies with the InSight instrumentation: surface layer scaling and exploration of the Monin-Obukhov model, aeolian surface changes and saltation / lifing studies, and monitoring of secular pressure changes. The InSight mission will be instrumental in broadening the knowledge of the Martian atmosphere, with a unique set of measurements from the surface of Mars

Numerical modeling of atmospheric waves due to Earth/Ocean/Atmosphere couplings and applications

Cette thèse se penche sur la propagation d’ondes au sein du système coupléTerre-océan-atmosphère. La compréhension de ces phénomènes a une importance majeure pour l’étude de perturbations sismiques et d’explosions atmosphériques notamment dans le cadre de missions spatiales planétaires. Les formes d’ondes issues du couplage fluide-solide permettent d’obtenir de précieuses informations sur la source du signal ou les propriétés des milieux de propagation. On développe donc deux outils numériques d’ordre élevé pour la propagation d’ondes acoustiques et de gravité. L'u en différences finies et se concentre sur le milieu atmosphérique et la propagation d’ondes linéaires dans un milieu stratifié visqueux et avec du vent. Cette méthode linéaire est validée par des solutions quasi-analytiques reposant sur les équations de dispersion dans une atmosphère stratifiée. Elle est aussi appliquée à deux cas d’études : la propagation d’ondes liée à l’impact d’une météorite à la surface de Mars (mission NASA INSIGHT), et la propagation d’ondes atmosphériques liées au tsunami de Sumatra en 2004. La seconde méthode résout la propagation non-linéaire d’ondes gravito-acoustiques dans une atmosphère couplée, avec topographie, à la propagation d’ondes élastiques dans un solide visco-élastique. Cette méthode repose sur sur le couplage d’une formulation en éléments finis discontinus, pour résoudre les équations de Navier-Stokes la partie fluide, par éléments finis continus pour résoudre les équations de l’élastodynamique dans la partie solide. Elle a été validée grâce à des solutions analytiques ainsi que par des comparaisons avec les résultats de la méthode par différences finies.This thesis deals with the wave propagation problem within the Earth-ocean-atmosphere coupled system. A good understanding of the these phenomena has a major importance for seismic and atmospheric explosion studies, especially for planetary missions. Atmospheric wave-forms generated by explosions or surface oscillations can bring valuable information about the source mechanism or the properties of the various propagation media. We develop two new numerical full-wave high-order modeling tools to model the propagation of acoustic and gravity waves in realistic atmospheres. The first one relies on a high-order staggered finite difference method and focus only on the atmosphere. It enables the simultaneous propagation of linear acoustic and gravity waves in stratified viscous and windy atmosphere. This method is validated against quasi-analytical solutions based on the dispersion equations for a stratified atmosphere. It has also been employed to investigate two cases : the atmospheric propagation generated by a meteor impact on Mars for the INSIGHT NASA mission and for the study of tsunami-induced acoutic and gravity waves following the 2004 Sumatra tsunami. The second numerical method resolves the non-linear acoustic and gravity wave propagation in a realistic atmosphere coupled, with topography, to the elastic wave propagation in a visco-elastic solid. This numerical tool relies on a discontinuous Galerkin method to solve the full Navier-Stokes equations in the fluid domain and a continuous Galerkin method to solve the elastodynamics equations in the solid domain. It is validated against analytical solutions and numerical results provided by the finite-difference method