Do objective neighbourhood characteristics relate to residents’ preferences for certain sports locations? A cross-sectional study using a discrete choice modelling approach

Background The number of sports facilities, sports clubs, or city parks in a residential neighbourhood may affect the likelihood that people participate in sports and their preferences for a certain sports location. This study aimed to assess whether objective physical and socio-spatial neighbourhood characteristics relate to sports participation and preferences for sports locations. Methods Data from Dutch adults (N = 1201) on sports participation, their most-used sports location, and socio-demographic characteristics were collected using an online survey. Objective land-use data and the number of sports facilities were gathered for each participant using a 2000-m buffer around their home locations, whereas socio-spatial neighbourhood characteristics (i.e., density, socio-economic status, and safety) were determined at the neighbourhood level. A discrete choice-modelling framework (multinomial probit model) was used to model the associations between neighbourhood characteristics and sports participation and location. Results Higher proportions of green space, blue space, and the number of sports facilities were positively associated with sports participation in public space, at sports clubs, and at other sports facilities. Higher degrees of urbanization were negatively associated with sports participation at public spaces, sports clubs, and other sports facilities. Conclusions Those with more green space, blue space or sports facilities in their residential neighbourhood were more likely to participate in sports, but these factors did not affect their preference for a certain sports location. Longitudinal study designs are necessary to assess causality: do active people choose to live in sports-facilitating neighbourhoods, or do neighbourhood characteristics affect sports participation?


Background
In many European countries, sports clubs have a long tradition in facilitating organized sports [1]. Research has shown that club-based sports participation is associated with improved psychological and social health, above and beyond the positive effects gained from other individual forms of physical activity [2,3]. However, societal developments, such as individualization, have had an important influence on leisure culture and lifestyle, and have led to a decline in participation in sports clubs. Hence, the importance of sports clubs as the main providers of sports activities is changing [4,5].
These changes have led to the development of new opportunities to practice sports. Although both traditional club-based sports participation (also referred to as 'heavy sports') and non-club sports participation (or 'light sports') have increased over time, non-club sports participation has increased more rapidly [4][5][6]. Over the past decades, informal, unorganized, non-competitive, individual forms of sports such as running, cycling, and working out in the gym have become more popular, which has resulted in a greater variety of locations used for sports activities, including public spaces [6,7]. According to the Eurobarometer survey 2014, the largest share of European citizens aged 15 years and older engage in sports or physical activity in informal settings, such as parks and outdoor spaces (i.e., 40%), 23% engage in health centres/gyms or other sports centres, and only 13% participate in sports or physical activity as member of a sports club [8]. In the Netherlands, data showed that in 2012, 43% of Dutch sports participants participated in sports as members of a sports club, whereas 63% (also) participated in sports individually, in non-organized settings [9]. These numbers indicate the growing importance of public spaces as facilitators of sports participation.
Increasing sports participation is a common policy goal in Western societies [10], and to stimulate participation in sports, it is important to understand which locations are used for sports participation (and by whom), as well as what determinants play a role in the use of such locations. Previous studies that investigated the determinants of sports participation have mostly investigated whether people participate in sports and how frequently they participate. These studies showed associations between sports participation and individual determinants, such as socio-demographical and psychological factors and household factors [11][12][13]. Additionally, characteristics of the social and physical environment (such as safety, neighbourhood socio-economic status, social network, social cohesion, access to sports facilities, urban density, and attractiveness of the physical environment) were associated with sports participation [14][15][16]. Other studies focused on the organization of sports as the outcome measure. For example, research has shown that women are more likely to engage in informal (light) sports in commercial or alternative settings [4,5,7], and adults of higher social classes and with higher incomes are more likely to engage in non-organized sports [11,17].
Although these studies provide useful insights, we know much less about the factors that relate to the preferred locations for sports participation, such as public spaces, sports facilities (e.g., health centres/gyms) or traditional sports clubs. Of special interest are objectively measured neighbourhood characteristics, as these provide opportunities to develop environmental interventions (e.g., increase the availability of walking, jogging, and cycling trails, or improve the accessibility of sports facilities) that aim to increase sports participation.
To what extent one's residential neighbourhood is conducive to sports participation (e.g., by providing sport facilities or attractive public open spaces) may not only affect whether one participates in sports but also at which location one prefers to participate, e.g., in a city park or at a sports facility, and whether they will do so as a member of a sports club. For example, adults who live near public natural spaces with good walking trails (e.g., city parks) may be more likely to use public spaces to go for a run than adults who live near a large business park. Additionally, adults living in areas with many sports facilities and sports clubs may be more likely to practice sports in a sports club or in a private health centre/gym or swimming pool than adults living in areas with only a few available sports facilities and sports clubs.
This study aims to examine whether objective neighbourhood characteristics (i.e., proportions of land use and the availability of sports facilities in the residential neighbourhood) and objective socio-spatial neighbourhood characteristics (i.e., urban density, neighbourhood socio-economic status, and neighbourhood safety) relate to sports participation and preferred sports locations. To do so, we distinguish between 1) no sports participation, 2) sports participation in public open spaces, 3) sports participation as a member of a sports club (e.g., soccer and hockey), and 4) sports participation in public or private sports facilities (e.g., gyms and swimming pools). A discrete choice-modelling approach, which is a common approach in the literature on marketing and travel behaviour [18,19], but relatively new in the literature on sports participation and health, is used to model associations between neighbourhood characteristics and people's propensity to fall into one of these four different groups.

