Next Article in Journal
An Optical Fibre Depth (Pressure) Sensor for Remote Operated Vehicles in Underwater Applications
Next Article in Special Issue
Improving the Accuracy of Urban Environmental Quality Assessment Using Geographically-Weighted Regression Techniques
Previous Article in Journal
A Biocompatible Colorimetric Triphenylamine- Dicyanovinyl Conjugated Fluorescent Probe for Selective and Sensitive Detection of Cyanide Ion in Aqueous Media and Living Cells
Previous Article in Special Issue
Urban Growth Modeling Using Cellular Automata with Multi-Temporal Remote Sensing Images Calibrated by the Artificial Bee Colony Optimization Algorithm
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Remote Sensing of Urban Microclimate Change in L’Aquila City (Italy) after Post-Earthquake Depopulation in an Open Source GIS Environment

1
Department of Civil, Constructional and Environmental Engineeering (DICEA), Sapienza University of Rome, I-00184 Rome, Italy
2
Fondazione Mach, Centre for Technology Transfer, I-38010 S. Michele all’Adige (TN), Italy
3
Department of Civil, Construction-Architectural and Environmental Engineering (DICEAA), University of L’Aquila, I-67100 L’Aquila, Italy
*
Author to whom correspondence should be addressed.
Sensors 2017, 17(2), 404; https://doi.org/10.3390/s17020404
Submission received: 11 October 2016 / Revised: 10 February 2017 / Accepted: 13 February 2017 / Published: 19 February 2017

Abstract

:
This work reports a first attempt to use Landsat satellite imagery to identify possible urban microclimate changes in a city center after a seismic event that affected L’Aquila City (Abruzzo Region, Italy), on 6 April 2009. After the main seismic event, the collapse of part of the buildings, and the damaging of most of them, with the consequence of an almost total depopulation of the historic city center, may have caused alterations to the microclimate. This work develops an inexpensive work flow—using Landsat Enhanced Thematic Mapper Plus (ETM+) scenes—to construct the evolution of urban land use after the catastrophic main seismic event that hit L’Aquila. We hypothesized, that, possibly, before the event, the temperature was higher in the city center due to the presence of inhabitants (and thus home heating); while the opposite case occurred in the surrounding areas, where new settlements of inhabitants grew over a period of a few months. We decided not to look to independent meteorological data in order to avoid being biased in their investigations; thus, only the smallest dataset of Landsat ETM+ scenes were considered as input data in order to describe the thermal evolution of the land surface after the earthquake. We managed to use the Landsat archive images to provide thermal change indications, useful for understanding the urban changes induced by catastrophic events, setting up an easy to implement, robust, reproducible, and fast procedure.

1. Introduction

The historic center of L’Aquila is still largely uninhabited after a major seismic event on 6 April 2009 (Figure 1); thus, it may have undergone a change in its microclimate due to the sudden absence of inhabitants, which may have caused a lower impact in heating during winter. It is likely that an equal but opposite change may have occurred in areas of new settlements, originating after the seismic event.
Broad literature describes procedures used, using commercial software, to evaluate the urban heat island effect and its changes [1,2,3]; however, nowadays, free and open source software offers the possibility of developing a complete workflow to describe such dynamics [4]. The choice to use freely-available remote-sensed images, and the free and open source tools for Geographic Information System (GIS) analysis of data could aid local governance in observing the spatial dynamics of urban areas, and to determine successful policies to manage such events in the long-term, and at low costs. On the other hand, there is not an extensive literature on the use of remote sensing in order to monitor depopulation: In Reference [5], depopulation is estimated using nighttime light images; in Reference [6], instead, land-cover maps using Landsat images were used; and, finally, in Reference [7], an exposure model at a global scale, suitable for different geo-hazards, using teledetected images, was implemented.
Satellite imagery was only used in the post-earthquake, L’Aquila City scenario, to detect damaged buildings for a quick damage assessment [8,9,10] using the geometric and thematic capabilities of Very High Resolution (VHR) images. To the best of the authors’ knowledge, no use of satellite thermal images on this area, in order to monitor depopulation, is known. Before the main shock, the inhabitants of L’Aquila numbered around 16,000, about 10,000 inhabitants, plus the university student population, which was approximately 6000 residents [11]. Some decommissioned industrial areas were present, neighboring the city, along the valley that crosses the city from NW to SE. The inhabitants of the center were relocated to areas neighboring the old city by two projects: C.A.S.E. (Project for sustainable and earthquake safe houses: 13,000 inhabitants) and M.A.P. (project for temporary housing modules: 2800 inhabitants). After the event, development of the city following the morphology of the NW to SE valley took place; thus, new residential and business buildings were relocated along the NW to SE direction, the center of the city almost empty. This trend can be observed, for example, in Figure 2a,b [11], where the public offices of national agencies (orange in Figure 2a,b), local agencies, such as prefecture and regional administration, courts (blue in Figure 2a,b), and university buildings (green in Figure 2a,b) are represented, before and after the event; the historic center is visible in the figures. A rigorous delineation of the urbanized area, after the event, is still not possible.
It can be easily noticed that, before the seismic event, the economic and social center of the city was concentrated in a very narrow area of the historic city center (Figure 2a).
After the earthquake, most of the public offices were relocated to a surrounding area where unused industrial buildings were available (Figure 2b), leaving the center almost deserted (see Figure 1).
In this study, we set the hypothesis that the sudden relocation of houses and offices changed the human-induced micro-climate, and we used remotely-sensed images in order to track the spatio-temporal evolution of the heat island effect.

