Cloud Top Height Research Articles

Identifying Hail Events Using Weather Radar and Satellite Data (Study Case in East Java on 23 December 2021)

Hail is an extreme weather phenomenon that rarely occurs in tropical regions such as Indonesia. However, if convective clouds are strong enough to develop higher, a hail event is likely to occur. Interestingly, on 23 December 2021, hail events occurred in four different nearby locations. Weather Radar data from the Juanda Meteorological Station and Himawari-8 Satellite data were used to analyse this hail phenomenon. The Radar measured reflectivity values, cloud top height, water vapour content, and the derivative product of the probability of hail occurrence. The Himawari-8 data were utilised to observe cloud top temperatures. From the results, it is known that the reflectivity value of hail is in the range of 45 dBZ to 55 dBZ, with a probability of hail above 80%, and the temperature at the top of the Cumulonimbus cloud is between -67.5°C to -82.5°C. These thresholds may be used in monitoring hail events, especially in the East Java region.

The effects of gravity waves on ozone over the Tibetan Plateau

The ozone valley over the Tibetan Plateau (TP) has an important effect on global weather and climate. As a significant source of atmospheric gravity waves (GWs), TP is also a key area of global stratosphere-troposphere exchange (STE), yet the exploration of the causes of low ozone values over the TP has been scarce from the perspective of GW. In this paper, we used the hyperspectral Atmospheric Infrared Sounder (AIRS) data to extract GWs and deep convection over the TP and analyzed the statistical characteristics of ozone for corresponding events. The results show that: (1) The GWs observed by AIRS mainly occur in the mid-latitude region, with a maximum frequency of >10%, and their intensity varies with longitude. Three large-value regions have been identified on the coast of Southeast Asia, Central America, and Central Africa. (2) Most deep convection occurs in the Intertropical Convergence Zone (ITCZ), with three large value centers in the East Asian monsoon region, Central America and Central Africa. Among them, the Indian monsoon region, the Northwest Pacific region to the east of the Philippines, and the western side of the South China Sea (SCS) are the large-value areas within the East Asian monsoon region, with the maximum occurrence frequency exceeding 15%. (3) TP is a region with a large frequency of convective gravity waves. In most areas of the southern foothills of the TP, the frequency of GWs related to convection is 30%–50%. Almost all deep convection over the TP and its adjacent areas are accompanied by GWs, with a frequency of 60–90%. (4) The composite analysis results show that the occurrence of GWs or deep convection is accompanied by negative anomalies of ozone in most areas of the southern foothills of the TP, indicating a decrease in the integrated column amount of ozone. (5) By analyzing the anomaly distribution of meteorological elements such as cloud fraction, cloud top height, tropopause height, and outgoing longwave radiation (OLR), the results show a dual effect in the presence of GWs or deep convection. The effect involves not only an increased generation of clouds, conducive to convection but also elevated cloud top height, promoting the formation of high clouds, some even beyond the tropopause. In addition, the occurrence of GW or deep convection can also affect the tropopause height, strengthen convective activity, and facilitate the upward transfer of low-concentration ozone from the lower layer to the upper layer, thereby causing a decrease in the integrated column amount of ozone.

Helios and Juliette: Two falsely acclaimed medicanes?

The present work analyzes the synoptic, thermodynamic, and microphysics characteristics of two Mediterranean cyclones that occurred in February–March 2023. The analysis is mainly carried out through the use of passive microwave (PMW) satellite measurements, which allow us to follow the cyclones' evolution and state whether Helios and Juliette can be considered as Mediterranean tropical-like cyclones (i.e., medicanes). Both cyclones show a very similar evolution, with a low-stratospheric warm air anomaly during the development phase, followed by the development of a warm anomaly in the low−/mid-troposphere. This feature is often observed in medicanes (e.g., Qendresa, Zorbas), except for few cases (i.e., Medicane Ianos, which shows a warm core (WC) development clearly driven by diabatic processes without a preliminary warming signal in the lower stratosphere and upper troposphere). The analysis carried out highlights that, while Helios maintains this setting through its whole lifetime, Juliette undergoes tropical transition in the final stage of its evolution. As opposed to most medicane cases, the PMW precipitation microphysics diagnostics shows the predominance of shallow clouds, with almost total absence of ice hydrometeors and deep convection in the proximity of the WC center (i.e., within 100 km radius) for Helios and during the initial stage of Juliette. PMW radiometry provides strong indication that diabatic heating plays a role in the WC development when the onset of deep convection features is identified in the proximity of the Juliette cyclone center. Moreover, the PMW cloud-top height product does not show a closed cloud-free eye for Helios, while it is observed for the final stage of Juliette as often happens during medicanes' mature phase. Therefore, we deem that while Helios can be labeled as a warm seclusion, Juliette can be included in the tropical-like cyclone category.

