A Review of Techniques for Diagnosing the Atmospheric Boundary Layer Height (ABLH) Using Aerosol Lidar Data

Retrieval of the planetary boundary layer height from lidar measurements by a deep-learning method based on the wavelet covariance transform

Understanding and characterization of the planetary boundary layer (PBL) are of great importance in terms of air pollution management, weather forecasting, modelling of climate change, etc. Although many lidar-based approaches have been proposed for the retrieval of the PBL height (PBLH) in case studies, development of a robust lidar-based algorithm without human intervention is still of great challenging. In this work, we have demonstrated a novel deep-learning method based on the wavelet covariance transform (WCT) for the PBLH evaluation from atmospheric lidar measurements. Lidar profiles are evaluated according to the WCT with a series of dilation values from 200 m to 505 m to generate 2-dimensional wavelet images. A large number of wavelet images and the corresponding PBLH-labelled images are created as the training set for a convolutional neural network (CNN), which is implemented based on a modified VGG16 (VGG - Visual Geometry Group) convolutional neural network. Wavelet images obtained from lidar profiles have also been prepared as the test set to investigate the performance of the CNN. The PBLH is finally retrieved by evaluating the predicted PBLH-labelled image and the wavelet coefficients. Comparison studies with radiosonde data and the Micro-Pulse-Lidar Network (MPLNET) PBLH product have successfully validated the promising performance of the deep-learning method for the PBLH retrieval in practical atmospheric sensing.

Simultaneous Extraction of Planetary Boundary-Layer Height and Aerosol Optical Properties from Coherent Doppler Wind Lidar

Planetary boundary-layer height is an important physical quantity for weather forecasting models and atmosphere environment assessment. A method of simultaneously extracting the surface-layer height (SLH), mixed-layer height (MLH), and aerosol optical properties, which include aerosol extinction coefficient (AEC) and aerosol optical depth (AOD), based on the signal-to-noise ratio (SNR) of the same coherent Doppler wind lidar (CDWL) is proposed. The method employs wavelet covariance transform to locate the SLH and MLH using the local maximum positions and an automatic algorithm of dilation operation. AEC and AOD are determined by the fitting curve using the SNR equation. Furthermore, the method demonstrates the influential mechanism of optical properties on the SLH and MLH. MLH is linearly correlated with AEC and AOD because of solar heating increasing. The results were verified by the data of an ocean island site in China.

Characterizing warm atmospheric boundary layer over land by combining Raman and Doppler lidar measurements

PBL plays a critical role in the atmosphere by transferring heat, moisture, and momentum. The warm PBL has a distinct diurnal cycle including daytime convective mixing layer (ML) and nighttime residual layer developments. Thus, for PBL characterization and process study, simultaneous determinations of PBL height (PBLH) and ML height (MLH) are necessary. Here, new approaches are developed to provide reliable PBLH and MLH to characterize warm PBL evolution. The approaches use Raman lidar (RL) water vapor mixing ratio (WVMR) and Doppler lidar (DL) vertical velocity measurements at the Southern Great Plains (SGP) atmospheric observatory, which was established by the Atmospheric Radiation Measurement (ARM) user facility. Compared with widely used lidar aerosol measurements for PBLH, WVMR is a better trace for PBL vertical mixing. For PBLH, the approach classifies PBL water vapor structures into a few general patterns, then uses a slope method and dynamic threshold method to determine PBLH. For MLH, wavelet analysis is used to re-construct 2-D variance from DL vertical wind velocity measurements according to the turbulence eddy size to minimize the impacts of gravity wave and eddy size on variance calculations; then, a dynamic threshold method is used to determine MLH. Remotely-sensed PBLHs and MLHs are compared with radiosonde measurements based on the Richardson number method. Good agreements between them confirm that the proposed new algorithms are reliable for PBLH and MLH characterization. The algorithms are applied to warm seasons' RL and ML measurements at the SGP site for five years to study warm season PBL structure and processes. The weekly composited diurnal evolutions of PBLHs and MLHs in warm climate were provided to illustrate diurnal and seasonal PBL evolutions. This reliable PBLH and MLH dataset will be valuable for PBL process study, model evolution, and PBL parameterization improvement.

