Next Article in Journal
Evaluation of HRCLDAS and ERA5 Datasets for Near-Surface Wind over Hainan Island and South China Sea
Next Article in Special Issue
Laser Beam Atmospheric Propagation Modelling for Aerospace LIDAR Applications
Previous Article in Journal
A Cluster Analysis of Forward Trajectory to Identify the Transport Pathway of Salt-Dust Particles from Dried Bottom of Aral Sea, Central Asia
Previous Article in Special Issue
Reduction of the VLF Signal Phase Noise Before Earthquakes
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Numerical Modeling of Coseismic Tropospheric Disturbances Arising from the Unstable Acoustic Gravity Wave Energetics

by
Esfhan A. Kherani
*,†,
Saul A. Sanchez
and
Eurico R. de Paula
National Institute for Space Research, Av. dos Astronautas, 1758, Sao Jose dos Campos 12227010, Brazil
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Atmosphere 2021, 12(6), 765; https://doi.org/10.3390/atmos12060765
Submission received: 5 March 2021 / Revised: 21 April 2021 / Accepted: 22 April 2021 / Published: 13 June 2021

Abstract

:
Numerous recent studies report the Coseismic Tropospheric Disturbances (CTDs) during large earthquakes. Their presence suggests the importance of atmospheric seismology in a possible earthquake forecasting scenario. The origin mechanism and associated energetics of CTDs are not well understood though the observations associate them with the atmospheric waves. We present the numerical modeling of coupled dynamics of Gravity waves (GWs) and convective instability (CI) in the dry troposphere that produces the CTDs, in the form of pressure disturbances, of observed magnitudes. The study reveals the altitude and epicentral distribution of CTDs and elaborates the relative role of GWs and CI in the generation and intensification of CTDs. The study finds that mega and strong earthquakes disturb the troposphere to a similar level as the severe meteorological weather.

1. Introduction

Troposphere hosts varieties of disturbances during near-surface meteorological, orographical, and lithospheric or seismic activities [1]. Among these activities, the former two forcings are omnipresent, the later forcing limits around the epicenter and depends on the magnitude of the earthquakes. Several studies in recent years report the detections of Coseismic Tropospheric Disturbances (CTDs) during strong earthquakes of magnitudes larger than 7 [2,3]. These CTDs are observed as tropospheric delay in the radio signal, as registered from the GNSS receivers scattered across the epicenter [4,5]. Their presence and associated delay in the radio signal influences the accuracy of the navigation and communication system. They are observed concurrently with the tropospheric pressure disturbances that are distinguished from the nonseismic components based on their association to the epicenter [6].
The origin mechanism and associated energetics of coseismic atmospheric disturbances at the ground are well understood. In the case of atmospheric disturbances of meteorological and orographical origins, Acoustic Gravity waves have received increasing attention, primarily due to their role in the generation of small to mesoscale structures. Mikumo (1968) [7] invoked energetics of Acoustic Gravity Waves (AGWs) to explain the detection of coseismic atmospheric disturbances at the ground during the great Alaskan earthquake of magnitude 9. He showed that most of the observed features of ground pressure disturbances are the manifestation of the AGWs of seismic origin. Based on infrasonic monitoring, numerous studies report the detection of coseismic pressure disturbances that may have their origin at various atmospheric heights, notably at stratospheric and thermospheric reflecting layers [8]. The presence of coseismic pressure disturbances at the long epicentral distance of about 6000 km is owing to the robust energetics of AGWs in the reflecting layers [2].
On the other hand, the detection of CTDs and their association to the energetics of AGWs are relatively recent topics in atmospheric seismology [6]. On the theoretical ground, the origin mechanism and origin altitudes of CTDs are the topics that remain unattended and unresolved to date. Recent observations suggest that the CTDs are the result of Lithosphere-Atmosphere coupling [4,5] that eventually leads to the Coseismic atmospheric Disturbances in the middle and upper atmosphere [9]. The coupling is established and energized by the AGWs in response to the ground uplift during the earthquakes. Owing to the decrease in density (or increase in buoyancy) with the height, the amplitude of AGWs from tiny ground uplift amplifies to many order at the thermospheric heights and forms the detectable atmospheric disturbances [9,10]. In the context of CTDs, the role of robust energetics of the AGWs is not examined. In particular, in the troposphere, where buoyancy increases rapidly with height owing to the steepest density gradient, the amplification may attain a factor of about 10–100 and therefore, CTDs of significant magnitudes are expected to be excited. Moreover, the convectively unstable troposphere is likely to intensify the CTDs. The present theoretical study examines the formation of CTDs from the lithosphere-atmosphere coupling energized by the gravity waves and convective instability. Moreover, the study examines the role of ambient conditions in the intensification and spatial distribution of CTDs.

