Boundary Layer Air Research Articles

Numerical analysis of the interaction between a droplet and an air boundary layer

The deformation and movement of droplets are widely utilized in many industrial applications. The present work investigates the evolution of a single droplet interacting with an air boundary layer numerically and validated by wind tunnel experiments. The volume of fluid method is employed to study the interaction from the micro-perspective. The influences of airflow velocity, droplet size, and depression angle on interactions are comprehensively discussed. The outcomes indicate that droplet diameter and airflow velocity significantly influence the interaction. Based on the morphological evolution of the droplets, the regimes of the interaction can be classified into three categories. It is shown that the airflow velocity, depression angle, and droplet diameter influence the droplet maximum streamwise spreading length. Furthermore, only the airflow velocity and droplet diameter influence the maximum height. The scaling law for the maximum streamwise spreading factor is revealed. Finally, the velocity profile of the boundary layer above the droplet maximum height is also analyzed, revealing a power-law relationship in its curve. These results provide valuable insight for further investigation on the droplet–air boundary layer interaction.

Simulating Mixed‐Phase Open Cellular Clouds Observed During COMBLE: Evaluation of Parameterized Turbulence Closure

AbstractMarine cold‐air outbreaks, or CAOs, are airmass transformations whereby relatively cold boundary layer (BL) air is transported over relatively warm water. To more deeply understand BL and mixed‐phase cloud properties during CAO conditions, the Cold‐Air Outbreaks in the Marine Boundary Layer Experiment (COMBLE) took place from late 2019 into early 2020. During COMBLE, the U.S. Department of Energy's first Atmospheric Radiation Measurement Mobile Facility (AMF1) was deployed to Andenes, Norway, far downstream (∼1,000 km) from the Arctic pack ice. This study examines the two most intense CAOs sampled at the AMF1 site. The observed BL structures are open cellular with high (∼3–5 km) and cold (−30 to −50 C) cloud tops, and they often have pockets of high liquid water paths (LWPs; up to ∼1,000 g ) associated with strong updrafts and enhanced turbulence. We use a high‐resolution mesoscale model to explore how well four turbulence closure methods represent open cellular clouds. After applying a radar simulator to model outputs for direct evaluation, cloud top properties agree well with AMF1 observations (within ∼10%), but radar reflectivity and LWP agreement is more variable. Results suggest that the turbulent Prandtl number may play an important role for the simulated BL and cloud properties. All simulations produce enhanced precipitation rates that are well‐correlated with a cloud transition. Finally, the eddy‐diffusivity/mass‐flux approach produces the deepest cloud layer and therefore the largest and most coherent cellular structures. We recommend the use of a non‐local turbulence closure approach to better capture turbulent processes in intense CAOs.

Frost crystal growth behavior on a hydrophilic surface over a wide range of cold surface temperature

The initial frosting phenomenon is a discontinuous phase nucleation process, the cold surface temperature and properties have a decisive influence on this phenomenon, especially in the initial frosting stage. With the development of aerospace and energy transportation technology, frost formation at low temperatures (-100 °C∼-30 °C) and ultra-low temperatures (-273 °C∼-100 °C) has gradually attracted the attention of researchers. In this paper, the initial frosting phenomena on hydrophilic surfaces with a contact angle of 10° (CA= 10°) and ordinary (CA= 95°) surfaces are studied experimentally in a wide range of cold surface temperatures (-190 °C∼-30 °C). Four modes are confirmed: cold surface condensation frosting, cold surface sublimation frosting, air boundary layer condensation frosting and air boundary layer sublimation frosting. It is also found that the four frosting modes do not appear in turn with the decrease of the cold surface temperature, but two or more frosting modes appear at the same time. And the surface contact angle has an important influence on the frosting mode. The initial frost crystal morphology mainly depends on the cold surface temperature and the corresponding frosting mode. Four different forms of frost crystals are observed: hexagonal prism (feather), branch (pine needle), cluster (shrub) and floc (grape), in which the cluster frost crystal is more sensitive to the surface contact angle and can appear in different temperature ranges due to different contact angles. Based on the statistics of the size, quantity, and distribution of the initial frost crystals, it is found that -70 °C is a major turning point for frost formation from the cold surface sublimation frosting to the air boundary layer sublimation frosting, and an important change has taken place near this point. Furthermore, it affects the shape and size distribution of frost crystals. These findings are of great significance for the study and understanding of frost crystal growth mechanism in the initial stage of frost formation at low and ultra-low temperatures.

