Convective Boundary Layer Height Research Articles

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.

The Influence of Terrain Smoothing on Simulated Convective Boundary-Layer Depths in Mountainous Terrain

Many applications rely on a correct estimation of the convective boundary layer (CBL) depth over mountainous terrain, but often these applications use numerical model simulations. Although models inevitably smooth terrain, the amount of smoothing depends on grid spacing. We investigate the behavior of the CBL in coarse- and fine-grid models applied to mountainous terrain by using output from an operational mesoscale modeling system and by performing quasi-idealized simulations. We investigate different areas in different climate zones using different CBL top derivation methods, grid spacing ratios, planetary boundary layer (PBL) schemes, and terrain smoothing. We find that when compared to fine-grid simulations, CBL depths are systematically larger in coarse domains over mountaintops, and to a lesser extent in valleys. On average, differences between coarse- and fine-domains over mountaintops could reach around 10%. In certain locations, differences could be as high as 25%. We attribute the result to terrain smoothing. Similarly, when using a coarse-grid CBL height (relative to mean sea level) interpolated using fine-grid terrain information, there is good agreement with fine-grid CBL depths over mountaintops and less agreement in valleys. Our results have implications for applications that use output from coarse model grids in mountainous terrain. These include inverse modeling studies (e.g., greenhouse gas budget estimations or integrated water vapor transport), PBL evaluation studies, climate research, air quality applications, planning and executing prescribed burns, and studies associated with precipitation over mountainous terrain.

Diagnosing Factors Leading to an Incorrect Supercell Thunderstorm Forecast

Abstract On 28 April 2019, hourly forecasts from the operational High-Resolution Rapid Refresh (HRRR) model consistently predicted an isolated supercell storm late in the day near Dodge City, Kansas, that subsequently was not observed. Two convection-allowing model (CAM) ensemble runs are created to explore the reasons for this forecast error and implications for severe weather forecasting. The 40-member CAM ensembles are run using the HRRR configuration of the WRF-ARW Model at 3-km horizontal grid spacing. The Gridpoint Statistical Interpolation (GSI)-based ensemble Kalman filter is used to assimilate observations every 15 min from 1500 to 1900 UTC, with resulting ensemble forecasts run out to 0000 UTC. One ensemble only assimilates conventional observations, and its forecasts strongly resemble the operational HRRR with all ensemble members predicting a supercell storm near Dodge City. In the second ensemble, conventional observations plus observations of WSR-88D radar clear-air radial velocities, WSR-88D diagnosed convective boundary layer height, and GOES-16 all-sky infrared brightness temperatures are assimilated to improve forecasts of the preconvective environment, and its forecasts have half of the members predicting supercells. Results further show that the magnitude of the low-level meridional water vapor flux in the moist tongue largely separates members with and without supercells, with water vapor flux differences of 12% leading to these different outcomes. Additional experiments that assimilate only radar or satellite observations show that both are important to predictions of the meridional water vapor flux. This analysis suggests that mesoscale environmental uncertainty remains a challenge that is difficult to overcome. Significance Statement Forecasts from operational numerical models are the foundation of weather forecasting. There are times when these models make forecasts that do not come true, such as 28 April 2019 when successive forecasts from the operational High-Resolution Rapid Refresh (HRRR) model predicted a supercell storm late in the day near Dodge City, Kansas, that subsequently was not observed. Reasons for this forecast error are explored using numerical experiments. Results suggest that relatively small changes to the prestorm environment led to large differences in the evolution of storms on this day. This result emphasizes the challenges to operational severe weather forecasting and the continued need for improved use of all available observations to better define the atmospheric state given to forecast models.

An ensemble method for improving the estimation of planetary boundary layer height from radiosonde data

