Search

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]

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]

Because of their high sensitivity to hydrometeors and their high vertical resolutions, space-borne radar observations are emerging as an undeniable asset for Numerical Weather Prediction (NWP) applications. The EUMETSAT (European Organisation for the Exploitation of Meteorological Satellites) NWP SAF (Satellite Application Facility) released an active sensor module within version 13 of the RTTOV (Radiative Transfer for TOVS) software with the goal of simulating both active and passive microwave instruments within a single framework using the same radiative transfer assumptions. This study provides an in-depth description of the radar simulator available within this software. In addition, this study proposes a revised version of the existing melting layer parametrization scheme of Bauer (2001) within the RTTOV-SCATT v13.1 model to provide a better fit to observations below the freezing level. Simulations are performed with and without melting layer schemes for the Dual precipitation radar (DPR) instrument onboard GPM using the ARPEGE (Action de Recherche Petite Echelle Grande Echelle) global NWP model running operationally at Météo-France for two different one-month periods (June, 2020 and January, 2021). Results for a case study over the Atlantic ocean show that the revised melting scheme produces more realistic simulations as compared to the default scheme both at Ku (13.5 GHz) and Ka (35.5 GHz) frequencies and these simulations are much closer to observations. A statistical assessment using more samples show significant improvement of the first-guess departure statistics with the revised scheme compared to the existing melting scheme. As a step further, this study showcases the use of melting layer simulations for the classification of precipitation (stratiform, convective and transition) using the Dual Frequency Ratio algorithm (DFR). The classification results also reveal a significant overestimation of the rain reflectivities in all hemispheres, which can either be due to a tendency of the ARPEGE model to produce a too large amount of convective precipitation, or to a mis-representation of the convective precipitation fraction within the forward operator. [ABSTRACT FROM AUTHOR]

Measurements with a 763 nm channel, located within the oxygen A-band and equipped on the Second-generation Global Imager (SGLI) onboard the JAXA's Global Change Observation Mission – Climate (GCOM-C) satellite, have the potential to retrieve cloud base height (CBH) and cloud geometric thickness (CGT) through passive remote sensing. This study implemented an algorithm to retrieve the CBH using the SGLI 763 nm channel in combination with several other SGLI channels in the visible, shortwave infrared, and thermal infrared regions. In addition to CBH, the algorithm can simultaneously retrieve other key cloud properties, including cloud optical thickness (COT), cloud effective radius, ice COT fraction as the cloud thermodynamic phase, cloud top height (CTH), and CGT. Moreover, the algorithm can be seamlessly applied to global clouds comprised of liquid, ice, and mixed phases. The SGLI-retrieved CBH exhibited quantitative consistency with CBH data obtained from the ground-based ceilometer network, ship-borne ceilometer, satellite-borne radar and lidar observations, as evidenced by sufficiently high correlations and small biases. These results provide practical evidence that the retrieval of CBH is indeed possible using the SGLI 763 nm channel. Moreover, the results lend credence to the future use of SGLI CBH data, including the estimation of the surface downward longwave radiative flux from clouds. Nevertheless, issues remain that must be addressed to enhance the value of SGLI-derived cloud retrieval products. These include the systematic bias of SGLI CTH related to cirrus clouds and the bias of SGLI CBH caused by multi-layer clouds. [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]

An experimental study aimed at identifying special rainfall regimes with the help of co-located disdrometers is performed. Eight potentially special events (i.e., four number-controlled events and four size-controlled events) are identified and examined. However, a detailed cross-check with additional, independent radar measurements reveals no clear evidence of special rainfall dynamics. The research underscores the difficulty of experimentally confirming seemingly straightforward questions about rainfall patterns and dynamics that have been theorized in the literature for several decades but never formally validated experimentally. The study also questions the reliability of previous claims and serves as a reminder to approach such problems with more caution, emphasizing the need for rigorous uncertainty analysis and multiple cross-checks between sensors to avoid misinterpretation. [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]

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]

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]

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]

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]

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]

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]

This paper introduces a machine-learning-driven approach for automated Nocturnal Low-Level Jet (NLLJ) identification using observations of wind profiles from a Radar Wind Profiler (RWP). The work discussed here is an effort to lay the groundwork for a systematic study of the Mid-Atlantic NLLJ's formation mechanisms and their influence on nocturnal and diurnal air quality in major urban regions by establishing a general framework of NLLJ features and characteristics with an identification algorithm. Leveraging a comprehensive wind profile dataset maintained by the Maryland Department of Environment's RWP network, our methodology employs supervised machine learning techniques to isolate the features of the south-westerly NLLJ. This methodology was developed to illuminate spatiotemporal patterns and nuanced characteristics of NLLJ events, unveiling their significant role in shaping the planetary boundary layer. This paper discusses the construction of this methodology, its performance against known NLLJs in the current literature, intended usage, and a preliminary statistical analysis. First light results from this analysis have identified a total of 90 south-westerly NLLJs from May–September of 2017–2021 as captured by the RWP stationed in Beltsville, MD (39.05°, -76.87°, 135 m ASL). A composite of these 90 jets is presented to better illustrate many of the bulk parameters, such as core height, duration, and maximum wind speed, associated with the onset and decay of the Mid-Atlantic NLLJ. We hope our study equips researchers and policymakers with further means to monitor, predict, and address these nocturnal dynamics phenomena that frequently influence boundary layer composition and air quality in the U.S. Mid-Atlantic and Northeastern regions. [ABSTRACT FROM AUTHOR]

