Abstract
Tourists’ activity patterns result from complex interactions between time–space constraints and cognitive, social, cultural, and emotional factors. Accordingly, tourists’ intradestination activity is studied today from multiple perspectives. Yet knowledge regarding the interrelationships between these factors is limited. The current research aims to contribute to the bridging of this gap, by studying tourists’ activity patterns and the time–space resource allocation decisions they reflect. Using a smartphone application, 384 tourists’ activity days in the cities of Tel-Aviv and Jerusalem were tracked and characterized in terms of resource allocation patterns, identifying three distinct behavioral patterns. The effects of various personal and external factors on group membership were then assessed. The results uncover complex relationships between various decision-making dimensions, where interdimensional interactions occur at multiple spatial scales and decision-making instances. These findings suggest that activity patterns emerge from a decision time–space in which effects are contingent upon each other across scales and behavioral dimensions.
Introduction
The rise of tourism as an important economic activity has led to the formation of tourism studies as a distinct research field (Archer, Cooper, and Ruhanen 2005; Butler 1980; Correia and Kozak 2017; Mayer et al. 2010). Tourism geographies had attracted much of the attention in this field because of this activity’s spatially explicit nature. The introduction of advanced tracking technologies, such as the global positioning system (GPS), over the last decade has made observing tourists’ behaviors with much detail and precision possible, thus promoting more research on intradestination behaviors (De Cantis et al. 2016; Grinberger, Shoval, and McKercher 2014; Hallo et al. 2012; McKercher et al. 2012, 2015; Shoval et al. 2011; Shoval and Isaacson 2007a, 2007b; Shoval, Schvimer, and Tamir 2017; Zhang, Huang, and Li 2017; Zillinger 2007; Zoltan and McKercher 2015). These studies are naturally related to the study of spatial behavior, a well-developed area of research within the field of geography (e.g., Walmsley and Jenkins 1992; Xiao-Ting and Bi-Hu 2012; Zhao et al. 2015; Zillinger 2007).
Behavioral geography, as a distinct research field, offers many insights regarding the motives and causes behind individuals’ daily spatiotemporal choices (Golledge 2008). These insights relate to the nature of the decision-making process (Golledge and Stimson 1997; Strauss 2008; Wolpert 1964), its cognitive foundations (Downs and Stea 1973; Kitchin 1994; Lynch 1960; Manley 2016; Manley, Addison, and Cheng 2015), the effects of external conditions and time–space constraints (Hägerstraand 1970; Long and Nelson 2013; Neutens, Schwanen, and Witlox 2011; Pred 1977), and also to the emotional, affectual, and cultural dimensions of behavior (Greenberg, Raanan, and Shoval 2014; McQuoid and Dijst 2012; Pile 2010). The relations between tourism geographies and time–space behavior are evident in a plethora of studies that incorporate behavioral-geographic ideas, such as the time–space budget (Xiao-Ting and Bi-Hu 2012; Zillinger 2007) and the cognitive map (Chang 2015; Guy, Curtis, and Crotts 1990; Pearce 1977; Walmsley and Jenkins 1992) into the study of intradestination behaviors. Yet, this reliance implies that the study of intradestination behavior suffers from the same constraints that plague behavioral studies in geography. One such issue is derived from overt spatial behavior being the result of a complex array of interrelationships between individual-specific and environmental factors (Chen et al. 2016; Golledge and Stimson 1997). Consequently, assessing the relative weight of each factor in the formation of behavioral patterns and identifying contingencies among effects present a great challenge within this type of studies (Chen et al. 2016).
In this context, the development of advanced tracking technologies holds a unique potential, allowing researchers to capture individuals’ behavioral patterns with much detail. Several studies over the last few years have shown how data collected via these methods can be used to uncover various effects within time–space behavior (in both tourism and behavioral studies; Grinberger, Shoval, and McKercher 2014; Manley 2016; Manley, Addison, and Cheng 2015; McKercher et al. 2012; Shoval et al. 2011). The research presented here joins these attempts, carried while considering multiple conceptual points of view on spatial behavior and using high-resolution data on tourists’ behaviors. As such, its main objective is to produce new insights regarding contingencies between the factors that affect tourists’ spatial behavior. This study uses a unique dataset that captures the behavior of tourists in Israel during their entire visit. A comparative analysis of this dataset identifies the similarities and differences in behaviors over two main urban destinations. As such, the analysis allows distinguishing between factors that exert a general influence, to those that are site- or individual-specific by nature (or combine both). Thus, this research contributes new knowledge regarding the hierarchical character of influences on time–space behavior, and contingencies that exist within it.
The article continues as follows. First, we present a review of related literature to identify the relevant theoretical background for this study. The second section introduces the database and discusses processing and analysis steps. The next section details the empirical results. The paper then concludes with a discussion of the findings, their theoretical significance, and the implication of these for the study of tourist and general time–space behaviors.
Literature Review
Time Geography
This article relies mainly on the research framework known as time–geography (TG). TG is a conceptual framework for the analysis of time–space behavior, developed by Torsten Hägerstraand and his students (Ellegård and Svedin 2012). The main argument within TG is that the study of spatial behavior must include a consideration of time, as each instance of such behavior (e.g., participation in activity, movement) requires the allocation of time resources (Hägerstraand 1970; Neutens, Schwanen, and Witlox 2011). Given that the human body is indivisible, allocation of time resources to one type of activity necessarily implies that less time would be available for other activities (Neutens, Schwanen, and Witlox 2011; Pred 1977). For example, for tourists, allocating more time for a stationary activity (e.g., visiting an attraction) means less time would be available to explore the destination. As this latter activity requires the consumption of a “space” resource, TG identifies time and space as substitute resources. Furthermore, the “amount” of time–space resources available to the individual is identified to be limited by several types of constraints (Hägerstraand 1970; Long and Nelson 2013; Pred 1977): capability constraints, which relate to the physical needs and abilities of the individual; coupling constraints, which result from the need to coordinate activities with the schedules of other individuals and objects; authority constraints, emerging from the set of laws and norms that apply within a specific time–space “domain.”
The implications of such constraints, when interacted with goals and purposes that guide spatial behavior, is that behavioral patterns are organized around a set of spatially and/or temporally fixed activities (known as “stations”; Neutens, Schwanen, and Witlox 2011; Pred 1977). Activity (mobile or stationary) within the time–spaces near these stations are constrained by a limited activity time budget, set by the temporal attributes of the stations. This conceptualization of behavior informs the development of several fundamental time–geographic concepts. One such concept is the time–space prism, a formal and graphical representation of accessible activity space at each time instance, as determined by the spatiotemporal attributes of two subsequent stations and other constraints on mobility, such as speed limits (Lenntorp 1976; Miller 1991; Neutens, Schwanen, and Witlox 2011). The prism encompasses all time–space locations for which visit probability is above zero for the individual. Its graphical representation is a three-dimensional polygon, where the X and Y axes represent space, and time being the third dimension (see Figure 1). When this polygon is projected onto a 2-dimensional space, the potential path area (PPA) is identified—a representation of an individual’s maximum geographical reach given a set of spatiotemporal constraints. The individual’s actual time–space path (see Figure 1), that is, the set of time–space locations visited during a specific time-period, must be traversed within the boundaries of this area.

