Letters to Nature
Nature 423, 528 - 531 (29 May 2003); doi:10.1038/nature01675

There is an corrigendum (04 September 2003) associated with this Letters to Nature. Please click here to view.

Impact of urbanization and land-use change on climate


University of Maryland, College Park, Maryland 20770-2425, USA

Correspondence and requests for materials should be addressed to E.K. (ekalnay@atmos.umd.edu).

The most important anthropogenic influences on climate are the emission of greenhouse gases1 and changes in land use, such as urbanization and agriculture2. But it has been difficult to separate these two influences because both tend to increase the daily mean surface temperature3, 4. The impact of urbanization has been estimated by comparing observations in cities with those in surrounding rural areas, but the results differ significantly depending on whether population data5 or satellite measurements of night light6-8 are used to classify urban and rural areas7, 8. Here we use the difference between trends in observed surface temperatures in the continental United States and the corresponding trends in a reconstruction of surface temperatures determined from a reanalysis of global weather over the past 50 years, which is insensitive to surface observations, to estimate the impact of land-use changes on surface warming. Our results suggest that half of the observed decrease in diurnal temperature range is due to urban and other land-use changes. Moreover, our estimate of 0.27 C mean surface warming per century due to land-use changes is at least twice as high as previous estimates based on urbanization alone7, 8.

Two methods used in the US to classify meteorological stations into urban and rural to 'correct' the observed surface temperature trends for urbanization effects are based on population data5 and satellite measurements of night-light6-8, respectively, and the corresponding estimates of the impact of urbanization differ in magnitude (0.06 and 0.15 C per century)7, 8. The finding that atmospheric temperatures as measured by satellites and weather balloons have smaller warming trends than surface observations has been the subject of much discussion9 centred mostly on the quality of the data, but it could be partially explained by a predominance of land-use effects over greenhouse warming near the surface.

We estimated the impact of urbanization and other land uses on climate change by comparing trends observed by surface stations with surface temperatures derived from the NCEP-NCAR 50-year Reanalysis (NNR)10. In the NNR (a statistical combination of 6-hour forecasts and observations), surface observations of temperature, moisture and wind over land are not used11. However, atmospheric vertical soundings of wind and temperature (rawinsondes and satellite soundings) strongly influence the NNR, and surface temperatures are estimated from the atmospheric values. As a result, the NNR should not be sensitive to urbanization or land-use effects, although it will show climate changes to the extent that they affect the observations above the surface.

As indicated by Fig. 1 and many other studies, the NNR captures well surface temperature variations caused by atmospheric storms, advection of warm/cold air, and variations in the frequency or track of major storms. In contrast to the actual surface observations, we find no statistically significant difference in the NNR estimation of urban and rural station trends (see Methods). These arguments suggest that we could attribute the differences between monthly or annually averaged surface-temperature trends derived from observations and from the NNR primarily to urbanization and other changes in land use.

Figure 1 Comparison of monthly mean station and NNR surface temperature anomalies with respect to their annual cycles for the city of Baltimore, Maryland, USA.   Full legend
High resolution image and legend (185k)

We compare the daily maximum and minimum temperatures of 1,982 surface stations located below 500 m in the 48 contiguous United States, and the daily surface maximum and minimum temperatures on a 2.5 gaussian grid from the NNR interpolated to the station locations, both for the period 1950–1999. We compute temperature anomalies with respect to the 50-year mean annual cycle for each site and each data set. Trends are computed as changes in decadal averages in the anomalies to reduce random errors. The NNR (1948 to the present) has been constructed with a model and data assimilation system kept unchanged, but it is affected by changes in the observing systems, especially the introduction of the satellite observing system in 1979. Therefore, in the computation of trends we exclude changes from the decade of the 1970s to the 1980s.

Figure 1 compares time series of 50 years of monthly mean temperature anomalies for Baltimore, a large city in Maryland, including the average decadal difference between observations and the NNR. There is good agreement in the interannual variability, with a correlation of over 0.9, but also a growing trend in the difference between the surface observations and NNR, increasing to 1.4 C during the 1990s, a difference we attribute to urbanization and other surface changes that do not affect the NNR. A similar analysis on all surface stations (Supplementary Fig. 1) indicates a 50-year correlation of about 0.9 everywhere except in mountainous regions, where it is between 0.4 and 0.7, which is why we only include stations located below 500 m. The correlation is also lower in the west coast, possibly owing to the proximity of mountains or to low data density in the Pacific Ocean. Decadal trends can be locally dominated by interannual and decadal variability of the temperature due to anomalies in the circulation rather than to land use change—effects that are excluded by taking the differences between surface and NNR temperatures.

The decadal trend averaged over the two separate 20-year periods (1980–1999, and 1960–1979) is computed for every station and averaged in boxes of 0.5 latitude by 0.5 longitude, with an overall average computed over all the boxes. Figures 2 and 3 show the average trends for the observations and from the NNR, and also the difference between these two trends, which is at least partially attributable to changes in use of the land surface.

