Cloud Base Height Research Articles

Analysis of Macro- and Microphysical Characteristics of Ice Clouds over the Tibetan Plateau Using CloudSat/CALIPSO Data

Utilizing CloudSat/CALIPSO satellite data and ERA5 reanalysis data from 2007 to 2016, this study analyzed the distributions of optical and physical characteristics and change characteristics of ice clouds over the Tibetan Plateau (TP). The results show that the frequency of ice clouds in the cold season (November to March) on the plateau is over 80%, while in the warm season (May to September) it is around 60%. The average cloud base height of ice clouds in the warm season is 3–5 km, and mostly around 2 km in the cold season. The average cloud top height in the warm season is around 5–8 km, while in the cold season it is mainly around 4.5 km. The average thickness of ice clouds in both seasons is around 2 km. The statistical results of microphysical characteristics show that the ice water content is around 10−1 to 103 mg/m3, and the effective radius of ice clouds is mainly in the range of 10–90 μm. Both have their highest frequency in the west of the TP and lowest in the northeast. A comprehensive analysis of the change in temperature, water vapor, and ice cloud occurrence frequency shows that the rate of increase in water vapor in the warm season is greater than that in the cold season, while the rates of increase in both surface temperature and ice cloud occurrence are smaller than in the cold season. The rate of increase in temperature in the warm season is around 0.038 °C/yr, and that in the cold season is around 0.095 °C/yr. The growth rate of thin ice clouds in the warm season is around 0.15% per year, while that in the cold season is as high as 1% per year. The results suggest that the surface temperature change may be related to the occurrence frequency of thin ice clouds, with the notable increase in temperature during the cold season possibly being associated with a significant increase in the occurrence frequency of thin ice clouds.

Cloud classification through machine learning and global horizontal irradiance data analysis

AbstractCloud observations and characterization are crucial owing to their influence on energy balance, climate, and weather. Their particular effects on radiation vary depending on different cloud parameters, such as cloud base or top height, water content, and cloud optical thickness, all of them closely related to the specific cloud type. Cloud classification therefore becomes a crucial task in meteorology, although it remains challenging for weather services worldwide owing to the intensive associated labor and cost. In this study we introduce a new low‐cost method for automating cloud classification based on a combination of ground‐based global horizontal irradiance (GHI) measurements, a clear‐sky model, and machine learning. Based on the hypothesis that different cloud types have their own GHI signatures, we trained different supervised learning algorithms using GHI data manually labeled by meteorological observers from time‐synchronized all‐sky images. Multiple time windows were extracted from each GHI series, with eight features defined in each case to characterize the sequence. The best outcome was achieved using an XGBoost model on features extracted on time windows of 33 min, obtaining an accuracy of 0.88 and a Cohen's kappa of 0.84 in a held‐out test set. The development presented in this study has the ability to provide low‐cost cloud classification from ground‐based observations, which is a challenge for weather services worldwide owing to intensive labor and cost.

Cloud base height variability observed using a Laser-Based Ceilometer over a tropical station Pune, India

ABSTRACT Clouds play a pivotal role in determining Earth’s climate, impacting several atmospheric processes influencing the weather. The study explores the temporal variations in cloud base height (CBH) over Pune, a tropical station, from 2017 to 2023 using a Laser-Based Ceilometer with a range of 0–7500 m. Our analysis shows seasonal trends in CBH, with the monsoon season exhibiting unique characteristics showing the lowest average CBH and minimal diurnal changes during the study period. Conversely, the pre-monsoon, post-monsoon, and winter season show a gradual rise of CBH in the morning, maintaining higher levels until late afternoon, followed by a decline in the evening. Analysis of the diurnal frequency of cloud occurrence at 500 m intervals reveals the monsoon season’s dominance in the 500 m to 1000 m layer. Throughout the year, low cloud bases (CBH <2000 m) are most frequent, outpacing the frequency of mid and high cloud bases. High cloud base (7500 < CBH > 6000 m) is the least common, showing a semi-diurnal pattern in winter with peaks around noon and midnight. These findings not only enhance our understanding of Pune’s atmospheric behaviour but also offer valuable information for weather prediction, climate modelling, and environmental monitoring.