2. Materials and Methods

We employ a nonhydrostatic and nonlinear simulation model of AGWs [10,11]. The model solves the following governing equations of wave amplitude (u), atmospheric density ( ρ ) and pressure (p):
2 u t 2 = 1 ρ ( γ p . u ) p ρ . u p ρ 2 u . ρ + 1 ρ ( u . ) p + Π ν + Π nl
ρ t + . ( ρ u ) = 0
p t + ( u . ) p + γ p . u = 0
Π ν = t ν 2 u + ( η + ν 3 ) ( . u ) ; Π nl = t ( u . u )
Here, Π ν and Π n l are the time derivatives of the viscous and nonlinear forces, ν = μ / ρ and η = η / ρ are the first and second kinematic viscosities, respectively, ( μ , η ) are the first and second dynamic viscosities respectively and γ is the specific heat ratio.
Much of the insight of AGWs can be obtained from a simplified form of the wave Equation (1), as presented in Appendix A.1 and Appendix A.2. The derivation considers strong vertical ( y ) gradient condition that prevails in the atmosphere and assumes nonlinear-nonlocal plane-wave solution of the following form:
u x ( x , y , t ) = u x ( y , t ) e x p ( i k x x + i k y y ) , u y ( x , y , t ) = u y ( y , t ) e x p ( i k x x + i k y y )
The derivation also employs the method of separation of variables, i.e.,
u y = u ys ( y ) u yt ( t ) , u x = u xs ( y ) u xt ( t )
We derive a set of simplified equations for the vertical amplitude of AGWs in the nonviscous atmosphere of the following form:
d 2 u yt dt 2 = Ω 0 2 u yt k x k y c 2 u xs u ys u xt , du ys dy = k 0 u ys
where,
Ω 0 2 = k y 2 c 2 + Ω b 2 , Ω b 2 = ( γ 1 ) k 0 2 k o c 2 d c 2 dy c 2 ζ ( γ a d + γ e )
γ e = 1 c 2 d c 2 d y , γ a d = ( γ 1 ) k 0
k 0 = ζ c 2 , ζ = 1 ρ d p d y , c 2 = γ p ρ
Here, Ω b is the nonisothermal nonhydrostatic Brunt-Vaisala frequency (Equation (6.7a) of [12]), γ e and γ a d are the lapse rate and adiabatic lapse rate respectively, γ is the ratio of the specific heats and c is the sound speed. The conditions Ω b 2 > 0 and Ω b 2 < 0 i.e., | γ e | < | γ a d | and | γ e | > | γ a d | , correspond to the stable and unstable conditions respectively. The condition Ω b 2 > 0 excites pure AGWs and the condition Ω b 2 < 0 excites CI. In the troposphere, both conditions may prevail at different altitudes. Therefore, the troposphere may host both AGWs and CI. The present study excludes the compressional dynamics and considers the energetics of gravity waves (GWs) and convective instability. In this case, the pressure disturbance is due to the air motion from advection and convection.
The numerical modeling solves Equations (1)–(3) in the Cartesian coordinate system where x , y and z respectively represents longitude (+ve towards west), altitude and latitude (+ve towards north). The simulation domain covers 0–600 km in altitude, x e p ± 10 Δ x e p km in longitude and z e p ± Δ z e p km in latitude with identical grid resolution Δ = 1 km. Here, x e p and z e p are the epicentral coordinates and Δ x e p and Δ z e p determines the size of the fault plane. The simulation begins 30 min before the earthquake onset time (t ep ) and span for 1 h with the the time resolution of 10 s.
Numerous studies have investigated the role of gravity waves in the formation of tropospheric disturbances of nonseismic origins. Three kinds of forcings are responsible for wave energetics: diabatic forcing from condensation/evaporation (e.g., Ref. [13]), obstacle effect associated with the wind shear layer within the troposphere that acts like a mountain and presents effective orographical forcing (e.g., Ref. [14]), and the mechanical oscillator mechanism (e.g., Refs. [15,16]). The present work opts for the mechanical oscillator mechanism in which ground vibration from seismic waves couples mechanically to the atmosphere without the loss of momentum. The continuity of vertical ground velocity across the Earth’s surface establishes the coupling. At the lower boundary i.e., at y = 0 km, the continuity of the vertical velocity of the ground associated with the earthquake acts as forcing for the generation of the AGWs, i.e.,
u y ( x , y = 0 , z , t ) = V SISM ( x , z , t )
Figure 1 demonstrates the V SISM that corresponds to a typical earthquake of magnitude ∼8 that attains maximum value of V SISM 1 cm/s, lasts for 10 min and has a fault plane size of about 40 km.
The numerical modeling derives atmosphere at t = 0 from the NRLMSISE model [17] for the epicenter location (31.06 N,103.37 E) of the Wenchuan earthquake of magnitude M w = 8.0 that occurred at 06:28:01.42 UTC on 12 May 2008, in Sichuan province, West China. Jin et al. (2011) [4] report the detection of CTDs from this earthquake event. Figure 2 demonstrates the ambient atmospheric density and temperature profiles. Figure 3 reveals that two temperature profiles represent convectively stable ( | γ e | < | γ a d | ) and unstable ( | γ e | > | γ a d | ) atmospheric conditions. We note that for the convectively unstable temperature profile, the lapse rate is greater than the the adiabatic lapse rate in the altitude range of 2.5–7.5 km.