Impacts of Anthropogenic Emissions and Open Biomass Burning in South Asia and Southeast Asia on Air Quality and Meteorology Over Southern China

AbstractThe emissions from South Asia and Southeast Asia significantly impact air quality and meteorological conditions in China. However, the individual or joint contributions of anthropogenic emissions and/or biomass burning from outside China to surface particulate matter with aerodynamic diameters less than 2.5 μm (PM2.5) and aerosol optical depth (AOD) over southern China have not been fully investigated. Here, five experiments were designed to investigate the impacts of these emissions in January (winter), March (pre‐monsoon), and October (post‐monsoon) for 2017. Aerosols from South Asia and Southeast Asia contributed less to southern China during winter and post‐monsoon seasons, whereas the wind patterns and emission intensity during the pre‐monsoon season were conducive to the transport of aerosols. During pre‐monsoon season, the total emissions contributed approximately 5.0 μg m−3 to surface PM2.5 in Xizang Province. Biomass burning in Southeast Asia increased PM2.5 in Yunnan Province by 37.9 μg m−3, while anthropogenic emissions increased it by 8.9 μg m−3. Transboundary aerosols can be transported to Xizang Province and Yunnan Province, primarily influencing PM2.5 below 2 km height. It mainly affected PM2.5 levels above the planetary boundary layer over southeast China. Aerosols from outside China can account for 79.5% and 54.8% AOD in Yunnan Province and southeast China, respectively. These aerosols reduced surface incident solar radiation by approximately 6%, leading to decreases in air temperature, wind speed, and boundary layer height. The findings are only applicable to the pre‐monsoon season.

An air quality and boundary layer dynamics analysis of the Los Angeles basin area during the Southwest Urban NOx and VOCs Experiment (SUNVEx)

Abstract. The NOAA Chemical Sciences Laboratory (CSL) conducted the Southwest Urban NOx and VOCs Experiment (SUNVEx) to study emissions and the role of boundary layer (BL) dynamics and sea-breeze (SB) transitions in the evolution of coastal air quality. The study presented utilizes remote sensing and in situ observations in Pasadena, California. Separate analyses are conducted on the synoptic conditions during ozone (O3) exceedance (>70 ppb) and non-exceedance (<70 ppb) days, as well as the fine-structure variability of in situ chemistry measurements during BL growth and SB transitions. Diurnal analyses spanning August 2021 revealed a markedly different wind direction during evenings preceding O3 exceedance (northerly) versus non-exceedance (easterly) days. Increased O3 occurred simultaneously with warmer and drier conditions, a reduction in winds, and an increase in volatile organic compounds (VOCs) and fine particulate matter (PM2.5). While the average BL height was lower and surface pressure was higher, the day-to-day variability of these quantities led to an overall weak statistical relationship. Investigations focused on the fine-structure variability of in situ chemistry measurements superimposed on background trends were conducted using a novel multivariate spectral coherence mapping (MSCM) technique that combined the spectral structure of two or more independent measurements through a wavelet analysis as reported by maximum-normalized scaleograms. A case study was chosen to illustrate the MSCM technique, where the dominant peaks in scaleograms were identified and compared to BL height during the growth phase. The temporal widths of peaks (τmax) derived from VOC and nitrogen oxide (NOx) scaleograms, as well as scaleograms combining VOCs, NOx, and variations in BL height, indicated a broadening with respect to time as the BL increased in depth. A separate section focused on comparisons between τmax and BL height during August 2021 revealed uncorrelated or weakly correlated scatter, except in the case of VOCs when really large τmax and relatively deep BL heights were ignored. Instances of large τmax and relatively deep BL heights occurred near sunrise and as onshore flow entered Pasadena, respectively. Wind transitions likely influenced both the dynamical evolution of the BL and tracer advection and thus offer additional challenges when separating factors contributing to the fine structure. Other insights gained from this work include observations of descending wind jets from the San Gabriel Mountains that were not resolved by the High-Resolution Rapid Refresh (HRRR) model and the derivation of intrinsic properties of oscillations observed in NOx and O3 during the interaction between an SB and enhanced winds above the BL that flowed in opposition to the SB.

Dynamic Wall Shear Stress Measurement Using Event-Based 3D Particle Tracking

