GOCE Mission Research Articles

Gravity Spectra from the Density Distribution of Earth’s Uppermost 435 km

The Earth masses reside in a near-hydrostatic equilibrium, while the deviations are, for example, manifested in the geoid, which is nowadays well determined by satellite gravimetry. Recent progress in estimating the density distribution of the Earth allows us to examine individual Earth layers and to directly see how the sum approaches the observed anomalous gravitational field. This study evaluates contributions from the crust and the upper mantle taken from the LITHO1.0 model and quantifies the gravitational spectra of the density structure to the depth of 435 km. This is done without isostatic adjustments to see what can be revealed with models like LITHO1.0 alone. At the resolution of 290 km (spherical harmonic degree 70), the crustal contribution starts to dominate over the upper mantle and at about 150 km (degree 130) the upper mantle contribution is nearly negligible. At the spatial resolution $$<150\,\hbox {km},$$ the spectra behavior is driven by the crust, the mantle lid and the asthenosphere. The LITHO1.0 model was furthermore referenced by adding deeper Earth layers from ak135, and the gravity signal of the merged model was then compared with the observed satellite-only model GOCO05s. The largest differences are found over the tectonothermal cold and old (such as cratonic), and over warm and young areas (such as oceanic ridges). The misfit encountered comes from the mantle lid where a velocity–density relation helped to reduce the RMS error by 40%. Global residuals are also provided in terms of the gravitational gradients as they provide better spatial localization than gravity, and there is strong observational support from ESA’s satellite gradiometry mission GOCE down to the spatial resolution of 80–90 km.

Applying the GOCE-based GGMs for the quasi-geoid modelling of Finland

AbstractThe gravity satellite mission GOCE made its final observations in the fall of 2013. Since the re-entry to the Earth’s atmosphere, the full cycle of the GOCE data has been published by ESA. At first, we evaluated all the GOCE-based global geoid models over Finland using terrestrial gravity and GNSS-levelling data. The most suitable model was selected as a global background model for the Finnish quasi-geoid calculations.Next, we combined the chosen model with terrestrial gravity data of Finland and surrounding areas. Quasi-geoid models with different modifications were calculated using the GOCE DIR5 model up to spherical harmonic degree and order (d/o) 240 and 300, and the high resolution EIGEN-6C4 (includes the complete GOCE data) model up to degree and order 1000 and 2190.The calculated quasi-geoid models were validated to the measurements on site with two independent GPS-levelling datasets. The best quasi-geoid models with GOCE gave standard deviations of 2.6 cm (FIN_DIR5 d/o 240) and 2.3 cm (FIN_DIR5 d/o 300) in Finland. For the high resolution model FIN_EIGEN-6C4, the results were 1.8 cm (d/o 1000) and 1.7 cm (d/o 2190). In addition, the results were compared with the latest geoid models available in Finland (FIN2005N00, NKG2004, NKG2015, EGG2008). The sub-2-centimetre (and near 2 cm, when using the GOCE-based models) accuracy is an improvement over the previous and current Finnish geoid models.

A re-evaluation of the relativistic redshift on frequency standards at NIST, Boulder, Colorado, USA**Contribution of US government, not subject to Copyright.

We re-evaluated the relativistic redshift correction applicable to the frequency standards at the National Institute of Standards and Technology (NIST) in Boulder, Colorado, USA, based on a precise GPS survey of three benchmarks on the roof of the building where these standards had been previously housed, and on global and regional geoid models supported by data from the GRACE and GOCE missions, including EGM2008, USGG2009, and USGG2012. We also evaluated the redshift offset based on the published NAVD88 geopotential number of the leveling benchmark Q407 located on the side of Building 1 at NIST, Boulder, Colorado, USA, after estimating the bias of the NAVD88 datum at our specific location. Based on these results, our current best estimate of the relativistic redshift correction, if frequency standards were located at the height of the leveling benchmark Q407 outside the second floor of Building 1, with respect to the EGM2008 geoid whose potential has been estimated to be , is equal to (−1798.50 ± 0.06) × 10−16. The corresponding value, with respect to an equipotential surface defined by the International Astronomical Union’s (IAU) adopted value of , is (−1798.53 ± 0.06) × 10−16. These values are comparable to the value of (−1798.70 ± 0.30) × 10−16, estimated by Pavlis and Weiss in 2003, with respect to an equipotential surface defined by . The minus sign implies that clocks run faster in the laboratory in Boulder than a corresponding clock located on the geoid.

