Changes in observed climate extremes in global urban areas

Climate extremes have profound implications for urban infrastructure and human society, but studies of observed changes in climate extremes over the global urban areas are few, even though more than half of the global population now resides in urban areas. Here, using observed station data for 217 urban areas across the globe, we show that these urban areas have experienced significant increases (p-value <0.05) in the number of heat waves during the period 1973–2012, while the frequency of cold waves has declined. Almost half of the urban areas experienced significant increases in the number of extreme hot days, while almost 2/3 showed significant increases in the frequency of extreme hot nights. Extreme windy days declined substantially during the last four decades with statistically significant declines in about 60% in the urban areas. Significant increases (p-value <0.05) in the frequency of daily precipitation extremes and in annual maximum precipitation occurred at smaller fractions (17 and 10% respectively) of the total urban areas, with about half as many urban areas showing statistically significant downtrends as uptrends. Changes in temperature and wind extremes, estimated as the result of a 40 year linear trend, differed for urban and non-urban pairs, while changes in indices of extreme precipitation showed no clear differentiation for urban and selected non-urban stations.


Introduction
For the first time, more than half of the global population resides in urban areas (Grimm et al 2008), a number that is expected to increase to 60% by 2030, and 70% by 2050 (World Health Organization 2014). Between 1900 and 2000, the number of urban areas with more than one million population increased from 17 to 388 (Millennium Ecosystem Assessment 2003). Urban areas are centers of wealth, human population, and built infrastructure and are considered by some to be 'first responders' to climate change (Rosenzweig et al 2010). Moreover, cities are the fundamental units for climate change mitigation and adaptation (Georgescu et al 2014).
Climate variability and change can exert profound stresses on urban environments, which are sensitive to heat waves, droughts, and changes in the frequency and magnitude of flash floods (Rosenzweig et al 2011).
Temperature-related climate extremes have been shown to be strongly related to human health (Patz et al 2005, Hayhoe et al 2010, Hondula et al 2014. Extreme precipitation events, which are projected to increase in a warming climate (Allen and Ingram 2002, O'Gorman and Schneider 2009, Min et al 2011a, cause disproportionate damage to urban transportation systems, and pose challenges to urban stormwater drainage systems (Schreider et al 2000, Rosenberg et al 2010, Mishra et al 2012a. In the recent past, flood disasters have affected many large urban areas including Bangkok (2011), Brisbane (2011), Guangdong (2007, Mumbai (2005), and Dresden (2002) (Liao 2012). In addition to temperature and precipitation extremes, extreme wind storms can cause enormous damage to urban areas (Munich Re 1999, Swiss Re 2001. Despite the potentially large and costly damages from climate extremes in urban areas, most Any further distribution of this work must maintain attribution to the author (s) and the title of the work, journal citation and DOI.
observation-based studies have focused on large spatial extents (i.e. national, regional, or continental) (Easterling et al 2000, Zhang et al 2005, Alexander et al 2006, Min et al 2011, Donat et al 2013. For instance, Donat et al (2013) analyzed trends in climatic extremes using a gridded dataset, but did not evaluate changes in extreme events in urban areas specifically. Studies that address climate extremes in global urban areas and disparities in urban and non-urban areas have been far fewer (Ashley et al 2005, Hallegatte et al 2007. This may be in part because of complications in obtaining high quality observations from station data or developing gridded datasets using stations within the vicinity of urban areas (Mishra and Lettenmaier 2011). Many studies have shown the influence of urbanization on meteorological forcing (Zhang et al 2009, Kishtawal et al 2010, Shepherd et al 2010a, 2010b, Chen et al 2011, which is not the focus of this work. Here we analyze observed changes in climate extremes in global urban areas during the last four decades  using daily data from selected Global Summary of the Day (GSOD) stations. We hypothesize that globally, observed climate changes in urban areas, notwithstanding major land use/ cover change over the last 40 years, are dominantly due to large scale changes, rather than local land cover. We report changes in daily temperature and precipitation extremes as well as changes in the frequency of extreme windy days in 217 urban areas across the globe. Moreover, we discuss disparities in changes in climate extremes in 142 paired urban and non-urban stations.

Analysis approach
We obtained daily observations for precipitation, air temperature (mean, maximum, and minimum), and mean wind speed from the GSOD data produced by the National Climatic Data Center (NCDC). GSOD data are quality controlled by NCDC through automated quality checks. Most random errors are removed and further corrections are applied to control for changes in instrumentation, station moves, and changes in time of observation (ftp://ftp.ncdc.noaa. gov/pub/data/gsod/readme.txt, accessed on 1 December 2013). The GSOD data are updated on a near realtime basis with a lag of 1-2 days. Historical data are available from 1929; however, we selected the period 1973-2012 for analysis as this period has the largest number of relatively complete and consistent records. Further details on the GSOD data can be obtained from the website http://www.climate.gov/globalsummary-day-gsod.
We identified all urban areas globally (about 650) with population greater than 250 000 using the ESRI Arc-GIS shapefile (ESRI 2008, http://www.baruch. cuny.edu/geoportal/data/esri/esri_intl.htm, assessed on 20 November 2013) and selected GSOD station(s) located within 25 km of these urban areas. GSOD records for these stations were checked for record length and missing values. Our selection of urban areas based on population was further verified with urban extent derived using the MODIS data (Schneider et al 2009). Schneider et al (2009) identified urban areas as built environments including non-vegetative and human constructed elements. They were identified using supervised classification of MODIS 1 km land cover data (see Schneider et al (2009) for details). This comparison shows that urban areas selected based on population fall within the urban extents derived using the satellite based dataset (figures S9-S14). We removed all stations that had one or more years of missing data or more than 10% missing values in any year during the period of 1973-2012. After identifying candidate stations based on the proximity criteria, and applying the quality checks, we were left with 217 stations (table S2) with relatively complete records for the period 1973-2012, most of which were located at airports close to urban areas (see supplemental table S1, available at stacks.iop.org/ERL/10/ 024005/mmedia). Furthermore, about 90% of the stations are within 15 km of the center of the urban area they represent.
We analyzed daily extremes for temperature, precipitation (24 h totals), and wind. For temperature extremes, we analyzed heat and cold waves as well as extreme hot days and nights. Heat waves were defined as periods during which the daily maximum temperature stayed above (below) the empirical 99th-percentile (estimated for the period 1973-2012) consecutively for six or more days. Similarly cold waves were identified using daily minimum temperatures. The frequency of extreme hot days was estimated using the 99th-percentile of daily maximum temperatures, also for the period 1973-2012. The frequency of extreme hot nights was identified using the 99th percentile of daily minimum temperatures. For precipitation extremes, annual maximum precipitation, frequency of precipitation extremes, and fraction of total precipitation occurring due to extreme events were estimated for each year during the period of 1973-2012. Annual maximum precipitation was defined as the maximum daily precipitation in each year. To estimate the frequency of extreme precipitation events, we used the 95th-percentile of daily precipitation for rain days (precipitation >1 mm d −1 ) during the 1973-2012 reference period. We selected the 95th percentile threshold for extreme precipitation events because the 99th percentile threshold leads to some years without extreme precipitation events, which makes detection of trends difficult. We carefully evaluated the influence of the selection of specific thresholds on trends of extreme precipitation events and found that the nature of the trends remains the same with a higher or lower threshold. The number of days with precipitation greater than or equal to the 95th-percentile was used as an index for extreme precipitation in each year. The fractional contribution from the extreme precipitation events to the total precipitation for each year was computed using the same 95th-percentile threshold. Extreme windy days were defined as days with a daily mean wind speed that exceeded the 99th-percentile threshold of daily mean wind speed.
We estimated changes in each of the statistics using the non-parametric Mann-Kendall (Mann 1945) trend test and Sen's slope method (Sen 1968). The Mann-Kendall method has been widely used for trend detection in hydrologic and climate data (Mishra and Lettenmaier 2011, Yue and Wang 2002). We also computed field significance (Livezey and Chen 1983) using the method described in Yue and Wang (2002). The statistical significance of changes in the mean and distribution was tested using the two sided Ranksum and Kolmogorov and Smirnov (KS) (Massey 1951) tests, respectively at 5% significance level. We divided the 217 locations into six regions based on the number of stations in each region (figure 1(a)): Africa (AF), East Asia (EA), Europe (EU), India (IN), North America (NA), and South America (SA) and applied the field significance test on a regional basis. Changes in climatic extremes were estimated by multiplying the linear trend obtained from the Mann-Kendall test with the length of the analysis period (i.e. 40 years).
To understand differences in precipitation and temperature extremes, we identified paired urban and non-urban stations based on population density as described in Mishra and Lettenmaier (2011). We obtained the gridded population density (1 km spatial resolution) from the Global Urban-Rural Mapping Project (GURMP v3) for the year 2000 from the CIE-SIN's web site (http://sedac.ciesin.columbia.edu/gpw/ global.jsp, accessed on 8 December 2013). Stations with population density less (more) than 200 (500) person/km 2 were identified as non-urban (urban) stations. We successfully verified urban and non-urban pairs based on urban extents derived from MODIS (figure S15). Zhou et al (2013) reported the influence of spatial clustering of stations on urban-heat island effect; however, our aim was to find out disparities in trends between stations located in urban and paired non-urban areas. We identified 142 urban and nonurban station pairs with median population densities 2540 and 99 persons/km 2 , respectively. For each urban station (out of 142), a non-urban station (based on the quality and completeness of the dataset) was selected using distance thresholds of 50, 75, 100, 125, 150, and 200 km from the urban station location. We excluded non-urban pairs that exceeded the 200 km threshold to avoid large changes in climate and topography. As much as possible, we selected pairs that are not affected by geographical disparities associated with elevation, location, and climate. However, to avoid these differences entirely, the station density in urban and non-urban areas needs to be greatly improved. Figure 1 shows changes in heat and cold waves globally during the reference period . Because heat and cold waves do not occur every year, we evaluated changes using the pooled data, where mean heat waves for each year were estimated for all the urban areas in a given region. Pooled time series of global urban areas showed statistically significant increases in the number of heat waves per urban area during the last four decades (change = 0.32 heat waves per urban area during the reference period, p < 0.05 on normalized data; figure 1(b)). During the last 40 years, the five years with the largest number of heat waves (aggregated over all regions) were 1998, 2009, 2010, 2011, and 2012. Taken over all 217 stations, the average number of heat waves by decade exhibited a consistent increase with the highest number of heat waves occurring during the most recent decade (2003-2012; figure 1(c)). Figure 1(d) shows that the frequency of cold waves generally declined. For instance, the pooled cold waves showed a statistically significant (p-value <0.05) decline of 0.16 cold waves per urban area during the reference period. The five years with the largest number of cold waves were 1973, 1974, 1976, 1981, and 1983. The average number of cold waves by decade declined until 1993-2002, followed by a slight increase in the number of cold waves during the most recent decade (2003)(2004)(2005)(2006)(2007)(2008)(2009)(2010)(2011)(2012). Median changes over all sites showed increases in heat waves and declines in cold waves (figures 1(f) and (g)). Increases in heat waves were field significant for AF, EA, EU, and NA. Increases in IN and SA were not field significant ( figure 1(f)). Field significant decreases in cold waves were found in NA (figure 1(g)), however, decreases in cold waves in EU were not field significant. No statistically significant trends in the number of cold waves were detected in the remaining regions.

Changes in temperature extremes
Changes in the frequency of extreme hot days (above 99th-percentile) are shown in figure 2. During the reference period, the number of extreme hot days increased significantly at many sites ( figure 2(a)). However, a few urban areas located in EA showed significant declining trends. About 48% of the sites showed statistically significant increases ( figure 2(b)). On the other hand, fewer than 2% of the total urban areas experienced significant declines in the frequency of extreme hot days. Median increases of eight days in the frequency of extreme hot days were found for urban areas with statistically significant uptrends ( figure 2(c)). Moreover, all regions showed field significance for increases in the frequency of extreme hot days (figure 2(d)) with median change (estimated using average linear trend multiplied by 40 years) ranging from three to nine days between 1973 and 2012. Similar to increases in the number of extreme hot days, the number of extreme hot nights (estimated from daily minimum temperature data) increased significantly in a majority of urban areas across the globe (supplemental figure S1(a)). About 63% of the urban areas experienced statistically significant increase in the frequency of hot nights ( figure S1(b)). Between 1973 and 2012, a median increase of ten hot nights was found in urban areas that showed significant increases (figure S1(c)). Increases in the number of hot nights were field significant in all the regions except for IN (figure S1(d)). Along with temperature extremes, a significant warming in mean temperature was found at a majority of urban areas across all seasons ( figure S6). Moreover, a majority of the selected urban areas experienced significant changes in both mean and distribution of mean annual temperature estimated using the two-sided Ranksum and KS tests, respectively (figure S7). The majority of urban areas showed changes in both the mean and distribution of mean annual air temperature rather than just the distribution (figure S7). Figure 3 shows changes (linear trend multiplied by the duration (40 years)) in the frequency of extreme windy days. The majority of sites experienced a significant (p < 0.05) decline in the number of extreme windy days during the last four decades ( figure 3(a)). However, in each region a few sites showed increases. About 75% of the total urban areas showed significant declines in number of extreme windy days ( figure 3(b)). On the other hand, statistically significant increases in extreme windy days occurred for less than 10% of the sites. Between 1973 and 2012, median declines in the number of extreme windy days were about 20 days for urban areas with statistically significant changes ( figure 3(c)). In all the regions except SA, results were field significant for declines in extreme windy days with the median ranging between 5 and 25 days during the period of 1973-2012. Our results also showed statistically significant declines in mean wind speed at many sites, which is consistent with trends in mean wind speed at the 925 mb level from the NCEP-NCAR reanalysis (figure S5). Moreover, these results are consistent with the findings of Vautard et al (2010) who reported declines in wind speed in continental areas with higher reductions in strong winds than in weaker winds.

Changes in extreme precipitation
We evaluated changes in the frequency of precipitation extremes for the period of 1973-2012 (figure S22). In contrast to temperature related extremes, changes in precipitation extremes exhibited greater spatial variability (figures S22(a), S2(a), S3(a)). A number of sites in IN and SA showed statistically significant increases in the frequency of precipitation extremes, but most sites elsewhere did not experience statistically significant changes in the number of extreme precipitation events. Taken over all sites, 17% showed a statistically significant increase and less than 5% showed a statistically significant decrease ( figure S22  (b)). Median changes in the number of extremes for sites with significant increases and declines were about the same (5% in both cases, figure S22(c)). Increases in the frequency of extreme precipitation events were field significant for only two regions (IN and SA; figure  S22(d)).
Many urban areas in EU experienced significant declines (estimated using linear trend multiplied by 40) in annual maximum precipitation over the reference period ( figure S2(a)). On the other hand, a few urban areas in IN showed increases in annual maximum precipitation, which is associated with large scale climate variability (Ali et al 2014). More sites had declines in annual maximum precipitation than increases ( figure S2(b)), while about the same number (10%) showed statistically significant increases and decreases. The median over sites in AF, IN, NA, and SA showed positive median changes, whereas EA and EU showed negative median changes in annual maximum precipitation ( figure S2(d)). EU showed a field significant decline in annual maximum precipitation, whereas changes were not field significant for any other region. These results showing declines or increases in annual maximum precipitation are based on a relatively small number of stations primarily located in the urban areas. Therefore, results based on gridded datasets over large areas of EU or IN for instance (e.g. Alaxender et al (2006) and Goswami et al (2006)) may be different from those reported here. Moreover, disparities in trends in different regions can be due to large scale climate variability. For instance, IN receives most of its annual precipitation during the monsoon (JJAS) season, which is strongly associated with sea surface temperature in the Pacific and Indian Ocean regions (Goswami et al 2006, Mishra et al 2012b. On the other hand, European climate is strongly associated with western circulation which transports moist air from Atlantic to the European Land mass (van Ulden and van Oldenborgh 2006).
Changes in the fraction of total precipitation contributed by extreme precipitation events were largely similar to those for annual maximum precipitation (figure S3). Field significant declines occurred for EU, while none of the other regions showed field significance. Precipitation extremes are more variable in space and time than are temperature extremes. For instance, Coumou and Rahmstorf (2012) reported that statistical detection of precipitation extremes remains a challenge due to their non-Gaussian behavior and localized spatial scales. Longer precipitation records in urban areas may provide more insights on changes in extreme precipitation events and their links with climate variability. Figure 4 shows changes in the number of hot days and extreme windy days estimated using linear trend multiplied by the duration of 40 years for urban and non-urban pairs. Both urban and non-urban pairs experienced significant increases in the number of hot days during the last four decades (figures 4(a) and (b)). However, stations located in non-urban areas showed somewhat lower increases in the number of hot days than those located in urban areas. Differences in changes in the number of hot days in urban and non-urban areas due to the urban heat island effect have been reported in many previous studies (Stone 2007, Stewart andOke 2012). Out of the 142 pairs, 72 (62) urban (non-urban) areas showed statistically significant increases in number of hot days (figure 4(e)). A majority of urban and non-urban pairs showed declines in extreme windy days; however, stations located in eastern Asia, EU, and eastern United States showed differences in direction and magnitude of changes (figures 4(c) and (d)). Out of the 142 pairs, statistically significant declines occurred at 88 (64) urban (non-urban) stations, respectively (figure 4(f)). Fortuniak et al (2006) studied urbanrural contrasts in meteorological variables in Poland and found that wind speed at urban stations was lower by 39 and 34% during day and night-time, respectively. Differences in wind speed in urban and nonurban pairs can be associated with land surface conditions as well as localized meteorological conditions as reported in Fortuniak et al (2006) and Vautard et al (2010). Changes in annual maximum precipitation and frequency of precipitation extremes for urban and non-urban pairs are compared in supplemental figure S8. Changes in extreme precipitation for urban and non-urban stations show no statistically significant difference. For instance, mixed changes were observed in annual maximum precipitation and frequency of extreme precipitation events for both urban and nonurban stations located in NA, EU, and EA (supplemental figures S8 (a) and (b)). Moreover, the number of urban/non-urban areas with significant increases/ declines in annual maximum precipitation and frequency of precipitation extremes were similar (supplemental figures S8 (e) and (f)).

Conclusions
Many urban areas across the globe experienced statistically significant increases in the number of heat waves during the reference period 1973-2012. Taken over all sites, the largest number of heat waves has occurred during the most recent decade (2003)(2004)(2005)(2006)(2007)(2008)(2009)(2010)(2011)(2012). Moreover, four of the five years with the largest number of heat waves are the last four years of the record (2009, 2010, 2011, and 2012). Many sites also experienced statistically significant declines in the number of cold waves. Statistically significant increases in the frequency of hot days and nights were detected at many sites, with almost half (48%) having statistically significant increases in the number of hot days, and almost 66% having statistically significant increases in the number of hot nights. Trends related to temperature related extremes are consistent with the findings of Donat et al (2013) who reported significant increases in heat waves and hot nights across the globe. Changes in temperatures extremes are largely driven by changes in both mean and distribution of air temperature (figure S7). Increasing and declining trends in temperature extremes may be associated with natural climate variability, anthropogenic climate warming, and land use/land cover (Kiktev et  Here we simply report the aggregate effect on urban extremes and have not attempted to quantify their separate contributions.
The frequency of extreme windy days declined significantly in a majority of the urban areas, with almost 60% of the sites having statistically significant declines. All regions except SA were field significant for declines in the frequency of extreme windy days. Results from station data for extreme windy days are consistent with the NCEP-NCAR reanalysis data (figure S5). Declining trends in wind speed have been noticed in many studies (Vautard et al 2010, Guo et al 2011, McVicar et al 2012, Troccoli et al 2012 and have been largely attributed to changes in atmospheric circulation and changes in land surface roughness driven by land cover changes. As compared with the prevalence of changes in indices related to extreme temperature and wind, precipitation extremes had statistically significant changes at far fewer sites, with only 17 and 10% of the sites experiencing statistically significant increases in the frequency of extreme precipitation events and annual maximum precipitation, respectively. Many previous studies demonstrated that climate warming will lead to more precipitation extremes (Christensen and Christensen 2003, Gutowski et al 2010, Min et al 2011c, however, consistent with the findings of Donat et al (2013), we found heterogeneous trends associated with extreme precipitation events. There were many urban areas with no trend in extreme precipitation events which highlights challenges in detection of trends related to extremes as described in Coumou and Rahmstorf (2012).
Urban and non-urban pairs showed disparate changes for temperature and wind related extremes (generally more increases in temperature-related extremes, and more decreases in wind-related extremes in urban as compared to non-urban stations), and hence appear to be counter to our overall hypothesis that large scale climate drivers dominate changes in climate extremes. However, precipitationrelated extremes were similar at urban and non-urban pairs, suggesting a prominent role of large scale climate variability, consistent with our overarching hypothesis. However, for a few regions (i.e. EU and EA), differences in land cover and local factors appear to be important contributors to observed changes in wind and temperature related extremes. As compared with results of a similar analysis for the conterminous US (Mishra and Lettenmaier 2011), larger disparities in temperature extremes in urban and non-urban areas can be attributed to record length of the data as well as lesser number of stations available in nonurban regions.
Our results have important implications for policy makers. For instance, increasing numbers of heat waves may lead to enhanced heat wave related mortality in urban areas Robinson 2013, Li et al 2013). Increased warming in urban areas will also have implications on residential heating and cooling demands (Frank 2005). Increases in precipitation extremes in urban areas will pose challenges for urban stormwater infrastructure. These implications argue for the importance of enhancing the density of climate stations in urban and surrounding non-urban areas to provide the baseline data that will be essential for climate change adaptation and decision making.