Variability in Human Mobility during the Third Wave of COVID-19 in Japan

: Understanding factors regarding individual variability in human mobility during the COVID-19 pandemic would help inform future political and medical decisions. Particularly, understanding environmental factors would provide effective evidence for future urban development and the construction of a mobility service. Given the limited related evidence thus far, we investigated inter-individual variability in human mobility regarding life space during COVID-19 in Japan. We conducted an online survey to assess human mobility, demographics, and personality, and city structure surveys for objective environmental factors. Human mobility was assessed using the Life Space Assessment (LSA) questionnaire. Result revealed that the total LSA score was signiﬁcantly higher in people who were men, middle-aged, working, living with their children, public transportation users, bicycle users, and car drivers, and those having a higher score for extraversion and ego resiliency and a lower score for conscientiousness. People living in the city with high trafﬁc on the roads between plains and mountains had a higher LSA score, and had a lower score where there were many plains with a consistent road density, revealing diverse individual and environmental factors associated with human mobility during the pandemic. Thus, political decisions for urban development should consider these characteristics, the pandemic, and individual convenience.


Introduction
A stretch of life space and sustainable transportation may contribute to the development of socio-economic activities, thus improving individuals' quality of life (QoL). However, after the outbreak of the coronavirus disease 2019 (COVID- 19), socio-economic activities involving outings were restricted based on political and medical decisions to reduce infection spread in Japan and other countries. Such restrictions on outings and human mobility can consequently vary a person's QoL and subjective happiness. Indeed, the lockdown imposed due to COVID-19 has resulted in negative consequences for individuals' mental health and happiness [1][2][3][4][5].
Reduced human mobility would have a large inter-individual variability depending on individual environments, such as the residential area, occupation, family structure, etc., as well as personality, in line with results of studies before the COVID-19 pandemic, suggesting that human mobility was affected by these factors [6]. The variability might also lead to a difference in the risk of infection and other diseases, including mental illness and disabilities due to reduced life space and physical activity. Thus, understanding factors regarding individual variability in human mobility during the COVID-19 pandemic would help inform future political and medical decisions; specifically, to seek and control high-risk populations vulnerable to infection and other diseases and support those requiring help.
Further, understanding both the types of environmental factors that unnecessarily increase human mobility and those that reduce human mobility and compensate for negative individual characteristics, would provide effective evidence for future urban development and the construction of a mobility service. However, evidence pertaining to these factors has been limited, at least in Japan. Hanibuchi et al. (2021) [7] reported that specific individual and social characteristics such as younger people, women, unemployed individuals, those using public transportation, those with chronic diseases, individuals living in metropolitan areas, and those living with more than two people were associated with spending less time outside their homes. However, the sample size at different time points for each analysis was not abundant, and there was a lack of information in terms of the impact on a range of environments in which the individuals move within, and not only on the time spent outside their homes and/or the frequency of outings. Furthermore, the previous study has been limited in investigation of the neighborhood city structure.
Therefore, this study aimed to investigate the factors of inter-individual variability in human mobility in each range of life space environments during the COVID-19 pandemic in Japan. In addition, we conducted a cross-sectional analysis of the association between human mobility and subjective happiness.

Smart Mobility Challenge
This study was conducted as part of the nationwide project, "Smart Mobility Challenge." The Smart Mobility Challenge, which is led by the Ministry of Economy, Trade and Industry and the Ministry of Land, Infrastructure, Transport and Tourism, and initiated in 2019, aims at the social implementation of novel mobility services using the Internet of Things and artificial intelligence.

Participants
The study included samples from both online and city structure surveys. The full online dataset contained data from 13,000 participants aged 18-89 years in 15 areas (cities) of Japan. The sample sizes in each area were constructed to be equal to the proportion of the difference in the area-wise population. However, the number of participants with both datasets (i.e., online and city structure surveys) was 11,886 out of 13,000. We then sampled the participants who were COVID-19 negative, with valid data on relevant variables for the study. In addition, the participants using walking aids and/or needing assistance from caregivers to move, as well as participants who did not leave their bedroom every day, were excluded to select healthy participants who could move freely, maintaining the reliability of the online data. The final dataset contained data from 7472 participants (Supplemental Figure S1).

