Introduction

After spreading around the world in a matter of months, the coronavirus (COVID-19) has become a leading cause of death in the United States. According to recent reports, COVID-19 “deaths have snowballed from a few isolated cases to thousands across the country each day” (Keating and Esteban 2020). Although the Centers for Disease Control and Prevention (CDC) (2020) have proposed several potential mitigation strategies, social distancing and staying at home have received the most national attention. CDC officials and frontline health care professionals advise that the best way to prevent exposure to the virus is to stay at home and to avoid close contact with people. As the official website of the Hopi Tribe (2020) explains, “the virus does not move, people move it…if people stop moving, the virus stops moving and dies.” If staying at home is so important for slowing the spread of the coronavirus, we must begin to consider the social patterning of mobility (i.e., how far people are traveling from home in their daily lives). The fundamental question is whether certain populations travel more or less than others during the pandemic.

In this paper, we consider the association between state religiosity and state mobility scores. Over the past several weeks, the institution of religion has received a great deal of media attention and public scrutiny. On March 18, Newsweek reported that “Pastor Holds Service with over 1000 Parishioners in Defiance of Large-Gathering Ban” (Slisco 2020). On March 19, The Washington Post advised that “Coronavirus Deniers and Hoaxers Persist despite Dire Warnings, Claiming ‘it’s Mass Hysteria’” (Gowen 2020). On March 27, The New York Times suggested that “The Religious Right’s Hostility to Science is Crippling our Coronavirus Response” (Stewart 2020). Although not representative of the religious spectrum in this country, all of these stories (and many more like them) seem to point to the same general conclusion: Religious populations and communities may be especially likely to acquire and spread the coronavirus.

In the pages that follow, we explore relevant research from the sociology of religion and recent religious rhetoric surrounding the coronavirus pandemic. We then model state mobility scores as a function of state religiosity and covariates. After summarizing our key results, we discuss the contributions and limitations of our study. We end with important directions for future research on the social patterning of mobility during the coronavirus pandemic.

Religion, Science, and the Pandemic

Our central argument is that more religious populations may be especially resistant to public health recommendations during the coronavirus pandemic (e.g., social distancing and staying at home) because they hold more negative views of science and scientists and strong religious beliefs concerning the pandemic itself. Several studies show that more religious populations tend to report less trust in science as a social institution and more anti-science attitudes (Evans 2013; Gauchat 2008, 2012). Of course, these positions are not representative of all religious groups. There is at least some evidence to suggest that more conservative-leaning Protestant denominations may be less literate in science and especially critical of the scientific community and the potential benefits of scientific progress (Ellison and Musick 1995; Evans 2013; Gauchat 2008; Sherkat 2011, 2017). For example, studies show that conservative Protestants are often less concerned with environmental degradation, less trusting of the findings of climate scientists, and more likely to endorse “a polluting creed” (Sherkat and Ellison 2007; Smiley 2019).

Many conservative Protestant denominations see the Bible as the ultimate source of authority and direction when interpreting and experiencing personal life and world events (Boone 1989; Ellison et al. 1996). In contrast to the positivist logic implied by the scientific method, so-called biblical literalists assess the legitimacy of scientific information by its apparent compatibility with scripture (Ellison and Musick 1995). Religious conservatives, guided by pastors and other religious elites, often draw on religious scripture to oppose scientific recommendations that are perceived as immoral or defined as encroaching on religious liberty or the will or grace of God. Moreover, tensions between religion and science are often rooted in fears concerning the profane influence of science on society (Evans 2013) and a “social conflict between institutions struggling for power” (Evans and Evans 2008:97).

Our perspective is that these religious belief systems are likely to serve as an ideological basis for resisting public health recommendations and initiatives during the pandemic. The specific flashpoints appear to center on (1) denying health information from health scientists (mistrust of science), (2) accepting health misinformation from religious and political leaders (religious authority), and (3) rejecting restrictions on in-person religious services (religious liberty). These themes are clearly represented in several recent cases of religious leadership.

Pastor Andrew of the USA Christian Church in California tells us that “our safety is at stake since national disobedience of God’s laws brings danger and diseases, such as coronavirus, but obeying God brings covenant protection. God protects the USA from danger as the country repents of LGBT, false gods, abortion and other sins” (Rosen 2020).

Reverend Curtis of Havre Assembly of God Church in Montana refers to COVID-19 as an acronym for “Christ Over Viruses and Infectious Diseases” (Curtis 2020). He also directs us to “Joshua 1:9—Be strong and courageous. Do not be frightened, and do not be dismayed, for the Lord your God is with you wherever you go.” Based on this passage, he asks, “Shall we deal with fright and dismay with the strength and courage that God prescribes? Shall we face an uncertain future, knowing that God is with us no matter what?”