Five-Year Analysis of Lightning Activities in Different Climatic Regions of Sichuan Province, China

Sichuan is a high-incidence area of thunderstorm activity in China. Based on the data of the total lightning location system from 2018 to 2022, the total lightning, cloud-to-ground (CG) lightning, and intracloud (IC) lightning activity regularity for the Sichuan province (SC) and its three climate subregions: Sichuan Basin (SB), Panxi district (PD), and West Sichuan Plateau (WSP) are analyzed, and the influences of different climate and topography conditions on lightning activities are also discussed. The results show that (1) for the whole province, the annual mean value of total lightning is about 850 thousand. The SB has the most lightning occurrences, and the WSP has the largest IC and +CG proportion. The southwest of PD, the north-center of PD, and the southeast of SB are the three high-value centers of lightning density. (2) For SB, the better thermodynamic and moisture conditions account for the most lightning occurrences. For PD, the lightning distribution is attributed to the joint effect of specific meteorological conditions and mountainous topography. For WSP, the convections are weak and shallow, which lead to high IC and +CG proportion. (3) The IC lightning mainly occurs below 12 km, and the IC height is much lower on WSP. The spatial and seasonal variation of IC height corresponds well to the cloud base height (CBH) and cloud top height (CTH). (4) The seasonal lightning frequency distribution in the three regions is similar, but the diurnal variation is quite different. The lightning activity mainly occurs between 1400 and 2200 LT on WSP but shows obvious nocturnal in SB. (5) Most CG intensity concentrates in the range below 50 kA, and IC concentrates in the range below 75 kA.

Cloud top heights and aerosol columnar properties from combined EarthCARE lidar and imager observations: the AM-CTH and AM-ACD products

Abstract. The Earth Cloud, Aerosol and Radiation Explorer (EarthCARE) is a combination of multiple active and passive instruments on a single platform. The Atmospheric Lidar (ATLID) provides vertical information of clouds and aerosol particles along the satellite track. In addition, the Multi-Spectral Imager (MSI) collects multi-spectral information from the visible to the infrared wavelengths over a swath width of 150 km across the track. The ATLID–MSI Column Products processor (AM-COL) described in this paper combines the high vertical resolution of the lidar along track and the horizontal resolution of the imager across track to better characterize a three-dimensional scene. ATLID Level 2a (L2a) data from the ATLID Layer Products processor (A-LAY), MSI L2a data from the MSI Cloud Products processor (M-CLD) and the MSI Aerosol Optical Thickness processor (M-AOT), and MSI Level 1c (L1c) data are used as input to produce the synergistic columnar products: the ATLID–MSI Cloud Top Height (AM-CTH) and the ATLID–MSI Aerosol Column Descriptor (AM-ACD). The coupling of ATLID (measuring at 355 nm) and MSI (at ≥670 nm) provides multi-spectral observations of the aerosol properties. In particular, the Ångström exponent from the spectral aerosol optical thickness (AOT 355/670 nm) adds valuable information for aerosol typing. The AOT across track, the Ångström exponent and the dominant aerosol type are stored in the AM-ACD product. The accurate detection of the cloud top height (CTH) with lidar is limited to the ATLID track. The difference in the CTH detected by ATLID and retrieved by MSI is calculated along track. The similarity of MSI pixels across track with those along track is used to transfer the calculated CTH difference to the entire MSI swath. In this way, the accuracy of the CTH is increased to achieve the EarthCARE mission's goal of deriving the radiative flux at the top of the atmosphere with an accuracy of 10 W m−2 for a 100 km2 snapshot view of the atmosphere. The synergistic CTH difference is stored in the AM-CTH product. The quality status is provided with the products. It depends, e.g., on day/night conditions and the presence of multiple cloud layers. The algorithm was successfully tested using the common EarthCARE test scenes. Two definitions of the CTH from the model truth cloud extinction fields are compared: an extinction-based threshold of 20 Mm−1 provides the geometric CTH, and a cloud optical thickness threshold of 0.25 describes the radiative CTH. The first CTH definition was detected with ATLID and the second one with MSI. The geometric CTH is always higher than or equal to the radiative CTH.

Research on Typhoon Multi-Stage Cloud Characteristics Based on Deep Learning