Cirrus cloud occurrence and its geometrical and optical properties during the transient monsoon conditions

Knowledge of variation in the percentage occurrence of the cirrus clouds (POC) during transient monsoon conditions is essential for understanding the role of the monsoon in transporting the water vapor into the lower stratosphere which is vital in quantifying the radiation budget of the earth–atmosphere system. In this paper, we present the spatial structure of the POC, the geometrical properties such as cloud top and base height (CTH & CBH), cloud thickness (CTH-CBH), optical properties such as optically thick, thin, and subvisible cirrus clouds during the active and break phases of the Asian summer monsoon using Cloud-Aerosol Lidar and Infrared Pathfinder Satellite Observation (CALIPSO) during July–August 2006–2018. The active and break phases are identified based on the central India rainfall from the India Meteorological Department dataset. The POC is found to be higher by 10–30 % over central India and some parts of the eastern Arabian Sea (AS) and Bay of Bengal (BoB) during the active phase compared to the break phase which shows enhancement in POC over the foothill of Himalaya and southern part of BoB and AS. The higher POC is attributed to the increase in the convective activities associated with the increase in tropical easterly jet during the active phase. The CTH and CBH range ~14–18 km and 12–16 km, respectively during the active phase compared to the break phase. The contributions of optically thick and thin clouds are higher in POC during both phases compared to subvisible clouds. Both convectively generated and in situ-formed cirrus clouds dominate during the active phase compared to the break phase indicating the dominance of the freeze-drying processes during the active phase.

Assessing Downburst Kinematics Using Video Footage Analysis

Measurements of downburst outflows using standard meteorological instruments (e.g., anemometers) are rare due to their transient and localized nature. However, video recordings of such events are becoming more frequent. This short communication (Technical Note) study presents a new approach to estimating the kinematics of a downburst event using video footage recordings of the event. The main geometric dimensions of the event, such as downdraft diameter, cloud base height, outflow depth, and the radius of the outflow at a given moment in time, are estimated by sizing them against reference structures of known dimensions that are present in the video footage. From this analysis, and knowing the frame rate of the video recording, one can estimate the characteristic velocities in the downburst event, such as the mean downdraft velocity and the mean velocity of the radial outflow propagation. The proposed method is tested on an August 2015 downburst event that hit Tucson, Arizona, United States. The diameter of the downburst outflow increased with the time from approximately 1.10 km to 3.35 km. This range of values indicates that the event was a microburst. The mean descending velocity of downburst downdraft was 8.9 m s−1 and the horizontal velocity of outflow propagation was 17.7 m s−1. The latter velocity is similar to the measured wind gust at the nearby weather station and Doppler radar. The outflow depth is estimated at 160 m, and the cloud base height was approximately 1.24 km. Estimating the kinematics of downbursts using video footage, while subject to certain limitations, does yield a useful estimation of the main downburst kinematics that contribute to a better quantification of these localized windstorms.

Improving Atmospheric Temperature and Relative Humidity Profiles Retrieval Based on Ground-Based Multichannel Microwave Radiometer and Millimeter-Wave Cloud Radar