Online Survey
The online survey was conducted by a subcontractor (VLC HOLDINGS CO., LTD., Tokyo, Japan) in the last week of January 2021 in 16 of the 52 areas where the Smart Mobility Challenge was implemented to understand the population's characteristics in those areas and the association between the convenience of transportation and human mobility, and subjective measures, including subjective satisfaction, QoL, and happiness. The subcontractor provided contracts to the survey participants, which included those relating to privacy protection, monetary compensation, and the consent to participate. We then obtained anonymous data under the contracts with the subcontractor. The online survey was conducted under an ethical manner of the human factors' ethical committee of the National Institute of Advanced Industrial Science and Technology.
The frequency of human mobility in each area of life space was assessed using the Japanese version of the Life Space Assessment (LSA) that was originally developed to assess the risk of long-term care in older adults [8]. For the LSA questionnaire, areas of Sustainability 2021, 13, 13131 3 of 17 human mobility ranged across five different levels ( Figure 1): leaving a bedroom, an outing outside the house, an outing within the neighborhood (0-800 m from the house), an outing outside the neighborhood (0.8-16 km from the house), and an outing outside the town (≥16 km from the house). In this study, we used the frequency of human mobility in each area (1-4 points: 1, less than 1 time/week; 2, 1-3 times/week; 3, 4-6 times/week; 4, everyday/week), as well as the total scores (the LSAtotal), without calculating the effect of walking aids and/or caregivers' assistance (the LSAtotal ranges between 0-60). Thus, the LSAtotal scores reflect the extent of human mobility and its frequency. each area (1-4 points: 1, less than 1 time/week; 2, 1-3 times/week; 3, 4-6 times/week; 4, everyday/week), as well as the total scores (the LSAtotal), without calculating the effect of walking aids and/or caregivers' assistance (the LSAtotal ranges between 0-60). Thus, the LSAtotal scores reflect the extent of human mobility and its frequency.
Personality was evaluated using the Japanese version of the Ten Item Personality Inventory (TIPI-J) [9,10] and the Ego Resiliency Scale (ER89) [11,12]. Anxiety about being infected with COVID-19 was assessed with the following question: How anxious do you feel about getting infected with COVID-19? Participants responded to this question using a 5-point Likert scale ranging from 1 (do not feel at all) to 5 (feel extremely). Subjective happiness was assessed using two tools: a Japanese version of the Subjective Happiness Scale (SHS) [13] that was originally developed by Lyubomirsky and Lepper in 1999 [14] and a sub-scale (general happiness, GH) of a Japanese version of the Subjective Well-being Inventory (SUBI) [15] which was originally developed by Nagpal and Sell in 1985 [16]. The SHS questionnaire contains four items scored on a 4-point Likert scale (ranged 1-4 points) [13]. The average SHS score of the individuals was calculated.

City Structure Survey
The city structure survey was conducted by a subcontractor (Nippon Koei Co., Ltd., Tokyo, Japan) in January 2021 in 15 areas, including 13 areas where the online-based survey was conducted. Statistical data on the living environment, such as the sufficiency rate for daily living services within walking distance, road network structures, and urban types (geography), were calculated using the National Land Numerical Information [17], the Digital Road Map database (Japan Digital Road Map Association, Tokyo, Japan), and traffic volume examinations, as listed in Table 2.

