Sunday, July 24, 2022

Coldest Days Trends

Changes in Cold Days

This is a companion post to the Hottest Days Trends analysis. 

Is the frequency of the coldest days of the year changing? The tl;dr answer is of course yes. There are any number of ways to measure this. In this blog post, I identify the frequency of cold days during a base period (1951-2000), and then evaluate the trends of those values. The methodology and results are presented below.

Methodology:

1) Find the temperature that occurs on average 20, 15, 10, and 5 days per year during a baseline period for all long-term climate stations in the U.S. (I chose 1951 to 2000 - 50 years). For a 50-year period, the 10-day per year temperature value is found by ordering all 18,263 low temperatures (assuming no missing obs) and identifying the 500th lowest value.

2) For each year between 1950 and 2021, count the number of days that met or exceeded the temperatures identified in Step 1. Also identify the first and last occurrence of these temperatures each year. For example, if 20F is the temperature that occurred an average of 10 days per year during the 1951-2000 baseline period, count the number of days each year that 20F occurred and find the first/last 20F day each year.

3) If a station contained 65 or more complete years of temperature data (no more than 10 missing days in a year), it was kept for the final analysis. There were 609 U.S. stations that successfully met the criteria, 15 Canadian stations, and 0 Mexican stations. 

4) Fit a linear regression line to the 72-year (1950-2021) time series of the number of annual occurrences (and the first/last dates) for each station. The different between the fitted (not observed) 2021 value and the fitted (not observed) 1950 value represents the change,

Example:

In Fig. 1, we see the change in days per year for the temperature that historically occurred 20 and 10 days per year for Boise, Idaho. Those temperatures are 17F and 1F respectively. A regression line is fitted to each time series and the begin and end points of the regression lines represent the change. In the Boise example, the 17F temperature, which occurred an average of 20 days per year during the 1951-2000 period, trended from 26 days per year to 12 days per year. An decrease of 24 days per year over the 72-year period.

Fig 1. Change in the number of days per year for the temperature that historically (1951-2000) occurred 10 and 20 days per year in Boise, Idaho.


Fig. 2. Change in the length of time (days) between the first/last occurrence of the temperature that historically (1951-2000) occurred 10 and 20 days per year in Boise, Idaho.

In Fig. 2 we see the same information for Boise, except it notes the length of time between the first and last occurrences of the temperature that historically occurred 20 and 10 days per year. A low temperature of 17F historically occurred 20 days per year in Boise. The trend in the length of time between the first and last 17F day of the year decreased from 85 to 57 days per year. We interpret this as a shortening of the coldest time of year by 28 days. Specifically, the typical first occurrence for Boise of 17F changed from November 30th to November 29th (1 day later) and the typical last occurrence of 17F changed from February 21st to January 25th (27 days earlier).

Maps: Modern Days per Year Count

Figs 3-6 below show the current number of days per year for the temperatures that historically occurred 20, 15, 10, and 5 days per year. If you refer back to Fig. 1, the value at the end of the trend (dotted) lines are what we are mapping below. 

Fig 3. Current number of days per year that the historically (1951-2020) coldest 20 days per year now occur. This is the end point of the 72-year regression line.

Fig 4. Current number of days per year that the historically (1951-2020) coldest 15 days per year now occur. This is the end point of the 72-year regression line.

Fig 5. Current number of days per year that the historically (1951-2020) coldest 10 days per year now occur. This is the end point of the 72-year regression line.

Fig 6. Current number of days per year that the historically (1951-2020) coldest 5 days per year now occur. This is the end point of the 72-year regression line.

Maps: Change in the Length of the Cold Portion of the Year

As Fig. 2 shows for Boise, we can track the length of time between the first/last occurrence of particularly hot days and see how that length changed over time. Figs 7-11 below show how the length of time of peak heating has changed. 

Fig 7. Change in the length of time that the historically (1951-2020) coldest 20 days per year now occur. This is the difference between end point and beginning point of the 72-year regression line.

