Ice Water Content Research Articles

Technique for comparison of backscatter coefficients derived from in situ cloud probe measurements with concurrent airborne lidar

Abstract. Jet engine power loss due to ice particle accumulation is a recognized aviation hazard occurring in cloud conditions difficult to forecast or visually recognize. High-altitude cirrus clouds can have ice particle concentrations high enough to be dangerous; therefore, pilots must be informed when aircraft enter such environments. One approach to determining ice particle concentration is an onboard lidar system. Concurrent lidar measurements are compared to backscatter coefficients derived from particle size distributions obtained from wing-mounted, in situ probes during four case studies consisting of sixty-second flight segments at different temperatures: +7 and +4 ∘C for water droplet analysis, and −33 and −46 ∘C for ice particle analysis. Backscatter coefficients derived from external cloud probes (ECP) are correlated (0.91) with measurements by an airborne lidar system known as the Optical Ice Detector (OID). Differences between OID and ECP backscatter coefficients range from less than 1 to over 3 standard deviations in terms of uncertainties. The backscatter coefficients are mostly in agreement for liquid clouds and are in disagreement for the −33 and −46 ∘C cases, with ECP-derived backscatter coefficients lower than the OID for three out of the four cases. Measurements over four 60 s research flight segments show that measured total water content is correlated (0.74) with the OID backscatter coefficient, which indicates that the OID is a useful instrument for determining ice particle concentrations over a broad range of environments, including at ice water contents as low as 0.02 g m−3. Additionally, concurrent measurements from cloud imaging probes and the OID provide improved knowledge of cloud conditions, which may help in understanding cloud processes.

Estimating Microphysics Properties in Ice-Dominated Clouds from Airborne Ka–W-band Dual-Wavelength Ratio Reflectivity Factor in Close Proximity to In Situ Probes

Abstract Properties of frozen hydrometeors in clouds remain difficult to sense remotely. Estimates of number concentration, distribution shape, ice particle density, and ice water content are essential for connecting cloud processes to surface precipitation. Progress has been made with dual-frequency radars, but validation has been difficult because of lack of particle imaging and sizing observations collocated with the radar measurements. Here, data are used from two airborne profiling (up and down) radars, the W-band Wyoming Cloud Radar and the Ka-band Profiling Radar, allowing for Ka–W-band dual-wavelength ratio (DWR) profiles. The aircraft (the University of Wyoming King Air) also carried a suite of in situ cloud and precipitation probes. This arrangement is optimal for relating the “flight-level” DWR (an average from radar gates below and above flight level) to ice particle size distributions measured by in situ optical array probes, as well as bulk properties such as minimum snow particle density and ice water content. This comparison reveals a strong relationship between DWR and the ice particle median-volume diameter. An optimal range of DWR values ensures the highest retrieval confidence, bounded by the radars’ relative calibration and DWR saturation, found here to be about 2.5–7.5 dB. The DWR-defined size distribution shape is used with a Mie scattering model and an experimental mass–diameter relationship to test retrievals of ice particle concentration and ice water content. Comparison with flight-level cloud-probe data indicate good performance, allowing microphysical interpretations for the rest of the vertical radar transects.

Responses of CIPS/AIM noctilucent clouds to the interplanetary magnetic field

Abstract. This study investigates the link between the interplanetary magnetic field (IMF) By component and the noctilucent clouds (NLCs) measured by the Cloud Imaging and Particle Size (CIPS) experiment onboard the Aeronomy of Ice in the Mesosphere (AIM) satellite. The mean ice particle radius in NLCs is found to be positively correlated with IMF By in the Southern Hemisphere (SH) and negatively correlated with IMF By in the Northern Hemisphere (NH), respectively, on a day-to-day timescale in most of the 20 summer seasons during the 2007–2017 period with a near 0 d lag time, and the response in the SH is stronger than that in the NH. Moreover, the albedo, ice water content and frequency of occurrence of NLCs present positive correlation with IMF By in the SH but no significant correlation in the NH. The superposed epoch analysis (SEA) further indicates the rm on average changes by about 0.73 nm after IMF By reversals, which is significant at the 90 % confidence level in Monte Carlo sensitivity tests. Our results suggest an IMF By-driven pathway: the influence of the solar wind on the polar ionospheric electric potential affects the nucleation processes in NLCs and consequently the ice particle radius and NLC brightness.

