Levels of domain-specific physical activity at work, in the household, for travel and for leisure among 327 789 adults from 104 countries

Objective To compare the country-level absolute and relative contributions of physical activity at work and in the household, for travel, and during leisure-time to total moderate-to-vigorous physical activity (MVPA). Methods We used data collected between 2002 and 2019 from 327 789 participants across 104 countries and territories (n=24 low, n=34 lower-middle, n=30 upper-middle, n=16 high-income) from all six World Health Organization (WHO) regions. We calculated mean min/week of work/household, travel and leisure MVPA and compared their relative contributions to total MVPA using Global Physical Activity Questionnaire data. We compared patterns by country, sex and age group (25–44 and 45–64 years). Results Mean MVPA in work/household, travel and leisure domains across the 104 countries was 950 (IQR 618–1198), 327 (190–405) and 104 (51–131) min/week, respectively. Corresponding relative contributions to total MVPA were 52% (IQR 44%–63%), 36% (25%–45%) and 12% (4%–15%), respectively. Work/household was the highest contributor in 80 countries; travel in 23; leisure in just one. In both absolute and relative terms, low-income countries tended to show higher work/household (1233 min/week, 57%) and lower leisure MVPA levels (72 min/week, 4%). Travel MVPA duration was higher in low-income countries but there was no obvious pattern in the relative contributions. Women tended to have relatively less work/household and more travel MVPA; age groups were generally similar. Conclusion In the largest domain-specific physical activity study to date, we found considerable country-level variation in how MVPA is accumulated. Such information is essential to inform national and global policy and future investments to provide opportunities to be active, accounting for country context.


INTRODUCTION
The physical and mental health benefits of physical activity are well established, 1-3 yet 27.5% of the global adult population do not reach the minimum recommended levels. 4 The World Health Organization (WHO) member states agreed to a 15% reduction in physical inactivity levels by 2030 2 but this is unlikely to be achieved, given current trajectories. 4 Physical activity is a complex behaviour; however, opportunities to be active exist in several domains in life: at work, household or at school, for travel or during leisure time. Different domains may contribute to health in different ways. 5 Understanding the composition of how and where people accumulate their activity has important implications for policy, clinical practice and future public health investments. Economic development and its associated technological, demographic and societal changes have the potential to influence the opportunities for activity across different settings. 6 7 One way of investigating these patterns is by comparing domain-specific physical activity behaviour between countries at different stages of economic development.
Occupational physical activity is influenced by the country's urbanisation level and occupational structure, with more urbanised, service-dominated countries showing lower levels. 8 Leisure activity levels appear to be higher in high compared with low-income countries, 9 although even between highincome countries, levels of urbanisation and wealth seem to be positively associated with leisure activity prevalence. 10 Active travel patterns are less clear by income group as high-income countries provide examples of both high and low levels. For instance, active travel accounts for half of all journeys in the Netherlands and Switzerland but only a quarter in the UK and California, 11 suggesting differences in available transport infrastructure, social norms or both. This illustrates the complexity of the global perspective on domain-specific physical activity levels, something that has not been systematically researched.
Older age and female sex are known correlates of lower overall physical activity levels. 12 13 While a number studies have shown differences in the absolute levels of domain-specific activity, [14][15][16] few have explored how the relative contributions by age and sex differ. 17 18 Understanding whether different at-risk subgroups rely on certain physical activity domains more than others is critical information for policy-makers to protect existing and promote new behaviours and to help reduce health inequalities.
This study uses the latest internationally comparable physical activity data to compare the absolute and relative contributions of physical activity at work and in the household, for travel and during leisure time to total moderate-to-vigorous physical activity (MVPA) for 104 countries across the spectrum of economic development. As a secondary aim, we also compared the domain-specific MVPA contributions by sex and age group.

Data sources
We obtained the most recent individual-level survey data from 104 countries and territories (hereafter called countries) in which the WHO Global