Fig 8. Change in the length of time that the historically (1951-2020) coldest 15 days per year now occur. This is the difference between end point and beginning point of the 72-year regression line.

Fig 9. Change in the length of time that the historically (1951-2020) coldest 10 days per year now occur. This is the difference between end point and beginning point of the 72-year regression line.

Fig 10. Change in the length of time that the historically (1951-2020) coldest 5 days per year now occur. This is the difference between end point and beginning point of the 72-year regression line.

Maps: Change in the Length of Consecutive Cold Days

This is where we will look at the changing length of coldwaves. As with before, we're using the temperature that historically occurred 20, 15, 10, and 5 days per year. Now we are looking at each year and identifying the maximum number of consecutive days that those temperatures occurred. If the 5-day per year temperature occurred 8 times in a specific year but there was at least 1 day between each occurrence, then the maximum consecutive streak would be 1 in that example. Figs. 11-14 show the change in the length of these streaks. The change is the end trend point minus the beginning trend point divided by the beginning trend point times 100. 

Fig 11. Change in the number of consecutive days at or below the temperature that historically (1951-2020) occurred 20 days per year. This is the difference between end point and beginning point of the 72-year regression line as a percentage.

Fig 12. Change in the number of consecutive days at or below the temperature that historically (1951-2020) occurred 15 days per year. This is the difference between end point and beginning point of the 72-year regression line as a percentage.

Fig 13. Change in the number of consecutive days at or below the temperature that historically (1951-2020) occurred 10 days per year. This is the difference between end point and beginning point of the 72-year regression line as a percentage.

Fig 14. Change in the number of consecutive days at or below the temperature that historically (1951-2020) occurred 5 days per year. This is the difference between end point and beginning point of the 72-year regression line as a percentage.

Mapping

The map in Fig. 15 shows the stations used in this analysis. The maps in Figs 3-10 were generated from station data using an inverse distance weighted surfacing algorithm with 50km grid cells. There is often a lot of variability in small geographical regions; therefore, a 5x5 smoothing filter was applied to all maps. No adjustment for topography was made. 

Fig 15. Map of stations used in the analysis presented in this blog post. There are 609 U.S. stations that successfully met the criteria, 15 Canadian stations, and 0 Mexican stations.

Conclusion

The frequency and duration of the coldest days of the year are dramatically decreasing in the U.S. Period. Unlike the heatwave analysis, there are no locations that have escaped this trend.


Saturday, July 16, 2022

Hottest Days Trends

** Update: there is a companion Coldest Days Trends blog post. **


There's a lot of discussion about the increase in heatwaves and hot days in general in a warming world. Unfortunately there is no standard definition of a "heatwave". Many places in the northeastern U.S. describe heatwaves as multiple days in a row with a high temperature above 90F. In the southern Great Plains, days above 100F are more relevant. From a public health perspective, nighttime low temperature may be most meaningful. 

In 2021, the EPA published a widely-shared blog post and map set of the trends in heatwaves from 1961-2019 for the 50 largest cities in the U.S. They defined heatwaves as: 

[A] period of two or more consecutive days when the daily minimum apparent temperature (the actual temperature adjusted for humidity) in a particular city exceeds the 85th percentile of historical July and August temperatures (1981–2010) for that city.

The 85th percentile equates to the 9 highest values per summer (July-August). I have no real objection to their criteria. If I were to critique it though, I would say that most people informally characterize heatwaves as a measure of the peak heat, not the overnight minimum (apparent) temperature. Again, from a public health point of view, nighttime recovery is crucially important. Their metric is more heat stress than heatwave. In my opinion, heatwaves are about high temperatures. Not low temperatures. Not apparent temperatures. In most cases, the very hottest temperatures occur with locally dry airmasses. But again, there's plenty of room to slice and dice this any number of ways.

Changes in Hot Days

In May 2021, I wrote a blog post on the change in days with heavy precipitation. The USA Today newspaper did a full-length write-up on the findings. For this analysis of change in hot days, I used a similar methodology.