3. Results

Based on stable and unstable temperature profiles, we carry out the following numerical experiments:
  • Stable+GWs: Simulation of GWs with the stable temperature profile,
  • Unstable+GWs: Simulation of GWs with the unstable temperature profile without the inclusion of CI,
  • Unstable+GWs+CI: Simulation of GWs with the unstable temperature profile and the inclusion of CI
Figure 4 demonstrates the frequency spectrum of the amplitude ( u y ) of GWs at an altitude range of 0–20 km. The black circles represent the Brunt-Vaisala period estimated from expression (5) for the stable+GWs case. The spectrum at the ground represents the seismic source spectrum that shows the principal spectral peak at about 2 min and secondary peaks at periods higher than 10 min. We note that the power spectrum at altitudes above ground also reveals spectral peaks at around 2 min and numerous spectral peaks at other periodicities higher than 2 min. Therefore, the GWs sustain the frequency components of the seismic source, as well as generate other frequency components with periods longer than the Brunt-Vaisala period in the troposphere.
The amplitudes u x and u y cause the disturbances in the atmospheric pressure and density, governed by Equations (2) and (3). We define the CTDs of the following form:
C T D δ p ( x , y , t ) = p ( x , y , t ) p ( x , y , t = 0 )
Figure 5 demonstrates the results from the stable+GWs and Unstable+GWs numerical experiments. In Figure 5A, we note the formation of CTDs in the stable troposphere within 20 min from the earthquake onset time. It attains a maximum of about 5 Pascals in the lower troposphere below the altitude of 10 km. At first, long period-long wavefronts with periods longer than 5 min and wavelength of about 10 km are launched in 0–10 km altitudes during the first 25 min from the earthquake onset. Afterward, short period-short wavefronts with periods of about 2 min and wavelengths lower than about 5 km are launched at altitudes higher than 4 km. Therefore, the phase velocities of these wavefronts are less than 40 m/s, a typical range of phase velocity of GWs in the troposphere [16]. The case demonstrates the energetics of GWs to gives rise to the CTDs. Figure 5B reveals the formation of CTDs in the unstable troposphere, similar to the stable+GWs case with few noticeable exceptions. We note the launch of long period-short wavefronts of periods longer than 5 min and wavelength of about 5 km in 0–10 km altitudes during 0–25 min, instead of long period-long wavefronts from the stable+GWs case. The formation of short wavefronts affirms the generation of small scale GWs in the altitude region of the steep unstable gradient. We also note the formation of descending layers in the lower troposphere, a known energetic of GWs. These layers are intense and long-lived in contrast to the weak and short-lived layers in stable+GWs case. In other words, the GWs create coherent descending layers in an unstable troposphere in comparison to incoherent layers in the stable troposphere.
Figure 6 demonstrates the results from the Unstable+GWs+CI numerical experiment. It reveals the formation of descending layers of CTDs, similar to the previous experiments. However, we note the intensification of layers in the altitude range of 2.5–7.5 km where the lapse rate is larger than the adiabatic rate. Therefore, the intensification is owing to the CI that has a positive growth in the unstable altitude range. The condition Ω b 2 < 0 represents the positive growth condition and the growth rate of CI equals to the imaginary part of Ω b . In the case of unstable+GWs without CI, these layers are less resolved and short-lived. Therefore, the unstable+GWs+CI experiment produces the CTDs that persist for a long duration after the mainshock onset time.
Figure 7 demonstrates the evolution of CTDs (pixmaps) and growth rate of CI (contours) in space at a few selected times from the unstable+GWs+CI experiment. The growth rate is positive in the altitude range of 2.5–7.5 km where the ambient lapse rate becomes steeper than the adiabatic lapse rate (Figure 3) The evolution of the growth rate is due to the varying ambient conditions during the production of CTDs. Thus, the evolution reflects the nonlinear dynamics of GWs and CI. We note the development of CTDs at off-epicenter such that they occupy broad epicentral distance, more than the limited ±20 km size of the fault plane. CTDs attain strong amplitudes above the epicenter and in the altitude range of 2.5–7.5 km where the growth rate of CI is positive. The broad epicentral distance coverage of CTDs is owing to the oblique propagating GWs that have their origin from two sources: the epicentral inhomogeneity in the seismic ground forcing, V SISM and secondary generation of GWs from the epicentral pressure gradients developed from the primary GWs.