Gravity field models derived from the second degree radial derivatives of the GOCE mission: a case study

&lt;p class="Default"&gt;The GOCE satellite mission is one of the main achievements of the satellite geodesy for the Earth’s gravitational field recovery. Three different approaches have been developed for the estimation of harmonic coefficients from gradiometry data measured on board of GOCE-satellite. In this paper a special version of the space-wise method based on the second method of Neumann for fast determination of the harmonic coefficients &lt;em&gt;C&lt;sub&gt;nm&lt;/sub&gt;, S&lt;sub&gt;nm&lt;/sub&gt;&lt;/em&gt; of the Earth’s gravitational potential is given based on the radial gravity gradients of the EGG_TRF_2 product, except of two polar gaps filled by radial gradients from the EGM2008 gravity model. In the pre-processing stage GOCE-based second degree radial derivatives were averaged to the regular grid through Kalman static filter with additional Gaussean smoothing of residual radial derivatives. All computations are made by iterations. As the first step the determination of the preliminary NULP-01S&lt;strong&gt; &lt;/strong&gt;model up to degree/order 220 derived from the Gaussean grid of the GOCE radial derivatives with respect to the WGS-84 reference field was developed based only one of the radial gradients EGG_TRF_2 in the EFRF-frame. In the second iteration the same algorithm is applied to build the NULP-02S gravity field model up to degree/order 250 using the same Gaussean grid with respect to the NULP-01S reference field. The NULP-02S model was verified by means of applying various approaches for the construction of the gridded gravity anomalies and geoid heights in the Black sea area using processing of datasets from six altimetry satellite missions. Comparison of different models with GNSS-levelling data in the USA area demonstrates the independent verification of achieved accuracy of the constructed NULP-02S Earth’s gravity field model.&lt;/p&gt;

Space-Wise approach for airborne gravity data modelling

Regional gravity field modelling by means of remove-compute-restore procedure is nowadays widely applied in different contexts: it is the most used technique for regional gravimetric geoid determination, and it is also used in exploration geophysics to predict grids of gravity anomalies (Bouguer, free-air, isostatic, etc.), which are useful to understand and map geological structures in a specific region. Considering this last application, due to the required accuracy and resolution, airborne gravity observations are usually adopted. However, due to the relatively high acquisition velocity, presence of atmospheric turbulence, aircraft vibration, instrumental drift, etc., airborne data are usually contaminated by a very high observation error. For this reason, a proper procedure to filter the raw observations in both the low and high frequencies should be applied to recover valuable information. In this work, a software to filter and grid raw airborne observations is presented: the proposed solution consists in a combination of an along-track Wiener filter and a classical Least Squares Collocation technique. Basically, the proposed procedure is an adaptation to airborne gravimetry of the Space-Wise approach, developed by Politecnico di Milano to process data coming from the ESA satellite mission GOCE. Among the main differences with respect to the satellite application of this approach, there is the fact that, while in processing GOCE data the stochastic characteristics of the observation error can be considered a-priori well known, in airborne gravimetry, due to the complex environment in which the observations are acquired, these characteristics are unknown and should be retrieved from the dataset itself. The presented solution is suited for airborne data analysis in order to be able to quickly filter and grid gravity observations in an easy way. Some innovative theoretical aspects focusing in particular on the theoretical covariance modelling are presented too. In the end, the goodness of the procedure is evaluated by means of a test on real data retrieving the gravitational signal with a predicted accuracy of about 0.4 mGal.

Analysis of the Illapel Mw = 8.3 Thrust Earthquake Rupture Zone Using GOCE-Derived Gradients