Reverend Spell of the Life Tabernacle Church in Louisiana is convinced that the virus is “not a concern” because it is “politically motivated” (Slisco 2020). Reverend Spell also expresses faith in the healing powers of his church: “Our church is a hospital where the sick can come and get healing. Cancers are healed here, people are healed of HIV in these services, and we believe that tonight, we’re also going to pass out anointed handkerchiefs to people who may have a fear, who may have a sickness, and we believe that when those anointed handkerchiefs go, that healing virtue is going to go on them as well.”

Pastor Howard-Browne of The River at Tampa Bay Church in Florida describes people who are concerned about the coronavirus as “pansies” and insists he would only close his church “when the rapture is taking place” (Stewart 2020). Days after making these comments, the pastor was arrested by the sheriff of Hillsborough County because “his reckless disregard for human life put hundreds of people in his congregation at risk and thousands of residents who may interact with them this week in danger” (Mazzei 2020).

Bishop Glenn of New Deliverance Evangelistic Church in Virginia told his congregation that “God is larger than this dreaded virus” and that “people are healed” in his church (Brown 2020). A few weeks later, the pastor died after being diagnosed with the coronavirus.

Even President Trump has attempted to connect with the rhetoric of religious leadership to push his precarious economic timelines, suggesting that “pews should be filled on Easter,” knowing that “his base will revel in the symbolism of resurrection” (Wise 2020). The general concern of health professionals is that churches may become “coronavirus hotspots” because such “meetings are counter to the prescribed public health policy of groups of people coming together” (Gattis 2020). Scholars in public health echo these sentiments, noting “the way people interact in churches, synagogues, mosques and other religious facilities—shaking hands, hugging, singing—appears conducive to what epidemiologists call ‘super-spreading events’” (Collier et al. 2020). Emboldened by religious leaders and the President, some churchgoers in Missouri have defied social distancing warnings with claims that they are protected from the coronavirus because they are “covered in Jesus’ blood” (Edwards 2020).

Hypotheses

In accordance with these arguments, we developed three hypotheses to guide subsequent analyses. Our first hypothesis (H1) is that more religious states will exhibit higher average mobility scores. In other words, more religious states will tend to travel more than less religious states during the pandemic. This hypothesis will be tested through the direct effect of state religiosity on state mobility scores. Our second hypothesis (H2) is that state mobility scores will decline at slower rates in more religious states. The idea is that religious states will be especially reluctant to change their behavior over the study period and that reductions in travel during the pandemic will be, on average, smaller in these states. This hypothesis will be tested through the moderation of state mobility trends (i.e., weekly variations in mobility scores) by state religiosity (i.e., week*religiosity). Our final hypothesis (H3) is that the inverse association between number of days with a state stay-at-home order and mobility scores will be less pronounced in more religious states. This hypothesis will be tested through the interaction of days with a state stay-at-home order and state religiosity (i.e., stay-at-home order*religiosity).

Data

We employ 8 weeks (February 24–April 13) of first-party mobility data from Cuebiq (COVID-19 Mobility Insights 2020), religion data from the 2010 U.S. Religion Census (Grammich et al. 2018) and the 2014 Religious Landscape Study (Pew Research Center 2015), demographic characteristics from the 2018 American Community Survey: 5-Year Estimates (American Community Survey 2018) and the 2019 Bureau of Labor Statistics (Bureau of Labor Statistics 2019), stay-at-home orders from the New York Times (Mervosh et al. 2020), and governor’s known political party. Data are limited to the contiguous United States because mobility estimates exclude Alaska and Hawaii. The District of Columbia is also omitted due to missing data on governor’s political party. All measures are state-level. Our final analytic sample size is 348 (48 states*8 weeks).

Measures

Mobility

Our outcome variable is mobility. Cuebiq, an offline intelligence and measurement marketing company, partners with 86 apps to collect first-party location data to understand population behavior during the coronavirus pandemic via Software Development Kit (SDK) technology. SDK technology improves accuracy and precision in location data collection by linking with mobile phone operating systems, global positioning systems, and Wi-Fi signals. Anonymous and privacy compliant geo-behavioral data are collected for opted-in users, including movement and stops to determine dwell time and visit frequency at locations. Each day, always-on data collection accumulates, on average, 100 data points for approximately 15,000,000 cell phone users. The mobility index measures the average median distance traveled by all devices for the weeks of February 24–April 13 for each state. This time frame was selected because it marks the first period of national coronavirus awareness and changing mobility across the United States. In preliminary analyses, we assessed the construct validity of the mobility index by testing whether these scores were associated with the number of days with a state stay-at-home order. Our analyses revealed a moderate inverse association (r = − 0.45, p < 0.001). In other words, states with longer periods of stay-at-home orders tended to exhibit lower mobility scores.

