Construction and Validation of an Urban Area Flow and Dispersion Model on Building Scales

This paper presents a numerical model that simulates the wind fields, turbulence fields, and dispersion of gaseous substances in urban areas on building to city block scales. A Computational Fluid Dynamics（CFD） approach using the steady-state, Reynolds-Averaged Navier-Stokes（RANS） equations with the standard k-ε turbulence model within control volumes of non-uniform cuboid shapes has been employed. Dispersion field is computed by solving an unsteady transport equation of passive scalar. Another approach based on Gaussian plume model is used to correct the turbulent Schmidt number of tracer, in order to improve the dispersion simulation. The experimental data from a wind tunnel under neutral conditions are used to validate the numerical results of velocity, turbulence, and dispersion fields. The numerical results show a reasonable agreement with the wind tunnel data. The deviation of concentration between the simulation with corrected turbulent Schmidt number and the wind tunnel experiments may arise from 1） imperfect point sources, 2） heterogeneous turbulent difusivity, and 3） the constant turbulent Schmidt assumption used in the model.     

Numerical study of flow and gas diffusion in the near-wake behind an isolated building

To assist validation of the experimental data of urban pollution dispersion, the effect of an isolated building on the flow and gaseous diffusion in the wake region have been investigated numerically in the neutrally stratified rough-walled turbulent boundary layer. Numerical studies were carried out using Computational Fluid Dynamics (CFD) models. The CFD models used for the simulation were based on the steady-state Reynolds-Average Navier-Stoke equations (RANS) with κ-ε turbulence models; standard κ-ε and RNG κ-ε models. Inlet conditions and boundary conditions were specified numerically to the best information available for each fluid modeling simulation. A gas pollutant was emitted from a point source within the recirculation cavity behind the building model. The accuracy of these simulations was examined by comparing the predicted results with wind tunnel experimental data. It was confirmed that simulation using the model accurately reproduces the velocity and concentration diffusion fields with a fine-mish resolution in the near wake region. Results indicated that there is a good agreement between the numerical simulation and the wind tunnel experiment for both wind flow and concentration diffusion. The results of this work can help to improve the understanding of mechanisms of and simulation of pollutant transport in an urban environment.     

Non-Gaussian Lagrangian Stochastic Model for Wind Field Simulation in the Surface Layer

Wind field simulation in the surface layer is often used to manage natural resources in terms of air quality,gene flow(through pollen drift),and plant disease transmission(spore dispersion).Although Lagrangian stochastic(LS)models describe stochastic wind behaviors,such models assume that wind velocities follow Gaussian distributions.However,measured surface-layer wind velocities show a strong skewness and kurtosis.This paper presents an improved model,a non-Gaussian LS model,which incorporates controllable non-Gaussian random variables to simulate the targeted non-Gaussian velocity distribution with more accurate skewness and kurtosis.Wind velocity statistics generated by the non-Gaussian model are evaluated by using the field data from the Cooperative Atmospheric Surface Exchange Study,October 1999 experimental dataset and comparing the data with statistics from the original Gaussian model.Results show that the non-Gaussian model improves the wind trajectory simulation by stably producing precise skewness and kurtosis in simulated wind velocities without sacrificing other features of the traditional Gaussian LS model,such as the accuracy in the mean and variance of simulated velocities.This improvement also leads to better accuracy in friction velocity(i.e.,a coupling of three-dimensional velocities).The model can also accommodate various non-Gaussian wind fields and a wide range of skewness–kurtosis combinations.Moreover,improved skewness and kurtosis in the simulated velocity will result in a significantly different dispersion for wind/particle simulations.Thus,the non-Gaussian model is worth applying to wind field simulation in the surface layer.     

Simulating urban flow and dispersion in Beijing by coupling a CFD model with the WRF model

The airflow and dispersion of a pollutant in a complex urban area of Beijing, China, were numerically examined by coupling a Computational Fluid Dynamics （CFD） model with a mesoscale weather model. The models used were Open Source Field Operation and Manipulation （OpenFOAM） software package and Weather Research and Forecasting （WRF） model. OpenFOAM was firstly validated against wind-tunnel experiment data. Then, the WRF model was integrated for 42 h starting from 0800 LST 08 September 2009, and the coupled model was used to compute the flow fields at 1000 LST and 1400 LST 09 September 2009. During the WRF-simulated period, a high pressure system was dominant over the Beijing area. The WRF-simulated local circulations were characterized by mountain valley winds, which matched well with observations. Results from the coupled model simulation demonstrated that the airflows around actual buildings were quite different from the ambient wind on the boundary provided by the WRF model, and the pollutant dispersion pattern was complicated under the influence of buildings. A higher concentration level of the pollutant near the surface was found in both the step-down and step-up notches, but the reason for this higher level in each configurations was different： in the former, it was caused by weaker vertical flow, while in the latter it was caused by a downward-shifted vortex. Overall, the results of this study suggest that the coupled WRF-OpenFOAM model is an important tool that can be used for studying and predicting urban flow and dispersions in densely built-up areas.     

