Search

Snow provides critical water resources for billions of people, making the remote sensing of snow water equivalent (SWE) a highly prioritized endeavor, particularly given ongoing climate change impacts. Synthetic aperture radar (SAR) is a promising method for remote sensing of SWE because radar penetrates snow, and SAR interferometry (InSAR) can be used to estimate changes in SWE (Δ SWE) between SAR acquisitions. We calculated Δ SWE retrievals from 10 NASA L-band (1–2 GHz, ∼25 cm wavelength) uninhabited aerial vehicle SAR (UAVSAR) acquisitions covering a ∼640 km 2 swath in northern Colorado during the winters of 2020 and 2021. UAVSAR acquisitions coincided with ∼117 mm of accumulation in 2020 and ∼282 mm of accumulation in 2021. Δ SWE retrievals were evaluated against measurements of SWE from repeat ground-penetrating radar (GPR) and terrestrial lidar scans (TLSs) collected during the NASA SnowEx time series campaigns at two field sites (total area =∼0.2 km 2) as well as SWE measurements from seven automated stations distributed throughout the UAVSAR swath. For single InSAR pairs, UAVSAR Δ SWE retrievals yielded an overall r of 0.72–0.79 and an RMSE of 19–22 mm when compared with TLS and GPR Δ SWE retrievals. UAVSAR Δ SWE showed some scatter with Δ SWE measured at automated stations for both study years, but cumulative UAVSAR SWE yielded a r of 0.92 and an RMSE of 42 mm when compared to total SWE measured by the stations. Further, UAVSAR Δ SWE RMSEs differed by <10 mm for coherences (i.e., the complex interferometric coherence) of 0.10 to 0.90, suggesting that coherence has only a small influence on the Δ SWE retrieval accuracy. Given the evaluations presented here and in other recent studies, the upcoming NASA-ISRO SAR (NISAR) satellite mission, with a 12 d revisit period, offers an exciting opportunity to apply this methodology globally. [ABSTRACT FROM AUTHOR]

Active-layer thickness (ALT) is estimated for a study area in northern Alaska's continuous-permafrost zone using satellite data from Sentinel-1 (radar) and ICESat-2 (lidar) for the period 2017 to 2022. Synthetic aperture radar (SAR) interferograms were generated using the Short Baseline Subset (SBAS) approach. Displacement time series over the thaw season (June–September) are fit well with a linear model (root mean square error (RMSE) scatter is less than 7 mm) and show maximum seasonal subsidence of 20–60 mm. ICESat-2 products were used to validate the interferometric synthetic aperture radar (InSAR) displacement time series. ALT was estimated from measured subsidence using a widely used model exploiting the volume difference between ice and water, reaching a maximum depth in our study area of 1.5 m. Estimated ALT is in good agreement with in situ and other remotely sensed data but is sensitive to assumed thaw season onset, indicating the need for reliable surface temperature data. Our results suggest the feasibility of long-term permafrost monitoring with satellite InSAR. However, the C-band (∼55 mm center wavelength) Sentinel radar is sensitive to vegetation cover and, in our studies, was not successful for similar monitoring in the heavily treed discontinuous-permafrost zone of central Alaska. [ABSTRACT FROM AUTHOR]

Remote hail detection and hail size estimation using weather radar observations has the advantage of wide spatial coverage and high spatial and temporal resolution. Switzerland's National Weather Service (MeteoSwiss) uses two radar-based hail metrics: the probability of hail on the ground (POH) to assess the presence of hail and the maximum expected severe hailstone size (MESHS) to estimate the largest hailstone diameter. However, radar-based metrics are not direct measurements of hail and have to be calibrated with and verified against ground-based observations of hail, such as crowdsourced hail reports. Switzerland benefits from a particularly rich and dense dataset of crowdsourced hail reports from the MeteoSwiss app. We combine a new spatiotemporal clustering method (Density-Based Spatial Clustering of Applications with Noise, ST-DBSCAN) with radar reflectivity to filter the reports and use the filtered reports to verify POH and MESHS in terms of the hit rate, false-alarm ratio (FAR), critical success index (CSI), and Heidke skill score (HSS). Using a 4 km × 4 km maximum upscaling approach, we find FAR values between 0.3 and 0.7 for POH and FAR > 0.6 for MESHS. For POH, the highest CSI (0.37) and HSS (0.52) are obtained using a 60 % threshold, while for MESHS the highest CSI (0.25) and HSS (0.4) are obtained using a 2 cm threshold. We find that the current calibration of POH does not correspond to a probability and suggest a recalibration based on the filtered reports. [ABSTRACT FROM AUTHOR]

Observations of glacier mass changes are key to understanding the response of glaciers to climate change and related impacts, such as regional runoff, ecosystem changes, and global sea level rise. Spaceborne optical and radar sensors make it possible to quantify glacier elevation changes, and thus multi-annual mass changes, on a regional and global scale. However, estimates from a growing number of studies show a wide range of results with differences often beyond uncertainty bounds. Here, we present the outcome of a community-based inter-comparison experiment using spaceborne optical stereo (ASTER) and synthetic aperture radar interferometry (TanDEM-X) data to estimate elevation changes for defined glaciers and target periods that pose different assessment challenges. Using provided or self-processed digital elevation models (DEMs) for five test sites, 12 research groups provided a total of 97 spaceborne elevation-change datasets using various processing approaches. Validation with airborne data showed that using an ensemble estimate is promising to reduce random errors from different instruments and processing methods but still requires a more comprehensive investigation and correction of systematic errors. We found that scene selection, DEM processing, and co-registration have the biggest impact on the results. Other processing steps, such as treating spatial data voids, differences in survey periods, or radar penetration, can still be important for individual cases. Future research should focus on testing different implementations of individual processing steps (e.g. co-registration) and addressing issues related to temporal corrections, radar penetration, glacier area changes, and density conversion. Finally, there is a clear need for our community to develop best practices, use open, reproducible software, and assess overall uncertainty to enhance inter-comparison and empower physical process insights across glacier elevation-change studies. [ABSTRACT FROM AUTHOR]

The firn layer on the Greenland Ice Sheet (GrIS) plays a crucial role in buffering surface meltwater runoff, which is constrained by the available firn pore space and impermeable ice layers that limit deeper meltwater percolation. Understanding these firn properties is essential for predicting current and future meltwater runoff and its contribution to global sea-level rise. While very-high-frequency (VHF) radars have been extensively used for surveying the GrIS, their lower bandwidth restricts direct firn stratigraphy extraction. In this study, we use concurrent VHF airborne radar and laser altimetry data collected as part of Operation IceBridge over the 2011–2019 period to investigate our hypothesis that vertical heterogeneities in firn (i.e. ice layers) cause vertical offsets in the radar surface reflection (dz). Our results, corroborated by modelling and firn core analyses, show that a dz larger than 1 m is strongly related to the vertical heterogeneity of a firn profile and effectively delineates between vertically homogeneous and vertically heterogeneous firn profiles over a depth range of ∼ 4 m. Temporal variations in dz align with climatic events and reveal an expansion of heterogeneous firn between 2011–2013 covering an area of ∼ 350 815 km 2 , followed by firn replenishment over the years 2014–2019 spanning an area of ∼ 667 725 km 2. Our approach reveals the firn evolution of key regions on the Greenland Ice Sheet, providing valuable insights for detecting potential alterations in meltwater runoff patterns. [ABSTRACT FROM AUTHOR]

We discuss robust estimations for the variance of normally distributed random variables in the presence of interference. The robust estimators are based on either ranking or the geometric mean. For the interference models used, estimators based on the geometric mean outperform the rank-based ones in both mitigating the effect of interference and reducing the statistical error when there is no interference. One reason for this is that estimators using the geometric mean do not suffer from the "heavy tail" phenomenon like the rank-based estimators do. The ratio of the standard deviation over the mean of the power random variable is sensitive to interference. It can thus be used as a criterion to combine the sample mean with a robust estimator to form a hybrid estimator. We apply the estimators to the Arecibo incoherent scatter radar signals to determine the total power and Doppler velocities in the ionospheric E-region altitudes. Although all the robust estimators selected deal with light contamination well, the hybrid estimator is most effective in all circumstances. It performs well in suppressing heavy contamination and is as efficient as the sample mean in reducing the statistical error. Accurate incoherent scatter radar measurements, especially at nighttime and at E-region altitudes, can improve studies of ionospheric dynamics and composition. [ABSTRACT FROM AUTHOR]

To better understand the interactions between C-band radar waves and snow, a tower-based experiment was set up in the Idaho Rocky Mountains for the period of 2021–2023. The experiment objective was to improve understanding of the sensitivity of Sentinel-1 C-band backscatter radar signals to snow. The data were collected in the time domain to measure the backscatter profile from the various snowpack and ground surface layers. The data show that scattering is present throughout the snow volume, although it is limited for low snow densities. Contrasting layer interfaces, ice features and metamorphic snow can have considerable impact on the backscatter signal. During snowmelt periods, wet snow absorbs the signal, and the soil backscatter becomes negligible. A comparison of the vertically integrated tower radar data with Sentinel-1 data shows that both systems have similar temporal behavior, and both feature an increase in backscatter during the dry-snow period in 2021–2022, even during weeks of nearly constant snow depth, likely due to morphological changes in the snowpack. The results demonstrate that C-band radar is sensitive to the dominant seasonal patterns in snow accumulation but that changes in microstructure, stratigraphy, melt–freeze cycles and snow wetness may complicate satellite-based snow depth retrievals. [ABSTRACT FROM AUTHOR]