Figure 2 Decadal trend of the maximum temperature averaged for every US station below an elevation of 500 m.   Full legend
High resolution image and legend (168k)

Figure 3 Decadal trend of the minimum temperature averaged for every US station below 500 m.   Full legend
High resolution image and legend (168k)

The maximum temperature (Fig. 2) shows a warming trend in the observations in the eastern and western US and a cooling trend in the Midwest, with a slightly negative overall average of -0.017 C per decade. The NNR is similar but smoother, with an average of +0.008 C per decade. The difference between the observed and NNR trends is somewhat negative in most of the country east of the Rockies, but is strongly positive in California and to a lesser extent, in Oregon and Washington, with an average difference of -0.025 C per decade.

The minimum temperature (Fig. 3) observations show a much stronger positive trend in most of the country, with an average of +0.193 C per decade. In the NNR, the minimum temperature increases everywhere except in the Midwest and California, with an average of +0.113 C per decade. The difference in minimum temperature trends between observed and NNR values is positive in most of the country, especially in California, with an average of 0.080 C per decade (40% of the observed trend).

Supplementary Fig. 2 shows the trend in the difference between maximum and minimum temperatures or diurnal temperature range (DTR). In the observations the DTR trend is strongly negative in most of the country with an average decrease of -0.210 C per decade. The NNR also shows a general decrease of DTR, with a national average of -0.105 C per decade so that about half of the decrease in DTR could be attributable to surface changes.

The daily mean temperature observation trends (Supplementary Fig. 3) obtained as the average of the maximum and minimum temperatures show an increase in most of the country, with an average trend of +0.088 C per decade. The NNR trends have an average of +0.061 C per decade. Of the two "urban correction" estimates7, 8, our estimate of 0.27 C per century attributable to land use is closer to the estimate based on the night-light urban effect (see centre bottom of plate 3C in ref. 7) than to the estimate based on population density (see centre bottom of plate 3B in ref. 7). It should be noted that our observed daily mean temperature trends (Supplementary Fig. 3a) are different from previous 50-year or 100-year trend estimates (see plates 3A and 7A in ref. 7) because in our computations we did not include (1) the decadal trends corresponding to the 1980s–1970s and especially the 1960s–1950s, and (2) urban and non-urban data adjustments. The non-urban adjustments tend to be strongly positive except over the Rockies (see plate 3B in ref. 7), so that if we had added them to the raw observations, our estimate of the land-use impact on the mean temperature trends would have been geographically similar but larger.

Although it is not possible definitively to attribute the differences between the observation and the NNR temperature trends solely to land use, including urbanization, agriculture and irrigation, our results are compatible with such an interpretation. The well-known 'urban heat island' effect actually takes place at night, when buildings and streets release the solar heating absorbed during the day. At the time of the maximum temperature the urban effect is one of slight cooling, owing to shading, aerosols, and to thermal inertia differences between city and country that are not currently well understood12.

The effect of agricultural development, increasing evaporation during the day, would also tend to decrease the maximum temperature: irrigation would increase the heat capacity of the soil, thus increasing the minimum temperature. Therefore, both urbanization and agriculture effects could be consistent with the general increase in the minimum temperature and slight decrease in the maximum temperature, and contribute to the reduction in the diurnal temperature range shown in our estimates east of the Rockies (Figs 2c and 3c).

This implies that the comparison of urban and rural stations without including agricultural effects would underestimate the total impact of land-use changes. More studies are needed, including a comparison of geographical distribution of NNR trends with other upper-air observations, such as rawinsondes and satellites, a more precise definition of the urban and rural observing stations, and the impact of other human activities such as contrails and aerosols that can also reduce the diurnal temperature range13.

Our method can incorporate updated observations as they become available, can be applied to land stations throughout the world, to other variables such as humidity and winds, detect seasonal trends, and signal changes in station locations that are otherwise difficult to identify.

Data For the surface observations, we use the daily surface maximum and minimum uncorrected surface station temperatures from the National Climate Data Center (NCDC) 'Cooperative Summary of the Day' data set over the 48 contiguous states of the United States for 1950–1999. For the NNR, we use the global daily surface maximum and minimum temperatures gridded on 2.5 gaussian boxes, also for the period 1950–1999.

Analysis We interpolate linearly the gridded NNR data to each observational site, and only consider the sites that have a total of at least 480 (whole) months of observations. In addition, because the NNR has surface heights different from those of the real locations, and extrapolations underground can introduce errors overwhelming the signal of the real trends (Supplementary Fig. 2), in the computation of the trends we only consider sites with elevations lower than 500 m. There are 1,982 US surface stations satisfying these two conditions. We obtain monthly means by averaging daily data; daily mean temperatures are obtained by averaging maximum and minimum temperatures, and daily temperature ranges by subtracting the minimum from the maximum temperature.