Methodology:

1) Find the temperature that occurs on average 20, 15, 10, and 5 days per year during a baseline period for all long-term climate stations in the U.S. (I chose 1951 to 2000 - 50 years). For a 50-year period, the 10-day per year temperature value is found by ordering all 18,263 high temperatures (assuming no missing obs) and identifying the 500th highest value.

2) For each year between 1950 and 2021, count the number of days that met or exceeded the temperatures identified in Step 1. Also identify the first and last occurrence of these temperatures each year. For example, if 90F is the temperature that occurred an average of 10 days per year during the 1951-2000 baseline period, count the number of days each year that 90F occurred and find the first/last 90F day each year.

3) If a station contained 65 or more complete years of temperature data (no more than 10 missing days in a year), it was kept for the final analysis. There were 706 U.S. stations that successfully met the criteria, 18 Canadian stations, and 0 Mexican stations. 

4) Fit a linear regression line to the 72-year (1950-2021) time series of the number of annual occurrences (and the first/last dates) for each station. The different between the fitted (not observed) 2021 value and the fitted (not observed) 1950 value represents the change,

Example:

In Fig. 1, we see the change in days per year for the temperature that historically occurred 20 and 10 days per year for Boise, Idaho. Those temperatures are 96F and 98F respectively. A regression line is fitted to each time series and the begin and end points of the regression lines represent the change. In the Boise example, the 96F temperature, which occurred an average of 20 days per year during the 1951-2000 period, trended from 14 days per year to 34 days per year. An increase of 20 days per year over the 72-year period.

Fig 1. Change in the number of days per year for the temperature that historically (1951-2000) occurred 10 and 20 days per year in Boise, Idaho.

In Fig. 2 we see the same information for Boise, except it notes the length of time between the first and last occurrences of the temperature that historically occurred 20 and 10 days per year. A high temperature of 96F historically occurred 20 days per year in Boise. The trend in the length of time between the first and last 96F day of the year increased from 65 to 89 days per year. We interpret this as a lengthening of the hottest time of year by 24 days. Specifically, the typical first occurrence for Boise of 96F changed from June 24th to June 13th (11 days earlier) and the typical last occurrence of 96F changed from August 26th to September 8th (13 days later).

Fig 2. Change in the length of time (days) between the first/last occurrence of the temperature that historically (1951-2000) occurred 10 and 20 days per year in Boise, Idaho.

Maps: Modern Days per Year Count

Figs 3-6 below show the current number of days per year for the temperatures that historically occurred 20, 15, 10, and 5 days per year. If you refer back to Fig. 1, the value at the end of the trend (dotted) lines are what we are mapping below. 

Analysis

Most readers are familiar with the great heatwaves of the 1930s during the Dust Bowl era. Poor agricultural practices led to a drying of the soils, which meant much more efficient solar heating (higher temperatures). In the 1950s, there were pronounced droughts in the central portion of the U.S. too. These droughts meant more sunshine and greater solar heating efficiency. I note this because the time period we begin with is 1950. This results is a muting of the trends for many regions. If we started in 1960, the change in hot days per year would be more dramatic. That said, starting in 1950 allowed for a maximum number of stations to be included in the analysis; and we should never arbitrarily pick a start date for the analysis to maximize (or minimize) the trends.

The obvious question that arises from looking at the maps are the blue areas in the Great Plains and Mississippi River valley. My assessment is that these regions had large values in the 1950s (flattens the curve) and also have seen a dramatic increase in precipitation during this same time period - which causes less solar energy to be used for warming the surface (more solar energy is used to evaporate water). Conversely, there has been a dramatic drying of the western U.S. This allows for more efficient solar heating.

Let's not lose sight of the fact that the net values are dramatically higher - even with the blue areas. The orange and blue areas do not cancel each other out. Looking at Fig. 6, we now see 7.8 days per year nationwide of the temperature that historically occurred 5 days per year. That is a greater than 50% increase!

