Observations of Lake-Breeze Events During the Toronto 2015 Pan-American Games

Enhanced meteorological observations were made during the 2015 Pan and Parapan American Games in Toronto in order to measure the vertical and horizontal structure of lake-breeze events. Two scanning Doppler lidars (one fixed and one mobile), a C-band radar, and a network including 53 surface meteorological stations (mesonet) provided pressure, temperature, humidity, and wind speed and direction measurements over Lake Ontario and urban areas. These observations captured the full evolution (prior, during, and after) of 27 lake-breeze events (73% of observation days) in order to characterize the convective and dynamic processes driving lake breezes at the local scale and mesoscale. The dominant signal of a passing lake-breeze front (LBF) was an increase in dew-point temperature of \(2.3 \pm 0.3 \,^{\circ }\hbox {C}\), coinciding with a \(180^{\circ }\) shift in wind direction and a decrease in air temperature of \(2.1 \pm 0.2 \,^{\circ }\hbox {C}\). Doppler lidar observations over the lake detected lake breezes 1 hour (on average) before detection by radar and mesonet. On days with the synoptic flow in the offshore direction, the lidars observed wedge-shaped LBFs with shallow depths, which inhibited the radar’s ability to detect the lake breeze. The LBF’s ground speed and inland penetration distance were found to be well-correlated (\(r = 0.78\)), with larger inland penetration distances occurring on days with non-opposing (non-offshore) synoptic flow. The observed enhanced vertical motion \(({>} 1\hbox { m s}^{-1})\) at the LBF, observed by the lidar on 54% of lake-breeze days, was greater (at times \({>} 2.5\hbox { m s}^{-1}\)) than that observed in previous studies and longer-lasting over the lake than over land. The weaker and less pronounced lake-breeze structure over land is illustrated in two case studies highlighting the lifetime of the lake-breeze circulation and the impact of propagation distance on lake-breeze intensity.     

The Impact of the Afternoon Planetary Boundary-Layer Height on the Diurnal Cycle of CO and $$\hbox {CO}_{2}$$ Mixing Ratios at a Low-Altitude Mountaintop

Mountaintop trace-gas mixing ratios are often assumed to represent free atmospheric values, but are affected by valley planetary boundary-layer (PBL) air at certain times. We hypothesize that the afternoon valley–PBL height relative to the ridgetop is important in the diurnal cycle of mountaintop trace-gas mixing ratios. To investigate this, we use, (1) 4-years (1 January 2009–31 December 2012) of CO and \(\hbox {CO}_{2}\) mixing-ratio measurements and supporting meteorological observations from Pinnacles (\(38.61^{\circ }\hbox {N}\), \(78.35^{\circ }\hbox {W}\), 1017 m a.s.l.), which is a monitoring site in the Appalachian Mountains, (2) regional \(\hbox {O}_{3}\) mixing-ratio measurements, and (3) PBL heights determined from a nearby sounding station. Results reveal that the amplitudes of the diurnal cycles of CO and \(\hbox {CO}_{2}\) mixing ratios vary as a function of the daytime maximum valley–PBL height relative to the ridgetop. The mean diurnal cycle for the subset of days when the afternoon valley–PBL height is at least 400 m below the ridgetop shows a daytime CO mixing-ratio increase, implying the transport of PBL air from the valley to the mountaintop. During the daytime, on days when the PBL heights exceed the mountaintop, PBL dilution and entrainment cause CO mixing ratios to decrease. This decrease in CO mixing ratio, especially on days when PBL heights are at least 400 m above the ridgetop, suggests that measurements from these days can be used as with afternoon measurements from flat terrain in applications requiring regionally-representative measurements.     

Recent changes in the precipitation-driving processes over the southern tropical Andes/western Amazon