Obtaining temperature and humidity profiles with high vertical resolution is essential for describing and predicting atmospheric motion, and, in particular, for understanding the evolution of medium- and small-scale weather processes, making short-range and near-term weather forecasting, and implementing weather modifications (artificial rainfall, artificial rain elimination, etc.). Ground-based microwave radiometers can acquire vertical tropospheric atmospheric data with high temporal and spatial resolution. However, the accuracy of temperature and relative humidity retrieval is still not as accurate as that of radiosonde data, especially in cloudy conditions. Therefore, improving the observation and retrieval accuracy is a major challenge in current research. The focus of this study was to further improve the accuracy of atmospheric temperature and humidity profile retrieval and investigate the specific effects of cloud information (cloud-base height and cloud thickness) on temperature and humidity profile retrieval. The observation data from the ground-based multichannel microwave radiometer (GMR) and the millimeter-wave cloud radar (MWCR) were incorporated into the retrieval process of the atmospheric temperature and relative humidity profiles. The retrieval was performed using the backpropagation neural network (BPNN). The retrieval results were quantified using the mean absolute error (MAE) and root mean square error (RMSE). The statistical results showed that the temperature profiles were less affected by the cloud information compared with the relative humidity profiles. Cloud thickness was the main factor affecting the retrieval of relative humidity profiles, and the retrieval with cloud information was the best retrieval method. Compared with the retrieval profiles without cloud information, the MAE and RMSE values of most of the altitude layers were reduced to different degrees after adding cloud information, and the relative humidity (RH) errors of some altitude layers were reduced by approximately 50%. The maximum reduction in the RMSE and MAE values for the retrieval of temperature profiles with cloud information was about 1.0 °C around 7.75 km, and the maximum reduction in RMSE and MAE values for the relative humidity profiles was about 10%, which was obtained around 2 km.

Deep-learning-driven simulations of boundary layer clouds over the Southern Great Plains

Abstract. Based on long-term observations at the Southern Great Plains site by the Atmospheric Radiation Measurement (ARM) program for training and validation, a deep-learning model is developed to simulate the daytime evolution of boundary layer clouds (BLCs) from the perspective of land–atmosphere coupling. The model takes ARM measurements (including early-morning soundings and diurnally varying surface meteorological conditions and heat fluxes) as inputs and predicts hourly estimates (including cloud occurrence, the positions of cloud boundaries, and the vertical profile of the cloud fraction) as outputs. The deep-learning model offers good agreement with the observed cloud fields, especially in the accuracy with which cloud occurrence and base height are reproduced. When the inputs are substituted by reanalysis data from ERA5 and MERRA-2, the outputs of the deep-learning model provide a better agreement with observation than the cloud fields extracted from ERA5 and MERRA-2 themselves. Thus, the deep-learning model shows great potential to serve as a diagnostic tool for the performance of physics-based models in simulating stratiform and cumulus clouds. By quantifying biases in clouds and attributing them to the simulated atmospheric state variables versus the model-parameterized cloud processes, this observation-based deep-learning model may offer insights into the directions needed to improve the simulation of BLCs in physics-based models for weather forecasting and climate prediction.

Ampycloud: an open-source algorithm to determine cloud base heights and sky coverage fractions from ceilometer data

Abstract. Ceilometers are used routinely at aerodromes worldwide to derive the height and sky coverage fraction of cloud layers. This information, possibly combined with direct observations by human observers, contributes to the production of meteorological aerodrome reports (METARs). Here, we present ampycloud, a new algorithm, and its associated Python package for automatic processing of ceilometer data with the aim of determining the sky coverage fraction and base height of cloud layers above aerodromes. The ampycloud algorithm was developed at the Swiss Federal Office of Meteorology and Climatology (MeteoSwiss) as part of the AMAROC (AutoMETAR/AutoReport rOund the Clock) program to help in the fully automatic production of METARs at Swiss civil aerodromes. ampycloud is designed to work with no direct human supervision. The algorithm consists of three distinct, sequential steps that rely on agglomerative clustering methods and Gaussian mixture models to identify distinct cloud layers from individual cloud base hits reported by ceilometers. The robustness of the ampycloud algorithm stems from the first processing step, which is simple and reliable. It constrains the two subsequent processing steps that are more sensitive but also better suited to handling complex cloud distributions. The software implementation of the ampycloud algorithm takes the form of an eponymous, pip-installable Python package developed on GitHub and made publicly accessible.

A note of double-layer cloud detection method based on cloud base height and brightness

