An urban neighborhood temperature and energy study from the CAPITOUL experiment with the Solene model

A methodology is proposed to analyze the radiative and thermal exchanges between a small urban neighborhood and the atmosphere based on the use of the thermoradiative model SOLENE and radiometric measurements to optimize the effective values of its constant parameters. Applied to the center of Toulouse, France, the optimization data are building surface temperatures measured with handheld radiometers and a downward-facing pyranometer during one of the intense observation periods of the CAnopy and Particles Interactions in TOulouse Urban Layer (CAPITOUL) experimental campaign. The quality of the simulations is assessed by comparing, without any other model adjustment, the model outputs for two diurnal cycles (1?day in summer and 1?day in winter) against two independent experimental datasets from fixed permanent radiometers and from sensors measuring the solar, infrared and sensible heat fluxes to the atmosphere at the top of a mast 30?m above the roofs. These simulations allow us to further analyze the separate contributions of the different surface classes, roofs, facades and pavement to these fluxes and to compare them with their counterparts observed over a neighborhood of Marseilles city center during a previous experiment. The partition is remarkably similar in Toulouse and Marseilles: 6–7% for the solar radiation, 73–76% for the infrared radiation and 17–20% for the sensible heat flux. The contribution of roofs to the infrared flux appears proportional to their plan area proportion (built density) but not the contributions to the other two fluxes. The contributions of facades to all three fluxes are roughly proportional to their fraction of the total surface area.     

Urban Morphology Influence on Urban Albedo: A Revisit with the Solene Model

This heuristic study of the urban morphology influence on urban albedo is based on some 3,500 simulations with the Solene model. The studied configurations include square blocks in regular and staggered rows, rectangular blocks with different street widths, cross-shaped blocks, infinite street canyons and several actual districts in Marseilles, Toulouse and Nantes, France. The scanned variables are plan density, facade density, building height, layout orientation, latitude, date and time of the day. The sky-view factors of the ground and canopy surfaces are also considered. This study demonstrates the significance of the facade density, in addition to the built plan density, as the explanatory geometrical factor to characterize the urban morphology, rather than building height. On the basis of these albedo calculations the puzzling results of Kondo et al. (Boundary-Layer Meteorol 100:225–242, 2001) for the influence of building height are explained, and the plan density influence is quantitatively assessed. It is shown that the albedo relationship with plan and facade densities obtained with the regular square plot configuration may be considered as a reference for all other configurations, with the exception of the infinite street canyon that shows systematic differences for the lower plan densities. The curves representing this empirical relationship may be used as a sort of abacus for all other geometries while an approximate simple mathematical model is proposed, as well as relationships between the albedo and sky-view factors.     

A Hybrid Wind-Farm Parametrization for Mesoscale and Climate Models

We investigate the area-averaged sensible heat flux (\(Q_{H}\)) obtained with a scintillometer along a 3.1-km path length over the city centre of ?ód?, Central Poland. The annual cycle of \(Q_{H}\) peaks in June but is lower by the middle of summer. In winter, due to a large amount of anthropogenic heat input, \(Q_{H}\) remains positive all day long, with positive night-time fluxes also found during months with frequent cold advection, e.g., June 2010. In the diurnal cycle of this flux, several features specific to urban areas are seen: the peak shifts 1–2 h after noon, the heat flux turns from positive to negative 1–2 h after sunset. In ?ód? \(Q_{H}\) was observed during inflow from the north and north-west, i.e. from the city centre. As this area is mostly covered with impervious materials, most of the heat exchanged between the ground and the overlying air is in the form of sensible heat flux. Under the conditions of inflow from the east and south-east, the maximum heat flux is approximately \(100\,\hbox {W}\,\hbox {m}^{-2}\) lower than during the inflow from the city centre, since more vegetation exists to the east and south-east of the scintillometer path. Cold and warm advection are found to be a vital factor in the observed heat-flux variability in the centre of ?ód?.     

Improving a tree-ring reconstruction from west-central Scandinavia: 900?years of warm-season temperatures