Other Variables
The total number of individuals who were infected with COVID-19 for a month up to the time of the online survey was calculated using prefecture's websites providing COVID-19 information for each relevant prefecture. We also used the variable of whether Personality was evaluated using the Japanese version of the Ten Item Personality Inventory (TIPI-J) [9,10] and the Ego Resiliency Scale (ER89) [11,12]. Anxiety about being infected with COVID-19 was assessed with the following question: How anxious do you feel about getting infected with COVID-19? Participants responded to this question using a 5-point Likert scale ranging from 1 (do not feel at all) to 5 (feel extremely). Subjective happiness was assessed using two tools: a Japanese version of the Subjective Happiness Scale (SHS) [13] that was originally developed by Lyubomirsky and Lepper in 1999 [14] and a sub-scale (general happiness, GH) of a Japanese version of the Subjective Well-being Inventory (SUBI) [15] which was originally developed by Nagpal and Sell in 1985 [16]. The SHS questionnaire contains four items scored on a 4-point Likert scale (ranged 1-4 points) [13]. The average SHS score of the individuals was calculated.

City Structure Survey
The city structure survey was conducted by a subcontractor (Nippon Koei Co., Ltd., Tokyo, Japan) in January 2021 in 15 areas, including 13 areas where the online-based survey was conducted. Statistical data on the living environment, such as the sufficiency rate for daily living services within walking distance, road network structures, and urban types (geography), were calculated using the National Land Numerical Information [17], the Digital Road Map database (Japan Digital Road Map Association, Tokyo, Japan), and traffic volume examinations, as listed in Table 2.

Other Variables
The total number of individuals who were infected with COVID-19 for a month up to the time of the online survey was calculated using prefecture's websites providing COVID-19 information for each relevant prefecture. We also used the variable of whether the areas were in a state of emergency (lockdown). Under the state of emergency, people were prohibited from unnecessary outings and crossing of prefectural boundaries.

Statistical Analyses
Prior to the analyses, all variables were tested for normal distribution based on kurtosis and skewness. Non-normally distributed variables were log-transformed. Gender differences, including gender interactions, were tested using a chi-square test for categorical variables and a t-test for numeric variables. Associations of measures of LSA with potential factors, such as personality, environmental, and demographic factors, were analyzed using general linear models. General linear models were also used to investigate the associations of measures of LSA with SHS and GH, with adjustments for potential covariates. For all models, the inclusion of variables was based on the hypotheses of the present study. All statistical analyses were performed using SAS Enterprise Guide 8.3 (SAS Institute Inc., Cary, NC, USA). The p-value was considered significant at 0.05. Data were presented as mean ± SD, unless otherwise indicated.

Participants' Characteristics and Geographic Characteristics of Cities
The participants' characteristics are displayed in Table 1. Except for marital and car status (shared family cars), there was a significant difference in the mean values of variables between gender and the proportion of gender in the other variables. The LSA measures were higher among men than women (all p < 0.0001), whereas both happiness scales were higher among women than men (both p < 0.0001).
A summary of the geographic characteristics is presented in Table 2. The population of each relevant area ranged from 18,297 to 799,598. The population density of each relevant area ranged from 27.4 people/km 2 to 5996.2 people/km 2 ; in brief, our study included data ranging from a small city in a rural area to a big city in a metropolitan suburban area. During the online survey, three areas were in a state of emergency.  Note. Sufficiency rate for daily living services within walking distance (≤0.8 km), the proportion of people who can access all the daily living services (including medical, welfare, and commercial facilities, and core public transportation routes) by walking; coverage of population within walking distance (≤0.8 km) from medical facilities, the proportion of people who can access medical facilities (a hospital for internal medicine or surgery, a clinic) by walking; coverage of population within walking distance (≤0.8 km) from welfare facilities, the proportion of people who can access welfare facilities by walking; coverage of population within walking distance (≤0.8 km) from commercial facilities, the proportion of people who can access commercial facilities (a supermarket, a department store) by walking; coverage of population within walking distance (≤0.8 km) from core public transportation routes (≥30 available /daily), the proportion of people who can access core public transportation stations (railway or bus with more than 30 operations per day) by walking; insufficiency rate for medical facilities within walking distance (≤0.8 km) from housing occupied by the elderly population (≥65 years), the proportion of housing that is not within walking distance for elderly people (≤0.8 km) from the medical facilities; coverage of the elderly population (≥65 years) within 1 km from elderly welfare facilities, the proportion of elderly people who live within 1 km from the welfare facilities.