Religiosity

Our focal predictor variable is religiosity. We measure religiosity as a mean index of six variables (α = 0.98), including (1) percent evangelical, (2) percent who report attending worship services at least weekly, (3) percent who identify as highly religious, (4) percent who say religion is very important in their lives, (5) percent who say they pray daily, and (6) the percent who say they believe in God with absolute certainty. This index assesses how religious a state is by assessing self-reported religious identities, beliefs, and practices. Percent evangelical was collected through the 2010 U.S. Religion Census: Religious Congregations and Membership Study. The remaining variables were collected through the 2014 Religious Landscape Study.

Background Variables

Our analyses include a range of state-level background variables that are at least theoretically related to mobility and religiosity, including (1) median age, (2) percent black, (3) the unemployment rate, (4) population density, (5) mobility lag (the mean of mobility scores for the first 2 weeks of the study period), (6) governor’s political party, and (7) number of days with a state stay-at-home order. Median age, percent black, and population density are 5-year estimates from the 2018 American Community Survey. The unemployment rate for December 2019 was obtained from the Bureau of Labor Statistics. Governor’s party identification is publicly available. Data on state stay-at-home orders as of April 6, 2020 were collected by the New York Times.

Analysis

Our analytic strategy proceeds in three steps. In Table 1, we present descriptive statistics for all study variables, including variable ranges, means, and standard deviations. In Table 2, we fit a model of state mobility scores using robust regression with the M estimator and iterated re-weighted least squares (IWLS) to down-weight the influence of outliers (Jorgenson 2007; Venables and Ripley 2002). Model 1 regresses state mobility scores on dummy variables for week. This model tests whether mobility scores vary across the 8-week study period. We use dummy variables for week to capture thresholds in the association between week and mobility. Model 2 adjusts for state religiosity, median age, percent black, the unemployment rate, population density, and the mobility lag. This model tests whether state religiosity is associated with state mobility scores net of state demographic characteristics and baseline mobility scores. The question here is whether more religious states tend to exhibit higher average mobility scores. Population density and the mobility lag capture factors related to state mobility, including the concentration of populations and the average median distance traveled for the first 2 weeks of the study period. Model 3 adjusts for governor’s party identification and number of days with a state stay-at-home order to assess the association between religiosity and mobility net of state political orientation and virus policy environments.

Table 1 Descriptive statistics for selected study variables
Table 2 Robust regression of mobility scores (direct effects)

Table 3 and Figs. 1, 2 and 3 present our moderation models. Model 1 of Table 3 adds week by religiosity interaction terms to the regression equation to formally test whether state mobility trends (weekly variations in mobility scores) vary by state religiosity. The question here is whether state mobility scores declined at a slower rate in more religious states. Model 2 adds the interaction term for number of days with a state stay-at-home order by religiosity to test whether the association between state stay-at-home order and mobility varies by state religiosity. The question here is whether the effects of stay-at-home orders are weaker for more religious states. Figure 1 provides a graphical illustration of Model 1 (Table 3) state mobility trends by state religiosity. Figure 2 presents overall mobility reductions (changes between the first and last week of the study period) for three exemplar states, one of low, moderate, and high religiosity. Finally, Fig. 3 provides a graphical illustration of the stay-at-home order*religiosity interaction terms in Model 2 (Table 3).

Table 3 Robust regression of mobility scores (moderation effects)
Fig. 1
figure 1

Adjusted state mobility trends by state religiosity

Fig. 2
figure 2

Percent change in mobility (2/24–4/13) by religious exemplar states

Fig. 3
figure 3

Adjusted state mobility by days under stay-at-home order and state religiosity

Results

Descriptive Statistics

In Table 1, we see an overall trend toward declining state mobility scores over the study period. State mobility was highest during the weeks of February 24 and March 2. From March 2 to March 23, we observed the first trend toward declining state mobility. After a short period of slightly increasing state mobility from March 23 to April 6, we recorded our lowest mobility scores yet for the week of April 13.

Direct Effects

In Table 2, Model 1 compares mobility scores for the first week (February 24) of our study period to each of the 7 weeks from March 2 to April 13. We observe that mobility scores for the week of February 24 are comparable to those for the weeks of March 2 and March 9. We then find lower mobility scores for the weeks of March 16 to April 13 than for the week of February 24.

Model 2 reveals higher state mobility scores for states that are more religious, older, and more mobile at the beginning of the study period (February 24–March 2). We also see lower state mobility scores for states with generally higher unemployment rates and greater population density. Mobility scores do not appear to vary by percent black in this model.