Analyzing the development and evolution characteristics of typhoons are conducive to improving typhoon monitoring and optimizing early warning models. Based on the deep learning model YOLOv5 and Himawari-8 data products, this study analyzes the movement path and cloud evolution of typhoon “Infa”. The specific conclusions of this study are as follows. (1) Based on the YOLOv5 model and brightness temperature perturbation algorithm, the central positioning of the typhoon is realized, where the Himawari-8 bright temperature image is used as the input of the model and the output of the model is the typhoon range boundary. The results show that this method was 90% accurate for monitoring ocular typhoons and 83% accurate for blind typhoons. The typhoon center location determined by the brightness temperature perturbation algorithm closely matched the CMA best-path dataset (CMA) (goodness of fit ≈0.99). (2) This study observed that as typhoons developed, cloud parameters evolved with the cloud cluster becoming denser. However, as the typhoon neared land, the cloud structure collapsed and cloud parameters decreased rapidly. (3) Changes in the typhoon cloud system were linked to topography and surface temperature. Changes in cloud optical thickness (COT) were influenced by the digital elevation model (correlation −0.18), while changes in cloud top temperature (CTT) and cloud top height (CTH) were primarily affected by surface temperature changes (correlation values: CTT −0.69, CTH −0.37). This suggests that the ocean environment supports the vertical development of typhoon clouds and precipitation. In summary, this study optimized the positioning simulation of typhoon movement paths and cloud change trends, and this is helpful for improving the early warning and response-ability of typhoons in coastal cities and for reducing the threat of typhoons to the daily life of residents in coastal areas.

Aerosol Effects on Water Cloud Properties in Different Atmospheric Regimes

AbstractAerosol‐cloud interaction remains one of the least understood processes in climate science arena, despite its profound impacts in radiation and water budget perturbations. The aerosol effects on clouds largely depend on aerosol characteristics. Here, we implemented 17‐year (2003–2020) data set of aerosol, cloud, and meteorological factors collected over East Asia—a highly polluted region with recent decreasing trend of air pollution due to control measures—to elucidate atmospheric regime‐dependent aerosol effects on water cloud properties by simultaneously accessing the response of air pollution control measures in cloud field. The study found a very close relationship between aerosol loading and cloud properties modifications in the continental region of East Asia with a significant response of air pollution control measures in the cloud field. The study further revealed that aerosols of the polluted continental atmosphere affected cloud micro‐ and macro‐physics differently than aerosols of the clean maritime atmosphere: in the former, increased aerosol loading increased the stability under cloud base and then enhanced cloud droplet collision‐coalescence process, resulting to increase cloud droplet size by decreasing cloud top height; whereas in the latter, increased aerosol loading decreased cloud droplet size without notable influence in the atmosphere thermodynamics and cloud top height. This study further showed a complex aerosol‐cloud interaction process in the polluted maritime atmosphere due to the mixed effect of polluted continental and clean maritime atmospheres. In all atmospheric regimes, cloud fraction was found to increase with the increase of aerosol loading.

Characterizing the near-global cloud vertical structures over land using high-resolution radiosonde measurements

Abstract. Cloud remains one of the largest uncertainties in weather and climate research due to the lack of fine-resolution observations of cloud vertical structure (CVS) on a large scale. In this study, near-global CVS is characterized by high-vertical-resolution twice-daily radiosonde observations from 374 stations over land, which are distributed in Europe, North America, East Asia, Australia, the Pacific Ocean, and Antarctica. To this end, we initially develop a novel method to determine CVS, by combining both the vertical gradients of air temperature and relative humidity (RH) and the altitude-dependent thresholds of RH. It is found that the cloud base heights (CBHs) from radiosondes have a higher correlation coefficient (R= 0.91) with the CBHs from a millimeter-wave cloud radar than those from the ERA5 reanalysis (R= 0.49). Overall, cloudy skies occur 65.3 % (69.5 %) of the time, of which 55.4 % (53.8 %) are one-layer clouds at 00:00 (12:00) UTC. Most multi-layer clouds are two-layer clouds, accounting for 62.2 % (61.1 %) among multi-layer clouds at 00:00 (12:00) UTC. Geographically, one-layer clouds tend to occur over arid regions, whereas two-layer clouds do not show any clear spatial preference. The cloud bases and tops over arid regions are higher compared with humid regions albeit with smaller cloud thickness (CT). Clouds tend to have lower bases and thinner layer thicknesses as the number of cloud layer increases. The global-mean CT, CBH, and cloud top height (CTH) are 4.89 ± 1.36 (5.37 ± 1.58), 3.15 ± 1.15 (3.07 ± 1.06), and 8.04 ± 1.60 (8.44 ± 1.52) km above ground level (a.g.l.) at 00:00 (12:00) UTC, respectively. The occurrence frequency of clouds is bimodal, with lower peaks between 0.5 and 3 km a.g.l. and upper peaks between 6 and 10 km a.g.l. The CBH, CTH, and CT undergo almost the same seasonality; namely, their magnitudes in boreal summer are greater than in boreal winter. As expected, the occurrence frequencies of clouds exhibit pronounced diurnal cycles in different seasons. In boreal summer, clouds tend to form as the sun rises and the occurrence frequencies increase from morning to late afternoon, with the peak in the early afternoon at the altitude of 6–12 km a.g.l., while in boreal winter, clouds have peak occurrence frequencies in the morning. The relations between surface meteorological variables and moisture with CBH are investigated as well, showing that CBHs are generally more significantly correlated with 2 m relative humidity (RH2 m) and 2 m air temperature (T2 m) than with surface pressure and 10 m wind speed. Larger T2 m and smaller RH2 m always correspond to higher CBH. In most cases CBHs are negatively correlated to soil water content. The near-global CVS obtained from high-vertical-resolution radiosondes in this study can provide key data support for improving the accuracy of cloud radiative forcing simulation in climate models.