Characteristics of two different types of thunderstorms in summer over the Nagqu area in China

To investigate further the characteristics of thunderstorms on the Qinghai–Tibet plateau, a plateau vortex (PV) thunderstorm and a local thermoconvective (TC) thunderstorm over the Nagqu area are analyzed using cloud-radar, microwave-radiometer, and raindrop-spectrometer data, and their macrophysical and microphysical evolution characteristics are discussed in terms of thermodynamic processes, microphysical processes, and lightning activities. The results show the following. 1) The cloud layer was deeper in the PV thunderstorm, but the TC thunderstorm had a stronger short-time updraft with a radial velocity exceeding 10 m/s, and the warming from the strong updraft action and latent heat release from the hydromorphic phase change was about twice that of the PV thunderstorm. 2) The water vapor density increased significantly when the thunderstorm cloud passed, and the liquid water content in the middle and lower layers exceeded 4 g/m3. The maximum ice water content in the TC thunderstorm was twice bigger than that in the PV thunderstorm. The trends of raindrop number concentration and rain intensity of the PV thunderstorm were similar, and the average particle size of raindrops was smaller than that of the TC thunderstorm. 3) Both types of thunderstorms accounted for more than 90% of negative cloud-to-ground (CG) lightning, and the regions with black-body temperature (TBB) less than −40 °C and a larger TBB gradient were more favorable for the occurrence of CG lightning. 4) Thermodynamic fields bring water vapor and lift for microphysical processes, and microphysical changes release latent heat to enhance the dynamic effects, which together promote the development of lightning activities. The peak radial velocity and ice-phase particle concentration were more than 10 min ahead of the active peak of the CG lightning. This study reveals the macrophysical and microphysical evolution characteristics of different types of thunderstorms and provides a certain scientific basis for disaster prevention and mitigation regarding thunderstorms over the Nagqu area.

Contrasts in the Evolution and Microphysical Features of Two Convective Systems during a Heavy Rainfall Event along the Coast of South China

On 1 June 2021, a heavy rainstorm hit the coast of South China (148.6 mm in 1 h, 361 mm over 12 h). The storm process was successively affected by two convective systems (CSs). The initial convection of the two CSs occurred at a similar location; however, they subsequently showed different evolution characteristics. Based on multi-source data, including dual-polarimetric radars, wind profiling radars, sounding, and automatic weather stations, we explored the differences in the key characteristics of these two CSs. It was found that the convection was initially triggered at a similar location for both CSs, closely related to the mesoscale boundary and the hilly terrain. After formation, CS1 moved eastward to the regions with lower surface temperature and weaker lower-level convergence but similar humidity, which means the environmental conditions for sustaining the CS became less favorable. As a result, CS1 dissipated rapidly and only lasted for about 90 min, resulting in 5% of the total precipitation of the overall storm. In contrast, during the lifespan of CS2, the southerly wind over the South China Sea became stronger. This caused an intense lower-level convergence zone along the coastal region of Guangdong Province, which provided favorable dynamic conditions for maintaining CS2. Favored by the strong coastal convergence and abundant moisture, new convective cells (CCs) were generated continuously and merged with CS2, acting as another favorable condition for its sustainment. Overall, CS2 lasted for 8 h, and its precipitation accounted for 95% of the total rainfall. In CS1, CCs showed a notable evaporation process below 4 km, manifested by the large raindrops. However, in CS2, the CCs had a higher concentration of small raindrops and higher ice and liquid water content. Since CS2 was close to the coastal region, the warm local environment promoted convection, leading to intense precipitation. In addition, the riming and melting processes were active, leading to a high precipitation efficiency and strong local precipitation during a short period of time.

