Search

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 planetary boundary layer (PBL) over the Tibetan Plateau (TP) exerts a significant influence on regional and global climate, while its vertical structures of turbulence and evolution features remain poorly understood, largely due to the scarcity of observation. This study examines the vertical profile and daytime variation of turbulence dissipation rate (ε) in the PBL over the TP using the high-resolution (6 min and 120 m) measurements from the radar wind profiler (RWP) network, combined with the hourly data from the ERA5 reanalysis. Observational analyses show that the magnitude of ε below 3 km under all-sky conditions exhibits large spatial discrepancy over the six RWP sites over the TP. Particularly, the values of ε at Minfeng and Jiuquan over the northern TP and Dingri over the southern TP are roughly an order of magnitude greater than those at Lijiang, Ganzi and Hongyuan over the eastern TP. This could be partially attributed to the difference of land cover across the six RWP sites. In terms of the diurnal variation, ε rapidly intensifies from 0900 local standard time (LST) to 1400 LST, and then gradually levels off in the late afternoon. Under clear-sky conditions, both ε and planetary boundary layer height (zi) are greater, compared with cloudy-sky conditions. This reveals that clouds would suppress the turbulence development and deduce zi. In the lower PBL (0.2< z / zi <0.5, where z is the height above ground level), the dominant influential factor for the development of turbulence is the surface-air temperature difference (Ts – Ta). By comparison, in the upper PBL (0.6< z / zi <1.0), both the and vertical wind shear (VWS) affect the development of turbulence. Above the PBL (1.0< z / zi <2.0), the shear production resulting from VWS dominates the variation of turbulence. Under cloudy-sky conditions, clouds are found to decrease the surface total solar radiation, thereby reducing Ts – Ta and surface sensible heat flux. This weakened sensible heat flux tends to inhibit the turbulent motion within PBL especially in the lower PBL and decrease the growth rate of zi. On the other hand, the strong VWS induced by clouds enhances the turbulence above the PBL. The findings obtained here underscore the importance of RWP network in revealing the fine-scale structures of the PBL over the TP and gaining new insight into the PBL evolution. [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]

Long-term observations of deep convective cloud (DCC) vertical velocity and mass flux were collected during the Observations and Modelling of the Green Ocean Amazon (GoAmazon2014/5) experiment. Precipitation echoes from a surveillance weather radar near Manaus, Brazil, are tracked to identify and evaluate the isolated DCC lifecycle evolution during the dry and wet seasons. A radar wind profiler (RWP) provides precipitation and air motion profiles to estimate the vertical velocity, mass flux, and mass transport rates within overpassing DCC cores as a function of the tracked cell lifecycle stage. The average radar reflectivity factor (Z), DCC area (A), and surface rainfall rate (R) increased with DCC lifetime as convective cells were developing, reached a peak as the cells matured, and decreased thereafter as cells dissipated. As the convective cells mature, cumulative DCC properties exhibit stronger updraft behaviors with higher upward mass flux and transport rates above the melting layer (compared with initial and later lifecycle stages). In comparison, developing DCCs have the lowest Z associated with weak updrafts, as well as negative mass flux and transport rates above the melting layer. Over the DCC lifetime, the height of the maximum downward mass flux decreased, whereas the height of the maximum net mass flux increased. During the dry season, the tracked DCCs had higher Z , propagation speed, and DCC area, and were more isolated spatially compared with the wet season. Dry season DCCs exhibit higher Z , mass flux, and mass transport rate while developing, whereas wet season DCCs exhibit higher Z , mass flux, and mass transport rates at later stages. [ABSTRACT FROM AUTHOR]

