A time-series method to identify and correct range sidelobes in meteorological radar data

The use of pulse compression techniques to improve the sensitivity of meteorological radars has become increasingly common in recent years. An unavoidable side-effect of such techniques is the formation of ‘range sidelobes’ which lead to spreading of information across several range gates. These artefacts are particularly troublesome in regions where there is a sharp gradient in the power backscattered to the antenna as a function of range. In this article we present a simple method for identifying and correcting range sidelobe artefacts. We make use of the fact that meteorological targets produce an echo which fluctuates at random, and that this echo, like a fingerprint, is unique to each range gate. By cross-correlating the echo time series from pairs of gates therefore we can identify whether information from one gate has spread into another, and hence flag regions of contamination. In addition we show that the correlation coefficients contain quantitative information about the fraction of power leaked from one range gate to another, and we propose a simple algorithm to correct the corrupted reflectivity profile

High-precision measurements of the co-polar correlation coefficient: non-Gaussian errors and retrieval of the dispersion parameter µ in rainfall

The co-polar correlation coefficient (ρhv) has many applications, including hydrometeor classification, ground clutter and melting layer identification, interpretation of ice microphysics and the retrieval of rain drop size distributions (DSDs). However, we currently lack the quantitative error estimates that are necessary if these applications are to be fully exploited. Previous error estimates of ρhv rely on knowledge of the unknown "true" ρhv and implicitly assume a Gaussian probability distribution function of ρhv samples. We show that frequency distributions of ρhv estimates are in fact highly negatively skewed. A new variable: L = -log10(1 - ρhv) is defined, which does have Gaussian error statistics, and a standard deviation depending only on the number of independent radar pulses. This is verified using observations of spherical drizzle drops, allowing, for the first time, the construction of rigorous confidence intervals in estimates of ρhv. In addition, we demonstrate how the imperfect co-location of the horizontal and vertical polarisation sample volumes may be accounted for. The possibility of using L to estimate the dispersion parameter (µ) in the gamma drop size distribution is investigated. We find that including drop oscillations is essential for this application, otherwise there could be biases in retrieved µ of up to ~8. Preliminary results in rainfall are presented. In a convective rain case study, our estimates show µ to be substantially larger than 0 (an exponential DSD). In this particular rain event, rain rate would be overestimated by up to 50% if a simple exponential DSD is assumed

A robust automated technique for operational calibration of ceilometers using the integrated backscatter from totally attenuating liquid clouds

A simple and robust method for calibrating ceilometers has been tested in an operational environment demonstrating that the calibrations are stable to better than ± 5 % over a period of a year. The method relies on using the integrated backscatter (B) from liquid clouds that totally extinguish the ceilometer signal; B is inversely proportional to the lidar ratio (S) of the backscatter to the extinction for cloud droplets. The calibration technique involves scaling the observed backscatter so that B matches the predicted value for S of 18.8 ± 0.8 sr for cloud droplets, at ceilometer wavelengths. For accurate calibration, care must be taken to exclude any profiles having targets with different values of S, such as drizzle drops and aerosol particles, profiles that do not totally extinguish the ceilometer signal, profiles with low cloud bases that saturate the receiver, and any profiles where the window transmission or the lidar pulse energy is low. A range dependent multiple scattering correction that depends on the ceilometer optics should be applied to the profile. A simple correction for water vapour attenuation for ceilometers operating at around 910 nm wavelength is applied to the signal using the vapour profiles from a forecast analysis. For a generic ceilometer in the UK the 90-day running mean of the calibration coefficient over a period of 20 months is constant to within 3 % with no detectable annual cycle, thus confirming the validity of the humidity and multiple scattering correction. For Gibraltar, where cloud cover is less prevalent than in the UK, the 90-day running mean calibration coefficient was constant to within 4 %. The more sensitive ceilometer model operating at 1064 nm is unaffected by water vapour attenuation but is more prone to saturation in liquid clouds. We show that reliable calibration is still possible, provided the clouds used are above a certain altitude. The threshold is instrument dependent but is typically around 2 km. We also identify a characteristic signature of saturation, and remove any profiles with this signature. Despite the more restricted sample of cloud profiles, a robust calibration is readily achieved, and, in the UK, the running mean 90-day calibration coefficients varied by about 4 % over a period of one year. The consistency of profiles observed by nine pairs of co-located ceilometers in the UK Met Office network operating at around 910 nm and 1064 nm provided independent validation of the calibration technique. EUMETNET is currently networking 700 European ceilometers so they can provide ceilometer profiles in near real time to European weather forecast centres and has adopted the cloud calibration technique described in this paper for ceilometers with a wavelength of around 910 nm