Dendroclimatological sampling of Scots pine (Pinus sylvestris L.) has been made in the province of J?mtland, in the west-central Scandinavian mountains, since the 1970s. The tree-ring width (TRW) chronology spans several thousand years and has been used to reconstruct June?CAugust temperatures back to 1632 bc. A maximum latewood density (MXD) dataset, covering the period ad 1107?C1827 (with gap 1292?C1315) was presented in the 1980s by Fritz Schweingruber. Here we combine these historical MXD data with recently collected MXD data covering ad 1292?C2006 into a single reconstruction of April?CSeptember temperatures for the period ad 1107?C2006. Regional curve standardization (RCS) provides more low-frequency variability than ??non-RCS?? and stronger correlation with local seasonal temperatures (51% variance explained). The MXD chronology shows a stronger relationship with temperatures than the TRW data, but the two chronologies show similar multi-decadal variations back to ad 1500. According to the MXD chronology, the period since ad 1930 and around ad 1150?C1200 were the warmest during the last 900?years. Due to large uncertainties in the early part of the combined MXD chronology, it is not possible to conclude which period was the warmest. More sampling of trees growing near the tree-line is needed to further improve the MXD chronology.     

Circa a.d. 626 volcanic eruption, climatic cooling, and the collapse of the Eastern Turkic Empire

During the late sixth century and early seventh century, the Eastern Turkic Empire (i.e., the Eastern Turkic Khanate) was the most powerful country in the Northeast Asia. It collapsed suddenly in a.d. 630, and historians concluded that the combination of social, political and economic factors, as well as the invasion of the Tang Empire, would be the root cause. Here we suggest that a climatic cooling event ca. a.d. 627–629 could be the direct cause. In a.d. 627–629, the Eastern Turkic Empire experienced severe disasters of snow and frost. Many of the sheep and horses died. People suffered great famine and massive deaths. The Empire fell into severe national crisis and collapsed in a.d. 630. Simultaneously, the Tang Empire also experienced three successive years of frost disasters. Climatic cooling possibly also occurred in other regions. Our investigation of the ca. a.d. 627–629 climatic cooling event also improved our understanding of another problem: was the climatic event due to the impact of a ca. a.d. 626 volcanic eruption?     

Palaeolimnological studies and ancient maps confirm secular climate fluctuations in Amazonia

Secular Amazon discharge oscillations were investigated comparing information from ancient cartography, satellite images, palaeovegetation, sediments, anthropology and climate, focusing on the evolution of Marajó Island in the Amazon estuary. Four phases were identified. (1) 800–1200 a.d.: Amazon discharge increased gradually and eastern Marajó was a vegetation-free tidal plain with energy oscillations. (2) 1200–1350 a.d. was a dry period with lower river discharges, lower energy, increasing water salinities and maximum mangrove extension, coinciding with the extinction of chiefdoms in Marajó. (3) 1350–1540 a.d.: Amazon discharge increased rapidly, inducing a dominance of freshwater vegetation, inundation of east Marajó and likely a regional rise of relative sea-level (RSL). (4) 1540–1750 a.d.: RSL decreased, coinciding with an El Niño-related intense drought ca. 1600 a.d. Emergence of eastern Marajó was rapid and completed in the eighteenth century, under RSL stabilization to current values. The approach used provided evidence of intermittent large changes in Amazon climate, and can facilitate the prediction of future regional dynamics.     

Long-term summer temperature reconstruction inferred from tree-ring records from the Eastern Carpathians

The first 1,000 year long Carpathian tree-ring width chronology was established based on living and subfossil stone pine (Pinus cembra L.) samples from an upper timberline forest located in Calimani Mts. (Romania). Tree-ring data were standardized using the regional curve standardization method in order to preserve the low and medium frequency climate signals. The de-trended index strongly correlates with summer mean temperature both at annual and decadal scales. The Calimani summer mean temperature anomalies were reconstructed for the period ad 1163-2005 applying the rescaling method. This new climate proxy from the Carpathians shows similar fluctuations to other North Hemispheric temperature reconstructions, but with periods of distinct differences. The fingerprint of Little Ice Age in the Calimani area is visible between ad 1370 and 1630 followed by lagged cold decades in ad 1820 and 1840. The recent warming is evident only after the 1980s in our reconstruction.     