Variation and Influencing Factors of Cloud Characteristics over Qinghai Lake from 2006 to 2019

Clouds are an indispensable part of climate change, and the occurrence and development of clouds in the Qinghai Lake area (QHL) have great significance for the regional energy budget and precipitation system. To a certain extent, clouds will affect the water resources, agriculture, animal husbandry, and photovoltaic power industry in this region. In this study, we used CloudSat satellite data, combined with meteorological elements and atmospheric circulation, to analyze the cloud occurrence frequency and cloud water content in QHL. The results demonstrate that the frequency of cloud occurrences in QHL is 33% with a decreasing trend from 2006 to 2019. Altostratus and Nimbostratus are the main types of cloud systems in QHL. The cloud ice water content is 62.21 mg/m3 and the cloud liquid water content is 261.66 mg/m3. The highest value of the vertical cloud fraction occurs from March to June, at a height of 7–11 km in QHL. The height of the mixed-phase clouds is approximately 4–8 km and the ice clouds are above 8 km. The vertical distribution of ice particles is relatively dispersed, while the vertical distribution of liquid particles is relatively concentrated. The time and height of high particle effective radius and high particle concentration are consistent with the high value of cloud water content. The decrease in total cloud occurrence frequency in QHL is caused by the increase in temperature. This study helps to clarify the detailed structure of clouds and the distribution of cloud water resources, which has an important reference value for the study of climate change impact and the sustainable development of lake resources in QHL.

The impacts of secondary ice production on microphysics and dynamics in tropical convection

Abstract. Secondary ice production (SIP) is an important physical phenomenon that results in an increase in the ice particle concentration and can therefore have a significant impact on the evolution of clouds. In this study, idealized simulations of a mesoscale convective system (MCS) were conducted using a high-resolution (250 m horizontal grid spacing) mesoscale model and a detailed bulk microphysics scheme in order to examine the impacts of SIP on the microphysics and dynamics of a simulated tropical MCS. The simulations were compared to airborne in situ and remote sensing observations collected during the “High Altitude Ice Crystals – High Ice Water Content” (HAIC-HIWC) field campaign in 2015. It was found that the observed high ice number concentration can only be simulated by models that include SIP processes. The inclusion of SIP processes in the microphysics scheme is crucial for the production and maintenance of the high ice water content observed in tropical convection. It was shown that SIP can enhance the strength of the existing convective updrafts and result in the initiation of new updrafts above the melting layer. Agreement between the simulations and observations highlights the impacts of SIP on the maintenance of tropical MCSs in nature and the importance of including SIP parameterizations in models.

A Modified SWAT Model to Simulate Soil Water Content and Soil Temperature in Cold Regions: A Case Study of the South Saskatchewan River Basin in Canada

Soil water content (SWC) and soil temperature are important hydrologic state variables. Accurate model simulation is critical in hydrologic regimes in cold regions dominated by spring snowmelt. In this study, we developed a combined physically-based soil temperature and energy-balance rain-on-snow (ROS) module for the Soil and Water Assessment Tool (SWAT) model and applied it to the South Saskatchewan River Basin (SSRB). We calibrated the SWAT base (SWAT-B) model and the SWAT modified (SWAT–M) model using daily measured soil temperature and SWC by hydrological response unit (HRU) for the years 2015 to 2020. The results of sensitivity analysis using the SUFI-2 technique in SWAT-CUP indicated that eight parameters have the most significant (p < 0.5) effect on streamflow, soil moisture, and snowmelt. Statistics for the SWAT-B and SWAT-M streamflow models revealed that the new module improved the Nash-Sutcliffe efficiency (NSE) from 0.39 to 0.71 and 0.42 to 0.76 for calibration and validation, respectively. The statistics for SWAT-simulated daily SWC showed that the measured data were a better fit with SWAT-M versus the SWAT-B output. Furthermore, SWAT-B values exceeded SWAT-M output and field measurements, and thus, the range of SWAT-M results was a better fit with observations. SWAT-B tended to underestimate soil temperature in the cold season, while SWAT-M significantly improved soil temperature simulation for winter. This new SWAT module simulated freeze-thaw cycles and captured the influence of snow cover on surface soil ice-water content. Spatial analysis of SWC and soil temperature across the SSRB showed that the SWAT-M model predicted more SWC and lower soil temperature in the western part of SSRB than SWAT-B, with higher soil temperature and lower SWC in the eastern region.