The Eastern Pacific Cloud Aerosol Precipitation Experiment (EPCAPE) was a year-round campaign conducted by the US Department of Energy at the Scripps Institution of Oceanography in La Jolla, CA, USA, with a focus on characterizing atmospheric processes at a coastal location. The ground-based prototype of a new Ka-, W-, and G-band (35.75, 94.88, and 238.8 GHz) profiling atmospheric radar, named CloudCube, which was developed at the Jet Propulsion Laboratory, took part in the experiment during 6 weeks in March and April 2023. This article describes the unique data sets that were obtained during the field campaign from a variety of marine clouds and light precipitation. These are, to the best of the authors' knowledge, the first observations of atmospheric clouds using simultaneous multifrequency measurements including 238.8 GHz. These data sets therefore provide an exceptional opportunity to study and analyze hydrometeors with diameters in the millimeter- and submillimeter size range that can be used to better understand cloud and precipitation structure, formation, and evolution. The data sets referenced in this article are intended to provide a complete, extensive, and high-quality collection of G-band data in the form of Doppler spectra and Doppler moments. In addition, Ka-band and W-band reflectivity and Ka-, W-, and G-band reflectivity ratio profiles are included for several cases of interest on 6 different days. The data sets can be found at https://doi.org/10.5281/zenodo.10076227 (Socuellamos et al., 2024). [ABSTRACT FROM AUTHOR]

Global marine gravity anomaly models are predominantly recovered from along-track radar altimeter data. While remarkable advancements has been achieved in gravity anomaly modelling, the quality of gravity anomaly model remains constrained by the absence of across-track geoid gradients and the reduction of radar altimeter data, particularly in coastal and high-latitudes regions. ICESat-2 laser altimetry operates three-pair laser beams with a small footprint and near-polar orbit, enabling the determination of across-track geoid gradients and providing more valid observations in certain regions. The ICESat-2 altimeter data processing method is presented including the determination of across-track geoid gradients and the combination of along/across-track geoid gradients. A new global marine gravity model, SDUST2022GRA, is recovered from radar and laser altimeter data using different method for determining each altimeter data error. The accuracy and spatial resolution of SDUST2022GRA is assessed by published global gravity anomaly models (DTU17, V32.1, NSOAS22) and available shipborne gravity measurements. The accuracy of SDUST2022GRA is 4.43 mGal on a global scale, which is at least 0.22 mGal better than that of others models. Moreover, in local coastal and high-latitude regions, SDUST2022GRA achieves an accuracy improvement of 0.16–0.24 mGal compared to others models. The spatial resolution of SDUST2022GRA is approximately 20 km in a certain region, slightly better superior others models. These assessments suggests that SDUST2022GRA is a reliable global marine gravity anomaly model. By comparing SDUST2022GRA with incorporating ICESat-2 and SDUST2021GRA without ICESat-2, the percentage contribution of ICESat-2 to the improvement of gravity anomaly model accuracy is 13 % in the global ocean region, and it is increasing with an proportion of ICESat-2 altimeter data in high-latitude and coastal regions. The SDUST2022GRA are freely available at the site of https://doi.org/10.5281/zenodo.8337387 (Li et al., 2023). [ABSTRACT FROM AUTHOR]

This paper presents an open-access data set of reanalysed radar reflectivities and rainfall rates at sub-kilometre spatial and minute temporal scales. Variability at these scales is a blind spot for both operational rain gauge networks and operational radar networks. In the urban area of Hamburg, precipitation measurements of a single-polarized X-band weather radar operating at high temporal (30s), range (60m), and azimuthal sampling (1°) resolutions are made available for a period of more than 8 years. We describe in detail the reanalysis of the raw radar data, outline the radar performance for the years 2013 to 2021, and discuss open issues and limitations of the data set. Several sources of radar-based errors were adjusted gradually, affecting the radar reflectivity and rainfall measurements, e.g. noise, alignment, non-meteorological echoes, radar calibration, and attenuation. The deployment of additional vertically pointing micro rain radars yields drop size distributions at the radar beam height, which effectively reduces errors concerning the radar calibration and attenuation correction and monitors the radar data quality. A statistical evaluation revealed that X-band radar reflectivities and rainfall rates are in very good agreement with the micro rain radar measurements. Moreover, the analyses of rainfall patterns shown for an event and accumulated rainfall of several months prove the quality of the data set. The provided radar reflectivities facilitate studies on attenuation correction and the derivation of further weather radar products, like an improved rainfall rate. The rainfall rates themselves can be used for studies on the spatial and temporal scales of precipitation and hydrological research, e.g. input data for high-resolution modelling, in an urban area. The radar reflectivities and rainfall rates are available at https://doi.org/10.26050/WDCC/LAWR_UHH_HHG_v2. [ABSTRACT FROM AUTHOR]

We present the geophysical data set acquired in summer 2022 close to Ny-Ålesund (western Svalbard, Brøggerhalvøya Peninsula, Norway) as part of the project ICEtoFLUX. The aim of the investigation is to characterize the role of groundwater flow through the active layer as well as through and/or below the permafrost. The data set is composed of electrical resistivity tomography (ERT) and ground-penetrating radar (GPR) surveys, which are well-known geophysical techniques for the characterization of glacial and hydrological processes and features. Overall, 18 ERT profiles and 10 GPR lines were acquired, for a total surveyed length of 9.3 km. The data have been organized in a consistent repository that includes both raw and processed (filtered) data. Some representative examples of 2D models of the subsurface are provided, that is, 2D sections of electrical resistivity (from ERT) and 2D radargrams (from GPR). The resistivity models revealed deep resistive structures, probably related to the heterogeneous permafrost, which are often interrupted by electrically conductive regions that may relate to aquifers and/or faults. The interpretation of these data can support the identification of the active layer, the occurrence of spatial variation in soil conditions at depth, and the presence of groundwater flow through the permafrost. To a large extent, the data set can provide new insight into the hydrological dynamics and polar and climate change studies of the Ny-Ålesund area. The data set is of major relevance because there are few geophysical data published about the Ny-Ålesund area. Moreover, these geophysical data can foster multidisciplinary scientific collaborations in the fields of hydrology, glaciology, climate, geology, and geomorphology, etc. The geophysical data are provided in a free repository and can be accessed at 10.5281/zenodo.10260056 (Pace et al., 2023). [ABSTRACT FROM AUTHOR]

Accurate estimations of extreme precipitation return levels are critical for many hydrological applications. Extreme precipitation is highly variable in both space and time; therefore, to better understand and manage the related risks, knowledge of their probability at different spatial–temporal scales is crucial. We employ a novel non-asymptotic framework to estimate extreme return levels (up to 100 years) at multiple spatial–temporal scales from weather radar precipitation estimates. The approach reduces uncertainties and enables the use of relatively short archives typical of weather radar data (12 years in this case). We focus on the eastern Mediterranean, an area of high interest due to its sharp climatic gradient, containing Mediterranean, semi-arid, and arid areas across a few tens of kilometres, and its susceptibility to flash flood. At-site intensity–duration–area–frequency relations are derived from radar precipitation data at various scales (10 min–24 h, 0.25–500 km 2) across the study area, using ellipses of varying axes and orientations to account for the spatial component of storms. We evaluate our analysis using daily rain gauge data over areas for which sufficiently dense gauge networks are available. We show that extreme return levels derived from radar precipitation data for 24 h and 100 km 2 are generally comparable to those derived from averaging daily rain gauge data over a similar areal scale. We then analyse differences in multi-scale extreme precipitation over coastal, mountainous, and desert regions. Our study reveals that the power-law scaling relationship between precipitation and duration (simple scaling) weakens for increasing area sizes. This finding has implications for temporal downscaling. Additionally, precipitation intensity varies significantly for different area sizes at short durations but becomes more similar at long durations, suggesting that, in the region, areal reduction factors may not be necessary for computing return levels over long durations. Furthermore, the reverse orographic effect, which causes decreased precipitation for hourly and sub-hourly durations, diminishes for larger areas. Finally, we discuss the effects of orography and coastline proximity on extreme precipitation intensity over different spatial–temporal scales. [ABSTRACT FROM AUTHOR]

Stratiform rain and the overlying ice play crucial roles in Earth's climate system. From a microphysics standpoint, water mass flux primarily depends on two variables: particles' concentration and their mass. The Dual-frequency Precipitation Radar (DPR) on the Global Precipitation Measurement mission core satellite is a spaceborne instrument capable of estimating these two quantities through dual-wavelength measurements. In this study, we evaluate bulk statistics on the ice particle properties derived from dual-wavelength radar data in relation to the properties of rain underneath. Specifically, we focus on DPR observations over stratiform precipitation, characterized by columns exhibiting a prominent bright band, where the melting layer can be easily detected. Our analysis reveals a large increase in the retrieved mass flux as we transition from the ice to the rain phase in the official DPR product. This observation is in disagreement with our expectation that mass flux should remain relatively stable across the bright band in cold-rain conditions. To address these discrepancies, we propose an alternative retrieval algorithm that ensures a gradual transition of Dm (mean mass-weighted particle melted-equivalent diameter) and the precipitation rate across the melting zone. This approach also helps in estimating bulk ice density above the melting level. These findings demonstrate that DPR observations can not only quantify ice particle content and their size above stratiform rain regions but also estimate bulk density, provided uniform conditions that minimize uncertainties related to partial beam filling. [ABSTRACT FROM AUTHOR]