Reliable estimation of precipitation fields at high resolution is a key issue for snow cover modelling in mountainous areas, where the density of precipitation networks is far too low to capture their complex variability with topography. Adequate quantification of the remaining uncertainty in precipitation estimates is also necessary for further assimilation of complementary snow observations in snow models. Radar observations provide spatialised estimates of precipitation with high spatial and temporal resolution, and are often combined with rain gauge observations to improve the accuracy of the estimate. However, radar measurements suffer from significant shortcomings in mountainous areas (in particular, unrealistic spatial patterns due to ground clutter). Precipitation fields simulated by high-resolution numerical weather prediction (NWP) models provide an alternative estimate, but suffer from systematic biases and positioning errors. Even though these uncertainties can be partially described by ensemble NWP systems and systematic errors can be reduced by statistical post-processing, NWP precipitation estimates are still not reliable enough for the requirements of high resolution snow cover modelling. In this study, better precipitation estimates are obtained through a specific analysis based on a combination of all these available products. First, a pre-processing step is proposed to mitigate the main deficiencies of radar and gauges precipitation estimation products, focusing on reducing unrealistic spatial patterns. This method also provides a spatialised estimate of the associated error in mountainous areas, based on a climatological analysis of both radar and NWP-estimated precipitation. Three ensemble daily precipitation analysis methods are then proposed, first using only the modified precipitation estimates and associated errors, then combining them with ensemble NWP simulations based on the Particle Filter and Ensemble Kalman Filter data assimilation algorithms. The performance of the different precipitation analysis methods is evaluated at a local scale using independent ski resort precipitation observations. The evaluation of the pre-processing step shows its ability to remove the main spatial artefacts coming from the radar measurements and to improve the precipitation estimates at the local scale. The local scale evaluations of the ensemble analyses do not demonstrate an additional benefit of ensemble NWP forecasts, but their contrasted spatial patterns are challenging to evaluate with the available data. [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 National Weather Service (MeteoSwiss) uses two radar-based hail metrics: the probability of hail at 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 (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 Alarms 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 for a 60 % threshold, while for MESHS the highest CSI (0.25) and HSS (0.4) are obtained for 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]

Precipitation is arguably one of the most relevant surface variables impacting human lives on the planet, but global-coverage, high-resolution and good-quality observations are not readily available. In particular, gridded observational datasets are much needed for model development and forecast quality assessment. Here we compare the quality of four types of gridded precipitation products over Europe, namely: a rain-gauge interpolation; a satellite-derived product; a radar composite; and a reanalysis. Each product has its own strengths and weaknesses, and since each precipitation estimate uses different measuring techniques, we can employ a triangulation method to estimate the error variance of each product with respect to the unknown true values. Results show that: a) the satellite product has limited quality over Europe and may be problematic to use in quantitative forecast evaluation and diagnostics; b) the radar composite has spurious features that need to be considered when used in verification; c) all products struggle in topographically complex areas; d) the rain-gauge interpolation is not free of errors, despite rain gauges being often treated as ground truth in the literature; and e) the reanalysis dataset produces in some cases the best available estimates, in particular over the European near-coastal waters. [ABSTRACT FROM AUTHOR]

Downhill thunderstorms frequently occur in Beijing during the rainy seasons, leading to substantial precipitation. The accurate intensity prediction of these events remains a challenge, partly attributed to insufficient observational studies that unveil the thermodynamic and dynamic structures along the vertical direction. This study provides a comprehensive methodology for identifying both enhanced and dissipated downhill thunderstorms. In addition, a radar wind profiler (RWP) mesonet has been built in Beijing to characterize the pre-storm environment downstream to the thunderstorms at the mountain foot. This involves deriving vertical distributions of high-resolution horizontal divergence and vertical motion from the horizontal wind profiles measured by the RWP mesonet. A case of enhanced downhill thunderstorm on 28 September 2018 is carried out, which support the idea that the enhanced southerly flow and corresponding convergence detected by the RWP mesonet could favorably support the development of thunderstorms in the afternoon. The results also indicate that low-level convergence is an effective signal in accounting for convective maintenance. Statistical analysis based on radar reflectivity from April to September 2018–2021 have shown that a total of 63 thunderstorm events tend to be enhanced after moving into downhill and urban areas, accounting for about 66 % of the total number of downhill thunderstorm events. A critical region for intensified thunderstorms lies on the downslope side of the mountains west to Beijing. The evolution of the downhill storm is associated with the dynamic conditions over the plain compared to its initial morphology. The lifting induced by stronger westerly winds and vertical shear in the low and midlevel troposphere exerts a critical influence on the enhancement of convection. The findings underscore the significant role of RWP network in elucidating the evolution of downhill storm. [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]

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]

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]

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]

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]

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]

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 that is rescaled from reflectivity and incorporates an under and over estimate 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 datasets and applications. [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]

