skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: Four-dimensional mesospheric and lower thermospheric wind fields using Gaussian process regression on multistatic specular meteor radar observations
Abstract. Mesoscale dynamics in the mesosphere and lower thermosphere (MLT) region have been difficult to study from either ground- or satellite-based observations. For understanding of atmospheric coupling processes, important spatial scales at these altitudes range between tens and hundreds of kilometers in the horizontal plane. To date, this scale size is challenging observationally, so structures are usually parameterized in global circulation models. The advent of multistatic specular meteor radar networks allows exploration of MLT mesoscale dynamics on these scales using an increased number of detections and a diversity of viewing angles inherent to multistatic networks. In this work, we introduce a four-dimensional wind field inversion method that makes use of Gaussian process regression (GPR), which is a nonparametric and Bayesian approach. The method takes measured projected wind velocities and prior distributions of the wind velocity as a function of space and time, specified by the user or estimated from the data, and produces posterior distributions for the wind velocity. Computation of the predictive posterior distribution is performed on sampled points of interest and is not necessarily regularly sampled. The main benefits of the GPR method include this non-gridded sampling, the built-in statistical uncertainty estimates, and the ability to horizontally resolve winds on relatively small scales. The performance of the GPR implementation has been evaluated on Monte Carlo simulations with known distributions using the same spatial and temporal sampling as 1 d of real meteor measurements. Based on the simulation results we find that the GPR implementation is robust, providing wind fields that are statistically unbiased with statistical variances that depend on the geometry and are proportional to the prior velocity variances. A conservative and fast approach can be straightforwardly implemented by employing overestimated prior variances and distances, while a more robust but computationally intensive approach can be implemented by employing training and fitting of model hyperparameters. The latter GPR approach has been applied to a 24 h dataset and shown to compare well to previously used homogeneous and gradient methods. Small-scale features have reasonably low statistical uncertainties, implying geophysical wind field horizontal structures as low as 20–50 km. We suggest that this GPR approach forms a suitable method for MLT regional and weather studies.  more » « less
Award ID(s):
1933005
PAR ID:
10391199
Author(s) / Creator(s):
; ; ; ; ;
Date Published:
Journal Name:
Atmospheric Measurement Techniques
Volume:
14
Issue:
11
ISSN:
1867-8548
Page Range / eLocation ID:
7199 to 7219
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract On 15 January 2022, the Hunga volcano produced a massive explosion that generated perturbations in the entire atmosphere. Nonetheless, signatures in the mesosphere and lower thermosphere (MLT) have been challenging to identify. We report MLT horizontal wind perturbations using three multistatic specular meteor radars on the west side of South America (spanning more than 3,000 km). The most notorious signal is an exceptional solitary wave with a large vertical wavelength observed around 18 UT at all three sites, with an amplitude of ∼50 m/s mainly in the westward direction. Using a customized analysis, the wave is characterized as traveling at ∼200 m/s, with a period of ∼2 hr and a horizontal wavelength of ∼1,440 km in the longitudinal direction, away from the source. The perturbation is consistent with anL1Lamb wave mode. The signal's timing coincides with the arrival time of the tsunami triggered by the eruption. 
    more » « less
  2. Abstract. The Hunga Tonga–Hunga Ha′apai volcano eruption was a unique event that caused many atmospheric phenomena around the globe. In this study, we investigate the atmospheric gravity waves in the mesosphere/lower-thermosphere (MLT) launched by the volcanic explosion in the Pacific, leveraging multistatic meteor radar observations from the Chilean Observation Network De Meteor Radars (CONDOR) and the Nordic Meteor Radar Cluster in Fennoscandia. MLT winds are computed using a recently developed 3DVAR+DIV algorithm. We found eastward- and westward-traveling gravity waves in the CONDOR zonal and meridional wind measurements, which arrived 12 and 48 h after the eruption, and we found one in the Nordic Meteor Radar Cluster that arrived 27.5 h after the volcanic detonation. We obtained observed phase speeds for the eastward great circle path at both locations of about 250 m s−1, and they were 170–150 m s−1 for the opposite propagation direction. The intrinsic phase speed was estimated to be 200–212 m s−1. Furthermore, we identified a potential lamb wave signature in the MLT winds using 5 min resolved 3DVAR+DIV retrievals. 
    more » « less
  3. Abstract Utilizing multistatic specular meteor radar (MSMR) observations, this study delves into global aspects of wind perturbations in the mesosphere and lower thermosphere (MLT) from the unprecedented 2022 eruption of the Hunga Tonga‐Hunga Ha'apai (HTHH) submarine volcano. The combination of MSMR observations from different viewing angles over South America and Europe, and the decomposition of the horizontal wind in components along and transversal to the HTHH eruption's epicenter direction allow an unambiguous detection and identification of MLT perturbations related to the eruption. The performance of this decomposition is evaluated using Whole Atmosphere Community Climate Model with thermosphere/ionosphere extension (WACCM‐X) simulations of the event. The approach shows that indeed the HTHH eruption signals are clearly identified, and other signals can be easily discarded. The winds in this decomposition display dominant Eastward soliton‐like perturbations observed as far as 25,000 km from HTHH, and propagating at 242 m/s. A weaker perturbation observed only over Europe propagates faster (but slower than 300 m/s) in the Westward direction. These results suggest that we might be observing the so‐called Pekeris mode, also consistent with theL1pseudomode, reproduced by WACCM‐X simulations at MLT altitudes. They also rule out the previous hypothesis connecting the observations in South America to the Tsunami associated with the eruption because these perturbations are observed over Europe as well. Despite the progress, theL0pseudomode in the MLT reproduced by WACCM‐X remains elusive to observations. 
    more » « less
  4. Idealized large-eddy simulations of shallow convection often utilize horizontally periodic computational domains. The development of precipitation in shallow cumulus convection changes the spatial structure of convection and creates large-scale organization. However, the limited periodic domain constrains the horizontal variability of the atmospheric boundary layer. Small computational domains cannot capture the mesoscale boundary layer organization and artificially constrain the horizontal convection structure. The effects of the horizontal domain size on large-eddy simulations of shallow precipitating cumulus convection are investigated using four computational domains, ranging from 40×40km2 to 320×320km2 and fine grid resolution (40 m). The horizontal variability of the boundary layer is captured in computational domains of 160×160km2. Small LES domains (≤40 km) cannot reproduce the mesoscale flow features, which are about 100km long, but the boundary layer mean profiles are similar to those of the larger domains. Turbulent fluxes, temperature and moisture variances, and horizontal length scales are converged with respect to domain size for domains equal to or larger than 160×160km2. Vertical velocity flow statistics, such as variance and spectra, are essentially identical in all domains and show minor dependence on domain size. Characteristic horizontal length scales (i.e., those relating to the mesoscale organization) of horizontal wind components, temperature and moisture reach an equilibrium after about hour 30. 
    more » « less
  5. The development of a multistatic radar system for the Colorado Zephyr Meteor Radar Network is described in this article. This system relies on recent developments in the field of meteor radar, including advancements in software-defined radio-based radar receivers, multistatic wind retrieval, coded constant- wave transmit signals, and transmit-side interferometry. We present the current status of a prototype multistatic radar transmitter deployed in Platteville, Colorado, and a forward look toward how it can inform the design of a large-scale radar network. 
    more » « less