Individual Factors of Human Mobility during the Pandemic
The results of the individual factors, such as age, gender, and occupation, of the LSAtotal are shown in Table 3. Regarding the demographic factors, male gender, public transportation use, bicycle riding, and car driving contributed to a higher LSAtotal score (Table 3). Further, people aged 40-59 years had a higher LSAtotal score compared to those aged 18-39 years even after a multiple comparison correction with the Tukey-Kramer method (all p < 0.05, Figure 2A). People who were working had the highest LSAtotal score (all p < 0.05, Figure 2B) and LSA at levels 2-5 among the occupations. People who were living with their children had a higher LSAtotal score (p = 0.0056, Figure 2C) and LSA at levels 2-4 (p < 0.05) regarding the child status even after a multiple comparison correction with the Tukey-Kramer method. People aged 18-29 years had the lowest scores, whereas those aged 70-79 years had the highest LSA scores at levels 2-4. However, at level 5, people aged 50-59 years had the highest LSA score, and there was a significant difference between individuals aged 30-39 years and those aged 50-59 (p = 0.0086).  LSAtotal, the total score of the life space assessment questionnaire. The TIPI-J is a ten-item personality inventory; the ER89 is an ego resiliency questionnaire. The statistics of all these variables were considered to be interrelated using the general linear model. 0.05, Figure 2B) and LSA at levels 2-5 among the occupations. People who were living with their children had a higher LSAtotal score (p = 0.0056, Figure 2C) and LSA at levels 2-4 (p < 0.05) regarding the child status even after a multiple comparison correction with the Tukey-Kramer method. People aged 18-29 years had the lowest scores, whereas those aged 70-79 years had the highest LSA scores at levels 2-4. However, at level 5, people aged 50-59 years had the highest LSA score, and there was a significant difference between individuals aged 30-39 years and those aged 50-59 (p = 0.0086). In addition, at levels 2-4 of LSA, vehicle use, including public transportation use, bicycle riding, and car driving contributed to a higher LSA score (all p < 0.05); however, bicycle riding's contribution disappeared at level 5, although the contribution of the other variables still appeared at level 5. At levels 3-5 of LSA, men had a higher LSA score (all p In addition, at levels 2-4 of LSA, vehicle use, including public transportation use, bicycle riding, and car driving contributed to a higher LSA score (all p < 0.05); however, bicycle riding's contribution disappeared at level 5, although the contribution of the other variables still appeared at level 5. At levels 3-5 of LSA, men had a higher LSA score (all p < 0.05). At level 5 of LSA, owning cars significantly contributed to a higher LSA score (all p < 0.0001).
Regarding personality factors, extraversion and ego resiliency scores were positively associated with scores of LSAtotal and LSA at levels 2 and 5, whereas the conscientiousness score was negatively associated with LSAtotal and LSA at levels 3 and 5 (all p < 0.05). At level 4 of the LSA, only extraversion had a positive effect on the LSA score (p = 0.0013).