Micrometeorological Methods for the Indirect Estimation of Odorous Emissions

Odors are typically released into the atmosphere as diffuse emissions from area and volume sources, whose detailed quantification in terms of odor emission rate is often hardly achievable by direct source sampling. Indirect methods, involving the use of micrometeorological methods in order to correlate downwind concentrations to the emission rates, are already mentioned in literature, but rarely found in real applications for the quantification of odor emissions. The instrumentation needed for the development of micrometeorological methods has nowadays become accessible in terms of prices and reliability, thus making the implementation of such methods to industrial applications more and more interesting. For this reason, this work aims to provide an overview of micrometeorological methods and investigate their effective applicability to odors, thereby providing a short description of the physics related to such methods and analyzing the relevant scientific literature. The theoretical basis of these methods is presented, and their advantages and disadvantages are discussed. Moreover, their applicability to the estimation of odor emissions is discussed by providing some suggestions about the suitable ways to evaluate the most critical parameters needed for the calculation of the odor emission rate.

Estimating Boundary Layer Height from LiDAR Data under Complex Atmospheric Conditions Using Machine Learning

Reliable estimation of the atmospheric boundary layer height (ABLH) is critical for a range of meteorological applications, including air quality assessment and weather forecasting. Several algorithms have been proposed to detect ABLH from aerosol LiDAR backscatter data. However, most of these focus on cloud-free conditions or use other ancillary instruments due to strong interference from clouds or residual layer aerosols. In this paper, a machine learning method named the Mahalanobis transform K-near-means (MKnm) algorithm is first proposed to derive ABLH under complex atmospheric conditions using only LiDAR-based instruments. It was applied to the micro pulse LiDAR data obtained at the Southern Great Plains site of the Atmospheric Radiation Measurement (ARM) program. The diurnal cycles of ABLH from cloudy weather were detected by using the gradient method (GM), wavelet covariance transform method (WM), K-means, and MKnm. Meanwhile, the ABLH obtained by these four methods under cloud or residual layer conditions based on micropulse LiDAR data were compared with the reference height retrieved from radiosonde data. The results show that MKnm was good at tracking the diurnal variation of ABLH, and the ABLHs obtained by it have remarkable correlation coefficients and smaller mean absolute error and mean deviation with the radiosonde-derived ABLHs than those measured by other three methods. We conclude that MKnm is a promising algorithm to estimate ABLH under cloud or residual layer conditions.

Analyzing Atmospheric Circulation Patterns Using Mass Fluxes Calculated from Weather Balloon Measurements: North Atlantic Region as a Case Study

In recent decades, efforts to investigate atmospheric circulation patterns have predominantly relied on either semi-empirical datasets (i.e., reanalyses) or modeled output (i.e., global climate models, GCMs). While both approaches can provide important insights, there is a need for more empirical data to supplement these approaches. In this paper, we demonstrate how the application of relatively simple calculations to the basic measurements from a standard weather balloon radiosonde can provide a vertical profile of the horizontal atmospheric mass fluxes. These mass fluxes can be resolved into their meridional (north/south) and zonal (east/west) components. This provides a new useful empirical tool for analyzing atmospheric circulations. As a case study, we analyze the results for a selected five stations along a fairly constant meridian in the North Atlantic sector from 2015–2019. For each station, we find the atmospheric mass flux profiles from the lower troposphere to mid-stratosphere are surprisingly coherent, suggesting stronger interconnection between the troposphere and stratosphere than previously thought. Although our five stations span a region nominally covered by the classical polar, Ferrel and Hadley meridional circulation cells, the results are inconsistent with those expected for polar and Ferrel cells and only partially consistent with that of a Hadley cell. However, the region is marked by very strong prevailing westerly (west to east) mass fluxes for most of the atmosphere except for the equatorial surface easterlies (“trade winds”). We suggest that the extension of the techniques of this case study to other stations and time periods could improve our understanding of atmospheric circulation patterns and their time variations.