4. Discussion

Figure 8 shows the pressure disturbances at the ground and mean pressure disturbance from the lower troposphere. It reveals the ground pressure disturbance of ∼1 Pa from the ground velocity V SISM 0.01 m/s. Moreover, the waveform of the ground disturbance is similar to the waveform of V SISM . Watada et al. (2006) [3] reported the ground pressure disturbances of about 1–5 Pa magnitudes in the vicinity of the epicenter of the Tokachi-Oki earthquake of magnitude 8.3 that had corresponding V SISM 0.005 m/s. In addition, he found similar waveforms of ground forcing and ground pressure disturbations for the first 20 min from the earthquake onset. Thus, the observed characteristics of ground pressure disturbances reported by Watada et al. (2006) [3] and simulated characteristics in the present study are similar.
We also note that the mean CTDs attain the amplitude of ∼2 Pa in the case of unstable+GWs+CI, in comparison to the lower amplitude in the case of Unstable+GWs case. The contribution from CI begins immediately after the mainshock and continues until 100 min from the mainshock onset. The most significant contribution comes during 40–80 min when the altitude region of the steep gradient is significantly perturbed (Figure 6 and Figure 7), i.e., when the nonlinear dynamics prevail in this altitude region. Therefore, the coupled energetics of GWs and CI can produce intense and long-lasting CTDs.
The presence of CTDs alters the refractivity ( N ) of the troposphere that takes the following form [4]:
N = 77.8 p T
where p and T are in units of millibars and kelvin, respectively. The refractivity introduces the time and phase delays in the radio signal that lead to the following Tropospheric Range Delay or TRD ( Λ ) :
Λ = C τ C 1 v d y = d y + 10 6 N d y
where we use the definition of wave velocity v = C / n and the relation between refractivity and refractive index n = 1 + 10 6 N . Therefore, the residual-refractivity ( δ N) and residual-TRD ( δ Λ ) can be written as follows:
δ N = N ( x , y , t ) N ( x , y , t = 0 ) 77.8 p T p 0 T 0 ; δ Λ = 10 6 δ N d y
Here, p 0 and T 0 are the pressure and temperature at time t = 0. Besides CTDs, the residual-TRD is another variable that can be monitored during the earthquakes [4]. The present study estimates the residual-TRD by estimating residual-refractivity from the pressure and temperature disturbances associated with the CTDs. In the observations, the residual-TRD is estimated from the phase and range delays of dual-frequency radio waves of GNSS satellites [6].
Figure 9 shows the time variation of residual-refractivity and residual-TRD, derived from Equation (7). The altitude range of integration of δ N in Equation (7) is 0–20 km. It reveals that similar to the long-lasting CTDs the residual-refractivity and residual-TRD last long after the mainshock. The residual-TRD remains larger than 0.5 mm during the evolution and attains the maximum amplitude of about 2 mm at 20 min from the mainshock onset. Recent studies report the coseismic residual-TRD of 1–3 mm during the Haida Gwaii earthquake of magnitude 7.8 and other similar magnitude earthquakes [4,5]. The residual-TRD of comparable amplitudes from these observations and present simulation affirm that the CTDs of 1–5 Pascals can introduce residual-TRD of a few millimeters in the radio signal. Therefore, the mechanism of unstable+GWs+CI, described in the present study for the generation of CTDs, can explain the observed residual-TRD.
Severe and extreme weathers, such as tropical cyclones, produce ground pressure disturbances of about 25 Pa [18]. Mega earthquakes have produced ground pressure disturbances of about 5 Pa [2]. The present numerical study reveals the CTDs of about 5–10 Pa at the tropospheric heights. Therefore, mega and strong earthquakes disturb the troposphere to a similar level as the severe meteorological weather.
The formation of CTDs is a part of atmospheric seismology that reveals their energizations by the GWs and CI. The energizations manifest as the amplification of tiny ground displacement by an order of magnitude at the tropospheric heights and formations of CTDs at a wider epicentral distance than at the ground. Such energizations favor CTDs detection at tropospheric heights during moderate and weak earthquakes that occur before the strong earthquakes. In this regard, the atmospheric seismology of CTDs contributes to the forecasting scenario of strong earthquakes.