Cloud Macro- and Microphysical Properties in Extreme Rainfall Induced by Landfalling Typhoons over China

Extreme rainfall induced by landfalling typhoon (ERLTC) can cause destructive natural disasters throughout China. Cloud properties in ERLTC are not yet well understood and parameterized, which limits the forecast accuracy of ERLTC to some extent. The 99th percentile intensity of daily rainfall associated with LTC is objectively defined as ERLTC and using the CloudSat tropical cyclone (CSTC) dataset from 2006 to 2018, cloud macro- and microphysical characteristics are statistically investigated. Results show that the proportion of single-layer (double-layered) clouds increases (decreases) significantly on the occurrence day of ERLTC. In the TC inner core region, the proportion of deep convective cloud at 2–10 km is the highest, reaching 50%. In the TC envelop region, deep convective cloud at the height of 3–8 km and cirrus at the height of 12–14 km account for the highest proportions. For the TC outer region, cirrus around 13 km has the highest proportion. During the ERLTC period, the ice-water content is mainly distributed in 5–18 km, and is mostly distributed in the TC inner core, followed by the envelop region. A large number of smaller ice particles are gathering in the upper troposphere at 13–18 km, while a small number of larger ones is gathering in the middle levels around 8–10 km. These results are useful for evaluating the ERLTC simulations and are expected to provide new forecasting factors for ERLTC in cloud macro- and microphysical perspectives.

Spatial and Temporal Distribution of Cloud Water in the Yellow River Basin, China

The Yellow River Basin is essential to China’s economic and social development and ecological security. In order to assess the temporal and spatial distribution of cloud water in the Yellow River Basin, we analyzed the distribution characteristics of water vapor content and cloud water content using ERA5 monthly averaged data from 1980 to 2019. The results showed that the high-value area of the annual average atmospheric water vapor content distribution was concentrated above the North China Plain in the eastern part of the basin, and the value was mostly between 21 and 24 mm. The low-value areas were mainly centered above the high mountain areas in the western part of the basin, and the value mostly fell between 3 and 6 mm. The seasonal distribution characteristics of the annual average water vapor content were relatively consistent with the annual average distribution characteristics. The high-value cloud ice water content area was in the northeastern part of the Qinghai-Tibet Plateau (4.5 g·m−2), while the low-value area appeared on the Loess Plateau (2 g·m−2). The high-value area of cloud liquid water content was on the north side of the West Qinling Mountains (12 g·m−2). The low-value area appeared on the Loess Plateau and the northern part of the Qinghai-Tibet Plateau (3 g·m−2). The cloud water content was higher in the eastern region than in the western region in the overall spatial distribution, and the content of cloud liquid water was higher than that of cloud ice water. The average annual atmospheric water vapor content was increasing, and the annual average content of cloud ice water and cloud liquid water was declining. The change in the total amount and spatial distribution of cloud water was not obvious in the Yellow River Basin.

Arctic Climate Feedbacks in ERA5 Reanalysis: Seasonal and Spatial Variations and the Impact of Sea‐Ice Loss

