Megacity-Induced Mesoclimatic Effects in the Lower Atmosphere : A Modeling Study for Multiple Summers over Moscow , Russia

Urbanization leads to distinct meteorological features of urban environments, and one the best-known is the urban heat island (UHI) effect. For megacities, these features become mesoscale phenomena (scale ≥ 10 km) that are amplified by the tropospheric feedbacks, and have substantial implications on human well-being. For the first time, a three-dimensional statistical description of the megacity-induced meteorological effects extending towards the lower troposphere for summer is acquired on a quasi-climatological timescale (a decade) based on high-resolution (1 km) simulations for Moscow with the COSMO-CLM model with and without its urban canopy model TERRA_URB. Our results confirm the features from previous observational and modeling studies, including the UHI itself, the cooling effect above established by the cross-over effect, the urban dry/moist islands and the urban breeze circulation. Particularly, the UHI shows a strong diurnal variation in terms of intensity and vertical extent between daytime (≈0.5 K/≈1.5 km) and nighttime (>3 K/≈150 m). We have discovered a systematic veering in the downwind shift of the UHI spatial pattern established by the Coriolis effect, and an enhanced stable stratification of the rural surroundings established by the urban plumes further downwind. Finally, extending the analysis to multiple summers demonstrates a substantial increase in summer precipitation (up to +25%) over the city center and its leeward side. These urban-caused mesoclimatic effects need to be taken into account in weather and climate services, including the design of future megacities.


Introduction
Essential features of urban climate are deviations of the temperature and humidity from surrounding rural environment, known as urban heat island (UHI) and urban dry/moisture islands (UDI/UMI).Investigation of these effects and taking them into account within the framework of different weather and climate services is especially important due to their significant environmental and social impacts.In general, such impacts are considered as negative due to urban-caused acceleration of climate warming in the cities [1][2][3] and amplification of the heat waves [4,5] and related heat stress [6,7] affecting the human health and even mortality rates [8][9][10].The drivers of urban-rural temperature and humidity differences and their near-surface climatology are well-studied for a big variety of cities all over the world (see classic studies [11,12] and review papers [13,14]).UDI and UMI effects are described, e.g., in [15][16][17][18].
The micro-climatological features established by the local urban environment (street canyons, parks, etc.) are generally limited to the urban canopy itself spanning in the air up to the roof level.However, for medium-sized urban areas and megacities, influence of surface parameters that differ from surrounding rural areas (in terms of albedo, roughness, thermal structure, evaporative cooling, etc.) is accumulated over the whole urban area reaching spatial extent of tens of kilometers or even more.In turn, these distinct surface features result in meteorological features of temperature, humidity and wind speed at similar spatial scales.According to the classification of atmospheric processes by scale (e.g., [19]), such features can be considered as mesoscale phenomena.
These "urban-caused mesoclimatic phenomena" involve interactions with the troposphere aloft by changing the size and properties of the atmospheric boundary layer (ABL) [20] particularly leading to the "boundary-layer urban heat island" (e.g., [21,22]) and urban-induced mesoscale circulations [23][24][25].They further affect other tropospheric features such as the presence of cumulonimbus clouds and convective precipitation [26][27][28] through the changes in the triggering of moist-convection [29].The urban-caused mesoclimatic phenomena interacting with the troposphere play the key role in the climate of the entire urban agglomerations: they do not only involve the "bottom-up" urban forcing affecting the troposphere that scales with the city extent as described above, but also the "top-down" tropospheric feedbacks that further amplify the distinct climatic conditions in the urban canopy.Examples of the latter are the advection of hot air between urban neighborhoods resulting in an UHI intensification with city sizes [7, 30,31], the down-mixing of hot air due to the cross-over effect [21,32,33], and cascade effects such as the changes in radiation and rain hitting the cities that lead to an additional modulation of the surface energy and water balance.
The number of studies devoted to urban-caused climate phenomena that consider the lower troposphere is much smaller than those focusing only on the near-surface climate features.This could be explained by three reasons: (1) interest in layer in which we live; (2) difficulties of making continuous observations above the roof level; and (3) significantly higher complexity of three-dimensional atmospheric processes in ABL and lower troposphere, especially for coastal and mountain cities.
As such, most existing studies considering the urban-tropospheric interactions listed above are based only on episodic observations (the case studies), which have the insufficient timespan for making conclusions on urban-climate statistics.Moreover, various urban-induced tropospheric effects are studied by different methods and for cities located in different climatic zones and terrain types, which makes the understanding of the overall picture especially difficult.Therefore, quantitative and even qualitative estimates of such urban effects vary strongly.For example, estimations of the vertical extent of the UHI over big cities vary from 100-300 m [21,22,34] to 1-2 km [1,33,35,36] in dependence from the day time and season.The estimates of climatological or seasonal means also vary significantly.Understanding of urban effects on precipitation [28], cloudiness and wind regime is even more complicated.
Recent developments in field of modeling of atmospheric dynamics and computation techniques have opened new opportunities to study these urban-induced mesoscale atmospheric processes and feedbacks with a high detail in a systematic way towards the climatological timescale.Modern limited-area atmospheric models are possible to simulate meteorological regime over wide areas (from hundreds to thousands of kilometers) with horizontal grid spacing of few kilometers and even less, which allows directly resolving most mesoscale processes.In addition, models allow running various sensitivity experiments for in depth investigation of urban-atmospheric interactions.
In this study, we consider the result of application of such modeling approach, based on numerical experiments with COSMO-CLM regional climate model, to study the urban-caused mesoclimatic features of Moscow megacity, which is the biggest monocentric agglomeration in Europe and a promising place for urban climate studies (see Section 2.1).
The aim of this paper is to investigate the influence of a big urban agglomeration on the meteorological regime of the lower troposphere in terms of the temperature, humidity, wind and precipitation within similar methodology and compare results with existing studies based on various observations.In contrast to several similar modeling studies, based on simulations for short periods (few days, e.g.[24,32]), we consider simulations conducted at a quasi-climatological time scale (a decade).Study is focused on summer due to: (1) higher probability of cases with anticyclonic weather and intensive UHI; (2) higher relevance of urban climate studies firstly because of urban-induced heat stress during heat waves [4][5][6]; and (3) existing problems for simulations winter-time meteorological regime within considered modeling framework [37].

The Study Area
Moscow megacity is the biggest monocentric agglomeration in Europe.The population of the administrative unit is about 12 million people, while population of whole urban agglomeration is 16-17 million people [38].The population of Moscow is rapidly increasing: since 1970, it has grown by more than 60%.The area of Moscow (the administrative unit) is ≈2500 km 2 , however it includes a wide area known as "New Moscow", which was joined to the city in 2011, but remains practically unbuilt.The actual area of the city (excluding the suburbs and satellite cities) is about 1000 km 2 .Unlike most European and American cities, Moscow is densely built by the midrise and high-rise block-houses, while private low-rise houses are practically absent.According to the classification of Stewart and Oke [39], the prevailing types of local climate zones in Moscow are 2 (compact midrise), 4 (open high-rise) and 5 (open midrise) [40].
The Köppen Climate Classification subtype for Moscow is "Dfb" (Warm Summer Continental Climate).The mean annual rural temperature (according to observations at nine rural stations, see Figure 1) for recent WMO standard period (1981-2010) is 5.1 • C, mean summer (JJA) temperature is 17.2 • C and mean winter (DJF) temperature is −7.2 • C. The mean annual precipitation amount is about 700 mm [41] with a typical summer maximum.Due to the size and listed climate and building features, Moscow forms an intensive UHI with a mean annual intensity value equal to 2 K (for Balchug weather station at the city center, see Figure 1) and its maximum values up to 12-13 K [41].Such values of mean annual UHI intensity are higher than those for a number of other big megacities including London (≈1 K [42,43]), Beijing (1.5-2 K [44][45][46]), and Shanghai (≈1 K [47]).
The city is located within a flat and homogeneous terrain of the East-European plain, far away from big water bodies and mountain ranges (Figure 1).In addition to the compact and relatively symmetric shape of the city, such terrain features enable us to investigate the pure urban-atmosphere interactions by minimizing the interference of other mesoscale features induced by non-urban surface properties.Finally, Moscow has a unique record of long-term urban-climate observations [1,41] and new networks of detailed intra-urban observations [30,48].All of this makes Moscow a promising and interesting megacity for urban climate studies.