5. Summary

The present study investigates the coupled dynamics of gravity waves and convective instability and their role in the formation of coseismic tropospheric disturbances (CTDs) and Tropospheric Range Delay (TRD). The study presents a number of numerical experiments that employ numerical code of energetics of GWs and CI in response to the ground uplift associated with an earthquake of magnitude ∼8. These experiments reveal the formation of CTDs of magnitude ∼5 Pascals, within about 20 min from the mainshock onset. It is found that the CTDs of 1–5 Pascals can introduce residual-TRD of a few millimeters in the radio signal, similar in magnitude from the recent observations. The CTDs are intense in the region of steep temperature gradients in the lower troposphere. In this altitude region, GWs attain large amplitude and the growth rate of CI is positive. The ground CTDs acquire characteristics similar to the ground pressure disturbances from earlier reports. We argue that the coseismic residual-delay in radio waves from previous reports are owing to the pressure disturbances in the altitude region of steep gradient above the epicenter. The study finds that the GWs alone can produce the CTDs on the conditionally stable troposphere. However, CTDs become stronger in the conditionally unstable troposphere owing to the convectively unstable wave energetics. The present study also reveals the formation of CTDs at off-epicenter distances which is due to the epicentral inhomogeneity of the ground vibration and subsequent nonlinear feedback that results in the atmospheric inhomogeneity and the excitation of secondary GWs.

Author Contributions

Conceptualization, E.A.K.; methodology, E.A.K.; formal analysis, E.A.K., S.A.S.; investigation, E.A.K., S.A.S.; resources, E.A.K.; writing—original draft preparation, E.A.K.; writing—review and editing, E.A.K., S.A.S., E.R.d.P.; visualization, E.A.K.; supervision, E.A.K., E.R.d.P.; project administration, E.A.K.; funding acquisition, E.A.K. All authors have read and agreed to the published version of the manuscript.

Funding

Eurico R. de Paula. acknowledges the support from Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq) under process number 202531/2019-0. Saul A. Sanchez is grateful to Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES) process 88887.351785/2019-00. Esfhan A. Kherani wishes to acknowledge the support from Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq) under process number 307496/2015-5.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Simulation data are too large (about 300 MB) in size and we are unable to upload it along with the manuscript. However, if MDPI and/or reviewers demand, we can make them available on googledrive or dropbox.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

Abbreviations

The following abbreviations are used in this manuscript:
CTDsCo-seismic Tropospheric Disturbances
GNSSGlobal Navigation Satellite System
AGWsAcoustic Gravity Waves
GWsGravity Waves
CIConvective Instability
TRDTropospheric Range Delay

Appendix A. Wave Equation of the Acoustic Gravity Waves

Campos et al. (1983) [19] and Kherani et al. (2011) [11] had derived the wave equation of AGWs of the following form:
2 u t 2 = c 2 ( . u ) + ( γ 1 ) p ρ . u p ρ u . log ρ + 1 ρ ( u . ) p
In the atmosphere, condition of strong vertical ( y ) gradient in comparison to the weak horizontal ( x ) gradient prevails. Therefore, the derivation considers ambient variations only in the vertical direction while the amplitudes ( u x , u y ) of the AGWs have variations in both directions. Under these conditions, the wave equation takes the following form:
2 u y t 2 = c 2 2 u y y 2 + c 2 2 u x y x + ( γ 1 ) ζ u y y + ( γ 1 ) ζ u x x + ζ y u y + ζ u y y
and
2 u x t 2 = c 2 2 u x x 2 + c 2 2 u y x y + ζ u y x
Here
ζ = 1 ρ p y , c 2 = γ p ρ
We consider a solution of the following form:
u y = u y ( y , t ) e i k y y + i k x x , u x = u x ( y , t ) e i k y y + i k x x
Which will reduce the derivatives into the following form:
x i k x , 2 x 2 k x 2
y i k y + y , 2 y 2 k y 2 + 2 y 2 + 2 i k y y
2 y x i k x y = i k x i k y + y = k x k y + i k x y
2 u y t 2 = d ζ dy u y k y 2 c 2 u y + c 2 2 u y y 2 + γ ζ u y y k x k y c 2 u x + i 2 k y c 2 u y y + k x c 2 u x y + γ ζ k y u y + ( γ 1 ) ζ k x u x
2 u x t 2 = k x 2 c 2 u x k x k y c 2 u y + i k x c 2 u y y + ζ u y
Equations (A5) and (A6) can be solved using the method of separation of variable such that
u y = u ys ( y ) u yt ( t ) , u x = u xs ( y ) u xt ( t )