AbstractRadiative climate feedbacks in the Arctic have been extensively studied, but their spatial and seasonal variations have not been thoroughly examined. Using ERA5 reanalysis data, we examine seasonal variations in Arctic climate feedbacks and their relationship to sea‐ice loss based on changes from 1950–1979 to 1990–2019. The spring and summer seasons experienced large sea‐ice loss, strong surface albedo feedback, and large oceanic heat uptake. Arctic clouds exerted small net cooling in May‐June‐July but moderate warming during the cold season, especially over areas with large sea‐ice loss where cloud liquid and ice water content increased. Arctic water vapor feedback peaked in summer but was weak and uncorrelated with sea‐ice loss. Arctic positive lapse rate feedback (LRF) was strongest in winter over areas with large sea‐ice loss and weak inversion but uncorrelated with atmospheric stability, suggesting that oceanic heating from sea‐ice loss led to enhanced surface warming and the positive LRF.

Retrieving ice-nucleating particle concentration and ice multiplication factors using active remote sensing validated by in situ observations

Abstract. Understanding the evolution of the ice phase within mixed-phase clouds (MPCs) is necessary to reduce uncertainties related to the cloud radiative feedback in climate projections and precipitation initiation. Both primary ice formation via ice-nucleating particles (INPs) and secondary ice production (SIP) within MPCs are unconstrained, not least because of the lack of atmospheric observations. In the past decades, advanced remote sensing methods have emerged which provide high-resolution data of aerosol and cloud properties and could be key in understanding microphysical processes on a global scale. In this study, we retrieved INP concentrations and ice multiplication factors (IMFs) in wintertime orographic clouds using active remote sensing and in situ observations obtained during the RACLETS campaign in the Swiss Alps. INP concentrations in air masses dominated by Saharan dust and continental aerosol were retrieved from a polarization Raman lidar and validated with aerosol and INP in situ observations on a mountaintop. A calibration factor of 0.0204 for the global INP parameterization by DeMott et al. (2010) is derived by comparing in situ aerosol and INP measurements, improving the INP concentration retrieval for continental aerosols. Based on combined lidar and radar measurements, the ice crystal number concentration and ice water content were retrieved and validated with balloon-borne in situ observations, which agreed with the balloon-borne in situ observations within an order of magnitude. For seven cloud cases the ice multiplication factors (IMFs), defined as the quotient of the ice crystal number concentration to the INP concentration, were calculated. The median IMF was around 80, and SIP was active (defined as IMFs > 1) nearly 85 % of the time. SIP was found to be active at all observed temperatures (−30 to −5 ∘C), with the highest IMFs between −20 and −5 ∘C. The introduced methodology could be extended to larger datasets to better understand the impact of SIP not only over the Alps but also at other locations and for other cloud types.

Contrail formation within cirrus: ICON-LEM simulations of the impact of cirrus cloud properties on contrail formation

Abstract. Contrail formation within natural cirrus introduces large perturbations in cirrus ice crystal number concentrations, leading to modifications in cirrus microphysical and optical properties. The number of contrail ice crystals formed in an aircraft plume depends on the atmospheric state as well as aircraft and fuel properties. Our aim is to study contrail formation within cirrus and, in particular, the impact of pre-existing cirrus on the contrail formation processes. We implement a parameterization for contrail ice nucleation and the survival of contrail ice crystals within the vortex phase within the high-resolution ICON-LEM at a horizontal resolution of 625 m over Germany and modify the parameterizations in order to consider the impact of the pre-existing cirrus. We then simulate contrail formation and analyze the change in ice nucleation and survival due to the presence of cirrus ice crystals. We have selected two different synoptic situations to sample a large range of cirrus cloud properties from those representative of average thick in situ formed cirrus and liquid origin cirrus that are connected to a frontal system down to very thin cirrus within a high-pressure system. We find that contrail formation within cirrus mostly leads to increases in cirrus ice crystal numbers by a few orders of magnitude. Pre-existing cirrus has a discernible impact on the contrail formation threshold temperature and on contrail ice crystal number concentrations only if the cirrus ice water content is high and ice supersaturation is low. The cirrus ice crystals sucked into and subsequently sublimated within the aircraft's engine and the ice crystals mixed into the aircraft plume experiencing sublimation and later deposition lead to an increase in the contrail formation threshold temperature most of the time. Often this increase is negligible, but in the above atmospheric conditions the increase in the contrail formation threshold can be large, reaching maximum values of 2 K. Consequently, significant increases in ice nucleation are very seldom, but in areas of high cirrus ice water content and ice crystal number concentrations, contrail ice nucleation rates can be significantly increased, in particular at lower flight levels. Cirrus ice crystals can also lead to a reduction in ice nucleation, but negative changes are significantly smaller and less common. After nucleation cirrus and contrail ice crystals compete for water vapor, while the combined contrail and cirrus ice water mass grows. Once the aircraft plume gets trapped within the wake vortices and descends the plume relative humidity decreases so that eventually both cirrus and contrail ice crystals sublimate. We find that the impact of cirrus ice crystals on the contrail ice crystal loss within the descending vortices is negligible. Only for nucleation rates several orders of magnitude lower than the soot number emission index can the impact be noticeable, but the relevance of this impact is limited.