Regional Climate Model
This study is based on numerical simulations with the regional climate model COSMO-CLM (Climate Limited-area Modeling community) [49,50].It is a special version of COSMO model (COnsortium for Small scale MOdeling), a non-hydrostatic mesoscale numerical weather prediction system, which is adapted for long-term numerical experiments.This model was applied for a three-step dynamic downscaling of ERA-Interim reanalysis [51] data (with longitude and latitude grid steps of 0.75 • ).For the first step of downscaling, 6-hourly reanalysis data were used as the initial and boundary conditions for running the model for a "basic" domain D1 covering the central part of European Russia with a size of approximately 1700 km × 1700 km and a 12-km horizontal grid step, 140 × 140 grid cells (Figure 1a).In addition, to ensure a more reliable binding of the internal model mode to the real atmospheric dynamics, the spectral nudging technique [52] was applied, which is important for such studies [53,54].Modeling results for D1 were used as the initial and boundary conditions for intermediate domain D2 with a size of 600 km × 600 km (200 × 200 grid cells) and a 3-km grid spacing covering the Moscow region and the neighboring regions of Russia.Then, the same procedure was repeated for domain D3 with a size of 180 km × 180 km with a 1-km grid spacing.The final 1-km step was selected to resolve urban-atmosphere interactions and the consequential urban-caused mesoscale features.Configuration of the model was optimized to minimize the known problems of overestimation of nocturnal temperatures by COSMO model [48,[55][56][57] and problems of temperature overestimation during summer heat waves reported in [58].Such optimization included the tuning of the turbulent diffusion scheme as suggested in [55], using the new vegetation-skin temperature formulation [59] and other modifications, see the study of Varentsov et al. (2017) [30] for more details.

Urban Canopy Model
The bulk urban canopy model (UCM) called TERRA_URB [60,61] coupled to COSMO-CLM was used to simulate the urban climate.To take into account the urban physics in terms of surface energy exchanges including the influence of street-canyon geometry, TERRA_URB provides corrections of the surface parameters (roughness length, albedo, emissivity, heat capacity, etc.) using the semiempirical urban canopy dependencies [61].Instead of explicitly calculating the urban physical processes (e.g., the reduced-albedo effect from multiple-reflections inside a street canyon), they are implicitly taken into account in the land-surface parameters in accordance to detailed observational studies, modeling experiments and available parameter inventories.This way, TERRA_URB bridges the gap between previous bulk schemes (e.g., URBCLIM [62], Noah/Urban [63] or IUM [64]) and explicit-canyon schemes (e.g., CLM-U [65], BEP [66], TEB [67] or SLUCM [63]), hence provides canopy-dependent urban physics with a low computational cost.
It further parameterizes the moisture balance of impervious surfaces considering a distribution of water puddles, which provides a better representation of the magnitude and duration of evaporation after rainfall for different cities compared to previous water-storage parameterizations [60,68].The anthropogenic heat flux (AHF) is parameterized according to Flanner (2009) [69] taking the annual-mean value into account together with its typical diurnal and annual variation and the dependence of its annual amplitude on the latitude.Modeling results for D1 were used as the initial and boundary conditions for intermediate domain D2 with a size of 600 km × 600 km (200 × 200 grid cells) and a 3-km grid spacing covering the Moscow region and the neighboring regions of Russia.Then, the same procedure was repeated for domain D3 with a size of 180 km × 180 km with a 1-km grid spacing.The final 1-km step was selected to resolve urban-atmosphere interactions and the consequential urban-caused mesoscale features.Configuration of the model was optimized to minimize the known problems of overestimation of nocturnal temperatures by COSMO model [48,[55][56][57] and problems of temperature overestimation during summer heat waves reported in [58].Such optimization included the tuning of the turbulent diffusion scheme as suggested in [55], using the new vegetation-skin temperature formulation [59] and other modifications, see the study of Varentsov et al. (2017) [30] for more details.

Urban Canopy Model
The bulk urban canopy model (UCM) called TERRA_URB [60,61] coupled to COSMO-CLM was used to simulate the urban climate.To take into account the urban physics in terms of surface energy exchanges including the influence of street-canyon geometry, TERRA_URB provides corrections of the surface parameters (roughness length, albedo, emissivity, heat capacity, etc.) using the semi-empirical urban canopy dependencies [61].Instead of explicitly calculating the urban physical processes (e.g., the reduced-albedo effect from multiple-reflections inside a street canyon), they are implicitly taken into account in the land-surface parameters in accordance to detailed observational studies, modeling experiments and available parameter inventories.This way, TERRA_URB bridges the gap between previous bulk schemes (e.g., URBCLIM [62], Noah/Urban [63] or IUM [64]) and explicit-canyon schemes (e.g., CLM-U [65], BEP [66], TEB [67] or SLUCM [63]), hence provides canopy-dependent urban physics with a low computational cost.
It further parameterizes the moisture balance of impervious surfaces considering a distribution of water puddles, which provides a better representation of the magnitude and duration of evaporation after rainfall for different cities compared to previous water-storage parameterizations [60,68].The anthropogenic heat flux (AHF) is parameterized according to Flanner (2009) [69] taking the annual-mean value into account together with its typical diurnal and annual variation and the dependence of its annual amplitude on the latitude.
To describe the heterogeneity of the urban surface, TERRA_URB implements the so-called tile approach for which urbanized and natural parts can coexist in each model grid cell.Calculation of the resulting fluxes from the model cells to the atmosphere accounts the urban fraction parameter (UF), which could be defined individually for each model cell.In addition, TERRA_URB allows defining a heterogeneous distribution of the AHF and urban canopy parameters (building height H, roof area fraction RF and street canyon aspect ratio H/W), which are used for resulting semi-empirical corrections of the surface parameters.
The definition of the urban canopy parameters is an important methodological problem of mesoscale modeling for urban areas.A popular approach, used, e.g., in well-known WRF model [70], is based on the classification of the urban area to a number of land-use classes, each characterized by specific values for urban canopy parameters.In our study, we used another approach.All of the required parameters (UF, AHF, H, RF, and H/W) were calculated individually for each model cell with the original technology of GIS-processing of open-access OpenStreetMaps data [71].This technology converts the elementary map units (buildings, roads segments, land-use polygons, etc.) to specific parameters of the urban canopy, averaged over the given model cells.Partially, it identifies street canyons between the real buildings and calculates their statistical parameters.
For the annual mean AHF, the existing estimate of its average value over urban areas as 75 W/m 2 made by Stewart and Kennedy [72] was used.Following the annual cycle of the AHF considered in TERRA_URB, this gives about 30 W/m 2 for the summer months.Spatial distribution of the AHF was obtained considering the values of UF, H and RF (for more details, see [30]).The spatial distribution of resulting values (UF, AHF, H, RF, and H/W) used in the study is shown in Figure S1; the UF is also shown in Figure 1b.

Configuration of the Numerical Experiments and Timespan of the Study
The TERRA_URB scheme was used only for D3 domain in the "URB" model runs.In addition, "noURB" runs with TERRA_URB switched off were conducted.All numerical experiments were initialized on 1 May and were run until 31 August.The first month (May) was considered as a spin-up period required for the adaptation of the model to the meteorological regime.Such experiments were performed for ten summer seasons 2007-2016.Despite obvious changes of the land-use and urban canopy parameters during this period, their values (calculated for conditions of 2015) were considered unchanged, because the simulations for multiple summers were aimed only to obtain more accurate statistical features of urban-caused meteorological effects over the 10 summer seasons.Hereby, it is considered that the additional urban expansion during the 10-year period is small compared to the urban extent of Moscow itself.Temporal resolution of model output is 1 h.
The detailed analysis of three-dimensional structure of urban-caused meteorological effects is shown for summer 2014, which is characterized by the higher probability of anticyclonic weather with intensive UHI.To enforce the signal, we consider results for the selection of days with pronounced UHI (with daily maximum of canopy-layer UHI intensity >4 K) which gives 78% of all summer days, or for the smaller sub-selections described in the text.The urban effects on cloudiness and precipitation are studied according to simulations for all the 10 summer seasons.

The Definition of Urban-Induced Meteorological Effects and Their Magnitude in This Study
The key strength of the modeling approach to study urban climate features is the opportunity to consider them as the model responses to land-use changes (rural to urban), modulated by switching on and off the UCM.This approach has been used in urban climatology only recently (e.g., [24,32,33]).Accordingly, in this study (except in Section 3.1), we analyze the urban-induced meteorological effects as model responses to switching on the TERRA_URB UCM (hereafter, called just responses), defined as differences between the results of "URB" and "noURB" simulations.For example, to analyze the UHI intensity, we consider temperature response ∆T, defined for any given point in space in time as: Similar values are used to characterize other urban-induced effects in fields of the lapse rate (∆γ), specific humidity (∆q), wind speed and direction (∆V), precipitation amounts (∆P), and cloudiness (∆N).For the conditions of the relatively flat and homogenous landscapes of Moscow region, the observed urban climate features are mostly driven by the influence of the megacity.This is further supported by the good agreement between various observations and modeling results (see below) that could only be reproduced when including the UCM.This justifies that such model responses could represent the urban-caused meteorological effects at least qualitatively.

Observations Used for Model Verification
For the purpose of model verification (Section 3.1), we used the data of meteorological stations located in the Moscow region (Figure 1b).Balchug station (WMO id 27605) is located in the city center, 600 m from the Kremlin [41], and characterizes the conditions of densely build urban environment.Its location is considered as geometric city center for analysis of modeling results in the following sections (from Section 3.2).Another urban station used for verification is meteorological observatory of Lomonosov Moscow State University (MSU), located in the park zone 8 km to the southwest from the city center [73].For rural conditions, we use the data from nine stations surrounding the city (Figure 1b).Even though all of them are marked near the smaller towns, all of them are actually located far away from the built environment and could be considered as representative for rural environment.More details about these station are shown in [30,36,41,48].

Model Verification Based on Canopy-Layer Observations
Verification of the model ability to simulate the key features of investigated phenomena is an essential part of modeling studies.Since the canopy-layer UHI and UDI/UMI intensities are considered as important for most of the studied features of urban climate, model-to-observations comparison was performed for this variables in addition to rural values of near-surface (2 m) air temperature and specific humidity.The observed canopy-layer UHI and UDI/UMI intensities were defined as deviations from the mean temperature and humidity over nine rural stations (for humidity, positive deviation indicates UMI while negative deviation indicates UDI).Such values were calculated separately for the city center (Balchug station) and the urban park (MSU).Modeled values for such comparison were calculated in the same way based on the "URB" simulations for the corresponding model cells, selected among the nine nearest cells around each of the weather stations.For Balchug station, we used the temperature and humidity values that were modeled for the urban tile in the nearest densely built urban cell.For MSU station, we used such values for the natural tile of the nearest semi-urban cell (UF ≈ 50%), and for the rural stations the values for the nearest non-urban cells were used.
Figure 2 shows the model-to-observation comparison for mean rural temperature and UHI intensity for two urban sites, while Figure 3 shows the same for the rural specific humidity and UDI/UMI intensity.Such comparison indicates a good agreement between modeled and observed values in terms of temporal variation (on diurnal and day-to-day timescales), average values and typical extremes, including high UHI intensity values reaching up to 10 K.The only significant bias is a systematic underestimation of the urban-rural moisture difference (an overestimation of the daytime UDI intensity and underestimating of the nocturnal UMI intensity), which is discussed in Section 3.4.
We would like to highlight that such a good quality of replication of observations by the model was achieved using the optimized model configuration described in [30].In addition, that study reports the good agreement between the modeling results and observations for Moscow for other summer seasons in terms of the temperature, UHI intensity and summer precipitation amounts.Additionally, it proves the adequacy of the simulations of the UHI spatial structure based on comparison between modeling results and data of the newest detailed observational networks.
summer seasons in terms of the temperature, UHI intensity and summer precipitation amounts.Additionally, it proves the adequacy of the simulations of the UHI spatial structure based on comparison between modeling results and data of the newest detailed observational networks.

Vertical Structure of the UHI
The vertical cross-sections trough the city center in field of temperature response to switching on the UCM (Figure 4), calculated as described in the Section 2.5, shows that daily-average UHI vertical extent reaches 250-400 m.Daytime (12 UTC/15 MSK) UHI extends up to 1-1.5 km, but its intensity even near the surface does not exceed 0.5 K.In the evening (16 UTC/19 MSK), UHI has similar vertical extent and strengthens within lower 400 m.

Vertical Structure of the UHI
The vertical cross-sections trough the city center in field of temperature response to switching on the UCM (Figure 4), calculated as described in the Section 2.5, shows that daily-average UHI vertical extent reaches 250-400 m.Daytime (12 UTC/15 MSK) UHI extends up to 1-1.5 km, but its intensity even near the surface does not exceed 0.5 K.In the evening (16 UTC/19 MSK), UHI has similar vertical extent and strengthens within lower 400 m.The averaging is made for the days with intensive urban heat island during June-August 2014.
The contour lines for the values higher than 1 K are indicated at each 1 K. Zero at horizontal axis is the city center (Balchug weather station, see Section 2.2).Bold black line shows the urbanized area (urban fraction higher than 25%).The vertical coordinate is elevation above the sea level (a.s.l.), gray shading represents the relief.
Intensive nocturnal (0 UTC/3 MSK) UHI exists only within thin (100-150 m) layer and is overlaid by a layer with negative temperature response that indicates a so-called cross-over effect [21] or cold lens [34].Its appearance could be explained by more intensive vertical mixing and upward motions over the city in conditions of stable stratification [74].The other explanation is the divergence of longwave radiation flux, typical for more moist and polluted urban air over the city [34,75].However, these factors do not play a significant role in the model because urban features of atmospheric composition are not considered, while urban humidity is underestimated by the model (see Sections 3.1 and 3.4).
The existing measurement campaigns for Moscow region makes a qualitative comparison possible between the modeling and the observational results.According to the recent study based on the data of contact measurements (analyzed on annual basis only) [22], daily average value of UHI vertical extent is ≈300 m and nocturnal crossover effect is identified higher than 100 m over the surface.Recent studies, based on the satellite remote sensing data [1,35,36], report that UHI vertical extent at midday and evening hours during the summer is 1.5-2 km, which agrees with modeling results for this time of the day.
Another estimation of UHI vertical extent could be obtained from the measurements of microwave temperature profilers MTP-5 [76,77], operating in Moscow region [34,78].The most relevant in the scope of this study are recent results [79], based on synchronous observations at urban park located close to city center (MSU site) and two suburb sites (Dolgoprudny and Kosino) during the warm season (April-July) of 2015.The average day-time temperature profiles over these sites are very similar in terms of vertical gradients.At the same time, the mean profiles over selection of nocturnal hours with pronounced canopy-layer UHI clearly show its vertical extent up to 200 m (Figure 5), which agrees with modeling results.Nocturnal crossover effect is not identified by these observations, which could be explained by its low magnitude or its spread over the nearest suburbs.
The qualitative agreement between modeling results, listed observations for Moscow megacity and also simulations for Paris within different modeling framework [32] gives us confidence that the COSMO-CLM model coupled to TERRA_URB UCM adequately simulates the UHI in the ABL and the precursor urban-atmospheric interactions.Thus, the additional features of the urban troposphere are investigated below.
city center (Balchug weather station, see Section 2.2).Bold black line shows the urbanized area (urban fraction higher than 25%).The vertical coordinate is elevation above the sea level (a.s.l.), gray shading represents the relief.
Intensive nocturnal (0 UTC/3 MSK) UHI exists only within thin (100-150 m) layer and is overlaid by a layer with negative temperature response that indicates a so-called cross-over effect [21] or cold lens [34].Its appearance could be explained by more intensive vertical mixing and upward motions over the city in conditions of stable stratification [74].The other explanation is the divergence of longwave radiation flux, typical for more moist and polluted urban air over the city [34,75].However, these factors do not play a significant role in the model because urban features of atmospheric composition are not considered, while urban humidity is underestimated by the model (see Sections 3.1 and 3.4).
The existing measurement campaigns for Moscow region makes a qualitative comparison possible between the modeling and the observational results.According to the recent study based on the data of contact measurements (analyzed on annual basis only) [22], daily average value of UHI vertical extent is ≈300 m and nocturnal crossover effect is identified higher than 100 m over the surface.Recent studies, based on the satellite remote sensing data [1,35,36], report that UHI vertical extent at midday and evening hours during the summer is 1.5-2 km, which agrees with modeling results for this time of the day.
Another estimation of UHI vertical extent could be obtained from the measurements of microwave temperature profilers MTP-5 [76,77], operating in Moscow region [34,78].The most relevant in the scope of this study are recent results [79], based on synchronous observations at urban park located close to city center (MSU site) and two suburb sites (Dolgoprudny and Kosino) during the warm season (April-July) of 2015.The average day-time temperature profiles over these sites are very similar in terms of vertical gradients.At the same time, the mean profiles over selection of nocturnal hours with pronounced canopy-layer UHI clearly show its vertical extent up to 200 m (Figure 5), which agrees with modeling results.Nocturnal crossover effect is not identified by these observations, which could be explained by its low magnitude or its spread over the nearest suburbs.
The qualitative agreement between modeling results, listed observations for Moscow megacity and also simulations for Paris within different modeling framework [32] gives us confidence that the COSMO-CLM model coupled to TERRA_URB UCM adequately simulates the UHI in the ABL and the precursor urban-atmospheric interactions.Thus, the additional features of the urban troposphere are investigated below.

Urban Heat Plumes
The modeling results clearly demonstrate urban heat plumes (UHP) [80] downwind of Moscow.Their existence follows from general features of ABL dynamic over inhomogeneous landscapes: warm urban air is carried by the wind, interacts with the colder and more stable suburban internal boundary layer and is pushed upward, which forms the elevated urban plume.While the influence of wind direction on the shape of the near-surface UHI is reported in the number of studies (e.g., [79,81,82]), three-dimensional structure of UHPs is practically unstudied.Probably, the only experimental proofs of their existence are based on the data of helicopter and airplane sounding over American cities for certain case-studies [23,80].
The modeling results allow us to obtain the more detailed vertical structure of the UHI and those UHP under different wind directions.The classification of the cases by wind direction is performed based on the mean wind direction within the lowest 150 m above the ground over the whole D3 domain (four lowest model levels).The results are grouped for prevailing northern (315-45 • ), eastern (45-135 • ), southern (135-225 • ) and western (225-315 • ) wind direction and for each hour.Vertical cross-sections of the temperature response, averages for such selections, shows that the UHI is shifted by an incoming stream to the lee side of the city for the entire day.This effect is most evidently pronounced at night and in the early morning.
At the same time, the formation of elevated UHPs downwind to the city is observed.They are most pronounced between 0 and 4 UTC (3-7 MSK), stretching for tens of km (Figure 6).The spatial structure of the UHI wind shear and UHPs also features the Ekman wind rotation with altitude (Figure S2).For example, the UHI at the lowest model level is shifted to northwest for the selection of nocturnal cases with prevailing northern wind and to southeast for the selection of cases with southern wind.With increasing altitude, the shear direction turns right following the wind rotation.

Urban Heat Plumes
The modeling results clearly demonstrate urban heat plumes (UHP) [80] downwind of Moscow.Their existence follows from general features of ABL dynamic over inhomogeneous landscapes: warm urban air is carried by the wind, interacts with the colder and more stable suburban internal boundary layer and is pushed upward, which forms the elevated urban plume.While the influence of wind direction on the shape of the near-surface UHI is reported in the number of studies (e.g., [79,81,82]), three-dimensional structure of UHPs is practically unstudied.Probably, the only experimental proofs of their existence are based on the data of helicopter and airplane sounding over American cities for certain case-studies [23,80].
The modeling results allow us to obtain the more detailed vertical structure of the UHI and those UHP under different wind directions.The classification of the cases by wind direction is performed based on the mean wind direction within the lowest 150 m above the ground over the whole D3 domain (four lowest model levels).The results are grouped for prevailing northern (315-45°), eastern (45-135°), southern (135-225°) and western (225-315°) wind direction and for each hour.Vertical cross-sections of the temperature response, averages for such selections, shows that the UHI is shifted by an incoming stream to the lee side of the city for the entire day.This effect is most evidently pronounced at night and in the early morning.
At the same time, the formation of elevated UHPs downwind to the city is observed.They are most pronounced between 0 and 4 UTC (3-7 MSK), stretching for tens of km (Figure 6).The spatial structure of the UHI wind shear and UHPs also features the Ekman wind rotation with altitude (Figure S2).For example, the UHI at the lowest model level is shifted to northwest for the selection of nocturnal cases with prevailing northern wind and to southeast for the selection of cases with southern wind.With increasing altitude, the shear direction turns right following the wind rotation.
In addition, the non-trivial effect of the UHPs on thermal stratification over rural areas was revealed.It is well-known that thermal stratification over the cities is in general more unstable than over rural areas due to higher buoyancy of warmer urban air.The modeling results reproduce this pattern: the response of the lapse rate ( = − ⁄ ), calculated between two lowest model levels (10 and 35 m) is positive over urbanized model cells (Figure 7a,b).However, over the nearby rural areas negative responses are observed.In other words, influence of the big city makes the atmosphere more stable over surrounding rural areas.Such model behavior is related to the UHPs.This is clearly illustrated by mean values of lapse rate response over selection of cases with prevailing northern wind (Figure 7c).In addition, the non-trivial effect of the UHPs on thermal stratification over rural areas was revealed.It is well-known that thermal stratification over the cities is in general more unstable than over rural areas due to higher buoyancy of warmer urban air.The modeling results reproduce this pattern: the response of the lapse rate (γ = −dT/dz), calculated between two lowest model levels (10 and 35 m) is positive over urbanized model cells (Figure 7a,b).However, over the nearby rural areas negative responses are observed.In other words, influence of the big city makes the atmosphere more stable over surrounding rural areas.Such model behavior is related to the UHPs.This is clearly illustrated by mean values of lapse rate response over selection of cases with prevailing northern wind (Figure 7c).

Vertical Structure of the UDI/UMI
The spatiotemporal variation of the UDI is partially similar to the UHI following general features of ABL diurnal dynamics.However, it is affected by negative phase shift: the maximum of canopylayer UDI intensity is reached at the evening instead night, which is consistent with classic studies [11, 12,15] and observations for Moscow [36].The negative daily-mean specific humidity response also has a domed shape and extends up to 200-250 m over the city center (Figure 8).
During the day-time, UDI grows up to 1-1.5 km, becomes the most intensive near the surface in the evening and weakens at night.This agrees with existing studies reporting near-zero or even positive (UMI effect) nocturnal urban-rural moisture differences, caused by later and less often reaching the dew point by warmer urban air [15].Nocturnal UMI effect is often detected in Moscow [36], however near-surface urban moisture seems to be underestimated by the model (Figure 3).This could originate from the fact that model does not take into account the vegetation in street canyons and anthropogenic moisture release.Consequently, near-surface UMI could not be detected in the mean seasonal values.However, it is detected above the surface, at the same height as elevated negative temperature response stemming from the cross-over effect.
The only observational study devoted to vertical structure of UDI/UMI over Moscow [36] is based on satellite remote sensing data and reports similar patterns of spatiotemporal variations of summer-time urban-rural moisture differences, but with more pronounced humidity cross-over effect at the evening.Significantly weaker midday humidity cross-over effect, also detected according to remote sensing data, is not reproduced by the model.This is consistent with the fact that the overall underestimation of the modeled urban humidity mentioned above.It should be noted that much more pronounced day-time humidity crossover effect is detected according to WRF simulations for London [33], however even the qualitative comparison in this case is incorrect because of different climate features of Moscow and London (Dfb and Cfb types according to Köppen classification, respectively) and because of the absence of any observational evidence for London.

Vertical Structure of the UDI/UMI
The spatiotemporal variation of the UDI is partially similar to the UHI following general features of ABL diurnal dynamics.However, it is affected by negative phase shift: the maximum of canopy-layer UDI intensity is reached at the evening instead night, which is consistent with classic studies [11, 12,15] and observations for Moscow [36].The negative daily-mean specific humidity response also has a domed shape and extends up to 200-250 m over the city center (Figure 8).
During the day-time, UDI grows up to 1-1.5 km, becomes the most intensive near the surface in the evening and weakens at night.This agrees with existing studies reporting near-zero or even positive (UMI effect) nocturnal urban-rural moisture differences, caused by later and less often reaching the dew point by warmer urban air [15].Nocturnal UMI effect is often detected in Moscow [36], however near-surface urban moisture seems to be underestimated by the model (Figure 3).This could originate from the fact that model does not take into account the vegetation in street canyons and anthropogenic moisture release.Consequently, near-surface UMI could not be detected in the mean seasonal values.However, it is detected above the surface, at the same height as elevated negative temperature response stemming from the cross-over effect.
The only observational study devoted to vertical structure of UDI/UMI over Moscow [36] is based on satellite remote sensing data and reports similar patterns of spatiotemporal variations of summer-time urban-rural moisture differences, but with more pronounced humidity cross-over effect at the evening.Significantly weaker midday humidity cross-over effect, also detected according to remote sensing data, is not reproduced by the model.This is consistent with the fact that the overall underestimation of the modeled urban humidity mentioned above.It should be noted that much more pronounced day-time humidity crossover effect is detected according to WRF simulations for London [33], however even the qualitative comparison in this case is incorrect because of different climate features of Moscow and London (Dfb and Cfb types according to Köppen classification, respectively) and because of the absence of any observational evidence for London.

Urban Effects on Wind Speed and Mesoscale Circulation
Existing studies have noted decrease of the mean near-surface wind speed for urban areas due to higher roughness [11,12,83], which, in general, should extend upward to the ABL.Our modeling results confirm this by negative wind speed response over the city, which is the most pronounced at 150-200 m above the surface (Figure 9).At night, this is more pronounced (within 100-300 m above the ground), but no systematic response could be found at midday.During the evening, the negative response is also observed at 100 m above the ground onwards.Surprisingly, it turns positive between 0 and 100 m during the evening (up to 1.5 m/s on average) and also slightly during the night.
Such effect could be caused by another factor altering the wind regime over urban areas, related to urban-induced mesoscale circulations known as the urban breeze.Such circulations appear

Urban Effects on Wind Speed and Mesoscale Circulation
Existing studies have noted decrease of the mean near-surface wind speed for urban areas due to higher roughness [11,12,83], which, in general, should extend upward to the ABL.Our modeling results confirm this by negative wind speed response over the city, which is the most pronounced at 150-200 m above the surface (Figure 9).At night, this is more pronounced (within 100-300 m above the ground), but no systematic response could be found at midday.During the evening, the negative response is also observed at 100 m above the ground onwards.Surprisingly, it turns positive between 0 and 100 m during the evening (up to 1.5 m/s on average) and also slightly during the night.
Such effect could be caused by another factor altering the wind regime over urban areas, related to urban-induced mesoscale circulations known as the urban breeze.Such circulations appear following the urban effects on stratification and upward motions.However, their existence is proven only by limited number of observational studies (e.g., [12,23,25]), idealistic laboratory and numerical experiments [84,85] or mesoscale simulations for certain case studies [24].Its influence on wind climatology is practically unstudied, and it has only been reported that the wind speed within urban canopy could be higher that over rural areas at calm evenings and nights [11,12,86].
Atmosphere 2017, 8, 50 13 of 25 following the urban effects on stratification and upward motions.However, their existence is proven only by limited number of observational studies (e.g., [12,23,25]), idealistic laboratory and numerical experiments [84,85] or mesoscale simulations for certain case studies [24].Its influence on wind climatology is practically unstudied, and it has only been reported that the wind speed within urban canopy could be higher that over rural areas at calm evenings and nights [11,12,86].To investigate whether the modeled positive wind speed response is related with urban breeze, radial (directed to the city center) wind speed component: To investigate whether the modeled positive wind speed response is related with urban breeze, radial (directed to the city center) wind speed component: where V is absolute wind speed for each model cell, D is the wind direction, and D c is the direction to city center.
In addition to the main selection of cases (days with pronounced UHI during summer 2014), cases with calm weather (mean wind speed within the lowest 150 m over whole D3 domain <3 m/s) were sampled out of the main selection separately for each hour of the day.This results 7-30 cases for each hour.
Temporal and vertical variations of the modeled responses of the wind speed and its radial component, averaged over the Moscow region (within the area extending by 50 and 40 km in latitude and longitude directions correspondingly), show that a positive wind speed response is accompanied by a positive response of its radial component, which is most pronounced in the evening, especially for calm weather conditions (Figure S3).Negative response of radial wind speed component, corresponding to "anti-breeze" part of urban-induced circulation, is detected in the evening higher than 1 km.
Such patterns are confirmed by spatial structure of vector field of the wind response to switching on the UCM.The amplification of wind convergence to the city center at lowest model levels is identified over whole urban area for evening cases with low wind speed (Figure 10a,b).Such pattern is in good agreement with results of idealized studies [84,85].Wind divergence is amplified at the height over 1.5 km (Figure 10c), but spatial structure of the response is more complicated, probably because of interactions with much stronger background flow.
Atmosphere 2017, 8, 50 14 of 25 where is absolute wind speed for each model cell, is the wind direction, and is the direction to city center.
In addition to the main selection of cases (days with pronounced UHI during summer 2014), cases with calm weather (mean wind speed within the lowest 150 m over whole D3 domain <3 m/s) were sampled out of the main selection separately for each hour of the day.This results 7-30 cases for each hour.
Temporal and vertical variations of the modeled responses of the wind speed and its radial component, averaged over the Moscow region (within the area extending by 50 and 40 km in latitude and longitude directions correspondingly), show that a positive wind speed response is accompanied by a positive response of its radial component, which is most pronounced in the evening, especially for calm weather conditions (Figure S3).Negative response of radial wind speed component, corresponding to "anti-breeze" part of urban-induced circulation, is detected in the evening higher than 1 km.
Such patterns are confirmed by spatial structure of vector field of the wind response to switching on the UCM.The amplification of wind convergence to the city center at lowest model levels is identified over whole urban area for evening cases with low wind speed (Figure 10a,b).Such pattern is in good agreement with results of idealized studies [84,85].Wind divergence is amplified at the height over 1.5 km (Figure 10c), but spatial structure of the response is more complicated, probably because of interactions with much stronger background flow.For nocturnal hours with maximum surface-layer UHI intensity, existence of pronounced urban breeze could be also suggested.However, the spatial patterns of the nocturnal wind response are more complicated (Figure 10d,e).They are much stronger affected by the inhibition of incoming air flow, which is characterized by northwestern direction.The magnitude of the wind response vector at southeastern part of the city exceeded the speed of undisturbed flow according to "noURB" simulations.Another interesting feature of nocturnal wind response is the amplification of cyclonic vorticity over the city at the height about 120 m (Figure 10d).Nocturnal "anti-breeze" effect is not pronounced (Figure 10f).Understanding the drivers of such response requires a separate, more detailed analysis.
It should be noted that similar spatial patterns of vector field of the wind response (with quite lower magnitude) are identified in the fields of mean values over much bigger selection of cases without respect to the background wind speed (Figure S4).As such, urban breeze effect is found to be a systematic urban mesoclimatic feature of Moscow.
Quantitative or even qualitative comparison between modeling results and observations in terms of urban breeze features is impossible for Moscow because of the absence of any previous urban breeze studies.Modeling study for Paris megacity [24] reports similar patterns of the late-afternoon urban breeze and its absence during the night with intensive UHI.Urban-induced wind speed amplification reported for Paris (5-7 m/s) is significantly stronger than revealed in this study.However, results for Paris characterize the specially selected case with favorable conditions instead of averaged over several formally selected cases for Moscow, and extreme values of the wind speed response for Moscow are comparable with those reported for Paris (not shown).

Urban Effects on Precipitation and Cloudiness
The modeled summer precipitation response to switching on the UCM, accumulated over a certain summer season, is characterized by a chaotic spatial structure (Figure 11a), which could be explained by the stochastic nature of summer convective showers.However, the noise pattern decreases when increasing the averaging period.Averaging over ten summer seasons allows us to clearly identify a positive precipitation response over the city and its surrounding (Figure 11b).Switching on the UCM increases the mean summer precipitation amount within R1 area (20 km from the city center, see Figure 11b) by 10% in relation to the results of "noURB" simulations, with local maxima up to 25%.Within the distance from 20 to 40 km from the center (R2 area), the mean response value is only ≈4%.However, local maxima up to 25% are also found to the east and northeast of the city, which is its leeward side relative to prevailing western winds.Values of such responses Switching on the UCM increases the mean summer precipitation amount within R1 area (20 km from the city center, see Figure 11b) by 10% in relation to the results of "noURB" simulations, with local maxima up to 25%.Within the distance from 20 to 40 km from the center (R2 area), the mean response value is only ≈4%.However, local maxima up to 25% are also found to the east and northeast of the city, which is its leeward side relative to prevailing western winds.Values of such responses (means over R1 and R2 areas) have a significant inter-annual variability as described in more detail in [48], however they are positive for each among the ten considered summer seasons.
Presented modeling results agrees with a number of studies reporting the urban-induced amplification of convective precipitation due to additional triggering of the moist convection (e.g., [26][27][28][29]).Midday and evening maximum in the daily course of the precipitation response (Figure S5) additionally supports its connection with urban-induced intensification of convective processes.Additionally, such hypothesis is supported by the presence of the positive response of day-time and evening cloudiness, also shifted to the eastern leeward side of the city (Figure 12).This pattern is detected for total, low-level, mid-level and upper-level cloudiness (accordingly, CLCT and CLCL model variables).Switching on the UCM increases the mean summer precipitation amount within R1 area (20 km from the city center, see Figure 11b) by 10% in relation to the results of "noURB" simulations, with local maxima up to 25%.Within the distance from 20 to 40 km from the center (R2 area), the mean response value is only ≈4%.However, local maxima up to 25% are also found to the east and northeast of the city, which is its leeward side relative to prevailing western winds.Values of such responses (means over R1 and R2 areas) have a significant inter-annual variability as described in more detail in [48], however they are positive for each among the ten considered summer seasons.
Presented modeling results agrees with a number of studies reporting the urban-induced amplification of convective precipitation due to additional triggering of the moist convection (e.g., [26][27][28][29]).Midday and evening maximum in the daily course of the precipitation response (Figure S5) additionally supports its connection with urban-induced intensification of convective processes.Additionally, such hypothesis is supported by the presence of the positive response of day-time and evening cloudiness, also shifted to the eastern leeward side of the city (Figure 12).This pattern is detected for total, low-level, mid-level and upper-level cloudiness (accordingly, CLCT and CLCL model variables).Surprisingly, a negative response of the nocturnal low-level cloudiness is simulated by the model over the city (not shown).It is probably caused by the fact that fog and haze (diagnosed as low-level clouds), frequently present over the rural areas during summer nights and occur less frequent due to the UHI that inhibits condensation.However, these processes are not fully understood and need to be analyzed in more detail in future research.This is one reason why dailymean and nocturnal cloudiness is not analyzed in details.
Existing studies for Moscow also report the urban-caused enhancement of the summer precipitation, however its quantities estimates vary significantly.According to the data of the dense network of the weather stations, which existed in Moscow city during the 1950s-1960s, the summer Surprisingly, a negative response of the nocturnal low-level cloudiness is simulated by the model over the city (not shown).It is probably caused by the fact that fog and haze (diagnosed as low-level clouds), frequently present over the rural areas during summer nights and occur less frequent due to the UHI that inhibits condensation.However, these processes are not fully understood and need to be analyzed in more detail in future research.This is one reason why daily-mean and nocturnal cloudiness is not analyzed in details.
Existing studies for Moscow also report the urban-caused enhancement of the summer precipitation, however its quantities estimates vary significantly.According to the data of the dense network of the weather stations, which existed in Moscow city during the 1950s-1960s, the summer precipitation amounts in the city exceed those for surrounding rural areas by 20% [87], while the more recent studies bases on smaller number of weather stations report increase of only 8% [88] or even close to zero [89].However, such variation of the estimates looks reasonable due to very high spatial and temporal variations of the urban-rural precipitation differences, as reported in [88] and confirmed by presented modeling results.
Modeled responses of day-time cloudiness are in good quantitative and qualitative agreement with remote sensing data for Moscow [90] (increase up to 0.4 tenth for daytime total cloudiness over the city in the model and ≈0.6 tenth according to observations).

Discussion
Urban-caused meteorological phenomena, considered within the current study (UHI, UDI and UMI, urban heat plumes and urban breeze, urban enhancement of precipitation and cloudiness) were previously known according to theoretical assumptions, several observational and modeling studies, and idealistic numerical and laboratory experiments.This is the first time that they were analyzed together within similar modeling and data-processing framework based on continuous high-resolution mesoscale simulations, conducted for multiple summer seasons.
In contrast to a wide part of previous studies, presented results shows that the megacity-induced meteorological effects in the lower troposphere persist, not only for certain case studies with favorable conditions, but also for the seasonal means or means over specific selection of the cases, performed in formal ways (e.g., according to wind direction or speed).
Although the detailed analysis of three-dimensional variables (temperature, humidity and wind components) was made only for one summer season (2014) due to technical limitations, preliminary results for other summer seasons shows very similar patterns, which supports the significance of presented estimates and allows to consider them as "quasi-climatologic".However, the urban effect on precipitation and cloudiness cannot be adequately diagnosed from one summer due to the stochastic nature of convective systems.Only the averaging for 10 seasons allows us to detect the urban effects on the precipitation amounts and cloudiness.This demonstrates that the consideration of (quasi-)climatological timescales is required for providing conclusive results about such effects.
In further studies, the developed framework of modeling and data analysis will allow us to make a more robust statistical analysis of the revealed mesoclimatic effects based on the data for longer period, including other seasons.This will allow us to go from "quasi-climatic" estimates of urban-induced meteorological effects to their actual climatology.Some of the presented results (the UHI and its vertical extent, urban breeze circulation, and urban enhancement of precipitation and cloudiness) are supported by earlier observational and modeling studies for Moscow and other big cities, hence confirming their relevance.This opens a wide perspective for more detailed quantitative comparison between modeling results and observational data of existing and new measurement campaigns for Moscow megacity, which seems to be a useful testbed for the verification of urban climate models.
Key scientific value of presented results, in our opinion, is related to the clarification of the various feedbacks and interrelations in the ABL and lower troposphere, interacting with urban climate features.In contrast to a number of studies, focusing on micrometeorological regime canopy layer and considering its features as a function only from land-use and building parameters (e.g., [91][92][93]), we show that the megacity climate features are expressed as mesoscale phenomena and affected by tropospheric feedbacks, including advection, diffusion, urban-induced mesoscale circulations and perturbations of atmospheric convection.These findings are in good agreement with previous modeling studies reporting the presence of non-local mesoscale feedbacks to urban expansion [7,30].The revealed Coriolis-forced veering of the UHI downwind shift by the Ekman layer wind rotation and amplification of cyclonic vorticity seems to be a new finding, not shown in existing studies, which requires more detailed investigation.
Special attention should be paid to our urban modeling results that revealed a systematic urban-induced enhancement of precipitation and daytime cloudiness over the 10 summers.In the number of previous studies about the urban features of atmospheric composition, the higher aerosol concentrations were considered as important factors of urban-induced precipitation increase [28,87,94].However, the current modeling experiments were performed without considering aerosol effects.This shows that urban precipitation enhancement is at least partially caused by pure physical (dynamical) mechanisms induced by surface temperature, humidity and roughness features.Such effect found in our results is consistent with the well-known "negative spatial soil moisture-precipitation coupling" (see, e.g., [95][96][97][98]) indicating that precipitation occurs preferentially over those patches-in our case, Moscow megacity-that are relatively drier than the surroundings.The mechanisms of urban-induced precipitation evidently require more detailed investigation in further studies, however their linkage with urban breeze circulation and urban effect on the lapse rate could be suggested.
From a practical point of view, presented results show the importance of taking the revealed urban-induced meteorological effects and urban climate features into account in various weather, air quality and climate services, not only for the cities, but also for their rural outskirts.Particularly, the revealed enhancement of stable stratification over downwind urban areas could affect the features of pollutants dispersion.Urban breeze circulation is another factor probably enhancing the ventilation of the city, while at the same time promoting the pollutants transfer from industrial areas, located in suburbs, to the city center.Thus, urban-induced processes should be essentially considered in air quality forecasts for urban areas, e.g., in COSMO-ART model [99], which is now used for this purpose for Moscow region [100,101].It is also important to take them into consideration in high-resolution weather forecast for urban areas, biometeorological applications and in downscaling of climate change scenarios as shown in [7,102].
The complex analysis of the presented urban-induced mesoclimatic effects, performed for Moscow megacity, could be a demonstrative illustration for future mega-city climate modeling and data processing frameworks, which shows: (1) the variety of the urban-atmosphere feedbacks affecting the urban climate features; (2) the importance of further research on in-depth understanding of these feedbacks, which would foster better forecasts and more effective urban climate adaption; and (3) the significance of taking them into account in various weather and climate services and applications, including weather and air-quality forecasts, modeling studies in field of biometeorology and the climate change assessments for urban areas.

Conclusions
Using the COSMO-CLM mesoscale model coupled to the urban canopy scheme (UCM) TERRA_URB at the horizontal resolution of 1 km, the urban-induced mesoclimatic effects in the lower atmosphere over Moscow megacity were analyzed for the summer season on continuous time-scales (July-August 2014 for three-dimensional fields of the temperature, humidity and wind, and 10 summer seasons for precipitation and cloudiness).Such effects were investigated as differences between the modeling results for present-day urban surface characteristics calculated according to OpenStreetMaps data and the results of the additional model runs with the UCM switched off.
The model-to-observation comparison shows that temperature and humidity and their respective urban-rural differences are simulated very well, except for a negative bias for the specific humidity in the city.Moreover, the vertical extent of the UHI also agrees qualitatively well with existing studies including the recent results based on the observations of three similar microwave temperature profilers in the city and its suburbs.
Our simulations revealed that the urban effects on the temperature, humidity and wind extends to the lower troposphere at least for 1-2 km and are affected by various atmospheric feedbacks, which defines the complicated spatiotemporal dynamics of such effects.Particularly, the day-time UHI is characterized by a weak intensity (≈0.5 K) and pronounced vertical extent (≈1.5 km).The nocturnal UHI is more intense (>3 K) and confined close to the ground (vertical extent ≈150 m), while cooling above (at the height ≈200-700 m above the ground) is established by the cross-over effect.The urban effects on moisture are characterized by similar vertical extent and negative phase shift in comparison to the UHI.The pronounced daytime and evening UDI extends up to 1-1.5 km, while the nocturnal elevated UMI is formed at same height as the evaluated cooling effect.The spatial structure of the UHI is shifted downwind of the city, and such shift is modulated by a double Coriolis effect resulting from the Eckman-layer forcing and additionally from cyclonic-vorticity forcing.Moreover, pronounced elevated heat plumes are formed downwind of the city during the night.They extend for tens of kilometers and significantly enhanced the stable stratification over rural surroundings.
The urban effects on wind also extend vertically up to 1-2 km above the ground and are characterized by: (1) a lower wind speed due to higher surface roughness; and (2) an amplification of the near-surface wind convergence around the thermal low known as the urban-breeze effect.The latter is most pronounced during the evening.On calm days, such effect results in a significant increase of the wind speed near the surface, which is also identified in terms of seasonal mean values.The urban-induced circulations are not so pronounced during the night and have more complicated spatial patterns, which are not fully understood.Particularly, they result in the amplification of the cyclonic vorticity around the city detected at the height of ≈ 100-150 m over the surface.
By addressing the quasi-climatological scales considering 10 summer seasons, our analysis has revealed a significant influence of the city on the dynamic processes in the lower troposphere by the enhanced precipitation and cloudiness over the city and its leeward side.Herein, the mean summer precipitation amount is increased by 10% over the city with local extremes up to 25%.Similar patterns are identified for the day-time total and low cloudiness.

Supplementary Materials:
The following are available online at http://www.mdpi.com/2073-4433/9/2/50/s1, Figure S1: Spatial distribution of: the urban fraction (a); mean annual value of anthropogenic heat flux (b); building fraction within urban fraction (c); mean building height (d); and street canyon aspect ratio (e) used in numerical simulations.Partially adopted from [43].Black lines represent primary road network in Moscow region according to OpenStreetMaps data.Figure S2: The modeled temperature response to switching on the UCM (∆T) at the: 1st (a,d); 2nd (b,e); and 3rd (c,f) model levels with corresponding heights above the surface equal to 10, 34 and 71 m, averaged over selection of nocturnal cases (0-1 UTC/3-4 MSK) with prevailing: northern (a-c); and southern wind (d-f).Considered cases (same as for Figure 6) are sampled out of selection of days with intensive UHI during June-August 2014.Blue arrows show the wind speed and direction at corresponding model levels according to "URB" simulations.

Data Availability:
The datasets generated during this study are available from the corresponding author on reasonable request.

Figure 1 .
Figure 1.The configuration of the nested domains, used for simulations with the COSMO-CLM model, with surface elevation shown by the color scale, and water surfaces shown by the light-blue color (a).A detailed map of the finest D3 domain, used for the urban climate simulations with the COSMO-CLM model coupled to TERRA_URB, with urban fraction (UF) shown by an additional color scale (b).The location of the weather stations used in the study (see Section 2.6) is also shown (b).The yellow cross indicates Balchug station (considered as the city center), the green cross indicates MSU (Moscow State University) station, and blue crosses indicate the nine rural reference stations around Moscow.The black lines in (b) represent the primary road network.

Figure 1 .
Figure 1.The configuration of the nested domains, used for simulations with the COSMO-CLM model, with surface elevation shown by the color scale, and water surfaces shown by the light-blue color (a).A detailed map of the finest D3 domain, used for the urban climate simulations with the COSMO-CLM model coupled to TERRA_URB, with urban fraction (UF) shown by an additional color scale (b).The location of the weather stations used in the study (see Section 2.6) is also shown (b).The yellow cross indicates Balchug station (considered as the city center), the green cross indicates MSU (Moscow State University) station, and blue crosses indicate the nine rural reference stations around Moscow.The black lines in (b) represent the primary road network.

Figure 2 .
Figure 2. Temporal variation of: the mean rural temperature (a); and urban heat island (UHI) intensity (ΔT) for two urban stations (b,c) according to observations and modeling results during July-August 2014.Mean rural temperature are averaged for each moment over nine rural stations around Moscow; UHI intensity is defined as the difference between the temperature at certain urban station and the mean rural value.

Figure 2 .Figure 2 .Figure 3 .
Figure 2. Temporal variation of: the mean rural temperature (a); and urban heat island (UHI) intensity (∆T) for two urban stations (b,c) according to observations and modeling results during July-August 2014.Mean rural temperature are averaged for each moment over nine rural stations around Moscow; UHI intensity is defined as the difference between the temperature at certain urban station and the mean rural value.

Figure 4 .
Figure 4. Vertical meridional cross-sections through the center of Moscow in the field of the modeled temperature response to switching on the urban canopy model (∆ ): daily-mean (a); mean daytime (13 UTC/16 MSK) (b); evening (16 UTC/19 MSK) (c); and nocturnal (0 UTC/3 MSK) (d) values.The averaging is made for the days with intensive urban heat island during June-August 2014.The contour lines for the values higher than 1 K are indicated at each 1 K. Zero at horizontal axis is the

Figure 4 .
Figure 4. Vertical meridional cross-sections through the center of Moscow in the field of the modeled temperature response to switching on the urban canopy model (∆T): daily-mean (a); mean daytime (13 UTC/16 MSK) (b); evening (16 UTC/19 MSK) (c); and nocturnal (0 UTC/3 MSK) (d) values.The averaging is made for the days with intensive urban heat island during June-August 2014.The contour lines for the values higher than 1 K are indicated at each 1 K. Zero at horizontal axis is the city center (Balchug weather station, see Section 2.2).Bold black line shows the urbanized area (urban fraction higher than 25%).The vertical coordinate is elevation above the sea level (a.s.l.), gray shading represents the relief.

Figure 5 .
Figure 5. Vertical profiles of mean temperature deviation from the top level of measurements (ΔT), calculated according to MTP-5 measurements for the selections of nocturnal cases with pronounced urban heat island during 17 April-3 July 2015.Red shading shows the confidence interval for the urban site (MSU).Adopted from [79].

Figure 5 .
Figure 5. Vertical profiles of mean temperature deviation from the top level of measurements (∆T), calculated according to MTP-5 measurements for the selections of nocturnal cases with pronounced urban heat island during 17 April-3 July 2015.Red shading shows the confidence interval for the urban site (MSU).Adopted from [79].

Figure 6 .
Figure 6.Vertical meridional cross-sections through the center of Moscow in the field of the modeled temperature response to switching on the urban canopy model (∆T), averaged over the selection of nocturnal (0-1 UTC/3-4 MSK) cases with prevailing: southern wind (a); and northern wind (b).Considered cases are sampled out of selection of days with intensive urban heat island during June-August 2014.Designations are similar to Figure 4.

Atmosphere 2017, 8 , 50 11 of 25 Figure 6 .Figure 7 .
Figure 6.Vertical meridional cross-sections through the center of Moscow in the field of the modeled temperature response to switching on the urban canopy model (∆ ), averaged over the selection of nocturnal (0-1 UTC/3-4 MSK) cases with prevailing: southern wind (a); and northern wind (b).Considered cases are sampled out of selection of days with intensive urban heat island during June-August 2014.Designations are similar to Figure 4.

Figure 7 .
Figure 7.The modeled lapse rate (calculated between two lowest model levels) response to switching on the urban canopy model (∆γ) for: daily-mean values (a); mean nocturnal (0-1 UTC/3-4 MSK) values (b); and mean values over nocturnal cases with prevailing northern wind (c).Considered cases are sampled out of selection of days with intensive urban heat island during June-August 2014.Black lines represent primary road network.

Figure 8 .
Figure 8. Vertical meridional cross-sections through the center of Moscow in the field of the modeled specific humidity response to switching on the urban canopy model (∆ ): daily-mean (a); mean daytime (13 UTC/16 MSK) (b); evening (16 UTC/19 MSK) (c); and nocturnal (23 UTC/02 MSK) (d) values averaged for the days with intensive urban heat island during June-August 2014.Designations are similar to Figure 4.

Figure 8 .
Figure 8. Vertical meridional cross-sections through the center of Moscow in the field of the modeled specific humidity response to switching on the urban canopy model (∆q): daily-mean (a); mean daytime (13 UTC/16 MSK) (b); evening (16 UTC/19 MSK) (c); and nocturnal (23 UTC/02 MSK) (d) values averaged for the days with intensive urban heat island during June-August 2014.Designations are similar to Figure 4.

Figure 9 .
Figure 9. Vertical meridional cross-sections through the center of Moscow in the field of the modeled wind speed response to switching on the urban canopy model (∆| |): daily-mean (a); mean daytime (13 UTC/16 MSK) (b); evening (16 UTC/19 MSK) (c); and nocturnal (0 UTC/03 MSK) (d) values.Averaging is made for the days with intensive urban heat island during June-August 2014.Designations are similar to Figure 4.

Figure 9 .
Figure 9. Vertical meridional cross-sections through the center of Moscow in the field of the modeled wind speed response to switching on the urban canopy model (∆|V|): daily-mean (a); mean daytime (13 UTC/16 MSK) (b); evening (16 UTC/19 MSK) (c); and nocturnal (0 UTC/03 MSK) (d) values.Averaging is made for the days with intensive urban heat island during June-August 2014.Designations are similar to Figure 4.

Figure 10 .Figure 10 .
Figure 10.The vector field of the modeled wind response to switching on the urban canopy model (∆ , shown by black arrows) and response of the radial wind speed component (∆ , shown by color) at: 2nd (a,d); 4th (b,f); and 16th (c,e) model levels with corresponding heights above the surface equal to 34 m, 120 m and ≈1.7 km, averaged over selection of: evening (15-16 UTC/18-19 MSK) (a-c); Figure 10.The vector field of the modeled wind response to switching on the urban canopy model (∆V, shown by black arrows) and response of the radial wind speed component (∆V rad , shown by color) at: 2nd (a,d); 4th (b,f); and 16th (c,e) model levels with corresponding heights above the surface equal to 34 m, 120 m and ≈1.7 km, averaged over selection of: evening (15-16 UTC/18-19 MSK) (a-c); and nocturnal (0-1 UTC/3-4 MSK) (d-f) cases with low background wind speed, sampled out from the selection of days with intensive urban heat island during June-August 2014.Green arrows represent wind speed and direction at corresponding model levels according to "noURB" simulations, and the scale is similar to black arrows.

Figure 11 .
Figure 11.Modeled response of summer precipitation amount to switching on the urban canopy model (ΔP): for 2014 year (a); and averaged over 10 summer seasons (2007-2016) (b).Percentage values are calculated in relation to the modeled precipitation amount for "noURB" experiments.Red circles in (b) represent the R1 and R2 areas used in analysis.

Figure 11 .
Figure 11.Modeled response of summer precipitation amount to switching on the urban canopy model (∆P): for 2014 year (a); and averaged over 10 summer seasons (2007-2016) (b).Percentage values are calculated in relation to the modeled precipitation amount for "noURB" experiments.Red circles in (b) represent the R1 and R2 areas used in analysis.

Figure 11 .
Figure 11.Modeled response of summer precipitation amount to switching on the urban canopy model (ΔP): for 2014 year (a); and averaged over 10 summer seasons (2007-2016) (b).Percentage values are calculated in relation to the modeled precipitation amount for "noURB" experiments.Red circles in (b) represent the R1 and R2 areas used in analysis.
Figure S3: The dependence of the modeled response of: the wind speed (∆|V|) (a,c); and its radial component (∆V rad ) (b,d) to switching on the UCM from the height and day time built: for the basic selection of days with pronounced UHI during summer 2014 (a,b); and for the selection of cases with low wind speed (c,d).
Figure S4: The vector field of the modeled wind response to switching on the UCM (∆V, shown by black arrows) and response of the radial wind speed component (∆V rad , shown by color) at 4th model level (20 m), averaged over: all evening (15-16 UTC/18-19 MSK) (a); and nocturnal cases (0-1 UTC/3-4 MSK) (b) for the days with pronounced UHI during the summer of 2014.Designations are similar to Figure 10, but a different color scale is used.
Figure S5: The diurnal course of the: relative (a); and absolute (b) values of the modeled summer precipitation response to switching on the UCM (∆P), averaged over R1 and R2 areas (red and black lines correspondingly) and over the 10 summer seasons.Value for each hour represents the accumulated amount during the previous hour.The green shading represents the diurnal course of hourly precipitation amounts, averaged over whole D3 domain according "noURB" simulations.