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.

Changes in Wind Speed and Extremes in Beijing during 1960–2008 Based on Homogenized Observations

Daily observations of wind speed at 12 stations in the Greater Beijing Area during 1960–2008 were homogenized using the Multiple Analysis of Series for Homogenization method. The linear trends in the regional mean annual and seasonal （winter, spring, summer and autumn） wind speed series were-0.26,-0.39,-0.30,-0.12 and-0.22 m s-1 （10 yr）-1 , respectively. Winter showed the greatest magnitude in declining wind speed, followed by spring, autumn and summer. The annual and seasonal frequencies of wind speed extremes （days） also decreased, more prominently for winter than for the other seasons. The declining trends in wind speed and extremes were formed mainly by some rapid declines during the 1970s and 1980s. The maximum declining trend in wind speed occurred at Chaoyang （CY）, a station within the central business district （CBD） of Beijing with the highest level of urbanization. The declining trends were in general smaller in magnitude away from the city center, except for the winter case in which the maximum declining trend shifted northeastward to rural Miyun （MY）. The influence of urbanization on the annual wind speed was estimated to be about-0.05 m s-1 （10 yr）-1 during 1960–2008, accounting for around one fifth of the regional mean declining trend. The annual and seasonal geostrophic wind speeds around Beijing, based on daily mean sea level pressure （MSLP） from the ERA-40 reanalysis dataset, also exhibited decreasing trends, coincident with the results from site observations. A comparative analysis of the MSLP fields between 1966–1975 and 1992–2001 suggested that the influences of both the winter and summer monsoons on Beijing were weaker in the more recent of the two decades. It is suggested that the bulk of wind in Beijing is influenced considerably by urbanization, while changes in strong winds or wind speed extremes are prone to large-scale climate change in the region.

Evolution of the Total Lightning Activity in a Leading-Line and Trailing Stratiform Mesoscale Convective System over Beijing

Data from the Beijing SAFIR 3000 lightning detection system and Doppler radar provided some insights into the three-dimensional lightning structure and evolution of a leading-line and trailing-stratiform （LLTS） mesoscale convective system （MCS） over Beijing on 31 July 2007. Most of the lightning in the LLTS-MCS was intracloud （IC） lightning, while the mean ratio of positive cloud-to-ground （＋CG） lightning to –CG lightning was 1：4, which was higher than the average value from previous studies. The majority of CG lightning occurred in the convective region of the radar echo, particularly at the leading edge of the front. Little IC lightning and little ＋CG lightning occurred in the stratiform region. The distribution of the CG lightning indicated that the storm had a tilted dipole structure given the wind shear or the tripole charge structure. During the storm’s development, most of the IC lightning occurred at an altitude of ～9.5 km; the lightning rate reached its maximum at 10.5 km, the altitude of IC lightning in the mature stage of the storm. When the thunderstorm began to dissipate, the altitude of the IC lightning decreased gradually. The spatial distribution of lightning was well correlated with the rainfall on the ground, although the peak value of rainfall appeared 75 min later than the peak lightning rate.

Diurnal Variations of Summer Precipitation in the Beijing Area and the Possible Effect of Topography and Urbanization

The present study examined the diurnal variations of summer precipitation in the Beijing area by usingsubdaily precipitation and wind observations. A combined effect of topography and urbanization on thecharacteristics of diurnal variations was suggested. It was shown that stations located in the plain areaexhibited typical night rain peaks, whereas those in the mountainous area exhibited clear afternoon peaks ofprecipitation diurnal variations. The precipitation peaks were associated with wind fields around the Beijingarea, which were found to be highly modulated by mountain-valley circulation and urban-country circulation.The lower-tropospheric wind exhibited a clear diurnal shift in its direction from north at 0800 LST to southat 2000 LST, which reflected mountain-valley circulation. The transitions from valley to mountain windand the opposite generally happened after sunset and sunrise, respectively, and both occurred earlier for thestations located closer to mountains. By comparing the diurnal variations of precipitation at stations in anortheast suburb, an urban area, and a southwest suburb, it was revealed that the northeast suburb grouphad the highest normalized rainfall frequency, but the southwest group had the lowest from late afternoon tolate evening. On the contrary, in the early morning from about 0200 to 1000 LST, the southwest group andurban group had the highest normalized rainfall frequency. This pattern might originate from the combinedeffects of mountain-valley topography and urbanization.

Changes of Precipitation and Extremes and the Possible Effect of Urbanization in the Beijing Metropolitan Region during 1960–2012 Based on Homogenized Observations