Novel Parameterization of Ice Cloud Effective Diameter from Collocated CALIOP-IIR and CloudSat Retrievals

Abstract Satellite-based measurements of global ice cloud microphysical properties are sampled to develop a novel set of physical parameterizations, relating to cloud layer temperature and effective diameter De, that can be implemented for two separate applications: in numerical weather prediction models and lidar-based cloud radiative forcing studies. Ice cloud optical properties (i.e., spectral scattering and absorption) are estimated based on the effective size and habit mixture of the cloud particles. Historically, the ice cloud De has been parameterized from aircraft in situ measurements. However, aircraft-based parameterizations are opportunistic in that they only represent specific types of clouds (e.g., convective anvil, tropopause-topped cirrus) in the regions in which they were sampled and, in some cases, are limited in fully resolving the entire vertical cloud layer. Breaking away from the aircraft-based parameterization paradigm, this study is the first of its kind to attempt a parameterization of De as a function of temperature, ice water content (IWC), and lidar-derived extinction from satellite-based global oceanic measurements of ice clouds. Data from both active and passive remote sensing sensors from two of NASA’s A-Train satellites, CloudSat and CALIPSO, are collected to guide development of globally robust parameterizations of all ice cloud types and one exclusively for cirrus clouds. Significance Statement We derived unique parameterizations of ice crystal effective size from global satellite measurements in an effort to more robustly and consistently represent ice clouds in numerical models for weather forecasting and climate energy balance studies. Based on our results, effective ice crystal size is easily solved based on temperature and visible cloud translucence. By knowing the size of the ice crystals, we can then estimate cloud scattering and absorption. In comparison with aircraft-based parameterizations, the satellite data reveal that ice crystal effective sizes are much smaller, on global average, for ice clouds occurring in relatively warm layers (>230 K), indicating that many ice clouds are more reflective than previously believed.

Analysing spatio-temporal patterns in wintertime rainfall across Iran’s deserts using GPM DPR data

Although known as the warmest deserts of the world, rainfall and integrated water of Lut and Kavir Deserts is still unknown due to insufficient weather stations. The Dual-Frequency Precipitation Radar (DPR) onboard the Global Precipitation Measurement (GPM) mission satellite created an opportunity for use to study the rainfall and vertically integrated liquid water content (LWC) and integrated non-liquid (ice) water content (IWC) using statistical and distance analyses over the two deserts during winter months (December to March) of 2015–2020. The results showed good similarity between the GPM DPR and station data with a Pearson correlation coefficient of 0.81. March has the highest contribution in rain amount, which is about 37% in both deserts. March also has the largest contribution of the rainy days, being about 41.4 and 37.6% in the Lut and Kavir Deserts, respectively. LWC has the highest amount in the Lut Desert, while IWC is the largest in the Kavir Desert. Distance analysis showed that there is a significant increasing trend of rainfall from west to east in the Lut Desert. Elevation does not affect rainfall distribution strongly, but rainfall is highly influenced by the atmospheric-driven large-scale parameters.

Diurnal Variation in Cloud and Precipitation Characteristics in Summer over the Tibetan Plateau and Sichuan Basin