2. Materials and Methods

The image set consists of 11 Landsat ETM+ scenes, all in the same path and row (Table 1). Each image covers the area of the central part of Italy: at the center of each image lies the urban area of Aquila (42°20′55″ N 13°23′54″ E) and the neighboring villages. The elevation of the entire area ranges between 559 and 2508 m above the sea level, with different combinations of urban and agricultural landscapes, different types of forest coverage, and bare soils (Figure 3).
The Digital Terrain Model (DTM) was obtained from the Shuttle Radar Topography Mission (SRTM) dataset, freely available online, and has a spatial resolution of 90 m [12].
Table 1 reports the attributes of the image set used in the analyses (Figure 4). The time of acquisition (Time) is expressed in Coordinated Universal Time. Latitude and Longitude refer to the center of each image and are expressed in geographic coordinates (decimal degrees) in the WGS84 datum.
It can be observed that all the images were acquired, almost at the same time, around 8:37 a.m. local time (GMT + 1), when domestic and office heating were surely being used.
Starting from the Landsat data, a procedure to estimate the Land Surface Temperature (LST) is presented.
Among the many free and open source tools for data analysis, we decided to choose the GRASS GIS software (GRASS stands for Geographic Resources Analysis Support System) and the R statistical software. These two tools can be easily coupled together to run geographic and statistical analyses on data in one seamless environment [13].
In order to filter the pixels where the Scan Line Corrector failed (SLC-off), the corresponding SLC-off mask was applied to each of the 9 bands of the Landsat ETM+ image.
Radiometric correction is applied to all bands of every Landsat scene, converting the Digital Numbers (DN) to at-sensor (or top-of-atmosphere) spectral radiance (L), using Equation (1). The obtained values must be corrected for solar variability, caused by annual changes in the Earth–Sun (Esun) distance (d), producing unitless, top-of-atmosphere reflectance (ρAS, Equation (2)) using data on solar intensity. Esun is the band-specific exo-atmospheric solar constant and is the solar elevation angle.
The thermal bands are first converted from DN to ‘at-sensor radiance’, and then to effective at-sensor temperature in Kelvin. The values involved in these transformations are obtained from the metadata available for each scene [14].
L = ( L m a x L m i n D N m a x D N m i n ) ( D N D N m i n ) + L m i n ,
ρ A S = π d 2 L E s u n cos ( 90 θ s )
The i.landsat.toar GRASS module [15] has been used to perform the radiometric correction of the images, and was used for the calculation of the top-of-atmosphere radiance and top-of-atmosphere reflectance. In order to calculate the cloud-cover contamination, obtained from both the reflective and thermal proprieties of the scene, we used the i.landsat.acca GRASS module [16,17]. The cloud-cover signature was then used to mask each Landsat scene.
For the application developed in this work, the values of ground reflectance are needed, thus, atmospheric correction is required. We chose an absolute atmospheric method to correct the image, relying on a mechanistic understanding of atmospheric effects, deducing the values for the atmospheric parameters from the information contained within the scene itself, rather than using externally-measured data and individually adjusting each image. In brief, the conversion of the at-sensor radiance to atmospherically-corrected surface reflectance is obtained using Equation (3), where the values of Tz, Tv, Edown, and Lhaze must be determined.
ρ = π d 2 ( L L h a z e ) T v ( E s u n cos ( 90 θ s ) T z + E d o w n )
Among the methods available, we chose the Dark Object Subtraction model (DOS), which performs a re-calibration of the image, assuming that, for pixels with low reflectance values, any apparent reflectance should be due to atmospheric scattering effects [18,19]. The DOS (Equation (4)) model searches for the lowest DN value found in a fixed number of pixels, and this selected DN, the Starting Haze Value (SHV), is converted into radiance. It is unlikely that images contain pixels that are completely black, thus, correction considers a 1% reflectance (L1%) of these areas [14].
L 1 % = 0.01 E s u n cos ( 90 θ s ) π d 2 L h a z e = S H V r a d L 1 % ,
The modified Dark Object Subtraction model was applied in order to incorporate the effects of atmospheric aerosols in the atmosphere [20]. This method works iteratively, following Equation (5).
τ = cos ( θ z ) l n ( 1 L m i n 0.01 T v π ( E d o w n + cos ( θ z ) d 2 T z E s u n ) cos ( θ z ) d 2 E s u n ) T v = e τ T z = e τ cos ( θ z ) ,
In order to obtain the LST maps, the Land Surface Emissivity (LSE) maps were calculated using a semi-empirical method (SEM), based on the Normalized Difference Vegetation Index (NDVI, Equation (6)): Among the different approaches found in literature, such as those described in References [21,22,23], the one based on NDVI Thresholds and Vegetation Proportion (PV, Equation (7)) has been chosen [24]. A detailed description of the method can be found in References [2,25]. This method is preferred among all others because it automates the land-cover detection, and should adjust for seasonal variability in the forested landscape.
N D V I = N I R V I S N I R + V I S ,
P v = ( N D V I N D V I m i n N D V I m a x N D V I m i n ) ,
Eventually, the LST maps were estimated using the procedure found in References [2,26], and using the radiative transfer equation (RTE, Equation (8)), where B(Ts) is the black-body radiance given by Plank’s law, Ts is the LST estimated for the pixel, and ε is the land surface emissivity.
L = ( ε B ( T s ) + ( 1 ε ) L ) τ + L ,
The down-welling (L↓) and up-welling (L↑) atmospheric radiance, and the transmittivity between the surface and the sensor, needed for the LST computation, were retrieved using Atmospheric Correction Parameter Calculator software, available online, which is based on a model (MODTRAN, MODerate resolution atmospheric TRANSmission) that simulates the emission and adsorption of electromagnetic radiation from 0.2 to 100 μm spectral range (from middle ultraviolet to far infrared) [27,28].
After processing, each LST map was cropped to the study area. The complete flowchart of the procedure is shown in Figure 5.
As we are not interested in the estimation of surface temperature, but, rather, on the description of urban exploitation, the LST of each scene was re-scaled in the range of [–1, 1] in order to enhance the visual contrast of the spatial gradient of the LST among the different land uses, and to reduce the inter-annual and intra-annual variability of LST due to meteorological conditions.
This post-processing was performed using R software, extended with the raster, ggplot2, and ggmap packages [29,30,31,32].