Double-layer cloud is a special cloud condition, and clouds in two layers have different effects on solar radiation reaching ground. Getting the shape of double-layer clouds is helpful to study the change of solar irradiance. The proposed double-layer clouds detection method is based on the cloud base height (CBH) and brightness of double-layer clouds. The CBH is calculated based on binocular parallax, which requires two images filmed synchronously by two cameras, and the brightness is offered by the brightness channel in HSV color space. The accuracy of this method is 98.5% for single-layer clouds and 89.4% for double-layer clouds. The reasons for misidentifications are (1) the brightness contrast decreases when the high-layer clouds are thick, (2) the distant low clouds have different brightness in top and bottom.

Deforestation amplifies climate change effects on warming and cloud level rise in African montane forests

Tropical montane forest ecosystems are pivotal for sustaining biodiversity and essential terrestrial ecosystem services, including the provision of high-quality fresh water. Nonetheless, the impact of montane deforestation and climate change on the capacity of forests to deliver ecosystem services is yet to be fully understood. In this study, we offer observational evidence demonstrating the response of air temperature and cloud base height to deforestation in African montane forests over the last two decades. Our findings reveal that approximately 18% (7.4 ± 0.5 million hectares) of Africa’s montane forests were lost between 2003 and 2022. This deforestation has led to a notable increase in maximum air temperature (1.37 ± 0.58 °C) and cloud base height (236 ± 87 metres), surpassing shifts attributed solely to climate change. Our results call for urgent attention to montane deforestation, as it poses serious threats to biodiversity, water supply, and ecosystem services in the tropics.

Evaluation and Improvement of a CALIPSO-Based Algorithm for Cloud Base Height in China

Clouds are crucial in regulating the Earth’s energy budget. Global cloud top heights have been easily retrieved from satellite measurements, but there are few methods for determining cloud base height (CBH) from satellite measurements. The Cloud Base Altitude Spatial Extrapolator (CBASE) algorithm was proposed to derive the height of the lower-troposphere liquid cloud base by using the Cloud-Aerosol Lidar with Orthogonal polarization cloud aerosol LiDAR (CALIOP) profiles and weather observations at airports from aviation routine and special weather report (METARs and SPECIs, called METAR) observation data in the United States. A modification to the CBASE algorithm over China (CNMETAR-CBASE) is presented in this paper. In this paper, the ability of the CBASE algorithm to calculate CBH in China is evaluated, and METAR observations over China (CNMETAR) were then used to modify the CBASE algorithm. The results including CNMETAR observation data in China can better retrieve CBH over China compared with the results using the original CBASE algorithm, and the accuracy of the global CBH results has been improved. Overestimations of CBH with the original algorithm range from 500 to 800 m in China, which have been reduced to about 300 m with an improved algorithm. The deviations calculated by the algorithm also have a significant reduction, from 480 m (CBASE) to 420 m (CNMETAR-CBASE). In conclusion, the modified CBASE algorithm not only calculates the CBH more accurately in China but also improves the results of the global CBH retrieved from satellites.

Multi-source observations on the effect of atmospheric blocking on air quality in İstanbul: a study case.

Air pollution is affected by the atmospheric dynamics. This study aims to determine that air pollution concentration values in İstanbul increased significantly and reached peak values due to atmospheric blocking between the 30th of December 2022 and the 5th of January 2023. In this study, hourly pollutant data was obtained from 16 air quality monitoring stations (AQMS), the exact reanalysis data was extracted from ERA5 database, and inversion levels and meteorological and synoptic analyses were used to determine the effects of atmospheric blocking on air pollution. Also, cloud base heights and vertical visibility measurements were taken with a ceilometer. Statistical calculations and data visualizations were performed using the R and Grads program. Omega-type blocking, which started in İstanbul on December 30, 2022, had a significant impact on the 1st and 2nd of January 2023, and PM10 and PM2.5 concentration values reached their peak values at 572.8 and 254.20µg/m3, respectively. In addition, it was found that the average concentration values in the examined period in almost all stations were higher than the averages for January and February. As a result, air quality in İstanbul was determined as "poor" between these calendar dates. It was found that the blocking did not affect the ozone (µg/m3) concentration. It was also found that the concentrations of particulate matter (PM) 10µm or less in diameter (PM10) and PM 2.5µm or less in diameter (PM2.5) were increased by the blocking effect in the İstanbul area. Finally, according to the data obtained using the ceilometer, cloud base heights decreased to 30m and vertical visibility to 10m.

