Diurnal patterns of sedentary behavior and changes in physical function over time among older women: a prospective cohort study
International Journal of Behavioral Nutrition and Physical Activity volume 17, Article number: 88 (2020)
Sedentary behavior (SB) is linked to negative health outcomes in older adults. Most studies use summary values, e.g., total sedentary minutes/day. Diurnal timing of SB accumulation may further elucidate SB-health associations.
Six thousand two hundred four US women (mean age = 79 ± 7; 50% White, 34% African-American) wore accelerometers for 7-days at baseline, yielding 41,356 person-days with > 600 min/day of data. Annual follow-up assessments of health, including physical functioning, were collected from participants for 6 years. A novel two-phase clustering procedure discriminated participants’ diurnal SB patterns: phase I grouped day-level SB trajectories using longitudinal k-means; phase II determined diurnal SB patterns based on proportion of phase I trajectories using hierarchical clustering. Mixed models tested associations between SB patterns and longitudinal physical functioning, adjusted for covariates including total sedentary time. Effect modification by moderate-vigorous-physical activity (MVPA) was tested.
Four diurnal SB patterns were identified: p1 = high-SB-throughout-the-day; p2 = moderate-SB-with-lower-morning-SB; p3 = moderate-SB-with-higher-morning-SB; p4 = low-SB-throughout-the-day. High MVPA mitigated physical functioning decline and correlated with better baseline and 6-year trajectory of physical functioning across patterns. In low MVPA, p2 had worse 6-year physical functioning decline compared to p1 and p4. In high MVPA, p2 had similar 6-year physical functioning decline compared to p1, p3, and p4.
In a large cohort of older women, diurnal SB patterns were associated with rates of physical functioning decline, independent of total sedentary time. In particular, we identified a specific diurnal SB subtype defined by less SB earlier and more SB later in the day, which had the steepest decline in physical functioning among participants with low baseline MVPA. Thus, diurnal timing of SB, complementary to total sedentary time and MVPA, may offer additional insights into associations between SB and physical health, and provide physicians with early warning of patients at high-risk of physical function decline.
With accelerated growth of the older adult population , identifying ways to prevent disease and promote healthy aging is becoming a major public health focus. Older women have proportionately higher morbidity and disability than older men , despite having a longevity advantage. Because of this, women have more hospitalizations and outpatient visits, longer utilization of long-term care services, and greater health care spending [3,4,5]. These disparities make the identification of modifiable risk factors for declining health especially important for older women.
A growing body of recent research suggests that a promising risk factor to target is sedentary behavior (SB). In the United States, older adults are estimated to spend an average of 8 to 9.5 h per day in SB [6,7,8]. Older adults accumulate the greatest volume of sedentary time of any age group [9,10,11,12], and older age is associated with more sedentary time even among older adults . In older adults, SB has been linked with higher risk of cardiovascular disease, metabolic syndrome, diabetes, and all-cause mortality [14,15,16].
Accelerometers have been increasing our understanding of SB and how it relates to health. Early studies used self-reported measurements of SB, which can be unreliable or inaccurate, especially in older age groups [17,18,19]. Accelerometer data are subject to less measurement error, and analyses of these data have revealed stronger associations between SB and health than were previously observed using self-reports [20,21,22].
Accelerometers also enable researchers to study diurnal patterns of sedentary time . Studies analyzing objectively measured sedentary time data usually reduce them to a single measure—total minutes per day spent sedentary—overlooking diurnal accumulation patterns that may yield important information about health outcomes. While it is recognized that high levels of SB are negatively related to health, the association may vary depending on when sedentary time is accumulated during the waking day.
While there is a paucity of evidence relating diurnal patterns of SB with health, emerging evidence suggests that diurnal patterns of physical activity (PA) are associated with fatigability ; early-stage Alzheimer’s disease ; insulin and C-reactive protein levels, and quality of life ; sleep efficiency, cognition, and all-cause mortality . There have also been reports that diurnal patterns of SB vary by sociodemographic variables; for example, older adults accumulate more SB during the day than younger adults, with smaller age group differences in evening SB accumulation . These data support the need to study diurnal patterns of SB to determine possible associations with health-related outcomes.
Using a novel, two-stage clustering procedure, we investigated diurnal patterns of SB and how these were related to trajectories of physical functioning. Data were drawn from a study in which 6489 women wore ActiGraph accelerometers on their waist for up to 7 days at baseline, and then were followed for up to 6 years with annual assessments of physical functioning . We expected diurnal patterns of SB to be associated with trajectories of physical functioning among this cohort of community-living older women.
Sample and design
Data were used from participants in the Objective Physical Activity and Cardiovascular Health Study (OPACH), a subset of Women’s Health Initiative participants who had enrolled in the Long Life Study. Women’s Health Initiative women in the Long Life Study consented to periodic in-home examinations to provide blood samples, updated health information, and physical measurements. More details of the ancillary OPACH study, which was designed specifically for collection of objective measures of physical behavior, have been previously published . In brief, the 7048 ambulatory, community-dwelling women who consented to participate in OPACH were given ActiGraph GT3X+ accelerometers on a waist belt and asked to wear them 24 h per day (except when bathing or swimming) for 7 consecutive days. They concurrently recorded in-bed and out-of-bed times each day in sleep logs. Accelerometers were returned by 95.4% (6721) of participating women and 92.1% (6489) contained evidence of human wear . Each participant was followed for up to 6 years with annual medical updates collecting information on physical functioning.
Outcome: physical functioning
Physical functioning was assessed yearly using the RAND 36 Health Survey physical function subscale, a 10-item, well-validated self-report measure [30, 31]. The physical functioning scale assesses current health limitations on physical functioning during daily activities. Item scores are averaged, producing a score from 0 to 100, with higher values indicating superior physical functioning. Median number of physical functioning assessments available (25th %, 75th %) was 5 (4, 5) time points.
Raw accelerometer data (30 hertz) were converted using ActiLife software (Version 6) to vector magnitude counts per 15-s epochs. To remove accelerometer non-wear time, we applied the Choi algorithm (90-min window, 30-min stream frame, and 2-min tolerance) to the minute-level vector magnitude counts  and identified in-bed time using data reported in the sleep logs. Missing bed times were imputed using the mean in-bed and out-of-bed time of each person when available; for the 434 participants who had in-bed and/or out-of-bed data missing in all days, the mean in-bed (10:45 pm) and/or out-of-bed (7:22 am) time for the population was used. We defined adherent days as calendar days with 10 h or more of awake wear time , and only analyzed adherent days. Because our metrics of sedentary time and sedentary accumulation patterns were intended to estimate typical behavior during a week, at least 4 adherent days were required for inclusion of a participant’s data in analyses . Out of the 6489 women with wear time, 95.6% (6204) had at least 4 adherent days.
Accelerometer data were processed using intensity-specific cutpoints that had been developed specifically for the OPACH cohort in a laboratory calibration study of 200 women aged 60–91 years who wore ActiGraph GT3X+ accelerometers on their hip . Sedentary time was classified as any 15-s epoch at or below 18 vector magnitude counts and MVPA was any 15-s epoch above 519 vector magnitude counts.
Questionnaire data collected at Women’s Health Initiative baseline included age, race/ethnicity (Black, White, or Hispanic) and education (high school graduate/general education development (GED) or less, some college, college graduate/advanced degree). Participants completed questionnaires at OPACH baseline that included measures of self-reported health (excellent/very good, good, fair/poor), frequency of alcohol consumption (non-drinker, less than one drink/week, one or more drinks/week, unknown), and current smoking status (smoker vs nonsmoker). Measures of weight and height were taken during Long Life Study in-home visits using a calibrated analog scale and a tape measure, and used to compute body mass index (BMI) as weight (kilograms)/height (meters2). Morbidity was measured at OPACH baseline as the number of the following chronic health conditions: cardiovascular disease; cancer; osteoarthritis; cognitive impairment; depression; history of falls over the previous 12 months; cerebrovascular disease; chronic obstructive pulmonary disease; diabetes; vision impairment; hearing loss .
Sedentary minutes were summed for each one-hour interval of the day (e.g., 6:00 am-6:59 am, 7:00 am – 7:59 am) to describe hourly behavior. Inclusive of the first hour that participants were awake and wore the accelerometer for 60 min, daily trajectories of sedentary behavior were obtained using data from 14 consecutive hours. This consistent duration across days was required to implement the longitudinal clustering, and 14 h was selected based on the degree of missingness (i.e., 23.6% of data from the 14th hour was missing while 49.5% of data from the 15th hour was missing). For hours with missing data, sedentary time was imputed based on linear interpolation with added stochastic variation based on population trajectories within missing time points (a method known as “Copy Mean”) . An alternative method for computing trajectories would be to use a fixed range of clock time (e.g. 8 am to 10 pm) instead of beginning with participants’ wake time. That alternative method, while a reasonable approach to obtaining SB diurnal patterns, does not conform with our understanding that participants’ diurnal patterns commence with their highly variable time of awakening rather than at a fixed hour of the day. For our primary analyses, we adopted the method of beginning the day at the participant’s hour of awakening, but conducted a concordance analysis using the alternative method (i.e., using clock-time to define start-of-day).
A novel two-phase clustering procedure was then used to group OPACH participants according to their diurnal and day-to-day variability in SB patterns. In Phase I, k-means for longitudinal data (kml) was conducted on all available days of data from the entire sample (ndays = 41,356) to determine clusters of days that optimally differentiated diurnal trajectories. To select the number of clusters in Phase I, we used the Calinski-Harabasz criterion . To address cluster stability, the kml algorithm was repeated five times with different starting conditions, and the partition that led to the highest Calinski-Harabasz criterion was selected. In Phase II, hierarchical clustering with complete linkage  was conducted on the proportion of a participant’s measurement days assigned to each of the identified day clusters, to determine the diurnal SB patterns of participants. The proportion of days was used because the number of days available for each participant varied from 4 to 7 days. To determine the number of SB pattern clusters, we selected the solution that maximized average silhouette width across observations . Thus, Phase I clustered days based on similarities in within-day timing of SB accumulations, whereas Phase II clustered individuals based on similarities in their between-day SB patterns, e.g., how they distribute days within a week into various SB accumulation clusters.
Socio-demographic, health-related, and activity-related variables were summarized for each of the groups of participants clustered according to their diurnal SB patterns, using means and standard deviations (SD) for continuous variables and percentages for categorical variables with differences tested using F-tests and chi-square tests, respectively.
Multivariable linear mixed-effects regression tested whether diurnal SB patterns were associated with trajectories of physical functioning. This enabled investigation of differences in baseline physical functioning by diurnal SB pattern (aka, the y-intercept of the trajectory), and how the slopes of the physical functioning trajectories changed over time in relation to these patterns. The general model was physical functioningit = patterni + timeit + patterni*timeit + covariatesi, with subjects i = 1,…,N, and year of follow-up t = 1,…,T. A person-level random intercept was also included in the model to account for within-person correlations. Model covariates included age, race-ethnicity, BMI, education, smoking status, alcohol use, number of morbidities, and self-rated health. To test whether associations of diurnal SB patterns with trajectories of physical functioning reflected daily sedentary time, we additionally adjusted for total sedentary time. Based on previous literature, we examined whether associations of diurnal SB patterns and physical functioning varied by MVPA (which, as in our previous studies [15, 40], has been modeled as a binary variable using a median split to increase statistical power within each MVPA subgroup and because there is no older-adult-specific threshold that has been established using accelerometer data) . This was accomplished by adding a pattern*time*MVPA interaction and using a likelihood ratio test to assess statistical significance. The potential influence of using complete case analysis was assessed by using multivariate imputation by chained equations  with predictive mean matching and Bayesian polytomous regression for continuous and categorical variables, respectively.
All analyses were conducted in R (version 3.5.3), including the use of packages kml , cluster , nlme , and mice. All tests were two-tailed with alpha = 0.05 unless otherwise specified. Sample R code for our cluster analysis is available at https://github.com/Aging-and-Behavioral-Epidemiology/two_stage_clustering.
Study sample and baseline descriptive characteristics
Our analytic sample comprised 6204 OPACH participants, with mean (SD) age of 79 (7) years at study entry, 49.6% White, 33.6% Black, 41.0% completed college or higher. Mean (SD) BMI at baseline was 28.1 (5.7). 2.6% of participants were current smokers, and 25.8% had 1 or more alcoholic drinks per week in past 3 months before baseline. 89.2% of participants rated their general health as excellent, very good, or good. 17.5, 34, 26.8, and 21.6% had no, 1, 2, or 3 or more morbidities at baseline, respectively. Participants wore accelerometers for a mean (SD) 6.9 (0.3) days, yielding 41,356 days with mean (SD) 14.9 (1.3) hours/day of data. At study entry, mean (SD) daily sedentary time and MVPA were respectively 597 (103), and 50 (34) minutes. Mean (SD) physical functioning scores at baseline were 69 (26).
Diurnal SB patterns defined by the two-phase clustering procedure
In Phase I, k-means for longitudinal data (kml) was conducted on all available days of data (ndays = 41,356) to select the optimal clusters of days. We found a solution of four clusters that had a sufficiently high Calinski-Harabasz criterion , yet with enough variability among clusters to clearly differentiate diurnal trajectories; the chosen partition also had the lowest AIC and BIC. Based on visual inspection, the four day clusters can be described as: (A) an overall high level of SB throughout the day, (B) a moderate level of SB, with a lower level of SB at start of day, relative to the end of day, (C) a moderate level of SB, with a higher level of SB at the start of day, relative to the end of day, and (D) an overall low level of SB throughout the day (Fig. 1a). Of note, clusters A, B, and D showed increasing SB trends over the course of the day, whereas cluster C had higher SB earlier followed by ~ 6 h of a decreasing trend, and then increasing for the rest of the day.
Then in Phase II, hierarchical clustering was conducted on the proportion of each participant’s days that fell within each of the four Phase I day clusters. We found that a 4-cluster solution maximized average silhouette width across observations . Diurnal SB patterns 1, 2, 3, and 4 can be described according to which day cluster A, B, C, and D, respectively, their days predominantly fell into (see Fig. 1 b, c). For example, for individuals in diurnal SB pattern 1 the median (25th%, 75th%) proportions of cluster A (high SB) days were 66.7% (50, 85.7%) vs 0% (0, 0%) proportion of cluster D (low SB) days. Conversely, for SB pattern 4 the median (25th%, 75th%) proportions were 0% (0, 0%) vs 60.0% (50.0, 85.7%) for cluster A and cluster D days, respectively. SB patterns 2 and 3 had similar distributions of high SB (cluster A) and low SB (cluster D) days, but differed according to the proportions of cluster B vs C days: for SB pattern 2 the median (25th%, 75th%) proportions for cluster B vs C days were 57.1% (42.9, 71.4%) vs 28.6% (14.3, 33.3%), respectively, whereas SB pattern 3 had 14.3% (0, 28.6%) vs 57.1% (50.0, 71.4%), respectively. Thus, on average, SB patterns 1–4 were distinguished by the corresponding day clusters A–D, of which each pattern was primarily composed.
Concordance analysis for the clustering procedure
Of two alternative accelerometer processing procedures on which to base our clustering analyses, we chose the procedure using data aligned to the per subject start of day rather than data aligned to clock time of day. Reassuringly, when we conducted a concordance analysis to compare the procedures, findings were similar: the two procedures resulted in 81.4% of the days being classified in the same day cluster, and 73.2% of the participants being classified in the same diurnal SB pattern.
Descriptive statistics by diurnal SB pattern
Baseline participant characteristics by SB pattern are in Table 1. The low SB pattern  and high SB pattern (1) had characteristics that correlated with better and worse health, respectively. For example, those with the high SB pattern had on average higher age, BMI, proportion of current smokers, number of morbidities, total sedentary time; they also had lower self-rated health, MVPA time, and physical functioning. Both moderate SB patterns (2 and 3) were on average similar on most characteristics, with values closer to the sample averages.
Mixed-effects regression: associations between physical functioning and diurnal SB patterns
In the study sample, mean physical functioning decline was 2.2 points/year (standard error = 0.07). MVPA was a significant modifier of the association between diurnal SB pattern and physical functioning (likelihood ratio p < 0.0001), and therefore results are presented stratified by high and low MVPA (Fig. 2). Baseline physical functioning levels were significantly higher for women with high MVPA (ranging from 78.6 to 81.5) than for women with low MVPA (range from 70.6 to 75.1). Women in all SB patterns experienced decline in physical functioning over time. The decline in physical functioning over time was significantly slower in SB pattern 2 and 3 women with high MVPA (1.9 and 1.6 points/year, respectively) than the physical functioning decline among women in the same SB pattern with low MVPA (3.0 and 2.5 points/year, respectively).
Among women with low MVPA, those in the high SB pattern  had lower baseline physical functioning (70.6 ± 1.0) than all other patterns, which was significantly lower than both of the moderate SB patterns (pattern 2 physical functioning score = 75.1 ± 1.2; pattern 3 physical functioning score = 74.4 ± 1.3). The pattern having lower SB at the start of the day relative to the end of day (pattern 2) had the steepest negative slope of physical functioning decline (− 3.0 ± 0.2), which was significantly steeper than that of the consistently low SB pattern  and consistently high SB pattern . There were no significant differences in baseline physical functioning levels or changes in physical functioning over time between the two moderate diurnal SB patterns.
Among women with high MVPA, the women in SB pattern 2 had significantly higher baseline physical functioning (81.5 ± 1.0) than those in SB patterns 1 and 3 (78.9 ± 1.3, and 78.6 ± 1.1, respectively). The women in the consistently high SB pattern  had the most rapid decline in physical functioning over time, though not significantly different than that of the other diurnal SB pattern groups. The rate of physical functioning decline among women with low morning and high evening SB pattern  was not significantly different than that in the most rapidly declining group (pattern 1; p = 0.25). Women with a consistently low diurnal SB pattern  and women with a high morning and low evening SB pattern  had slightly slower declines in physical functioning (1.7 ± 0.1 and 1.6 ± 0.2 physical functioning points per year) than those in the persistently high SB pattern (1; 2.1 ± 0.2 physical functioning points per year; respective p-values = 0.10 and 0.08). Results were similar with and without adjustment for total sedentary time (Additional File 1).
Similar results were obtained with and without multiple imputation, likely due to the relatively small amount of missingness (records removed due to listwise deletion = 1796/27717, or 6.5%). As a result, all analyses were conducted using complete cases.
Summary measures of sedentary behavior, e.g., daily sedentary minutes, are linked to negative health in older adults. By leveraging minute-level accelerometry, we derived clusters of diurnal timing of sedentary behavior to further elucidate associations with health. We used a two-phase clustering approach that first elicited daily (hour-by-hour) trajectories of SB time accumulation, and then clustered individuals based on their daily patterns across multiple days. This approach allowed us to group participants according to the amount of daily SB, as well as the timing of their SB accumulation. E.g., we identified two clusters that had similar (moderate) total daily SB yet differed in the phase of the day (morning vs evening) wherein most sedentary time was accumulated. Such delineation is not possible if only total daily SB is used to characterize an individual’s SB.
We next investigated whether the derived diurnal SB patterns were associated with cross-sectional and longitudinal markers of health. Unsurprisingly, we found that higher sedentary behavior throughout the day was generally associated with poorer health indicators (high BMI, lower physical functioning, more morbidities) compared to moderate or lower sedentary behavior. However, a novel and intriguing finding was that among participants who had below median MVPA with moderate daily sedentary behavior, those who were less sedentary earlier in their day (vs later) had significantly higher physical functioning scores at study entry but also steeper decline in longitudinal physical health compared to those who had high levels of sedentary behavior throughout the day. Thus, we seem to have identified a latent subgroup that is at high risk of physical functioning decline, a finding that would have been masked had we only used total daily SB as the exposure of interest. Interestingly, these effects were not as apparent in women who had above-median levels of MVPA. In particular, although patterns 1 and 2 exhibited higher rates of decline than patterns 3 and 4, the effect-sizes were smaller, suggesting that MVPA may partially mitigate the negative effects of diurnal sedentary behavior patterns. This finding is similar to recently-published results showing that self-reported sitting time is more strongly related to all-cause mortality among adults with low levels of MVPA compared to adults with higher levels of MVPA [44, 45].
While to our knowledge this is the first study to evaluate diurnal patterns of SB in relation to physical function, Schrack et al. observed trends in physical activity accumulation across the day, showing that morning PA was on average similar across age groups for the 611 adults (aged 32–93) they studied, but that older adults had significantly less activity as the day progressed than younger adults . They hypothesized that, “Older adults tended to reach their peak activity level much earlier in the day, which may indicate that daily tasks such as bathing, dressing, running errands, and volitional physical exercise are performed in the morning hours with little activity later in the day.” Our data expand on this work by showing that even among a population consisting solely of older adults and after accounting for their chronological age, there is noticeable diurnal variation in accumulation of SB. Furthermore, compared to consistently high SB throughout the day (i.e., cluster p1), we found that a pattern with lower SB earlier and higher SB later in the day (i.e., cluster p2) is related to accelerated functional decline, among participants who engaged in low MVPA. The extent to which this novel observation represents a causal relationship that evening accumulation of SB is more deleterious to physical function than morning SB, requires further study. It is possible that there is an underlying biological ageing process that affects both diurnal SB patterns and trajectories of physical function. Or, sedentary behavior later in the day may disrupt or change circadian rhythms of biological processes linked to ageing and chronic disease. It is important to note that the high-risk cluster (p2) in fact generally had better health (i.e., fewer morbidities, lower BMI) at baseline relative to the high-SB (p1) group. If confirmed in other studies, knowing that specific diurnal activity patterns are associated with increased risk for accelerated decline might identify high-risk patients, who might not otherwise be considered high risk, and alert clinicians early to these patients. In the future, this might trigger a suite of interventions to prevent such decline.
In our older adult population, multiple health conditions are prevalent, which could confound the observed associations. We attempted to adjust for these health conditions by including a “number of morbidities” covariate. However, it is possible that specific morbidities may exert stronger confounding effects. To this end, we conducted a sensitivity analysis in which we controlled for specific conditions, such as stroke, diabetes, and arthritis, which would be strongly related to SB and physical function. The results of this sensitivity analysis were similar to our original findings (Additional File 2). Another important confounding factor is vision loss and the likely higher fear of falling among participants with vision loss . We emphasize that while our analysis adjusted for a host of measured confounders, we cannot rule out residual confounding and acknowledge this as a limitation.
We used physical functioning trajectories across 6 years as a measure of overall functioning and health in our population of older women. The validated RAND-36 physical functioning subscale is reported to be a comprehensive marker of mobility-disability among older adults . It is recognized that aging-related decline in health is a multidimensional phenomenon [49, 50], and the physical functioning score used herein may be regarded as an overarching health indicator that captures multiple facets and morbidities of health decline in older adults. A natural next step would be to identfy associations between diurnal SB patterns and specific biological and self-reported outcomes that may underlie the current finding. We leave this to future research.
Many investigators have developed new analytic approaches for accelerometer data in an effort to not just measure total accumulation of PA/SB time but to detect patterns in the data. Several studies have investigated patterns based on lengths of uninterrupted intervals (“bouts”) and/or degree of fragmentation of PA/SB [15, 22, 51,52,53,54,55,56,57]. More relevant to the current investigation, studies have also described variations in the timing of sedentary time and PA time accumulation with respect to hour of the day and/or day of the week, including differences in timing by sociodemographic factors [23, 28, 58,59,60,61,62]. None of these studies about the timing of PA/SB attempted to determine whether the identified patterns were associated with health status. We identified a few studies that examined diurnal patterns of PA/SB, (i.e., patterns based on the time of day that active/sedentary time was accumulated) and health. Using functional principal components analysis, Xu et al. reported that higher evening (vs mid-day) activity was associated with worse mental quality-of-life . Zeitzer et al. found principal components-derived patterns that were predictive of changes in sleep and cognition, as well as cardiovascular-related mortality and all-cause mortality . However, both studies modeled accelerometer/ActiGraph counts as the exposure measure of activity and did not examine SB directly. Also, the study design used by Xu et al. was cross-sectional, and they did not investigate longitudinal trajectories of quality-of-life. In another study, using k-means analysis of METs data (metabolic equivalent of tasks), Fukuoka et al. discriminated three groups of participants according to the times of day at which their peak levels of MVPA occurred . The MVPA evening peak group had significantly higher BMI, waist circumference, and hip circumference than the MVPA noon peak group. An x-means cluster analysis study of diurnal patterns, similar to our own, identified 4 distinct PA clusters on the basis of intensity and temporal patterns of activity, and found that risk of cardiovascular disease varied across clusters in both males and females . The functional principal components analysis approach discerns the major modes of variation in PA/SB patterns, while our and other clustering methods identify distinct homogeneous subgroups of participants. Fukuoka et al. and Niemala et al. derived clusters using diurnal patterns of overall activity measured via METs, whereas our approach focuses specifically on SB. For behavioral intervention studies, eliciting timing of specific behaviors (such as SB or MVPA) may be most useful for tailoring interventions to the needs of the individual. For instance, given findings from our study, a behavioral intervention aiming to reduce SB would offer different strategies to individuals who have one versus the other of the two moderate SB patterns.
Our study has many strengths. Our sample comprised a large diverse and well-characterized cohort of > 6000 older participants from the Women’s Health Initiative. The rich Women’s Health Initiative database contains detailed data on a large number of variables, allowing us to develop robust models and adjust for a number of potential confounders. We used objective accelerometers to measure SB and PA, and leveraged minute-level assessments to ascertain diurnal variation in SB patterns. Our 2-stage clustering approach uncovered latent subgroups based on within-day and between-day similarity in SB patterns, which could be useful for developing targeted interventions to reduce SB, e.g., when during the day to intervene (morning vs afternoon) and which day (e.g., a high SB-day). Due to the availability of multiple longitudinal measures of physical functioning, we were able to examine how SB pattern impacts not only current physical functioning but also future rate of decline. There are also limitations. Clustering analysis is an unsupervised, inherently exploratory technique. While we conducted a variety of sensitivity analyses to assess the stability of our results, additional independent studies are needed to replicate (or refute) our findings. Also, the derived clusters in our study provide broad groupings of similar SB patterns. Although finer-grained groupings are possible, we elected to maintain a coarser classification, so as not to overfit the data at hand. If our findings are replicated, further subdivisions of clusters may offer additional insights. In our study, we treated weekday and weekend days equally. While studies among similarly aged older adults indicate that total sedentary time does not significantly vary by weekday vs weekend day [65, 66], there is evidence that diurnal accumulation patterns may differ and future studies should further explore this to guide development of personalized interventions. Our physical functioning measure, albeit validated and used extensively in “aging” research, is self-reported. Replicating our results using longitudinal objective outcomes of physical health would strengthen our findings. Even though we were able to adjust for a large number of covariates, unmeasured confounding could still be present (e.g., by variables such as occupation that were not measured in OPACH), and must be acknowledged when interpreting our results. E.g., the finding that the moderate-SB-with-lower-morning-SB pattern  had steeper physical functioning decline could signify a latent subgroup at high risk of decline, but could also be due to an unmeasured baseline factor. Similarly, in our observational study it is not possible to infer the direction of associations between MVPA (or other baseline factors) and baseline physical functioning. Even so, availability of longitudinal measures of physical functioning allows us to infer effects of baseline MVPA on future physical functioning decline, thus partially addressing this concern. Further investigation of the sequence of cause and effects is warranted.
In summary, in a large cohort of older, ambulatory, community-living women, we found that diurnal SB patterns were associated with differential rates of physical functioning decline. Diurnal timing of SB, complementary to total sedentary time and MVPA, may offer additional insights into associations between SB and physical health.
Availability of data and materials
The data that support the findings of this study are available from Dr. Andrea LaCroix on reasonable request in accordance with the WHI (Women’s Health Initiative) publications and presentations policy.
Body mass index
General educational development
Metabolic equivalent of task
Moderate-to-vigorous physical activity
Objective Physical Activity and Cardiovascular Health
United Nations Department of Economic and Social Affairs/Population Division. World Population Prospects: The 2017 Revision, Key Findings and Advance Tables [Internet]. Available from: https://esa.un.org/unpd/wpp/Publications/Files/WPP2017_KeyFindings.pdf. Accessed 1 July 2020.
Crimmins EM, Beltrán-Sánchez H. Mortality and morbidity trends: Is there compression of morbidity? J Gerontol - Ser B Psychol Sci Soc Sci. 2011;66 B(1):75–86.
Harris-Kojetin L, Sengupta M, Lendon JP, Rome V, Valverde R, Caffrey C. Long-Term Care Providers and Services Users in the United States: Data From the National Study of Long-Term Care Providers, 2015–2016. Vital Heal Stat. 2019;3(43):1–78..
Rice DP. Older women’s health and access to care. Women’s Heal Issues. 2000;10(2):42–6.
Spillman BC, Lubitz J. The effect of longevity on spending for acute and long-term care. N Engl J Med. 2000;342(19):1409–15.
Shiroma EJ, Freedson PS, Trost SG, Lee I-M. Patterns of accelerometer-assessed sedentary behavior in older women. J Am Med Assoc. 2013;310(23):2562–3. https://doi.org/10.1001/jama.2013.278896.
Evenson KR, Buchner DM, Morland KB. Objective measurement of physical activity and sedentary behavior among US adults aged 60 years or older. Prev Chronic Dis. 2012;9(1):2–11..
LaCroix AZ, Rillamas-Sun E, Buchner D, Evenson KR, Di C, Lee I-MI-M, et al. The Objective Physical Activity and Cardiovascular Disease Health in Older Women (OPACH) Study. BMC Public Health. 2017;17(192) https://doi.org/10.1186/s12889-017-4065-6.
Matthews CE, Chen KY, Freedson PS, Buchowski MS, Beech BM, Pate RR, et al. Amount of time spent in sedentary behaviors in the United States, 2003-2004. Am J Epidemiol. 2008;167(7):875–81.
Sparling PB, Howard BJ, Dunstan DW, Owen N. Recommendations for physical activity in older adults. BMJ. 2015;350:1–5.
Colley RC, Garriguet D, Janssen I, Craig CL, Clarke J, Tremblay MS. Physical activity of Canadian children and youth: Accelerometer results from the 2007 to 2009 Canadian health measures survey. Heal Rep. 2011;22(1):10.
Colley RC, Garriguet D, Janssen I, Craig CL, Clarke J, Tremblay MS. Physical activity of Canadian adults: Accelerometer results from the 2007 to 2009 Canadian health measures survey. Heal Rep. 2011;22(1):9.
Diaz KM, Howard VJ, Hutto B, Colabianchi N, Vena JE, Blair SN, et al. Patterns of sedentary behavior in US middle-age and older adults: the REGARDS study. Med Sci Sports Exerc. 2016;48(3):430–8.
Patterson R, McNamara E, Tainio M, de Sá TH, Smith AD, Sharp SJ, et al. Sedentary behaviour and risk of all-cause, cardiovascular and cancer mortality, and incident type 2 diabetes: a systematic review and dose response meta-analysis. Eur J Epidemiol. 2018;33(9):811–29.
Bellettiere J, LaMonte MJ, Evenson KR, Rillamas-Sun E, Kerr J, Lee I-M, et al. Sedentary behavior and cardiovascular disease in older women: the OPACH study. Circulation. 2019;139(8):1036–46. https://doi.org/10.1161/CIRCULATIONAHA.118.035312.
de Rezende LFM, Rey-López JP, VKR M, do Luiz O C. Sedentary behavior and health outcomes among older adults: a systematic review. BMC Public Health. 2014;14(333):1–9.
Van Uffelen JGZ, Heesch KC, Hill RL, Brown WJ. A qualitative study of older adults’ responses to sitting-time questions: do we get the information we want? BMC Public Health. 2011;11:1–11.
Jefferis BJ, Sartini C, Ash S, Lennon LT, Wannamethee SG, Whincup PH. Validity of questionnaire-based assessment of sedentary behaviour and physical activity in a population-based cohort of older men; comparisons with objectively measured physical activity data. Int J Behav Nutr Phys Act. 2016;13(1):14 Available from: http://www.scopus.com/inward/record.url?eid=2-s2.0-84959368086&partnerID=tZOtx3y1.
LaMonte MJ, Lee I-M, Rillamas-Sun E, Bellettiere J, Evenson KR, Buchner DM, et al. Comparison of questionnaire and device measures of physical activity and sedentary behavior in a multi-ethnic cohort of older women. J Meas Phys Behav. 2019;2(2):82–93.
Celis-Morales CA, Perez-Bravo F, Ibañez L, Salas C, Bailey MES, Gill JMR. Objective vs. self-reported physical activity and sedentary time: Effects of measurement method on relationships with risk biomarkers. PLoS One. 2012;7(5):1–5.
Ekelund U, Tarp J, Steene-Johannessen J, Hansen BH, Jefferis B, Fagerland MW, et al. Dose-response associations between accelerometry measured physical activity and sedentary time and all cause mortality: systematic review and harmonised meta-analysis. Bmj. 2019;366:1–10.
Chang YJ, Bellettiere J, Godbole S, Keshavarz S, Maestas JP, Unkart JT, et al. Total sitting time and sitting pattern in postmenopausal women differ by Hispanic ethnicity and are associated with Cardiometabolic risk biomarkers. J Am Heart Assoc. 2020;9(4):e013403.
Bellettiere J, Carlson JA, Rosenberg D, Singhania A, Natarajan L, Berardi V, et al. Gender and age differences in hourly and daily patterns of sedentary time in older adults living in retirement communities. PLoS One. 2015;10(8):1–15.
Wanigatunga AA, Simonsick EM, Zipunnikov V, Spira AP, Studenski S, Ferrucci L, et al. Perceived fatigability and objective physical activity in mid- to late-life. J Gerontol - Ser A Biol Sci Med Sci. 2018;73(5):630–5.
Varma VR, Watts A. Daily physical activity patterns during the early stage of Alzheimer’s disease. J Alzheimers Dis. 2017;55(2):659–67.
Xu SY, Nelson S, Kerr J, Godbole S, Johnson E, Patterson RE, et al. Modeling temporal variation in physical activity using functional principal components analysis. Stat Biosci. 2019;11(2):403–21.
Zeitzer JM, Blackwell T, Hoffman AR, Cummings S, Ancoli-Israel S, Stone K. Daily patterns of accelerometer activity predict changes in sleep, cognition, and mortality in older men. J Gerontol - Ser A Biol Sci Med Sci. 2018;73(5):682–7.
Van Cauwenberg J, Van Holle V, De Bourdeaudhuij I, Owen N, Deforche B. Diurnal patterns and correlates of older adults’ sedentary behavior. PLoS One. 2015;10(8):1–15.
Rillamas-Sun E, Buchner DM, Di C, Evenson KR, AZ LC. Development and application of an automated algorithm to identify a window of consecutive days of accelerometer wear for large-scale studies. BMC Res Notes. 2015;8(1):1–8 Available from: http://www.biomedcentral.com/1756-0500/8/270. [cited 2020 Feb 24].
Ware JE. SF-36 health survey update. Spine (Phila Pa 1976). 2000;25(24):3130–9.
Ware CJE, Sherbourn D. The MOS 36-item short-form health survey (SF-36): I. conceptual framework and item selection. Med Care. 1992;30(6):473–83.
Choi L, Liu Z, Matthews CE, Buchowski MS. Validation of accelerometer wear and nonwear time classification algorithm. Med Sci Sports Exerc. 2011;43(2):357–64.
Migueles JH, Cadenas-Sanchez C, Ekelund U, Nyström CD, Mora-Gonzalez J, Löf M, et al. Accelerometer data collection and processing criteria to assess physical activity and other outcomes: a systematic review and practical considerations. Sport Med. 2017;47(9):1821–45.
Evenson KR, Wen F, Herring AH, Di C, LaMonte MJ, Tinker LF, et al. Calibrating physical activity intensity for hip-worn accelerometry in women age 60 to 91years: The Women’s Health Initiative OPACH Calibration Study. Prev Med Reports. 2015;2:750–6. https://doi.org/10.1016/j.pmedr.2015.08.021.
Rillamas-Sun E, LaCroix AZ, Bell CL, Ryckman K, Ockene JK, Wallace RB. The impact of multimorbidity and coronary disease comorbidity on physical function in women aged 80 years and older: the Women’s health initiative. J Gerontol Ser A Biol Sci Med Sci. 2016;71(Suppl 1):S54–61.
Genolini C, Écochard R, Jacqmin-Gadda H. Copy mean: a new method to impute intermittent missing values in longitudinal studies. Open J Stat. 2013;03(04):26–40.
Calinski T, Harabasz J. Communications in Statistics a dendrite method for cluster analysis. Commun Stat. 1974;3(1):1–27.
James G, Witten D, Hastie T, Tibshirani R. An introduction to statistical learning. New york: Springer; 2013.
Maechler M, Rousseeuw P, Struyf A, Hubert M, Hornik K. cluster: Cluster Analysis Basics and Extensions. R package version 2.1.0. 2019. Available from: https://cran.r-project.org/web/packages/cluster/citation.html.
LaCroix AZ, Bellettiere J, Rillamas-Sun E, Di C, Evenson KR, Lewis CE, et al. Association of light physical activity measured by accelerometry and the incidence of coronary heart disease and cardiovascular disease in older women. JAMA Open. 2019;2(3):e190419.
van Buuren S, Groothuis-Oudshoorn K. mice: Multivariate Imputation by Chained Equations in R. J Stat Software, Artic. 2011;45(31):1–67 Available from: https://www.jstatsoft.org/v045/i03.
Genolini C, Alacoque X, Sentenac M, Arnaud C. Kml and kml3d: R packages to cluster longitudinal data. J Stat Softw. 2015;65(4):1–34.
Pinheiro J, Bates D, DebRoy S, Sarkar D, R Core Team. nlme: Linear and Nonlinear Mixed Effects Models. R package version 3.1–141. 2019. Available from: https://cran.r-project.org/web/packages/nlme/citation.html.
Stamatakis E, Gale J, Bauman A, Ekelund U, Hamer M, Ding D. Sitting time, physical activity, and risk of mortality in adults. J Am Coll Cardiol. 2019;73(16):2062–72.
Ekelund U, Steene-Johannessen J, Brown WJ, Fagerland MW, Owen N, Powell KE, et al. Does physical activity attenuate, or even eliminate, the detrimental association of sitting time with mortality? A harmonised meta-analysis of data from more than 1 million men and women. Lancet [Internet]. 2016;388(10051):1302–10. https://doi.org/10.1016/S0140-6736(16)30370-1.
Schrack JA, Zipunnikov V, Goldsmith J, Bai J, Simonsick EM, Crainiceanu C, et al. Assessing the physical cliff: detailed quantification of age-related differences in daily patterns of physical activity. J Gerontol - Ser A Biol Sci Med Sci. 2014;69(8):973–9.
Jian-Yu E, Mihailovic A, Kuo PL, West SK, Friedman DS, Gitlin LN, et al. Characterizing the impact of fear of falling on activity and falls in older adults with glaucoma. J Am Geriatr Soc. 2020:1–5.
Syddall HE, Martin HJ, Harwood RH, Cooper C, Sayer AA. The SF-36: a simple, effective measure of mobility-disability for epidemiological studies. J Nutr Heal Aging. 2009;13(1):57–62.
Calderón-Larrañaga A, Vetrano DL, Ferrucci L, Mercer SW, Marengoni A, Onder G, et al. Multimorbidity and functional impairment–bidirectional interplay, synergistic effects and common pathways. J Intern Med. 2019;285(3):1–17.
Kritchevsky SB, Forman DE, Callahan KE, Ely EW, High KP, McFarland F, et al. Pathways, contributors, and correlates of functional limitation across specialties: workshop summary. J Gerontol - Ser A Biol Sci Med Sci. 2019;74(4):534–43.
Davis M, Fox K, Stathi A, Trayers T, Thompson J, Cooper A. Objectively measured sedentary time and its association with physical function in older adults. J Aging Phys Act. 2015;22(4):474–81.
Sardinha LB, Ekelund U, dos Santos L, Cyrino ES, Silva AM, Santos DA. Breaking-up sedentary time is associated with impairment in activities of daily living. Exp Gerontol. 2015;72:57–62.
Bellettiere J, Healy GN, LaMonte MJ, Kerr J, Evenson KR, Rillamas-Sun E, et al. Sedentary behavior and prevalent diabetes in 6,166 older women: the objective physical activity and cardiovascular health stmudy. Journals Gerontol - Ser a biol Sci med Sci. 2019;74(3):387–95 Available from: http://www.ncbi.nlm.nih.gov/pubmed/29726906.
Diaz KM, Howard VJ, Hutto B, Colabianchi N, Vena JE, Safford MM, et al. Patterns of sedentary behavior and mortality in U.S. middle-aged and older adults a national cohort study. Ann Intern Med. 2017;167(7):465–75.
Schrack JA, Kuo PL, Wanigatunga AA, Di J, Simonsick EM, Spira AP, et al. Active-to-sedentary behavior transitions, fatigability, and physical functioning in older adults. J Gerontol - Ser A Biol Sci Med Sci. 2019;74(4):560–7.
Wanigatunga AA, Di J, Zipunnikov V, Urbanek JK, Kuo P, Simonsick EM. Association of Total Daily Physical Activity and Fragmented Physical Activity with Mortality in older adults. JAMA Netw Open. 2019;2(10):1–11.
Hartman SJ, Marinac CR, Cadmus-Bertram L, Kerr J, Natarajan L, Godbole S, et al. Sedentary behaviors and biomarkers among breast cancer survivors. J Phys Act Health. 2018;15(1):1–6.
Chinapaw MJ, Wang X, Andersen LB, Altenburg TM. From Total volume to sequence maps: sophisticated accelerometer data analysis. Med Sci Sports Exerc. 2019;51(4):814–20.
Yerrakalva D, Cooper AJ, Westgate K, Khaw KT, Wareham NJ, Brage S, et al. The descriptive epidemiology of the diurnal profile of bouts and breaks in sedentary time in older English adults. Int J Epidemiol. 2017;46(6):1871–81.
Sartini C, Wannamethee SG, Iliffe S, Morris RW, Ash S, Lennon L, et al. Diurnal patterns of objectively measured physical activity and sedentary behaviour in older men. BMC Public Health. 2015;15(1):1–13.
Kurita S, Shibata A, Ishii K, Koohsari MJ, Owen N, Oka K. Patterns of objectively assessed sedentary time and physical activity among Japanese workers: a cross-sectional observational study. BMJ Open. 2019;9(2):1–9.
McVeigh JA, Winkler EAH, Howie EK, Tremblay MS, Smith A, Abbott RA, et al. Objectively measured patterns of sedentary time and physical activity in young adults of the Raine study cohort. Int J Behav Nutr Phys Act [Internet]. 2016;13(1):1–12. https://doi.org/10.1186/s12966-016-0363-0.
Fukuoka Y, Zhou M, Vittinghoff E, Haskell W, Goldberg K, Aswani A. Objectively measured baseline physical activity patterns in women in the mped trial: cluster analysis. J Med Internet Res. 2018;20(2):1–12.
Niemelä M, Kangas M, Farrahi V, Kiviniemi A, Leinonen AM, Ahola R, et al. Intensity and temporal patterns of physical activity and cardiovascular disease risk in midlife. Prev med (Baltim) [internet]. 2019;124(April):33–41. Available from. https://doi.org/10.1016/j.ypmed.2019.04.023.
Marshall S, Kerr J, Carlson J, Cadmus-Bertram L, Patterson R, Wasilenko K, et al. Patterns of weekday and weekend sedentary behavior among older adults. J Aging Phys Act. 2015;23(4):534–41.
Schlaff RA, Baruth M, Boggs A, Hutto B. Patterns of sedentary behavior in older adults. Am J Health Behav. 2017;41(4):411–8.
We thank the Women’s Health Initiative participants, staff, and investigators. The full list of Women’s Health Initiative Investigators can be found at the following site: www.whi.org/researchers/Documents%20%20Write%20a%20Paper/WHI%20Investigator%20Long%20List.pdf
The National Heart, Lung, and Blood Institute provided funding for the OPACH study (grant number RO1 HL105065 to AZL). Funding also came from a training grant provided by the National Institutes of Health (grant numbers T32HL079891 to JB). The Women’s Health Initiative program was funded by the National Heart, Lung, and Blood Institute, National Institutes of Health, US Department of Health and Human Services (contract numbers HHSN268201600018C, HHSN268201600001C, HHSN268201600002C, HHSN268201600003C, and HHSN268201600004C). The current analysis was partially funded by a program project grant supported by the National Institute on Aging (1P01AG052352-01A1 PI AZL). All funding agencies had no role in the design of the study and collection, analysis, or interpretation of data and and no role in writing the manuscript.
Ethics approval and consent to participate
The study protocol was approved by the ethics review board at the Fred Hutchinson Cancer Research Center, approval number 140654XL, and all participants gave informed consent. Women’s Health Initiative women in the Long Life Study consented to periodic in-home examinations to provide blood samples, updated health information, and physical measurements.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Reuter, C., Bellettiere, J., Liles, S. et al. Diurnal patterns of sedentary behavior and changes in physical function over time among older women: a prospective cohort study. Int J Behav Nutr Phys Act 17, 88 (2020). https://doi.org/10.1186/s12966-020-00992-x