We describe the implementation of a 3d Lagrangian particle tracking (LPT) system based on event-based vision (EBV) sensors and demonstrate its application for the near-wall characterization of a turbulent boundary layer (TBL) in air. The viscous sublayer of the TBL is illuminated by a thin light sheet that grazes the surface of a glass window inserted into the flat wall of the wind tunnel wall. The data simultaneously captured by three synchronized EBV-cameras is used to reconstruct the 3d particle tracks within 0.4 mm of the wall on a field of view of 12.0 mm by 7.5 mm. The velocity and position of particles within the viscous sub-layer permit the estimation of the local, unsteady wall shear stress vector under the assumption of linearity between particle velocity and wall shear stress. Thereby, 2d distributions and higher order statistics of the unsteady wall shear stress are obtained. The systematic underestimation of the spanwise shear stress fluctuation can be explained on the basis of DNS data; a methodology for its mitigation is proposed. The employed EBV hardware coupled with suited LPT algorithms provide data quality on par with currently used, considerably more expensive, high-speed framing cameras.

Experimental investigation of interactions between a water droplet and an airflow boundary layer

The deformation and movement of droplets is widely relevant in many fields of research. The present work experimentally investigates the evolution of a single droplet interacting with an air boundary layer. A series of experiments are carried out using a high-speed photography technique to determine the effects of the airflow velocity, drop height, and droplet size. The morphological characteristics can be classified into three types according to the experiments. The outcomes indicate that both the drop height and the airflow velocity significantly influence the maximum streamwise spreading length, but only the drop height has an impact on the maximum lateral spreading width. The maximum streamwise spreading factor follows a power function relationship with WeRe−0.5. In addition, the crater maximum streamwise and lateral spreading diameters are mainly influenced by the drop height. An energy conversion model is established by considering the effects of the aerodynamic drag force, surface tension, and viscous force. This study provides experimental reference data for the scenario of a droplet interacting with an air boundary layer.

An improved Baldwin–Lomax algebraic wall model for high-speed canonical turbulent boundary layers using established scalings

In this work, we employ well-established relations for compressible turbulent mean flows, including the velocity transformation and algebraic temperature–velocity (TV) relation, to systematically improve the algebraic Baldwin–Lomax (BL) wall model for high-speed zero-pressure-gradient air boundary layers. Any new functions or coefficients fitted by ourselves are avoided. Twelve published direct numerical simulation (DNS) datasets are employed for a priori inspiration and a posteriori examination, with Mach numbers up to 14 under adiabatic, cold and heated walls. The baseline BL model is the widely used one with semilocal scalings. Three targeted modifications are made. First, we employ a total-stress-based transformation (Griffin et al., Proc. Natl Acad. Sci. USA, vol. 118, issue 34, 2021, e2111144118) to the inner-layer eddy viscosity for improved scaling up to the logarithmic region. Second, we utilize the van Driest transformation in the outer layer based on the compressible defect velocity scaling. Third, considering the difficulty in modelling the rapidly varying and singular turbulent Prandtl number near the temperature peak in cold-wall cases, we design a two-layer strategy and use the TV relation to formulate the inner-layer temperature. Numerical results prove that the modifications take effect as designed. The prediction accuracy for mean streamwise velocity is notably improved for diabatic cases, especially in the logarithmic region. Moreover, a significant improvement in mean temperature is realized for both adiabatic and diabatic cases. The mean relative errors of temperature to DNS for all cases are down to 0.4 % in the logarithmic wall-normal coordinate and 3.4 % in the outer coordinate, around one-third of those in the baseline model.

Elevated Mixed Layers during Great Lake Lake-Effect Events: An Investigation and Case Study from OWLeS

Abstract Data from rawinsondes launched during intensive observation periods (IOPs) of the Ontario Winter Lake-Effect Systems (OWLeS) field project reveal that elevated mixed layers (EMLs) in the lower troposphere were relatively common near Lake Ontario during OWLeS lake-effect events. Conservatively, EMLs exist in 193 of the 290 OWLeS IOP soundings. The distribution of EML base pressure derived from the OWLeS IOP soundings reveals two classes of EML, one that has a relatively low-elevation base (900–750 hPa) and one that has a relatively high-elevation base (750–500 hPa). It is hypothesized that the former class of EML, which is the focus of this research, is, at times, the result of mesoscale processes related to individual Great Lakes. WRF reanalysis fields from a case study during the OWLeS field project provide evidence of two means by which low-elevation base EMLs can originate from the lake-effect boundary layer convection and associated mesoscale circulations. First, such EMLs can form within the upper-level outflow branches of mesoscale solenoidal circulations. Evacuated Great Lakes–modified convective boundary layer air aloft then lies above ambient air of a greater static stability, forming EMLs. Second, such EMLs can form in the absence of a mesoscale solenoidal circulation when Great Lake–modified convective boundary layers overrun ambient air of a greater density. The reanalysis fields show that EMLs and layers of reduced static stability tied to Great Lakes–modified convective boundary layers can extend downwind for hundreds of kilometers from their areas of formation. Operational implications and avenues for future research are discussed.