Abstract. The planetary boundary layer (PBL) height (PBLH) is an important parameter for weather, climate, and air quality models. Radiosonde is one of the most commonly used instruments for PBLH determination and is generally accepted as a standard for other methods. However, mainstream approaches for the estimation of PBLH from radiosonde present some uncertainties and even show disadvantages under some circumstances, and the results need to be visually verified, especially during the transition period of different PBL regimes. To avoid the limitations of individual methods and provide a benchmark estimation of PBLH, we propose an ensemble method based on high-resolution radiosonde data collected in Beijing in 2017. Seven existing methods including four gradient-based methods are combined along with statistical modification. The ensemble method is verified at 08:00, 14:00, and 20:00 Beijing time (BJT = UTC+8), respectively. The overestimation of PBLH can be effectively eliminated by setting thresholds for gradient-based methods, and the inconsistency between individual methods can be reduced by clustering. Based on the statistics of a 1-year observational analysis, the effectiveness (E) of the ensemble method reaches up to 62.6 %, an increase of 6.5 %–53.0 % compared to the existing methods. Nevertheless, the ensemble method suffers to some extent from uncertainties caused by the consistent overestimation of PBLH, the profiles with a multi-layer structure, and the intermittent turbulence in the stable boundary layer (SBL). Finally, this method has been applied to characterize the diurnal and seasonal variations of different PBL regimes. Particularly, the average convective boundary layer (CBL) height is found to be the highest in summer, and the SBL is lowest in summer with about 200 m. The average PBLH at the transition stage lies around 1100 m except in winter. These findings imply that the ensemble method is reliable and effective.

Thermodynamic structure of the atmospheric boundary layer over a coastal station in India for contrasting sky conditions during different seasons

The thermodynamic structure of the atmospheric boundary layer for contrasting sky conditions over Chennai, a coastal station in the Indian subcontinent, is investigated through conserved variable analysis of equivalent potential temperature and specific humidity. Simultaneous radiosonde and micropulse lidar measurements undertaken in 2018 constitute the primary database in this investigation. One of the most prominent features of this analysis is a consistent occurrence of a double mixed layer structure during the clear-sky and cloudy conditions throughout the year. The occurrence frequency of the double mixing lines is higher during the pre-monsoon season compared with the winter and northeast (NE) monsoon seasons. The advection mainly dominates the formation of double mixing lines during the winter and pre-monsoon seasons. In contrast, convection and advection dominate during the southwest (SW) and NE monsoon seasons. The frequent double mixing lines over Chennai occur mainly from the restratification of the convective boundary layer (CBL) due to the sea-breeze onset and the cloud layer. Occasionally, a triple mixing line structure is also observed during the fair-weather boundary layer (FWBL) of the pre-monsoon and SW monsoon seasons. Among the 355 total observations collected during 2018, the first, second, and third mixing lines occurred 100%, 70%, and 14%, respectively. The thermal internal boundary layer (TIBL), CBL, and FWBL occur ∼50%, ∼97%, and ∼ 30%, respectively. The first mixing line is represented by both the TIBL and CBL, and CBL and FWBL represent the second mixing line, whereas the third mixing line is represented solely by the FWBL. The first and second mixing line shows strong seasonal variations with lower altitudes during the pre-monsoon season and higher altitudes in the SW monsoon season, almost in the same phase as the CBL variation but in the opposite phase of the TIBL variations. The CBL height attains a minimum during the winter season and maximum during the SW monsoon season, while TIBL becomes minimum during the pre-monsoon and SW monsoon seasons and maximum during winter and NE monsoon seasons.

Boundary Layer Height Characteristics in Mexico City from Two Remote Sensing Techniques