Original research
Physical Activity Questionnaire (GPAQ) was used. 19 These included 94 surveys using the WHO STEPwise approach to noncommunicable disease surveillance (STEPS), 20 6 from the WHO Study on global AGEing and adult health (SAGE), 21 and 4 from public archives (see Acknowledgments and Data Availability statements for details).
The countries covered were from all six WHO regions: Africa (n=37), Americas (n=15), Eastern Mediterranean (n=12), Europe (n=8), South-East Asia (n=8) and Western Pacific (n=24). According to World Bank Income Classification (2020), 22 24 were low, 24 were lower-middle, 30 were uppermiddle, and 16 were high-income countries (online supplemental file 1). These countries make up 71% of the global population in 2020. 23 In addition, an earlier STEPS survey was identified for four countries in the African region (Benin, Botswana, Malawi, Seychelles) for an exploratory analysis of domain-specific trends. These were countries within the same region, spanning all income classification groups, and with national data from two points in time with a minimum time gap of 5 years between them and the most recent time point within the last 5 years.

Representativeness of data
Survey sampling was designed to obtain a nationally representative (n=90 + 4 trend datasets) or a subnational (n=14) sample. In the majority of cases, this was a stratified multistage clustered design but some of the subnational surveys or smaller islands used simple random, quota or interval sampling, or a census (online supplemental file 2). Sample weights were provided for 91 (+3 trend) surveys to adjust for unequal selection probabilities. These also adjusted for the oversampling of those ≥50 years in the SAGE surveys, and normalised to population distributions of age, sex and in some cases region and ethnicity in the publicly available surveys.
Eight of the subnational surveys were urban only. For these, we imputed rural data and weighted them by national urbanrural prevalence to derive national estimates. In the urban-rural mixed subnational survey, we performed this weighting step with the data collected. The remaining five subnational surveys were analysed with no adjustment. Full methods are provided in online supplemental file 3.
The analysis was limited to the common age-range across all surveys (25-64 years) for comparability reasons so the estimates presented in this paper are not official national statistics.

Global Physical Activity Questionnaire
We only included surveys which had used the WHO GPAQ. This instrument captures MVPA undertaken in a typical week at work (paid or unpaid, including household chores), for travel to and from places by walking or cycling, and during leisure time (including sports and fitness-enhancing activities). 19 Work/household and leisure domains are split into moderate and vigorous intensity. Only MVPA undertaken in bouts of 10 min or more is reported, and participants are reminded not to report activities already included in an earlier answer. The GPAQ was translated into national languages, and back translated for quality assurance, and show-cards were used to provide culturally relevant examples of physical activities in each domain.
The GPAQ has shown comparable levels of validity to other self-report instruments across a variety of cultural settings. [24][25][26][27][28] In the absence of domain-specific validity evidence, we estimated the Pearson's correlation coefficients for the associations between the country-level domain-specific MVPA (median min/ week) and indicators selected a priori (online supplemental file 4). These were work/household and the proportion of the workforce in the agricultural sector (males: r=0.43, females: r=0.46) and other occupations involving a high degree of manual labour (males: r=0.47, females: r=0.48); travel and the number of vehicle registrations per adult in the population (r=−0.49); and leisure and the GINI measure of economic inequality across the population (r=0.03), and the Human Development Index (r=0.27).

Data processing
The standard GPAQ data processing protocol was followed. 19 Individuals reporting implausible (>16 hours/day in one of the domains) or inconsistent (eg, no days but valid duration) were excluded. We considered moderate and vigorous intensity work/household and leisure physical activity separately for this step. We further excluded those with missing data in either the frequency or duration variables in any domain (online supplemental file 5). We calculated total weekly duration in minutes for each individual by multiplying the number of days and daily time variables for each domain and summing the totals. For the main analysis, we combined the moderate and vigorous components of work/household and leisure as reported. In a sensitivity analysis, we doubled the reported duration of vigorous intensity activity to reflect the physical activity guidelines that suggest 75 min of vigorous intensity activity provides comparable benefits to 150 min of moderate intensity activity. 3 We estimated the absolute levels of domain-specific MVPA (min/week) at a country level, and by sex and age group (25)(26)(27)(28)(29)(30)(31)(32)(33)(34)(35)(36)(37)(38)(39)(40)(41)(42)(43)(44) and 45-64 years) strata based on all those meeting the inclusion criteria. We calculated relative contributions at an individual level, and group-level arithmetic means were presented. Those who reported no MVPA across all domains could not provide any information to the relative contributions (online supplemental file 5). We described differences with reference to the World Bank Income Classification 2020. 22 We estimated Pearson's correlation coefficients for the associations between the mean domain-specific MVPA and total MVPA at a country level. We performed the analyses in Stata V.16 (StataCorp) and produced the ternary plots in R Studio using the ggtern 29 and tricolore packages. 30