A multi-scale urban atmospheric dispersion model for emergency management

To assist emergency management planning and prevention in case of hazardous chemical release into the atmosphere,especially in densely built-up regions with large populations,a multi-scale urban atmospheric dispersion model was established.Three numerical dispersion experiments,at horizontal resolutions of 10 m,50 m and 3000 m,were performed to estimate the adverse effects of toxic chemical release in densely built-up areas.The multi-scale atmospheric dispersion model is composed of the Weather Forecasting and Research （WRF） model,the Open Source Field Operation and Manipulation software package,and a Lagrangian dispersion model.Quantification of the adverse health effects of these chemical release events are given by referring to the U.S.Environmental Protection Agency＇s Acute Exposure Guideline Levels.The wind fields of the urban-scale case,with 3 km horizontal resolution,were simulated by the Beijing Rapid Update Cycle system,which were utilized by the WRF model.The sub-domain-scale cases took advantage of the computational fluid dynamics method to explicitly consider the effects of buildings.It was found that the multi-scale atmospheric dispersion model is capable of simulating the flow pattern and concentration distribution on different scales,ranging from several meters to kilometers,and can therefore be used to improve the planning of prevention and response programs.     

An Application of the RAMS/FLUENT System on the Multi-Scale Numerical Simulation of the Urban Surface Layer A Preliminary Study

The Regional Atmospheric Modeling System （RAMS） and the computational fluid dynamics （CFD） codes known as FLUENT are combinatorially applied in a multi-scale numerical simulation of the urban surface layer （USL）. RAMS and FLUENT are combined as a multi-scale numerical modeling system, in which the RAMS simulated data are delivered to the computational model for FLUENT simulation in an offline way. Numerical simulations are performed to present and preliminarily validate the capability of the multi-scale modeling system, and the results show that the modeling system can reasonably provide information on the meteorological elements in an urban area from the urban scale to the city-block scale, especially the details of the turbulent flows within the USL.     

Study on Effects of Building Morphology on Urban Boundary Layer Using an Urban Canopy Model

An urban canopy model is incorporated into the Nanjing University Regional Boundary Layer Model. Temperature simulated by the urban canopy model is in better agreement with the observation, especially in the night time, than that simulated by the traditional slab model. The coupled model is used to study the effects of building morphology on urban boundary layer and meteorological environment by changing urban area, building height, and building density. It is found that when the urban area is expanded, the urban boundary layer heat flux, thermal turbu- lence, and the turbulent momentum flux and kinetic energy all increase or enhance, causing the surface air temperature to rise up. The stability of urban atmospheric stratiˉcation is a?ected to diffierent extent at diffierent times of the day. When the building height goes up, the aerodynamic roughness height, zero plane displacement height of urban area, and ratio of building height to street width all increase. Therefore, the increase in building height results in the decrease of the surface heat flux, urban surface temperature, mean wind speed, and turbulent kinetic energy in daytime. While at night, as more heat storage is released by higher buildings,thermal turbulence is more active and surface heat flux increases, leading to a higher urban temperature. As the building density increases, the aerodynamic roughness height of urban area decreases, and the effect of urban canopy on radiation strengthens. The increase of building density results in the decrease in urban surface heat flux, momentum flux, and air temperature, the increase in mean wind speed, and the weakening of turbulence in the daytime. While at night, the urban temperature increases due to the release of more heat storage.     

Observation and Simulation of Near-Surface Wind and Its Variation with Topography in Urumqi,West China

