Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Predicting climate-change induced heat-related illness risk in Grand Canyon National Park visitors

  • Danielle E. Buttke ,

    Roles Conceptualization, Data curation, Formal analysis, Methodology, Writing – original draft, Writing – review & editing

    danielle_buttke@nps.gov

    Affiliation Natural Resource Stewardship and Science Directorate, National Park Service, Fort Collins, Colorado, United States of America

  • Brinkley Raynor,

    Roles Formal analysis, Validation, Visualization, Writing – review & editing

    Affiliation Department of Biostatistics, Epidemiology, and Informatics, University of Pennsylvania, Philadelphia, Pennsylvania, United States of America

  • Gregor W. Schuurman

    Roles Conceptualization, Data curation, Project administration, Writing – original draft, Writing – review & editing

    Affiliation Natural Resource Stewardship and Science Directorate, National Park Service, Fort Collins, Colorado, United States of America

Abstract

Background

The climate crisis is the greatest public health threat of the 21st century. Excessive heat is responsible for more deaths than any other extreme weather event, and the frequency, intensity, and duration of extreme heat events are increasing globally due to climate change. Exposure to excessive heat can result in heat related illnesses (HRIs) and long-term poor health outcomes. Physical exertion, sudden exposure to excessive heat, and the lack of physical or behavioral adaptation resources are all associated with greater HRI risk, which is expected to increase for visitors to Grand Canyon National Park (GCNP) and other public lands as climate change worsens.

Objectives

Our objectives were to understand 1) the relationship between weather and HRI in GCNP visitors, 2) how future HRI rates may change, and 3) how land management agencies can update risk mitigation strategies to match changing risk and better manage an increased HRI burden.

Methods

We utilized previously published data on HRI in GCNP visitors, and records of daily visitation, temperatures, and maximum and minimum daily humidity from the same study period to develop a model estimate for HRI risk. We then used future climate projections from the World Climate Research Programme’s Coupled Model Intercomparison Project phase 5 multi-model dataset to model future HRI risk under different climate scenarios.

Results

The incidence of HRI was significantly associated with maximum daily temperature and minimum relative humidity, and was more common in the shoulder season months. We estimated that HRI will increase 29%-137% over 2004–2009 levels through 2100, assuming no change in visitation.

Discussion

Climate change will continue to increase HRI risk for GCNP visitors and poses risks to public land managers’ mission to provide for safe recreation experiences for the benefit of this and future generations in places like GCNP. Excessive risk during the shoulder season months presents an opportunity to increase preventative search and rescue and education efforts to mitigate increased risk.

Introduction

Climate change is the greatest public health threat of the 21st century [1]. Negative public health implications of climate change include reduced air quality and access to clean water, increasing malnutrition, vector-borne, and diarrheal diseases, increasing civil conflict and migration, and negatively impacted mental health [2, 3]. Excessive heat is responsible for more deaths than any other extreme weather event, and the frequency, intensity, and duration of extreme heat events are increasing globally due to climate change [2, 4, 5]. Heat-related illness (HRI) includes acute conditions of heat cramps, heat exhaustion, and heat stroke, but can also exacerbate underlying health conditions such as respiratory or heart disease, diabetes, or cerebrovascular disease [68]. High ambient temperatures have significant negative impacts on mental health, including a significant increase in suicide risk and violent crime, and can even hinder learning ability, productivity, and cognitive performance [912]. Those who recover from HRI face increased risk of cardiovascular and chronic disease and early death [1316]. Like many other aspects of climate change, the burden of excessive heat falls disproportionately on disadvantaged and minority communities, who also have fewer adaptation and recovery resources [5, 17, 18].

National parks and other public lands are important public health resources where people increasingly seek the physical-, mental-, and social-health benefits of nature [19]. However, the nature of people’s use of public lands also poses distinct climate change-related health risks because 1) visitors often engage in more strenuous exertion and experience greater exposure to the elements with fewer resources than experienced in daily life, and 2) these areas are increasingly experiencing unprecedented climate and weather conditions that may surprise experienced visitors and even park staff [17, 20, 21]. Warming temperatures, including increasing frequency of very hot days, may be the most consistent expression of ongoing climate change across the U.S. National Park System—a vast array of over 400 management units centered in North America but spanning the globe from the south and west Pacific to the Arctic and the Caribbean. Most U.S. national parks are at or beyond the extreme warm edge of their historical ranges of variability [21]. Unanticipated heat during historically cooler months and in park units with historically milder climates increases the risk of HRI for unprepared visitors. Parks such as Grand Canyon National Park (GCNP) with a history of HRI are therefore increasingly modifying their approaches in light of recent increases in rates, and parks such as Yosemite National Park in California and New River Gorge National Park and Preserve in West Virginia that have not historically encountered the issue are developing/adopting mitigation approaches (e.g., [22]).