Basic time–geographical entities.
TG presents a valuable foundation for explorations in many mobility-related fields, and studies using it have managed to produce important insights regarding issues such as the effects of gender on behavior (Kwan 1999), the relations between transport and activity (Ellegård and Svedin 2012), and even wildlife movement (Long and Nelson 2015). Yet, the reliance of TG on a very constrained view of time–space behavior arguably limits its relevance to the study of tourist behavior, given tourism’s unique nature. Yet, tourists’ behaviors are still very much constrained by many time–space elements, such as their personal capabilities, attractions and facilities operation hours, the capabilities of group members (when traveling in group), and the time available for activities (number of days, hours per day). The applicability of TG to the case of tourist activity is evident in a growing number of studies that use it to investigate the geographies of tourism (e.g., Grinberger, Shoval, and McKercher 2014; Shoval and Isaacson 2007a, 2007b; Xiao-Ting and Bi-Hu 2012; Zillinger 2007). The uniqueness of the tourist case within the TG framework lies in the enhanced freedom of choice that tourists enjoy within the boundaries of the prism. In other words, the tension between potential behavior (represented by the prism) and its specific realization requires more attention in TG-based tourism studies (Grinberger, Shoval, and McKercher 2014).
This realization is dependent on many elements; not all of them are spatiotemporal in nature. For example, multidestination travel is a phenomenon that was identified to be highly contingent on the spatial configuration of the destination (Kim and Fesenmaier 1990; Lue, Crompton, and Fesenmaier 1993). Nevertheless, elements such as preferences, decision-making procedures, perceptions, travel group composition, and cultural background also exert significant influences on travel sequencing, influences that exist alongside this spatiotemporal effect and interact with it (Bowden 2003; Hwang and Fesenmaier 2003; Hwang, Gertzel, and Fesenmaier 2006; Tideswell and Faulkner 1999). These effects are not entirely captured by the time–geographic framework, and are better considered by other conceptual points of view.
Cognitive and Emotional Dimensions of Time–Space and Tourist Behavior
When time–space conditions create a choice space, conscious decision-making behaviors, as well as cognitive biases and precognitive effects, determine spatiotemporal patterns. The classical model used to describe individuals’ decision-making processes is that of the rational decision maker—an all-knowing, utility-maximizing individual (Chen et al. 2016; Strauss 2008). Stemming from the understanding that this model does not reflect real-life choice-making procedures (Kahneman and Tversky 1979; Simon 1955, 1972; Thaler 1980, 1985), the subfield of behavioral geography has emerged (Golledge 2008; McCormack and Schwanen 2011). Studies in this subfield use disaggregated individual-level data to study actual behaviors, producing more descriptive (instead of normative) models of spatial behavior (Golledge 2008). As such, the field has quickly turned to studying individuals’ cognitive representations of the external world and their accuracy as a critical factor that informs behavioral choices (Couclelis et al. 1987; Downs and Stea 1973; Kitchin 1994; Lynch 1960). Such representations are constructed from knowledge gathered via personal experience or secondhand information (e.g., maps, images, and textual information; Richardson, Montello, and Hegarty 1999). Individuals’ cognitive maps are known to be distorted, thus biasing locational and travel distance/duration estimations (Brunyé, Mahoney, and Taylor 2015; Tversky 1992; biases may vary by source of information—see Richardson, Montello, and Hegarty 1999). The amount of existing knowledge and the extent of its distortion influence behavior, at times leading to incongruence between intentions and actual behavior (Cadwallader 1975).
Yet, some behaviors are more intuitive or even automatic in nature, rather than being the result of well-crafted decision-making processes. In such cases, the cognitive element may be less important and precognitive elements, such as the emotional and the affectual dimensions, may prove to be the critical ones (Pile 2010). Furthermore, as spatial behavior is embedded within a social context, it may be guided by social structures, group decision-making processes, and power relations as well (McCormack and Schwanen 2011). While each of these factors cannot contradict the spatial reality within which behavior is carried, their effects may at times be strong enough to outweigh the constraining role of spatiotemporal conditions. This may be especially true when the choice margin is very large, such as in the case of tourism activities.
Factors from environmental psychology were suggested to be of relevance for tourism studies as early as the 1980s (Fridgen 1984) and indeed, the tourism literature shows some consideration of the roles of these factors in the formation of activity patterns. Tourists are known to develop a cognitive representation of the destination quite quickly, with the nature of this representation depending on familiarity with the destination, individual characteristics, and on the type of navigation aid used (Chang 2015; Fridgen 1987; Guy, Curtis, and Crotts 1990; Pearce 1977; Walmsley and Jenkins 1992). In addition, knowledge differences between first-time and repeat visitors, or between well-prepared and unprepared visitors, were found to affect both cognitive representation (Young 1999) and observed behavioral patterns (Bauder and Freytag 2015; McKercher et al. 2012). Yet, tourism is largely an activity motivated by precognitive and emotional effects that may define the entire tourist experience (Massara and Severino 2013). This is evident in the formation of destination images, relying on both “cognitive” and “affective” conceptions of the destination (San Martín and Rodríguez del Bosque 2008), and also in the final outcomes such as visitors’ satisfaction and destination loyalty (Bigné and Andreu 2004; Yuskel, Yuskel, and Bilim 2010). In terms of intradestination behavior, cultural background has been shown to affect choices such as accommodation, mode of transport, and movement patterns (Dejbakhsh, Arrowsmith, and Jackson 2011) while behavior was claimed to change with activity context (Sarmento 2017). While these findings present new knowledge regarding the formation of behavioral patterns, they are limited in the sense that most of them relate to a subset of the multiple aspects motivating spatial behavior. While relating to all aspects may be beyond the scope of a single article, studying behavior from a more integrative point of view may prove to be beneficial in this context.
Tourist Behaviors: Tracking and Trajectory Processing
Many of the studies mentioned above use high-resolution activity data collected by means such as GPS loggers. This phenomenon is not new, as studies of tourism geographies have been frequently relying on data acquired via advanced tracking technologies for more than a decade now (Shoval and Ahas 2016). These methods replace more traditional observation methods, such as the travel diary and participating observations (Shoval and Isaacson 2007a). These novel methods, which also include passive tracking via cellular networks and social network data, have given rise to a new generation of tourism studies that have produced new insights regarding fundamental questions in the field (Shoval and Ahas 2016). This is not the case only in the field of tourism, and tracking data were already used to produce findings regarding the cognitive (Manley 2016; Manley, Addison, and Cheng 2015), emotional (Shoval, Schvimer, and Tamir 2017, 2018), and social (Slim and Ahas 2016) dimensions of human spatial behavior.
The growing availability of spatiotemporal data calls for the development of approaches for processing and analyzing these data. This need is answered within the realm of trajectory pattern analysis, a field of research focused on the development of methods for analyzing time–space trajectories (the collection of traces depicting the movement of objects during a specific time-period; Nanni and Pedreschi 2006). As sampling conditions are rarely optimal, such data can be incomplete or inaccurate at times. Consequently, raw trajectory data need to be preprocessed, with steps like outlier identification, track smoothing, and map-matching commonly employed (Yan et al. 2013). Furthermore, trajectories in their raw form are rarely useful and most analyses require further processing steps, such as segmenting the data into stops and movements (Renso et al. 2013) or semantically enriching trajectories with information from some other data repositories (Gong et al. 2014; Yan et al. 2013).
Upon the completion of the preprocessing stage, the analysis of trajectories is usually carried with the objective of measuring the similarity between different trajectories and clustering those together, thus identifying general behavioral patterns (i.e., trajectory pattern mining; D’Urso and Massari 2013). This approach was implemented in the context of tourism using social media data (Bermingham and Lee 2014; Yin et al. 2011) and GPS-based trajectories (Shoval and Isaacson 2007b; Tchetchik, Fleischer, and Shoval 2009). Yet, this descriptive approach is not enough to achieve a fundamental understanding of the reasons behind tourists’ observed space–time behaviors. Therefore, while this study uses an approach that generally follows the trajectory pattern mining approach, more analysis steps are carried to expose the intertwined contributions of the cognitive, emotional, and cultural dimensions to the formation of identified patterns.
Methodology
Analysis Approach
The approach utilized in this research follows the one presented by Grinberger, Shoval, and McKercher (2014). This approach, rooted in time–geographic thinking, aims at clustering tourists’ trajectories according to the time–space resource allocation behaviors they reflect. For example, a tendency to cluster stops together would lead to the minimization of travel times, thus reflecting a bias toward maximizing the consumption of the time resource. Yet, a movement path not following the shortest or fastest route between these stops would result in more spatially explorative behaviors, thus indicating a preference towards the space resource.
The approach consists of three steps. First, raw trajectories are transformed into segmented (or structured/syntactic) trajectories via a stop-and-movement identification procedure, a process in which a movement trajectory is broken down into stationary (“stop”) and mobile behavior episodes. Second, two network-based time–space measures are computed for each trajectory: the average distance between pairs of subsequent stops (“average stop distance”; computed using fastest/shortest route distances) and the ratio between the length of the shortest route along the sequence of stops (in their original order) to the distance actually traveled (“movement efficiency”). The first measure indicates the extent to which stops cluster in space while the latter relates to the extent to which chosen paths show travel-time efficiency. Finally, trajectories are clustered according to these measures’ values using the K-means algorithm. These identified clusters represent and characterize general behavioral patterns displayed by tourists within a destination.
This study includes an additional analysis step that studies the effects of different factors on the emergence of identified patterns. This step analyzes the relations between individual-, site-, and activity-specific variables and the emergence of different behavioral patterns by performing a multinomial logistic regression.
Database
The database used in this research includes GPS traces depicting the movement of tourists in Israel, all through their visit, collected as part of a tourist travel survey carried by the Jerusalem Transport Master-Plan Team (JTMT). Tourists were approached by representatives of the JTMT upon arriving at Ben Gurion airport in Israel and were offered to participate in the research in exchange for a smartphone, a communication package offering unlimited local text messaging and 3G surfing, and a $40-worth coupon applicable in the airport’s Duty-Free store (distributed upon the successful completion of the survey). Tourists who agreed to participate filled a preliminary questionnaire, used to screen individuals visiting for more than 15 days or visiting as part of an organized group, and to control the sampling process so as to avoid oversampling specific populations. The questionnaire collected general sociodemographic and tourist profile data, such as country of origin, number of people in the visiting group, visit purposes, and number of previous visits (Table 1 presents the general characteristics of the sample). The surveyor then activated a location tracking application installed on the phone, which was set to sample location every 3 seconds. The application continuously uploaded the data to a main server, where imprecise traces were automatically filtered and trajectories were segmented into stop and movement episodes. Each evening, a representative of the JTMT contacted the tourist by phone and administered a short prompted-recall survey, based on the processed tracking data. As part of the survey, the surveyor asked the participant to indicate which activity (if any) took place during each of the identified stops and the mode of transport used to travel to its location. The participant also provided information regarding stops not captured in the processed trajectory (deu to the participant not carrying the phone with him or her, to a technical malfunction, or to the identification procedure simply not identifying it). In such a case, the surveyor added new stops to the data and labeled them as either “on path” (when they correlated with the tracking data) or “out of path.” The participation ended on returning to the airport, where the participant filled a short satisfaction survey.
Sample Characteristics.
Up to two different purposes were indicated by each participant.
Case Studies
While the spatial extent of the database covers the entire state of Israel (as well as several locations in Jordan and Egypt), a national-level analysis would not make much sense in the current context, as the presence of interdestination movement would bias the analysis. In addition, multiple destinations imply multiple spatiotemporal behavioral contexts, thus making the interpretation of results an extremely complex task. Consequently, the analysis focused on two specific destinations—the cities of Tel-Aviv and Jerusalem (as well as municipalities in their close surroundings, such as Bethlehem for Jerusalem and Herzelia for Tel-Aviv). These cities, while being two of the main tourist destinations in Israel, are very different in both character and spatial structure. Jerusalem is a city sacred to the three largest monotheistic religions, including sites such as the Temple Mount, the Wailing Wall, the Church of the Holy Sepulture in the Old City, and the Mount of Olives, alongside more modern sites such as the Yad Vashem Museum and the Mahne Yehuda Market (see Figure 2). The spatial distribution of attractions in Jerusalem is related to its topographically fragmented structure that includes many hills and valleys. Mobility is further constrained by the cultural context, as participation in activities in the Palestinian neighborhoods of Eastern Jerusalem or the city of Bethlehem may be unavailable to Jewish populations, or simply be (cognitively) unwelcoming.

