SEARCH · Search NASA
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.
Daily files containing radar reflectivity and mean Doppler velocity collected from vertically pointing W-band radar during the Macquarie Island Cloud & Radiation Experiment (MICRE): April 2016 to March 2017.
The retrieval of hail kinetic energy with weather radars or its simulation in numerical models is challenging because of the shape complexity and variable density of hailstones. We combine 3D scans of individual hailstones with measurements of the particle size distributions (PSD) and T-matrix calculations to understand how hail reflectivity Z changes when approximating hailstones as spheroids, as compared to the realistic shapes obtained by 3D scanning technology. Additionally, recent terminal velocity relations are used to compare Z to the hail kinetic energy flux E ˙ . We parameterize the hail backscattering cross sections at L, S, C, and X bands as a function of size between 0.5 and 5.0 cm, matching the range of the observed PSDs. The scattering calculations use the T-matrix method for size parameters below 1.0 and the discrete dipole approximation (DDA) method otherwise. The DDA calculations are done for 48 digital models of realistic hailstones of sizes between 1 and 5 cm. The DDA cross sections are calculated for multiple orientations and averaged assuming a fully random orientation distribution to provide a single value per hailstone. The T-matrix reflectivity assuming solid ice spheres presents negligible differences to DDA results for size parameters below 1.0. Therefore, T matrix was used to fill in the gaps left by the DDA calculations. The results are mapped to the same size bins of the observed PSDs, allowing the calculation of the radar reflectivity. This is then correlated to E ˙ , allowing a potential improvement of past retrieval methods of E ˙ from Z in multiple wavelengths.
The Echo Classification from COnvectivity (ECCO) algorithm identifies convective and stratiform types of radar echo in three dimensions. It is based on the calculation of reflectivity texture—a combination of the intensity and the heterogeneity of the radar echoes on each horizontal plane in a 3D Cartesian volume. Reflectivity texture is translated into convectivity, which is designed to be a quantitative measure of the convective nature of each 3D radar grid point. It ranges from 0 (100% stratiform) to 1 (100% convective). By thresholding convectivity, a more traditional qualitative categorization is obtained, which classifies radar echoes as convective, mixed, or stratiform. In contrast to previous algorithms, these echo-type classifications are provided on the full 3D grid of the reflectivity field. The vertically resolved classifications, in combination with temperature data, allow for subclassifications into shallow, mid-, deep, and elevated convective features, and low, mid-, and high stratiform regions—again in three dimensions. The algorithm was validated using datasets collected over the U.S. Great Plains during the PECAN field campaign. An analysis of lightning counts shows ~90% of lightning occurring in regions classified as convective by ECCO. A statistical comparison of ECCO echo types with the well-established GPM radar precipitation-type categories show 84% (88%) of GPM stratiform (convective) echo being classified as stratiform (convective) or mixed by ECCO. ECCO was applied to radar grids for the continental United States, the United Arab Emirates, Australia, and Europe, illustrating its robustness and adaptability to different radar grid characteristics and climatic regions.
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.
Abstract. The detection of the early growth of drizzle particles in marine stratocumulus clouds is important for studying the transition from cloud water to rainwater. Radar reflectivity is commonly used to detect drizzle; however, its utility is limited to larger drizzle particles. Alternatively, radar Doppler spectrum skewness has proven to be a more sensitive quantity for the detection of drizzle embryos. Here, a machine learning (ML)-based technique that uses radar reflectivity and skewness for detecting small drizzle particles is presented. Aircraft in situ measurements are used to develop and validate the ML algorithm. The drizzle detection algorithm is applied to three Atmospheric Radiation Measurement (ARM) observational campaigns to investigate the drizzle occurrence in marine boundary layer clouds. It is found that drizzle is far more ubiquitous than previously thought; the traditional radar-reflectivity-based approach significantly underestimates the drizzle occurrence, especially in thin clouds with liquid water paths lower than 50 g m−2. Furthermore, the drizzle occurrence in marine boundary layer clouds differs among the three ARM campaigns, indicating that the drizzle formation, which is controlled by the microphysical process, is regime dependent. A complete understanding of the drizzle distribution climatology in marine stratocumulus clouds calls for more observational campaigns and continuing investigations.
In this study, we investigate how the regional climate model HIRHAM5 reproduces the spatial and temporal distribution of Arctic snowfall when compared to CloudSat satellite observations during the examined period of 2007–2010. For this purpose, both approaches, i.e., the assessments of the surface snowfall rate (observation-to-model) and the radar reflectivity factor profiles (model-to-observation), are carried out considering spatial and temporal sampling differences. The HIRHAM5 model, which is constrained in its synoptic representation by nudging to ERA-Interim, represents the snowfall in the Arctic region well in comparison to CloudSat products. The spatial distribution of the snowfall patterns is similar in both identifying the southeastern coast of Greenland and the North Atlantic corridor as regions gaining more than twice as much snowfall as the Arctic average, defined here for latitudes between 66 and 81°N. Excellent agreement (difference less than 1%) in the Arctic-averaged annual snowfall rate between HIRHAM5 and CloudSat is found, whereas ERA-Interim reanalysis shows an underestimation of 45% and significant deficits in the representation of the snowfall rate distribution. From the spatial analysis, it can be seen that the largest differences in the mean annual snowfall rates are an overestimation near the coastlines of Greenland and other regions with large orographic variations as well as an underestimation in the northern North Atlantic Ocean. To a large extent, the differences can be explained by clutter contamination, blind zone or higher resolution of CloudSat measurements, but clearly HIRHAM5 overestimates the orographic-driven precipitation. The underestimation of HIRHAM5 within the North Atlantic corridor south of Svalbard is likely connected to a poor description of the marine cold air outbreaks which could be identified by separating snowfall into different circulation weather type regimes. By simulating the radar reflectivity factor profiles from HIRHAM5 utilizing the Passive and Active Microwave TRAnsfer (PAMTRA) forward-modeling operator, the contribution of individual hydrometeor types can be assessed. Looking at a latitude band at 72–73°N, snow can be identified as the hydrometeor type dominating radar reflectivity factor values across all seasons. The largest differences between the observed and simulated reflectivity factor values are related to the contribution of cloud ice particles, which is underestimated in the model, most likely due to the small sizes of the particles. The model-to-observation approach offers a promising diagnostic when improving cloud schemes, as illustrated by comparison of different schemes available for HIRHAM5.
We characterize convective clouds associated with sea‐breeze circulations (SBC) using multi‐agency observations and multi‐case ensemble model simulations. The focus is on assessing convective cloud lifecycle properties and their merging behavior, as well as the environmental conditions they are embedded in, particularly SBC features. In total, 46 SBC days over the Houston‐Galveston region are selected and simulated using the Weather Research and Forecasting (WRF) model at a gray zone scale with a forecast‐like parameterization setup. Advanced techniques, including change‐point detection, a Lagrangian cloud tracking method, and a newly developed cell merging and splitting detection algorithm, are applied and/or developed for this study. Our findings indicate that the WRF model at 1 km grid spacing well represents the thermodynamic conditions over the region, as well as SBC timing and intensity. However, for the associated convective cells, WRF overestimates the 30‐dBZ echo top height, cell area, and maximum radar reflectivity compared to radar observations. This overestimation is potentially due to under‐resolved entrainment processes, an overestimated merging frequency, and the overestimation of updraft intensity. Furthermore, the model exhibits a deficiency in simulating congestus clouds, showing a more rapid transition from shallow to deep convection compared to observed behavior. Moreover, observations indicate stronger, deeper, and wider clouds when merging happens. Conversely, in simulations, the merging process does not necessarily lead to higher or longer‐lived cells, as many cases experience rapid and frequent merging and splitting which may result in more variance in convective updraft velocity during the convection lifetime.
Vertically pointing radars (VPRs) operating at millimeter wavelengths measure the power return from raindrops enabling precipitation retrievals as a function of height. However, as the rain rate increases, there are combinations of rain rate and rain path length that produce sufficient attenuation to prevent the radar from detecting raindrops all the way through rain shafts. This study explores the question: Which rain rate and path length combinations completely extinguish radar return signals for VPRs operating between 3 and 200 GHz? An important step in these simulations is converting attenuated radar reflectivity factor into radar received signal-to-noise ratio (SNR) in order to determine the range where the SNR drops below the receiver detection threshold. Configuring the simulations to mimic a U.S. Department of Energy Atmospheric Radiation Mission (ARM) W-band (95 GHz) radar deployed in Brazil, the simulation results indicate that a W-band radar could observe raindrops above 3.5 km only when the rain rate was less than approximately 4 mm h -1 . The deployed W-band radar measurements confirm the simulation results with maximum observed heights ranging between 3 and 4.5 km when a surface disdrometer measured 4 mm h -1 rain rate (based on 25-to-75 percentiles from over 25,000 W-band radar profiles). In summary, this study contributes to our understanding of how rain and atmospheric gas attenuation impacts the performance of millimeter-wave VPRs and will help with the design and configuration of multi-frequency VPRs deployed in future field campaigns.
The profile of the liquid water content (LWC) in clouds provides fundamental information for understanding the internal structure of clouds, their radiative effects, propensity to precipitate, and degree of entrainment and mixing with the surrounding environment. In principle, differential absorption techniques based on coincident dual-frequency radar reflectivity observations have the potential to provide the LWC profile. Previous differential frequency radar reflectivity (DFR) efforts were challenged by the fact that the measurable differential attenuation for small quantities of LWC is usually comparable to the system measurement error. This typically renders the retrieval impractical, as the uncertainty can become many times greater than the retrieved value itself. Theoretically, this drawback can be mitigated following two interconnected approaches: (1) increasing the frequency separation between the dual-frequency radar system to measure greater differential attenuation and (2) increasing the radar operating frequency to reduce the instrument measurement random error. Our recently developed 239 GHz radar was deployed during the Eastern Pacific Cloud Aerosol Precipitation Experiment (EPCAPE) along with a variety of collocated remote sensing and in situ instruments. We have combined Ka-band (35 GHz) and G-band (239 GHz) observations to retrieve the LWC from more than 100 vertical profiles of shallow clouds with typical amounts of LWC smaller than 1 g m -3 . We theoretically and experimentally demonstrate that the Ka-band and G-band pair of frequencies offers at least a 65 % relative improvement in the LWC retrieval sensitivity compared to previous works reported in the literature using lower-frequency radars. This new technique provides a missing capability to determine the LWC in the challenging low liquid water path (LWP) range (< 200 g m -2 ) and suggests a way forward to characterize microphysical and dynamical processes more precisely in shallow clouds.
Over the years rainfall data for the Savannah River Site has been obtained from ground level measurements made by rain gauges. These instruments have inherent errors or biases that can impact the measured rainfall totals but are assumed as ground truth for most climatological and weather applications. With the development of weather radar technologies, various methods to derive rainfall totals from radar reflectivity values have been developed and have continued to improve. The Z-R relationship, which uses an exponential relationship to estimate rainfall rate based on radar reflectivity values, provides estimates of rainfall amounts (Z-R Level III) for locations within the radar’s detection range. The recently developed Multi-Radar Multi-Sensor (MRMS) dataset combines reflectivity-based estimates using the Z-R relationship with a network of gauges and other rainfall estimates to produce a refined set of precipitation estimates for each grid point within its domain. Comparisons were done between gauge observations and radar estimates for various SRS locations to assess whether radar derived estimates are representative of rainfall measurements at the site. Results obtained show good agreement between radar derived amounts and ground measurements, with MRMS showing stronger correlations and lower spread than Z-R Level III estimates. Outliers and errors observed appear to be related to hydrometeor classification schemes. In general, MRMS proved to provide sufficiently representative estimates of daily rainfall amounts to replace existing SRS rain gauge measurements.
Following methods outlined in Remillard et al. (2012), we classify seven cloud types using radar reflectivity, best-estimated cloud base, and cloud-layer product from the Active Remotely Sensed Cloud Locations (ARSCL) product (Kollias et al. 2007). A cloud mask is created based on the detectable radar reflectivity (>-40 dBZ) combined with the best-estimated cloud base height. Each cloud object is analyzed individually as contiguous cloudy pixels, and its type is determined based on the cloud’s boundaries and duration. Focusing on marine boundary-layer clouds, low clouds are further classified into four types: shallow cumulus, broken stratocumulus (Sc) or cumulus clouds, single-layer Sc, and multi-layer Sc or Sc coupled with cumulus. The remaining three categories are middle clouds, high clouds, and deep convective clouds. For detailed definition of each cloud type, please refer to Zheng et al. (2024).
Shallow cloud decks residing in or near the boundary layer cover a large fraction of the Southern Ocean (SO) and play a major role in determining the amount of shortwave radiation reflected back to space from this region. In this article, we examine the macrophysical characteristics and thermodynamic phase of low clouds (tops <3 km) and precipitation using ground-based ceilometer, depolarization lidar and vertically-pointing W-band radar measurements collected during the Macquarie Island Cloud and Radiation Experiment (MICRE) from April 2016 to March 2017. During MICRE, low clouds occurred ~65% of the time on average (slightly more often in austral winter than summer). About 2/3 of low clouds were cold-topped (temperatures ≤0°C). These were thicker and had higher bases on average than warm-topped clouds. 83%–88% of cold-topped low clouds were liquid phase at cloud base (depending on the season). The majority of low clouds had precipitation in the vertical range 150–250 m below cloud base, a significant fraction of which did not reach the surface. Phase characterization is limited to the period between April 2016 and November 2016. Small-particle (low-radar-reflectivity) precipitation (which dominates precipitation occurrence) was mostly liquid below-cloud, while large-particle precipitation (which dominates total accumulation) was predominantly mixed/ambiguous or ice phase. Approximately 40% of cold-topped clouds had mixed/ambiguous or ice phase precipitation below (with predominantly liquid phase cloud droplets at cloud base). Below-cloud precipitation with radar reflectivity factors below about -10 dBZ were predominantly liquid, while reflectivity factors above about 0 dBZ were predominantly ice.
In this paper, a novel methodology for modeling ice-ice aggregation is presented. This methodology combines a modified hydrodynamic collection algorithm with bulk aggregate characteristic information from an offline simulator that collects ice particles, namely, the Ice Particle and Aggregate Simulator, and has been implemented into the Adaptive Habit Microphysics scheme in the Weather Research and Forecasting Model. Aggregates, or snow, are formed via collection of cloud ice particles, where initial ice characteristics and the resulting geometry determine aggregate characteristics. Upon implementation, idealized squall-line simulations are performed to examine the new methodology in comparison with commonly used bulk microphysics schemes. It is found that the adaptive habit aggregation parameterization develops snow and reduces ice mass and number concentrations compared to other schemes. The development of aggregates through the new methodology cascades into other interesting effects, including enhancements in ice and snow growth, as well as homogeneous freezing. Further microphysical analyses reveal varying sensitivities, where snow processes are most sensitive to the new parameterization, followed by ice, then cloud, rain, and graupel processes. Further, the new scheme results in enhancements in surface precipitation due to the persistence of snow at lower altitudes. This persistence is a result of shape-dependent melting and sublimation, increasing the residence time. Moreover, these low-level enhancements are reflected in increases in radar reflectivity at the surface and its spatial distribution. Finally, the ability to predict snow shape and density allows for the simulation of polarimetric radar quantities, resulting in signature enhancements compared to schemes that do not consider spatial and temporal variations in snow shape and density.
Explore the source record for details and available documents.
The SGP Central Facility (C1) radar wind profiler (RWP) was calibrated using nearby surface disdrometer observations. Between 2011 and 2019, the SGP C1 RWP operated in two modes. The vertically pointing mode (named the precipitation mode) transmitted a short and long pulse length to have two different range resolutions and the beam-swinging mode (named the wind mode) transmitted one pulse length into three different beam directions. The precipitation-mode observations were available and calibrated from April 2011 through mid-August 2019. The wind-mode observations were available and calibrated between April 2014 and March 2019. The RWP spectra were processed to account for Nyquist velocity aliasing and coherent integration filtering effects before calculating the spectrum moments. During intense precipitation events, the calculated signal-to-noise ratio (SNR) is biased low due to signal power being distributed across the velocity spectrum such that some signal power is erroneously included in the noise level estimate, causing the noise level to be biased high. To correct for the low SNR bias, a new noise level is estimated using observations without precipitation and the SNR is increased accordingly. The adjusted SNR was converted to radar reflectivity factor and then calibrated against a nearby surface disdrometer. The calibration methodology is fully described in: Williams, CR, J Barrio, PE Johnston, P Muradyan, and S Giangrande. 2023. “Calibrating radar wind profiler reflectivity factor using surface disdrometer observations.” Atmospheric Measurement Techniques, https://egusphere.copernicus.org/preprints/2023/egusphere-2022-1405/
The instrument is a polarimetric W-Band simultaneous transmission simultaneous reception (STSR) Doppler Cloud Radar manufactured by Radiometer Physics GmbH (RPG, instrument type RPG-FMCW-94-DP). It was deployed about 500m away from the ARM M1 facility in Gothic (CO) near the Ore House between Nov 15, 2022 and Jun 5, 2023 at an altitude of 2913 m. The Doppler cloud radar was mounted on a prototype of a cold temperature scanner manufactured by RPG. The dataset contains common radar variables like equivalent radar reflectivity factor, mean Doppler velocity, Doppler spectrum width as well as common polarimetric radar variables like differential phase shift (PhiDP), specific differential phase shift (KDP), correlation coefficient (RHV) and differential reflectivity (ZDR). Between Nov 15, 2022, and Feb 7, 2023 the radar was operated at 40° constant elevation and 151° azimuth (towards the KAZR at the AMF1). Between Feb 7 and Mar 15, 2023 the following scan cycle was repeated in 15min intervals: zenith pointing measurements, RHI scans from 0° to 90° and 90° to 0° elevation at 151° azimuth, constant elevation measurements at 40° elevation and 151° azimuth and PPI scans at 85° elevation for wind retrievals (one PPI per hour). Between Mar 15 18:18 UTC and Apr 8, 2023 11:04 UTC, the scanner was stuck at 72° elevation, which resulted in corrupted files for this period. For that period, the files containing data are uploaded but should be treated with care. Between Apr 8 and Apr 14, 2023 the scanning pattern described above was performed again. Afterwards, except for May 17 and May 18, 2023 where scanning patterns as described above were performed, measurements continued at 40° elevation and 151° azimuth until May 22, 2023 09:00 UTC. The scanning pattern was then started again, but without hourly PPI measurements until May 31, 2023. Between May 31 and Jun 5, 2023 solely zenith-pointing observations were performed. The data provided here are the RPG Level 1 (L1) files, Doppler spectra Level 0 (L0) files are available upon request.
ARM's scanning cloud radars are fully coherent dual-frequency, dual-polarization Doppler radars mounted on a common scanning pedestal. Each pedestal includes a Ka-band radar (2kW peak power) and the deployment location determines whether the second radar is a W-band (WSACR; 1.7 kW peak power) or X-band (XSACR; 20 kW peak power). Beamwidths at Ka and W bands are roughly matched at 0.3 degrees. Due to the narrow antenna beamwidth, ARM’s scanning cloud radars use scanning strategies that are unlike typical weather radars. Rather than focusing on plan position indicator, or PPI, scans, the Ka-SACR uses range height indicator, or RHI, scans at numerous azimuths to obtain cloud volume data. Measurements collected with the W-SACR are copolar and cross-polar radar reflectivity, Doppler velocity, spectra width and spectra when not scanning, and linear depolarization ration.
ARM's scanning cloud radars are fully coherent dual-frequency, dual-polarization Doppler radars mounted on a common scanning pedestal. Each pedestal includes a Ka-band radar (2kW peak power) and the deployment location determines whether the second radar is a W-band (WSACR; 1.7 kW peak power) or X-band (XSACR; 20 kW peak power). Beamwidths at Ka and W bands are roughly matched at 0.3 degrees. Due to the narrow antenna beamwidth, ARM’s scanning cloud radars use scanning strategies that are unlike typical weather radars. Rather than focusing on plan position indicator, or PPI, scans, the Ka-SACR uses range height indicator, or RHI, scans at numerous azimuths to obtain cloud volume data. Measurements collected with the W-SACR are copolar and cross-polar radar reflectivity, Doppler velocity, spectra width and spectra when not scanning, and linear depolarization ration.