Statistics of convective cloud turbulence from a comprehensive turbulence retrieval method for radar observations

Turbulent mixing processes are important in determining the evolution of convective clouds,and the production of convective precipitation. However, the exact nature of these impacts remains uncertain due to limited observations. Model simulations show that assumptions made in parametrizing turbulence can have a marked effect on the characteristics of simulated clouds. This leads to significant uncertainty in forecasts from convection‐permitting numerical weather prediction (NWP) models. This contribution presents a comprehensive method to retrieve turbulence using Doppler weather radar to investigate turbulence in observed clouds. This method involves isolating the turbulent component of the Doppler velocity spectrum width, expressing turbulence intensity as an eddy dissipation rate, ϵ. By applying this method throughout large datasets of observations collected over the southern United Kingdom using the (0.28° beam‐width) Chilbolton Advanced Meteorological Radar (CAMRa), statistics of convective cloud turbulence are presented. Two contrasting case days are examined: a shallow “shower” case, and a “deep convection” case, exhibiting stronger and deeper updraughts. In our observations, ϵ generally ranges from 10−3 to 10−1 m2/s3, with the largest values found within, around and above convective updraughts. Vertical profiles of ϵ suggest that turbulence is much stronger in deep convection; 95th percentile values increase with height from 0.03 to 0.1 m2/s3, compared to approximately constant values of 0.02–0.03 m2/s3 throughout the depth of shower cloud. In updraught regions on both days, the 95th percentile of ϵ has significant (p < 10−3) positive correlations with the updraught velocity, and the horizontal shear in the updraught velocity, with weaker positive correlations with updraught dimensions. The ϵ‐retrieval method presented considers a very broad range of conditions, providing a reliable framework for turbulence retrieval using high‐resolution Doppler weather radar. In applying this method across many observations, the derived turbulence statistics will form the basis for evaluating the parametrization of turbulence in NWP models

An accurate and computationally cheap microwave scattering method for ice aggregates: the Independent Monomer Approximation

The Discrete Dipole Approximation (DDA) is widely used to simulate scattering of microwaves by snowflakes, by discretising the snowflake into small “dipoles” which oscillate in response to (i) the incident wave and (ii) scattered waves from all the other dipoles in the particle. It is this coupling between all dipole pairs which makes solving the DDA system computationally expensive, and that cost grows non‐linearly as the number of crystals n within an aggregate is increased. Motivated by this, many studies have ignored the dipole coupling (the Rayleigh‐Gans Approximation, RGA). However, use of RGA leads to systematic underestimation of both scattering and absorption, and an inability to predict polarimetric properties. To address this, we present a new approach (the Independent Monomer Approximation, IMA) which solves the DDA system for each crystal “monomer” separately, then combines them to construct the full solution. By including intra‐monomer coupling, but neglecting inter‐monomer coupling, we save a factor of n in computation time over DDA. Benchmarking IMA against DDA solutions indicates that its accuracy is greatly superior to RGA, and provides ensemble scattering cross sections which closely agree with their more expensive DDA counterparts, particularly at size parameters smaller than ∼5. Addition of rime to the aggregates does not significantly degrade the results, despite the increased density. The use of IMA for radar remote sensing is evaluated, and we show that multi‐wavelength and multi‐polarisation parameters are successfully captured to within a few tenths of a dB for aggregates probed with frequencies between 3 and 200GHz, in contrast to RGA where errors of up to 2.5dB are observed. Finally we explore the realism of the IMA solutions in greater detail by analysing internal electric fields, and discuss some broader insights that IMA provides into the physical features of aggregates that are important for microwave scattering

Use of ZDR columns for early detection of severe convection within the operational radar network of the United Kingdom