Participants
For this cross-sectional study, data were collected in six municipalities in the Netherlands (Amsterdam, Utrecht, Alphen aan den Rijn, Heerlen, Berkelland, and Roerdalen) in September 2014. These municipalities were selected based on their differences in size, population density, and geographical location in the Netherlands (i.e., more central vs. more peripheral) to have sufficient variation in presence, type, and accessibility of sports facilities.
Eighteen thousand adults (3000 adults per municipality), aged 18-80 years old, were randomly selected from municipal population registers. An information letter was sent to the home addresses of these adults, in which they were invited to participate in an online survey on sports participation and the use of sports locations. The online survey was used to obtain data on sports participation characteristics, principal sports location, and the personal characteristics of the respondents. Adults were asked to fill in their main type of sport, that is, the sport in which they participated most frequently during the 12 months prior to the online survey. Subsequently, they were asked where that sports activity mostly occurred (e.g., a public space, a sports club, or a registered sports facility) and if they participated as member of a sports club individually or in an informal group.
In total, 1663 respondents completed the survey (9.2% response rate). Data of respondents who met the following characteristics were excluded for analyses: participation in an inactive form of sports (e.g., bridge) (N = 20), participation in sports at home (N = 40) or at non-official sports facilities (e.g., community centres) (N = 64), unknown or incomplete data with regard to the postal code of their home address (N = 236), and other socio demographical characteristics (N = 69). Adults who were unable to participate in sports due to disabilities (N = 21) or health constraints (N = 12) were also excluded. Complete data were available for 1201 adults, and these respondents were included in further analyses.

Sports participation at specific locations
Based on survey questions about sports participation (at least once per month vs. less than once a month), the sports location that was used most often for sports participation over the past month (i.e., a public open space, a sports club, or a sports facility) and sports club membership (yes or no), the independent variable 'sports participation at specific locations' was categorized in four groups: 1) no sports participation (i.e., no sports participation at all or less than once per month), 2) sports participation in public spaces, 3) sports participation as a member of a sports club, using sports club facilities, and 4) sports participation at indoor (private or public) sports facilities, without club membership, e.g., health centres/gyms and swimming pools. Sports participation in public spaces included both unorganized sports (e.g., individually, with a friend, or in a small informal group), and organized sports (e.g., in a running group but without club membership). The four different types of sports participation are further referred to as no sports participation, sports participation in public spaces, sports participation at sport clubs, and sports participation at sports facilities.

