In this work, some preliminary numerical simulations with the COSMO model including TERRA_URB parameterization have been performed. In particular, this work concerns simulations over a small domain located in southern Italy, in order to test the capabilities of the model in reproducing the main climate features of Urban Heat Islands over this area. Model evaluation has been performed in terms of 2 m temperature in an urban area and in a rural one, in order to highlight the behaviour of the parameterization in different contexts. Numerical results encourage further investigation and development of urban parameterization in very high-resolution configuration of limited area models and specifically of COSMO, to improve the representation of the maximum daily values of temperature and diurnal cycle especially in urban, but also in rural areas. Furthermore a better parameter tuning is still required for this specific test case.
The Urban Heat Island (UHI) effect is one of the most relevant meteorological features of urban environment, since it has substantial implications on human well-being. The increase of built surfaces constitutes the main reason for the formation of UHIs (Milelli, 2016). While natural soil with vegetation uses most of the absorbed radiation in evapotranspiration processes with release of water vapour cooling the surrounding air, paved terrains and buildings tend to absorb a lot of the incident radiation, which is then released as heat. The main characteristic of UHIs is that during the warmest hours of the day, there are small differences between urban and suburban areas; moreover, at sunset the thermal inertia of the city is higher than elsewhere, so there the temperature decreases less than in rural areas, leading to the maximum temperature differences during the night.
In the last years, Limited Area Models have employed very high spatial resolution to improve local forecasts, especially for dangerous convective phenomena and to capture better the urban features. In the same period, a large number of urban land-surface schemes have been developed. They enable the convection permitting atmospheric models to resolve the heterogeneity of cities with applications for heat stress assessment and the development of urban climate adaptation and mitigation strategies.
The main aim of this work is to test the capabilities of the COSMO-LM model
with urban parameterization TERRA_URB in reproducing the main
climate features of UHIs in southern Italy. TERRA_URB is a
bulk parameterization scheme including prescribed Anthropogenic Heat Flux
and Impervious Surface Area datasets (Wouters et al., 2015), offering
intrinsic representation of the urban physics with modifications of input
data, soil module and land atmospheric interactions. Numerical simulations
have been performed over a domain centred over Campania region, employing a
spatial resolution of 0.009
COSMO (Steppeler et al., 2003) is the operational weather forecast model of the Consortium for Small-scale Modeling (COSMO) originally developed by the German Weather Service (DWD). COSMO is based on the governing equations for a compressible flow. The unresolved scale phenomena are taken into account in a statistical manner through parameterizations. In particular, cities are represented by natural land surface with an increased surface roughness length and a reduced vegetation cover. This representation could not reliably capture the urban physics and associated climatic effects including UHIs. The model version used in the present work includes the urban land use scheme TERRA_URB, which parametrizes the effects of buildings, streets and other man-made impervious surfaces. TERRA_URB implements the Semi-empirical Urban canopy parametrization (SURY) (Wouters et al., 2016), which translates canopy parameters (with 3D information) into bulk parameters while preserving the low computational cost. Bulk schemes are able to simulate the general urban climate features at convection permitting scales. They include the interactions between the excess conversion of incoming radiation into sensible heat, heat accumulation, reduced wind speeds and the anthropogenic heating, resulting in the urban heat island effect. However, they do not explicitly resolve the complex processes depending on the local characteristics and the three-dimensional structure of the canopy, which further modulate the urban climate. The model includes two schemes for simulating the bare soil evaporation: a Biosphere-Atmosphere Transfer Scheme (BATS) (Dickinson, 1984) and a resistance-based formulation (RB) (Schulz and Vogel, 2020).
The domain considered, shown in Fig. 1, is centred over Campania and Lazio
regions in southern Italy. The two test cases simulated as preliminary
investigation of the COSMO version including urban parameterization scheme
concerned the following periods: the first one is 8–10 August 2017,
since in these days the city of Naples experienced extreme temperature and
uncomfortable conditions for the population. The second one is 6–8 June 2019, since in these days a sudden increase of the daily maximum
temperature was recorded, from about 27 to about 33
The computational domain with the position of the 2 stations used for the verification (Naples and Grazzanise).
For each day considered, three simulations have been performed:
TERRA_URB switched off (NOURB), TERRA_URB
switched on and BATS (URB_bats), TERRA_URB
switched on and RB (URB_rb.). A spin-up of 6 h is
assumed. Model evaluation has been performed in terms of 2 m temperature
(T_2m) in an urban area (Naples) and in a rural one
(Grazzanise), in order to highlight the behaviour of TERRA_URB in different contexts. The following observational data (identified as
OBS in the figures) have been considered: hourly time series of
T_2m provided by the University of Naples, Fuorigrotta
quarter (urban area 14.19
Figure 2 shows the time series of T_2m over the period
8–10 August 2017 at (a) the urban area of Naples and (c) the rural area
of Grazzanise. Panel (b) shows the spatial T_2m distribution
for 10 August 2017 at 13.00 in a series of points close to Naples
location (Naples is located in point 8). AHF in the figure represents the
Anthropogenic Heat Flux Distribution, which is used as input by
TERRA_URB. In Naples, URB_bats scheme improves
the representation of T_2m with respect to NOURB, especially
the daily minimum values. Anyway, a general underestimation is still evident
since the simulated bare soil evaporation by BATS may be too high, thus
creating a moist and cold bias, particularly during daytime.
URB_rb allows a better behaviour, increasing the average
daily T_2m of 0.2
Time series of T_2m over the period 8–10 August 2017 at
Time series of T_2m over the period 6–8 June 2019 at
The main aim of the present work was to show some preliminary results obtained testing the capabilities of the COSMO model with urban parameterization TERRA_URB in reproducing the main climate features of UHIs in specific areas of southern Italy. The numerical simulations performed on selected test cases have shown that the usage of TERRA_URB produces a remarkable increase of temperature in urban areas, reducing the model bias, especially when the resistance based scheme is adopted for bare soil evaporation. In rural areas TERRA_URB, as expected, does not produce relevant effects so that the model underestimates the diurnal cycle and is affected by a cold bias at day time. Preliminary tests conducted using a “skin temperature formulation” (Schulz and Vogel, 2020) have shown a reduction of the nocturnal temperatures in rural areas and an increase the UHI intensity (urban-rural temperature difference). It is worth mentioning that the proximity of the sea has an impact on the results, since UHIs are less evident close to the sea. Further investigations with a more accurate tuning of this scheme are required and this will be the topic of future work.
The data for this paper are available at the CIRA supercomputing center, Capua (Italy). For data requests please contact Edoardo Bucchignani (e.bucchignani@cira.it).
EB and PM conceived of the presented idea and developed the theory. EB performed the computations. EB wrote the manuscript with support from PM.
The authors declare that they have no conflict of interest.
This article is part of the special issue “19th EMS Annual Meeting: European Conference for Applied Meteorology and Climatology 2019”. It is a result of the EMS Annual Meeting: European Conference for Applied Meteorology and Climatology 2019, Lyngby, Denmark, 9–13 September 2019.
The authors acknowledge the University of Naples “Federico II” for providing the observational datasets. The authors want to thank the colleagues of the COSMO Consortium participating to the PT_AEVUS activity. Their collaboration and support has been fundamental for the present work. For more information, see the dedicated page on the COSMO web page.
This paper was edited by Daniel Reinert and reviewed by two anonymous referees.