High Resolution Simulation of the Variability of Surface Energy Balance Fluxes Across Central London with Urban Zones for Energy Partitioning

The parameterization of surface heat-flux variability in urban areas relies on adequate representation of surface characteristics. Given the horizontal resolutions (e.g. $\approx $ 0.1–1 km) currently used in numerical weather prediction (NWP) models, properties of the urban surface (e.g. vegetated/built surfaces, street-canyon geometries) often have large spatial variability. Here, a new approach based on Urban Zones to characterize Energy partitioning (UZE) is tested within a NWP model (Weather Research and Forecasting model; WRF v3.2.1) for Greater London. The urban land-surface scheme is the Noah/Single-Layer Urban Canopy Model (SLUCM). Detailed surface information (horizontal resolution 1 km) in central London shows that the UZE offers better characterization of surface properties and their variability compared to default WRF-SLUCM input parameters. In situ observations of the surface energy fluxes and near-surface meteorological variables are used to select the radiation and turbulence parameterization schemes and to evaluate the land-surface scheme and choice of surface parameters. For radiative fluxes, improved performance (e.g. $>$ 25 W m $^{-2}$ root-mean-square error reduction for the net radiation) is attained with UZE parameters compared to the WRF v3.2.1 default for all three methods from the simplest to the most detailed. The UZE-based spatial fluxes reproduce a priori expectations of greater energy storage and less evaporation in the dense city centre compared to the residential surroundings. Problems in Noah/SLUCM partitioning of energy between the daytime turbulent fluxes are identified with the overestimation of the turbulent sensible heat and underestimation of the turbulent latent heat fluxes.     

Surface energy balance measurements over a banana plantation in South China

The land surface energy exchange depends highly on the surface properties. Little is known of the energy balance over a typical banana plantation of humid tropics. In this study, we examine the characteristics of surface energy exchange over a typical banana field in South China during the period of May 2010 to April 2011 by using the eddy covariance and micrometeorological tower. The results showed that the diurnal and seasonal variations in surface latent heat flux were larger compared with those over the nearby grassland. The dominant energy partitioning varies with season. The latent heat flux was the main consumer of net radiation in summer, whereas the sensible heat flux was the main consumer in winter. The increasing cloud coverage and rain appear to control the surface energy balance with the development of the monsoon. Due to increased afternoon convective cloud systems in the monsoon active period, downward shortwave radiation was dramatically diminished around 14:00?pm. The annual mean Bowen ratio was 0.69, which fell within the range of other vegetated surfaces. The observed surface energy components were not closed, and the ratio of turbulent fluxes to the available energy was about 77 % in October–January and about 85 % in the other months after considering soil heat and air heat storage.     

Large-Eddy Simulation of Heat Transfer from a Single Cube Mounted on a Very Rough Wall

The local thermal effects in the wake of a single cube with a strong heated rear face, representing a large building in an urban area, are studied using large-eddy simulations (LES) for various degrees of heating, which are characterized by the local Richardson number, $Ri$ . New wall models are implemented for momentum and temperature and comparison of the flow and thermal fields with the wind-tunnel data of Richards et al. (J Wind Eng Ind Aerodyn 94, 621–636, 2006) shows fair agreement. Buoyancy effects are quite evident at low $Ri$ and a significant increase in the turbulence levels is observed for such flows. Apart from the comparisons with experiments, further analysis included the estimation of the thermal boundary-layer thickness and heat transfer coefficient for all $Ri$ . For sufficiently strong heating, the heat transfer coefficient at the leeward face is found to be higher than the roof surface. This suggests that, beyond a certain $Ri$ value, buoyancy forces from the former surface dominate the strong streamwise convection of the latter. Quadrant analysis along the shear layer behind the cube showed that the strength of sweeps that contribute to momentum flux is considerably enhanced by heating. The contribution of different quadrants to the heat flux is found to be very different to that of the momentum flux for lower $Ri$ .     