Edge Detection Method for Determining Boundary Layer Height Based on Doppler Lidar

The top of the boundary layer, referred to as the planetary boundary layer height (BLH), is an important physical parameter in atmospheric numerical models, which has a critical role in atmospheric simulation, air pollution prevention, and climate prediction. The traditional methods for determining BLHs using Doppler lidar vertical velocity variance (σw2) can be classified into the variance and peak methods, which depend on atmospheric conditions due to their use of a single threshold, hence limiting their ability to estimate diurnal BLHs. Edge detection (ED) was later introduced in BLH estimation due to its ability to identify the 2D gradient of an image. A key step in ED is automatically identifying the edge of BLHs based on the peaks of the profile, hence avoiding the influence of extreme atmospheric conditions. Two cases in the diurnal cycle on 4 March 2019 and 8 July 2019 reveal that ED outperforms both the variance and peak methods in nighttime and extreme atmospheric conditions. The retrieved BLHs from 2018 to 2020 were compared with radiosonde (RS) measurements for the same time at the neutral, stable, and convective boundary layers. The correlation coefficient (R: 0.4 vs. 0.05, 0.14; 0.26 vs. −0.10, −0.16; 0.35 vs. 0.01, 0.16) and root mean square error (RMSE (km): 0.58 vs. 0.82, 0.90; 0.37 vs. 1.01, 0.50; 0.66 vs. 0.98, 0.82) obtained by the ED method were higher and lower than those obtained by the variance and peak methods, respectively. The mean absolute error (MAE) of the ED method under the NBL, SBL, and CBL conditions are lower than the variance and peak methods (MAE (km): 0.44, 0.14, 0.50 vs. 0.62, 0.34, 0.64; 0.59, 0.75, 0.74), respectively. The mean relative error (MRE) of the ED method is lower than the variance and peak methods under the NBL condition (MRE: −8.88% vs. −18.39%, 13.91%). Under the SBL, the MRE of the ED method is lower than the variance method and higher than the peak method (−38.64%, vs. −152.23%; 14.02%). Under the CBL, the MRE of the ED method is lower than the variance method and higher than the peak method (−15.07% vs. 2.24%; 5.64%). In addition, the comparison between ED and wavelet covariance transform (WCT) method and RS measurements showed that the ED method has a similar performance with the WCT method and is even better. In the long-term analysis, the hourly and monthly BLHs in the diurnal and annual cycles, respectively, as obtained by ED, were highly consistent with the RS measurements and obtained the lowest standard error. In the annual cycle, the retrieved BLHs in summer and autumn were higher than those retrieved in spring and winter.

Atmospheric Boundary Layer Wind Profile Estimation Using Neural Networks Applied to Lidar Measurements

This paper introduces a new methodology for estimating the wind profile within the ABL (Atmospheric Boundary Layer) using a neural network and a single-point near-ground measurement. An important advantage of this solution when compared with others available in the literature is that it only requires near surface measurements for the prognosis once the neural network is trained. Another advantage is that it can be used to study the wind profile temporal evolution. This work uses data collected by a lidar sensor located at the Universidad de León (Spain). The neural network best configuration was determined using sensibility analyses. The result is a multilayer perceptron with three layers for each altitude: the input layer has six nodes for the last three measurements, the second has 128 nodes and the third consists of two nodes that provide u and v. The proposed method has better performance than traditional methods. The obtained wind profile information obtained is useful for multiple applications, such as preliminary calculations of the wind resource or CFD models.

A Comparison of Wintertime Atmospheric Boundary Layer Heights Determined by Tethered Balloon Soundings and Lidar at the Site of SACOL