Near-surface wind measurements obtained with five 100-m meteorology towers, 39 regional automatic stations, and simulations by the Weather Research and Forecasting(WRF) model were used to investigate the spatial structure of topography-driven flows in the complex urban terrain of Urumqi, China. The results showed that the wind directions were mainly northerly and southerly within the reach of 100 m above ground in the southern suburbs, urban area, and northern suburbs, which were consistent with the form of the Urumqi gorge. Strong winds were observed in southern suburbs, whereas the winds in the urban,northern suburbs, and northern rural areas were weak. Static wind occurred more frequently in the urban and northern rural areas than in the southern suburbs. In the southern suburbs, wind speed was relatively high throughout the year and did not show significant seasonal variations. The average annual wind speed in this region varied among 1.9–5.5, 1.1–3.6, 1.2–4.3, 1.2–4.3, and 1.1–3.5 m s~(-1) with in the reach of 100 m above ground at Yannanlijiao, Shuitashan, Liyushan, Hongguangshan, and Midong, respectively. The flow characteristics comprised more airflows around the mountain, where the convergence and divergence were dominated by the terrain in eastern and southwestern Urumqi. Further analysis showed that there was a significant mountain–valley wind in spring, summer, and autumn, which occurred more frequently in spring and summer for 10–11 h in urban and northern suburbs. During daytime, there was a northerly valley wind,whereas at night there was a southerly mountain wind. The conversion time from the mountain wind to the valley wind was during 0800–1000 LST(Local Standard Time), while the conversion from the valley wind to the mountain wind was during 1900–2100 LST. The influence of the mountain–valley wind in Urumqi City was most obvious at 850 h Pa, according to the WRF model.     

The multi-scale numerical modeling system for research on the relationship between urban planning and meteorological environment

Considering the urban characteristics, a customized multi-scale numerical modeling system is established to simulate the urban meteorological environment.THE system mainly involves three spatial scales;the urban scale, urban sub-domain scale, and single to few buildings scale. In it, different underlying surface types are employed, the building drag factor is used to replace its roughness in the influence on the urban wind field, the effects of building distribution, azimuth and screening of shortwave radiation are added, and the influence of anthropogenic heating is also taken into account. All the numerical tests indicate that the simulated results are reasonably in agreement with the observational data, so the system can be used to simulate the urban meteorological environment. Making use of it, the characteristics of the meteorological environment from the urban to urban sub-domain scales, even the among-buildings scale,can be recognized. As long as the urban planning scheme is given, the corresponding simulated results can be obtained so as to meet the need of optimizing urban planning.     

ANALYSIS OF THE EFFECT OF 3DVAR AND ENSRF DIRECT ASSIMILATION OF RADAR DATA ON THE FORECAST OF A HEAVY RAINFALL EVENT

The present study designs experiments on the direct assimilation of radial velocity and reflectivity data collected by an S-band Doppler weather radar (CINRAD WSR-98D) at the Hefei Station and the reanalysis data produced by the United States National Centers for Environmental Prediction using the Weather Research and Forecasting (WRF) model, the WRF model with a three-dimensional variational (3DVAR) data assimilation system and the WRF model with an ensemble square root filter (EnSRF) data assimilation system. In addition, the present study analyzes a Meiyu front heavy rainfall process that occurred in the Yangtze -Huaihe River Basin from July 4 to July 5, 2003, through numerical simulation. The results show the following. (1) The assimilation of the radar radial velocity data can increase the perturbations in the low-altitude atmosphere over the heavy rainfall region, enhance the convective activities and reduce excessive simulated precipitation. (2) The 3DVAR assimilation method significantly adjusts the horizontal wind field. The assimilation of the reflectivity data improves the microphysical quantities and dynamic fields in the model. In addition, the assimilation of the radial velocity and reflectivity data can better adjust the wind fields and improve the intensity and location of the simulated radar echo bands. (3) The EnSRF assimilation method can assimilate more small-scale wind field information into the model. The assimilation of the reflectivity data alone can relatively accurately forecast the rainfall centers. In addition, the assimilation of the radial velocity and reflectivity data can improve the location of the simulated radar echo bands. (4) The use of the 3DVAR and EnSRF assimilation methods to assimilate the radar radial velocity and reflectivity data can improve the forecast of precipitation, rain-band areal coverage and the center location and intensity of precipitation.     

Representation of the Madden–Julian Oscillation in CAMS-CSM