CLAAS-3: the third edition of the CM SAF cloud data record based on SEVIRI observations

Abstract. CLAAS-3, the third edition of the Cloud property dAtAset using SEVIRI (Spinning Enhanced Visible and InfraRed Imager), was released in December 2022. It is based on observations from SEVIRI, on board geostationary satellites Meteosat-8, 9, 10 and 11, which are operated by the European Organisation for the Exploitation of Meteorological Satellites (EUMETSAT). CLAAS-3 was produced and released by the EUMETSAT Satellite Application Facility on Climate Monitoring (CM SAF), which aims to provide high-quality satellite-based data records suitable for climate monitoring applications. Compared to previous CLAAS releases, CLAAS-3 is expanded in terms of both temporal extent and cloud properties included, and it is based on partly updated retrieval algorithms. The available data span the period from 2004 to present, covering Europe; Africa; the Atlantic Ocean; and parts of South America, the Middle East and the Indian Ocean. They include cloud fractional coverage, cloud-top height, phase (liquid or ice) and optical and microphysical properties (water path, optical thickness, effective radius and droplet number concentration), from instantaneous data (every 15 min) to monthly averages. In this study we present an extensive evaluation of CLAAS-3 cloud properties, based on independent reference data sets. These include satellite-based retrievals from active and passive sensors, ground-based observations and in situ measurements from flight campaigns. Overall results show very good agreement, with small biases attributable to different sensor characteristics, retrieval/sampling approaches and viewing/illumination conditions. These findings demonstrate the fitness of CLAAS-3 to support the intended applications, which include evaluation of climate models, cloud characterisation and process studies focusing especially on the diurnal cycle and cloud filtering for other applications. The CLAAS-3 data record is publicly available via the CM SAF website at https://doi.org/10.5676/EUM_SAF_CM/CLAAS/V003 (Meirink et al., 2022).

Drivers of Cloud Condensation Nuclei in the Eastern North Atlantic as Observed at the ARM Site

AbstractSeven years of data collected at the Atmospheric Radiation Measurement's Eastern North Atlantic (ENA) site are analyzed to understand the controls of Cloud Condensation Nuclei (CCN) concentrations in the region. Day‐night differences in the aerosol data as segregated by wind direction demonstrate the aerosol observations to be impacted by local emissions when the wind direction (wdir) is between 90° and 310° (measured clockwise from the North where air is coming from). Data collected during marine conditions (wdir <90° or wdir >310°) show the CCN concentrations to be higher in the summer months as compared to the winter months. CCN budget analysis revealed advection and precipitation scavenging being primarily responsible for modulating the CCN concentrations at the site on monthly timescales, with rain rates driving the precipitation scavenging term. High (greater than 75th percentile) and low (lower than 25th percentile) CCN events were identified for each month to characterize the sub‐monthly variability of CCN concentrations. Low CCN events had thicker clouds, stronger rain rates, and lower reanalysis reported free‐tropospheric aerosol pseudo number concentration at the ENA site as compared to the high CCN events. Analysis of satellite data of air‐parcels 48 hr prior to their arrival at the ENA site demonstrated the air parcels during low CCN events to encounter higher cloudiness, stronger rain rates, and higher cloud top heights as compared to the high CCN events. The results presented herein provide key constraints for model evaluation studies and climatological studies conducted at the ENA site.

Characteristics and Variability of Precipitation Across Different Sectors of an Extra‐Tropical Cyclone: A Case Study Over the High‐Latitudes of the Southern Ocean