Hailstorms, while fascinating from a meteorological perspective, pose significant risks to communities, agriculture, and infrastructure. In regions such as the Central Andes of Peru, the characteristics and frequency of these extreme weather events remain largely uncharted. This study fills this gap by investigating the historical frequency and vertical structure of hailstorms in this region. We analyzed historical hailstorm records dating back to 1958 alongside 4 years of observations (2017–2021) from the Parsivel2 disdrometer and a cloud-profiling radar MIRA35c. Our findings indicate a trend of decreasing hail frequency (- 0.5 events per decade). However, the p value of 0.07 suggests the need for further investigation, particularly in relation to environmental changes and reporting methods. The results show that hailstorms predominantly occur during the austral summer months, with peak frequency in December, and are most common during the afternoon and early evening hours. The analysis of radar variables such as reflectivity, radial velocity, spectral width, and linear depolarization ratio (LDR) reveals distinct vertical profiles for hail events. Two case studies highlight the diversity in the radar measurements of hailstorms, underscoring the complexity of accurate hail detection. This study suggests the need for refining the Parsivel2 algorithm and further understanding its classification of hydrometeors. Additionally, the limitations of conventional radar variables for hail detection are discussed, recommending the use of LDR and Doppler spectrum analysis for future research. Our findings lay the groundwork for the development of more efficient hail detection algorithms and improved understanding of hailstorms in the Central Andes of Peru. [ABSTRACT FROM AUTHOR]

The structure and distribution of sub-glacial water directly influences Antarctic ice mass loss by reducing or enhancing basal shear stress and accelerating grounding line retreat. A common technique for detecting sub-glacial water involves analyzing the spatial variation in reflectivity from an airborne radar echo sounding (RES) survey. Basic RES analysis exploits the high dielectric contrast between water and most other substrate materials, where a reflectivity increase ≥ 15 dB is frequently correlated with the presence of sub-glacial water. There are surprisingly few additional tools to further characterize the size, shape, or extent of hydrological systems beneath large ice masses. We adapted an existing radar backscattering simulator to model RES reflections from sub-glacial water structures using the University of Texas Institute for Geophysics (UTIG) Multifrequency Airborne Radar Sounder with Full-phase Assessment (MARFA) instrument. Our series of hypothetical simulation cases modeled water structures from 5 to 50 m wide, surrounded by bed materials of varying roughness. We compared the relative reflectivity from rounded Röthlisberger channels and specular flat canals, showing both types of channels exhibit a positive correlation between size and reflectivity. Large (> 20 m), flat canals can increase reflectivity by more than 20 dB , while equivalent Röthlisberger channels show only modest reflectivity gains of 8–13 dB. Changes in substrate roughness may also alter observed reflectivity by 3–6 dB. All of these results indicate that a sophisticated approach to RES interpretation can be useful in constraining the size and shape of sub-glacial water features. However, a highly nuanced treatment of the geometric context is necessary. Finally, we compared simulated outputs to actual reflectivity from a single RES flight line collected over Thwaites Glacier in 2022. The flight line crosses a previously proposed Röthlisberger channel route, with an obvious bright bed reflection in the radargram. Through multiple simulations comparing various water system geometries, such as canals and sub-glacial lakes, we demonstrated the important role that topography and water geometry can play in observed RES reflectivity. From the scenarios that we tested, we concluded the bright reflector from our RES flight line cannot be a Röthlisberger channel but could be consistent with a series of flat canals or a sub-glacial lake. However, we note our simulations were not exhaustive of all possible sub-glacial water configurations. The approach outlined here has broad applicability for studying the basal environment of large glaciers. We expect to apply this technique when constraining the geometry and extent of many sub-glacial hydrologic structures in the future. Further research may also include comprehensive investigations of the impact of sub-glacial roughness, substrate heterogeneity, and computational efficiencies enabling more complex and complete simulations. [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]

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]

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]

Monitoring of lake ice is important to maintain transportation routes, but in recent decades the number of in situ observations have declined. Remote sensing has worked to fill this gap in observations, with active microwave sensors, particularly synthetic-aperture radar (SAR), being a crucial technology. However, the impact of wet conditions on radar and how interactions change under these conditions have been largely ignored. It is important to understand these interactions as warming conditions are likely to lead to an increase in the occurrence of slush layers. This study works to address this gap using the Snow Microwave Radiative Transfer (SMRT) model to conduct forward-modelling experiments of backscatter for Lake Oulujärvi in Finland. Experiments were conducted under dry conditions, under moderate wet conditions, and under saturated conditions. These experiments reflected field observations during the 2020–2021 ice season. Results of the dry-snow experiments support the dominance of surface scattering from the ice–water interface. However, conditions where layers of wet snow are introduced show that the primary scattering interface changes depending on the location of the wet layer. The addition of a saturated layer at the ice surface results in the highest backscatter values due to the larger dielectric contrast created between the overlying dry snow and the slush layer. Improving the representation of these conditions in SMRT can also aid in more accurate retrievals of lake ice properties such as roughness, which is key for inversion modelling of other properties such as ice thickness. [ABSTRACT FROM AUTHOR]