High-precision and -resolution atmospheric boundary layer height (BLH) has received increasing attention in air pollution research in recent years. The low time resolution of sounding data is the main challenge to validate BLH retrieval from lidar observations. To resolve this issue, we conducted simultaneous tethered balloon sounding and lidar observations at the Semi-Arid Climate and Environment Observatory of Lanzhou University (SACOL) during winter 2019–2020. The BLHs derived from the tethered balloon sounding data were 170, 210, 393, 676, 423, and 190 m at 02:00, 08:00, 11:00, 14:00, 17:00, and 20:00 (Beijing time), respectively. The diurnal evolution of BLH was reasonably captured by lidar observation-based wavelet covariance transform and ideal profile fitting methods, which exhibited correlation coefficients of 0.91 and 0.89, respectively, with the BLHs determined from tethered balloon sounding data. The lidar results slightly overestimated the BLHs, though all results were acceptable when considering both the absolute and relative errors with respect to BLHs from the tethered balloon data. Our results revealed high-precision and -resolution diurnal variations in BLH at SACOL in Northwest China and suggest the importance of validating lidar-based BLHs using simultaneous sounding data.

Vertical Profiles of Atmospheric Species Concentrations and Nighttime Boundary Layer Structure in the Dry Season over an Urban Environment in Central Amazon Collected by an Unmanned Aerial Vehicle

Nighttime vertical profiles of ozone, PM2.5 and PM10 particulate matter, carbon monoxide, temperature, and humidity were collected by a copter-type unmanned aerial vehicle (UAV) over the city of Manaus, Brazil, in central Amazon during the dry season of 2018. The vertical profiles were analyzed to understand the structure of the urban nighttime boundary layer (NBL) and pollution within it. The ozone concentration, temperature, and humidity had an inflection between 225 and 350 m on most nights, representing the top of the urban NBL. The profile of carbon monoxide concentration correlated well with the local evening vehicular congestion of a modern transportation fleet, providing insight into the surface-atmosphere dynamics. In contrast, events of elevated PM2.5 and PM10 concentrations were not explained well by local urban emissions, but rather by back trajectories that intersected regional biomass burning. These results highlight the potential of the emerging technologies of sensor payloads on UAVs to provide new constraints and insights for understanding the pollution dynamics in nighttime boundary layers in urban regions.

Study of the Vertical Structure of the Coastal Boundary Layer Integrating Surface Measurements and Ground-Based Remote Sensing

The understanding of the atmospheric processes in coastal areas requires the availability of quality datasets describing the vertical and horizontal spatial structure of the Atmospheric Boundary Layer (ABL) on either side of the coastline. High-resolution Numerical Weather Prediction (NWP) models can provide this information and the main ingredients for good simulations are: an accurate description of the coastline and a correct subgrid process parametrization permitting coastline discontinuities to be caught. To provide an as comprehensive as possible dataset on Mediterranean coastal area, an intensive experimental campaign was realized at a near-shore Italian site, using optical and acoustic ground-based remote sensing and surface instruments, under different weather characteristic and stability conditions; the campaign is also fully simulated by a NWP model. Integrating information from instruments responding to different atmospheric properties allowed for an explanation of the development of various patterns in the vertical structure of the atmosphere. Wind LiDAR measurements provided information of the internal boundary layer from the value of maximum height reached by the wind profile; a height between 80 and 130 m is often detected as an interface between two different layers. The NWP model was able to simulate the vertical wind profiles and the eight of the ABL.

Tailored Algorithms for the Detection of the Atmospheric Boundary Layer Height from Common Automatic Lidars and Ceilometers (ALC)