'Train effect' is defined as a cloud system in which several convective cells pass over the same place in a short time. Trains produce large amount of rainfall, frequently leading to flash floods, reported mainly over North America during spring and summer. Thirty train events were identified, using radar images, calibrated by rain-gauges, for four winters, all associated with Cyprus Lows (CL). The dynamic factors responsible for their formation in Israel were examined, utilizing the ECMWF Integrated Forecast System of 0.1° resolution. Seventeen out of the 30 events share common features. Each one was found at the cold sector in the southern periphery of a CL at its occluded stage, and located in the left flank of a maximum wind belt, where cyclonic shear vorticity exists. The trains cross the Israeli coast near 32.2°N, with a mean length of 45 km, last 1-3 hrs, and yield ~35 mm rainfall. The maximum wind belts right of the trains were found to delineate the limit of the precipitative region of the CLs. Unlike classical trains, activated by thermal or frontal forcing, the EM trains that develop in cold air-mass, can be entitled 'cold trains', rather than the classical 'warm trains'. [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]

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]

This work presents the importance of polarimetric variables as an additional data source for nowcasting thunderstorm hazards using an existing neural network architecture with recurrent-convolutional layers. The model can be trained to predict different target variables, which enables nowcasting of hail, lightning, and heavy rainfall for lead times up to 60 min with a 5 min resolution, in particular. The exceedance probabilities of Swiss thunderstorm warning thresholds are predicted. This study is based on observations from the Swiss operational radar network, which consists of five operational polarimetric C-band radars. The study area of the Alpine region is topographically complex and has a comparatively very high thunderstorm activity. Different model runs using combinations of single- and dual-polarimetric radar observations and radar quality indices are compared to the reference run using only single-polarimetric observations. Two case studies illustrate the performance difference when using all predictors compared to the reference model. The importance of the predictors is quantified by investigating the final training loss of the model, with skill scores such as critical success index (CSI), precision, recall, precision–recall area under the curve, and the Shapley value. Results indicate that single-polarization radar data are the most important data source. Adding polarimetric observations improves the model performance compared to reference model in term of the training loss for all three target variables. Adding quality indices does so, too. Including both polarimetric variables and quality indices at the same time improves the accuracy of nowcasting heavy precipitation and lightning, with the largest improvement found for heavy precipitation. No improvement could be achieved for nowcasting of the probability of hail in this way. [ABSTRACT FROM AUTHOR]

In this work, we study a Mediterranean cyclone, Helios, which took place during 9–11 February 2023 in the southeastern part of Sicily and Malta, by a multiparametric approach combining microseism results with sea state and meteorological data provided by wavemeter buoy, HF radar, hindcast maps and satellite SEVIRI images. The sub-tropical system Helios caused heavy rainfall, strong wind gusts and violent storm surges with significant wave heights greater than 5 m. We deal with the relationships between such a system and the features of microseism (the most continuous and ubiquitous seismic signal on Earth) in terms of spectral content, space–time variation of the amplitude and source locations tracked by means of two methods (amplitude-based grid search and array techniques). By comparing the location of the microseism sources and the area affected by significant storm surges derived from sea state data, we note that the microseism location results are in agreement with the real position of the storm surges. In addition, we are able to obtain the seismic signature of Helios using a method that exploits the coherence of continuous seismic noise. Hence, we show how an innovative monitoring system of the Mediterranean cyclones can be designed by integrating microseism information with other techniques routinely used to study meteorological phenomena. [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 snow melt 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 a 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]

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 and comprise information about past accumulation, ice deformation and allow to link 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 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 (LO-SAR) processing scheme to this setup. We showcase the system capability at Colle Gnifetti, Switzerland/Italy, and detect significantly deeper and older IRHs compared to previously deployed pulsed radar systems. Continuous IRHs are now apparent down to the base of the glacier. Corresponding reflection mechanisms for this glacier are linked to a stratified acidic impurity which was deposited at a higher rate due to increased industrial activity in the area. Possible improvements of 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]