NASA's Vapor In-cloud Profiling Radar (VIPR) is a tunable G-band radar designed for in-cloud and precipitation humidity remote sensing. VIPR estimates humidity using the differential absorption radar (DAR) technique. This technique exploits the difference between atmospheric attenuation at different frequencies ("on" and "off" an absorption line) and combines it with the ranging capabilities of the radar to estimate the absorbing gas concentration along the radar path. We analyze the VIPR humidity measurements during two NASA field campaigns: (1) the Investigation of Microphysics and Precipitation for Atlantic Coast-Threatening Snowstorms (IMPACTS) campaign, with the objective of studying wintertime snowstorms focusing on east coast cyclones; and (2) the Synergies Of Active optical and Active microwave Remote Sensing Experiment (SOA 2 RSE) campaign, which studied the synergy between DAR (VIPR) and differential absorption lidar (DIAL, the High altitude Lidar Observatory – HALO) measurements. We discuss a comparison with dropsondes launched during these campaigns as well as an intercomparison against the ERA5 reanalysis fields. Thus, this study serves as an additional evaluation of ERA5 lower tropospheric humidity fields. Overall, in-cloud and in-snowstorm comparisons suggest that ERA5 and VIPR agree within 20 % or better against the dropsondes. The exception is during SOA 2 RSE (i.e., in fair weather), where ERA5 exhibits up to a 50 % underestimation above 4 km. We also show a smooth transition in water vapor profiles between the in-cloud and clear-sky measurements obtained from VIPR and HALO respectively, which highlights the complementary nature of these two measurement techniques for future airborne and space-based missions. [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]

Large hail events are typically infrequent, with significant time gaps between occurrences at specific locations. However, when these events do happen, they can cause rapid and substantial economic losses within a matter of minutes. Therefore, it is crucial to have the ability to accurately observe and understand hail phenomena to improve the mitigation of this impact. While in situ observations are accurate, they are limited in number for an individual storm. Weather radars, on the other hand, provide a larger observation footprint, but current radar-derived hail size estimates exhibit low accuracy due to horizontal advection of hailstones as they fall, the variability of hail size distributions (HSDs), complex scattering and attenuation, and mixed hydrometeor types. In this paper, we propose a new radar-derived hail product developed using a large dataset of hail damage insurance claims and radar observations. We use these datasets coupled with environmental information to calculate a hail damage estimate (HDE) using a deep neural network approach aiming to quantify hail impact, with a critical success index of 0.88 and a coefficient of determination against observed damage of 0.79. Furthermore, we compared HDE to a popular hail size product (MESH), allowing us to identify meteorological conditions that are associated with biases on MESH. Environments with relatively low specific humidity, high CAPE and CIN, low wind speeds aloft, and southerly winds at the ground are associated with a negative MESH bias, potentially due to differences in HSD, hail hardness, or mixed hydrometeors. In contrast, environments with low CAPE, high CIN, and relatively high specific humidity aloft are associated with a positive MESH bias. [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]

Polar mesosphere summer echoes (PMSE) have been observed for more than 30 years with 50 MHz VHF radars at various locations in the Northern Hemisphere. Continuous observation of PMSE is conducted on the northern Norwegian island of Andøya (69.3° N) using the ALWIN radar (1999–2008) and MAARSY (since 2010). The same kind of PMSE measurements began in 2004 in the Southern Hemisphere with the Australian Antarctic Division's VHF radar at Davis Station in Antarctica (68.6° S), which is at an opposite latitude to Andøya. Since the radars at both sites are calibrated, the received echo strength of PMSE from more than 1 decade of mesospheric observations on both hemispheres could be converted to absolute signal power, allowing for direct comparison of the measurements. Comparison of PMSE observations obtained at both radar sites during a period of 23 boreal summers (Andøya) and 15 austral summers (Davis) shows that their PMSE signal strengths are of the same order of magnitude, but significantly fewer PMSE are observed in the Southern Hemisphere than in the Northern Hemisphere. Compared to Andøya, the PMSE season over Davis starts about 7 d later on average and ends 9 d earlier, making it 16 d shorter. PMSE over Davis occur less frequently but with greater variability in seasonal, diurnal, and altitudinal occurrence. For example, PMSE over Davis reach maximum altitudes about 1.5 km higher than those over Andøya. [ABSTRACT FROM AUTHOR]

Solar tides play a major role in the dynamics of the upper mesosphere and lower thermosphere (MLT). Hence, a comprehensive understanding of these tides is important for successful modelling of the MLT region. Most ground-based observations of tidal variations in the MLT have been from meteor radar measurements with a temporal resolution of 1 h. Here, we take a different perspective on these tidal variations using high-resolution 1 min neutral-wind measurements from the Falkland Islands SuperDARN (Super Dual Auroral Radar Network) radar. This analysis shows that these higher-resolution data can be used to identify higher frequency tidal components than are typically observed by meteor radars (up to a heptadiurnal component). It also shows evidence of significant power in these higher frequency components, particularly in the quaddiurnal component, which may be particularly suitable for a global analysis using high-resolution SuperDARN neutral-wind measurements. The high-resolution analysis also shows evidence of fluctuations with a frequency of 1.5 cycles per day, as well as higher-frequency fluctuations, accompanying a quasi-2 d (two day) wave. We discuss the limitations of this high-resolution analysis method and the new opportunities that it may provide. We conclude that higher-resolution SuperDARN neutral-wind measurements need to be better exploited in the future as they provide a complementary way of studying tides and waves in the MLT. [ABSTRACT FROM AUTHOR]