Climate Dynamics - Analyzing December–February (DJF) precipitation in the southern tropical Andes—STA ($$12^{\circ }\,\hbox {S}$$–$$20^{\circ }\,\hbox {S}$$; > 3000...     

Wake Response to an Ocean-Feedback Mechanism: Madeira Island Case Study

We focus on an island wake episode that occurred in the Madeira Archipelago region of the north-east Atlantic at $32.5^{\circ }\mathrm{N}, 17^{\circ }\mathrm{W}$ . The Weather Research and Forecasting numerical model was used in a (one-way) downscaling mode, considering initial and boundary conditions from the European Centre for Medium-range Weather Forecasts system. The current literature emphasizes adiabatic effects on the dynamical aspects of atmospheric wakes. Changes in mountain height and consequently its relation to the atmospheric inversion layer should explain the shift in wake regimes, from a ‘strong-wake’ to ‘weak-wake’ scenario. Nevertheless, changes in sea-surface temperature variability in the lee of an island can induce similar regime shifts because of exposure to stronger solar radiation. Increase in evaporation contributes to the enhancement of convection and thus to the uplift of the stratified atmospheric layer above the critical height, with subsequent internal gravity wave activity.     

The Effects of Heat Advection on UK Weather and Climate Observations in the Vicinity of Small Urbanized Areas

Weather and climate networks traditionally follow rigorous siting guidelines, with individual stations located away from frost hollows, trees or urban areas. However, the diverse nature of the UK landscape suggests that the feasibility of siting stations that are truly representative of regional climate and free from distorting local effects is increasingly difficult. Whilst the urban heat island is a well-studied phenomenon and usually accounted for, the effect of warm urban air advected downwind is rarely considered, particularly at rural stations adjacent to urban areas. Until recently, urban heat advection (UHA) was viewed as an urban boundary-layer process through the formation of an urban plume that rises above the surface as it is advected. However, these dynamic UHA effects are shown to also have an impact on surface observations. Results show a significant difference in temperatures anomalies (\(p\,< \,0.001\)) between observations taken downwind of urban and rural areas. For example, urban heat advection from small urbanized areas (\(\sim \)1\(\,\hbox {km}^{2}\)) under low cloud cover and wind speeds of 2–3\(\,\hbox {m}\,\hbox {s}^{-1}\) is found to increase mean nocturnal air temperatures by 0.6\(\,^{\circ }\hbox {C}\) at a horizontal distance of 0.5 km. Fundamentally, these UHA results highlight the importance of careful interpretation of long-term temperature data taken near small urban areas.     

Effect of Urbanization on Land-Surface Temperature at an Urban Climate Station in North China

While the land-surface temperature (LST) observed at meteorological stations has significantly increased over the previous few decades, it is still unclear to what extent urbanization has affected these positive trends. Based on the LST data recorded at an urban station in Shijiazhuang in North China, and two rural meteorological stations, the effect of urbanization at the Shijiazhuang station for the period 1965–2012 is examined. We find, (1) a statistically-significant linear trend in annual mean urban–rural LST difference of \(0.27\,^{\circ }\hbox {C}\) \(\hbox {(10 year)}^{-1}\), with an urbanization contribution of 100% indicating that the increase in the annual mean LST at the urban station is entirely caused by urbanization. The urbanization effects in spring, summer and autumn on the trends of mean LST are also significant; (2) the urbanization effect is small for time series of the annual mean minimum LST, and statistically marginal for the trend in annual mean maximum LST [\(0.19\,^{\circ }\hbox {C}\,\hbox {(10 year)}^{-1}\)]; (3) the urbanization effect on the annual mean diurnal LST range (\(\Delta {LST}\)) at the urban station is a strongly significant trend of \(0.23\,^{\circ }\hbox {C (10\,year)}^{-1}\), with an urbanization contribution of 21%. The urbanization effects on trends in the spring and autumn mean \(\Delta {LST}\) are also larger and more significant than for the other seasons; (4) the urbanization effects on the long-term LST trends are remarkably different from those on the near-surface air temperature at the same urban station. Nonetheless, the significant warming of the urban boundary layer is expected to affect the urban environment and ecosystems. However, the problem of data representativeness at an urban station for the monitoring and investigation of large-scale climate change remains.     

Simulations of a Cold-Air Pool in Utah’s Salt Lake Valley: Sensitivity to Land Use and Snow Cover

Obtaining realistic land-surface states for initial and boundary conditions is important for the numerical weather prediction of many atmospheric phenomena. Here we investigate model sensitivity to land use and snow cover for a persistent wintertime cold-air pool in northern Utah during 1–8 January 2011. A Weather Research and Forecast model simulation using the 1993 United States Geological Survey land-use and North American Mesoscale model reanalysis snow-cover datasets is compared to an improved configuration using the modified 2011 National Land Cover Database and a more realistic representation of snow cover. The improved surface specification results in an increase (decrease) in urban land cover (Great Salt Lake surface area), and changes to the snow-cover initialization, depth, extent, and albedo. The results obtained from the model simulations are compared to observations collected during the Persistent Cold-Air Pool Study. The changes in land use and snow cover and the resulting impacts on the surface albedo and surface heat fluxes contributed to near-surface air temperature increases of 1–\(2\,^{\circ }\hbox {C}\) in urban areas and decreases of 2–\(4\,^{\circ }\hbox {C}\) in areas surrounding the Great Salt Lake. Although wind speeds in the boundary layer were overestimated in both simulations, shallow thermally-driven and terrain-forced flows were generally lessened in intensity and breadth in response to the decreased areal extent of the Great Salt Lake and increases in the urban footprint.     

Estimate of Boundary-Layer Depth Over Beijing,China, Using Doppler Lidar Data During SURF-2015

Planetary boundary-layer (PBL) structure was investigated using observations from a Doppler lidar and the 325-m Institute of Atmospheric Physics (IAP) meteorological tower in the centre of Beijing during the summer 2015 Study of Urban-impacts on Rainfall and Fog/haze (SURF-2015) field campaign. Using six fair-weather days of lidar and tower data under clear to cloudy skies, we evaluate the ability of the Doppler lidar to probe the urban boundary-layer structure, and then propose a composite method for estimating the diurnal cycle of the PBL depth using the Doppler lidar. For the convective boundary layer (CBL), a threshold method using vertical velocity variance \((\sigma _w^2 >0.1\,\hbox {m}^{2}\hbox {s}^{-2})\) is used, since it provides more reliable CBL depths than a conventional maximum wind-shear method. The nocturnal boundary-layer (NBL) depth is defined as the height at which \(\sigma _w^2\) decreases to 10 % of its near-surface maximum minus a background variance. The PBL depths determined by combining these methods have average values ranging from \(\approx \)270 to \(\approx \)1500 m for the six days, with the greatest maximum depths associated with clear skies. Release of stored and anthropogenic heat contributes to the maintenance of turbulence until late evening, keeping the NBL near-neutral and deeper at night than would be expected over a natural surface. The NBL typically becomes more shallow with time, but grows in the presence of low-level nocturnal jets. While current results are promising, data over a broader range of conditions are needed to fully develop our PBL-depth algorithms.     

Quantifying and Modelling the Effect of Cloud Shadows on the Surface Irradiance at Tropical and Midlatitude Forests

Cloud shadows lead to alternating light and dark periods at the surface, with the most abrupt changes occurring in the presence of low-level forced cumulus clouds. We examine multiyear irradiance time series observed at a research tower in a midlatitude mixed deciduous forest (Harvard Forest, Massachusetts, USA: \(42.53{^{\circ }}\hbox {N}\), \(72.17{^{\circ }}\hbox {W}\)) and one made at a similar tower in a tropical rain forest (Tapajós National Forest, Pará, Brazil: \(2.86{^{\circ }}\hbox {S}\), \(54.96{^{\circ }}\hbox {W}\)). We link the durations of these periods statistically to conventional meteorological reports of sky type and cloud height at the two forests and present a method to synthesize the surface irradiance time series from sky-type information. Four classes of events describing distinct sequential irradiance changes at the transition from cloud shadow and direct sunlight are identified: sharp-to-sharp, slow-to-slow, sharp-to-slow, and slow-to-sharp. Lognormal and the Weibull statistical distributions distinguish among cloudy-sky types. Observers’ qualitative reports of ‘scattered’ and ‘broken’ clouds are quantitatively distinguished by a threshold value of the ratio of mean clear to cloudy period durations. Generated synthetic time series based on these statistics adequately simulate the temporal “radiative forcing” linked to sky type. Our results offer a quantitative way to connect the conventional meteorological sky type to the time series of irradiance experienced at the surface.     

The Effect of Breaking Waves on $$\hbox {CO}_2$$ Air–Sea Fluxes in the Coastal Zone

The influence of wave-associated parameters controlling turbulent \(\hbox {CO}_2\) fluxes through the air–sea interface is investigated in a coastal region. A full year of high-quality data of direct estimates of air–sea \(\hbox {CO}_2\) fluxes based on eddy-covariance measurements is presented. The study area located in Todos Santos Bay, Baja California, Mexico, is a net sink of \(\hbox {CO}_2\) with a mean flux of \(-1.3\, \upmu \hbox {mol m}^{-2}\hbox {s}^{-1}\) (\(-41.6\hbox { mol m}^{-2}\hbox {yr}^{-1}\)). The results of a quantile-regression analysis computed between the \(\hbox {CO}_2\) flux and, (1) wind speed, (2) significant wave height, (3) wave steepness, and (4) water temperature, suggest that the significant wave height is the most correlated parameter with the magnitude of the flux but the behaviour of the relation varies along the probability distribution function, with the slopes of the regression lines presenting both positive and negative values. These results imply that the presence of surface waves in coastal areas is the key factor that promotes the increase of the flux from and into the ocean. Further analysis suggests that the local characteristics of the aqueous and atmospheric layers might determine the direction of the flux.     

One Year of Surface-Based Temperature Inversions at Dome C, Antarctica

In 2005 the Study of Stable Boundary Layer Environment at Dome C (STABLEDC) experimental campaign was conducted at the plateau station of Concordia at Dome C, Antarctica. Temperature profiles measured with a microwave radiometer were used to study the characteristics of surface-based temperature inversions over the course of a year. Statistics of temperature profiles for every month are discussed; the difference between daytime and nocturnal cases observed during the summer months disappears during winter. Surface-based temperature inversions occurred in 70 % of the time during summer, and almost all of the time during winter. During winter the occurrence of warming events leads to a decrease in the temperature difference between the top and the base of the inversion (i.e. the inversion strength). The inversion strength maxima ranged between $3\,^{\circ }\mathrm{C}$ (December) and $35\,^{\circ }\mathrm{C}$ (August) corresponding to gradients of 0.1 and $0.3\,^{\circ }\mathrm{C}\, \mathrm{m}^{-1}$ , respectively. The average surface-based inversion height presents a daily cycle during the summer months with values up to 200 m in the morning hours, while it affects a layer always deeper than 100 m during the winter months. The relationships between inversion strength and the downward longwave radiative flux, absolute temperature, and wind speed are examined. The inversion strength decreases as the longwave radiation increases. A clear anti-correlation between inversion strength and near-surface temperature is evident throughout the year. During the winter, the largest inversion strength values were observed under low wind-speed conditions; in contrast, a clear dependence was not found during the summer.     

From Near-Neutral to Strongly Stratified: Adequately Modelling the Clear-Sky Nocturnal Boundary Layer at Cabauw

The performance of an atmospheric single-column model (SCM) is studied systematically for stably-stratified conditions. To this end, 11 years (2005–2015) of daily SCM simulations were compared to observations from the Cabauw observatory, The Netherlands. Each individual clear-sky night was classified in terms of the ambient geostrophic wind speed with a \(1\hbox { m} \hbox { s}^{-1}\) bin-width. Nights with overcast conditions were filtered out by selecting only those nights with an average net radiation of less than \(-\,30\hbox { W }\hbox {m}^{-2}\). A similar procedure was applied to the observational dataset. A comparison of observed and modelled ensemble-averaged profiles of wind speed and potential temperature and time series of turbulent fluxes showed that the model represents the dynamics of the nocturnal boundary layer (NBL) at Cabauw very well for a broad range of mechanical forcing conditions. No obvious difference in model performance was found between near-neutral and strongly-stratified conditions. Furthermore, observed NBL regime transitions are represented in a natural way. The reference model version performs much better than a model version that applies excessive vertical mixing as is done in several (global) operational models. Model sensitivity runs showed that for weak-wind conditions the inversion strength depends much more on details of the land-atmosphere coupling than on the turbulent mixing. The presented results indicate that in principle the physical parametrizations of large-scale atmospheric models are sufficiently equipped for modelling stably-stratified conditions for a wide range of forcing conditions.     

An Investigation on Temperature Variance Scaling in the Atmospheric Surface Layer

We present surface-layer measurements of temperature fluctuation variance from a site characterized by small-scale inhomogeneities. Periods of marked radiative forcing are selected. The data characterized by diabatic conditions and vertical heat flux larger than some threshold (here, chosen to be 0.01 K ms⁻¹) agree quite well with the convective scaling in unstable cases, and with the z-less parameterisation (with a large scatter) in stable cases. For near-neutral cases, the similarity function diverges because of the loss of significance of the temperature scale. Departures from similarity are highlighted in cases with smaller thermal fluxes, because horizontal heterogeneity and unsteadiness become important as production terms.     

The Growth of the Planetary Boundary Layer at a Coastal Site: a Case Study

A lidar system is used to determine the diurnal evolution of the planetary boundary layer (PBL) height on a summer day characterized by anticyclonic conditions. The site is located some 15 km distant from the sea, on a peninsula in south-east Italy. Contrary to expectations, the PBL height, after an initial growth consequent to sunrise, ceases to increase about 2 h before noon and then decreases and stabilizes in the afternoon. An interpretation of such anomalous behaviour is provided in terms of trajectories of air parcels towards the lidar site, which are influenced by the sea breeze, leading to a transition from a continental boundary layer to a coastal internal boundary layer. The results are analyzed using mesoscale numerical model simulations and a simple model that allows for a more direct interpretation of experimental results.     

Land-Surface Heterogeneity Effects in the Planetary Boundary Layer

We investigate the cumulative added value of assimilating temperature, moisture, and wind observations in the three-dimensional non-hydrostatic Fifth-Generation Pennsylvania State University/National Center for Atmospheric Research Mesoscale Model MM5 and use these forecasts to analyze the relationship between surface forcing and planetary boundary-layer (PBL) depth. A data assimilation methodology focused on the surface and the PBL, previously tested in a one-dimensional version of MM5, is applied to 29 May, 6 June, and 7 June 2002 during the International $\hbox {H}_{2}\hbox {O}$ Project over the Southern Great Plains. Model-predicted PBL depth is evaluated against PBL depth diagnosed from data across 4,800 km of airborne lidar data (flight tracks 100–300 km long). The forecast with data assimilation verifies better against observations and is thus used to investigate the environmental conditions that govern PBL depth. The spatial structure in PBL depth is found to be most affected by spatial variations in surface buoyancy flux and capping inversion strength. The spatial scales of surface flux forcing reflected in the PBL depth are found through Fourier analysis and multiresolution decomposition. Correlations are ${<}0.50$ at scales of 64 km or less and increase at larger scales for 29 May and 6 June, but on 7 June low correlations are found at all scales, possibly due to greater within-PBL wind speeds, a stronger capping inversion on this day, and clouds. The results suggest a minimum scale, a function of wind speed, below which heterogeneity in surface buoyancy fluxes is not reflected directly in PBL depth.     

Investigation of Turbulence Parametrization Schemes with Reference to the Atmospheric Boundary Layer Over the Aegean Sea During Etesian Winds

The spatial structure of the marine atmospheric boundary layer (MABL) over the Aegean Sea is investigated using the Weather Research and Forecasting (WRF) mesoscale model. Two ‘first-order’ non-local and five ‘1.5-order’ local planetary boundary-layer (PBL) parametrization schemes are used. The predictions from the WRF model are evaluated against airborne observations obtained by the UK Facility for Airborne Atmospheric Measurements BAe-14 research aircraft during the Aegean-GAME field campaign. Statistical analysis shows good agreement between measurements and simulations especially at low altitude. Despite the differences between the predicted and measured wind speeds, they reach an agreement index of 0.76. The simulated wind-speed fields close to the surface differ substantially among the schemes (maximum values range from 13 to \(18\hbox { m s}^{-1}\) at 150-m height), but the differences become marginal at higher levels. In contrast, all schemes show similar spatial variation patterns in potential temperature fields. A warmer (1–2 K) and drier (2–3\(\hbox { g kg}^{-1})\) layer than is observed, is predicted by almost all schemes under stable conditions (eastern Aegean Sea), whereas a cooler (up to 2 K) and moister (1–2\(\hbox { g kg}^{-1})\) layer is simulated under near-neutral to nearly unstable conditions (western Aegean Sea). Almost all schemes reproduce the vertical structure of the PBL and the shallow MABL (up to 300 m) well, including the low-level jet in the eastern Aegean Sea, with non-local schemes being closer to observations. The simulated PBL depths diverge (up to 500 m) due to the different criteria applied by the schemes for their calculation. Under stable conditions, the observed MABL depth corresponds to the height above the sea surface where the simulated eddy viscosity reaches a minimum; under neutral to slightly unstable conditions this is close to the top of the simulated entrainment layer. The observed sensible heat fluxes vary from ?40 to \(25\hbox { W m}^{-2}\), while the simulated fluxes range from ?40 to \(40\hbox { W m}^{-2}\); however, all of the schemes’ predictions are close to the observations under unstable conditions. Finally, all schemes overestimate the friction velocity, although the simulated range (from 0.2 to \(0.5\hbox { m s}^{-1})\) is narrower than that observed (from 0.1 to \(0.7\hbox { m s}^{-1})\).     

Characteristics and Mechanisms of Low-Level Jets in the Yangtze River Delta of China

A dataset obtained using a wind-profile radar located at the Yangtze River Delta in China ( $31.14^{\circ }$ N, $121.81^{\circ }$ E) in 2009 was used to investigate the characteristics and evolution of low-level jets (LLJs) along the east China coast. The study investigated the daily and seasonal structures of LLJs as well as several possible causes. A total of 1,407 1-h LLJ periods were detected based on an adaptive definition that enabled determination of four LLJ categories. The majority (77 %) of LLJs were found to have speeds $<$ 14.0 m s $^{-1}$ (maximum of 34.6 m s $^{-1})$ and occur at an average altitude below 600 m (76 % of the observed LLJs). The dominant direction of the LLJs was from the south-south-west, which accounted for nearly 32 %, with the second most common wind direction ranging from $040^{\circ }$ to $100^{\circ }$ , albeit with a number of stronger LLJs from the west-south-west. A comparison of LLJs and South-west Jets revealed that the frequencies of occurrence in summer are totally different. Results also revealed that in spring and summer, most LLJs originate from the south-south-west, whereas in autumn and winter, north-east is the dominant direction of origin. The peak heights of LLJs tended to be higher in winter than in other seasons. The horizontal wind speed and peak height of the LLJs displayed pronounced diurnal cycles. The Hilbert–Huang transform technique was applied to demonstrate that the intrinsic mode functions with a cycle of nearly 23 h at levels below 800 m, and the instantaneous amplitudes of inertial events (0.0417–0.0476 h $^{-1}$ frequencies) have large values at 300–600 m. The variations in the occurrences of LLJs suggested connections between the formation mechanisms of LLJs and the South-west Jet stream, steady occupation of synoptic-scale pressure system, and land–sea temperature contrasts.     

Modeling of urban heat island and its impacts on thermal circulations in the Beijing–Tianjin–Hebei region,China

Through regulating the land–atmosphere energy balance, urbanization plays an important role in modifying local circulations and cross-border transport of air pollutants. The Beijing–Tianjin–Hebei (BTH) metropolitan area in northern China is frequently influenced by complex atmospheric thermal circulations due to its special topography and geographic position. In this study, the Weather Research and Forecasting (WRF) model combined with remote sensing is used to explore the urbanization impacts on local circulations in the BTH region. The urban heat island (UHI) effect generated around Beijing and Tianjin shows complex interactions with local thermal circulations. Due to the combined effects of UHI and topography, the UHI circulation around Beijing and valley breeze at the southern slopes of Yan Mountain are coupled together to reinforce each other. At the coastal cities, the increased land/sea temperature gradient considerably accelerates the sea breeze along Bohai Bay and moves the sea breeze front further inland to reach as far as Beijing. This study may lay a foundation for the better understanding of air pollutant dispersion on complex terrain.

     

Impact of Megacity Shanghai on the Urban Heat-Island Effects over the Downstream City Kunshan

The impact of upstream urbanization on the enhanced urban heat-island (UHI) effects between Shanghai and Kunshan is investigated by analyzing seven years of surface observations and results from mesoscale model simulations. The observational analysis indicates that, under easterly and westerly winds, the temperature difference between Shanghai and Kunshan increases with wind speed when the wind speed \(<\) 5 m s \(^{-1}\) . The Weather Research and Forecasting (WRF) numerical model, coupled with a one-layer urban canopy model (UCM), is used to examine the UHI structure and upstream effects by replacing the urban surface of Shanghai and/or Kunshan with cropland. The WRF/UCM modelling system is capable of reproducing the surface temperature and wind field reasonably well. The simulated urban canopy wind speed is a better representation of the near-surface wind speed than is the 10-m wind speed at the centre of Shanghai. Without the urban landscape of Shanghai, the surface air temperature over downstream Kunshan would decrease by 0.2–0.4  \(^{\circ }\) C in the afternoon and 0.4–0.6  \(^{\circ }\) C in the evening. In the simulation with the urban landscape of Shanghai, a shallow cold layer is found above the UHI, with a minimum temperature of about \(-0.2\) to \(-\) 0.5  \(^{\circ }\) C during the afternoon hours. Strong horizontal divergence is found in this cold layer. The easterly breeze over Shanghai is strengthened at the surface by strong UHI effects, but weakened at upper levels. With the appearance of the urban landscape specific humidity decreases by 0.5–1 g kg \(^{-1}\) within the urban area because of the waterproof property of an urban surface. On the other hand, the upper-level specific humidity is increased because of water vapour transferred by the strong upward vertical motions.     

A puff model revised by Monte-Carlo method on mesoscale range

A puff model is developed in this study, which simultaneously considers the Monte-Carlo technique, the time and space changes of atmospheric parameters, multiple continuity pollutant sources, linear chemical trans-formation and removal of pollutants, and the effect of complex terrain. The continuously observed turbulent statistical quantities, Lagrangian time scales, mesoscale flow field, and mixing layer depth in the PBL in the Dianchi area in China are directly put into the model, and the diurnal variations of air pollution are forecasted, which are dominated by such mesoscale local circulations as mountain and valley breeze, land and lake breeze, and city heat island (Kunming City). The results show that in the case of inputting the same data, they are in good agreement with the experimental data, as well as with the results of the three-dimensional advection-diffusion model (TD-ADM); the diurnal variation of mesoscale local circulation results in the obvious diurnal variation of mesoscale concentration distribution patterns; the Dianchi lake (appr. 300 km2) has a considerable effect on the distribution of air pollution in the area.