Objective physical and socio-spatial neighbourhood characteristics
The independent variables were objectively measured neighbourhood characteristics and included land-use data, number of sports facilities, and socio-spatial data. Land-use data of respondents' home environments were obtained using ArcMAP 10.3.1. The coordinates of the 6-digit postal codes of respondents' home addresses were uploaded in ArcMAP. Mean coordinates of the 6digit postal codes (i.e., polygon features) were calculated, and Euclidean buffers of different sizes (i.e., 400, 800, 1600, and 2000 m) were drawn around these coordinates. The buffers were used to calculate the proportions of different types of land use (available from Statistics Netherlands, 2012) and the number of sports facilities (available from the Dutch Facility Monitor Sport (FMS), see [20]). The following types of land use were distinguished, as it is plausible that these may be related to sports behaviour: roads, facilities (e.g., churches, hospitals, shops, restaurants, and educational institutes), green space (e.g., parks, allotments, forest, and moorland), and blue space (e.g., rivers, lakes, and sea). These land use types were chosen based on associations shown in previous literature with physical activity [21,22] and sports participation [23], as well as based on their potential relation with sports behaviour, for instance due to the exercise friendly design of public spaces. As there is no consensus on buffer size for assessing associations between environmental characteristics and sports participation, we assessed models with various Eucledian buffers (i.e., 400, 800, 1600, and 2000 m). Based on the following reasons we decided to use the 2000 m buffers. First, the model with the 2000-m buffer had the best model fit compared to the models we have tested with other buffer sizes (McFadden R 2 , see Table 2). Second, the 2000-m buffer size corresponded best to our assumptions. We assumed that sports participants using the public space (for sports such as running et cetera) usually go further than their immediately neighbourhood of 400 or 800 m around their homes. For instance, a previous study found that runners not only use the public space in their neighbourhood, however, they also go outside the neighbourhood and out of town [24]. Moreover, our data showed that sports participants using sports clubs or private sports facilities on average travelled 3082 m (SD = 3.843 m) to their sports activities (see also [16] based on the same data), and those who use the public space for sports such as running will probably use an even wider area around their homes.
The socio-spatial data included urban density, neighbourhood socio-economic status (SES), and safety. Urban density was estimated as the average number of addresses within a radius of one square kilometre (available from Statistics Netherlands [25]). Three categories of address density were distinguished: rural (< 500 addresses per km 2 ), hardly to moderately urbanized (500-1500 addresses per km 2 ), and strongly to extremely urbanized (> 1500 per km 2 ). Objectively measured neighbourhood safety (on 4digit postal code level) was obtained from the 'Leefbaarometer 2.0' [26]. This measure includes items such as reported demolitions, crime and theft. Neighbourhood safety (mean = −0.002, SD = 0.13) was defined relative to the Dutch average score that had a standardized score of zero and was classified into three categories: safety level below the national average (score ≤ −0.05), approximately equal to the national average (score − 0.049-0.049) and above the national average (score ≥ 0.05). Neighbourhood SES (mean = −0.043, SD = 1.20), on 4-digit postal code level, was obtained from The Netherlands Institute for Social Research [27]. The SES scores were based on an aggregated indicator consisting of the following variables derived from Statistics Netherlands: average neighbourhood income, proportion of residents with a low income, proportion of residents with a low education level, and proportion of unemployed residents. We categorized the SES scores relative to the Dutch average score into three categories: neighbourhoods with an SES below (< −-1), approximately equal to (−-0.99-0.99) and above (> 1) the national average.

Confounders
We controlled for the following demographical characteristics: age, gender, education, having children who live at home (yes or no), and employment (yes or no). Education was classified into three levels based on the self-reported highest level of completed education: 1) lower education (i.e., no education, primary education, and lower professional education), 2) middle education (i.e., intermediate and higher general education), and 3) higher education (i.e., higher professional education and university).
Statistical analysis SPSS 23.0 was used to provide descriptive statistics on respondents' personal characteristics and objective neighbourhood characteristics (i.e., socio-spatial characteristics and proportions of land use). The influence of the discussed determinants on the use of different locations for sports participation (i.e., participants belonging to one of the four distinguished type of sports participant categories) was analysed through the application of a discrete choice modelling approach. In this approach, type of sports participant is considered to be a choice out of four alternatives available: non-participants, public space participants, sports club participants, and other sports facility participants. In discrete choice modelling, individuals are assumed to choose the alternatives that provide the highest utility [18,19]. The utility of alternative j (j = 1,…,J) for individual n can be represented by the following function: Here, is a vector of the observed characteristics (objective physical and socio-spatial neighbourhood characteristics as well as individuals' socio-demographic characteristics), which is the deterministic part of the utility function and in the context of this study only includes individualspecific variables, and an error term , which is the stochastic component of the utility function. The probability that individual n will choose alternative i is the probability that the utility derived from alternative i exceeds the utility of the other alternatives, which can be represented by: Under the assumption of independently and identically distributed (IID) error terms, the logit probabilities underlying the popular multinomial logit (MNL) model become: Central to the MNL model is the independence of irrelevant alternatives (IIA) property, which implies that the preference for an alternative is not affected by the inclusion or exclusion of other alternatives in the choice set. This property allows the use of independent standard normally distributed error terms and thus is fundamentally related to the IID assumption. However, many choice situations do not comply with IIA, as alternatives often share certain attributes that are unobserved by the researcher and therefore lead to correlations in the error terms of these alternatives. Additionally, in the case of sports location choice, such correlations can potentially be present. For example, being a 'public space participant' or a 'other sports facility participant' might be a shared preference of persons motivated to participate in individual sports but with a dislike for joining formal sports clubs. The Hausman-McFadden test [28] offers a procedure to test the IIA hypothesis for an MNL model, and applying this test, our data showed that the IIA property was violated for the estimated MNL model. In the presence of only individual-specific variables, the multinomial probit model (MNP) offers an attractive alternative model specification that can handle dependence across alternatives. The MNP model assumes that errors follow a multivariate normal distribution with mean 0 and covariance matrix ∑ [18]. The probabilities can be written as: where f(.) is the probability density function of the multivariate normal distribution.
For the estimation of our MNP model, the software platform 'R' , with the 'mlogit' package, has been used [29,30]. Table 1 shows socio-demographic characteristics of the four groups of non-sports participants and sports participants using different locations and environmental characteristics of their residential neighbourhoods. Table 2 shows the results of the raw (model 1) and adjusted (model 2) multinomial probit analyses. Green space, blue space, and number of sports facilities showed significant positive associations with sports participation in public space, in sports clubs, and in other sports facilities. Of these three variables, blue space showed the   High education was positively associated with sports participation in all three types of sports locations. Significant effects of neighbourhood characteristics (i.e., regarding the direction and size) were similar for each of the sports locations.   The significant correlation between the error terms of public space and other facilities shows that the multinomial probit model is an appropriate model to use to analyse the data. The correlation between public space and other sports facility participants indicates a shared preference for sports participants in public space and sports participants in other sports facilities based on non-measured factors, such as attitudes towards club membership. A potential explanation is that those who do not want to participate in organized sports because of obligations and fixed times prefer more flexible forms of sports, which may occur either in public spaces or other facilities such as gyms and swimming pools.