Environmental Factors of Human Mobility during the Pandemic
The LSAtotal of the environmental factors, such as the area characteristics for the living environment and urban types, are shown in Table 4. Each factor (i.e., area characteristics for the living environment and urban types) was independently analyzed as solely determining the area characteristics to prevent multicollinearity between variables.
In terms of the area characteristics for a living environment, each sub-variable was also analyzed separately. The sufficiency rate for daily living services within walking distance (Standardized β = −0.045, p < 0.0001), the population coverage within walking distance from welfare facilities (Standardized β = −0.039, p = 0.0003), and the population coverage within walking distance from core public transportation routes (Standardized β = −0.053, p < 0.0001) were negatively associated with the LSAtotal, while the insufficiency rate for medical facilities within walking distance from housing occupied by the elderly population (Standardized β = 0.045, p < 0.0001) was positively associated with the LSAtotal. For the LSA scores at each level, the sufficiency rate for daily living services within walking distance was negatively associated with the LSA scores at levels 4 (Standardized β = −0.027, p = 0.0176) and 5 (Standardized β = −0.057, p < 0.0001), the population coverage within walking distance from welfare facilities was negatively associated with the LSA score at level 5 (Standardized β = −0.060, p < 0.0001), and the population coverage within walking distance from commercial facilities was positively associated with the LSA score at level 5 (Standardized β = 0.024, p = 0.0371). Similarly, the population coverage within walking distance from core public transportation routes was negatively associated with the LSA scores at levels 4 (Standardized β = −0.025, p = 0.0303) and 5 (Standardized β = −0.075, p < 0.0001); the population density in areas around train stations or bus stops for public transportation lines was negatively associated with the LSA scores at levels 2 (Standardized β = −0.025, p = 0.0346), 3 (Standardized β = −0.041, p = 0.0005), and 4 (Standardized β = −0.042, p = 0.0003); and the insufficiency rate for medical facilities within walking distance from housing occupied by the elderly was positively associated with the LSA scores at level 5 (Standardized β = 0.064, p < 0.0001). These variables were adjusted for the variables in Table 3, except for the area-related variables. LSAtotal, the total score of the life space assessment questionnaire.
For the urban types, after multiple comparison correction with the Tukey-Kramer method (Figure 3), people living in the area where there were many plains and where road density was almost consistent had the lowest LSAtotal score, whereas people living in the area where traffic on the roads between plains and mountains was high had the highest LSAtotal score. For the LSA score at level 5, the relationships among urban types were almost consistent with those of the LSAtotal. However, for levels 2 and 3, the LSA score of people living in the area where there were plains and mountains, and where road links were concentrated in the plains, was higher than in the area where hubs with high road density were dispersed.
in the area where traffic on the roads between plains and mountains was high had the highest LSAtotal score. For the LSA score at level 5, the relationships among urban types were almost consistent with those of the LSAtotal. However, for levels 2 and 3, the LSA score of people living in the area where there were plains and mountains, and where road links were concentrated in the plains, was higher than in the area where hubs with high road density were dispersed.