The Atmospheric Boundary Layer (ABL) height is a key parameter in air quality research as well as for numerical simulations and forecasts. The identification of thermally stable layers, often with radiosondes, has been a common approach for estimating ABL height, though with limited temporal coverage. Remote sensing techniques offer essentially continuous measurements. Nevertheless, ABL height retrievals from different methods can vary greatly when compared, which is particularly notable for topographically complex terrains, such as that surrounding Mexico City. This study, employing one year of data in Mexico City, reveals that the daytime convective boundary layer height (retrieved from Doppler lidar data) is typically lower than the aerosol layer height (retrieved from ceilometer data). Although both estimated heights evolved diurnally, the more elevated aerosol layer decays more slowly, suggesting that the mechanisms that elevate aerosols are not limited to convective motions. Additionally, both diurnal and seasonal variability are investigated, comparing remotely sensed-retrieved heights with thermally stable layers estimated from radiosonde data. Multiple stable layers often develop, those at higher levels have similar values to the ceilometer-retrieved heights, while stable layers at lower heights are similar to Doppler lidar height retrievals. The present research constitutes the first detailed analysis of ceilometer backscatter and Doppler lidar thresholding methods for estimating ABL height over Mexico City, and our results illustrate the complexity of mixing mechanisms on the ABL in this region of complex orography.

Performance Evaluation of Planetary Boundary Layer Schemes in Simulating Structures of Wintertime Lower Troposphere in Seoul Using One-Hour Interval Radiosonde Observation

We investigated the structures of the wintertime lower troposphere in Seoul, South Korea on 17 and 18 January 2017 by performing 1 h interval radiosonde observation and numerical simulations. In the daytime on 17 January, the height of the convective boundary layer (CBL) sharply and quickly increased when the residual layer became a part of the CBL. From the afternoon on 17 January, moist air with clouds began to substantially intrude in the lower troposphere in Seoul, and radiative heating/cooling weakened. As a result, the mixing of air in the lower troposphere was inhibited and the vertical gradients of potential temperature and water vapor mixing ratio changed little on 18 January. We evaluated the performance of four planetary boundary layer (PBL) parameterization schemes (the Yonsei University (YSU), Mellor–Yamada–Janjić (MYJ), Mellor–Yamada–Nakanishi–Niino (MYNN), and Asymmetric Convective Model version 2 (ACM2) schemes) coupled with the Weather Research and Forecasting model in simulating the structures of the lower troposphere against 1 h interval radiosonde observation. The general tendencies of the air temperature and wind speed in the lower troposphere were well-reproduced in the four simulations. However, the sharp increase in the CBL height did not appear in the four simulations, implying that the process of the residual layer becoming a part of the CBL in the daytime is not well-parameterized. Additionally, the simulated water vapor mixing ratio near the surface was smaller compared with the observation. We found that small-scale turbulence in the CBL, which mixes advected air and pre-existing air, was not reproduced well by the PBL parameterization schemes. Compared with the other simulations, the most accurate air temperature and wind speed were reproduced in the simulation with the MYJ scheme, while the CBL development and moisture advection were reproduced relatively well in the simulation with the MYNN scheme.

Evaluating convective planetary boundary layer height estimations resolved by both active and passive remote sensing instruments during the CHEESEHEAD19 field campaign

Abstract. During the Chequamegon Heterogeneous Ecosystem Energy-balance Study Enabled by a High-density Extensive Array of Detectors 2019 (CHEESEHEAD19) field campaign, held in the summer of 2019 in northern Wisconsin, USA, active and passive ground-based remote sensing instruments were deployed to understand the response of the planetary boundary layer to heterogeneous land surface forcing. These instruments include radar wind profilers, microwave radiometers, atmospheric emitted radiance interferometers, ceilometers, high spectral resolution lidars, Doppler lidars, and collaborative lower-atmospheric mobile profiling systems that combine several of these instruments. In this study, these ground-based remote sensing instruments are used to estimate the height of the daytime planetary boundary layer, and their performance is compared against independent boundary layer depth estimates obtained from radiosondes launched as part of the field campaign. The impact of clouds (in particular boundary layer clouds) on boundary layer depth estimations is also investigated. We found that while all instruments are overall able to provide reasonable boundary layer depth estimates, each of them shows strengths and weaknesses under certain conditions. For example, radar wind profilers perform well during cloud-free conditions, and microwave radiometers and atmospheric emitted radiance interferometers have a very good agreement during all conditions but are limited by the smoothness of the retrieved thermodynamic profiles. The estimates from ceilometers and high spectral resolution lidars can be hindered by the presence of elevated aerosol layers or clouds, and the multi-instrument retrieval from the collaborative lower atmospheric mobile profiling systems can be constricted to a limited height range in low-aerosol conditions.