Discussion
This study adds to the existing literature on relationships between neighbourhood environments and sports participation by examining whether objectively measured residential neighbourhood characteristics (e.g., availability of sports facilities, green spaces) not only relates to sports participation but also to the preferred sports location of sports participants (e.g., an indoor sports facility or a city park). The results indicate that those with more green space, blue space or sports facilities in their residential neighbourhood were more likely to participate in sports, but it did not affect their preference for a certain sports location. Neighbourhood blue space in particular was associated equally strongly with sports participation at public spaces, sports clubs, and other sports facilities.
That natural environmental characteristics (i.e., green and blue space) are related to sports participation is in line with previous research that showed that attractive and liveable environments stimulate and invite people to exercise or participate in sports [31,32]. Moreover, such environments are associated with less experienced constraints to participating in sports more frequently [16]. The strong effect of green, and particularly blue space, on sports participation in all types of locations might be related to different types of sports that different natural areas are suitable for. For instance, previous research has shown that different green and blue areas (with different sizes) are related to different modalities and intensities of physical activity [21]. However, it would be of great interest to explore the importance of blue space for sports participation in more detail, as this may provide policy makers and urban designers with more accurate information on how to design both physical activity-and sports-friendly environments..
Moreover, higher proportions of green and blue space in the residential neighbourhood were also associated with a higher chance of sports participation in sports clubs and other public or private sports facilities. Additionally, sports facilities within a 2,000-metrem buffer around the home were positively associated with sports participation at sports facilities as well as with increased sports participation in public open spaces. Somehow, a sports-promoting neighbourhood environment increases the likelihood that residents participate in sports, but it does not affect their preference for a certain location. One explanation could be processes of 'modelling' and perceived social norms [33,34]. In a sports-facilitating neighbourhood, more people are out on the streets running, biking, walking, or travelling in their active wear to their sports. Viewing these active people may make residents more likely to participate in sports themselves ('modelling') or make them feel that being active is the social norm, which in turn increases their likelihood to become active. Another explanation may be that people with active lifestyles choose to live in (relatively 'blue' or 'green') neighbourhoods that facilitate such lifestyles; they may not necessarily use their neighbourhood facilities to practice sports but may use facilities in other parts of the city. Due to the cross-sectional design of this study, it remains unclear whether more green or blue space and sports facilities within the residential environment contribute to an increase in sports participation at public spaces, sports clubs, and other sports facilities or if adults who participate in sports choose to live in an environment that has a certain amount of sports facilities and green or blue space available. Hence, to better inform policy and intervention development, longitudinal or retrospective research is needed.
Regarding socio-spatial neighbourhood characteristics, we only found an effect of urban density. With higher urbanity levels, the chance of participating in sports declines in all types of sports locations. We found a similar effect for hardly to moderate urbanity levels, which was negatively associated with participation in public space and other (private or public) sports facilities. In other words, the higher the urbanity degree, the more likely it is that adults do not participate in sports. This is supported by existing evidence that sports participation rates are higher in rural areas in the Netherlands [15]. Literature has explained such urbanity differences in sports participation by determinants related to sociospatial neighbourhood characteristics, such as neighbourhood SES and safety [15,35]. However, these associations between safety, neighbourhood SES, and sports participation were not found in the current study. This inconsistency with previous research regarding the role of safety may be due to differences in measures of safety (i.e., objective vs. subjective measures). For example, it may be that perception of safety (which is mostly how safety was measured in previous studies) plays a different role in sports participation than objectively measured safety (which is how safety was measured in this study). Although neighbourhood SES was not associated with sports participation in any of the three types of sports locations, education on the individual level was associated with sports participation in all types of sports locations. This finding is in accordance with many other studies [11,36] and implies that the effects of green space, blue space, and sports facilities occur irrespective of education.
Strengths of this study include the novel approach within sports participation literature, focusing on an outcome variable that allows for a distinction between nonparticipants and different types of sports participants related to their sports location. Moreover, we used a strong analytical discrete choice approach and took various objectively measured environmental measures into account. We have checked for correlations between the independent variables and found no issues with collinearity.
The low response rate (9.2%) is a limitation of this study. The response rate was lowest in the largest municipalities of Amsterdam and Utrecht (Table 1), which is probably related to underrepresentation of adults of non-native Dutch origin in the sample (10.8% compared to 21.4% nationally in 2014) [37], who more frequently live in the larger cities. Furthermore, lower educated respondents were underrepresented as they represented 15.9% of our sample whereas their share of the total population is 33%, according to 2014 data. [25]. It is unclear how these issues affected the representativeness of the results. However, the high share of higher-educated and native Dutch adults did not lead to an overrepresentation of sports participants in this study, as according to the national Dutch figures for sports participation, 30% of adults do not participate at all or participate less than once per month in sports [9], which is similar to the percentage of non-participants in the current study. This suggests that a selection bias toward more sports minded respondents is unlikely. In addition, the fact that we correct for education level in our multivariate analysis implies that the effects of e.g. environmental factors represent the general effect across education levels in a reliable way. Another limitation of this study was that we were unable to include various buffers to compare with the 2,000-metrem buffer, due to data loss when using smaller buffers.
Future research should consider taking the use of multiple sports locations into account for adults who participate in more than one type of sport because this might have an impact on sports location choice. For instance, adults who participate in multiple sports may join a sports training programme at a sports club but may also engage in sports in a public space. With regard to buffer sizes used to calculate the proportions of land use available around one's home, future research could investigate if different buffer sizes can be used for different land use variables to better reflect the characteristics of these land-uses for sports participation behaviour. For example, it could be that the presence of (larger) green and blue environments lead to effects on sports participation that manifest themselves at a larger distance than the presence of roads and public facilities. In addition, to gain more insight into the importance of public space for sports participation, research should investigate which, how frequently, by whom and for what reasons specific locations in public space are used for sports participation. Moreover, more insight into the motivations of participants might add to the understanding of the use of different locations for sports. Longitudinal or retrospective research is needed to gain insight into the causality of the relationships between environmental characteristics (including the difference between green and blue spaces) and sports participation (in different locations).