Modeling regional coral reef responses to global warming and changes in ocean chemistry: Caribbean case study

Climatic change threatens the future of coral reefs in the Caribbean and the important ecosystem services they provide. We used a simulation model [Combo (??COral Mortality and Bleaching Output??)] to estimate future coral cover in the part of the eastern Caribbean impacted by a massive coral bleaching event in 2005. Combo calculates impacts of future climate change on coral reefs by combining impacts from long-term changes in average sea surface temperature (SST) and ocean acidification with impacts from episodic high temperature mortality (bleaching) events. We used mortality and heat dose data from the 2005 bleaching event to select historic temperature datasets, to use as a baseline for running Combo under different future climate scenarios and sets of assumptions. Results suggest a bleak future for coral reefs in the eastern Caribbean. For three different emissions scenarios from the Intergovernmental Panel on Climate Change (IPCC; B1, A1B, and A1FI), coral cover on most Caribbean reefs is projected to drop below 5% by the year 2035, if future mortality rates are equivalent to some of those observed in the 2005 event (50%). For a scenario where corals gain an additional 1?C1.5°C of heat tolerance through a shift in the algae that live in the coral tissue, coral cover above 5% is prolonged until 2065. Additional impacts such as storms or anthropogenic damage could result in declines in coral cover even faster than those projected here. These results suggest the need to identify and preserve the locations that are likely to have a higher resiliency to bleaching to save as many remnant populations of corals as possible in the face of projected wide-spread coral loss.     

Observing Entrainment Processes Using a Small Unmanned Aerial Vehicle: A Feasibility Study

Measurement flights with the meteorological mini aerial vehicle (M $^2$ AV) were performed in spring 2011 to assess the capability of an unmanned aerial vehicle (UAV) to measure the structure of the transition zone between the convective boundary layer and the stably stratified free atmosphere. The campaign took place at the Meteorological Observatory Lindenberg/Richard-Aßmann-Observatory of the German Meteorological Service. Besides the M $^2$ AV flights, observations were made from a 12-m and a 99-m tower, a sodar, two ceilometers, radiosondes, and a tethered balloon with sensor packages at six different levels. M $^2$ AV measurements were intentionally combined with remote sensing systems. The height range of the entrainment zone as well as its diurnal cycle were provided by the remote sensing instruments. The UAV provided the high-resolution in situ data of temperature and wind for the study of turbulent processes. It is shown that the M $^2$ AV is able to maintain constant altitude with very small deviations—a pre-requisite to study processes inside the often quite thin entrainment zone and that M $^2$ AV high-resolution wind and temperature measurements allow for very detailed studies of the fine structure of the atmosphere and thus for the identification of quite local and/or short-duration processes such as overshooting thermals or downward intrusions of warm air. Spatial series measured by the M $^2$ AV during horizontal flights show turbulent exchange of heat in short turbulent bursts at heights close to and within the entrainment zone. Scaled vertical profiles of vertical velocity, potential temperature variance, and sensible heat flux confirm the general shape found by previous measurements and numerical studies.     

Large-Eddy Simulation of Mesoscale Circulations Forced by Inhomogeneous Urban Heat Island

The large-eddy simulation mode of the Weather Research and Forecasting model is employed to simulate the planetary boundary-layer characteristics and mesoscale circulations forced by an ideal urban heat island (UHI). In our simulations, the horizontal heterogeneity of the UHI intensity distribution in urban areas is considered and idealized as a cosine function. Results indicate that the UHI heating rate and the UHI intensity heterogeneity affect directly the spatial distribution of the wind field; a stronger UHI intensity produces a maximum horizontal wind speed closer to the urban centre. The strong advection of warm air from the urban area to the rural area in the upper part of the planetary boundary-layer causes a more stable atmospheric stratification over both the urban and rural areas. The mesoscale sensible heat flux caused by the UHI circulation increases with UHI intensity but vanishes when the background wind speed is sufficiently high $(>$ 3.0  $\mathrm{{m\,s}}^{-1})$ .     