3. Results

As the combined effects of cloud cover and SLC-off greatly degraded the informative power of seven out of 11 LST maps, only four scenes (Figure 6 and Figure 7) were used to draw a qualitative analysis of the dynamics of the urban developing of Aquila and its neighbors, by using the urban micro-climate evolution as an index of the evolution (Figure 8).
A noticeable amount of qualitative information can be drawn from the interpretation of the LST-rescaled maps (LSTr, Figure 9). Although the algorithm based on NDVI values is one of the simplest methods available for obtaining Land Surface Emissivity, it proved to be reliable in identifying the different landscapes of the study area, and the overall dynamics in urban landscape evolution after the earthquake.
The first LSTr map, which refers to the first Landsat scene available before the earthquake, is used as reference for the other three scenes, which recorded the events that followed. In this first LSTr map, urban structures (towns and roads) are effectively matched by pixels with higher temperature values. One month after the earthquake, and as of then, these structures changed remarkably: The Aquila town center is still visible, but the values of LSTr are more similar to the values of the surrounding landscapes.
Particularly interesting are the high values of LSTr recorded in the pixels near the area of Coppito, the industrial district of Bazzano, and the industrial district Pile (Figure 9). For these three areas, the values clearly remained high, and are distinct from those of surrounding pixels, before and after the disaster. For the Coppito area, this evidence could be explained by the forced displacement of survivors into existing military structures in the aftermath of the earthquake.
The explanation for the constant high values of LSTr in the two industrial areas is linked to the more resistant structures, built using anti-seismic criteria. These earthquake-resistant building were used by the local university and other public agencies after the disaster (Figure 2b).
On the other hand, the destructive aftermath of the earthquake can be quickly identified by matching the first LSTr map with the other three, focusing on the towns of Paganica and Tempera (northeast of L’Aquila (Figure 9)): The urban structures of the two villages disappear from the LSTr time series. Indeed, it is well-known that Paganica and Tempera were severely hit by the disaster and remain uninhabited.