Appendix A.1. Governing Equations in Space

Imaginary parts of (A5) and (A6) lead to the governing equations in the space of the following form:
du ys dy = k 0 u ys , k 0 = ζ c 2
and
du xs dy + ( γ 1 ) k 0 u xs + ( γ 2 ) k 0 k y k x u yt u xt u ys = 0

Appendix A.2. Governing Equations in Time

Real parts of (A5) and (A6) lead to the governing equations in the time of the following form:
d 2 u yt dt 2 = d ζ dy u yt k y 2 c 2 u yt + c 2 u yt u ys d 2 u ys dy 2 + γ ζ u yt u ys du ys dy k x k y c 2 u xs u ys u xt
and
d 2 u xt dt 2 = k x 2 c 2 u xt k x k y c 2 u ys u xs u yt
Since from (A8),
du ys dy = k 0 u ys , d 2 u ys dy 2 = k 0 2 u ys d k 0 dy u ys
Therefore,
d 2 u yt dt 2 = Ω 0 2 u yt k x k y c 2 u xs u ys u xt , d 2 u xt dt 2 = k x 2 c 2 u xt k x k y c 2 u ys u xs u yt
where
Ω 0 2 = d ζ dy + k y 2 c 2 + ( γ 1 ) k 0 2 + d k 0 dy c 2 k y 2 c 2 + Ω b 2
and
Ω b 2 = d ζ dy + ( γ 1 ) k 0 2 + d k 0 dy c 2 ( γ 1 ) k 0 2 k o c 2 d c 2 dy c 2
Here Ω b is the nonisothermal nonhydrostatic Brunt-Vaisala frequency (Equation (6.7a) of [12]).