Torneträsk tree-ring width and density ad 500–2004: a test of climatic sensitivity and a new 1500-year reconstruction of north Fennoscandian summers

This paper presents updated tree-ring width (TRW) and maximum density (MXD) from Torneträsk in northern Sweden, now covering the period ad 500–2004. By including data from relatively young trees for the most recent period, a previously noted decline in recent MXD is eliminated. Non-climatological growth trends in the data are removed using Regional Curve Standardization (RCS), thus producing TRW and MXD chronologies with preserved low-frequency variability. The chronologies are calibrated using local and regional instrumental climate records. A bootstrapped response function analysis using regional climate data shows that tree growth is forced by April–August temperatures and that the regression weights for MXD are much stronger than for TRW. The robustness of the reconstruction equation is verified by independent temperature data and shows that 63–64% of the instrumental inter-annual variation is captured by the tree-ring data. This is a significant improvement compared to previously published reconstructions based on tree-ring data from Torneträsk. A divergence phenomenon around ad 1800, expressed as an increase in TRW that is not paralleled by temperature and MXD, is most likely an effect of major changes in the density of the pine population at this northern tree-line site. The bias introduced by this TRW phenomenon is assessed by producing a summer temperature reconstruction based on MXD exclusively. The new data show generally higher temperature estimates than previous reconstructions based on Torneträsk tree-ring data. The late-twentieth century, however, is not exceptionally warm in the new record: On decadal-to-centennial timescales, periods around ad 750, 1000, 1400, and 1750 were equally warm, or warmer. The 200-year long warm period centered on ad 1000 was significantly warmer than the late-twentieth century (p < 0.05) and is supported by other local and regional paleoclimate data. The new tree-ring evidence from Torneträsk suggests that this “Medieval Warm Period” in northern Fennoscandia was much warmer than previously recognized.     

Dynamical changes in the ENSO system in the last 11,000 years

A thorough analysis of a proxy El Nino/Southern Oscillation (ENSO) record indicates that a bifurcation occurred in the ENSO system sometime around 5,000 years b.p. As a result of this bifurcation the attractor became higher dimensional and a new mechanism of instability was introduced. As a consequence of these changes the system switched from a dynamics where the normal condition (La Nina) was dominant to a dynamics characterized by more frequent and stronger El Nino events.     

Turbulent Transport of Momentum and Scalars Above an Urban Canopy

Turbulent transport of momentum and scalars over an urban canopy is investigated using the quadrant analysis technique. High-frequency measurements are available at three levels above the urban canopy (47, 140 and 280 m). The characteristics of coherent ejection–sweep motions (flux contributions and time fractions) at the three levels are analyzed, particularly focusing on the difference between ejections and sweeps, the dissimilarity between momentum and scalars, and the dissimilarity between the different scalars (i.e., temperature, water vapour and $\hbox {CO}_{2})$ . It is found that ejections dominate momentum and scalar transfer at all three levels under unstable conditions, while sweeps are the dominant eddy motions for transporting momentum and scalars in the urban roughness sublayer under neutral and stable conditions. The flux contributions and time fractions of ejections and sweeps can be adequately captured by assuming a Gaussian joint probability density function for flow variables. However, the inequality of flux contributions from ejections and sweeps is more accurately reproduced by the third-order cumulant expansion method (CEM). The incomplete cumulant expansion method (ICEM) also works well except for $\hbox {CO}_{2}$ at 47 m where the skewness of $\hbox {CO}_{2}$ fluctuations is significantly larger than that for vertical velocity. The dissimilarity between momentum and scalar transfers is linked to the dissimilarity in the characteristics of ejection–sweep motions and is further quantified by measures of transport efficiencies. Atmospheric stability is the controlling factor for the transport efficiencies of momentum and heat, and fitted functions from the literature describe their behaviour fairly accurately. However, transport efficiencies of water vapour and $\hbox {CO}_{2}$ are less affected by the atmospheric stability. The dissimilarity among the three scalars examined in this study is linked to the active role of temperature and to the surface heterogeneity effect.     