Rain cells are the most common units in the natural precipitation system. Enhancing the understanding of these rain cell characteristics can significantly improve the cognition of the precipitation system. Previous studies have mostly analyzed rain cells from a single radar data. In this study, we merged the precipitation parameters measured by the Tropical Rainfall Measuring Mission (TRMM) precipitation radar (PR) with the muti-channel cloud-top radiance measured by the visible and infrared scanner (VIRS) and the muti-channel brightness temperature measured by the TRMM microwave imager (TMI). The rain cells were identified within the PR orbit, and the swath truncation effect was eliminated. We used two methods for rain cell identification: the minimum bounding rectangle (MBR) method and the best fit ellipse (BFE) method, and compared the differences between these two methods in describing the rain cell characteristics. The results indicate that both methods can better reflect the geometric characteristics of rain cells. Compared with the MBR method, the BFE method can obtain a smaller rain cell area, and the filling ratio is better. However, the MBR method can simplify the data storage volume. Consequently, we employed the MBR method to analyze the precipitation structure of two typical rain cell precipitation cases. The results show that the new rain cell dataset can be used for the analysis of rain cell precipitation parameters and visible/infrared and microwave signals, which provides valuable data for comprehensive studies on the rain cell structural characteristics and furthers the understanding of precipitation mechanisms. The data which were used in this paper are freely available at https://doi.org/10.5281/zenodo.8352587 (Wu et al., 2023). [ABSTRACT FROM AUTHOR]

A physical–statistical framework to estimate snow water equivalent (SWE) and snow depth from synthetic aperture radar (SAR) measurements is presented and applied to four SnowSAR flight-line data sets collected during the SnowEx'2017 field campaign in Grand Mesa, Colorado, USA. The physical (radar) model is used to describe the relationship between snowpack conditions and volume backscatter. The statistical model is a Bayesian inference model that seeks to estimate the joint probability distribution of volume backscatter measurements, snow density and snow depth, and physical model parameters. Prior distributions are derived from multilayer snow hydrology predictions driven by downscaled numerical weather prediction (NWP) forecasts. To reduce the signal-to-noise ratio, SnowSAR measurements at 1 m resolution were upscaled by simple averaging to 30 and 90 m resolution. To reduce the number of physical parameters, the multilayer snowpack is transformed for Bayesian inference into an equivalent one- or two-layer snowpack with the same snow mass and volume backscatter. Successful retrievals meeting NASEM (2018) science requirements are defined by absolute convergence backscatter errors ≤1.2 dB and local SnowSAR incidence angles between 30 and 45 ∘ for X- and Ku-band VV-pol backscatter measurements and were achieved for 75 % to 87 % of all grassland pixels with SWE up to 0.7 m and snow depth up to 2 m. SWE retrievals compare well with snow pit observations, showing strong skill in deep snow with average absolute SWE residuals of 5 %–7 % (15 %–18 %) for the two-layer (one-layer) retrieval algorithm. Furthermore, the spatial distributions of snow depth retrievals vis-à-vis lidar estimates have Bhattacharya coefficients above 94 % (90 %) for homogeneous grassland pixels at 30 m (90 m resolution), and values up to 76 % in mixed forest and grassland areas, indicating that the retrievals closely capture snowpack spatial variability. Because NWP forecasts are available everywhere, the proposed approach could be applied to SWE and snow depth retrievals from a dedicated global snow mission. [ABSTRACT FROM AUTHOR]

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. 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 the evaluation of kilometer-scale numerical weather prediction models with EarthCARE CPR observations. [ABSTRACT FROM AUTHOR]

A new dual-frequency stratospheric–tropospheric (ST) and meteor radar has been built and installed at the Langfang Observatory in northern China. It utilizes a new two-frequency system design that allows interleaved operation at 53.8 MHz for ST mode and at 35.0 MHz for meteor mode, thus optimizing performance for both ST wind retrieval and meteor trail detection. In dedicated meteor mode, the daily meteor count rate reaches over 40 000 and allows wind estimation at finer time resolutions than the 1 h typical of most meteor radars. The root mean square uncertainty of the ST wind measurements is better than 2 m s -1 when estimating the line of best fit with radiosonde winds. Preliminary observation results for 1 month of winter gravity wave (GW) momentum fluxes in the mesosphere, lower stratosphere and troposphere are also presented. A case of waves generated by the passage of a cold front is found. [ABSTRACT FROM AUTHOR]

Mixed-phase clouds are not well represented in climate and weather forecasting models, due to a lack of the key processes controlling their life cycle. Developing methods to study these clouds is therefore essential, despite the complexity of mixed-phase cloud processes and the difficulty of observing two cloud phases simultaneously. We propose in this paper a new method to retrieve the microphysical properties of mixed-phase clouds, ice clouds and supercooled water clouds using airborne or satellite radar and lidar measurements, called VarPy-mix. This new approach extends an existing variational method developed for ice cloud retrieval using lidar, radar and passive radiometers. We assume that the lidar attenuated backscatter β at 532 nm is more sensitive to particle concentration and is consequently mainly sensitive to the presence of supercooled water. In addition, radar reflectivity Z at 95 GHz is sensitive to the size of hydrometeors and hence more sensitive to the presence of ice particles. Consequently, in the mixed phase the supercooled droplets are retrieved with the lidar signal and the ice particles with the radar signal, meaning that the retrievals rely strongly on a priori and error values. This method retrieves simultaneously the visible extinction for ice αice and liquid αliq particles, the ice and liquid water contents IWC and LWC, the effective radius of ice re,ice and liquid re,liq particles, and the ice and liquid number concentrations Nice and Nliq. Moreover, total extinction αtot , total water content (TWC) and total number concentration Ntot can also be estimated. As the retrieval of ice and liquid is different, it is necessary to correctly identify each phase of the cloud. To this end, a cloud-phase classification is used as input to the algorithm and has been adapted for mixed-phase retrieval. The data used in this study are from DARDAR-MASK v2.23 products, based on the Cloud-Aerosol Lidar with Orthogonal Polarization (CALIOP) and Cloud Profiling Radar (CPR) observations from the CALIPSO and CloudSat satellites, respectively, belonging to the A-Train constellation launched in 2006. Airborne in situ measurements performed on 7 April 2007 during the Arctic Study of Tropospheric Aerosol, Clouds and Radiation (ASTAR) campaign and collected under the track of CloudSat–CALIPSO are compared with the retrievals of the new algorithm to validate its performance. Visible extinctions, water contents, effective radii and number concentrations derived from in situ measurements and the retrievals showed similar trends and are globally in good agreement. The mean percent error between the retrievals and in situ measurements is 39 % for αliq , 398 % for αice , 49 % for LWC and 75 % for IWC. It is also important to note that temporal and spatial collocations are not perfect, with a maximum spatial shift of 1.68 km and a maximum temporal shift of about 10 min between the two platforms. In addition, the sensitivity of remote sensing and that of in situ measurements is not the same, and in situ measurement uncertainties are between 25 % and 60 %. [ABSTRACT FROM AUTHOR]

Global measurements of horizontal winds in cloud and precipitation systems represent a gap in the global observation system. The Wind Velocity Radar Nephoscope (WIVERN) mission, one of the two candidates to be the ESA's Earth Explorer 11 mission, aims at filling this gap based on a conically scanning W-band Doppler radar instrument. The determination of the antenna boresight mispointing angles and the impact of their uncertainty on the line of sight Doppler velocities is critical to achieve the mission requirements. While substantial industrial efforts are on their way to achieving accurate determination of the pointing, alternative (external) calibration approaches are currently under scrutiny. The correction of the line of sight Doppler velocity error introduced by the mispointing only needs knowledge of such mispointing angles and does not need the correction of the mispointing itself. Thus, this work discusses four methods applicable to the WIVERN radar that can be used at different timescales to characterize the antenna mispointing both in the azimuthal and in the elevation directions and to correct the error in the Doppler velocity induced by such mispointing. Results show that elevation mispointing is well corrected at very short timescales by monitoring the range at which the surface peak occurs. Azimuthal mispointing is harder but can be tackled by using the expected profiles of the non-moving surface Doppler velocity. Biases in pointing at longer timescales can be monitored by using a well-established reference database (e.g. ECMWF reanalysis) or ad-hoc ground-based calibrators. Although tailored to the WIVERN mission, the proposed methodologies can be extended to other Doppler concepts featuring conically scanning or slant viewing Doppler systems. [ABSTRACT FROM AUTHOR]