Hail remains a major threat to agriculture in Switzerland and beyond and assessments of current and future hail risk are of paramount importance for decision-making in the insurance industry and the agricultural sector. However, relating observational information on hail with crop-specific damages is challenging. Here, we build and systematically assess a model to predict hail damage footprints for field crops (wheat, maize, barley, rapeseed) and grapevine from the operational radar product Maximum Expected Severe Hail Size (MESHS) at different spatial resolutions. To this end, we combine the radar information with detailed geospatial information on agricultural land use and geo-referenced damage data from a crop insurer for 12 recent hail events in Switzerland. We find that for field crops, model skill gradually increases when the spatial resolution is reduced from 1 km down to 8 km. For even lower resolutions, the skill is diminished again. On the contrary, for grapevine, a lower model resolution tends to reduce skill, which is attributed to the different spatial distribution of field crops and grapevine in the landscape. It is shown that identifying a suitable MESHS thresholds to model damage footprints always involves trade-offs. For the lowest possible MESHS threshold (20 mm) the model predicts damage about two times too often (high frequency bias and number of false alarms) but also has a high probability of detection (80 %). The frequency bias decreases for larger thresholds and reaches an optimal value close to 1 for MESHS thresholds of 30–40 mm. However, this comes at the cost of a substantially lower probability of detection (around 50 %) while overall model skill remains largely unchanged. We argue that, ultimately, the best threshold selection therefore depends on the user need and the costs of a false alarm or a missed event. Finally, the frequency of false alarms can be substantially reduced when only areas with high cropland density are considered. Results from this simple, open-source model show that modelling of hail damage footprints to crops from single-polarization radar in Switzerland is skillful and is best done at 8 km resolution for field crops and 1 km for grapevine. They further allow different users of such models to identify the suitable threshold for their application, taking into account associated trade-offs. [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]

This study assimilates for the first time polarimetric C-band radar observations from the German meteorological service (DWD) into DWD's convective-scale model ICON-D2 using DWD's ensemble-based KENDA assimilation framework. We compare the assimilation of conventional observations (CNV) with the additional assimilation of radar reflectivity Z (CNV + Z), with the additional assimilation of liquid or ice water content (LWC or IWC) estimates below or above the melting layer instead of Z where available (CNV + LWC /Z or CNV + IWC /Z respectively). Hourly quantitative precipitation forecasts (QPF) are evaluated for two stratiform and one convective rainfall events in the summers of 2017 and 2021. With optimized data assimilation settings (e.g., observation errors), the assimilation of LWC mostly improves first-guess QPF compared with the assimilation of Z alone (CNV + Z), whereas the assimilation of IWC does not, especially for convective cases, probably because of the lower quality of the IWC retrieval in these situations. Improvements are, however, notable for stratiform rainfall in 2021, for which the IWC estimator profits from better specific differential phase estimates owing to a higher radial radar resolution than the other cases. The assimilation of all radar data sets together (CNV + LWC + IWC + Z) yields the best first guesses. All assimilation configurations with radar information consistently improve deterministic 9 h QPF compared with the assimilation of only conventional data (CNV). Forecasts based on the assimilation of LWC and IWC retrievals on average slightly improve Fraction Skill Score (FSS) and Frequency Bias (FBI) compared with the assimilation of Z alone (CNV + Z), especially when LWC is assimilated for the 2017 convective case and when IWC is assimilated for the high-resolution 2021 stratiform case. However, IWC assimilation again degrades forecast FSS for the convective cases. Forecasts initiated using all radar data sets together (CNV + LWC + IWC + Z) yield the best FSS. The development of IWC retrievals that are more adequate for convection constitutes one next step to further improving the exploitation of ice microphysical retrievals for radar data assimilation. [ABSTRACT FROM AUTHOR]

Long term observations of deep convective cloud (DCC) vertical velocity and mass flux were collected during the GoAmazon2014/5 experiment. Precipitation echoes from a surveillance weather radar near Manaus, Brazil are tracked to identify and evaluate the isolated DCC lifecycle evolution during the dry and wet seasons. A Radar Wind Profiler (RWP) provides precipitation and air motion profiles to estimate the vertical velocity, mass flux, and mass transport rates within overpassing DCC cores as a function of the tracked cell lifecycle stage. The average radar reflectivity factor (Z), DCC area (A), and surface rainfall rate (R) increased with DCC lifetime as convective cells were developing, reached a peak as the cells matured, and decreased thereafter as cells dissipated. As the convective cells mature, cumulative DCC properties exhibit stronger updraft behaviors with higher upward mass flux and transport rates above the melting layer (compared to initial and later lifecycle stages). In comparison, developing DCCs have the lowest Z associated with weak updrafts, and negative mass flux and transport rates above the melting layer. Over the DCC lifetime, the height of the maximum downward mass flux decreased whereas the height of maximum net mass flux increased. During the dry season, the tracked DCCs had higher Z, propagation speed, and DCC area, and were more isolated spatially compared to the wet season. Dry season DCCs exhibit higher Z, mass flux, and mass transport rate while developing whereas wet season DCCs exhibit higher Z, mass flux, and mass transport rates at later stages. [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]