A detailed understanding of atmospheric boundary layer (ABL) processes is key to improve forecasting of pollution dispersion and cloud dynamics in the context of future climate scenarios. International networks of automatic lidars and ceilometers (ALC) are gathering valuable data that allow for the height of the ABL and its sublayers to be derived in near real time. A new generation of advanced methods to automatically detect the ABL heights now exist. However, diversity in ALC models means these algorithms need to be tailored to instrument-specific capabilities. Here, the advanced algorithm STRATfinder is presented for application to high signal-to-noise ratio (SNR) ALC observations, and results are compared to an automatic algorithm designed for low-SNR measurements (CABAM). The two algorithms are evaluated for application in an operational network setting. Results indicate that the ABL heights derived from low-SNR ALC have increased uncertainty during daytime deep convection, while high-SNR observations can have slightly reduced capabilities in detecting shallow nocturnal layers. Agreement between the ALC-based methods is similar when either is compared to the ABL heights derived from temperature profile data. The two independent methods describe very similar average diurnal and seasonal variations. Hence, high-quality products of ABL heights may soon become possible at national and continental scales.

An automated common algorithm for planetary boundary layer retrievals using aerosol lidars in support of the U.S. EPA Photochemical Assessment Monitoring Stations Program

A unique automated planetary boundary layer (PBL) retrieval algorithm is proposed as a common cross-platform method for use with commercially available ceilometers for implementation under the redesigned U.S. Environmental Protection Agency Photochemical Assessment Monitoring Stations program. This algorithm addresses instrument signal quality and screens for precipitation and cloud layers before the implementation of the retrieval methodology using the Haar wavelet covariance transform method. Layer attribution for the PBL height is supported with the use of continuation and time-tracking parameters, and uncertainties are calculated for individual PBL height retrievals. Commercial ceilometer retrievals are tested against radiosonde PBL height and cloud-base height during morning and late afternoon transition times, critical to air quality model prediction and when retrieval algorithms struggle to identify PBL heights. A total of 58 radiosonde profiles were used and retrievals for nocturnal stable layers, residual layers and mixing layers were assessed. Overall good agreement was found for all comparisons with one system showing limitations for the cases of nighttime surface stable layers and daytime mixing layer. It is recommended that nighttime shallow stable layer retrievals be performed with a recommended minimum height or with additional verification. Retrievals of residual layer heights and mixing layer comparisons revealed overall good correlations to radiosonde heights (correlation coefficients, r², ranging from 0.89 - 0.96 and bias ranging from ~ -131 to +63 m, and r² from 0.88 - 0.97 and bias from -119 to +101 m, respectively).

Determination of Planetary Boundary Layer height with Lidar Signals Using Maximum Limited Height Initialization and Range Restriction (MLHI-RR)

The planetary boundary layer height (PBLH) is a vital parameter to characterize the surface convection, which determines the diffusion of air pollutants. The accurate inversion of PBLH is extremely important for the study of aerosol concentrations, in order to predict air quality and provide weather forecast. Aerosol lidar, a powerful remote sensing instrument for detecting the characteristics of atmospheric temporal and spatial evolution, can continuously retrieve the planetary boundary layer (PBL) and obtain high resolution measurements. However, multi-layer conditions, including one or more layers of aerosol, or cloud above the PBL, can seriously interfere the accuracy of PBLH determined by lidar. A new technique of maximum limited height initialization and range restriction (MLHI-RR) is proposed to eliminate the impact of multi-layer conditions on PBLH determination. Four widely used methods for deriving PBLH are utilized, in addition to the MLHI-RR constraint. Comparisons demonstrate that the proposed technique can determine the PBLH in multi-layer conditions with higher accuracy. The proposed technique requires no affiliate information besides lidar signals, which provide a convenient method for PBLH determination under complicated conditions.

Variability of the Boundary Layer Over an Urban Continental Site Based on 10 Years of Active Remote Sensing Observations in Warsaw