Selected sites in case study cities. Jerusalem: (1) Ein Karem, (2) Yad Vashem and Mount Herzel, (3) Israel Museum, (4) Mahne Yehuda market, (5) Mea She’arim, (6) City Center, (7) Church of the Holy Sepulcher, (8) Temple Mount, (9) Mount of Olives, (10) Rachel’s Tomb, and (11) Church of Nativity; Tel-Aviv: (12) Jaffa Port, (13) Hatakhana Campus, (14) Neve Tsedek, (15) HaCarmel market, (16) Tel-Aviv Museum of Art, (17) Tel-Aviv Port, and (18) Erets Israel Museum.
Tel-Aviv, on the other hand, is the economic capital of Israel, connected to the global economic system. As such, it offers a variety of cultural, dining, and nightlife activities. None of its attractions is as salient as the main sites in Jerusalem and the dominant element in Tel-Aviv seems to be the beach of the Mediterranean Sea along which the city is located. The area is relatively flat, meaning that topography does not exert a significant effect on mobility and that the city is generally accessible. The beach acts as a linear boundary in relation to which the city is organized. In accordance, tourist attractions in Tel-Aviv almost form a continuum from the old city of Jaffa in the south to the old Port in the north (see Figure 2). Cultural and ethnic tensions are also rare in the city.
The choice of such contradicting pair of destinations was made deliberately. These contrasts allow differentiating between components of the behavioral decision regarding resource allocation related to the specific conditions in one site, to ones transcending the effects of specific time–space behavioral contexts.
Preprocessing Steps
As noted above, GPS data, although being much more precise than data produced via other tracking methods (such as cellular positioning), may still contain many errors due to interferences in reception, sensor malfunction, or misuse. Furthermore, data are valid only at times when the participant carried the logger. Finally, as the focus of this article is on the behavior of individuals, only nonorganized tourism-related patterns were deemed as relevant for the analysis. In accordance with the above, we first carried several preprocessing stages, filtering invalid days from the data set and improving the quality of the data.
These steps (summarized in Table 2) included first a phase of removing outlier samples (traces showing speed and acceleration values that are strikingly abnormal relative to traces recorded before and after them) from the database (see Table 2, step 1, and Appendix A). Then, invalid activity days were removed from the data set of each participant, based on the number of activity hours sampled in each day (step 2, Table 2), the existence of unsampled movement periods (i.e., “gaps”) within the GPS data (step 3), or according to the amount of time the participant flagged as “out-of-path” (step 4). Finally, after more processing steps (steps 5 and 6; see the “analysis specifications” section), two more filtering steps were preformed (see Table 2), based on the types of activities reported in each day (focusing on unorganized tourism activities; step 7), and on the degree to which they match the GPS data (step 8). Matching included finding the most similar GPS-based stops, in terms of location and timing, for each of the reported activities (see Appendix B for the full description of the procedure).
Preprocessing Steps and Derived Data Set Size.
This somewhat strict multistage filtering process, while required for increasing the validity of the analysis, greatly reduced the sample size (see Table 2). The choice of two specific case studies, and not the entire country, further aggravated the situation. To test whether the filtering process produced a biased dataset, the characteristics of the valid subsample were compared with those of the nonvalid sample 1 (see Table 3). While differences were relatively small for most variables, some significant incongruences existed. Some of these—the overrepresentation of “touring” visitors, the underrepresentation of business visitors, and the increase in number of nights—are the natural product of the filtering process (e.g., the longer the visit, the greater the probability that at least one day would be validly sampled). Yet in terms of representation of tourists visiting for religious purposes, of group size, and of the number of visitors traveling with children, the subsample displays some level of bias.
Final Subsample Characteristics and Differences from Nonvalid Sample.
95% confidence level significance.
However, these biases did not seem to affect spatial coverage: a comparison of tourism activity density-surfaces 2 for both the entire dataset (since the nonvalid dataset includes many inaccuracies) and the subsample of valid days shows that the coverage of the main touristic sites and areas did not change by much, despite a general decrease in coverage (Figure 3).

