Chapter 3 Participation

Does ethnic density influence community participation in mass participation physical activity events?: a case of parkrun

Smith R1, Schneider PP1, Bullas A2, Haake SSJ2, Goyder E1


1School of Health and Related Research, University of Sheffield
2Sheffield Hallam University Advanced Wellbeing Centre


Keywords: parkrun, Physical Activity, Deprivation, Ethnic density.

Intended Journal: Research in Sport, Exercise and Health, Sports Medicine Open,Leisure Sciences, BMC Public Health.


3.1 Introduction

parkrun is a collection of free mass participation 5km running events that takes place every Saturday morning. There are currently over 500 locations in England, with a combined weekly attendance of over 100,000. parkrun has been identified as being successful at engaging with individuals who may not otherwise have taken part in organised physical activity (S. Haake 2018b; Stevinson and Hickson 2013b), and there is some evidence that it has increased overall physical activity levels in participants (Stevinson and Hickson 2018). Overall, there is a consensus that parkrun has huge public health potential (Reece et al. 2019b).

However, qualitative research from Sheffield (Goyder et al. 2018) and more broadly the United Kingdom (Fullagar et al. 2019) identified that parkruns located in more deprived areas have lower attendances, and that ethnic diversity in parkrun was limited. This leads to concern that as with many public health interventions, parkrun is “likely to be responsible for significant intervention generated inequalities in uptake of opportunities for physically active recreation” (Goyder et al. 2018).

Undertaking quantitative analysis of the determinants of participation in parkrun is therefore long overdue. Aside from a single previous study from Australia (Cleland et al. 2019), with substantial limitations including, as noted by the authors, that “The sample was limited to a non-random sample of parkrun participants in one State of Australia and may not be generalizable to other parkrun populations.” (p.21), no other studies have attempted to identify the determinants of participation in parkrun.

Our previous work revealed that there is substantial heterogeneity in parkrun participation across different communities in England: after controlling for geographical distance to nearest event, deprived communities have significantly lower participation rates (P. P. Schneider et al. 2019b). The analysis was able to quantify, for the first time, how participation in parkrun varied in different communities in England. However, the analysis was interested only in the relationship between participation, access and deprivation and did not consider ethnic density as a potential determinant of participation in parkrun. Yet, evidence from survey data shows that non-White-British individuals in England are less likely to be physically active, and to engage in sport in general (Rowe and Champion 2000). We thus hypothesised that at the community level, all else being equal, areas with higher ethnic density have lower levels of participation in parkrun participation.

3.2 Methods

We undertook an ecological analysis of parkrun participation in England in 2018. Data was obtained from multiple sources at the Lower layer Super Output Level (LSOA). There are 32,844 LSOAs in England, each of which is a geographical area containing around 1,000-3,000 people.

parkrunUK provided data on the number of parkrun finishers from each LSOA in England between the 1st January and 10th December 2018, we use the number of finishers as a proxy for parkrun participation, althoough we appreciate that people participate in parkrun in other ways (e.g. volunteering). We also used parkrun event location data, which are publicly available on the parkrunUK website.

The rest of the data, including Index of Multiple Deprivation (IMD) Score, Ethnic Density, Rural-Urban Classification, Population Density, Percentage Working Age and LSOA centroids were obtained from the Office of National Statistics (ONS). Full sources are listed in the table below, and all ONS data is provided open source on the author’s GitHub page.

Table 1. Variables used in the Analysis

Variable Description Source
run_count number of finishers from each LSOA in England between 1st January and 10th December 2018 parkrunUK
imd Index of Multiple Deprivation scores for each LSOA ONS
total_pop total number of individuals in each LSOA ONS.
pop_density population density for each LSOA ONS
rural_urban Rural-Urban Classification ONS
perc_bme Ethnic Density: percentage of population non-white-british ONS
mn_dstn distance from LSOA centroid to nearest parkrun derived from ONS
perc_non_working_age derived from ONS data on age-groups in each LSOA ONS
run_rate derived from run_count and LSAO populations derived

After merging these datasets we had detailed data on 32,844 LSOAs, including participation (number of finishers) and several characteristics of the LSOAs which we hypothesised may influence participation. Since previous work (P. P. Schneider et al. 2019b) has found correlations between participation and deprivation, distance to nearest event, and population density we included all these variables. We also extended the analysis to include ethnic density (we use the percentage of the population that reported being non-White-British as a proxy for ethnic density) and the percentage of the population of working age. We are interested in ethnic density as we hypothesised that areas with higher ethnic density would have lower participation rates, all else being equal. We included the percentage of the population that is working age as a control to limit for the effect of populations heavily skewed toward older people (e.g. care homes), or young people (e.g near schools). Since participation in parkrunUK is dominated by those aged 20-60 S. Haake (2018b) we felt this was justified.