Atmospheric boundary layer height (ABLH) was observed by the CHM15k ceilometer (January 2008 to October 2013) and the PollyXT lidar (July 2013 to December 2018) over the European Aerosol Research LIdar NETwork to Establish an Aerosol Climatology (EARLINET) site at the Remote Sensing Laboratory (RS-Lab) in Warsaw, Poland. Out of a maximum number of 4017 observational days within this period, a subset of quasi-continuous measurements conducted with these instruments at the same wavelength (1064 nm) was carefully chosen. This provided a data sample of 1841 diurnal cycle ABLH observations. The ABLHs were derived from ceilometer and lidar signals using the wavelet covariance transform method (WCT), gradient method (GDT), and standard deviation method (STD). For comparisons, the rawinsondes of the World Meteorological Organization (WMO 12374 site in Legionowo, 25 km distance to the RS-Lab) were used. The ABLHs derived from rawinsondes by the skew-T-log-p method and the bulk Richardson (bulk-Ri) method had a linear correlation coefficient (R2) of 0.9 and standard deviation (SD) of 0.32 km. A comparison of the ABLHs obtained for different methods and instruments indicated a relatively good agreement. The ABLHs estimated from the rawinsondes with the bulk-Ri method had the highest correlations, R2 of 0.80 and 0.70 with the ABLHs determined using the WCT method on ceilometer and lidar signals, respectively. The three methods applied to the simultaneous, collocated lidar, and ceilometer observations (July to October 2013) showed good agreement, especially for the WCT method (R2 of 0.94, SD of 0.19 km). A scaling threshold-based algorithm was proposed to homogenize ceilometer and lidar datasets, which were applied on the lidar data, and significantly improved the coherence of the results (R2 of 0.98, SD of 0.11 km). The difference of ABLH between clear-sky and cloudy conditions was on average below 230 m for the ceilometer and below 70 m for the lidar retrievals. The statistical analysis of the long-term observations indicated that the monthly mean ABLHs varied throughout the year between 0.6 and 1.8 km. The seasonal mean ABLH was of 1.16 ± 0.16 km in spring, 1.34 ± 0.15 km in summer, 0.99 ± 0.11 km in autumn, and 0.73 ± 0.08 km in winter. In spring and summer, the daytime and nighttime ABLHs appeared mainly in a frequency distribution range of 0.6 to 1.0 km. In winter, the distribution was common between 0.2 and 0.6 km. In autumn, it was relatively balanced between 0.2 and 1.2 km. The annual mean ABLHs maintained between 0.77 and 1.16 km, whereby the mean heights of the well-mixed, residual, and nocturnal layer were 1.14 ± 0.11, 1.27 ± 0.09, and 0.71 ± 0.06 km, respectively (for clear-sky conditions). For the whole observation period, the ABLHs below 1 km constituted more than 60% of the retrievals. A strong seasonal change of the monthly mean ABLH diurnal cycle was evident; a mild weakly defined autumn diurnal cycle, followed by a somewhat flat winter diurnal cycle, then a sharp transition to a spring diurnal cycle, and a high bell-like summer diurnal cycle. A prolonged summertime was manifested by the September cycle being more similar to the summer than autumn cycles.

Vertical Profiles of Ozone Concentration Collected by an Unmanned Aerial Vehicle and the Mixing of the Nighttime Boundary Layer over an Amazonian Urban Area

The nighttime boundary layer was studied in an urban area surrounded by tropical forest by use of a copter-type unmanned aerial vehicle (UAV) in central Amazonia during the wet season. Fifty-seven vertical profiles of ozone concentration, potential temperature, and specific humidity were collected from surface to 500 m above ground level (a.g.l.) at high vertical and temporal resolutions by use of embedded sensors on the UAV. Abrupt changes in ozone concentration with altitude served as a proxy of nighttime boundary layer (NBL) height for the case of a normal, undisturbed, stratified nighttime atmosphere, corresponding to 40% of the cases. The median height of the boundary layer was 300 m. A turbulent mixing NBL constituted 28% of the profiles, while the median height of the boundary layer was 290 m. The remaining 32% of profiles corresponded to complex atmospheres without clear boundary layer heights. The occurrence of the three different cases correlated well with relative cloud cover. The results show that the standard nighttime model widely implemented in chemical transport models holds just 40% of the time, suggesting new challenges in modeling of regional nighttime chemistry. The boundary layer heights were also somewhat higher than observed previously over forested and pasture areas in Amazonia, indicating the important effect of the urban heat island.