Dust Characteristics Observed by Unmanned Aerial Vehicle over the Taklimakan Desert

Based on observations from the Unmanned Aerial Vehicle (UAV) together with an environmental particulate matter analyzer (Grimm-180) and Global Positioning System (GPS) sounding balloons, the vertical structure of dust with different particle sizes was explored over the Taklimakan Desert (TD) during an intensive observation from 1 July 2021 to 31 July 2021. The power functions were fitted between the particle counts and particle sizes, indicating negative correlations with an R2 higher than 0.99 under different dust pollution conditions in Tazhong (TZ). The dust concentrations show a sharp vertical increase over the TD during dust pollution; however, more particles with larger sizes are entrained into the air in TZ compared with Minfeng (MF). The total solar radiation during dust pollution days is significantly weakened, accompanied by major modifications in the temperature stratification, which were characterized by low-level cooling (with −2.71 K mean intensity) and high-level heating (with +0.70 K mean intensity). On clear days, the average convective boundary layer (CBL) heights at the TZ and MF are approximately 3.94 and 2.84 km, respectively, and the average stable boundary layer (SBL) height at the TZ and MF are approximately 0.19 and 0.14 km, respectively. With the increasing dust pollution level, the CBL height decreases rapidly while the SBL height shows the opposite trend. The unique ultra-high atmospheric boundary layer structure in daytime provides beneficial conditions for the suspension and vertical transportation of dust over TD. Moreover, a negative correlation between the CBL height and near-surface PM10/PM2.5/PM1.0 concentration in TD is revealed by power function fittings.

Intercomparison of planetary boundary-layer height in Mexico City as retrieved by microwave radiometer, micro-pulse lidar and radiosondes

The planetary boundary layer (PBL) height plays a major role in air quality and weather forecast studies. However, it cannot be directly measured as it can only be determined from the measured profiles of atmospheric parameters such as., temperature, moisture, and aerosol backscatter based on different retrieval mechanisms. This paper presents the PBL features such as., stable boundary layer (SBL), convective boundary layer (CBL), and residual layer (RL) detected from the microwave radiometer (MWR), mini Micro Pulse Lidar (MPL), and radiosondes (RS) and an intercomparison was done between the instruments during the period of 16 February – 31 May 2019 in Mexico City. RS were launched thrice a day during this period at 06 LST, 12 LST, and 18 LST respectively. The PBL heights comparison was classified into 3 categories: MWR-RS, mini MPL-RS, and mini MPL-MWR. The daytime CBL heights for the MWR-RS comparison correlate well (r = 0.96 at 12 LST, r = 0.88 at 18 LST) as well as RL height (r = 0.94 at 06 LST). The CBL and RL heights for the mini MPL-RS and mini MPL-MWR comparison also agree very well (0.89 < r < 0.99). Overall, the corresponding SBL height comparisons at 06 LST yielded lower agreements (r < 0.66) with the lowest correlation values for the mini MPL-MWR comparison (r = 0.39), which is likely due to the inability of the mini MPL to detect aerosol backscatter below ~120 m above ground level (agl) as well as the stability of the atmosphere, which prevents uniform mixing of the aerosol particles.

Characteristics of the atmospheric boundary layer's structure and heating (cooling) rate in summer over the Northern Tibetan Plateau

