Overall Warming with Reduced Seasonality: Temperature Change in New England, USA, 1900–2020

: The ecology, economy, and cultural heritage of New England is grounded in its seasonal climate, and this seasonality is now changing as the world warms due to human activity. This research uses temperature data from the U.S. Historical Climatology Network (USHCN) to analyze annual and seasonal temperature changes in the New England region of the United States from 1900 to 2020 at the regional and state levels. Results show four broad trends: (1) New England and each of the states (annually and seasonally) have warmed considerably between 1900 and 2020; (2) all of the states and the region as a whole show three general periods of change (warming, cooling, and then warming again); (3) the winter season is experiencing the greatest warming; and (4) the minimum temperatures are generally warming more than the average and maximum temperatures, especially since the 1980s. The average annual temperature (analyzed at the 10-year and the ﬁve-year average levels) for every state, and New England as a whole, has increased greater than 1.5 ◦ C from 1900 to 2020. This warming is diminishing the distinctive four-season climate of New England, resulting in changes to the region’s ecology and threatening the rural economies throughout the region.


Introduction
It is now well documented that our world is warming [1][2][3], and much of this warming is due to the emission of greenhouse gases from human activity [4]. Warming from anthropogenic emissions since the start of the industrial revolution, and especially in the past 60 years, will persist for centuries to millennia, and will cause additional long-term changes to our climate system [5]. Rising global temperatures will lead to a variety of detrimental consequences, including the melting of snow, glaciers and sea ice, warming of the oceans, bleaching of coral reefs, rising sea levels, changes in precipitation patterns, increased drought and flooding, heat waves, increased fire activity, increased storm activity, changes in species' physiology, phenology, and distribution to list a few [6,7]. Models show that we will see continuing changes potentially for centuries to come [8,9], based on the amount of greenhouse gases that we have already put into the atmosphere.
To avoid the most catastrophic climate changes in the future, we need to limit future global temperatures. There is now considerable ongoing work trying to determine a threshold of acceptable temperature change. For many years, the long-term goal has been to keep global temperatures from rising 2 • C or more since the start of the Industrial Revolution. The 2009 Copenhagen Accord (COP 15) created an aspirational goal of limiting global temperature increase to 2 • C to avoid dangerous anthropogenic interference with the climate system [10]. Research after the Copenhagen Accord indicates that 2 • C may be too high [11]. The 2015 Paris climate agreement sought to curb greenhouse-gas emissions and limit global temperature increases to between 1.5 and 2 • C [12]. However, climate models now project considerable differences in climate characteristics between the present and showed that all but two stations had an increase in temperature. Through an autocorrelation analysis, they determined that southeastern New England was warming the fastest [29]. Some researchers have combined multiple data sets to analyze temperature change. Burakowski et al. (2008) analyzed temperature data along with snowfall and snow depth data for northeastern U.S. (1965 to 2005) revealing that temperatures (mean, minimum and maximum) increased during this period with the greatest warming occurring in the coldest months of the winter (January and February) [30]. They also discovered decreases in snow cover days and snowfall for the region, reinforcing the warming that was observed. Some researchers have analyzed historical station data to test simulation models of potential future climate change. Hayhoe et al. (2007) used USHCN station data for the northeastern U.S., showing that annual temperatures rose an average of +0.08 ± 0.01 • C per decade in the 20th century and +0.25 ± 0.01 • C per decade for the last three decades of the century [31]. Temperature extremes have also been examined with Brown et al. (2010) measuring climate extremes for the northeastern U.S. from 1870 to 2005 and discovering that temperature indices showed strong warming with an increased frequency of warm events and a decrease in cold events [32]. Thibeault and Seth (2014) examined temperature extremes in the northeastern U.S. and found winters warming about three times faster than summers [33]. There are also reviews of New England's climate literature with the most recent thorough review concluding that New England's seasonality is diminishing as winter temperatures increase [18]. NOAA's National Centers for Environmental Information also provides State Climate Summaries which review a variety of broad temperature changes for each state from 1900 to 2014, showing that in New England, Rhode Island had the greatest increase in annual temperature from the 1900-1904 period to the 2010-2014 period [34].
This research uses USHCN temperature data to analyze how seasonal and annual temperatures have been changing in each of the six New England states as well as at the New England regional level. Our analysis investigates temperature change at the decadal (10-year) and half-decadal (five-year) averages, which past research has not addressed, as well as employing the univariate differencing method which has also not been employed with such data for New England.
This research has two main objectives: 1.
To determine how the annual and seasonal minimum, average, and maximum temperatures have changed during this time period for New England and each of the states; 2.
To determine if New England, and any of the states, have passed the threshold of average temperature change beyond the 1.5 • C or 2 • C thresholds during the time period of 1900-2020.