Conclusions
More and more, sports participants find their way to alternatives to traditional sports clubs. This study is among the first to investigate the extent to which characteristics of the physical and socio-spatial environment are associated with the use of different locations for sports participation and whether these are traditional sports clubs, private or public sports facilities (such as gyms and swimming pools) or the more informal, flexible public spaces. We found that neighbourhood characteristics were similarly associated (i.e., their direction and size of the effect) with sports participation in public space, sports clubs and sports facilities. The more green space, sports facilities, and (especially) blue space available around the home environment, the greater the chance one participates in sports in each of the three types of sports locations. In other words, physical neighbourhood characteristics do not affect one's preferences for a certain sports location to use most often. Furthermore, higher urbanity levels were associated with lower chances of participating in sports in all three types of sports locations. For hardly to moderate urban areas this applies especially to lower probabilities participating in sports in public space and other (private or public) sports facilities. The associations we found indicate the possibly important role that the environment may have in facilitating and stimulating sports participation. For policy makers in the sports and health domains, this underlines the importance of facilitating and investing in exercise-friendly public spaces, whether participants will use these for sports activities, or join a (municipal subsidized) sports club or other sports facility. However, longitudinal study designs are necessary to assess the causality of these associations in order to adequately inform policy makers and urban designers for intervention development.