Satellite gravimetry has proven to be a useful tool to identify mass anomalies along a subduction interface, interpreted as heterogeneities related to the rupture process during megathrust earthquakes. In the last years, different works, reinforced with data derived from satellite gravity missions as GRACE and now GOCE, have analyzed not only the static component of the Earth gravity field, but also its temporal variations and relation to the seismic cycle.In particular, during the last decade, the Chilean margin has been affected by three megathrust earthquakes (withMw>8): Maule 2010 Mw = 8.8, Pisagua 2014 Mw = 8.2 and recently the Mw = 8.3 Illapel event. Then, the recently completed GOCE mission (November 2009 to November 2013) offered a unique opportunity to study the Maule February 2010 and Pisagua April 2014 events by means of gravity gradients, directly measured at satellite height altitudes, which allowed mapping density heterogeneities with greater detail than the gravity anomaly which has been used in most studies up to now. In the present work, we use the last GOCE model (GO_CONS_GCF_2_DIR_R5), the one of higher spatial resolution (N = 300, λ/2 ≈ 66 km) derived from satellite-only data. The methodology used is the same as that to study the previous events, with the addition that now we derived a relation between the associated depths of a causative mass with a determined degree of the spherical harmonic expansion. This allowed to ‘‘decompose’’ the gravimetric signal, by cutting off the degree/order of the harmonic expansion, as depth increases. From this analysis, we found that prominent oceanic features such as the Challenger fracture zone and the Juan Fernandez ridge played a key role in latitudinal seismic segmentation for the Illapel earthquake rupture zone, acting as barriers/attenuators to the seismic energy release. We compared the slip model from Tilmann et al. (Geophysical Research Letters 43: 574–583. doi: 10.1002/2015GL066963, 2016) for the Illapel earthquake with vertical gravity gradient with and without sediment correction, and at different degree/order of the harmonic expansion. From this analysis, we inferred that prominent oceanic features over the subducting Nazca plate play a key role in seismic segmentation not only at heavily sedimented trenches, but also at sediment-starved segments.

Comparison of Satellite-Only Gravity Field Models Constructed with All and Parts of the GOCE Gravity Gradient Dataset

ABSTRACTThe impact of GOCE Satellite Gravity Gradiometer data on gravity field models was tested. All models were constructed with the same Laser Geodynamics Satellite (LAGEOS) and Gravity Recovery and Climate Experiment (GRACE) data, which were combined with one or two of the diagonal gravity gradient components for the entire GOCE mission (November 2009 to October 2013). The Stokes coefficients were estimated by solving large normal equation (NE) systems (i.e., the direct numerical approach). The models were evaluated through comparisons with the European Space Agency's (ESA) gravity field model DIR-R5, by GPS/Leveling, GOCE orbit determination, and geostrophic current evaluations. Among the single gradient models, only the model constructed with the vertical ZZ gradients gave good results that were in agreement with the formal errors. The model based only on XX gradients is the least accurate. The orbit results for all models are very close and confirm this finding. All models constructed with two diagonal gradient components are more accurate than the ZZ-only model due to doubling the amount of data and having two complementary observation directions. This translates also to a slower increase of model errors with spatial resolution. The different evaluation methods cannot unambiguously identify the most accurate two-component model. They do not always agree, emphasizing the importance of evaluating models using many different methods. The XZ gravity gradient gives a small positive contribution to model accuracy.

Spherical Harmonic Analysis of Gravitational Curvatures and Its Implications for Future Satellite Missions

In this study we assume that a gravitational curvature tensor, i.e. a tensor of third-order directional derivatives of the Earth’s gravitational potential, is observable at satellite altitudes. Such a tensor is composed of ten different components, i.e. gravitational curvatures, which may be combined into vertical–vertical–vertical, vertical–vertical–horizontal, vertical–horizontal–horizontal and horizontal–horizontal-horizontal gravitational curvatures. Firstly, we study spectral properties of the gravitational curvatures. Secondly, we derive new quadrature formulas for the spherical harmonic analysis of the four gravitational curvatures and provide their corresponding analytical error models. Thirdly, requirements for an instrument that would eventually observe gravitational curvatures by differential accelerometry are investigated. The results reveal that measuring third-order directional derivatives of the gravitational potential imposes very high requirements on the accuracy of deployed accelerometers which are beyond the limits of currently available sensors. For example, for orbital parameters and performance similar to those of the GOCE mission, observing third-order directional derivatives requires accelerometers with the noise level of \({\sim}10^{-17}\,\hbox {m}\,\hbox {s}^{-2}\) Hz\(^{-1/2}\).

Swarm kinematic orbits and gravity fields from 18 months of GPS data