The atmospheric boundary layer (ABL) is one of the most fundamental yet complex components of the Earth's atmosphere. Hence, studying the ABL has important theoretical value and practical significance. In this paper, the structural characteristics and heating (cooling) rate of the ABL in summer over the Northern Tibetan Plateau (NTP) were analysed using radiosonde observation data from the Amdo and Nagqu regions. The results indicate that the summertime ABL height over the NTP exhibited obvious diurnal variations, with the ABL height during the dry season being greater than that during the rainy season. The maximum convective boundary layer (CBL) height during the daytime reached 3200 m and 2500 m in the dry and rainy seasons, respectively, and the mean maximum CBL height was approximately 2500 m; the maximum stable boundary layer (SBL) height at night reached 900 m, and the mean maximum SBL height was approximately 500 m. The wind speed dominated by westerly wind in the dry season was greater than that dominated by easterly wind in the rainy season, and the zonal (meridional) wind speed (shear) on sunny days was greater than that on cloudy days. The inverse humidity phenomenon occurred in both Amdo and Nagqu, and the strong humidity inversion occurred mainly at midnight on sunny days and at noon on cloudy days. The heating (cooling) rate of the ABL displayed obvious diurnal variations, with the rates being greater on sunny days and lower on cloudy and rainy days. Furthermore, the mean values of the daytime heating rate and nighttime cooling rate of the ABL were relatively equal, indicating that the atmospheric energy budget was, for the most part, balanced.

Estimating the urban atmospheric boundary layer height from remote sensing applying machine learning techniques

This study proposes a new methodology to estimate the Atmospheric Boundary Layer Height (ABLH), discriminating between Convective Boundary Layer and Stable Boundary Layer heights, based on the machine learning algorithm known as Gradient Boosting Regression Tree. The algorithm proposed here uses a first estimation of the ABLH derived applying the gradient method to a ceilometer signal and several meteorological variables to obtain ABLH values comparable to those derived from a microwave radiometer. A deep analysis of the model configuration and its inputs has been performed in order to avoid the model overfitting and ensure its applicability. The hourly and seasonal values and variability of the ABLH values obtained with the new algorithm have been analyzed and compared with the initial estimations obtained using only the ceilometer signal. Mean Relative Errors (MRE) between the ABLH estimated with the new algorithm and microwave radiometer show a daily pattern with their highest values during the night-time (stable situations) and their lowest values along the day-time (convective situations). This pattern has been observed for all the seasons with MRE ranging between −5% and 35%. This result notably improves those ABLH values derived by applying the gradient method to ceilometer data during convective situations and enables the Stable Boundary Layer height detection at night and early morning, instead of only Residual Layer top height. Finally, the model performance has been directly validated in three particular cases: clear-sky day, presence of low-clouds and dust outbreak event. In these three particular situations, ABLH values obtained with the new algorithm follow the pattern obtained with the microwave radiometer presenting very similar values, thus confirming the good model performance. In this way it is feasible by the combination of the proposed method with gradient method, to estimate Convective, Stable and Residual Boundary Layer height from ceilometer data and surface meteorological data in extended network that include ceilometer profiling.

Evaluation of convective boundary layer height estimates using radars operating at different frequency bands

Abstract. Knowledge of the atmospheric boundary layer state and evolution is important for understanding air pollution and low-level cloud development, among other things. There are a number of instruments and methods that are currently used to estimate boundary layer height (BLH). However, no single instrument is capable of providing BLH measurements in all weather conditions. We proposed a method to derive a daytime convective BLH using clear air echoes in radar observations and investigated the consistency of these retrievals between different radar frequencies. We utilized data from three vertically pointing radars that are available at the SMEAR II station in Finland, i.e. the C band (5 GHz), Ka band (35 GHz) and W band (94 GHz). The Ka- or W-band cloud radars are an integral part of cloud profiling stations of pan-European Aerosol, Clouds and Trace Gases Research Infrastructure (ACTRIS). Our method will be utilized at ACTRIS stations to serve as an additional estimate of the BLH during summer months. During this period, insects and Bragg scatter are often responsible for clear air echoes recorded by weather and cloud radars. To retrieve a BLH, we suggested a mechanism to separate passive and independently flying insects that works for all analysed frequency bands. At the lower frequency (the C band) insect scattering has been separated from Bragg scattering using a combination of the radar reflectivity factor and linear depolarization ratio. Retrieved values of the BLH from all radars are in a good agreement when compared to the BLH obtained with the co-located HALO Doppler lidar and ERA5 reanalysis data set. Our method showed some underestimation of the BLH after nighttime heavy precipitation yet demonstrated a potential to serve as a reliable method to obtain a BLH during clear-sky days. Additionally, the entrainment zone was observed by the C-band radar above the CBL in the form of a Bragg scatter layer. Aircraft observations of vertical profiles of potential temperature and water vapour concentration, collected in the vicinity of the radar, demonstrated some agreement with the Bragg scatter layer.

