skip to main content

Explore Research Products in the NSF-PAR

  • NSF Public Access
  • Search Results
  • Four-dimensional mesospheric and lower thermospheric wind fields using Gaussian process regression on multistatic specular meteor radar observations
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
NSF-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. ; ; ; ; ; ; ; ( , Earth, Planets and Space)
    Abstract

    In recent years, multistatic specular meteor radars (SMRs) have been introduced to study the Mesosphere and Lower Thermosphere (MLT) dynamics with increasing spatial and temporal resolution. SMRs, compared to other ground-based observations, have the advantage of continuously measuring the region between 80 and 100 km independent of weather, season, or time of day. In this paper, frequency spectra of MLT horizontal winds are explored through observations from a campaign using the SIMONe (Spread-spectrum Interferometric Multistatic meteor radar Observing Network) approach conducted in northern Germany in 2018 (hereafter SIMONe 2018). The 7-day SIMONe 2018 comprised of fourteen multistatic SMR links and allows us to build a substantial database of specular meteor trail events, collecting more than one hundred thousand detections per day within a geographic area of$$\sim$$∼500 km$$\times$$×500 km. We have implemented two methods to obtain the frequency spectra of the horizontal wind components: (1) Mean Wind Estimation (MWE) and (2) Wind field Correlation Function Inversion (WCFI), which utilizes the mean and the covariances of the line of sight velocities, respectively. Monte Carlo simulations of a gravity wave spectral model were implemented to validate and compare both methods. The simulation analyses suggest that the WCFI helps us to capture the energy of smaller scale wind fluctuations than those capture with MWE. Characterization of the spectral slope of the horizontal wind at different MLT altitudes has been conducted on the SIMONe 2018, and it provides evidence that gravity waves with periods smaller than 7 h and greater than 2 h dominate with horizontal structures significantly larger than 500 km. In the future, these analyses can be extended to understand the significance of small-scale fluctuations in the MLT, which were not possible with conventional MWE methods.

    Graphical Abstract

     
    more » « less
  2. ; ; ; ; ; ; ; ( , Earth and Space Science)

    There are few observational techniques for measuring the distribution of kinetic energy within the mesosphere with a wide range of spatial and temporal scales. This study describes a method for estimating the three‐dimensional mesospheric wind field correlation function from specular meteor trail echoes. Each radar echo provides a measurement of a one‐dimensional projection of the wind velocity vector at a randomly sampled point in space and time. The method relies on using pairs of such measurements to estimate the correlation function of the wind with different spatial and temporal lags. The method is demonstrated using a multistatic meteor radar data set that includes ≈105meteor echoes observed during a 24‐hr time period. The new method is found to be in good agreement with the well‐established technique for estimating horizontal mean winds. High‐resolution correlation functions with temporal, horizontal, and vertical lags are also estimated from the data. The temporal correlation function is used to retrieve the kinetic energy spectrum, which includes the semidiurnal mode and a 3‐hr period wave. The horizontal and vertical correlation functions of the wind are then used to derive second‐order structure functions, which are found to be compatible with the Kolmogorov prediction for spectral distribution of kinetic energy in the turbulent inertial range. The presented method can be used to extend the capabilities of specular meteor radars. It is relatively flexible and has a multitude of applications beyond what has been shown in this study.

     
    more » « less
  3. ; ; ; ; ; ; ; ; ( , Earth and Space Science)
    Abstract

    The mesosphere and lower thermosphere (MLT) region is dominated globally by dynamics at various scales: planetary waves, tides, gravity waves, and stratified turbulence. The latter two can coexist and be significant at horizontal scales less than 500 km, scales that are difficult to measure. This study presents a recently deployed multistatic specular meteor radar system, SIMONe Peru, which can be used to observe these scales. The radars are positioned at and around the Jicamarca Radio Observatory, which is located at the magnetic equator. Besides presenting preliminary results of typically reported large‐scale features, like the dominant diurnal tide at low latitudes, we show results on selected days of spatially and temporally resolved winds obtained with two methods based on: (a) estimation of mean wind and their gradients (gradient method), and (b) an inverse theory with Tikhonov regularization (regularized wind field inversion method). The gradient method allows improved MLT vertical velocities and, for the first time, low‐latitude wind field parameters such as horizontal divergence and relative vorticity. The regularized wind field inversion method allows the estimation of spatial structure within the observed area and has the potential to outperform the gradient method, in particular when more detections are available or when fine adaptive tuning of the regularization factor is done. SIMONe Peru adds important information at low latitudes to currently scarce MLT continuous observing capabilities. Results contribute to studies of the MLT dynamics at different scales inherently connected to lower atmospheric forcing and E‐region dynamo related ionospheric variability.

     
    more » « less
  4. ; ; ; ; ; ; ; ; ; ; et al ( , Atmospheric Measurement Techniques)
    Abstract. Ground-based remote sensing of atmospheric parameters is often limited to single station observations by vertical profiles at a certain geographic location. This is a limiting factor for investigating gravity wave dynamics as the spatial information is often missing, e.g., horizontal wavelength, propagation direction or intrinsic frequency. In this study, we present a new retrieval algorithm for multistatic meteor radar networks to obtain tomographic 3-D wind fields within a pre-defined domain area. The algorithm is part of the Agile Software for Gravity wAve Regional Dynamics (ASGARD) and called 3D-Var, and based on the optimal estimation technique and Bayesian statistics. The performance of the 3D-Var retrieval is demonstrated using two meteor radar networks: the Nordic Meteor Radar Cluster and the Chilean Observation Network De Meteor Radars (CONDOR). The optimal estimation implementation provide statistically sound solutions and diagnostics from the averaging kernels and measurement response. We present initial scientific results such as body forces of breaking gravity waves leading to two counter-rotating vortices and horizontal wavelength spectra indicating a transition between the rotational k−3 and divergent k-5/3 mode at scales of 80–120 km. In addition, we performed a keogram analysis over extended periods to reflect the latitudinal and temporal impact of a minor sudden stratospheric warming in December 2019. Finally, we demonstrate the applicability of the 3D-Var algorithm to perform large-scale retrievals to derive meteorological wind maps covering a latitude region from Svalbard, north of the European Arctic mainland, to central Norway. 
    more » « less
  5. ; ; ; ; ; ; ; ; ; ; et al ( , Atmospheric Measurement Techniques)
    Abstract. Meteor radars have become widely used instruments to study atmospheric dynamics, particularly in the 70 to 110 km altitude region. Thesesystems have been proven to provide reliable and continuous measurements of horizontal winds in the mesosphere and lower thermosphere. Recently,there have been many attempts to utilize specular and/or transverse scatter meteor measurements to estimate vertical winds and vertical windvariability. In this study we investigate potential biases in vertical wind estimation that are intrinsic to the meteor radar observation geometryand scattering mechanism, and we introduce a mathematical debiasing process to mitigate them. This process makes use of a spatiotemporal Laplacefilter, which is based on a generalized Tikhonov regularization. Vertical winds obtained from this retrieval algorithm are compared to UA-ICON modeldata. This comparison reveals good agreement in the statistical moments of the vertical velocity distributions. Furthermore, we present the firstobservational indications of a forward scatter wind bias. It appears to be caused by the scattering center's apparent motion along the meteortrajectory when the meteoric plasma column is drifted by the wind. The hypothesis is tested by a radiant mapping of two meteor showers. Finally, weintroduce a new retrieval algorithm providing a physically and mathematically sound solution to derive vertical winds and wind variability frommultistatic meteor radar networks such as the Nordic Meteor Radar Cluster (NORDIC) and the Chilean Observation Network De meteOr Radars(CONDOR). The new retrieval is called 3DVAR+DIV and includes additional diagnostics such as the horizontal divergence and relative vorticity toensure a physically consistent solution for all 3D winds in spatially resolved domains. Based on this new algorithm we obtained vertical velocitiesin the range of w = ± 1–2 m s−1 for most of the analyzed data during 2 years of collection, which is consistent with the values reportedfrom general circulation models (GCMs) for this timescale and spatial resolution. 
    more » « less
  • Citation Formats
  • MLA
  • APA
  • Chicago
  • BibTeX
  • Save / Share this Record
  • Send to Email