Materials and Methods
This research analyzed temperature change (minimum, mean, and maximum) between 1900 and 2020 for all New England states and the entire region as a whole for annual and seasonal temperatures. Although some station data start as far back as 1885, 1900 was chosen as a starting date so that most stations would have starting data. Over 77% of the stations in New England have data for 1900 and beyond. Because annual data vary from year to year, and one extreme year can highly skew the data, the data were averaged into five-year and 10-year units. The five-year and 10-year data units were graphed over time as well as differenced. Single-year data were not differenced because of the greater variability of single year of data. Figure 1 outlines the methodology of this research.
Air temperature change is analyzed using the USHCN data set (Version 2.5). Three air temperature data sets (monthly mean minimum, monthly mean average, and monthly mean maximum) were downloaded from the National Centers for Environmental Information (formerly the National Climate Data Center) website (https://www.ncdc.noaa. gov/ushcn/data-access last accessed 6 September 2021). The USHCN stations have been adjusted over time to take into account the validity of extreme outliners, time of observation bias, changes in instrumentations, random relocations of stations, and urban warming biases [35]. Menne et al. (2000) [36] analyzed poorly-sited USHCN stations with good-sited stations and found that adjustments applied to USHCN Version 2 data largely account for the impact of instrument and siting changes. Menne et al. (2010) do note that the adjusted data should not be considered to be completely free from error [36]. However, for New England, Hayhoe et al. (2007) state that the USHCN station data are the best available data for investigating changes since 1900, as the stations have been selected based on the quality and length of data collection and have undergone numerous quality adjustments and validations [31]. In our analysis the quality of the New England USHCN station data were reviewed using the Mann-Kendall Test [37][38][39] and the standard normal homogeneity test [40,41] (details below). Air temperature change is analyzed using the USHCN data set (Version 2.5). Three air temperature data sets (monthly mean minimum, monthly mean average, and monthly mean maximum) were downloaded from the National Centers for Environmental Infor-  There are 44 USHCN stations in New England: Connecticut (CT) (four stations), Maine (ME) (12 stations), Massachusetts (MA) (12 stations), New Hampshire (NH) (five stations), Rhode Island (RI) (three stations), and Vermont (VT) (eight stations) ( Figure 2). Not all of the stations have a complete set of data from 1900 to 2020. The total percentage of missing years for the whole region is 6.7% of the total potential years. For monthly data, 10.8% of data are missing one to five days of observations per month and 2.8% lack six or more days of observation per month.
al. (2010) do note that the adjusted data should not be considered to be completely free from error [36]. However, for New England, Hayhoe et al. (2007) state that the USHCN station data are the best available data for investigating changes since 1900, as the stations have been selected based on the quality and length of data collection and have undergone numerous quality adjustments and validations [31]. In our analysis the quality of the New England USHCN station data were reviewed using the Mann-Kendall Test [37][38][39] and the standard normal homogeneity test [40,41] (details below).
There are 44 USHCN stations in New England: Connecticut (CT) (four stations), Maine (ME) (12 stations), Massachusetts (MA) (12 stations), New Hampshire (NH) (five stations), Rhode Island (RI) (three stations), and Vermont (VT) (eight stations) ( Figure 2). Not all of the stations have a complete set of data from 1900 to 2020. The total percentage of missing years for the whole region is 6.7% of the total potential years. For monthly data, 10.8% of data are missing one to five days of observations per month and 2.8% lack six or more days of observation per month. The three air temperature data sets (minimum, average, and maximum) for the six New England states were downloaded and imported into Microsoft Excel. The data were in monthly mean format for each year. The climate data had four issues that required preprocessing (some stations had missing years of data, some years were missing monthly data, some monthly means were missing days of observations, and data were in integer The three air temperature data sets (minimum, average, and maximum) for the six New England states were downloaded and imported into Microsoft Excel. The data were in monthly mean format for each year. The climate data had four issues that required preprocessing (some stations had missing years of data, some years were missing monthly data, some monthly means were missing days of observations, and data were in integer format in hundredths of degrees Celsius). If a station had a missing year of data, we skipped that year and did not process anything for that year (6.7% of years were missing from the entire data set). For some months, there were one to five days of missing observations for the monthly means (10.8% of total months of the entire data set). In this study, we tolerated one to five missing days of observations per month and used these values with the other monthly means. If there were more than five days of missing observations in a month, then the missing data value (−9999) was a placeholder in the USHCN data. We filled these months with averages for that month determined from two years before and two years after the missing data. If one of the two months before or after were missing, we moved on to the nearest year for that month. Moreover, 2.8% of total months of the entire data set had missing months that were replaced with averages. With so few months and years missing, our averaging of replacement months and skipping of missing years should have a minor influence on the analysis. No state had a disproportionate amount of averaged Climate 2021, 9, 176 6 of 24 or missing data. Temperature values were in hundredths of a degree Celsius and were in integer format. The data were converted into decimal format by multiplying each number by 0.01. All data preprocessing was done using the python programming language.
Once the data were preprocessed, annual, and seasonal (spring, summer, fall, winter) averages and anomalies were created from the following months: Various researchers have organized USHCN data and other temperature data sets in annual and seasonal formats specifically to reduce noise and provide a clearer signal [42][43][44]. We further reduced the noise by creating five-year and 10-year averages for the annual and seasonal data as others have done [34,45,46].
The annual and seasonal anomalies were created using the 30-year base period of 1951-1980 [47], which is near the middle of the analyzed time period. For every station's data set, a long-term average (based on the 30-year period: 1951-1980) was created for the annual and seasonal data. The yearly anomalies were created by subtracting the long-term average from each of the yearly annual and seasonal data. The annual and seasonal anomalies from each state were averaged, providing data for all of New England. The annual and seasonal anomalies from every station for every year were then averaged into five-year (half-decade) and 10-year (decadal) data sets. The annual and seasonal data at the five-year and 10-year levels were then graphed and analyzed. Change over time for the annual and seasonal data were analyzed at the five-year and 10-year levels using a univariate differencing method [48] by subtracting the first five years (1900-1904) from the last five years (2016-2020) and the first 10 years (1900-1909) from the last 10 years (2011-2020). Decadal temperature trends in India have been analyzed with a similar univariate differencing method [46]. The differencing was done at the state and New England regional levels. After the differencing, a t-test (one tailed distribution, sample with equal variance) was run for every result to determine if it was significant at the 95th (p < 0.05) or 99th (p < 0.01) percentile. Annual change (2020 minus 1900) was not analyzed because one year can vary greatly and skew the results. For example, for all of New England at the annual average level, 2020 minus 1900 = +1.87 • C, while 2019 minus 1900 = +0.56 • C. The five-year and 10-year averages reduce the annual variability and the t-tests indicate significance. The standard error was determined and graphed for all analyses and in Appendix A, the standard errors of all decades used in the univariate differencing are presented.
To analyze the quality of the USHCN data, the Mann-Kendall Test (significance level set at 5%) was run on each of the 44 USHCN stations' minimum, average, and maximum annual data as well as on each state's, and New England's, minimum, average, and maximum annual and seasonal data. The Mann-Kendall Test is a non-parametric test, which has no prerequisite conditions on the data to be normally distributed [39]. For the test, a null hypothesis (Ho) indicates that there is no trend; the data are independent and random while the alternative hypothesis (Ha) supposes that there is a trend. The standard normal homogeneity test (significance level set at 5%) was also run on the 44 USCHN stations' annual data. For this test the null hypothesis (Ho) states that the data are homogeneous while the alternative hypothesis (Ha) indicates that there is a change in the data. When the p-value is lower than the 5% significance level, one should reject Ho and accept Ha [40,41,49]. The standard normal homogeneity test can determine where there are non-homogenous breaks in the time series. Concerning sudden changes in the USHCN station data there are a number of possible climatic and non-climatic reasons. One issue about abrupt changes in the USHCN data concerns the replacement of liquid-in-glass thermometers with electronic thermometers at the USHCN sites. To evaluate this potential external force, the changing of thermometers at all 44 stations was investigated and if the timing of the changes occurred at the breaking point indicated by the standard normal homogeneity test, then this potential external error was noted.
The decadal air temperature amplitude of New England's annual minimum and maximum temperature change was created by subtracting the 10-year maximum anomalies from the 10-year minimum anomalies. If values increase, the minimum temperatures are rising fastest and if the trend declines then the maximum values are increasing faster. The GISS Global Surface Temperature data (in monthly and annual means) were downloaded for the years 1900 to 2020 from the NASA Goddard Institute for Space Science (https: //data.giss.nasa.gov/gistemp/ last accessed 6 September 2021). The data came in anomaly format, which used the 1951-1980 base period, the same base period for which we processed the USHCN data [47]. The GISS data were graphed along with the USHCN data to see how New England's temperature is changing compared with the global temperature change.

Results
The results show four broad trends (Tables 1 and 2, Figures 3 and 4). First, New England as a whole, and each of the states (annually and seasonally), have warmed considerably between 1900 and 2020, with a strong warming trend since the 1960s. Second, for all of the states and the region as a whole, there are three general periods of change: warming from 1900 to the 1950s, cooling into the 1960s, and then warming through to the 2010s. Third, the season experiencing the greatest warming is the winter season, with some states warming more than 3 • C (p < 0.01) based on the 10-year analysis, and some states warming more than 4 • C (p < 0.01) based on the five-year analysis. Fourth, the minimum temperatures are generally warming more than the average and maximum temperatures, especially since the 1980s.      For the 10-year (decadal) temperature change (Table 1, Figure 3), 97% of the differencing results were significant at the 95th percentile (p < 0.05), and 82% were significant at the 99th percentile (p < 0.01). The average and minimum annual temperatures for every state, and New England have increased greater than 1.5 • C with all values being significant (p < 0.01), and most of the minimums increased greater than 2 • C. Concerning the seasons, winter warmed the greatest followed by the summer, fall and spring. For the winter, the average winter temperature change was greater than 2 • C for New England and all six states with all values being significant (p < 0.01). Connecticut's winter warmed the most with more than a 3 • C increase. The summer and fall seasons warmed in a similar fashion with spring warming the least. For the summer season, the average temperature changed greater than 1.5 • C for New England and all states, except Maine and Vermont, with all values significant (p < 0.01). For every data point (maximum, average, minimum) all spring temperature changes were less than those of the other seasons, except for Vermont and New Hampshire's maximum spring temperature change. The significance of the spring values were generally lower than the other seasons with the three insignificant figures for the decadal data all occurring in the spring. For the change analysis, there were a total of 105 data points (four seasons and the annual average for six states plus the New England region at the minimum, average, and maximum levels of temperature change, 5 × 7 × 3 = 105). At the decadal level, there was not a single data point that showed a decline in temperature, and 69 out of 105 data points showed an increase of temperature greater than 1.5 • C (66% of points greater than 1.5 • C) and 28 data points greater than 2 • C (27%). All data points greater than 1.5 • C were significant (p < 0.05). For most data points, especially annually and winter, the minimum temperatures increased faster than the average and maximum temperatures, except for Rhode Island, where the maximum temperatures increased the fastest. For the New England region, since the 1990s the minimum temperatures have been increasing significantly faster than the maximum values ( Figure 3). The air temperature amplitude of New England's annual minimum and maximum temperature change shows that since the 1980s the minimum temperature anomalies have been increasing rapidly compared to the maximum temperature anomalies ( Figure 5). This shows that annually in New England the minimum temperatures are rising faster than the maximum temperatures. the spring values were generally lower than the other seasons with the three insignificant figures for the decadal data all occurring in the spring. For the change analysis, there were a total of 105 data points (four seasons and the annual average for six states plus the New England region at the minimum, average, and maximum levels of temperature change, 5 × 7 × 3 = 105). At the decadal level, there was not a single data point that showed a decline in temperature, and 69 out of 105 data points showed an increase of temperature greater than 1.5 °C (66% of points greater than 1.5 °C) and 28 data points greater than 2 °C (27%). All data points greater than 1.5 °C were significant (p < 0.05). For most data points, especially annually and winter, the minimum temperatures increased faster than the average and maximum temperatures, except for Rhode Island, where the maximum temperatures increased the fastest. For the New England region, since the 1990s the minimum temperatures have been increasing significantly faster than the maximum values ( Figure 3). The air temperature amplitude of New England's annual minimum and maximum temperature change shows that since the 1980s the minimum temperature anomalies have been increasing rapidly compared to the maximum temperature anomalies ( Figure 5). This shows that annually in New England the minimum temperatures are rising faster than the maximum temperatures. The five-year data analysis (Table 2, Figure 4) showed similar results as the 10-year analysis. However, the significance level for the results of the five-year analysis was lower than for the 10-year analysis. Moreover, 74% of the results were significant at the 95th percentile and 57% at the 99th percentile. All of the annual and winter values were signif- The five-year data analysis (Table 2, Figure 4) showed similar results as the 10-year analysis. However, the significance level for the results of the five-year analysis was lower than for the 10-year analysis. Moreover, 74% of the results were significant at the 95th percentile and 57% at the 99th percentile. All of the annual and winter values were significant (p < 0.01). The average annual temperature for New England and every state has increased greater than 1.5 • C, and for the southern three states (Connecticut, Massachusetts, and Rhode Island), average annual temperatures rose more than 2 • C. The annual average temperature change for New England was also almost 2 • C (1.99 • C). Like the 10-year data, the winter season warmed the most followed by the summer, fall, and spring. For the winter, the average temperature change was greater than 3 • C for all states and New England, with Connecticut warming the most at 4.73 • C. All of the spring values except for Connecticut's minimum were insignificant, making up 74% of the insignificant values. Like the 10-year results, for most data points, the minimum values were greater than the average and maximum, except for Rhode Island where the maximum values were the highest for all categories.
When looking at the temperature change (minimum, average, and maximum) for New England over the course of the 1900-2020 time period (Figures 3 and 4), there are three general periods of change: warming (1900-1950s), cooling (1950s-1960s), and warming again (1960s-2020). The cooling period appears in the middle of the entire time frame (1900-2020) and divides the changes into two periods of warming. Decadal average annual anomalies reached a peak in the 1950s (0.33) then dropped to a low in the 1960s (-0.33) before rising again (Figure 3). Concerning the five-year average annual anomalies, temperatures reached a peak in the 1950-1954 period (0.60) and declined to a low in the 1965-1969 period (-0.27) before rising again (Figure 4). The minimum, average, and maximum temperature changes at the state levels reflect those of the New England region ( Figure 6).
Concerning data for the entire New England region, we used this inflection point of the 1950s to divide the decadal time period into two 60-year warming periods, from the low of 1900 to the warming peak of 1950s and from the low point of the 1960s to the peak of 2010s (Table 3). These two warming periods were similar, as the winter season warmed the most and the spring season the least. However, the second warming period warmed more than the first. As a percentage of the complete period (1900-2020), the second warming period (1960-2020) accounted for more than 70% of the warming for every season and the annual period except for the summer maximum, which was 49%. Only 8 out of 15 values for the first warming period were significant at the 95th percentile, while for the second warming period, 14 out of 15 values were significant at the 95th percentile, and 13 were significant at the 99th percentile.  with Italic font were at 1.5 °C or above. Numbers with Italic font and Bold were at 2 °C or above. * Significant at the 95th percentile (p < 0.05). ** Significant at the 99th percentile (p < 0.01). Concerning the quality of the USHCN data, the results of the Mann-Kendall Test and the standard normal homogeneity test shows that there are some potential errors in the data and that they may be reducing the amount of warming shown in the results. The Mann-Kendall Test for the 44 stations shows that 25% of the minimum, 14% of the average, and 29% of the maximum station data are not significant (5% level) ( Table 4). When Concerning the quality of the USHCN data, the results of the Mann-Kendall Test and the standard normal homogeneity test shows that there are some potential errors in the data and that they may be reducing the amount of warming shown in the results. The Mann-Kendall Test for the 44 stations shows that 25% of the minimum, 14% of the average, and 29% of the maximum station data are not significant (5% level) ( Table 4). When the insignificant station data are removed from the 10-year and five-year annual data analysis, values for some states increased and others decreased, and for the entire New England region, maximum values increased by 8 to 9% and minimum values rose by 9 to 11% while average values remained almost the same (Table 5). Using all 44 station data, creating annual and seasonal data sets for all states, and New England at the yearly level, the Mann-Kendal analysis only found one data set, Rhode Island's minimum fall data to be insignificant (Table 6). Thus, the combination of all station data did not create multiple insignificant data sets.
x indicates station data that were not significant at the 95th percentile.   The standard normal homogeneity test run on all 44 station's minimum and maximum data sets showed that 84 out of 88 data sets were non-homogeneous with specific breaking points ( Table 7). Of the 84 data sets with breaking points, 6 had breaking points where new thermometers were replaced ( Table 7). Eight of the forty-four stations did not have information about the dates of thermometer changes. Removing these six potentially problematic data sets from the 10-year and five-year annual change analysis created some changes at the state level, but at the New England regional level, changes were less than 2% (Table 8).  When new thermometer equipment occurred on dates where there were breaking points, these dates were flagged as potential problem data. It is clear from this analysis that the New England region is warming in a way similar to much of the world, with an early warming period followed by a brief cooling period and then a return to a warming pattern (Figure 7). New England experienced a greater early rise and a sharper decline than that experienced across the globe. The graph also clearly shows that since the 1960s, New England is warming faster than the world as a whole (Figure 7). The pronounced cooling from the early 1950s through the 1960s is due to the cooling effect of atmospheric aerosols at the global scale, primarily due to humancaused air pollution [50]. The US experienced a greater cooling than the globe as a whole where between the 1930s and 1970s the US cooled 0.5 • C compared to the global cooling of 0.1 • C (48). This mid-century cooling period due to human-caused air pollution has extended longer in parts of Asia such as in northeast India, which had cooling periods due to air pollution into the 2000s [46]. It is clear from this analysis that the New England region is warming in a way similar to much of the world, with an early warming period followed by a brief cooling period and then a return to a warming pattern (Figure 7). New England experienced a greater early rise and a sharper decline than that experienced across the globe. The graph also clearly shows that since the 1960s, New England is warming faster than the world as a whole (Figure 7). The pronounced cooling from the early 1950s through the 1960s is due to the cooling effect of atmospheric aerosols at the global scale, primarily due to humancaused air pollution [50]. The US experienced a greater cooling than the globe as a whole where between the 1930s and 1970s the US cooled 0.5 °C compared to the global cooling of 0.1 °C (48). This mid-century cooling period due to human-caused air pollution has extended longer in parts of Asia such as in northeast India, which had cooling periods due to air pollution into the 2000s [46].  Decadal temperature anomalies for New England and the World, 1900-2020 (New England data from the USHCN data set. World data from the NASA Goddard Institute for Space Science GISS data set. Both data sets created anomalies from the 1951-1980 base period. Error bars based on standard error. For comparative purposes, starting points for both data sets were calibrated to zero by raising each line by the 1900s anomaly (0.313 for the World and 0.645 for New England)).

Discussion
New England appears to be warming faster than the world as a whole. It is clear from the research that New England has warmed past the 1.5 • C level, which the IPCC has set as a do-not-pass threshold for the world [5], and New England is close to passing the 2 • C level. Regions in the higher latitudes, such as New England, are generally warming faster than the world as a whole. It is also clear from the research that, over the past few decades, the colder temperatures (minimum temperatures and winter temperatures) are warming the fastest. This might be a reason why people in New England are not as aware, or are not as concerned, about the warming temperatures, as they would be if the hottest (maximum and summer season) temperatures were warming the fastest. The warming of the cold temperatures reduces negative effects on human thermal comfort whereas if the hottest temperatures warmed faster, it would create more stress on people's thermal comfort [51]. The warming of the winter season found in this research reinforces earlier works, which found similar results [18,[24][25][26][27][28]30,32,33]. In addition to the air temperatures in New England warming faster than much of the world, the neighboring Gulf of Maine is warming faster than most other ocean water bodies [52].
New England's rising temperatures are diminishing the distinct seasonality of the region by vastly reducing winter's cold temperatures as well as increasing temperatures in all of the other seasons. The differences between the four seasons is decreasing. Rising temperatures have resulted in a change in snow cover in the winter. Every decade between 1965 and 2005, New England has lost nine snow-cover days due to less precipitation falling as snow and from the snow melting faster [30]. The loss of snow cover days is occurring throughout much of northeastern North America [53]. Shorter, milder winters will present a variety of challenges for rural industries such as logging, maple syrup harvesting, and the ski industry. Poor road conditions and washouts along with unfrozen forest floors have the potential to limit future logging operations, which need frozen or snow-covered soils [54]. Decreasing snow cover days and warmer winter temperatures are creating the closure of smaller ski areas and driving up the costs at large ski areas [55].
In addition to a loss of snow cover days, warmer winters with rising minimum temperatures will allow insect pests to survive in areas where they could not in the past. Tree pests such as the hemlock woolly adelgid and emerald ash borer have expanded their ranges in New England due to the warmer winters [56,57]. Concerning the spring season, an increase in temperature will shift spring's seasonal time period (beginning and ending earlier) and affect its seasonal environment. Changes to spring temperatures will have a ripple effect on different species. Biologically, the start of spring is measured by the dates of when leaves and blooms are present on diverse species, which is connected to temperature. As spring temperatures start earlier, so blooms, leaves, and insects will arrive earlier as well. This will affect the insect-eating birds, whose migration is connected to the length of day and not as much on temperature. If the insect population has already peaked when the birds arrive, the birds will not have adequate food sources [58,59]. The late winter, early spring maple syrup season is now in flux with changing temperatures with a threat to the viability of the industry in parts of New England [60].
Summer's temperatures are increasing and they are expected to continue to increase with more days above 32.2 • C (90 • F) [18]. As temperatures increase, the summer season will become lengthier in New England. Summer in New Hampshire is projected to start three weeks earlier and last three more weeks into the fall season [61]. In New England, the number of summer droughts are expected to rise due to earlier snowmelt and more runoff from heavier summer rainfall combined with increased evaporation [62]. Concerning the fall season, with both summer and winter temperatures warming, the fall season will arrive later and extend further. Increasing temperatures will potentially cause New England trees to become drought-stressed, which will change their autumn colors. This would have a negative effect on the fall tourists industry [63].
One outlier in the New England analysis concerns Rhode Island, which is different from all of the other states when considering changes in minimum, average, and maximum temperatures. For New England as a whole and all states except for Rhode Island, the minimum temperatures are rising faster than the average temperatures, which are then rising faster than the maximum temperatures. For Rhode Island, the maximum temperatures are rising faster than the average and the minimum. There are a few possible reasons for this. One is that Rhode Island only had three stations and so the state's data may be strongly influenced by one or two stations. The USHCN station on Block Island had a par-ticularly strong increase in maximum temperatures relative to the average and minimum temperatures. The three Rhode Island stations are near the ocean and could be influenced by the Long Island Sound. Future research should look at coastal stations in New England outside of Rhode Island to see if there is a similar pattern. It would also be interesting to see if the coastal stations of Maine have warmed due to the warming of the Gulf of Maine.
Although the USHCN is the best data set to analyze long-term temperature change in New England [31], there are some concerns about the quality of the temperature data. The USHCN data set has many stations across the country that have been collecting temperature data for a long time, some for more than 130 years. Over that time span some stations moved their instruments, changed the types of instruments used, changed the times that they made observations and for some places the land cover around the station changed, such as increasing urban environments. The USHCN data has been processed to remove as many of the inconsistencies as possible [35,36]. One inconsistency for the past 70 years or so is a residual maximum temperature bias where USHCN stations underestimate maximum (and mean) temperature trends [64]. So some of the disparity between the minimum and maximum temperatures found in this research might be due to this bias apparently still in the USHCN data. Another bias of concern is from built-environment encroachment near stations, which can artificially increase temperature readings [65]. In this paper the USHCN data for New England was analyzed with the Mann-Kendall Test and the standard normal homogeneity test with results indicating that there may be some issues with a few of the station time series. The potentially problematic data sets were removed from the analysis and compared with the use of all station data. These comparisons indicate that the problematic station data might be lowering the values showing change. Future research needs to be undertaken to analyze the USHCN data for New England in more detail. This analysis of the region used all USHCN station data so that it can be comparable with other USHCN data studies and to be more conservative with the results, as removing stations tended to increase the warming values. Although there may be some issues with a few of the USHCN stations in New England, the overall USHCN observations in New England fall in line with the sheer number of other indicators, such as earlier lake ice-out dates, decreasing snow cover days, the blooming of plants earlier and earlier, or the increasing number of extreme heat events and heat waves with decreasing extreme cold events and cold waves.
Temperatures in New England will not return to where they were a century ago, or even stay where they are now. One of the fundamentals of our current climate is that the temperatures are no longer stable over long periods of time (decades), but they are changing in a warming trend as the data here show. Our world will continue to warm for some time to come. The hope has been to prevent global temperatures from increasing more than 1.5 • C since the start of the Industrial Revolution [5], to avoid the devastating effects of higher seas, more powerful storms, more droughts, increased flooding, more frequent heat waves, and more fires to name some of the consequences that we are now starting to witness due to our warmer world. A warming world will bring a number of positive (or reinforcing) feedbacks, such as melting more snow and ice reducing the reflection of sunlight and increasing the absorption of solar energy making the world even warmer, which then melts more snow and continues the warming. This is now happening in New England and Northeast North America [53]. This, and other positive feedback scenarios, may potentially accelerate warming across the globe in addition to the warming caused by the continuing, increasing concentration of greenhouse gasses such as carbon dioxide [3]. The increasing global greenhouse gasses and positive feedback loops will continue to warm New England and diminish New England's distinctive fourseason climate.
The results of this research have shown areas in need of future research. The influence of the ocean on the temperature data needs to be investigated by examining near-ocean stations and in-land stations to see if there are significant differences between these two sets of stations on the seasonal and annual data. The amplitude of summer vs. winter should be investigated to determine how these two seasons are changing. The amplitude of seasonal maximum and minimum temperatures at the state level needs to be examined. The reasons for various shifts in temperature change at the regional and state levels should be studied, such as the influence of the North Atlantic Oscillation (NAO) on temperatures in New England. Finally, the quality of the USHCN station data for New England needs to be investigated more deeply than we did here.

Conclusions
This research clearly shows that New England is warming faster than the world average temperature change, with every season experiencing a warming trend. Winters are warming the fastest followed by the summer, fall, and spring seasons. In addition to the winters warming the fastest, the minimum temperatures across the region are warming faster than the average and maximum temperatures, except for Rhode Island where the maximum temperatures are warming the fastest. The result of this warming is diminishing New England's distinctive four-season climate, which will alter the ecology of the region as well as having a negative economic impact, especially for outdoor sports, tourism, and biological natural resource-based economies. The changing climate not only alters the ecology and upends the livelihood of many in New England, it also threatens various infrastructures in New England as well the health and well-being of people through more extreme weather, warmer temperatures, floods and droughts, degradation of air and water quality, expansion of diseases, and sea-level rise.
Based on the data presented here, and the continuing increase of greenhouse gases, it is clear that humanity does not have its hand on the rudder of climate control. It is clear from both the 10-year and five-year analysis that New England's temperature does not show signs of bending towards equilibrium. The science is clear: we are in a climate crisis and we need to take concerted steps to reduce our production of greenhouse gasses as soon as possible. The temperature changes that are currently happening in New England, as outlined in this paper, threaten to disrupt the seasonality of New England, which will disrupt the ecosystems and the economy of New England. Data Availability Statement: All data are available in the public domain. Data can be downloaded from the National Centers for Environmental Information (formerly the National Climate Data Center) website (https://www.ncdc.noaa.gov/ushcn/data-access last accessed 6 September 2021). Processed data may be provided upon reasonable request from the authors. GISS data are available from the NASA Goddard Institute for Space Science (https://data.giss.nasa.gov/gistemp/ last accessed 6 September 2021).