Assessment of Trends and Uncertainties in the Atmospheric Boundary Layer Height Estimated Using Radiosounding Observations over Europe

Trends in atmospheric boundary layer height may represent an indication of climate changes. The related modified interaction between the surface and free atmosphere affects both thermodynamics variables and dilution of chemical constituents. Boundary layer is also a major player in various feedback mechanisms of interest for climate models. This paper investigates trends in the nocturnal and convective boundary layer height at mid-latitudes in Europe using radiosounding profiles from the Integrated Global Radiosounding Archive (IGRA). Atmospheric data from the European Centre for Medium-Range Weather Forecasts (ECMWF) ReAnalysis v5 (ERA5) and from the GCOS Reference Upper-Air Network (GRUAN) Lindenberg station are used as intercomparison datasets for the study of structural and parametric uncertainties in the trend analysis. Trends are calculated after the removal of the lag-1 autocorrelation term for each time series. The study confirms the large differences reported in literature between the boundary layer height estimates obtained with the two different algorithms used for IGRA and ERA5 data: ERA5 shows a density distribution with median values of 350 m and 1150 m for the night and the daytime data, respectively, while the corresponding IGRA median values are of 1150 m and 1750 m. An overall good agreement between the estimated trends is found for nighttime data, while daytime ERA5 boundary layer height estimates over Europe are characterized by a lower spatial homogeneity than IGRA. Parametric uncertainties due to missing data in both the time and space domain are also investigated: the former is not exceeding 1.5 m, while the latter are within 10 m during night and 17 m during the day. Recommendations on dataset filtering based on time series completeness are provided. Finally, the comparison between the Lindenberg data as processed at high-resolution by GRUAN and as provided to IGRA at a lower resolution, shows the significant impact of using high-resolution data in the determination of the boundary layer height, with differences from about 200 m to 450 m for both night and day, as well as a large deviation in the estimated trend.

The impacts of the atmospheric boundary layer on regional haze in North China

An intensive atmospheric boundary layer (ABL) experiment was conducted simultaneously at six stations arranged in a cross shape on the North China Plain (NCP) from 26 November to 26 December 2019. The impacts of the regional ABL structure on heavy haze pollution and the relationship between the ABL height and aerosol accumulation layer (AAL) depth were discussed. Bouts of downdrafts generate a persistent descending elevated inversion layer, helping the maintenance and exacerbation of haze pollution. Continuous weak wind layers contribute to the pollutants accumulation, and low-level jets promote the removal of air pollutants. The unique landform conditions of the NCP are reflected in its regional ABL structure and further affect the spatial distribution of haze pollution. Due to the drainage flow and strong downdrafts, the western stations near the mountains have a colder surface and warmer upper air masses, resulting in a more stable stratification and worse diffusion conditions; these stations also experience a thicker weak wind layer caused by increased friction. Thus, the spatial distribution of haze is heavier in the west and lighter in the east. The convective boundary layer (CBL) height declines evidently during haze episodes, usually lower than 560 m. Furthermore, as the vertical distribution of aerosols is mainly influenced by daytime thermal turbulence and maintained at night, it is appropriate to determine the CBL height using the AAL depth. However, the AAL depth is not consistent with the stable boundary layer height due to the influence of the residual layer at night.