The Madden–Julian Oscillation(MJO) has a significant impact on global weather and climate and can be used as a predictability resource in extended-term forecasting. We evaluate the ability of the Chinese Academy of Meteorological Sciences Climate System Model(CAMS-CSM) to represent the MJO by using the diagnostic method proposed by the US Climate Variability and Predictability Program(CLIVAR) MJO Working Group(MJOWG). In general,the model simulates some major characteristics of MJO well, such as the seasonality characteristics and geographical dependence, the intensity of intraseasonal variability(ISV), dominant periodicity, propagation characteristics, coherence between outgoing longwave radiation(OLR) and wind, and life cycle of MJO signals. However, there are a few biases in the model when compared with observational/reanalyzed data. These include an overestimate of precipitation in the convergence zone of the North and South Pacific, a slightly weaker eastward propagation, and a shift in the dominant periodicity toward lower frequencies with slower speeds of eastward propagation. The model gives a poor simulation of the northward propagation of MJO in summer and shows less coherence between the MJO convection and wind. The role of moistening in the planetary boundary layer(PBL) in the eastward/northward propagation of MJO was also explored. An accurate representation of the vertical titling structure of moisture anomalies in CAMS-CSM leads to moistening of the PBL ahead of convection, which accounts for the eastward/northward propagation of MJO. Poor simulation of the vertical structure of the wind and moisture anomalies in the western Pacific leads to a poor simulation of the northward propagation of MJO in this area. Budget analysis of the PBL integral moisture anomalies shows that the model gives a good simulation of the moisture charging process ahead of MJO convection and that the zonal advection of moisture convergence term has a primary role in the detour of MJO over the Maritime Continent.     

Effects of Wind Fences on the Wind Environment around Jang Bogo Antarctic Research Station

This study investigated the flow characteristics altered by Jang Bogo Antarctic Research Station using computational fluid dynamics(CFD) modeling. The topography and buildings around Jang Bogo Station were constructed with computeraided-design data in the CFD model domain. We simulated 16 cases with different inflow directions, and compared the flow characteristics with and without Jang Bogo Station for each inflow direction. The wind data recorded by the site's automatic weather station(AWS) were used for comparison. Wind rose analysis showed that the wind speed and direction after the construction of Jang Bogo Station were quite different from those before construction. We also investigated how virtual wind fences would modify the flow patterns, changing the distance of the fence from the station as well as the porosity of the fence. For westerly inflows, when the AWS was downwind of Jang Bogo Station, the decrease in wind speed was maximized(-81% for west-northwesterly). The wind speed reduction was also greater as the distance of the fence was closer to Jang Bogo Station. With the same distance, the fence with medium porosity(25%–33%) maximized the wind speed reduction.These results suggest that the location and material of the wind fence should be selected carefully, or AWS data should be interpreted cautiously, for particular prevailing wind directions.     

Numerical and Experimental Studies on Flow and Pollutant Dispersion in Urban Street Canyons

In this study numerical simulations and water tank experiments were used to investigate the flow and pollutant dispersion in an urban street canyon. Two types of canyon geometry were tested. The studies indicate that in a step-up notch canyon （higher buildings on the downstream side of the canyon）, the height and shape of the upstream lower buildings plays an important role in flow pattern and pollutant dispersion, while in a step-down notch canyon （lower buildings on the downstream side）, the downstream lower buildings have little influence. The studies also show that the substitution of tall towers for parailelepiped buildings on one side of the canyon may enhance the street ventilation and decrease the pollutant concentration emitted by motor vehicles.     

Air,Water and Soil Environmental Pollution and Regulation Mechanisms in Beijing and Its Ambient Area

1 Physical and Chemical Structural Characteristics of Urban ABLObtained preliminarily the atmospheric dynamic structure of the urban ABL: Urban heat island and urban wind, the zero-plane displacement and roughness of urban building blocks, the distribution and dynamic and thermal impacts of urban fog. In the aspect of city-scale simulation research, the effect of the topography around the outskirts of Beijing on the urban wind field and the effect of buildings on the urban atmospheric circulation are studied. The possible influencing mechanisms of the U-shaped mountains to the west and north of Beijing on the daily variation of the urban wind field, and the formation, transport and cumulation of urban pollutants are analyzed and some distribution charts are obtained.     

Move a Tropical Cyclone with 4D-Var and Vortex Dynamical Initialization in WRF Model

Previous studies showed that 4 D-Var technique used for data assimilation could be modified for weather control. This study demonstrates the ability of 4 D-Var to influence the future path of a tropical cyclone by calculating perturbations in WRF simulation. Given the background error covariance matrix, the initial field is improved by the vortex dynamic initialization technique. Our results show that 4 D-Var can be applied to control the trajectory of simulated tropical cyclones by producing "optimal" perturbations. In the numerical simulation experiment of Typhoon Mitag in 2019, after this kind of weather control similar to data assimilation, the tropical cyclone moved obviously,and the damaging wind over the coastline weakened. The prediction results after the initial field modified by 4 D-Var have a great change, and the position of the tropical cyclone moved about 0.5° southeastward after assimilation,which misses the southeast coast of China. Moreover, the damaging wind is also weakened. Since the 4 D-Var is premised on the assumption that the model is perfect and does not consider the model error, then the research plan to consider model error and introduce new methods is discussed in the paper.     