4. Discussion and Further Developments

This work aims at the development of a parsimonious work-flow for processing Landsat Enhanced Thematic Mapper Plus scenes (ETM+) in order to track changes to the urban micro-climate of L’Aquila City and neighboring towns after an earthquake event on 6 April 2009.
We know that the historic center of L’Aquila, which is still largely unpopulated, was suddenly abandoned after the event, as were some neighboring towns. However, some new towns have been built in order to help people recover, thus reshaping the urban structure of the area.
The idea underlying this work is that the forced relocation of inhabitants may have caused a lower impact of heating in winter, and of air conditioning in summer, and that these changes could potentially be reflected in the thermal signature, remotely sensed, in the images. In addition, an opposite change should have occurred in areas of new settlements that originated after the event.
We have found that the Land Surface Temperature (LST) estimation, when rescaled in order to enhance spatial gradients, effectively tracks changes to the urban micro-climate, and can be used to describe the evolution of urban land use after the catastrophic earthquake event was taken into consideration, even when the quality of the image is sub-optimal due to failure of the Scan Line corrector (SCL-off).
Although an assessment of the quantitative effects of the aftermaths was not the main focus of this work (the procedure here uses a very parsimonious model, in terms of input data, in order to precisely estimate land surface temperature), the procedure presented here can be used to perform a qualitative, but correct, description of the dynamics of people’s relocation.
With this work, we do not intend to propose a novel method for LST estimation, but, rather, we set up a pragmatic, robust, and easy-to-use work-flow, which can be replicated by non-technical users, ensuring the conceptual correctness and reproducibility of image processing. One important point to underline here is that the ability to apply this work-flow to Landsat scenes is limited to the availability of the Atmospheric Correcting Parameters of thermal bands (average atmospheric transmission, effective bandpass upwelling radiance, and effective bandpass down-welling radiance), which can be easily obtained using free web services, or can be directly implemented into the workflow using the Application Programming Interface (API) available with the MODTRAN6 model, with enhanced physics features, including a line-by-line (LBL) algorithm [33,34,35,36].
Clearly, the overall precision and performance of the work-flow could be improved by using alternative and more sophisticated algorithms for the estimation of LST; however, this step can only be performed with a procedure of validation of the output maps with near-sensed measures, when available. These data were not available to us at the time of writing, and this step is envisaged for inclusion in future work.
The authors are aware that it is impossible to attribute all surface temperature variability of an affected or unaffected landscape to earthquake effects. Certainly, there are effects of differential solar zenith angles, and other ambient weather variables, due to the non-identical temporal scale of images on potential surface temperature variability. For instance, the results presented in Figure 9 show over a month difference between March 2009, and May 2010, which may have accounted for some of the spatial variations in the surface temperature of the study area. However, as noted in the Abstract, the authors managed to use ¥ Landsat archive images to provide relative thermal change indications in the study area, which is useful to understand the relative urban changes (more than in finding absolute values of temperatures). Additionally, as stated, all images were acquired at almost the same time (around 8:37 a.m., local time, GMT + 1) when heating was surely being used in houses and offices, for all the four images, and sunshine had just started to warm the ground. If this work captures the interest of both local and national governance, the procedure presented here could be encapsulated in a general-purpose GIS software package as a Quantum-GIS (QGIS) plug-in, an R package, or a GRASS module, in order to provide an additional tool, such as a black-box, for decision-makers.
It is worth mentioning that this black-box, if the GNU (General Public License) is granted with the software, could be improved in a “community driven way”; the more this software is used, the more weak points of the work-flow can be tackled, and new approaches to LST estimation can be integrated into the software, as well as alternative algorithms for the atmospheric correction of Landsat scenes and for cloud cover assessment.
Our proposed work-flow could be used retrospectively, as in this paper, or online, after the occurrence of similar disastrous events. In this case a specialization of the algorithm for the Landsat 8 product should also be considered. The usage of a pool of free and open source software and freely-available data sources greatly lowers the total cost for the adoption of the procedure presented in this paper, and could help in “search and rescue” (SAR) campaigns, especially in third-world countries or in remote regions that are difficult to reach.
Finally, the use of demographic and economic data, when available, could provide an additional source of information for a more precise description of urban micro-climate evolution after a disastrous event.

