The extensive use of fireworks generates large amounts of pollutants, deteriorating air quality and potentially causing adverse health impacts. In Medellín and its metropolitan area, although fireworks are banned during December, their use is widespread during the Christmas season, particularly during the midnight of November 30 (La Alborada) and New Year’s Eve (NYE). It is therefore essential to assess the effects of these celebrations on air quality in the region. Air-quality data from the official monitoring network and a low-cost particulate matter (PM) citizen science project, backscattering intensity (BI) retrievals from a ceilometer network, potential temperature from a microwave radiometer, and information from a radar wind profiler provide an excellent platform to study the spatio-temporal distribution of contaminants resulting from the La Alborada and NYE celebrations. Substantial increases in PM2.5 and PM10 mass concentrations due to La Alborada and NYE, ranging in some cases from 50 to 100 μgm⁻³, are observed in the Aburrá Valley and particularly in the densely populated communes of Medellín, with most concentration changes corresponding to ultrafine and fine particles. The PM increments resulting from fireworks show almost no increase in the net amount of black carbon in the atmosphere. Ceilometer BI profiles show a substantial change immediately after the La Alborada and NYE midnights, confined to the atmospheric boundary layer (ABL). Strong thermal inversions lead to fairly homogeneous increments in BI within the ABL, lasting until the onset of the convective boundary layer. In contrast, weak thermal inversions lead to rapid dispersion of aerosols, allowing them to episodically escape above the ABL.

Aerosol-cloud interactions (ACIs) are vital for regulating Earth’s climate by influencing energy and water cycles. Yet, effects of ACI bear large uncertainties, evidenced by systematic discrepancies between observed and modeled estimates. This study quantifies a major bias in ACI determinations, stemming from conventional surface or space measurements that fail to capture aerosol at the cloud level unless the cloud is coupled with land surface. We introduce an advanced approach to determine radiative forcing of ACI by accounting for cloud-surface coupling. By integrating field observations, satellite data, and model simulations, this approach reveals a drastic alteration in aerosol vertical transport and ACI effects caused by cloud coupling. In coupled regimes, aerosols enhance cloud droplet number concentration across the boundary layer more homogeneously than in decoupled conditions, under which aerosols from the free atmosphere predominantly affect cloud properties, leading to marked cooling effects. Our findings spotlight cloud-surface coupling as a key factor for ACI quantification, hinting at potential underassessments in traditional estimates.

Knowledge of the planetary boundary layer height (PBLH) is crucial for various applications in atmospheric and environmental sciences. Lidar measurements are frequently used to monitor the evolution of the PBLH, providing more frequent observations than traditional radiosonde‐based methods. However, lidar‐derived PBLH estimates have substantial uncertainties, contingent upon the retrieval algorithm used. In addressing this, we applied the Different Thermo‐Dynamic Stabilities (DTDS) algorithm to establish a PBLH data set at five separate Department of Energy's Atmospheric Radiation Measurement sites across the globe. Both the PBLH methodology and the products are subject to rigorous assessments in terms of their uncertainties and constraints, juxtaposing them with other products. The DTDS‐derived product consistently aligns with radiosonde PBLH estimates, with correlation coefficients exceeding 0.77 across all sites. This study delves into a detailed examination of the strengths and limitations of PBLH data sets with respect to both radiosonde‐derived and other lidar‐based estimates of the PBLH by exploring their respective errors and uncertainties. It is found that varying techniques and definitions can lead to diverse PBLH retrievals due to the inherent intricacy and variability of the boundary layer. Our DTDS‐derived PBLH data set outperforms existing products derived from ceilometer data, offering a more precise representation of the PBLH. This extensive data set paves the way for advanced studies and an improved understanding of boundary‐layer dynamics, with valuable applications in weather forecasting, climate modeling, and environmental studies.

The Planetary Boundary Layer Height (PBLH) significantly impacts weather, climate, and air quality. Understanding the global diurnal variation of the PBLH is particularly challenging due to the necessity of extensive observations and suitable retrieval algorithms that can adapt to diverse thermodynamic and dynamic conditions. This study utilized data from the Cloud-Aerosol Transport System (CATS) to analyze the diurnal variation of PBLH in both continental and marine regions. By leveraging CATS data and a modified version of the Different Thermo-Dynamics Stability (DTDS) algorithm, along with machine learning denoising, the study determined the diurnal variation of the PBLH in continental mid-latitude and marine regions. The CATS DTDS-PBLH closely matches ground-based lidar and radiosonde measurements at the continental sites, with correlation coefficients above 0.6 and well-aligned diurnal variability, although slightly overestimated at nighttime. In contrast, PBLH at the marine site was consistently overestimated due to the viewing geometry of CATS and complex cloud structures. The study emphasizes the importance of integrating meteorological data with lidar signals for accurate and robust PBLH estimations, which are essential for effective boundary layer assessment from satellite observations.

Remote sensing measurements have been widely used to estimate the planetary boundary layer height (PBLHT). Each remote sensing approach offers unique strengths and faces different limitations. In this study, we use machine learning (ML) methods to produce a best-estimate PBLHT (PBLHT-BE-ML) by integrating four PBLHT estimates derived from remote sensing measurements at the Department of Energy (DOE) Atmospheric Radiation Measurement (ARM) Southern Great Plains (SGP) observatory. Three ML models – random forest (RF) classifier, RF regressor, and light gradient-boosting machine (LightGBM) – were trained on a dataset from 2017 to 2023 that included radiosonde, various remote sensing PBLHT estimates, and atmospheric meteorological conditions. Evaluations indicated that PBLHT-BE-ML from all three models improved alignment with the PBLHT derived from radiosonde data (PBLHT-SONDE), with LightGBM demonstrating the highest accuracy under both stable and unstable boundary layer conditions. Feature analysis revealed that the most influential input features at the SGP site were the PBLHT estimates derived from (a) potential temperature profiles retrieved using Raman lidar (RL) and atmospheric emitted radiance interferometer (AERI) measurements (PBLHT-THERMO), (b) vertical velocity variance profiles from Doppler lidar (PBLHT-DL), and (c) aerosol backscatter profiles from micropulse lidar (PBLHT-MPL). The trained models were then used to predict PBLHT-BE-ML at a temporal resolution of 10 min, effectively capturing the diurnal evolution of PBLHT and its significant seasonal variations, with the largest diurnal variation observed over summer at the SGP site. We applied these trained models to data from the ARM Eastern Pacific Cloud Aerosol Precipitation Experiment (EPCAPE) field campaign (EPC), where the PBLHT-BE-ML, particularly with the LightGBM model, demonstrated improved accuracy against PBLHT-SONDE. Analyses of model performance at both the SGP and EPC sites suggest that expanding the training dataset to include various surface types, such as ocean and ice-covered areas, could further enhance ML model performance for PBLHT estimation across varied geographic regions.