Fig 3. Current number of days per year that the historically (1951-2020) hottest 20 days per year now occur. This is the end point of the 72-year regression line.

Fig 4. Current number of days per year that the historically (1951-2020) hottest 15 days per year now occur. This is the end point of the 72-year regression line.

Fig 5. Current number of days per year that the historically (1951-2020) hottest 10 days per year now occur. This is the end point of the 72-year regression line.

Fig 6. Current number of days per year that the historically (1951-2020) hottest 5 days per year now occur. This is the end point of the 72-year regression line.

Maps: Change in the Length of the Hot Portion of the Year

Ultimately I was interested in how the length of summer has changed with respect to peak heating. As Fig. 2 shows for Boise, we can track the length of time between the first/last occurrence of particularly hot days and see how that length changed over time. Figs 7-11 below show how the length of time of peak heating has changed. 

Analysis

In general, the patterns are consistent with the maps showing the change in the number of days (Figs 3-6). I would say that the primary difference is a sharper peak in max summer heating. Using the example of the historically hottest 10 days per year, Fig. 7 shows a nationwide average increase of 44% (now 14.4 days per year vs 10 days per year). The length of time that the historically warmest 10 days per year occurs increased from 59.6 days to 67.7 days - a length increase of 13.6%.

Fig 7. Change in the length of time that the historically (1951-2020) hottest 20 days per year now occur. This is the difference between end point and beginning point of the 72-year regression line.

Fig 8. Change in the length of time that the historically (1951-2020) hottest 15 days per year now occur. This is the difference between end point and beginning point of the 72-year regression line.

Fig 9. Change in the length of time that the historically (1951-2020) hottest 10 days per year now occur. This is the difference between end point and beginning point of the 72-year regression line.

Fig 10. Change in the length of time that the historically (1951-2020) hottest 5 days per year now occur. This is the difference between end point and beginning point of the 72-year regression line.

Maps: Change in the Length of Consecutive Hot Days

This is where we will look at the changing length of heatwaves. As with before, we're using the temperature that historically occurred 20, 15, 10, and 5 days per year. Now we are looking at each year and identifying the maximum number of consecutive days that those temperatures occurred. If the 5-day per year temperature occurred 8 times in a specific year but there was at least 1 day between each occurrence, then the maximum consecutive streak would be 1 in that example. Figs. 11-14 show the change in the length of these streaks. The change is the end trend point minus the beginning trend point divided by the beginning trend point times 100. As with the other maps, there is a decline in the Mississippi River valley, but a marked increase everywhere else. A quick summary is that heat waves are longer in most places, except for the Mississippi River valley.



Fig 11. Change in the number of consecutive days at or above the temperature that historically (1951-2020) occurred 20 days per year. This is the difference between end point and beginning point of the 72-year regression line as a percentage.

Fig 12. Change in the number of consecutive days at or above the temperature that historically (1951-2020) occurred 15 days per year. This is the difference between end point and beginning point of the 72-year regression line as a percentage.

Fig 13. Change in the number of consecutive days at or above the temperature that historically (1951-2020) occurred 10 days per year. This is the difference between end point and beginning point of the 72-year regression line as a percentage.

Fig 14. Change in the number of consecutive days at or above the temperature that historically (1951-2020) occurred 5 days per year. This is the difference between end point and beginning point of the 72-year regression line as a percentage.

Mapping

The map in Fig. 15 shows the stations used in this analysis. The maps in Figs 3-10 were generated from station data using an inverse distance weighted surfacing algorithm with 50km grid cells. There is often a lot of variability in small geographical regions; therefore, a 5x5 smoothing filter was applied to all maps. No adjustment for topography was made. 

Fig 15. Map of stations used in the analysis presented in this blog post. There are 706 U.S. stations that successfully met the criteria, 18 Canadian stations, and 0 Mexican stations.

IPCC Comparison

The maps presented above are broadly consistent with the IPCC analysis of trends in hot extremes (see Fig. 16). They also looked at trends since the 1950s and found parts of central North America were the only region in the world not experiencing an increase in extreme hot temperatures.