Evaluation of Weather Research and Forecasting Model Parameterizations under Sea-Breeze Conditions in a North Sea Coastal Environment

Three atmospheric boundary layer(ABL) schemes and two land surface models that are used in the Weather Research and Forecasting(WRF) model, version 3.4.1, were evaluated with numerical simulations by using data from the north coast of France(Dunkerque). The ABL schemes YSU(Yonsei University),ACM2(Asymmetric Convective Model version 2), and MYJ(Mellor–Yamada–Janjic) were combined with two land surface models, Noah and RUC(Rapid Update Cycle), in order to determine the performances under sea-breeze conditions. Particular attention is given in the determination of the thermal internal boundary layer(TIBL), which is very important in air pollution scenarios. The other physics parameterizations used in the model were consistent for all simulations. The predictions of the sea-breeze dynamics output from the WRF model were compared with observations taken from sonic detection and ranging, light detection and ranging systems and a meteorological surface station to verify that the model had reasonable accuracy in predicting the behavior of local circulations. The temporal comparisons of the vertical and horizontal wind speeds and wind directions predicted by the WRF model showed that all runs detected the passage of the sea-breeze front. However, except for the combination of MYJ and Noah, all runs had a time delay compared with the frontal passage measured by the instruments. The proposed study shows that the synoptic wind attenuated the intensity and penetration of the sea breeze. This provided changes in the vertical mixing in a short period of time and on soil temperature that could not be detected by the WRF model simulations with the computational grid used. Additionally, among the tested schemes, the combination of the localclosure MYJ scheme with the land surface Noah scheme was able to produce the most accurate ABL height compared with observations, and it was also able to capture the TIBL.     

Assimilating surface observations in a four-dimensional variational Doppler radar data assimilation system to improve the analysis and forecast of a squall line case

This paper examines how assimilating surface observations can improve the analysis and forecast ability of a fourdimensional Variational Doppler Radar Analysis System(VDRAS).Observed surface temperature and winds are assimilated together with radar radial velocity and reflectivity into a convection-permitting model using the VDRAS four-dimensional variational(4DVAR) data assimilation system.A squall-line case observed during a field campaign is selected to investigate the performance of the technique.A single observation experiment shows that assimilating surface observations can influence the analyzed fields in both the horizontal and vertical directions.The surface-based cold pool,divergence and gust front of the squall line are all strengthened through the assimilation of the single surface observation.Three experiments—assimilating radar data only,assimilating radar data with surface data blended in a mesoscale background,and assimilating both radar and surface observations with a 4DVAR cost function—are conducted to examine the impact of the surface data assimilation.Independent surface and wind profiler observations are used for verification.The result shows that the analysis and forecast are improved when surface observations are assimilated in addition to radar observations.It is also shown that the additional surface data can help improve the analysis and forecast at low levels.Surface and low-level features of the squall line—including the surface warm inflow,cold pool,gust front,and low-level wind—are much closer to the observations after assimilating the surface data in VDRAS.     

Summertime Thermally-Induced Circulations over the Lake Nam Co Region of the Tibetan Plateau

Performance of the fifth-generation Pennsylvania State University/National Center for Atmospheric Research Mesoscale Model(MM5) over the Lake Nam Co region of the Tibetan Plateau was evaluated based on the data from five surface observation sites in 2006. The interaction between two thermally-induced circulations(lake breezes and mountain-valley winds) was also investigated. The results show that MM5 could be used to simulate 2-m air temperature; however, MM5 needs improvement in wind field simulation.Two numerical simulations were conducted to study the effect of the lake on the local weather and wind system. The original land cover of the model was used in the control experiment, and the lake was replaced with grassland resembling the area surrounding the lake in the sensitive experiment. The results of the simulations indicate that the lake enhanced the north slope mountain-valley wind and the mountain changed the offshore flow direction at the north shore. During the day, a clear convergent zone and a strong upflow were observed over the north slope of the Nyainq?entanglha Range, which may cause frequent precipitation over the north slope. During the night, the entire area was controlled by a south flow.