Daily precipitation series at 15 stations in the Beijing metropolitan region（BMR） during 1960–2012 were homogenized using the multiple analysis of series for homogenization method,with additional adjustments based on analysis of empirical cumulative density function（ECDF） regarding climate extremes.The cumulative density functions of daily precipitation series,the trends of annual and seasonal precipitation,and summer extreme events during 1960–2012 in the original and final adjusted series at Beijing station were comparatively analyzed to show the necessity and efficiency of the new method.Results indicate that the ECDF adjustments can improve the homogeneity of high-order moments of daily series and the estimation of climate trends in extremes.The linear trends of the regional-mean annual and seasonal（spring,summer,autumn,and winter）precipitation series are-10.16,4.97,-20.04,5.02,and-0.11 mm（10 yr）-1,respectively.The trends over the BMR increase consistently for spring/autumn and decrease for the whole year/summer; however,the trends for winter decrease in southern parts and increase in northern parts.Urbanization affects local trends of precipitation amount,frequency,and intensity and their geographical patterns.For the urban-influenced sites,urbanization tends to slow down the magnitude of decrease in the precipitation and extreme amount series by approximately-10.4% and-6.0%,respectively; enhance the magnitude of decrease in precipitation frequency series by approximately 5.7%; reduce that of extremes by approximately-8.9%; and promote the decreasing trends in the summer intensity series of both precipitation and extremes by approximately 6.8% and51.5%,respectively.

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 Numerical Study of the Urban Intensity Effect on Fog Evolution in the Beijing-Tianjin-Hebei Region

The influence of urban intensity on fog evolution in the Beijing-Tianjin-Hebei (BTH) region (China) is investigated numerically with the the Weather Research and Forecasting (WRF) model coupled with the urban canopy parameterization-building energy model (UCP- BEM) urban physics scheme. The experiments were designed with a focus on the influence of different urban intensities, which are represented by a different fractional coverage of natural land, buildings, and energy consumption inside buildings in an urban environment. The results of this study indicate that urban areas notably influence fog evolution when natural land is reduced to a small fraction (e.g., less than 10%). Developed land changes fog evolution through urban effects. Higher urban intensity (HUI) generally results in warmer temperatures and lower wind speeds throughout the day, while inhibiting morning specific humidity loss and afternoon specific humidity gain because of the HUI effect on surface heat flux, surface roughness, and surface moisture flux. HUI leads to later and weaker liquid water content formation, with a higher liquid water content base, primarily due to its effect on near surface temperatures. This finding implies that HUI may inhibit the conditions for fog formation. In addition, urban areas with equal natural and developed land coverage seem to greatly enhance the upward surface moisture flux, which is attributed to the combination of a relatively large potential evaporation on developed land and an ample moisture supply from natural land. As a result, the specific humidity increases in the afternoon.

Effect of different cold air intensities and their lagged effects on outpatient visits for respiratory illnesses in Handan in different seasons

本文利用2016年到2019年邯郸市气象要素和呼吸系统疾病门诊数据,分析了不同季节不同强度的冷空气过程及其对呼吸系统疾病的影响,结果显示:尽管呼吸系统疾病在冬季高发,夏季最低,但冷空气对呼吸系统疾病的影响在夏,春季最大,就诊人数分别在冷空气日后两天和五天增加18.4%和13.3%,而冬季就诊人数在冷空气日后三天仅增加3.2%.冷空气对疾病影响的滞后时间在夏,秋和冬季随冷空气强度的增加而减少,而春季的滞后时间总是很长.这些发现可为科学应对气候异常导致的人群健康风险提供针对性依据.

Regional Climate Damage Quantification and Its Impacts on Future Emission Pathways Using the RICE Model

This study quantified the regional damages resulting from temperature and sea level changes using the Regional Integrated of Climate and Economy(RICE)model,as well as the effects of enabling and disabling the climate impact module on future emission pathways.Results highlight varied damages depending on regional economic development and locations.Specifically,China and Africa could suffer the most serious comprehensive damages caused by temperature change and sea level rise,followed by India,other developing Asian countries(OthAsia),and other high-income countries(OHI).The comprehensive damage fractions for China and Africa are projected to be 15.1%and 12.5%of gross domestic product(GDP)in 2195,with corresponding cumulative damages of 124.0 trillion and 87.3 trillion United States dollars(USD)from 2005 to 2195,respectively.Meanwhile,the comprehensive damage fractions in Japan,Eurasia,and Russia are smaller and projected to be lower than 5.6%of GDP in 2195,with cumulative damages of 6.8 trillion,4.2 trillion,and 3.3 trillion USD,respectively.Additionally,coastal regions like Africa,the European Union(EU),and OHI show comparable damages for sea level rise and temperature change.In China,however,sea level-induced damages are projected to exceed those from temperature changes.Moreover,this study indicates that switching the damage modules on or off affects the regional and global emission trajectories,but the magnitude is relatively small.By 2195,global emissions under the experiments with all of the damage modules switched off,only the sea level damage module switched on,and only the temperature damage module switched on,were 3.5%,2.3%and 1.2%higher than those with all of the damage modules switched on,respectively.