Because the NNR can have systematic differences with observations, especially near the surface, owing to deficiencies in the model forecast or the method of assimilation, we remove the 50-year monthly mean annual cycle for each site from both the observations and the NNR. We are thus comparing anomalies with respect to the 50-year mean annual cycle. In the results we present both comparisons of the 50-year time series and trends. The trends are computed as changes in decadal averages in order to reduce random errors. We only consider two decadal trends: the decade 1990–1999 minus 1980–1989, and 1970–1979 minus 1960–1969. We do not include in the trends the difference between the decades 1960–1969 and 1950–1959, because the observing system during the 1950s was considerably less reliable than in later decades, and it underwent significant scheduling changes during 1958 (ref. 11).

In addition, we have to address changes in the observing systems, especially the introduction of the satellite observing system (of which the most important is the TIROS-N Operational Vertical Sounder, TOVS) starting in 1979. These two major changes are the main reason why trends in the NNR need to be carefully estimated. We therefore do not include the changes 1980–1989 minus 1970–1979. The two decadal changes that we keep correspond to the 1990s minus 1980s (20 years with satellite data), and 1970s minus 1960s (20 years essentially without satellite data). Thus, when we average them we obtain decadal trends from two independent and largely homogeneous 20-year periods.

We compared the 1990s versus 1980s trend of 775 stations classified as urban versus 167 stations classified as rural. The mean surface temperature increased by 0.31 C for the urban stations and 0.13 C for the rural stations, with standard deviations of about 0.5 C each. The difference between urban and rural warming, 0.18 C, is significant at a 99% level of significance. The trends for the reanalysis station estimates are 0.26 C for urban and 0.25 C for rural, with standard deviations of about 0.22 C, and the difference 0.01 C between urban and rural is insignificant, showing that the NNR is insensitive to surface effects.

In the time series we compute the 1950–1959 average temperature difference between the NNR and the surface station at each station and subtract it from the NNR. This forces the two time series to have the same 10-year time average during the 1950s and is done for display but does not affect the computation of the trends or correlations.

Supplementary information accompanies this paper.

Received 18 December 2002;accepted 23 April 2003


1. IPCC Climate Change 2001: The Scientific Basis (Cambridge Univ. Press, Cambridge, UK, 2001)
2. Pielke, R. A. Sr et al. The influence of land-use change and landscape dynamics on the climate system: Relevance to climate-change policy beyond the radiative effects of greenhouse gases. Phil. Trans. R. Soc. Lond. A 360, 1-15 (2002)
3. Gallo, K. P. & Owen, T. W. Satellite-based adjustments for urban heat island temperature bias. J. Appl. Meteorol. 38, 806-813 (1999) | Article | ISI |
4. Owen, T. W., Gallo, K. P., Elvidge, C. D. & Baugh, K. E. Using DMSP-OLS light frequency data to categorize urban environments associated with US climate observing stations. Int. J. Remote Sensing 19, 3451-3456 (1998) | Article | ISI |
5. Easterling, D. R. et al. Maximum and minimum temperature trends for the globe. Science 277, 364-367 (1997) | Article | ISI | ChemPort |
6. Gallo, K. P., Easterling, D. R. & Peterson, T. C. The influence of land use/land cover on climatological values of the diurnal temperature range. J. Clim. 9(11), 2941-2944 (1996) | Article | ISI |
7. Hansen, J. E. et al. A closer look at United States and global surface temperature change. J. Geophys. Res. 106, 23947-23963 (2001) | ISI |
8. Gallo, K. P., Owen, T. W., Easterling, D. R. & Jamason, P. F. Temperature trends of the US historical climatology network based on satellite-designated land use/land cover. J. Clim. 12(5), 1344-1348 (1999) | Article | ISI |
9. NRC Board on Atmosphere Sciences and Climate Committee Panel, John M. Wallace, chair, pages 1-71. Reconciling Observations of Global Temperature Change (National Research Council, Washington DC, 2000).
10. Kalnay, E. et al. The NCEP/NCAR 40-year reanalysis project. Bull. Am. Meteorol. Soc. 77, 437-431 (1996) | Article | ISI |
11. Kistler, R. et al. The NCEP/NCAR 50-year reanalysis: Monthly means CD-ROM and documentation. Bull. Am. Meteorol. Soc. 82, 247-267 (2000) | Article | ISI |
12. Runnalls, K. E. & Oke, T. R. Dynamics and controls of the near-surface heat island of Vancouver, BC. Phys. Geogr. 21, 283-304 (2000) | ISI |
13. Travis, D. J., Carleton, A. M. & Lauritsen, R. G. Contrails reduce daily temperature range. Nature 418, 601 (2002) | Article | PubMed | ISI | ChemPort |

Acknowledgements. This study was partially supported by a Risk Prediction Initiative grant. We are grateful to A. Senserini, who performed most of the computations, to J. E. Janowiak and W. Ebisuzaki who provided the data, and to R. Murnane, T. Oke, J. Hansen, E. Rassmusson, R. Pielke Sr, T. vonder Haar and Z. Li for discussions.

Competing interests statement. The authors declare that they have no competing financial interests.

2003 Nature Publishing Group
Privacy Policy