June 30, 2026
The rapid advance of smart cities increasingly depends on trajectory data mining, yet underrepresented demographic groups, particularly the elderly, are often sparsely represented in public mobility datasets. This underrepresentation can introduce systematic bias into mobility modeling and downstream urban planning. Using the 2016–2020 Jersey City subset of the Citi Bike System Data, this study quantitatively examines how the absence of underrepresented subgroups’ mobility signatures affects mobility modeling, using synthetic trajectory generation as a case study. The analysis reveals that elderly riders exhibit a structurally distinct mobility signature, including localized activity spaces (958 m vs.,189 m for young riders), lower mobility entropy (1.82 vs.), and asymmetric off-peak temporal patterns. To demonstrate that relying on majority-dominated training data yields biased synthetic outcomes, we further evaluate both a first-order Markov chain and a Qwen3-4B model fine-tuned with QLoRA across three demographic training settings: the full population, young riders only, and elderly riders only. Results show that models trained on majority-dominated populations systematically misrepresent elderly mobility behavior, particularly for spatial mobility metrics. The Markov model trained on the full population overestimates elderly step length by 4.5% and dwell time by 8.9%, whereas the elderly-specific model achieves substantially lower errors across most metrics. Comparisons between the Markov and LLM-based frameworks further show that higher-capability models do not necessarily improve subgroup-level fidelity under limited demographic data. These findings underscore the importance of demographic representation in mobility modeling and its downstream applications for underrepresented populations.
The integration of smart technologies into urban infrastructure has reshaped transportation planning and public health, with trajectory data mining at the core of this transformation [1]. As global populations age, ensuring that smart-city systems serve all demographic groups equitably has become a pressing concern in social computing [2]. Yet the datasets underpinning these systems rarely carry granular demographic annotations, rendering the mobility patterns of older adults largely invisible to standard analytical pipelines [3]–[5]. Recent large-scale analyses confirm systematic demographic disparities in human mobility data [6], [7], and a growing body of work calls for algorithmic fairness in urban computing [8], [9].
In age-friendly urban planning [10], capturing the distinct travel patterns of older adults is critical: their mobility is shaped by physical accessibility, safety concerns, and different daily activity schedules [11], [12]. Mobility models trained on broad, unannotated datasets that are disproportionately representative of younger, working-age populations inherently encode demographic biases [13].
When synthetic mobility data derived from majority populations informs downstream decisions such as facility allocation or public transit routing, the mobility patterns of elderly riders may become systematically underrepresented, potentially undermining initiatives such as the 15-minute city vision of localized and accessible neighborhoods [14].
A critical barrier to studying this problem is the scarcity of mobility datasets with reliable demographic labels. The Citi Bike bike-sharing system historically collected riders’ birth year and gender at account registration, providing one of the few large-scale mobility sources with age metadata [15]. However, citing privacy concerns, Citi Bike removed demographic attributes from public trip data beginning in 2021.
To examine how demographic underrepresentation propagates bias into downstream mobility modeling, this work uses synthetic trajectory generation as a controlled testbed. We construct an experimental framework comparing a first-order Markov chain and a Large Language Model(LLM) under identical demographic training conditions. We focus on the 2016-2020 JC Citi Bike dataset and construct a controlled experimental design that isolates the effect of demographic composition. The main contributions are as follows:
Age-Specific Mobility Characterization: Through comparative analysis of elderly (age \(\ge 65\)) and young (age 18–35) rider groups in the 2016–2020 Jersey City Citi Bike, we show that the elderly exhibit a structurally distinct mobility signature: spatially confined (radius of gyration of 958 m versus 1,189 m), behaviorally routine (entropy 1.82 versus 4.15), and temporally offset from commuting peaks.
Controlled Multi-Model Experiment: We evaluate both a first-order Markov chain and a Qwen3-4B model fine-tuned with QLoRA under three demographic training settings (full population, young-only, and elderly-only). This controlled design isolates the effect of demographic composition on downstream mobility fidelity while enabling direct comparison between structured probabilistic and higher-capacity sequence modeling approaches.
Effects of Demographic Composition on Downstream Mobility Fidelity: Our results quantify how demographic imbalance in training data propagates bias into downstream mobility modeling. Models trained on majority-dominated populations systematically overestimate elderly mobility characteristics, while elderly-specific training achieves substantially closer agreement with empirical elderly mobility patterns.
Reproducible Evaluation Framework: We propose a reproducible evaluation pipeline spanning data preprocessing, mobility metric construction, synthetic trajectory generation, and subgroup-specific evaluation. The framework provides a benchmark for evaluating demographic-specific fidelity in downstream mobility modeling.
Understanding human mobility is foundational to urban computing [1]. Early seminal work by González et al. [16] and Song et al. [17] utilized mobile phone data to demonstrate that human movement is highly regular and predictable, exhibiting bounded spatial variance. Alessandretti et al. [18] further revealed that human mobility operates across characteristic spatial and temporal scales, with individual trajectories exhibiting stable, idiosyncratic patterns.
Building upon these foundations, trajectory generation has evolved from statistical modeling to sophisticated deep learning techniques. Traditional approaches, such as grid-based Markov models, simulate movement via transition matrices, while data-driven routine models synthesize trajectories by learning historical motifs [19]. Recently, deep generative models, including Generative Adversarial Networks (GANs) [20], [21], sequence-to-sequence architectures, and spatial-temporal graph neural networks [22], have been employed to capture complex spatial-temporal dependencies. The emergence of diffusion probabilistic models has further advanced the field: DiffTraj [23] and related methods [24] leverage iterative denoising to generate high-fidelity GPS trajectories with improved sample quality and diversity. Most recently, LLM-based approaches such as MobilityGPT [25] have been proposed, treating trajectories as token sequences for autoregressive generation. Several comprehensive surveys have catalogued this rapidly expanding landscape [26], [27]. However, the majority of these models are optimized for general population accuracy, often overlooking the nuanced mobility patterns of underrepresented demographic subgroups [13].
Evaluation of synthetic trajectory quality presents a multifaceted challenge. Classical approaches rely on statistical distribution matching (e.g., comparing step length, speed, radius of gyration), while sequence alignment metrics such as Dynamic Time Warping [28] and Fréchet distance [29] assess structural fidelity. Recent work by Lucas et al. [30] proposed a comprehensive quality framework spanning statistical, spatial, and temporal dimensions, while Wang et al. [31] systematically reviewed trajectory similarity metrics in the era of deep learning.
Alongside quality concerns, privacy-preserving trajectory generation has received growing attention. Cunningham et al. [32] established foundational methods for synthetic location data with formal privacy guarantees. Chen et al. [33] advanced differentially private trajectory generation that balances utility and privacy, while Van de Ven et al. [34] critically examined the limitations of generative models for trajectory privacy. Due to privacy concerns, demographic attributes such as birth year and gender were no longer publicly available in Citi Bike trip data beginning in 2018, leaving fewer demographic signals available for evaluating algorithmic fairness on this platform. Our work quantifies how the absence of such signals can reduce the fidelity of mobility models and their downstream outputs for older adults.
The intersection of social computing and urban gerontology highlights the necessity of age-friendly urban environments [35]. Research indicates that older adults often face mobility constraints related to physical capabilities and access to public transportation, contributing to a strong reliance on localized, neighborhood-centric travel [3], [36]. Recent large-scale studies leveraging mobile phone data [11] and GPS tracking [12] have provided quantitative evidence that older adults exhibit significantly smaller activity spaces, lower mobility entropy, and distinct temporal patterns compared to younger demographics.
Despite the growing emphasis on inclusive cities [10] and the 15-minute city paradigm [14], privacy concern [4] and the limited availability of age annotations in public mobility datasets have constrained data-driven research on older adults’ mobility. Recent analyses of demographic disparities in mobility data [6], [7], together with emerging work on algorithmic fairness in urban computing [8], underscore the importance of addressing representation gaps. Our work bridges this gap by leveraging the 2016–2020 Jersey City subset of the Citi Bike System Data, which contains publicly available age-annotated trip records, to systematically evaluate how training-data composition affects the fidelity of mobility models for underrepresented populations.
Before evaluating how demographic composition biases downstream mobility modeling, we first establish whether elderly and young riders exhibit meaningfully different mobility patterns in the first place. This section processes the Citi Bike trip records into mobility trajectories, defines the spatial and temporal metrics used throughout the paper, and applies them to compare the two groups. The resulting characterization confirms that elderly riders follow a structurally distinct mobility signature, setting the stage for the controlled experiments in Section 4.
Citi Bike is a bike-sharing system operating in Jersey City and New York City. It provides trip-level records containing origin and destination stations, station coordinates, trip duration, user type, and, for records prior to 2021, rider birth year and gender. Using rider birth year, we classify 1,638,153 Jersey City trips recorded between 2016 and 2020 into age groups over a 60-month period.
After removing trips with Haversine distance greater than 10 km or computed speed exceeding 10 m/s, 93.3% of records were retained, yielding a full-population dataset containing 1,528,325 trips from 6,102 bikes. Using rider birth year, we define elderly (65+) and young (18–35) subsets from the filtered data. The elderly subset contains 14,394 trips from 2,007 unique bikes, representing 0.9% of all filtered trips, whereas the young subset contains 783,558 trips from 5,832 bikes, representing 51.3%. Trips associated with the same bike are ordered chronologically to construct mobility trajectories and compute dwell times.
Table 1 summarizes the post-filtering descriptive statistics of the JC Citi Bike 2016–2020 dataset stratified by age group.
| Group | Trips | Trajectories | Trips/Traj. | Unique Stations | Age |
|---|---|---|---|---|---|
| Elderly (65+) | 14,394 | 2,007 | 7.2 | 50 | 65–80 |
| Young (18–35) | 783,558 | 5,832 | 134.4 | 55 | 18–35 |
| All Riders | 1,528,325 | 6,102 | 250.5 | 55 | 16–80 |
The elderly group is substantially smaller than the young group, accounting for 0.9% of all filtered trips compared to 51.3% for young riders. Elderly riders also exhibit shorter trajectories on average (7.2 vs. trips per trajectory) and interact with fewer unique stations (50 vs.), indicating lower overall mobility intensity and spatial coverage.
Gender and subscription patterns are broadly similar across the two groups. Elderly riders are composed of 70% male and 28% female users (2% unspecified), compared to 75% male and 25% female among young riders. Both groups consist almost entirely of annual Subscribers (99%).
We evaluate synthetic mobility fidelity using spatial and temporal metrics that characterize mobility range, visitation regularity, and travel dynamics. Let a trajectory consist of \(n\) trips for a given bike, where each trip \(i\) contains origin and destination coordinates together with departure and arrival timestamps.
Step Length: The Haversine distance between trip origin and destination stations, averaged across all trips within a group.
Radius of Gyration (RoG): A measure of spatial dispersion computed relative to the trajectory center of mass \((\overline{\text{lat}}, \overline{\text{lon}})\): \[r_g = \sqrt{\frac{1}{n} \sum_{i=1}^{n} d\big((\text{lat}_i^e,\text{lon}_i^e),\;(\overline{\text{lat}}, \overline{\text{lon}})\big)^2}\]
Mobility Entropy: The Shannon entropy of station visitation frequencies: \[S = - \sum_{j} p_j \log_2(p_j)\] where \(p_j\) denotes the fraction of trips ending at station \(j\). Lower entropy indicates more regular and predictable mobility patterns.
Speed: The ratio of trip distance to trip duration.
Dwell Time: The idle interval \(t_{i+1}^s - t_i^e\) between two consecutive trips on the same bike. Dwell times are restricted to \([0,86400)\) seconds to exclude overnight or long-term inactivity.
Intra-day Temporal Distribution: The normalized hourly distribution of trip start times aggregated over 24-hour bins.
Figure 1 summarizes demographic differences in mobility structure and temporal activity patterns between elderly and young riders.
Spatial Mobility Differences. Figure 1 (a) compares mobility metrics between elderly and young riders. Step lengths and travel speeds remain broadly similar across the two groups (1,004 m vs.,102 m; 2.41 m/s vs. m/s). In contrast, elderly riders exhibit lower spatial dispersion and mobility diversity, with a smaller radius of gyration (958 m vs.,189 m) and substantially lower mobility entropy (1.82 vs.).
The aggregate all-rider entropy (4.36) exceeds that of the young group alone, indicating that elderly riders contribute additional station visitation diversity beyond the dominant young-rider mobility pattern. This difference is further reflected in the station-set Jaccard similarity between the two groups, which is 0.521.
Temporal Activity Patterns. Figure 1 (b) compares intra-day mobility activity between the two groups. Young riders exhibit pronounced morning and evening activity peaks around 08:00 and 17:00–18:00, whereas elderly riders maintain more evenly distributed daytime activity concentrated between 09:00 and 17:00. Elderly riders also exhibit shorter dwell intervals between trips (18,329 s vs.,928 s), indicating distinct temporal usage patterns between the two groups.
Section 3 established that elderly and young riders follow structurally distinct mobility patterns. These differences carry direct implications for downstream urban computing tasks, such as facility allocation and transit routing, that depend on faithful mobility representations. A foundational component shared by these tasks is trajectory generation: synthetic mobility data are used to augment sparse observations, power urban simulations, and enable privacy-preserving data sharing. As LLMs emerge as the emerging paradigm for trajectory-based mobility analysis, understanding whether demographic bias propagates through generative pipelines becomes increasingly urgent. This section examines the question through a controlled experiment in which a first-order Markov chain and a Qwen3-4B model (QLoRA fine-tuned) are each trained under elderly-only, young-only, and full-population settings, and evaluated against real elderly trajectories using the metrics defined in Section 3.3.
To evaluate how training-data demographics affect synthetic elderly mobility fidelity, we construct a controlled generative experiment using three training populations: elderly-only (65+), young-only (18–35), and the full rider population. For both the Markov and LLM paradigms, separate models are trained on each population:
\(\mathcal{M}_{\text{elderly}}\) / \(\mathcal{L}_{\text{elderly}}\): Trained exclusively on elderly trips (65+).
\(\mathcal{M}_{\text{young}}\) / \(\mathcal{L}_{\text{young}}\): Trained exclusively on young trips (18–35).
\(\mathcal{M}_{\text{all}}\) / \(\mathcal{L}_{\text{all}}\): Trained on the full filtered dataset (all age groups).
Each model generates 2,007 synthetic trajectories to match the empirical elderly trajectory count. Synthetic trajectory lengths are sampled from the real elderly trajectory-length distribution (mean 7.2 trips per trajectory). All generated datasets are evaluated against the real elderly trajectories using the identical mobility evaluation pipeline.
Comparisons across training populations assess how demographic composition affects the fidelity of synthetic elderly mobility trajectories, while comparisons between the Markov and LLM paradigms examine the role of generative framework choice under severe minority-data sparsity.
We implement a first-order Markov chain model to generate synthetic mobility trajectories from empirical trip sequences. The model estimates station-to-station transition probabilities \(P(s_j \mid s_i)\), start-station probabilities \(P_{\text{start}}(s_i)\), and empirical distributions for dwell time, departure hour, and trajectory length. Trip durations for each station pair \((i,j)\) are modeled using Gaussian distributions parameterized as \(\mathcal{N}(\mu_{ij}, \sigma_{ij})\).
Synthetic trajectories are generated by first sampling a start station and trajectory length, followed by iterative sampling of destination stations conditioned on the current station. Each transition is assigned a duration sampled from the corresponding transition-specific Gaussian distribution, truncated to a minimum of 60 seconds. Dwell intervals between consecutive trips are sampled from the empirical dwell-time distribution and capped at 24 hours, while departure hours are sampled from the empirical hour-of-day distribution with randomized minutes and seconds. For stations without outgoing transitions in the training data, destination sampling falls back to the start-station distribution \(P_{\text{start}}(s_i)\). Synthetic datasets are generated by repeating this procedure until the target trajectory count is reached.
We also implement a Qwen3-4B [37] LLM fine-tuned using 4-bit NormalFloat QLoRA [38] to evaluate a higher-capacity generative framework under sparse demographic data. LoRA adapters are applied to all linear layers with rank \(r=16\), \(\alpha=32\), and dropout 0.05, yielding approximately 100M trainable parameters.
Each trajectory is represented as a whitespace-separated sequence of station IDs, such as 3195 3205 3186 …. Trajectories are tokenized using the Qwen3 tokenizer and truncated to 768 tokens. Separate LoRA adapters are fine-tuned for 3 epochs
using the elderly-only, young-only, and full-population training corpora described in Section 4.1, with batch size 4 and learning rate \(2\times10^{-4}\). Training corpus
sizes are summarized in Table 2.
| Training Corpus | Trajectories | Trips |
|---|---|---|
| Elderly (65+) | 1,631 | 13,928 |
| Young (18–35) | 2,897 | 469,801 |
| All Riders | 2,984 | 667,013 |
The fine-tuned model generates trajectories via batched autoregressive decoding (batch size 32) with temperature 0.7, top-\(p\) 0.9, and top-\(k\) 50, producing up to 128 new tokens per trajectory. Decoding is seeded with a start-of-sequence token. The generated token stream is parsed back into station ID integers by concatenating consecutive digit tokens; non-numeric tokens and station IDs not present in the training vocabulary are discarded. Each model generates 2,007 synthetic trajectories to match the real elderly trajectory count.
Since the LLM operates on station sequences only and does not model continuous temporal quantities, we synthesize temporal attributes post hoc using fixed heuristics. Trip duration is set to trip distance divided by 4.17 m/s (the real elderly mean speed), and inter-trip dwell time is sampled from an exponential distribution with mean 600 s. These heuristics are identical across all three LLM-trained models.
We evaluate how training-data composition affects the fidelity of synthetic elderly mobility trajectories using three Markov models trained on different demographic populations: \(\mathcal{M}_{\text{full}}\), \(\mathcal{M}_{\text{young}}\), and \(\mathcal{M}_{\text{elderly}}\). Each model generates 2,007 synthetic trajectories matching the empirical elderly trajectory count and length distribution. Figure 2 and Table 3 summarize the resulting mobility metrics relative to the real elderly baseline.
| Metric | Real Elderly | \(\mathcal{M}_{\text{full}}\) | \(\mathcal{M}_{\text{young}}\) | \(\mathcal{M}_{\text{elderly}}\) |
|---|---|---|---|---|
| Step Length (m) | 1,004.3 | 1,049.5 (+4.5%) | 1,067.7 (+6.3%) | 984.2 (\(-\)2.0%) |
| Speed (m/s) | 2.41 | 2.00 (\(-\)17.0%) | 2.27 (\(-\)5.9%) | 2.56 (+6.2%) |
| Dwell Time (s) | 19,338 | 21,064 (+8.9%) | 27,893 (+44.3%) | 18,392 (\(-\)4.9%) |
| RoG (m) | 957.8 | 1,101.1 (+15.0%) | 1,106.9 (+15.6%) | 738.8 (\(-\)22.9%) |
| Entropy | 1.82 | 4.41 (+141.9%) | 4.21 (+130.8%) | 1.79 (\(-\)1.8%) |
\(\mathcal{M}_{\text{young}}\) exhibits substantial deviations from the real elderly baseline, particularly for Dwell Time (+44.3%), Entropy (+130.8%), and RoG (+15.6%), indicating that mobility patterns learned exclusively from young riders do not transfer well to elderly mobility behavior.
Including all demographic groups in the training data improves temporal fidelity. Compared with \(\mathcal{M}_{\text{young}}\), \(\mathcal{M}_{\text{full}}\) reduces the Dwell Time error from 44.3% to 8.9%. However, spatial metrics remain poorly reproduced, with Entropy and RoG errors remaining comparable to the young-only model. This suggests that aggregate transition patterns remain dominated by the majority population.
Among the three variants, \(\mathcal{M}_{\text{elderly}}\) achieves the closest agreement with the real elderly baseline for Step Length, Dwell Time, and Entropy, with errors below 5%. However, it underestimates RoG by 22.9%, indicating reduced spatial coverage in the generated trajectories. This contraction likely reflects the sparsity of elderly transition observations in the training data.
Despite the substantially larger training corpus used in the full-population setting, the spatial fidelity gap remains unresolved. Together, these results indicate that demographic inclusion alone is insufficient to recover minority-group spatial mobility structure when observational coverage is sparse.
We repeat the synthetic trajectory experiment using Qwen3-4B with QLoRA fine-tuning. Three adapters (\(\mathcal{L}_{\text{elderly}}\), \(\mathcal{L}_{\text{young}}\), \(\mathcal{L}_{\text{all}}\)) are trained on the same demographic splits used in the Markov experiments, with each model generating 2,007 synthetic trajectories. Table 4 summarizes the evaluation against the real elderly baseline.
| Metric | Real Elderly | \(\mathcal{L}_{\text{elderly}}\) | \(\mathcal{L}_{\text{young}}\) | \(\mathcal{L}_{\text{all}}\) |
|---|---|---|---|---|
| Step Length (m) | 1,004.3 | 1,189.2 (+18.4%) | 967.3 (\(-\)3.7%) | 878.8 (\(-\)12.5%) |
| Speed (m/s) | 2.41 | 4.15 (+72.3%) | 4.15 (+72.2%) | 4.16 (+72.4%) |
| Dwell Time (s) | 19,338 | 597 (\(-\)96.9%) | 601 (\(-\)96.9%) | 605 (\(-\)96.9%) |
| RoG (m) | 957.8 | 903.6 (\(-\)5.7%) | 846.9 (\(-\)11.6%) | 773.5 (\(-\)19.2%) |
| Entropy | 1.82 | 2.72 (+49.1%) | 3.06 (+67.6%) | 2.86 (+57.0%) |
Across the three LLM variants, \(\mathcal{L}_{\text{elderly}}\) achieves the lowest errors for the spatial metrics RoG and Entropy. However, all LLM variants exhibit substantial errors for the temporal metrics Speed and Dwell Time. Because temporal attributes are synthesized post hoc rather than directly generated by the model, temporal patterns remain poorly aligned with the empirical elderly baseline.
Compared with the Markov models, the LLM-based models generally achieve lower fidelity across most metrics. \(\mathcal{L}_{\text{elderly}}\) underperforms \(\mathcal{M}_{\text{elderly}}\) on Step Length, Speed, Dwell Time, and Entropy. The primary exception is RoG, where the LLM achieves a smaller error (\(-\)5.7%) than the Markov model (\(-\)22.9%). This suggests that autoregressive sequence generation partially alleviates the spatial contraction caused by sparse transition observations. However, this improvement in spatial coverage comes at the cost of reduced behavioral regularity, as the LLM models tend to overestimate mobility diversity, particularly for Entropy.
Overall, the results suggest that higher-capacity sequence models do not necessarily improve minority-group mobility fidelity when demographic training data remain limited.
This study examined how demographic underrepresentation biases mobility models, using synthetic trajectory generation as a controlled testbed. Elderly and young riders exhibit structurally distinct mobility patterns across spatial and temporal dimensions. Controlled experiments show that models trained on majority-dominated populations systematically misrepresent elderly mobility behavior, while elderly-specific training achieves closer agreement with empirical patterns, although spatial coverage remains constrained by data sparsity. Comparisons between Markov and LLM-based generation further reveal that higher-capacity models do not necessarily improve subgroup-level fidelity under limited demographic data.
The results suggest that aggregate mobility metrics may mask substantial behavioral differences across demographic groups. Models trained primarily on majority populations can misrepresent minority mobility behavior even when overall performance appears acceptable. Incorporating demographic-specific evaluation into mobility analysis may therefore improve the reliability of data-driven transportation planning and policy decisions.
The empirical mobility patterns observed among elderly riders also indicate that bike-sharing systems may serve older populations differently from younger users. Elderly riders exhibited more localized spatial activity, lower mobility diversity, and distinct temporal usage patterns, suggesting that age-inclusive mobility infrastructure and improved station accessibility may help broaden participation among older adults.
More broadly, as data-driven mobility models become increasingly integrated into urban analytics and smart-city systems, evaluation based solely on aggregate performance metrics may overlook systematic demographic-specific errors. Whether the downstream task is generation, prediction, or resource allocation, demographic-aware validation procedures can help improve the robustness and representativeness of mobility models used in urban decision-making.
This study is subject to several limitations. First, the analysis is restricted to the Jersey City Citi Bike system, and future work should evaluate whether similar subgroup-specific mobility patterns emerge across larger and more diverse urban settings. Second, demographic analysis is limited to age group. Extending subgroup-level mobility evaluation to additional demographic dimensions and downstream tasks remains an important direction for future research. Finally, this study captures mobility patterns within a specific historical period (2016–2020), while urban mobility systems and demographic travel behaviors may evolve gradually over longer time horizons. Future research should examine whether subgroup-specific mobility disparities persist under long-term changes in urban mobility systems.
To appear in 7th International Conference on Social Computing (ICSC)↩︎