The Goddard Profiling algorithm (GPROF) converts radiometer observations from Global Precipitation Measurement (GPM) constellation satellites into precipitation estimates. Typically, high-quality ground-based estimates serve as reference to evaluate GPROF's performance. To provide a fair comparison, the ground-based estimates are often spatially aligned to GPROF. However, GPROF combines observations from various sensors and channels, each associated with a distinct footprint. Consequently, uncertainties related to the representativeness of the sampled areas are introduced in addition to the uncertainty when converting brightness temperatures into precipitation intensities. The exact contribution of resampling precipitation estimates, required to spatially and temporally align different resolutions when combining or comparing precipitation observations, to the overall uncertainty remains unknown. Here, we analyze the current performance of GPROF over the Netherlands during a 4-year period (2017–2020) while investigating the uncertainty related to sampling. The latter is done by simulating the reference precipitation as satellite footprints that vary in size, geometry, and applied weighting technique. Only GPROF estimates based on observations from the conical-scanning radiometers of the GPM constellation are used. The reference estimates are gauge-adjusted radar precipitation estimates from two ground-based weather radars from the Royal Netherlands Meteorological Institute (KNMI). Echo top heights (ETHs) retrieved from the same radars are used to classify the precipitation as shallow, medium, or deep. Spatial averaging methods (Gaussian weighting vs. arithmetic mean) minimally affect the magnitude of the precipitation estimates. Footprint size has a higher impact but cannot explain all discrepancies between the ground- and satellite-based estimates. Additionally, the discrepancies between GPROF and the reference are largest for low ETHs, while the relative bias between the different footprint sizes and implemented weighting methods increase with increasing ETHs. Lastly, our results do not show a clear difference between coastal and land simulations. We conclude that the uncertainty introduced by merging different channels and sensors cannot fully explain the discrepancies between satellite- and ground-based precipitation estimates. Hence, uncertainties related to the retrieval algorithm and environmental conditions are found to be more prominent than resampling uncertainties, in particular for shallow and light precipitation. [ABSTRACT FROM AUTHOR]

Observations of the vertical wind profile in Chongqing, a typical mountainous city in China, are important, but they are sparse and have low resolution. To obtain more wind profile data, this study matched the Aeolus track with ground-based wind observation sites in Chongqing in 2021. Based on the obtained results, verification and quality control studies were conducted on the wind observations of a wind profile radar (WPR) with radiosonde (RS) data, and a comparison of the Aeolus Mie-cloudy and Rayleigh-clear wind products (Aeolus winds measured in cloudy and aerosol-rich atmospheric conditions from Mie-channel-collected data and winds measured in clear-air conditions from Rayleigh-collected data) with WPR data was then performed. The conclusions can be summarized as follows: (1) a clear correlation between the wind observations of WPR and RS was found, with a correlation coefficient (R) of 0.71. Their root mean square deviation increased with height but decreased at heights between 3 and 4 km. (2) After quality control using Gaussian filtering (GF) and empirical orthogonal function construction (EOFc; G=87.23 %) of the WPR data, the R between the WPR and RS reached 0.83 and 0.95, respectively. The vertical distribution showed that GF could better retain the characteristics of WPR wind observations but with limited improvement in decreasing deviations, whereas EOFc performed better in decreasing deviations but considerably modified the original characteristics of the wind field, especially regarding intensive vertical wind shear in strong convective weather processes. (3) In terms of the differences between the Aeolus and WPR data, 56.0 % and 67.8 % deviations were observed within ±5 m s -1 for Rayleigh-clear and Mie-cloudy winds (Aeolus winds measured in cloudy and aerosol-rich atmospheric conditions from Mie-channel-collected data and winds measured in clear-air conditions from Rayleigh-collected data) vs WPR winds, respectively. Vertically, large mean differences of both Rayleigh-clean and Mie-cloudy winds versus WPR winds appeared below 1.5 km, which is attributed to the prevailing quiet and small winds within the boundary layer in Chongqing; in this case the movement of molecules and aerosols is mostly affected by irregular turbulence. Additionally, large mean differences at a height range between 4 and 8 km for Mie-cloudy versus WPR winds may be related to the high content of cloud liquid water in the middle troposphere of Chongqing. (4) The differences in both Rayleigh-clear and Mie-cloudy versus WPR winds had changed. Deviations of 58.9 % and 59.6 % were concentrated within ±5 m s -1 for Rayleigh-clear versus WPR winds with GF and EOFc quality control, respectively. In contrast, 69.1 % and 70.2 % of deviations appeared within ±5 m s -1 for Rayleigh-clear versus WPR and EOFc WPR winds, respectively. These results shed light on the comprehensive applications of multi-source wind profile data in mountainous cities or areas with sparse ground-based wind observations. [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]

Precipitation nowcasting has important implications for urban operation and flood prevention. Radar echo extrapolation is a common method in precipitation nowcasting. Using deep learning models to extrapolate radar echo data has great potential. The increase of lead time leads to a weaker correlation between the real rainfall evolution and the generated images. The evolution information is easily lost during extrapolation, which is reflected as echo attenuation. Existing models, including generative adversarial network (GAN)-based models, have difficulty curbing attenuation, resulting in insufficient accuracy in rainfall prediction. To solve this issue, a spatiotemporal process enhancement network (GAN-argcPredNet v2.0) based on GAN-argcPredNet v1.0 has been designed. GAN-argcPredNet v2.0 curbs attenuation by avoiding blurring or maintaining the intensity. A spatiotemporal image correlation (STIC) prediction network is designed as the generator. By suppressing the blurring effect of rain distribution and reducing the negative bias by STIC attention, the generator generates more accurate images. Furthermore, the discriminator is a channel–spatial (CS) convolution network. The discriminator enhances the discrimination of echo information and provides better guidance to the generator in image generation by CS attention. The experiments are based on the radar dataset of southern China. The results show that GAN-argcPredNet v2.0 performs better than other models. In heavy rainfall prediction, compared with the baseline, the probability of detection (POD), the critical success index (CSI), the Heidke skill score (HSS) and bias score increase by 18.8 %, 17.0 %, 17.2 % and 26.3 %, respectively. The false alarm ratio (FAR) decreases by 3.0 %. [ABSTRACT FROM AUTHOR]