Providing safe and healthy visitor experiences is an important part of the U.S. National Park Service (NPS) mission. Understanding conditions associated with increased HRI is critical for anticipating and managing this risk, both in the near term and in longer climate change adaptation planning and management of visitor services. Common NPS visitor HRI symptoms include heat exhaustion, heat stroke, and hyponatremia (low blood sodium), and these events result in visitor hospitalizations and deaths every year [20, 23]. Although the creation of a preventative search and rescue (SAR) program in GCNP in 1997 substantially improved health outcomes and reduced emergency medical services (EMS) HRI responses, staffing for this program is limited to the peak visitation season, which may not mirror peak HRI risk [24, 25]. Additionally, most public lands lack these resources and do not have a data-informed model to appropriately target limited resources to best affect public health outcomes. As temperatures rise and patterns of extreme heat become more variable, preventative SAR resources may not be ideally structured to address changing risk patterns nor increasing demand, and GCNP park EMS and preventative SAR services are again at risk of becoming overwhelmed by demand. At the same time, public lands lacking these HRI mitigation resources would benefit from greater understanding of how these programs might be developed to both understand and address peak risk intervals.

By integrating historical NPS HRI data and current NPS mitigation protocols and climate projections, we aim to better understand 1) the relationship between weather, human behavior, and HRI, 2) how future HRI rates may change, and 3) how land management agencies can update risk mitigation strategies to match changing risk and better manage an increased HRI burden. This information will allow the development of adaptive SAR and education programming targeting the highest risk conditions to improve public health outcomes as well as improve public land managers’ understanding of how climate change threatens their mission to provide for safe recreation experiences for the benefit of this and future generations.

Materials and methods

Study site

Our study site, GCNP provides a unique opportunity to understand HRI in a national-park context because its well documented history of individual HRI incidents and reliable information about visitor numbers provide a unique opportunity to understand HRI in a national-park context. In 2018, visitor numbers to GCNP reached a peak of 6.38 million people with most visits occurring during the peak season months [26]. Every year, park visitors unprepared for the intense exertion and high temperatures within the canyon suffer from HRI events; particularly because conditions at the top of the canyon where visitors begin their hikes are often much cooler than conditions at the bottom of the canyon, where visitors also transition from hiking downslope to having to hike up a considerable incline (Fig 1; [20, 23, 24]).

thumbnail
Fig 1. Map of the study site with Grand Canyon National Park boundaries depicted in white.

Park boundary data are publicly available through NPS IRMA [26] and base map data is publicly available [27].

https://doi.org/10.1371/journal.pone.0288812.g001

Data

All data used on our study are publicly available either online or provided by the study authors of a previous publication with permission granted for use and publication. A previously anonymized and published dataset of GCNP HRI case reports from 2004–2009 [23] was provided by the study authors and used to build the outcome of interest (Table 1). Both paper and electronic medical records were reviewed by a trained public health officer and records meeting the case definition were extracted for the study. HRI events were defined as:

“Persons whom EMS treated (or found if deceased) were eligible for inclusion into the study if they were a visitor during the study period and had a documented outdoor heat exposure of more than 1 hour within the previous 24 hours. To be a case there need to be one of the following documented indicator(s) of HRI: 1) paramedic assessment of HRI; 2) self-report or phrases that indicated “heat” (e.g., “overheated, dehydration”) as a contributing factor; or 3) at least 2 signs and symptoms consistent with HRI not associated with another etiology (e.g., nausea and vomiting while hiking)” [23].

thumbnail
Table 1. Data type, purpose, and source used in a predictive model of heat illness risk in visitors to Grand Canyon National Park, 2004–2009.

https://doi.org/10.1371/journal.pone.0288812.t001

Only records with complete data on the date of the event, paramedic assessment or self-reported symptoms, location, and reported activity were included in the study.

We estimated daily visitor numbers by first summing estimates of automobile-based visitor numbers from the park’s three automobile entry points—the South Rim, North Rim and Desert View visitor stations [28]—and then adding to this sum an estimate of daily visitors arriving by other means. The vast majority of park visitors (>75%, NPS [unpublished]) arrive in automobiles, for which the NPS provides daily visitor estimates. However, small proportions arrive by other means, including by bus (<20%), train (<5%), and boat (<<1%), for which only monthly counts exist. For this study, we distributed these monthly counts evenly across the days in the month, adding them to the automobile-based visitor-count estimate for each day to obtain a daily total park visitor estimate. NPS social science program staff suggested this approach because most such visitors are international travelers, a group whose visitation is not strongly influenced by day of the week or month.

Historical climate data, including daily minimum temperature (Tmin), maximum temperature (Tmax), minimum relative humidity (RHmin), and maximum relative humidity (RH max), were extracted from the Abatzoglou gridded daily surface meteorological (gridMET) dataset, for the period from 1979 to 2018 [29]. We used data from the 4-km gridMET grid cell that includes Grand Canyon Village (Lat 36.054390/Lon -112.140143) in the park’s heavily visited South Rim District. We calculated heat index using the standard methods described by the U.S. National Oceanographic and Atmospheric Administration [30]. We chose daily Tmax, daily RHmin, and time of year (month) as our predictor variables based on a strong correlation between these indices and HRI and death [6].