The Goddard Profiling Algorithm (GPROF) is used operationally for the retrieval of surface precipitation and hydrometeor profiles from the passive microwave (PMW) observations of the Global Precipitation Measurement (GPM) mission. Recent updates have led to GPROF V7, which has entered operational use in May 2022. In parallel, development is underway to improve the retrieval by transitioning to a neural-network-based algorithm called GPROF-NN. This study validates retrievals of liquid precipitation over snow-free and non-mountainous surfaces from GPROF V7 and multiple configurations of GPROF-NN against ground-based radar measurements over the conterminous United States (CONUS) and the tropical Pacific. GPROF retrievals from the GPM Microwave Imager (GMI) are validated over several years, and their ability to reproduce regional precipitation characteristics and effective resolution is assessed. Moreover, the retrieval accuracy for several other sensors of the constellation is evaluated. The validation of GPROF V7 indicates that the retrieval produces reliable estimates of liquid precipitation over the CONUS. During all four assessed years, annual mean precipitation is within 8 % of gauge-corrected radar measurements. Although biases of up to 25 % are observed over sub-regions of the CONUS and the tropical Pacific, the retrieval reliably reproduces each region's diurnal and seasonal precipitation characteristics. The effective resolution of GPROF V7 is found to be 51 km over the CONUS and 18 km over the tropical Pacific. GPROF V7 also produces robust precipitation estimates for the other sensors of the GPM constellation. The evaluation further shows that the GPROF-NN retrievals have the potential to significantly improve the GPM PMW precipitation retrievals. GPROF-NN 1D, the most basic neural network implementation of GPROF, improves the mean-squared error, mean absolute error, correlation and symmetric mean absolute percentage error of instantaneous precipitation estimates by about 20 % for GPROF GMI while the effective resolution is improved to 31 km over land and 15 km over oceans. The two GPROF-NN retrievals that are based on convolutional neural networks can further improve the accuracy up to the level of the combined radar–radiometer retrievals from the GPM core observatory. However, these retrievals are found to overfit on the viewing geometry at the center of the swath, reducing their overall accuracy to that of GPROF-NN 1D. For the other sensors of the constellation, the GPROF-NN retrievals produce larger biases than GPROF V7 and only GPROF-NN 3D achieves consistent improvements compared to GPROF V7 in terms of the other assessed error metrics. This points to shortcomings in the hydrometeor profiles or radiative transfer simulations used to generate the training data for the other sensors of the GPM constellation as a critical limitation for improving GPM PMW retrievals. [ABSTRACT FROM AUTHOR]

Measurements of radar backscatter are sensitive to snow water equivalent (SWE) across a wide range of frequencies, motivating proposals for satellite missions to measure global distributions of SWE. However, radar backscatter measurements are also sensitive to snow stratigraphy, to microstructure, and to ground surface roughness, complicating SWE retrieval. A number of recent advances have created new tools and datasets with which to address the retrieval problem, including a parameterized relationship between SWE, microstructure, and radar backscatter, and methods to characterize ground surface scattering. Although many algorithms also introduce external (prior) information on SWE or snow microstructure, the precision of the prior datasets used must be high in some cases in order to achieve accurate SWE retrieval. We hypothesize that a time series of radar measurements can be used to solve this problem and demonstrate that SWE retrieval with acceptable error characteristics is achievable by using previous retrievals as priors for subsequent retrievals. We demonstrate the accuracy of three configurations of prior information: using a global SWE model, using the previously retrieved SWE, and using a weighted average of the model and the previous retrieval. We assess the robustness of the approach by quantifying the sensitivity of the SWE retrieval accuracy to SWE biases artificially introduced in the prior. We find that the retrieval with the weighted averaged prior demonstrates SWE accuracy better than 20 % and an error increase of only 3 % relative RMSE per 10 % change in prior bias; the algorithm is thus both accurate and robust. This finding strengthens the case for future radar-based satellite missions to map SWE globally. [ABSTRACT FROM AUTHOR]

Abstract. In this study, we present a new forward polarimetric radar operator called the Accurate and Efficient Radar Operator designed by ZheJiang University (ZJU-AERO). This operator was designed to interface with the numerical weather prediction NWP) model of the global forecast system/regional mesoscale model of the China Meteorology Administration (CMA- GFS/MESO). The main objective of developing this observation operator was to assimilate observations from the Precipitation Measurement Radar (PMR) and ground-based radar's polarimetric radar variables, excluding the Doppler variables. To calculate the hydrometeor optical properties of ZJU-AERO, we utilized the invariant-imbedding T-matrix (IITM) method, which can handle non-spherical and inhomogeneous hydrometeor particles in the atmosphere. The optical database of ZJUAERO was designed with a multi-layered architecture to ensure the flexibility in hydrometeor morphology and orientation specifications, while maintaining operational efficiency. In this work, we elaborate on the design concepts, physical basis, and hydrometeor specifications of ZJU-AERO. Additionally, we present a case study demonstrating the application of ZJU-AERO in simulating radar observations of Typhoon Haishen. [ABSTRACT FROM AUTHOR]

Ground-based radar precipitation products typically need adjustment with rain gauge accumulations to achieve a reasonable accuracy. This is certainly the case for the pan-European radar precipitation products. The density of (near) real-time rain gauge accumulations from official networks is often relatively low. Crowdsourced rain gauge networks have a much higher density than conventional ones and are a potentially interesting (complementary) source to merge with radar precipitation accumulations. Here, a 1-year personal weather station (PWS) rain gauge dataset of ∼ 5 min accumulations is obtained from the private company Netatmo over the period 1 September 2019–31 August 2020, which is subjected to quality control using neighbouring PWSs and, after aggregating to 1 h accumulations, using unadjusted radar data. The PWS 1 h gauge accumulations are employed to spatially adjust OPERA radar accumulations, covering 78 % of geographical Europe. The performance of the merged dataset is evaluated against daily and disaggregated 1 h gauge accumulations from weather stations in the European Climate Assessment & Dataset (ECA&D). Results are contrasted to those from an unadjusted OPERA-based radar dataset and from EURADCLIM. The severe average underestimation for daily precipitation of ∼ 28 % from the unadjusted radar dataset diminishes to ∼ 3 % for the merged radar–PWS dataset. A station-based spatial verification shows that the relative bias in 1 h precipitation is still quite variable and suggests stronger underestimations for colder climates. A dedicated evaluation with scatter density plots reveals that the performance is indeed less good for lower temperatures, which points to limitations in observing solid precipitation by PWS gauges. The outcome of this study confirms the potential of crowdsourcing to improve radar precipitation products in (near) real time. [ABSTRACT FROM AUTHOR]

Some of the highest specific mass change rates in Antarctica are reported for the Antarctic Peninsula. However, the existing estimates for the northern Antarctic Peninsula (<70 ∘ S) are either spatially limited or are affected by considerable uncertainties. The complex topography, frequent cloud cover, limitations in ice thickness information, boundary effects, and uncertain glacial–isostatic adjustment estimates affect the ice sheet mass change estimates using altimetry, gravimetry, or the input-output method. Within this study, the first assessment of the geodetic mass balance throughout the ice sheet of the northern Antarctic Peninsula is carried out employing bi-static synthetic aperture radar (SAR) data from the TanDEM-X satellite mission. Repeat coverages from the austral winters of 2013 and 2017 are employed. Overall, coverage of 96.4 % of the study area by surface elevation change measurements and a total mass budget of -24.1±2.8 Gt a -1 are revealed. The spatial distribution of the surface elevation and mass changes points out that the former ice shelf tributary glaciers of the Prince Gustav Channel, Larsen A and B, and Wordie ice shelves are the hotspots of ice loss in the study area and highlights the long-lasting dynamic glacier adjustments after the ice shelf break-up events. The highest mass change rate is revealed for the Airy–Seller–Fleming glacier system at -4.9±0.6 Gt a -1 , and the highest average surface elevation change rate of -2.30±0.03 m a -1 is observed at Drygalski Glacier. The comparison of the ice mass budget with anomalies in the climatic mass balance indicates, that for wide parts of the southern section of the study area, the mass changes can be partly attributed to changes in the climatic mass balance. However, imbalanced high ice discharge drives the overall ice loss. The previously reported connection between mid-ocean warming along the southern section of the west coast and increased frontal glacier recession does not repeat in the pattern of the observed glacier mass losses, excluding in Wordie Bay. The obtained results provide information on ice surface elevation and mass changes for the entire northern Antarctic Peninsula on unprecedented spatially detailed scales and with high precision and will be beneficial for subsequent analysis and modeling. [ABSTRACT FROM AUTHOR]

Upper-atmosphere winds from a meteor radar and a field-widened Michelson interferometer, co-located at the Polar Environment Atmospheric Research Laboratory in Eureka, Nunavut, Canada (80° N, 86° W) are compared. The two instruments implement different wind-measuring techniques at similar heights and have very different temporal and spatial observational footprints. The meteor radar provides winds averaged over a ∼ 300 km horizontal area in 3 km vertical bins between 82 and 97 km on a 1 h cadence. The E-Region Wind Interferometer II (ERWIN) provides airglow-weighted winds (averaged over volumes of ∼ 8 km in height by ∼ 5 km radius) from three nightglow emissions (O(1S) , oxygen green line, 557.7 nm , 97 km ; an O2 line, 866 nm , 94 km ; and an OH line, 843 nm , 87 km) on a ∼ 5 min cadence. ERWIN's higher precision (1–2 ms-1 for the O(1S) and OH emissions and ∼ 4 ms-1 for the O2 emissions) and higher cadence allows more substantive comparisons between winds measured by meteor radar and Doppler shifts in airglow emissions than previously possible for similar meteor radar/airglow Doppler shift comparisons using Fabry–Perot interferometers. The best correlation is achieved using Gaussian weighting of meteor radar winds with peak height and vertical width being optimally determined. Peak heights agree well with co-located SABER airglow observations. Offsets between the two instruments are ∼ 1–2 ms-1 for the O2 and O(1S) emissions and less than 0.3 ms-1 for the OH emission. Wind directions are highly correlated with a ∼ 1:1 correspondence. On average, meteor radar wind magnitudes are ∼ 40 % larger than those from ERWIN. Gravity wave airglow brightness weighting of observations is discussed. Non-quadrature phase offsets between the airglow weighting and gravity wave associated wind and temperature perturbations will result in enhanced or reduced layer-weighted wind amplitudes. [ABSTRACT FROM AUTHOR]