We present a comprehensive quality-controlled 15-month dataset of remote sensing observations of low-level mixed-phase clouds (LLMPCs) taken at the high Arctic site of Ny-Ålesund, Svalbard, Norway. LLMPCs occur frequently in the Arctic region and extensively affect the energy budget. However, our understanding of the ice microphysical processes taking place in these clouds is incomplete. The dual-wavelength and polarimetric Doppler cloud radar observations, which are the cornerstones of the dataset, provide valuable fingerprints of ice microphysical processes, and the high number of cases included allows for the compiling of robust statistics for process studies. The radar data are complemented with thermodynamic retrievals from a microwave radiometer, liquid base height from a ceilometer, and wind fields from large-eddy simulations. All data are quality controlled, especially the cloud radar data, which are accurately calibrated, matched, and corrected for gas and liquid-hydrometeor attenuation, ground clutter, and range folding. We finally present an analysis of the temperature dependence of Doppler, dual-wavelength, and polarimetric radar variables, to illustrate how the dataset can be used for cloud microphysical studies. The dataset has been published in and is freely available at: 10.5281/zenodo.7803064. [ABSTRACT FROM AUTHOR]

A new dual-frequency Stratospheric Tropospheric (ST) / Meteor radar has been built and installed at the Langfang Observatory in northern China. It utilizes a novel 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-hour typical of most meteor radars. The root mean square uncertainty of the ST wind measurements is better than 2 m/s when estimating the line of best fit with radiosonde winds. Preliminary observation results for one 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]

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]

Weather radars are increasingly being used to study the interaction between wildfires and the atmosphere, owing to the enhanced spatio-temporal resolution of radar data compared to conventional measurements, such as satellite imagery and in situ sensing. An important requirement for the continued proliferation of radar data for this application is the automatic identification of fire-generated particle returns (pyrometeors) from a scene containing a diverse range of echo sources, including clear air, ground and sea clutter, and precipitation. The classification of such particles is a challenging problem for common image segmentation approaches (e.g. fuzzy logic or unsupervised machine learning) due to the strong overlap in radar variable distributions between each echo type. Here, we propose the following two-step method to address these challenges: (1) the introduction of secondary, texture-based fields, calculated using statistical properties of gray-level co-occurrence matrices (GLCMs), and (2) a Gaussian mixture model (GMM), used to classify echo sources by combining radar variables with texture-based fields from (1). Importantly, we retain all information from the original measurements by performing calculations in the radar's native spherical coordinate system and introduce a range-varying-window methodology for our GLCM calculations to avoid range-dependent biases. We show that our method can accurately classify pyrometeors' plumes, clear air, sea clutter, and precipitation using radar data from recent wildfire events in Australia and find that the contrast of the radar correlation coefficient is the most skilful variable for the classification. The technique we propose enables the automated detection of pyrometeors' plumes from operational weather radar networks, which may be used by fire agencies for emergency management purposes or by scientists for case study analyses or historical-event identification. [ABSTRACT FROM AUTHOR]

Wind turbines negatively affect the performance of weather radars, especially when located in the proximity of a radar site. In March 2019, MeteoSwiss performed a measurement campaign by deploying a mobile X-band radar in Schaffhausen. It proved to be useful for mapping and characterizing the maximum power returns by three wind turbines observed using standard scanning strategies. In March 2020, the campaign was repeated using a more sophisticated scan strategy: ∼ 100 min special sessions of fixed pointing an antenna towards the nacelle of the closest wind turbine (WT) located within a range of 7766 m from the radar, interleaved every 2 h by a scanning protocol identical to that of the March 2019 campaign. Polarimetric radar signatures were derived every 64 ms using 128 radar pulses transmitted every 0.5 ms (pulse repetition frequency (PRF) = 2000 Hz). A thorough overview of the polarimetric signatures of the WT in still or quasi-still conditions has been obtained based on 30 000 polarimetric measurables acquired over 32 min on the first day of the campaign (4 March 2020). During the first 2 min with zero rotor speed, the co-polar correlation coefficient between the orthogonal polarization states, ρHV , was persistently equal to 1, similarly to the signature of a bright scatterer observed by a non-rotating antenna. The changes between two consecutive values of the differential reflectivity and radar reflectivity factor were either 0 dBz or ± 0.5 dBz. Due to the absence of precipitation, one could assume that the standard deviation of the differential phase shift, which was as small as 3.0 ∘ , can be entirely attributed to the variability of the differential backscattering phase shift. There were two 10 min periods during which the rotor moved less than 1 revolution. It is worth noting that this slow movement could be associated with a change in the blade pitch angle and the nacelle orientation, which caused extreme changes in the radar reflectivity factor. For instance, two pairs of 64 ms consecutive values reached 78.5 dBz, which is the absolute maximum reached in the whole campaign (4–21 March 2020). [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]

