- Open Access
Using an agent-based model to simulate children’s active travel to school
© Yang and Diez-Roux; licensee BioMed Central Ltd. 2013
- Received: 10 December 2012
- Accepted: 22 May 2013
- Published: 26 May 2013
Despite the multiple advantages of active travel to school, only a small percentage of US children and adolescents walk or bicycle to school. Intervention studies are in a relatively early stage and evidence of their effectiveness over long periods is limited. The purpose of this study was to illustrate the utility of agent-based models in exploring how various policies may influence children’s active travel to school.
An agent-based model was developed to simulate children’s school travel behavior within a hypothetical city. The model was used to explore the plausible implications of policies targeting two established barriers to active school travel: long distance to school and traffic safety. The percent of children who walk to school was compared for various scenarios.
To maximize the percent of children who walk to school the school locations should be evenly distributed over space and children should be assigned to the closest school. In the case of interventions to improve traffic safety, targeting a smaller area around the school with greater intensity may be more effective than targeting a larger area with less intensity.
Despite the challenges they present, agent based models are a useful complement to other analytical strategies in studying the plausible impact of various policies on active travel to school.
- School Location
- Baseline Scenario
- Traffic Safety
- Active Travel
- Travel Mode
Active travel to school (ATS) constitutes a substantial portion of children’s overall physical activity  and ATS is associated with higher overall physical activity [2–10]. ATS can decrease traffic and pollution, reduce children’s dependence on parents, improve social interactions, and promote healthy life styles which may be maintained into adulthood . Despite the multiple merits associated with ATS, today fewer than 15% US children and adolescents walk or bicycle to school, compared to about 50% in 1969 . Recently, a number of intervention programs [13–15] have been applied to encourage ATS. However, intervention studies are in a relatively early stage and evidence of their effectiveness over long periods is limited .
Research has identified a number of factors that are associated with ATS including the characteristics of children and families [5, 10, 17, 18] (such as age and fitness of the child, neighbourhood safety and traffic safety, and household transportation options), features of schools [5, 10, 19, 20] (such as distance from the household, and school bus policy), and features of neighbourhoods [5, 10, 17, 21, 22] (such as urban form and social norms). Generally, boys, Hispanic and African American children, children from lower SES families, and children attending public schools, living close to schools, and living in neighbourhoods with walking or bicycling friendly environments are more like to actively travel to school .
Long distance to school is the most common barrier to ATS among US children . Distance to school is a function of the number of schools, their location, and the definition of catchment areas. A longer distance to school has been found to be associated with lower ATS, especially when the distance is less than one mile [25, 26]. It has been estimated that 47% of the decline in walking to school that occurred between 1969 and 2001 is attributable to increases in the distances between households and schools  resulting at least in part from suburbanization and from the creation of larger school districts. Communities across the USA have been shifting away from small, neighbourhood schools to larger schools located in areas where housing densities are low and road networks lack connectivity . Studies [29–31] have reported that high population density, and small enrollment size are associated with higher percentage of ATS. However, understanding of how the definition of school catchment areas affects ATS is limited. It is also unclear how school catchment area definition, school size, and population density may jointly affect ATS.
Traffic-related danger is another major barrier to ATS . However evidence of the impact of changing traffic safety levels is mixed [10, 32]. The Safe Routes to School  program is a nationwide initiative to improve traffic safety. One way to increase traffic safety is to modify road conditions such as creating crossings, improving sidewalks, and implementing traffic calming measures. An important decision is where interventions should be targeted in order to get the maximum benefit with limited resources, for example, should the investment focus on a more intense intervention on a smaller area, or a less intense intervention on a bigger area?
The majority of research on factors related to children’s school commuting has relied on statistical models applied to observational data. Most analyses have been cross-sectional, have investigated only a limited number of characteristics, and have not fully explored interactions between factors. In general, existing work fails to capture the dynamic interactions between individuals (e.g. the effects of social networks) or the dynamic interactions between individual and environmental factors (e.g. feedback loops between the walking behavior of individuals and the safety of an area). In addition, features of built and social environments may dynamically interact (e.g. mixed land use and land use density may affect levels of safety). Existing research may provide an incomplete basis for policy because of its inability to draw inferences regarding the plausible impact of interventions in the context of dynamic systems.
Agent-based modeling (ABM) has been increasingly applied to population health problems [33, 34] because of its utility in studying numerous agents specified at various scales, capturing the complex dynamic relationships, and accounting for feedbacks between individuals and between individuals and their environments. By focusing on explicit processes and considering how the outcomes of the functioning of the systems would change if features were modified, ABMs can be used to conduct various “what-if” experiments relevant to policy or intervention . To our knowledge, applications of ABMs to the study of how the social and built environments shape people’s active travel behavior remains limited [36–38], and has not yet focused on children.
We developed an ABM to simulate children’s school travel behavior within a hypothetical city. The model was calibrated to existing data. We illustrate the utility of this approach by applying the model to two outstanding questions in ATS: (1) the impact of school location, catchment area definition, school size, and population density; and (2) the impact of changing traffic safety levels.
The model was developed in Java and Repast . It includes a hypothetical city with a road network, schools, households and children. The city has a size of 8 km × 8 km, and has a grid road network with 41 vertical and 41 horizontal roads and evenly sized blocks. Each road has 40 segments separated by road intersections. Each segment (composed of 10 cells) denotes 200 m in reality. Each cell (denoting a location with the size of 20 m × 20 m in reality) can be occupied by a household, a school or be empty. The city includes four middle schools and 3000 households. Households are randomly distributed across the city. Each household includes one child and the child is enrolled in one of the four middle schools. Each household has a level of concern towards traffic safety (denoted by C t ) with respect to children walking to school. C t is a threshold value: if traffic safety is perceived to be above that value, the household may allow the child to walk to school. Each child has an attitude towards walking (denoted by A w ): a higher value of A w indicates a higher probability of walking to school (assuming other conditions are met, as discussed below). Both C t, and A w are assigned at model initialization using a random draw from a uniform distribution between 0 and 1. For the purposes of this illustration, C t, remains constant but A w is updated over time based on various feedbacks. The model is time-discrete with each step being one day.
Where d is the distance from household to school and β is the distance decay parameter. The distance decay parameter was calibrated by comparing model output to empirical distributions as described below. Condition 2 implies that the child may walk to school if the child’s attitude towards walking is high and/or if the distance to the school is short (resulting in a higher P d ). The use of the sum allows us to account for the fact that children with a high attitude towards walking may walk even if the distance to school is long, and vice versa, when the distance is very short even children with a low attitude may walk to school.
Where W is the total number of people who walk by that cell each day, which includes the children walking to school and the background walkers (the number of background walkers is assumed to be constant across the city). Prior research has shown that W-0.6 is a reasonable estimate of the probability that a pedestrian is struck by a car, given W walkers in the area . We therefore use 1-W-0.6 as a measure of traffic safety.
Where and A w are the attitudes towards walking today and yesterday, respectively. N is the total number of children enrolled in the school, and and N w are the total number of children enrolled in the school who travelled to school by walking today and yesterday, respectively. According to this formula, if within a given school, more children walk to school today than yesterday, then the attitudes towards walking of all the children within the school will increase.
To investigate our research questions, we contrasted the percent of children who walk to school under various scenarios. For each scenario the model was run until the number of children who walk to school remains constant (approximately 5–10 days). For each scenario, we report the mean percentage of children walking to school and the distribution of children by the distance from their households to school over 20 simulations.
Percent of children who walk to school, by distance, from model predictions and 2009 NHTS data
Summary statistics generated by the model (baseline scenario, with 20 simulations)
Values from NHTS data
> = 2
The impact of school location, catchment area definition, school size and population density
Features and simulated results for four scenarios with different school location and catchment area policies
Definition of catchment area
Children walking to school (%)
Percent distribution of children by distance to school (in miles)
> = 2
Children attend nearest school
Randomly located across the whole city
Same as baseline
The city is divided to four school zones of equal size, and the school is randomly located within the zone
Attend school within the same zone
Same as S2
Same as baseline
Changing the traffic safety level
To investigate this question we used our model to contrast three strategies: (1) increase the safety of all roads within 0.5 miles from the school with value of a; (2) increase the safety of all roads within 0.71 miles from the school with value of a/2 [we use a/2 because the size of area within 0.71 miles is twice the size of the area within 0.5 miles]; and (3) increase the safety of all roads within one mile from the school with value of a/4 [we use a/4 because the size of area within one mile is four times the size of area within 0.5 miles]. We adjust the intensity of the intervention so that the intensity-area product is equivalent across scenarios (the intervention per unit area is more intense when a smaller area is targeted than when a larger area is targeted). In order to investigate effects across a range of values, three levels of a (0.1, 0.2 and 0.3) were examined. Because the value of the safety ranges between 0 and 1 by definition, these three values cover a reasonable range of intensity of the intervention.
The impact of school location, catchment area definition, school size and population density
Simulated percent of children walking to school for scenarios with different numbers of schools and households in the city
Number of schools
Number of households
Changing the traffic safety level
The main contribution of our study is to illustrate the utility of agent-based models to gain insight into fundamental dynamics driving patterns of active travel to school. We also show how the model can be used to identify the plausible impact of various policies or interventions. Our results indicate that to maximize the percent of children who walk to school the school locations should be evenly distributed over space and children should be assigned to the closest school. We also illustrate the beneficial impact of smaller catchment areas and higher population density. In the case of interventions to improve traffic safety, targeting a smaller area around the school with greater intensity may be more effective than targeting a larger area with less intensity.
The impact of school location, catchment area definition, school size and population density
We used our model to illustrate how ABMs can help shed light on the impact of policies affecting school location and school catchment area which jointly determine the distribution of the distances to schools. The baseline scenario has the highest percentage of children walking to school because in this scenario, schools are distributed evenly across the city, and this minimizes the distances from households to schools across the whole sample. The percent of children walking to school is higher in scenario S3 than in scenario S1, because both scenarios allow the children to attend the school which is nearest to their household, but in S3, schools were located randomly within the zone while in S1, schools were located randomly within the whole city. Thus, on average the probability that schools are too close to each other (resulting in large distances for some children) is smaller in S3 than S1. Scenario S3 is also better than S2, because although the rules for the school location are the same in both scenarios, in S3 the catchment area is based on distance to the school resulting in a distribution of distance to school that is more skewed towards shorter distances. Walking levels in scenario S1 are very similar to S2. The implication is that the disadvantages in terms of ATS of the school location strategy for S1 (compared to S2) can be compensated by an advantageous school catchment policy that requires students to attend the closest school.
The school location and catchment polices in our scenarios are clearly unrealistic. However, a major advantage of applying ABMs to these questions is the ability to predict the joint impact of policies affecting multiple different domains (such as school location and catchment areas) in quantitative terms. The exploration of various more complicated “what if “scenarios can allow identification of the optimized school location and catchment area strategies  appropriate for different contexts, for example, across different spatial patterns of land use mix or in cities with varying spatial distributions of crime or traffic safety.
It is not surprising that with a fixed total number of households, a larger number of schools will increase the percent of children who engage in ATS. This is because more schools leads to smaller catchment areas and shorter distances to schools. This can also account for the diagonal pattern observed in Table 3: when the school size is constant (750 children each school), a higher population density means smaller catchment areas. An interesting insight arising from our model is that with a fixed number of schools (that is fixed school catchment areas), a higher population density also results in higher percent of children using ATS. This is because, when catchment areas are held constant, a higher population density will result in higher absolute number of children walking to school. This increase in walkers increases the traffic safety level which in turn results in more walking, leading to a reinforcing feedback loop.
Using the model, we provided a possible explanation of the mechanism by which high population density could increase ATS. Of course, this explanation is based on the assumptions encoded in our model including the assumption that the traffic level is only affected by the number of walkers which may be incorrect. For example, it is possible that higher population density will also result in higher absolute number of parents who drive their children to school which could have a detrimental effect of traffic safety. Nevertheless, our results illustrate the potential utility of ABMs, more specifically the ability of these models to identify dynamic mechanism and explain the synergistic effects of various features through modelling of simple feedbacks.
Changing the traffic safety level
Our results show that strategies that target a smaller area around the school with greater intensity are more effective in increasing ATS than strategies that target a broader area with lower intensity. The implication of this finding is that investments to increase safety should begin in the areas closer to school and then extend to remote areas, in order to maximize the benefit. At least three different processes could account for this result: (1) increasing traffic safety close to the school with greater intensity results in a greater improvement because it directly and intensely influences children living shorter distances from school who are more likely to respond to changes in safety (children living further from the school will not be affected as much because distance becomes a barrier to walking regardless of safety). (2) increasing traffic safety closer to the school is beneficial to all children because all children using ATS must walk through the area around the school; (3) the strategy also has an important indirect effect because it changes the “social norms” regarding ATS within the school, which in turn affects the attitude towards walking for all children within the school regardless of whether they live in the area targeted for intervention. This means that children living out of the area targeted for safety improvements can benefit as well.
Among the three processes, the first two could be inferred through reasoning, but identification of the third process requires modelling. By including a feedback mechanism in the model through which each child’s attitude towards walking is affected by school social norms, we allow children to influence each other. Another example of a dynamic process in our model is the dynamic interaction among travel modes especially between walking and driving. Some of the traffic around the schools is contributed by the parents themselves . When many parents drive their children, other parents may feel obliged to do the same to avoid traffic related risks to their children if they walk or cycle [43, 44]. On the other hand, if more children walk or bicycle to school, the traffic around the school may be alleviated and the level of traffic safety may be increased, which may in turn encourage more children to join ATS. The explicit implementation of these dynamic processes in the model allows us to test intervention strategies that could trigger these kinds of feedback loops, such as an education campaign to increase the children’s attitude towards active travel. Modelling approaches can allow us to identify the existence of a “tipping point” beyond which the feedback triggers a trend in the system in the direction we desire. We can also investigate ways to maintain the desired tendency. The implication is that we should take advantage of these “self-reinforcing” feedbacks to encourage ATS and maximize the benefits from limited investment.
The cases presented in this paper are very simple illustrations, but already begin to shed light on some fundamental dynamics that may be operating. The use of the model to answer more specific research questions will require refinements. For example, bicycling and using a school bus (or other public transportation) may be added as additional travel modes. It is also possible to modify the model by allowing for greater heterogeneity across children and households. For example, marked variations in ATS among children by age, gender and income level have been reported . These refinements may be especially important if the research questions are related to differences in the distribution of ATS across various populations. Another important improvement is to add more dynamic processes to the model. Currently, two dynamic processes were implemented for illustrative purposes. More dynamic interactions may be identified as important to ATS including for example, the interaction between the travel modes of various family members. If there is more than one child within a family enrolled in the same school, the probability of ATS may increase because children can walk together to be safe from traffic and crime . Also, a parent who walks to his/her workplace may be more likely to chaperone a child who is walking if spatially convenient. In addition, the impact of variations across features of the built and social environment such as land use mix and density, safety from crime and weather could be explored.
A major challenge in developing these models is limited knowledge on the specific processes that are involved in decision-making regarding mode choices for travel to school. Another challenge pertains to making these models empirically grounded given that data to support some parameters may be unavailable or may be impossible to obtain. Despite these challenges, our examples illustrate the potentialities of ABMs to gain insights that cannot be obtained from observational data and to generate hypotheses that can be tested with other standard approaches including usual observational and experimental studies.
This research was funded by grants from the NIH, National Heart, Lung and Blood Institute (Grants R21-HL106467 and R01-HL071759, and Office of Behavioral and Social Sciences Research (Grant #HHSN276200800013C).
- Tudor-Locke C: Omission of active commuting to school and the prevalence of children’s health-related physical activity levels: the Russian longitudinal monitoring study. Child Care Health Dev. 2002, 28: 507-512. 10.1046/j.1365-2214.2002.00295.x.View ArticleGoogle Scholar
- Sirard J: Physical activity and active commuting to elementary school. Medicine & Science in Sports & Exercise. 2005, 37 (12): 2062-2069. 10.1249/01.mss.0000179102.17183.6b.View ArticleGoogle Scholar
- Cooper AR: Physical Activity Levels of Children Who Walk, Cycle, or Are Driven to School. Am J Prev Med. 2005, 29 (3): 179-184. 10.1016/j.amepre.2005.05.009.View ArticleGoogle Scholar
- Heelan K: Active commuting to and from School and BMI in elementary school children - preliminary Data. Child: Care, Health and Devevelopment. 2005, 31 (3): 341-349. 10.1111/j.1365-2214.2005.00513.x.Google Scholar
- Lee M, Orenstein M, Richardson M: Systematic review of active commuting to school and children’s physical activity and weight. J Phys Act Health. 2008, 5 (6): 930-949.Google Scholar
- Saksvig B: Travel by walking before and after school and physical activity among adolescent girls. Arch Pediatr Adolesc Med. 2007, 161 (2): 153-158. 10.1001/archpedi.161.2.153.View ArticleGoogle Scholar
- Alexander L: The broader impact of walking to school among adolescents: seven-day accelerometry based study. Br Med J. 2005, 331: 1061-1062. 10.1136/bmj.38567.382731.AE.View ArticleGoogle Scholar
- Cooper A: Active travel to school and cardiovascular fitness in Danish children and adolescents. Med Sci Sports Exerc. 2006, 38: 1724-1731. 10.1249/01.mss.0000229570.02037.1d.View ArticleGoogle Scholar
- Landsberg B: Associations between active commuting to school, fat mass and lifestyle factors in adolescents: the Kiel Obesity Prevention Study (KOPS). Eur J Clin Nutr. 2008, 62 (6): 739-747. 10.1038/sj.ejcn.1602781.View ArticleGoogle Scholar
- Davison K, Werder J, Lawson C: Children’s active commuting to school: current knowledge and future directions. Prev Chronic Dis. 2008, 5: A100.Google Scholar
- US Department of Health and Human Services: Physical activity and health: a report of the surgeon general. 1996, Atlanta, GA: U.S. Department of Health and Human Services, Centers for Disease Control and Prevention, National Center for Chronic Disease Prevention and Health PromotionGoogle Scholar
- McDonald NC, School Travel US: An Assessment of Trends. Am J Prev Med. 2011, 41 (2): 146-151. 10.1016/j.amepre.2011.04.006.View ArticleGoogle Scholar
- CDC: KidsWalk-to-School. 2010, http://www.cdc.gov/nccdphp/dnpa/kidswalk/.Google Scholar
- National Center for Safe Routes to School: Safe Routes to School. 2010, http://www.saferoutesinfo.org.Google Scholar
- Walking School Bus: Walking School Bus. 2010, http://www.walkingschoolbus.org/.Google Scholar
- Chillón P: A systematic review of interventions for promoting active transportation to school. Int J Behav Nutr Phys Act. 2011, 8 (1): 10-10.1186/1479-5868-8-10.View ArticleGoogle Scholar
- Trapp GSA: Increasing Children’s Physical Activity: Individual, Social, and Environmental Factors Associated With Walking to and From School. Health Educ Behav. 2012, 39 (2): 172-182. 10.1177/1090198111423272.View ArticleGoogle Scholar
- Mammen G: Understanding the drive to escort: a cross-sectional analysis examining parental attitudes towards children’s school travel and independent mobility. BMC Publ Health. 2012, 12 (1): 862-10.1186/1471-2458-12-862.View ArticleGoogle Scholar
- McMillan TE: Urban Form and a Child’s Trip to School: The Current Literature and a Framework for Future Research. J Plan Lit. 2005, 19 (4): 440-456. 10.1177/0885412204274173.View ArticleGoogle Scholar
- McMillan TE: The relative influence of urban form on a child’s travel mode to school. Transportation Research Part A: Policy and Practice. 2007, 41 (1): 69-79. 10.1016/j.tra.2006.05.011.Google Scholar
- Mitra R, Buliung RN: Built environment correlates of active school transportation: neighborhood and the modifiable areal unit problem. J Transp Geogr. 2012, 20 (1): 51-61. 10.1016/j.jtrangeo.2011.07.009.View ArticleGoogle Scholar
- Zhu X: Context-Specific Correlates of Walking Behaviors to and From School: Do They Vary Across Neighborhoods and Populations?. J Phys Act Health. 2011, 8 (S1): S59-S71.Google Scholar
- Babey SH: Sociodemographic, Family, and Environmental Factors Associated with Active Commuting to School among US Adolescents. J Public Health Pol. 2009, 30 (S1): S203-S220.View ArticleGoogle Scholar
- CDC: Barriers to Children Walking to or from School - United States. Morb Mortal Wkly Rep. 2005, 54 (38): 949-952.Google Scholar
- Wong BY-M, Faulkner G, Buliung R: GIS measured environmental correlates of active school transport: A systematic review of 14 studies. Int J Behav Nutr Phys Act. 2011, 8: 39-10.1186/1479-5868-8-39.View ArticleGoogle Scholar
- McDonald N: Travel and the social environment: Evidence from Alameda County, California. Transp Res Part D: Transp Environ. 2007, 12: 53-63. 10.1016/j.trd.2006.11.002.View ArticleGoogle Scholar
- McDonald NC: Active Transportation to School: Trends Among U.S. Schoolchildren, 1969–2001. Am J Prev Med. 2007, 32 (6): 509-516. 10.1016/j.amepre.2007.02.022.View ArticleGoogle Scholar
- Ewing R, Greene W: Travel and environmentalimplications of school siting. 2003, Washington, DC: US Environmental Protection AgencyGoogle Scholar
- Falb MD: Estimating the Proportion of Children Who Can Walk to School. Am J Prev Med. 2007, 33 (4): 269-275. 10.1016/j.amepre.2007.05.005.View ArticleGoogle Scholar
- Nelson N: Active commuting to school: How far is too far?. Int J Behav Nutr Phys Act. 2008, 5 (1): 1-10.1186/1479-5868-5-1.View ArticleGoogle Scholar
- McDonald NC: Critical Factors for Active Transportation to School Among Low-Income and Minority Students: Evidence from the 2001 National Household Travel Survey. Am J Prev Med. 2008, 34 (4): 341-344. 10.1016/j.amepre.2008.01.004.View ArticleGoogle Scholar
- Pont K: Environmental correlates of children’s active transportation: A systematic literature review. Health Place. 2009, 15 (3): 849-862. 10.1016/j.healthplace.2009.02.002.View ArticleGoogle Scholar
- Maglio PP, Mabry PL: Agent-Based Models and Systems Science Approaches to Public Health. Am J Prev Med. 2011, 40 (3): 392-394. 10.1016/j.amepre.2010.11.010.View ArticleGoogle Scholar
- Luke DA, Stamatakis KA, Systems Science Methods in Public Health: Systems Science Methods in Public Health: Dynamics, Networks, and Agents. Annu Rev Public Health. 2012, 33 (1): 357-376. 10.1146/annurev-publhealth-031210-101222. p. nullView ArticleGoogle Scholar
- Bonabeau E: Agent-based modeling: Methods and techniques for simulating human systems. Proc Natl Acad Sci. 2002, 99: 7280-7287. 10.1073/pnas.082080899.View ArticleGoogle Scholar
- Lu Y, Kawamura K, Zellner M: Exploring the Influence of Urban Form on Work Travel Behavior with Agent-Based Modeling. Transportation Research Record: Journal of the Transportation Research Board. 2082, 2008: 132-140.Google Scholar
- Yang Y: A spatial agent-based model for the simulation of adults’ daily walking within a city. Am J Prev Med. 2011, 40 (3): 353-361. 10.1016/j.amepre.2010.11.017.View ArticleGoogle Scholar
- Yang Y: Exploring walking differences by socioeconomic status using a spatial agent-based model. Health Place. 2012, 18 (1): 96-99. 10.1016/j.healthplace.2011.08.010.View ArticleGoogle Scholar
- Repast: Repast. 2010, http://repast.sourceforge.net/repast_3/.Google Scholar
- Jacobsen PL: Safety in numbers: more walkers and bicyclists, safer walking and bicycling. Inj Prev. 2003, 9 (3): 205-209. 10.1136/ip.9.3.205.View ArticleGoogle Scholar
- McDonald NC: School Siting : contested visions of the community school. J Am Plann Assoc. 2010, 76 (2): 184-198. 10.1080/01944361003595991.View ArticleGoogle Scholar
- Dumbaugh E, Frank LD: Traffic safety and safe routes to schools: Synthesizing the empirical evidence. Transportation Research Record: Journal of the Transportation Research Board. 2009, 2007: 89-97.Google Scholar
- Fyhri A: Children’s active travel and independent mobility in four countries: Development, social contributing trends and measures. Transport Policy. 2011, 18 (5): 703-710. 10.1016/j.tranpol.2011.01.005.View ArticleGoogle Scholar
- Bradshaw R: School children’s travel. Geography. 2001, 86 (1): 77-78.Google Scholar
- Chillón P: Six-Year Trend in Active Commuting to School in Spanish Adolescents. Int J Behav Med. 2012, Oct 5 [Epub ahead of print]Google Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.