Public involvement
The experiences of data collectors informed the development of the GPAQ, 24 and assisted in the cultural and language adaptations made for each country.

RESULTS
We included data from 327 789 individuals aged 25-64 years from 104 countries collected between 2002 and 2019. The prevalence of meeting the WHO physical activity guidelines ranged between 36% and 98%, and the proportion reporting zero min/ week of total MVPA ranged between 1% and 51% (table 1).

Absolute levels of domain-specific activity
The distributions of the reported weekly durations of work/ household, travel and leisure MVPA were highly skewed for all countries; the mode was almost always zero min/week (online supplemental file 6). The country-level means for the domains were 950 (IQR 618-1198), 327 (190-405) and 104 (51-131) min/week, for work/household, travel and leisure time, respectively (table 2, figure 1, online supplementary file 6). Mean work/ household MVPA in low-income countries was almost twice that

Original research
of high-income countries (1233 vs 668 min/week). However, the reverse was seen for mean leisure MVPA (72 vs 143 min/week). Travel MVPA was approximately three times higher in lowincome versus high-income countries (499 vs 158 min/week). We found a very strong correlation between work/household MVPA and total MVPA (r=0.95) at the country level. The correlation was 0.58 for travel and 0.04 for leisure MVPA.

Relative contributions of domain-specific activity to total MVPA
The 47 946 individuals who reported 0 min/week of total MVPA did not contribute to the relative analyses (denominator of zero; online supplemental file 5). Based on the remaining 279 843, the mean contribution of work/household MVPA to total MVPA was 52% (IQR 44%-63%); the respective means for travel and leisure were 36% (25%-45%) and 12% (4%-15%; table 2). Work/household domain was the largest contributor to total MVPA in 80 of the 104 countries, and 70 of these had travel as the second-largest contributor. Travel was the largest contributor in 23 countries, and leisure was dominant in only one. There was a tendency towards higher contributions of leisure MVPA in high-income countries (mean of 28% compared with 4%, 8%, and 13% in low, lower-middle, and uppermiddle-income countries, respectively; table 2, figure 2, online supplementary files 7 and 8). The mean contributions of work/ household MVPA in low and lower-middle-income countries were 57%, but were 44%-47% for upper-middle and highincome countries. The contribution of travel to total MVPA did not follow a clear pattern by World Bank income group.
Doubling the duration of vigorous intensity activity (ie, 1 min of vigorous physical activity=2 min of moderate intensity physical activity) in the work/household and leisure domain obviously influenced the absolute durations but did not affect the relative contributions in a meaningful way. This was because those individuals reporting high levels of work/household vigorous intensity activity were already obtaining close to 100% of their MVPA in this domain and so further increases in the absolute durations did not shift the balance of the domain contributions. This conversion therefore had limited influence on the countrylevel mean relative contributions which varied by a maximum of 2 percentage points for all domains and the highest contributing domain only changed in five countries.
There was a tendency towards greater travel and lower work/ household contributions among women compared with men,

Exploratory analysis of domain-specific trends
We conducted exploratory analysis of trends using a subsample of data from four African countries (n=25 749), spanning the income classifications ( figure 4). The results show considerable heterogeneity in the directions and magnitudes of the changes in domain-specific MVPA over time, even within income classifications.

Main findings
This study is the largest comparison of country-level domainspecific physical activity to date, including 372 789 individuals from 104 countries. We found the work and household domain to be the highest contributor to total MVPA levels in threequarters of countries. Travel was the second-largest contributor for all but one of the remaining countries, and was rarely the lowest contributor. There was a trend towards greater contributions of MVPA through work/household in lower income countries, and higher contributions of physical activity through leisure in high-income countries. This pattern was suspected but is shown quantitatively for the first time in a large sample of representative surveys across multiple regions. There were differences in the relative contributions of the domains by sex and age group within countries, with a tendency towards greater differences in upper-middle and high-income countries.

Comparison with previous results
Previous national-level domain-specific analyses have been focused on specific regions, for example, Africa, 31 South America, 16 Asia-Pacific, 14 while other cross-national studies have investigated the proportions reporting a certain threshold of domain-specific activity (eg, 0, 10, 150 min/week) which limit comparisons with the present results. 15 16 32 Our results extend the range of comparisons to be cross continental allowing greater comparisons between different World Bank income classifications, and show the relative importance of the domains in addition to the absolute values. Although high-income countries are relatively under-represented in our analyses, the 16 that are included act as helpful comparators to contrast the dominance of work/household-related activity and the low prevalence of leisure activity among many low and middle-income countries. Comparisons in total and domain-specific physical activity by country-level and individual-level income groups is an important yet under-researched area of global physical activity surveillance. However, it is clear that factors beyond income level and economic development, such as cultural and historical aspects, likely play a role in explaining differences between countries.
Our exploratory analysis of changes over time within domains is an indicator of this complexity. Some but not all countries showed a decline in the work/household domain, a feature that was expected as the epidemiological transition is accompanied by a shift away from manual labour and increased access to labour-saving technology in the home. 8 The variations seen in the work/household and travel domain trajectories is likely to be explained by differences in countries' socioeconomic

Original research
circumstances, level of urbanisation and rates of development. Urbanisation and uptake of digital technologies as well as other cultural, environmental and social factors influence the opportunities for increasing levels of participation in the different domains in different ways. These results indicate the importance of further research on this topic with a larger set of comparable data. It is also a reminder that while global comparisons can identify broad patterns, this needs to be coupled with a local understanding of the context in order to develop appropriate policies and interventions.
Building on previous research that has shown lower levels of overall MVPA among women compared with men across the majority of countries worldwide, 4 12 we showed that the domain composition also varies. The profile for women, at least in high-income and upper-middle-income countries, had a slightly higher contribution of travel and lower contribution of work/ household. These results provide a more nuanced understanding of how physical activity varies by sex.
We did not identify any notable differences in the relative contributions of the activity domains between 25-44 and 45-64 year olds, although where differences were found, they tended to be in upper-middle and high-income countries. Research from UK samples indicated that the greater shifts in domain-related contributions occur after 65 years, that is, around the age of retirement. 17 18

Implications of the findings
These results can inform the development of interventions and policies aimed at increasing overall physical activity levels as called for by the 2018 Global Action Plan on Physical Activity. 2 The high contribution of work/household activity for many low and middleincome countries is important for a number of reasons. First, it identifies a reliance on this domain for meeting activity recommendations that may not persist as countries develop economically. 7 8 Policy action is necessary to ensure that alternative types of activity are available and affordable to those who might otherwise decrease their activity levels as consequence of such societal and economic changes. The travel domain is an obvious area of focus as our results show that it already is an important contributor to overall levels across all countries. Policies and infrastructure development that promote walking and cycling would also support efforts to reduce air pollution and carbon emissions to combat climate change as set out in the Sustainable Development Goals. 33 Second, these data stress the importance of assessing all domains of physical activity, since solely concentrating on the leisure domain (as has often been the case 34 ) would grossly misrepresent, and for many populations underestimate, overall physical activity prevalence. 35 Third, if domain-specific activity were to show any differential health benefits, as recently debated, 36 mapping absolute and relative profiles of domain-specific activity is critical for estimations of global disease burden.