Fig 16. IPCC 6 figure (SPM.3) showing the trend in extreme

Conclusion

The frequency and duration of the hottest days of the year are dramatically increasing in the U.S. Period. The rate and magnitude of increase differs across space. Some areas have even seen a decrease. Why does this spatial variation exist? My guess is that the 1950s were an especially warm and dry decade in the central U.S. and a dramatic increase in precipitation in this region in recent decades make achieving hot temperatures especially difficult. Conversely, the western CONUS has dried substantially in recent decades. This has worked to magnify the increase in hot days. Remember that the orange and blue areas on the maps do not cancel each other out.

The blue colors in Figs 3-14 (the analysis maps) unfortunately draw attention away from the strong warming signal in the rest of the maps. It provides "deniers" with a false narrative to latch on to. Such is life. Scientists should always publish data and analysis and accept that bad faith actors will willfully misinterpret. It is incumbent upon us to explain when the bad faith arguments are incorrect and why the people pushing those arguments are bad actors. 

Wednesday, July 6, 2022

Chugach 5,000' Peaks

The Chugach Mountains outside of Anchorage, Alaska, are a vast, unpopulated wilderness full of wildlife, glaciers, and more mountains than you can hike/climb in a lifetime. The portion of Chugach State Park closest to Anchorage boasts a handful of peaks above 5,000'. Specifically, the region south and west of Ship Creek and west of Indian Creek has a dozen (12) peaks that are identified on USGS topographic maps as having a summit more than 5,000' above sea level. Many people have identified these peaks as goals to summit in their hiking excursions. A now famous endurance challenge involves summiting all 12 peaks in a single outing (see: https://www.stockalpine.com/posts/chugach-front-linkup.html ). The routes (and there are several routes) have an elevation gain of ~20,000' and a distance of ~40 miles.

But are there really 12 peaks above 5,000'? Again, these are what the USGS topographic maps indicate. The maps were generated from stereo aerial photographs taken in the late 1950s and early 1960s, and were published in the early to mid 1960s.


The image above is a screenshot of the Anchorage A-7 USGS 1:63,360 topographic map (1960) that shows several of the peaks above 5,000'. It is important to note that the elevations were determined from photogrammetric estimates - not cadastral surveys. The operative work is "estimate". Few, if any, of the highest peaks in this region were actually surveyed. I have been to all of these summits (identified as 5,000'+) and do not recall seeing any survey monuments. The Bureau of Land Management's land survey website indicates no formal survey activities in the area except for section corner and witness points, and the Powerline utility easement.

How accurate are USGS elevation estimates for mountain peaks? My assessment is that they are pretty good, but not perfect. Recall that a few years ago the quite famous Mt. Marathon elevation of 3,022' was determined to be 2,974' when a surveyor decided to test the USGS topographic map elevation estimate. This is a discrepancy of 48'. Let me be clear, USGS topographic maps are the most amazing feat of mapping in the history of humanity. Every American should be in awe of this map set. Full stop. But they are not intended to be a perfect, infallible representation of the landscape. 

Can we do better than the USGS estimates? Yes. We could survey the summits. This is a costly and time consuming process. I imagine that like Mt. Marathon, perhaps someone might do it out of personal interest. What's the next best thing? Aerial or satellite systems that send pulses of energy to the ground and measure the return time. The most common systems are LIDAR (LIght Detection and Ranging) and InSAR (Synthetic Aperture Radar Interferometry). 

