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.

Determinación de longitudes efectivas de lluvia en trayectos de enlaces vía satélite a partir de los datos de reflectividad del Radar meteorológico de Kapildui

[ES]En este proyecto se presenta un estudio sobre la estimación de la longitud efectiva de lluvia derivada de los escaneos de elevación obtenidos por el radar meteorológico de Kapildui, en Álava. Se estudia la altura y la longitud de la lluvia para distintos eventos: para lluvia estratiforme y para lluvia convectiva. Se analizará la variabilidad espacial y temporal para diferentes ángulos de elevación del radar. Finalmente, se presentará una versión del algoritmo implementado para el cálculo de longitudes efectivas de lluvia y se realizará un estudio estadístico de la variabilidad de ésta para diferentes direcciones y con diferentes eventos de lluvia.

Anàlisi d'observacions de superfície en relació al diagnòstic de tipus de precipitació observat amb radar

Màster en Meteorologia

Convective updraught evaluation in high-resolution NWP simulations using single-Doppler radar measurements

This study presents an evaluation of the size and strength of convective updraughts in high-resolution simulations by the UK Met Office Unified Model (UM). Updraught velocities have been estimated from range–height indicator (RHI) Doppler velocity measurements using the Chilbolton advanced meteorological radar, as part of the Dynamical and Microphysical Evolution of Convective Storms (DYMECS) project. Based on mass continuity and the vertical integration of the observed radial convergence, vertical velocities tend to be underestimated for convective clouds due to the undetected cross-radial convergence. Velocity fields from the UM at a resolution corresponding to the radar observations are used to scale such estimates to mitigate the inherent biases. The analysis of more than 100 observed and simulated storms indicates that the horizontal scale of updraughts in simulations tend to decrease with grid length; the 200 m grid length agreed most closely with the observations. Typical updraught mass fluxes in the 500 m grid length simulations were up to an order of magnitude greater than observed, and greater still in the 1.5 km grid length simulations. The effect of increasing the mixing length in the sub-grid turbulence scheme depends on the grid length. For the 1.5 km simulations, updraughts were weakened though their horizontal scale remained largely unchanged. Progressively more so for the sub-kilometre grid lengths, updraughts were broadened and intensified; horizontal scale was now determined by the mixing length rather than the grid length. In general, simulated updraughts were found to weaken too quickly with height. The findings were supported by the analysis of the widths of reflectivity patterns in both the simulations and observations.

Análise de superfície de tendência aplicada à chuva, medida por radar meteorológico, nas regiões de Assis e Piracicaba, SP

This paper presents the results from lineal trend surface analysis technique application. The purpose was to detect positive and negative anomalies in the rain measure distribution obtained by the meteorological radar Doppler, band S, located in Bauru, during the period of 21 of October/2004 to 29 of April/2005 in the areas of Assis and Piracicaba. Using three Z-R radar relations for rain quantification was chosen the specific equation Z = 32R1,65, as the best one. The results showed that the applied methodology was able to indicate the space distribution of the rain accumulated, identifying and locating the regions where there was rainy excess and rainy lack during each analyzed period. Such results indicate areas with larger pluvial impact and consequently more favorable for environmental damages.

The three-dimensional morphology of simulated and observed convective storms over southern England

A set of high-resolution radar observations of convective storms has been collected to evaluate such storms in the UK Met Office Unified Model during the DYMECS project (Dynamical and Microphysical Evolution of Convective Storms). The 3-GHz Chilbolton Advanced Meteorological Radar was set up with a scan-scheduling algorithm to automatically track convective storms identified in real-time from the operational rainfall radar network. More than 1,000 storm observations gathered over fifteen days in 2011 and 2012 are used to evaluate the model under various synoptic conditions supporting convection. In terms of the detailed three-dimensional morphology, storms in the 1500-m grid-length simulations are shown to produce horizontal structures a factor 1.5–2 wider compared to radar observations. A set of nested model runs at grid lengths down to 100m show that the models converge in terms of storm width, but the storm structures in the simulations with the smallest grid lengths are too narrow and too intense compared to the radar observations. The modelled storms were surrounded by a region of drizzle without ice reflectivities above 0 dBZ aloft, which was related to the dominance of ice crystals and was improved by allowing only aggregates as an ice particle habit. Simulations with graupel outperformed the standard configuration for heavy-rain profiles, but the storm structures were a factor 2 too wide and the convective cores 2 km too deep.