Exploring hail and lightning diagnostics over the Alpine-Adriatic region in a km-scale climate model

Abstract. The north and south of the Alps, as well as the eastern shores of the Adriatic Sea, are hot spots of severe convective storms, including hail and lightning associated with deep convection. With advancements in computing power, it has become feasible to simulate deep convection explicitly in climate models by decreasing the horizontal grid spacing to less than 4 km. These kilometer-scale models improve the representation of orography and reduce uncertainties associated with the use of deep convection parameterizations. In this study, we perform km-scale simulations for eight observed cases of severe convective storms (seven with and one without observed hail) over the Alpine-Adriatic region. The simulations are performed with the climate version of the regional model Consortium for Small-scale Modeling (COSMO) that runs on graphics processing units (GPUs) at a horizontal grid spacing of 2.2 km. To analyze hail and lightning we have explored the hail growth model (HAILCAST) and lightning potential index (LPI) diagnostics integrated with the COSMO-crCLIM model. Comparison with available high-resolution observations reveals good performance of the model in simulating total precipitation, hail, and lightning. By performing a detailed analysis of three of the case studies, we identified the importance of significant meteorological factors for heavy thunderstorms that were reproduced by the model. Among these are the moist unstable boundary layer and dry mid-level air, the topographic barrier, as well as an approaching upper-level trough and cold front. Although COSMO HAILCAST tends to underestimate the hail size on the ground, the results indicate that both HAILCAST and LPI are promising candidates for future climate research.

Experimental investigation and simulation of an ion anemometer and studying the effect of the air boundary layer on its performance

Experimental investigation and simulation of an ion anemometer and studying the effect of the air boundary layer on its performance

Thermal performance augmentation of solar air heater with curved path

Poor heat convection from absorber to airstream limits the benefits of solar air heater (SAH) utilization. This study aims to maximize the extracted heat from the absorber by increasing air flow mixing rate and boundary layers separation process. Thereby, SAHs with smooth curvilinear paths of radiuses of 5 cm (SAH5), 15 cm (SAH15), and 25 cm (SAH25) are introduced. Moreover, CFD is employed to conduct the present study, and model predictions are found to yield good agreement with respect to well-known correlations from literature. The findings reveal that thermal performance of curved SAH does not vary significantly after path's radius of 15 cm. Additionally, data of Nusselt number of SAH15 are found 3–7%, and 55–83% greater than those of SAH5 and conventional SAH, respectively. Nevertheless, data of friction factor of SAH15 are found 3.2–9.6% lower than those of SAH5. Overall, SAH15 contains smaller stagnant areas, yields better flow distribution, and owns higher overall thermal performance with respect to SAH5.

Integrated Numerical Investigation on the Aerodynamics Characteristic and Vortex Development of Airfoil using Spalart-Allmaras Model

In this paper, the numerical simulation successfully obtained good results in analyzing fluid flow over the airfoil. Detailed explanations of simulation steps were also presented. The flow characteristics over three airfoil models were numerically simulated in this work: NACA0021, NACA2409, NACA2409+Fowler flap. The reliable Spalart-Allmaras (S-A) turbulent model was used and validated using reported data from experimental in terms of lift and drag coefficients. In this regard, the discrepancies of less than 10% were obtained for both coefficients, respectively. The boundary layer separation, vortex development, and air separation were clearly captured. The results of symmetric airfoil showed that the vortex shedding regimes occurred at α = 8o, and the stall critical-angle was about 14o. The value was higher for the NACA2409, where the airflow exhibited a relatively more stable behavior. Moreover, it is evident that flap addition altered lift-drag characteristics. The value of the lift-to-drag ratio increased due to the increase of Cl and the reduction of Cd. The parametric study was done on the α and flap deflection angle to attain the desirable airfoil configuration. The maximum result of airfoil configuration was obtained on the NACA2409 at α = 12o with 100 flap deflection angle while it enhanced the lift coefficient by about 54%. This result strengthens the robustness of the S-A turbulence model and projects the use of the S-A model for dealing with the aerodynamics analysis. This study is beneficial for initial aircraft design on the aerodynamics aspect of a wing.