We first studied the bivariate Pearson correlations between the dependant variable (weekly parkrun participation per 1000 population) and all indepdentant variables described. Results are visually illustrated using the correlation plots (Wei and Simko 2017) and stratified (by rural/urban classification) heat maps. We then fitted a multivariable regression model to investigate the independent effects of the predictors on parkrun participation. A Poisson distribution with a log-link and with the LSOA total population as an offset variable was used to model parkrun participation as a rate (runners per LSOA population). Model fit was assessed using Pseudo R^2, based on quasi-likelihood functions (Zhang 2017). All analyses were conducted in R (R Core Team 2018).

3.3 Results

3.3.1 Descriptive Statistics

Participation in parkrun varies across LSOAs. Around half of all communities (LSOA) average less than 1 finisher per week per 1000 people. Approximately a quarter average between 1 and 2 finishers, and around an eighth between 2 and 3 finishers. There is considerable variation in ethnic density, with most LSOA having a large majority of White-British residents, and few areas having over 50% non-White-British residents. Deprivation is positively skewed, meaning that most areas are not deprived, with a few very deprived areas. Finally, around 70% of LSOAs are within 5km, a parkrun, of a parkrun. Again, this is positively skewed with most LSAO being within 3-4km of their nearest event.

Table 2. Descriptive Statistics

3.3.2 Correlation Matrix

There is a negative correlation between particpation (run_count) and: deprivation (imd), distance to nearest parkrun (mn_distance), population density (pop_density) and ethnic density (perc_bme). Ethnic density is strongly postively correlated with population density, negatively correlated with percentage non-working age, and moderately positvely correlated with IMD suggesting that areas with higher ethnic density are more densely populated overall, more deprived and have higher percentage working age people.

The colour plots below show the participation levels for LSOA by deprivation and ethnic density for Urban and Rural areas (as defined in Bibby and Brindley (2013)). Yellow, green and blue indicate high, moderate and low levels of participation respectively.

The plot shows that participation is generally greatest in areas that have low levels of deprivation and low levels of ethnic density (bottom left), and lowest in areas with high levels of deprivation and high ethnic density (top-right). Areas with either high deprivation, or high ethnic density, tended to have low participation, suggesting that both are important independently. The relationship was robust to Urban Major areas and Urban Minor areas but did not hold in Rural areas where data was more limited. It is important to note that we do not control for other factors, such as the age of residents or the population density and there are therefore many confounding factors.

3.3.3 Poisson Model

The results of three Poisson regression models are shown in Table 3 below. All models include the control variables: population density, distance to nearest event and percentage of the population of non-working age. Model 1 includes IMD Score, Model 2 includes Ethnic Density and Model 3 includes both IMD and Ethnic Density. All coefficients are significant at the p<0.01 level.

Model 1 shows that, controlling for population density, distance to nearest event and age of population, areas with higher IMD (more deprived) have lower participation.

Model 2 shows that, with the same controls, areas with higher ethnic density (% non-White-British) have lower participation.

Model 3 shows that when both independent variables are included the coefficients decrease, suggesting that some of the effect previously attributed to deprivation is indeed due to lower participation in areas with higher ethnic density.

3.4 Discussion

Our findings show that more deprived areas and areas with higher ethnic density have lower participation rates. This effect persists after controlling for other area characteristics such as deprivation, access to events and population density. While our previous analysis (P. P. Schneider et al. 2019b) showed that participation in parkrun is lower in more deprived communities the present results suggest that a small part of the negative effect on participation previously attributed to deprivation can actually be attributed to ethnic density.

parkrun’s vision of creating a “healthier and happier planet by continually breaking down barriers to participation and bringing people together from all walks of life whenever they want to come along” (p.5) Cutforth (2017) has potential to improve population physical activity and therefore public health. Understanding the determinants of participation at the community level is a useful first step, but qualitative work to understand why this is the case is an obvious next step. Replicating this study in several years will enable parkrun to monitor trends in participation from different groups in society, and therefore the effectiveness of efforts to reach minority communities and those living in deprived areas.