In dry southwestern South Tyrol, Italy, rock glaciers are dominant landforms of the high-mountain cryosphere. Their spatial distribution and degree of activity hold critical information on the current state of discontinuous permafrost and consequently on the response potential to climate warming. Traditional geomorphologic mapping, however, owing to the qualitative expert-based nature, typically displays a high degree of uncertainty and variability among operators with respect to the dynamic classification of intact (permafrost-bearing) and relict (permafrost-devoid) rock glaciers. This limits the reliability of geomorphologic rock glacier inventories for basic and applied purposes. To address this limitation, (i) we conduct a systematic evaluation of the improvements that synthetic aperture radar interferometry (InSAR) can afford to the detection and dynamic classification of rock glaciers and (ii) build an integrated inventory that combines the strengths of geomorphologic- and InSAR-based approaches. To exploit fully InSAR-based information towards a better understanding of the topo-climatic conditions that sustain creeping permafrost, we further explore how velocity and the spatial distribution of moving areas (MAs) within rock glaciers may vary as a function of simple topographic variables known to exert first-order controls on incoming solar radiation, such as elevation and aspect. Starting from a geomorphologic inventory (n=789), we characterize the kinematics of InSAR-based MAs and the relevant hosting rock glaciers on 36 Sentinel-1 interferograms in the 2018–2019 period. With respect to the original inventory, InSAR analysis allowed us to identify 14 previously undetected rock glaciers. Further, it confirmed that 246 (76 %) landforms, originally interpreted as intact, do exhibit detectable movement (i.e., ≥1 cm yr -1) and that 270 (60 %) of the relict labeled counterparts do not, whereas 144 (18 %) were kinematically undefined due to decorrelation. Most importantly, InSAR proved critical for reclassifying 121 (15 %) rock glaciers, clarifying that 41 (13 %) of those interpreted as being intact do not exhibit detectable movement and that 80 (17 %) of the original relict ones do move. Reclassification (i) allowed us to identify a cluster of intact rock glaciers below 2000 m a.s.l. associated with positive mean annual air temperature (MAAT), and (ii) by increasing the altitudinal overlap between intact and relict rock glaciers, it depicts a broad transition belt in the aspect–elevation space, which varies from 50 m on west-facing slopes to 500 m on easterly ones. This finding deteriorates the significance of elevation and aspect as topographic proxies for modeling permafrost occurrence and highlights the importance of using InSAR to inform such models. From a process-oriented standpoint, InSAR information proves fundamental for imaging how this altitudinal transition manifests through changing rates and styles of rock glacier surface deformation. Specifically, we find that, as rock glaciers move faster, an increasingly larger proportion of their surface becomes kinematically involved (i.e., percent MA cover) and that this proportion increases with elevation up to 2600–2800 m, beyond which an inflection occurs and consistent average values are attained. Considering that the inflection falls between the -1 and -2 °C MAAT – the lower boundary for discontinuous permafrost – and is independent of slope gradient, we conclude that this altitudinal pattern represents a geomorphic signature: the dynamic expression of increasing permafrost distribution, from sporadic to discontinuous. [ABSTRACT FROM AUTHOR]

In recent years, the application of compact and cost-effective deployable X-band polarimetric radars has gained in popularity, particularly in regions with complex terrain. The deployable radars generally use a radome constructed by joining multiple panels using metallic threads to facilitate easy transportation. As a part of the ESPOIRS project, Laboratoire de l'Atmosphεave;re et des Cyclones has acquired an X-band meteorological radar with four panel radome configuration. In this study, we investigated the effect of the radome on the measured polarimetric variables, particularly differential reflectivity and differential phase. Our observations reveal that the metallic threads connecting the radome panels introduce power loss at vertical polarization, leading to a positive bias in the differential reflectivity values. To address the spatial variability bias observed in differential reflectivity and differential phase, we have developed a novel algorithm based on the Discrete Fourier Transform. The algorithm's performance was tested during an intense heavy rainfall event caused by the Batsirai cyclone on Reunion Island. The comparative and joint histogram analysis demonstrates the algorithm's effectiveness in correcting the spatial bias in the polarimetric variables. [ABSTRACT FROM AUTHOR]

Cloud radar Doppler spectra are of particular interest for investigating cloud microphysical processes, such as ice formation, riming and ice multiplication. When hydrometeor types within a cloud radar observation volume have sufficiently different terminal fall velocities, they produce individual Doppler spectrum peaks, convoluted by dynamical effects. If these (sub-)peaks can be separated, properties of the underlying hydrometeor populations can potentially be estimated, such as their fall velocity, number, size and to some extent their shape. However, this task is complex and dependent on the cloud radar operation settings, atmospheric dynamics and hydrometeor characteristics. As a consequence, there is a need for adjustable tools that are able to detect peaks in cloud radar Doppler spectra to extract the valuable information contained in them. This paper presents the synergistic use of two cloud radar Doppler spectra peak analysis algorithms, PEAKO and peakTree. PEAKO is a supervised machine learning tool that can be trained to obtain the optimal parameters for peak detection in cloud radar Doppler spectra for specific cloud radar instrument settings. The learned Doppler spectrum peak detection parameters can then be applied by peakTree, which is used to detect, structure and interpret Doppler spectrum peaks. The application of the improved PEAKO-peakTree toolkit is demonstrated in two case studies. The interpretation is supported by forward simulated cloud radar Doppler spectra by the Passive and Active Microwave TRAnsfer tool (PAMTRA), which are also used to explore the limitations of the algorithm toolkit posed by turbulence and the number of spectral averages chosen in the radar settings. From the PAMTRA simulations, we can conclude that a minimum number of 20–40 spectral averages is desirable for Doppler spectrum peak discrimination. Furthermore, liquid peaks can only be reliably separated for turbulence eddy dissipation rate values up to approximately 0.0002 m2 s−3. The first case study demonstrates that the methods work for different radar systems and settings by comparing the results for two cloud radar systems which were operated simultaneously at a site in Punta Arenas, Chile. Detected peaks which can be attributed to liquid droplets agree well between the two systems, as well as with an independent liquid-predicting neural network. The second case study compares PEAKO-peakTree-detected cloud radar Doppler spectra peaks to in situ observations collected by a balloon-based holographic imager during a campaign in Ny-Ålesund, Svalbard. This case showcases the Doppler spectrum peak detection algorithms' ability to identify different hydrometeor types, but also reveals the limitations of the algorithm toolkit posed by strong turbulence and a low number of spectral averages. Despite these challenges, the algorithm toolkit offers a powerful means of extracting comprehensive information from cloud radar observations. In the future, we envision PEAKO-peakTree application on the one hand for interpreting cloud microphysics in case studies. The identification of liquid cloud peaks emerges as a valuable asset e.g. in studies on cloud radiative effects, seeder-feeder processes, or for tracing vertical air motions. Furthermore, the computation of the moments for each sub-peak enables the tracking of hydrometeor populations and the observation of growth processes along fall streaks. On the other hand, PEAKO-peakTree application could be extended to statistical evaluations of longer data sets. Both algorithms are openly available on GitHub, offering accessibility for the scientific community. [ABSTRACT FROM AUTHOR]

In this study, we present a new forward polarimetric radar operator called the Accurate and Efficient Radar Operator designed by ZheJiang University (ZJU-AERO). This operator was designed to interface with the numerical weather prediction (NWP) model of the global forecast system/regional mesoscale model of the China Meteorological Administration (CMA-GFS/MESO). The main objective of developing this observation operator was to assimilate observations from the precipitation measurement radar (PMR). It is also capable of simulating the ground-based radar's polarimetric radar variables, excluding the Doppler variables such as radial velocity and spectrum width. To calculate the hydrometeor optical properties of ZJU-AERO, we utilize the invariant-imbedding T-matrix (IITM) method, which can handle non-spherical and inhomogeneous hydrometeor particles in the atmosphere. The optical database of ZJU-AERO was designed with a multi-layered architecture to ensure the flexibility in hydrometeor morphology and orientation specifications while maintaining operational efficiency. Specifically, three levels of databases are created that store the single-scattering properties for different shapes at discrete sizes for various fixed orientations, integrated single-scattering properties over shapes and orientations, and bulk-scattering properties incorporating the size average, respectively. In this work, we elaborate on the design concepts, physical basis, and hydrometeor specifications of ZJU-AERO. Additionally, we present a case study demonstrating the application of ZJU-AERO in simulating the radar observations of Typhoon Haishen. [ABSTRACT FROM AUTHOR]

The Cloud Profiling Radar (CPR) of the Earth Cloud, Aerosol and Radiation Explorer (EarthCARE) has a new capability to observe the Doppler velocity related to the vertical air motion of the terminal velocity of hydrometeors. The new observation from space will be used to evaluate and improve the model. Before the launch of EarthCARE, we need to develop a methodology for using the CPR data for model evaluations. In this study, we evaluated simulated data by a stretched version of the global non-hydrostatic model over Japan with a ground-based CPR using an instrument design similar to the EarthCARE CPR. We chose two cases with different precipitation events in September 2019 using two cloud microphysics schemes. We introduced the categorization method for evaluating microphysics using Doppler velocity. The results show that the liquid and solid phases of hydrometeors are divided in Doppler velocity, and the model's terminal velocities of rain, snow, and graupel categories can be evaluated with the observation. The results also show that the choice of microphysics scheme has a more significant impact than the dependence on precipitation cases. We discussed the application of the EarthCARE-like simulation results using a satellite simulator. [ABSTRACT FROM AUTHOR]