To characterize how HRI incidence might change in the park through the rest of the 21st century, we obtained daily Tmax and RHmin data (projections) for the period 2018–2098 from 14 global climate models (GCMs) in the World Climate Research Programme’s CMIP5 (Coupled Model Intercomparison Project phase 5) multi-model dataset [31], statistically downscaled using the MACA (Multivariate Adaptive Constructed Analogs) method [32]. The MACA method uses training data from gridMET and therefore MACA-based projections are compatible with gridMET records. We obtained these data for the visitation season of interest (April 1-Sept 30) from the 4-km MACA grid cell that includes Grand Canyon Village (Lat 36.054390/Lon -112.140143) in the South Rim District. We worked with climate uncertainty by considering a range of plausible climate futures [33, 34] represented by 28 projections. These 28 projections consisted of two simulations each of 14 downscaled CMIP5 global climate models, in which one simulation used a moderate greenhouse gas emissions pathway that assumes lower future emissions rates due to technological advancements and policy change (Representative Concentration Pathway [RCP] 4.5) and the other used a high emissions pathway (RCP 8.5) [35]. Pathways of GHG emissions have been standardized as a set of Representative Concentration Pathways (RCPs; [35]) that reflect the climate consequences of a broad range of socio-economic futures and have been accepted by the international community. A large set of GCMs, reflecting different quantitative representations of atmospheric dynamics and different sets of Earth system processes, have used these standardized RCPs to generate climate projections through the twenty-first century.

Data cleaning

We compared the original dataset to preventative SAR reports from the GCNP Visitor and Resource Protection Office and compared deaths in our data set with public records to ensure completeness. We imputed missing visitor numbers from any individual visitor station on a given day using the MICE package in R and the CART method due to the high degree of correlation between visitor numbers from all other stations and date, with daily visitor numbers from all other visitor stations and day of the year as predictors [36].

Our primary outcome was the number of HRIs per week per estimated 100,000 visitors during the peak season. The peak season here is defined as April 1 through September 30. Therefore, we binned the data (structured by date) into weeks, with weeks starting on Sunday and ending on Saturday. Recorded HRIs and recorded visitation were summed per week, and temperature and humidity data were averaged across the week. Month was assigned per week based on which month the majority of the week fell into.

Regression model development

We built general linear regressions predicting the number of heat events per estimated 100,000 visitors per week based on climate variables including Tmax, RHmin, heat index, year, time of year, and number of visitors. Following methods described in Hosmer and Lemeshow [37], we used a forward selection method to build a suite of reasonable negative binomial and Poisson models (S1 File). Our final model selection was based on a minimum Akaike Information Criterion [38] and significantly improved performance via a likelihood ratio test [39].

We applied the regression predictions to climate projections from 14 different global climate models and then fitted a linear model to the projections relating HRI events to year. (The term ‘projection’ refers to a climate model [e.g., HadGEM2-CC365] driven by a specific radiative forcing [e.g., RCP 4.5] Fig 3, S2 File).

Statistical analysis

We used R for all statistical analysis [40]. We used the “stats” package generalized linear model function to construct Poisson models and the negative binomial generalized linear model function in the “MASS” package to construct negative binomial regressions [41]. We made the map using the “leaflet” package and all figures using the “ggplot2” package [42].

Results

For the years 2004–2009, estimated average weekly peak season visitation to Grand Canyon National Park during the months of April 1 to September 30 was 110,596 individuals, with a range of 13,486 to 153,560. Peak visitation occurred in the month of July. A total of 483 HRI events was recorded, including 6 deaths. The average weekly number of HRI events per 100,000 visitors per week during our study period was 2.68 events, with the highest number of HRIs (12 cases) occurring in July 2005 and no HRIs recorded in 23 out of 120 (19%) weeks over the study period.

The best-fitting model was a Poisson regression with the variables of RHmin, Tmax, month (nominal), and month-RHmin interaction (Table 2). We plotted predicted heat events against recorded heat events in Fig 2 to test model performance. Fig 3 shows the projected increase in mean daily heat events for RCP 4.5 and RCP 8.5 for 14 global climate models.

thumbnail
Fig 2. Modeled heat events per 100,000 visitors in Grand Canyon National Park from 1979 through 2100.

A Poisson regression was constructed using observed heat events in GRCA from 2004–2009. Represented here is the daily average per each year of both the modeled (black line) and each year’s daily average observed events (purple points) and its 95% confidence interval (purple error bars). Additionally, the regression model was applied to 14 global climate models and projected heat events are shown in blue (RCP 4.5) and red (RCP 8.5) for each model with the 10th and 90th percentile of the ensemble marked with a dotted line.

https://doi.org/10.1371/journal.pone.0288812.g002