Winds from two wind-profiling radars, ESRAD (ESrange atmospheric RADar) in Arctic Sweden and MARA (Moveable Atmospheric Radar for Antarctica) on the coast of Antarctica, are compared with collocated (within 100 km) winds measured by the Doppler lidar on board the Aeolus satellite for the time period July 2019–May 2021 (baseline 2B11). Data are considered as a whole and subdivided into summer and winter as well as ascending (afternoon) and descending (morning) passes. Mean differences (bias) and random differences are categorized (standard deviation and scaled median absolute deviation) and the effects of different quality criteria applied to the data are assessed, including the introduction of the "modified Z score" to eliminate gross errors. This last criterion has a substantial effect on the standard deviation, particularly for Mie winds. Significant bias is found in two cases, for Rayleigh winds for the descending satellite passes. at MARA (-1.4 (+ 0.7) m s -1) and for all Mie winds at ESRAD (+ 1.0 (+ 0.3) m s -1). For the Rayleigh winds at MARA, there is no obvious explanation for the bias in the data distribution. The Mie wind error with respect to the wind data measured at ESRAD shows a skewed distribution toward positive values (Aeolus horizontal line-of-sight wind > ESRAD wind). Random differences (scaled median absolute deviation) for all data together are 5.9 and 5.3 m s -1 for Rayleigh winds at MARA and ESRAD, respectively, and 4.9 and 3.9 m s -1 for Mie winds. When the comparison is restricted to Aeolus measurements with a mean location within 25 km from the radars, there is no change to the random differences for Rayleigh winds, but for Mie winds they are reduced to 3.3 and 3.6 m s -1. These represent an upper bound for Aeolus wind random errors since they are due to a combination of spatial differences and random errors in both radar winds and Aeolus winds. The random errors in radar winds are < 2 m s -1 and therefore contribute little, but spatial variability clearly makes a significant contribution for Mie winds, especially at MARA. [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]

This paper explores the possibility of using multi-source precipitation estimates for climatological applications. A data-processing algorithm (RainGRS Clim) has been developed to work on precipitation accumulations such as daily or monthly totals, which are significantly longer than operational accumulations (generally between 5 min and 1 h). The algorithm makes the most of additional opportunities, such as the possibility of complementing data with delayed data, access to high-quality data that are not operationally available, and the greater efficiency of the algorithms for data quality control and merging with longer accumulations. Verification of the developed algorithms was carried out using monthly accumulations through comparison with precipitation from manual rain gauges. As a result, monthly accumulations estimated by RainGRS Clim were found to be significantly more reliable than accumulations generated operationally. This improvement is particularly noticeable for the winter months, when precipitation estimation is much more difficult due to less reliable radar estimates. [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]