Radio-echo sounding is a standard technique for imaging the englacial stratigraphy of glaciers and ice sheets. In most cases, internal reflection horizons (IRHs) represent former glacier surfaces, comprise information about past accumulation and ice deformation, and enable the linking of ice core chronologies. IRHs in the lower third of the ice column are often difficult to detect or coherently trace. In the polar ice sheets, progress in IRH detection has been made by using multistatic, phase-coherent radars, enabling focused synthetic aperture radar (SAR) processing. However, these radar systems are often not suitable for deployment on mountain glaciers. We present a proof-of-concept study for a lightweight, phase-coherent and ground-based radar system, based on the phase-sensitive radio-echo sounder (pRES). To improve the detectability of IRHs we additionally adapted a layer-optimized SAR processing scheme to this setup. We showcase the system capability at Colle Gnifetti, Switzerland and Italy, where specular reflections are now apparent down to the base of the glacier. Compared to previously deployed impulse radar systems, with the mobile pRES the age of the oldest continuously traceable IRH could be increased from 78±12 to 288±35 a. Corresponding reflection mechanisms for this glacier are linked to stratified acidic impurities which in the upper part were deposited at a higher rate due to increased industrial activity in the area. Possible improvements to the system are discussed. If successfully implemented, these may provide a new way to map the deep internal structure of Colle Gnifetti and other mountain glaciers more extensively in future deployments. [ABSTRACT FROM AUTHOR]

Weather radar is a crucial tool in rainfall estimation, providing high-resolution estimates in both space and time. Despite this, radar rainfall estimates are subject to many error sources – including attenuation, ground clutter, beam blockage and the drop-size distribution – with the true rainfall field unknown. A flexible stochastic model for simulating errors relating to the radar rainfall estimation process is implemented, inverting standard weather radar processing methods, imposing path-integrated attenuation effects, a stochastic drop-size distribution field, along with sampling and random errors. This can provide realistic weather radar images, of which we know the true rainfall field, and the corrected 'best guess' rainfall field which would be obtained if they were observed in the real-world case. The structure of these errors is then investigated, with a focus on frequency and behaviour of 'rainfall shadows'. Half of simulated weather radar images have at least 3 % of significant rainfall rates shadowed and 25 % had at least 45 km2 containing rainfall shadows, resulting in underestimation of potential impacts of flooding. A model framework for investigating the behaviour of errors relating to the radar rainfall estimation process is demonstrated, with the flexible and efficient tool performing well at generating realistic weather radar images visually, for a large range of event types. [ABSTRACT FROM AUTHOR]

Radar observations of winter storms often exhibit locally enhanced linear features in reflectivity, sometimes labeled as snow bands. We have developed a new, objective method for detecting locally enhanced echo features in radar data from winter storms. In comparison to convective cells in warm season precipitation, these features are usually less distinct from the background echo and often have more fuzzy or feathered edges. This technique identifies both prominent, strong features and more subtle, faint features. A key difference from previous radar reflectivity feature detection algorithms is the combined use of two adaptive differential thresholds, one that decreases with increasing background values and one that increases with increasing background values. The algorithm detects features within a snow rate field rather than reflectivity and incorporates an underestimate and overestimate of feature areas to account for uncertainties in the detection. We demonstrate the technique on several examples from the US National Weather Service operational radar network. The feature detection algorithm is highly customizable and can be tuned for a variety of data sets and applications. [ABSTRACT FROM AUTHOR]

The Goddard Profiling Algorithm (GPROF) is used operationally for the retrieval of surface precipitation and hydrometeor profiles from the passive microwave (PMW) observations of the Global Precipitation Measurement (GPM) mission. Recent updates have led to GPROF V7, which has entered operational use in May 2022. In parallel, development is underway to improve the retrieval by transitioning to a neural-network-based algorithm called GPROF-NN. This study validates GPROF V7 and multiple configurations of the GPROF-NN retrieval against ground-based radar mea- surements over the conterminous United States (CONUS) and the tropical Pacific. GPROF retrievals from the GPMMicrowave Imager (GMI) are validated over several years and their ability to reproduce regional precipitation characteristics and effective resolution is assessed. Moreover, the retrieval accuracy for several other sensors of the constellation is evaluated. The validation of GPROF V7 indicates that the retrieval produces reliable precipitation estimates over CONUS. During all four assessed years, annual mean precipitation is within 8% of gauge-corrected radar measurements. Although biases of up to 25% are observed over sub-regions of the CONUS and the tropical Pacific, the retrieval reproduces the principal precipitation characteristics of each region. The effective resolution of GPROF V7 is found to be 51km over CONUS and 18km over the tropical Pacific. GPROF V7 produces robust precipitation estimates also for the other sensors of the GPM constellation. The evaluation further shows that the GPROF-NN retrievals have the potential to significantly improve the GPROF precipi- tation retrievals. GPROF-NN 1D, the most basic neural network implementation of GPROF, improves the mean-squared error, mean absolute error, correlation and symmetric mean absolute percentage error by about twenty percent for GPROF GMI while the effective resolution is improved to 31km over land and 15km over oceans. The two GPROF-NN retrievals that are based on convolutional neural networks can further improve the accuracy up to the level of the combined radar/radiometer retrievals from the GPM core observatory. However, these retrievals are found to overfit on the viewing geometry at the center of the swath, reducing their overall accuracy to that of GPROF-NN 1D. For the other sensors of the constellation, the GPROF-NN retrievals produce larger biases than GPROF V7 and only GPROF-NN 3D achieves consistent improvements compared to GPROF V7 in terms of the other assessed error metrics. This points to shortcomings in the hydrometeor profiles or radiative transfer simulations used in the training of the retrievals for the other sensors of the GPM constellation as a critical limitation for improving GPM PMW retrievals. [ABSTRACT FROM AUTHOR]

The potential ability of VHF or UHF Doppler radars to measure turbulence kinetic energy (TKE) dissipation rate ε in the atmosphere is a major asset of these instruments because of the possibility to continuously monitor turbulence in the atmospheric column above the radars. Several models have been proposed over the past decades to relate ε to half the Doppler spectral width σ , corrected for non-turbulent contributions, but their relevance remains unclear. Recently, Luce et al. (2023) tested the performance of a new model expected to be valid for weakly stratified or strongly sheared conditions, i.e., for low Richardson (Ri) numbers. Its simplest expression is εS=CSσ2S , where CS∼0.64 and S=|dV/dz| is the vertical shear of the horizontal wind V. We assessed the relevance of this model with a UHF (1.357 GHz) wind profiler called WPR LQ-7, which is routinely operated at the Shigaraki Middle and Upper Atmosphere (MU) observatory (34.85 ∘ N, 136.10 ∘ E) in Japan. For this purpose, we selected turbulence events associated with Kelvin–Helmholtz (KH) billows, because their formation necessarily requires Ri<0.25 somewhere in the flow, a condition a priori favorable to the application of the model. Eleven years of WPR LQ-7 data were used for this objective. The assessment of εS was first based on its consistency with an empirical model εLout=σ3/Lout , where Lout has the dimension of an outer scale of turbulence. It was found to compare well in a KH layer with direct estimates of ε from in situ measurements for Lout≈70 m. Some degree of equivalence between εS and εLout was confirmed by statistical analysis of 192 KH layers found in the height range 0.3–5.0 km indicating that Lout≈LH/0.64 , where LH=σ/S is the Hunt scale defined for neutral turbulence. The degree of equivalence is even significantly improved if Lout is not treated as a constant but depends on the depth D of the layer. We found Lout≈0.0875D or equivalently LH∼0.056D , which also means that σ is proportional to the apparent variation in the horizontal velocity (S×D) over the depth of the KH layer. Consequently, εS=0.64σ2S and εLout=σ3/0.0875D would express the same model for KH layers when Ri remains low. For such a condition, we provide a physical interpretation of εLout , which would be qualitatively identical to that for neutral boundary layers. [ABSTRACT FROM AUTHOR]

Accurate and precise KDP estimates are essential for radar-based applications, especially in quantitative precipitation estimation and radar data quality control routines. The accuracy of these estimates largely depends on the post-processing of the radar's measured Φ DP , which aims to reduce noise and backscattering effects while preserving fine-scale precipitation features. In this study, we evaluate the performance of several publicly available KDP estimation methods implemented in open-source libraries such as PyArt and Wradlib, and the method used in the Vaisala weather radars. To benchmark these methods, we employ a polarimetric self-consistency approach that relates KDP to reflectivity and differential reflectivity in rain, providing a reference self-consistency KDP (K DPSC) for comparison. This approach allows for the construction of the reference KDP observations that can be used to assess the accuracy and robustness of the studied KDP estimation methods. We assess each method by quantifying uncertainties using C-band weather radar observations where the reflectivity values ranged between 20 and 50 dBZ. Using the proposed evaluation framework we could define optimized parameter settings for the methods that have user-configurable parameters. Most of such methods showed significant reduction in the estimation errors after the optimization with respect to the default settings. We have found significant differences in the performances of the studied methods, where the best performing methods showed smaller normalized biases in the high reflectivity values (i.e., ≥ 40 dBZ) and overall smaller normalized root mean squared errors across the range of reflectivity values. [ABSTRACT FROM AUTHOR]