Acknowledgments

This work was partially financed by a grant funded by Italian Ministry for Research and University (MIUR) in the framework of a “Progetto di rilevante interesse nazionale” (PRIN) project N. 20107XJHCC_004, “Tecniche geomatiche innovative ed emergenti di rilievo, telerilevamento (da aereo, satellite, uav) e webgis per la mappatura del rischio in tempo reale e la prevenzione del danno ambientale”, coordinated by Professor Raffaele Santamaria of Parthenope University, Naples, Italy.

Author Contributions

Study conception and design, Acquisition of Data, Analysis and Interpreation of Data, Drafting of the Article and Critical Revision are based on Valerio Baiocchi, Donatella Dominici and Fabio Zottele‘s collaborative effort.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Voogt, J.A.; Oke, T.R. Thermal remote sensing of urban climates. Remote Sens. Environ. 2003, 86, 370–384. [Google Scholar] [CrossRef]
  2. Sobrino, J.A.; Jiménez-Muñoz, J.C.; Paolini, L. Land surface temperature retrieval from LANDSAT TM 5. Remote Sens. Environ. 2004, 90, 434–440. [Google Scholar] [CrossRef]
  3. Imhoff, M.L.; Zhang, P.; Wolfe, R.E.; Bounoua, L. Remote sensing of the urban heat island effect across biomes in the continental USA. Remote Sens. Environ. 2010, 114, 504–513. [Google Scholar] [CrossRef]
  4. Costantino, D.; Angelini, M.G. Thermal monitoring using an ASTER image. J. Appl. Remote Sens. 2016, 10, 046031. [Google Scholar] [CrossRef]
  5. Uchida, S.; Takeuchi, W.; Hatoyama, K.; Mazurov, Y. Socio-economic impact of Trans-Siberian railway after the collapse of Soviet Union by integrated spatial data analysis. In Proceedings of the 8th IGRSM International Conference and Exhibition on Geospatial & Remote Sensing (IGRSM 2016), Kuala Lumpur, Malaysia, 13–14 April 2016.
  6. Bagan, H.; Yamagata, Y. Landsat analysis of urban growth: How Tokyo became the world’s largest megacity during the last 40 years. Remote Sens. Environ. 2012, 127, 210–222. [Google Scholar] [CrossRef]
  7. Pittore, M.; Wieland, M.; Fleming, K. Perspectives on global dynamic exposure modelling for geo-risk assessment. Nat. Hazards 2016, 1–24. [Google Scholar] [CrossRef]
  8. Baiocchi, V.; Brigante, R.; Dominici, D.; Milone, M.V.; Mormile, M.; Radicioni, F. Automatic three-dimensional features extraction: The case study of L’Aquila for collapse identification after April 06, 2009 earthquake. Eur. J. Remote Sens. 2014, 47, 413–435. [Google Scholar] [CrossRef]
  9. Baiocchi, V.; Dominici, D.; Giannone, F.; Zucconi, M. Rapid building damage assessment using EROS B data: the case study of L’Aquila earthquake. Eur. J. Remote Sens. 2012, 44, 153–165. [Google Scholar] [CrossRef]
  10. Contreras, D.; Blaschke, T.; Tiede, D.; Jilge, M. Monitoring recovery after earthquakes through the integration of remote sensing, GIS, and ground observations: The case of L’Aquila (Italy). Cartogr. Geogr. Inf. Sci. 2016, 43, 115–133. [Google Scholar] [CrossRef]
  11. Linee di Indirizzo Strategico per la Ricostruzione ed il Rilancio Della Citta. Available online: http://www.comune.laquila.gov.it/moduli/output_immagine.php?id=3030 (accessed on 28 January 2017).
  12. Hole-Filled SRTM for the Globe Version 4. Available online: http://srtm.csi.cgiar.org (accessed on 28 January 2017).
  13. Geographic Resources Analysis Support System (GRASS) Software. Available online: http://grass.osgeo.org (accessed on 28 January 2017).
  14. Goslee, S.C. Analyzing Remote Sensing Data in R: The landsat package. J. Stat. Softw. 2011, 43. [Google Scholar] [CrossRef]
  15. Irish, R.R. Landsat 7 Automatic Cloud Cover Assessment. In Algorithms for Multispectral, Hyperspectral, and Ultraspectral Imagery VI, Proceedings of the SPIE, Orlando, FL, USA, 24 April 2000; Shen, S.S., Descour, M.R., Eds.; SPIE: Bellingham, WA, USA, 2004; Volume 4049, pp. 348–355. [Google Scholar]
  16. Irish, R.R. Landsat 7. Science Data Users Handbook; Goddard Space Flight Center: Greenbelt, MD, USA, 1998. Available online: http://landsathandbook.gsfc.nasa.gov/orbit_coverage/ (accessed on 28 January 2017).
  17. Irish, R.R.; Barker, J.L.; Goward, S.N.; Arvidson, T. Characterization of the Landsat-7 ETM+ Automated Cloud-Cover Assessment (ACCA) Algorithm. Photogramm. Eng. Remote Sens. 2006, 72, 1179–1188. [Google Scholar] [CrossRef]
  18. Chavez, P.S., Jr. Radiometric calibration of Landsat Thematic Mapper multispectral images. Photogramm. Eng. Remote Sens. 1989, 55, 1285–1294. [Google Scholar]
  19. Chavez, P.S., Jr. Image-based atmospheric corrections revisited and improved. Photogramm. Eng. Remote Sens. 1996, 62, 1025–1036. [Google Scholar]
  20. Song, C.; Woodcock, C.E.; Seto, K.C.; Lenney, M.P.; Macombe, S.A. Classification and Change Detection Using Landsat TM Data, When and How to Correct Atmospheric Effects? Remote Sens. Environ. 2001, 75, 230–243. [Google Scholar] [CrossRef]
  21. Sobrino, J.A.; Raissouni, N. Towards remote sensing methods for land cover dynamic monitoring: Application to Morocco. Int. J. Remote Sens. 2000, 21, 353–366. [Google Scholar] [CrossRef]
  22. Valor, E.; Caselles, V. Mapping Land Surface Emissivity from NDVI: Application to European, African, and South American Areas. Remote Sens. Environ. 1996, 57, 167–184. [Google Scholar] [CrossRef]
  23. Van de Griend, A.A.; Owe, M. On the Relationship between Thermal Emissivity and the Normalized Difference Vegetation Index for Natural Surfaces. Int. J. Remote Sens. 1993, 14, 1119–1131. [Google Scholar] [CrossRef]
  24. Carlson, T.N.; Ripley, D.A. On the Relation between NDVI, Fractional Vegetation Cover, and Leaf Area Index. Remote Sens. Environ. 1997, 62, 241–252. [Google Scholar] [CrossRef]
  25. Alipour, T.; Sarajian, M.R.; Esmaeily, A. Land Surface Temperature Estimation from Thermal Band of Landsat Sensor, Case Study: Alashtar City. The International Archieves of the Photogrammetry, Remote Sensing and Spatial Information Sciences. XXXVIII-4/C7. Available online: https://www.researchgate.net/publication/215444122_Land_Surface_Temprature_Estimation_Base_on_Thermal_Band_of_ETMSensor (accessed on 28 January 2017).
  26. Coll, C.J.; Galve, J.M.; Sanchez, M.; Caselles, V. Validation of Landsat-7/ETM+Thermal-Band Calibration and Atmospheric Correction with Ground-Based Measurements. IEEE Trans. Geosci. Remote Sens. 2014, 48, 547–555. [Google Scholar] [CrossRef]
  27. Barsi, J.A.; Barker, J.L.; Schott, J.R. An Atmospheric Correction Parameter Calculator for a Single Thermal Band Earth-Sensing Instrument. In Proceedings of the IGARSS03, Toulouse, France, 21–25 July 2003.
  28. Barsi, J.A.; Schott, J.R.; Palluconi, F.D.; Hook, S.J. Validation of a Web-Based Atmospheric Correction Tool for Single Thermal Band Instruments. In Earth Observing Systems, Proceedings of X SPIE, San Diego, CA, USA; SPIE: Bellingham, WA, USA, 2005; Volume 5882. [Google Scholar]
  29. Hijmans, R.J. Raster: Raster: Geographic Data Analysis and Modeling. R Package Version 2.3-12. Available online: http://CRAN.R-project.org/package=raster (accessed on 28 January 2017).
  30. Wickham, H. Ggplot2: Elegant Graphics for Data Analysis; Springer: New York, NY, USA, 2009. [Google Scholar]
  31. Kahle, D.; Wickham, H. Ggmap: A Package for Spatial Visualization with Google Maps and OpenStreetMap. R Package Version 2.3. Available online: http://CRAN.R-project.org/package=ggmap (accessed on 28 January 2017).
  32. R Core Team. R: A Language and Environment for Statistical Computing; R Foundation for Statistical Computing: Vienna, Austria, 2011; Available online: http://www.R-project.org/ (accessed on 28 January 2017).
  33. Berk, A.; Conforti, P.; Kennett, R.; Perkins, T.; Hawes, F.; Van den Bosh, J. MODTRAN6: A major upgrade of the MODTRAN radiative transfer code. Proc. SPIE 2014, 9088. [Google Scholar] [CrossRef]
  34. Mlawer, E.J.; Payne, V.H.; Moncet, J.L.; Delamere, J.S.; Alvarado, M.J.; Tobin, D.C. Development and recent evaluation of the MT_CKD model of continuum absorption. Philos. Trans. R. Soc. A 2012, 370, 1–37. [Google Scholar] [CrossRef] [PubMed]
  35. Clough, S.A.; Shephard, M.W.; Mlawer, E.J.; Delamere, J.S.; Iacono, M.J.; Cady-Pereira, K.; Boukabara, S.; Brown, P.D. Atmospheric radiative transfer modeling: A summary of the AER codes. J. Quant. Spectrosc. Radiat. Transf. 2005, 91, 233–244. [Google Scholar] [CrossRef]
  36. Mishchenko, M.I.; Travis, L.D.; Lacis, A.A. Multiple Scattering of Light by Particles: Radiative Transfer and Coherent Backscattering; Cambridge University Press: Cambridge, UK, 2006. [Google Scholar]