thumbnail
Fig 3. Projected increase in mean weekly heat related illness events per 100,000 visitors to Grand Canyon National Park for representative concentration pathway (RCP) 4.5 (blue) and RCP 8.5 (red) for 14 global climate models from the World Climate Research Programme’s Coupled Model Intercomparison Project phase 5 multi-model dataset.

https://doi.org/10.1371/journal.pone.0288812.g003

thumbnail
Table 2. Variables, model estimates, standard error, and p-values of parameters evaluated for inclusion in the final regression model predicting the rate of heat related illness events in Grand Canyon National Park visitors, May through October, 2004–2009.

https://doi.org/10.1371/journal.pone.0288812.t002

Our model predicted an increase in the weekly average peak season HRI events per 100,000 visitors for all projection models examined, with a range of 0.0035 (RCP 4.5 CanESM2) to 0.032 (RCP 8.5 MROC-ESM-CHEM). This range equates to a total of 0.10 to 0.96 more weekly HRI events each peak season, assuming visitation remains at 2004–2009 levels. This corresponds to an average increase of 29%-137% over 2004–2009 levels, assuming no change in visitation. However, average annual visitation to GCNP is increasing—the estimate for the six-year period on which this study focused (2004–2009) was 4,365,707.5 visitors/year whereas the estimate across six years a decade later (2014–2019) was 5,809,410, an increase of over 33%. If we assume 2014–2019 visitation levels in 2100, 138–254 total peak season HRI events would be expected for RCP 8.5 models, and that assumption may also ultimately prove to be an underestimate because visitor numbers may increase further [43].

Discussion

Heat projection

We found that increased temperatures and decreased relative humidity increased HRI risk. The expected number of HRIs increases by a multiplicative factor of 0.03 for every 1-degree-Fahrenheit increase in maximum temperature with all other factors being equal, and decreases by a multiplicative factor of 0.9 for every 1-unit increase in minimum relative humidity with all other factors being equal. We also found increased risk associated with perceived cooler peak-season months (April and May) when visitors may not be as prepared for HRI. This finding is consistent with findings of [23, 44], who attributed this effect to lack of behavioral adaptation to unexpected heat, which is more likely to occur during the shoulder seasons and is a critical focus for reducing morbidity and mortality from HRI [45]. Finally, we found a significant interaction between month and RHmin, suggesting that the effect of very dry conditions (minimum relative humidity) varies across months, again perhaps due to visitors’ expectations, water packing, and exertion. Prior studies have also found a close association between physical exertion and HRI as well as a change in the climate zone and wearing clothing that is not suitable for the climate, which is common with air travel to national parks, as significant risk factors for HRI [46]. Behavioral adaptation is key to reducing mortality from extreme heat events but requires planning and resources that national park visitors may not be prepared to use without additional awareness and education efforts [47, 48].

HRIs increased significantly under all climate projections. Recent work has found RCP 8.5 models most closely model historical cumulative emissions as well as match most closely with projected emissions based on current policy through the mid-21st century [49]. These results suggest the RCP 8.5 models may more closely predict the actual increase in HRI most likely to be seen at GCNP and predict a concerning increase in disease burden among visitors.

Management implications

The Center for Disease Control and Prevention’s Climate and Health Technical Report identifies health communication as a primary intervention to decrease HRI risk [17]. At GCNP, alerts and warnings are published at NPS.gov/GRCA, on signs throughout the park, and on social media and virtual communication channels. These communications will likely remain a major part of HRI prevention activities for the NPS. Park managers have also made structural changes in recent years, such as building more shade structures. As temperatures increase in coming decades, more such structures may be needed, especially along popular inner-canyon trails (North and South Kebob trails and Bright Angel trails). Another important resource along these trails is potable water. There are several water stations along these inner-canyon trails and a team of rangers called the preventative SAR team provides water and assistance along popular trails, but water availability is inherently limited in the natural, back-country, desert environment found in the Grand Canyon [24]. Similarly, air-conditioned cooling centers are not available in these remote locations nor is there sufficient space in air-conditioned spaces in developed areas of the park to serve all visitors during the peak season months [17]. Occupational heat exposure poses an even greater risk for outdoor workers with physical exertion during heat exposure, which includes many NPS employees such as preventative SAR rangers [50]. As temperatures continue to rise and resources in these back-country locations remain limited, certain previously-used parts of national parks and public lands may not be safe for recreation or stewardship operations during periods of extreme heat. As the climate continues to warm, these periods of extreme heat will become longer, more intense, and more frequent. The NPS is strongly committed to considering climate change adaptation in all its planning processes, but significant changes in important places and loss of access to many important resources may occur as the climate crisis worsens and recreation becomes less safe in these remote settings.