Radar based microphysical retrievals of cloud and droplet properties are vital for informing model parameterisations of clouds and precipitation but these retrievals often do not capture the details of small droplets in light rain or drizzle. A state-of-the-art G-band radar is used here to demonstrate improvements to microphysical retrievals in a case study featuring light rain. Improvements are seen, as compared to W-band radar, in the retrieval of vertical wind speed, due to the location of Mie minima at smaller droplet sizes with the G-band radar. This, in turn, has an impact on the retrieval of the drop size distribution, allowing for better accuracy in the retrieval of the characteristic drop diameter and improvements in the retrieval of the number of concentration of small droplet sizes. The differential Doppler velocity between Ka- and G-bands shows increased dynamic range compared to the Ka-W pairing, particularly for instances presenting small characteristic drop diameters. The increased attenuation experienced at G-band enables improved retrievals of liquid water content and precipitation rate when paired with W-band or Ka-band as compared to the W-band and Ka- band pairing. This is particularly noticeable in periods of light rain where the W-band and Ka-band radars receive negligible attenuation while the attenuation at G-band is much greater. [ABSTRACT FROM AUTHOR]

Operational rain gauge networks provide reliable local precipitation measurements, but are unable to represent rainfall variability for large domains at small temporal scales. In the urban area of Hamburg, precipitation measurements of a single-polarized X-band weather radar operating at a high temporal (30s), range (60m), and azimuthal sampling (1°) resolution are available for more than eight years. These measurements refine observations of the German nationwide C-band radars within the 20km scan radius. Studies on short time periods (several months) prove the performance of this low-cost local area weather radar. A case study on a tornado in a rain event demonstrates its refined resolution compared to the German nationwide C-band radars. The deployment of additional vertically pointing micro rain radars yields drop size distributions at relevant heights, which reduces errors effectively concerning the radar calibration and required statistical relations (k-Z and Z-R relation), and monitors the radar data quality. This paper describes an open-access data set covering reanalysed radar reflectivities and rainfall estimates measured by an X-band radar at high spatio-temporal resolution in the urban environment of Hamburg between 2013 and 2021. Additionally, this paper outlines the performance of this X-band radar for long time periods, and discuss open issues and limitations of the data set. The provided radar reflectivities facilitate studies on attenuation correction and the derivation of further weather radar products, like an improved rainfall rate. The rainfall rates will be used for studies on the spatial and temporal scale of precipitation and hydrological research, e.g. input data for high-resolution modelling, in an urban area. [ABSTRACT FROM AUTHOR]

WindRAD (Wind Radar) is a dual-frequency rotating fan-beam scatterometer instrument on the FY-3E (FengYun-3E) satellite. Scatterometers are generally calibrated using the linear NOC (NWP Ocean Calibration) method to control the main gain factor of the radar. While WindRAD is stable, the complex geometry, the design of the instrument, and the rotating antenna make the backscatter (σ∘) distributions persistently non-linear; hence NOC is insufficient. Therefore, a higher-order calibration (HOC) method is proposed. The CDF (cumulative distribution function) matching technique is employed to match the CDF of measured σ∘ instances to simulated σ∘ instances. HOC removes the non-linearities for each incidence angle. However, it is not constructed to remove the anomalous harmonic azimuth dependencies caused by the antenna rotation. These azimuth dependencies are reduced by NOCant (NOC as a function of incidence angle and relative antenna azimuth angle). Therefore, the combination of HOC and NOCant is implemented to correct both anomalous σ∘ amplitude and azimuth variations. The wind retrieval performance is evaluated with NOCant, HOC, and HOC and NOCant combined. The wind statistics and the cone distance metric both show that HOC&NOCant achieves the optimal winds for C-band and Ku-band. The calibrations have been tested on two operational input data versions; HOC works well on both data versions and HOC&NOCant can achieve the optimal wind performance for both data versions. This confirms the usefulness of HOC in the case of non-linear instrument gain anomalies. [ABSTRACT FROM AUTHOR]

The Dome Fuji (DF) region in Antarctica is a potential site for an ice core with a record of over 1 Myr. Here, we combine large-scale internal airborne radar stratigraphy with a 1-D model to estimate the age of basal ice in the DF region. The radar data used in the study were collected in a survey during the 2016–2017 Antarctic season. We transfer the latest age–depth scales from the DF ice core to isochrones traced in radargrams in the surrounding 500 km × 550 km region. At each point of the survey the 1-D model uses the ages of isochrones to construct the age–depth scale at depths where dated isochrones do not exist, the surface accumulation rate and the basal thermal condition, including melt rate and the thickness of stagnant ice. Our resulting age distribution and age density suggest that several promising sites with ice older than 1.5 Myr in the DF region might exist. The deduced melt rates and presence of stagnant ice provide more constraints for locating sites with a cold base. The accumulation rates range from 0.015 to 0.038 m a -1 ice equivalent. Based on sensitivity studies we find that the number and depth of picked isochrones and the timescale of the ice core severely affect the model results. Our study demonstrates that constraints from deep radar isochrones and a trustworthy timescale could improve the model estimation to find old ice in the DF region. [ABSTRACT FROM AUTHOR]

Stratiform rain and the overlying ice play crucial roles in the Earth's climate system. From a microphysics standpoint, water mass flux primarily depends on two variables: particles concentration and their mass. The Dual-frequency Precipitation Radar (DPR) on the Global Precipitation Measurement mission core satellite is a space-borne instrument capable of estimating these two quantities through dual-wavelength measurements. In this study, we evaluate bulk statistics on the ice particle properties derived from dual-wavelength radar data in relation to the properties of rain underneath. Specifically, we focus on DPR observations over stratiform precipitation, characterized by columns exhibiting a prominent bright band, where the melting layer can be easily detected. Our analysis reveals a significant increase in the retrieved mass flux as we transition from the ice to the rain phase in the official DPR product. This observation is in disagreement with our expectation that mass flux should remain relatively stable across the bright band in cold rain conditions. To address these discrepancies, we propose an alternative retrieval algorithm that ensures a gradual transition of Dm (mean mass-weighted particle melted-equivalent diameter) and precipitation rate across the melting zone. This approach also helps in estimating bulk ice density above the melting level. These findings demonstrate that DPR observations can not only quantify ice particle content and their size above stratiform rain regions but also estimate bulk density, provided uniform conditions that minimize uncertainties related to partial beam filling. [ABSTRACT FROM AUTHOR]

Most Arctic clouds occur below 2 km altitude, as revealed by CloudSat satellite observations. However, recent studies suggest that the relatively coarse spatial resolution, low sensitivity, and blind zone of the radar installed on CloudSat may not enable it to comprehensively document low-level clouds. We investigate the impact of these limitations on the Arctic low-level cloud fraction, which is the number of cloudy points with respect to all points as a function of height, derived from CloudSat radar observations. For this purpose, we leverage highly resolved vertical profiles of low-level cloud fraction derived from down-looking Microwave Radar/radiometer for Arctic Clouds (MiRAC) radar reflectivity measurements. MiRAC was operated during four aircraft campaigns that took place in the vicinity of Svalbard during different times of the year, covering more than 25 000 km. This allows us to study the dependence of CloudSat limitations on different synoptic and surface conditions. A forward simulator converts MiRAC measurements to synthetic CloudSat radar reflectivities. These forward simulations are compared with the original CloudSat observations for four satellite underflights to prove the suitability of our forward-simulation approach. Above CloudSat's blind zone of 1 km and below 2.5 km , the forward simulations reveal that CloudSat would overestimate the MiRAC cloud fraction over all campaigns by about 6 percentage points (pp) due to its horizontal resolution and by 12 pp due to its range resolution and underestimate it by 10 pp due to its sensitivity. Especially during cold-air outbreaks over open water, high-reflectivity clouds appear below 1.5 km , which are stretched by CloudSat's pulse length causing the forward-simulated cloud fraction to be 16 pp higher than that observed by MiRAC. The pulse length merges multilayer clouds, whereas thin low-reflectivity clouds remain undetected. Consequently, 48 % of clouds observed by MiRAC belong to multilayer clouds, which reduces by a factor of 4 for the forward-simulated CloudSat counterpart. Despite the overestimation between 1 and 2.5 km , the overall low-level cloud fraction is strongly reduced due to CloudSat's blind zone that misses a cloud fraction of 32 % and half of the total (mainly light) precipitation amount. [ABSTRACT FROM AUTHOR]

Solid precipitation measurements in Antarctica are crucial as snowfall represents the main water input term for the polar cap, and its probable increase in the coming century can mitigate sea-level rise caused by global warming. This paper presents 7 years of Micro Rain Radar (Metek MRR-2) data at the Dumont d'Urville station in coastal Adélie Land, East Antarctica. Statistics are calculated on 3 radar variables (equivalent reflectivity, mean Doppler velocity and signal-to-noise ratio) to outline the main characteristics of the radar dataset. Seasonal and interannual variabilities are also investigated, but no significant temporal trends are detected except for the seasonal mean Doppler velocity which is higher in summer and lower in winter. We then use the snowfall rate (S) data from a colocated snow-gauge to estimate the MRR precipitation profile from the radar equivalent reflectivity (Ze) through a locally derived Ze-S relation. We find the relation Ze = 43.3 S0.88. The processing method used to obtain this relation, data quality and uncertainty considerations are discussed in the manuscript. In order to give an example of application of the dataset, a brief statistical comparison of the MRR precipitation rate along the vertical with model data from the ERA5 reanalysis and the LMDZ climate model is performed, and notably shows that models underestimate heavy precipitation events. [ABSTRACT FROM AUTHOR]