The diurnal variation in precipitation and cloud parameters and their influencing factors during summer over the Tibetan Plateau (TP) and Sichuan Basin (SB) were investigated using the Hydro-Estimator satellite rainfall estimates, ground observations, and ERA5 dataset. The precipitation and cloud parameters show diurnal propagation over the SB during the mei-yu period in contrast to such parameters over the TP. The diurnal maximum precipitation from the Hydro-Estimator satellite and cloud ice and liquid water content (cloud LWC and IWC) from the ERA5 dataset are concentrated in the early evening, while their diurnal minimums manifest in the morning. Cloud LWC accounts for more than 60% of the total water during almost the entire diurnal cycle over the inner TP and SB during the mei-yu period. The IWC accounts for more than 60% of the total water in the late afternoon over the edge of the SB and TP. The cloud base height (CBH) above ground level (AGL), the lifting condensation level (LCL) AGL, and the zero degree level AGL are almost equal over the TP during the summer period. The zero degree level AGL over the SB is higher than that over the TP because the air temperature lapse rate over the TP is larger. The thickness of liquid water cloud over the SB is larger than that over the TP. The correlation analysis shows that the CBH AGL and LCL AGL over the TP are related to the dewpoint spread, but less so over the SB because of the stronger turbulence and lower air density over the TP than the SB. Convective available potential energy has a larger impact on precipitation over the TP than the SB. The cloud LWC makes a larger contribution to the precipitation over the SB than over the TP, which is related to the mean zonal wind and diurnal cycle of low-level winds. The precipitation at the edge of the TP and SB (i.e., the steep downstream slope) is largely influenced by the ice water contained within clouds owing to the convergence rising motion over the slopes.

Dual-Frequency Radar Retrievals of Snowfall Using Random Forest

The microphysical parameters of snowfall directly impact hydrological and atmospheric models. During the International Collaborative Experiment hosted at the Pyeongchang 2018 Olympic and Paralympic Winter Games (ICE-POP 2018), dual-frequency radar retrievals of particle size distribution (PSD) parameters were produced and assessed over complex terrain. The NASA Dual-frequency Dual-polarized Doppler Radar (D3R) and a collection of second-generation Particle Size and Velocity (PARSIVEL2) disdrometer observations were used to develop retrievals. The conventional look-up table method (LUT) and random forest method (RF) were applied to the disdrometer data to develop retrievals for the volume-weighted mean diameter (Dm), the shape factor (mu), the normalized intercept parameter (Nw), the ice water content (IWC), and the snowfall rate (S). Evaluations were performed between the D3R radar and disdrometer observations using these two methods. The mean errors of the retrievals based on the RF method were small compared with those of the LUT method. The results indicate that the RF method is a promising way of retrieving microphysical parameters, because this method does not require any assumptions about the PSD of snowfall.

High Ice Water Content Conditions Associated with Wintertime Elevated Convection in the Midwest

AbstractAircraft flying through areas of high ice water content (HIWC) can experience engine damage and/or failure. HIWC is typically associated with convection and the microphysical properties of tropical oceanic and coastal convection are well documented as a result of several field campaigns in the past decade. HIWC appears to be less common in extratropical convection, but instances of HIWC-related aircraft issues have been recorded in extratropical weather, even during winter. The present study documents the microphysical properties of HIWC between −25° and 0°C and the meteorological and thermodynamic conditions around that HIWC from five flights from the In-Cloud Icing and Large-Drop Experiment (ICICLE) in the midwestern United States in February 2019. All five cases contained elevated convection above a strong low-level temperature inversion. Values for top-of-inversion mixing ratios and total column precipitable water were about 5 g kg−1and 20 mm, respectively, according to soundings near each case. A maximum ice water content of 2.1 g m−3was observed over a length scale of about 500 m, and ice particle size distributions had mean volume equivalent diameters around 1000μm. Supercooled drizzle droplets were also observed in the vicinity of the HIWC, raising questions about the possible role of secondary ice production via the freezing and shattering of supercooled large droplets in HIWC formation. The generalizability of these results is limited by the small number of cases, but they provide some of the first in situ observations of extratropical winter HIWC and highlight the need for continued research on these conditions.Significance StatementHigh ice water content (HIWC) conditions can cause engine damage, stall, and failure in aircraft and can cause air data probes to report erroneous values. Most research on HIWC has focused on tropical convection. This paper serves to draw more attention to the aviation hazard posed by extratropical winter HIWC and the need for additional research into these environments. The cases examined also contained supercooled precipitation-sized liquid droplets. Combined with other recent research, these observations may help to motivate laboratory experiments investigating the role of secondary ice production via droplet freezing and shattering in the formation of HIWC.