Tourism activity density surfaces for the entire sample and the valid days sub-samples. Tel-Aviv: (a) subsample, (b) entire sample; Jerusalem: (c) subsample, (d) entire sample.
The analysis used only activity days from the final database, which included at least two tourism-related activities in one of the case study cities, resulting in a sample size of 384 activity days (188 for Jerusalem, 196 for Tel-Aviv) representing 285 participants.
Analysis Specifications
As noted above, the analysis procedure included three stages: the segmentation of raw trajectories into stop and movement episodes, the calculation of network-based time–space resource allocation measures, and the clustering of trajectories based on these measures.
Segmentation of trajectories into stops and movement episodes was carried via a stop-identification procedure, which required three parameters: an intrastop maximal distance threshold, beyond which two traces were not considered to be part of the same episode (70 m here); a minimal duration threshold used to decide whether a set of subsequent points that satisfy the distance threshold indeed make a stop or just a momentary pause (300 seconds here); and an interstop minimal distance threshold used to merge subsequent stops that are located nearby (100 m here). The procedure iterated through all samples, and for each sample identified the sequence of subsequent points that satisfied the maximal distance threshold. The algorithm identified this sequence as a stop only if it also satisfied the minimal duration threshold. The location of each stop was set according to the average coordinates over all samples within the sequence (where the weight of each trace in this average was set by its duration, i.e., “duration-weighted” mean). If the distance between centers of two subsequent stops was below the value of the third parameter, the stops were merged into one stop and the mean center was recalculated.
The calculation of resource allocation measures was based only on stops that were matched with a tourism-related activity (each activity could have been matched with more than one stop). The shortest distance between each pair of stops was calculated using Google Maps API, while considering the mode of transport: if during an activity day the participant reported using a private car or a taxi, it was assumed that private vehicle–based transport was available to this individual throughout the day; otherwise only transit transport mode was considered. The distance between each pair of stops was set as the lower one between the shortest walking and the shortest vehicle-based (private or public) route distances. Vehicle-based routes were computed for a mid-week day, six months into the future, using Google’s “best guess” estimation model. This meant that the estimation would rely mostly on Google Map’s historical traffic data, thus avoiding a bias toward momentary conditions in the analysis. The average stop-distance measure was derived by dividing the sum of distances between each pair of stops by the number of pairs minus 1. Movement efficiency was calculated by dividing this sum by the total length of measured movements between stops (aerial distances from the start point to the first movement trace and from the end point to the last trace were added to the length of observed routes). As the two measures use different units (meters per stop for average distance, percentage for movement efficiency), their values were normalized.
Finally, clustering was carried using the K-means clustering algorithm (McQueen 1967), using the relevant command in Stata 13.0. Following the findings of Grinberger, Shoval, and McKercher (2014), K (the number of clusters parameter) was set to three.
Results
Time–Space Resource Allocation
As noted above, the two case studies chosen for the paper are strikingly different, presenting varying cultural and behavioral time–space environments. Therefore, the first issue that required consideration is whether these differences are large enough to affect the distribution of resource allocation behaviors within the population. Graphically (Figure 4), it seems that is not so, as the distribution of the data is quite similar between the two subsamples of visit days. Statistically, however, one difference was found as the movement efficiency values in Jerusalem (M = 0.65, SD = 0.02) were significantly lower than in Tel-Aviv (M = 0.73, SD = 0.02; t = −2.64; p = 0.009). Average distance between stops did not present such a significant effect (estimated using log values, Jerusalem: M = 2.92, SD = 0.03; Tel-Aviv: M = 2.94, SD = 0.03; t = −0.38, p = 0.70). In light of these partial differences, we chose to cluster the data as one sample and later test for differences.

Distribution of time–space resource allocation measures’ values.
Following Grinberger, Shoval, and McKercher (2014), clusters were characterized by comparing cluster-level and dataset–level average values. Each cluster was represented by a two-letter combination, the first letter representing relations between cluster-level and dataset’s movement efficiency average values, and the second relating to the average value of average distances. Cases in which the cluster’s average was higher than the dataset’s were represented by the letter H (standing for higher than average), and L was used otherwise. For example, the combination HL indicated a cluster that presented, in relation to the average values over the entire dataset, a higher average efficiency ratio along with a lower average value of average distances. The three identified clusters (Figure 5A) represented distinct behavioral strategies, with the HL cluster representing a disposition toward maximizing activity time over space, and the other two reflecting substitution effects: in the HH cluster, tourists consume more space by increasing distances between stops, while in the LL cluster the same is achieved via inefficient movement routes. Interestingly, these results echoed Grinberger, Shoval, and McKercher’s (2014) findings regarding the behavior of tourists in Hong Kong (Figure 5b). Differences between the two cases were evident in the share of each cluster out of the total sample—while in Hong Kong the largest cluster was the LL group (with an average size of 38.93, 3 57.25% of the sample), in the current case clusters were more evenly distributed, with HL being the largest group (169 days, 44.59%). In both cases, HH was identified as the smallest group (76 days, 20.05%, in Israel; 10.8, 15.88%, in Hong Kong).

(A) Clustering results for Tel-Aviv and Jerusalem. (B) Clustering results for Hong Kong (source: Grinberger, Shoval, and McKercher 2014, Taylor and Francis Ltd: www-tandfonline-com.web.bisu.edu.cn).
Analysis of Allocation Patterns
The larger sample size and wealth of information allowed an in-depth exploration of the factors that lead to the emergence of the identified clusters, an analysis not carried for the Hong Kong case study. Thus, several variables that could affect the choice of time–space pattern were selected for the analysis (Table 4). These variables were used in a multinomial regression model studying their effects on cluster membership. While this analysis did not rely on the common conceptualization of such models as a method representing utility-maximizing choices, we do borrow here some of its notation in order to examine the role of each variable. The multinomial logistic model relies on a definition of a linear function of the form Y = βx+α. Traditionally, a differentiation between individual-specific and alternative-specific variables is made. Following the different perspectives on elements motivating tourists’ behavior, we suggest here a different, two-level, classification: according to the behavioral dimension (external elements, time–space, cognitive, emotional–cultural) and to the relevant temporal unit (visit day, entire visit). Thus, it possible to identify eight types of variables, the variables used in the current analysis representing six of those (Table 4).
Analysis Variables, by Source of Information and Type.
The multinomial regression used the time-maximizing behavioral pattern (HL) as the baseline category, thus allowing us to explore which variables induced a more space-oriented behavioral pattern. Since religion data was not available for every participant, thus reducing sample sizes, two models were estimated—one not using this variable and one with it, so as to estimate the robustness of effects. The main interesting result from these models (Table 5, models 1 and 2) is that the specific destination has a consistent effect on behavioral patterns. A second stage of analysis was thus carried, where the two models were estimated for each city separately (with and without religion; Jerusalem: models 3 and 4, Tel-Aviv: 5 and 6) (see Table 5).
Multinomial Regression Analysis Results (Unstandardized Coefficients).
Note: Models (1) and (2): entire sample; models (3) and (4): activity days in Jerusalem; models (5) and (6): activity days in Tel-Aviv. Significant values are in bold.
p < 0.10, bp < 0.05, cp < 0.01.
Interpretation of Results
Several interesting findings emerge in Table 5. First, the only variable that remains significant regardless of the type of model is transport mode, where private vehicle–based transport seems to encourage the use of the more spatially explorative behavioral strategies. Other than that, all other effects that were identified in models 1 and 2 either proved to be destination-specific or became insignificant. Some of the effects make intuitive sense:
The effect of transport mode: it is possible to attribute the tendency away from the HL pattern, evident in the positive coefficient values over all models and groups, to the enhanced flexibility that tourists using private vehicle–based transport enjoy, thus adopting less time-focused strategies.
The effects of elderly members of the visiting group: the negative coefficient values for the HH pattern in Jerusalem (models 3 and 4), which represent these groups’ tendency to avoid that pattern in Jerusalem, could be related to the topographical structure of the city and its effects on mobility, along with the more dispersed distribution of activities.
The effects of visiting friends and relatives: it is possible to attribute the tendency toward more focused visits (HL) in Jerusalem, visible in the significantly negative coefficient values for the LL pattern in models 3 and 4, to these tourists having a local guiding them through the city.
The “furthest” hotel effect: allocating time to movement between destinations reduces the amount of time available for activity, thus encouraging tourists to adopt more time-efficient strategies. In Tel-Aviv, where sites are more adjacent to each other and mobility between them is easier, this leads to a significant decrease in the implementation of the more spatially constrained LL strategy, as evident in the negative coefficient values in models 5 and 6. In Jerusalem, where movement between sites is carried over longer distances, the coefficients for the HH pattern are significantly negative (models 3 and 4), meaning its application becomes less probable.
Weekend effect in Tel-Aviv: in Israel, many venues close down during the Sabbath (the time period between Friday afternoon and Saturday evening). While in Tel-Aviv this is less common, main sites like the Carmel and Flea markets, stores in Neve Tzedek, or museums are closed or limit their opening hours during this time. This breaks the north–south continuum of sites in which distances between sites are relatively small, thus driving tourists toward a more dispersed HH pattern (seen in the positive coefficient values in models 5 and 6 for this behavioral pattern).
In contrast to these effects, the interpretation of other effects requires a more in-depth analysis of the behavior of subpopulations. For instance, the differentiating effect of previous visits between the cities, where the coefficients in Jerusalem are almost always positive and significant while in Tel-Aviv no significant effect is registered, can be better understood only when the population is divided by the number of previous visits, and each subpopulation’s spatial patterns are explored (Figure 6A–D). Previous research has shown that repeat visitors present more concentrated visit patterns than first-time visitors (McKercher et al. 2012). This is the picture that rises in Jerusalem where the activities of repeat visitors are concentrated on the route between the city center and the Old City, while first-time visitors show more dispersed patterns, with the Old City being the sole anchor. In Tel-Aviv, however, no real difference is evident. It is possible to attribute these patterns to the organizing role of the Old City in Jerusalem—being the most dominant “must-see” area of the city—which leads returning visitors in the city to focus on it and on the path leading to it, thus creating HL kind of patterns. In Tel-Aviv, no such dominant sites exist (at least not as salient as the Old City), thus enabling more explorative behaviors even for repeat visitors.

