MU Radar Research Articles

Turbulence kinetic energy dissipation rate: assessment of radar models from comparisons between 1.3 GHz wind profiler radar (WPR) and DataHawk UAV measurements

Abstract. The WPR-LQ-7 is a UHF (1.3575 GHz) wind profiler radar used for routine measurements of the lower troposphere at Shigaraki MU Observatory (34.85∘ N, 136.10∘ E; Japan) at a vertical resolution of 100 m and a time resolution of 10 min. Following studies carried out with the 46.5 MHz middle and upper atmosphere (MU) radar (Luce et al., 2018), we tested models used to estimate the rate of turbulence kinetic energy (TKE) dissipation ε from the Doppler spectral width in the altitude range ∼ 0.7 to 4.0 km above sea level (a.s.l.). For this purpose, we compared LQ-7-derived ε using processed data available online (http://www.rish.kyoto-u.ac.jp/radar-group/blr/shigaraki/data/, last access: 24 July 2023) with direct estimates of ε (εU) from DataHawk UAVs. The statistical results reveal the same trends as reported by Luce et al. (2018) with the MU radar, namely (1) the simple formulation based on dimensional analysis εLout=σ3/Lout, with Lout∼70 m, provides the best statistical agreement with εU. (2) The model εN predicting a σ2N law (N is Brunt–Vaïsälä frequency) for stably stratified conditions tends to overestimate for εU≲5×10-4 m2 s−3 and to underestimate for εU≳5×10-4 m2 s−3. We also tested a model εS predicting a σ2S law (S is the vertical shear of horizontal wind) supposed to be valid for low Richardson numbers (Ri=N2/S2). From the case study of a turbulent layer produced by a Kelvin–Helmholtz (K–H) instability, we found that εS and εLout are both very consistent with εU, while εN underestimates εU in the core of the turbulent layer where N is minimum. We also applied the Thorpe method from data collected from a nearly simultaneous radiosonde and tested an alternative interpretation of the Thorpe length in terms of the Corrsin length scale defined for weakly stratified turbulence. A statistical analysis showed that εS also provides better statistical agreement with εU and is much less biased than εN. Combining estimates of N and shear from DataHawk and radar data, respectively, a rough estimate of the Richardson number at a vertical resolution of 100 m (Ri100) was obtained. We performed a statistical analysis on the Ri dependence of the models. The main outcome is that εS compares well with εU for low Ri100 (Ri100≲1), while εN fails. εLout varies as εS with Ri100, so that εLout remains the best (and simplest) model in the absence of information on Ri. Also, σ appears to vary as Ri100-1/2 when Ri100≳0.4 and shows a degree of dependence on S100 (vertical shear at a vertical resolution of 100 m) otherwise.

Determination of Meteor Vector Velocity Using MU Interferometry Measurements of Head Echoes

A new method for measuring the vector velocity of meteoroids using meteor head echoes is proposed in this study. The lateral velocity is determined by utilizing the phase interference measurement between channels, while the radial velocity is obtained using a conventional Doppler frequency shift measurement. Compared to previous studies, this method does not require multi-site observations and can calculate the vector velocity of meteors in real-time. This paper provides the complete process for the inversion of the meteor vector velocity, detailing the analyzing process using MU radar head echo data. First, the MUSIC algorithm was used to estimate the DOA of the meteor target, which is a parameter required for lateral velocity measurement. Channel calibration is required before this estimation. Next, delay-Doppler matched filter processing was performed on each receiving channel’s data to determine the distance and radial velocity of the meteor target. Subsequently, the lateral velocity component was synthesized using the least squares method from the phase difference rate extracted from the matched filter output results of multiple channel pairs. Then, the vector velocity and trajectory of the meteor could be determined. The method was verified using MU radar head echo data. Different groups of channel pairs were selected for calculating the lateral velocity, and the results were found to be close, demonstrating the self-consistency of the method. Additionally, the calculated vector velocity is consistent with the direction and magnitude of the meteor’s motion trajectory, confirming the feasibility of the proposed approach. The method allows for the observation of more prominent characteristics of meteoroid motion, providing a more detailed observation capability of velocity variations in other directions than previous methods.

DDMA-MIMO Observations With the MU Radar: Validation by Measuring a Beam Broadening Effect

The phased-array radar, basically developed with the defense systems, has mainly utilized for the atmospheric radar and the wind profiling radar in the field of meteorological remote sensing, and recently it has also applied to the weather radar for research purposes. As a further development of phased-array technology, the “Multiple-Input Multiple-Output (MIMO) technique,” which has been developed in the field of communication systems, has been applied to radars. With the MIMO radar, it is possible to create a virtual antenna aperture plane beyond the actual antenna and to reduce the actual antenna size compared to that of the conventional antenna while maintaining the angular resolution. This effect is expected to reduce costs, which is one of the major hurdles in expanding phased array radars instead of parabolic antenna systems. In order to confirm the effect, an experimental observation was performed using the MU radar which is a VHF-band phased array atmospheric radar with multi-channel receivers. The MIMO technique requires orthogonal waveforms on each transmitter to identify the transmitted signals with multiple receivers, and various methods are known to realize orthogonality. In this paper, we focus on the “ Doppler Division Multiple Access (DDMA)” MIMO technique, with which slightly different frequencies are selected as transmit waveforms to separate in each receiver in the Doppler frequency domain. The observation results by measuring a beam broadening effect with the MU radar indicate that it will be a key technique for the atmospheric radar in the near future.

High-altitude meteors detected by the interferometric MU radar

ABSTRACTWe have re-analysed part of the middle and upper atmosphere (MU) radar meteor head echo data set collected during 2009–2010 and confirmed the existence of a rare high-altitude radar meteor population reaching up to ∼150 km altitude. The number of detections decreases significantly as a function of initial altitude. Out of the total amount of 106 000 events, 74 had an initial altitude >130 km while four of those had an initial altitude >145 km. High-altitude radar meteor observations have been reported before, e.g. using the EISCAT VHF radar and the Jicamarca Radio Observatory. The main novelty of this study is that the observations were performed using methods that render the final data set unambiguous in direction of arrival together with rigorously tested analysis routines that were validated by noisy raw data simulations. Due to our experimental set-up the maximum detectable range was limited to 148 km. Hence, we cannot confirm or deny the existence of radar meteors above that altitude.

Radar analysis algorithm for determining meteor head echo parameter probability distributions

ABSTRACT We present an automated radar data analysis algorithm developed to calculate probability distributions of meteor- and meteoroid parameters for head echoes detected with the Middle and Upper atmosphere (MU) radar in Shigaraki, Japan. The algorithm utilizes direct Monte Carlo simulations of uncertainties, with Bayesian Markov-chain Monte Carlo estimation of meteor model parameters and N-body propagation of distributions to perform orbit determination. The implementation has been validated using raw data simulations and a comparison with previous analysis methods. The concepts are applicable on a wide range of possible head echo measurements with other radar systems. The generated probability distributions provide quantitative reliability, which enables improved statistical studies and investigating the origins of detected meteoroids. The methodology section is highly detailed in order for the methods to be reproducible and provide a solid reference foundation for future studies. One such study is presented in a companion paper called ‘High-altitude meteors detected by the interferometric MU radar’.

Measurement of Aspect Angle of Field-Aligned Plasma Irregularities in Mid-Latitude E Region Using VHF Atmospheric Radar Imaging and Interferometry Techniques

Multireceiver and multifrequency radar imaging were carried out with the 46.5 MHz MU radar in Japan (34.85°N and 136.10°E) to examine the aspect sensitivity of field-aligned plasma irregularities (FAIs) in the mid-latitude ionosphere E region. A radar beam was directed to geographic north and at 51° zenith angle, which was normal to the geomagnetic field line around 100 km height. Nineteen receivers and five carrier frequencies were used for radar imaging to retrieve the power distribution in the radar volume, and then the aspect angle along the geomagnetic field line was calculated according to the angular power distribution. Retrieval algorithms such as Fourier, Capon, and norm-constrained Capon (NC-Capon) were employed, in which the NC-Capon was applied to FAIs for the first time and found to be more suitable for the present study. The aspect angles estimated by the NC-Capon ranged between 0.1° and 0.4° mostly, and averaged around 0.2°, which were the same order to the previous measurements with radar interferometry (RI), made for equatorial electrojet irregularities and the lower mid-latitude sporadic E region. For comparison, RI-estimated aspect angles were also investigated and found to be close to that of NC-Capon, but distributed over a wider extent of angles.

Probabilistic analysis of ambiguities in radar echo direction of arrival from meteors

Abstract. Meteors and hard targets produce coherent radar echoes. If measured with an interferometric radar system, these echoes can be used to determine the position of the target through finding the direction of arrival (DOA) of the incoming echo onto the radar. Depending on the spatial configuration of radar-receiving antennas and their individual gain patterns, there may be an ambiguity problem when determining the DOA of an echo. Radars that are theoretically ambiguity-free are known to still have ambiguities that depend on the total radar signal-to-noise ratio (SNR). In this study, we investigate robust methods which are easy to implement to determine the effect of ambiguities on any hard target DOA determination by interferometric radar systems. We apply these methods specifically to simulate four different radar systems measuring meteor head and trail echoes, using the multiple signal classification (MUSIC) DOA determination algorithm. The four radar systems are the Middle And Upper Atmosphere (MU) radar in Japan, a generic Jones 2.5λ specular meteor trail radar configuration, the Middle Atmosphere Alomar Radar System (MAARSY) radar in Norway and the Program of the Antarctic Syowa Mesosphere Stratosphere Troposphere Incoherent Scatter (PANSY) radar in the Antarctic. We also examined a slightly perturbed Jones 2.5λ configuration used as a meteor trail echo receiver for the PANSY radar. All the results are derived from simulations, and their purpose is to grant understanding of the behaviour of DOA determination. General results are as follows: there may be a region of SNRs where ambiguities are relevant; Monte Carlo simulation determines this region and if it exists; the MUSIC function peak value is directly correlated with the ambiguous region; a Bayesian method is presented that may be able to analyse echoes from this region; the DOA of echoes with SNRs larger than this region are perfectly determined; the DOA of echoes with SNRs smaller than this region completely fail to be determined; the location of this region is shifted based on the total SNR versus the channel SNR in the direction of the target; and asymmetric subgroups can cause ambiguities, even for ambiguity-free radars. For a DOA located at the zenith, the end of the ambiguous region is located at 17 dB SNR for the MU radar and 3 dB SNR for the PANSY radar. The Jones radars are usually used to measure specular trail echoes far from zenith. The ambiguous region for a DOA at 75.5∘ elevation and 0∘ azimuth ends at 12 dB SNR. Using the Bayesian method, it may be possible to analyse echoes down to 4 dB SNR for the Jones configuration when given enough data points from the same target. The PANSY meteor trail echo receiver did not deviate significantly from the generic Jones configuration. The MAARSY radar could not resolve arbitrary DOAs sufficiently well enough to determine a stable region. However, if the DOA search is restricted to 70∘ elevation or above by assumption, stable DOA determination occurs above 15 dB SNR.

Relationship between radar cross section and optical magnitude based on radar and optical simultaneous observations of faint meteors

Radar and optical simultaneous observations of meteors are important to understand the size distribution of the interplanetary dust. However, faint meteors detected by high power large aperture radar observations, which are typically as faint as 10 mag. In optical, have not been detected until recently in optical observations, mainly due to insufficient sensitivity of the optical observations. In this paper, two radar and optical simultaneous observations were organized. The first observation was carried out in 2009–2010 using Middle and Upper Atmosphere Radar (MU radar) and an image-intensified CCD camera. The second observation was carried out in 2018 using the MU radar and a mosaic CMOS camera, Tomo-e Gozen, mounted on the 1.05-m Kiso Schmidt Telescope. In total, 331 simultaneous meteors were detected. The relationship between radar cross sections and optical V-band magnitudes was well approximated by a linear function. A transformation function from the radar cross section to the V-band magnitude was derived for sporadic meteors. The transformation function was applied to about 150,000 meteors detected by the MU radar in 2009–2015, large part of which are sporadic, and a luminosity function was derived in the magnitude range of −1.5–9.5 mag. The luminosity function was well approximated by a single power-law function with the population index of r=3.52±0.12. The present observation indicates that the MU radar has capability to detect interplanetary dust of 10−5–100g in mass as meteors.

Atmospheric structures in the troposphere as revealed by high-resolution backscatter images from MU radar operating in range-imaging mode

VHF band stratosphere/troposphere (ST) radars around the globe are seldom operated in range-imaging mode. As such, the typical range resolution of their backscatter images is about 150 m. The only exception is the Kyoto University’s Middle and Upper Atmosphere (MU) radar in Shigaraki, Japan. Range imaging using frequency diversity was implemented there in 2005 and has often been used since then. During the Shigaraki UAV Radar Experiment (ShUREX) campaigns in the spring/summers of 2015, 2016, and 2017, the MU radar was operated in range-imaging mode to provide a range resolution of typically 20 m, for good signal to noise (SNR) ratios. The resulting Capon backscatter images revealed a variety of atmospheric structures in the moist troposphere in great detail. They were also quite useful in deploying in situ sensors on board unmanned aerial vehicles (UAVs) to probe such structures in near real time guided by the images. The goal of this paper is to present and discuss some such structures of interest to atmospheric dynamics collectively, to provide an overarching view. They include Kelvin-Helmholtz (KH) billows generated by shear instability, mid-level cloud-base turbulence (MCT) layers generated by convective instability in a moist troposphere, convective boundary layer (CBL), and sheet and layer (S&L) structures in a stably stratified atmospheric column. Videos of radar images collected during the 2015 and 2016 campaigns are included as Additional files 1, 2 and 3 with a brief explanatory text as Additional file 4 to demonstrate the fascinating, ever-changing evolution of atmospheric structures over the MU radar.

Midlevel Cloud‐Base Turbulence: Radar Observations and Models

AbstractMidlevel Cloud‐base Turbulence (MCT) caused by ice/snow precipitation from midlevel clouds falling into dry air below the cloud base and sublimating is investigated. MCT phenomena in the lower troposphere as revealed by the Middle and Upper atmosphere (MU) radar in Shigaraki, Japan, during the Shigaraki unmanned aerial vehicle and radar Experiment campaigns in the spring‐summer of 2015 and 2016 are described. The MU radar was operated in a high‐resolution (~20‐m) range imaging mode and hence revealed the structure of these MCTs in great detail. These MCT layers grew to hundreds of meters in thickness, often reaching nearly 2,000 m in depth, lasting often for as much as 24 hr. A simple analytical model and a second‐moment closure‐based turbulent mixing model are used to estimate the levels of turbulence kinetic energy and its dissipation rate in the MCT layer and compare them to those measured by the MU radar. The study shows that MCT can give rise to moderate levels of turbulence of potential interest to aviation, confirming conclusions reached by earlier pioneering studies. However, to our knowledge, this important process has not received the attention it deserves by the atmospheric community so far, and therefore, it is the principal goal of this paper to not only present some observations and preliminary modeling results, but also redraw the attention of atmospheric scientists to this important process.

Turbulence kinetic energy dissipation rates estimated from concurrent UAV and MU radar measurements

We tested models commonly used for estimating turbulence kinetic energy dissipation rates varepsilon from very high frequency stratosphere–troposphere radar data. These models relate the root-mean-square value sigma of radial velocity fluctuations assessed from radar Doppler spectra to varepsilon. For this purpose, we used data collected from the middle and upper atmosphere (MU) radar during the Shigaraki unmanned aerial vehicle (UAV)—radar experiment campaigns carried out at the Shigaraki MU Observatory, Japan, in June 2016 and 2017. On these occasions, UAVs equipped with fast-response and low-noise Pitot tube sensors for turbulence measurements were operated in the immediate vicinity of the MU radar. Radar-derived dissipation rates varepsilon estimated from the various models at a range resolution of 150 m from the altitude of 1.345 km up to the altitude of ~ 4.0 km, a (half width half power) beam aperture of 1.32° and a time resolution of 24.6 s, were compared to dissipation rates (varepsilon_{U}) directly obtained from relative wind speed spectra inferred from UAV measurements. Firstly, statistical analysis results revealed a very close relationship between enhancements of sigma and varepsilon_{U} for varepsilon_{U} ,{ gtrsim },10^{ - 5} ,{text{m}}^{2} ,{text{s}}^{ - 3}, indicating that both instruments detected the same turbulent events with varepsilon_{U} above this threshold. Secondly, varepsilon_{U} was found to be statistically proportional to sigma^{3}, whereas a sigma^{2} dependence is expected when the size of the largest turbulent eddies is smaller than the longitudinal and transverse dimensions of the radar sampling volume. The sigma^{3} dependence was found even after excluding convectively generated turbulence in the planetary boundary layer and below clouds. The best agreement between varepsilon_{U} and radar-derived varepsilon was obtained with the simple formulation based on dimensional analysis varepsilon = sigma^{3} /L_{c} where LC ≈ 50–70 m. This empirical expression constitutes a simple way to estimate dissipation rates in the lower troposphere from MU radar data whatever the sources of turbulence be, in clear air or cloudy conditions, consistent with UAV estimates.

On a numerical model for extracting TKE dissipation rate from very high frequency (VHF) radar spectral width

The standard model (e.g., Hocking in Earth Planets Space 51:525–541, 1999), varepsilon = c_{0} sigma^{2} N, (where sigma is the radar spectral width assumed to be equal to vertical turbulence velocity fluctuation sqrt {overline{{w^{2} }} }, N is the buoyancy frequency, and c0 is a constant), derived from Weinstock (J Atmos Sci 35:1022–1027, 1978; J Atmos Sci 38:880–883, 1981) formulation, has been used extensively for estimating the turbulence kinetic energy (TKE) dissipation rate varepsilon under stable stratification from VHF radar Doppler spectral width sigma. The Weinstock model can be derived by simply integrating the TKE spectrum in the wavenumber space from the buoyancy wavenumber k_{text{B}} = frac{N}{sigma } to infty. However, it ignores the radar volume dimensions and hence its spatial weighting characteristics. Labitt (Some basic relations concerning the radar measurements of air turbulence, MIT Lincoln Laboratory, ATC Working Paper NO 46WP-5001, 1979) and White et al. (J Atmos Ocean Technol 16:1967–1972, 1999) formulations do take into account the radar spatial weighting characteristics, but assume that the wavenumber range in the integration of TKE spectrum extends from 0 to infty. The White et al. model accounts for wind speed effects, whereas the other two do not. More importantly, all three formulations make the assumption that k−5/3 spectral shape of TKE spectrum extends across the entire wavenumber range of integration. It is traditional to use Weinstock formulation for k_{text{B}}^{ - 1} < 2a,2b (where a and b are radar volume dimensions in the horizontal and vertical directions) and White et al. formulation (without wind advection) for k_{text{B}}^{ - 1} > 2a,2b. However, there is no need to invoke these asymptotic limits. We present here a numerical model, which is valid for all values of buoyancy wavenumber k_{text{B}} and transitions from varepsilon simsigma^{2} behavior at lower values of sigma in accordance with Weinstock’s model, to varepsilon simsigma^{3} at higher values of sigma, in agreement with Chen (J Atmos Sci 31:2222–2225, 1974) and Bertin et al. (Radio Sci 32:791–804, 1997). It can also account for the effects of wind speed, as well the beam width and altitude. Following Hocking (J Atmos Terr Phys 48:655–670, 1986, Earth Planets Space 51:525–541, 1999), the model also takes into account contributions of velocity fluctuations beyond the inertial subrange. The model has universal applicability and can also be applied to convective turbulence in the atmospheric column. It can also be used to explore the parameter space and hence the influence of various parameters and assumptions on the extracted varepsilon values. In this note, we demonstrate the utility of the numerical model and make available a MATLAB code of the model for potential use by the radar community. The model results are also compared against in situ turbulence measurements using an unmanned aerial vehicle (UAV) flown in the vicinity of the MU radar in Shigaraki, Japan, during the ShUREX 2016 campaign.

Atmospheric Kelvin\u2013Helmholtz billows captured by the MU radar, lidars and a fish-eye camera

On June 11, 2015, a train of large-amplitude Kelvin–Helmholtz (KH) billows was monitored by the Middle and Upper Atmosphere (MU) radar (Shigaraki MU Observatory, Japan) at the altitude of ~ 6.5 km. Four to five KH billows in formation and decay stages were observed for about 20 min at the height of a strong speed shear (> ~ 30 m s−1km−1), just a few hundred meters above a mid-level cloud base. The turbulent billows had a spacing of about 3.5–4.0 km (3.71 km in average) and an aspect ratio (depth/spacing) of ~ 0.3. The turbulence kinetic energy dissipation rate estimated was of the order of 10–50 {text{mWkg}}^{ - 1}, corresponding to moderate turbulence according to ICAO (2010) classification. By chance, an upward-looking fish-eye camera producing pictures once every minute detected smooth protuberances at the cloud base caused by the KH billows so that comparisons of their characteristics could be made for the first time between the radar observations and the pictures. The main characteristics of the KH wave (horizontal wavelength, phase front direction and phase speed) obtained from the analysis of the pictures were fully consistent with those found from radar data. The pictures indicated that the billows were advected by the wind observed at the height of the critical level. They also revealed a very small transverse extent (about twice the horizontal spacing) suggesting that the large-amplitude KH billows were generated by a very localized source. Micro-pulse lidar and Raman–Rayleigh–Mie lidar data also collected during the event permitted us to confirm some of the characteristics of the billows.

Vertical structure of the lower troposphere derived from MU radar, unmanned aerial vehicle, and balloon measurements during ShUREX 2015

The ShUREX (Shigaraki UAV Radar Experiment) 2015 campaign carried out at the Shigaraki Middle and Upper atmosphere (MU) observatory (Japan) in June 2015 provided a unique opportunity to compare vertical profiles of atmospheric parameters estimated from unmanned aerial vehicle (UAV), balloon, and radar data in the lower troposphere. The present work is intended primarily as a demonstration of the potential offered by combination of these three instruments for studying the small-scale structure and dynamics in the lower troposphere. Here, we focus on data collected almost simultaneously by two instrumented UAVs and two meteorological balloons, near the MU radar operated continuously during the campaign. The UAVs flew along helical ascending and descending paths at a nearly constant horizontal distance from the radar (~ 1.0 km), while the balloons launched from the MU radar site drifted up to ~ 3–5 km in the altitude range of comparisons (~ 0.5 to 4.0 km) due to wind advection. Vertical profiles of squared Brünt-Väisälä frequency N2 and squared vertical gradient of generalized potential refractive index M2 were estimated at a vertical resolution of 20 m from pressure, temperature, and humidity data collected by UAVs and radiosondes. Profiles of M2 were also estimated from MU radar echo power at vertical incidence at a vertical sampling of 20 m and various time resolutions (1–4 min). The balloons and the MU radar provided vertical profiles of wind and wind shear S so that two independent estimates of the gradient Richardson number (Ri = N2/S2) could be obtained at a range resolution of 150 m. The two estimates of Ri profiles also showed remarkable agreement at all altitudes. We show that all three instruments detected the same prominent temperature and humidity gradients, down to decameter scales in stratified conditions. These gradients extended horizontally over a few kilometers at least and persisted for hours without significant changes, indicating that the turbulent diffusion was weak. Large discrepancies between N2and M2 profiles derived from the balloon, UAV, and radar data were found in a turbulent layer generated by a Kelvin-Helmholtz (KH) shear flow instability in the height range from 1.80 to 2.15 km. The cause of these discrepancies appears to depend on the stage of the KH billows.

Doppler Profile Tracing Using MPCF on MU Radar and Sodar: Performance Analysis

Wind profilers (WPs) are ground-based pulsed Doppler radars operating in the UHF/very high frequency (VHF) bands. They use backscatter from clear-air turbulence as a tracer of background wind in the troposphere and up to mesosphere. Wind profiling systems range from large research radar systems to smaller operational units. In recent times, multiparameter cost function (MPCF) has emerged as a computationally efficient Doppler profile estimation method. It has been claimed that the MPCF method can be easily migrated to any vertically sounding clear-air wind profiling system that works in Doppler beam swinging mode. In order to investigate this claim, MPCF was applied to the following wind profiling systems: the middle and upper atmosphere (MU) radar, located at Shigaraki, Japan, which is an active phased array VHF radar system with 475 transceivers, and Doppler sodar, located at Pune, India, which uses acoustic frequency to obtain echoes in the planetary boundary layer. Both these systems are complementary in the sense that they cover ranges from ground to about 80 km. During the experimentation, the MPCF algorithm did not need any change of parameters except the matching of the data reading formats. The results of MPCF on the MU radar were validated with Radiosonde data. These results indicate that the MPCF works seamlessly on all types of WP systems irrespective of the carrier, range, and radar type.

On the Performance of the Range Imaging Technique Estimated Using Unmanned Aerial Vehicles During the ShUREX 2015 Campaign

For the first time, the performance of the range imaging technique of the very high frequency band middle and upper atmosphere (MU) radar (Shigaraki observatory, Japan), when using frequency diversity, is assessed. This is done by the detection of unmanned aerial vehicles (UAVs) operated near the radar during the Shigaraki UAV radar experiment campaign, carried out from June 1, 2015 to June 14, 2015. During this campaign devoted to the measurements of fine-scale turbulence and stability in the lower troposphere using the turbulence sensors on the DataHawk UAV and the MU radar, the detection of the small UAV by the MU radar provided an excellent opportunity for taking stock of the range imaging technique in the presence of a single hard target. It was found that range imaging reproduces a faithful image of the aircraft position and its displacements with an excellent accuracy (of the order of ~10 m), giving extra credence to the thin echo layers and their vertical displacements generally observed from the range imaging technique in stably stratified conditions.

Shigaraki UAV-Radar Experiment (ShUREX): overview of the campaign with some preliminary results

The Shigaraki unmanned aerial vehicle (UAV)-Radar Experiment (ShUREX) is an international (USA-Japan-France) observational campaign, whose overarching goal is to demonstrate the utility of small, lightweight, inexpensive, autonomous UAVs in probing and monitoring the lower troposphere and to promote synergistic use of UAVs and very high frequency (VHF) radars. The 2-week campaign lasting from June 1 to June 14, 2015, was carried out at the Middle and Upper Atmosphere (MU) Observatory in Shigaraki, Japan. During the campaign, the DataHawk UAV, developed at the University of Colorado, Boulder, and equipped with high-frequency response cold wire and pitot tube sensors (as well as an iMET radiosonde), was flown near and over the VHF-band MU radar. Measurements in the atmospheric column in the immediate vicinity of the radar were obtained. Simultaneous and continuous operation of the radar in range imaging mode enabled fine-scale structures in the atmosphere to be visualized by the radar. It also permitted the UAV to be commanded to sample interesting structures, guided in near real time by the radar images. This overview provides a description of the ShUREX campaign and some interesting but preliminary results of the very first simultaneous and intensive probing of turbulent structures by UAVs and the MU radar. The campaign demonstrated the validity and utility of the radar range imaging technique in obtaining very high vertical resolution (~20 m) images of echo power in the atmospheric column, which display evolving fine-scale atmospheric structures in unprecedented detail. The campaign also permitted for the very first time the evaluation of the consistency of turbulent kinetic energy dissipation rates in turbulent structures inferred from the spectral broadening of the backscattered radar signal and direct, in situ measurements by the high-frequency response velocity sensor on the UAV. The data also enabled other turbulence parameters such as the temperature structure function parameter {C}_T^2 and refractive index structure function parameter {C}_n^2 to be measured by sensors on the UAV, along with radar-inferred refractive index structure function parameter {C}_{n,mathrm{radar}}^2 . The comprehensive dataset collected during the campaign (from the radar, the UAV, the boundary layer lidar, the ceilometer, and radiosondes) is expected to help obtain a better understanding of turbulent atmospheric structures, as well as arrive at a better interpretation of the radar data.

Comparisons between high-resolution profiles of squared refractive index gradient &amp;lt;i&amp;gt;M&amp;lt;/i&amp;gt;&amp;lt;sup&amp;gt;2&amp;lt;/sup&amp;gt; measured by the Middle and Upper Atmosphere Radar and unmanned aerial vehicles (UAVs) during the Shigaraki UAV-Radar Experiment 2015 campaign

Abstract. New comparisons between the square of the generalized potential refractive index gradient M2, estimated from the very high-frequency (VHF) Middle and Upper Atmosphere (MU) Radar, located at Shigaraki, Japan, and unmanned aerial vehicle (UAV) measurements are presented. These comparisons were performed at unprecedented temporal and range resolutions (1–4 min and ∼ 20 m, respectively) in the altitude range ∼ 1.27–4.5 km from simultaneous and nearly collocated measurements made during the ShUREX (Shigaraki UAV-Radar Experiment) 2015 campaign. Seven consecutive UAV flights made during daytime on 7 June 2015 were used for this purpose. The MU Radar was operated in range imaging mode for improving the range resolution at vertical incidence (typically a few tens of meters). The proportionality of the radar echo power to M2 is reported for the first time at such high time and range resolutions for stratified conditions for which Fresnel scatter or a reflection mechanism is expected. In more complex features obtained for a range of turbulent layers generated by shear instabilities or associated with convective cloud cells, M2 estimated from UAV data does not reproduce observed radar echo power profiles. Proposed interpretations of this discrepancy are presented.

A scanning Raman lidar for observing the spatio-temporal distribution of water vapor

We have constructed a scanning Raman lidar to observe the cross-sectional distribution of the water vapor mixing ratio and aerosols near the Earth's surface, which are difficult to observe when a conventional Raman lidar system is used. The Raman lidar is designed for a nighttime operating system by employing a ultra-violet (UV) laser source and can measure the water vapor mixing ratio at an altitude up to 7km using vertically pointing observations. The scanning mirror system consists of reflective flat mirrors and a rotational stage. By using a program-controlled rotational stage, a vertical scan can be operated with a speed of 1.5°/s. The beam was pointed at 33 angles over range of 0–48° for the elevation angle with a constant step width of 1.5°. The range-height cross sections of the water vapor and aerosol within a 400m range can be obtained for 25min. The lidar signals at each direction were individually smoothed with the moving average to spread proportionally with the distance from the laser-emitting point. The averaged range at a distance of 200m (400m) from the lidar was 30.0m (67.5m) along the lidar signal in a specific direction.The experimental observations using the scanning lidar were conducted at night in the Shigaraki MU radar observatory located on a plateau with undulating topography and surrounded by forests. The root mean square error (RMSE) between the temporal variations of the water vapor mixing ratio by the scanning Raman lidar and by an in-situ weather sensor equipped with a tethered balloon was 0.17g/kg at an altitude of 100m. In cross-sectional measurements taken at altitudes and horizontal distances up to 400m from the observatory, we found that the water vapor mixing ratio above and within the surface layer varied vertically and horizontally. The spatio-temporal variability of water vapor near the surface seemed to be sensitive to topographic variations as well as the wind field and the temperature gradient over the site. From the wide-range cross-sectional observations of the water vapor mixing ratio and the backscatter ratio of aerosols within a 2000m range, we can detect small-scale water vapor structures on a horizontal scale of several hundred meters in the atmospheric boundary layer.

Evaluation of momentum flux with radar

The statistics of gravity wave momentum flux estimation are investigated using data from the MU radar at Shigariki, Japan (136°E, 35°N). The radar has been operating during campaign periods since 1986. The first part of the paper focuses on a multi-day campaign during October 13–31, 1986. The second part of the paper investigates data after 2006 when the radar was operated in a meteor scatter mode. Momentum fluxes are derived from both the turbulent scatter and the meteor scatter measurements, but the techniques are quite different. Probability Distribution Functions are formed using turbulent scatter data. These show that wave packets sometimes have momentum flux magnitudes in excess of 100m2s−2. The technique for meteor radars, introduced by Hocking (2005), has been widely adopted by the radar community in recent years. The momentum flux estimated using this technique is found to be anti-correlated with the background tidal winds. A validation investigation is carried out for periods with a high meteor echo data rate. The conclusion was that the method can be used to calculate the sign of momentum flux, but does not accurately specify the magnitude.