Abstract
One important challenge facing the urbanization and global environmental change community is to understand the relation between urban form, energy use and carbon emissions. Missing from the current literature are scientific assessments that evaluate the impacts of different urban spatial units on energy fluxes; yet, this type of analysis is needed by urban planners, who recognize that local scale zoning affects energy consumption and local climate. Satellite-based estimation of urban energy fluxes at neighbourhood scale is still a challenge. Here we show the potential of the current satellite missions to retrieve urban energy budget fluxes, supported by meteorological observations and evaluated by direct flux measurements. We found an agreement within 5% between satellite and in-situ derived net all-wave radiation; and identified that wall facet fraction and urban materials type are the most important parameters for estimating heat storage of the urban canopy. The satellite approaches were found to underestimate measured turbulent heat fluxes, with sensible heat flux being most sensitive to surface temperature variation (−64.1, +69.3 W m−2 for ±2 K perturbation). They also underestimate anthropogenic heat fluxes. However, reasonable spatial patterns are obtained for the latter allowing hot-spots to be identified, therefore supporting both urban planning and urban climate modelling.
Similar content being viewed by others
Introduction
The Urban Energy Balance (UEB) needs to account for the 3D nature of cities, quantifying the fluxes into, out of and the storage change within the control volume (Fig. 1). In the last 15 years, significant advances in understanding urban processes have benefited from enhanced computational capacity, improved resolution of satellite sensors and increased ability to couple advance urban surface parameterization schemes with atmospheric models1,2,3,4,5,6,7,8,9. Earth Observation (EO) has been widely used to study the Urban Heat Island (UHI) phenomenon, but to a lesser extent to quantify heat fluxes10,11,12,13,14,15. Recent studies16,17 have investigated the potential of EO to derive turbulent heat fluxes and identify and analyse the associated uncertainties.
Both Earth system science and urban planning communities need spatially disaggregated UEB data at local scale18,19,20 (neighbourhood, e.g., order (100 m × 100 m) or larger). However, such information is practically impossible to derive for extensive areas by in-situ flux measurements and EO-based estimation of spatio-temporal patterns of different UEB components is challenging. Therefore, the question arises whether EO can provide reliable estimates of UEB at the times of satellite acquisitions. Here we address this by investigating the potential of the current satellite missions to retrieve UEB fluxes at local scale, supported by meteorological observations. Recognising the range of city forms, size and settings, we explored three locations: London (UK), Basel (Switzerland) and Heraklion (Greece).
Results and Discussion
EO data from various sources were used to extract urban surface morphology and characteristics. Surface cover and material type were derived using advanced machine learning and fusion techniques and detailed spectral un-mixing approaches21 (e.g., Fig. 2a). Using EO-derived high resolution Digital Surface Models (DSM), surface roughness parameters (such as plan area index, frontal area index, roughness length and zero-displacement height) were calculated by morphometric analysis22. Examples of the morphometric analysis results are given in Supplementary Fig. S1. Satellite-derived thermal radiance at 1 km × 1 km was downscaled23 to 100 m × 100 m. Although uncertainties of the downscaling method exist, as Mitraka et al.24 noted, the error of downscaled temperatures, evaluated using in-situ surface temperature observations from the micrometeorological flux towers, is within 2.4 K (Supplementary Fig. S2).
The study period (January 2016–May 2017) allowed the seasonal cycle to be examined (Supplementary Fig. S3, London). For all cities this was appropriate with net all-wave radiation (Q*), net change in heat storage (ΔQS) and turbulent sensible heat flux (QH) peaking in summer. QH and ΔQS peak in August, whereas turbulent latent heat flux (QE) is relatively minor throughout the year in city centres.
Deriving Q* maps from EO is not straightforward as urban reflectance and thermal emission present anisotropic behaviour, caused by illumination geometry, 3D urban geometry and the distribution of urban material optical and thermal properties2. Here, the DART (Direct Anisotropic Radiative Transfer) model6 takes into account this anisotropy. It simulates reflected and emitted radiation using information derived from EO satellite images (see more details in Methods). For example, Q* in Basel (Fig. 3a) range between 496 and 633 W m−2. The highest fluxes occur near the river and in a vegetated region in the upper right of the domain. The mean absolute error (MAE, observed - modelled (EO-derived)) in Q* for 2016 was 22.5 W m−2 (Fig. 3d), with agreement within 5%. This corresponds to an average root mean square error (RMSE) of 13 W m−2 and R² > 0.99 (Supplementary Fig. S4). Uncertainty arises from simplifying the complex urban surface structure and the surface elements (e.g., walls) unseen by satellites, except through multiple scattering, although they contribute to Q*.
ΔQS of an urban canopy is approximately 2–6 times larger than for non-urban canopies. Here, ΔQS is estimated using the element surface temperature method (ESTM)3, which reduces the 3D urban structure to one-dimensional (1D) elements (see details in Methods). The sensitivity analysis identified that the fraction of wall facet and materials are the two most important variables in the estimation of ΔQS. Daytime peak (11:00 local time) ΔQS is largest (up to 400 W m−2) in densely built areas where tall buildings dominate (e.g., City of London and Canary Wharf in London Fig. 2b). Areas with low building density have, as expected, smaller ΔQS values. Uncertainty is assessed using two alternative approaches implemented within the SUEWS7 modelling platform: OHM25, and AnOHM26. All approaches produced similar results, but different spatial patterns are identified, mostly depending on the forcing data. To evaluate ESTM performance, observations from less complex sites with uniform land cover are used, where ΔQS can be directly measured with ground heat flux plates. MAE between 1.1 and 21.7 W m−2 is found.
The turbulent heat fluxes are strongly modified by large roughness elements and a complex mix of sources/sinks of heat and water in cities. Here, the Aerodynamic Resistance Method (ARM) is used to derive QH16,17 and QE16 at local scale. To implement this, EO-derived downscaled land surface temperature (LST) and surface roughness parameters are used combined with air temperature and humidity from meteorological stations (see details in Methods). For example, for 23 August 2016 in Basel (Fig. 3b), the highest QH values (up to 270 W m−2) are found in the most densely built-up and industrial areas. QH values have similar spatial patterns and magnitudes to ΔQS, whereas, QE (Fig. 3c) is much smaller in the built-up areas (less than 150 W m−2) and has largest values in vegetated areas. The eddy covariance flux footprint27 is used to compare measured with modelled turbulent heat fluxes. For the Basel BAES (Aeschenplatz) flux tower, the QH MAE for 2016 is 41.2 W m−2 (Fig. 3e), with modelled QH underestimated compared to measured QH. Similarly, the EO-derived QE is underestimated compared to flux tower observations (Fig. 3f). This relates to the uncertainties associated with ARM16,17 and the required input variables. Supplementary Figs S5 and S6 compare the EO-based turbulent sensible and latent heat fluxes with the ground truth at the three sites. A sensitivity analysis17 indicates that, overall, the EO-derived QH is most sensitive to LST variation (+69.3, −64.1 W m−2 for ±2 K perturbation).
Direct measurements of the anthropogenic heat flux QF are extremely difficult28. Here, QF is obtained from the residual of the UEB, and thus includes the cumulative uncertainties of all the terms. Energy balance closure typically compares available energy (Q* - ΔQs) with turbulent heat fluxes (QH + QE). Basel’s two flux towers, have total EC turbulent heat fluxes that tend to be slightly greater than the available energy during winter months (Fig. 4a,b), when QF is expected to be highest due to building space heating. The relative spatial pattern of EO-derived QF for central Basel (winter 2016, Fig. 4c) appears reasonable with clear correlation with building density. However, the absolute QF values for individual pixels include negatives (physically unrealistic). This suggests that individual UEB components are incorrect most likely due to the underestimation of QH. Highest relative QF is found south of the river in the most densely built-up areas (Fig. 4c). Independent QF models based on inventories29 also confirm that buildings are the dominant source of QF. This positive correlation between built-up areas indicates that EO can be used to identify hot-spots with relatively high QF.
Conclusions
Based on the analysis of individual UEB components, we conclude that current satellite missions have the potential to provide information about spatial patterns of urban energy exchanges, if supported by suitable meteorological measurements. Synergistic analysis of specific satellite observations is able to identify and map spatial distributions of heat fluxes at local scale, several times per day (if clear weather). The EO-based methods developed are easily transferable to any city and have the potential to support sustainable planning strategies, since knowledge of UEB patterns at neighbourhood level is needed in urban planning (e.g., to reduce or prevent QH and QF hot spots), health studies (e.g., to estimate impact on thermal comfort) and future proofing (e.g., to plan and implement interventions to reduce heat emissions). Finally, the frequency of heat waves is expected to increase30 with UHI and other urban characteristics exacerbating the respective warming, resulting in increased energy demand for cooling systems in low and mid-latitude cities, which in turn adds to heat emissions and raises temperatures further31. Our satellite-based approach, strengthened by uncertainty minimization in future applications, is expected to advance the current knowledge of the impacts of heat fluxes on energy consumption in cities, leading to the development of tools and strategies to mitigate these impacts, improving thermal comfort and energy efficiency.
Methods
Case studies and datasets
Three cities were selected: highly urbanized mega-city (London, UK); typical central European medium size city, that requires substantial winter heating (Basel, Switzerland); and a smaller, low latitude Mediterranean city that requires substantial summer cooling (Heraklion, Greece). Meteorological data from wireless automatic weather station networks (AWSN) and flux towers (net all-wave radiation and turbulent heat fluxes) in all cities were supplemented by Large Aperture Scintillometers28 in London. High spatial resolution satellite imagery (between 10 m × 10 m and 30 m × 30 m) from Copernicus Sentinel 2 and Landsat 5/7/8, as well as very high resolution (VHR) data (<10 m × 10 m) from Copernicus Contributing Missions, such as SPOT (Satellite Pour l’Observation de la Terre) and WorldView-2, RapidEye, TerraSAR-X and TanDEM-X are used. The EO thermal infrared data are from MODIS (Moderate Resolution Imaging Spectroradiometer), Sentinel 3 (1 km × 1 km) and ASTER (Advanced Spaceborne Thermal Emission and Reflection Radiometer) observations (at 90 m × 90 m)21. Synergies between Sentinels 2 and 3 proved useful to identify fractional surface cover needed in local scale surface emissivity and LST23 estimates. All UEB fluxes are estimated at 100 m × 100 m resolution.
EO data analysis
Surface characteristics include morphology, cover and biophysical parameters. Morphology requires the location and dimensions of urban objects. For this purpose, DSMs are derived from VHR optical stereo imagery (Heraklion) or airborne LiDAR observations (Basel21 and London32). Surface structure parameters (i.e., building volume, sky view factor, plan area index, frontal area index, roughness length and zero-plane displacement height) are derived by morphometric analysis22, using the Urban Multi-scale Environmental Predictor (UMEP)33. Surface cover (plus its spatial and temporal changes) are obtained from EO using advanced machine learning techniques and detailed spectral mixture models34. Multi-temporal acquisitions are used to update the surface cover fraction maps for the classes of interest. By fusing VHR data with Sentinel-2 observations it is possible to categorize different roof materials35 and build spectral libraries36. Other required parameters (e.g., emissivity, albedo and vegetation indices) are derived as outlined in Marconcini et al.21. As daily high spatial resolution thermal imagery is unavailable, the thermal radiance retrieved at 1 km × 1 km resolution (MODIS, Sentinel 3) are downscaled to 100 m × 100 m LST and emissivity using Mitraka et al.’s23 novel technique. Uncertainties in this approach are quantified by Mitraka et al.24. The downscaled LST values are compared to in-situ surface temperature observations from the micrometeorological flux towers.
Estimation of net all-wave radiation
Q* is the difference of the incoming solar (shortwave) radiation (K↓) and atmosphere thermal (longwave) emission (L↓), minus the outgoing reflected solar radiation (K↑) and urban surface thermal emission (L↑). Satellites observe in few spectral bands from restricted viewing directions, whereas Q* is an integral over a wide spectral domain and the whole hemisphere. Here, the DART model9,37 is used to estimate local scale Q*. DART simulates radiative transfer of the urban surface-atmosphere system. It operates on scenes simulated as 3D arrays of rectangular cells. Here a processing chain is developed that calibrates DART with satellite images: the iterative comparison of the model output with images allows the material optical properties (OP) per urban element at the spatial resolution of the satellite data to be obtained. The iterative procedure accounts for the multiple scattering mechanisms and has six steps: 1) Urban morphology and material type are given as input to DART. 2) DART simulates at high spatial resolution (e.g., 2.5 m) all spectral bands of the satellite sensor being used for comparison. The simulation accounts for the specific atmosphere and illumination conditions of each input satellite scene. In the first iteration, spatially constant OP are set per urban element type. 3) DART-simulated spectral image output are georeferenced and spatially resampled to the input satellite scene. 4) DART-simulated and satellite images are compared on a pixel basis to improve the OP map per urban element, using the area of the urban elements within each satellite pixel (derived from the urban morphology and material database). If the DART-simulated and satellite images do not match, the procedure re-iterates from step 2. 5) DART simulates the angular spectral radiative flux along all directions sampling the upper hemisphere. 6) Shortwave K* (K↓ - K↑) and longwave L* (L↓ - L↑) exitance maps are computed as a double integral over the spectrum and the upper hemisphere. Finally, Q* (K* + L*) is resampled to 100 m × 100 m spatial resolution and compared against observations from the micrometeorological towers.
Estimation of the net change in heat storage
ΔQS is the net flow of heat stored in the urban volume (i.e., air, trees, buildings, ground, etc.). In urban areas, the net heat stored in the canopy is a relatively large fraction of Q* and directly evaluating ΔQS in the urban canopy is difficult38. To determine ΔQS we use the Element Surface Temperature Method (ESTM), which reduces the 3D urban volume to four elements, i.e., building roofs, walls, and internal mass and ground (road, vegetation, etc.)6:
where ΔΤi/Δt is the rate of temperature change over the period for each element i, ρci is the volumetric heat capacity, Δxi is the element thickness and fi is the fraction of each surface type. Each element type has sublayers (e.g., a wall can be built up by brick, insulation and wood). Without internal element temperature data the average is determined:
where Q is the conductive heat flux through the surface and k the thermal conductivity. For the internal surfaces (i.e., roof, exterior and interior walls, and floors) the surface temperature of element i is determined by setting the conductive heat transfer out of (in to) the surface equal to the radiative and convective heat losses (gains), as described by Offerle et al.6.
The morphology of the urban surface was derived from high resolution DSM, including 3D information of vegetation32,33. As land cover thermal properties differ, detailed land cover information was required. The seven land cover classes (Supplementary Fig. S1) have three surface materials types for impervious covers and five building categories. These are derived from Urban Atlas39 land use, street view photographs (ground level provided by Google Maps) and local knowledge. The meteorological forcing from the AWSN with the EO downscaled LST are used.
Estimation of the turbulent heat fluxes
The methodology uses the ARM approach40 to estimate QH:
where ρ is the density of air, cp the specific heat of air at constant pressure (1005 J kg−1 K−1), Ta is the air temperature provided by the AWSN and ra is the aerodynamic resistance (s m−1). Analogously, QE is expressed as:
where es* is the saturation water vapour pressure (hPa) at surface air temperature, ea is the atmospheric water vapour pressure (hPa), γ is the psychrometric constant (0.67 hPa K−1) and rs is the stomatal resistance (s m−1). Stomatal resistance is calculated following Kato et al.41 using a simplified equation from Nishida et al.42:
where PAR is the photosynthetic active radiation, rsMIN is the minimum stomatal resistance and rcuticle is the canopy resistance related to the diffusion through the cuticle layer of leaves (105 s m−1). Functions f1 and f2 are from Nishida et al.42 and rsMIN is for each vegetation type41. QE is calculated by land cover type and weighted by the pervious rsMIN type in every pixel. The aerodynamic resistance ra is40:
with u* the friction velocity (m s−1):
where κ is von Kármán’s constant (0.4), zref is height of wind measurements, zd is the zero-plane displacement height, L is the Obukhov length, z0m (z0h) is the roughness length and ψm (ψh) is the stability function for momentum (heat). z0h is often reported within the dimensionless number kβ −1, defined as
EC and scintillometry28 measurements can be used to determine kβ −1 for their flux footprint. Here, flux tower measurements are used as reference values for the magnitude of the fluxes of momentum, sensible and latent heat during the satellite overpass. z0h is calculated from the roughness Reynold’s number (Re*) as43:
where α = 1.29 and Re* (=z0mu*/ν) uses the kinematic molecular viscosity ν (1.461 × 10−5 m s−1).
The displacement height (zd) and z0m are calculated with morphometic parameterization of Kanda et al.44, using urban characteristics (mean and maximum building height, standard deviation of building height, plan area index and frontal area index) derived from 1 m × 1 m DSM, using UMEP22,33.
To evaluate the EO-derived turbulent heat fluxes, EC observations are compared to modeled fluxes in the flux tower source area. Here, the Kormann and Meixner27 footprint model is used.
Estimation of the anthropogenic heat flux
With the other UEB components known, QF is estimated as a residual and hence contains the net errors of all the components. The spatial Energy Balance Closure (EBC) can be determined from regression between (QH + QE) and (Q* - ΔQS) to give an estimate of QF and uncertainty. This requires three assumptions: 1) All sources of anthropogenic heat are released into the environment. 2) Advection is negligible at the scale of interest or cancels between scales: at the microscale (e.g., QH between shadowed and sunlit patches, QE between wet and dry patches); at the local scale (e.g., between parks, water bodies, and built-up areas of different density); and at the meso-scale (e.g., between city and surrounding rural environment: urban breeze; coastal cities sea breezes; topographic induced anabatic/katabatic effects). Moisture advection enhances latent heat flux and can be of similar size to QH, but opposite in sign, thereby essentially offsetting each other45. For this local scale study the advection error is within the QF estimate from energy balance closure. 3) Any unmeasured terms are incorporated in the error of the QF estimate. The resulting QF is evaluated by comparison to QF estimates from alternative approaches, based on inventories29,33,46.
Data availability
This research was conducted under the framework of the URBANFLUXES project that received funding from the European Union’s Horizon 2020 Research and Innovation Programme and joined the H2020 Pilot on Open Research Data. Therefore, all datasets related to the present study are included in the URBANFLUXES Data Repository as Open Data and can be accessed at: http://urbanfluxes.eu/data/.
References
Arnfield, A. J. Two decades of urban climate research: a review of turbulence, exchanges of energy and water, and the urban heat island. Int. J. Climatol. 23, 1–26 (2003).
Voogt, J. A. & Oke, T. R. Thermal remote sensing of urban climates. Remote Sens. Environ. 86, 370–384 (2003).
Offerle, B., Grimmond, C. S. B. & Fortuniak, K. Heat Storage and Anthropogenic Heat Flux in Relation to the Energy Balance of a Central European City Centre. Int. J. Climatol. 25, 1405–1419 (2005).
Weng, Q. Thermal infrared remote sensing for urban climate and environmental studies: Methods, applications, and trends. ISPRS J. Photogramm. Remote Sens. 64, 335–344 (2009).
Grimmond, C. S. B., Roth, M. & Oke, T. R. et al. Climate and More Sustainable Cities: Climate Information for Improved Planning and Management of Cities (Producers/Capabilities Perspective). Procedia Environ. Sci. 1, 247–274 (2010).
Gastellu-Etchegorry, J. P., Yin, T. & Lauret, N. et al. Discrete Anisotropic Radiative Transfer (DART 5) for Modeling Airborne and Satellite Spectroradiometer and LIDAR Acquisitions of Natural and Urban Landscapes. Remote Sens. 7, 1667–1701 (2015).
Ward, H. C., Kotthaus, S., Järvi, L. & Grimmond, C. S. B. Surface Urban Energy and Water Balance (SUEWS): development and evaluation at two UK sites. Urban Climate 18, 1–32 (2016).
Ando, T. & Ueyama, M. Surface energy exchange in a dense urban built-up area based on two-year eddy covariance measurements in Sakai, Japan. Urban Climate 19, 155–169 (2017).
Oke, T. R., Mills, G., Christen, A. & Voogt, J. A. Urban Climates. Cambridge University Press (2017).
Chrysoulakis, N. Estimation of the all-wave Net Radiation Balance in Urban Environment with the combined use of Terra/ASTER multispectral imagery and in-situ spatial data. J. Geophys. Res. 108(D18), 4582 (2003).
Kato, S. & Yamaguchi, Y. Analysis of urban heat-island effect using ASTER and ETM+ Data: Separation of anthropogenic heat discharge and natural heat from sensible heat flux. Remote Sens. Environ. 99, 44–54 (2005).
Kato, S. & Yamaguchi, Y. Estimation of storage heat flux in an urban area using ASTER data. Remote Sens. Environ. 110, 1–17 (2007).
Xu, W., Wooster, M. J. & Grimmond, C. S. B. Modelling of urban sensible heat flux at multiple spatial scales: A demonstration using airborne hyperspectral imagery of Shanghai and a temperature–emissivity separation approach. Remote Sens. Environ. 112, 3493–3510 (2008).
Weng, Q., Hu, X., Quattrochi, D. & Liu, H. Assessing intra-urban surface energy fluxes using remotely sensed aster imagery and routine meteorological data: A case study in Indianapolis, USA. IEEE J. Sel. Top. Appl. Earth Obs. Remote Sens. 7, 4046–4057 (2014).
Chen, S. & Hu, D. Parameterizing Anthropogenic Heat Flux with an Energy-Consumption Inventory and Multi-Source Remote Sensing Data. Remote Sens. 9, 1165 (2017).
Feigenwinter, C., Parlow, E. & Vogt, R. et al. Spatial Distribution of Sensible and Latent Heat Flux in the City of Basel (Switzerland). IEEE J. Sel. Top. Appl. Earth Obs. Remote Sens. https://doi.org/10.1109/JSTARS.2018.2807815 (2018).
Crawford, B. et al. Variability of urban surface temperatures and implications for aerodynamic energy exchange in cities. Q. J. R. Meteorol. Soc. in press (2018).
Seto, K. C., Dhakal, S. & Bigio, A. et al. Human Settlements, Infrastructure and Spatial Planning, Cli- mate Change 2014: Mitigation of Climate Change. Contribution of Working Group III to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change (2014).
Creutzig, F., Baiocchi, G. & Bierkandt, R. et al. Global typology of urban energy use and potentials for an urbanization mitigation wedge. Proc. Natl. Acad. Sci. USA 112, 6283–6288 (2015).
Seto, K. C. & Christensen, P. Remote sensing science to inform urban climate change mitigation strategies. Urban Climate 3, 1–6 (2013).
Marconcini, M., Heldens, W. & Del Frate, F. et al. EO-based Products in support of urban heat fluxes estimation. Proceeding of Joint Urban Remote Sensing Event JURSE 2017, Dubai, UAE, March 30 – April 1, https://doi.org/10.1109/JURSE.2017.7924592 (2017).
Kent, C. W. et al. Evaluation of urban local-scale aerodynamic parameters: Implications for the vertical profile of wind speed and for source areas. Boundary-Layer Meteorol. 164, 183–213 (2017).
Mitraka, Z., Chrysoulakis, N. & Doxani, G. et al. Urban surface temperature time series estimation at the local scale by spatial-spectral unmixing of satellite observations. Remote Sens. 7, 4139–4156 (2015).
Mitraka, Z., Doxani, G., Del Frate, F. & Chrysoulakis, N. Uncertainty Estimation of local-scale land surface temperature products over urban areas using Monte Carlo simulations. IEEE Geosci. Remote Sens. Lett. 13, 917–921 (2016).
Grimmond, C. S. B. & Oke, T. R. Heat storage in urban areas: Local-scale observations and evaluation of a simple model. J. Appl. Meteorol. 38, 922–940 (1999).
Sun, T., Wang, Z. H., Oechel, W. C. & Grimmond, S. The Analytical Objective Hysteresis Model (AnOHMv1.0): methodology to determine bulk storage heat flux coefficients. Geosci. Model Dev. 10, 2875–2890 (2017).
Kormann, R. & Meixner, F. X. An analytical footprint model for non-neutral stratification. Bound.-Layer Meteorol. 99, 207–224 (2001).
Crawford, B., Grimmond, C. S. B. & Ward, H. C. et al. Spatial and temporal patterns of surface–atmosphere energy exchange in a dense urban environment using scintillometry. Q. J. R. Meteorol. Soc. 143, 817–833 (2017).
Gabey, A. M., Grimmond, C. S. B. & Capel-Timms, I. Anthropogenic Heat Flux: advisable spatial resolutions when input data are scarce. Theor. Appl. Climatol. https://doi.org/10.1007/s00704-018-2367-y (2018).
Russo, S., Sillmann, J. & Fischer, E. M. Top ten European heatwaves since 1950 and their occurrence in the coming decades. Environ. Res. Lett. 10, 124003 (2015).
Founda, D. & Santamouris, M. Synergies between Urban Heat Island and Heat Waves in Athens (Greece), during an extremely hot summer (2012). Sci. Rep. 7, 10973 (2017).
Lindberg, F. & Grimmond, C. S. B. Nature of vegetation and building morphology characteristics across a city: influence on shadow patterns and mean radiant temperatures in London. Urban Ecosyst. 14, 617–634 (2011).
Lindberg, F., Grimmond, C. S. B. & Gabey, A. et al. Urban Multi-scale Environmental Predictor (UMEP): an integrated tool for city-based climate services. Environ. Modell. Softw. 99, 70–87 (2018).
Mitraka, Z., Del Frate, F. & Carbone, F. Nonlinear Spectral Unmixing of Landsat Imagery for Urban Surface Cover Mapping. IEEE J. Sel. Top. Appl. Earth Obs. Remote Sens. 9, 3340–3350 (2016b).
Wang, Q., Shi, W., Li, Z. & Atkinson, P. M. Fusion of Sentinel-2 images. Remote Sens. Environ. 187, 241–252 (2016).
Kotthaus, S., Smith, T. E., Wooster, M. J. & Grimmond, C. S. B. Derivation of an urban materials spectral library through emittance and reflectance spectroscopy. ISPRS J. Photogramm. Remore Sens. 94, 194–212 (2014).
Gastellu-Etchegorry, J.-P., Demarez, V., Pinel, V. & Zagolski, F. Modeling radiative transfer in heterogeneous 3-D vegetation canopies. Remote Sens. Environ. 58, 131–156 (1996).
Meyn, S. K. & Oke, T. R. Heat fluxes through roofs and their relevance to estimates of urban heat storage. Energy and Build. 41, 745–752 (2009).
European Environment Agency. Urban Atlas, https://www.eea.europa.eu/data-and-maps/data/urban-atlas (2012).
Voogt, J. A. & Grimmond, C. S. B. Modeling surface sensible heat flux using surface radiative temperatures in a simple urban area. J. Appl. Meteor. 39, 1679–1699 (2000).
Kato, S., Yamaguchi, Y., Liu, C. C. & Sun, C. Y. Surface Heat Balance Analysis of Tainan City on March 6, 2001 Using ASTER and Formosat-2 Data. Sensors 8, 6026–6044 (2008).
Nishida, K., Nemani, R. R., Running, S. W. & Glassy, J. M. An operational remote sensing algorithm of land surface evaporation. J. Geophys. Res. 108, 4720 (2003).
Kanda, M. & Moriizumi, T. Momentum and heat transfer over urban-like surfaces. Boundary- Layer Meteorol. 131, 385–401 (2009).
Kanda, M., Inagaki, A., Miyamoto, T., Gryschka, M. & Raasch, S. A new aerodynamic parametrization for real urban surfaces. Boundary- Layer Meteorol. 148, 357–377 (2013).
Pigeon, G., Lemonsu, A., Grimmond, C. S. B., Durand, P. & Masson, V. Divergence of turbulent fluxes in the surface layer: case of a coastal city. Bound.-Layer Meteorol. 124, 269–290 (2007).
Lindberg, F., Grimmond, C. S. B., Yogeswaran, N., Kotthaus, S. & Allen, L. Impact of city changes and weather on anthropogenic heat flux inEurope 1995-2015. Urban Climate 4, 1–15 (2013).
Acknowledgements
The research was conducted in the framework of the project URBANFLUXES (http://urbanfluxes.eu), that has received funding from the European Union’s Horizon 2020 Research and Innovation Programme, under Grant Agreement No 637519. We would like to thank all URBANFLUXES team members for their contributions to the project and especially Judith Klostermann from ALTERRA; Fabio Del Frate and Danile Latini from GEOK; Andy Gabey, Ting Sun and Christoph Kent from University of Reading; Roland Vogt and Andreas Wicki from University of Basel; Thomas Esch and Wieke Heldens from DLR; Ahmad Al Bitar and Nicolas Lauret from CESBIO; and Giannis Lantzanakis, Nektarios Spyridakis, Giorgos Kogxylakis and Nikos Manioudakis from FORTH. We would also like to thank Michael Abrams from Jet Propulsion Laboratory, California Institute of Technology, for the provision of time series of ASTER images. Other support was received by NERC (TRUC, AirPro).
Author information
Authors and Affiliations
Contributions
Individual contribution of each co-author to the reported research: N.C. led the Heraklion case study, contributed to satellite data analysis, EC observations analysis, interpretation of results and wrote the first version of the paper. S.G. was responsible for the anthropogenic heat flux estimation, contributed to interpretation of results and led the London case study. C.F. was responsible for the turbulent heat fluxes estimation, contributed to interpretation of results and led the Basel case study. F.L. was responsible for the net change in heat storage estimation and contributed to interpretation of results. J.P.G.E. was responsible for the net all-wave radiation estimation and contributed to interpretation of results. M.M. was responsible for satellite data analysis and contributed to interpretation of results. Z.M. developed the surface temperature downscaling method, contributed to the satellite data analysis, to the interpretation of results and to the preparation of the Figures. S.S. contributed to the EC observations analysis, to the satellite data analysis, to the interpretation of results and to the preparation of the Figures. B.C. performed the Scinitillometry and ARM analyses in London and error analyses and QF analyses. F.O. implemented the ESTM approach for net change in heat storage estimation. L.L. performed the DART simulations for net all-wave radiation estimation. W.M. performed the surface temperature observations in London and contributed to the interpretation of results. E.P. supported the turbulent heat fluxes estimation and contributed to the interpretation of results.
Corresponding author
Ethics declarations
Competing Interests
The authors declare no competing interests.
Additional information
Publisher's note: Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Electronic supplementary material
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Chrysoulakis, N., Grimmond, S., Feigenwinter, C. et al. Urban energy exchanges monitoring from space. Sci Rep 8, 11498 (2018). https://doi.org/10.1038/s41598-018-29873-x
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/s41598-018-29873-x
This article is cited by
-
Unprivileged groups are less served by green cooling services in major European urban areas
Nature Cities (2024)
-
Copernicus for urban resilience in Europe
Scientific Reports (2023)
-
Impact of urbanization on surface energy balance components over metropolitan cities of India during 2000–2018 winter seasons
Theoretical and Applied Climatology (2022)
-
Evaluation of the SPARTACUS-Urban Radiation Model for Vertically Resolved Shortwave Radiation in Urban Areas
Boundary-Layer Meteorology (2022)
-
The role of urban trees in reducing land surface temperatures in European cities
Nature Communications (2021)