The Swarm mission consists of three satellites orbiting the Earth at low orbital altitudes. The onboard GPS receivers, star cameras, and laser retro-reflectors make the Swarm mission an interesting candidate to explore the contribution of Swarm GPS data to the recovery of both the static and time-variable gravity fields. We use 1.5years of Swarm GPS and attitude data to generate kinematic positions of high quality to perform gravity field determination using the Celestial Mechanics Approach. The generated gravity fields reveal severe systematic errors along the geomagnetic equator. Their size is correlated with the ionospheric density and thus strongly varying over the analyzed time period. Similar to the findings of the GOCE mission, the systematic errors are related to the Swarm GPS carrier phase data and may be reduced by rejecting GPS data affected by large ionospheric changes. Such a measure yields a strong reduction of the systematic errors along the geomagnetic equator in the gravity field recovery. Long wavelength signatures of the gravity field may then be recovered with a similar quality as achieved with GRACE GPS data, which makes the Swarm mission well suited to bridge a potential gap between the current GRACE and the future GRACE Follow-On mission.

A new generation of ultra-sensitive electrostatic accelerometers for GRACE Follow-on and towards the next generation gravity missions

The sensor core configuration of the electrostatic accelerometers of the CHAMP, GRACE and GOCE missions has been especially designed for space applications and so optimized in regard to the weak level of acceleration to sustain and measure in orbit. The return of experience of these three gravity space missions which demonstrate the robustness and the performances of this family of space instruments, allows a better optimization of the design of the accelerometer in terms of thermal stability and operation reliability for the near future missions as GRACE follow-on. The paper presents the improvement of the GRACE-FO accelerometer with respect to the still in-orbit previous models and a status of its development.Nevertheless because the next generation of low-low satellite to satellite tracking missions will take advantage of interferometer laser ranging methods to improve their performance, the noise level of the accelerometer has also to be lowered, especially in the low frequency bandwidth. In addition to the measurement of the surface forces exerted on the spacecraft by the atmospheric drag, by the Sun radiation and by the Earths albedo and infrared pressures, the accelerometer instrument becomes a major part of the attitude and orbit control system by acting as a drag free sensor and by accurately measuring the angular accelerations. Onera proposes a new configuration much more compact with performances in a better adequacy with a next generation of small but drag compensated micro-satellites or geodesy missions.

Evaluation of gravity field model EIGEN-6C4 by means of various functions of gravity potential, and by GNSS/levelling

&lt;p&gt;The combined gravity field model EIGEN-6C4 (Förste et al., 2014) is the latest combined global gravity field model of GFZ Potsdam and GRGS Toulouse. EIGEN-6C4 has been generated including the satellite gravity gradiometry data of the entire GOCE mission (November 2009 till October 2013) and is of maximum spherical degree and order 2190. In this study EIGEN-6C4 has been compared with EGM2008 to its maximum degree and order via gravity disturbances and T&lt;sub&gt;zz&lt;/sub&gt; part of the Marussi tensor of the second derivatives of the disturbing potential. The emphasis is put on such areas where GOCE data (complete set of gradiometry measurements after reductions) in EIGEN-6C4 obviously contributes to an improvement of the gravity field description. &lt;/p&gt;&lt;p&gt;GNSS/levelling geoid heights are independent data source for the evaluation of gravity field models. Therefore, we use the GNSS/levelling data sets over the territories of Europe, Czech Republic and Slovakia for the evaluation of EIGEN-6C4 w.r.t. EGM2008.&lt;/p&gt;

Moho topography, ranges and folds of Tibet by analysis of global gravity models and GOCE data.

The determination of the crustal structure is essential in geophysics, as it gives insight into the geohistory, tectonic environment, geohazard mitigation, etc. Here we present the latest advance on three-dimensional modeling representing the Tibetan Mohorovičić discontinuity (topography and ranges) and its deformation (fold), revealed by analyzing gravity data from GOCE mission. Our study shows noticeable advances in estimated Tibetan Moho model which is superior to the results using the earlier gravity models prior to GOCE. The higher quality gravity field of GOCE is reflected in the Moho solution: we find that the Moho is deeper than 65 km, which is twice the normal continental crust beneath most of the Qinghai-Tibetan plateau, while the deepest Moho, up to 82 km, is located in western Tibet. The amplitude of the Moho fold is estimated to be ranging from −9 km to 9 km with a standard deviation of ~2 km. The improved GOCE gravity derived Moho signals reveal a clear directionality of the Moho ranges and Moho fold structure, orthogonal to deformation rates observed by GPS. This geophysical feature, clearly more evident than the ones estimated using earlier gravity models, reveals that it is the result of the large compressional tectonic process.

Технология определения параметров гравитационного поля Земли по градиентометрическим измерениям. 2. Системы координат и времени в спутниковой градиентометрии