The role of evaporating showers in the transfer of sting-jet momentum to the surface

The most damaging winds in a severe extratropical cyclone often occur just ahead of the evaporating ends of cloud filaments emanating from the so-called cloud head. These winds are associated with low-level jets (LLJs), sometimes occurring just above the boundary layer. The question then arises as to how the high momentum is transferred to the surface. An opportunity to address this question arose when the severe ‘St Jude's Day’ windstorm travelled across southern England on 28 October 2013. We have carried out a mesoanalysis of a network of 1 min resolution automatic weather stations and high-resolution Doppler radar scans from the sensitive S-band Chilbolton Advanced Meteorological Radar (CAMRa), along with satellite and radar network imagery and numerical weather prediction products. We show that, although the damaging winds occurred in a relatively dry region of the cyclone, there was evidence within the LLJ of abundant precipitation residues from shallow convective clouds that were evaporating in a localized region of descent. We find that pockets of high momentum were transported towards the surface by the few remaining actively precipitating convective clouds within the LLJ and also by precipitation-free convection in the boundary layer that was able to entrain evaporatively cooled air from the LLJ. The boundary-layer convection was organized in along-wind rolls separated by 500 to about 3000 m, the spacing varying according to the vertical extent of the convection. The spacing was greatest where the strongest winds penetrated to the surface. A run with a medium-resolution version of the Weather Research and Forecasting (WRF) model was able to reproduce the properties of the observed LLJ. It confirmed the LLJ to be a sting jet, which descended over the leading edge of a weaker cold-conveyor-belt jet.

Lidar observation campaign of sugar cane fires and industrial emissions in the State of São Paulo, Brazil

Brazil has an important role in the biomass burning, with the detection of approximately 100,000 burning spots in a single year (2007). Most of these spots occur in the southern part of the Amazon basin during the dry season (from August to november) and these emissions reach the southeast of the country, a highly populated region and with serious urban air pollution problems. With the growing demand on biofuels, sugarcane is considerably expanding in the state of São Paulo, being a strong contributor to the bad air quality in this region. In the state of São Paulo, the main land use are pasture and sugarcane crop, that covers around 50% and 10% of the total area, respectively. Despite the aerosol from sugarcane burning having reduced atmospheric residence time, from a few days to some weeks, they might get together with those aerosol which spread over long distances (hundreds to thousands of kilometers). In the period of June through February 2010 a LIDAR observation campaign was carried in the state of São Paulo, Brazil, in order to observe and characterize optically the aerosols from two distinct sources, namely, sugar cane biomass burning and industrial emissions. For this purpose 2 LIDAR systems were available, one mobile and the other placed in a laboratory, both working in the visible (532 nm) and additionally the mobile system had a Raman channel available (607 nm). Also this campaign counted with a SODAR, a meteorological RADAR specially set up to detect aerosol echoes and gas-particle analyzers. To guarantee a good regional coverage 4 distinct sites were available to deploy the instruments, 2 in the near field of biomass burning activities (Rio Claro and Bauru), one for industrial emissions (Cubatão) and others from urban sources (São Paulo). The whole campaign provide the equivalent of 30 days of measurements which allowed us to get aerosol optical properties such as backscattering/extinction coefficients, scatter and LIDAR ratios, those were used to correlate with air quality and meteorological indicators and quantities. In this paper we should focus on the preliminary results of the Raman LIDAR system and its derived aerosol optical quantities. © 2010 Copyright SPIE - The International Society for Optical Engineering.

Um filtro gaussiano adaptativo para análise comparativa de imagens de radares e satélites meteorológicos