AbstractShipborne observations from the CAPRICORN‐2018 field campaign were used to investigate the characteristics and variability of precipitation across different sectors of an extra‐tropical cyclone on 16 February 2018, over the Southern Ocean (SO). Three distinct time periods—frontal, post‐frontal, and cyclone—were identified during the day. The frontal passage recorded a total accumulation of 1.9 mm, where the precipitation phases were primarily composed of rain (96%), while the cyclone period recorded the largest precipitation (4.0 mm), where the precipitation phases varied with snow (10%), mixed‐phase (40%), and rain (50%). The BASTA radar suggests the freezing level was shallow (∼500 m) with snow present above. The cloud top heights, observed by a C‐band radar, were shallower in the cyclone period, although deeper cloud depths of ∼6 km were sporadically recorded. Increased surface fluxes and a southerly wind direction indicate that cold air advection, was likely the main cause of high precipitation during the cyclone period. A non‐precipitating multi‐layer cloud structure with a geometrically thin (200 m) homogeneous layer of supercooled liquid water (SLW) overlaying shallow boundary layer convection was seen during the post‐frontal period. The ship‐borne observations were used to evaluate Weather Research & Forecasting (WRF) simulations with different microphysics settings. We found the frontal precipitation intensity is well reproduced, but it is underestimated during the cyclone period. This study represents a unique set of observations and highlights the need for understanding how ice processes and potentially horizontal advection contribute to the development of precipitation and convection over the SO.

The evolution of deep convective systems and their associated cirrus outflows

Abstract. Tropical deep convective clouds, particularly their large cirrus outflows, play an important role in modulating the energy balance of the Earth's atmosphere. Whilst the cores of these deep convective clouds have a significant short-wave (SW) cooling effect, they dissipate quickly. Conversely, the thin cirrus that flow from these cores can persist for days after the core has dissipated, reaching hundreds of kilometres in extent. These thin cirrus have a potential for large warming in the tropics. Understanding the evolution of air parcels from deep convection, clouds along these trajectories, and how they change in response to anthropogenic emissions is therefore important to understand past and future climate change. This work uses a novel approach to investigate the evolution of tropical convective clouds by introducing the concept of “time since convection” (TSC). This is used to build a composite picture of the lifecycle of air parcels from deep convection. Cloud properties are a strong function of TSC, showing decreases in the optical thickness, cloud-top height, and cloud fraction over time, thereby driving the latitudinal structure of cloudiness. After an initial dissipation of the convective core, changes in thin cirrus cloud amount were seen beyond 200 h from convection. Changes in cloud are shown to be a strong function of TSC and not simply reflective of latitudinal changes as air moves from the tropics to the extratropics. Finally, in the initial stages of convection there was a large net negative cloud radiative effect (CRE). However, once the convective core had dissipated, the sign of the CRE flipped and there was a sustained net warming CRE beyond 120 h from the convective event. Changes are present in the cloud properties long after the main convective activities have dissipated, signalling the need to continue further analysis at longer timescales than previously studied.

Historical (1960–2014) lightning and LNOx trends and their controlling factors in a chemistry–climate model

Abstract. Lightning can cause natural hazards that result in human and animal injuries and fatalities, infrastructure destruction, and wildfire ignition. Lightning-produced NOx (LNOx), a major NOx (NOx=NO+NO2) source, plays a vital role in atmospheric chemistry and global climate. The Earth has experienced marked global warming and changes in aerosol and aerosol precursor emissions (AeroPEs) since the 1960s. Investigating long-term historical (1960–2014) lightning and LNOx trends can provide important indicators for all lightning-related phenomena and for LNOx effects on atmospheric chemistry and global climate. Understanding how global warming and changes in AeroPEs influence historical lightning and LNOx trends can be helpful in providing a scientific basis for assessing future lightning and LNOx trends. Moreover, global lightning activities' responses to large volcanic eruptions such as the 1991 Pinatubo eruption are not well elucidated and are worth exploring. This study employed the widely used cloud top height lightning scheme (CTH scheme) and the newly developed ice-based ECMWF-McCAUL lightning scheme to investigate historical (1960–2014) lightning and LNOx trends and variations as well as their influencing factors (global warming, increases in AeroPEs, and the Pinatubo eruption) in the framework of the CHASER (MIROC) chemistry–climate model. The results of the sensitivity experiments indicate that both lightning schemes simulated almost flat global mean lightning flash rate anomaly trends during 1960–2014 in CHASER (the Mann–Kendall trend test (significance inferred as 5 %) shows no trend for the ECMWF-McCAUL scheme, but a 0.03 % yr−1 significant increasing trend is detected for the CTH scheme). Moreover, both lightning schemes suggest that past global warming enhances historical trends for global mean lightning density and global LNOx emissions in a positive direction (around 0.03 % yr−1 or 3 % K−1). However, past increases in AeroPEs exert an opposite effect on the lightning and LNOx trends (−0.07 % to −0.04 % yr−1 for lightning and −0.08 % to −0.03 % yr−1 for LNOx) when one considers only the aerosol radiative effects in the cumulus convection scheme. Additionally, effects of past global warming and increases in AeroPEs in lightning trends were found to be heterogeneous across different regions when analyzing lightning trends on the global map. Lastly, this paper is the first of study results suggesting that global lightning activities were markedly suppressed during the first year after the Pinatubo eruption as shown in both lightning schemes (global lightning activities decreased by as much as 18.10 % as simulated by the ECMWF-McCAUL scheme). Based on the simulated suppressed lightning activities after the Pinatubo eruption, the findings also indicate that global LNOx emissions decreased after the 2- to 3-year Pinatubo eruption (1.99 %–8.47 % for the annual percentage reduction). Model intercomparisons of lightning flash rate trends and variations between our study (CHASER) and other Coupled Model Intercomparison Project Phase 6 (CMIP6) models indicate great uncertainties in historical (1960–2014) global lightning trend simulations. Such uncertainties must be investigated further.