Climate and cultural history of the Northeastern Yucatan Peninsula, Quintana Roo, Mexico

We inferred the Holocene paleoclimate history of the northeastern Yucatan Peninsula, Mexico, by studying stratigraphic variations in stable isotopes (δ ¹⁸O and δ ¹³C) and lithologic properties (organic matter and carbonate content) in sediment cores taken in 6.3 and 16.2 m of water from Lake Punta Laguna. We present a simple model to explain the lithologic and isotopic variations, and discuss the inferred paleoclimate history in terms of its relation to ancient Maya cultural development. We find evidence for lower lake level and drier climate at about the same time as each major discontinuity in Maya cultural history: Preclassic Abandonment (150–250 A.D.), Maya Hiatus (534 to 593 A.D.), Terminal Classic Collapse (750–1050 A.D.), and Postclassic Abandonment (mid-fifteenth century). Although these broad temporal correlations suggest climate played a role in Maya cultural evolution, chronological uncertainties preclude a detailed analysis of climate changes and archaeologically documented cultural transformations.     

Temperature variability in China in an ensemble simulation for the last 1,200?years

Regional temperature anomalies in China during 800?C2005 ad in an ensemble simulation with the atmosphere?Cocean general circulation model ECHAM5/MPIOM subject to anthropogenic and natural forcings are compared to reconstructions. In a mutual assessment of three reconstructed data sets and two ensemble simulations with different solar forcings, a reconstructed data set and a simulated ensemble for weak solar variability are selected for further comparison. Temperature variability in the selected simulated and reconstructed data shows a continuous power spectrum with weak long-term memory. The simulation reveals weak long-term anomaly periods known as the Medieval Warm Period (MWP), the Little Ice Age (LIA), and the Modern Warming (MW) in the three considered regions: Northeast, Southeast, and West China. The ensemble spread yields an uncertainty of ±0.5°C in all regions. The simulated temperature varies nearly synchronously in all three regions, whereas reconstructed data hint to increased decadal variability in the West and centennial variability in the Northeast. Cold periods are found in 1200?C1300 and in 1600?C1900 ad in all regions. The coldest anomalies which are caused by volcanic eruptions in the beginnings of the thirteenth and the nineteenth centuries are only partly consistent with reconstructed data. After 1800, the annual cycle reduces in the Northeast and on the Tibetan plateau, whereas the eastern Pacific shows an enhanced summer?Cwinter contrast.     

Thermal Stratification Effects on Flow Over a Generic Urban Canopy

The influence of local surface heating and cooling on flow over urban-like roughness is investigated using large-eddy simulations. By adjusting the incoming or outgoing heat flux from the ground surface, various degrees of local thermal stratification, represented by a Richardson number \((Ri_\tau )\) , were attained. Drag and heat transfer coefficients, turbulence structure, integral length scales, and the strength of quadrant events that contribute to momentum and heat fluxes were obtained and are compared with locally stable, neutral and unstable flows. With increasing \(Ri_\tau \) , or equivalently as the flow characteristics change from local thermal instability to stability, a gradual decline in the drag and heat transfer coefficients is observed. These values are found to be fairly independent of the type of thermal boundary condition (constant heat flux or constant temperature) and domain size. The maps of anisotropy invariants showed that for the values of \(Ri_\tau \) considered, turbulence structures are almost the same in shape for neutral and unstable cases but differ slightly from those in the stable case. The degree of anisotropy is found to decrease as \(Ri_\tau \) increases from \(-2\) to 2.5. Compared to the neutral case, the integral length scales are shortened in the streamwise and vertical direction by ground cooling, but enhanced in the vertical direction with ground heating. Quadrant analysis showed that an increase in floor heating increases the strength of ejections above the canopy. However, the contributions of updrafts or downdrafts to the heat flux are found not to be significantly influenced by the type of local thermal stratification for the values of \(Ri_\tau \) considered. From the octant analysis, the transport mechanisms of momentum and heat above the canopy are found to be very similar in both locally unstable and stable flows.     

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.