References

  1. Oberheide, J.; Shiokawa, K.; Gurubaran, S.; Ward, W.E.; Fujiwara, H.; Kosch, M.J.; Makela, J.J.; Takahashi, H. The geospace response to variable inputs from the lower atmosphere: A review of the progress made by Task Group 4 of CAWSES-II. Prog. Earth Planet. Sci. 2015, 2, 1–31. [Google Scholar] [CrossRef] [Green Version]
  2. Mikumo, T.; Garces, M.; Shibutani, T.; Morii, W.; Okawa, T.; Ishihara, Y. Acoustic-gravity waves from the source region of the 2011 great Tohoku earthquake (M w = 9.0). J. Geophys. Res. Solid Earth 2013, 118, 1534–1545. [Google Scholar] [CrossRef] [Green Version]
  3. Watada, S.; Kunugi, T.; Hirata, K.; Sugioka, H.; Nishida, K.; Sekiguchi, S.; Oikawa, J.; Tsuji, Y.; Kanamori, H. Atmospheric pressure change associated with the 2003 Tokachi-Oki earthquake. Geophys. Res. Lett. 2006, 33, L24306. [Google Scholar] [CrossRef] [Green Version]
  4. Jin, S.; Han, L.; Cho, J. Lower atmospheric anomalies following the 2008 Wenchuan Earthquake observed by GPS measurements. J. Atmos. Sol. Terr. Phys. 2011, 73, 810–814. [Google Scholar] [CrossRef]
  5. Yao, Y.B.; Lei, X.X.; Liu, Q.; He, C.Y.; Zhang, B.; Zhang, L. Anomalies of zenith tropospheric delay following the M w 7.8 Haida Gwaii earthquake. Nat. Hazards Earth Syst. Sci. Discuss. 2014, 2, 3533–3559. [Google Scholar] [CrossRef] [Green Version]
  6. Jin, S.; Jin, R.; Liu, X. GNSS Atmospheric Seismology: Theory, Observations and Modeling; eBook; Springer Nature Singapore Pte Ltd.: Singapore, 2019. [Google Scholar] [CrossRef]
  7. Mikumo, T. Atmospheric pressure waves and tectonic deformation associated with the Alaskan Earthquake of March 28, 1964. J. Geophys. Res. 1968, 73, 2009–2025. [Google Scholar] [CrossRef] [Green Version]
  8. Le Pichon, A.; Mialle, P.; Guilbert, J.; Vergoz, J. Multistation infrasonic observations of the Chilean earthquake of 2005 June 13. Geophys. J. Int. 2006, 167, 838–844. [Google Scholar] [CrossRef] [Green Version]
  9. Lognonné, P. Seismic waves from atmospheric sources and atmospheric/ionospheric signatures of seismic waves. In Infrasound Monitoring for Atmospheric Studies; Springer: Berlin/Heidelberg, Germany, 2009; pp. 281–308. [Google Scholar]
  10. Kherani, E.A.; Lognonné, P.; Hébert, H.; Rolland, L.; Astafyeva, E.; Occhipinti, G.; Coïsson, P.; Walwer, D.; de Paula, E.R. Modelling of the total electronic content and magnetic field anomalies generated by the 2011 Tohoku-Oki tsunami and associated acoustic-gravity waves. Geophys. J. Int. 2012, 191, 1049–1066. [Google Scholar] [CrossRef] [Green Version]
  11. Kherani, E.A.; Abdu, M.A.; Fritts, D.C.; de Paula, E.R. The acoustic gravity wave induced disturbances in the equatorial Ionosphere. In Aeronomy of the Earth’s Atmosphere and Ionosphere; Springer: Dordrecht, The Netherlands, 2011; pp. 141–162. [Google Scholar] [CrossRef]
  12. Kelley, M.C. The Earth’s Ionosphere: Plasma Physics and Electrodynamics, 2nd ed.; Academic Press: San Diego, CA, USA, 2009. [Google Scholar]
  13. Chun, H.; Baik, J. Momentum flux by thermally induced internal gravity waves and its approximation for large-scale models. J. Atmos. Sci. 1998, 55, 3299–3310. [Google Scholar] [CrossRef]
  14. Beres, J.H.; Alexander, M.J.; Holton, J.R. Effects of tropospheric wind shear on the spectrum of convectively generated gravity waves. J. Atmos. Sci. 2002, 59, 1805–1824. [Google Scholar] [CrossRef] [Green Version]
  15. Fovell, R.; Durran, D.; Holton, J.R. Numerical simulations of convectively generated stratospheric gravity waves. J. Atmos. Sci. 1992, 49, 1427–1442. [Google Scholar] [CrossRef] [Green Version]
  16. Lane, T.P.; Reeder, M.J.; Clark, T.L. Numerical modeling of gravity-wave generation by deep tropical convection. J. Atmos. Sci. 2001, 58, 1249–1274. [Google Scholar] [CrossRef]
  17. Picone, J.; Hedin, A.; Drob, D.; Aiken, A. NRLMSISE-00 empirical model of the atmosphere: Statistical comparisons and scientific issues. J. Geophys. Res. 2002, 107, 1468. [Google Scholar] [CrossRef]
  18. Shi, L.; Olabarrieta, M.; Nolan, D.S.; Warner, J.C. Tropical cyclone rainbands can trigger meteotsunamis. Nat. Commun. 2020, 11, 678. [Google Scholar] [CrossRef] [PubMed]
  19. Campos, L.M.B.C. On three dimnsional Acoustic-gravity waves in model non-isothermal atmospheres. Wave Motion 1983, 5, 1–14. [Google Scholar] [CrossRef]