A Parallax Shift Effect Correction Based on Cloud Top Height for FY-4A Lightning Mapping Imager (LMI)

The Lightning Mapping Imager (LMI) onboard the Fengyun-4A (FY-4A) satellite is the first independently developed satellite-borne lightning imager in China. It enables continuous lightning detection in China and surrounding areas, regardless of weather conditions. The FY-4A LMI uses a Charge-Coupled Device (CCD) array for lightning detection, and the accuracy of lightning positioning is influenced by cloud top height (CTH). In this study, we proposed an ellipsoid CTH parallax correction (ECPC) model for lightning positioning applicable to FY-4A LMI. The model utilizes CTH data from the Advanced Geosynchronous Radiation Imager (AGRI) on FY-4A to correct the lightning positioning data. According to the model, when the CTH is 12 km, the maximum deviation in lightning positioning caused by CTH in Beijing is approximately 0.1177° in the east–west direction and 0.0530° in the north–south direction, corresponding to a horizontal deviation of 13.1558 km, which exceeds the size of a single ground detection unit of the geostationary satellite lightning imager. Therefore, it is necessary to be corrected. A comparison with data from the Beijing Broadband Lightning Network (BLNET) and radar data shows that the corrected LMI data exhibit spatial distribution that is closer to the simultaneous BLNET lightning positioning data. The coordinate differences between the two datasets are significantly reduced, indicating higher consistency with radar data. The correction algorithm decreases the LMI lightning location deviation caused by CTH, thereby improving the accuracy and reliability of satellite lightning positioning data. The proposed ECPC model can be used for the real-time correction of lightning data when CTH is obtained at the same time, and it can be also used for the post-correction of space-based lightning detection with other cloud top height data.

Estimation of downwelling surface longwave radiation for cloudy skies by considering the radiation effect from the entire cloud layers

Cloud base height (CBH) is a key parameter to characterize the cloud radiation effect. However, the CBH used in downwelling surface longwave radiation (DSLR) estimation is generally obtained indirectly through cloud top parameters retrieved by passive optical remote sensing instruments, which is of high uncertainty. At the same time, it is unreasonable to replace the effective radiation of the entire cloud layers by only using the cloud base radiation with single layer cloud model. This study proposes a new method to estimate cloudy-sky DSLR, which considers the radiation effect of the entire cloud layers from the cloud base to top. First, the CBH estimation model is established by the genetic algorithm-artificial neural network (GA-ANN) algorithm. The cloud top height and cloud attribute parameters (cloud optical depth, cloud water path, and cloud phase) from the passive remote sensing are used as the input features, and meanwhile CBH data from the active remote sensing are output features in the training and testing process of the model. Then, the cloud base temperature (CBT) is estimated based on the CBH combined with the temperature profile data in the EAR5 reanalysis data. Finally, the effective temperature of the entire cloud layers is calculated by using CBT and cloud top temperature. The verification results of CBH estimation showed that R2 is 0.83, the bias and root mean square error (RMSE) are 0.02 km and 1.56 km, respectively, which indicates a comparable accuracy and higher stability compared with the previous studies. The ground-based measurements in the SURFRAD network are used to validate the newly proposed DSLR estimation method, and the results showed that the bias and RMSE are 5.27 W/m2 and 28.48 W/m2, respectively. Additionally, this study found that although the effective temperature of the entire cloud layers has a weaker linear correlation with DSLR, the radiation contribution generated by cloud still occupies a certain weight, and the maximum ratio of cloud radiation in DSLR estimation can account for 30%. Therefore, the cloud radiation effect must be taken into account in the estimation of cloudy-sky DSLR.