On the Ability of the Regional Climate Model RIEMS to Simulate the Present Climate over Asia

A continuous 10-year simulation in Asia for the period of 1 July 1988 to 31 December 1998 was conducted using the Regional Integrated Environmental Model System （RIEMS） with NCEP Reanalysis II data as the driving fields. The model processes include surface physics state package （BATS le）, a Holtslag explicit planetary boundary layer formulation, a Grell cumulus parameterization, and a modified radiation package （CCM3）. Model-produced surface temperature and precipitation axe compared with observations from 1001 meteorology stations distributed over Asia and with the 0.5° × 0.5° CRU gridded dataset. The analysis results show that： （1） RIEMS reproduces well the spatial pattern and the seasonal cycle of surface temperature and precipitation; （2） When regionally averaged, the seasonal mean temperature biases are within 1-2℃. For precipitation, the model tends to give better simulation in winter than in summer, and seasonal precipitation biases are mostly in the range of - 12%-50%; （3） Spatial correlation coefficients between observed and simulated seasonal precipitation are higher in north of the Yangtze River than in the south and higher in winter than in summer; （4） RIEMS can well reproduce the spatial pattern of seasonal mean sea level pressure. In winter, the model-simulated Siberian high is stronger than the observed. In summer, the simulated subtropical high is shifted northwestwards; （5） The temporal evolution of the East Asia summer monsoon rain belt, with steady phases separated by more rapid transitions, is reproduced.

Aerosol Optical Properties Affected by a Strong Dust Storm over Central and Northern China

Aerosol observational data at 8 ground-based observation sites in the Chinese Sun Hazemeter Network （CSHNET） were analyzed to characterize the optical properties of aerosol particles during the strong dust storm of 16-21 April 2005. The observational aerosol optical depth （AOD） increased significantly during this dust storm at sites in Beijing city （86%）, Beijing forest （84%）, Xianghe （13%）, Shapotou （27%）, Shenyang （47%）, Shanghai （23%）, and Jiaozhou Bay （24%）. The API （air pollution index） in Beijing and Tianjin also had a similar rise during the dust storm, while the Angstrhm exponent （a） declined evidently at sites in Beijing city （21%）, Beijing forest （39%）, Xianghe （19%）, Ordos （77%）, Shapotou （50%）, Shanghai （12%）, and Jiaozhou Bay （21%）, respectively. Furthermore, The observational AOD and a demonstrated contrary trends during M1 storm stages （pre-dust storm, dust storm, and post-dust storm）, with the AOD indicating an obvious ＂Valley Peak-Valley＂ pattern of variation, while a demonstrated a ＂Peak-Valley- Peak＂ pattern. In addition, the dust module in a regional climate model （RegCM3） simulated the dust storm occurrence and track accurately and RegCM3 was able to basically simulate the trends in AOD. The simulation results for the North China stations were the best, and the simulation for dust-source stations was on the high side, while the simulation was on the low side for coastal sites.

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.

Variation in Summer Rainfall in North China during the Period 1956-2007 and Links with Atmospheric Circulation

Using gauge precipitation data and NCEP-NCAR reanalysis data, the interdecadal changes in summer precipitation during the period 1956-2007 in North China and the link with atmospheric circulation change over Eurasia are studied. Results show that precipitation amounts decreased by 16.2 mm per decade, which was attributable to a significant reduction in precipitation frequency. Contrary wave trains were found in the subtropical westerly jet （wave guide） over Eurasia for the wet and dry years of North China. When the wave trains had a ridge （trough） around the Korean Peninsula, conditions favored （disfavored） the westward and northward extension of the West Pacific subtropical high. The westward and northward extension of the West Pacific subtropical high is, and was, beneficial to rainfall in North China. The downstream propagation of Rossby waves was found to favor the maintenance of these wave trains. Sensible heating in the south of Lake Baikal and latent heating from the Korean Peninsula to the south of Japan increased during the period 1980-2007, as compared to that during 1957-1979. the wet period. These changes had positive influences on the maintenance of Anticyclonic-Cyclonic anomaly centers in the wave trains. Furthermore, northerly winds were prevalent in the lower troposphere during the dry period （1980-2007）. which prohibited the transportation of water vapor to North China from the seas and thereby led to a decrease in rainfall in North China. The weakening of the Indian Monsoon during the dry period might be one of reasons for the reduction in water vapor transportation.

