Our analysis calculates the (aliased) range, (aliased) Doppler, acceleration, azimuth and elevation of detections. For known RSO’s (resident space objects) these parameters can be compared to those obtained from North American Aerospace Defense Command (NORAD) two line ephemeris (TLE) data, which is forwards propagated in time using a simplified perturbations model. Non-catalogued objects are arguably more important information to Defence, and we will show examples of such detections.
A major aim of this work is to simultaneously detect ST winds and satellites, thus enabling all wind profiling radars to play a role in space domain awareness while maintaining their primary purpose. We will present results of operational detections, in conjunction with ST winds from the BPST. We will discuss the advantages and challenges of multi-purpose radar, and consider the future role these class of radars will play in Australia.
References
[1] Holdsworth, David A., Spargo, Andrew J., Reid, Iain M. and Adami, Chris, (2020a), Low Earth Orbit object observations using the Buckland Park VHF radar, Radio Sci., 55, e2019RS006873, https://doi.org/10.1029/2019RS006873
[2] Heading, E., S. T. Nguyen, D. Holdsworth, I.M. Reid, Micro-Doppler Signature Analysis for Space Domain Awareness using VHF radar, Remote Sens. 2024, 16, 1354. https://doi.org/10.3390/rs16081354
[3] Holdsworth, D.A.; Spargo, A.J.; Reid, I.M.; Adami, C.L. Space Domain Awareness Observations Using the Buckland Park VHF Radar. Remote Sens. 2024, 16, 1252. https://doi.org/10.3390/rs16071252
To delve even further into the study of head echoes, we harnessed the power of a specially trained Sequence-to-One Regression convolutional neural network. This advanced tool allowed us to categorize echoes into six distinct types, each revealing unique signal patterns and behaviors within radar systems. Type I echoes, the most common, exhibit a low signal-to-noise ratio and lack distinct features, constituting approximately 75% of all detections. In contrast, Type II echoes demonstrate sudden power changes within a 0.4-microsecond Inter-Pulse Period, equivalent to 60 km. Type III echoes, considered as 'good' meteors, consistently align with the radar beam pattern and display a smooth bell curve. Types IV and V echoes are characterized by strong, irregular, and regular interference patterns, respectively. Finally, Type VI echoes stand out for their remarkable ability to traverse multiple radar lobes.
This meticulous categorization significantly amplifies the meteor community's capability to analyze meteoroid entry, atmospheric dynamics, and radar scattering properties. It represents a significant leap forward in the field, offering not only enhanced detection methods but also profound insights into atmospheric science. These findings have the potential to shape future research endeavors and practical applications in meteor research and related disciplines.
To better understand the effect of GWs on the background flow, more observations in the altitude range above 80 km are needed. One form of observation is provided by airglow imagers, which detect, for example, the infrared emission from OH at an altitude of ~86 km. GWs modulate the temperature of the atmosphere at this altitude, which in turn determines among other things the intensity of the detected radiation. Airglow imagers such as the Advanced Mesospheric Temperature Mapper (AMTM) developed by Mike Taylor and Dominique Pautet observe GWs over a horizontal range of approximately 200 km x 200 km throughout the night. Instead of linear monochromatic waves, the observations reveal a chaotic superposition of waves with many different properties.
These GW properties such as amplitude, wavelength, and frequency must be measured precisely in order to systematically determine the GW momentum flux. Since GWs occur not only intermittently in localized regions but also can occupy quasi-stationary very large spatial areas and their spectral components are functions of space and time, a mathematical tool is needed that allows for the examination of nonlinear non-stationary signals in multiple dimensions. This tool is provided by the multidimensional continuous wavelet transformation (CWT).
Applying the CWT to a 1D dataset, such as a time series, results in a 2D wavelet spectrum that possesses spectral power as a function of time and frequency. Since the AMTM detects wave signatures in two spatial and one time dimension, the result of the CWT is a 6D space. The question arises of how to most effectively explore this multidimensional wavelet spectrum and how to identify individual wave packets. In my approach, I use hierarchical density-based spatial clustering of applications with noise (HDBSCAN). The clusters identified in spectral power are used to derive the spectral properties of individual wave packets and to reconstruct the wave packets themselves.
The methodology is first applied to artificial datasets where multiple waves with different properties are superimposed. This allows for the determination of the methodology's constraints. Finally, the methodology is applied to real observations from the AMTM in Río Grande, Argentina.
Ultimately, the addition of co-located temperature and wind measurements by a middle atmosphere Rayleigh lidar and a meteor radar in Río Grande, enable a precise determination of GW momentum flux in the mesosphere / lower thermosphere.
A strong variation (~60 m/s) in the seasonal mean diurnal variation is observed, suggesting the effect of tides in the MLT region. Background mean zonal winds show a clear semi-annual oscillation in the mesosphere (MSAO), with a peak amplitude of 20 m/s at 80 km, whereas meridional winds are characterized by an annual oscillation as expected. Weak magnitudes in MSAO are observed during the disturbed stratospheric quasi-biannual oscillation (SQBO) in 2015 and 2020-2021. Interestingly, these years also coincide with the strongest cold and warm phases of ENSO, respectively, suggesting the role of lower atmospheric processes in the middle atmospheric structure and dynamics. Long-term trends in mean winds constructed using Rocket, Satellite (HRDI), MST radar, and Meteor radar data during 1971-2024 in the overlapping regions are also presented and discussed in relation to anthropogenic influences on middle atmospheric structure and dynamics in a changing climate.
To understand the role of neutral dynamics on the October effect seen in very low frequency (VLF) observations, we have studied the interannual variation of tides (particularly semi-diurnal-S2) and PWs amplitude during the summer and fall transition period at the high and middle latitudes. To estimate the tidal and PW amplitudes, we used 1 hour (1 km) temporal (vertical) resolution Saura and Juliusruh partial reflection radar (PRR) zonal and meridional winds at high and mid-latitudes, respectively. To extract the tidal and PW amplitudes, we applied a fast Fourier transform of 42 days window with 1 hour sliding. This analysis has been carried out for the altitude range between 75 and 90 km for Saura and Juliusruh radar winds from the year 2005 to 2023. To investigate the interannual variations, we looked into the maximum and minimum amplitudes during the summer and fall-to-winter transition periods. Furthermore, we also looked into the zonal wind reversal from summer to fall/winter. We found that the VLF October effect is more connected to the reversal in the zonal wind than the S2 amplitude minima both in the high and mid-latitudes. The day of the S2 maximum and VLF October effect coincide during the high and moderate solar activity years, however, during the low solar activity condition the S2 maximum was observed earlier than the VLF October effect. Detailed results and physical reasoning will be presented in this presentation.
Detailed analysis of anomalies from the 7-year climatology unveiled the successive interaction between QTDWs and GWs that plays a pivotal role in the IHC: When IHC occurs associated with a warming in the winter (NH) stratosphere, a positive (eastward) anomaly of GW forcing appears in a height range of z= 60–70 km at ~10° S–30° S. This positive forcing anomaly is due to GW filtering by anomalously strong westward wind in thermal-wind balance with a cold anomaly in the equatorial stratosphere. Positive wave forcing is equivalent to northward flux of potential vorticity (PV). Thus, the positive GW forcing anomaly acts to steepen a negative latitudinal PV gradient in the summer low-latitude mesosphere. In the region of the strong negative PV gradient, anomaly of EP flux divergence associated with QTDWs becomes positive. This fact indicates that the GW-induced negative anomaly of PV gradient enhances QTDW generation through barotropic and/or baroclinic instability. The QTDWs break and exert anomalously strong negative forcing in z= 85–95 km at 40°–70° S. Downwelling induced by the QTDW forcing yields a positive temperature anomaly in z=80–90 km at ~60° S. To satisfy the thermal-wind balance, an eastward anomaly appears in zonal winds in the polar summer upper MLT. Above the wind anomaly, a negative anomaly of GW forcing occurs. The GW forcing seems to be the direct cause of downwelling and thus a warming in the polar summer upper MLT.
We also compared the resolved GW forcing with parameterized GW forcing in reanalysis data produced by a data assimilation system JAGUAR-DAS using a medium-resolution JAGUAR. The GW parameterizations underestimate the GW forcing anomaly enhancing QTDW generation, whose magnitude is ~1/10 to 1/5 the magnitude of the resolved forcing. This result suggests that GWs not satisfying the assumptions in GW parameterizations make a considerable contribution to the enhancement of QTDW activity and the resulting IHC.
We have processed and analyzed new two-dimensional (2D) observations of F-region echoes made by a 14-panel version of the AMISR system installed at the JRO. Semi-routine (~200 nights per year) observations have been made since July 2021. AMISR-14 scans the magnetic equatorial plane and produces “images” of ESF echoes over a zonal distance of +/- 200 km with respect to the JRO. Therefore, the images allow us to determine, unambiguously, whether radar ESF was generated locally or not. Here, we use these images to better understand the development of post-midnight ESF events.
We present and discuss results of the analyses of 2D observations of post-midnight ESF made by AMISR-14 between July 2021 and August 2023. The results include occurrence rates of post-midnight ESF events detected by the system and the seasonal and geomagnetic conditions under which these events developed. We also present our findings about the local versus non-local development of these post-midnight ESF events.
Zhan, W., Rodrigues, F., & Milla, M. (2018). On the genesis of postmidnight equatorial spread F: Results for the American/Peruvian sector. Geophysical Research Letters, 45, 7354–7361. https://doi.org/10.1029/2018GL078822.
Acknowledgment: This work was supported by NSF award AGS-1916055 and by a NDSEG fellowship.
Detection of such echoes are important for global climate-change studies altitude, as their detection highlights regions of cold temperatures at ~90 km altitude. Typically the existence of such scattering entities relates to ice-crystals at temperatures below 140K. The importance of choice of frequency is also discussed, and it is emphasized that frequencies less than 30 MHz offer nothing of value in regard to temperatures within the scattering echoes; only detection at frequencies > 30 MHz are valuable in this regard.
Observational results of this campaign will be presented and discussed.
References
Hocking, W.K. and V.L. Pinnegar, "Mid-Latitude Detection of High Schmidt-Number Turbulent Echoes, and Comparison to PMSE and Geomagnetic Variations". Atmosphere, 13, 396, https://doi.org/10.3390/atmos13030396, 2022.
http://www.mdpi.com/2073-4433/13/3/396.
During a Solar Energetic Particle event in February 2024 resulting in proton precipitation, we observed Polar Mesosphere Winter Echoes (PMWE) from Svalbard for the first time. This is a record high latitude for the phenomena.
This shows the improved capability of the radar after the upgrade, and is promising since we are planning further improvements to the system in the future.
During a 5-minute measurement cycle, the radar measures 6 vertical dwells and four dwells at 6° off-vertical in four orthogonal azimuths. The raw spectral width of each spectrum was obtained by first subtracting the noise, then identifying the strongest peak and the spectral region around that peak where the signal exceeded the noise, and finally deriving the Gaussian-equivalent width σw from the first and second moments of the identified signal region. Individual dwells from all four 6° beams were averaged over six measurement cycles (33 minutes) to produce a raw spectral width measurement. These raw values were corrected for beam and shear broadening using the methods described by Hocking et al (2016) to give estimates of the broadening due to turbulence. Because the spectral resolution of the radar is 0.154 ms-1, values of corrected spectral width need to be considerably larger than this to be practically useful; for this study a threshold of 0.4 ms-1 was adopted. By the formula presented by Luce et al (2018), based on the MU radar, the turbulent dissipation rate ε ~ σw3/L where L~ 60 m, so 0.4 ms-1 corresponds to ε~10-3 m2s-3, which is strong turbulence (e.g. Ko et al (2019) quote the mean ε for turbulent layers to be 1.84x10-4 and 1.37x10-4 m2s-3 in the troposphere and stratosphere respectively).
A preliminary study was conducted for September-October 2016 to establish a method for analysing the turbulence. For each day, image-analysis software (available in IDL) was used to identify and extract turbulent patches meeting the threshold criterion. Erosion and dilatation operations, together with a minimum size (in terms of time-height pixels) was used to retain only the most significant patches. In the resulting set of 221 patches, 217 showed a mean aspect sensitivity (vertical – off-vertical signal power) < 5 dB, consistent with near-isotropic scattering as expected from turbulence. The influence of vertical wind shear on the strength of turbulence was clear from the positive correlation (r=0.62) between the maximum spectral width and maximum shear in the patch. However, mountain waves and convection were also found to be influential in generating turbulence, especially in the troposphere.
References:
W. K. Hocking et al., 2016. ISBN 978-1-107-14746-1
H.-C. Ko et al., 2019. doi: https://doi.org/10.1029/2019JD030287.
H. Luce et al, 2018. doi: 10.1186/s40623-018-0979-1
With this overview, over the past decades, MST/ST radars operating at VHF frequencies have emerged as ideal tools for studying these waves due to their optimal vertical and temporal resolution. There have been numerous observational case studies on convectively generated gravity waves over the tropical Indian region using MST radars at Gadanki (13.5oN, 79.2oE) and at Indonesia (0.2oS, 100.3oE). However, no detailed investigations have been done over the subtropical Himalayan region, where the topography and the convective system can simultaneously affect the generation of gravity waves.
In this context, we have utilised the 206.5 MHz ST Radar installed at Nainital (29.4oN, 79.5oE;1793m amsl) in the Himalayan foothills to examine the gravity waves generation during deep convective systems induced by western disturbance as well as south-west monsoon. A preliminary analysis during deep convective storms revealed upward propagating gravity waves, usually originating from mid-troposphere travelling into the lower stratosphere region. A combination of spectral and hodograph analysis of zonal, meridional and vertical wind perturbations have been utilised to estimate wave parameters. The interaction of the wave with the highly variable wind leads to the distribution of the intrinsic periods varying from ~34 min – 1100 min and vertical wavelength of ~5.6 – 8.4 km. Furthermore, momentum flux has been estimated to be about 10 times with respect to the pre-convective period and reversing its direction in the Upper Troposphere and the Lower Stratosphere (UTLS) region.
This study will offer a comprehensive outlook on the parameterisation of convective gravity waves produced by western disturbances and the monsoon over the Himalayan site. Additionally, it will delve into the mechanisms behind their generation, propagation characteristics, and the budget of momentum flux in the UTLS region. Details will be presented at the upcoming conference.
In this work, we use multi-year observations of winds in the lower and middle atmosphere from the Middle Atmosphere Alomar Radar System (MAARSY) located in Norway (69.30°N and 16.04°E) to infer temporal high-order SFs, which are systematically converted to spatial SFs employing the Taylor hypothesis. We describe the characteristics of second- and third-order structure functions in detail. The limitations of the method, as well as comparisons with measurements from other sources, are also discussed.
We present a detailed analysis of the probabilities of convective and dynamic instabilities in the 85-100 km altitude region based on over 2000 hours of high-quality sodium lidar temperature and horizontal wind measurements made at the Andes Lidar Observatory (ALO) in Cerro Pachón, Chile (30.25S, 70.73W). An important aspect of the analysis is the removal of probability bias caused by measurement errors. On average, probabilities vary from 1-4% for convective and 6-9% for dynamic instabilities in the 85-100 km altitude region. Short-period (<1 hour) waves contribute most to the occurrence of instabilities but the combination of waves at different scales is the most important process that creates instabilities. It is also found that turbulence heat flux is much stronger in unstable regions than in stable regions.
In this work, we have applied velocity filtering techniques to both multistatic specular meteor radar measurements and global circulation model simulations to analyze horizontal wind frequency spectra over southern Patagonia. We decompose the frequency spectra into divergent and rotational parts and find that their energy contributions are equipartitioned at high frequencies. This mesoscale energy equipartition is a feature of stratified turbulence and illustrates the complexity of the mesoscale dynamics in the summer mesopause region.
We also analyze the power spectra of observed and simulated mesoscale zonal and meridional winds at middle latitudes in the Southern Hemisphere and show that stratified turbulence (in the statistical sense) plays a role in the mesopause region during summer.
We showcase the findings from our investigation on the shape, movement, and microphysical characteristics of clouds using lidar, radar, and camera data. To provide insight into their relationship to the background environment, observations will be analysed on time periods ranging from seconds to decades and horizontal scales ranging from a few meters to hundreds of kilometres. The knowledge gained from NLC research is used for studying other phenomena from the lower Stratosphere to the Edge of Space.
VHF radars have become indispensable for providing vertical profiles of three-dimensional winds from the surface to the lower stratosphere with high height (~75 m) and temporal resolution (~1 min). In this context, Stratosphere-Troposphere (ST) radar located at Nainital (29.4oN, 79.5oE, 1793 m AMSL) in the central Himalayan foothills has been used to estimate the momentum flux during MCS and in normal conditions. We have utilized the modified expression of the coplanar beam method pioneered by Vincent and Reid (1983) to estimate the vertical structure of the horizontal momentum flux components Zonal, u'w' and Meridional, v'w'.
u'w'=(vE'-vw')/2sinθ [(vE'+vW')+(vN'+vS')]/4cosθ
v'w'=(vN'-vS')/2sinθ [(vE'+vW')+(vN'+vS')]/4cosθ
where u'w' and v'w' are the vertical flux of zonal and meridional components of momentum flux, respectively. vE', vW', vN', vS' are the radial perturbation velocities in the east, west, north and south beams, respectively. The MCS are identified utilizing X-band weather radar and the outgoing longwave radiation obtained from the geostationary satellite. Therefore, a statistical description of the momentum flux estimate will be presented for the MCS. These results are expected to serve as observational evidence for testing the numerical models and their simulation of tracer transport in the UTLS region.
Reference
Vincent, R.A. Vincent, I.M. Reid (1983), HF doppler measurements of mesospheric gravity wave momentum fluxes. J. Atmos. Sci., 40 (5), 1321-1333
Ground based techniques can be used to measure zonal winds, such as the ISR technique which is usually performed to estimate zonal wind profiles from measured ion drifts at mid- and high-latitude but cannot be used at low-latitudes due to clutter from coherent echoes at E-layer (EEJ). Additionally, the Spread Spectrum Interferometric Multistatic meteor radar Observing Network (SIMONe) system, deployed in September 2019, can estimate neutral winds from specular meteor trail echoes in the mesosphere and lower thermosphere (MLT) region (70 - 110 km) which includes the EEJ layer, however the time resolution is 1 hour and it probes a larger volume than our oblique radar, then estimated winds will be a result of a contribution from different regions. Space-based missions have also been used to measure these winds such as the Michelson Interferometer for Global High-Resolution Thermospheric Imaging (MIGHTI) on board the ICON satellite, the Wind Imaging Interferometer (WINDII) on the Upper Atmosphere Research Satellite (UARS) among others, but their measurements are not continuous over a specific latitude and longitude as they orbit around the planet. On the other hand, a technique to estimate zonal neutral winds from oblique EEJ type II Doppler shifts, sampled every minute approximately, was proposed by Shume et al. [2005]. This method predicts Doppler shifts from neutral winds, compares the Doppler with the measured ones and updates the winds until the best data-model agreement of the Doppler shift (RMSE < 2.0). Although wind estimations are limited to the EEJ echoes detection and SNR intensity, this method provides high time resolution wind profiles from 97 to 107 km approximately and might be complemented by other techniques. This work presents the first results of neutral wind estimations over IGP-JRO at the EEJ region with a resolution of a few minutes and 700 meters over 10 years of oblique spectra data fitted by a skewed Gaussian distribution.
In this study, we investigate the characteristics of interannual oscillations using the meteor radar network winds observed over Juliusruh and Collm from 2004 to 2021. Furthermore, we also utilized the WACCM-X winds and temperatures for 40-60°N and 0-20°E latitude and longitude grids. To identify the interannual oscillations we applied the Lomb-Scargle periodogram (LSP) analysis on the monthly and annual mean observation and model simulation data. We found the 2.4 and 3.2-year oscillations in both annual mean zonal and meridional winds, 4.2 and 7-9-year oscillations in the zonal winds, and 5.5-year oscillations in the annual mean meridional winds. To explore the plausible causative mechanism, we segregate the data into four seasons namely, winter (November- February mean), summer (May-August mean), spring (March- April), and fall (September-October), and carried out the LSP analysis. The obtained results and the potential physical mechanism that causes the interannual oscillation in the mid-latitude MLT region will be detailed in this presentation.
Buhler et al. (2014) have presented one method that is suited to instruments mounted on a sensor moving through the fluid, but it is not well suited to radar applications. In this poster, we present another method which is based on gravity-waves spectra deduced by Fritts and VanZandt, (1987). In particular, the slope of the spectrum at wave periods of a few hours is particularly sensitive to the nature of the underlying cause of the oscillations.
Studies have been made at 10 sites in Canada, situated at sites from the poles to mid-latitudes, over a 10-year period of observation. In this poster, results are presented which show that the relative contributions are seasonally- and site-dependent. Further results can be found in Hocking et al., (2021).
References
Bühler, O., J. Calliesand R. Ferrari,(2014), "Wave-vortex decomposition of one-dimensional ship-track data", J. Fluid Mech, vol. 756,1007-1026.
Fritts, D.C. and T. E. VanZandt, (1987), "Effects of Doppler Shifting on the Frequency Spectra of Atmospheric Gravity Waves", J. Geophys. Res., 92(D8), 9723-9732.
Hocking, W.K., S. Dempsey, M. C. Wright, P.A. Taylor and F. Fabry, "Studies of Relative Contributions of Internal Gravity Waves and 2-D Turbulence to Tropospheric and Lower Stratospheric Temporal Wind Spectra measured by a Network of VHF Windprofiler Radars using a Decade-long Data-set in Canada", Q. J. R. Meteorol Soc., Vol. 147(740), 3735-3758, https://rmets.onlinelibrary.wiley.com/doi/10.1002/qj.4152, 2021.
In this study, we investigate the relationship between CAT and gravity waves, with a specific focus on tracking the formation of these waves around regions of inertial instability. Previously, [4] showed the emission of inertia–gravity waves following the release of inertial instability using idealised model simulations. Here, we use the WRF model to consider some real-world examples of where regions of low potential vorticity (PV) in the vicinity of the jet stream are associated with inertia–gravity waves. We track the waves as they propagate and investigate whether the causal link found by Thompson and Schultz can be observed in more realistic simulations.
We present results from several case studies exhibiting this behaviour, identifying the sources of the gravity waves observed in simulations. The characteristics of these waves will be compared to those in the idealised model simulations, and gravity-wave parameters will be calculated. Finally, we widen our analysis by examining the broader upstream pattern that contributes to the development of the initial inertial instabilities and explore the different regimes under which these phenomena occur.
References:
[1] Gultepe, I. et al. (2019), "A review of high impact weather for aviation meteorology." Pure and Applied Geophysics, 176, pp.1869–1921.
[2] Williams, J. K. (2014), "Using random forests to diagnose aviation turbulence. " Machine Learning, 95, pp.51-70.
[3] Meneguz, E., Wells, H. and Turp, D. (2016), "An automated system to quantify aircraft encounters with convectively induced turbulence over Europe and the Northeast Atlantic." Journal of Applied Meteorology and Climatology, 55(5), pp.1077–1089.
[4] Thompson, C. F. and Schultz, D. M. (2021), "The release of inertial instability near an idealized zonal jet." Geophysical Research Letters, 48(14), e2021GL092649.
Two different mesospheric phenomena are most frequently observed, which are the summer echoes and the so-called winter echoes (MSE and MWE, and PMSE, PMWE for polar locations, respectively). While the summer echoes are related to reduced electron diffusivity due the presence of charged ice particles around 84km altitude, the generally much fainter winter echoes typically occur between 50 and 85km height.
Specifically in the last two decades PMWE gained interest and case studies as well as statistical analysis of their occurrence have been published.
Winter echoes are most frequently observed around the equinoxes and show a clear diurnal pattern, which is related to the available electron density controlled by solar and geomagnetic activity.
Given the remaining open questions on the cause, their formation and propagation dedicated sounding rocket campaigns have been performed in the last few years.The dependence of the PMWE visibility with the VHF radar MAARSY on the ambient electron density has been investigated in depth with the co-located partial reflection radar Saura for the solar minimum conditions.
Besides a minimum required electron density neutral dynamics are assumed to play a crucial role for the formation, i.e. creation of turbulent structures of suitable scales and their transport. In this study we will analyze radar observations of PMWE in northern Norway using MAARSY multibeam experiments and connect them to background dynamics using the Saura partial reflection radar as well as ECMWF renanalysis data.
As part of this unattended mode, low range resolution mesospheric echoes were detected, and we were able to estimate the zonal and meridional winds at the detected ranges.
In the present study, we will present the preliminary results of the climatology of the winds estimates obtained from the mesospheric winds and compared with other systems installed at the same location.
Several non-specular meteors studies have been conducted with the high-power large-aperture (HPLA) radar at the Instituto Geofísico del Perú's Jicamarca Radio Observatory. Studies of interest are the ones that use the high-power meteor trails to estimate Mesosphere and Lower Thermosphere (MLT) winds (80 - 120 km) by combining 3 or more received signals using interferometry. The main limitation with this technique is the strong presence of the Equatorial ElectroJet (EEJ) between 90 - 110 km.
A novel method to extract non-specular meteors trails based on their shape through a computer vision model called RetinaNet in presence of the EEJ is presented in this work. The obtained MLT winds obtained with the HPLA is compared with the routinely winds obtained with the specular meteor radar SIMONe collocated with the Jicamarca radar.
Hourly values of the peak electron density of the ionospheric F2-layer (NmF2) from the ionosonde at Juliusruh station (54.6°N, 13.4°E) are analyzed. The analysis period is between 1957 and 2023. Geomagnetic perturbations are removed based on kp index with a threshold of 3. Linear regression and 3rd order polynomial fit models dependent on different solar EUV flux indices (MgII, F30, and F10.7) are generated for each solar cycle and for winter (January), summer (July), and two equinox months (March and October) separately.
The coefficient R^2 is used to measure the quality of the NmF2 model dependent on the solar EUV flux proxies. It shows a well-pronounced local time dependence in January with R^2 value being maximum around noon hours. During the other months, R^2 is rather constant throughout the day. R^2 is highest for F30 and MgII in general and especially during winter noon conditions. Comparing the slopes of the model curves for the different solar cycles, a clear decrease in the slope with each solar cycle becomes visible. We conclude that the best solar EUV flux proxies for describing the NmF2 dependence at all LT hours are MgII and F30. In addition, the response of NmF2 to solar EUV flux shows a clear long-term change as the model slope decreases with time for each solar cycle.
The tropopause, a critical atmospheric boundary, separates the dynamically active troposphere from the more stable stratosphere. Its characteristics influence atmospheric processes, including stratosphere-troposphere exchange (STE). The tropopause plays a significant role in regulating the transport of water vapor into the stratosphere, thus impacting Earth's radiative balance and climate. Accurately determining tropopause height and its variability is crucial for climate research. Radiosondes remain a primary tool for tropopause measurements, but it faces face inherent limitations in providing sufficient spatial and temporal information. Very High Frequency (VHF) radar offers an alternative. The radar technique involves assessing the rate of change of Signal-to-Noise Ratio (SNR) values with altitude. As the tropopause is characterized by distinct shifts in SNR, detecting the gradient helps identify the boundary between the troposphere and the stratosphere. However, the SNR gradient method can be less accurate in the presence of multiple tropopauses or strong convective activity. We propose a novel approach by using the Change Point Detection (CPD) method for the first time to estimate tropopause height from an SNR profile. Our novel approach employs one of the CPD approaches, Jenks natural breaks classification to detect the tropopause based on the SNR profile. This data-driven method pinpoints abrupt SNR shifts, demonstrating adaptability under various atmospheric conditions. Our technique shows a strong correlation (0.89) and low root mean square error (RMSE) of 0.29 km compared to radiosonde-derived tropopause heights. This presents a promising new technique for tropopause identification using radar. The paper also details other CPD approaches we tested.
Using the Equatorial Atmosphere Radar (EAR) in Kototabang, west Sumatra, we have applied a method for estimating the specific humidity, q, which was originally developed with the MU radar. We employed a relation between the refractive index gradient squared, M2, and the turbulence echo power (the volume reflectivity). The key technique here is the Radio Acoustic Sounding System (RASS) for simultaneously monitoring profiles of the temperature, T, as well as the Brunt–Väisälä frequency squared N2. This method, however, determined only the magnitude of |M|, therefore, delineation of its sign is the critical problem. The other uncertainty is an accurate estimation of the filling factor of turbulence within a radar range volume. Simultaneous q profile derived from radiosonde data was used as an initial condition. From the continuous measurements of the precipitable water vapor (PWV) with a GNSS receiver, we constrained the integrated amount of water vapor so as to delineate a better estimate of the sign of |M|.
The other unique method of the q estimate has been recently developed by using the two-cascaded steps of several machine learning (ML) models, which processed the basic moment values of the Doppler spectra, such as the power (first moment), wind velocity (second moment), and the spectrum width (third moment), from WPR. Radiosonde data is used as a reference.
This study combines the above two methods to estimate a q profile by applying various ML algorithms to the basic parameters observed with EAR-RASS measurements, such as T , echo power and spectral width . These parameters have a relation with turbulence energy dissipation rate, volume reflectivity of the turbulence echo, and the volume filling factor of turbulence layers, which are essential in calculating specific humidity .
We processed the EAR-RASS data collected between 28 August until 3 September 2016. First, we removed clutter echoes and outliers, and interpolation was done when necessary, then, the data is converted into an appropriate scale. As a reference and label for ML inputs, we use q from the ERA5 global reanalysis and radiosonde data launched at the EAR observatory. We examined performance of several ML models, such as Random Forest, Gradient Boosting, XGBoost, and Artificial Neural Network (ANN). Exploring ML models for the atmospheric parameter prediction has been becoming very active. Our study shows that the ML estimation of q by involving RASS temperature data improves reliability and performance compared to the earlier methods, while maintaining comparable accuracy.
The Zephyr network is similar to other multiple-in, multiple-out (MIMO) continuous coded wave radar networks, but with additional fully polarimetric capability to enable the collection of a complete scattering matrix at each timestep of observation. The network currently consists of a 6-channel transmitter, 3 deployed receivers, and 3 more receivers under deployment. With a center frequency of 31.25 MHz, each transmitter element radiates a unique pseudorandom code at 500 W. The network is scoped to grow to 3 transmit sites and at least 24 receivers scattered throughout the front range of Colorado and beyond.
While efforts over the last five years have focused on hardware development and preliminary science results, a recently submitted NSF DASI Track II proposal outlines plans to expand the network. Expansion includes deploying additional transmitters and receivers, creating scalable digital infrastructure for wider data assimilation and easier public access, and developing retrieval techniques that leverage additional information inherent to a large MIMO network. Additional information includes assimilation of observations across many links, higher observation densities and counts, and utilizing the increased geometric diversity. Additional nodes covering the same volume allows for the specular observation of meteors with trail orientations that would otherwise be missed, and the resulting increase in measurement density would greatly improve wind field estimations. This poster details the accomplishments, current state, and future plans of the Zephyr meteor radar network in Colorado.
During the usage of these instruments we are also able to observe a lot more different types of echoes besides the specular-meteor echoes that are currently not used on a routine bases but are interesting targets for further studies.
Some studies on such echoes were already done in a single event based analysis (e.g. Huyghebaert et al. 2023 https://doi.org/10.1029/2023ja031682).
These types of echoes (mainly non-specular meteors, E-region, Aurora, PMSE and EEJ echoes) typically last much longer than the specular-meteor echoes, which typically last far less than a second.
As these echoes typically have a relatively small spectral-width, it is more efficient to detect and store them in the spectral-domain.
We are proposing a method that allows to efficiently detect this echoes and store them much more space efficient than just saving the raw data making a further analysis and saving the data on a regular basis much easier.
This method might also be applicable to other radars observing similar echoes.
In South America, the most notorious signal was 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. The wave was characterized as traveling at ∼240 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 an L1 Lamb wave pseudomode and the so-called Pekeris mode.
The extent of such a signal is investigated using the European systems. To do this, the signal characteristics are studied separately as it propagates westward and eastward. The winds display dominant eastward perturbations reaching distances of up to 25,000 km, with high altitudinal coherence, consistent with the South American signal. Additionally, weaker westward perturbations, propagating faster, are identified. Throughout the study, the longitudinal and transverse horizontal wind decomposition approach is systematically used. This approach shows benefits over the traditional zonal and meridional decomposition for studying wind perturbations originating from localized sources, such as extreme geophysical events.
Here, we study the impact of four SSW/ES events on the VLF signal amplitude between the high latitude Tx-Rx link NRK-NyÅlesund, to gain further knowledge about interactions between the ionospheric D-region and the atmosphere during these significant atmospheric phenomena.
For three of four SSW/ES events occurring during the winter periods of 2009, 2010, and 2018, a very similar VLF signal amplitude variation is observed, which is characterized by a significant increase during the wind reversal in the signal amplitude followed by a decrease during the ES. This study aims to reveal a possible mechanism driving these similar VLF signal amplitude variations, involving an unexpectedly strong correlation with the mesospheric water vapor concentration. Furthermore, causes for the different VLF signal amplitude variations for the SSW/ES event in winter 2019 are discussed.
In this presentation, we give a retrospective overview of the Aeolus mission. We introduce the satellite’s measurement geometry and wind detection method, explore the mission’s primary aims and objectives, and demonstrate the potential of Aeolus and its successors for future research within atmospheric dynamics. Notably, we present the first satellite observations of atmospheric GWs using a Doppler Wind Lidar (DWL). We also display Aeolus’ ability to observe the QBO, including the 2019/2020 QBO disruption, and show the benefit of the high vertical resolution of its observations when compared to atmospheric reanalysis. Finally, we consider a range of additional novel applications of Aeolus’ measurement data, and look ahead to its follow-on mission, EPS-Aeolus, due to be launched by EUMETSAT in 2031.
Here we present Rayleigh-Mie wind lidar, measuring Doppler shift between the emitted and backscattered light by means of a Fabry-Perot interferometer. The wind lidars are deployed at Observatory of Haute-Provence (OHP, Southern France) and at Maido observatory in Reunion Island (tropical Indian Ocean). The instruments are capable of measuring horizontal wind components from the middle troposphere up to 75 km altitude with vertical resolution as high as 100 m and temporal resolution down to 5 minutes. Both lidars have been operated quasi-continuously since 2019 thus providing the first multi-year vertically-resolved time series of wind velocity in the northern midlatitudes and southern tropics.
In this study, we present and discuss various scientific applications of the wind lidar observations. We report examples of a weekly and an hourly observation series, reflecting various dynamical events in the middle atmosphere, such as a Sudden Stratospheric Warming events and occurrence of orographic and non-orographic gravity waves. Wind velocity profiles are analysed in conjunction with temperature lidar profiling, radiosoundings, satellite observations as well as meteorological analysis/reanalysis in order to characterize gravity waves and their vertical propagation and filtering.
Finally, we present a summary of results of the ESA Aeolus space-borne wind lidar validation using its ground-based predecessor in the frame of a series of dedicated campaigns at both observatories.
The European Space Agency's Aeolus satellite, equipped with the Atmospheric LAser Doppler INstrument (ALADIN), provides comprehensive global wind profiling, offering valuable insights into convection-induced gravity waves (GWs) in the upper troposphere and lower stratosphere (UTLS). This study analyzes Aeolus-derived wind data alongside ECMWF ERA5 reanalysis, revealing a migrating hotspot of enhanced GW activity. This activity transitions from the Indian Ocean during the Boreal Summer to the Maritime Continent in Boreal Winter, correlating strongly with convective activity as evidenced by variability in outgoing longwave radiation. The results underscore the importance of Aeolus in improving the representation of convective-induced GWs in ERA5, particularly in the under-sampled Indian Ocean region, thus enhancing our understanding of GW phenomena in the UTLS.
While the benefits of Aeolus global wind profiling are straightforward, important limitations exist, particularly when studying extratropical UTLS dynamics. We identify and characterize these limitations in our study, focusing on the Range Bin Setting (RBS) of Aeolus. This RBS adjusts dynamically based on seasonal and geographical variations, affecting our ability to measure and interpret gravity wave dynamics accurately.
We also examine oscillatory perturbations (OPs) identified in Aeolus data, particularly between 12 and 19 km altitudes, displaying a unique oscillatory signature. These OPs emerge prominently above high-albedo regions such as Greenland and Antarctica. The nature of these oscillations remains unexplained although the dark current anomalies on the instruments (hot pixels) may amplify these oscillations without explaining their inherent characteristics. The absence of these phenomena in GNSS Radio Occultation data suggests that the observed perturbations in Aeolus data are not related to gravity waves but are likely linked to the satellite's specific observational characteristics.
These variations introduce substantial uncertainty into our models and highlight the necessity for developing more precise metrics that can reliably detect and interpret such phenomena.
Our results highlight the system's effectiveness in detecting small-scale wind asymmetries, emphasizing its importance for atmospheric science, particularly in measuring vertical wind velocities. The introduction of the Multi FOV upgrade, featuring an innovative optical arrangement that allows for rapid switching between different viewing angles, represents a major enhancement in our observational capabilities and allows for an in-depth examination of atmospheric dynamics across various scales.
Furthermore, a comparison to winds from ECMWF indicates a significant underestimation of vertical winds by ECMWF, underscoring the VAHCOLI system's accuracy and reliability in capturing complex three-dimensional wind patterns. A validation against Aeolus data, demonstrates a minimal bias. This underlines the potential of ground-based lidar systems to complement satellite observations and advance our understanding of atmospheric phenomena. Our work lays the groundwork for future atmospheric research, providing a valuable tool for validating successors to Aeolus and other spaceborne lidar systems.
In this presentation, we present a novel new way to determine layer tilts using VHF windprofiler radar. The method not only determines layer tilts from horizontal, but also reveals the azimuthal alignment of the layers. Compensation for possible "tilts" in the nominally vertical beam of the radar are also considered.
The angular tilts (both in zenith and azimuth) are determined by cross-correlating vertical winds with radial velocities deduced using the nominally vertical radar beam. The correlations can be ascribed to either (i) tilts in the scatterers within the layers, (ii) tilts in the layers, and/or (iii) geophysical correlations between true vertical velocities and horizontal winds, such as in gravtiy waves. These different interpretations are discussed. Results show typical seasonal and annual variations in tilts, and emphasize the dependence on geographical location.
References
Hocking, W.K., J. Röttger, R.D. Palmer, T. Sato and P.B. Chilson, " Atmospheric Radar: Application and Science of MST Radars in
the Earth's Mesosphere, Stratosphere, Troposphere, and weakly ionized regions", Cambridge University Press, 2016. ISBN
9781316556115, DOI: https://doi.org/10.1017/9781316556115
We have used non-standard pulse codes with a coherent meteor radar network to separate the detections from 2 monostatic GPS locked radar transmitters operating at the same frequency and with the same parameters, and we show representative
results from this meteor network. The technique can also be applied to provide additional forward scatter detection paths within the network, and the future extension of the technique to more than 2 monostatic radars is discussed.
In 1997, Gerrard et al. proposed a novel lidar system targeting metastable helium, which exists in measurable quantities in a layer extending from 250 km to 1000 km or above. Recent advances in detector technology have since made this proposed system feasible, and, starting in 2021, a prototype experiment has been under development at the German Aerospace Center in Oberpfaffenhofen. In early 2022, the first profiles of metastable helium were measured, reaching an altitude of 700 km. Subsequent measurements, with increased system performance, suggest that the helium layer extends well beyond these initial results. Most recently, we have expanded the measurement capabilities to allow for spectroscopic measurements of the helium spectrum.
We will present an overview of the instrumental setup, report on the recent measurements in more detail, and discuss the future potential of the instrument to measure thermospheric wind and temperature profiles.
Lidars are capable to monitor the atmosphere from the troposphere to the thermosphere. In general, such instruments focus onto a limited altitude range for a certain scientific task, such as wind, temperature, metals, ions or aerosols. Moreover, current systems focus on a single method, such as Rayleigh, Mie or resonance scattering applying different technologies. In case of resonance scattering or molecule scattering, such as water vapor, the required wavelength is only addressable by a specially build laser. For a future network of capable lidars, a single technology should cover all applications, regardless of altitude, scattering or scientific goal as far as possible. Since many systems must be build and maintained over long periods of time standardization of the whole chain from the hardware over data analysis to long-term data archives is required. Since multiple systems are required for a network, fast assembling, cost efficiency, and robustness are major goals. Each system must be installed world-wide, operated and repaired. Therefore, each system should be compact, transporatble and self-contained operating in automatic mode after installation.
General-purpose ground-based spectral Doppler-lidars are nowadays capable to obtain multiple parameters, such as wind, temperature, aerosols and metals throughout the whole atmosphere out of compact and robust systems with a size of ~1 m3 in automatic 24/7 operation. The project VAHCOLI (Vertical And Horizontal COverage by Lidar) has shown that a networks of such lidars nowadays is capable of applying a bunch of new measurements methods at different altitudes not even possible with existing lidars specialized for a certain application. VAHCOLI is a software-defined, general purpose Doppler-lidar applying advanced general sub-MHz spectroscopy beyond the capability of existing lidars.
VAHCOLI bases on a tunable alexandrite ring laser with a line width of 3.3 MHz developed in close cooperation of IAP and ILT. Within the EU-project EULIAA the technology is currently extended to the 386 nm line of Fe for a future UV lidar network including the demonstration of measurements at hard-to reach environments and real-time data processing for data bases such as Copernicus. Due to the deep Fraunhofer-line daylight observations at higher altitudes can be achieved for Rayleigh scattering as has been demonstrated in the past by the IAP for Doppler resonance measurements at this wavelength in the mesosphere. With intra-cavity single harmonic generation of the alexandrite laser the required UV-laser does not need amplification stages. The whole system reduces further in size, electric power and allows eye safe operation from the troposphere to the thermosphere. Whereas EULIAA will be the next generation system, IAP and ILT currently transfer the technology of VAHCOLI already within the project LidarCUBE to the industry. Both activities focus currently onto future lidar networks for the lower atmosphere (Doppler Rayleigh & Doppler Mie) as a demonstration of the feasibility of such networks with nowadays technology.
The fully steerable, tri-static, phased-array incoherent scatter radar is located in Skibotn (inland from Tromsø, Norway), Karesuvanto (Finland, north of Kiruna), and Kaiseniemi (Sweden, west of Kiruna). The transmit-receive array at Skibotn consists of about 10,000 aerials and ten 91-aerial outrigger receivers in the immediate vicinity. The receive-only arrays of Kaiseniemi and Karesuvanto consist of about 5,000 aerials each.
Construction of the facility began after the project kick-off in September 2017. First test measurements have been carried out with a single "antenna unit", which consists of 91 crossed-dipole aerials. The next step will be the so-called "PET-7" system, which consists of seven antenna units with full transmit and receive capabilities. PET-7, with 637 aerials with 1 kW power each, will allow for mono-static incoherent scatter radar measurements. Thereafter, EISCAT_3D will gradually expand to fully tri-static operations.
EISCAT_3D will eventually replace the EISCAT mainland radars, i.e. the 930-MHz UHF radar and the 224-MHz VHF radar at Tromsø, even though some overlap is foreseen not least for cross-calibration and comparison purposes. The remote receivers at Sodankylä (Finland) and Kiruna (Sweden) are already inoperable due to failures not feasible to repair.
The EISCAT Svalbard Radar (ESR) and the Ionospheric Heating facility at Tromsø will not be affected by EISCAT_3D becoming operational.
Here we give an overview of the current status and the road ahead.
This work presents the development of two spaced receiving stations, which in combination with the high power-large aperture 50 MHz radar of IGP-ROJ constitutes the J-ARGUS multi-static system, which will provide wavenumber diversity to radar experiments and allow for better estimates of vector Doppler velocities and also enhance the research areas at the ionospheric equator. The first receiving station will be located at the IGP’s Huancayo Observatory, about 170 km east of IGP-ROJ, and the second station will be located at Santa Maria campus of Pontificia Universidad Católica del Perú about 50 km south of IGP-ROJ. Each station will operate over a wide band (~10-80 MHz) and will allow for perpendicular-to-B and oblique radar observations.
Our system is called atc2met (Air-Traffic-Control radio communication to METeorological data) and comprises small and cheap Raspberry Pi-based receivers and a data processing server connected via the internet. Currently, we have three receiver sites in Japan and three in Indonesia. The operation started in 2021, and data sharing began in 2023.
In this presentation, we briefly explain the methodologies and structures of the developed system and show results, public data products, and plans for future development.
This technique is being applied to MF radars at the northern hemisphere, Saura (69N, 16E) and Juliusruh (54N, 13E), routinely operated by Leibniz Institute of Atmospheric Physics. These radars are well equipped with an interferometer capability with 9 and 6 receiver channels, respectively. We have found that existing archived data of these radar systems can be applicable to the meteor echo analyses as those we have done with the Syowa system. The same techinique is now further being planned to apply to the MF system operated by Arctic University of Norway at Tromso (70N, 19E), by adding digital receiver and antenna system to enable all-sky interferometer measurements. The Saura and Tromso systems are ideally separated to share a common illuminating volume, and are expected to resolve fine time and spatial structures of wind fields, even with much better resolutions than those of recent VHF meteor radar based network measurements.
In this work, we present the design of two antenna arrays that will be part of a frequency-agile multistatic radio system for geo-space imaging. Each of these arrays, comprising 256 antennas, will be installed in two cities in Peru (Lima and Huancayo). Alongside the Jicamarca Radio Observatory (JRO) radar, these arrays will serve as reception stations for a multistatic radar system to be used for ionospheric studies.
The objective is to identify the optimal configuration for the elements of both arrays. Firstly, a mathematical model for the antenna arrays is introduced, considering the requirements for scientific objectives as well as any constraints. Subsequently, various techniques for sparse array design are explored, primarily relying on deterministic methods (such as Kogan), stochastic approaches (Simulated Annealing), and convex optimization algorithms. The aim is to compare their solutions with the desired performance of the array. The analysis focuses on obtaining a globally optimal solution, surpassing local optima, thus evaluating different cost functions and initial conditions. For instance, an average peak-to-side-lobe level (PSLL) of -28 dB was achieved with cosine-type elements (inverted V antennas) for a 100-meter diameter available area. Finally, since the arrays are expected to operate within the 10-88 MHz range for radio astronomy purposes, the stability of the solution was verified for specific frequencies and steering angles.