COVID-19-Related Factors of Human Mobility during the Pandemic
Anxiety about getting infected with COVID-19 was associated with a lower LSA score only at level 3 (β = −0.0167, Standardized β = −0.0247, p = 0.0287), after the adjustment of potential covariates, including the variable "areas." Relevant analyses were performed for the independent variables below by replacing them with variable "areas." The state of emergency was associated with higher LSAtotal scores (β = 0.7108, p = 0.0198, Figure 4A) and LSA at level 5 (β = 0.2753, p < 0.0001, Figure  4E), whereas the state of emergency was associated with a lower LSA score at levels 2 (β

COVID-19-Related Factors of Human Mobility during the Pandemic
Anxiety about getting infected with COVID-19 was associated with a lower LSA score only at level 3 (β = −0.0167, Standardized β = −0.0247, p = 0.0287), after the adjustment of potential covariates, including the variable "areas." Relevant analyses were performed for the independent variables below by replacing them with variable "areas." The state of emergency was associated with higher LSAtotal scores (β = 0.7108, p = 0.0198, Figure 4A) and LSA at level 5 (β = 0.2753, p < 0.0001, Figure 4E), whereas the state of emergency was associated with a lower LSA score at levels 2 (β = −0.0354, p = 0.0123, Figure 4B), 3 (β = −0.0918, p < 0.0001, Figure 4C), and 4 (β = −0.0799, p = 0.0054, Figure 4D The rapid increase in new COVID-19 cases in each prefecture where the relevant area was located (i.e., ratio of new COVID-19 cases in the past month to total cases since the reporting of the first case in each prefecture) was associated with a lower LSAtotal score (β = −0.1087, Standardized β = −0.0463, p < 0.0001) and LSA at level 5 (β = −0.0224, Standardized β = −0.0818, p < 0.0001) but was not significantly associated with LSA at levels 2-4.  The rapid increase in new COVID-19 cases in each prefecture where the relevant area was located (i.e., ratio of new COVID-19 cases in the past month to total cases since the reporting of the first case in each prefecture) was associated with a lower LSAtotal score (β = −0.1087, Standardized β = −0.0463, p < 0.0001) and LSA at level 5 (β = −0.0224, Standardized β = −0.0818, p < 0.0001) but was not significantly associated with LSA at levels 2-4.

Associations of Human Mobility with Subjective Happiness
The LSAtotal score was significantly and positively associated with the GH score (β = 0.0012 ± 0.004, r = 0.0376, p = 0.0014; Table 5 Model 1) but not with the SHS score (β = 0.0011 ± 0.006, r = 0.0214, p = 0.0646, Table 5 Model 1) after adjustment for gender, age, household income, occupation, residential area, marital status, child status, and anxiety about getting infected with COVID-19. However, the association between LSAtotal score and GH score disappeared after adjustment for personality, such as the measures of the TIPI-J and ER89. In addition, there were no associations between the SHS or GH scores and any level of LSA scores after adjustment for potential covariates including personality.

Discussion
The COVID-19 pandemic has led to a reduction in human mobility in Japan [18][19][20][21]. The government and other institutions have had to make swift political and medical decisions without established evidence and past experience. Therefore, the present nationwide study investigated the determinants, including personality, demographic, and environmental factors, of the inter-individual difference in human mobility during the pandemic in Japan to provide evidence for the next phase of COVID-19. This study was performed during the third wave (January 2021) of the pandemic in Japan. During this period, people's lifestyles began to change from the acute phase of the pandemic to the next stage, where people have begun to develop a new lifestyle to prevent the COVID-19 infection. For example, almost all leading companies have introduced work-from-home jobs depending on the occupational category. The state of emergency announcements and restrictions of human mobility have been based on several criteria established by the government and/or each local government. Therefore, our results can serve as evidence for a guide to a new lifestyle.
Our study revealed several features that might contribute to individual differences in human mobility. Notably, gender, occupation, vehicle use types, personality, and several types of living environments strongly contributed to the inter-individual variation in the LSAtotal score ( Figure 5A).
In addition, the features likely varied depending on the level of life space ( Figure 5B). To further interpret these associations, it would be necessary to consider that the features varied depending on the level of life space because the LSAtotal tends to be affected by the score of a wider range of life space due to its calculation method.

Effects of Individual Factors on Human Mobility
Once the demographic factors of human mobility during the pandemic are clarified, this should provide significant evidence on population groups requiring help. This study indicated that people who were men, middle-aged, and working had higher scores of LSAtotal and scores at almost all levels of LSA, implying that specific types occupations, such as an administrative position and manual labor, as well as essential workers, might contribute to greater human mobility. This result might suggest that even though workfrom-home was imposed in society, these individuals had to go to their workplace. The impact of gender and occupation on human mobility was consistent with the results during the early phase of the pandemic in Japan, as indicated by Hanibuchi et al. [7].
People living with their children also had higher LSAtotal scores and scores at almost all levels of LSA. This may be because they had to visit the grocery store more often and pick their children up from school as well as from leisure activities. In addition, younger people aged 20-39 years had relatively lower LSA scores. This might indicate that this generation can enjoy indoor life because of the popularization of video games or the Internet.

Effects of Individual Factors on Human Mobility
Once the demographic factors of human mobility during the pandemic are clarified, this should provide significant evidence on population groups requiring help. This study indicated that people who were men, middle-aged, and working had higher scores of LSAtotal and scores at almost all levels of LSA, implying that specific types occupations, such as an administrative position and manual labor, as well as essential workers, might contribute to greater human mobility. This result might suggest that even though workfrom-home was imposed in society, these individuals had to go to their workplace. The impact of gender and occupation on human mobility was consistent with the results during the early phase of the pandemic in Japan, as indicated by Hanibuchi et al. [7].
People living with their children also had higher LSAtotal scores and scores at almost all levels of LSA. This may be because they had to visit the grocery store more often and pick their children up from school as well as from leisure activities. In addition, younger people aged 20-39 years had relatively lower LSA scores. This might indicate that this generation can enjoy indoor life because of the popularization of video games or the Internet.
Not surprisingly, people who used public transportation and rode bicycles had greater LSA scores, although a previous study indicated that there was no association between the use of public transportation and outings [7]. The association between bicycle riding and the LSA score disappearing at level 5 was also understandable. Although, as indicated by these results, vehicle use definitely contributes to the expansion of human mobility, the risk of infection by vehicle use requires further investigation.
As indicated by previous studies that personality factors, such as the Big Five, contribute to inter-individual variability in physical activity [22][23][24], our study revealed that the extraversion score and conscientiousness score was positively and was negatively associated with human mobility, respectively. In addition, our study found that ego resiliency as a factor of personality, which has received attention in the field of mental health, Not surprisingly, people who used public transportation and rode bicycles had greater LSA scores, although a previous study indicated that there was no association between the use of public transportation and outings [7]. The association between bicycle riding and the LSA score disappearing at level 5 was also understandable. Although, as indicated by these results, vehicle use definitely contributes to the expansion of human mobility, the risk of infection by vehicle use requires further investigation.
As indicated by previous studies that personality factors, such as the Big Five, contribute to inter-individual variability in physical activity [22][23][24], our study revealed that the extraversion score and conscientiousness score was positively and was negatively associated with human mobility, respectively. In addition, our study found that ego resiliency as a factor of personality, which has received attention in the field of mental health, also had an association with human mobility. This is the first study to describe the associations between ego resiliency and human mobility, as well as personality and human mobility during the COVID-19 pandemic.

Effects of Environmental Factors on Human Mobility
Our study also revealed several environmental factors affecting human mobility during the pandemic. Notably, the variables of environmental factors around one's home had an impact on human mobility in areas far from the home. For example, insufficient daily living services within walking distance was associated with greater human mobility for the areas of LSA at levels 4-5. In addition, human mobility in the area where traffic on the roads between plains and mountains was high would be greater than in other areas. These results provide useful information for future urban development projects.

Effects of COVID-19 on Human Mobility
This study indicated that the impact of COVID-19 on human mobility depended on the area and the range of life space. Particularly, associations of the state of emergency with human mobility were varied at levels of life space in that, for the areas near the home, a state of emergency might have an impact on human mobility restriction, with no impact for the areas far from home. This might be explained by the fact that the areas under the state of emergency consisted of relatively larger cities. It is expected that people would commute a longer way than those living in areas without a state of emergency.
Although anxiety about being infected with COVID-19 had no impact on human mobility since the rapid increase in new COVID-19 cases was associated with lower human mobility, some subjective feelings might influence the related restriction.

Effects of Human Mobility on Subjective Happiness
Although several studies have pointed out that the pandemic [25,26] and the lockdown associated with it [1][2][3][4][5] have worsened mental health and lowered happiness, our results did not indicate a clear association between human mobility and measures of subjective happiness after adjustment for personality. This implies that subjective happiness reflects a large inter-individual variability in personality and the impact of the pandemic on subjective happiness did not outweigh the impact of personality. The results of this study are understandable because of its cross-sectional analysis approach.
In addition, the impact on mental health was affected by the time period from the start of the pandemic. A previous study indicated that the effect on mental health gradually became apparent [2]. In Japan, it has been reported that the suicide rate was lower during the first phase of the pandemic but then increased [4]. Furthermore, for some people, restricted human mobility might improve subjective happiness by improving relationships among family members and QoL at home.
Therefore, further longitudinal studies are needed to clarify the association between reduced human mobility, related restrictions, and subjective happiness.

Limitation
Several limitations should be considered when interpreting the results of this study. First, because of its cross-sectional nature, we could not identify causality. Second, because our survey included online data, the assessment of human mobility was obtained using a subjective method. This might cause some misevaluation of the associations. Third, although our dataset had a relatively larger sample size from each area, it could not fully avoid sampling bias due to randomized anonymous sampling. This bias might lead to a misunderstanding of the association between environmental factors and human mobility. To support our results, further studies using objective methods for the assessment of human mobility are required.

Summary of Factors That Increase/Decrease Human Mobility by Each Life Space
For human mobility in the areas within the home and neighborhood, men, older adults aged 70-79 years, working individuals, those living with their children, public transportation users, bicycle users, and car drivers, including those who had a higher score of extraversion and ego resiliency and lower score of conscientiousness, and were living in an area where traffic on the roads between plains and mountains was high had a larger amount of human mobility. A higher population density in areas around train stations or bus stops for public transportation lines also had a negative impact on human mobility for areas ranging from the neighborhood to those far from home. In addition, the state of emergency lowered human mobility in the neighborhood. Transportation from one's home to the neighborhood is considered important for pursuing daily activities such as grocery shopping. To restrain the human traffic flow, the presence of environmental factors that could suitably support the individual factors would be ideal, but few meaningful relationships were found in this regard. By dispersing urban structures consisting of hubs with high road density, it may be possible to reduce the human traffic flow. In addition, although the population density in areas around train stations or bus stops for public transportation lines was negatively associated with the LSA scores at level 2, this association might be the result of an inherent fear of the COVID-19 infection and restraint from contact with too many people.
For the areas far from home, men, middle-aged individuals between 50-59 years, working individuals, public transportation users, and car drivers but not bicycle users, and those who had a higher score of extraversion and ego resiliency and a lower score of conscientiousness had a larger amount of human mobility. Furthermore, insufficient daily living services within walking distance and higher population coverage within walking distance from commercial facilities were associated with greater human mobility in areas far from home. In contrast, higher population coverage within walking distance from welfare facilities and core public transportation routes was associated with lower human mobility for areas far from home. In addition, the area where traffic on the roads between plains and mountains was high contributed to a larger amount of human mobility, whereas the area where there were many plains and road density was almost consistent was associated with lower human mobility for the areas far from home.

Conclusions
This study identified several factors that are determinants of human mobility during the COVID-19 pandemic in Japan. With regard to individual factors, including personality, our study revealed that gender, occupation, child status, vehicle use types, and specific personality types would be strong determinants of inter-individual differences. In particular, multi-aspect support is needed for people living with children during the pandemic. With respect to environmental factors, our study revealed that there is a lack of daily living services, including medical, welfare, and commercial facilities, and core public transportation routes, which increases human mobility in areas far from home. Thus, it is important to complement systems and/or services, such as home delivery services, to reduce the extent of human mobility in rural and suburban areas during the pandemic, in addition to urban reconstruction. However, even when these systems are developed, workers responsible for logistics and delivery are constantly needed in real time. It will also be necessary to investigate the construction of a system and urban design that will enable the minimum amount of human mobility and optimal logistics. Together, we suggest the necessity of making political decisions based on the characteristics of each area and individual, while also considering resilience measures against the pandemic and individual convenience for urban development.  Institutional Review Board Statement: Ethical review and approval were waived for this study, due to the anonymous online survey by the external institution.

Informed Consent Statement: Not applicable.
Data Availability Statement: The data are not publicly available; however, the data that support the findings of this study may be available from the corresponding author, T.A., upon reasonable request.