Coastal areas can be tremendously biodiverse and host a substantial part of the world's population and critical infrastructure. However, there are often fragile environments that face various hazards such as flooding, coastal erosion, land salinization or pollution, earthquake-induced land motion, or anthropogenic processes. In this article, we investigate the stability of the Nice Côte d'Azur Airport, which has been built on reclaimed land in the Var River delta (French Riviera, France). This infrastructure, as well as the ongoing subsidence of the airport runways, has been a permanent concern since the partial collapse of the platform in 1979. Here, we used the full archive of ESA SAR (synthetic-aperture radar) images from 1992 to 2020 to comprehensively monitor the dynamics of the airport subsidence. We found that the maximum downward motion rate has been slowing down from 16 mm yr -1 in the 1990s to 8 mm yr -1 today. However, sediment compaction is still active, and an acceleration phase of the continuous creep leading to a potential failure of a part of the platform cannot be excluded. Our study demonstrates the importance of remotely monitoring of the platform to better understand the motion of coastal land, which will ultimately help evaluate and reduce associated hazards. [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]

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]

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]

The Goddard Profiling algorithm (GPROF) converts radiometer observations aboard Global Precipitation Measurement (GPM) constellation satellites to precipitation estimates. Analyzing the accuracy of GPROF's estimates is vital to further improve the algorithm. Such analyses often use high-quality ground-based estimates as reference with a different spatial resolution. Often, the reference is resampled to match the satellite's resolution. However, the implemented sampling method to simulate the satellite's resolution varies amongst studies, which limits the transferability of conclusions. additionally, GPROF combines observations from various sensors and frequency channels, each with its own footprint size. Hence, uncertainties related to sampling are added on top of the uncertainty introduced when converting brightness temperatures to precipitation intensities. The contribution of sampling to the total amount of uncertainty remains unknown. Here, we quantify the uncertainty related to sampling while analyzing the current performance of GPROF over the Netherlands during a four year period (2017–2020). In this area, shallow and light precipitation frequently occur. Both precipitation types are often subject to research, as both types are difficult to detect with space-borne sensors. Only GPROF estimates based on observations from the conical-scanning radiometers of the GPM constellation are used. We investigate the uncertainty related to sampling by simulating the reference precipitation as satellite footprints that vary in size, geometry, and applied weighting technique. 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 (ETH) retrieved from the same radars are used to classify the precipitation as shallow, medium, or deep. The method used to spatially average the reference into a satellite footprint, i.e. using Gaussian weighting or the arithmetic mean, is found to exhibit a minimal influence on the retrieved estimate. The size of the sampled area is found to be the most influential. Still, the effect of using different footprint sizes cannot explain all the differences between the ground- and satellite based precipitation estimates. Additionally, the discrepancies between GPROF and the reference are largest for low ETH, while the relative bias between the different footprint sizes and implemented weighting methods increase with increasing ETH. Lastly, our results do not show a clear difference between coastal simulations and simulations over land. We conclude that the uncertainty introduced by merging different channels and sensors cannot fully explain the errors introduced by the retrieval algorithm. Hence, retrieval errors are found to be more prominent than sampling uncertainties, in particular for shallow and light precipitation. [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]

In this research, we delve into the influence of cloud condensation nuclei (CCN) and ice-nucleating particles (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 INP 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 are the observation that increased INP concentrations yield smaller ice particles, while a surge in CCN concentrations leads to a subtle growth in their dimensions. Consistent with existing literature, our results spotlight plate-like crystals as dominant between temperatures of −20 to −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 INP 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]

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 (SOA2RSE) 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. In addition, we 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]

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 km2–500 km2) 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 km2 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 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]

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]