Seismic physics-based characterization of permafrost sites using surface waves

Abstract. The adverse effects of climate warming on the built environment in (sub-)arctic regions are unprecedented and accelerating. The planning and design of climate-resilient northern infrastructure, as well as predicting deterioration of permafrost from climate model simulations, require characterizing permafrost sites accurately and efficiently. Here, we propose a novel algorithm for the analysis of surface waves to quantitatively estimate the physical and mechanical properties of a permafrost site. We show the existence of two types of Rayleigh waves (R1 and R2; R1 travels faster than R2). The R2 wave velocity is highly sensitive to the physical properties (e.g., unfrozen water content, ice content, and porosity) of active and frozen permafrost layers, while it is less sensitive to their mechanical properties (e.g., shear modulus and bulk modulus). The R1 wave velocity, on the other hand, depends strongly on the soil type and mechanical properties of permafrost or soil layers. In situ surface wave measurements revealed the experimental dispersion relations of both types of Rayleigh waves from which relevant properties of a permafrost site can be derived by means of our proposed hybrid inverse and multiphase poromechanical approach. Our study demonstrates the potential of surface wave techniques coupled with our proposed data-processing algorithm to characterize a permafrost site more accurately. Our proposed technique can be used in early detection and warning systems to monitor infrastructure impacted by permafrost-related geohazards and to detect the presence of layers vulnerable to permafrost carbon feedback and emission of greenhouse gases into the atmosphere.

ENTICE Satellite Orbital Simulator to Study Ice Clouds

AbstractClouds play a significant role in the Earth's energy balance and hydrologic cycle through their effects on radiation and precipitation, and therefore are crucial for life on Earth. Earth's NexT‐generation ICE mission (ENTICE) is proposed to measure diurnally resolved global vertical profiles of cloud ice particle size, ice water content, and in‐cloud humidity and temperature combining a 94 GHz cloud radar and multi‐frequency sub‐millimeter (sub‐mm) microwave radiometers from space. The scientific objective of ENTICE is to identify the important processes by which anvil clouds evolve and interact with ambient thermodynamic conditions to advance our fundamental understanding of clouds and reduce uncertainties in cloud climate feedback. Whether such a science objective could be achieved depends on the orbital sampling characteristics of the mission. In this study, ENTICE sampling statistics are simulated with more than 1 billion (109) profiles and using five different scanning methods in a 400 km altitude precession orbit with an inclination of 65°: nadir, forward pointing, side scanning, and conical scanning for the radiometers, and nadir pointing for the radar. These statistics can then be used to help with future mission planning efforts. Using the GEOS‐5 model atmosphere (produced by nature‐run at 7‐km and 30‐min resolution), the simulator calculates sampling statistics related to cloud types and overpass times with enhancement from radar. The wide swath of ENTICE radiometers by conical and side scanning methods ensures ample high cloud samples (1.3 × 109 samples) gathered by ENTICE over its 2‐year mission for different types of clouds with sufficient sampling over the diurnal cycles. Sampling differences between radar and radiometers at nadir shows the combination of radar and radiometers will allow for measurements of cloud vertical profiles. Therefore, our results show that the designed orbit sampling of ENTICE is sufficient to fulfill the mission science goals.