Observed Interactions Among Haze, Fog and Atmospheric Boundary Layer during a Haze-fog Episode in the Yangtze River Delta Region, Eastern China

A severe haze-fog episode occurred in the Yangtze River Delta region of eastern China during 22–30 November, 2018. In this period, the PM2.5 mass concentration and meteorological parameters at the surface were collected at the Station for Observing Regional Processes of the Earth System site in Nanjing. The vertical distributions of PM2.5, humidity and potential temperature below 500 m were observed simultaneously by an unmanned aerial vehicle, and the profile of potential temperature at 1400 local standard time on each day was also observed by radiosonde at the same site. During the first four days, the PM2.5 mass concentration increased, the maximum convective planetary boundary layer height (CBLH) decreased, and the air humidity increased. These are favorable conditions for fog formation. In the latter five days, fog formed on four days, with a lowering of the CBLH and a further increase in PM2.5 mass concentration. We found that the fog top cooling induced a potential temperature jump (i.e., sharp increase of potential temperature) with much warmer temperatures above the cloud top cooling and that this particular thermal structure was maintained until the end of the fog period, which significantly suppressed the daytime development of the planetary boundary layer after fog dissipation. The fog-induced reduction of the CBLH further increased the PM2.5 mass concentration. We also found that the wet deposition of fog on PM2.5 was negligible. The scavenging effect of fog on aerosols only acts during a fog period. When the fog dissipates, the aerosols are liberated from the fog droplets to the atmosphere.

Large eddy simulations of the Martian convective boundary layer: Towards developing a new planetary boundary layer scheme

The grid resolutions typically employed in atmospheric global circulation and mesoscale models are not sufficient to explicitly resolve the turbulence processes within the planetary boundary layer (PBL). turbulent fluxes are, therefore, fully parameterized in those models, based on empirical relationships for the mixing length scale using PBL schemes. However, microscale models use the large-eddy simulation (LES) technique to resolve turbulence processes. Here, we perform LES computations for the daytime Martian planetary boundary layer, ranging from weakly to strong convective conditions, using the Mars implementation of planetWRF, the MarsWRF model. In this study, our main focus is to investigate the turbulence statistics and turbulent spectrum utilizing our LES results. Then, using the computed turbulence kinetic energy and its dissipation rate, a generic formulation for the mixing length scale variation in the Martian convective boundary layer is proposed. This mixing length formulation is used to derive a Mars-specific PBL scheme and its performance is compared to the PBL scheme currently in use by the MarsWRF model, the MRF scheme. The proposed scheme is tested both in global and mesoscale simulations, which are used to evaluate the convective boundary layer height and near-surface meteorology conditions at the InSight landing site. The presently proposed PBL scheme results in an improved prediction of convective boundary layer height that agree better with the observational estimations acquired by radio occultations of Mars Express in comparison to the MRF scheme. Also, the prediction of near-surface winds at the InSight landing site is slightly improved.

Anatomy of the March 2016 severe ozone smog episode in Mexico-City

The diurnal evolution of the planetary boundary layer (PBL) is crucial to air quality studies as it impacts the exchange and distribution of pollutants close to the surface. This paper reports continuous detection of the daytime convective boundary layer height, the stable boundary layer height, and the residual layer height as estimated from the vertical profiles of virtual potential temperature, and moisture retrieved from a microwave radiometer (MWR) in Mexico City for the period 6–18 March 2016. This period included a severe smog episode. We analyzed the anatomy of this episode utilizing continuous air quality measurements recently deployed at elevated locations surrounding the basin of Mexico City, which were used to determine the impact of the background or residual pollutants during the severe smog episode in combination with back trajectory analysis and radar wind profiles data to track transport processes within the Mexico City basin. The first few days prior to the smog episode were impacted by the passage of a deep upper tropospheric trough and strong advection. Shortly before the smog episode, daytime maximum PBL height still reached ~2.5 km above ground level but then dropped to ~1.2–1.7 km above ground level for the most severe pollution days. During the first days with strong advection, the pollutant concentrations were flushed out from the basin and/or could not accumulate (maximum hourly ozone and carbon monoxide mixing ratios of ~50 ppbv and ~0.5 ppmv, respectively). At the departure of the storm, the winds became weaker, and a strong near surface temperature inversion was observed at nighttime increasing the nighttime mixing ratio of carbon monoxide to ~2.5 ppmv and daytime ozone mixing ratio to ~200 ppbv, which resulted in one of the most severe smog episodes in Mexico City over the last decade. Our results point to strong photochemical processes confined to the PBL within the Mexico City basin, whose maximum daytime convective boundary layer heights hardly surpassed the surrounding average mountain top heights.