Abstract: This paper is the second of paper series devoted to the technology description of the gravity field determination using satellite gradiometry data. The objective of this paper is the description of reference frames and time systems used for mathematical processing gradiometry data of mission GOCE.

Iterative Spherical Downward Continuation Applied to Magnetic and Gravitational Data from Satellite

In the last few decades, satellites have acquired various potential data sets hundreds of kilometers above the Earth’s surface. Conventionally, these global magnetic and gravitational data sets are approximated by using spherical harmonics that allow straightforward work with both fields outside the Earth’s mass. In this article, we present an alternative approach for working with potential data in mass-free space given over a regular coordinate grid on a spherical surface. The algorithm is based on an iterative scheme and the Poisson integral equation for the sphere. With help from the Fourier transform, global potential (magnetic or gravitational) data can efficiently be continued from a mean orbital sphere down to a reference surface without using the spherical harmonics. This is illustrated both with simulated magnetic field data and with real data from the satellite gradiometry mission GOCE. In the case of simulated magnetic data and the downward continuation for 450 km, we have achieved a root mean square at the level of 0.05 nT, while it was <1 E (eotvos) for real GOCE data continued for 250 km. The crucial point is to apply the algorithm twice as a large part of noise can be removed from the input data.

The construction of gravitational field models on the basis of satellite measurements of gravitational potential derivatives

In [1] expressions were constructed for the derivatives of all the orders of a planet’s gravitational potential with respect to the rectangular coordinates related to the gravity center of a planet. These expressions are series of spherical functions. The coefficients of the series of first-order derivatives depend on two Stokes constants, whereas the coefficients of next-order derivatives are linear combinations of the coefficients of preceding-order derivatives. In the present paper the derived expressions for the first and second potential derivatives are transformed into the form that is most convenient for solving the inverse problem, i.e., evaluating Stokes constants from satellite measurements of these derivatives. Each term of the new series for a derivative depends on a sum of two Stokes constants multiplied by linear combinations of several spherical functions. The new form of the expansions for the potential derivatives makes it possible to calculate Stokes constants by simultaneously applying satellite data either for all three first-order potential derivatives, or for all six second-order derivatives. The constructed series may be applied for modeling the Earth’s gravitational field from the satellite data obtained in the international CHAMP, GRACE, and GOCE missions.

The computation of the geoid model in the state of São Paulo using two methodologies and GOCE models