Sensitivity of thermodynamic profiles retrieved from ground-based microwave and infrared observations to additional input data from active remote sensing instruments and numerical weather prediction models

Abstract. Accurate and continuous estimates of the thermodynamic structure of the lower atmosphere are highly beneficial to meteorological process understanding and its applications, such as weather forecasting. In this study, the Tropospheric Remotely Observed Profiling via Optimal Estimation (TROPoe) physical retrieval is used to retrieve temperature and humidity profiles from various combinations of input data collected by passive and active remote sensing instruments, in situ surface platforms, and numerical weather prediction models. Among the employed instruments are microwave radiometers (MWRs), infrared spectrometers (IRSs), radio acoustic sounding systems (RASSs), ceilometers, and surface sensors. TROPoe uses brightness temperatures and/or radiances from MWRs and IRSs, as well as other observational inputs (virtual temperature from the RASS, cloud-base height from the ceilometer, pressure, temperature, and humidity from the surface sensors) in a physical iterative retrieval approach. This starts from a climatologically reasonable profile of temperature and water vapor, with the radiative transfer model iteratively adjusting the assumed temperature and humidity profiles until the derived brightness temperatures and radiances match those observed by the MWR and/or IRS instruments within a specified uncertainty, as well as within the uncertainties of the other observations, if used as input. In this study, due to the uniqueness of the dataset that includes all the abovementioned sensors, TROPoe is tested with different observational input combinations, some of which also include information higher than 4 km above ground level (a.g.l.) from the operational Rapid Refresh numerical weather prediction model. These temperature and humidity retrievals are assessed against independent collocated radiosonde profiles under non-cloudy conditions to assess the sensitivity of the TROPoe retrievals to different input combinations.

Characteristics and seasonal variations of cirrus clouds from ground-based lidar and satellite observations over Shouxian Area, China

This study investigates the macroscopic and optical properties of cirrus clouds in the 32N region from July 2016 to May 2017, leveraging data from ground-based lidar observations and CALIOP to overcome the inconsistencies in detected cirrus cloud samples. Through extensive data analysis, statistical characteristics of cirrus clouds were discerned, revealing lidar ratio values of 28.5 ± 10.8 from ground-based lidar and 27.4 ± 11.2 from CALIOP. Validation with a decade of CALIOP data (2008-2018) confirmed these findings, presenting a consistent lidar ratio of 27.4 ± 12.0. A significant outcome of the analysis was the identification of a positive correlation between the lidar ratio and cloud centroid temperature, indicating a gradual decrease in the lidar ratio as temperatures dropped. The study established a fundamental consistency in their macroscopic properties, including cloud base height, cloud top height, cloud thickness, cloud centroid height, and cloud centroid temperature. The results for ground-based lidar (CALIOP) are: 10.0 ± 2.1 km (10.0 ± 2.2 km), 11.8 ± 2.1 km (11.5 ± 2.3 km), 1.87 ± 0.83 km (1.52 ± 0.71 km), and 10.5 ± 2.2 km, -46.9 ± 9.7°C (-47.1 ± 10.0°C).These properties exhibited seasonal variations, with cirrus clouds reaching higher altitudes in summer and lower in winter, influenced by the height of the tropopause. The optical properties of cirrus clouds were also analyzed, showing an annual average optical depth of 0.31 ± 0.35 for ground-based lidar and 0.32 ± 0.44 for CALIOP. The study highlighted the distribution of subvisible, thin, and thick cirrus clouds, with a notable prevalence of subvisible clouds during summer, suggesting their frequent formation above 14 km. Furthermore, the study observed linear growth in geometric thickness and optical depth up to 2.5 km from CALIOP and 2.9 km from ground-based lidar. Maximum optical depth was observed at cloud centroid temperatures of -35°C for CALIOP and -40°C for ground-based lidar, with optical depth decreasing as temperatures fell. This suggests that fully glaciated cirrus clouds exhibit the highest optical depth at warmer temperatures, within the complete glaciation temperature range of -35°C to -40°C.