Figure 1. The city center after the event.
Figure 1. The city center after the event.
Sensors 17 00404 g001
Figure 2. L’Aquila, before (a) and after (b) the main seismic event: Overlapped grid is Universal Transverse Mercator Projection (UTM) WGS84/ETRF00. Scale bar (upper right) is 4 km, and the north-pointing arrow (lower left) indicates the geographic north.
Figure 2. L’Aquila, before (a) and after (b) the main seismic event: Overlapped grid is Universal Transverse Mercator Projection (UTM) WGS84/ETRF00. Scale bar (upper right) is 4 km, and the north-pointing arrow (lower left) indicates the geographic north.
Sensors 17 00404 g002
Figure 3. Study area of L’Aquila and its neighboring areas. Scenes copyrighted by Google Imagery, Mapdata, and Terrametrics. The image was reprojected from Spherical Mercator to WGS84; image is around 20 km in width.
Figure 3. Study area of L’Aquila and its neighboring areas. Scenes copyrighted by Google Imagery, Mapdata, and Terrametrics. The image was reprojected from Spherical Mercator to WGS84; image is around 20 km in width.
Sensors 17 00404 g003
Figure 4. The complete image set of the Landsat ETM+ scenes where the three bands (red, green blue) acquired in the visible spectrum were blended together (Red, Blue and Green (RGB) composite), projected to UTM33N WGS84 datum, images are around 23 km in width.
Figure 4. The complete image set of the Landsat ETM+ scenes where the three bands (red, green blue) acquired in the visible spectrum were blended together (Red, Blue and Green (RGB) composite), projected to UTM33N WGS84 datum, images are around 23 km in width.
Sensors 17 00404 g004
Figure 5. The workflow of the procedure adopted for image processing.
Figure 5. The workflow of the procedure adopted for image processing.
Sensors 17 00404 g005
Figure 6. RGB composite of the four image-sets used for Land Surface Temperature estimation Scan Line corrector (SCL)-off mask applied, no cloud cover assessment); images are around 12 km in width.
Figure 6. RGB composite of the four image-sets used for Land Surface Temperature estimation Scan Line corrector (SCL)-off mask applied, no cloud cover assessment); images are around 12 km in width.
Sensors 17 00404 g006
Figure 7. RGB composite of the four image-sets used for Land Surface Temperature estimation, with both SCL-off and cloud cover assessment masks. Images are around 12 km in width.
Figure 7. RGB composite of the four image-sets used for Land Surface Temperature estimation, with both SCL-off and cloud cover assessment masks. Images are around 12 km in width.
Sensors 17 00404 g007
Figure 8. Land Surface Temperature Maps (LST) without rescale. The gray pixels correspond to the SCL-off and cloud cover assessment masks. Images are around 12 km in width.
Figure 8. Land Surface Temperature Maps (LST) without rescale. The gray pixels correspond to the SCL-off and cloud cover assessment masks. Images are around 12 km in width.
Sensors 17 00404 g008
Figure 9. Land Surface Temperature Maps with rescale applied (LSTr). The gray pixels correspond to the SCL-off and cloud cover assessment masks. Images are around 12 km in width.
Figure 9. Land Surface Temperature Maps with rescale applied (LSTr). The gray pixels correspond to the SCL-off and cloud cover assessment masks. Images are around 12 km in width.
Sensors 17 00404 g009
Table 1. Descriptive attributes of the Landsat Enhanced Thematic Mapper Plus (ETM+) images.
Table 1. Descriptive attributes of the Landsat Enhanced Thematic Mapper Plus (ETM+) images.
IdDateGreenwich Medium Time (GMT)Lat (Deg.)Lon (Deg.)LST Successfully Processed
LE71900312008084ASN0024/03/200809:37:2141.7666513.61165
LE71900312008132ASN0011/05/200809:37:1341.7650413.59003
LE71900312008308ASN0003/11/200809:36:1141.7628813.63519
LE71900312008324ASN0019/11/200809:36:2541.7602413.64066
LE71900312009086ASN0027/03/200909:37:2641.7639313.6099*
LE71900312009134ASN0014/05/200909:37:3841.7623413.59009*
LE71900312009278ASN0005/10/200909:37:1841.7651713.60086*
LE71900312009326ASN0022/11/200909:37:5441.7631213.65504
LE71900312010089ASN0030/03/201009:39:2441.7544813.58893
LE71900312010121ASN0001/05/201009:39:2341.7554213.62054*
LE71900312010249ASN0006/09/201009:39:2941.7560113.62678
* Successfully processed.

Share and Cite

MDPI and ACS Style

Baiocchi, V.; Zottele, F.; Dominici, D. Remote Sensing of Urban Microclimate Change in L’Aquila City (Italy) after Post-Earthquake Depopulation in an Open Source GIS Environment. Sensors 2017, 17, 404. https://doi.org/10.3390/s17020404

AMA Style

Baiocchi V, Zottele F, Dominici D. Remote Sensing of Urban Microclimate Change in L’Aquila City (Italy) after Post-Earthquake Depopulation in an Open Source GIS Environment. Sensors. 2017; 17(2):404. https://doi.org/10.3390/s17020404

Chicago/Turabian Style

Baiocchi, Valerio, Fabio Zottele, and Donatella Dominici. 2017. "Remote Sensing of Urban Microclimate Change in L’Aquila City (Italy) after Post-Earthquake Depopulation in an Open Source GIS Environment" Sensors 17, no. 2: 404. https://doi.org/10.3390/s17020404

APA Style

Baiocchi, V., Zottele, F., & Dominici, D. (2017). Remote Sensing of Urban Microclimate Change in L’Aquila City (Italy) after Post-Earthquake Depopulation in an Open Source GIS Environment. Sensors, 17(2), 404. https://doi.org/10.3390/s17020404

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop