title: Multiscale Dynamic Human Mobility Flow Dataset in the U.S. during the COVID-19 Epidemic
authors: Yuhao Kang Song Gao Yunlei Liang Mingxiao Li Jinmeng Rao Jake Kruse
text: 
  	
  	
  	
  	
  	
  	
  	
  	Multiscale dynamic human mobility flow dataset in the U.S. during the COVID-19 epidemic
      Yuhao Kang1, Song Gao1*, Yunlei Liang1, Mingxiao Li1,2,3, Jinmeng Rao1, Jake Kruse1
      August, 2020
  =======================================================================================
  
  fancy
  	1. GeoDS Lab, Department of Geography, University of Wisconsin-Madison, WI 53706, United States; 2. Institute of Geographic Sciences and Natural Resources Research, Chinese Academy of Sciences, Beijing 100101, China; 3. School of Architecture and Urban Planning, Shenzhen University, Shenzhen 518061, China
  
  	*Corresponding Author: Song Gao (song.gao@wisc.edu)
  	
  	
  	 
  		
  		[A preprint version. Final corrected proof is available on Scientific Data: <https://www.nature.com/articles/s41597-020-00734-5>]Understanding dynamic human mobility changes and spatial interaction patterns at different geographic scales is crucial for assessing the impacts of non-pharmaceutical interventions (such as stay-at-home orders) during the COVID-19 pandemic. 
  		In this data descriptor, we introduce a regularly-updated multiscale dynamic human mobility flow dataset across the United States, with data starting from March 1st, 2020.
  		By analysing millions of anonymous mobile phone users’ visits to various places provided by SafeGraph, the daily and weekly dynamic origin-to-destination (O-D) population flows are computed, aggregated, and inferred at three geographic scales: census tract, county, and state.
  		There is high correlation between our mobility flow dataset and openly available data sources, which shows the reliability of the produced data. 
  		Such a high spatiotemporal resolution human mobility flow dataset at different geographic scales over time may help monitor epidemic spreading dynamics, inform public health policy, and deepen our understanding of human behaviour changes under the unprecedented public health crisis. 
  		This up-to-date O-D flow open data can support many other social sensing and transportation applications.
  	
  	
  	
  
  § BACKGROUND & SUMMARY
  
  	
  	The outbreak of the novel coronavirus disease SARS-CoV-2 (also known as COVID-19) in December 2019 has become a global threat to public health and human societies.
  	Thus far, more than 25 million people have been infected by the virus with more than eight hundred thousand death cases globally <cit.>.
  	To contain the transmission of the COVID-19, social distancing has been proved as the most effective non-pharmaceutical intervention <cit.>, and governments have applied various policies to reduce human mobility and restrict large gatherings, such as regional lockdowns <cit.>, stay-at-home orders <cit.>, and travel restrictions <cit.>.
  	Tracking dynamic human mobility changes and spatial interaction patterns is therefore a prerequisite for measuring the effects of human mobility and interventions on predicting the virus spread <cit.>. Several recent works have employed human movement flow matrices in understanding spatial interaction changes and social impact, and enabling network-based epidemic models to project the numbers of COVID-19 infected population in different countries such as China, Japan, Italy, France, Chile, and UK <cit.>, which requires up-to-date inbound and outbound human movement flow information. 
  	However, there is no such openly and timely updated human movement origin-to-destination (O-D) flow matrix data at a fine spatiotemporal resolution available in many other countries where researchers can only use historical O-D survey data and other proxies as a compromise <cit.>.
  	
  	Human mobility has been widely studied in multiple disciplines such as geography, transportation, urban planning, physics, computer sciences, and public health <cit.>. It reflects patterns about how people move from place to place and serves as an indicator of human behaviour and underlying socioeconomic environments.
  	With the rapid development of information and communication technologies (ICT) and GPS embedded devices, large-scale mobile phone data provides an unprecedented opportunity in tracking human trajectories, which benefits research about human mobility patterns.
  	Existing studies have used such data to investigate basic laws governing human movements <cit.>, model regional transportation connectedness and economy <cit.>, describe daily commuting flows <cit.>, compute urban vibrancy <cit.>, inform public health policy <cit.>, and understand spatial interaction patterns <cit.>.
  	
  	As pointed out by several recent studies <cit.>, human mobility data plays a key role and serves as data foundation in the fight against the COVID-19 pandemic.
  	Although companies such as Descartes Labs (<https://www.descarteslabs.com/mobility/>) <cit.>, Apple (<https://covid19.apple.com/mobility>), Google (<https://www.google.com/covid19/mobility/>) <cit.>, and Facebook (<https://dataforgood.fb.com/docs/covid19/>) <cit.>, have released a set of near real-time mobility-related open datasets for monitoring human mobility changes and social distancing behaviour during the COVID-19 period, these datasets are lacking in three respects.
  	First, human mobility flow matrices, which describe movement patterns from origin geographic units to destination regions, are often unavailable, even though such O-D paired flow matrices are incredibly valuable for epidemic transmission modeling and spatial-social interaction measurements <cit.>. For instance, only aggregated mobility indices (such as median travel distance and foot-traffic) at a specific region is provided by the Descartes Labs, Apple, and Google mobility datasets.
  	Second, there is a lack of fine-resolution datasets regarding the privacy-accuracy trade off. 
  	Most openly available datasets aggregate mobility patterns to the state, county or city scale, while higher spatial resolution (e.g., census tract) datasets, which provide more detailed human mobility patterns, are necessary to more accurately characterize heterogeneous human mobility changes within cities or intra-counties during the pandemic <cit.>.
  	Third, the mobile phone or other sensor data-driven mobility patterns often only provide a sample (e.g., 10%) of the entire population. There is no entire-population-level estimated flow data.
  	
  	
  	To address the limitations of existing mobility databases, we introduce an openly available dataset that provides an estimation of dynamic population flows at multiple spatial scales (at the census tract, county, and state) and temporal resolutions (daily and weekly) across the U.S. during the COVID-19 pandemic considering the findability, accessibility, interoperability, and reusability of data <cit.>.
  	A similar mobility dataset has been released by researchers from Italy <cit.>, while we focus on the mobility in the U.S.. 
  	The O-D format dataset is generated by tracking millions of anonymous mobile phone users' trajectories collected by SafeGraph <cit.>.
  	When producing the dataset, great efforts were taken to protect personal privacy by aggregating to various geographic scales so that individual information cannot be traced.
  	Other public datasets, such as the the American Community Survey (ACS) commuting flows and the Descartes Lab COVID-19 mobility dataset are then compared to illustrate the reliability of the produced dataset.
  	Such an up-to-date dataset can be a useful supplement in human mobility observation. It can be used not only in the fight against the COVID-19 pandemic, but also to benefit other researches and applications such as emergency response <cit.>, urban planning <cit.>, and population migration <cit.>.
  	
  	
  	
  	
  
  § METHODS
  
  	Figure <ref> illustrates detailed processing steps for how this human mobility flow dataset is generated.
  	The dynamic population O-D flow matrices are estimated using mobile phone location data provided by SafeGraph and demographic data retrieved from the ACS.
  	Based on millions of anonymous mobile phone user visits to various places tracked by SafeGraph, two types of visitor flows, namely daily census block group (CBG) to CBG visitors and weekly CBG to point of interest (POI) visitors are computed, respectively.
  	After spatially joining the place visitors to the administration regions, the visitor O-D flows are computed at three different spatial scales: census tract, county, and state, which are used to provide multi-scale views of human mobility and spatial interactions patterns between different places.
  	Since the number of mobile phone users detected by SafeGraph is about 10% sample of the entire population <cit.>, we further employ the ACS population data with mobile phone data samples to infer the population level of dynamic O-D flows.
  	
  	
  
   §.§ Track Place Visits
  
  	The place visitor patterns are retrieved from the SafeGraph COVID-19 Data Consortium <cit.>.
  	Millions of anonymous GPS pings collected from numerous mobile applications are tracked and then cleaned to remove noise. 
  	Then, users' home places are estimated and aggregated (e.g., at the level of a CBG), and those users' visits from home places to POIs are tracked.
  	POIs are the primary venue for tracking place foot-traffic by SafeGraph, while CBG is one of the fine-resolution geographical units the United States Census Bureau used for publishing demographic and socioeconomic data.
  	A home place of a user refers to his/her most common nighttime location during the last six weeks. 
  	For each day, GPS pings of each device are clustered and only those clusters during nighttime hours (6pm - 7am local time) are kept.
  	The CBG with the most clusters in that day is recorded.
  	Based on this, the most frequent CBG over the last six weeks that reflects the primary nighttime location is used as the “home location” for each user.
  	By aggregating home places to CBGs, user privacy can be protected as no individual records can be traced and accessed <cit.>.
  	
  	Active users' visits to POIs are produced with a similar strategy.
  	Using several clustering methods such as density-based spatial clustering for applications with noise (DBSCAN)<cit.>, GPS pings are grouped together in which each cluster contains a set of potential POIs and associates with CBGs.
  	The best place for a given cluster is classified by performing machine learning methods involving several entangled features.
  	Thereby each user's visits from home place to various POIs and CBGs are identified <cit.>.
  	
  	In total, there are more than 5 million POIs stored in the database, as well as more than 220 thousand CBGs retrieved from the ACS.
  	The spatial density distribution of POIs across the Contiguous United States is mapped in Figure <ref>, which shows that places cluster in major cities and are generally located along streets.
  	The more places, the brighter the region in the map.
  	
  	
  
   §.§ Compute Visitor Flows
  
  	Two major human mobility flow metrics are employed in data production, and are denoted as daily CBG to CBG visitor flows and weekly CBG to POI visitor flows.
  	In the daily CBG to CBG visitor flows metric, each row contains an origin CBG and a destination CBG, as well as the number of mobile phone-based visitor flows from the origin CBG to the destination CBG.
  	Every day, the number of unique mobile phone users who live in the origin CBG and visits to the destination CBG are recorded.
  	More specifically, GPS pings of each user are clustered first.
  	Only those clusters (i.e., not a single trajectory point) with a duration of at least one minute are counted as a “visit” <cit.>.
  	By doing so, the daily mobile phone-based visitor flows between CBG and CBG are grouped and summed up.
  	A sample record of the daily CBG to CBG visitor flows is included in Table <ref>.
  	
  	For the weekly CBG to POI visitor flows metric, different from the daily CBG to CBG visitor flows metric which aggregates visitors between origin CBG and destination CBG directly, it provides a mapping of CBGs to POIs.
  	In other words, the number of unique visitors who live inside the origin CBG and visit the destination POI in one week are counted.
  	A sample record of the weekly CBG to POI visitor flows is also included in Table <ref>.
  	
  	
  
   §.§ Multiscale Aggregation
  
  	The two mobile phone-based visitor flows metrics (from CBG to CBG and from CBG to POI) are both processed at the CBG scale. 
  	After obtaining these two metrics, all data are further aggregated into three different spatial scales: census tract, county, and state.
  	The motivations for providing data products at multiple spatial scales are discussed as follows.
  	First, using coarser/finer analysis unit may lead to different outputs, which is known as the scale effect.
  	As an important and fundamental concept in geography, the scale effect exists in almost all geographic phenomena <cit.>, including human mobility patterns <cit.>.
  	Providing a multiscale flow dataset allows us to have a more comprehensive view of human mobility and spatial interaction patterns.
  	Second, various research projects may acquire data at different spatial scales, depending on the usage.
  	For example, for research focusing on a macro view of spatial interactions, state or county scale data might be more suitable as it reflects general regional mobility patterns, while census tract scale data can be used for describing a micro view of human movement patterns such as within cities.
  	Also, considering the data size-accuracy trade-off (i.e., the higher the spatial resolution, the higher the accuracy but the larger the data size is), providing a multiscale dataset enables users to download the data that fit their own needs.
  	Third, although the daily CBG to CBG visitor flows and weekly CBG to POI visitor flows have been computed at the CBG scale, which is a finer spatial scale, aggregating them to coarser-level can preserve the data privacy better.
  	Therefore, the O-D flow dataset is generated at three geographical scales, respectively.
  	To do so, we assign census tract, county, and state's geographically unique identifier to each origin CBG and destination CBG that they belong to, and group all records according to the O-D pairs.
  	In total, there are 74,001 census tracts, 3,219 counties, and 50 states & Washington D.C. & Puerto Rico at each spatial scale for aggregation in the U.S.
  	
  	The aggregated mobile phone-based weekly CBG to POI visitor flows and the daily CBG to CBG visitor flows are termed as weekly visitor flows and daily visitor flows hereafter, respectively.
  	
  	
  
   §.§ Infer Dynamic Population Flows
  
  	The above-mentioned visitor flows at the three spatial scales are based on mobile phone users detected by SafeGraph, not on the entire population.
  	These users account for about 10% to the entire population in the U.S., and the sampling ratio of unique mobile devices to population varies across CBGs to CBGs <cit.>.
  	Existing studies have shown that a good representative sample of the entire population can reflect general human mobility patterns <cit.>. Also, short-term mobility and long-term migration are affected by multiple factors such as physical travel costs (distance, time, and money), economic and health benefits, social and political frictions  <cit.>. Given that SafeGraph' samples are highly correlated with the true Census populations regarding several socio-economic attributes <cit.>, we aim to infer the short-term population-level dynamic mobility flows during the COVID-19 pandemic as it is important to accurately estimate meta-population infection cases and in other human mobility applications <cit.>.
  	To do so, by utilizing the official ACS population data with mobile phone visitor patterns, the dynamic population flows are inferred using the following equation:
  
  	
      pop_flows(o, d) = visitor_flows(o, d) ×pop(o)/num_devices(o)
  
  	where pop_flows(o, d) is the estimated dynamic population flows from geographic unit o to geographic unit d, visitor_flows(o, d) is the computed mobile phone-based visitor flow from o to d, pop(o) indicates the population at the geographic unit o extracted from the ACS,  and num_devices(o) refers to the number of unique mobile devices residing in o. In addition, we also compare the estimation results with a gravity model and a radiation model (see more details in the “Technical Validation” section).
  	
  	
  
  § DATA RECORDS
  
  	We have produced two data products: weekly flow data and daily flow data, both of which are provided at the census tract, county, and state scales, starting from March 1st, 2020.
  	A static copy of the dataset has been uploaded on Figshare <cit.>, while the live version of the dataset will be kept up-to-date with new data stream and can be downloaded from an open data repository on Github: <https://github.com/GeoDS/COVID19USFlows>.
  	Data provided in this repository are separated into two folders daily_flows and weekly_flows to store daily flow data and weekly flow data, respectively.
  	The two folders are organized according to the geographic scale, where ct2ct indicates flows between census tract to census tract, county2county refers to flows between county to county, and state2state contains flow data that originate from one state to others.
  	All files are stored in a comma-separated values (CSV) format, which has been widely used for storing, transferring, and sharing data publicly.
  	File names are formatted as {data_type}_{spatial_scale}_{date}.csv, e.g., weekly_ct2ct_03_02.csv and daily_state2state_04_19.csv.
  	Specifically, for weekly flow data, the dates in file name refers to the date of the Monday in that week but summarize all mobility flows in that week from Monday to Sunday.
  	Since the file size of flow data at the census tract scale exceeds the GitHub disk limit, each flow data file is split into 20 subfiles (that can be merged after downloading).
  	The daily and weekly aggregated directional flows are coded by pairwise origin to destination units using geo-identifiers (GEOIDs). 
  	The GEOIDs are numeric codes that uniquely identify different administrative levels (e.g., census tract, county, and state) in the U.S. Census Bureau data portal <cit.>. 
  	For the GEOIDs at each scale, census tract is using an 11-digit number, county is a 5-digit number, and state is a 2-digit number.
  	The coordinates of origins and destinations can be used for creating spatial interaction flow maps. 
  	Reference shapefiles at each scale are available from TIGER/Line Shapefiles (<https://www.census.gov/geographies/mapping-files/time-series/geo/tiger-line-file.html>).
  	External demographic and socioeconomic statistical information at different spatial scales can be accessed directly from the U.S. Census Bureau and joined to each origin and destination using GEOIDs (<https://www.census.gov/data.html>).  
  	A description of all attributes in the dataset is shown below.
  	
  	
  	
  
    §.§.§ Weekly Flow Data
  
  	geoid_o - Unique identifier of the origin geographic unit (census tract, county, and state). Type: string. 
  
  	geoid_d - Unique identifier of the destination geographic unit (census tract, county, and state). Type: string.
  
  	lat_o - Latitude of the geometric centroid of the origin unit. Type: float.
  
  	lng_o - Longitude of the geometric centroid of the origin unit. Type: float.
  
  	lat_d - Latitude of the geometric centroid of the destination unit. Type: float.
  
  	lng_d - Longitude of the geometric centroid of the destination unit. Type: float.
  
  	date_range - Date range of the records. Type: string.
  
  	visitor_flows - Estimated number of visitors detected by SafeGraph between the two geographic units (from geoid_o to geoid_d), computed and aggregated from weeekly CBG to POI flows. Type: float.
  
  	pop_flows - Estimated entire population flows between the two geographic units (from geoid_o to geoid_d), inferred from visitor_flows as described in the subsection `Infer Dynamic Population Flows'. Type: float.
  
  	
  	
  	
  
    §.§.§ Daily Flow Data
  
  	geoid_o -  Unique identifier of the origin geographic unit (census tract, county, and state). Type: string. 
  
  	geoid_d - Unique identifier of the destination geographic unit (census tract, county, and state). Type: string.  
  
  	lat_o - Latitude of the geometric centroid of the origin unit. Type: float.
  
  	lng_o - Longitude of the geometric centroid of the origin unit. Type: float.
  
  	lat_d - Latitude of the geometric centroid of the destination unit. Type: float.
  
  	lng_d - Longitude of the geometric centroid of the destination unit. Type: float.
  
  	date - Date of the records. Type: string.
  
  	visitor_flows - Estimated number of visitors between the two geographic units (from geoid_o to geoid_d), computed and aggregated from daily CBG to CBG flows. Type: float.
  
  	pop_flows - Estimated entire population flows between the two geographic units (from geoid_o to geoid_d), inferred from visitor_flows as described in the subsection 'Infer Dynamic Population Flows'. Type: float.
  
  	
  	
  	
  
  § TECHNICAL VALIDATION
  
  	To check the data distribution and ensure the reliability of the produced mobility O-D flow dataset, three complementary methods are employed for data validation.
  	We first compare the probability distributions of the computed visitor flows and the estimated entire population flows to check if the data distributions stay consistent during the production process.
  	Then, we compare the estimated population flows using our inference approach with results estimated from a gravity model and a radiation model. Lastly, the released mobility flow dataset is compared with two other openly available data sources: ACS commuting flows and the Descartes Labs mobility changes.
  	The hypothesis is that mobility patterns at the same geographic scale should be consistent across multiple data sources.
  	
  	
  
   §.§ Checking Distributions of Visitor Flows and Population Flows
  
  	Following the methods described above, the two types (daily and weekly) of mobile phone-based visitor flows are directly aggregated at different geographic scales, while the demographic data are involved in estimating the entire population flows.
  	Visitor flows and population flows are supposed to have similar distributions, and thus to make sure the inferring process keeps the distributions of the mobility flows unchanged, we made the Q–Q (quantile-quantile) plots to compare their distributions.
  	Visitor flows and population flows are first normalized, as they have different value ranges.
  	If the distributions of the two metrics being compared are linearly related, the scatter points should locate following a line y = kx, where x,y are percentiles of the two datasets and k is the coefficient.
  	Figure <ref> shows the Q-Q plots of the visitor flows and population flows at three spatial scales (a: census tract, b: county, and c: state) based on the weekly flow data in the week of 03/02/2020 to 03/08/2020.
  	Though we only plot the two distributions using one week data as an example, the associations between visitor flows and population flows are similar for other dates.
  	As is shown in Figure <ref>, scatter points are distributed along a straight line at both county scale and state scale. Even though the line is not y = x, the inferred entire population flows are linearly related to the mobile phone-based visitor flows (R^2 is 0.958 at county scale and 0.953 at state scale) and keep consistent distributions (Figure <ref>b and <ref>c).
  	The distribution was slightly different at census tract scale. 
  	Though most scatter points are distributed following the fitting line (R^2 is 0.865 at census tract scale), those points with relatively high visitor flows and populations flows are located below the fitting line (Figure <ref>a).
  	The reason is that most paired census tracts have only a few visits and scatter points aggregate near the coordinate origin (and these points dominate the fitting line slope parameter).
  	Therefore, the distributions of visitor flows and population flows keep consistency especially for those paired regions with a small number of flows but have larger uncertainty for regions with large number of flows.
  	In sum, it can be concluded that the computed mobile phone-based visitor flows and the inferred entire population flows generally have linear associations across different spatial scales. Since the overall flow distributions are similar and generally consistent across geographic scales, the population flows inference process is reliable.
  	
  	
  	
  
   §.§ Comparison with Gravity Model and Radiation Model Estimations
  
  	
  	The gravity-style and radiation-style models that have been widely used in modeling spatial interaction patterns and human mobility O-D flows <cit.>. Thus, we also employ a classic gravity model and a radiation model to estimate the visitor/population flows between counties and then examine the correlation between the model outputs and our dataset estimates.
  	The comparison experiments are performed at the county scale using the weekly flow data considering that most of COVID-19 open data are available at the county scale in the U.S. <cit.>.
  	
  	The gravity model assumes that the magnitude of the flow between two places is related to their nodal attraction and in reverse proportion to the distance between them as shown by the following equation <cit.>. 
  	
  	
      F_ij = k× P_i× P_j/d_ij^β
  
  	
  	where F_ij represents the flow between county i and j; P_i and P_j are the total number of mobile devices (or the population) in the county i, j respectively; d_ij is the distance between the two counties. k and β are parameters that are used to adjust the model to fit the observations. For each week's data, we use the Particle Swarm Optimization technique <cit.> to find the best parameter combination of k and β that minimize the RMSE (Root Mean Square Error) between the estimated flow and the directly computed mobile phone-based visitor flow for that week. To make sure the estimated flow and the visitor flows are at the same scale, the range of k is set to be between 0 to 0.001. The range of β is between 0 and 2. The gravity model parameter fitting results for each week between 03/02/2020 and 05/31/2020 are shown in Table <ref>. The value of β ranges from 0.85 to 1.00. After the statewide lockdowns in the U.S. <cit.>, the distance decay coefficient β increased slightly, which showed that there were fewer longer-distance trips compared to that before the pandemic.
  	Then we applied the parameters in the census population to estimate the population flow of each week. 
  	The Pearson's correlation coefficient between the gravity model output and the population flow estimates (pop_flows) is about 0.56 to 0.64.
  	
  	One problem with the gravity model is that it can only capture undirected spatial interactions as the origin and the destination are not specified.
  	So we also use the radiation model that considers the direction of the flow for estimation. The original radiation model is used to estimate the commuting flows based on the U.S. census population distribution <cit.>. Here we applied this model to estimate the visitor flows based on the number of mobile phone users that travel across counties, the flow T_ij from county i to county j can be computed as follows:
  	
  	
      T_ij = T_im_i × n_j/(m_i + s_ij)×(m_i + n_j + s_ij)
  
  	where m_i is number of mobile phone devices in the source county i; n_j is number of mobile phone devices in the destination county j; s_ij is the total number of mobile phone devices in the counties that are within the circle of radius r_ij (distance between i, j) centered at county i (exclude i, j). T_i represents all the visitor flows that come from county i, and the model assumes that this is proportional to the total number of movement population of the source county <cit.>. Therefore, we have
  	
      T_i = m_i(N_c/N)
  
  	
  	where N_c is considered as the total number of devices that move (having at least one trip) during the study period, N is the total number of devices. In our daily data, we have the attribute that indicates the number of devices that are completely at home during that day. Therefore, we used the subtraction of the total number of devices and the number of `completely at home' devices as the N_c. The weekly N_c is the average of the 7-day's N_c values. For each week, we estimated the flow T_ij for all counties and computed the Pearson's correlation between the estimated flow and the `pop_flows' in our dataset (in Table <ref>). The correlation between the radiation model estimates and the `pop_flows' is about 0.75, which is better than that of the gravity model estimates.
  	
  	As there is no ground-truth, such a comparison with the gravity model and the radiation model provides a cross-referencing aspect and can inform the potential data users about their differences. 
  	
  	
  	
  
   §.§ Comparison with Other Data Sources
  
  	To illustrate that the data quality is high, we then compare two openly available data sources to the produced dataset from two different aspects, namely O-D flow patterns and temporal patterns.
  	Correlation analysis is conducted to check if mobility patterns of the compared two datasets differ.
  	Though there is no ground-truth data which characterizes the real dynamic population flows between two geographic regions as of yet, the comparison results are still useful for evaluating the credibility of the produced mobility flow dataset.
  	
  	In terms of the O-D flow patterns between two regions, we take the ACS commuting flows as baseline to see if the mobility patterns before the pandemic detected in our data products have a high correlation with the ACS commuting flows patterns <cit.>.
  	ACS commuting flows are generated by asking participants about their residence locations and primary workplace locations, and such an O-D flow dataset informs the understanding of interconnectedness between communities.
  	ACS commuting flows provide data at county (and county subdivision) scales <cit.> and thereby our comparison experiments are conducted only at the county scale.
  	The ACS commuting flows reflect the general patterns of commuting flows in the U.S., i.e., normal flows not affected by the COVID-19 pandemic.
  	Due to the fact that a national emergency concerning the COVID-19 pandemic was declared in the U.S. on March 13, 2020 and followed by statewide stay-at-home orders <cit.>, three time slices are picked up for comparison as they are supposed to reflect the mobility patterns before, during, and after the stay-at-home orders, respectively.
  	For daily flow data, we chose March 2nd, April 6th, May 11th, and May 25th as examples for comparison, while for weekly flow data we chose the weeks of March 2nd-8th (before the orders), April 6th-12th (during the orders), May 11th-17th (several states reopened after lifting the orders), and May 25th-31st  (business reopened in most states after the Memorial day <cit.>) for comparison.
  	
  	Comparison results between our produced mobility flow dataset and the ACS commuting flows data are illustrated in Table <ref>.
  	Regarding the number of records, the ACS commuting flows data contains 137,806 rows, while the generated mobility flow dataset at the county scale contains more origin to destination pairs than the ACS commuting flow data.
  	Our mobility flow dataset has a higher quantity of spatial interactions between county to county as more types of flows (including not just the home-to-job commuting flow) are captured.
  	The total number of records before the stay-at-home orders is greater than the number of records during and after stay-at-home orders, which is an intuitive consequence of people reducing their movements during the pandemic.
  	After joining the two datasets and removing records with no O-D matches, as is shown in Table <ref>, both weekly flow data and daily flow data had a high agreement (with greater than 0.93 Pearson's correlation) with the ACS commuting flow data.
  	In particular, the inferred entire population flows have higher correlation coefficients with ACS commuting flows than the directly computed mobile phone-based visitor flows.
  	With respect to the temporal changes, the flow patterns before stay-at-home orders have higher correlation coefficients, which matches expectations as people reduced their mobility when the stay-at-home orders in place <cit.>, leading to a decrease in the correlation coefficient values.
  	The comparison results illustrate that our proposed dataset may be complementary to the official ACS commuting flow data.
  	In addition, different from the ACS commuting flow data which only provides every 5-year updates at county scale and at minor civil division scale in major metropolitan areas, our dataset enables researchers to explore human mobility patterns at multiple spatial scales and in higher temporal resolutions (i.e., daily and weekly time windows).
  	
  	Additionally, the temporal patterns of the introduced mobility flow dataset and the Descartes Lab's COVID-19 mobility changes dataset is compared.
  	Descartes Lab's mobility changes dataset provides a mobility index which shows the median of the max-distance mobility of all users detected in one specific region <cit.>.
  	Such an index has been widely used for monitoring the daily mobility changes in the U.S. <cit.>. These two datasets are assumed to have similar temporal trends.
  	To do so, the total number of mobile phone-based visitor flows and entire population flows in the daily mobility flow dataset and their mobility changes are matched in five U.S. metropolitan areas: New York, Los Angeles, Chicago, Houston, as they are the top four metropolitan areas with largest population in the U.S., and Seattle, where the first confirmed COVID-19 case was reported <cit.>.
  	Correlation analysis is performed to check if these two datasets capture similar mobility changes during the epidemic since March 1st, 2020.
  	
  	As is shown in the Table <ref>, the produced mobility flow datasets have high correlation coefficients (at least 0.92) with the Descartes Lab's mobility changes dataset in all five metropolitan areas.
  	It shows that the generated datasets capture similar mobility temporal patterns with other open data sources.
  	As mentioned above, there is no ground-truth in measuring such high resolution dynamic human mobility patterns for the entire population, and as such it is therefore impossible to know which one can characterize human mobility more accurately.
  	However, by cross-referencing other sources, high correlation coefficients indeed show the reliability of our generated mobility flow dataset. 
  	
  	
  	
  
  § USAGE NOTES
  
  	Different datasets have their own pros and cons. 
  	In this data descriptor, we introduce two types of datasets in characterizing human dynamic O-D flows: daily flows and weekly flows, each at three geographic scales (i.e., census tract, county, and state).
  	A detailed description of the statistics and distributions of the datasets may benefit different applications.
  	Here, we discuss the characteristics and limitations of our mobility flow dataset to guide potential usages.
  	
  	Figure <ref> shows the temporal changes of the total number of mobility flows in five metropolitan areas: New York, Los Angeles, Chicago, Houston, and Seattle.
  	Accordingly, the daily flow data provides a more detailed temporal pattern description of human mobility changes compared to the weekly flow data. 
  	As illustrated in Figure <ref>, the temporally-changing curves have more fluctuations over time as human mobility patterns might be influenced by various factors and individual events (e.g., statewide lockdowns, presidential primary election day, Memorial day, or street protests), while weekly flow data reflects more general mobility patterns;
  	the temporal curves for the weekly flow data are more smooth and stable in comparison to the daily flow data.
  	For example, people may not visit the supermarket everyday and therefore the visitor volumes may vary by day of a week, while the sum of volumes in one week is more stable.
  	
  	Figure <ref> shows the temporal variations of daily and weekly active user sample size and the ratio of users having at least one trip (N_c/N) respectively. We find that the number of active users and the N_c/N values change over time and have a similar changing trend to the total number of mobility flows. One would expect fewer active users to be detected after the stay-at-home mandates as more people stayed at home more frequently.
  	The weekly number of active devices detected in the dataset reached over 22 million at the beginning and decreased to a minimum of 13 million active users in the week of April 13-19, 2020.
  	Similarly, as the outbreak progressed, the daily number of active users identified in the dataset decreased to a minimum of 16 million residing users on April 18.
  	After that, more active users were detected due to reduced compliance with stay-at-home mandates and the eventual lifting of stay-at-home orders.
  	
  	Figure <ref> shows the spatial patterns of O-D flow changes during the COVID-19 pandemic based on the weekly flow mobility data.
  	Figure <ref>a and <ref>b show the spatial interaction patterns of population flows across the Contiguous U.S. at the state and county scales, respectively. 
  	Figure <ref>c shows the population flow patterns at the census tract scale in the New York metropolitan area (New York was selected as it had the most confirmed COVID-19 cases during the study period).
  	We take four weekly flow data to represent the spatial patterns of mobility flows before (March 2nd to March 8th), during (April 6th to April 12th), and after (May 11th to May 17th and May 25th to May 31st) the stay-at-home orders.
  	At all three spatial scales, the movement flows decrease significantly from March to April due to the stay-at-home orders, with certain increases in May with the start of state partial reopenings. The decrease and subsequent increase in human movement flows clearly show how the outbreak of the COVID-19 pandemic and social distancing-related policies affect human mobility changes.
  	In particular, according to Figure <ref>b, during the stay-at-home order period long-range spatial interactions decrease to a small quantity, while most human movements appear as short-range movements to the adjacent counties.
  	When cities begin to reopen in May, long-term spatial interactions gradually rebound at both state scale and county scale.
  	
  	
  	
  	
  	While we made great efforts to guarantee the reliability of our produced mobility flow datasets and reduce the data uncertainty, a few limitations should be acknowledged.
  	
  	First, since we acquire anonymous mobile phone data from SafeGraph, the privacy policies applied should be considered by the end-user, as they may influence data uncertainty.
  	The weekly flow data is derived from CBG to POI visitor flow metric.
  	Please note that when a CBG has no visitors or has only one visitor who originate from that CBG to one POI, the visit between CBG to POI will not be recorded.
  	If the CBG has two to four visitors who originate from that CBG to another POI, the visitor count will be recorded but shown as four to enhance differential privacy.
  	Only when the CBG has more than five visitors whose home places inside the CBG to one POI will the real visitor count be displayed correctly.
  	Consequently, the least number of visitors from a CBG to POI as well as the weekly flow data computed in one week is four.
  	In addition, only when a visitor stays at a place (CBG or POI) for more than one minute will  a “visit” be counted. Such a relatively short time period may also affect the foot-traffic counting and cause data uncertainty.
  	
  	
  	Second, in terms of foot-traffic counting criteria, it should be noted that, for the daily flow data, it records the unique visitors to one CBG in one day, while for the weekly flow data it records the unique visitors to one POI in one week.
  	Due to the different counting methods, the sum of daily flow data in one week (i.e., seven days) are not equal to the weekly flow data. 
  	For the daily flow data, which is computed directly based on the daily CBG to CBG visitor flows, it may underestimate the mobility flows to POIs within one CBG.
  	That is, a visitor may visit more than one POI inside that CBG (e.g., different buildings in a campus) during the same day, it would still only be recorded once from the origin CBG to the destination CBG.
  	For the weekly flow data, which is based on the weekly CBG to POI visitor flows, mobility flows to POIs in one week may be underestimated as a visitor may visit one POI multiple times (e.g., work places) in one week, which would only be counted once as we compute the unique visitor devices rather than raw visits.
  	
  	In addition, visitor duplication may also exist when aggregating the flows from lower level spatial scales to upper level spatial scales (e.g., from CBG level to census tract scale, county scale, or state scale), and it may lead to the inflation of the inferred dynamic population flows.
  	For example, a user may originate from his/her home CBG and visit two POIs (e.g., one coffee shop and one grocery store) which are located inside different CBGs but at the same census tract.
  	Such a visitor should be counted as one unique visitor from the home census tract to the POI census tract. 
  	However, when aggregating the visits from CBG to CBG mobility flow to census tract to census tract mobility flow, because the individual travel behaviour cannot be traced, the mobility flows between these two census tract are counted as two and thereby inflate the real mobility flows.
  	
  	
  	
  	
  	Last but not least, data bias is a common issue for large-scale mobile phone data that may influence the representativeness of our produced dataset.
  	While dynamic mobility flows are inferred from mobile phone applications by users, not everyone in the population has a mobile phone, and not everyone uses smartphone applications, especially elderly people and children. 
  	Given these differences in mobile phone usage, age groups and demographic composition might influence the estimated entire population mobility flows. With these caveats in mind, this data still provides important up-to-date mobility flow information at scale.
  	
  	In conclusion, such a timely-produced dynamic O-D human flow dataset at different geographic scales can help deepen our understanding of human dynamics under this public health crisis, inform public health policy making, and support many other social sensing and transportation applications. 
  	
  	
  	
  	
  	
  
  § CODE AVAILABILITY
  
  	Data processing and data analysis were performed on a Linux server using the Python version 3.7.
  	All codes used for analysis are available in the public GitHub repository that hosts the data: <https://github.com/GeoDS/COVID19USFlows>. 
  	
  	
  
  § ACKNOWLEDGEMENTS
  
  	We would like to thank the SafeGraph Inc. for providing the anonymous and aggregated place visit data. S.G. acknowledges the funding support provided by the National Science Foundation (Award No. BCS-2027375). Any opinions, findings, and conclusions or recommendations expressed in this material are those of the author(s) and do not necessarily reflect the views of the National Science Foundation. Support for this research was partly provided by the University of Wisconsin - Madison Office of the Vice Chancellor for Research and Graduate Education with funding from the Wisconsin Alumni Research Foundation.
  	
  	
  	
  
  § AUTHOR CONTRIBUTIONS
  
  	Research design and conceptualization: S.G., Y.K.; Data collection and processing: Y.K., S.G.; Result analysis: Y.K., S.G., Y.L.; Visualization: M.L., J.R., Y.K., S.G.; Writing: all authors. 
  	
  	
  	
  	
  
  § COMPETING INTERESTS
  
  	
  	The authors declare no competing interests.
  	
  	
  
  § FIGURES
  
  	
  	
  	
  	
  	
  	
  	
  	
  	
  	
  	
  	
  	
  	
  	
  	
  	
  	
  	
  	
  	
  	
  	
  	
  	
  	
  	
  	
  	
  	
  	
  	
  	
  
  § TABLES
  
  	
  	
  	
  	
  	
  	
  	
  	
  	
  	
  	
  	
  	
  	
  	
  	
  	
  	
  	
  	
  	
  	
  	
  	
  	
  	
  	
  	
  	
  	
  	
  	
  	
  	
  	
  	
  	
  	
  	
  	
  	
  	
  	acm
  	
  	
  	
  	
  	
  