Numerical simulation of the flow within and over an intersection model with Reynolds-averaged Navier-Stokes method

In this study, the Reynolds-averaged Navier-Stokes （RANS） method is employed to simulate the flow within and over an intersection model with three kinds of k-ε turbulence closure schemes, namely, standard model, renormalization group （RNG） model and realizable k-ε model. The comparison between the simulated and observed flow fields shows that the RANS simulation with all the three turbulence models cannot completely and accurately reproduce the observed flow field in all details. A detailed comparison between the predicted profiles of wind velocities and the measured data shows that the realizble k-ε model is the best one among the three turbulence closure models in general. However, the extent to which the improvement is achieved by the realizable k-ε model is still not enough to completely and accurately describe the turbulent flow in a relatively complex environment.

Updated Homogenized Chinese Temperature Series with Physical Consistency

Most methods of homogenization of climate data are applied to time series of a single variable, such as daily maximum temperature(Tmax) or daily minimum temperature(Tmin). Consequently, the physical relationship among different variables, e.g., Tmax>Tmin, may be distorted after homogenization of climate series of individual variables. The authors develop a solution to improve consistency among diurnal temperature records, while using the Multiple Analysis of Series for Homogenization(MASH) method to homogenize the observation series of daily mean temperature(Tm), Tmin, and Tmax at 545 stations in China for the period 1960–2011, called CHTM2.0. In the previous version of this homogenized dataset based on MASH(CHTM1.0) for the period 1960–2008, there are a few records(0.039% of the total) that are physically inconsistent. For developing CHTM2.0, the authors apply additional adjustments for each day with inconsistent temperature records, in order to hold Tmax>Tm>Tmin. Although the additional adjustments are barely influential for estimating long-term climate trends in China as a whole(because very few records are additionally adjusted), the newly introduced solution improves the physical consistency throughout the dataset. It is also helpful for developing more reasonable homogenized climate datasets with regard to physical consistency among multiple variables. Based on CHTM2.0, the annual Tmax/Tm/Tmin series averaged over China for the period 1960–2011 show significant warming trends of about 0.19/0.25/0.34°C per decade, respectively. Large warming trends of up to 0.425/0.596/ 0.704°C per decade occur in northeastern and northwestern China.

Simulation of Quasi-Linear Mesoscale Convective Systems in Northern China：Lightning Activities and Storm Structure

Two intense quasi-linear mesoscale convective systems（QLMCSs） in northern China were simulated using the WRF（Weather Research and Forecasting） model and the 3D-Var（three-dimensional variational） analysis system of the ARPS（Advanced Regional Prediction System） model.A new method in which the lightning density is calculated using both the precipitation and non-precipitation ice mass was developed to reveal the relationship between the lightning activities and QLMCS structures.Results indicate that,compared with calculating the results using two previous methods,the lightning density calculated using the new method presented in this study is in better accordance with observations.Based on the calculated lightning densities using the new method,it was found that most lightning activity was initiated on the right side and at the front of the QLMCSs,where the surface wind field converged intensely.The CAPE was much stronger ahead of the southeastward progressing QLMCS than to the back it,and their lightning events mainly occurred in regions with a large gradient of CAPE.Comparisons between lightning and non-lightning regions indicated that lightning regions featured more intense ascending motion than non-lightning regions;the vertical ranges of maximum reflectivity between lightning and non-lightning regions were very different;and the ice mixing ratio featured no significant differences between the lightning and non-lightning regions.

The Characteristics and Parameterization of Aerodynamic Roughness Length over Heterogeneous Surfaces

Aerodynamic roughness length （Z0m） is a key factor in surface flux estimations with remote sensing algorithms and/or land surface models. This paper calculates ZOrn over several land surfaces, with 3 years of experimental data from Xiaotangshan. The results show that Z0m is direction-dependent, mainly due to the heterogeneity of the size and spatial distribution of the roughness elements inside the source area along different wind directions. Furthermore, a heuristic parameterization of the aerodynamic roughness length for heterogeneous surfaces is proposed. Individual Z0m over each surface component （patch） is calculated firstly with the characteristic parameters of the roughness elements （vegetation height, leaf area index, etc.）, then Z0m over the whole experimental field is aggregated, using the footprint weighting method.