HRI is perhaps the most obvious direct impact of warming temperatures [5, 17]. Some parts of the globe are expected to cease be habitable to humans due to elevated ambient temperatures, so heat is also viewed as a direct driver of climate-induced migration and civil conflict [5, 17, 51]. However, other significant indirect health impacts may also be detrimental to public health [5, 17]. Increased frequency of extreme weather events and wildfires take a dramatic toll on human life and mental health, and declines in agricultural productivity due to climate change are projected to threaten global food security in all regions. Rising temperatures and increased variability in precipitation can dramatically increase the potential for vector-borne and waterborne disease transmission. National parks are highly visited and loved resources where visitors are more likely to be open to new information, behaviors, and transformative experiences [52, 53]. Additionally, exposure to extreme heat or prior HRI events is associated with higher levels of behavioral adaptation and concern about heat risks [54]. Becoming aware of and managing the risk of extreme heat impacts during a recreational visit to a park may therefore have significant potential to stimulate climate change adaptation and mitigation behaviors in visitors elsewhere after their visit, and therefore national parks may have an opportunity for sustained impact on public ethics and behaviors by highlighting climate change in HRI prevention messaging.

Previous studies have found that the greatest burden, most severe disease, and highest overall expenses for treatment of HRI fall disproportionately on minorities and older populations, worsening existing inequities and disproportionately burdening publicly funded systems [55]. Populations without access to air conditioning or with inadequate access to water or health services are more vulnerable to adverse health outcomes from exposure to excessive heat, conditions that can occur from a lack of development or when temperatures rise faster than communities are able to adapt, resulting in a significant disproportionate burden on disadvantaged and minority communities impacted by structural racism [5, 17]. Visitors to the U.S. national parks are disproportionately white, educated, and affluent compared to the general U.S. population, and therefore less likely to be as directly impacted by the earliest impacts of climate change. By highlighting the impacts of climate change to national park visitors, the NPS therefore may engage broader audiences in the conversation about climate change impacts on human health and the need for adaptation and mitigation than through approaches that focus on personal impacts alone.

Limitations

This study is subject to several limitations. It uses data from 2004–2009, after a preventative SAR program was developed by the park and after significant warming had already occurred, which may bias our projections downwards. Other areas where such harm-reduction measures have not been put in place are likely to have a more significant burden of disease than that predicted by our model. The study also likely underestimates HRI events because not all affected individuals sought health care and some victims may have sought care at health care facilities outside of the park, such as at local health care facilities in gateway communities, and/or waited to seek care until they returned home. We also cannot account for adaptation, neither behavioral nor physiological. Although physiological adaptation may be less feasible in non-local visitors, behavioral adaptation is likely to occur, particularly as extreme heat becomes more common. GCNP established a preventative SAR program in 1996, which provides support staff with food and water to prevent adverse health outcomes in visitors, and this program significantly reduced heat-illness related SAR responses and highlights the success of behavioral adaptation in preventing HRI [24].

The range of climate projections used in this study yielded a 10-fold difference in predicted HRI burden. We cannot know today which projection will ultimately turn out to be the most accurate because we cannot know either precisely 1) what future emissions choices humans will make or 2) which GCM among the set we used for our climate projections will most skillfully characterize how climate in this region responds to those future emissions. Therefore, we present the full set of projections, so that managers can understand and address the full range of plausible outcomes [34]. However, as noted above, RCP 8.5-based models characterize historical emissions most closely and are consistent with current emissions-related policy, and thus likely provide the best approximation of future conditions [49] unless humans reduce greenhouse gas emission rates substantially and soon. Our model also cannot address the interaction of heat and air quality. Heat impacts air quality by increasing production of both ground-level ozone and particulate matter (through increased evapotranspiration), and health impacts from poor air quality can mimic and worsen some signs and symptoms of HRI [56]. Additionally, this model is based only on heat events occurring during the current main season at GCNP from April through September and therefore cannot predict any change in events occurring outside of this season. Our model also does not evaluate the impact of rising minimum temperatures, which can play an important role in HRI as physiologic night-time cooling is impaired [57]. We are currently working to expand this work into additional months and contemporary years to better understand HRI risk. Finally, our study focuses exclusively on GCNP visitors and may not be useful for predicting future HRI in other populations whose behaviors may differ significantly from those of the recreating public or for public lands with fewer resources for prevention and response. Specifically, other public lands may see more significant increases in HRI due to the presence of fewer resources for prevention and response, particularly the lack of a preventative SAR program focused on HRI prevention.

Conclusions

Rates of HRI will increase significantly in GCNP visitors in coming decades, with shoulder-season extreme heat increasing heat illness risk more than mid-season extreme heat, as well as extreme heat occurring during low humidity. Because RCP 8.5 climate models most accurately track historical and current warming trends, the highest heat illness risk projections are the most likely and make heat illness risk a significant need for public land management planning and staffing efforts. Like many other aspects of climate change, extreme heat is a threat to human health and wellbeing and poses a significant challenge for public land managers to carry out the mission to provide for safe visitor experiences in our national parks. This threat will worsen with time regardless of how quickly society acts to mitigate continued warming. Land managers and the visiting public alike should be aware of the significant risk climate change poses to human health and wellbeing and work to mitigate these risks in their daily work and lives.

Supporting information

S1 Fig. Predicted heat events vs recorded heat events.