Model 3 shows that state mobility scores are lower in states with longer periods of state stay-at-home orders. State mobility scores do not appear to vary by Governor’s party identification in this model. In short, state mobility scores are comparable for states with Republican and Democratic Governors.

Moderation Effects

In Table 3, the interaction terms in Model 1 suggest that, compared to the week of February 24, state mobility scores declined at slower rates in more religious states for five of the 8 weeks under study (March 16–April 13). Model 2 of Table 3 indicates that the inverse association between number of days with a state stay-at-home order and mobility is less pronounced in more religious states.

Figure 1 represents how mobility for three levels of religiosity changed over the 8-week study period, controlling for all covariates. From March 16 to April 13, states with high religiosity exhibited the highest state mobility scores. Figure 2 presents the percent reduction in overall state mobility scores between the weeks of February 24 and April 13 for three exemplar states. From the first week to the last week of the study period, state mobility scores declined by 22% for Alabama (high religiosity), by 32% for Indiana (moderate religiosity), and by 64% for Massachusetts (low religiosity). Finally, Fig. 3 clearly demonstrates that state stay-at-home orders had a weaker impact on mobility in more religious states.

Supplemental Analyses

In supplemental analyses, we estimated our regression models using ordinary least squares regression with cluster robust standard errors for state (Robitzsch and Grund 2020) to assess potential violations of regression assumptions (serial correlation and heteroscedasticity). The results were substantively identical to our robust regressions models. In other sensitivity analyses, we controlled for total population size, percent Hispanic, percent with a college degree, per capital income, income inequality (Gini coefficient), and percent obese. In each case, the additional variable was either unrelated to mobility or failed to alter any of our substantive conclusions.

Discussion

In this paper, we considered the association between state religiosity and state mobility scores during the coronavirus (COVID-19) pandemic. Our first hypothesis (H1) specified that more religious states would tend to exhibit higher average mobility scores. In support of this hypothesis, our analyses showed that more religious states tended to travel more than less religious states during the pandemic. This general pattern persisted with adjustments for median age, percent black, the unemployment rate, population density, a mobility lag, governor’s political party identification, and number of days with a state stay-at-home order. Our second hypothesis (H2) indicated that state mobility scores would decline at slower rates in more religious states. Consistent with this hypothesis, our analyses revealed that reductions in travel during the pandemic were, on average, smaller in more religious states. Our final hypothesis (H3) stated that that the inverse association between number of days with a state stay-at-home order and mobility scores would be less pronounced in more religious states. In accordance with this hypothesis, we observed that religious states were more resistant to stay-at-home orders.

To our knowledge, this is the first empirical study of religious contextual variations in coronavirus pandemic behavior. Our findings are important because they seem to confirm the suspicion that religious populations and communities may be especially likely to acquire and spread the coronavirus. In the context of a pandemic, the risk may be even more pervasive. Religious states include people who are religious and people who are not. Because the coronavirus is so highly contagious, any ideological resistance to public health recommendations could present an existential threat to society. We need to begin to think about ways of overcoming religious cultural barriers to critical pandemic responses. Potential strategies or interventions must systematically address obstacles related to the mistrust of science, religious authority, and religious liberty.

We acknowledge that our analyses are limited in three key respects. Although our data suggest that more religious states tend to be more mobile during the coronavirus pandemic, we cannot conclude that religious individuals are more mobile without individual-level data. Because we are currently in the early stages of the pandemic, we are unable to assess concurrent changes in our state-level predictors across years. Finally, our mobility data are also limited in the sense that they are based on geo-behavioral data for opted-in users, not probability samples in states. We are nevertheless encouraged by the fact that our mobility scores are predictably associated with a range of variables, including median age, the unemployment rate, population density, a mobility lag, and number of days with a state stay-at-home order.

Despite these limitations, we provided the first empirical study of state religiosity and mobility during the coronavirus pandemic. Our analyses consistently showed that, in the early weeks of the pandemic, religious states tended to travel more and were more resistant to changing their movement patterns. Our analyses are important because they contribute to our understanding of the social patterning of pandemic mobility, which is ultimately relevant to slowing the spread of the coronavirus. More research is needed to replicate our findings using longer longitudinal designs and data collected at different levels of analyses, including the county and individual levels. As more valid and reliable epidemiological data become available, we will need to assess whether state infection and mortality rates also vary according to state religiosity. Future work should continue to consider the social patterning of pandemic mobility more broadly. Understanding the role of state political environments is perhaps most pressing. Research along these lines would provide a more thorough understanding of the impact of social and ideological forces on pandemic behavior.