Locating Parent Lightning Strokes of Sprites Observed over a Mesoscale Convective System in Shandong Province,China

In this paper, we report the location results for the parent lightning strokes of more than 30 red sprites observed over an asymmetric mesoscale convective system(MCS) on 30 July 2015 in Shandong Province, China, with a long-baseline lightning location network of very-low-frequency/low-frequency magnetic field sensors. The results show that almost all of these cloud-to-ground(CG) strokes are produced during the mature stage of the MCS, and are predominantly located in the trailing stratiform region, which is similar to analyses of sprite-productive MCSs in North America and Europe. Comparison between the location results for the sprite-producing CG strokes and those provided by the World Wide Lightning Location Network(WWLLN) indicates that the location accuracy of WWLLN for intense CG strokes in Shandong Province is typically within 10 km, which is consistent with the result based on analysis of 2838 sprite-producing CG strokes in the continental United States. Also, we analyze two cases where some minor lightning discharges in the parent flash of sprites can also be located, providing an approach to confine the thundercloud region tapped by the sprite-producing CG strokes.

Impacts of Two Types of El Nio on Atmospheric Circulation in the Southern Hemisphere

Based on NCEP/NCAR （National Centers for Environmental Prediction/National Center for Atmo- spheric Research） reanalysis data from 1979 to 2010, the impacts of two types of E1 Nino on atmospheric circulation in the Southern Hemisphere （SH） are analyzed. It is shown thaL when a warming event occurs in the equatorial eastern Pacific （EP E1 Nino）, there is a negative sea level pressure （SLP） anomaly in the east- ern Pacific and a positive one in the western Pacific. Besides, there exists a negative anomaly between 40°S and 60°S and a positive anomaly to the south of 60°S. When a warming event in the central Pacific （CP E1 Nino） occurs, there appears a negative SLP anomaly in the central Pacific and a positive SLP anomaly in the eastern and western Pacific, but the SLP anomalies are not so evident in the SH extratropics. In particular, the Pacific-South America （PSA） pattern induced by the CP E1 Nino is located more northwestward, with a weaker anomaly compared with the EP E1 Nino. This difference is directly related with the different position of heating centers associated with the two types of E1 Nino events. Because the SST anomaly associated with CP E1 Nino is located more westward than that associated with EP El Nino, the related heating center tends to move westward and the response of SH atmospheric circulation to the tropical heating changes accordingly, thus exciting a different position of the PSA pattern. It is also noted that the local meridional cell plays a role in the SH high latitudes during EP E1 Nino. The anomalous ascending motion due to the enhancement of convection over the eastern Pacific leads to an enhancement of the local Hadley cell and the meridional cell in the middle and high latitudes, which in turn induces an anomalous descending motion and the related positive anomaly of geopotential height over the Amundsen-Bellingshausen Sea.

THE NORTH-SOUTH ANTI-PHASE DISTRIBUTION OF RAINFALL IN MEIYU PERIODS AND ITS RELATIONSHIP WITH QUASI-BIWEEKLY OSCILLATION IN THE ATMOSPHERE

Based on the daily rainfall datasets from 740 stations in China from 1954 to 2005 and the NCEP/NCAR reanalysis data, the relationship between the north-south anti-phase distribution(APD) of rainfall during Meiyu periods and the Quasi-Biweekly Oscillation(QBWO) in the atmosphere was analyzed. Diagnostic results are as follows:(1) there was significant north-south oscillation of Meiyu rainfall during the 16 years from 1954 to 2005. Since the 1990 s, the APD enhanced significantly and showed 2- and 4-6-year period. In the region with more rainfall, the QBWO was always more active.(2) The APD of Meiyu and north-south movements of precipitation in eastern China belong to the same phase.(3) The 10-25 day filtered water vapor flux could spread to the area north of 30°N in 1991. The divergence of the water vapor flux which propagated from middle- and higher- latitudes to the of Yangtze-Huaihe River Basins(YHRB) was significant in 1991, but the latitudes that the water vapor flux could reach were further southward and there was no southward propagation of divergence in 1993.(4) The locations of Western Pacific Subtropical High(WPSH) and 10-25 day anti-cyclone, which modulated WPSH's advancement in and out of the South China Sea, were relatively northward in 1991. Furthermore, the vertical circulation showed north-south deviation between 1991 and 1993, just as other elements of the circulation did.