Figure of recorded heat events (black) and predicted heat events (red line) at Grand Canyon National Park from 2004–2009; the pink shaded areas representing 95% prediction confidence intervals from the final negative binomial regression model.

https://doi.org/10.1371/journal.pone.0288812.s001

(TIF)

S1 File. Binned HRI data.

Deidentified binned data of visitation estimates and recorded HRIs used in regressions.

https://doi.org/10.1371/journal.pone.0288812.s002

(CSV)

S2 File. Poisson model AICs table.

Table of AIC of different Poisson models tested.

https://doi.org/10.1371/journal.pone.0288812.s003

(CSV)

S3 File. Predicted HRI increase table.

Table of predicted mean increase in cases per year for each climate model, estimated using lm(Events ~ year). Events increased per year taken as coefficient of slope [95% CI].

https://doi.org/10.1371/journal.pone.0288812.s004

(CSV)

S4 File. R code of regression model fitting.

Rmd file with code and output for regression model fitting.

https://doi.org/10.1371/journal.pone.0288812.s005

(HTML)

S5 File. R code of regression model prediction.

Rmd file with code and output for regression model prediction.

https://doi.org/10.1371/journal.pone.0288812.s006

(HTML)

Acknowledgments

We would like to thank Pam Ziesler (NPS) for assistance with GCNP visitation data; Amber Runyon (NPS) for assistance with climate model data and manuscript review; Joel Reynolds (NPS) for careful review and comments on statistical methods and interpretation; Nick Fudala (NPS) for assistance with interpreting daily visitation data for GCNP; David Lawrence (NPS) for review of data visualization and manuscript; Jennifer Proctor (NPS) for access to the historical Noe et al. 2013 dataset. We also extend our gratitude to the GCNP Preventative and regular Search and Rescue Rangers who protect visitor health and safety.