Observations and modelling of the winter thunderstorm on 4 February 2022 at the Milešovka meteorological observatory

AbstractThe study analyses a winter thunderstorm that passed over the Milešovka meteorological observatory on 4 February 2022, between 2300 and 2330 UTC. Lightning was recorded directly over the observatory by both the observer and the EUCLID lightning network at 2320 UTC. To analyse the state of the atmosphere at the time when the lightning occurred, we used data from the X‐band Doppler polarimetric radar and the Ka‐band Doppler polarimetric vertical profiler, both located at the observatory. We also applied data from the Meteosat Second Generation satellite, and data from standard meteorological instruments located at the observatory. In addition, we run our cloud electrification model to simulate cloud electrification of the winter thunderstorm to find out whether the model develops conditions suitable for the occurrence of lightning and if so, under what circumstances. Our results show that the lightning appeared at the very end of the storm passage defined by high radar reflectivity. At the same time, it is clear from the radar observations that before lightning occurred, the cloud contained hydrometeors (graupel, cloud or rain water, and ice or snow) which are commonly associated with charge separation by collisions. Our analysis of the radar data also suggests that in at least several parts of the cloud the electric field was strong. Although the cloud top height was very low compared to summer storms, the model results indicate conditions suitable for lightning occurrence. However, uncertainty remains on how to properly formulate the initial conditions for model simulations for this type of storm which was shallow and occurs rarely in winter.

Cloud top heights and aerosol layer properties from EarthCARE lidar observations: the A-CTH and A-ALD products

Abstract. The high-spectral-resolution Atmospheric Lidar (ATLID) on the Earth Cloud, Aerosol and Radiation Explorer (EarthCARE) provides vertically resolved information on aerosols and clouds with unprecedented accuracy. Together with the Cloud Profiling Radar (CPR), the Multi-Spectral Imager (MSI), and the Broad-Band Radiometer (BBR) on the same platform, it allows for a new synergistic view on atmospheric processes related to the interaction of aerosols, clouds, precipitation, and radiation at the global scale. This paper describes the algorithms for the determination of cloud top height and aerosol layer information from ATLID Level 1b (L1b) and Level 2a (L2a) input data. The ATLID L2a Cloud Top Height (A-CTH) and Aerosol Layer Descriptor (A-ALD) products are developed to ensure the provision of atmospheric layer products in continuation of the heritage from the Cloud–Aerosol Lidar and Infrared Pathfinder Satellite Observations (CALIPSO). Moreover, the products serve as input for synergistic algorithms that make use of data from ATLID and MSI. Therefore, the products are provided on the EarthCARE joint standard grid (JSG). A wavelet covariance transform (WCT) method with flexible thresholds is applied to determine layer boundaries from the ATLID Mie co-polar signal. Strong features detected with a horizontal resolution of 1 JSG pixel (approximately 1 km) or 11 JSG pixels are classified as thick or thin clouds, respectively. The top height of the uppermost cloud layer together with information on cloud layering are stored in the A-CTH product for further use in the generation of the ATLID-MSI Cloud Top Height (AM-CTH) synergy product. Aerosol layers are detected as weaker features at a resolution of 11 JSG pixels. Layer-mean optical properties are calculated from the ATLID L2a Extinction, Backscatter and Depolarization (A-EBD) product and stored in the A-ALD product, which also contains the aerosol optical thickness (AOT) of each layer, the stratospheric AOT, and the AOT of the entire atmospheric column. The latter parameter is used to produce the synergistic ATLID-MSI Aerosol Column Descriptor (AM-ACD) later in the processing chain. Several quality criteria are applied in the generation of A-CTH and A-ALD, and respective information is stored in the products. The functionality and performance of the algorithms are demonstrated by applying them to common EarthCARE test scenes. Conclusions are drawn for the application to real-world data and the validation of the products after the launch of EarthCARE.

Earth-system-model evaluation of cloud and precipitation occurrence for supercooled and warm clouds over the Southern Ocean's Macquarie Island

