Assimilating atmospheric infrasound data to constrain atmospheric winds in a two‐dimensional grid

Abstract

AbstractInfrasound waves travelling through atmospheric channels are affected by the conditions they encounter along their path. The shift in the back azimuth angle of a wavefront detected at the reception site depends on the cross‐winds it encountered. Estimating the original field from this integrated measurement is an (ill‐posed) inverse problem. By using a prior, this can be converted into a Bayesian estimation problem. In this work, we use the (ensemble) Kalman filter (EnKF) to tackle this problem. In particular, we provide an illustration of the setup and solution of the problem in a two‐dimensional grid, depending on both across‐track distance and height, which has not been done in previous works. We use a synthetic setup to discuss the details of the method. We show that one of the effects of along‐track averaging (done in previous studies to simplify the problem) is to overestimate the magnitude of the analysed values, and propose that this will be a source of model error. We also illustrate the process with real data corresponding to nine controlled ammunition explosions that took place in the summer of 2018. In these cases, the real infrasound waves we study seldom reach higher than km in height. However, the use of covariance‐based methods (e.g., the EnKF) allows for updates in higher regions where the wave did not travel and where traditional observations are sparse. In fact, the larger impacts from observations in these cases are in the region of 40–60 km, in agreement with previous works. This study contributes to paving the way towards the ultimate goal of assimilating information derived from infrasound waves into operational numerical weather forecasting. More studies in quality control of the observations and proper validation of the results are urgently needed.