Impact of Aerosols on the Macrophysical and Microphysical Characteristics of Ice-Phase and Mixed-Phase Clouds over the Tibetan Plateau

Using CloudSat/CALIPSO satellite data and ERA5 reanalysis data from 2006 to 2010, the effects of aerosols on ice- and mixed-phase, single-layer, non-precipitating clouds over the Tibetan Plateau during nighttime in the MAM (March to May), JJA (June to August), SON (September to November), and DJF (December to February) seasons were examined. The results indicated the following: (1) The macrophysical and microphysical characteristics of ice- and mixed-phase clouds exhibit a nonlinear trend with increasing aerosol optical depth (AOD). When the logarithm of AOD (lnAOD) was ≤−4.0, with increasing AOD during MAM and JJA nights, the cloud thickness and ice particle effective radius of ice-phase clouds and mixed-phase clouds, the ice water path and ice particle number concentration of ice-phase clouds, and the liquid water path and cloud fraction of mixed-phase clouds all decreased; during SON and DJF nights, the cloud thickness of ice-phase clouds, cloud top height, liquid droplet number concentration, and liquid water path of mixed-phase clouds all decreased. When the lnAOD was &gt;−4.0, with increasing AOD during MAM and JJA nights, the cloud top height, cloud base height, cloud fraction, and ice particle number concentration of ice-phase clouds, and the ice water path of mixed-phase clouds all increased; during SON and DJF nights, the cloud fraction of mixed-phase clouds and the ice water path of ice-phase clouds all increased. (2) Under the condition of excluding meteorological factors, including the U-component of wind, V-component of wind, pressure vertical velocity, temperature, and relative humidity at the atmospheric pressure heights near the average cloud top height, within the cloud, and the average cloud base height, as well as precipitable water vapor, convective available potential energy, and surface pressure. During MAM and JJA nights. When the lnAOD was ≤−4.0, an increase in aerosols may have led to a decrease in the thickness of ice and mixed-phase cloud layers, as well as a reduction in cloud water path values. In contrast, when the lnAOD was &gt;−4.0, an increase in aerosols may contribute to elevated cloud base and cloud top heights for ice-phase clouds. During SON and DJF nights, changes in various cloud characteristics may be influenced by both aerosols and meteorological factors.

Impact of Urbanization on Cloud Characteristics over Sofia, Bulgaria

Urban artificial surfaces and structures induce modifications in land–atmosphere interactions, affecting the exchange of energy, momentum, and substances. These modifications stimulate urban climate formation by altering the values and dynamics of atmospheric parameters, including cloud-related features. This study evaluates the presence and quantifies the extent of such changes over Sofia, Bulgaria. The findings reveal that estimations of low-level cloud base height (CBH) derived from lifting condensation level (LCL) calculations may produce unexpected outcomes due to microclimate influence. Ceilometer data indicate that the CBH of low-level clouds over urban areas exceeds that of surrounding regions by approximately 200 m during warm months and afternoon hours. Moreover, urban clouds exhibit reduced persistence relative to rural counterparts, particularly pronounced in May, June, and July afternoons. Reanalysis-derived low-level cloud cover (LCC) shows no significant disparities between urban and rural areas, although increased LCC is observed above the western and northern city boundaries. Satellite-derived cloud products reveal that the optically thinnest low-level clouds over urban areas exhibit slightly higher cloud tops, but the optically thickest clouds are more prevalent during warm months. These findings suggest an influence of urbanization on cloudiness, albeit nuanced and potentially influenced by the city size and surrounding physical and geographical features.