Meteorological satellite and radar data comparative analysis allows to correlate the precipitation structures observed in both images. Such analysis would make feasible the extension of the range of ground-based meteorological radars. In addition to the different spatial and temporal resolution of these images this comparative analysis presents difficulties due to the effects of rotation and distortion, besides the different formats, projections, and coordinate systems. This work employed an approach based on a Gaussian adaptive filter in order to compare such images. The statistical results obtained from the comparison of the images are matched to those produced by other methods.

Initial analysis from a Lidar observation campaign of sugar cane fires in the central and western portion of the S̃ao Paulo State - Brazil

The central and western portion of the S̃ao Paulo State has large areas of sugar cane plantations, and due to the growing demand for biofuels, the production is increasing every year. During the harvest period some plantation areas are burnt a few hours before the manual cutting, causing significant quantities of biomass burning aerosol to be injected into the atmosphere. During August 2010, a field campaign has been carried out in Ourinhos, situated in the south-western region of S̃ao Paulo State. A 2-channel Raman Lidar system and two meteorological S-Band Doppler Radars are used to indentify and quantify the biomass burning plumes. In addiction, CALIPSO Satellite observations were used to compare the aerosol optical properties detected in that region with those retrieved by Raman Lidar system. Although the campaign yielded 30 days of measurements, this paper will be focusing only one case study, when aerosols released from nearby sugar cane fires were detected by the Lidar system during a CALIPSO overpass. The meteorological radar, installed in Bauru, approximately 110 km northeast from the experimental site, had recorded echoes (dense smoke comprising aerosols) from several fires occurring close to the Raman Lidar system, which also detected an intense load of aerosol in the atmosphere. HYSPLIT model forward trajectories presented a strong indication that both instruments have measured the same air masss parcels, corroborated with the Lidar Ratio values from the 532 nm elastic and 607 nm Raman N2 channel analyses and data retrieved from CALIPSO have indicated the predominance of aerosol from biomass burning sources. © 2011 SPIE.

Assimilation of data from conventional and non conventional networks through a LETKF scheme in the COSMO model

The present work studies a km-scale data assimilation scheme based on a LETKF developed for the COSMO model. The aim is to evaluate the impact of the assimilation of two different types of data: temperature, humidity, pressure and wind data from conventional networks (SYNOP, TEMP, AIREP reports) and 3d reflectivity from radar volume. A 3-hourly continuous assimilation cycle has been implemented over an Italian domain, based on a 20 member ensemble, with boundary conditions provided from ECMWF ENS. Three different experiments have been run for evaluating the performance of the assimilation on one week in October 2014 during which Genova flood and Parma flood took place: a control run of the data assimilation cycle with assimilation of data from conventional networks only, a second run in which the SPPT scheme is activated into the COSMO model, a third run in which also reflectivity volumes from meteorological radar are assimilated. Objective evaluation of the experiments has been carried out both on case studies and on the entire week: check of the analysis increments, computing the Desroziers statistics for SYNOP, TEMP, AIREP and RADAR, over the Italian domain, verification of the analyses against data not assimilated (temperature at the lowest model level objectively verified against SYNOP data), and objective verification of the deterministic forecasts initialised with the KENDA analyses for each of the three experiments.

Aircraft hazard assessment from a clear-air radar and meteorological tower study of gravity wave events /

"June 1983."

Towards the Use of Radar Winds in Numerical Weather Prediction