Limitations
Methodological limitations include the unavailability of sampling weights for 13 countries and the necessity of using imputation methods to account for urban-rural variations when including subnational samples. Our aim was to be as inclusive as possible with the available data. Nonetheless, imputation likely increases the error regarding the representativeness of these estimates, over and above that of the national sampling. To maintain comparability between samples, we restricted the age range to 25-64 year olds. For these reasons, the descriptive statistics of overall activity and the proportions meeting the guidelines do not replace any official national statistics and are presented for context only. The novel estimates of the domain-specific relative contributions are valid in the context of this paper, but further work is necessary to generate official statistics that include adults of all ages, and the multitude of other national survey data collected using alternative questionnaires.
The validity of domain-specific GPAQ estimates has not been evaluated against an appropriate criterion, such as a domain log. 37 Some have raised concerns about the implausible work/household values. 38 Inaccurate recall, social desirability and activity misclassification are possible sources of error. In lieu of such a validity study, the correlations presented between the work/household and travel domains and the selected country-level indicators provide tentative support for domain-specific analysis using GPAQ data. We believe the weaker results for the leisure domain are in part due to the number of countries with a median MVPA of zero min/week and the lack of specificity in the macro-economic indicators. In addition, the GPAQ states a minimum 10 min bout length, a requirement that is no longer included in the WHO physical activity guidelines. 1 The inclusion of shorter bouts might alter the duration of reported activity across all domains, but it is difficult to predict which domains would be more affected as the magnitude of these increments is likely context specific. This would influence the relative contributions. Also, the GPAQ does not contain questions about light-intensity activity which is a major contributor to overall activity. 39 It is likely that its inclusion would change the relative domain-specific results as well the absolute.
Finally, in the relative contribution analysis, we presented the arithmetic mean values (of individuals for each country, and of countries within income classifications). Although the scaled geometric mean is often used in compositional data, this was not an option when so many individuals reported 0 min/week in at least one domain. However, our analysis method is somewhat protected against individuals reporting very high values of domain-specific activity as the relative contributions were calculated at an individual level, before being summarised at a country level. This means a high value in one domain is limited to a maximum of 100% relative contribution. Of course, this does not protect against issues of differential validity between the domains, which is a plausible concern as even the ordering of questions can affect reporting. 40 In the GPAQ, the work/ household domain is asked about first and all surveys used the same ordering. There is also the potential for cultural differences in the reporting as there are different interpretations of work/household, travel and leisure activity across populations. We were also unable to account for seasonal differences as surveys are often conducted within a short period of time that does not span all seasons.

Future directions
Filling in the data gaps will require wider adoption of the GPAQ in national surveys or harmonisation between different existing questionnaires. This would be pertinent for the European region as many countries have nationally representative data using alternative questionnaires. Nonetheless, it is rare in physical activity studies for the majority of data to be from low and middle-income countries and so this present work does provide an important perspective on the topic.
Future work could extend these methods to investigate subnational regional differences that others have shown to exist. 41 42 Exploring variation by socioeconomic differences would also have utility, and require suitable harmonisation of indicators of socioeconomic position across nations.
We have also highlighted the need for analysis of trends in domain-specific data. As countries repeat their STEPS surveys, more comparable data will become available. In the meantime, efforts should be made to harmonise existing data.
Lastly, although the present study reflects the optimal analysis of currently available data, the advent of combined location and physical activity sensors may allow objective assessment of domain-specific activity in future research.

Conclusion
This study is the largest analysis of domain-specific physical activity to date comparing data from 104 countries. Activity at work/household is the main contributor to total MVPA levels, particularly in lower income countries. Leisure activity was the smallest contributor but was highest in high-income countries. Achieving the 15% decrease in global levels of physical inactivity by 2030 as outlined in the Global Action Plan on Physical Activity, 2 and agreed by 194 Member States of WHO 43 will require detailed understanding of the context in which people accumulate their physical activity, especially for countries currently undergoing rapid economic development and urbanisation. Across all nations, policy actions will need to align efforts to support and enable opportunities to be active in different domains according to the social, economic and demographic changes, population needs and local context.

What are the new findings?
► Physical activity in the domains of work and household, travel and leisure vary broadly across countries in both their absolute levels and relative contributions to total moderateto-vigorous physical activity (MVPA). ► Work/household MVPA provided the largest contribution in 80 of the 104 countries included in this study. ► Women tended to have relatively less work/household and more travel MVPA. ► Differences between 25-44 and 45-64 year olds were less apparent than differences between the sexes, but were most evident in high and upper-middle-income countries.

How might it impact on clinical practice in the near future?
► Clinicians and policy-makers alike should be aware that MVPA is accumulated across a variety of domains, and this varies across countries. ► Work/household is a dominant source of MVPA for many; this has implications for continuing to meet guidelines at the point of changing jobs or retirement. ► Women tended to have relatively lower contributions of work/household and more travel MVPA and this should be considered when tailoring behavioural advice.