3.5 Limitations

This analysis is ecological and therefore it is not possible to make conclusions at an individual level without risking an ecological inference fallacy. We have been careful throught to make conclusions at the level of the LSOA, rather than te individual. Nevertheless, given that the evidence at the individual level points to lower participation in organised sport by those from ethnic minority backgrounds (insert REF), we think it is likely that the same effect exists at the individual level.

Our dependent variable is the number of finishers by residents of each LSOA. This is a count variable where each walk or run finished is treated equally (e.g. 10 finishes by one person is equal to 10 people completing one event). We cannot draw inferences on the number of people who took part within each LSOA at some point in the year, but instead focus on the total finisher count. We do not expect that this will affect the core finding of the paper.

We use percent non-White-British as a crude proxy for ethnic density, and do not estimate participation by ethnic groups seperately. It is possible that there are significant differences between participation rates of different minority ethnic groups. Future analysis could look into which groups are most/least engaged in order to target efforts most efficiently. Furthermore we controlled for several variables which we thought would influence participation, it is possible that there are other confounding factors which have not been included.

3.6 Conclusion

parkrun is already in the process of increasing the number of events in deprived areas of England to encourage participation from disadvantaged groups. Our findings show, however, that in addition to deprivation and access, ethnic density is another important determinant of participation. Breaking down barriers to engagement in parkrun has the potential to improve overall population physical activity and therefore improve overall health and reduce health inequalities.

References

Bibby, Peter, and P Brindley. 2013. “The 2011 Rural-Urban Classification for Small Area Geographies: A User Guide and Frequently Asked Questions (V1. 0).” Government Statistical Service.

Cleland, Verity, Meredith Nash, Melanie J Sharman, and Suzi Claflin. 2019. “Exploring the Health-Promoting Potential of the ‘Parkrun’ Phenomenon: What Factors Are Associated with Higher Levels of Participation?” American Journal of Health Promotion 33 (1): 13–23.

Cutforth, Christopher. 2017. “So Much More Than a Run in the Park.” The Leisure Review, no. 84.

Fullagar, Simone, Sandra Petris, Julia Sargent, Stephanie Allen, M Akhtar, and Gozde Ozakinci. 2019. “Action Research with Parkrun Uk Volunteer Organisers to Develop Inclusive Strategies.” Health Promotion International, In–Press.

Goyder, EC, C Edmonds, A Sabey, D Lawrence, A Bullas, M Taylor, and S Potter. 2018. “P2 What Factors Predict Participation in a Mass Community Physical Activity Programme? The Case of the Five Sheffield ‘Parkruns’.” BMJ Publishing Group Ltd.

Haake, Steve. 2018b. “Parkrun: A New Model of Physical Activity for Large Populations?” The Sport and Exercise Scientist 57: 18–19.

R Core Team. 2018. R: A Language and Environment for Statistical Computing. Vienna, Austria: R Foundation for Statistical Computing. https://www.R-project.org/.

Reece, Lindsey J, Helen Quirk, Chrissie Wellington, Steve J Haake, and Fiona Wilson. 2019b. “Bright Spots, Physical Activity Investments That Work: Parkrun; a Global Initiative Striving for Healthier and Happier Communities.” Br J Sports Med 53 (6): 326–27.

Rowe, Nick, and Ross Champion. 2000. “Sport Participation and Ethnicity in England.” Headline Findings. Sport England, London.

Schneider, Paul P, Robert A Smith, Alice M Bullas, Thomas Bayley, Steve SJ Haake, Alan Brennan, and Elizabeth Goyder. 2019b. “Where Should New Parkrun Events Be Located? Modelling the Potential Impact of 200 New Events on Socio-Economic Inequalities in Access and Participation.” medRxiv, 19004143.

Stevinson, Clare, and Mary Hickson. 2013b. “Exploring the Public Health Potential of a Mass Community Participation Event.” Journal of Public Health 36 (2): 268–74.

Stevinson, Clare, and Mary Hickson. 2018. “Changes in Physical Activity, Weight and Wellbeing Outcomes Among Attendees of a Weekly Mass Participation Event: A Prospective 12-Month Study.” Journal of Public Health.

Wei, Taiyun, and Viliam Simko. 2017. R Package "Corrplot": Visualization of a Correlation Matrix. https://github.com/taiyun/corrplot.

Zhang, Dabao. 2017. “Rsq. Partial: Partial R-Squared for Generalized Linear Models.”