Variability of Index of Coalescence Activity (ICA) over a rain-shadow region during monsoon and its role in cloud seeding programs in India

The Index of Coalescence Activity (ICA) is an indicator of the microphysical state of the cloud. ICA is related to cloud base temperature and buoyancy of cloud at 500 hPa level. It is a useful parameter for making a decision on the type of seeding to be carried out. During Cloud Aerosol Interaction and Precipitation Enhancement EXperiment (CAIPEEX) over the Indian region, aircraft operations for in situ measurements of aerosols and clouds were conducted from the base station Hyderabad to cover parts of north peninsular India falling under the rain-shadow region. Special Radiosonde flights were conducted in the afternoon hours prior to aircraft operations from the locations Hyderabad and Mahabubnagar in addition to regular India Meteorological Department (IMD) Radiosonde flights at 00 UTC from Hyderabad. The parameters like cloud base height and buoyancy at 500 hPa level obtained from Radiosonde flights have been used to estimate ICA values during the period of the experiment (June to October 2009–2011). The variability of ICA has been studied using morning and afternoon sounding profiles over the rain-shadow region. Morning hours showed highly negative ICA values, as compared to noon hours. During noon time, positive ICA values are observed before monsoon onset (in 2009, monsoon onset was 21st June over the study area) and on few occasions during post-monsoon (October). Positive ICA values represent that collision-coalescence (CC) process is not fully developed and hence, hygroscopic seeding is the proper decision. While, negative ICA values indicate active CC process, suggesting glaciogenic seeding. The soundings conducted just before the cloud seeding operations provide information about the actual thermodynamic and microphysical state of the cloud. Hence, ICA information is useful for providing guidelines to the operational cloud seeding programs.

CALOTRITON: a convective boundary layer height estimation algorithm from ultra-high-frequency (UHF) wind profiler data

Abstract. Long time series of observations of atmospheric dynamics and composition are collected at the French Pyrenean Platform for Observation of the Atmosphere (P2OA). Planetary boundary layer depth is a key variable of the climate system, but it remains difficult to estimate and analyse statistically. In order to obtain reliable estimates of the convective boundary layer height (Zi) and to allow long-term series analyses, a new restitution algorithm, named CALOTRITON, has been developed. It is based on the observations of an ultra-high-frequency (UHF) radar wind profiler (RWP) from P2OA with the help of other instruments for evaluation. Estimates of Zi are based on the principle that the top of the convective boundary layer is associated with both a marked inversion and a decrease in turbulence. Those two criteria are respectively manifested by larger RWP reflectivity and smaller vertical-velocity Doppler spectral width. With this in mind, we introduce a new UHF-deduced dimensionless parameter which weighs the air refractive index structure coefficient with the inverse of vertical velocity standard deviation to the power of x. We then search for the most appropriate local maxima of this parameter for Zi estimates with defined criteria and constraints such as temporal continuity. Given that Zi should correspond to fair-weather cloud base height, we use ceilometer data to optimize our choice of the power x and find that x=3 provides the best comparisons. The estimates of Zi by CALOTRITON are evaluated using different Zi estimates deduced from radiosounding according to different definitions. The comparison shows excellent results with a regression coefficient of up to 0.96 and a root-mean-square error of 71 m, which is close to the vertical resolution of the UHF RWP of 75 m, when conditions are optimal. In more complex situations, that is when the atmospheric vertical structure is itself particularly ambiguous, secondary retrievals allow us to identify potential thermal internal boundary layers or residual layers and help to qualify the Zi estimations. Frequent estimate errors are observed nevertheless; for example, when Zi is below the UHF RWP first reliable gate or when the boundary layer begins its transition to a stable nocturnal boundary layer.

Characteristics and mechanisms of near-surface negative atmospheric electric field anomalies preceding the 5 September 2022, Ms 6.8 Luding earthquake in China