The purpose of this manuscript is to compute and to evaluate the geoid model in the State of São Paulo from two methodologies (Stokes' integral through the Fast Fourier Transform - FFT and Least Squares Collocation - LSC). Another objective of this study is to verify the potentiality of GOCE-based. A special attention is given to GOCE mission. The theory related to Stokes' integral and Least Squares Collocation is also discussed in this work. The spectral decomposition was employed in the geoid models computation and the long wavelength component was represented by EGM2008 up to degree and order 150 and 360 and GOCE-based models up to 150. The models were compared in terms of geoid height residual and absolute and relative comparisons from GPS/leveling and the results show consistency between them. In addition, a comparison in the mountain regions was carried out to verify the methodologies behavior in this area; the results showed that LSC is less consistent than FFT.

Spherical integral formulas for upward/downward continuation of gravitational gradients onto gravitational gradients

New integral formulas for upward/downward continuation of gravitational gradients onto gravitational gradients are derived in this article. They provide more options for continuation of gravitational gradient combinations and extend available mathematical apparatus formulated for this purpose up to now. The starting point represents the analytical solution of the spherical gradiometric boundary value problem in the spatial domain. Applying corresponding differential operators on the analytical solution of the spherical gradiometric boundary value problem, a total of 18 integral formulas are provided. Spatial and spectral forms of isotropic kernels are given and their behaviour for parameters of a GOCE-like satellite is investigated. Correctness of the new integral formulas and the isotropic kernels is tested in a closed-loop simulation. The derived integral formulas and the isotropic kernels form a theoretical basis for validation purposes and geophysical applications of satellite gradiometric data as provided currently by the GOCE mission. They also extend the well-known Meissl scheme.

The static gravity field model DGM-1S from GRACE and GOCE data: computation, validation and an analysis of GOCE mission’s added value

We present a global static model of the Earth’s gravity field entitled DGM-1S based on GRACE and GOCE data. The collection of used data sets includes nearly 7 years of GRACE KBR data and 10 months of GOCE gravity gradient data. The KBR data are transformed with a 3-point differentiation into quantities that are approximately inter-satellite accelerations. Gravity gradients are processed in the instrumental frame. Noise is handled with a frequency-dependent data weighting. DGM-1S is complete to spherical harmonic degree 250 with a Kaula regularization being applied above degree 179. Its performance is compared with a number of other satellite-only GRACE/GOCE models by confronting them with (i) an independent model of the oceanic mean dynamic topography, and (ii) independent KBR and gravity gradient data. The tests reveal a competitive quality for DGM-1S. Importantly, we study added value of GOCE data by comparing the performance of satellite-only GRACE/GOCE models with models produced without GOCE data: either ITG-Grace2010s or EGM2008 depending on which of the two performs better in a given region. The test executed based on independent gravity gradients quantifies this added value as 25–38 % in the continental areas poorly covered with terrestrial gravimetry data (Equatorial Africa, Himalayas, and South America), 7–17 % in those with a good coverage with these data (Australia, North America, and North Eurasia), and 14 % in the oceans. This added value is shown to be almost entirely related to coefficients below degree 200. It is shown that this gain must be entirely attributed to gravity gradients acquired by the mission. The test executed based on an independent model of the mean dynamic topography suggests that problems still seem to exist in satellite-only GRACE/GOCE models over the Pacific ocean, where noticeable deviations between these models and EGM2008 are detected, too.

Gravitational Gradients at Satellite Altitudes in Global Geophysical Studies

Due to the ESA’s satellite mission GOCE launched in March 2009, gravitational gradients sampled along the orbital trajectory approximately 250 km above the Earth’s surface have become available. Since 2010, gravitational gradients have routinely been applied in geodesy for the derivation of global Earth’s gravitational models provided in terms of fully normalized coefficients in a spherical harmonic series representation of the Earth’s gravitational potential. However, in geophysics, gravitational gradients observed by spaceborne instruments have still been applied relatively seldom. This contribution describes their possible geophysical applications in structural studies where gravitational gradients observed at satellite altitudes are compared with those derived by a spectral forward modeling technique using available global models of selected Earth’s mass components as input data. In particular, GOCE gravitational gradients are interpreted in terms of a superposition principle of gravitation as combined gravitational effects generated by a homogeneous reference ellipsoid of revolution, mean topographic and ice mass density distributions, depth-dependent mass density contrasts within bathymetry and lateral mass density anomalies with sediments and crustal layers. Respective gravitational effects are one by one removed from gravitational gradients observed at approximately 250 km elevation above ground. Removing respective gravitational gradients from observed gravitational gradients gradually reveals problematic geographic areas with model deficiencies. For the full interpretation of observed gravitational gradients, deficiencies of CRUST2.0 must be corrected and effects of deeper laying mass anomalies not included in the study considered. These findings are confirmed by parameters describing spectral properties of the gravitational gradients. The methodology can be applied for validating Earth’s gravitational models and for constraining crustal models in the development phase.

GOCE: The first seismometer in orbit around the Earth

The first in situ sounding of a post‐seismic infrasound wavefront is presented, using data from the GOCE mission. The atmospheric infrasounds following the great Tohoku earthquake (on 11 March 2011) induce variations of air density and vertical acceleration of the GOCE platform. These signals are detected at two positions along the GOCE orbit corresponding to a crossing and a doubling of the infrasonic wavefront created by seismic surface waves. Perturbations up to 11% of air density and 1.35 × 10 − 7 m/s2 of vertical acceleration are observed and modeled with two different solid‐atmosphere coupling codes. These perturbations are a due to acoustic waves creating vertical velocities up to 130 m/s. Amplitudes and arrival times of these perturbations are reproduced respectively within a factor 2, and within a 60 s time window. Waveforms present a good agreement with observed data. The vertical acceleration to air density perturbation ratio is higher for these acoustic waves than for gravity waves. Combining these two pieces of information offers a new way to distinguish between these two wave types. This new type of data is a benchmark for the models of solid‐atmosphere coupling. Amplitude and frequency content constrain the infrasound attenuation related to atmosphere viscosity and thermal conductivity. Observed time shifts between data and synthetics are ascribed to lateral variations of the seismic and atmospheric sound velocities and to the influence of atmospheric winds. These effects should be included in future modeling. This validation of our modeling tools allows to specify more precisely future observation projects.