Seasonal Temperature Identification 1. Indicator Description This indicator describes changes in average seasonal air temperature for the United States from 1896 to 2020. Temperature data for each season (winter, spring, summer, fall) are presented as trends in anomalies. Air temperature is an important indicator of climate, and changes in temperature can have wide-ranging direct and indirect effects on the environment and society. Evidence of a warming climate often begins with average surface temperature (USGCRP, 2017). Components of this indicator include: • Changes in average air temperature by season in the contiguous 48 states over time (Figure 1). • Maps showing the total change in average air temperature by state and by season, across the contiguous 48 states (Figure 2). • A comparison of long-term trends in the contiguous 48 states for all four seasons (Figure 3). 2. Revision History April 2021: Indicator published. Data Sources 3. Data Sources This indicator is based on temperature anomaly data provided by the National Oceanic and Atmospheric Administration's (NOAA's) National Centers for Environmental Information (NCEI), formerly the National Climatic Data Center (NCDC). Specifically, all three figures in this indicator use data from NCEI's nClimDiv data set, which is based on data from the daily version of the Global Historical Climatology Network (GHCN-Daily). GHCN-Daily data undergo more extensive processing by NCEI on a monthly basis for inclusion in nClimDiv. 4. Data Availability The nClimDiv data used in each figure can be accessed online, along with descriptions and metadata. Surface temperature time series data for the contiguous 48 states (Figure 1 and Figure 3) and by state (Figure 2) are based on nClimDiv data that were obtained from NCEI's "Climate at a Glance" web interface (www.ncdc.noaa.gov/cag). For access to underlying nClimDiv data and documentation, see: www.ncdc.noaa.gov/monitoring-references/maps/us-climate-divisions.php. EPA downloaded data in early February 2021, at which point data for all four seasons were available through 2020 (where winter 2020 is the season running from December 2019 through February 2020). Technical Documentation: Seasonal Temperature 1 ------- Methodology 5. Data Collection This indicator is based on surface temperature measurements that were compiled from thousands of weather stations throughout the United States using standard meteorological instruments. All of the stations are overseen by NOAA, and their methods of site selection and quality control have been extensively peer reviewed. Data from these stations for the contiguous 48 states were compiled in the nClimDiv data set, which represents the most complete long-term instrumental data set for analyzing recent climate trends in the contiguous 48 states. The nClimDiv divisional data set incorporates temperature data from GHCN-Daily stations in the contiguous 48 states and Alaska. This data set includes stations that were previously part of the U.S. Historical Climatology Network (USHCN), as well as additional stations that were able to be added to nClimDiv as a result of quality-control adjustments and digitization of paper records. Altogether, nClimDiv incorporates data from more than 10,000 stations. These stations are spread among 357 climate divisions in the contiguous 48 states and Alaska. In addition to incorporating more stations, the nClimDiv data set differs from the USHCN because it incorporates a grid-based computational approach known as climatologically aided interpolation (Willmott and Robeson, 1995), which helps to address topographic variability. Data from individual stations are combined in a grid that covers the entire contiguous 48 states and Alaska with 5-kilometer resolution. These improvements have led to a new data set that maintains the strengths of its predecessor data sets while providing more robust estimates of area averages and long-term trends. The nClimDiv data set is NOAA's official temperature data set for the contiguous 48 states and Alaska, replacing USHCN. The nClimDiv data set starts in 1895. This indicator starts in 1896 for consistency, due to the way winter is defined (see more details in Section 6). To learn more about nClimDiv, see: www.ncdc.noaa.gov/news/ncdc-introduces-national-temperature-index-page and: www.ncdc.noaa.gov/monitoring-references/maps/us-climate-divisions.php. 6. Indicator Derivation General Data Processing Steps NOAA calculated monthly temperature means for each station represented in the nClimDiv data set. In populating the data set, NOAA adjusted the data to remove biases introduced by differences in the time of observation. NOAA also employed a homogenization algorithm to identify and correct for substantial shifts in local-scale data that might reflect changes in instrumentation, station moves, or urbanization effects. These adjustments were performed according to published, peer-reviewed methods. For more information on these quality assurance and error correction procedures, see Section 7. To achieve uniform spatial coverage (i.e., not biased toward areas with a higher concentration of measuring stations), NOAA calculated area-weighted averages of grid-point estimates interpolated from station data. The surface time series for the contiguous 48 states (Figure 1) is based on the nClimDiv gridded data set, which reflects a high-resolution (5-kilometer) interpolated grid that accounts for station density and topography. Results were initially averaged within each climate division, then Technical Documentation: Seasonal Temperature 2 ------- aggregated by state and for the contiguous 48 states as a whole based on an area-weighted average of climate divisions. See: wwwl.ncdc.noaa.gov/pub/data/cirs/climdiv/divisional-readme.txt for more information. In this indicator, the four seasons correspond with the meteorological seasons. Winter corresponds to the months December-February, spring is March-May, summer is June-August, and fall is September- November. NOAA averaged monthly means across each three-month period to derive seasonal temperature data points for each year. Winter data are nominally assigned to the year in which the winter ended; for example, "1896" refers to the consecutive three-month period of December 1895, January 1896, and February 1896. Because the nClimDiv data set does not include data for December 1894, the seasonal temperature time series in this indicator start in the year 1896 to allow for consistency across all four seasons. Figure 1. Average Seasonal Temperatures in the Contiguous 48 States, 1896-2020 Figure 1 presents temperature data across the contiguous 48 states in terms of anomalies for each season. An anomaly represents the difference between an observed value and the corresponding value from a baseline period. This indicator uses a baseline period of 1901 to 2000. The choice of a sufficiently long-term average or baseline period (e.g., at least 30 years) will not substantially affect the shape or the statistical significance of the overall trend in anomalies. For absolute anomalies in degrees, it only moves the trend up or down on the graph in relation to the point defined as "zero." Figure 2. Change in Seasonal Temperatures by State, 1896-2020 The maps in Figure 2 show the long-term change in temperature by state for each season for the period from 1896 to 2020. The slope of each temperature trend was calculated from the three-month climate division anomalies by ordinary least-squares regression and then multiplied by the total number of years (124) to estimate the total change over the period of record. Figure 3. Temperature Change by Season in the Contiguous 48 States, 1896-2020 Figure 3 presents long-term temperature changes in all four seasons. The graph shows the total change for each season, derived from the ordinary least-squares regression slope (°F per year) multiplied by 124 years. Indicator Development NOAA continually refines historical data points in the GHCN and nClimDiv, often as a result of improved methods to reduce bias and exclude erroneous measurements. As EPA updates this indicator to reflect these upgrades, slight changes to some historical data points may become apparent. No attempt has been made to portray data beyond the time and space in which measurements were made. 7. Quality Assurance and Quality Control NCEI's databases have undergone extensive quality assurance procedures to identify errors and biases in the data and either remove these stations from the time series or apply correction factors. The nClimDiv data set follows the USHCN's methods to detect and correct station biases brought on by changes to the station network over time. The transition to a grid-based calculation did not significantly change Technical Documentation: Seasonal Temperature 3 ------- national averages and totals, but it has led to improved historical temperature values in certain regions, particularly regions with extensive topography above the average station elevation—topography that is now being more thoroughly accounted for. An assessment of the major impacts of the transition to nClimDiv can be found at: www.ncdc.noaa.gov/monitoring-references/docs/GrDD-Transition.pdf. Analysis 8. Comparability Over Time and Space Both nClimDiv and the GHCN have undergone extensive testing to identify errors and biases in the data and either remove these stations from the time series or apply scientifically appropriate correction factors to improve the utility of the data. In particular, these corrections address changes in the time-of- day of observation, advances in instrumentation, and station location changes. See Section 7 for documentation. All GHCN-Daily stations are routinely processed through a suite of logical, serial, and spatial quality assurance reviews to identify erroneous observations. For nClimDiv, all such observations were set to "missing" before computing monthly values, which in turn were subjected to additional serial and spatial checks to eliminate residual outliers. Stations having at least 10 years of valid monthly data since 1950 were used in nClimDiv. For temperature, bias adjustments were computed to account for historical changes in observation time, station location, temperature instrumentation, and siting conditions. As with USHCN, the method of Karl et al. (1986) was applied to remove the observation time bias from the COOP network, and the pairwise method of Menne and Williams (2009) was used to address changes in station location and instrumentation in all networks. Because the pairwise method also largely accounts for local, unrepresentative trends that arise from changes in siting conditions, nClimDiv contains no separate adjustment in that regard. For more documentation about nClimDiv, see: wwwl.ncdc.noaa.gov/pub/data/cirs/climdiv/divisional- readme.txt. 9. Data Limitations Factors that may impact the confidence, application, or conclusions drawn from this indicator are as follows: 1. Biases in surface measurements may have occurred as a result of changes over time in instrumentation, measuring procedures (e.g., time of day), and the exposure and location of the instruments. Where possible, data have been adjusted to account for changes in these variables. For more information on these corrections, see Section 8. Some scientists believe that the empirical debiasing models used to adjust the data might themselves introduce non-climatic biases (e.g., Pielke et al., 2007). 2. Uncertainties in surface temperature data increase as one goes back in time, as there are fewer stations early in the record. These uncertainties are not sufficient, however, to mislead the user about fundamental trends in the data (Vose and Menne, 2004). Technical Documentation: Seasonal Temperature 4 ------- 3. Alaska, Hawaii, and U.S. territories are not included in this indicator. Alaska data are included in the nClimDiv data set, but sufficient data are not available before 1925. Hawaii and U.S. territories also have limited data, and nClimDiv spatial interpolation methods have not been developed and validated for these regions. 10. Sources of Uncertainty Error estimates are not readily available for U.S. temperature. Uncertainties in temperature data increase as one goes back in time, as there are fewer stations early in the record. These uncertainties are not sufficient, however, to undermine the fundamental trends in the data. Vose and Menne (2004) suggest that the station density in the U.S. climate network is sufficient to produce a robust spatial average. 11. Sources of Variability Seasonal temperature anomalies naturally vary from location to location and from year to year as a result of normal variations in weather patterns, multi-year climate cycles such as the El Nino-Southern Oscillation and Pacific Decadal Oscillation, and other factors. This indicator minimizes the effect of these factors by presenting a long-term record (more than a century of data) and averaging consistently over time and space. 12. Statistical/Trend Analysis This indicator uses ordinary least-squares regression to calculate the slope of the observed trends in temperature for each season. This slope is the basis for presenting total long-term change in Figures 2 and 3. Table TD-1 shows the resulting regression coefficients and the corresponding p-values (via simple t-test) for the contiguous 48 states, which are significant at the 95 percent confidence level (p < 0.05) in all four seasons. Table TD-1 also summarizes the results of state-level significance testing. All significant state trends showed increasing temperatures. Table TD-1. Linear Regression of Seasonal Average Temperature Data, 1896-2020 Season Contiguous 48 states State-specific trends Slope (°F/year) P-value Number of states with significant change (p < 0.05) Number of states with no significant change Winter 0.023 <0.0001 39 9 Spring 0.015 <0.0001 36 12 Summer 0.011 <0.0001 29 19 Fall 0.011 <0.0001 29 19 Use of linear regression and a simple t-test inevitably requires some assumptions about the shape and underlying characteristics of the data. These tests are used to provide initial approximations about the statistical attributes of this indicator. Technical Documentation: Seasonal Temperature 5 ------- References Karl, T.R., C.N. Williams Jr., P.J. Young, and W.M. Wendland. 1986. A model to estimate the time of observation bias associated with monthly mean maximum, minimum, and mean temperature for the United States. J. Climate Appl. Meteor. 25:145-160. Menne, M.J., and C.N. Williams, Jr. 2009. Homogenization of temperature series via pairwise comparisons. J. Climate 22(7):1700-1717. Pielke, R., J. Nielsen-Gammon, C. Davey, J. Angel, O. Bliss, N. Doesken, M. Cai, S. Fall, D. Niyogi, K. Gallo, R. Hale, K.G. Hubbard, X. Lin, H. Li, and S. Raman. 2007. Documentation of uncertainties and biases associated with surface temperature measurement sites for climate change assessment. B. Am. Meteorol. Soc. 88:913-928. USGCRP (U.S. Global Change Research Program). 2017. Climate science special report: Fourth National Climate Assessment, volume I. Wuebbles, D.J., D.W. Fahey, K.A. Hibbard, D.J. Dokken, B.C. Stewart, and T.K. Maycock (eds.). https://science2017.globalchange.gov. doi:10.7930/J0J964J6 Vose, R.S., and M.J. Menne. 2004. A method to determine station density requirements for climate observing networks. J. Climate 17(15):2961-2971. Willmott, C.J., and S.M. Robeson. 1995. Climatologically aided interpolation (CAI) of terrestrial air temperature. Int. J. Climatol. 15(2):221-229. Technical Documentation: Seasonal Temperature 6 ------- |