Differential reflectivity () columns were observed using a Met Office three‐dimensional radar composite. An algorithm for automatic detection of columns was developed, based on dB and dBZ. Across three case days, detected columns were found to precede severe convection in tracked convective cells with a range of lead times from 0 to 20 min depending on the case day. Requiring maxima above 1.4 dB and 30 dBZ of and respectively was an appropriate second condition for all three cases although the skill in the early detection of severe convection varied across case days. Despite the high probability of detections, the high false alarm rate accompanied by low critical success index and data latency limit performance based on the three cases considered in this study. Nevertheless, the ability to detect columns in operational radar data with a useful lead time prior to severe convection in certain conditions is a promising development towards advancing nowcasting of severe convection in the United Kingdom

Evolutionary Diversification of SPANX-N Sperm Protein Gene Structure and Expression

The sperm protein associated with nucleus in the X chromosome (SPANX) genes cluster at Xq27 in two subfamilies, SPANX-A/D and SPANX-N. SPANX-A/D is specific for hominoids and is fairly well characterized. The SPANX-N gave rise to SPANX-A/D in the hominoid lineage ∼7 MYA. Given the proposed role of SPANX genes in spermatogenesis, we have extended studies to SPANX-N gene evolution, variation, regulation of expression, and intra-sperm localization. By immunofluorescence analysis, SPANX-N proteins are localized in post-meiotic spermatids exclusively, like SPANX-A/D. But in contrast to SPANX-A/D, SPANX-N are found in all ejaculated spermatozoa rather than only in a subpopulation, are localized in the acrosome rather than in the nuclear envelope, and are expressed at a low level in several nongametogenic adult tissues as well as many cancers. Presence of a binding site for CTCF and its testis-specific paralogue BORIS in the SPANX promoters suggests, by analogy to MAGE-A1 and NY-ESO-1, that their activation in spermatogenesis is mediated by the programmed replacement of CTCF by BORIS. Based on the relative density of CpG, the more extended expression of SPANX-N compared to SPANX-A/D in nongametogenic tissues is likely attributed to differences in promoter methylation. Our findings suggest that the recent duplication of SPANX genes in hominoids was accompanied by different localization of SPANX-N proteins in post-meiotic sperm and additional expression in several nongonadal tissues. This suggests a corresponding functional diversification of SPANX gene families in hominoids. SPANX proteins thus provide unique targets to investigate their roles in the function of spermatozoa, selected malignancies, and for SPANX-N, in other tissues as well

In-home evaluation of efficacy and titration of a mandibular advancement device for obstructive sleep apnea

There is increasing evidence that mandibular advancement devices (MADs) can be an effective treatment for some patients with obstructive sleep apnea, a highly prevalent chronic disease. In this study, the objectives were to objectively assess the effectiveness of MAD therapy using a limited channel recorder, and to develop a model for identifying patients who may be appropriate for MAD therapy as the initial treatment option. Thirty patients were prospectively recruited and studied at two independent dentist offices and the participants’ homes. Subjects wore the ARES Unicorder for two nights before insertion of the MAD, and again when the dentist determined that the patient had reached the titration endpoint. Self-reported measures of depression, sleepiness, and quality of life were obtained pre- and posttreatment. The reviewer was blinded to the study status while the physiological signals were being visually inspected. Significant reductions in the apnea/hypopnea index (AHI), hypoxemia measures, and snoring level were observed posttreatment. Twenty-seven of the 30 (90%) patients had a posttreatment AHI (using a 4% desaturation for hypopneas) below a clinical cut-off of 10. All but one patient (97%) exhibited at least a 50% decrease in AHI or had a posttreatment AHI ≤ 10. Significant differences in body mass index, weight, and neck circumference in patients with posttreatment AHIs above and below a clinical cut-off of five were identified. The linear regression used to predict the posttreatment AHI using pretreatment data resulted in an R2 of 0.68. The model correctly predicted two patients who were unable to obtain a posttreatment AHI of 10 or less. This study was designed to demonstrate two models of collaboration between a dental sleep medicine specialist and a sleep medicine physician in the monitoring of a patient treated with a MAD. The outcome data suggest that the limited channel recording system can be used as an alternative to laboratory polysomnography to reduce the cost of MAD treatment, and to improve the quality and consistency of posttreatment patient care