References

  1. 1. Atwoli L, Baqui AH, Benfield T, Bosurgi R, Godlee F, Hancocks A, et al. Call for Emergency Action to Limit Global Temperature Increases, Restore Biodiversity, and Protect Health. N Engl J Med. 2021, 385:1134–1137. pmid:34491006
  2. 2. Pachauri RK, Mayer L, Intergovernmental Panel on Climate Change (Eds.). Climate change 2014: Synthesis report. Intergovernmental Panel on Climate Change: 2015.
  3. 3. WHO. WHO calls for urgent action to protect health from climate change–Sign the call. WHO; World Health Organization. 2020. http://www.who.int/globalchange/global-campaign/cop21/en/
  4. 4. Greene S, Kalkstein LS, Mills DM, Samenow J. An Examination of Climate Change on Extreme Heat Events and Climate–Mortality Relationships in Large U.S. Cities. Weather Clim Soc. 2011, 3(4): 281–292.
  5. 5. Smith KR, Woodward A, Campbell-Lendrum D, Chadee DD, Honda Y, Liu Q, et al. Human health: Impacts, adaptation, and co-benefits. In Climate Change 2014: Impacts, Adaptation, and Vulnerability: 2014.
  6. 6. Schwartz JD, Lee M, Kinney PL, Yang S, Mills D, Sarofim MC, et al. Projections of temperature-attributable premature deaths in 209 U.S. cities using a cluster-based Poisson approach. Environ Hlth. 2015, 14(1):85. Available from: pmid:26537962
  7. 7. Stonehill RB. Heat illnesses. J Am Coll Emerg Physicians. 1972, 1(6):21–23.
  8. 8. Schulte F, Röösli M, Ragettli MS. Heat-related cardiovascular morbidity and mortality in Switzerland: a clinical perspective. Swiss Med Wkly. 2021,151:w30013. pmid:34519460
  9. 9. Park RJ, Goodman J, Behrer AP. Learning is inhibited by heat exposure, both internationally and within the United States. Nat Hum Behav. 2020, 1–9. Available from: pmid:33020588
  10. 10. Thompson R, Hornigold R, Page L, Waite T. Associations between high ambient temperatures and heat waves with mental health outcomes: A systematic review. Public Health. 2018, 161:171–191. Available from: pmid:30007545
  11. 11. Tiihonen J, Halonen P, Tiihonen L, Kautiainen H, Storvik M, Callaway J. The Association of Ambient Temperature and Violent Crime. Sci Reps. 2017, 7(1):6543. Available from: pmid:28754972
  12. 12. Wang JC, Chien WC, Chu P, Chung CH, Lin CY, Tsai SH. The association between heat stroke and subsequent cardiovascular diseases. PloS One. 2019, 14(2):e0211386. Available from: pmid:30759128
  13. 13. Glaser J, Lemery J, Rajagopalan B, Diaz HF, García-Trabanino R, Taduri G, et al. Climate Change and the Emergent Epidemic of CKD from Heat Stress in Rural Communities: The Case for Heat Stress Nephropathy. Clin J Am Soc Nephrol. 2016, 11(8):1472–1483. pmid:27151892
  14. 14. Nerbass FB, Pecoits-Filho R, Clark WF, Sontrop JM, McIntyre CW, Moist L. Occupational Heat Stress and Kidney Health: From Farms to Factories. Kidney Int Rep. 2017, 2(6):998–1008. pmid:29270511
  15. 15. Tseng MF, Chou CL, Chung CH, Chen YK, Chien WC, Feng CH, et al. Risk of chronic kidney disease in patients with heat injury: A nationwide longitudinal cohort study in Taiwan. PLoS One. 2020, 15(7): e0235607. pmid:32614909
  16. 16. Wallace RF, Kriebel D, Punnett L, Wegman DH, Amoroso PJ. Prior heat illness hospitalization and risk of early death. Environ Res. 2007, 104(2):290–295. Available from: pmid:17306249
  17. 17. Abbinett, J, Schramm PJ, Widerynski S, Saha S, Beavers S, Eaglin M, et al. Climate and Health Technical Report Series. 2020. 60.
  18. 18. Bekkar B, Pacheco S, Basu R, DeNicola N. Association of Air Pollution and Heat Exposure With Preterm Birth, Low Birth Weight, and Stillbirth in the US. JAMA Network Open. 2020, 3(6). pmid:32556259
  19. 19. Buttke D, Allen D, Higgins C. Benefits of Biodiversity to Human Health and Well-being (U.S. National Park Service). 2013. October 27, 2020, https://www.nps.gov/articles/parksciencev31-n1_buttke_etal-htm.htm.
  20. 20. NPS. Park Statistics—Grand Canyon National Park (U.S. National Park Service). 2018. https://www.nps.gov/grca/learn/management/statistics.htm.
  21. 21. Monahan WB, Fisichelli NA. Climate Exposure of US National Parks in a New Era of Change. PLOS ONE. 2014, 9(7):e101302. pmid:24988483
  22. 22. Yosemite National Park. A wave of heat-related illnesses. 2002. Accessed 2 May 2022. https://www.nps.gov/yose/blogs/a-wave-of-heat-related-incidents.htm?msclkid=4596d24bca3911ecb13292d670c33b3b.
  23. 23. Noe RS, Choudhary E, Cheng-Dobson J, Wolkin AF, Newman SB. Exertional heat-related illnesses at the Grand Canyon National Park, 2004–2009. Wilderness & Environmental Medicine. 2013, 24(4):422–428. pmid:24119571
  24. 24. Malcolm CJ, Heinrich H. Preventive success! Grand Canyon’s response to search-and-rescue overload. Park Science. 2017, 33(1):99–107.
  25. 25. Pearce EA, Jelínková L, Fullerton L, Malcolm CJ, Heinrich HL, Norwil EJ, et al. Observational Study of Grand Canyon Rim-to-Rim Day Hikers: Determining Behavior Patterns to Aid in Preventive Search and Rescue Efforts. Wilderness Environ Med. 2019 Mar;30(1):4–11. pmid:30391110
  26. 26. NPS. Administrative Boundaries of National Park System Units—National Geospatial Data Asset (NGDA) NPS National Parks Dataset. NPS—Land Resources Division. Accessed 2 February 2023. Accessible: https://irma.nps.gov/DataStore/Reference/Profile/2224545.
  27. 27. ESRI. Maxar, Earthstar Geographics, and the GIS User Community. World_Imagery (Map Server). Accessed 2 February 2023. Accessible: World_Imagery (MapServer) https://services.arcgisonline.com/ArcGIS/rest/services/World_Imagery/MapServerarcgisonline.com.
  28. 28. NPS. Grand Canyon National Park Total Recreational Visitors. Accessed 4 May 2022. Accessible: Stats Report Viewer (nps.gov).
  29. 29. Abatzoglou JT. Development of gridded surface meteorological data for ecological applications and modelling. Int J Climatol. 2013. 33(1):121–131.
  30. 30. NOAA. Heat Index Equation. 2014. https://www.wpc.ncep.noaa.gov/html/heatindex_equation.shtml.
  31. 31. Abatzoglou JT, Brown TJ. A comparison of statistical downscaling methods suited for wildfire applications. Int J Climatol. 2012. 32(5):772–780.
  32. 32. Taylor KE, Stouffer RJ, Meehl GA. An overview of CMIP5 and the experimental design. Bulletin of the American Meteorological Society. 2012, 93(4):485–498.
  33. 33. Runyon AN, Carlson AR, Gross J, Lawrence DJ, Schuurman GW. Repeatable approaches to work with scientific uncertainty and advance climate change adaptation in US national parks. Parks Stewardship Forum. 2020, 36(1). Available from:
  34. 34. Lawrence DJ, Runyon AN, Gross JE, Schuurman GW, Miller BW. Divergent, plausible, and relevant climate futures for near- and long-term resource planning. Clim Change. 2021, 167:1–20.
  35. 35. Moss RH, Edmonds JA, Hibbard KA, Manning MR, Rose SK, van Vuuren DP, et al. The next generation of scenarios for climate change research and assessment. Nature. 2010, 463(7282):747–756. pmid:20148028
  36. 36. van Buuren S, Groothuis-Oudshoorn K. mice: Multivariate Imputation by Chained Equations in R. J Stat Softw. 2011, 45(1):1–67.
  37. 37. Hosmer DW, Lemeshow S. Applied Logistic Regression:1992. https://books.google.com/books/about/Applied_Logistic_Regression.html?id=Po0RLQ7USIMC.
  38. 38. Akaike H. A new look at the statistical model identification. IEEE Trans Automat Contr. 1974, 19(6):716–723.
  39. 39. Buse A. The Likelihood Ratio, Wald, and Lagrange Multiplier Tests: An Expository Note. Am Stat. 1982, 36(3):153.
  40. 40. The R core team. R: The R Project for Statistical Computing. 2019. https://www.r-project.org/
  41. 41. Venables WN. Modern applied statistics with S /: 2002.
  42. 42. Wickham H, Navarro D, Lin T. ggplot2: Elegant Graphics for Data Analysis. 2016. https://ggplot2-book.org/.
  43. 43. Fisichelli NA, Schuurman GW, Monahan WB, Ziesler PS. Protected Area Tourism in a Changing Climate: Will Visitation at US National Parks Warm Up or Overheat? PLOS ONE. 2015, 10(6): e0128226. pmid:26083361
  44. 44. Liss A, Naumova EN. Heatwaves and hospitalizations due to hyperthermia in defined climate regions in the conterminous USA. Environ Monit Assess. 2019, 191(Suppl 2):394. pmid:31254102
  45. 45. Ebi KL, Capon A, Berry P, Broderick C, de Dear R, Havenith G, et al. Hot weather and heat extremes: health risks. Lancet. 2021, 398(10301):698–708. pmid:34419205
  46. 46. Leyk D, Hoitz J, Becker C, Glitz KJ, Nestler K, Piekarski C. Health Risks and Interventions in Exertional Heat Stress. Dtsch Arztebl Int. 2019,116(31–32):537–544. pmid:31554541
  47. 47. Arbuthnott K, Hajat S, Heaviside C, Vardoulakis S. Changes in population susceptibility to heat and cold over time: assessing adaptation to climate change. Environ Health. 2016, 15 Suppl 1(Suppl 1):33. pmid:26961541
  48. 48. Howe PD, Marlon JR, Wang X, Leiserowitz A. Public perceptions of the health risks of extreme heat across US states, counties, and neighborhoods. Proc Natl Acad Sci U S A. 2019, 116(14):6743–6748. pmid:30862729
  49. 49. Schwalm CR, Glendon S, Duffy PB. RCP8.5 tracks cumulative CO2 emissions. PNAS. 2020, 117(33):19656–19657. Available from: pmid:32747549
  50. 50. Chapman CL, Hess HW, Lucas RAI, Glaser J, Saran R, Bragg-Gresham J, et al. Occupational heat exposure and the risk of chronic kidney disease of nontraditional origin in the United States. Am J Physiol Regul Integr Comp Physiol. 2021, 321(2):R141–R151. pmid:34161738
  51. 51. Watts N, Amann M, Arnell N, Ayeb-Karlsson S, Belesova K, Berry H, et al. The 2018 report of the Lancet Countdown on health and climate change: Shaping the health of nations for centuries to come. Lancet. 2018, 392(10163):2479–2514. Available from: pmid:30503045
  52. 52. Robledo MA, Batle J. Transformational tourism as a hero’s journey. Current Issues in Tourism. 2017, 20:1736–1748.
  53. 53. Morse WC, Stern M, Blahna D, Stein T. Recreation as a transformative experience: Synthesizing the literature on outdoor recreation and recreation ecosystem services into a systems framework, Journal of Outdoor Recreation and Tourism. 2022, 38:100492,
  54. 54. Liu T, Xu YJ, Zhang YH, Yan QH, Song XL, Xie HY,et al. Associations between risk perception, spontaneous adaptation behavior to heat waves and heatstroke in Guangdong province, China. BMC Public Health. 2013, 13:913. pmid:24088302
  55. 55. Schmeltz MT, Petkova EP, Gamble JL. Economic Burden of Hospitalizations for Heat-Related Illnesses in the United States, 2001–2010. International Journal of Environmental Research and Public Health. 2016, 13(9). Available from: pmid:27618079
  56. 56. Abed Al Ahad M, Sullivan F, Demšar U, Melhem M, Kulu H. The effect of air-pollution and weather exposure on mortality and hospital admission and implications for further research: A systematic scoping review. PLoS ONE. 2020, 15(10):e0241415. pmid:33119678
  57. 57. He C, Kim H, Hashizume M, Lee W, Honda Y, Kim SE, et al. The effects of night-time warming on mortality burden under future climate change scenarios: a modelling study. Lancet Planet Health. 2021, 6(8):e648–e657. pmid:35932785