The Possibilities of Using the Minimax Method to Diagnose the State of the Atmosphere

The article is devoted to the discussion of the possibilities of approbation of one of the probabilistic methods of verification of evaluation works - the minimax method or the method of establishing the minimum risk of making erroneous diagnoses of the instability of the planetary boundary layer of air. Within the framework of this study, the task of probabilistic forecasting of diagnostic parameters and their combinations, leading in their totality to the formation of an unstable state of the planetary boundary layer of the atmosphere, was carried out. It is this state that, as shown by previous studies, a priori contribution to the development of a number of weather phenomena dangerous for society (squalls, hail, heavy rains, etc.). The results of applying the minimax method made it possible to identify a number of parameters, such as the intensity of circulation, the activity of the Earth’s magnetosphere, and the components of the geostrophic wind velocity, the combination of which led to the development of instability. In the future, it is possible to further expand the number of diagnosed parameters to identify more sensitive elements. In this sense, the minimax method, the usefulness of which is shown in this study, can be considered as one of the preparatory steps for the subsequent more detailed method for forecasting individual hazardous weather phenomena.

A modelling framework for local thermal comfort assessment related to bicycle helmet use

Thermal discomfort due to accumulated sweat increasing head skin wettedness may contribute to low wearing rates of bicycle helmets. Using curated data on human head sweating and helmet thermal properties, a modelling framework for the thermal comfort assessment of bicycle helmet use is proposed. Local sweat rates (LSR) at the head were predicted as the ratio to the gross sweat rate (GSR) of the whole body or by sudomotor sensitivity (SUD), the change in LSR per change in body core temperature (Δtre). Combining those local models with Δtre and GSR output from thermoregulation models, we simulated head sweating depending on the characteristics of the thermal environment, clothing, activity, and exposure duration. Local thermal comfort thresholds for head skin wettedness were derived in relation to thermal properties of bicycle helmets. The modelling framework was supplemented by regression equations predicting the wind-related reductions in thermal insulation and evaporative resistance of the headgear and boundary air layer, respectively. Comparing the predictions of local models coupled with different thermoregulation models to LSR measured at the frontal, lateral and medial head under bicycle helmet use revealed a large spread in LSR predictions predominantly determined by the local models and the considered head region. SUD tended to overestimate frontal LSR but performed better for lateral and medial head regions, whereas predictions by LSR/GSR ratios were lower and agreed better with measured frontal LSR. However, even for the best models root mean squared prediction errors exceeded experimental SD by 18–30%. From the high correlation (R > 0.9) of skin wettedness comfort thresholds with local sweating sensitivity reported for different body regions, we derived a threshold value of 0.37 for head skin wettedness. We illustrate the application of the modelling framework using a commuter-cycling scenario, and discuss its potential as well as the needs for further research.

Experimental investigation of turbulent flow in a rectangular bonneted slide gate and eliminating random fluctuating loads

Bonneted slide gates are widely used in dam bottom outlets for flow regulation. High kinetic energy flow generates turbulence, which induces vibrations in bonneted slide gates under partially open conditions. Significant vibrations can indicate problems and cause damage or expedited deterioration of the gate if left unchecked. This experimental study focuses on two aspects: turbulent flow formation due to rectangular bonneted slide gates and elimination of random fluctuating loads on the gate using a turbulence inhibitor. Five reservoir heads were combined with nine gate openings (10%–90%) to produce 45 different discharges of 12–182 l/s. The different types of the gate flow create different turbulence and random loads based on the water–air mixture and the kinetic energy of the flow. Flow analysis indicates considerable random static pressure fluctuations under the slide gate. The results show the formation of swirling secondary flows in the gate's side guide slots, and their movement and collision with the gate are the leading cause of the turbulent multiphase flow and random fluctuations in static pressure. Also, the use of turbulence inhibitors on the gate can prevent the formation of secondary flows, which results in a 98.92% reduction in the amplitude of random static pressure fluctuations and the elimination of the random fluctuating loads on the slide gate. After removing the swirling secondary flows, a stable air boundary layer was formed under the gate. Finally, the gate flow changed from a turbulent two-phase flow to a steady single-phase jet of water.

Evaluation of the boundary layer mixing height and air pollution in Arak, Iran

Evaluation of the boundary layer mixing height and air pollution in Arak, Iran

A Multimodel Investigation of Asian Summer Monsoon UTLS Transport Over the Western Pacific