Figure 1. The ground vertical velocity V SISM associated with the earthquake of magnitude 8 (Here the seismic waveform of 2008 Wenchuan earthquake is considered).
Figure 1. The ground vertical velocity V SISM associated with the earthquake of magnitude 8 (Here the seismic waveform of 2008 Wenchuan earthquake is considered).
Atmosphere 12 00765 g001
Figure 2. Ambient conditions derived from the NRLMSISE model. In (A,B), atmospheric density and temperature are shown. Two temperature profiles represent covectively stable and unstable troposphere.
Figure 2. Ambient conditions derived from the NRLMSISE model. In (A,B), atmospheric density and temperature are shown. Two temperature profiles represent covectively stable and unstable troposphere.
Atmosphere 12 00765 g002
Figure 3. Temperature lapse rate ( γ e ) in red and adiabatic lapse rate ( γ a d ) in green, are shown. The dashed and solid curves correspond to the stable and unstable temperature profiles.
Figure 3. Temperature lapse rate ( γ e ) in red and adiabatic lapse rate ( γ a d ) in green, are shown. The dashed and solid curves correspond to the stable and unstable temperature profiles.
Atmosphere 12 00765 g003
Figure 4. Frequency spectrum of the amplitude u y of GWs at each 2 km in 0–20 km altitude range. The spectral amplitudes are normalized to the maximum spectral amplitude at the ground. For better visualization, the line plots are multiplied by 0.5. The black circles represent the Brunt-Vaisala periods.
Figure 4. Frequency spectrum of the amplitude u y of GWs at each 2 km in 0–20 km altitude range. The spectral amplitudes are normalized to the maximum spectral amplitude at the ground. For better visualization, the line plots are multiplied by 0.5. The black circles represent the Brunt-Vaisala periods.
Atmosphere 12 00765 g004
Figure 5. Vertical Travel Time diagram above the epicenter: It shows the CTDs from the stable GWs case and unstable GWs case in (A,B) respectively. CTDs defined as expression (6) is plotted in two types of format: as line plot at each height and as color pixmap. For better visualization, CTDs in line plot is multiplied by factor 0.2.
Figure 5. Vertical Travel Time diagram above the epicenter: It shows the CTDs from the stable GWs case and unstable GWs case in (A,B) respectively. CTDs defined as expression (6) is plotted in two types of format: as line plot at each height and as color pixmap. For better visualization, CTDs in line plot is multiplied by factor 0.2.
Atmosphere 12 00765 g005
Figure 6. Vertical Travel Time diagram above the epicenter: it shows the CTDs from the unstable GWs+CI case. CTDs defined as expression (6) is plotted in two types of format: as line plot at each height and as color pixmap. For better visualization, CTDs in line plot is multiplied by factor 0.2.
Figure 6. Vertical Travel Time diagram above the epicenter: it shows the CTDs from the unstable GWs+CI case. CTDs defined as expression (6) is plotted in two types of format: as line plot at each height and as color pixmap. For better visualization, CTDs in line plot is multiplied by factor 0.2.
Atmosphere 12 00765 g006
Figure 7. It presents snapshots of spatial distribution of CTDs as pixmaps and growth rate of CI as contours for the Unstable GWs+CI case. The growth rate increases from yellow to red colors.
Figure 7. It presents snapshots of spatial distribution of CTDs as pixmaps and growth rate of CI as contours for the Unstable GWs+CI case. The growth rate increases from yellow to red colors.
Atmosphere 12 00765 g007
Figure 8. It presents the time vs. pressure disturbances over the epicenter. The curves correspond to the ground pressure disturbance and mean pressure disturbance from lower tropospheric altitudes between 1 and 15 km.
Figure 8. It presents the time vs. pressure disturbances over the epicenter. The curves correspond to the ground pressure disturbance and mean pressure disturbance from lower tropospheric altitudes between 1 and 15 km.
Atmosphere 12 00765 g008
Figure 9. It demonstrates the time vs. Residual-Refractivity and time vs. Residual-TRD over the epicenter.
Figure 9. It demonstrates the time vs. Residual-Refractivity and time vs. Residual-TRD over the epicenter.
Atmosphere 12 00765 g009
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Kherani, E.A.; Sanchez, S.A.; de Paula, E.R. Numerical Modeling of Coseismic Tropospheric Disturbances Arising from the Unstable Acoustic Gravity Wave Energetics. Atmosphere 2021, 12, 765. https://doi.org/10.3390/atmos12060765

AMA Style

Kherani EA, Sanchez SA, de Paula ER. Numerical Modeling of Coseismic Tropospheric Disturbances Arising from the Unstable Acoustic Gravity Wave Energetics. Atmosphere. 2021; 12(6):765. https://doi.org/10.3390/atmos12060765

Chicago/Turabian Style

Kherani, Esfhan A., Saul A. Sanchez, and Eurico R. de Paula. 2021. "Numerical Modeling of Coseismic Tropospheric Disturbances Arising from the Unstable Acoustic Gravity Wave Energetics" Atmosphere 12, no. 6: 765. https://doi.org/10.3390/atmos12060765

APA Style

Kherani, E. A., Sanchez, S. A., & de Paula, E. R. (2021). Numerical Modeling of Coseismic Tropospheric Disturbances Arising from the Unstable Acoustic Gravity Wave Energetics. Atmosphere, 12(6), 765. https://doi.org/10.3390/atmos12060765

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop