Estimating radar reflectivity - Snowfall rate relationships and their uncertainties over Antarctica by combining disdrometer and radar observations
Not provided.
SEARCH · Engineering Papers
Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.
Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.
Not provided.
The Earth Cloud, Aerosol and Radiation Explorer (EarthCARE) satellite developed by the European Space Agency (ESA) and the Japan Aerospace Exploration Agency (JAXA) launched in May 2024 carries a novel 94 GHz cloud profiling radar (CPR) with Doppler capability. This work describes the open-source instrument simulator Orbital-Radar, which transforms high-resolution radar data from field observations or forward simulations of numerical models to CPR primary measurements and uncertainties. The transformation accounts for sampling geometry and surface effects. We demonstrate Orbital-Radar's ability to provide realistic CPR views of typical cloud and precipitation scenes. The presented case studies show small-scale convection, marine stratus clouds, and Arctic mixed-phase cloud cases. These results provide valuable insights into the capabilities and challenges of the EarthCARE CPR mission and its advantages over the CloudSat CPR. Finally, Orbital-Radar allows for evaluating kilometre-scale numerical weather prediction models with EarthCARE CPR observations. So, Orbital-Radar can generate calibration and validation (Cal/Val) data sets already pre-launch. Nevertheless, an evaluation of synthetic CPR output data to accurate EarthCARE CPR data is missing.
1. Provide cleaned and filtered data stream from multiple ARM disdrometers; 2. Estimate commonly-used physical parameters of the Drop Size Distribution (DSD). 3. Calculate radar equivalent quantities such as the Radar Reflectivity Factor Z.
Abstract. This study extends the relative calibration adjustment technique for calibration of weather radars to higher-frequency radars as well as range–height indicator (RHI) scans. The calibration of weather radars represents one of the most dominant sources of error for their use in a variety of fields including quantitative precipitation estimation and model comparisons. While most weather radars are routinely calibrated, the frequency of calibration is often less than required, resulting in miscalibrated time periods. While full absolute calibration techniques often require the radar to be taken offline for a period of time, there have been online calibration techniques discussed in the literature. The relative calibration adjustment (RCA) technique uses the statistics of the ground clutter surrounding the radar as a monitoring source for the stability of calibration but has only been demonstrated to work at S- and C-band for plan-position indicator (PPI) scans at a constant elevation. In this work the RCA technique is modified to work with higher-frequency radars, including Ka-band cloud radars. At higher frequencies the properties of clutter can be much more variable. This work introduces an extended clutter selection procedure that incorporates the temporal stability of clutter and helps to improve the operational stability of RCA for relatively higher-frequency radars. The technique is also extended to utilize range–height scans from radars where the elevation is varied rather than the azimuth. These types of scans are often utilized with research radars to examine the vertical structure of clouds. The newly extended technique (eRCA) is applied to four Department of Energy Atmospheric Radiation Measurement (DOE ARM) weather radars ranging in frequency from C- to Ka-band. Cross comparisons of three co-located radars with frequencies C, X, and Ka at the ARM Cloud, Aerosol, and Complex Terrain Interactions (CACTI) site show that the technique can determine changes in calibration. Using an X-band radar at the ARM Eastern North Atlantic (ENA) site, we show how the technique can be modified to be more resilient to clutter fields that show increased variability, in this case due to sea clutter. The results show that this technique is promising for a posteriori data calibration and monitoring.
The U.S. Department of Energy (DOE) Atmospheric Radiation Measurement (ARM) Climate Research Facility’s radar organization was formed to provide an integrated approach to meet the scientific objectives of ARM. The radar organization consists of two groups: Radar Science and Radar Engineering and Operations. The two groups work closely for planning, coordination, and communications to produce high-quality observations and products from the ARM radar network, which consists of complex radars deployed in very challenging environments. In addition, scanning cloud radars are relatively new in the community, with very limited products that are mature enough to be applied for modeling studies. The Radar Science Group and Radar Engineering and Operations Group meet regularly along with researchers and users from the community in a workshop. The overarching objective of the workshop is to maximize the usage of ARM radar data for improving our understanding of processes and, ultimately improving representation of clouds and precipitation in models. The 6th DOE ARM Radar Engineering and Radar Science Workshop was held October 24-26, 2017, at Stony Brook University. Approximately 40 participants were present; a detailed list of participants is listed in Section 5. This report documents recommendations made by members of the science community who attended the workshop. Their recommendations appear italicized in the report. ARM management will consider these recommendations in their future planning and prioritization activities. The topics addressed during the workshop are contained in the agenda (Section 4). Significant action items, recommendations, and comments that came out of the meeting are contained in Section 3. The content in Section 3 presents a path through important discussions that were pursued during the Workshop.
The potential of millimeter-wavelength radar-based ice water content (IWC) estimation is demonstrated using a Ka-band Scanning Polarimetric Radar (KASPR) for the U.S. northeast coast winter storms. Two IWC relations for Ka-band polarimetric radar measurements are proposed: one that uses a combination of the radar reflectivity Z and the estimated total number concentration of snow particles N t and the other based on the joint use of Z, specific differential phase K DP , and the degree of riming f rim . A key element of the algorithms is to obtain the “Rayleigh-equivalent” value of Z measured at the Ka band, i.e., the corresponding Z at a longer radar wavelength for which Rayleigh scattering takes place. This is achieved via polarimetric retrieval of the mean volume diameter D m and incorporating the relationship between the dual-wavelength ratio DWR S/Ka and D m . Those techniques allow for retrievals from single millimeter-wavelength radar measurements and do not necessarily require the dual-wavelength ratio (DWR) measurements, if the DWR–D m relation and Rayleigh assumption for Ka-band K DP are valid. Comparison between the quasivertical profile product obtained from KASPR and the columnar vertical profile product generated from the nearby WSR-88D S-band radar measurements demonstrates that the DWR S/Ka can be estimated from the two close radars without the need for collocated radar beams and synchronized antenna scanning and can be used for determining the Rayleigh-equivalent value of Z. As a result, the performance of the suggested techniques is evaluated for seven winter storms using surface disdrometer and snow accumulation measurements.
The altitude distribution of meteors detected by a radar is sensitive to the instrument’s response function and can thus provide insight into the physical processes involved in radar measurements. This, in turn, can be used to determine the rate of ablation and ionization of the meteoroids and ultimately the input flux on Earth. In this work, we model the radar meteor head echo altitude distribution for three High Power and Large Aperture radar systems, by considering meteoroid populations from the main cometary family sources. In this simulation, we first use the results of a dynamical model of small meteoroids impacting Earth’s upper atmosphere to model the incoming mass, velocity, and entry angular distributions. We then combine these with the Chemical Ablation Model and establish the meteoroid ionization rates as a function of mass, velocity, and entry angle in order to determine the altitude at which these radars should detect the produced meteors and the portion of produced meteors from each population that are detected by these radars. We explore different sizes of head plasma as well as the possible effects on radar scattering of the head echo aspect sensitivity. We find that the modeled altitude distributions are generally in good agreement with measurements, particularly for ultra-high-frequency radars. In addition, our results indicate that the number of particles from Jupiter Family Comets (JFCs) required to fit the observations is lower than predicted by astronomical models. It is not clear yet if this discrepancy is due to the overprediction of JFC meteoroids by dynamical models or due to unaccounted physical processes in the treatment of ablation, ionization, and detections of meteoroids as they pass through Earth’s atmosphere.
The utilization of polarimetric weather radars for optimizing cloud models is a next frontier of research. It is widely understood that inadequacies in microphysical parameterization schemes in numerical weather prediction (NWP) models is a primary cause of forecast uncertainties. Due to its ability to distinguish between hydrometeors with different microphysical habits and to identify “polarimetric fingerprints” of various microphysical processes, polarimetric radar emerges as a primary source of needed information. There are two approaches to leverage this information for NWP models: (1) radar microphysical and thermodynamic retrievals and (2) forward radar operators for converting the model outputs into the fields of polarimetric radar variables. In this paper, we will provide an overview of both. Polarimetric measurements can be combined with cloud models of varying complexity, including ones with bulk and spectral bin microphysics, as well as simplified Lagrangian models focused on a particular microphysical process. Combining polarimetric measurements with cloud modeling can reveal the impact of important microphysical agents such as aerosols or supercooled cloud water invisible to the radar on cloud and precipitation formation. Some pertinent results obtained from models with spectral bin microphysics, including the Hebrew University cloud model (HUCM) and 1D models of melting hail and snow coupled with the NSSL forward radar operator, are illustrated in the paper.
Detailed observations from DOE ARM field campaigns were used to evaluate and improve microphysical and dynamic fields simulated by 3D Cloud Resolving Models (CRMs), employing both bulk and size-resolved bin microphysics. One of the main challenges in bridging the gap between observations and models is placing the analysis in a common framework. To address this, a polarimetric radar simulator was built, consisting of a forward operator to calculate model-consistent polarimetric radar fields, and an analysis package to calculate microphysical and dynamic fields from radar data in a self-consistent manner. The assumptions required for the transformation from model output to radar parameters were carefully considered and tested. It was found that although rain is fairly well captured by the forward model, ice characteristics are difficult to constrain. This polarimetric radar simulator, POLArimetric Radar Retrieval and Instrument Simulator (POLARRIS), was then applied to WRF simulations and polarimetric C-band radar observations from Darwin, Australia and the SGP site in Oklahoma. These two locations contrast maritime and mid-latitude continental convection respectively. POLARRIS was used with NASA Unified WRF simulations of the same storms in order to compare the observed land and ocean contrasts in microphysics and dynamics. An improved model framework for dynamic downscaling of large-scale aerosol products was also developed and applied to the NASA Unified WRF (NU-WRF) cloud resolving model. This improved framework was used to explore the effects of thermodynamic and aerosol invigoration in deep land and ocean convection through a series of sensitivity studies by perturbating aerosol (CCN) and convective available potential energy (CAPE). It was found that increases in aerosol concentrations yielded increases in low-density graupel in the maritime case and high-density hail in the continental case through riming by increased supercooled water. This result implies that native differences in aerosol concentrations may enhance the inherent contrast in riming effects between deep convective systems in the maritime and continental conditions. Formal publications were generated describing this research.
Abstract. A relative calibration technique has been developed for the US Department of Energy's (DOE's) Atmospheric Radiation Measurement (ARM) user facility Ka-band ARM Zenith Radars (KAZRs). This method uses the signal attenuation caused by water on the radome to estimate reflectivity factor (Ze) offsets. The wet-radome attenuation (WRA) is assumed to follow a log-linear relationship with rainfall rate during light and moderate rain, as measured by a collocated surface disdrometer. The technique has an uncertainty of approximately 3 dB, due to factors such as disdrometer measurement error, rain variability between radar and disdrometer sample volumes, and the fitting function's uncertainty for the WRA behavior. A practical advantage of this WRA-based approach to shorter-wavelength radar monitoring is that, while it requires a reference disdrometer, it proves feasible for a wider range of collocated disdrometer measurements compared to traditional direct disdrometer comparison at the onset of light rain. This technique thus offers a cost-effective monitoring tool for remote or long-term radar deployments. This calibration technique was applied during the ARM Tracking Aerosol Convection Interactions Experiment (TRACER) from October 2021 through September 2022. The estimated Ze offsets were compared against traditional radar calibration and monitoring methods using available datasets from this campaign. Results show that the WRA-based offsets align closely with mean offsets found between cloud radars and from direct disdrometer comparison near the onset of rain, while also reflecting similar offset and campaign-long trends when compared to collocated, independently calibrated radar wind profilers. Nevertheless, overall, the KAZR Ze offsets estimated during TRACER remained stable at approximately 2 dB lower than the disdrometer estimates from the campaign start until the end of June 2022; afterward, the offsets increased to around 7 dB by the campaign's end. This increase is linked to a drop of about 1 dB in transmitter power toward the end of the project.
Traditional fuzzy logic hydrometeor classification algorithm is a common way to classify precipitation type from dual polarization doppler radar. We propose a deep learning-based method to estimate hydrometeors efficiently using observed radar variables such as horizontal reflectivity (Z H ), differential reflectivity (Z DR ), correlation coefficient (ρ HV ) and specific differential phase (K DP ) from National Weather Service NEXRAD collected at Vance AFB facility at the first elevation angle from January 1st, 2015 to July 31th, 2019. We stack matrices of values from these four polarimetric variables as one 3D array. Samples are preprocessed and divided into training, validation and test set with four target hydrometeor categories (Ice Crystals (IC), Dry Snow (DS), Light and/or Moderate Rain (RA) and Big Drops (rain) (BD)). We developed and optimized five Convolutional Neural Networks (CNNs) architectures and achieved an accuracy of 87.23% and 93.736% respectively using modified ResNet with two different input data selection strategies for hydrometeor classification. Training data selection strategies were important to ensure use of available samples in training for robust performance evaluated by applying the models to novel time period beyond what was use to train the model. Seasonal variation in atmospheric conditions lead to seasonal patterns of liquid vs solid forms of precipitation, that poses challenge for classifier and offer insights into domain specific approaches required for problem of hydrometeor identification. Computationally efficient and scalable approach for classification of hydrometeors offer opportunties to effectively use the large volumes of rich time series of radar observations that are becoming increasingly available.
Clouds play an important role in Earth’s radiation budget and hydrological cycle. However, current global climate models (GCMs) have had difficulties in accurately simulating clouds and precipitation. To improve the representation of clouds in climate models, it is crucial to identify where simulated clouds differ from real world observations of them. This can be difficult, since significant differences exist between how a climate model represents clouds and what instruments observe, both in terms of spatial scale and the properties of the hydrometeors which are either modeled or observed. To address these issues and minimize impacts of instrument limitations, the concept of instrument “simulators”, which convert model variables into pseudo-instrument observations, has evolved with the goal to improve and to facilitate the comparison of modeled clouds with observations. Many simulators have (and continue to be developed) for a variety of instruments and purposes. A community satellite simulator package, the Cloud Feedback Model Intercomparison Project (CFMIP) Observation Simulator Package (COSP; Bodas-Salcedo et al. 2011), contains several independent satellite simulators and is being widely used in the global climate modeling community to exploit satellite observations for model cloud evaluation (e.g., Klein et al. 2013; Zhang et al. 2010). This article introduces a ground-based cloud radar simulator developed by the U.S. Department of Energy (DOE) Atmospheric Radiation Measurement (ARM) program for comparing climate model clouds with ARM observations from its vertically pointing 35-GHz radars. As compared to CloudSat radar observations, ARM radar measurements occur with higher temporal resolution and finer vertical resolution. This enables users to investigate more fully the detailed vertical structures within clouds, resolve thin clouds, and quantify the diurnal variability of clouds. Particularly, ARM radars are sensitive to low-level clouds, which are difficult for the CloudSat radar to detect due to surface contamination (Mace et al. 2007; Marchand et al. 2008). Therefore, the ARM ground-based cloud observations can provide important observations of clouds that complement measurements from space.
Abstract. The Energy Exascale Earth System Model (E3SM) developed by the Department of Energy has a goal of addressing challenges in understanding the global water cycle. Success depends on correct simulation of cloud and precipitation elements. However, lack of appropriate evaluation metrics has hindered the accurate representation of these elements in general circulation models. We derive metrics from the three-dimensional data of the ground-based Next-Generation Radar (NEXRAD) network over the US to evaluate both horizontal and vertical structures of precipitation elements. We coarsened the resolution of the radar observations to be consistent with the model resolution and improved the coupling of the Cloud Feedback Model Intercomparison Project Observation Simulator Package (COSP) and E3SM Atmospheric Model Version 1 (EAMv1) to obtain the best possible model output for comparison with the observations. Three warm seasons (2014–2016) of EAMv1 simulations of 3-D radar reflectivity features at an hourly scale are evaluated. A general agreement in domain-mean radar reflectivity intensity is found between EAMv1 and NEXRAD below 4 km altitude; however, the model underestimates reflectivity over the central US, which suggests that the model does not capture the mesoscale convective systems that produce much of the precipitation in that region. The shape of the model-estimated histogram of subgrid-scale reflectivity is improved by correcting the microphysical assumptions in COSP. Different from previous studies that evaluated modeled cloud top height, we find the model severely underestimates radar reflectivity at upper levels – the simulated echo top height is about 5 km lower than in observations – and this result is not changed by tuning any single physics parameter. For more accurate model evaluation, a higher-order consistency between the COSP and the host model is warranted in future studies.
This study explores the updraft and downdraft properties of mature stage mesoscale convective systems (MCSs) in terms of draft core width, shape, intensity, and mass flux characteristics. The observations use extended radar wind profiler (RWP) and surveillance radar datasets from the U.S. Department of Energy Atmospheric Radiation Measurement program for midlatitude (Oklahoma, USA) and tropical (Amazon, Brazil) sites. MCS drafts behave qualitatively similar to previous aircraft and RWP cloud summaries. The Oklahoma MCSs indicate larger and more intense convective up- and downdraft cores, and greater mass flux than Amazon MCS counterparts. However, similar size-intensity relationships and draft vertical profile behaviors are observed for both regions. Additional similarities include weak positive correlations between core intensity and core width (correlation coefficient r ~ 0.5), and increases in draft intensity with altitude. A model-observational intercomparison for draft properties (core width, intensity, mass flux) is also performed to illustrate the potential usefulness of statistical observed draft characterizations. Idealized simulations with the Weather Research and Forecasting model aligned with midlatitude MCS conditions are performed at model grid spacings (?x) that range from 4 km to 250 m. It is shown that the simulations performed at ?x = 250 m at similar mature MCS lifecycle stages are those that exhibit draft intensity, width, mass flux, and shape parameter performances best matching with observed properties.
The Next Generation Weather Radar (NEXRAD) system is a network of doppler radar operated jointly by the National Weather Service (NWS), the Federal Aviation Administration (FAA), and the U.S. Air Force. This dataset represents a line-of-sight for each radar station. Radar line-of-sight can become saturated with wind energy due to interference. For further details and citation, please refer to the publication linked below: Lopez, Anthony, Pavlo Pinchuk, Michael Gleason, Wesley Cole, Trieu Mai, Travis Williams, Owen Roberts, Marie Rivers, Mike Bannister, Sophie-Min Thomson, Gabe Zuckerman, and Brian Sergi. 2024. Solar Photovoltaics and Land-Based Wind Technical Potential and Supply Curves for the Contiguous United States: 2023 Edition. Golden, CO: National Renewable Energy Laboratory. NREL/TP-6A20-87843.
The Python ARM Radar Toolkit is a package for reading, visualizing, correcting and analysing data from weather radars. Development began to meet the needs of the Atmospheric Radiation Measurement Climate Research Facility and has since expanded to provide a general-purpose framework for working with data from weather radars in the Python programming language. The toolkit is built on top of libraries in the Scientific Python ecosystem including NumPy, SciPy, and matplotlib, and makes use of Cython for interfacing with existing radar libraries written in C and to speed up computationally demanding algorithms. As a result, the source code for the toolkit is available on GitHub and is distributed under a BSD license.
Corrected Moments in Antenna Coordinates (CMAC) calculates quantitative precipitation estimates (QPE) from empirical relationships based on equivalent radar reflectivity factor, specific differ- ential phase, and specific attenuation. To evaluate these empirical relationships, the Extracted Radar Columns and In-Situ Sensors (RADclss) product was developed. Utilizing Py-ART, RAD- clss extracts CMAC radar columns above various ARM and partner locations. These columns are then spatiotemporally synced with in-situ observations at the surface utilizing the Atmospheric data Community Toolkit (ACT; Theisen et al. 2025), allowing direct comparison of radar parameters with rain gauges and laser disdrometers for further investigation.
This data set contains multi-Doppler radar analysis from CSAPR-2, CSU-CHIVO, COW, and RMA-1 radars. These radars were collecting dual-polarization data during the CACTI/RELAMPAGO experiments in Argentina in 2018. Doppler analysis is systematically conducted for 31 days with convection. Three-dimensional wind fields are retrieved using the PyDDA algorithm. Dual-polarization information is also included in the data set.