Numerical weather prediction (NWP) models provide the basis for weather forecasting by simulating the evolution of the atmospheric state. A good forecast requires that the initial state of the atmosphere is known accurately, and that the NWP model is a realistic representation of the atmosphere. Data assimilation methods are used to produce initial conditions for NWP models. The NWP model background field, typically a short-range forecast, is updated with observations in a statistically optimal way. The objective in this thesis has been to develope methods in order to allow data assimilation of Doppler radar radial wind observations. The work has been carried out in the High Resolution Limited Area Model (HIRLAM) 3-dimensional variational data assimilation framework. Observation modelling is a key element in exploiting indirect observations of the model variables. In the radar radial wind observation modelling, the vertical model wind profile is interpolated to the observation location, and the projection of the model wind vector on the radar pulse path is calculated. The vertical broadening of the radar pulse volume, and the bending of the radar pulse path due to atmospheric conditions are taken into account. Radar radial wind observations are modelled within observation errors which consist of instrumental, modelling, and representativeness errors. Systematic and random modelling errors can be minimized by accurate observation modelling. The impact of the random part of the instrumental and representativeness errors can be decreased by calculating spatial averages from the raw observations. Model experiments indicate that the spatial averaging clearly improves the fit of the radial wind observations to the model in terms of observation minus model background (OmB) standard deviation. Monitoring the quality of the observations is an important aspect, especially when a new observation type is introduced into a data assimilation system. Calculating the bias for radial wind observations in a conventional way can result in zero even in case there are systematic differences in the wind speed and/or direction. A bias estimation method designed for this observation type is introduced in the thesis. Doppler radar radial wind observation modelling, together with the bias estimation method, enables the exploitation of the radial wind observations also for NWP model validation. The one-month model experiments performed with the HIRLAM model versions differing only in a surface stress parameterization detail indicate that the use of radar wind observations in NWP model validation is very beneficial.

On the use of weather radar for mesoscale applications in northern conditions

Mesoscale weather phenomena, such as the sea breeze circulation or lake effect snow bands, are typically too large to be observed at one point, yet too small to be caught in a traditional network of weather stations. Hence, the weather radar is one of the best tools for observing, analyzing and understanding their behavior and development. A weather radar network is a complex system, which has many structural and technical features to be tuned, from the location of each radar to the number of pulses averaged in the signal processing. These design parameters have no universal optimal values, but their selection depends on the nature of the weather phenomena to be monitored as well as on the applications for which the data will be used. The priorities and critical values are different for forest fire forecasting, aviation weather service or the planning of snow ploughing, to name a few radar-based applications. The main objective of the work performed within this thesis has been to combine knowledge of technical properties of the radar systems and our understanding of weather conditions in order to produce better applications able to efficiently support decision making in service duties for modern society related to weather and safety in northern conditions. When a new application is developed, it must be tested against ground truth . Two new verification approaches for radar-based hail estimates are introduced in this thesis. For mesoscale applications, finding the representative reference can be challenging since these phenomena are by definition difficult to catch with surface observations. Hence, almost any valuable information, which can be distilled from unconventional data sources such as newspapers and holiday shots is welcome. However, as important as getting data is to obtain estimates of data quality, and to judge to what extent the two disparate information sources can be compared. The presented new applications do not rely on radar data alone, but ingest information from auxiliary sources such as temperature fields. The author concludes that in the future the radar will continue to be a key source of data and information especially when used together in an effective way with other meteorological data.

Integration of speckle de-noising and image segmentation using synthetic aperture radar image for flood extent extraction

Flood is one of the detrimental hydro-meteorological threats to mankind. This compels very efficient flood assessment models. In this paper, we propose remote sensing based flood assessment using Synthetic Aperture Radar (SAR) image because of its imperviousness to unfavourable weather conditions. However, they suffer from the speckle noise. Hence, the processing of SAR image is applied in two stages: speckle removal filters and image segmentation methods for flood mapping. The speckle noise has been reduced with the help of Lee, Frost and Gamma MAP filters. A performance comparison of these speckle removal filters is presented. From the results obtained, we deduce that the Gamma MAP is reliable. The selected Gamma MAP filtered image is segmented using Gray Level Co-occurrence Matrix (GLCM) and Mean Shift Segmentation (MSS). The GLCM is a texture analysis method that separates the image pixels into water and non-water groups based on their spectral feature whereas MSS is a gradient ascent method, here segmentation is carried out using spectral and spatial information. As test case, Kosi river flood is considered in our study. From the segmentation result of both these methods are comprehensively analysed and concluded that the MSS is efficient for flood mapping.