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.
Seven years of autumnal ground-based observations (September–November 2002–2008) at the U.S. Department of Energy Atmospheric Radiation Measurement North Slope of Alaska site have been analyzed for addressing the occurrence frequency and macrophysical and microphysical properties of Arctic mixed-phase clouds (AMC), as well as the relationship between environmental parameters and AMC properties. In September and October, AMC occurrence frequency is 20–30% lower during a southerly wind when compared to the other wind directions; in November, the variation of AMC occurrence frequency with wind direction is small. The mean liquid water path in November is about half of that in October and September. When the surface is snow free, temperature ( T ) and specific humidity ( q ) profiles during a northerly wind are warmer and moister than those for the southerly wind. Northerly wind profiles have a higher relative humidity to ice (RH i ) and lower atmosphere stability. Furthermore, the AMC occurrence frequency has a positive relationship with RH i and a negative relationship with stability. These two points may explain the lower AMC occurrence frequency during a southerly wind. During a northerly wind, AMCs have larger radar reflectivity, wider spectrum width, and larger Doppler velocity signatures. The stronger precipitation for AMC during a northerly wind is possibly due to the cleaner air masses from the ocean (north). Lastly, with the same amount of q, the radar spectrum width has a higher frequency in the larger bins during a northerly wind. Both T, q , and radar reflectivity, radar spectrum width profiles show evidence of deposition in the sub-cloud layer in September and October.
Cloud radars are unique instruments for observing cloud processes, but uncertainties in radar calibration have frequently limited data quality. Thus far, no single robust method exists for assessing the calibration of past cloud radar data sets. Here, we investigate whether observations of microphysical processes in liquid clouds such as the transition of cloud droplets to drizzle drops can be used to calibrate cloud radars. Specifically, we study the relationships between the radar reflectivity factor and three variables not affected by absolute radar calibration: the skewness of the radar Doppler spectrum ( γ ), the radar mean Doppler velocity ( W ), and the liquid water path (LWP). For each relation, we evaluate the potential for radar calibration. For γ and W , we use box model simulations to determine typical radar reflectivity values for reference points. We apply the new methods to observations at the Atmospheric Radiation Measurement (ARM) sites North Slope of Alaska (NSA) and Oliktok Point (OLI) in 2016 using two 35 GHz K a -band ARM Zenith Radars (KAZR). For periods with a sufficient number of liquid cloud observations, we find that liquid cloud processes are robust enough for cloud radar calibration, with the LWP-based method performing best. We estimate that, in 2016, the radar reflectivity at NSA was about 1±1  dB too low but stable. For OLI, we identify serious problems with maintaining an accurate calibration including a sudden decrease of 5 to 7 dB in June 2016.
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.
Downdrafts and cold pool characteristics for strong mesoscale convective systems (MCSs) and isolated, unorganized deep precipitating convection are analyzed using multi-instrument data from the DOE Atmospheric Radiation Measurement (ARM) GoAmazon2014/5 campaign. Increases in column water vapor (CWV) are observed leading convection, with higher CWV preceding MCSs than for isolated cells. For both MCSs and isolated cells, increases in wind speed, decreases in surface moisture and temperature, and increases in relative humidity occur coincidentally with system passages. Composites of vertical velocity data and radar reflectivity from a radar wind profiler show that the downdrafts associated with the sharpest decreases in surface equivalent potential temperature ( θ e ) have a probability of occurrence that increases with decreasing height below the freezing level. Both MCSs and unorganized convection show similar mean downdraft magnitudes and probabilities with height. Mixing computations suggest that, on average, air originating at heights greater than 3 km must undergo substantial mixing, particularly in the case of isolated cells, to match the observed cold pool θ e , implying a low typical origin level. Precipitation conditionally averaged on decreases in surface equivalent potential temperature (Δ θ e ) exhibits a strong relationship because the most negative Δ θ e values are associated with a high probability of precipitation. The more physically motivated conditional average of Δ θ e on precipitation shows that decreases in θ e level off with increasing precipitation rate, bounded by the maximum difference between surface θ e and its minimum in the profile aloft. In conclusion, robustness of these statistics observed across scales and regions suggests their potential use as model diagnostic tools for the improvement of downdraft parameterizations in climate models.
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.
Two years of volumetric radar observations acquired during the GoAmazon2014/5 field campaign by the SIPAM S-band radar in Manaus, Brazil were used to create rain rate estimates from radar reflectivity by Texas A&M University. Estimates of rain rates were generated using radar reflectivity data from the 2.5 km level and a Z-R relationship of Z=174.8*R^1.56.
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.
Abstract Accurate measurements of snow amounts by radar are very difficult to achieve. The inherent uncertainty in radar snow estimates that are based on the radar reflectivity factor Z is caused by the variability of snow particle size distributions and snow particle density as well as the large diversity among snow growth habits. In this study, a novel method for snow quantification that is based on the joint use of radar reflectivity Z and specific differential phase KDP is introduced. An extensive dataset of 2D-video-disdrometer measurements of snow in central Oklahoma is used to derive polarimetric relations for liquid-equivalent snowfall rate S and ice water content IWC in the forms of bivariate power-law relations S = and along with similar relations for the intercept N0s and slope Λs of the exponential snow size distribution. The physical basis of these relations is explained. Their multipliers are sensitive to variations in the width of the canting angle distribution and to a lesser extent the particles’ aspect ratios and densities, whereas the exponents are practically invariant. This novel approach is tested against the S(Z) relation using snow disdrometer measurements in three geographical regions (Oklahoma, Colorado, and Canada). Significant improvement in snow estimates relative to the traditional Z-based methods is demonstrated.
Here, the enhancement of precipitation processes aloft over complex terrain is documented using reflectivity data from an S–band scanning radar (NPOL) that was deployed on the west coast of Washington State during the Olympic Mountains Experiment (OLYMPEX). From November 2015 through mid–January 2016, NPOL obtained high–resolution data within sectors over the ocean and over the windward slopes of the Olympic Mountains. Contoured Frequency by Altitude Diagrams of radar reflectivity highlight a higher frequency of occurrence of larger reflectivities for all heights between 2 and 8 km over land compared to ocean, with the largest difference in the 4– to 6–km range indicating a robust signature of enhancement aloft over the windward slopes. This enhancement pattern is found to some degree under all environmental conditions considered but is especially pronounced during periods of high vapor transport, high melting level height, southwest low–level winds, and neutral stability. These conditions are generally associated with warm sectors of midlatitude cyclones and atmospheric rivers. Past studies have postulated that a secondary enhancement in reflectivity aloft was an intrinsic part of atmospheric river type systems. However, these results show that further significant enhancement of this signature occurs as deep moist–neutral, high water vapor content flow is lifted when it encounters a mountain range. Reflectivity data from the dual–precipitation radar aboard the Global Precipitation Measurement satellite also documents this reflectivity increase aloft over the Olympic Mountains compared to the adjacent ocean, showing the potential for Global Precipitation Measurement to provide reliable estimates of precipitation structure over remote mountainous regions.