Abstract. A magnitude 6.8 strike-slip earthquake (EQ) struck Luding, Sichuan Province, China, on 5 September 2022, resulting in significant damage to nearby Ganzi Prefecture and the city of Ya'an. In this research, the near-surface atmospheric electric field (AEF) recorded at four sites 15 d before the Luding EQ was analyzed and differentiated, and multisource auxiliary data including precipitation, cloud base height, and low cloud cover were used at the same time. Nine possible seismic AEF anomalies at four sites were obtained preliminarily. Accordingly, microwave brightness temperature (MBT) data, which are very sensitive to the surface dielectrics and are closely related to the air ionization, together with surface soil moisture, lithology, and a 3D-simulated crustal stress field, were jointly analyzed to confirm the seismic relations of the obtained negative AEF anomalies. The geophysical environment for crustal high-stress concentration, positive charge carrier transfer, and surface accumulation was demonstrated to exist and to meet the conditions necessary to generate local negative AEF anomalies. Furthermore, to deal with the spatial disparities in sites and regions with potential atmospheric ionization, near-surface wind field data were employed to scrutinize the reliability of the AEF anomalies by comprehensively analyzing the spatial relationships among surface charges accumulation areas, wind direction and speed, and the AEF sites. Finally, four negative AEF anomalies were deemed to be closely related to the Luding EQ, and the remaining five possible anomalies were ruled out. A possible mechanism of negative AEF anomalies before the Luding EQ is proposed: positive charge carriers were generated from the underground high-stress concentration areas and then transferred to and accumulated on the ground surface to ionize the surface air, thus disturbing the AEF above the ground. This study presents a method for identifying and analyzing seismic AEF anomalies and is also beneficial for the examination of the pre-earthquake coupling process between the coversphere and the atmosphere.

Geometrical and optical properties of cirrus clouds in Barcelona, Spain: analysis with the two-way transmittance method of 4 years of lidar measurements

Abstract. In this paper a statistical study of cirrus geometrical and optical properties based on 4 years of continuous ground-based lidar measurements with the Barcelona (Spain) Micro Pulse Lidar (MPL) is analysed. First, a review of the literature on the two-way transmittance method is presented. This method is a well-known lidar inversion method used to retrieve the optical properties of an aerosol–cloud layer between two molecular (i.e. aerosol and cloud-free) regions below and above, without the need to make any a priori assumptions about their optical and/or microphysical properties. Second, a simple mathematical expression of the two-way transmittance method is proposed for both ground-based and spaceborne lidar systems. This approach of the method allows the retrieval of the cloud optical depth, the cloud column lidar ratio and the vertical profile of the cloud backscatter coefficient. The method is illustrated for a cirrus cloud using measurements from the ground-based MPL and from the spaceborne Cloud-Aerosol Lidar with Orthogonal Polarization (CALIOP). Third, the database is then filtered with a cirrus identification criterion based on (and compared to) the literature using only lidar and radiosonde data. During the period from November 2018 to September 2022, 367 high-altitude cirrus clouds were identified at 00:00 and 12:00 UTC, of which 203 were successfully inverted with the two-way transmittance method. The statistical results of these 203 high-altitude cirrus clouds show that the cloud thickness is 1.8 ± 1.1 km, the mid-cloud temperature is −51 ± 8 ∘C and the linear cloud depolarization ratio is 0.32 ± 0.13. The application of the transmittance method yields an average cloud optical depth (COD) of 0.36 ± 0.45 and a mean effective column lidar ratio of 30 ± 19 sr. Statistical results of the errors associated with the two-way transmittance method retrievals are also provided. The highest occurrence of cirrus is observed in spring and the majority of cirrus clouds (48 %) are visible (0.03 &lt; COD &lt; 0.3), followed by opaque (COD &gt; 0.3) with a percentage of 38 %. Together with results from other sites, possible latitudinal dependencies have been analysed together with correlations between cirrus cloud properties. For example, we noted that in Barcelona the COD correlates positively with the cloud base temperature, effective column lidar ratio and linear cloud depolarization ratio and negatively with the cloud base height.