Estimations of the Mexicali Valley (Mexico) Mixing Height

We report estimations of the Mexicali Valley (Mexico) mixing height for three seasons. Surface and upper air meteorological measurements were carried out nearby Cerro Prieto geothermal power plant during July 2010 (summer), January 2012 (winter), and October 2016 (autumn). Four different methods were applied to estimate the convective boundary layer (CBL) height from radiosonde (RS) profiles: the parcel method, the gradients method, the least-squares variational approach based on the slab model of the CBL structure, and a covariance method. For nocturnal conditions, we used diagnostic models based on friction velocity and Monin–Obukhov length. Under unstable conditions, we obtained (on average) mixing heights of 497 m at 06:00 LST, 1242 m at 12:00 LST, 1404 m at 15:00 LST, and 482 at 18:00 LST during summer; 754 m at 12:00 LST during winter; and 1195 m at 12:00 LST and 13:00 m between the 14:00 and 15:00 LST during the autumn. The results allowed adjusting a semiempirical model to evaluate mixing height from turbulent sensible heat flux and friction velocity data. Our results provide practical tools that could facilitate the application of regulatory dispersion models to assess air quality in the region.

Assessment of New Alternative Scaling Properties of the Convective Boundary Layer: Application to Velocity and Temperature Spectra

There are two models for the surface layer in the convective boundary layer (CBL). First is the standard Monin–Obukhov similarity theory, and second is the McNaughton–Laubach model (Laubach and McNaughton 2009, Boundary-Layer Meteorology 133:219–252; hereafter MNL model) developed based on the complex dynamical system approach to address characteristics of the unstable surface layer. The fundamental difference between the Monin–Obukhov similarity theory and the MNL model is the use of local and non-local parameters for analyzing surface-layer spectra in the CBL. However, there is a need to check applicability of this new model at various flow conditions before the model could be extensively used. Subsequently, applicability of the MNL model is tested in comparison to the standard model using CBL observations from three different regions with increasing terrain complexity (i.e. over flat-terrain, onslope, and ridge-top sites). The MNL model is tested by estimating and using the non-local scaling parameters to collapse the power spectra of velocity and temperature on the frequency–amplitude scale under the generalized hypothesis that convective surface layer depends on non-local outer variables. We find that the u and v spectra for all sites indicate run-to-run similarity of each spectra with MNL scaling irrespective of the height limiting role of local buoyancy on the shape of the spectra. Similarly, w spectra from all sites indicate transitions between the surface friction layer and the outer layer are governed by flow in the entire CBL. The temperature spectra collapse using $$(z \epsilon _o)^{2/3} H_0^{-2}$$ as amplitude scaling and $$kz_i^{1/2}z^{1/2}$$ as wavenumber scaling, is a new observation within the surface friction layer, where the streamwise wavenumber is k, measurement height is z, CBL height is $$z_i$$ , the dissipation rate of turbulence energy in the outer CBL is $$\epsilon _o$$ , and the surface heat flux is $$H_0$$ . These observations corroborate well with the MNL model conjecture that the convective temperature spectra do not depend only on local stability, and CBL parameters affect spectra when a subset of local factors remains constant.