Activity densities for different subpopulations: (A) repeat visitors in Jerusalem, (B) first-time visitors in Jerusalem, (C) repeat visitors in Tel-Aviv, (D) first-time visitors in Tel-Aviv, (E) Jewish (religious and other) visitors in Jerusalem, (F) Christian visitors in Jerusalem.
The saliency of the Old City also explains the positive coefficient values of Jewish religious affiliation for the HH pattern in Jerusalem. While the sites in the Old City are the most attractive for both Jewish and Christian tourists, Jewish visitors are more likely to find several other sites to be important or interesting (e.g., Mea She’arim, Yemin Moshe, the City of David). It is thus not surprising to see that their visit more frequently included these sites (Figure 6E), thus encouraging HH behaviors. As the Old City is the only organizing element for Christian visitors, their visit patterns outside of it were dispersed (Figure 6F), enabling the implementation of other strategies.
The role of the Old City affects not only visit patterns within Jerusalem but also the entire visit to Israel. This is the factor behind a seemingly odd finding—the change in behavioral patterns in Tel-Aviv found for tourists visiting for religious purposes (evident in significant coefficient values for this variable in models 5 and 6). Why should religiously oriented individuals change their behavior in a secular city? A possible explanation lies in the perception of Jerusalem as a more central destination within Israel, consequently attracting more visit days and affecting the amount of time available for visiting other locations, including Tel-Aviv. Religiously oriented visitors may thus decide either to skip Tel-Aviv entirely, or to limit their visiting time in it. The latter group may try to compensate for this lack of time by adopting more space-oriented strategies (HH, LL) in Tel-Aviv. To verify this argument, we carried a multivariate Tobit regression (Huang, Sloan, and Adamache 1987; estimated using the mvtobit model in Stata 13.0; Barslund 2015), allowing us to determine the effect of religious activity on the number of visit days in both Tel-Aviv and Jerusalem. This model was chosen since visit length is a censored variable (lower limit of zero) and since the values of dependent variables are determined simultaneously (increasing the number of days in Jerusalem is expected to decrease the duration of visit in Tel-Aviv and vice versa). To control for mediating influences, additional variables were included in the model. The results (Table 6) confirm the hypothesis, showing that the two variables are negatively correlated (correlation value of −0.68), and that visit for religious purposes negatively affects the length of visit in Tel-Aviv (coefficient value of −1.95). Other results, such as the shorter visit in Jerusalem for first-time visitors, the positive effect of total visit duration, and the effects of religious affiliation also make sense.
Multivariate Tobit Model Estimation Results (Unstandardized Coefficients).
Note: Bold values indicate statistical significance.
p < 0.10, bp < 0.05, cp < 0.01.
Discussion and Conclusion
This article presents an analysis of the factors behind the emergence of tourists’ observed visit patterns in the cities of Jerusalem and Tel-Aviv, Israel. The analysis included characterizing tourists’ daily activity trajectories, derived from high-resolution GPS traces data, according to preferences regarding the allocation of time–space resources they reflect. Trajectories were clustered into three unique behavioral patterns, each reflecting a different reaction to the inherent tension between the consumption of time and space resources. The factors leading to the emergence of these patterns were studied via the analysis of the effect of individual-, day-, and destination-specific variables along the emotional–cultural, time–geographic, and cognitive dimensions.
The results point to the complex interrelationships between these dimensions and the unique spatial conditions at the destination, as most of the identified effects were valid only within one specific spatial context. In fact, only two results were replicated regardless of which subset of the data was used to investigate behavior (Jerusalem, Tel-Aviv, or the entire data set): no LH cluster being identified (which also replicates previous results from Hong Kong; Grinberger, Shoval, and McKercher 2014), and an increased tendency to choose the LL or HH patterns when using private vehicle–based transport. The different clusters represent the participants’ choice patterns in reaction to the constraints of reality, a process that is cognitive in nature. Thus, it is possible to interpret the repeating tendency to avoid the space-maximizing LH pattern as resulting from a cognitive overvaluation of the time resource. In this context, it is also possible to explain the effect of transport mode as another cognitive reaction to space–time constraints: the increase in the size of accessible space derived from the use of private transport seems to encourage a less stressful approach toward the time resource, thus leading to the adoption of more space-oriented approaches (LL or HH) at the expense of the most time-efficient one (HL). Thus, these primary effects seem to reflect more generic cognitive biases than responses to the time–space realities of the choice situation.
All other identified effects displayed time–space dependencies related to the morphology of the destination and to the distribution of activities within it. Yet, except maybe for the “furthest hotel” effect (in which reduction in time available for visiting the city is interacted with spatial structure), these effects seem to be very much dependent on the specific cultural, social, cognitive, and precognitive decision-making contexts. For instance, the destinations’ morphology may be used to explain the findings regarding the effects of number of previous visits to Israel, religious affiliation, and day in the week (weekend or not) on behavior. According to Lynch’s seminal findings (1960) on how individuals perceive their urban environments, the distribution of activities is sufficient for creating a cognitive effect on behavior: in Jerusalem, the Old City may act as an organizing “district,” while in Tel-Aviv the north–south continuum parallel to the coast can be viewed as a linear “path.” Yet, as behavior varied in relation to the religious affiliation variable, an element that is cultural/emotional in nature, understanding this finding as representing a purely spatiotemporal–cognitive effect is incorrect. The precognitive (emotional) or metacognitive (cultural) salience, or lack of it, of different sites is the one that probably guided choices in this case.
This dimension seemed to operate at even deeper levels, as the purpose of visit determined the amount of attention (i.e., visit duration) given to each city. This cultural/affectual effect leads to a cognitive one in Tel-Aviv where the scarcity of the time resource (i.e., shorter visit durations) alters behavioral decisions. These decisions also seem to be socially distributed, as the company of the individual affects behavior, as evident in the effects of visiting groups with old-age members or for tourists visiting friends and family. The variance of effects in terms of significance and sign between the cities again suggests that these effects are dependent upon spatial context.
The main contribution of these findings lies in them uncovering the complexity of the procedure through which visit patterns emerge. The roles of different decision-making dimensions change with scale, place, and individual, while the time–space dimension is intertwined through all decision-making processes, when they are carried in situ. Accordingly, this dimension always interacts with other dimensions, where cognitive and precognitive dimensions both exert an in situ effect and guide choice stages that precede the actual activity (such as the effect of visit for religious purposes on visit length in each city). This contribution, beyond adding to the theoretical knowledge, also holds practical implications for design and marketing efforts. For example, it answers the need for a better understanding of the behavioral basis of tourism, as identified within the literature discussing the application of information technologies to assist and guide tourists during their visit (Gretzel 2011; Gretzel, Hwang, and Fesenmaier 2012). This discussion stresses the decision-making context as a critical factor in the design of mobile environments for tourists (Lamsfus et al. 2011, 2012). The results here show that context is an always-evolving concept rather than a static one, thus suggesting that context-aware information applications for tourism need to be designed with flexibility in mind, similar to the view suggested by Lamsfus and his colleagues (2015).
Therefore, it is not advisable to think of decision making as a linear, one-dimensional sequence of events that consists of distinct stages at which each behavioral component has a unique influence. Instead, as proposed by McCormack and Schwanen (2011), the decision is better described as a multidimensional space–time in which complex relationships between different behavioral components exist over multiple spatial scales and temporal instances, in ways that “complicate a neat linearity between past, present, and future” (McCormack and Schwanen 2011, p. 2808). The results here support this view, as well as some of the propositions these authors make regarding the nature of decision based on this conceptualization: the effects of personal characteristics show that the decision is differentiated and experienced differently by different individuals; the effects of group members and social–cultural components point to the decision as socially distributed; and the dependence of many of the effects on the time–space realities (physical and others) exposes the decision to be contextually related. As such, the results here unravel the space–times through which tourists’ behaviors unfold, thus enriching the conceptual and empirical base for the study of intra- and interdestination time–space behaviors.
A second contribution of this article is of a more methodological or conceptual nature. Tourism and spatial behavior studies increasingly use high-resolution mobility in recent years (Shoval and Ahas 2016). While such data is indeed rich in the spatiotemporal sense, it generally lacks when other behavioral dimensions are considered. Recently, several studies have shown how such data can be used to gain knowledge on individuals’ and populations’ cognitive images of space (Manley 2016; Manley, Addison, and Cheng 2015). Such works join previous tourism-related findings on the effect of time–space elements or knowledge on behavior (McKercher et al. 2012; Shoval et al. 2011). The analysis carried out here presents an example of how the intertwined effects of these and other dimensions can be exposed using large mobility data sets. It thus answers a need in this age of data abundance (Chen et al. 2016), suggesting that even without extensive exploration of each individual’s internal world, it is still possible to gather behavioral insights. This is not to say that such studies of the individual are redundant. In fact, we agree here with Chen and colleagues (2016) and suggest that future research aiming to expand our understanding of behavioral processes must rely on the development of large-scale analysis techniques along with in-depth investigations of individuals’ motives. Such joint efforts would serve to indicate how time–spaces, personal characteristics, and social–cultural contexts intermingle to generate spatiotemporal behavior patterns. The new knowledge gained in such a way would not only enrich theoretical thinking but would also lead the way for the design of better management tools for urban tourism.
Footnotes
Appendix A
Appendix B
Acknowledgements
The authors would like to thank the Jerusalem Transport Master Plan Team for providing the data, and to Mr. Muli Yehezkel from MPayer for his support.
Declaration of Conflicting Interests
The author(s) declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article.
Funding
The author(s) received no financial support for the research, authorship, and/or publication of this article.