Abstract. Over the remote Southern Ocean (SO), cloud feedbacks contribute substantially to Earth system model (ESM) radiative biases. The evolution of low Southern Ocean clouds (cloud-top heights < ∼ 3 km) is strongly modulated by precipitation and/or evaporation, which act as the primary sink of cloud condensate. Constraining precipitation processes in ESMs requires robust observations suitable for process-level evaluations. A year-long subset (April 2016–March 2017) of ground-based profiling instrumentation deployed during the Macquarie Island Cloud and Radiation Experiment (MICRE) field campaign (54.5∘ S, 158.9∘ E) combines a 95 GHz (W-band) Doppler cloud radar, two lidar ceilometers, and balloon-borne soundings to quantify the occurrence frequency of precipitation from the liquid-phase cloud base. Liquid-based clouds at Macquarie Island precipitate ∼ 70 % of the time, with deeper and colder clouds precipitating more frequently and at a higher intensity compared to thinner and warmer clouds. Supercooled cloud layers precipitate more readily than layers with cloud-top temperatures > 0 ∘C, regardless of the geometric thickness of the layer, and also evaporate more frequently. We further demonstrate an approach to employ these observational constraints for evaluation of a 9-year GISS-ModelE3 ESM simulation. Model output is processed through the Earth Model Column Collaboratory (EMC2) radar and lidar instrument simulator with the same instrument specifications as those deployed during MICRE, therefore accounting for instrument sensitivities and ensuring a coherent comparison. Relative to MICRE observations, the ESM produces a smaller cloud occurrence frequency, smaller precipitation occurrence frequency, and greater sub-cloud evaporation. The lower precipitation occurrence frequency by the ESM relative to MICRE contrasts with numerous studies that suggest a ubiquitous bias by ESMs to precipitate too frequently over the SO when compared with satellite-based observations, likely owing to sensitivity limitations of spaceborne instrumentation and different sampling methodologies for ground- versus space-based observations. Despite these deficiencies, the ESM reproduces the observed tendency for deeper and colder clouds to precipitate more frequently and at a higher intensity. The ESM also reproduces specific cloud regimes, including near-surface clouds that account for ∼ 25 % of liquid-based clouds during MICRE and optically thin, non-precipitating clouds that account for ∼ 27 % of clouds with bases higher than 250 m. We suggest that the demonstrated framework, which merges observations with appropriately constrained model output, is a valuable approach to evaluate processes responsible for cloud radiative feedbacks in ESMs.

Analysis of lightning flash rate with the occurrence of flash floods and hailstorms in Peninsular Malaysia

This paper will present five flash flood cases in Johor and two hailstorms from Kuala Lumpur and Melaka reported during lightning activities to determine the relationship between lightning flash rate and rainfall that occurred from 2013 to 2017. The flash floods were observed based on lightning data obtained from Tenaga Nasional Berhad Research Sdn. Bhd. (TNBR) and rainfall data from the Department of Irrigation and Drainage (DID). On the other hand, storm-producing hails were observed based on the lightning flash rate using electromagnetic field measurements and weather radar for their reflectivity value and cloud structure. The detection of a lightning flash, especially the negative narrow bipolar event (−NBE) and positive cloud-to-ground (+CG), could be a possible indicator of a deeper convective storm that could lead to the production of hails. In the 24-h observation of selected days when a flash flood occurred, lightning strikes correlation with rainfall had r-values ranging from 0.6 to 0.9, where the lightning data may enhance flood forecasting accuracy when combined with rainfall. Meanwhile, during the first hailstorms in Kuala Lumpur, the highest −NBEs flash rate detected during the hails appeared where 16 −NBEs were detected in a span of 5 min with reflectivity value between 53 and 55 dBz and cloud top height of 20 km. During the second hailstorms in Melaka, more + CG flashes rather than −NBE were detected during the appearance of the hail, with 64 +CG flashes and five −NBEs with reflectivity values between 60 and 66 dBz and cloud top height of 18.8 km. Based on the high convective strength, cloud top height, and presence of hailstones, −NBE and +CG flash could be possible indicators of storm severity, especially in Malaysia.

Parameterization of Entrainment Rate for Cumulus Clouds with WRF Simulation

By using Weather Research and Forecasting Model (WRF) to simulate a southwest vortex precipitation process, this work studies the correlations between entrainment rate (λ) and dynamical parameters in the cloud and further fit λ. We relate the probability density distribution (PDF) to the parameterization of λ and find that the greater the probability, the larger the slope of the logarithmic liner function. The slope of the log-linear fitting function in fitting decreases for developing and enhancing cumulus clouds, which is related to the increase in updraft motion and the decrease in λ. Then, we group clouds according to cloud top heights and calculate average λ and dynamic parameters, and the results indicate that when only one dynamic parameter is used, vertical wind velocity (w) is more suitable than buoyancy (B) to be used to fit λ. The fitting functions combing one single parameter and more parameters by principal components regression are compared with two traditional schemes, and we found that λ obtained by our fitting schemes are between the two traditional schemes. Because the principal component regression method takes into account the interaction between more dynamic factors and entrainment, the fitting function, including w and B, is suitable to be applied to fit λ in the parameterization scheme for cumulus clouds.