Radar Doppler spectra observations provide a wealth of information about cloud and precipitation microphysics and dynamics. The interpretation of these measurements depends on our ability to simulate these observations accurately using a forward model. The effect of small-scale turbulence on the radar Doppler spectra shape has been traditionally treated by implementing the convolution process on the hydrometeor reflectivity spectrum and environmental turbulence. This approach assumes that all the particles in the radar sampling volume respond the same to turbulent-scale velocity fluctuations and neglects the particle inertial effect. Here, we investigate the inertial effects of liquid-phase particles on the forward modeled radar Doppler spectra. A physics-based simulation (PBS) is developed to demonstrate that big droplets, with large inertia, are unable to follow the rapid change of the velocity field in a turbulent environment. These findings are incorporated into a new radar Doppler spectra simulator. Comparison between the traditional and newly formulated radar Doppler spectra simulators indicates that the conventional simulator leads to an unrealistic broadening of the spectrum, especially in a strong turbulent environment. This study provides clear evidence to illustrate the droplet inertial effect on radar Doppler spectrum and develops a physics-based simulator framework to accurately emulate the Doppler spectrum for a given droplet size distribution (DSD) in a turbulence field. The proposed simulator has various potential applications for the cloud and precipitation studies, and it provides a valuable tool to decode the cloud microphysical and dynamical properties from Doppler radar observation. [ABSTRACT FROM AUTHOR]

The national German weather radar network operates in C-band between 5.6 and 5.65 GHZ. In a radar network, individual transmit frequencies have to be chosen such that radar-radar induced interferences are avoided. In an unique experiment the Hohenpeissenberg research radar and five operational systems from the radar network were used to characterize radar-radar induced interferences as a function of the radar frequency. The results allow to assess the possibility to add additional C-band radars with magnetron transmitters into the existing network. Based on the experiment, at least a 15 MHz separation of the nominal radar frequency is needed to avoid a radar-radar interference. The most efficient mitigation of radar-radar interference is achieved by the 'Radar Tango'. Latter refers to the synchronized scanning of all radar systems in the network. Based on those results, additional C-band radar systems can be added to the German weather radar network if a further improvement of the radar coverage is needed. [ABSTRACT FROM AUTHOR]

Given that the Gaussianity of observation error distribution is the fundamental principle of some data assimilation and machine learning algorithms, the error structure of radar reflectivity becomes increasingly important with the development of high resolution forecasts and nowcasts of convective systems. This study examines the error distribution of radar reflectivity and discusses what give rise to the non-Gaussian error distribution by using 6 month observations minus backgrounds (OmBs) of composites of vertical maximum reflectivity (CVMRs) in mountainous and hilly areas. By following the symmetric error model in all-sky satellite radiance assimilation, we unveil the error structure of CVMRs as a function of symmetric rain rates, which is the average of observed and simulated rain rates. Unlike satellite radiance, the error structure of CVMRs shows a sharper slope in light precipitations than moderate precipitations. Thus, a three-piecewise fitting function is more suitable for CVMRs. The probability density functions of OmBs normalized by symmetric rain rates become more Gaussian in comparison with the probability density function normalized by the whole samples. Moreover, the possibility of using third-party predictor to construct the symmetric error model are also discussed in this study. The Gaussianity of OmBs can be further improved by using a more accurate precipitation observations. According to the Jensen-Shannon divergence, a more linear predictor, the logarithmic transformation of rain rate, can provide the most Gaussian error 20 distribution in comparison with other predictors. [ABSTRACT FROM AUTHOR]

Having knowledge of the drop size distribution (DSD) is of particular interest to researchers as it is widely applied to quantitative precipitation estimation (QPE) methods. Polarimetric radar measurements have previously been utilized to derive DSD curve characteristics frequently modeled as a gamma distribution. Likewise, approaches using dual-frequency measurements have shown positive results. Both cases have relied on the need to constrain the relationship between the DSD parameters based on location or assumed weather conditions. This paper presents a methodology for retrieving the DSD parameters using the dual-frequency and polarimetric nature of measurements from a unique data set taken at co-located S-band and C-band dual-polarization radars. Using the reflectivity and differential-phase measurements from each radar, an optimization routine employing particle swarm optimization (PSO) and T-matrix computation of radar parameters is able to accurately retrieve the gamma distribution parameters without the constraints required in previous methods. Retrieved results are compared to known truth data collected using a network of OTT Parsivel disdrometers in Taiwan in order to assess the success of this procedure. [ABSTRACT FROM AUTHOR]

The Southern Ocean (SO) provides a unique natural laboratory for studying cloud formation and cloud-aerosol interactions with minimal anthropogenic influence. The Southern Ocean Clouds, Radiation, Aerosol Transport Experimental Study (SOCRATES), was an aircraft-based campaign conducted from January 15 to February 28, 2018, off the coast of Hobart, Tasmania. During SOCRATES, the NSF/NCAR GV research aircraft, equipped with in-situ probes and remote sensors, observed aerosol, cloud, and precipitation properties, and provided detailed vertical structure of clouds over the SO, particularly for the low-level clouds (below 3 km). The HIAPER Cloud Radar (HCR) and in-situ cloud and drizzle probes (CDP and 2DS) measurements were used to provide comprehensive statistical and phase-relevant macrophysical properties for the low-level clouds sampled by the 15 research flights during SOCRATES. A new method based on HCR reflectivity and spectrum width gradient was developed to estimate cloud boundaries (cloud-base and -top heights) and classify cloud types based on their top and base heights. Low-level clouds were found to be the most prevalent, with an almost 90 % occurrence frequency. A new phase determination method was also developed to identify the single-layered low-level clouds as liquid, ice, and mixed phases, with occurrence frequencies of 45.4 %, 32.5 %, and 22.2 %, respectively. Low-level clouds over the SO have significantly higher SLW concentrations, with liquid being most prevalent at higher temperatures, ice phase dominating at lower temperatures, and mixed-phase being least common due to its thermodynamic instability. Regarding their vertical distributions, the liquid phase occurs most frequently in the lower mid-cloud range (from 500 m to 1 km), the mixed phase dominates at cloud bases lower than 1 km but is well distributed along the vertical cloud layer, while the ice phase is prevalent from the middle to upper cloud levels (1–3 km). The higher occurrence of the mixed phase at the cloud base could be attributed to large drizzle-sized drops and/or ice particles. [ABSTRACT FROM AUTHOR]

Given that the Gaussianity of the observation error distribution is the fundamental principle of some data assimilation and machine learning algorithms, the error structure of radar reflectivity has become increasingly important with the development of high-resolution forecasts and nowcasts of convective systems. This study examines the error distribution of radar reflectivity and discusses what causes the non-Gaussian error distribution using 6-month observations minus backgrounds (OmBs) of composites of vertical maximum reflectivity (CVMRs) in mountainous and hilly areas. By following the symmetric error model in all-sky satellite radiance assimilation, we reveal the error structure of CVMRs as a function of symmetric rain rates, which is the average of the observed and simulated rain rates. Unlike satellite radiance, the error structure of CVMRs shows a sharper slope for light precipitation than for moderate precipitation. Thus, a three-piecewise fitting function is more suitable for CVMRs. The probability density functions of OmBs normalized by symmetric rain rates become more Gaussian than the probability density functions normalized by all samples. Moreover, the possibility of using a third-party predictor to construct the symmetric error model is also discussed in this study. The result shows that the Gaussian distribution of OmBs can be further improved via more accurate precipitation observations. According to the Jensen–Shannon divergence, a more linear predictor, the logarithmic transformation of the rain rate, can provide the most Gaussian error distribution in comparison with other predictors. [ABSTRACT FROM AUTHOR]

In this research, we delve into the influence of cloud condensation nuclei (CCN) and ice-nucleating particle (INP) concentrations on the morphology and abundance of ice particles in mixed-phase clouds, emphasizing the consequential impact of ice particle shape, number, and size on cloud dynamics and microphysics. Leveraging the synergy of the Advanced Microphysics Prediction System (AMPS) and the Kinematic Driver (KiD) model, we conducted simulations to capture cloud microphysics across diverse CCN and INP concentrations. The Passive and Active Microwave radiative TRAnsfer (PAMTRA) radar forward simulator further augmented our study, offering insights into how the concentrations of CCN and INPs affect radar reflectivities. Our experimental framework encompassed CCN concentrations ranging from 10 to 5000 cm -3 and INP concentrations from 0.001 to 10 L -1. Central to our findings is the observation that higher INP concentrations yield smaller ice particles, while an increase in CCN concentrations leads to a subtle growth in their dimensions. Consistent with existing literature, our results spotlight oblate-like crystals as dominant between temperatures of −20 and −16 °C. Notably, high-INP scenarios unveiled a significant prevalence of irregular polycrystals. The aspect ratio (AR) of ice particles exhibited a decline with the rise in both CCN and INP concentrations, highlighting the nuanced interrelation between CCN levels and ice particle shape, especially its ramifications on the riming mechanism. The forward-simulated radar reflectivities, spanning from -11.83 dBZ (low INP, 0.001 L -1) to 4.65 dBZ (high INP, 10 L -1), elucidate the complex dynamics between CCN and INPs in determining mixed-phase cloud characteristics. Comparable differences in radar reflectivity were also reported from observational studies of stratiform mixed-phase clouds in contrasting aerosol environments. Our meticulous analysis of KiD-AMPS simulation outputs, coupled with insights into aerosol-driven microphysical changes, thus underscores the significance of this study in refining our ability to understand and interpret observations and climate projections. [ABSTRACT FROM AUTHOR]