In 2010-2017, USGS flew the entire state of Alaska to generate a highly detailed elevation data (see: https://eros.usgs.gov/doi-remote-sensing-activities/2017/usgs/alaska-insar-elevation-data-status ). The gridded data set has 5m horizontal resolution, and from what I can tell has an estimated vertical accuracy of 1.5m. Using this data set, we can interpolate the elevation for all peaks in our area of interest. No longer are we limited to 1950s photogrammetric technology. 

So what does the InSAR data show? It shows that the 12 peaks are all indeed over 5,000' and that Mt. Williwaw is the highest and S. Suicide Peak is the lowest. But it also shows something interesting. It shows that Ptarmigan Peak is over 5,000'. Instead of 12 peaks over 5,000', the InSAR elevation data shows 13 peaks over 5,000'. 


The USGS map shows Ptarmigan Peak is 4,910'. However, when standing on the summit, it clearly is very close to the same elevation of S. Suicide Peak. Several people have reported that the GPS on their phone shows the summit to be well over 5,000', myself included. GPS is certainly more accurate in the horizontal direction than the vertical direction - but the > 5,000' GPS reading is a strong indicator that the topographic map value is far too low in elevation. 

Is the InSAR data to be trusted? An independent measure to cross reference would be nice. It turns out that we have a second data set for comparison. In 2015, the Municipality of Anchorage commissioned a LIDAR elevation mapping survey (also airborne). Unfortunately their area of study did not cover the 5,000' peaks in question; however, there are areas of overlap for many of the high peaks near Eagle River and Girdwood. A sample of peaks in the 4,000' to 5,500' range shows the 2015 LIDAR elevation data is within 3m of the 2010 InSAR elevation data. This is a very strong, independent verification of the InSAR accuracy. Even if the InSAR elevation for Ptarmigan Peak is 5m to 10m too high, it still is safely over 5,000'. 

The National Geodetic Survey has thousands of survey markers around the country that are used as control points for land surveys. In the part of the Chugach Mountains that we are discussing here, there are three monuments above 3,000' that were surveyed in the mid-1980s. One is at Near Point, one is at Rusty Point, and a third is on the McHugh Peak ridge line. Looking at the National Geodetic Survey data for these three sites and comparing it to the InSAR data, the National Geodetic Survey points are 3.8m higher in elevation than the InSAR. This may be a result of different vertical datums used by the National Geodetic Survey (NVD29) and the InSAR data (NVD88). In any event, the InSAR data is actually a little lower than the National Geodetic Survey data, confirming the 5,000' elevation of the peaks in question.

Finally, in 2021, the USGS published a new 1:24,000 topographic map series using the InSAR data (see map below). Those maps clearly show Ptarmigan Peak is over 5,000' (contours are 80'). Unfortunately they do not have Peak elevations noted like the previous maps did. 


Since the underlying data are provided, we can extract the highest grid cell value (5,052'). We expect the true peak elevations to be up to a few feet higher since the horizontal resolution of the underlying grid cells are 5x5 meters and the grid cell average value will not capture the highest elevation. This is why the nice rounded western summit of Ptarmigan Peak shows as 1 foot higher than the narrow, sharp, pointed eastern summit (true summit).

There is now overwhelming evidence that shows  Ptarmigan Peak is over 5,000' in elevation and therefore there are 13 peaks over 5,000' in this portion of the Chugach Mountains.

Postscript: On August 3, 2022, I hiked up the standard route of Ptarmigan with my drone in hand. I know, you are not supposed to have drones in the Park, but this was for science. Drones are very good at keeping a constant elevation. From the main summit, I calibrated the drone to the highest location and then flew it 360 meters to the summit immediately to the east. What I found is that the eastern summit is a 7.3 meters (24 feet) higher and is therefore the true summit of Ptarmigan Peak. In addition, my inReach and cell phone GPS both had elevations comfortably over 5,000'. For consistency, I will not edit the table or the map of the InSAR elevations since those reflect their true data source.



Google Earth File

Here is a Google Earth file that I made in GIS using the InSAR DEM and identifying peaks. It only covers the area south/west of Ship Creek and Indian Creek. I only left the peaks that have a page on the peakbagger.com site that are over 3,000' in elevation in InSAR. 



Tuesday, June 7, 2022

10F Increment Low Temperature Days per Year

This is the companion to the 10°F high temperature post. It shows the number of days per year with low temperatures in each 10F temperature range.