AbstractThe Asian summer monsoon (ASM) as a chemical transport system is investigated using a suite of models in preparation for an airborne field campaign over the Western Pacific. Results show that the dynamical process of anticyclone eddy shedding in the upper troposphere rapidly transports convectively uplifted Asian boundary layer air masses to the upper troposphere and lower stratosphere over the Western Pacific. The models show that the transported air masses contain significantly enhanced aerosol loading and a complex chemical mixture of trace gases that are relevant to ozone chemistry. The chemical forecast models consistently predict the occurrence of the shedding events, but the predicted concentrations of transported trace gases and aerosols often differ between models. The airborne measurements to be obtained in the field campaign are expected to help reduce the model uncertainties. Furthermore, the large‐scale seasonal chemical structure of the monsoon system is obtained from modeled carbon monoxide, a tracer of the convective transport of pollutants, which provides a new perspective of the ASM circulation, complementing the dynamical characterization of the monsoon.

Superhydrophobic Strategy for Nature-Inspired Rotating Microfliers: Enhancing Spreading, Reducing Contact Time, and Weakening Impact Force of Raindrops

Wind-dispersal of seeds is a remarkable strategy in nature, enlightening the construction of microfliers for environmental monitoring. However, the flight of these microfliers is greatly affected by climatic conditions, especially in rainy days, they suffer serious raindrop impact. Here, a hierarchical superhydrophobic surface is fabricated and a novel strategy is demonstrated that the superhydrophobic coating can enhance spreading while reduce contact time and impact force of raindrops, all of which are beneficial for the rotating microfliers. When the surface rotating speed exceeds a critical value, the effect of centrifugal force becomes considerable so that the droplet spreading is enhanced. The rotating superhydrophobic surface can rotate an impacting droplet by the tangential drag force from the air boundary layer, and the rotation of the droplet generates a negative pressure zone inside it, reducing the contact time by more than 30%. The impact force by the droplet on the rotating superhydrophobic surface also has a remarkable reduction of 53% compared to that on unprocessed hydrophilic surfaces, which helps maintain the flight stability of the microfliers. This work pioneers in revealing the droplet impact effect on rotating microflier surfaces and demonstrates the effectiveness of protecting microfliers with superhydrophobic coatings, which shall guide the manufacture and flight of microfliers in rainy conditions.

The impacts of urban anthropogenic heat and surface albedo change on boundary layer meteorology and air pollutants in the Beijing-Tianjin-Hebei region

In this study, present and future anthropogenic heat flux (AH) and surface albedo (SA) in Beijing are estimated and incorporated into the WRF-Chem/UCM model. A series of nested model simulations are conducted at 9 km and 27 km horizontal grid resolution with 30 layers in the vertical to explore the sensitivity of model results to urban surface characteristics during the study period from 23 February to 13 March 2014, and the impacts of anthropogenic heat and surface albedo changes in the future on urban meteorology and air pollutants in Beijing. Model validation demonstrates that the model with the up-to-date urban surface parameters apparently improves predictions for both meteorological variables and chemical species compared with model results with default urban parameters. The combined effect of SA increase and AH decrease leads to consistent decreases in T2 (air temperature at 2 m above ground), WS10 (wind speed at 10 m) and PBLH (planetary boundary layer height) by up to 1.2 °C, 0.4 m s−1 and 200 m (33%), respectively, in Beijing averaged over the study period. In response to the meteorological changes, near surface PM2.5 (particulate matter with aerodynamic equivalent diameter below 2.5 μm) concentration increases by up to 15.6 μg m−3 (8.6%), whereas O3 concentration decreases by 2.6 ppb (25.1%), respectively. In daytime, near surface cooling due to SA increase and AH decrease causes a strong downdraft over Beijing urban areas, a wind divergence near the surface and an upslope wind over the western mountain, resulting in aerosol transport from urban district to western mountain slope. At nighttime, surface cooling is weaker than that in daytime, leading to a weaker downdraft over Beijing and a slight downslope wind in the mountain. The combined changes in SA and AH lead to a decrease in T2 by 1.2 °C, a decrease of O3 concentration by approximately 20%, and increases in PM2.5 and its components by 10–17% in daytime averaged over the urban areas of Beijing during the study period, and the above changes are weaker during haze days than those in clean days. The changes in urban surface parameters in Beijing in the future may be benefit for mitigating heat wave and O3 pollution risks, but could increase fine aerosol level and human health risk, which require an optimal strategy for co-benefits of climate change mitigation and air pollution control in future city planning.