The WPR-LQ-7 is a UHF (1.3575 GHz) wind profiler radar used for routine measurements of the lower troposphere at Shigaraki MU Observatory (34.85 ∘ N, 136.10 ∘ E; Japan) at a vertical resolution of 100 m and a time resolution of 10 min. Following studies carried out with the 46.5 MHz middle and upper atmosphere (MU) radar (Luce et al., 2018), we tested models used to estimate the rate of turbulence kinetic energy (TKE) dissipation ε from the Doppler spectral width in the altitude range ∼ 0.7 to 4.0 km above sea level (a.s.l.). For this purpose, we compared LQ-7-derived ε using processed data available online (http://www.rish.kyoto-u.ac.jp/radar-group/blr/shigaraki/data/ , last access: 24 July 2023) with direct estimates of ε (εU) from DataHawk UAVs. The statistical results reveal the same trends as reported by Luce et al. (2018) with the MU radar, namely (1) the simple formulation based on dimensional analysis εLout=σ3/Lout , with Lout∼70 m, provides the best statistical agreement with εU. (2) The model εN predicting a σ2N law (N is Brunt–Vaïsälä frequency) for stably stratified conditions tends to overestimate for εU≲5×10-4 m 2 s -3 and to underestimate for εU≳5×10-4 m 2 s -3. We also tested a model εS predicting a σ2S law (S is the vertical shear of horizontal wind) supposed to be valid for low Richardson numbers (Ri=N2/S2). From the case study of a turbulent layer produced by a Kelvin–Helmholtz (K–H) instability, we found that εS and εLout are both very consistent with εU , while εN underestimates εU in the core of the turbulent layer where N is minimum. We also applied the Thorpe method from data collected from a nearly simultaneous radiosonde and tested an alternative interpretation of the Thorpe length in terms of the Corrsin length scale defined for weakly stratified turbulence. A statistical analysis showed that εS also provides better statistical agreement with εU and is much less biased than εN. Combining estimates of N and shear from DataHawk and radar data, respectively, a rough estimate of the Richardson number at a vertical resolution of 100 m (Ri100) was obtained. We performed a statistical analysis on the Ri dependence of the models. The main outcome is that εS compares well with εU for low Ri100 (Ri100≲1), while εN fails. εLout varies as εS with Ri100 , so that εLout remains the best (and simplest) model in the absence of information on Ri. Also, σ appears to vary as Ri100-1/2 when Ri100≳0.4 and shows a degree of dependence on S100 (vertical shear at a vertical resolution of 100 m) otherwise. [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]

Large hail events are typically infrequent, with significant time gaps between occurrences at specific locations. However, when these events do happen, they can cause rapid and substantial economic losses within a matter of minutes. Therefore, it is crucial to have the ability to accurately observe and understand hail phenomena to improve the mitigation of this impact. While in-situ observations are accurate, they are limited in number for an individual storm. Weather radars, on the other hand, provide a larger observation footprint, but current radar-derived hail size estimates exhibit low accuracy due to horizontal advection of hailstones as they fall, the variability of hail size distributions (HSD), complex scattering and attenuation, and mixed hydrometeor types. In this paper, we propose a new radar-derived hail product that is developed using a large dataset of hail damage insurance claims and radar observations. We use these datasets coupled with environmental information to calculate a Hail Damage Estimate (HDE) using a deep neural network approach aiming to quantify hail impact, with a critical success index of 0.88 and a coefficient of determination against observed damage of 0.78. Furthermore, we compared HDE to a popular hail size product (MESH), allowing us to identify meteorological conditions that are associated with biases on MESH. Environments with relatively low specific humidity, high CAPE and CIN, low wind speeds aloft and southerly winds at ground are associated with a negative MESH bias, potentially due to differences in HSD or mixed hydrometeors. In contrast, environments with low CAPE, high CIN, and relatively high specific humidity aloft are associated with a positive MESH bias. [ABSTRACT FROM AUTHOR]

Observations of vertical wind in Chongqing, a typical mountainous city in China, are important, but sparse and have low resolution. To obtain more wind profile data, this study matched the Aeolus track with ground-based wind observation sites in Chongqing in 2021. Based on the obtained results, verification and quality control studies were conducted on the wind observations of a wind profile radar (WPR) with radiosonde (RS) data, and a comparison of the Aeolus Mie-cloudy and Rayleigh-clear wind products with WPR data was then performed. The conclusions can be summarized as follows: (1) A clear correlation between the wind observations of WPR and RS was found, with a correlation coefficient (R) of 69.92%. Their root-mean-square deviation increased with height but decreased by 3 - 4 km. (2) After quality control of Gaussian filtering (GF) and empirical orthogonal function construction (EOFc, G = 87.23%) of the WPR data, the R between the WPR and RS reached 76.00% and 95.44%, respectively. The vertical distribution showed that GF could better retain the characteristics of WPR wind observations, but with limited improvement in decreasing deviations, whereas EOFc performed better in decreasing deviations, but considerably modified the original characteristics of the wind field, especially regarding intensive vertical wind shear in strong convective weather processes. (3) In terms of the differences between the Aeolus and WPR data, 56.0% and 67.8% deviations were observed between ± 5 m/s for Rayleigh-clear and Mie-cloudy winds vs. WPR winds, respectively. Vertically, the mean differences of both Rayleigh-clean and Mie-cloudy winds versus WPR winds appeared below 1.5 km, which is attributed to the prevailing quiet and small winds within the boundary layer in Chongqing, such that the movement of molecules and aerosols is mostly affected by irregular turbulence. Additionally, large mean differences of 4-8 km for Mie-cloudy versus WPR winds may be related to the high content of cloud liquid water in the middle troposphere, influenced by the topography of the Tibetan Plateau. (4) The differences in both Rayleigh-clear and Mie-cloudy versus WPR winds had changed. Deviations of 58.9% and 59.6% were concentrated between ±5m/s for Rayleigh-clear versus WPR winds with GF and EOFc quality control, respectively. In contrast, 69.1% and 70.2% of deviations appeared between ± 5 m/s for Rayleigh-clear versus WPR and EOFc WPR winds, respectively. These results shed light on the comprehensive applications of multi-source wind profile data in mountainous cities or areas with sparse ground-based wind observations. [ABSTRACT FROM AUTHOR]

Sea ice volume's significant interannual variability requires long-term series of observations to identify trends in its evolution. Despite improvements in sea ice thickness estimations from altimetry during the past few years thanks to CryoSat-2 and ICESat-2, former ESA radar altimetry missions such as the Environmental Satellite (Envisat) and especially the European Remote-Sensing Satellite (ERS-1 and ERS-2) have remained under-exploited so far. Although solutions have already been proposed to ensure continuity of measurements between CryoSat-2 and Envisat, there is no time series integrating ERS. The purpose of this study is to extend the Arctic radar freeboard time series back to 1995. The difficulty in handling ERS measurements comes from a technical issue known as the pulse blurring effect, altering the radar echoes over sea ice and the resulting surface height estimates. Here we present and apply a correction for this pulse blurring effect. To ensure consistency of the CryoSat-2, Envisat and ERS-2 time series, a multiparameter neural-network-based method to calibrate Envisat against CryoSat-2 and ERS-2 against Envisat is presented. The calibration is trained on the discrepancies observed between the altimeter measurements during the mission-overlap periods and a set of parameters characterizing the sea ice state. Monthly radar freeboards are provided with uncertainty estimations based on a Monte Carlo approach to propagate the uncertainties all along the processing chain, including the neural network. Comparisons of corrected radar freeboards during overlap periods reveal good agreement between the missions, with a mean bias of 0.30 cm and a standard deviation of 9.7 cm for Envisat and CryoSat-2 and a 0.20 cm bias and a standard deviation of 3.8 cm for ERS-2 and Envisat. The monthly corrected radar freeboards obtained from Envisat and ERS-2 are then validated by comparison with several independent datasets such as airborne, mooring, direct-measurement and other altimeter products. Except for two datasets, comparisons lead to correlations ranging from 0.41 to 0.94 for Envisat and from 0.60 to 0.74 for ERS-2. The study finally provides radar freeboard estimation for winters from 1995 to 2021 (from the ERS-2 mission to CryoSat-2). [ABSTRACT FROM AUTHOR]

The cloud-profiling radar (CPR) on the Earth Clouds, Aerosol, and Radiation Explorer (EarthCARE) satellite (EC-CPR) is the first satellite-borne Doppler radar. In a previous study, we examined the effects of horizontal (along-track) integration and simple unfolding methods on the reduction of Doppler errors in the EC-CPR observations, and those effects were evaluated using two limited scenes in limited-latitude and low-pulse-repetition-frequency (PRF) settings. In this study, the amount of data used was significantly increased, and the area of the data used was extended globally. Not only low-PRF but also high-PRF settings were examined. We calculated the EC-CPR-observed Doppler velocity from pulse-pair covariances using the radar reflectivity factor and Doppler velocity obtained from a satellite data simulator and a global storm-resolving simulation. The global data were divided into five latitudinal zones, and each standard deviation of Doppler errors for 5 dB Ze after 10 km integration was calculated. In the case of the low-PRF setting, the error without unfolding correction for the tropics reached a maximum of 2.2 m s -1 and then decreased toward the poles (0.43 m s -1). The error with unfolding correction for the tropics became much smaller at 0.63 m s -1. In the case of the high-PRF setting, the error without unfolding correction for the tropics reached a maximum of 0.78 m s -1 and then decreased toward the poles (0.19 m s -1). The error with unfolding correction for the tropics was 0.29 m s -1 , less than half the value without the correction. The results of the analyses of the simulated data indicated that the zonal mean frequency of precipitation echoes was highest in the tropics and decreased toward the poles. Considering a limitation of the unfolding correction for discrimination between large upward velocity and large precipitation falling velocity, the latitudinal variation in the standard deviation of Doppler error can be explained by the precipitation echo distribution. [ABSTRACT FROM AUTHOR]

Convective wind gusts (CGs) are usually related to thunderstorms, and they may cause great structural damage and serious hazards, such as train derailment, service interruption, and building collapse. Due to the small-scale and nonstationary nature of CGs, reliable CG nowcasting with high spatial and temporal resolutions has remained unattainable. In this study, a novel nowcasting model based on deep learning – namely, CGsNet – is developed for 0–2 h lead times of quantitative CG nowcasting, achieving minute–kilometer-level forecasts. CGsNet is a physics-constrained model established by training on large corpora of average surface wind speed (ASWS) and radar observations; it can produce realistic and spatiotemporally consistent ASWS predictions in CG events. By combining the gust factor (1.77, the ratio of the observed peak wind gust speed (PWGS) to the ASWS) with the ASWS predictions, the PWGS forecasts are estimated with a spatial resolution of 0.01 ∘ × 0.01 ∘ and a 6 min temporal resolution. CGsNet is shown to be effective, and it has an essential advantage in learning the spatiotemporal features of CGs. In addition, quantitative evaluation experiments indicate that CGsNet exhibits higher generalization performance for CGs than the traditional nowcasting method based on numerical weather prediction models. CG-nowcasting technology can be applied to provide real-time quantitative CG forecasts. [ABSTRACT FROM AUTHOR]

We present a comprehensive quality-controlled 15-month dataset of remote sensing observations of low-level mixedphase clouds (LLMPCs) taken at the high Arctic site of Ny-Ålesund, Svalbard, Norway. LLMPCs occur frequently in the Arctic region, and extensively affect the energy budget. However, our understanding of the ice microphysical processes taking place in these clouds is incomplete. The dual-wavelength and polarimetric Doppler cloud radar observations, which are the cornerstones of the dataset, provide valuable fingerprints of icemicrophysical processes, and the high number of cases included allows for the compiling of robust statistics for process studies. The radar data are complemented with thermodynamic retrievals from a microwave radiometer, liquid base height from a ceilometer and wind fields from large-eddy simulations. All data are quality controlled, especially the cloud radar data, which are accurately calibrated, matched, and corrected for gas and liquid-hydrometeor attenuation, ground clutter and range folding. We finally present an analysis of the temperature dependence of Doppler, dual wavelength, and polarimetric radar variables, to illustrate how the dataset can be used for cloud microphysical studies. The dataset has been published in Chellini et al. (2023) and is freely available at: www.doi.org/10.5281/zenodo.7803064. [ABSTRACT FROM AUTHOR]

Given that the Gaussianity of observation error distribution is the fundamental principle of most current modern data assimilation methods, the error structure of radar reflectivity becomes increasingly important with the development of reflectivity assimilation in convection-allowing numerical weather prediction. This study examines the error distribution of radar reflectivity and discusses what give rise to the non-Gaussian error distribution by using 6 month reflectivity departures between observations and simulations in the Southwest China. By following the symmetric error model in all-sky satellite radiance assimilation, we unveil the error structure of radar reflectivity as a function of symmetric rainrates, which is the average of observed and simulated rainrates. Unlike satellite radiance, the reflectivity error shows a sharper slope in light precipitations than moderate precipitations. Thus, a three-piecewise fitting function is more suitable for radar reflectivity than a two-piecewise fitting function. The probability distribution functions of reflectivity departures normalized by symmetric rainrates become more Gaussian in comparison with the raw probability distribution function. Moreover, the possibility of using third-party predictor to construct the symmetric error model are also discussed in this study. According to the Jensen-Shannon divergence, a more linear predictor, the logarithmic transformation of rainrate, can provide the most Gaussian error distribution. [ABSTRACT FROM AUTHOR]